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

Membrane lateral organization from potential energy disconnectivity graph

Understanding the thermodynamic and kinetic properties of biomolecules requires elucidation of their complex energy landscape. A disconnectivity graph analysis of the energy landscape provides a framework for mapping the multi-dimensional landscape onto a two-dimensional representation while preserving the key features of the energy landscape. Several studies show that the structure or shape of the disconnectity graph is directly associated with the function of protein and nucleic acid molecules. In this review, we discuss how disconnectivity analysis of the potential energy surface can be extended to lipid molecules to glean important information about membrane organization. The shape of the disconnectivity graphs can be used to predict the lateral organization of multi-component lipid bilayer. We hope that this review encourages the use of disconnectivity graphs routinely by membrane biophysicists to predict the lateral organization of lipids.

Recent advances in de novo computational design and redesign of intrinsically disordered proteins and intrinsically disordered protein regions

In the early 2000s, the concept of "unstructured biology" has emerged to be an important field in protein science by generating various new research directions. Many novel strategies and methods have been developed that are focused on effectively identifying/predicting intrinsically disordered proteins (IDPs) and intrinsically disordered protein regions (IDPRs), identifying their potential functions, disorder based drug design etc. Due to the range of functions of IDPs/IDPRs and their involvement in various debilitating diseases they are of contemporary interest to the scientific community. Recent researches are focused on designing/redesigning specific IDPs/IDPRs de novo. These de novo design/redesigns of IDPs/IDPRs are carried out by altering compositional biases and specific sequence patterning parameters. The main focus of these researches is to influence specific molecular functions, phase behavior, cellular phenotypes etc. In this review, we first provide the differences of natively folded and natively unfolded or IDPs with respect to their potential energy landscapes. Here, we provide current understandings on the different computational design strategies and methods that have been utilized in de novo design and redesigns of IDPs and IDPRs. Finally, we conclude the review by discussing the challenges that have been faced during the computational design/design attempts of IDPs/IDPRs.

Unusual Paradigm for DNA-DNA Recognition and Binding: "Socket-Plug" Complementarity

DNA is the key informational polymer in biology by virtue of its precisely defined self-assembling properties. Watson-Crick complementarity, which underlies DNA's self-assembly, is required not only in biology but has also proved powerful in the field of nanoscience, where it has been utilized to assemble complex 2D and 3D architectures and nanodevices built from the DNA double-helix. Aside from Watson-Crick base-pairing, however, DNA also participates in alternative base pairing schemes, giving rise to DNA triplexes and G-quadruplexes. Herein, we describe "sticky-ended" DNA triplex-quadruplex composites that specifically recognize and bind to each other using a wholly different logic, "socket-plug" complementarity, a shape-sensing fitting of guanine "prongs" into guanine-lacking "cavities." A remarkable property of this kind of complementarity is the key role played in it by specific counter-cations: thus, exclusive "self" socket-plug recognition occurs over "other" in sodium salt solutions while precisely the reverse occurs in potassium salt solutions. We have used gel electrophoresis, Förster resonance energy transfer, alkylation protection, and structural modeling to study this remarkable fundamental property of DNA, that we anticipate will find wide practical application.

The Energy Landscape Perspective: Encoding Structure and Function for Biomolecules

The energy landscape perspective is outlined with particular reference to biomolecules that perform multiple functions. We associate these multifunctional molecules with multifunnel energy landscapes, illustrated by some selected examples, where understanding the organisation of the landscape has provided new insight into function. Conformational selection and induced fit may provide alternative routes to realisation of multifunctionality, exploiting the possibility of environmental control and distinct binding modes.

Computer-aided comprehensive explorations of RNA structural polymorphism through complementary simulation methods

While RNA folding was originally seen as a simple problem to solve, it has been shown that the promiscuous interactions of the nucleobases result in structural polymorphism, with several competing structures generally observed for non-coding RNA. This inherent complexity limits our understanding of these molecules from experiments alone, and computational methods are commonly used to study RNA. Here, we discuss three advanced sampling schemes, namely Hamiltonian-replica exchange molecular dynamics (MD), ratchet-and-pawl MD and discrete path sampling, as well as the HiRE-RNA coarse-graining scheme, and highlight how these approaches are complementary with reference to recent case studies. While all computational methods have their shortcomings, the plurality of simulation methods leads to a better understanding of experimental findings and can inform and guide experimental work on RNA polymorphism.

Insights into the free energy landscape and salt-controlled mechanism of the conformational conversions between human telomeric G-quadruplex structures

G-quadruplexes have become attractive drug targets in cancer therapy. However, due to the polymorphism of G-quadruplex structures, it is difficult to experimentally verify the relevant structures of multiple intermediates and transition states in dynamic equilibrium. Hence, understanding the mechanism by which structural conversions of G-quadruplexes occur is still challenging. We conducted targeted molecular dynamics simulation with umbrella sampling to investigate how salt affects the conformational conversion of human telomeric G-quadruplex. Our results explore a unique view into the structures and energy barrier of the intermediates and transition states in the interconversion process. The pathway of G-quadruplex conformational interconversion was mapped out by a free energy landscape, consisting of branched parallel pathways with multiple energy basins. We propose a salt-controlled mechanism that as the salt concentration increases, the conformational conversion mechanism switches from multi-pathway folding to sequential folding pathways. The hybrid-I and hybrid-II structures are intermediates in the basket-propeller transformation. In high-salt solutions, the conformational conversion upon K⁺ binding is more feasible than upon Na⁺ binding. The free energy barrier for conformational conversions ranges from 1.6 to 4.6 kcal/mol. Our work will be beneficial in developing anticancer agents.

Biasing RNA Coarse-Grained Folding Simulations with Small-Angle X-ray Scattering Data

RNA molecules can easily adopt alternative structures in response to different environmental conditions. As a result, a molecule's energy landscape is rough and can exhibit a multitude of deep basins. In the absence of a high-resolution structure, small-angle X-ray scattering data (SAXS) can narrow down the conformational space available to the molecule and be used in conjunction with physical modeling to obtain high-resolution putative structures to be further tested by experiments. Because of the low resolution of these data, it is natural to implement the integration of SAXS data into simulations using a coarse-grained representation of the molecule, allowing for much wider searches and faster evaluation of SAXS theoretical intensity curves than with atomistic models. We present here the theoretical framework and the implementation of a simulation approach based on our coarse-grained model HiRE-RNA combined with SAXS evaluations "on-the-fly" leading the simulation toward conformations agreeing with the scattering data, starting from partially folded structures as the ones that can easily be obtained from secondary structure prediction-based tools. We show on three benchmark systems how our approach can successfully achieve high-resolution structures with remarkable similarity with the native structure recovering not only the overall shape, as imposed by SAXS data, but also the details of initially missing base pairs.

RNA Modeling with the Computational Energy Landscape Framework

The recent advances in computational abilities, such as the enormous speed-ups provided by GPU computing, allow for large scale computational studies of RNA molecules at an atomic level of detail. As RNA molecules are known to adopt multiple conformations with comparable energies, but different two-dimensional structures, all-atom models are necessary to better describe the structural ensembles for RNA molecules. This point is important because different conformations can exhibit different functions, and their regulation or mis-regulation is linked to a number of diseases. Problematically, the energy barriers between different conformational ensembles are high, resulting in long time scales for interensemble transitions. The computational potential energy landscape framework was designed to overcome this problem of broken ergodicity by use of geometry optimization. Here, we describe the algorithms used in the energy landscape explorations with the OPTIM and PATHSAMPLE programs, and how they are used in biomolecular simulations. We present a recent case study of the 5'-hairpin of RNA 7SK to illustrate how the method can be applied to interpret experimental results, and to obtain a detailed description of molecular properties.

Is the H4 histone tail intrinsically disordered or intrinsically multifunctional?

The structural versatility of histone tails is one of the key elements in the organisation of chromatin, which allows for the compact storage of genomic information. However, this structural diversity also complicates experimental and computational studies. Here, the potential and free energy landscape for the isolated and bound H4 histone tail are explored. The landscapes exhibit a set of distinct structural ensembles separated by high energy barriers, with little difference between isolated and bound tails. This consistency is a desirable feature that facilitates the formation of transient interactions, which are required for the liquid-like chromatin organisation. The existence of multiple, distinct structures on a multifunnel energy landscape is likely to be associated with multifunctionality, i.e. a set of evolved, distinct functions. Contrasting it with previously reported results for other disordered peptides, this type of landscape may be associated with a conformational selection based binding mechanism. Given the similarity to other systems exhibiting similar multifunnel energy landscapes, the disorder in histone tails might be better described in context of multifunctionality.

An in silico investigation of the binding modes and pathway of APTO-253 on c-KIT G-quadruplex DNA

The stability of c-KIT G-quadruplex DNA via ligands has been a significant concern in the growing field of cancer therapy. Thus, it is very important to understand the mechanism behind the high binding affinity of the small drug molecules on the c-KIT G-quadruplex DNA. In this study, we have investigated the binding mode and pathway of the APTO-253 ligand on the c-KIT G-quadruplex DNA employing a total of 10 μs all atom molecular dynamics simulations and further 8.82 μs simulations via the umbrella sampling method using both OL15 and BSC1 latest force fields for DNA structures. From the cluster structure analysis, mainly three binding pathways i.e., top, bottom and side loop stacking modes are identified. Moreover, RMSD, RMSF and 2D-RMSD values indicate that the c-KIT G-quadruplex DNA and APTO-253 molecules are stable throughout the simulation run. Furthermore, the number of hydrogen bonds in each tetrad and the distance between the two central K⁺ cations confirm that the c-KIT G-quadruplex DNA maintains its conformation in the process of complex formation with the APTO-253 ligand. The binding free energies and the minimum values in the potential of mean forces suggest that the binding processes are energetically favorable. Furthermore, we have found that the bottom stacking mode is the most favorable binding mode among all the three modes for the OL15 force field. However, for the BSC1 force field, both the top and bottom binding modes of the APTO-253 ligand in c-KIT G-quadruplex DNA are comparable to each other. To investigate the driving force for the complex formation, we have noticed that the van der Waals (vdW) and π-π stacking interactions are mainly responsible. Our detailed studies provide useful information for the discovery of novel drugs in the field of stabilization of G-quadruplex DNAs.

Insights into the Kinetic Partitioning Folding Dynamics of the Human Telomeric G-Quadruplex from Molecular Simulations and Machine Learning

The human telomeric DNA G-quadruplex follows a kinetic partitioning folding mechanism. The underlying folding landscape potentially has many minima separated by high free-energy barriers. However, using current theoretical models to characterize this complex folding landscape has remained a challenging problem. In this study, by developing a hybrid atomistic structure-based model that merges structural information on the hybrid-1, hybrid-2, and chair-type G-quadruplex topologies, we investigated a kinetic partitioning folding process of human telomeric DNA involving three native folds. The model was validated as it reproduced the experimental observation that the hybrid-1 conformation is the major fold and the hybrid-2 conformation is kinetically more accessible. A three-step mechanism was revealed for the formation of the hybrid-1 conformation, while a two-step mechanism was demonstrated for the formation of hybrid-2 and chair-type conformations. Likewise, a class of state in which structures adopted inappropriate combinations of syn/anti guanine nucleotides was found to greatly slow down the folding process. In addition, by employing the XGBoost machine learning algorithm, three interatom distances and six dihedral angles were identified as essential internal coordinates to represent the low-dimensional folding landscape. The strategy of coupling the multibasin model and the machine learning algorithm may be useful to investigate the conformational dynamics of other multistate biomolecules.

Understanding The Role of Reline, a Natural DES, on Temperature-Induced Conformational Changes of C-Kit G-Quadruplex DNA: A Molecular Dynamics Study

The noncanonical guanine-rich DNAs have drawn particular attention to the scientific world due to their controllable diverse and polymorphic structures. Apart from biological and medical significance, G-quadruplex DNAs are widely used in various fields such as nanotechnology, nanomachine, biosensors, and biocatalyst. So far, the applications of the G-quadruplex DNA are mainly limited in the water medium. Recently, a new generation of solvent named deep eutectic solvent (DES) has become very popular and has been widely used as a reaction medium of biocatalytic reactions and long-term storage medium for nucleic acids, even at high temperature. Hence, it is essential to understand the role of DES on temperature-induced conformational changes of a G-quadruplex DNA. In this research work, we have explored the temperature-mediated conformational dynamics of c-kit oncogene promoter G-quadruplex DNA in reline medium in the temperature range of 300-500 K, using a total of 10 μs unbiased all-atom molecular dynamics simulation. Here, from RMSD, RMSF, R_g and principal component analyses, we notice that the c-kit G-quadruplex DNA is stable up to 450 K in reline medium. However, it unfolds in water medium at 450 K. It is found that the hydrogen bonding interactions between c-kit G-quadruplex DNA and reline play a key role in the stabilization of the G-quadruplex DNA even at high temperature. Furthermore, in this work we have observed a very interesting and distinctive phenomenon of the central cation of the G-quadruplex DNA. Its position was seen to fluctuate between the two tetrad cores, that is, the region between tetrad-1 and tetrad-2 and that between tetrad-2 and tetrad-3 and vice versa at 450 and 500 K in reline medium which is absent in water medium at 450 K. Moreover, the rate of its oscillation is increased when temperature is increased.

Structural transitions in the RNA 7SK 5' hairpin and their effect on HEXIM binding

7SK RNA, as part of the 7SK ribonucleoprotein complex, is crucial to the regulation of transcription by RNA-polymerase II, via its interaction with the positive transcription elongation factor P-TEFb. The interaction is induced by binding of the protein HEXIM to the 5′ hairpin (HP1) of 7SK RNA. Four distinct structural models have been obtained experimentally for HP1. Here, we employ computational methods to investigate the relative stability of these structures, transitions between them, and the effects of mutations on the observed structural ensembles. We further analyse the results with respect to mutational binding assays, and hypothesize a mechanism for HEXIM binding. Our results indicate that the dominant structure in the wild type exhibits a triplet involving the unpaired nucleotide U40 and the base pair A43-U66 in the GAUC/GAUC repeat. This conformation leads to an open major groove with enough potential binding sites for peptide recognition. Sequence mutations of the RNA change the relative stability of the different structural ensembles. Binding affinity is consequently lost if these changes alter the dominant structure.

Parallel G-triplexes and G-hairpins as potential transitory ensembles in the folding of parallel-stranded DNA G-Quadruplexes

Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.

Coarse-grained dynamic RNA titration simulations

Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule's structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local pK _a values for each nucleotide.

Radicals Generated in Tetramolecular Guanine Quadruplexes by Photoionization: Spectral and Dynamical Features

G-quadruplexes are four-stranded DNA structures playing a key role in many biological functions and are promising for applications in the field of nanoelectronics. Characterizing the generation and fate of radical cations (electron holes) within these systems is important in relation to the DNA oxidative damage and/or conductivity issues. This study focuses on guanine radicals in G-quadruplexes formed by association of four TGGGGT strands in the presence of Na⁺ cations, (TG4T)₄/Na⁺. Using nanosecond transient spectroscopy with 266 nm excitation, we quantitatively characterize hydrated ejected electrons and three types of guanine radicals. We show that, at an energy lower by 2.7 eV than the guanine ionization potential, one-photon ionization occurs with quantum yield of (3.5 ± 0.5) × 10^-3. Deprotonation of the radical cations is completed within 20 μs, leading to the formation of (G-H2)^• radicals, following a strongly nonexponential decay pattern. Within 10 ms, the latter undergoes tautomerization to deprotonated (G-H1)^• radicals. The dynamics of the various radicals determined for (TG4T)₄/Na⁺, in connection to those reported previously for telomeric G-quadruplexes TEL21/Na⁺, is correlated with energetic factors computed by quantum chemical methods. The faster deprotonation of radical cations in (TG4T)₄/Na⁺ compared to TEL21/Na⁺ explains that irradiation of the former does not generate 8-oxodGuo, which is readily detected by high-performance liquid chromatography/mass spectrometry in the case of TEL21/Na⁺.

Conformational deviation of Thrombin binding G-quadruplex aptamer (TBA) in presence of divalent cation Sr2+: A classical molecular dynamics simulation study

Thrombin binding TBA-G-quadruplex aptamer (TBA) plays a major role in blood coagulation cascade. The 15-mer TBA sequence tends to form four-stranded TBA-G-quadruplex structure. In this research work, a series of explicit solvent classical MD simulations of the TBA is carried out using different salt (SrCl₂) concentrations (0, 50, 100 and 200 mM). Here we have also testified the effect of salt concentration of divalent cation Sr²⁺ on the conformational change of quadruplex DNA. The structural deviations, fluctuations, torsional angles and the affinity of the ion are explored at different salt concentrations. It is found that the conformation of TBA-G-quadruplex at 0 mM and 50 mM salt concentrations, is very much different than the other salt concentrations (100 mM and 200 mM). Also observed are as follows: (i) no exchange of Sr²⁺ ion between inside and outside of the channel, (ii) an enhancement in the Sr²⁺ ion density around the phosphate region of the loop residues as salt concentration increases and (iii) the stacking of T3 and T4 residues of loop-1 that appears up to 50 mM concentration, vanishes as the salt concentration is increased further.

Structural dynamics of propeller loop: towards folding of RNA G-quadruplex

We have carried out an extended set of standard and enhanced-sampling MD simulations (for a cumulative simulation time of 620 μs) with the aim to study folding landscapes of the rGGGUUAGGG and rGGGAGGG parallel G-hairpins (PH) with propeller loop. We identify folding and unfolding pathways of the PH, which is bridged with the unfolded state via an ensemble of cross-like structures (CS) possessing mutually tilted or perpendicular G-strands interacting via guanine-guanine H-bonding. The oligonucleotides reach the PH conformation from the unfolded state via a conformational diffusion through the folding landscape, i.e. as a series of rearrangements of the H-bond interactions starting from compacted anti-parallel hairpin-like structures. Although isolated PHs do not appear to be thermodynamically stable we suggest that CS and PH-types of structures are sufficiently populated during RNA guanine quadruplex (GQ) folding within the context of complete GQ-forming sequences. These structures may participate in compact coil-like ensembles that involve all four G-strands and already some bound ions. Such ensembles can then rearrange into the fully folded parallel GQs via conformational diffusion. We propose that the basic atomistic folding mechanism of propeller loops suggested in this work may be common for their formation in RNA and DNA GQs.

Energy Landscapes of Mini-Dumbbell DNA Octanucleotides

Single-stranded DNA structures play a significant role in biological systems, in particular during replication, translation, and DNA repair. Tracts of simple repetitive DNA are associated with slipped-strand mispairing, which can lead to genetic diseases. Recent NMR studies of TTTA and CCTG repeats have shown that these sequences form mini-dumbbells (MDBs), leading to frameshift mutations. Here we explore the energy landscapes of (CCTG)₂ and (TTTA)₂, which are currently the smallest known molecules to form MDBs. While (CCTG)₂ MDBs are stable, (TTTA)₂ exhibits numerous other structures with lower energies. A key factor identified in the stabilization of MDB structures is the bonding strength between residues 1 and 4, and 5 and 8.