Resolving the Ligand-Binding Specificity in c-MYC G-Quadruplex DNA: Absolute Binding Free Energy Calculations and SPR Experiment

Applications of Surface Plasmon Resonance for Advanced Studies Involving Nucleic Acids

Surface plasmon resonance (SPR) is increasingly recognized as one of the most widely used techniques for studying nucleic acid interactions. The main advantage of SPR is its ability to measure the binding affinities and association/dissociation kinetics of complexes in real-time, in a label-free environment, and using relatively small quantities of materials. The method is based on the immobilization of one of the binding partners, ligand, on a dedicated sensor surface. Immobilization is followed by the injection of the other partner, analyte, over the surface containing the ligand. The binding is monitored by subsequent changes in the refractive index of the medium close to the sensor surface upon injection of the analyte. In the field of Nucleic Acid, SPR has been intensively used in the study of various artificial and naturally occurring RNA/DNA molecules interaction with large molecular weight mass proteins and small organic molecules because of its ability to detect highly dynamic complexes that are difficult to investigate using other techniques. This mini review aims to provide a short guideline for setting up SPR experiments to identify nucleic acid complexes and assess their binding affinity or kinetics. It covers protocols for (i) nucleic acid immobilization methods, including biotin-streptavidin, metal ion-based affinity, and amine coupling, (ii) analyte-binding analysis, (iii) affinity and kinetic measurements, and (iv) data interpretation. Determining the affinity and kinetics of nucleic acid interactions through SPR is essential for gaining insights into molecular-level binding mechanisms, thus supporting advancements in nucleic acid nanotechnology. The review also highlights the various sections of SPR applications in nucleic acid research, including nucleic acid-probe immobilization, interactions with biomolecules, aptamer studies, and small molecule binding, concluding with perspectives on future developments in the field.

Modeling and NMR Data Elucidate the Structure of a G-Quadruplex-Ligand Interaction for a Pu22T-Cyclometalated Iridium(III) System

Cyclometalated iridium(III) complexes are increasingly being developed for application in G-quadruplex (GQ) nucleic acid biosensors. We monitored the interactions of a GQ structure with an iridium(III) complex by nuclear magnetic resonance (NMR) titrations and subsequently compared the binding site inferred from NMR with binding positions modeled by molecular docking and molecular dynamics simulations. When titrated into a solution of G-quadruplex Pu22T, compound 1(PF6), [Ir(ppy)2(pizp)](PF6), where ppy is 2-phenylpyridine and pizp is 2-phenylimidazole[4,5f][1,10]phenanthroline, had the greatest impact on the hydrogen chemical shifts of G5, G8, G9, G13, and G17 residues of Pu22T, indicating end-stacking at the 5' tetrad. In blind cross-docking studies with Autodock 4, end-stacking at the 5' tetrad was found as the lowest energy binding position. AMBER molecular dynamics simulations resulted in a refined binding position at the 5' tetrad with improved pi stacking. For this model system, Pu22T-1, molecular docking and molecular dynamics simulations are tools that are able to predict the experimentally determined binding position.

Mechanistic Insights into the c-MYC G-Quadruplex and Berberine Binding inside an Aqueous Two-Phase System Mimicking Biomolecular Condensates

We investigated the binding between the c-MYC G-quadruplex (GQ) and berberine chloride (BCl) in an aqueous two-phase system (ATPS) with 12.3 wt % polyethylene glycol and 5.6 wt % dextran, mimicking the highly crowded intracellular biomolecular condensates formed via liquid-liquid phase separation. We found that in the ATPS, complex formation is significantly altered, leading to an increase in affinity and a change in the stoichiometry of the complex with respect to neat buffer conditions. Thermodynamic studies reveal that binding becomes more thermodynamically favorable in the ATPS due to entropic effects, as the strong excluded volume effect inside ATPS droplets reduces the entropic penalty associated with binding. Finally, the binding affinity of BCl for the c-MYC GQ is higher than those for other DNA structures, indicating potential specific interactions. Overall, these findings will be helpful in the design of potential drugs targeting the c-MYC GQ structures in cancer-related biocondensates.

Application and Method of Surface Plasmon Resonance Technology in the Preparation and Characterization of Biomedical Nanoparticle Materials

Surface Plasmon Resonance (SPR) technology, as a powerful analytical tool, plays a crucial role in the preparation, performance evaluation, and biomedical applications of nanoparticles due to its real-time, label-free, and highly sensitive detection capabilities. In the nanoparticle preparation process, SPR technology can monitor synthesis reactions and surface modifications in real-time, optimizing preparation techniques and conditions. SPR enables precise measurement of interactions between nanoparticles and biomolecules, including binding affinities and kinetic parameters, thereby assessing nanoparticle performance. In biomedical applications, SPR technology is extensively used in the study of drug delivery systems, biomarker detection for disease diagnosis, and nanoparticle-biomolecule interactions. This paper reviews the latest advancements in SPR technology for nanoparticle preparation, performance evaluation, and biomedical applications, discussing its advantages and challenges in biomedical applications, and forecasting future development directions.

Exploring the Dispersion and Electrostatic Components in Arene-Arene Interactions between Ligands and G4 DNA to Develop G4-Ligands

G-Quadruplex (G4) DNA structures are important regulatory elements in central biological processes. Small molecules that selectively bind and stabilize G4 structures have therapeutic potential, and there are currently >1000 known G4 ligands. Despite this, only two G4 ligands ever made it to clinical trials. In this work, we synthesized several heterocyclic G4 ligands and studied their interactions with G4s (e.g., G4s from the c-MYC, c-KIT, and BCL-2 promoters) using biochemical assays. We further studied the effect of selected compounds on cell viability, the effect on the number of G4s in cells, and their pharmacokinetic properties. This identified potent G4 ligands with suitable properties and further revealed that the dispersion component in arene-arene interactions in combination with electron-deficient electrostatics is central for the ligand to bind with the G4 efficiently. The presented design strategy can be applied in the further development of G4-ligands with suitable properties to explore G4s as therapeutic targets.

Towards Computational Modeling of Ligand Binding to the ILPR G-Quadruplex

Using a combination of unconstrained and constrained molecular dynamics simulations, we have evaluated the binding affinities between two porphyrin derivatives (TMPyP4 and TEGPy) and the G-quadruplex (G4) of a DNA fragment modeling the insulin-linked polymorphic region (ILPR). Refining a well-established potential of mean force (PMF) approach to selections of constraints based on root-mean-square fluctuations results in an excellent agreement between the calculated and observed absolute free binding energy of TMPyP4. The binding affinity of IPLR-G4 toward TEGPy is predicted to be higher than that toward TMPyP4 by 2.5 kcal/mol, which can be traced back to stabilization provided by the polyether side chains of TMPyP4 that can nestle into the grooves of the quadruplex and form hydrogen bonds through the ether oxygen atoms. Because our refined methodology can be applied to large ligands with high flexibility, the present research opens an avenue for further ligand design in this important area.

The effect of side chain variations on quinazoline-pyrimidine G-quadruplex DNA ligands

G-quadruplex (G4) DNA structures are involved in central biological processes such as DNA replication and transcription. These DNA structures are enriched in promotor regions of oncogenes and are thus promising as novel gene silencing therapeutic targets that can be used to regulate expression of oncoproteins and in particular those that has proven hard to drug with conventional strategies. G4 DNA structures in general have a well-defined and hydrophobic binding area that also is very flat and featureless and there are ample examples of G4 ligands but their further progression towards drug development is limited. In this study, we use synthetic organic chemistry to equip a drug-like and low molecular weight central fragment with different side chains and evaluate how this affect the compound's selectivity and ability to bind and stabilize G4 DNA. Furthermore, we study the binding interactions of the compounds and connect the experimental observations with the compound's structural conformations and electrostatic potentials to understand the basis for the observed improvements. Finally, we evaluate the top candidates' ability to selectively reduce cancer cell growth in a 3D co-culture model of pancreatic cancer which show that this is a powerful approach to generate highly active and selective low molecular weight G4 ligands with a promising therapeutic window.

Using Macrocyclic G-Quadruplex Ligands to Decipher the Interactions Between Small Molecules and G-Quadruplex DNA

This study aims to deepen the knowledge of the current state of rational G4-ligand design through the design and synthesis of a novel set of compounds based on indoles, quinolines, and benzofurans and their comparisons with well-known G4-ligands. This resulted in novel synthetic methods and G4-ligands that bind and stabilize G4 DNA with high selectivity. Furthermore, the study corroborates previous studies on the design of G4-ligands and adds deeper explanations to why a) macrocycles offer advantages in terms of G4-binding and -selectivity, b) molecular pre-organization is of key importance in the development of strong novel binders, c) an electron-deficient aromatic core is essential to engage in strong arene-arene interactions with the G4-surface, and d) aliphatic amines can strengthen interactions indirectly through changing the arene electrostatic nature of the compound. Finally, fundamental physicochemical properties of selected G4-binders are evaluated, underscoring the complexity of aligning the properties required for efficient G4 binding and stabilization with feasible pharmacokinetic properties.

DNA G-Quadruplex in Human Telomeres and Oncogene Promoters: Structures, Functions, and Small Molecule Targeting

DNA G-quadruplex secondary structures formed in guanine-rich human telomeres and oncogene promoters are functionally important and have emerged as a promising new class of cancer-specific drug targets. These globular intramolecular structures are stabilized by K+ or Na+ and form readily under physiological solution conditions. Moreover, G-quadruplexes are epigenetic features and can alter chromatin structure and function together with interactive proteins. Here, we discuss our efforts over the last two decades to understand the structures and functions of DNA G-quadruplexes formed in key oncogene promoters and human telomeres and their interactions with small molecules. Using high-field NMR spectroscopy, we determined the high-resolution structures of physiologically relevant telomeric G-quadruplexes in K+ solution with a major form (hybrid-2) and a minor form (hybrid-1), as well as a two-tetrad intermediate. The intrinsic structural polymorphism of telomeric DNA may be important for the biology of human telomeres, and we proposed a model for the interconversion. More recently, we have worked on G-quadruplexes of MYC, BCL2, PDGFR-β, VEGF, and k-RAS oncogene promoters. We determined the structure of the major G-quadruplex formed in the MYC promoter, a prototype for parallel G-quadruplexes. It is the first example of the parallel-stranded G3NG3 structure motif with a 1-nt loop, which is prevalent in promoter sequences and likely evolutionarily selected to initiate folding. Remarkably, the parallel MYC promoter G-quadruplexes are highly stable. Additionally, we determined the molecular structures of G-quadruplexes formed in human BCL2, VEGF, and PDGFR-β promoters, each adopting a unique structure. For example, the BCL2 promoter contains distinct interchangeable G-quadruplexes in two adjacent regions, suggesting precise regulation by different proteins. The PDGFR-β promoter adopts unique "broken-strand" and vacancy G-quadruplexes, which can be recognized by cellular guanine metabolites for a potential regulatory role.Structural information on G-quadruplexes in complex with small-molecules is critical for understanding specific recognition and structure-based rational drug design. Our studies show that many G-quadruplexes contain unique structural features such as capping and loop structures, allowing specific recognition by drugs and protein. This represents a paradigm shift in understanding DNA as a drug target: Rather than a uniform, nonselective binding site in duplex DNA, the G-quadruplex is being pursued as a new class of selectively targetable drug receptors. We focus on targeting the biologically relevant MYC promoter G-quadruplex (MycG4) with small molecules and have determined its first and additional drug complex structures. Very recently, we have discovered clinically tested indenoisoquinolines as strong MycG4 binders and potent MYC inhibitors. We have also discovered drugs targeting the unique dGMP-bound-vG4 formed in the PDGFR-β promoter. Moreover, we determined the complex structures of the first small molecules that specifically recognize the physiologically relevant human telomeric G-quadruplexes. Unlike the previously recognized dogma that the optimal G-quadruplex ligands are large aromatic or cyclic compounds, our results suggest that smaller asymmetric compounds with appropriate functional groups are better choices to specifically bind G-quadruplexes. This body of work lays a strong foundation for future work aimed at understanding the cellular functions of G-quadruplexes and G-quadruplex-targeted drug design.

Structure-based virtual screening workflow to identify antivirals targeting HIV-1 capsid

We have identified novel HIV-1 capsid inhibitors targeting the PF74 binding site. Acting as the building block of the HIV-1 capsid core, the HIV-1 capsid protein plays an important role in the viral life cycle and is an attractive target for antiviral development. A structure-based virtual screening workflow for hit identification was employed, which includes docking 1.6 million commercially-available drug-like compounds from the ZINC database to the capsid dimer, followed by applying two absolute binding free energy (ABFE) filters on the 500 top-ranked molecules from docking. The first employs the Binding Energy Distribution Analysis Method (BEDAM) in implicit solvent. The top-ranked compounds are then refined using the Double Decoupling method in explicit solvent. Both docking and BEDAM refinement were carried out on the IBM World Community Grid as part of the FightAIDS@Home project. Using this virtual screening workflow, we identified 24 molecules with calculated binding free energies between − 6 and − 12 kcal/mol. We performed thermal shift assays on these molecules to examine their potential effects on the stability of HIV-1 capsid hexamer and found that two compounds, ZINC520357473 and ZINC4119064 increased the melting point of the latter by 14.8 °C and 33 °C, respectively. These results support the conclusion that the two ZINC compounds are primary hits targeting the capsid dimer interface. Our simulations also suggest that the two hit molecules may bind at the capsid dimer interface by occupying a new sub-pocket that has not been exploited by existing CA inhibitors. The possible causes for why other top-scored compounds suggested by ABFE filters failed to show measurable activity are discussed.

Developing end-point methods for absolute binding free energy calculation using the Boltzmann-quasiharmonic model

Understanding the physical forces underlying receptor-ligand binding requires robust methods for analyzing the binding thermodynamics. In end-point binding free energy methods the binding free energy is naturally decomposable into physically intuitive contributions such as the solvation free energy and configurational entropy that can provide insights. Here we present a new end-point method called EE-BQH (Effective Energy-Boltzmann-Quasiharmonic) which combines the Boltzmann-Quasiharmonic model for configurational entropy with different solvation free energy methods, such as the continuum solvent PBSA model and the integral equation-based 3D-RISM, to estimate the absolute binding free energy. We compare EE-BQH with other treatments of configurational entropy such as Quasiharmonic models in internal coordinates (QHIC) and in Cartesian coordinates (QHCC), and Normal Mode analysis (NMA), by testing them on the octa acids host-guest complexes from the SAMPL8 blind challenge. The accuracies in the calculated absolute binding free energies strongly depend on the configurational entropy and solvation free energy methods used. QHIC and BQH yield the best agreements with the established potential of mean force (PMF) estimates, with R² of ∼0.7 and mean unsigned error of ∼1.7 kcal mol^-1. These results from the end-point calculations are also in similar agreement with experiments. While 3D-RISM in combination with QHIC or BQH lead to reasonable correlations with the PMF results and experiments, the calculated absolute binding free energies are underestimated by ∼5 kcal mol^-1. While the binding is accompanied by a significant reduction in the ligand translational/rotational entropy, the change in the torsional entropy in these host-guest systems is slightly positive. Compared with BQH, QHIC underestimates the reduction of configurational entropy because of the non-Gaussian probability distributions in the ligand rotation and a small number of torsions. The study highlights the crucial role of configurational entropy in determining binding and demonstrates the potential of using the new end-point method to provide insights in more complex protein-ligand systems.

Indoloquinoline Ligands Favor Intercalation at Quadruplex-Duplex Interfaces

Quadruplex‐duplex (Q‐D) junctions are increasingly considered promising targets for medicinal and technological applications. Here, a Q‐D hybrid with a hairpin‐type snapback loop coaxially stacked onto the quadruplex 3’‐outer tetrad was designed and employed as a target structure for the indoloquinoline ligand SYUIQ‐5. NMR spectral analysis demonstrated high‐affinity binding of the ligand at the quadruplex‐duplex interface with association constants determined by isothermal titration calorimetry of about 10⁷ M⁻¹ and large exothermicities ΔH° of −14 kcal/mol in a 120 mM K⁺ buffer at 40 °C. Determination of the ligand‐bound hybrid structure revealed intercalation of SYUIQ‐5 between 3’‐outer tetrad and the neighboring CG base pair, maximizing π–π stacking as well as electrostatic interactions with guanine carbonyl groups in close vicinity to the positively charged protonated quinoline nitrogen of the tetracyclic indoloquinoline. Exhibiting considerable flexibility, the SYUIQ‐5 sidechain resides in the duplex minor groove. Based on comparative binding studies with the non‐substituted N5‐methylated indoloquinoline cryptolepine, the sidechain is suggested to confer additional affinity and to fix the alignment of the intercalated indoloquinoline aromatic core. However, selectivity for the Q‐D junction mostly relies on the geometry and charge distribution of the indoloquinoline ring system. The presented results are expected to provide valuable guidelines for the design of ligands specifically targeting Q‐D interfaces.

Evaluating Molecular Docking Software for Small Molecule Binding to G-Quadruplex DNA

G-quadruplexes are four-stranded nucleic acid secondary structures of biological significance and have emerged as an attractive drug target. The G4 formed in the MYC promoter (MycG4) is one of the most studied small-molecule targets, and a model system for parallel structures that are prevalent in promoter DNA G4s and RNA G4s. Molecular docking has become an essential tool in structure-based drug discovery for protein targets, and is also increasingly applied to G4 DNA. However, DNA, and in particular G4, binding sites differ significantly from protein targets. Here we perform the first systematic evaluation of four commonly used docking programs (AutoDock Vina, DOCK 6, Glide, and RxDock) for G4 DNA-ligand binding pose prediction using four small molecules whose complex structures with the MycG4 have been experimentally determined in solution. The results indicate that there are considerable differences in the performance of the docking programs and that DOCK 6 with GB/SA rescoring performs better than the other programs. We found that docking accuracy is mainly limited by the scoring functions. The study shows that current docking programs should be used with caution to predict G4 DNA-small molecule binding modes.

Estimating the Roles of Protonation and Electronic Polarization in Absolute Binding Affinity Simulations

Accurate prediction of binding free energies is critical to streamlining the drug development and protein design process. With the advent of GPU acceleration, absolute alchemical methods, which simulate the removal of ligand electrostatics and van der Waals interactions with the protein, have become routinely accessible and provide a physically rigorous approach that enables full consideration of flexibility and solvent interaction. However, standard explicit solvent simulations are unable to model protonation or electronic polarization changes upon ligand transfer from water to the protein interior, leading to inaccurate prediction of binding affinities for charged molecules. Here, we perform extensive simulation totaling ∼540 μs to benchmark the impact of modeling conditions on predictive accuracy for absolute alchemical simulations. Binding to urokinase plasminogen activator (UPA), a protein frequently overexpressed in metastatic tumors, is evaluated for a set of 10 inhibitors with extended flexibility, highly charged character, and titratable properties. We demonstrate that the alchemical simulations can be adapted to utilize the MBAR/PBSA method to improve the accuracy upon incorporating electronic polarization, highlighting the importance of polarization in alchemical simulations of binding affinities. Comparison of binding energy prediction at various protonation states indicates that proper electrostatic setup is also crucial in binding affinity prediction of charged systems, prompting us to propose an alternative binding mode with protonated ligand phenol and Hid-46 at the binding site, a testable hypothesis for future experimental validation.

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.

Molecular Dynamics Study on the Binding of an Anticancer DNA G-Quadruplex Stabilizer, CX-5461, to Human Telomeric, c-KIT1, and c-Myc G-Quadruplexes and a DNA Duplex

DNA G-quadruplex (G4) stabilizer, CX-5461, is in phase I/II clinical trials for advanced cancers with BRCA1/2 deficiencies. A FRET-melting temperature increase assay measured the stabilizing effects of CX-5461 to a DNA duplex (∼10 K), and three G4 forming sequences negatively implicated in the cancers upon its binding: human telomeric (∼30 K), c-KIT1 (∼27 K), and c-Myc (∼25 K). Without experimentally solved structures of these CX-5461-G4 complexes, CX-5461's interactions remain elusive. In this study, we performed a total of 73.5 μs free ligand molecular dynamics binding simulations of CX-5461 to the DNA duplex and three G4s. Three binding modes (top, bottom, and side) were identified for each system and their thermodynamic, kinetic, and structural nature were deciphered. The molecular mechanics/Poisson Boltzmann surface area binding energies of CX-5461 were calculated for the human telomeric (-28.6 kcal/mol), c-KIT1 (-23.9 kcal/mol), c-Myc (-22.0 kcal/mol) G4s, and DNA duplex (-15.0 kcal/mol) systems. These energetic differences coupled with structural differences at the 3' site explained the different melting temperatures between the G4s, while CX-5461's lack of intercalation to the duplex explained the difference between the G4s and duplex. Based on the interaction insight, CX-5461 derivatives were designed and docked, showing higher selectivity to the G4s over the duplex.

Exploring the Free-Energy Landscape and Thermodynamics of Protein-Protein Association

We report the free-energy landscape and thermodynamics of the protein-protein association responsible for the drug-induced multimerization of HIV-1 integrase (IN). Allosteric HIV-1 integrase inhibitors promote aberrant IN multimerization by bridging IN-IN intermolecular interactions. However, the thermodynamic driving forces and kinetics of the multimerization remain largely unknown. Here, we explore the early steps in the IN multimerization by using umbrella sampling and unbiased molecular dynamics simulations in explicit solvent. In direct simulations, the two initially separated dimers spontaneously associate to form near-native complexes that resemble the crystal structure of the aberrant tetramer. Most strikingly, the effective interaction of the protein-protein association is very short-ranged: the two dimers associate rapidly within tens of nanoseconds when their binding surfaces are separated by d ≤ 4.3 Å (less than two water diameters). Beyond this distance, the oligomerization kinetics appears to be diffusion controlled with a much longer association time. The free-energy profile also captured the crucial role of allosteric IN inhibitors in promoting multimerization and explained why several C-terminal domain mutations are remarkably resistant to the drug-induced multimerization. The results also show that at small separation, the protein-protein binding process contains two consecutive phases with distinct thermodynamic signatures. First, interprotein water molecules are expelled to the bulk, resulting in a small increase in entropy, as the solvent entropy gain from the water release is nearly cancelled by the loss of side-chain entropies as the two proteins approach each other. At shorter distances, the two dry binding surfaces adapt to each other to optimize their interaction energy at the expense of further protein configurational entropy loss. Although the binding interfaces feature clusters of hydrophobic residues, overall, the protein-protein association in this system is driven by enthalpy and opposed by entropy.

Subtle structural alterations in G-quadruplex DNA regulate site specificity of fluorescence light-up probes

G-quadruplex (G4) DNA structures are linked to key biological processes and human diseases. Small molecules that target specific G4 DNA structures and signal their presence would therefore be of great value as chemical research tools with potential to further advance towards diagnostic and therapeutic developments. However, the development of these types of specific compounds remain as a great challenge. In here, we have developed a compound with ability to specifically signal a certain c-MYC G4 DNA structure through a fluorescence light-up mechanism. Despite the compound's two binding sites on the G4 DNA structure, only one of them result in the fluorescence light-up effect. This G-tetrad selectivity proved to originate from a difference in flexibility that affected the binding affinity and tilt the compound out of the planar conformation required for the fluorescence light-up mechanism. The intertwined relation between the presented factors is likely the reason for the lack of examples using rational design to develop compounds with turn-on emission that specifically target certain G4 DNA structures. However, this study shows that it is indeed possible to develop such compounds and present insights into the molecular details of specific G4 DNA recognition and signaling to advance future studies of G4 biology.

Combining Alchemical Transformation with a Physical Pathway to Accelerate Absolute Binding Free Energy Calculations of Charged Ligands to Enclosed Binding Sites

We present a new approach to more accurately and efficiently compute the absolute binding free energy for receptor-ligand complexes. Currently, the double decoupling method (DDM) and the potential of mean force method (PMF) are widely used to compute the absolute binding free energy of biomolecular complexes. DDM relies on alchemically decoupling the ligand from its environments, which can be computationally challenging for large ligands and charged ligands because of the large magnitude of the decoupling free energies involved. In contrast, the PMF method uses a physical pathway to directly transfer the ligand from solution to the receptor binding pocket and thus avoids some of the aforementioned problems in DDM. However, the PMF method has its own drawbacks: because of its reliance on a ligand binding/unbinding pathway that is free of steric obstructions from the receptor atoms, the method has difficulty treating ligands with buried atoms. To overcome the limitation in the standard PMF approach and enable buried ligands to be treated, here we develop a new method called AlchemPMF in which steric obstructions along the physical pathway for binding are alchemically removed. We have tested the new approach on two important drug targets involving charged ligands. One is HIV-1 integrase bound to an allosteric inhibitor; the other is the human telomeric DNA G-quadruplex in complex with a natural product protoberberine buried in the binding pocket. For both systems, the new approach leads to more reliable estimates of absolute binding free energies with smaller error bars and closer agreements with experiments compared with those obtained from the existing methods, demonstrating the effectiveness of the new method in overcoming the hysteresis often encountered in PMF binding free energy calculations of such systems. The new approach could also be used to improve the sampling of water equilibration and resolvation of the binding pocket as the ligand is extracted.

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.

To probe the binding pathway of a selective compound (D089-0563) to c-MYC Pu24 G-quadruplex using free ligand binding simulations and Markov state model analysis

D089-0563 is a highly promising anti-cancer compound that selectively binds the transcription-silencing G-quadruplex element (Pu27) at the promoter region of the human c-MYC oncogene; however, its binding mechanism remains elusive.

Why do G-quadruplexes dimerize through the 5'-ends? Driving forces for G4 DNA dimerization examined in atomic detail

G-quadruplexes (G4) are secondary structures formed by guanine-rich nucleic acid sequences and shown to exist in living cells where they participate in regulation of gene expression and chromosome maintenance. G-quadruplexes with solvent-exposed guanine tetrads show the tendency to associate together through cofacial stacking, which may be important for packaging of G4-forming sequences and allows for the design of higher-order G4 DNA structures. To understand the molecular driving forces for G4 association, here, we study the binding interaction between two parallel-stranded G-quadruplexes using all-atom molecular dynamics simulations. The predicted dimerization free energies show that direct binding through the 5’-G-tetrads is the most preferred of all possible end-to-end stacking orientations, consistently with all available experimental data. Decomposition of dimerization enthalpies in combination with simulations at varying ionic strength further indicate that the observed orientational preferences arise from a fine balance between the electrostatic repulsion of the sugar-phosphate backbones and favorable counterion binding at the dimeric interface. We also demonstrate how these molecular-scale findings can be used to devise means of controlling G4 dimerization equilibrium, e.g., by altering salt concentration and using G4-targeted ligands.

Native DNA usually folds to form the canonical double helix, however, under certain conditions, it can also fold into other secondary structures. Some of the most interesting ones are G-quadruplexes (G4)—compact DNA structures in which guanines assemble into multilayered tetrads, and whose formation has been reported at the ends of linear chromosomes (telomeres) and at different regulatory regions of the genome. Although structural and basic energetic properties, as well as some biological functions of G-quadruplexes are quite well understood, not much is known about their propensity to form agregated structures. A very high density of G-quadruplexes at telomeres along with their large exposed planar surfaces indeed favor G4 aggregation through end-to-end stacking, which might be important for the protection of telomeres and DNA packaging. In this research, using computer simulations, we provide insight into molecular origins of stability of the higher-order G-quadruplexes and explain in structural and energetic terms a strong preference for one particular end-to-end stacking orientation. Based on the recognized aggregation driving forces, we also suggest methods for controling the aggregation preferences openining up new opportunities for designing oligomeric G-quadruplexes.

Surface plasmon resonance based analysis of the binding of LYAR protein to the rs368698783 (G>A) polymorphic Aγ-globin gene sequences mutated in β-thalassemia

Recent studies have identified and characterized a novel putative transcriptional repressor site in a 5' untranslated region of the Aγ-globin gene that interacts with the Ly-1 antibody reactive clone (LYAR) protein. LYAR binds the 5'-GGTTAT-3' site of the Aγ-globin gene, and this molecular interaction causes repression of gene transcription. In β-thalassemia patients, a polymorphism has been demonstrated (the rs368698783 G>A polymorphism) within the 5'-GGTTAT-3' LYAR-binding site of the Aγ-globin gene. The major results gathered from surface plasmon resonance based biospecific interaction analysis (SPR-BIA) studies (using crude nuclear extracts, LYAR-enriched lysates, and recombinant LYAR) support the concept that the rs368698783 G>A polymorphism of the Aγ-globin gene attenuates the efficiency of LYAR binding to the LYAR-binding site. This conclusion was fully confirmed by a molecular docking analysis. This might lead to a very important difference in erythroid cells from β-thalassemia patients in respect to basal and induced levels of production of fetal hemoglobin. The novelty of the reported SPR-BIA method is that it allows the characterization and validation of the altered binding of a key nuclear factor (LYAR) to mutated LYAR-binding sites. These results, in addition to theoretical implications, should be considered of interest in applied pharmacology studies as a basis for the screening of drugs able to inhibit LYAR-DNA interactions. This might lead to the identification of molecules facilitating induced increase of γ-globin gene expression and fetal hemoglobin production in erythroid cells, which is associated with possible reduction of the clinical severity of the β-thalassemia phenotype. Graphical abstract.

Ligand Selectivity in the Recognition of Protoberberine Alkaloids by Hybrid-2 Human Telomeric G-Quadruplex: Binding Free Energy Calculation, Fluorescence Binding, and NMR Experiments

The human telomeric G-quadruplex (G4) is an attractive target for developing anticancer drugs. Natural products protoberberine alkaloids are known to bind human telomeric G4 and inhibit telomerase. Among several structurally similar protoberberine alkaloids, epiberberine (EPI) shows the greatest specificity in recognizing the human telomeric G4 over duplex DNA and other G4s. Recently, NMR study revealed that EPI recognizes specifically the hybrid-2 form human telomeric G4 by inducing large rearrangements in the 5′-flanking segment and loop regions to form a highly extensive four-layered binding pocket. Using the NMR structure of the EPI-human telomeric G4 complex, here we perform molecular dynamics free energy calculations to elucidate the ligand selectivity in the recognition of protoberberines by the human telomeric G4. The MM-PB(GB)SA (molecular mechanics-Poisson Boltzmann/Generalized Born) Surface Area) binding free energies calculated using the Amber force fields bsc0 and OL15 correlate well with the NMR titration and binding affinity measurements, with both calculations correctly identifying the EPI as the strongest binder to the hybrid-2 telomeric G4 wtTel26. The results demonstrated that accounting for the conformational flexibility of the DNA-ligand complexes is crucially important for explaining the ligand selectivity of the human telomeric G4. While the MD-simulated (molecular dynamics) structures of the G-quadruplex-alkaloid complexes help rationalize why the EPI-G4 interactions are optimal compared with the other protoberberines, structural deviations from the NMR structure near the binding site are observed in the MD simulations. We have also performed binding free energy calculation using the more rigorous double decoupling method (DDM); however, the results correlate less well with the experimental trend, likely due to the difficulty of adequately sampling the very large conformational reorganization in the G4 induced by the protoberberine binding.

Importance of Chiral Recognition in Designing Metal-Free Ligands for G-Quadruplex DNA

Four pairs of amino acid-functionalized naphthalenediimide enantiomers (d- and l-lysine derived NDIs) were screened toward G-quadruplex forming sequences in telomeres (h-TELO) and oncogene promoters: c-KIT1, c-KIT2, k-RAS and BCL-2. This is the first study to address the effect of point chirality toward G-quadruplex DNA stabilization using purely small organic molecules. Enantioselective behavior toward the majority of ligands was observed, particularly in the case of parallel conformations of c-KIT2 and k-RAS. Additionally, N_ε-Boc-l-Lys-NDI and N_ε-Boc-d-Lys-NDI discriminate between quadruplexes with parallel and hybrid topologies, which has not previously been observed with enantiomeric ligands.

Binding of BRACO19 to a Telomeric G-Quadruplex DNA Probed by All-Atom Molecular Dynamics Simulations with Explicit Solvent

Although BRACO19 is a potent G-quadruplex binder, its potential for clinical usage is hindered by its low selectivity towards DNA G-quadruplex over duplex. High-resolution structures of BRACO19 in complex with neither single-stranded telomeric DNA G-quadruplexes nor B-DNA duplex are available. In this study, the binding pathway of BRACO19 was probed by 27.5 µs molecular dynamics binding simulations with a free ligand (BRACO19) to a DNA duplex and three different topological folds of the human telomeric DNA G-quadruplex (parallel, anti-parallel and hybrid). The most stable binding modes were identified as end stacking and groove binding for the DNA G-quadruplexes and duplex, respectively. Among the three G-quadruplex topologies, the MM-GBSA binding energy analysis suggested that BRACO19′s binding to the parallel scaffold was most energetically favorable. The two lines of conflicting evidence plus our binding energy data suggest conformation-selection mechanism: the relative population shift of three scaffolds upon BRACO19 binding (i.e., an increase of population of parallel scaffold, a decrease of populations of antiparallel and/or hybrid scaffold). This hypothesis appears to be consistent with the fact that BRACO19 was specifically designed based on the structural requirements of the parallel scaffold and has since proven effective against a variety of cancer cell lines as well as toward a number of scaffolds. In addition, this binding mode is only slightly more favorable than BRACO19s binding to the duplex, explaining the low binding selectivity of BRACO19 to G-quadruplexes over duplex DNA. Our detailed analysis suggests that BRACO19′s groove binding mode may not be stable enough to maintain a prolonged binding event and that the groove binding mode may function as an intermediate state preceding a more energetically favorable end stacking pose; base flipping played an important role in enhancing binding interactions, an integral feature of an induced fit binding mechanism.

Using Molecular Dynamics Free Energy Simulation to Compute Binding Affinities of DNA G-Quadruplex Ligands

We provide a practical guide for using molecular dynamics simulation to compute the binding affinity of small molecules in complex with G-quadruplex DNA. Such calculations have a number of applications, such as rescoring docking results and validating docked poses, to inform the discovery of G-quadruplex binders with high affinity and selectivity. This chapter describes two binding free energy protocols: the double decoupling method (DDM) and the potential of mean force method (PMF). We illustrate the application of the two methods using a recent case study of the binding of quindoline with the c-MYC G-quadruplex DNA. For this system, the two methods yield absolute binding free energies within ~2 kcal/mol of the experimental value. We discuss the advantages and disadvantages of these binding free energy methods.

Binding of Telomestatin, TMPyP4, BSU6037, and BRACO19 to a Telomeric G-Quadruplex-Duplex Hybrid Probed by All-Atom Molecular Dynamics Simulations with Explicit Solvent

A promising anticancer therapeutic strategy is the stabilization of telomeric G-quadruplexes using G-quadruplex-binding small molecules. Although many G-quadruplex-specific ligands have been developed, their low potency and selectivity to G-quadruplexes over duplex remains unsolved. Recently, a crystal structure of a telomeric 3' quadruplex-duplex hybrid was reported and the quadruplex-duplex interface was suggested to a good target to address the issues. However, there are no high-resolution complex structures reported for G-quadruplex ligands except for a docked BSU6037. In this study, molecular dynamic (MD) binding simulations with a free ligand were used to study binding poses and dynamics of four representative ligands: telomestatin, TMPyP4, BSU6037, and BRACO19. The MD data showed that BSU6037 was able to fully intercalate into the interface whereas TMPyP4 and BRACO19 could only maintain partial intercalation into the interface and telomestatin only binds at the quadruplex and duplex ends. Both linear ligands, BSU6037 and BRACO19, were able to interact with the interface, yet they were not selective over duplex DNA. The DNA geometry, binding modes, and binding pathways were systematically characterized, and the binding energy was calculated and compared for each system. The interaction of the ligands to the interface was by the means of an induced-fit binding mechanism rather than a lock-key mechanism, consisting of the DNA unfolding at the interface to allow entrance of the drug and then the refolding and repacking of the DNA and the ligand to further stabilize the G-quadruplex. On the basis of the findings in this study, modifications were suggested to optimize the interface binding for TMPyp4 and telomestatin.

Chemical and structural studies provide a mechanistic basis for recognition of the MYC G-quadruplex

G-quadruplexes (G4s) are noncanonical DNA structures that frequently occur in the promoter regions of oncogenes, such as MYC, and regulate gene expression. Although G4s are attractive therapeutic targets, ligands capable of discriminating between different G4 structures are rare. Here, we describe DC-34, a small molecule that potently downregulates MYC transcription in cancer cells by a G4-dependent mechanism. Inhibition by DC-34 is significantly greater for MYC than other G4-driven genes. We use chemical, biophysical, biological, and structural studies to demonstrate a molecular rationale for the recognition of the MYC G4. We solve the structure of the MYC G4 in complex with DC-34 by NMR spectroscopy and illustrate specific contacts responsible for affinity and selectivity. Modification of DC-34 reveals features required for G4 affinity, biological activity, and validates the derived NMR structure. This work advances the design of quadruplex-interacting small molecules to control gene expression in therapeutic areas such as cancer.

Targeting noncoding nucleic acids with small molecules represents an important and significant challenge in chemical biology and drug discovery. Here the authors characterize DC-34, a small molecule that exhibits selective binding to specific G4 structures, and provide a structural basis for its selectivity

Fluorescence-based tools to probe G-quadruplexes in cell-free and cellular environments

Biophysical and biochemical investigations provide compelling evidence connecting the four-stranded G-quadruplex (GQ) structure with its role in regulating multiple cellular processes. Hence, modulating the function of GQs by using small molecule binders is being actively pursued as a strategy to develop new chemotherapeutic agents. However, sequence diversity and structural polymorphism of GQs have posed immense challenges in terms of understanding what conformation a G-rich sequence adopts inside the cell and how to specifically target a GQ motif amidst several other GQ-forming sequences. In this context, here we review recent developments in the applications of biophysical tools that use fluorescence readout to probe the GQ structure and recognition in cell-free and cellular environments. First, we provide a detailed discussion on the utility of covalently labeled environment-sensitive fluorescent nucleoside analogs in assessing the subtle difference in GQ structures and their ligand binding abilities. Furthermore, a detailed discussion on structure-specific antibodies and small molecule probes used to visualize and confirm the existence of DNA and RNA GQs in cells is provided. We also highlight the open challenges in the study of tetraplexes (GQ and i-motif structures) and how addressing these challenges by developing new tools and techniques will have a profound impact on tetraplex-directed therapeutic strategies.

Assessment of a Single Decoupling Alchemical Approach for the Calculation of the Absolute Binding Free Energies of Protein-Peptide Complexes

The computational modeling of peptide inhibitors to target protein-protein binding interfaces is growing in interest as these are often too large, too shallow, and too feature-less for conventional small molecule compounds. Here, we present a rare successful application of an alchemical binding free energy method for the calculation of converged absolute binding free energies of a series of protein-peptide complexes. Specifically, we report the binding free energies of a series of cyclic peptides derived from the LEDGF/p75 protein to the integrase receptor of the HIV1 virus. The simulations recapitulate the effect of mutations relative to the wild-type binding motif of LEDGF/p75, providing structural, energetic and dynamical interpretations of the observed trends. The equilibration and convergence of the calculations are carefully analyzed. Convergence is aided by the adoption of a single-decoupling alchemical approach with implicit solvation, which circumvents the convergence difficulties of conventional double-decoupling protocols. We hereby present the single-decoupling methodology and critically evaluate its advantages and limitations. We also discuss some of the challenges and potential pitfalls of binding free energy calculations for complex molecular systems which have generally limited their applicability to the quantitative study of protein-peptide binding equilibria.