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.

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.

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.

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.

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.