Nejvíce citovaný článek - PubMed ID 26588288
Base Pair Fraying in Molecular Dynamics Simulations of DNA and RNA
The transition from B-DNA to A-DNA occurs in many protein-DNA interactions or in DNA/RNA hybrid duplexes, and thus plays a role in many important biomolecular processes that convey the biological function of DNA. However, the stability of A-DNA is severely underestimated in current AMBER force fields such as OL15, OL21 or bsc1, potentially leading to unstable or deformed protein-DNA complexes. In this study, we refine the deoxyribose dihedral potential to increase the stability of the north (N) puckering present in A-DNA. The new parameters, termed OL24, model A/B equilibrium in B-DNA duplexes in water in good agreement with nuclear magnetic resonance (NMR) experiment. They also improve the description of DNA/RNA hybrids and the transition of the DNA duplex to the A-form in concentrated ethanol solutions. These refinements significantly improve the modeling of protein-DNA complexes, increasing their structural stability and A-form population, while maintaining accurate representation of canonical B-DNA duplexes. Overall, the new parameters should allow more reliable modeling of the thermodynamic equilibrium between A- and B-DNA forms and the interactions of DNA with proteins.
The ribosome, owing to its exceptional conservation, harbours a remarkable molecular fossil known as the protoribosome. It surrounds the peptidyl transferase center (PTC), responsible for peptide bond formation. While previous studies have demonstrated the PTC activity in RNA alone, our investigation reveals the intricate roles of the ribosomal protein fragments (rPeptides) within the ribosomal core. This research highlights the significance of rPeptides in stability and coacervation of two distinct protoribosomal evolutionary stages. The 617nt 'big' protoribosome model, which associates with rPeptides specifically, exhibits a structurally defined and rigid nature, further stabilized by the peptides. In contrast, the 136nt 'small' model, previously linked to peptidyltransferase activity, displays greater structural flexibility. While this construct interacts with rPeptides with lower specificity, they induce coacervation of the 'small' protoribosome across a wide concentration range, which is concomitantly dependent on the RNA sequence and structure. Moreover, these conditions protect RNA from degradation. This phenomenon suggests a significant evolutionary advantage in the RNA-protein interaction at the early stages of ribosome evolution. The distinct properties of the two protoribosomal stages suggest that rPeptides initially provided compartmentalization and prevented RNA degradation, preceding the emergence of specific RNA-protein interactions crucial for the ribosomal structural integrity.
- MeSH
- konformace nukleové kyseliny MeSH
- molekulární modely MeSH
- peptidy chemie metabolismus MeSH
- peptidyltransferasy metabolismus chemie MeSH
- ribozomální proteiny * metabolismus chemie MeSH
- ribozomy * metabolismus MeSH
- RNA metabolismus chemie MeSH
- stabilita RNA MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- peptidy MeSH
- peptidyltransferasy MeSH
- ribozomální proteiny * MeSH
- RNA MeSH
Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.
- MeSH
- DNA * chemie MeSH
- konformace nukleové kyseliny * MeSH
- párování bází MeSH
- RNA * chemie MeSH
- simulace molekulární dynamiky * MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- DNA * MeSH
- RNA * MeSH
Advances in molecular dynamics (MD) software alongside enhanced computational power and hardware have allowed for MD simulations to significantly expand our knowledge of biomolecular structure, dynamics, and interactions. Furthermore, it has allowed for the extension of conformational sampling times from nanoseconds to the microsecond level and beyond. This has not only made convergence of conformational ensembles through comprehensive sampling possible but consequently exposed deficiencies and allowed the community to overcome limitations in the available force fields. The reproducibility and accuracy of the force fields are imperative in order to produce biologically relevant data. The Amber nucleic acid force fields have been used widely since the mid-1980s, and improvement of these force fields has been a community effort with several artifacts revealed, corrected, and reevaluated by various research groups. Here, we focus on the Amber force fields for use with double-stranded DNA and present the assessment of two recently developed force field parameter sets (OL21 and Tumuc1). Extensive MD simulations were performed with six test systems and two different water models. We observe the improvement of OL21 and Tumuc1 compared to previous generations of the Amber DNA force. We did not detect any significant improvement in the performance of Tumuc1 compared to OL21 despite the reparameterization of bonded force field terms in the former; however, we did note discrepancies in Tumuc1 when modeling Z-DNA sequences.
- MeSH
- DNA * chemie MeSH
- molekulární konformace MeSH
- reprodukovatelnost výsledků MeSH
- simulace molekulární dynamiky MeSH
- Z-DNA * MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- DNA * MeSH
- Z-DNA * MeSH
Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.
- MeSH
- DNA * MeSH
- křížová struktura DNA * MeSH
- molekulární konformace MeSH
- oprava DNA MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- DNA * MeSH
- křížová struktura DNA * MeSH
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
- MeSH
- DNA chemie MeSH
- katalýza MeSH
- konformace nukleové kyseliny * MeSH
- počítačová simulace MeSH
- RNA chemie MeSH
- simulace molekulární dynamiky * MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- přehledy MeSH
- Research Support, N.I.H., Extramural MeSH
- Názvy látek
- DNA MeSH
- RNA MeSH
The utility of molecular dynamics (MD) simulations to model biomolecular structure, dynamics, and interactions has witnessed enormous advances in recent years due to the availability of optimized MD software and access to significant computational power, including GPU multicore computing engines and other specialized hardware. This has led researchers to routinely extend conformational sampling times to the microsecond level and beyond. The extended sampling time has allowed the community not only to converge conformational ensembles through complete sampling but also to discover deficiencies and overcome problems with the force fields. Accuracy of the force fields is a key component, along with sampling, toward being able to generate accurate and stable structures of biopolymers. The Amber force field for nucleic acids has been used extensively since the 1990s, and multiple artifacts have been discovered, corrected, and reassessed by different research groups. We present a direct comparison of two of the most recent and state-of-the-art Amber force field modifications, bsc1 and OL15, that focus on accurate modeling of double-stranded DNA. After extensive MD simulations with five test cases and two different water models, we conclude that both modifications are a remarkable improvement over the previous bsc0 force field. Both force field modifications show better agreement when compared to experimental structures. To ensure convergence, the Drew-Dickerson dodecamer (DDD) system was simulated using 100 independent MD simulations, each extended to at least 10 μs, and the independent MD simulations were concatenated into a single 1 ms long trajectory for each combination of force field and water model. This is significantly beyond the time scale needed to converge the conformational ensemble of the internal portions of a DNA helix absent internal base pair opening. Considering all of the simulations discussed in the current work, the MD simulations performed to assess and validate the current force fields and water models aggregate over 14 ms of simulation time. The results suggest that both the bsc1 and OL15 force fields render average structures that deviate significantly less than 1 Å from the average experimental structures. This can be compared to similar but less exhaustive simulations with the CHARMM 36 force field that aggregate to the ∼90 μs time scale and also perform well but do not produce structures as close to the DDD NMR average structures (with root-mean-square deviations of 1.3 Å) as the newer Amber force fields. On the basis of these analyses, any future research involving double-stranded DNA simulations using the Amber force fields should employ the bsc1 or OL15 modification.
- MeSH
- B-DNA chemie MeSH
- DNA chemie MeSH
- konformace nukleové kyseliny MeSH
- magnetická rezonanční spektroskopie MeSH
- párování bází MeSH
- simulace molekulární dynamiky MeSH
- voda chemie MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- B-DNA MeSH
- DNA MeSH
- voda MeSH
We present a systematic study of the long-timescale dynamics of the Drew-Dickerson dodecamer (DDD: d(CGCGAATTGCGC)2) a prototypical B-DNA duplex. Using our newly parameterized PARMBSC1 force field, we describe the conformational landscape of DDD in a variety of ionic environments from minimal salt to 2 M Na(+)Cl(-) or K(+)Cl(-) The sensitivity of the simulations to the use of different solvent and ion models is analyzed in detail using multi-microsecond simulations. Finally, an extended (10 μs) simulation is used to characterize slow and infrequent conformational changes in DDD, leading to the identification of previously uncharacterized conformational states of this duplex which can explain biologically relevant conformational transitions. With a total of more than 43 μs of unrestrained molecular dynamics simulation, this study is the most extensive investigation of the dynamics of the most prototypical DNA duplex.
- MeSH
- B-DNA chemie ultrastruktura MeSH
- chlorid draselný chemie MeSH
- chlorid sodný chemie MeSH
- konformace nukleové kyseliny * MeSH
- molekulární modely MeSH
- simulace molekulární dynamiky * MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- B-DNA MeSH
- chlorid draselný MeSH
- chlorid sodný MeSH
DNA G-hairpins are potential key structures participating in folding of human telomeric guanine quadruplexes (GQ). We examined their properties by standard MD simulations starting from the folded state and long T-REMD starting from the unfolded state, accumulating ∼130 μs of atomistic simulations. Antiparallel G-hairpins should spontaneously form in all stages of the folding to support lateral and diagonal loops, with sub-μs scale rearrangements between them. We found no clear predisposition for direct folding into specific GQ topologies with specific syn/anti patterns. Our key prediction stemming from the T-REMD is that an ideal unfolded ensemble of the full GQ sequence populates all 4096 syn/anti combinations of its four G-stretches. The simulations can propose idealized folding pathways but we explain that such few-state pathways may be misleading. In the context of the available experimental data, the simulations strongly suggest that the GQ folding could be best understood by the kinetic partitioning mechanism with a set of deep competing minima on the folding landscape, with only a small fraction of molecules directly folding to the native fold. The landscape should further include non-specific collapse processes where the molecules move via diffusion and consecutive random rare transitions, which could, e.g. structure the propeller loops.
The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.
- MeSH
- denaturace nukleových kyselin MeSH
- draslík chemie MeSH
- G-kvadruplexy * MeSH
- kationty MeSH
- párování bází MeSH
- promotorové oblasti (genetika) * MeSH
- protoonkogenní proteiny c-kit genetika MeSH
- simulace molekulární dynamiky MeSH
- sodík chemie MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Názvy látek
- draslík MeSH
- kationty MeSH
- protoonkogenní proteiny c-kit MeSH
- sodík MeSH
