The 3D structures of G-Quadruplexes of HIV-1 integrase inhibitors: molecular dynamics simulations in aqueous solution and in the gas phase

Enhancement of the thermal stability of G-quadruplex structures by urea

The solute urea has been used extensively as a denaturant in protein folding studies; double-stranded nucleic acid structures are also destabilized by urea, but comparatively less than proteins. In previous research, the solute has been shown to strongly destabilize folded G-quadruplex DNA structures. This contribution demonstrates the stabilizing effect of urea on the G-quadruplex formed by the oligodeoxyribonucleotide (ODN), G3T (d[5'-GGGTGGGTGGGTGGG-3']), and related sequences in the presence of sodium or potassium cations. Stabilization is observed up to 7 M urea, which was the highest concentration we investigated. The folded structure of G3T has three G-tetrads and three loops that consist of single thymine residues. ODNs related to G3T, in which the thymine residues in the loop are substituted by adenosine residues, also exhibit enhanced stability in the presence of molar concentrations of urea. The circular dichroism (CD) spectra of these ODNs in the presence of urea are consistent with that of a G-quadruplex. As the urea concentration increases, the spectral intensities of the peaks and troughs change, while their positions change very little. The heat-induced transition from the folded to unfolded state, T_m, was measured by monitoring the change in the UV absorption as a function of temperature. G-quadruplex structures with loops containing single bases exhibited large increases in T_m with increasing urea concentrations. These data imply that the loop region play a significant role in the thermal stability of tetra-helical DNA structures in the presence of the solute urea.

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.

Three-Dimensional Structure of RNA Monomeric G-Quadruplex Containing ALS and FTD Related G4C2 Repeat and Its Binding with TMPyP4 Probed by Homology Modeling based on Experimental Constraints and Molecular Dynamics Simulations

The G-quadruplex-forming hexanucleotide repeat expansion (HRE), d(G4C2)_n, within the human C9orf72 gene is the root cause for familial amyotrophic lateral sclerosis-frontotemporal dementia (ALS-FTD). A recent study has shown that TMPyP4 has good potential to work as a RNA G-quadruplex binder in treating ALS and FTD. Although the high-resolution structure of the monomeric DNA antiparallel G-quadruplex form of the monomeric hexanucleotide repeat was recently solved, the RNA parallel G-quadruplex structure and its complex with TMPyP4 are not available yet. In this study, we first constructed the homology model for the parallel monomeric RNA G-quadruplex of r(G4C2)₃G4 based on experimental constraints and the parallel monomeric G-quadruplex DNA crystal structure. Although the G-tetra core of the homology model was stable observed in 15 μs molecular dynamics (MD) simulations, we observed that the loops adopt additional conformations besides the initial crystal conformation, where TMPyP4 binding was found to reduce the loop fluctuation of the RNA monomeric G-quadruplex. Next, we probed the elusive binding behavior of TMPyP4 to the RNA monomeric G-quadruplex. Encouragingly, the binding modes observed are similar to the modes observed in two experimental complexes of a parallel DNA G-quadruplex with TMPyP4. We also constructed a Markov state model to provide insights into the binding pathways. Together, the findings from our study may assist future development of G-quadruplex-specific ligands in the treatment of neurodegenerative diseases like ALS and FTD.

Computational Approaches to Prioritize Cancer Driver Missense Mutations

Cancer is a complex disease that is driven by genetic alterations. There has been a rapid development of genome-wide techniques during the last decade along with a significant lowering of the cost of gene sequencing, which has generated widely available cancer genomic data. However, the interpretation of genomic data and the prediction of the association of genetic variations with cancer and disease phenotypes still requires significant improvement. Missense mutations, which can render proteins non-functional and provide a selective growth advantage to cancer cells, are frequently detected in cancer. Effects caused by missense mutations can be pinpointed by in silico modeling, which makes it more feasible to find a treatment and reverse the effect. Specific human phenotypes are largely determined by stability, activity, and interactions between proteins and other biomolecules that work together to execute specific cellular functions. Therefore, analysis of missense mutations’ effects on proteins and their complexes would provide important clues for identifying functionally important missense mutations, understanding the molecular mechanisms of cancer progression and facilitating treatment and prevention. Herein, we summarize the major computational approaches and tools that provide not only the classification of missense mutations as cancer drivers or passengers but also the molecular mechanisms induced by driver mutations. This review focuses on the discussion of annotation and prediction methods based on structural and biophysical data, analysis of somatic cancer missense mutations in 3D structures of proteins and their complexes, predictions of the effects of missense mutations on protein stability, protein-protein and protein-nucleic acid interactions, and assessment of conformational changes in protein conformations induced by mutations.

Anti-HIV Activities of Intramolecular G4 and Non-G4 Oligonucleotides

New natural and chemically modified DNA aptamers that inhibit HIV-1 activity at submicromolar concentrations (presumably via preventing viral entry into target cells) are reported. The new DNA aptamers were developed based on known intramolecular G-quadruplexes (G4s) that were functionally unrelated to HIV inhibition [the thrombin-binding aptamer and the fragment of the human oncogene promoter (Bcl2)]. The majority of previously described DNA inhibitors of HIV infection adopt intermolecular structures, and thus their folding variability represents an obvious disadvantage. Intramolecular architectures refold correctly after denaturation and are generally easier to handle. However, whether the G4 topology or other factors account for the anti-HIV activity of our aptamers is unknown. The impact of chemical modification (thiophosphoryl internucleotide linkages) on aptamer activity is discussed. The exact secondary structures of the active compounds and further elucidation of their mechanisms of action hopefully will be the subjects of future studies.

Molecular Dynamics Simulation Study of Parallel Telomeric DNA Quadruplexes at Different Ionic Strengths: Evaluation of Water and Ion Models

Most molecular dynamics (MD) simulations of DNA quadruplexes have been performed under minimal salt conditions using the Åqvist potential parameters for the cation with the TIP3P water model. Recently, this combination of parameters has been reported to be problematic for the stability of quadruplex DNA, especially caused by the ion interactions inside or near the quadruplex channel. Here, we verify how the choice of ion parameters and water model can affect the quadruplex structural stability and the interactions with the ions outside the channel. We have performed a series of MD simulations of the human full-parallel telomeric quadruplex by neutralizing its negative charge with K(+) ions. Three combinations of different cation potential parameters and water models have been used: (a) Åqvist ion parameters, TIP3P water model; (b) Joung and Cheatham ion parameters, TIP3P water model; and (c) Joung and Cheatham ion parameters, TIP4Pew water model. For the combinations (b) and (c), the effect of the ionic strength has been evaluated by adding increasing amounts of KCl salt (50, 100, and 200 mM). Two independent simulations using the Åqvist parameters with the TIP3P model show that this combination is clearly less suited for the studied quadruplex with K(+) as counterions. In both simulations, one ion escapes from the channel, followed by significant deformation of the structure, leading to deviating conformation compared to that in the reference crystallographic data. For the other combinations of ion and water potentials, no tendency is observed for the channel ions to escape from the quadruplex channel. In addition, the internal mobility of the three loops, torsion angles, and counterion affinity have been investigated at varied salt concentrations. In summary, the selection of ion and water models is crucial as it can affect both the structure and dynamics as well as the interactions of the quadruplex with its counterions. The results obtained with the TIP4Pew model are found to be closest to the experimental data at all of the studied ion concentrations.

Understanding the conformation transition in the activation pathway of β2 adrenergic receptor via a targeted molecular dynamics simulation

The conformation transition in the activation pathway of β2 adrenergic receptor was explored mainly using a target molecular dynamics simulation.

MAZ-binding G4-decoy with locked nucleic acid and twisted intercalating nucleic acid modifications suppresses KRAS in pancreatic cancer cells and delays tumor growth in mice

KRAS mutations are primary genetic lesions leading to pancreatic cancer. The promoter of human KRAS contains a nuclease-hypersensitive element (NHE) that can fold in G4-DNA structures binding to nuclear proteins, including MAZ (myc-associated zinc-finger). Here, we report that MAZ activates KRAS transcription. To knockdown oncogenic KRAS in pancreatic cancer cells, we designed oligonucleotides that mimic one of the G-quadruplexes formed by NHE (G4-decoys). To increase their nuclease resistance, two locked nucleic acid (LNA) modifications were introduced at the 3'-end, whereas to enhance the folding and stability, two polycyclic aromatic hydrocarbon units (TINA or AMANY) were inserted internally, to cap the quadruplex. The most active G4-decoy (2998), which had two para-TINAs, strongly suppressed KRAS expression in Panc-1 cells. It also repressed their metabolic activity (IC50 = 520 nM), and it inhibited cell growth and colony formation by activating apoptosis. We finally injected 2998 and control oligonucleotides 5153, 5154 (2 nmol/mouse) intratumorally in SCID mice bearing a Panc-1 xenograft. After three treatments, 2998 reduced tumor xenograft growth by 64% compared with control and increased the Kaplan-Meier median survival time by 70%. Together, our data show that MAZ-specific G4-decoys mimicking a KRAS quadruplex are promising for pancreatic cancer therapy.

d(GGGT) 4 and r(GGGU) 4 are both HIV-1 inhibitors and interleukin-6 receptor aptamers

Aptamers are oligonucleotides that bind targets with high specificity and affinity. They have become important tools for biosensing, target detection, drug delivery and therapy. We selected the quadruplex-forming 16-mer DNA aptamer AID-1 [d(GGGT) 4] with affinity for the interleukin-6 receptor (IL-6R) and identified single nucleotide variants that showed no significant loss of binding ability. The RNA counterpart of AID-1 [r(GGGU) 4] also bound IL-6R as quadruplex structure. AID-1 is identical to the well-known HIV inhibitor T30923, which inhibits both HIV infection and HIV-1 integrase. We also demonstrated that IL-6R specific RNA aptamers not only bind HIV-1 integrase and inhibit its 3' processing activity in vitro, but also are capable of preventing HIV de novo infection with the same efficacy as the established inhibitor T30175. All these aptamer target interactions are highly dependent on formation of quadruplex structure.

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.

Cation binding to 15-TBA quadruplex DNA is a multiple-pathway cation-dependent process

A combination of explicit solvent molecular dynamics simulation (30 simulations reaching 4 µs in total), hybrid quantum mechanics/molecular mechanics approach and isothermal titration calorimetry was used to investigate the atomistic picture of ion binding to 15-mer thrombin-binding quadruplex DNA (G-DNA) aptamer. Binding of ions to G-DNA is complex multiple pathway process, which is strongly affected by the type of the cation. The individual ion-binding events are substantially modulated by the connecting loops of the aptamer, which play several roles. They stabilize the molecule during time periods when the bound ions are not present, they modulate the route of the ion into the stem and they also stabilize the internal ions by closing the gates through which the ions enter the quadruplex. Using our extensive simulations, we for the first time observed full spontaneous exchange of internal cation between quadruplex molecule and bulk solvent at atomistic resolution. The simulation suggests that expulsion of the internally bound ion is correlated with initial binding of the incoming ion. The incoming ion then readily replaces the bound ion while minimizing any destabilization of the solute molecule during the exchange.

Stacking of G-quadruplexes: NMR structure of a G-rich oligonucleotide with potential anti-HIV and anticancer activity

G-rich oligonucleotides T30695 (or T30923), with the sequence of (GGGT)₄, and T40214, with the sequence of (GGGC)₄, have been reported to exhibit anti-HIV and anticancer activity. Here we report on the structure of a dimeric G-quadruplex adopted by a derivative of these sequences in K⁺ solution. It comprises two identical propeller-type parallel-stranded G-quadruplex subunits each containing three G-tetrad layers that are stacked via the 5′-5′ interface. We demonstrated control over the stacking of the two monomeric subunits by sequence modifications. Our analysis of possible structures at the stacking interface provides a general principle for stacking of G-quadruplexes, which could have implications for the assembly and recognition of higher-order G-quadruplex structures.

HIV-1 integrase inhibitor T30177 forms a stacked dimeric G-quadruplex structure containing bulges

T30177 is a G-rich oligonucleotide with the sequence (GTGGTGGGTGGGTGGGT) which inhibits the HIV-1 integrase activity at nanomolar concentrations. Here we show that this DNA sequence forms in K(+) solution a dimeric G-quadruplex structure comprising a total of six G-tetrad layers through the stacking of two propeller-type parallel-stranded G-quadruplex subunits at their 5'-end. All twelve guanines in the sequence participate in the G-tetrad formation, despite the interruption in the first G-tract by a thymine, which forms a bulge between two adjacent G-tetrads. In this work, we also propose a simple analytical approach to stoichiometry determination using concentration-dependent melting curves.