PubMed Health⌕ Search

Biomedical subjects

Shugo Nakamura

Publications and source records attributed to Shugo Nakamura.

15 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

Potential for assessing quality of protein structure based on contact number prediction.

We developed a novel knowledge-based residue environment potential for assessing the quality of protein structures in protein structure prediction. The potential uses the contact number of residues in a protein structure and the absolute contact number of residues predicted from its amino acid sequence using a new prediction method based on a support vector regression (SVR). The contact number of an amino acid residue in a protein structure is defined by the number of residues around a given residue. First, the contact number of each residue is predicted using SVR from an amino acid sequence of a target protein. Then, the potential of the protein structure is calculated from the probability distribution of the native contact numbers corresponding to the predicted ones. The performance of this potential is compared with other score functions using decoy structures to identify both native structure from other structures and near-native structures from nonnative structures. This potential improves not only the ability to identify native structures from other structures but also the ability to discriminate near-native structures from nonnative structures.

Algorithms↗
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↗
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

Structure, epitope mapping, and docking simulation of a gibberellin mimic peptide as a peptidyl mimotope for a hydrophobic ligand.

Using NMR spectroscopy and simulated annealing calculations, we determined the solution structure of the disulfide-linked cyclized decapeptide ACLPWSDGPC (SD), which is bound to an anti-(gibberellin A(4)) mAb 4-B8(8)/E9 and was found to be the first peptidyl mimotope for a hydrophobic ligand. The resulting structure of the peptide showed a beta-turn-like conformation in residues three to seven and the region converges well (average rmsd 0.54 A). The binding activity and the epitopes of the peptide to the antibody were assessed using saturation transfer difference (STD)-NMR experiments. We also conducted docking simulations between the peptide and the mAb to determine how the peptide is bound to the mAb. Resonances around the beta-turn-like conformation of peptide SD (residues 3-5) showed strong STD enhancement, which agreed well with results from docking simulation between peptide SD and the mAb. Together with the commonality of amino acid residues of the mAb involved in interactions with gibberellin A(4) (GA(4)) and peptide SD, we concluded that peptide SD is bound to the antigen-binding site of mAb 4-B8(8)/E9 as a GA(4) mimic, confirming evidence for the existence of peptide mimics even for hydrophobic ligands.

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↗
PubMed2004

Gene expression profiling of cranial sensory ganglia that transmit food intake stimuli.

Peripheral cranial sensory nerves projecting into the oral cavity receive food intake stimuli and transmit sensory signals to the central nervous system. They are derived from four cranial sensory ganglia, trigeminal, geniculate, petrosal, and nodose ganglia, each of which contains multiple kinds of sensory neurons with different cell morphologies and neuronal properties. We investigated the complex properties of these neurons from the viewpoint of gene expression using DNA microarrays. The 498 genes were selected from a total of 8,740 genes as showing tissue-dependent expression on the microarray by hierarchical cluster analysis, in which several genes known to be differentially expressed in cranial sensory ganglia are included. This suggests that DNA microarray cluster analysis revealed a number of characteristic genes for sensory neurons in these ganglia. Among the selected 498 genes, 44 genes are associated with neurotransmission, such as neuropeptides, their receptors, and vesicle transport, and 26 are ion channels regulating membrane potentials. The identification of a number of genes related directly to neural properties indicates that these sensory ganglia contain heterogeneous types of neurons with different neural properties.

Animals↗
PubMed2003

DNA microarray cluster analysis reveals tissue similarity and potential neuron-specific genes expressed in cranial sensory ganglia.

Each of four cranial sensory ganglia, trigeminal, geniculate, petrosal, and nodose ganglia, contains multiple kinds of sensory neurons with different cell morphologies and neuronal properties that transmit information about sensory stimuli received peripherally. Here we analyze the complex properties of these neurons from the viewpoint of gene expression using DNA microarrays by cluster analysis. From a total of 8,740 genes, 498 genes were selected as showing tissue-dependent expression on the microarray by hierarchical cluster analysis, and their profiles indicated that, among the four sensory ganglia, the petrosal and trigeminal ganglia are intimately related. Tissue trees of 37 subclusters containing the 498 genes showed that the profiles of gene expression and the subclusters were classified into a smaller number of groups (18 groups) when information on the amounts of expression was added. In situ hybridization analysis of 21 genes selected from 13 different groups was carried out, and the gene expression patterns were classified into eight categories. The putative profiles postulated from the microarray data were essentially consistent with the patterns of expression at the cellular level as shown by in situ hybridization. In conclusion, from the overall analyses of gene expression by DNA microarray, we can identify a number of candidate genes showing neuron type-specific expression in the peripheral ganglia.

Animals↗
PubMed2003

Detection of genes with tissue-specific expression patterns using Akaike's information criterion procedure.

We applied a method based on Akaike's information criterion (AIC) to detect genes whose expression profile is considerably different in some tissue(s) than in others. Such observations are detected as outliers, and the method we used was originally developed to detect outliers. The main advantage of the method is that objective decisions are possible because the procedure is independent of a significance level. We applied the method to 48 expression ratios corresponding to various tissues in each of 14,610 clones obtained from the RIKEN Expression Array Database (READ; http://read.gsc.riken.go.jp). As a result, for several tissues (e.g., muscle, heart, and tongue tissues that contain similar cell types) we objectively obtained specific clones without any "thresholding." Our study demonstrates the feasibility of the method for detecting tissue-specific gene expression patterns.

Algorithms↗
PubMed2003

Purification and characterization of meta-cleavage compound hydrolase from a carbazole degrader Pseudomonas resinovorans strain CA10.

2-Hydroxy-6-oxo-6-(2'-aminophenyl)-hexa-2,4dienoic acid [6-(2'-aminophenyl)-HODA] hydrolase, involved in carbazole degradation by Pseudomonas resinovorans strain CA10, was purified to near homogeneity from an overexpressing Escherichia coli strain. The enzyme was dimeric, and its optimum pH was 7.0-7.5. Phylogenetic analysis showed the close relationship of this enzyme to other hydrolases involved in the degradation of monocyclic aromatic compounds, and this enzyme was specific for 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid (6-phenyl-HODA), having little activity toward 2-hydroxy-6-oxohepta-2,4-dienoic acid and 2-hydroxymuconic semialdehyde. The enzyme had a Km of 2.51 microM and k(cat) of 2.14 (s(-1)) for 6-phenyl-HODA (50 mM sodium phosphate, pH 7.5, 25 degrees C). The effect of the presence of an amino group or hydroxyl group at the 2'-position of phenyl moiety of 6-phenyl-HODA on the enzyme activity was found to be small; the activity decreased only in the order of 6-(2'-aminophenyl)-HODA (2.44 U/mg) > 6-phenyl-HODA (1.99 U / mg) > 2-hydroxy-6-oxo-6-(2'-hydroxyphenyl)-hexa-2,4-dienoic acid (1.05 U/mg). The effects of 2'-substitution on the activity were in accordance with the predicted reactivity based on the calculated lowest unoccupied molecular orbital energy for these substrates.

Biotransformation↗
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↗
PubMed2002

Identification of genes regulating colorectal carcinogenesis by using the algorithm for diagnosing malignant state method.

We studied the expression profiles of various stages of colorectal tumors (adenoma (AD), seven samples; carcinoma (CA), 16 samples) by using cDNA microarrays and developed ADMS (algorithm for diagnosing malignant state) method, selecting 335 clones characteristic of CA state. We, then, applied ADMS to 12 additional samples (five from primary lesions with metastasis and seven metastases); all 16 CAs and 12 metastatic tumors were diagnosed correctly as cancerous states. Although three of the seven ADs were diagnosed as "cancerous," the large size of two of these tumors suggested their potential malignancy. Our strategy for selecting clones characteristic of the malignant state is widely applicable to diagnosis and for predicting the stage of progression during multistep carcinogenesis. Of the 335 clones we selected, 135 were known genes. Included in the 135 genes were tumor suppressor and growth factor-related genes and were consistent with the literature. ADMS is a reliable means for identifying genes useful for the diagnosis of cancer.

Adenoma↗
PubMed2002

New method for parallel computation of Hessian matrix of conformational energy function in internal coordinates.

A new algorithm for parallel calculation of the second derivatives (Hessian) of the conformational energy function of biomolecules in internal coordinates is proposed. The basic scheme of this algorithm is the division of the entire calculation of the Hessian matrix (called "task") into subtasks and the optimization of the assignment of processors to each subtask by considering both the load balancing and reduction of the communication cost. A genetic algorithm is used for this optimization considering the dependencies between subtasks. We applied this method to a glutaminyl transfer RNA (Gln-tRNA) molecule for which the scalability of our previously developed parallel algorithm was significantly decreased when the large number of processors was used. The speedup for the calculation was 32.6 times with 60 processors, which is considerably better than the speedup for our previously reported parallel algorithm. The elapsed time for the calculation of subtasks, data sending, and data receiving was analyzed, and the effect of the optimization using the genetic algorithm is discussed.

Algorithms↗