We present parmbsc1, a new force-field for DNA atomistic simulation, which has been parameterized from high-level quantum mechanical data and tested for nearly 100 systems (~140 μs) covering most of the DNA structural space. Parmbsc1 provides high quality results in diverse systems, solving problems of previous force-fields. Parmbsc1 aims to be a reference force-field for the study of DNA in the next decade. Parameters and trajectories are available at http://mmb.irbbarcelona.org/ParmBSC1/.

Histone tails and their epigenetic modifications play crucial roles in gene expression regulation by altering the architecture of chromatin. However, the structural mechanisms by which histone tails influence the interconversion between active and inactive chromatin remain unknown. Given the technical challenges in obtaining detailed experimental characterizations of the structure of chromatin, multiscale computations offer a promising alternative to model the effect of histone tails on chromatin folding. Here we combine multi-microsecond atomistic molecular dynamics simulations of dinucleosomes and histone tails in explicit solvent and ions, performed with three different state-of-the-art force fields and validated by experimental NMR measurements, with coarse-grained Monte Carlo simulations of 24-nucleosome arrays to describe the conformational landscape of histone tails, their roles in chromatin compaction, and the impact of lysine acetylation, a widespread epigenetic change, on both. We find that while the wild-type tails are highly flexible and disordered, the dramatic increase of secondary-structure order by lysine acetylation unfolds chromatin by decreasing tail availability for crucial fiber-compacting internucleosome interactions. This molecular level description of the effect of histone tails and their charge modifications on chromatin folding explains the sequence sensitivity and underscores the delicate connection between local and global structural and functional effects. Our approach also opens new avenues for multiscale processes of biomolecular complexes.

Transporters of the amino acid, polyamine and organocation (APC) superfamily play essential roles in cell redox balance, cancer, and aminoacidurias. The bacterial L-arginine/agmatine antiporter, AdiC, is the main APC structural paradigm and shares the “5 + 5 inverted repeat” fold found in other families like the Na + -coupled neurotransmitter transporters. The available AdiC crystal structures capture two states of its transport cycle: the open-to-out apo and the outward-facing Arg + -bound occluded. However, the role of Arg + during the transition between these two states remains unknown. Here, we report the crystal structure at 3.0 Å resolution of an Arg + -bound AdiC mutant (N101A) in the open-to-out conformation, completing the picture of the major conformational states during the transport cycle of the 5 + 5 inverted repeat fold-transporters. The N101A structure is an intermediate state between the previous known AdiC conformations. The Arg + -guanidinium group in the current structure presents high mobility and delocalization, hampering substrate occlusion and resulting in a low translocation rate. Further analysis supports that proper coordination of this group with residues Asn101 and Trp293 is required to transit to the occluded state, providing the first clues on the molecular mechanism of substrate-induced fit in a 5 + 5 inverted repeat fold-transporter. The pseudosymmetry found between repeats in AdiC, and in all fold-related transporters, restraints the conformational changes, in particular the transmembrane helices rearrangements, which occur during the transport cycle. In AdiC these movements take place away from the dimer interface, explaining the independent functioning of each subunit.

Last generation of force-fields are raising expectations on the quality of molecular dynamics (MD) simulations of DNA, as well as to the belief that theoretical models can substitute experimental ones in several cases. However these claims are based on limited benchmarks, where MD simulations have shown the ability to reproduce already existing ‘experimental models’, which in turn, have an unclear accuracy to represent DNA conformation in solution. In this work we explore the ability of different force-fields to predict the structure of two new B-DNA dodecamers, determined herein by means of 1H nuclear magnetic resonance (NMR). The study allowed us to check directly for experimental NMR observables on duplexes previously not solved, and also to assess the reliability of ‘experimental structures’. We observed that technical details in the annealing procedures can induce non-negligible local changes in the final structures. We also found that while not all theoretical simulations are equally reliable, those obtained using last generation of AMBER force-fields (BSC1 and BSC0OL15) show predictive power in the multi-microsecond timescale and can be safely used to reproduce global structure of DNA duplexes and fine sequence-dependent details.

Hybrids of RNA with arabinonucleic acids 2′F-ANA and ANA have very similar structures but strikingly different thermal stabilities. We now present a thorough study combining NMR and other biophysical methods together with state-of-the-art theoretical calculations on a fully modified 10-mer hybrid duplex. Comparison between the solution structure of 2′F-ANA•RNA and ANA•RNA hybrids indicates that the increased binding affinity of 2′F-ANA is related to several subtle differences, most importantly a favorable pseudohydrogen bond (2′F–purine H8) which contrasts with unfavorable 2′-OH–nucleobase steric interactions in the case of ANA. While both 2′F-ANA and ANA strands maintained conformations in the southern/eastern sugar pucker range, the 2′F-ANA strand’s structure was more compatible with the A-like structure of a hybrid duplex. No dramatic differences are found in terms of relative hydration for the two hybrids, but the ANA•RNA duplex showed lower uptake of counterions than its 2′F-ANA•RNA counterpart. Finally, while the two hybrid duplexes are of similar rigidities, 2′F-ANA single strands may be more suitably preorganized for duplex formation. Thus the dramatically increased stability of 2′F-ANA•RNA and ANA•RNA duplexes is caused by differences in at least four areas, of which structure and pseudohydrogen bonding are the most important.

The local aromaticity of a series of benzenoid systems was determined through the use of structurally (HOMA) and magnetically (NICS) based measures, and also by using a new electronically based indicator of aromaticity, the para-delocalization index (PDI). The results were compared with the predictions of Clar's aromatic -sextet rule. It is found that, for all analyzed benzenoid hydrocarbons having a single Clar structure, local aromaticity orderings of the different six-membered rings given by all descriptors of aromaticity tested are identical and in agreement with Clar's aromatic -sextet rule. For benzenoid species that are described by a superposition of Clar structures, HOMA and PDI indices provide local aromaticity values that are totally consistent with Clar's rule, whereas NICS results deviate somewhat from the estimations based on this model.

BackgroundRNA secondary structures in the 5′-untranslated regions (5′-UTR) of mRNAs are key to the post-transcriptional regulation of gene expression. While it is evident that non-canonical Hoogsteen-paired G-quadruplex (rG4) structures somehow contribute to the regulation of translation initiation, the nature and extent of human mRNAs that are regulated by rG4s is not known. Here, we provide new insights into a mechanism by which rG4 formation modulates translation.ResultsUsing transcriptome-wide ribosome profiling, we identify rG4-driven mRNAs in HeLa cells and reveal that rG4s in the 5′-UTRs of inefficiently translated mRNAs associate with high ribosome density and the translation of repressive upstream open reading frames (uORF). We demonstrate that depletion of the rG4-unwinding helicases DHX36 and DHX9 promotes translation of rG4-associated uORFs while reducing the translation of coding regions for transcripts that comprise proto-oncogenes, transcription factors and epigenetic regulators. Transcriptome-wide identification of DHX9 binding sites shows that reduced translation is mediated through direct physical interaction between the helicase and its rG4 substrate.ConclusionThis study identifies human mRNAs whose translation efficiency is modulated by the DHX36- and DHX9-dependent folding/unfolding of rG4s within their 5′-UTRs. We reveal a previously unknown mechanism for translation regulation in which unresolved rG4s within 5′-UTRs promote 80S ribosome formation on upstream start codons, causing inhibition of translation of the downstream main open reading frames. Our findings suggest that the interaction of helicases with rG4s could be targeted for future therapeutic intervention.Electronic supplementary materialThe online version of this article (10.1186/s13059-018-1602-2) contains supplementary material, which is available to authorized users.

Extensive (more than 90 microseconds) molecular dynamics simulations complemented with ion-mobility mass spectrometry experiments have been used to characterize the conformational ensemble of DNA triplexes in the gas phase. Our results suggest that the ensemble of DNA triplex structures in the gas phase is well-defined over the experimental time scale, with the three strands tightly bound, and for the most abundant charge states it samples conformations only slightly more compact than the solution structure. The degree of structural alteration is however very significant, mimicking that found in duplex and much larger than that suggested for G-quadruplexes. Our data strongly supports that the gas phase triplex maintains an excellent memory of the solution structure, well-preserved helicity, and a significant number of native contacts. Once again, a linear, flexible, and charged polymer as DNA surprises us for its ability to retain three-dimensional structure in the absence of solvent. Results argue against the generally assumed roles of the different physical interactions (solvent screening of phosphate repulsion, hydrophobic effect, and solvation of accessible polar groups) in modulating the stability of DNA structures.