Design, Synthesis, and Investigation of the Pharmacokinetics and Anticancer Activities of Indenoisoquinoline Derivatives That Stabilize the G-Quadruplex in the MYC Promoter and Inhibit Topoisomerase I.

G-quadruplexes are noncanonical four-stranded DNA secondary structures. MYC is a master oncogene and the G-quadruplex formed in the MYC promoter functions as a transcriptional silencer and can be stabilized by small molecules. We have previously revealed a novel mechanism of action for indenoisoquinoline anticancer drugs, dual-downregulation of MYC and inhibition of topoisomerase I. Herein, we report the design and synthesis of novel 7-aza-8,9-methylenedioxyindenoisoquinolines based on desirable substituents and π-π stacking interactions. These compounds stabilize the MYC promoter G-quadruplex, significantly lower MYC levels in cancer cells, and inhibit topoisomerase I. MYC targeting was demonstrated by differential activities in Raji vs CA-46 cells and cytotoxicity in MYC-dependent cell lines. Cytotoxicities in the NCI-60 panel of human cancer cell lines were investigated. Favorable pharmacokinetics were established, and in vivo anticancer activities were demonstrated in xenograft mouse models. Furthermore, favorable brain penetration, brain pharmacokinetics, and anticancer activity in an orthotopic glioblastoma mouse model were demonstrated.

Negative Electron Transfer Collision-Induced Dissociation of G-Quadruplexes: Uncovering the Guanine Radical Anion Loss Pathway.

G-quadruplex (G4) DNA can form highly stable secondary structures in the presence of metal cations, and research has shown its potential as a transcriptional regulator for oncogenes in the human genome. In order to explore the interactions of DNA with metal cations using mass spectrometry, employing complementary fragmentation methods can enhance structural information. This study explores the use of ion-ion reactions for sequential negative electron transfer collision-induced dissociation (nET-CID) as a complement to traditional ion-trap CID (IT-CID). The resulting nET-CID data for G4 anions with and without metal cations show an increase in fragment ion type diversity and yield of structurally informative ions relative to IT-CID. The nET-CID yields greater sequence coverage by virtue of fragmentation at the 3'-side of thymine residues, which is lacking with IT-CID. Potassium adductions to backbone fragments in IT-CID and nET-CID spectra were nearly identical. Of note is a prominent fragment resulting from a loss of a 149 Da anion seen in nET-CID of large, G-rich sequences, proposed to be radical anion guanine loss. Neutral loss of neutral guanine (151 Da) and deprotonated nucleobase loss (150 Da) have been previously reported, but this is the first report of radical anion guanine loss (149 Da). Confirmation of the identity of the 149 Da anion results from the examination of the homonucleobase sequence 5'-GGGGGGGG-3'. Loss of a charged adenine radical anion at much lower relative abundance was also noted for the sequence 5'-AAAAAAAA-3'. DFT modeling indicates that the loss of a nucleobase as a radical anion from odd-electron nucleic acid anions is a thermodynamically favorable fragmentation pathway for G.

Best method to determine DNA G-quadruplex folding: The 1H-13C HSQC NMR experiment.

NMR spectroscopy is the major method for G-quadruplex structure determination under physiologically relevant solution conditions. Unlike duplex B-DNA, in which all nucleotides adopt an anti glycosidic conformation, the core tetrad-guanines in a G-quadruplex can adopt anti or syn glycosidic conformation depending on the folding structure. An experimental method that can clearly and unambiguously determine syn and anti tetrad-Gs in a G-quadruplex is highly desirable and necessary. In the present study, we exploit the advantages of the 1H-13C HSQC experiment to determine tetrad-G's glycosidic conformation and thus folding topology of G-quadruplexes. We use several examples to demonstrate the clear and straightforward determination of the guanine glycosidic conformations and G-quadruplex folding structures. Moreover, 1H-13C HSQC data can readily identify adenine H2 resonances as well as determine unusual syn conformation in loop and flanking sequences, a challenging task by standard 2D NOESY.

Gas-Phase Fragmentation as a Probe of G-Quadruplex Formation.

G-quadruplex (G4) DNA is found in oncogene promoters and human telomeres and is an attractive anticancer target. Stable G4 structures form in guanine-rich sequences in the presence of metal cations and can stabilize further with specific ligand adduction. To explore the preservation and stability of this secondary structure with mass spectrometry, gas-phase collision-induced dissociation kinetics of G4-like and non-G4-like ion structures were determined in a linear quadrupole ion trap. This study focused on a sequence from the promoter of the MYC oncogene, MycG4, and a mutant non-G4-forming sequence, MycNonG4. At relatively high ion activation energies, the backbone fragmentation patterns of the MycG4 and MycNonG4 are similar, while potassium ion-stabilized G4-folded [MycG4 + 2K-7H]5- and counterpart [MycG4-5H]5- ions are essentially indistinguishable, indicating that high-energy fragmentation is not sensitive to the G4 structure. At low energies, the backbone fragmentation patterns of MycG4 and MycNonG4 are significantly different. For MycG4, fragmentation over time differed significantly between the potassium-bound and free structures, reflecting the preservation of the G4 structure in the gas phase. Kinetic measurements revealed the [MycG4 + 2K-7H]5- ions to fragment two to three times more slowly than the [MycG4-5H]5-. Results for the control MycNonG4 indicated that the phenomena noted for [MycG4 + 2K-7H]5- ions are specific to G4-folding. Therefore, our data show that gentle activation conditions can lead to fragmentation behavior that is sensitive to G-quadruplex structure, revealing differences in kinetic stabilities of isomeric structures as well as the regions of the sequence that are directly involved in forming these structures.

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.

Decoding regulatory associations of G-quadruplex with epigenetic and transcriptomic functional components.

G-quadruplex (G4) has been previously observed to be associated with gene expression. In this study, we performed integrative analysis on G4 multi-omics data from in-silicon prediction and ChIP-seq in human genome. Potential G4 sites were classified into three distinguished groups, such as one group of high-confidence G4-forming locations (G4-II) and groups only containing either ChIP-seq detected G4s (G4-I) or predicted G4 motif candidates (G4-III). We explored the associations of different-confidence G4 groups with other epigenetic regulatory elements, including CpG islands, chromatin status, enhancers, super-enhancers, G4 locations compared to the genes, and DNA methylation. Our elastic net regression model revealed that G4 structures could correlate with gene expression in two opposite ways depending on their locations to the genes as well as G4-forming DNA strand. Some transcription factors were identified to be over-represented with G4 emergence. The motif analysis discovered distinct consensus sequences enriched in the G4 feet, the flanking regions of two groups of G4s. We found high GC content in the feet of high-confidence G4s (G4-II) when compared to high TA content in solely predicted G4 feet of G4-III. Overall, we uncovered the comprehensive associations of G4 formations or predictions with other epigenetic and transcriptional elements which potentially coordinate gene transcription.

Subnanomolar Sensitive Stable Isotopic Determination in CO2 by Tunable Infrared Laser Absorption Spectroscopy.

Oxygen and carbon stable isotope ratios (18O/16O, 13C/12C, and 17O/16O) of CO2 have been crucial in helping us understand Earth and planetary systems. These ratios have also been used in medicine for the noninvasive diagnosis of diseases from exhaled breath and for quantifying biochemical or metabolic reactions and in determining the production area of agricultural products. The current method for measuring the stable isotope ratios of CO2 is primarily gas-source isotope ratio mass spectroscopy (IRMS). Due to the recent demand for isotopic microanalysis of carbonates and organic compounds, the sample size required for isotopic measurements has been reduced to approximately 2 nmol CO2 (equivalent to 0.2 µg CaCO3 and 24 ng carbon) by using high-precision IRMS. We report a novel method using tunable mid-infrared laser direct absorption spectroscopy (TILDAS) for sensitive measurements of 18O/16O and 13C/12C in subnanomolar CO2. This method can accurately measure 18O/16O and 13C/12C in CO2 with a repeatability of less than 0.03‰ (n = 28) in a range of 0.3 nmol (equivalent to 0.03 µg CaCO3 and 3.8 ng carbon) to 30 nmol. This is a sample size 1 order of magnitude smaller than currently available sensitive analytical techniques. In addition, the TILDAS system measures 17O/16O simultaneously with a repeatability of less than 0.06‰ (n = 28). Our method is a major advance in supersensitive CO2 stable isotopic analyses for various fields.

Oxidative Damage Induces a Vacancy G-Quadruplex That Binds Guanine Metabolites: Solution Structure of a cGMP Fill-in Vacancy G-Quadruplex in the Oxidized BLM Gene Promoter.

Guanine (G)-oxidation to 8-oxo-7,8-dihydroguanine (OG) by reactive oxygen species in genomic DNA has been implicated with various human diseases. G-quadruplex (G4)-forming sequences in gene promoters are highly susceptible to G-oxidation, which can subsequently cause gene activation. However, the underlying G4 structural changes that result from OG modifications remain poorly understood. Herein, we investigate the effect of G-oxidation on the BLM gene promoter G4. For the first time, we show that OG can induce a G-vacancy-containing G4 (vG4), which can be filled in and stabilized by guanine metabolites and derivatives. We determined the NMR solution structure of the cGMP-fill-in oxidized BLM promoter vG4. This is the first complex structure of an OG-induced vG4 from a human gene promoter sequence with a filled-in guanine metabolite. The high-resolution structure elucidates the structural features of the specific 5'-end cGMP-fill-in for the OG-induced vG4. Interestingly, the OG is removed from the G-core and becomes part of the 3'-end capping structure. A series of guanine metabolites and derivatives are evaluated for fill-in activity to the oxidation-induced vG4. Significantly, cellular guanine metabolites, such as cGMP and GTP, can bind and stabilize the OG-induced vG4, suggesting their potential regulatory role in response to oxidative damage in physiological and pathological processes. Our work thus provides exciting insights into how oxidative damage and cellular metabolites may work together through a G4-based epigenetic feature for gene regulation. Furthermore, the NMR structure can guide the rational design of small-molecule inhibitors that specifically target the oxidation-induced vG4s.

RNA helicase DDX5 enables STAT1 mRNA translation and interferon signalling in hepatitis B virus replicating hepatocytes.

OBJECTIVE: RNA helicase DDX5 is downregulated during HBV replication and poor prognosis HBV-related hepatocellular carcinoma (HCC). The objective of this study is to investigate the role of DDX5 in interferon (IFN) signalling. We provide evidence of a novel mechanism involving DDX5 that enables translation of transcription factor STAT1 mediating the IFN response. DESIGN AND RESULTS: Molecular, pharmacological and biophysical assays were used together with cellular models of HBV replication, HCC cell lines and liver tumours. We demonstrate that DDX5 regulates STAT1 mRNA translation by resolving a G-quadruplex (rG4) RNA structure, proximal to the 5' end of STAT1 5'UTR. We employed luciferase reporter assays comparing wild type (WT) versus mutant rG4 sequence, rG4-stabilising compounds, CRISPR/Cas9 editing of the STAT1-rG4 sequence and circular dichroism determination of the rG4 structure. STAT1-rG4 edited cell lines were resistant to the effect of rG4-stabilising compounds in response to IFN-α, while HCC cell lines expressing low DDX5 exhibited reduced IFN response. Ribonucleoprotein and electrophoretic mobility assays demonstrated direct and selective binding of RNA helicase-active DDX5 to the WT STAT1-rG4 sequence. Immunohistochemistry of normal liver and liver tumours demonstrated that absence of DDX5 corresponded to absence of STAT1. Significantly, knockdown of DDX5 in HBV infected HepaRG cells reduced the anti-viral effect of IFN-α. CONCLUSION: RNA helicase DDX5 resolves a G-quadruplex structure in 5'UTR of STAT1 mRNA, enabling STAT1 translation. We propose that DDX5 is a key regulator of the dynamic range of IFN response during innate immunity and adjuvant IFN-α therapy.

Berberine Molecular Recognition of the Parallel MYC G-Quadruplex in Solution.

The medicinal natural product berberine is one of the most actively studied and pursued G-quadruplex (G4)-ligands. The major G-quadruplex formed in the promoter region of the MYC oncogene (MycG4) is an attractive drug target and a prominent example and model structure for parallel G-quadruplexes. G4-targeted berberine derivatives have been actively developed; however, the analogue design was based on a previous crystal structure in which berberine binds as a dimer to a parallel G-quadruplex. Herein, we show that in solution, the binding mode and stoichiometry of berberine are substantially different from the crystal structure: berberine binds as a monomer to MycG4 using a base-recruitment mechanism with a reversed orientation in that the positively charged convex side is actually positioned above the tetrad center. Our structure provides a physiologically relevant basis for the future structure-based rational design of G4-targeted berberine derivatives, and this study demonstrates that it is crucial to validate the ligand-DNA interactions.

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.

Solution Structure of Ternary Complex of Berberine Bound to a dGMP-Fill-In Vacancy G-Quadruplex Formed in the PDGFR-ß Promoter.

The G-quadruplexes (G4s) formed in the PDGFR-ß gene promoter are transcriptional modulators and amenable to small-molecule targeting. Berberine (BER), a clinically important natural isoquinoline alkaloid, has gained increasing attention due to its potential as anticancer drug. We previously showed that the PDGFR-ß gene promoter forms a unique vacancy G4 (vG4) that can be filled in and stabilized by guanine metabolites, such as dGMP. Herein, we report the high-resolution NMR structure of a ternary complex of berberine bound to the dGMP-fill-in PDGFR-ß vG4 in potassium solution. This is the first small-molecule complex structure of a fill-in vG4. This ternary complex has a 2:1:1 binding stoichiometry with a berberine molecule bound at each the 5'- and 3'-end of the 5'-dGMP-fill-in PDGFR-ß vG4. Each berberine recruits the adjacent adenine residue from the 5'- or 3'-flanking sequence to form a "quasi-triad plane" that covers the external G-tetrad of the fill-in vG4, respectively. Significantly, berberine covers and stabilizes the fill-in dGMP. The binding of berberine involves both π-stacking and electrostatic interactions, and the fill-in dGMP is covered and well-protected by berberine. The NMR structure can guide rational design of berberine analogues that target the PDGFR-ß vG4 or dGMP-fill-in vG4. Moreover, our structure provides a molecular basis for designing small-molecule guanine conjugates to target vG4s.

Novel DNA Bis-Intercalator XR5944 as a Potent Anticancer Drug-Design and Mechanism of Action.

This review is dedicated to Professor William A. Denny's discovery of XR5944 (also known as MLN944). XR5944 is a DNA-targeted agent with exceptionally potent antitumor activity and a novel DNA binding mode, bis-intercalation and major groove binding, as well as a novel mechanism of action, transcription inhibition. This novel anticancer compound represents a remarkable accomplishment resulting from two decades of drug discovery by Professor Denny and coworkers. Here, we review our work on the structural study of the DNA binding mode of XR5944 and mechanistic study of XR5944 action.

Structural recognition of the MYC promoter G-quadruplex by a quinoline derivative: insights into molecular targeting of parallel G-quadruplexes.

DNA G-Quadruplexes (G4s) formed in oncogene promoters regulate transcription. The oncogene MYC promoter G4 (MycG4) is the most prevalent G4 in human cancers. However, the most studied MycG4 sequence bears a mutated 3'-residue crucial for ligand recognition. Here, we report a new drug-like small molecule PEQ without a large aromatic moiety that specifically binds MycG4. We determined the NMR solution structures of the wild-type MycG4 and its 2:1 PEQ complex, as well as the structure of the 2:1 PEQ complex of the widely used mutant MycG4. Comparison of the two complex structures demonstrates specific molecular recognition of MycG4 and shows the clear effect of the critical 3'-mutation on the drug binding interface. We performed a systematic analysis of the four available complex structures involving the same mutant MycG4, which can be considered a model system for parallel G4s, and revealed for the first time that the flexible flanking residues are recruited in a conserved and sequence-specific way, as well as unused potential for selective ligand-G4 hydrogen-bond interactions. Our results provide the true molecular basis for MycG4-targeting drugs and new critical insights into future rational design of drugs targeting MycG4 and parallel G4s that are prevalent in promoter and RNA G4s.

Custom G4 Microarrays Reveal Selective G-Quadruplex Recognition of Small Molecule BMVC: A Large-Scale Assessment of Ligand Binding Selectivity.

G-quadruplexes (G4) are considered new drug targets for human diseases such as cancer. More than 10,000 G4s have been discovered in human chromatin, posing challenges for assessing the selectivity of a G4-interactive ligand. 3,6-bis(1-Methyl-4-vinylpyridinium) carbazole diiodide (BMVC) is the first fluorescent small molecule for G4 detection in vivo. Our previous structural study shows that BMVC binds to the MYC promoter G4 (MycG4) with high specificity. Here, we utilize high-throughput, large-scale custom DNA G4 microarrays to analyze the G4-binding selectivity of BMVC. BMVC preferentially binds to the parallel MycG4 and selectively recognizes flanking sequences of parallel G4s, especially the 3'-flanking thymine. Importantly, the microarray results are confirmed by orthogonal NMR and fluorescence binding analyses. Our study demonstrates the potential of custom G4 microarrays as a platform to broadly and unbiasedly assess the binding selectivity of G4-interactive ligands, and to help understand the properties that govern molecular recognition.

Custom DNA Microarrays Reveal Diverse Binding Preferences of Proteins and Small Molecules to Thousands of G-Quadruplexes.

Single-stranded DNA (ssDNA) containing four guanine repeats can form G-quadruplex (G4) structures. While cellular proteins and small molecules can bind G4s, it has been difficult to broadly assess their DNA-binding specificity. Here, we use custom DNA microarrays to examine the binding specificities of proteins, small molecules, and antibodies across ∼15,000 potential G4 structures. Molecules used include fluorescently labeled pyridostatin (Cy5-PDS, a small molecule), BG4 (Cy5-BG4, a G4-specific antibody), and eight proteins (GST-tagged nucleolin, IGF2, CNBP, FANCJ, PIF1, BLM, DHX36, and WRN). Cy5-PDS and Cy5-BG4 selectively bind sequences known to form G4s, confirming their formation on the microarrays. Cy5-PDS binding decreased when G4 formation was inhibited using lithium or when ssDNA features on the microarray were made double-stranded. Similar conditions inhibited the binding of all other molecules except for CNBP and PIF1. We report that proteins have different G4-binding preferences suggesting unique cellular functions. Finally, competition experiments are used to assess the binding specificity of an unlabeled small molecule, revealing the structural features in the G4 required to achieve selectivity. These data demonstrate that the microarray platform can be used to assess the binding preferences of molecules to G4s on a broad scale, helping to understand the properties that govern molecular recognition.

PDGFR-ß Promoter Forms a Vacancy G-Quadruplex that Can Be Filled in by dGMP: Solution Structure and Molecular Recognition of Guanine Metabolites and Drugs.

Aberrant expression of PDGFR-ß is associated with a number of diseases. The G-quadruplexes (G4s) formed in PDGFR-ß gene promoter are transcriptional modulators and amenable to small molecule targeting. The major G4 formed in the PDGFR-ß gene promoter was previously shown to have a broken G-strand. Herein, we report that the PDGFR-ß gene promoter sequence forms a vacancy G-quadruplex (vG4) which can be filled in and stabilized by physiologically relevant guanine metabolites, such as dGMP, GMP, and cGMP, as well as guanine-derivative drugs. We determined the NMR structure of the dGMP-fill-in PDGFR-ß vG4 in K+ solution. This is the first structure of a guanine-metabolite-fill-in vG4 based on a human gene promoter sequence. Our structure and systematic analysis elucidate the contributions of Hoogsten hydrogen bonds, sugar, and phosphate moieties to the specific G-vacancy fill-in. Intriguingly, an equilibrium of 3'- and 5'-end vG4s is present in the PDGFR-ß promoter sequence, and dGMP favors the 5'-end fill-in. Guanine metabolites and drugs were tested and showed a conserved selectivity for the 5'-vacancy, except for cGMP. cGMP binds both the 3'- and 5'-end vG4s and forms two fill-in G4s with similar population. Significantly, guanine metabolites are involved in many physiological and pathological processes in human cells; thus, our results provide a structural basis to understand their potential regulatory functions by interaction with promoter vG4s. Moreover, the NMR structure can guide rational design of ligands that target the PDGFR-ß vG4.

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.

Structures of 1:1 and 2:1 complexes of BMVC and MYC promoter G-quadruplex reveal a mechanism of ligand conformation adjustment for G4-recognition.

BMVC is the first fluorescent probe designed to detect G-quadruplexes (G4s) in vivo. The MYC oncogene promoter forms a G4 (MycG4) which acts as a transcription silencer. Here, we report the high-affinity and specific binding of BMVC to MycG4 with unusual slow-exchange rates on the NMR timescale. We also show that BMVC represses MYC in cancer cells. We determined the solution structures of the 1:1 and 2:1 BMVC-MycG4 complexes. BMVC first binds the 5'-end of MycG4 to form a 1:1 complex with a well-defined structure. At higher ratio, BMVC also binds the 3'-end to form a second complex. In both complexes, the crescent-shaped BMVC recruits a flanking DNA residue to form a BMVC-base plane stacking over the external G-tetrad. Remarkably, BMVC adjusts its conformation to a contracted form to match the G-tetrad for an optimal stacking interaction. This is the first structural example showing the importance of ligand conformational adjustment in G4 recognition. BMVC binds the more accessible 5'-end with higher affinity, whereas sequence specificity is present at the weaker-binding 3'-site. Our structures provide insights into specific recognition of MycG4 by BMVC and useful information for design of G4-targeted anticancer drugs and fluorescent probes.