Free-energy landscape of reverse tRNA translocation through the ribosome analyzed by electron microscopy density maps and molecular dynamics simulations

Free Energy Landscape of H2A-H2B Displacement From Nucleosome

Nucleosome reconstitution plays an important role in many cellular functions. As an initial step, H2A-H2B dimer displacement, which is accompanied by disruption of many of the interactions within the nucleosome, should occur. To understand how H2A-H2B dimer displacement occurs, an adaptively biased molecular dynamics (ABMD) simulation was carried out to generate a variety of displacements of the H2A-H2B dimer from the fully wrapped to partially unwrapped nucleosome structures. With regards to these structures, the free energy landscape of the dimer displacement was investigated using umbrella sampling simulations. We found that the main contributors to the free energy were the docking domain of H2A and the C-terminal of H4. There were various paths for the dimer displacement which were dependent on the extent of nucleosomal DNA wrapping, suggesting that modulation of the intra-nucleosomal interaction by external factors such as histone chaperons could control the path for the H2A-H2B dimer displacement. Key residues which contributed to the free energy have also been reported to be involved in the mutations and posttranslational modifications (PTMs) which are important for assembling and/or reassembling the nucleosome at the molecular level and are found in cancer cells at the phenotypic level. Our results give insight into how the H2A-H2B dimer displacement proceeds along various paths according to different interactions within the nucleosome.

Probing Interplays between Human XBP1u Translational Arrest Peptide and 80S Ribosome

The ribosome stalling mechanism is a crucial biological process, yet its atomistic underpinning is still elusive. In this framework, the human XBP1u translational arrest peptide (AP) plays a central role in regulating the unfolded protein response (UPR) in eukaryotic cells. Here, we report multimicrosecond all-atom molecular dynamics simulations designed to probe the interactions between the XBP1u AP and the mammalian ribosome exit tunnel, both for the wild type AP and for four mutant variants of different arrest potencies. Enhanced sampling simulations allow investigating the AP release process of the different variants, shedding light on this complex mechanism. The present outcomes are in qualitative/quantitative agreement with available experimental data. In conclusion, we provide an unprecedented atomistic picture of this biological process and clear-cut insights into the key AP–ribosome interactions.

Torsional stress can regulate the unwrapping of two outer half superhelical turns of nucleosomal DNA

Torsional stress has a significant impact on the structure and stability of the nucleosome. RNA polymerase imposes torsional stress on the DNA in chromatin and unwraps the DNA from the nucleosome to access the genetic information encoded in the DNA. To understand how the torsional stress affects the stability of the nucleosome, we examined the unwrapping of two half superhelical turns of nucleosomal DNA from either end of the DNA under torsional stress with all-atom molecular dynamics simulations. The free energies for unwrapping the DNA indicate that positive stress that overtwists DNA facilitates a large-scale asymmetric unwrapping of the DNA without a large extension of the DNA. During the unwrapping, one end of the DNA was dissociated from H3 and H2A-H2B, while the other end of the DNA stably remained wrapped. The detailed analysis indicates that this asymmetric dissociation is facilitated by the geometry and bendability of the DNA under positive stress. The geometry stabilized the interaction between the major groove of the twisted DNA and the H3 αN-helix, and the straightened DNA destabilized the interaction with H2A-H2B. Under negative stress, the DNA became more bendable and flexible, which facilitated the binding of the unwrapped DNA to the octamer in a stable state. Consequently, we conclude that the torsional stress has a significant impact on the affinity of the DNA and the octamer through the inherent nature of the DNA and can change the accessibility of regulatory proteins.

Efficient Flexible Fitting Refinement with Automatic Error Fixing for De Novo Structure Modeling from Cryo-EM Density Maps

Structural modeling of proteins from cryo-electron microscopy (cryo-EM) density maps is one of the challenging issues in structural biology. De novo modeling combined with flexible fitting refinement (FFR) has been widely used to build a structure of new proteins. In de novo prediction, artificial conformations containing local structural errors such as chirality errors, cis peptide bonds, and ring penetrations are frequently generated and cannot be easily removed in the subsequent FFR. Moreover, refinement can be significantly suppressed due to the low mobility of atoms inside the protein. To overcome these problems, we propose an efficient scheme for FFR, in which the local structural errors are fixed first, followed by FFR using an iterative simulated annealing (SA) molecular dynamics protocol with the united atom (UA) model in an implicit solvent model; we call this scheme "SAUA-FFR". The best model is selected from multiple flexible fitting runs with various biasing force constants to reduce overfitting. We apply our scheme to the decoys obtained from MAINMAST and demonstrate an improvement of the best model of eight selected proteins in terms of the root-mean-square deviation, MolProbity score, and RWplus score compared to the original scheme of MAINMAST. Fixing the local structural errors can enhance the formation of secondary structures, and the UA model enables progressive refinement compared to the all-atom model owing to its high mobility in the implicit solvent. The SAUA-FFR scheme realizes efficient and accurate protein structure modeling from medium-resolution maps with less overfitting.

Free energy profile for unwrapping outer superhelical turn of CENP-A nucleosome

Eukaryotic genome is packaged in a nucleus in the form of chromatin. The fundamental structural unit of the chromatin is the protein-DNA complex, nucleosome, where DNA of about 150 bp is wrapped around a histone core almost twice. In cellular processes such as gene expression, DNA repair and duplication, the nucleosomal DNA has to be unwrapped. Histone proteins have their variants, indicating there are a variety of constitutions of nucleosomes. These different constitutions are observed in different cellular processes. To investigate differences among nucleosomes, we calculated free energy profiles for unwrapping the outer superhelical turn of CENP-A nucleosome and compared them with those of the canonical nucleosome. The free energy profiles for CENP-A nucleosome suggest that CENP-A nucleosome is the most stable when 16 to 22 bps are unwrapped in total whereas the canonical nucleosome is the most stable when it is fully wrapped. This indicates that the flexible conformation of CENP-A nucleosome is ready to provide binding sites for the structural integrity of the centromere.

Acceleration of cryo-EM Flexible Fitting for Large Biomolecular Systems by Efficient Space Partitioning

Flexible fitting is a powerful technique to build the 3D structures of biomolecules from cryoelectron microscopy (cryo-EM) density maps. One popular method is a cross-correlation coefficient-based approach, where the molecular dynamics (MD) simulation is carried out with the biasing potential that includes the cross-correlation coefficient between the experimental and simulated density maps. Here, we propose efficient parallelization schemes for the calculation of the cross-correlation coefficient to accelerate flexible fitting. Our schemes are tested for small, medium, and large biomolecules using CPU and hybrid CPU + GPU architectures. The scheme for the atomic decomposition MD is suitable for small proteins such as Ca²⁺-ATPase with the all-atom Go model, while that for the domain decomposition MD is better for larger systems such as ribosome with the all-atom Go or the all-atom explicit solvent models. Our methods allow flexible fitting for various biomolecules with reasonable computational cost. This approach also connects high-resolution structure refinements with investigation of protein structure-function relationship.

Free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA

The eukaryotic genome is packaged into a nucleus in the form of chromatin. The fundamental structural unit of chromatin is a protein-DNA complex, the nucleosome, where 146 or 147 base pairs of DNA wrap 1.75 times around a histone core. To function in cellular processes, however, nucleosomal DNA must be unwrapped. Although this unwrapping has been experimentally investigated, details of the process at an atomic level are not yet well understood. Here, we used molecular dynamics simulation with an enhanced sampling method to calculate the free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA. A free energy change of about 11.5 kcal/mol for the unwrapping agrees well with values obtained in single molecule experiments. This simulation revealed a variety of conformational states, indicating there are many potential paths to outer superhelicdal turn unwrapping, but the dominant path is likely asymmetric. At one end of the DNA, the first five bps unwrap, after which a second five bps unwrap at the same end with no increase in free energy. The unwrapping then starts at the other end of the DNA, where 10 bps are unwrapped. During further unwrapping of 15 bps, the unwrapping advances at one of the ends, after which the other end of the DNA unwraps to complete the unwrapping of the outer superhelical turn. These results provide insight into the construction, disruption, and repositioning of nucleosomes, which are continuously ongoing during cellular processes.

Molecular simulations of the ribosome and associated translation factors

The ribosome is a macromolecular complex which is responsible for protein synthesis in all living cells according to their transcribed genetic information. Using X-ray crystallography and, more recently, cryo-electron microscopy (cryo-EM), the structure of the ribosome was resolved at atomic resolution in many functional and conformational states. Molecular dynamics simulations have added information on dynamics and energetics to the available structural information, thereby have bridged the gap to the kinetics obtained from single-molecule and bulk experiments. Here, we review recent computational studies that brought notable insights into ribosomal structure and function.

H4 Tails Potentially Produce the Diversity in the Orientation of Two Nucleosomes

Histone tails play an important role in internucleosomal interaction and chromatin compaction. To understand how the H4 tails are involved in the internucleosomal interaction, an adaptively biased molecular dynamics simulation of 63 models of two stacked nucleosomes, each with the H4 tails in different locations, was carried out. This simulation generated a variety of orientations of the separated nucleosomes depending on the formation of the H4 tail bridge between the H4 tails and the DNA of the neighboring nucleosomes. For the models that showed distinctive orientations of the two nucleosomes, the free energies of the separation of the nucleosomes were further investigated using umbrella sampling simulations. The attractive force between the nucleosomes was estimated from the free energies; the force when two H4 tail bridges formed varied from 36 to 63 pN, depending on the formation of the H4 tail-bridge and the interfacial interaction, whereas the force reduced to 15-18 pN after either one of the H4 tail bridges had broken, regardless of the conformation of the H4 tail. Additional simulations of the nucleosomes show that when the H4 tail was truncated, the force between the nucleosomes became repulsive (from-3 to -7 pN). We concluded that the H4 tails potentially produce the diversity in the orientation of the two nucleosomes, which would contribute to the polymorphism of the chromatin structure.

Investigation of Ribosomes Using Molecular Dynamics Simulation Methods

The ribosome as a complex molecular machine undergoes significant conformational changes while synthesizing a protein molecule. Molecular dynamics simulations have been used as complementary approaches to X-ray crystallography and cryoelectron microscopy, as well as biochemical methods, to answer many questions that modern structural methods leave unsolved. In this review, we demonstrate that all-atom modeling of ribosome molecular dynamics is particularly useful in describing the process of tRNA translocation, atomic details of behavior of nascent peptides, antibiotics, and other small molecules in the ribosomal tunnel, and the putative mechanism of allosteric signal transmission to functional sites of the ribosome.

myPresto/omegagene: a GPU-accelerated molecular dynamics simulator tailored for enhanced conformational sampling methods with a non-Ewald electrostatic scheme

Molecular dynamics (MD) is a promising computational approach to investigate dynamical behavior of molecular systems at the atomic level. Here, we present a new MD simulation engine named "myPresto/omegagene" that is tailored for enhanced conformational sampling methods with a non-Ewald electrostatic potential scheme. Our enhanced conformational sampling methods, e.g., the virtual-system-coupled multi-canonical MD (V-McMD) method, replace a multi-process parallelized run with multiple independent runs to avoid inter-node communication overhead. In addition, adopting the non-Ewald-based zero-multipole summation method (ZMM) makes it possible to eliminate the Fourier space calculations altogether. The combination of these state-of-the-art techniques realizes efficient and accurate calculations of the conformational ensemble at an equilibrium state. By taking these advantages, myPresto/omegagene is specialized for the single process execution with Graphics Processing Unit (GPU). We performed benchmark simulations for the 20-mer peptide, Trp-cage, with explicit solvent. One of the most thermodynamically stable conformations generated by the V-McMD simulation is very similar to an experimentally solved native conformation. Furthermore, the computation speed is four-times faster than that of our previous simulation engine, myPresto/psygene-G. The new simulator, myPresto/omegagene, is freely available at the following URLs: http://www.protein.osaka-u.ac.jp/rcsfp/pi/omegagene/ and http://presto.protein.osaka-u.ac.jp/myPresto4/.

Modeling of ribosome dynamics on a ds-mRNA under an external load

Protein molecules in cells are synthesized by macromolecular machines called ribosomes. According to the recent experimental data, we reduce the complexity of the ribosome and propose a model to express its activity in six main states. Using our model, we study the translation rate in different biological relevant situations in the presence of external force and the translation through the RNA double stranded region in the absence or presence of the external force. In the present study, we give a quantitative theory for translation rate and show that the ribosome behaves more like a Brownian Ratchet motor. Our findings could shed some light on understanding behaviors of the ribosome in biological conditions.

Mechanism for verification of mismatched and homoduplex DNAs by nucleotides-bound MutS analyzed by molecular dynamics simulations

In order to understand how MutS recognizes mismatched DNA and induces the reaction of DNA repair using ATP, the dynamics of the complexes of MutS (bound to the ADP and ATP nucleotides, or not) and DNA (with mismatched and matched base-pairs) were investigated using molecular dynamics simulations. As for DNA, the structure of the base-pairs of the homoduplex DNA which interacted with the DNA recognition site of MutS was intermittently disturbed, indicating that the homoduplex DNA was unstable. As for MutS, the disordered loops in the ATPase domains, which are considered to be necessary for the induction of DNA repair, were close to (away from) the nucleotide-binding sites in the ATPase domains when the nucleotides were (not) bound to MutS. This indicates that the ATPase domains changed their structural stability upon ATP binding using the disordered loop. Conformational analysis by principal component analysis showed that the nucleotide binding changed modes which have structurally solid ATPase domains and the large bending motion of the DNA from higher to lower frequencies. In the MutS-mismatched DNA complex bound to two nucleotides, the bending motion of the DNA at low frequency modes may play a role in triggering the formation of the sliding clamp for the following DNA-repair reaction step. Moreover, MM-PBSA/GBSA showed that the MutS-homoduplex DNA complex bound to two nucleotides was unstable because of the unfavorable interactions between MutS and DNA. This would trigger the ATP hydrolysis or separation of MutS and DNA to continue searching for mismatch base-pairs. Proteins 2016; 84:1287-1303. © 2016 Wiley Periodicals, Inc.

Selective Preference of Parallel DNA Triplexes Is Due to the Disruption of Hoogsteen Hydrogen Bonds Caused by the Severe Nonisostericity between the G*GC and T*AT Triplets

Implications of DNA, RNA and RNA.DNA hybrid triplexes in diverse biological functions, diseases and therapeutic applications call for a thorough understanding of their structure-function relationships. Despite exhaustive studies mechanistic rationale for the discriminatory preference of parallel DNA triplexes with G*GC & T*AT triplets still remains elusive. Here, we show that the highest nonisostericity between the G*GC & T*AT triplets imposes extensive stereochemical rearrangements contributing to context dependent triplex destabilisation through selective disruption of Hoogsteen scheme of hydrogen bonds. MD simulations of nineteen DNA triplexes with an assortment of sequence milieu reveal for the first time fresh insights into the nature and extent of destabilization from a single (non-overlapping), double (overlapping) and multiple pairs of nonisosteric base triplets (NIBTs). It is found that a solitary pair of NIBTs, feasible either at a G*GC/T*AT or T*AT/G*GC triplex junction, does not impinge significantly on triplex stability. But two overlapping pairs of NIBTs resulting from either a T*AT or a G*GC interruption disrupt Hoogsteen pair to a noncanonical mismatch destabilizing the triplex by ~10 to 14 kcal/mol, implying that their frequent incidence in multiples, especially, in short sequences could even hinder triplex formation. The results provide (i) an unambiguous and generalised mechanistic rationale for the discriminatory trait of parallel triplexes, including those studied experimentally (ii) clarity for the prevalence of antiparallel triplexes and (iii) comprehensive perspectives on the sequence dependent influence of nonisosteric base triplets useful in the rational design of TFO's against potential triplex target sites.

Using Molecular Simulation to Model High-Resolution Cryo-EM Reconstructions

An explosion of new data from high-resolution cryo-electron microscopy (cryo-EM) studies has produced a large number of data sets for many species of ribosomes in various functional states over the past few years. While many methods exist to produce structural models for lower resolution cryo-EM reconstructions, high-resolution reconstructions are often modeled using crystallographic techniques and extensive manual intervention. Here, we present an automated fitting technique for high-resolution cryo-EM data sets that produces all-atom models highly consistent with the EM density. Using a molecular dynamics approach, atomic positions are optimized with a potential that includes the cross-correlation coefficient between the structural model and the cryo-EM electron density, as well as a biasing potential preserving the stereochemistry and secondary structure of the biomolecule. Specifically, we use a hybrid structure-based/ab initio molecular dynamics potential to extend molecular dynamics fitting. In addition, we find that simulated annealing integration, as opposed to straightforward molecular dynamics integration, significantly improves performance. We obtain atomistic models of the human ribosome consistent with high-resolution cryo-EM reconstructions of the human ribosome. Automated methods such as these have the potential to produce atomistic models for a large number of ribosome complexes simultaneously that can be subsequently refined manually.