Search PubMed⌕ Search

Biomedical subjects

Haruki Nakamura

Publications and source records attributed to Haruki Nakamura.

16 recordsLinked to original sources

Functional and structural analyses of cryptochrome. Vertebrate CRY regions responsible for interaction with the CLOCK:BMAL1 heterodimer and its nuclear localization.

Mouse mCRY1 and zebrafish zCRY1a and zCRY3 belong to the DNA photolyase/Cryptochrome family. mCRY1 and zCRY1a repress CLOCK:BMAL1-mediated transcription, whereas zCRY3 does not. Reciprocal chimeras between zCRY1a and zCRY3 were generated to determine the zCRY1a regions responsible for nuclear translocation, interaction with the CLOCK:BMAL1 heterodimer, and repression of CLOCK:BMAL1-mediated transcription. Three regions, RD-2a-(126-196), RD-1-(197-263), and RD-2b-(264-293), were identified. Proteins in this family consist of an N-terminal alpha/beta domain and a C-terminal helical domain connected by an interdomain loop. RD-2a is within this loop, RD-1 is at the N-terminal 50 amino acids, and RD-2b at the following 31 amino acid residues of the helical domain. Either RD-2a or RD-1 is required for interaction with the CLOCK: BMAL1 heterodimer, and either RD-1 or RD-2b is required for the nuclear translocation of CRY. Both of these functions are prerequisites for the transcriptional repressor activity. The functional nuclear localizing signal in the RD-2b region also was identified. The sequence is well conserved among repressor-type CRYs, including mCRY1. Mutations in the nuclear localizing signal of mCRY1 reduce the extent of its nuclear localization. These findings show that both nuclear localization and interaction with the CLOCK:BMAL heterodimer are essential for transcriptional repression by CRY.

Amino Acid Sequence↗

Mucopolysaccharidosis IVA: characterization of a common mutation found in Finnish patients with attenuated phenotype.

Mucopolysaccharidosis IVA (MPS IVA) is caused by the deficiency of the lysosomal enzyme N-acetylgalactosamine-6-sulfate sulfatase encoded by the GALNS gene on chromosome 16. We describe, in detail, the clinical phenotype of five patients from three unrelated Finnish families and have characterized the disease-causing mutations in GALNS. Genotypes of the patients are D60N/A291T, D60N/W230X, and D60N/1374delT. Mutation 1374delT introduces premature termination of GALNS. Cells over-expressing the novel mutation W230X and A291T had no residual GALNS activity, whereas D60N gave 12.2% residual activity compared with the wild type. Co-transfection of D60N/A291T and D60N/W230X showed 5.5% and 6.7% of wild type activity, respectively. The precursor proteins of D60N and A291T were observed at 55 kDa and 57 kDa, respectively, whereas there was no detectable band in cells over-expressing W230X. At 55 degrees C, the mutant protein showed lower thermostability than the wild type protein at pH 3.8 and 7.0. The tertiary structural model of the GALNS protein revealed that aspartic acid at position 60 is located on the surface of the molecule, away from the active site. This makes it unlikely that the enzymatic function of the protein with D60N is severely impaired. On the other hand, A291 and W230 are localized near the active site. The molecular characteristics of the D60N mutation explain the attenuated clinical phenotype of the patients.

Child↗

Dynamical origin of uniform sampling in multicanonical ensemble.

The stochastic model describing the sampling process in multicanonical ensemble has been derived by considering the sampling process as an overdamped Brownian motion on the free energy surface. The essential dynamics of the multicanonical sampling has been characterized by a Langevin equation in a piecewise multivalleyed free energy landscape, modulated by a temperature-dependent curvature. Based on the stochastic model we showed that the multicanonical weight can be determined by interpolating maximum probability energy points of the canonical samplings at different temperatures.

Journal Article↗

beta-Hairpins, alpha-helices, and the intermediates among the secondary structures in the energy landscape of a peptide from a distal beta-hairpin of SH3 domain.

Energy landscape of a peptide, extracted from a distal beta-hairpin of src SH3 domain, in explicit water was obtained with the multicanonical molecular dynamics. A variety of beta-hairpins with various strand-strand hydrogen bonds were found in the energy landscape at 300 K. There was no energy barrier between random-coil and hairpins. Thus, the peptide conformation can easily change from the random-coil to the hairpins in the thermal fluctuations at 300 K. The landscape also included two clusters of alpha-helices, among which an energy barrier existed, and besides, these helix clusters were separated from the other conformations. Thus, the free-energy barrier exists among the helices and the other conformations. Intermediate clusters were found between the helix and the hairpin clusters. The current study showed that the isolated state of this peptide in water fluctuates among random-coil, beta-hairpin, and alpha-helix. In SH3 domain, which has a topology of mainly beta-protein, the whole-protein folding may proceed when the segment is folded in the beta-hairpin and the other parts of the protein are coupled with the beta-hairpin in an energetically or kinetically favorite way.

Algorithms↗

Protein informatics towards function identification.

The study of structural genomics and structural proteomics has determined the tertiary structures of many hypothetical proteins, whose molecular functions could not be understood using conventional methods. In order to infer the geometrical location of the functional site, the biochemical function and the biological function of the hypothetical protein, much effort has been made in protein informatics. The importance of heterogeneous databases and various descriptors of amino acid sequences, tertiary structures and pathways on the proteome scale has been emphasised.

Binding Sites↗

Structural genomics of membrane proteins.

A program on the structural genomics of membrane proteins has started at the BIRC, AIST, involving other academic institutions and industrial companies. Emphasis is being put on the development of techniques for the structural determination of membrane proteins of biological importance and ligand-receptor interactions by means of electron microscopy, X-ray diffraction, NMR, and computer simulation. Most efforts at the present stage, however, are being directed to finding suitable expression and purification systems and crystallization conditions for such proteins. The program is expected to be linked with the human full-length cDNA project and should lead to medical and industrial uses.

Bacteriorhodopsins↗

Crystal structure of the conserved protein TT1542 from Thermus thermophilus HB8.

The TT1542 protein from Thermus thermophilus HB8 is annotated as a conserved hypothetical protein, and belongs to the DUF158 family in the Pfam database. A BLAST search revealed that homologs of TT1542 are present in a wide range of organisms. The TT1542 homologs in eukaryotes, PIG-L in mammals, and GPI12 in yeast and protozoa, have N-acetylglucosaminylphosphatidylinositol (GlcNAc-PI) de-N-acetylase activity. Although most of the homologs in prokaryotes are hypothetical and have no known function, Rv1082 and Rv1170 from Mycobacterium tuberculosis are enzymes involved in the mycothiol detoxification pathway. Here we report the crystal structure of the TT1542 protein at 2.0 A resolution, which represents the first structure for this superfamily of proteins. The structure of the TT1542 monomer consists of a twisted beta-sheet composed of six parallel beta-strands and one antiparallel beta-strand (with the strand order 3-2-1-4-5-7-6) sandwiched between six alpha-helices. The N-terminal five beta-strands and four alpha-helices form an incomplete Rossmann fold-like structure. The structure shares some similarity to the sugar-processing enzymes with Rossmann fold-like domains, especially those of the GPGTF (glycogen phosphorylase/glycosyl transferase) superfamily, and also to the NAD(P)-binding Rossmann fold domains. TT1542 is a homohexamer in the crystal and in solution, the six monomers forming a cylindrical structure. Putative active sites are suggested by the structure and conserved amino acid residues.

Amino Acid Sequence↗

Identification of protein biochemical functions by similarity search using the molecular surface database eF-site.

The identification of protein biochemical functions based on their three-dimensional structures is strongly required in the post-genome-sequencing era. We have developed a new method to identify and predict protein biochemical functions using the similarity information of molecular surface geometries and electrostatic potentials on the surfaces. Our prediction system consists of a similarity search method based on a clique search algorithm and the molecular surface database eF-site (electrostatic surface of functional-site in proteins). Using this system, functional sites similar to those of phosphoenoylpyruvate carboxy kinase were detected in several mononucleotide-binding proteins, which have different folds. We also applied our method to a hypothetical protein, MJ0226 from Methanococcus jannaschii, and detected the mononucleotide binding site from the similarity to other proteins having different folds.

Algorithms↗

Intramolecular interaction of SUR2 subtypes for intracellular ADP-Induced differential control of K(ATP) channels.

ATP-sensitive K+ (K(ATP)) channels are composed of sulfonylurea receptors (SURs) and inwardly rectifying Kir6.2-channels. The C-terminal 42 amino acid residues (C42) of SURs are responsible for ADP-induced differential activation of K(ATP) channels in SUR-subtypes. By examining ADP-effect on K(ATP) channels containing various chimeras of SUR2A and SUR2B, we identified a segment of 7 residues at central portion of C42 critical for this phenomenon. A 3-D structure model of the region containing the second nucleotide-binding domain (NBD2) of SUR and C42 was developed based on the structure of HisP, a nucleotide-binding protein forming the bacterial Histidine transporter complex. In the model, the polar and charged residues in the critical segment located within a distance that allows their electrostatic interaction with Arg1344 at the Walker-A loop of NBD2. Therefore, the interaction might be involved in the control of ADP-induced differential activation of SUR2-subtype K(ATP) channels.

ATP-Binding Cassette Transporters↗

Tsallis dynamics using the Nosé-Hoover approach.

On the basis of the Nosé-Hoover method, we developed a deterministic algorithm that produces an arbitrary probability density. An ordinary differential equation in the algorithm can realize the Tsallis distribution density. The Tsallis distribution has been considered a candidate of a distribution that represents a physical system in a heat bath. The Tsallis distribution density employed in this algorithm is defined using a full energy function form E(x,p), along with the Tsallis index q > or = 1. Using the current equation, numerical simulations were performed for simple systems and the Tsallis distributions were observed.

Journal Article↗

Calibration of force-field dependency in free energy landscapes of peptide conformations by quantum chemical calculations.

The free energy landscapes of peptide conformations were calibrated by ab initio quantum chemical calculations, after the enhanced conformational diversity search using the multicanonical molecular dynamics simulations. Three different potentials of mean force for an isolated dipeptide were individually obtained by the multicanonical molecular dynamics simulations using the conventional force fields, AMBER parm94, AMBER parm96, and CHARMm22. Each potential of mean force was then calibrated based upon the umbrella sampling algorithm from the adiabatic energy map that was calculated separately by the ab initio molecular orbital method, and all of the calibrated potentials of mean force coincided well. The calibration method was also applied to the simulations of a peptide dimer in explicit water models, and it was shown that the calibrated free energy landscapes did not depend on the force field used in the classical simulations, as far as the conformational space was sampled well. The current calibration method fuses the classical free energy calculation with the quantum chemical calculation, and it should generally make simulations for biomolecular systems much more reliable when combining with enhanced conformational sampling.

Models, Molecular↗

Characterization of six mutations in five Spanish patients with mitochondrial acetoacetyl-CoA thiolase deficiency: effects of amino acid substitutions on tertiary structure.

Mitochondrial acetoacetyl-CoA thiolase (T2) deficiency is an inborn error of ketone body and isoleucine metabolism. We identified and characterized 6 mutations, DelE85, K124R, A127V, Q145E, G152A, and E345V in 5 Spanish T2-deficient patients. Transient expression of mutant cDNAs was done at 37 and at 30 degrees C. Expression of the Q145E mutant cDNA resulted in about 12.5% normal amount at 37 degrees C and it retained 15% residual T2, indicating that specific activity of Q145E mutant protein was almost normal. This mutation reduced the heat stability of T2 activity. Although no significant residual activity was detected in either the G152A and A127V substitution, mutant proteins were detected, at 12.5% the normal amount at 37 degrees C and one-half normal at 30 degrees C for A127V, and 25 % only at 30 degrees C for G152A. Mutant proteins with Q145E, G152A, or A127V accumulated at 30 degrees C expression were stable for 48 h at 37 degrees C after cycloheximide treatment. Expression of DelE85, K124R, and E345V cDNAs gave neither residual T2 protein nor T2 activity. We constructed an improved tertiary structural model of T2 based on the X-ray crystal structure of acetoacetyl-CoA thiolase of Zoogloea ramigera. On the basis of this model, K124, A127, and G152 are located near the active site, mutations of which might affect catalytic function whereas Q145E, De185E, and E345V are distant from the active site with mutants being expected to destabilize the tertiary structure, especially during protein folding and dimerization.

Acetyl-CoA C-Acyltransferase↗

Identification of protein functions from a molecular surface database, eF-site.

A bioinformatics method was developed to identify the protein surface around the functional site and to estimate the biochemical function, using a newly constructed molecular surface database named the eF-site (electrostatic surface of Functional site. Molecular surfaces of protein molecules were computed based on the atom coordinates, and the eF-site database was prepared by adding the physical properties on the constructed molecular surfaces. The electrostatic potential on each molecular surface was individually calculated solving the Poisson-Boltzmann equation numerically for the precise continuum model, and the hydrophobicity information of each residue was also included. The eF-site database is accessed by the internet (http://pi.protein.osaka-u.ac.jp/eF-site/). We have prepared four different databases, eF-site/antibody, eF-site/prosite, eF-site/P-site, and eF-site/ActiveSite, corresponding to the antigen binding sites of antibodies with the same orientations, the molecular surfaces for the individual motifs in PROSITE database, the phosphate binding sites, and the active site surfaces for the representatives of the individual protein family, respectively. An algorithm using the clique detection method as an applied graph theory was developed to search of the eF-site database, so as to recognize and discriminate the characteristic molecular surfaces of the proteins. The method identifies the active site having the similar function to those of the known proteins.

Antibodies↗

Conformational transition states of a beta-hairpin peptide between the ordered and disordered conformations in explicit water.

The conformational transition states of a beta-hairpin peptide in explicit water were identified from the free energy landscapes obtained from the multicanonical ensemble, using an enhanced conformational sampling calculation. The beta-hairpin conformations were significant at 300 K in the landscape, and the typical nuclear Overhauser effect signals were reproduced, consistent with the previously reported experiment. In contrast, the disordered conformations were predominant at higher temperatures. Among the stable conformations at 300 K, there were several free energy barriers, which were not visible in the landscapes formed with the conventional parameters. We identified the transition states around the saddle points along the putative folding and unfolding paths between the beta-hairpin and the disordered conformations in the landscape. The characteristic features of these transition states are the predominant hydrophobic contacts and the several hydrogen bonds among the side-chains, as well as some of the backbone hydrogen bonds. The unfolding simulations at high temperatures, 400 K and 500 K, and their principal component analyses also provided estimates for the transition state conformations, which agreed well with those at 400 K and 500 K deduced from the current free energy landscapes at 400 K and 500 K, respectively. However, the transition states at high temperatures were much more widely distributed on the landscape than those at 300 K, and their conformations were different.

Computer Simulation↗