Toward Force Fields with Improved Base Stacking Descriptions

Kinetics and dynamics of oligonucleotide hybridization

The hybridization of short nucleic acid strands is a remarkable spontaneous process that is foundational to biotechnology and nanotechnology and plays a crucial role in gene expression, editing and DNA repair. Decades of research into the mechanism of hybridization have resulted in a deep understanding of its thermodynamics, but many questions remain regarding its kinetics and dynamics. Recent advances in experiments and molecular dynamics simulations of nucleic acids are enabling more direct insight into the structural dynamics of hybridization, which can test long-standing assumptions regarding its mechanism. In this Review, we summarize the current state of knowledge of hybridization kinetics, discuss the barriers to a molecular description of hybridization dynamics, and highlight the new approaches that have begun uncovering the dynamics of hybridization and the duplex ensemble. The kinetics and dynamics of hybridization are highly sensitive to the composition of nucleic acids, and we emphasize recent discoveries and open questions on the role of nucleobase sequence and chemical modifications.

On the specificity of the recognition of m6A-RNA by YTH reader domains

Most processes of life are the result of polyvalent interactions between macromolecules, often of heterogeneous types and sizes. Frequently, the times associated with these interactions are prohibitively long for interrogation using atomistic simulations. Here, we study the recognition of N6-methylated adenine (m6A) in RNA by the reader domain YTHDC1, a prototypical, cognate pair that challenges simulations through its composition and required timescales. Simulations of RNA pentanucleotides in water reveal that the unbound state can impact (un)binding kinetics in a manner that is both model- and sequence-dependent. This is important because there are two contributions to the specificity of the recognition of the Gm6AC motif: from the sequence adjacent to the central adenine and from its methylation. Next, we establish a reductionist model consisting of an RNA trinucleotide binding to the isolated reader domain in high salt. An adaptive sampling protocol allows us to quantitatively study the dissociation of this complex. Through joint analysis of a data set including both the cognate and control sequences (GAC, Am6AA, and AAA), we derive that both contributions to specificity, sequence, and methylation, are significant and in good agreement with experimental numbers. Analysis of the kinetics suggests that flexibility in both the RNA and the YTHDC1 recognition loop leads to many low-populated unbinding pathways. This multiple-pathway mechanism might be dominant for the binding of unstructured polymers, including RNA and peptides, to proteins when their association is driven by polyvalent, electrostatic interactions.

The impact of sequence periodicity on DNA mechanics: investigating the origin of A-tract's curvature

Periodic sequences in phase with DNA helical shape are prevalent in genomes due to their capacity to modulate DNA elasticity on a global scale. However, how this occurs is not well understood. We use all-atom molecular dynamics simulations on 40 bp DNA fragments to assess the effect of periodicity on bending, twisting, and stretch elasticity. We observe that DNA static curvature is the mechanical parameter most influenced by periodicity, with A-tract sequences having the greatest effect. A-tracts generate global curvature by bending in distinct directions (minor groove and backbones) that complement the bending of the rest of DNA, which predominantly is towards the major groove. Even if A-tracts are rigid at the local scale, these small bends integrate with the greater bends from the sequences between, producing an amplifying effect. As a result, our findings support a 'delocalized bend' model in which the A-tract operates as an 'adaptable mechanical part'. By understanding how global curvature emerges from local fluctuations, we reconcile previous contradictory theories and open an avenue for manipulating DNA mechanics through sequence design.

Molecular dynamics simulations reveal the parallel stranded d(GGGA)3GGG DNA quadruplex folds via multiple paths from a coil-like ensemble

G-quadruplexes (G4s) are non-canonical nucleic acid structures that fold through complex processes. Characterization of the G4 folding landscape may help to elucidate biological roles of G4s but is challenging both experimentally and computationally. Here, we achieved complete folding of a three-quartet parallel DNA G4 with (GGGA)3GGG sequence using all-atom explicit-solvent enhanced-sampling molecular dynamics (MD) simulations. The simulations suggested early formation of guanine stacks in the G-tracts, which behave as semi-rigid blocks in the folding process. The folding continues via the formation of a collapsed compact coil-like ensemble. Structuring of the G4 from the coil then proceeds via various cross-like, hairpin, slip-stranded and two-quartet ensembles and can bypass the G-triplex structure. Folding of the parallel G4 does not appear to involve any salient intermediates and is a multi-pathway process. We also carried out an extended set of simulations of parallel G-hairpins. While parallel G-hairpins are extremely unstable when isolated, they are more stable inside the coil structure. On the methodology side, we show that the AMBER DNA force field predicts the folded G4 to be less stable than the unfolded ensemble, uncovering substantial force-field issues. Overall, we provide unique atomistic insights into the folding landscape of parallel-stranded G4 but also reveal limitations of current state-of-the-art MD techniques.

Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA

Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.