Polarizable Molecular Dynamics Simulations of Two c-kit Oncogene Promoter G-Quadruplexes: Effect of Primary and Secondary Structure on Loop and Ion Sampling

Can We Ever Develop an Ideal RNA Force Field? Lessons Learned from Simulations of the UUCG RNA Tetraloop and Other Systems

Molecular dynamics (MD) simulations are an important and well-established tool for investigating RNA structural dynamics, but their accuracy relies heavily on the quality of the employed force field (ff). In this work, we present a comprehensive evaluation of widely used pair-additive and polarizable RNA ffs using the challenging UUCG tetraloop (TL) benchmark system. Extensive standard MD simulations, initiated from the NMR structure of the 14-mer UUCG TL, revealed that most ffs did not maintain the native state, instead favoring alternative loop conformations. Notably, three very recent variants of pair-additive ffs, OL3CP-gHBfix21, DES-Amber, and OL3R2.7, successfully preserved the native structure over a 10 × 20 μs time scale. To further assess these ffs, we performed enhanced sampling folding simulations of the shorter 8-mer UUCG TL, starting from the single-stranded conformation. Estimated folding free energies (ΔG°fold) varied significantly among these three ffs, with values of 0.0 ± 0.6, 2.4 ± 0.8, and 7.4 ± 0.2 kcal/mol for OL3CP-gHBfix21, DES-Amber, and OL3R2.7, respectively. The ΔG°fold value predicted by the OL3CP-gHBfix21 ff was closest to experimental estimates, ranging from -1.6 to -0.7 kcal/mol. In contrast, the higher ΔG°fold values obtained using DES-Amber and OL3R2.7 were unexpected, suggesting that key interactions are inaccurately described in the folded, unfolded, or misfolded ensembles. These discrepancies led us to further test DES-Amber and OL3R2.7 ffs on additional RNA and DNA systems, where further performance issues were observed. Our results emphasize the complexity of accurately modeling RNA dynamics and suggest that creating an RNA ff capable of reliably performing across a wide range of RNA systems remains extremely challenging. In conclusion, our study provides valuable insights into the capabilities of current RNA ffs and highlights key areas for future ff development.

Insights into the Molecular Structure, Stability, and Biological Significance of Non-Canonical DNA Forms, with a Focus on G-Quadruplexes and i-Motifs

This article provides a comprehensive examination of non-canonical DNA structures, particularly focusing on G-quadruplexes (G4s) and i-motifs. G-quadruplexes, four-stranded structures formed by guanine-rich sequences, are stabilized by Hoogsteen hydrogen bonds and monovalent cations like potassium. These structures exhibit diverse topologies and are implicated in critical genomic regions such as telomeres and promoter regions of oncogenes, playing significant roles in gene expression regulation, genome stability, and cellular aging. I-motifs, formed by cytosine-rich sequences under acidic conditions and stabilized by hemiprotonated cytosine-cytosine (C:C+) base pairs, also contribute to gene regulation despite being less prevalent than G4s. This review highlights the factors influencing the stability and dynamics of these structures, including sequence composition, ionic conditions, and environmental pH. Molecular dynamics simulations and high-resolution structural techniques have been pivotal in advancing our understanding of their folding and unfolding mechanisms. Additionally, the article discusses the therapeutic potential of small molecules designed to selectively bind and stabilize G4s and i-motifs, with promising implications for cancer treatment. Furthermore, the structural properties of these DNA forms are explored for applications in nanotechnology and molecular devices. Despite significant progress, challenges remain in observing these structures in vivo and fully elucidating their biological functions. The review underscores the importance of continued research to uncover new insights into the genomic roles of G4s and i-motifs and their potential applications in medicine and technology. This ongoing research promises exciting developments in both basic science and applied fields, emphasizing the relevance and future prospects of these intriguing DNA structures.

Quantifying Induced Dipole Effects in Small Molecule Permeation in a Model Phospholipid Bilayer

The cell membrane functions as a semipermeable barrier that governs the transport of materials into and out of cells. The bilayer features a distinct dielectric gradient due to the amphiphilic nature of its lipid components. This gradient influences various aspects of small molecule permeation and the folding and functioning of membrane proteins. Here, we employ polarizable molecular dynamics simulations to elucidate the impact of the electronic environment on the permeation process. We simulated eight distinct amino-acid side chain analogs within a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer using the Drude polarizable force field (FF). Our approach includes both unbiased and umbrella sampling simulations. By using a polarizable FF, we sought to investigate explicit dipole responses in relation to local electric fields along the membrane normal. We evaluate molecular dipole moments, which exhibit variation based on their localization within the membrane, and compare the outcomes with analogous simulations using the nonpolarizable CHARMM36 FF. This comparative analysis aims to discern characteristic differences in the free energy surfaces of permeation for the various amino-acid analogs. Our results provide the first systematic quantification of the impact of employing an explicitly polarizable FF in this context compared to the fixed-charge convention inherent to nonpolarizable FFs, which may not fully capture the influence of the membrane dielectric gradient.

Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junction

G-quadruplexes (GQs) play key regulatory roles within the human genome and have also been identified to play similar roles in other eukaryotes, bacteria, archaea, and viruses. Human immunodeficiency virus 1, the etiological agent of acquired immunodeficiency syndrome, can form two GQs in its long terminal repeat (LTR) promoter region, each of which act to regulate viral gene expression in opposing manners. The major LTR GQ, called LTR-III, is a distinct hybrid GQ containing a 12-nucleotide duplex loop attached to the quadruplex motif. The resulting quadruplex:duplex junction (QDJ) has been hypothesized to serve as a selective drug targeting site. To better understand the dynamics of this QDJ, we performed conventional and enhanced-sampling molecular dynamics simulations using the Drude-2017 force field. We observed unbiased and reversible formation of additional base pairs in the QDJ, between Ade4:Thy14 and Gua3:Thy14. Both base pairs were electrostatically favored, but geometric constraints within the junction may drive the formation of, and preference for, the Ade4:Thy14 base pair. Finally, we demonstrated that the base pairs are separated only by small energy barriers that may enable transitions between both base-paired states. Together, these simulations provide new insights into the dynamics, electrostatics, and thermodynamics of the LTR-III QDJ.

RNA G-quadruplex folding is a multi-pathway process driven by conformational entropy

The kinetics of folding is crucial for the function of many regulatory RNAs including RNA G-quadruplexes (rG4s). Here, we characterize the folding pathways of a G-quadruplex from the telomeric repeat-containing RNA by combining all-atom molecular dynamics and coarse-grained simulations with circular dichroism experiments. The quadruplex fold is stabilized by cations and thus, the ion atmosphere forming a double layer surrounding the highly charged quadruplex guides the folding process. To capture the ionic double layer in implicit solvent coarse-grained simulations correctly, we develop a matching procedure based on all-atom simulations in explicit water. The procedure yields quantitative agreement between simulations and experiments as judged by the populations of folded and unfolded states at different salt concentrations and temperatures. Subsequently, we show that coarse-grained simulations with a resolution of three interaction sites per nucleotide are well suited to resolve the folding pathways and their intermediate states. The results reveal that the folding progresses from unpaired chain via hairpin, triplex and double-hairpin constellations to the final folded structure. The two- and three-strand intermediates are stabilized by transient Hoogsteen interactions. Each pathway passes through two on-pathway intermediates. We hypothesize that conformational entropy is a hallmark of rG4 folding. Conformational entropy leads to the observed branched multi-pathway folding process for TERRA25. We corroborate this hypothesis by presenting the free energy landscapes and folding pathways of four rG4 systems with varying loop length.

Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes

The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication and maintaining the genome, among others; thus, changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and telomeric repeat-containing RNA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on the GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in the electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle but important differences between DNA and RNA GQs with the same fold.

Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes

The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication; and maintaining the genome, among others, thus changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and TERRA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle, but important, differences between DNA and RNA GQs with the same fold.

Computer Aided Development of Nucleic Acid Applications in Nanotechnologies

Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.

Ion-Dependent Conformational Plasticity of Telomeric G-Hairpins and G-Quadruplexes

Telomeric DNA is guanine-rich and can adopt structures such as G-quadruplexes (GQs) and G-hairpins. Telomeric GQs influence genome stability and telomerase activity, making understanding of enzyme-GQ interactions and dynamics important for potential drug design. GQs have a characteristic tetrad core, which is connected by loop regions. Within this architecture are G-hairpins, fold-back motifs that are thought to represent the first intermediate in GQ folding. To better understand the relationship between G-hairpin motifs and GQs, we performed polarizable simulations of a two-tetrad telomeric GQ and an isolated SC11 telomeric G-hairpin. The telomeric GQ contains a G-triad, which functions as part of the tetrad core or linker regions, depending on local conformational change. This triad and another motif below the tetrad core frequently bound ions and may represent druggable sites. Further, we observed the unbiased formation of a G-triad and a G-tetrad in simulations of the SC11 G-hairpin and found that cations can be partially hydrated while facilitating the formation of these motifs. Finally, we demonstrated that K⁺ ions form specific interactions with guanine bases, while Na⁺ ions interact nonspecifically with bases in the structure. Together, these simulations provide new insights into the influence of ions on GQs, G-hairpins, and G-triad motifs.

Refining the RNA Force Field with Small-Angle X-ray Scattering of Helix-Junction-Helix RNA

The growing recognition of the functional and therapeutic roles played by RNA and the difficulties in gaining atomic-level insights by experiments are paving the way for all-atom simulations of RNA. One of the main impediments to the use of all-atom simulations is the imbalance between the energy terms of the RNA force fields. Through exhaustive sampling of an RNA helix-junction-helix (HJH) model using enhanced sampling, we critically assessed the select Amber force fields against small-angle X-ray scattering (SAXS) experiments. The tested AMBER99SB, DES-AMBER, and CUFIX force fields show deviations from measured profiles. First, we identified parameters leading to inconsistencies. Then, as a way to balance the forces governing RNA folding, we adopted strategies to refine hydrogen bonding, backbone, and base-stacking parameters. We validated the modified force field (HB-CUFIX) against SAXS data of the HJH model in different ionic strengths. Moreover, we tested a set of independent RNA systems to cross-validate the force field. Overall, HB-CUFIX demonstrates improved performance in studying thermodynamics and structural properties of realistic RNA motifs.

Polymorphic and Higher-Order G-Quadruplexes as Possible Transcription Regulators: Novel Perspectives for Future Anticancer Therapeutic Applications

In the past two decades, significant efforts have been put into designing small molecules to target selected genomic sites where DNA conformational rearrangements control gene expression. G-rich sequences at oncogene promoters are considered good points of intervention since, under specific environmental conditions, they can fold into non-canonical tetrahelical structures known as G-quadruplexes. However, emerging evidence points to a frequent lack of correlation between small molecule targeting of G-quadruplexes at gene promoters and the expression of the associated protein, which hampers pharmaceutical applications. The wide genomic localization of G-quadruplexes along with their highly polymorphic behavior may account for this scenario, suggesting the need for more focused drug design strategies. Here, we will summarize the G4 structural features that can be considered to fulfill this goal. In particular, by comparing a telomeric sequence with the well-characterized G-rich domain of the KIT promoter, we will address how multiple secondary structures might cooperate to control genome architecture at a higher level. If this holds true, the link between drug–DNA complex formation and the associated cellular effects will need to be revisited.

Computational and Experimental Characterization of rDNA and rRNA G-Quadruplexes

DNA G-quadruplexes in human telomeres and gene promoters are being extensively studied for their role in controlling the growth of cancer cells. G-quadruplexes have been unambiguously shown to exist both in vitro and in vivo, including in the guanine (G)-rich DNA genes encoding pre-ribosomal RNA (pre-rRNA), which is transcribed in the cell's nucleolus. Recent studies strongly suggest that these DNA sequences ("rDNA"), and the transcribed rRNA, are a potential anticancer target through the inhibition of RNA polymerase I (Pol I) in ribosome biogenesis, but the structures of ribosomal G-quadruplexes at atomic resolution are unknown and very little biophysical characterization has been performed on them to date. In the present study, circular dichroism (CD) spectroscopy is used to show that two putative rDNA G-quadruplex sequences, NUC 19P and NUC 23P and their counterpart rRNAs, predominantly adopt parallel topologies, reminiscent of the analogous telomeric quadruplex structures. Based on this information, we modeled parallel topology atomistic structures of the putative ribosomal G-quadruplexes. We then validated and refined the modeled ribosomal G-quadruplex structures using all-atom molecular dynamics (MD) simulations with the CHARMM36 force field in the presence and absence of stabilizing K⁺. Motivated by preliminary MD simulations of the telomeric parallel G-quadruplex (TEL 24P) in which the K⁺ ion is expelled, we used updated CHARMM36 force field K⁺ parameters that were optimized, targeting the data from quantum mechanical calculations and the polarizable Drude model force field. In subsequent MD simulations with optimized CHARMM36 parameters, the K⁺ ions are predominantly in the G-quadruplex channel and the rDNA G-quadruplexes have more well-defined, predominantly parallel-topology structures as compared to rRNA. In addition, NUC 19P is more structured than NUC 23P, which contains extended loops. Results from this study set the structural foundation for understanding G-quadruplex functions and the design of novel chemotherapeutics against these nucleolar targets and can be readily extended to other DNA and RNA G-quadruplexes.

Double-stranded flanking ends affect the folding kinetics and conformational equilibrium of G-quadruplexes forming sequences within the promoter of KIT oncogene

G-quadruplexes embedded within promoters play a crucial role in regulating the gene expression. KIT is a widely studied oncogene, whose promoter contains three G-quadruplex forming sequences, c-kit1, c-kit2 and c-kit*. For these sequences available studies cover ensemble and single-molecule analyses, although for kit* the latter were limited to a study on a promoter domain comprising all of them. Recently, c-kit2 has been reported to fold according to a multi-step process involving folding intermediates. Here, by exploiting fluorescence resonance energy transfer, both in ensemble and at the single molecule level, we investigated the folding of expressly designed constructs in which, alike in the physiological context, either c-kit2 or c-kit* are flanked by double stranded DNA segments. To assess whether the presence of flanking ends at the borders of the G-quadruplex affects the folding, we studied under the same protocols oligonucleotides corresponding to the minimal G-quadruplex forming sequences. Data suggest that addition of flanking ends results in biasing both the final equilibrium state and the folding kinetics. A previously unconsidered aspect is thereby unravelled, which ought to be taken into account to achieve a deeper insight of the complex relationships underlying the fine tuning of the gene-regulatory properties of these fascinating DNA structures.

Combining Electrospray Mass Spectrometry (ESI-MS) and Computational Techniques in the Assessment of G-Quadruplex Ligands: A Hybrid Approach to Optimize Hit Discovery

Guanine-rich sequences forming G-quadruplexes (GQs) are present in several genomes, ranging from viral to human. Given their peculiar localization, the induction of GQ formation or GQ stabilization with small molecules represents a strategy for interfering with crucial biological functions. Investigating the recognition event at the molecular level, with the aim of fully understanding the triggered pharmacological effects, is challenging. Native electrospray ionization mass spectrometry (ESI-MS) is being optimized to study these noncovalent assemblies. Quantitative parameters retrieved from ESI-MS studies, such as binding affinity, the equilibrium binding constant, and sequence selectivity, will be overviewed. Computational experiments supporting the ESI-MS investigation and boosting its efficiency in the search for GQ ligands will also be discussed with practical examples. The combination of ESI-MS and in silico techniques in a hybrid high-throughput-screening workflow represents a valuable tool for the medicinal chemist, providing data on the quantitative and structural aspects of ligand–GQ interactions.

Benchmark Force Fields for the Molecular Dynamic Simulation of G-Quadruplexes

G-quadruplexes have drawn widespread attention for serving as a potential anti-cancer target and their application in material science. Molecular dynamics (MD) simulation is the key theoretical tool in the study of GQ's structure-function relationship. In this article, we systematically benchmarked the five force fields of parmbsc0, parmbsc1, OL15, AMOEBA, and Drude2017 on the MD simulation of G-quadruplex from four aspects: structural stability, central ion channel stability, description of Hoogsteen hydrogen bond network, and description of the main chain dihedral angle. The results show that the overall performance of the Drude force field is the best. Although there may be a certain over-polarization effect, it is still the best choice for the MD simulation of G-quadruplexes.

Atomistic Polarizable Embeddings: Energy, Dynamics, Spectroscopy, and Reactivity

The computational modeling of realistic extended systems, relevant in, e.g., Chemistry and Biophysics, is a fundamental problem of paramount importance in contemporary research. Enzymatic catalysis and photoinduced processes in pigment-protein complexes are typical problems targeted by computer-aided approaches, to complement experiments as interpretative tools at a molecular scale. The daunting complexity of this task lies in between the opposite stringent requirements of results' reliability for structural/dynamical properties and related intermolecular interactions, and a mandatory principle of realism in the modeling strategy. Therefore, in practice, a truly realistic computational model of a biologically relevant system can easily fail to meet the accuracy requirement, in order to balance the excessive computational cost necessary to reach the desired precision.To address such an "accuracy vs reality" dualistic requirement, mixed quantum mechanics/classical mechanics approaches within Atomistic (i.e., preserving the discrete particle configuration) Polarizable Embeddings (QM/APEs) methods have been proposed over the years. In this Account, we review recent developments in the design and application of general QM/APE methods, targeting situations where a local intrinsically quantum behavior is coupled to a large molecular system (i.e., an environment), often involving processes with different dynamical time scales, in order to avoid brute-force, unpractical quantum chemistry calculations on the complete system.In the first place, our interest is devoted to the available APEs models presently implemented in computational software, highlighting the quantum chemistry methods that can be used to treat the QM subsystem. We review the coupling strategy between the QM subsystem and the APE, which requires to examine the way the QM/MM mutual interactions are accounted for and how the polarization of the classical environment is considered with respect to (wrt) the quantum variables. Because of the need of reliable molecular and macromolecular structures, a pivotal aspect to address here is the handling of the system dynamics (i.e., gradients wrt nuclear positions are required), especially for large molecular assemblies composed by an overwhelming number of atoms, exploring many conformations on a complex energy landscape.Alongside, we highlight our views on the necessary steps to take toward more accurate general-purposes and transferable explicit embeddings. The main objective to achieve here is to design a more physically grounded multiscale approach. To do so, one should apply advanced new generation classical models to account for refined induction effects that are able to (i) improve the quality of QM/MM interaction energies; (ii) enhance transferability by avoiding the compulsory partial (or total) reparameterization of the classical model. Moreover, the extension of recent developments originating from the field of advanced classical molecular dynamics (MD) to the realm of QM/APE methods is a key direction to improve both speed and efficiency for the phase space exploration of systems of growing size and complexity.Lastly, we point out specific research topics where an advanced QM/APE dynamics can certainly shed some light. For example, we discuss chemical reactions in "harsh" environments and the case of spectroscopic theoretical modeling where the inclusion of refined environment effects is often mandatory.

Nanomechanics of G-quadruplexes within the promoter of the KIT oncogene

G-quadruplexes (G4s) are tetrahelical DNA structures stabilized by four guanines paired via Hoogsteen hydrogen bonds into quartets. While their presence within eukaryotic DNA is known to play a key role in regulatory processes, their functional mechanisms are still under investigation. In the present work, we analysed the nanomechanical properties of three G4s present within the promoter of the KIT proto-oncogene from a single-molecule point of view through the use of magnetic tweezers (MTs). The study of DNA extension fluctuations under negative supercoiling allowed us to identify a characteristic fingerprint of G4 folding. We further analysed the energetic contribution of G4 to the double-strand denaturation process in the presence of negative supercoiling, and we observed a reduction in the energy required for strands separation.

Multimolecular complexes of the phosphodiester anion with Zn(II) or Mg(II) and water molecules-Preliminary validations of a polarizable potential by ab initio quantum chemistry

Dimethyl phosphate (DMP^- ) is a model for the phosphodiester backbone of DNA, RNA, and phospholipids. It is central for the binding of divalent cations and water along the backbone of nucleic acids. Significant polarization and charge-transfer contributions and nonadditivity come into play in the multimolecular complexes organized around phosphate. Prior to large-scale molecular dynamics (MD) with advanced polarizable potentials, it is essential to evaluate how well the values and trends of intermolecular interaction energies (ΔE) from ab initio quantum chemistry (QC) and their individual contributions are reproduced in a diversity of such complexes. These differ by the starting binding modes of a divalent cation, Zn(II), namely direct, bi- or mono-dentate to anionic and/or ester oxygens, versus through-water binding. We present first the results from automated refinements of the individual contributions of the SIBFA potential with respect to their QC counterparts using a Zn(II) or a water probe. This is followed by validations on eight relaxed multimolecular complexes of DMP^- with Zn(II) or Mg(II) and seven waters, then on sixteen complexes of DMP^- with Zn(II) and eight waters in arrangements extracted from MD or energy-minimization on a droplet of sixty-four waters. This monitors the compared evolutions of SIBFA and QC ΔE and their individual contributions in the competing arrangements. Some waters, bridging Zn(II) and DMP^- , were found to have exceptionally large dipole moments, of up to 3.8 Debye. The perspectives of extension to a flexible phosphodiester backbone are discussed in the context of the SIBFA potential for DNA and RNA.

Cation competition and recruitment around the c-kit1 G-quadruplex using polarizable simulations

Nucleic acid-ion interactions are fundamentally important to the physical, energetic, and conformational properties of DNA and RNA. These interactions help fold and stabilize highly ordered secondary and tertiary structures, such as G-quadruplexes (GQs), which are functionally relevant in telomeres, replication initiation sites, and promoter sequences. The c-kit proto-oncogene encodes for a receptor tyrosine kinase and is linked to gastrointestinal stromal tumors, mast cell disease, and leukemia. This gene contains three unique GQ-forming sequences that have proposed antagonistic effects on gene expression. The dominant GQ, denoted c-kit1, has been shown to decrease expression of c-kit transcripts, making the c-kit1 GQ a promising drug target. Toward disease intervention, more information is needed regarding its conformational dynamics and ion binding properties. Therefore, we performed molecular dynamics simulations of the c-kit1 GQ with K⁺, Na⁺, Li⁺, and mixed salt solutions using the Drude-2017 polarizable force field. We evaluated GQ structure, ion sampling, core energetics, ion dehydration and binding, and ion competition and found that each analysis supported the known GQ-ion specificity trend (K⁺ > Na⁺ > Li⁺). We also found that K⁺ ions coordinate in the tetrad core antiprismatically, whereas Na⁺ and Li⁺ align coplanar to guanine tetrads, partially because of their attraction to surrounding water. Further, we showed that K⁺ occupancy is higher around the c-kit1 GQ and its nucleobases than Na⁺ and Li⁺, which tend to interact with backbone and sugar moieties. Finally, we showed that K⁺ binding to the c-kit1 GQ is faster and more frequent than Na⁺ and Li⁺. Such descriptions of GQ-ion dynamics suggest the rate of dehydration as the dominant factor for preference of K⁺ by DNA GQs and provide insight into noncanonical nucleic acids for which little experimental data exist.

G-Quadruplex Modulation of SP1 Functional Binding Sites at the KIT Proximal Promoter

The regulation of conformational arrangements of gene promoters is a physiological mechanism that has been associated with the fine control of gene expression. Indeed, it can drive the time and the location for the selective recruitment of proteins of the transcriptional machinery. Here, we address this issue at the KIT proximal promoter where three G-quadruplex forming sites are present (kit1, kit2 and kit*). On this model, we focused on the interplay between G-quadruplex (G4) formation and SP1 recruitment. By site directed mutagenesis, we prepared a library of plasmids containing mutated sequences of the WT KIT promoter that systematically exploited different G4 formation attitudes and SP1 binding properties. Our transfection data showed that the three different G4 sites of the KIT promoter impact on SP1 binding and protein expression at different levels. Notably, kit2 and kit* structural features represent an on-off system for KIT expression through the recruitment of transcription factors. The use of two G4 binders further helps to address kit2-kit* as a reliable target for pharmacological intervention.

Ion Binding Properties and Dynamics of the bcl-2 G-Quadruplex Using a Polarizable Force Field

G-quadruplexes (GQs) are topologically diverse, highly thermostable noncanonical nucleic acid structures that form in guanine-rich sequences in DNA and RNA. GQs are implicated in transcriptional and translational regulation and genome maintenance, and deleterious alterations to their structures contribute to diseases such as cancer. The expression of the B-cell lymphoma 2 (Bcl-2) antiapoptotic protein, for example, is under transcriptional control of a GQ in the promoter of the bcl-2 gene. Modulation of the bcl-2 GQ by small molecules is of interest for chemotherapeutic development but doing so requires knowledge of the factors driving GQ folding and stabilization. To develop a greater understanding of the electrostatic properties of the bcl-2 promoter GQ, we performed molecular dynamics simulations using the Drude-2017 polarizable force field and compared relevant outcomes to the nonpolarizable CHARMM36 force field. Our simulation outcomes highlight the importance of dipole-dipole interactions in the bcl-2 GQ, particularly during the recruitment of a bulk K⁺ ion to the solvent-exposed face of the tetrad stem. We also predict and characterize an "electronegative pocket" at the tetrad-long loop junction that induces local backbone conformational change and may induce local conformational changes at cellular concentrations of K⁺. These outcomes suggest that moieties within the bcl-2 GQ can be targeted by small molecules to modulate bcl-2 GQ stability.

Dynamics and Barrier of Movements of Sodium and Potassium Ions Across the Oxytricha nova G-Quadruplex Core

G-quadruplexes (GQs) are highly stable noncanonical forms of nucleic acids that are present in important genomic regions. The central core of the GQ is lined up by four closely spaced carbonyl groups from the G-quartets, and the resulting electrostatic repulsion is neutralized by the coordinating cations. In spite of several reports on GQ structure and cation-GQ interactions, the atomic- to molecular-level understanding of the ion dynamics and ion exchange in the GQ core is quite poor. Here, we attempt to elucidate the mechanism of Na⁺ and K⁺ binding to the GQ core and trace the exchange of these ions with the ions in bulk by means of all-atomic molecular dynamics (MD) simulations. One of the most studied GQs, Oxytricha nova telomeric G-quadruplex (OxyGQ), is taken as the representative GQ. Subsequently, umbrella sampling MD simulations were performed to elucidate the energetics of ion translocation from one end to the other end of the GQ central core. Our study highlights the importance of ion hydration for the uptake and correct positioning of the cations in the core. The free-energy landscape of ion transport has shown favorable in-plane binding of Na⁺ ions with GQ quartets, which matches very well with the crystal structure. The binding of K⁺ ions, on the other hand, was out-of-plane and its translocation required a larger barrier to cross.

Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics

Nucleic acids play critical roles in carrying genetic information, participating in catalysis, and preserving chromosomal structure. Despite over a century of study, efforts to understand the dynamics and structure-function relationships of DNA and RNA at the atomic level are still ongoing. Molecular dynamics (MD) simulations augment experiments by providing atomistic resolution and quantitative relationships between structure and conformational energy. Steady advancements in computer hardware, software, and atomistic force fields (FFs) over 40 years have facilitated new discoveries. Here, we review nucleic acid FF development with emphasis on recent refinements that have improved descriptions of important nucleic acid properties. We then discuss several key examples of successes and challenges in modeling nucleic acid structure and dynamics using the latest FFs.