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.

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.

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.

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.

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.

The interiors of all living cells are highly crowded with macromolecules, which results in a considerable difference between the thermodynamics and kinetics of biological reactions in vivo from that in vitro. 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 E. coli cytoplasm modelcomprised of 15 different macromolecule types at physiological concentrations, where macromolecules were represented by spheres with their Stokes radii, BD simulations were performed with and without hydrodynamic interactions (HI). Without HI, the calculated diffusion constant of green fluorescent protein (GFP) is much larger than experiment. On the other hand, when HI were considered, the in vivo experimental GFP diffusion constant is almost reproduced without adjustable parameters. In addition, HI give rise to significant, size independent intermolecular dynamic correlations. These results suggest that HI play an important role on macromolecular dynamics in vivo.

Some organs in animals display left-right (LR) asymmetry. To better understand LR asymmetric morphogenesis in Drosophila, we studied LR directional rotation of the hindgut epithelial tube. Hindgut epithelial cells adopt a LR asymmetric (chiral) cell shape within their plane, and we refer to this cell behavior as planar cell-shape chirality (PCC). Drosophila E-cadherin (DE-Cad) is distributed to cell boundaries with LR asymmetry, which is responsible for the PCC formation. Myosin ID switches the LR polarity found in PCC and in DE-Cad distribution, which coincides with the direction of rotation. An in silico simulation showed that PCC is sufficient to induce the directional rotation of this tissue. Thus, the intrinsic chirality of epithelial cells in vivo is an underlying mechanism for LR asymmetric tissue morphogenesis.

The interiors of all living cells are highly crowded with macromolecules, which differs considerably the thermodynamics and kinetics of biological reactions between in vivo and in vitro. For example, the diffusion of green fluorescent protein (GFP) in E. coli is ~10-fold slower than in dilute conditions. In this study, we performed Brownian dynamics (BD) simulations of rigid macromolecules in a crowded environment mimicking the cytosol of E. coli to study the motions of macromolecules. The simulation systems contained 35 70S ribosomes, 750 glycolytic enzymes, 75 GFPs, and 392 tRNAs in a 100 nm × 100 nm × 100 nm simulation box, where the macromolecules were represented by rigid-objects of one bead per amino acid or four beads per nucleotide models. Diffusion tensors of these molecules in dilute solutions were estimated by using a hydrodynamic theory to take into account the diffusion anisotropy of arbitrary shaped objects in the BD simulations. BD simulations of the system where each macromolecule is represented by its Stokes radius were also performed for comparison. Excluded volume effects greatly reduce the mobility of molecules in crowded environments for both molecular-shaped and equivalent sphere systems. Additionally, there were no significant differences in the reduction of diffusivity over the entire range of molecular size between two systems. However, the reduction in diffusion of GFP in these systems was still 4–5 times larger than for the in vivo experiment. We will discuss other plausible factors that might cause the large reduction in diffusion in vivo.

Self-organization is one of the most fundamental phenomena in biology. Elucidating the mechanisms of the self-assembly of biological molecules, such as protein folding/binding and formation of lipid bilayer membranes, have been a fascinating topic over decades in biological field. Hydrodynamic interactions (HIs) give rise to collective motions between molecules. Therefore, they may play an important role on dynamics of self-assembly. Although there are significant advances in theoretical and computational studies on thermodynamics and kinetics of self-assembling systems recent years, analysis on effects of HIs to the reactions are limited due to their high complexity. Here, by using a Brownian dynamics (BD) simulation technique, we evaluate the importance of HIs on kinetic of self-assembly in a set of examples, formation of lipid membrane and actin polymerization reaction. We first built coarse-grained models of lipid and actin monomer. A cluster of lipid molecules was represented by two particles, one is hydrophilic head group and the other is hydrophobic tail group, and attractive interaction between tail groups was applied. For actin polymerization system, each monomer was represented by two particles and an angle-dependent attractive interaction was applied to form long filaments. In BD simulations, the Rotne-Prager-Yamakawa tensor was used to account for HIs between particles. In the early stage of self-assembly, HIs between inter-monomer particles decelerated the association rate of monomers to form oligomers compared with that in the simulations without inter-monomer HIs in both systems. On the other hand, the HIs accelerate the reaction of oligomers to form a fully-assembled state or much longer polymers. This phenomenon could be explained by the Kirkwood-Riseman theory. Those results clearly suggest that HIs greatly affect kinetics of self-assembly reactions and considering the interactions is crucial for studying dynamics of biological systems.