Search

tRNATyr of Nanoarchaeum equitans has a remarkable feature with an extra guanosine residue at the 5ʹ-terminus. However, the N. equitans tRNATyr mutant without extra guanosine at the 5ʹ-end was tyrosylated by tyrosyl-tRNA synthase (TyrRS). We solved the crystal structure of N. equitans TyrRS at 2.80 A resolution. By comparing the present solved structure with the complex structures TyrRS with tRNATyr of Thermus thermophilus and Methanocaldococcus jannaschii, an arginine substitution mutant of N. equitans TyrRS at Ile200 (I200R), which is the putative closest candidate to the 5ʹ-phosphate of C1 of N. equitans tRNATyr, was prepared. The I200R mutant tyrosylated not only wild-type tRNATyr but also the tRNA without the G-1 residue. Further tyrosylation analysis revealed that the second base of the anticodon (U35), discriminator base (A73), and C1:G72 base pair are strong recognition sites.

The nucleosome is the basic organizational unit of the genome. The folding structure of nucleosomes is closely related to genome functions, and has been reported to be in dynamic interplay with binding of various nuclear proteins to genomic loci. Here, we describe our high-throughput chromosome conformation capture with nucleosome orientation (Hi-CO) technology to derive 3D nucleosome positions with their orientations at every genomic locus in the nucleus. This technology consists of an experimental procedure for nucleosome proximity analysis and a computational procedure for 3D modeling. The experimental procedure is based on an improved method of high-throughput chromosome conformation capture (Hi-C) analysis. Whereas conventional Hi-C allows spatial proximity analysis among genomic loci with 1-10 kbp resolution, our Hi-CO allows proximity analysis among DNA entry or exit points at every nucleosome locus. This analysis is realized by carrying out ligations among the entry/exit points in every nucleosome in a micrococcal-nuclease-fragmented genome, and by quantifying frequencies of ligation products with next-generation sequencing. Our protocol has enabled this analysis by cleanly excluding unwanted non-ligation products that are abundant owing to the frequent genome fragmentation by micrococcal nuclease. The computational procedure is based on simulated annealing-molecular dynamics, which allows determination of optimized 3D positions and orientations of every nucleosome that satisfies the proximity ligation data sufficiently well. Typically, examination of the Saccharomyces cerevisiae genome with 130 million sequencing reads facilitates analysis of a total of 66,360 nucleosome loci with 6.8 nm resolution. The technique requires 2-3 weeks for sequencing library preparation and 2 weeks for simulation.

Vasohibin-1 and small vasohibin-binding protein (SVBP) form an intermolecular heterotetramer in the crystal. The heterotetramer was stabilized by exchange of the conserved N-terminal region., Vasohibins regulate angiogenesis, tumor growth, metastasis and neuronal differentiation. They form a complex with small vasohibin-binding protein (SVBP) and show tubulin tyrosine carboxypeptidase activity. Recent crystal structure determinations of vasohibin–SVBP complexes have provided a molecular basis for complex formation, substrate binding and catalytic activity. However, the regulatory mechanism and dynamics of the complex remain elusive. Here, the crystal structure of the VASH1–SVBP complex and a molecular-dynamics simulation study are reported. The overall structure of the complex was similar to previously reported structures. Importantly, however, the structure revealed a domain-swapped heterotetramer that was formed between twofold symmetry-related molecules. This heterotetramerization was stabilized by the mutual exchange of ten conserved N-terminal residues from the VASH1 structural core, which was intramolecular in other structures. Interestingly, a comparison of this region with previously reported structures revealed that the patterns of hydrogen bonding and hydrophobic interactions vary. In the molecular-dynamics simulations, differences were found between the heterotetramer and heterodimer, where the fluctuation of the N-terminal region in the heterotetramer was suppressed. Thus, heterotetramer formation and flexibility of the N-terminal region may be important for enzyme activity and regulation.

tRNA

Highlights • Generalized epilepsy with tremor phenotype ss a clinical marker of 6q22.1 microdeletion. • Deep brain stimulation may be applicable choince for intractable myoclonic tremor. • Over activation of basal ganglia may be related to this patient's symptoms., We report a patient with a 6q22.1 deletion, who presented with a rare syndrome of generalized epilepsy, myoclonic tremor, and intellectual disability. There was no clinical progression after follow-up for more than 10 years. Our report presents the genetic basis for a phenotype involving a non-progressive generalized epilepsy with tremor. The efficacy of valproic acid for seizure control and the partial efficacy of deep brain stimulation with propranolol for myoclonic tremor is detailed.

Nanoarchaeum equitans is a species of hyperthermophilic archaea with the smallest genome size. Its alanyl-tRNA synthetase genes are split into AlaRS-α and AlaRS-β, encoding the respective subunits. In the current report, we surveyed N. equitans AlaRS-dependent alanylation of RNA minihelices, composed only of the acceptor stem and the T-arm of tRNAAla. Combination of AlaRS-α and AlaRS-β showed a strong alanylation activity specific to a single G3:U70 base pair, known to mark a specific tRNA for charging with alanine. However, AlaRS-α alone had a weak but appreciable alanylation activity that was independent of the G3:U70 base pair. The shorter 16-mer RNA tetraloop substrate mimicking only the first four base pairs of the acceptor stem of tRNAAla, but with C3:G70 base pair, was also successfully aminoacylated by AlaRS-α. The end of the acceptor stem, including CCA-3' terminus and the discriminator A73, was able to function as a minimal structure for the recognition by the enzyme. Our findings imply that aminoacylation by N. equitans AlaRS-α may represent a vestige of a primitive aminoacylation system, before the appearance of the G3:U70 pair as an identity element for alanine.

GENeralized-Ensemble SImulation System (GENESIS) is a software package for molecular dynamics (MD) simulation of biological systems. It is designed to extend limitations in system size and accessible time scale by adopting highly parallelized schemes and enhanced conformational sampling algorithms. In this new version, GENESIS 1.1, new functions and advanced algorithms have been added. The all-atom and coarse-grained potential energy functions used in AMBER and GROMACS packages now become available in addition to CHARMM energy functions. The performance of MD simulations has been greatly improved by further optimization, multiple time-step integration, and hybrid (CPU + GPU) computing. The string method and replica-exchange umbrella sampling with flexible collective variable choice are used for finding the minimum free-energy pathway and obtaining free-energy profiles for conformational changes of a macromolecule. These new features increase the usefulness and power of GENESIS for modeling and simulation in biological research. © 2017 Wiley Periodicals, Inc.

This study reports the X-ray crystallographic structure of the glycyl-tRNA synthetase (GlyRS) of Nanoarchaeum equitans — a hyperthermophilic archaeal species. This is the first archaeal GlyRS crystal structure elucidated. The GlyRS comprises an N-terminal catalytic domain and a C-terminal anticodon-binding domain with a long β-sheet inserted between these domains. An unmodified transcript of the wild-type N. equitans tRNAGly was successfully glycylated using GlyRS. Substitution of the discriminator base A73 of tRNAGly with any other nucleotide caused a significant decrease in glycylation activity. Mutational analysis of the second base-pair C2G71 of the acceptor stem of tRNAGly elucidated the importance of the base-pair, especially G71, as an identity element for recognition by GlyRS. Glycylation assays using tRNAGly G71 substitution mutants and a GlyRS mutant where Arg223 is mutated to alanine strengthen the possibility that the carbonyl oxygen at position 6 of G71 would hydrogen-bond with the guanidine nitrogen of Arg223 in N. equitans GlyRS.

The origin of homochirality in L-amino acid in proteins is one of the mysteries of the evolution of life. Experimental studies show that a non-enzymatic aminoacylation reaction of an RNA minihelix has a preference for L-amino acid over D-amino acid. The reaction initiates by approaching of a 3′-oxygen of the RNA minihelix to the carbonyl carbon of an aminoacyl phosphate oligonucleotide. Here, employing molecular dynamics simulations, we examined the possible mechanisms that determine this chiral selectivity. The simulation system adopted a geometry required for the chemical reaction to occur more frequently with L-alanine than that with D-alanine. For L-alanine, the structure with this geometry was formed by a combination of stable dihedral angles along alanyl phosphate backbone with a canonical RNA structure, where the methyl group of alanine was placed on the opposite side of the approaching 3′-hydroxyl group with respect to the carbonyl plane. For D-alanine, the methyl group and the 3′-hydroxyl group were placed on the same side with respect to the carbonyl plane, which significantly decreased its ability to approach 3′-oxygen close to the carbonyl carbon compared to L-alanine. The mechanism suggested herein can explain experimentally observed chiral preferences.

Many animals show left-right (LR) asymmetric morphology. The mechanisms of LR asymmetric development are evolutionarily divergent, and they remain elusive in invertebrates. Various organs in Drosophila melanogaster show stereotypic LR asymmetry, including the embryonic gut. The Drosophila embryonic hindgut twists 90° left-handedly, thereby generating directional LR asymmetry. We recently revealed that the hindgut epithelial cell is chiral in shape and other properties; this is termed planar cell chirality (PCC). We previously showed by computer modeling that PCC is sufficient to induce the hindgut rotation. In addition, both the PCC and the direction of hindgut twisting are reversed in Myosin31DF (Myo31DF) mutants. Myo31DF encodes Drosophila MyosinID, an actin-based motor protein, whose molecular functions in LR asymmetric development are largely unknown. Here, to understand how PCC directs the asymmetric cell-shape, we analyzed PCC in genetic mosaics composed of cells homozygous for mutant Myo31DF, some of which also overexpressed wild-type Myo31DF. Wild-type cell-shape chirality only formed in the Myo31DF-overexpressing cells, suggesting that cell-shape chirality was established in each cell and reflects intrinsic PCC. A computer model recapitulating the development of this genetic mosaic suggested that mechanical interactions between cells are required for the cell-shape behavior seen in vivo. Our mosaic analysis also suggested that during hindgut rotation in vivo, wild-type Myo31DF suppresses the elongation of cell boundaries, supporting the idea that cell-shape chirality is an intrinsic property determined in each cell. However, the amount and distribution of F-actin and Myosin II, which are known to help generate the contraction force on cell boundaries, did not show differences between Myo31DF mutant cells and wild-type cells, suggesting that the static amount and distribution of these proteins are not involved in the suppression of cell-boundary elongation. Taken together, our results suggest that cell-shape chirality is intrinsically formed in each cell, and that mechanical force from intercellular interactions contributes to its formation and/or maintenance.

The cytoplasm of a cell is crowded with many different kinds of macromolecules. The macromolecular crowding affects the thermodynamics and kinetics of biological reactions in a living cell, such as protein folding, association, and diffusion. Theoretical and simulation studies using simplified models focus on the essential features of the crowding effects and provide a basis for analyzing experimental data. In most of the previous studies on the crowding effects, a uniform crowder size is assumed, which is in contrast to the inhomogeneous size distribution of macromolecules in a living cell. Here, we evaluate the free energy changes upon macromolecular association in a cell-like inhomogeneous crowding system via a theory of hard-sphere fluids and free energy calculations using Brownian dynamics trajectories. The inhomogeneous crowding model based on 41 different types of macromolecules represented by spheres with different radii mimics the physiological concentrations of macromolecules in the cytoplasm of Mycoplasma genitalium. The free energy changes of macromolecular association evaluated by the theory and simulations were in good agreement with each other. The crowder size distribution affects both specific and non-specific molecular associations, suggesting that not only the volume fraction but also the size distribution of macromolecules are important factors for evaluating in vivo crowding effects. This study relates in vitro experiments on macromolecular crowding to in vivo crowding effects by using the theory of hard-sphere fluids with crowder-size heterogeneity.

Biological macromolecules function in highly crowded cellular environments. The structure and dynamics of proteins and nucleic acids are well characterized in vitro, but in vivo crowding effects remain unclear. Using molecular dynamics simulations of a comprehensive atomistic model cytoplasm we found that protein-protein interactions may destabilize native protein structures, whereas metabolite interactions may induce more compact states due to electrostatic screening. Protein-protein interactions also resulted in significant variations in reduced macromolecular diffusion under crowded conditions, while metabolites exhibited significant two-dimensional surface diffusion and altered protein-ligand binding that may reduce the effective concentration of metabolites and ligands in vivo. Metabolic enzymes showed weak non-specific association in cellular environments attributed to solvation and entropic effects. These effects are expected to have broad implications for the in vivo functioning of biomolecules. This work is a first step towards physically realistic in silico whole-cell models that connect molecular with cellular biology. DOI: http://dx.doi.org/10.7554/eLife.19274.001, eLife digest Much of the work that has been done to understand how cells work has involved studying parts of a cell in isolation. This is particularly true of studies that have examined the arrangement of atoms in large molecules with elaborate structures like proteins or DNA. However, cells are densely packed with many different molecules and there is little proof that proteins keep the same structures inside cells that they have when they are studied alone. To really understand how cells work, new ways to understand how molecules behave inside cells are needed. While this cannot be achieved directly, technology has now reached the stage where we can, to some extent, study living cells by recreating them virtually. Simulated cells can copy the atomic details of all the molecules in a cell and can estimate how different molecules might behave together. Yu et al. have now developed a computer simulation of part of a cell from the bacterium, Mycoplasma genitalium, one of the simplest forms of life on Earth. This model suggested new possible interactions between molecules inside cells that cannot currently be studied in real cells. The model shows that some proteins have a much less rigid structure in cells than they do in isolation, whilst others are able to work together more closely to carry out certain tasks. Finally, the model predicted that small molecules such as food, water and drugs would move more slowly through cells as they become stuck or trapped by larger molecules. These results could be particularly important in helping to improve drug design. Currently the simulations are limited, and can only model parts of simple cells for less than a thousandth of a second. However, in future it should be possible to recreate larger and more complex cells, including human cells, for longer periods of time. These could be used to better study human diseases and help to design new treatments. The ultimate goal is to simulate a whole cell in full detail by combining all the available experimental data. DOI: http://dx.doi.org/10.7554/eLife.19274.002

Elucidating the global and local rules that govern genome-wide, hierarchical chromatin architecture remains a critical challenge. Current high-throughput chromosome conformation capture (Hi-C) technologies have identified large-scale chromatin structural motifs, such as topologically associating domains and looping. However, structural rules at the smallest or nucleosome scale remain poorly understood. Here, we coupled nucleosome-resolved Hi-C technology with simulated annealing-molecular dynamics (SA-MD) simulation to reveal 3D spatial distributions of nucleosomes and their genome-wide orientation in chromatin. Our method, called Hi-CO, revealed distinct nucleosome folding motifs across the yeast genome. Our results uncovered two types of basic secondary structural motifs in nucleosome folding: α-tetrahedron and β-rhombus analogous to α helix and β sheet motifs in protein folding. Using mutants and cell-cycle-synchronized cells, we further uncovered motifs with specific nucleosome positioning and orientation coupled to epigenetic features at individual loci. By illuminating molecular-level structure-function relationships in eukaryotic chromatin, our findings establish organizational principles of nucleosome folding.

Hydrodynamic interactions (HI) give rise to collective motions between molecules, which are known to be important in the dynamics of random coil polymers and colloids. However, their role in the biological self-assembly of many molecule systems has not been investigated. Here, using Brownian dynamics simulations, we evaluate the importance of HI on the kinetics of self-assembly of lipid membranes. One-thousand coarse-grained lipid molecules in periodic simulation boxes were allowed to assemble into stable bilayers in the presence and absence of intermolecular HI. Hydrodynamic interactions reduce the monomer-monomer association rate by 50%. In contrast, the rate of association of lipid clusters is much faster in the presence of intermolecular HI. In fact, with intermolecular HI, the membrane self-assembly rate is 3–10 times faster than that without intermolecular HI. We introduce an analytical model to describe the size dependence of the diffusive encounter rate of particle clusters, which can qualitatively explain our simulation results for the early stage of the membrane self-assembly process. These results clearly suggest that HI greatly affects the kinetics of self-assembly and that simulations without HI will significantly underestimate the kinetic parameters of such processes.

The effect of macromolecular crowding on protein diffusion has been extensively studied in order to understand the protein dynamics in cellular environments. Previous studies have generally focused on proteins as crowding agents as they are known to occupy 20-30% of the cell volume. The volume fraction of the nucleoid in bacterial cells is also significant, around 10-20%, but the effect of chromosomal DNA as a crowding agent on protein diffusion has drawn relatively limited attention. Using recently build high-resolution models of bacterial nucleoids, Brownian and Stokesian dynamics simulations of coarse-grained model systems containing chromosomal DNA and proteins were carried out to investigate the effect of nucleoid crowding on protein diffusion. The coarse-grained model for the nucleoid is an experimentally-derived 3D model of the Caulobacter crescentus chromosome at base-pair resolution built by combining topological information of chromosomal DNA with distance restraints obtained from high-throughput Chromosome Conformation Capture (HiC) experiments, while proteins in the systems are modeled as hard-spheres. Results show that diffusion of proteins are slower around the nucleoid compared to in protein crowding and that there is size-dependent exclusion of proteins from the chromosome interior. Further analysis was carried to examine the degree of anomalous diffusion at different time and length scales.

To evaluate the enthalpic and entropic contributions of polar residues in alanine-based peptides to α-helical propensities of peptides, we applied a Brownian dynamics simulation, together with the umbrella sampling, to two alanine-based 21-residue peptides: one composed of alanine only (AAA), the other composed of 18 alanines and 3 arginines (ARA). Higher α-helical propensity of ARA than that of AAA was obtained. However, they showed similar conformational stability in enthalpy, considering the contribution of the solvation energy and the potential energy. As an evaluation of entropic effects, the fluctuation of a dihedral angle, ψ, was investigated. Arginine residue showed smaller fluctuation than alanine in elongated states. Higher α-helical propensity might originate from entropic effects. Furthermore, arginine seems to affect the α-helical propensity of alanines interacting with arginine.

To begin to elucidate the principles of intermolecular dynamics in the crowded environment of cells, employing Brownian dynamics (BD) simulations, we examined possible mechanism(s) responsible for the great reduction in diffusion constants of macromolecules in vivo from that at infinite dilution. In an Escherichia coli cytoplasm model comprised of 15 different macromolecule types at physiological concentrations, BD simulations of molecular-shaped and equivalent sphere representations were performed with a soft repulsive potential. At cellular concentrations, the calculated diffusion constant of GFP is much larger than experiment, with no significant shape dependence. Next, using the equivalent sphere system, hydrodynamic interactions (HI) were considered. Without adjustable parameters, the in vivo experimental GFP diffusion constant was reproduced. Finally, the effects of nonspecific attractive interactions were examined. The reduction in diffusivity is very sensitive to macromolecular radius with the motion of the largest macromolecules dramatically slowed down; this is not seen if HI dominate. In addition, long-lived clusters involving the largest macromolecules form if attractions dominate, whereas HI give rise to significant, size independent intermolecular dynamic correlations. These qualitative differences provide a testable means of differentiating the importance of HI vs. nonspecific attractive interactions on macromolecular motion in cells.

To evaluate the relationship between the amount of available Silicon (Si) in paddy soils and their mineral properties on the Shounai Plain in Japan, which is formed from several parent materials, we evaluated the amount of available Si, the particle size distribution, the oxide composition of crystalline minerals and the amount of oxalate-extractable Si (Sio), iron (Feo) and aluminum (Alo) in the soil. The amount of available Si in the soil and the oxide content of the crystalline minerals differed among four soil groups that were distinguished by their clay mineral composition. There was no difference in the particle size distribution among the soil groups. The amount of available Si was positively related to the SiO2/Al2O3 ratio of clay, the CaO concentration of silt and fine sand, and the amounts of Sio, Feo and Alo in the soil. The amount of available Si in the soils was negative correlated with the Na2O and K2O concentrations of silt, the K2O concentration of fine sand, and the coarse sand c...

ブラウン動力学（Brownian Dynamics；BD）法は溶媒を明示的に扱わないため計算を高速化することができる。従来のBD法は水環境でのタンパク質シミュレーションを対象にしていたため、膜タンパク質には適用できなかった。本研究では膜環境をも非明示的に再現することによってBD法を膜タンパク質のシミュレーションに適用できるように拡張した。αヘリックス構造のポリアラニンペプチドやパピロマーウイルス由来のE5タンパク質、蜂毒のメリチンをモデルとしてシミュレーションした結果、疎水性であるポリアラニンとE5膜タンパク質は膜中で安定に存在した。また、両親媒性であるメリチンペプチドは膜表面に安定に結合していた。これらの結果から、本膜モデルを用いたBD法は膜タンパク質のシミュレーションに有効であると考えられる。

Chromatin dynamics and arrangement are involved in many biological processes in nuclei of eukaryotes including plants. Plants have to respond rapidly to various environmental stimuli to achieve growth and development because they cannot move. It is assumed that the alteration of chromatin dynamics and arrangement support the response to these stimuli; however, there is little information in plants. In this study, we investigated the chromatin dynamics and arrangement with DNA damage in Arabidopsis thaliana by live-cell imaging with the lacO/LacI-EGFP system and simulation analysis. It was revealed that homologous loci kept a constant distance in nuclei of A. thaliana roots in general growth. We also found that DNA double-strand breaks (DSBs) induce the approach of the homologous loci with γ-irradiation. Furthermore, AtRAD54, which performs an important role in the homologous recombination repair pathway, was involved in the pairing of homologous loci with γ-irradiation. These results suggest that homologous loci approach each other to repair DSBs and AtRAD54 mediates these phenomena.

Interlaminar fracture toughness and delamination fatigue crack growth behavior were investigated for carbon fiber (CF)/epoxy laminates with the self-same epoxy interleaf. The matrix epoxy with a thickness of 50 μm was chosen as the interleaf material in order to clarify the effect of resin-rich layer thickness on the delamination fatigue crack growth behavior. Tests under mode I loading were carried out using double cantilever beam specimens. For tests under mode II loading, three-point end notched flexure specimens were used for interlaminar fracture toughness tests, while four-point end notched flexure specimens were used for delamination fatigue tests. The mode I properties (interlaminar fracture toughness and delamination fatigue threshold) of these epoxy-interleaved CFRP laminates were almost identical to those of the laminates without interleaf (base CFRP laminates). The effect of epoxy interleaf was completely different under mode II loading. The mode II interlaminar fracture toughness for the epoxy-interleaved laminates was 1.6 (initial value) and 3.4 (propagation value) times higher than that for the base CFRP laminates. The mode II delamination fatigue threshold of the epoxy-interleaved laminates was 2–2.3 times higher than those of the base CFRP laminates. While the toughness of the interleaf is the key factor under mode I, the thickness of interlayer is the key factor under mode II. The difference in the effect of the self-same epoxy interlayer on the interlaminar fracture properties under modes I and II loadings was discussed on the bases of the fractographic observations and mechanism considerations.

DNA binding proteins efficiently search for their cognitive sites on long genomic DNA by combining 3D diffusion and 1D diffusion (sliding) along the DNA. Recent experimental results and theoretical analyses revealed that the proteins show a rotation-coupled sliding along DNA helical pitch. Here, we performed Brownian dynamics simulations using newly developed coarse-grained protein and DNA models for evaluating how hydrodynamic interactions between the protein and DNA molecules, binding affinity of the protein to DNA, and DNA fluctuations affect the one dimensional diffusion of the protein on the DNA. Our results indicate that intermolecular hydrodynamic interactions reduce 1D diffusivity by 30%. On the other hand, structural fluctuations of DNA give rise to steric collisions between the CG-proteins and DNA, resulting in faster 1D sliding of the protein. Proteins with low binding affinities consistent with experimental estimates of non-specific DNA binding show hopping along the CG-DNA. This hopping significantly increases sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome., Author Summary DNA binding proteins efficiently search for their cognitive sites on long genomic DNA in cells to control biological activities. Recent experimental studies have revealed that the proteins use not only three-dimensional diffusion, but also one-dimensional diffusion (sliding) on DNA for this search process. For a better understanding of this biological process, we need to elucidate the mechanism of sliding. We report here molecular simulations using newly developed coarse-grained protein and DNA models for elucidating the nature of the sliding motions. Our simulation results show that: 1) hydrodynamic interactions between protein and DNA reduce sliding rate by 30%, 2) structural fluctuations of DNA give rise to steric collisions between proteins and DNA, which facilitate sliding motions, and 3) proteins with low binding affinities to DNA can hop along the DNA, resulting in a significant increase in sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome.

We applied an atomistic Brownian dynamics (BD) simulation with multiple time step method for the folding simulation of a 13-mer α-helical peptide and a 12-mer β-hairpin peptide, giving successful folding simulations. In this model, the driving energy contribution towards folding came from both electrostatic and van der Waals interactions for the α-helical peptide and from van der Waals interactions for the β-hairpin peptide. Although, many non-native structures having the same or lower energy than that of native structure were observed, the folded states formed the most populated cluster when the structures obtained by the BD simulations were subjected to the cluster analysis based on distance-based root mean square deviation of side-chains between different structures. This result indicates that we can predict the native structures from conformations sampled by BD simulation.

Structural genomics projects are beginning to produce protein structures with unknown functions, thereby creating a need for high-throughput methods to predict functions. Although sequence-based function prediction methods have been used extensively, structure-based prediction is believed to provide higher specificity and sensitivity because functions are closely related to the three-dimensional structures of functional sites, which are more strongly conserved during evolution than sequence. We have developed FCANAL, a method to predict functions using a score matrix obtained from the distances between Cα atoms and frequencies of appearance [1]. The previous report used key residues predicted from sequence comparisons (motifs). In this report, we have expanded the method to include enzymes and binding proteins with key residues predicted on the basis of three-dimensional structures. Using FCANAL, we constructed score matrices for 31 enzymes. When we applied them to all of the structure entries deposited in the Protein Data Bank, FCANAL could detect functional sites with high accuracy. This suggests that FCANAL will help identify the functions of newly determined structures and pinpoint their functionally important regions.

Object. Levodopa-induced dyskinesia (LID) in patients with Parkinson disease (PD) mimics acute dystonic reactions induced by antipsychotic agents, possibly mediated by σ-receptors; however, there are few reports in which the relationship between σ-receptors and LID in advanced PD is investigated. The binding potential of cerebellar σ-receptors before and after a pallidal surgery for dyskinesia in patients with advanced PD is assessed. Methods. Six patients with advanced PD (male/female ratio 3:3, age 56.7 ± 9.8 years) underwent stereotactic pallidal surgery (two posteroventral pallidotomy procedures and four deep brain stimulation of the globus pallidus internus, including one bilateral case). Clinical features of patients with PD were assessed using Hoehn and Yahr (H & Y) stages, the Unified Parkinson's Disease Rating Scale (UPDRS), and the Schwab and England Activities of Daily Life Scale (S & E). The LID was evaluated by LID severity score. The binding potential of cerebellar σ-receptors was determined before and after the surgery by 11C-nemonapride positron emission tomoraphy, a specific radioligand for σ-receptors in the cerebellum. All clinical scores, especially the LID severity score, were dramatically improved after the surgery (p < 0.05). Preoperatively, contralateral cerebellar binding potential was significantly elevated (p < 0.01), and it was reduced after the surgery, but it was still higher than that of healthy volunteers (p < 0.05). The ipsilateral cerebellar binding potential remained unchanged after the surgery. The level of binding potential did not correlate with H & Y stage, UPDRS, or S & E score, but a strong positive correlation was seen between the binding potential and the preoperative LID severity score when the patients were receiving medication (r = 0.893, p < 0.05). Conclusions. Cerebellar σ-receptors may potentially involve the genesis of LID in advanced PD.

A new simple implicit solvent model, effective charge (EC) model, was introduced into the Brownian dynamics algorithm based on AMBER united-atom force field. In the EC model, an atomic charge was decreased as a function of solvent-accessible surface area of the atom. We carried out the Brownian dynamics simulations of a 28-mer ββα fold peptide using four implicit solvent models: a generalized Born/solvent-accessible surface area (GB/SA) model, a solvent-accessible surface area (SA) based solvent model, a SA in combination with distance-dependent dielectric (DD/SA) and the EC combined with DD/SA (DD/SA/EC) model; and the calculated results on structure and dynamics of the peptide were compared with those of molecular dynamics simulation using explicit solvent model. Several artifacts were observed in the simulation using the GB/SA model. On the other hand, simulation using the DD/SA and DD/SA/EC implicit solvent models were free from such artifacts. Especially BD with the DD/SA/EC model gave the most stable trajectory as judged by root mean square deviations from the initial structure without large computational cost.

We report a multiple time step algorithm applied to an atomistic Brownian dynamics simulation for simulating the long time scale dynamics of biomolecules. The algorithm was based on the original multiple time step method; a short time step was used to keep faster motions in local equilibrium. When applied to a 28-mer # # ! folded peptide, the simulation gave stable trajectories and the computation time was reduced by a factor of 160 compared to a conventional molecular dynamics simulation using explicit water molecules. We applied it for the folding simulation of a 13-mer ! -helical peptide, giving a successful folding simulation. These results indicate that the Brownian dynamics with the multiple time step algorithm is useful for studies of biomolecular motions by long time simulation.

Decadal variability of subsurface temperature in the North Pacific has been investigated. Two dominant regions were found; the central subarctic region (CSa) and the north-eastern subtropical region (NESt). In CSa, cooling (warming) of wintertime subsurface temperature corresponds to the large (small) temperature gradient and southward (northward) shift of subsurface temperature front, associated with the increase (decrease) of positive wind stress curl and the southward (northward) shift of curl τ zero line with 2 years delay. It is suggested that the relocation of subtropical-subarctic boundary plays an important role. In NESt, importance of heat flux through the sea surface and heat divergence in the Ekman layer is also discussed.