Biomedical subjects

Tohru Terada

Publications and source records attributed to Tohru Terada.

10 recordsLinked to original sources

PubMed2006

Mechanism of the difference in the binding affinity of E. coli tRNAGln to glutaminyl-tRNA synthetase caused by noninterface nucleotides in variable loop.

Aminoacyl-tRNA synthetases (ARSs) distinguish their cognate tRNAs from many other kinds of tRNAs, despite the very similar tertiary structures of tRNAs. Many researchers have supported the view that this recognition is achieved by intermolecular interactions between tRNA and ARS. However, one of the aptamers of Escherichia coli glutamine specific tRNA, var-AGGU, has a higher affinity to ARS than the wild-type, although the sequence difference only lies in the variable loop located on the opposite side of the binding interface with ARS. To understand the reason for the difference in affinity, we did molecular dynamics simulations on tRNAs and their complexes with ARS. We calculated the enthalpic and entropic contributions to the binding free energy with the molecular mechanics-Poisson-Boltzmann/surface area method and found that the entropic difference plays an important role in the difference in binding free energies. During the molecular dynamics simulations, dynamic rearrangements of hydrogen bonds occurred in the tertiary core region of the wild-type tRNA, whereas they were not observed in the free var-AGGU simulation. Since the internal mobility was suppressed upon complex formation with ARS, the entropy loss in the wild-type was larger than that of the aptamer. We therefore concluded that the sequence difference in the variable loop caused the difference in the internal mobility of the tertiary core region tRNAs and led to the difference in the affinity to ARS through the entropy term.

Amino Acyl-tRNA Synthetases↗

PubMed2006

ROKU: a novel method for identification of tissue-specific genes.

BACKGROUND: One of the important goals of microarray research is the identification of genes whose expression is considerably higher or lower in some tissues than in others. We would like to have ways of identifying such tissue-specific genes. RESULTS: We describe a method, ROKU, which selects tissue-specific patterns from gene expression data for many tissues and thousands of genes. ROKU ranks genes according to their overall tissue specificity using Shannon entropy and detects tissues specific to each gene if any exist using an outlier detection method. We evaluated the capacity for the detection of various specific expression patterns using synthetic and real data. We observed that ROKU was superior to a conventional entropy-based method in its ability to rank genes according to overall tissue specificity and to detect genes whose expression pattern are specific only to objective tissues. CONCLUSION: ROKU is useful for the detection of various tissue-specific expression patterns. The framework is also directly applicable to the selection of diagnostic markers for molecular classification of multiple classes.

Algorithms↗

PubMed2006

Evaluation of efficiency of reference potential spatial warping algorithm in conformational sampling of peptides.

Molecular dynamics (MD) simulation using the reference potential spatial warping algorithm (REPSWA) is a promising method for computing the conformational equilibrium of a system with a rugged energy surface. Its effectiveness has previously been demonstrated using only a simplified model system of a hydrocarbon chain omitting nonbonded interaction terms from the potential energy function. To evaluate the applicability of REPSWA MD simulation to more realistic problems, we applied it to small peptide systems in an aqueous environment. Difficulties were encountered, however, forcing us to devise several modifications. We investigated their effectiveness in comparison to conventional constant-temperature and multicanonical MD simulations. We found that the sampling efficiency of the modified REPSWA MD after a careful optimization of its parameters was better than that of the constant-temperature MD and comparable to that of the multicanonical MD in several cases.

Alanine↗

PubMed2006

Folding free-energy landscape of a 10-residue mini-protein, chignolin.

Chignolin is an artificial mini-protein composed of 10 residues (GYDPETGTWG) that has been shown to cooperatively fold into a beta-hairpin structure in water. We extensively explored the conformational space of chignolin using a 180-ns multicanonical molecular dynamics (MD) simulation and analyzed its folding free-energy landscape. In the MD trajectory, we found structures that satisfy 99% of the experimental restraints and are quite close to the experimentally determined structures with C(alpha) root-mean-square-deviations of less than 0.5 Angstroms. These structures formed a large cluster in the conformational space with the largest probability of existence, agreeing well with the experiment.

Cluster Analysis↗

PubMed2006

Crystal structure of neoculin: insights into its sweetness and taste-modifying activity.

Although the majority of sweet compounds are of low molecular mass, several proteins are known to elicit sweet taste responses in humans. The fruit of Curculigo latifolia contains a heterodimeric protein, neoculin, which has both sweetness and a taste-modifying activity that converts sourness to sweetness. Here, we report the crystal structure of neoculin at 2.76A resolution. This is the first well-defined tertiary structure of a taste-modifying protein of this kind. The overall structure is quite similar to those of monocot mannose-binding lectins. However, crucial topological differences are observed in the C-terminal regions of both subunits. In both subunits of neoculin, the C-terminal tails turn up to form loops fixed by inter-subunit disulfide bonds that are not observed in the lectins. Indeed, the corresponding regions of the lectins stretch straight over the surface of another subunit. Such a C-terminal structural feature as is observed in neoculin results in a decrease in subunit-subunit interactions. Moreover, distribution of electrostatic potential on the surface of neoculin is unique and significantly different from those of the lectins, particularly in the basic subunit (NBS). We have found that there is a large cluster composed of six basic residues on the surface of NBS, and speculate that it might be involved in the elicitation of sweetness and/or taste-modifying activity of neoculin. Molecular dynamics simulation based on the crystallography results suggests that neoculin may adopt a widely "open" conformation at acidic pH, while unprotonated neoculin at neutral pH is in a "closed" conformation. Based on these simulations and the generation of a docking model between neoculin and the sweet-taste receptor, T1R2-T1R3, we propose the hypothesis that neoculin is in dynamic equilibrium between open and closed states, and that the addition of an acid shifts the equilibrium to the open state, allowing ligand-receptor interaction.

Amino Acid Sequence↗

PubMed2005

Structure of the terminal oxygenase component of angular dioxygenase, carbazole 1,9a-dioxygenase.

Carbazole 1,9a-dioxygenase (CARDO) catalyzes the dihydroxylation of carbazole by angular position (C9a) carbon bonding to the imino nitrogen and its adjacent C1 carbon. This reaction is an initial degradation reaction of the carbazole degradation pathway by various bacterial strains. Only a limited number of Rieske non-heme iron oxygenase systems (ROSs) can catalyze this novel reaction, termed angular dioxygenation. Angular dioxygenation is also involved in the degradation pathways of carbazole-related compounds, dioxin, and CARDO can catalyze the angular dioxygenation for dioxin. CARDO consists of a terminal oxygenase component (CARDO-O), and the electron transport components, ferredoxin (CARDO-F) and ferredoxin reductase (CARDO-R). CARDO-O has a homotrimeric structure, and governs the substrate specificity of CARDO. Here, we have determined the crystal structure of CARDO-O of Janthinobacterium sp. strain J3 at a resolution of 1.95A. The alpha3 trimeric overall structure of the CARDO-O molecule roughly corresponds to the alpha3 partial structures of other terminal oxygenase components of ROSs that have the alpha3beta3 configuration. The CARDO-O structure is a first example of the terminal oxygenase components of ROSs that have the alpha3 configuration, and revealed the presence of the specific loops that interact with a neighboring subunit, which is proposed to be indispensable for stable alpha3 interactions without structural beta subunits. The shape of the substrate-binding pocket of CARDO-O is markedly different from those of other oxygenase components involved in naphthalene and biphenyl degradation pathways. Docking simulations suggested that carbazole binds to the substrate-binding pocket in a manner suitable for catalysis of angular dioxygenation.

Bacterial Proteins↗

PubMed2005

Characterization of the interaction between interleukin-13 and interleukin-13 receptors.

Interleukin-13 (IL-13) possesses two types of receptor: the heterodimer, composed of the IL-13Ralpha1 chain (IL-13Ralpha1) and the IL-4Ralpha chain (IL-4Ralpha), transducing the IL-13 signals; and the IL-13Ralpha2 chain (IL-13Ralpha2), acting as a nonsignaling "decoy" receptor. Extracellular portions of both IL-13Ralpha1 and IL-13Ralpha2 are composed of three fibronectin type III domains, D1, D2, and D3, of which the last two comprise the cytokine receptor homology modules (CRHs), a common structure of the class I cytokine receptor superfamily. Thus far, there has been no information about the critical amino acids of the CRHs or the role of the D1 domains of IL-13Ralpha1 and IL-13Ralpha2 in binding to IL-13. In this study, we first built the homology modeling of the IL-13.hIL-13 receptor complexes and then predicted the amino acids involved in binding to IL-13. By incorporating mutations into these amino acids, we identified Tyr-207, Asp-271, Tyr-315, and Asp-318 in the CRH of human IL-13Ralpha2, and Leu-319 and Tyr-321 in the CRH of human IL-13Ralpha1, as critical residues for binding to IL-13. Tyr-315 in IL-13Ralpha2 and Leu-319 in IL-13Ralpha1 are positionally conserved hydrophobic amino acid residues. Furthermore, by using D1 domain-deleted mutants, we found that the D1 domain is needed for the expression of IL-13Ralpha2, but not IL-13Ralpha1, and that the D1 domain of IL-13Ralpha1 is important for binding to IL-13, but not to IL-4. These results provide the basis for a precise understanding of the interaction between IL-13 and its receptors.

Amino Acid Sequence↗

PubMed2005

Crystal structure of the terminal oxygenase component of cumene dioxygenase from Pseudomonas fluorescens IP01.

The crystal structure of the terminal component of the cumene dioxygenase multicomponent enzyme system of Pseudomonas fluorescens IP01 (CumDO) was determined at a resolution of 2.2 A by means of molecular replacement by using the crystal structure of the terminal oxygenase component of naphthalene dioxygenase from Pseudomonas sp. strain NCIB 9816-4 (NphDO). The ligation of the two catalytic centers of CumDO (i.e., the nonheme iron and Rieske [2Fe-2S] centers) and the bridging between them in neighboring catalytic subunits by hydrogen bonds through a single amino acid residue, Asp231, are similar to those of NphDO. An unidentified external ligand, possibly dioxygen, was bound at the active site nonheme iron. The entrance to the active site of CumDO is different from the entrance to the active site of NphDO, as the two loops forming the lid exhibit great deviation. On the basis of the complex structure of NphDO, a biphenyl substrate was modeled in the substrate-binding pocket of CumDO. The residues surrounding the modeled biphenyl molecule include residues that have already been shown to be important for its substrate specificity by a number of engineering studies of biphenyl dioxygenases.

Amino Acid Sequence↗

PubMed2005

A fast protein-protein docking algorithm using series expansion in terms of spherical basis functions.

We describe a fast protein-protein docking algorithm using a series expansion in terms of newly designed bases to efficiently search the entire six-dimensional conformational space of rigid body molecules. This algorithm is an ab initio docking algorithm designed to list candidates of putative conformations from a global conformational space for unbound docking. In our algorithm, a scoring function is constructed from terms that are the inner products of two scalar fields expressing individual molecules. The mapping from a molecule to a scalar field can be arbitrarily defined to express an energy term. Since this scoring scheme has the same expressiveness as that of a method using a fast Fourier transform (FFT), it has the flexibility to introduce various physicochemical energies. Currently, we are using scalar fields that approximate desolvation free energy and steric hindrance energy. Fast calculation of the scoring function for each conformation of the six-dimensional search space is realized by expansion of the fields in terms of basis functions which are combinations of spherical harmonics and modified Legendre polynomials, and the use of only low-order terms, which carry most of the information on the scalar field. We have implemented this algorithm and evaluated the computation time and precision by using actual protein structure data of complexes and their monomers. This paper presents the results for six unbound cases and in all the cases we obtained at least one conformation close to the native structures (interface RMSD < 3.0 A) within the top 1000 candidates with about 40 seconds of computation time using a single Pentium4 2.4 GHz CPU.

Algorithms↗

PubMed2003

Development of an ab initio protein structure prediction system ABLE.

An ab initio protein structure prediction system called ABLE is described. It is based on the fragment assembly method, which consists of two steps: dividing a target sequence into overlapping subsequences (fragments) of short length and assigning a local structure to each fragment; and generating models by assembling the local structures and selecting the models with low potential energy. One of the most important problems in conventional fragment assembly methods is the difficulty of selecting native-like structures by energy minimization only. ABLE thus employs a structural clustering method to select the native-like models from among the generated models. By applying the unit-vector root mean square distance (URMS) as a measure of structure similarity, we achieve more robust, effective structural clustering. When no enough clusters of good quality are obtained, ABLE runs the energy minimization procedure again by incorporating structural restraint conditions obtained from the consensus substructures in the previously generated models. This approach is based on our observation that there is a high probability that the consensus substructures of the generated models have native-like structures. Another feature of ABLE is that in assigning local structures to fragments, it assigns mainchain dihedral angles (phi, psi) to the central residue of each fragment according to a probability distribution map built from candidate sequences similar to each fragment. This enables the system to generate appropriate local structures that may not already exist in a protein structure database. We applied our system to 25 small proteins and obtain near-native folds for more than half of them. We also demonstrate the performance of our structural clustering method, which can be applied to other protein structure prediction systems.

Cluster Analysis↗