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.