Studying the Dynamics of a Complex G-Quadruplex System: Insights into the Comparison of MD and NMR Data

AMOEBA Polarizable Molecular Dynamics Simulations of Guanine Quadruplexes: From the c-Kit Proto-Oncogene to HIV-1

Long oligomer sequences, rich in guanine and cytosine, such as c-kit1 and the HIV-1 LTR-III sequence, are prevalent in oncogenes and retroviruses and play crucial roles in cancer. Understanding the conformational dynamics of such guanine quadruplexes and identifying druggable regions are therefore essential for developing new inhibition strategies. In this study, we used extensive AMOEBA polarizable force field molecular dynamics simulations combined with data-driven adaptive sampling and clustering algorithms, reaching a cumulative simulation time of 7.5 μs for c-kit1. Such simulations identified novel structural motives and showcased the flexible loop dynamics, as well as the role of polarizable water in transient stabilization of the G-quadruplex. They also identified two druggable pockets in c-kit1. The 400 ns simulation of the HIV-1 LTR-III sequence confirmed its quadruplex stability and uncovered a potentially druggable cryptic pocket.

Unveiling a novel RNA G-triplex structure: its function and potential in CRISPR-based diagnostics

We report the discovery of a novel higher-order RNA structure, RNA G-triplex (rG3), formed by the TERRA sequence. Through CD spectroscopy, NMR analysis, and molecular modeling, we confirmed its stable, parallel conformation. rG3 exhibits strong binding to thioflavin T (ThT), N-methyl mesoporphyrin IX (NMM), and hemin, showcasing its potential as a biosensing element. Additionally, CRISPR-Cas13a trans-cleaves rG3, demonstrating its utility as a sensitive reporter in diagnostic applications. These findings expand the structural diversity of RNA and suggest new avenues for RNA-based biosensors and CRISPR diagnostics.

Exploring the interaction between a fluorescent Ag(I)-biscarbene complex and non-canonical DNA structures: a multi-technique investigation

Silver compounds are mainly studied as antimicrobial agents, but they also have anticancer properties, with the latter, in some cases, being better than their gold counterparts. Herein, we analyse the first example of a new Ag(I)-biscarbene that can bind non-canonical structures of DNA, more precisely G-quadruplexes (G4), with different binding signatures depending on the type of G4. Moreover, we show that this Ag-based carbene binds the i-motif DNA structure. Alternatively, its Au(I) counterpart, which was investigated for comparison, stabilises mitochondrial G4. Theoretical in silico studies elucidated the details of different binding modes depending on the geometry of G4. The two complexes showed increased cytotoxic activity compared to cisplatin, overcoming its resistance in ovarian cancer. The binding of these new drug candidates with other relevant biosubstrates was studied to afford a more complete picture of their possible targets. In particular, the Ag(I) complex preferentially binds DNA structures over RNA structures, with higher binding constants for the non-canonical nucleic acids with respect to natural calf thymus DNA. Regarding possible protein targets, its interaction with the albumin model protein BSA was also tested.

From RNA sequence to its three-dimensional structure: geometrical structure, stability and dynamics of selected fragments of SARS-CoV-2 RNA

In this computational study, we explore the folding of a particular sequence using various computational tools to produce two-dimensional structures, which are then transformed into three-dimensional structures. We then study the geometry, energetics and dynamics of these structures using full electron quantum-chemical and classical molecular dynamics calculations. Our study focuses on the SARS-CoV-2 RNA fragment GGaGGaGGuguugcaGG and its various structures, including a G-quadruplex and five different hairpins. We examine the impact of two types of counterions (K+ and Na+) and flanking nucleotides on their geometrical characteristics, relative stability and dynamic properties. Our results show that the G-quadruplex structure is the most stable among the constructed hairpins. We confirm its topological stability through molecular dynamics simulations. Furthermore, we observe that the nucleotide loop consisting of seven nucleotides is the most flexible part of the RNA fragment. Additionally, we find that RNA networks of intermolecular hydrogen bonds are highly sensitive to the surrounding environment. Our findings reveal the loss of 79 old hydrogen bonds and the formation of 91 new ones in the case when the G-quadruplex containing flanking nucleotides is additionally stabilized by Na+ counterions.

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.

Unraveling the structural aspects of the G-quadruplex in SMO promoter and elucidating its contribution in transcriptional regulation

We located a 25 nt G-rich sequence in the promoter region of SMO oncogene. We performed an array of biophysical and biochemical assays and confirmed the formation of a parallel G quadruplex (SMO1-GQ) by the identified sequence. SMO1-GQ is highly conserved in primates. For a comprehensive characterization of the SMO quadruplex structure, we have performed spectroscopic and in silico analysis with established GQ binder small molecules TMPyP4 and BRACO-19. We observed comparatively higher stable interaction of BRACO-19 with SMO1-GQ. Structure-based, rational drug design against SMO1-GQ to target SMO oncogene requires a detailed molecular anatomy of the G-quadruplex. We structurally characterised the SMO1-GQ using DMS footprinting assay and molecular modelling, docking, and MD simulation to identify the probable atomic regions that interact with either of the small molecules. We further investigated SMO1-GQ in vivo by performing chromatin immunoprecipitation (ChIP) assay. ChIP data revealed that this gene element functions as a scaffold for a number of transcription factors: specificity protein (Sp1), nucleolin (NCL), non-metastatic cell 2 (NM23-H2), cellular nucleic acid binding protein (CNBP), and heterogeneous nuclear ribonucleoprotein K (hnRNPK) which reflects the SMO1-P1 G-quadruplex to be the master regulator of SMO1 transcriptional activity. The strong binding interaction detected between SMO1-GQ and BRACO-19 contemplates the potential of the G quadruplex as a promising anti-cancer druggable target to downregulate SMO1 oncogene driven cancers.Communicated by Ramaswamy H. Sarma.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Polarizable Molecular Block Model: Toward the Development of an Induced Dipole Force Field for DNA

For flexible and highly ionized macromolecules such as DNA, it is important to correctly evaluate the intramolecular polarization in an induced dipole force field. In a proposed polarizable molecular block (PMB) model, a large molecule is divided into several molecular blocks. The atomic charges of the blocks are optimized by using the respective electrostatic potentials (ESPs) on the molecular surface. By using the capped hydrogen removal operation, the total charge of the blocks is controlled exactly to have an integer charge. The atomic polarizabilities of the blocks are optimized by using the respective polarized one-electron potentials that are the differences between ESPs with and without an external test charge. Induced dipole-charge interactions between the blocks are all included, but those interactions within the blocks are strictly excluded. All dipole-dipole interactions are included, but the damping functions are applied to the close dipole-dipole pairs. Several types of damping (simple scaling, exponential, linear, and Gaussian) are evaluated. The validity of the PMB model was verified by using trinucleotide duplexes which have A-, B-, and Z-DNA forms. The reference energies of trinucleotide duplexes including counterions (GGT3Na-ACC3Na, GAC3Na-GTC3Na, and GCG3Na-CGC3Na) are calculated using ωB97XD/aug-cc-pVDZ. All damping types reproduced well the reference interaction energies, dipole moments, and ESPs. Among them, the simple scaling with strong attenuation to 1-2 atomic pairs showed the highest stability against the polarization catastrophe. This study shows that it is possible to develop a high-quality polarizable force field by treating the intramolecular polarization on a block-by-block basis.

Enhancement of DNAzymatic activity using iterative in silico maturation

Enhancement of DNZymatic activity using a combined iterative in silico and in vitro method as a cheaper and more stable alternative to antibodies or enzymes.