Debayan Chakraborty

Energy landscapes, folding mechanisms, and kinetics of RNA tetraloop hairpins.

RNA hairpins play a pivotal role in a diverse range of cellular functions, and are integral components of ribozymes, mRNA, and riboswitches. However, the mechanistic and kinetic details of RNA hairpin folding, which are key determinants of most of its biological functions, are poorly understood. In this work, we use the discrete path sampling (DPS) approach to explore the energy landscapes of two RNA tetraloop hairpins, and provide insights into their folding mechanisms and kinetics in atomistic detail. Our results show that the potential energy landscapes have a distinct funnel-like bias toward the folded hairpin state, consistent with efficient structure-seeking properties. Mechanistic and kinetic information is analyzed in terms of kinetic transition networks. We find microsecond folding times, consistent with temperature jump experiments, for hairpin folding initiated from relatively compact unfolded states. This process is essentially driven by an initial collapse, followed by rapid zippering of the helix stem in the final phase. Much lower folding rates are predicted when the folding is initiated from extended chains, which undergo longer excursions on the energy landscape before nucleation events can occur. Our work therefore explains recent experiments and coarse-grained simulations, where the folding kinetics exhibit precisely this dependency on the initial conditions.

Computational Investigation of RNA CUG Repeats Responsible for Myotonic Dystrophy 1

Myotonic Dystrophy 1 (DM1) is a genetic disease caused by expansion of CTG repeats in DNA. Once transcribed, these repeats form RNA hairpins with repeating 1×1 nucleotide UU internal loop motifs, r(CUG)n, which attract muscleblind-like 1 (MBNL1) protein leading to the disease. In DM1 CUG can be repeated thousands of times, so these structures are intractable to characterization using structural biology. However, inhibition of MBNL1-r(CUG)n binding requires a detailed analysis of the 1×1 UU internal loops. In this contribution we employ regular and umbrella sampling molecular dynamics (MD) simulations to describe the structural and thermodynamic properties of 1×1 UU internal loops. Calculations were run on a reported crystal structure and a designed system, which mimics an infinitely long RNA molecule with continuous CUG repeats. Two-dimensional (2D) potential of mean force (PMF) surfaces were created by umbrella sampling, and the discrete path sampling (DPS) method was utilized to investigate the energy landscape of 1×1 UU RNA internal loops, revealing that 1×1 UU base pairs are dynamic and strongly prefer the anti–anti conformation. Two 2D PMF surfaces were calculated for the 1×1 UU base pairs, revealing several local minima and three syn–anti ↔ anti–anti transformation pathways. Although at room temperature the syn–anti ↔ anti–anti transformation is not observed on the MD time scale, one of these pathways dominates the dynamics of the 1×1 UU base pairs in temperature jump MD simulations. This mechanism has now been treated successfully using the DPS approach. Our results suggest that local minima predicted by umbrella sampling calculations could be stabilized by small molecules, which is of great interest for future drug design. Furthermore, distorted GC/CG conformations may be important in understanding how MBNL1 binds to RNA CUG repeats. Hence we provide new insight into the dynamic roles of RNA loops and their contributions to presently incurable diseases.

Structure and properties of DNA in apolar solvents.

The study of nucleic acids in low-polarity environments paves the way for novel biotechnological applications of DNA. Here, we use a repertoire of atomistic molecular simulation tools to study the nature of DNA when placed in a highly apolar environment and when transferred from aqueous to apolar solvent. Our results show that DNA becomes stiffer in apolar solvents and suggest that highly negatively charged states, which are the most prevalent in water, are strongly disfavored in apolar solvents and neutral states with conformations not far from the aqueous ones are the dominant forms. Transfer from water to an apolar solvent such as CCl4 is unlikely to occur, but our results suggest that if forced, the DNA would migrate surrounded by a small shell of water (the higher the DNA charge, the larger the number of water molecules in this shell). Even the neutral form (predicted to be the dominant one in apolar solvents) would surround itself by a small number of highly stable water molecules when moved from water to a highly apolar environment. Neutralization of DNA charges seems a crucial requirement for transfer of DNA to apolar media, and the most likely mechanism to achieve good transfer properties.

Energy Landscape and Pathways for Transitions between Watson–Crick and Hoogsteen Base Pairing in DNA

The recent discovery that Hoogsteen (HG) base pairs are widespread in DNA across diverse sequences and positional contexts could have important implications for understanding DNA replication and DNA-protein recognition. While evidence is emerging that the Hoogsteen conformation could be a thermodynamically accessible conformation of the DNA duplex and provide a means to expand its functionality, relatively little is known about the molecular mechanism underlying the Watson-Crick (WC) to HG transition. In this Perspective, we describe pathways and kinetics for this transition at an atomic level of detail, using the energy landscape perspective. We show that competition between the duplex conformations results in a double funnel landscape, which explains some recent experimental observations. The interconversion pathways feature a number of intermediates, with a variable number of WC and HG base pairs. The relatively slow kinetics, with possible deviations from two-state behavior, suggest that this conformational switch is likely to be a challenging target for both simulation and experiment.

Multifunctional energy landscape for a DNA G-quadruplex: An evolved molecular switch

We explore the energy landscape for a four-fold telomere repeat, obtaining interconversion pathways between six experimentally characterised G-quadruplex topologies. The results reveal a multi-funnel system, with a variety of intermediate configurations and misfolded states. This organisation is identified with the intrinsically multi-functional nature of the system, suggesting a new paradigm for the classification of such biomolecules and clarifying issues regarding apparently conflicting experimental results.

Conformational Energy Landscape of the Ritonavir Molecule

Conformational polymorphism of ritonavir, a well-known pharmaceutical drug, is intricately linked to its efficacy in the treatment of acquired immunodeficiency syndrome (AIDS). Polymorphic transition from the crystalline form I to form II leads to the loss of bioactivity. The constituent ritonavir molecules adopt a trans configuration about the carbamate torsion angle in the form I crystal, and a cis configuration in the form II crystal. Investigating the energetics and mechanistic features of conformational transitions at the single molecule level is a key step toward decoding the complex features of the solid state polymorphism. In this work, we employ the energy landscape framework to investigate the conformational transitions of an isolated ritonavir molecule. The landscape is explored using discrete path sampling (DPS) and visualized in terms of disconnectivity graphs. We identify two distinct funnels corresponding to the two molecular forms that are identified by crystallography. The two regions can be reliably distinguished using the carbamate torsion angle, and the corresponding interconversion rates are predicted to follow Arrhenius behavior. The results provide mechanistic insight into pathways for cis ↔ trans interconversion at the molecular level and may also help in elucidating the polymorphic transitions in the crystal state.