Structure of a 28.5 kDa duplex-embedded G-quadruplex system resolved to 7.4 Å resolution with cryo-EM.

Genomic regions with high guanine content can fold into non-B form DNA four-stranded structures known as G-quadruplexes (G4s). Extensive in vivo investigations have revealed that promoter G4s are transcriptional regulators. Little structural information exists for these G4s embedded within duplexes, their presumed genomic environment. Here, we report the 7.4 Å resolution structure and dynamics of a 28.5 kDa duplex-G4-duplex (DGD) model system using cryo-EM, molecular dynamics, and small-angle X-ray scattering (SAXS) studies. The DGD cryo-EM refined model features a 53° bend induced by a stacked duplex-G4 interaction at the 5' G-tetrad interface with a persistently unstacked 3' duplex. The surrogate complement poly dT loop preferably stacks onto the 3' G-tetrad interface resulting in occlusion of both 5' and 3' tetrad interfaces. Structural analysis shows that the DGD model is quantifiably more druggable than the monomeric G4 structure alone and represents a new structural drug target. Our results illustrate how the integration of cryo-EM, MD, and SAXS can reveal complementary detailed static and dynamic structural information on DNA G4 systems.

G-quadruplex DNA: A Longer Story.

G-quadruplexes (G4s) are distinctive four-stranded DNA or RNA structures found within cells that are thought to play functional roles in gene regulation and transcription, translation, recombination, and DNA damage/repair. While G4 structures can be uni-, bi-, or tetramolecular with respect to strands, folded unimolecular conformations are most significant in vivo. Unimolecular G4 can potentially form in sequences with runs of guanines interspersed with what will become loops in the folded structure: 5'GxLyGxLyGxLyGx, where x is typically 2-4 and y is highly variable. Such sequences are highly conserved and specifically located in genomes. In the folded structure, guanines from each run combine to form planar tetrads with four hydrogen-bonded guanine bases; these tetrads stack on one another to produce four strand segments aligned in specific parallel or antiparallel orientations, connected by the loop sequences. Three types of loops (lateral, diagonal, or "propeller") have been identified. The stacked tetrads form a central cavity that features strong coordination sites for monovalent cations that stabilize the G4 structure, with potassium or sodium preferred. A single monomeric G4 typically forms from a sequence containing roughly 20-30 nucleotides. Such short sequences have been the primary focus of X-ray crystallographic or NMR studies that have produced high-resolution structures of a variety of monomeric G4 conformations. These structures are often used as the basis for drug design efforts to modulate G4 function.We believe that the focus on monomeric G4 structures formed by such short sequences is perhaps myopic. Such short sequences for structural studies are often arbitrarily selected and removed from their native genomic sequence context, and then are often changed from their native sequences by base substitutions or deletions intended to optimize the formation of a homogeneous G4 conformation. We believe instead that G-quadruplexes prefer company and that in a longer natural sequence context multiple adjacent G4 units can form to combine into more complex multimeric G4 structures with richer topographies than simple monomeric forms. Bioinformatic searches of the human genome show that longer sequences with the potential for forming multiple G4 units are common. Telomeric DNA, for example, has a single-stranded overhang of hundreds of nucleotides with the requisite repetitive sequence with the potential for formation of multiple G4s. Numerous extended promoter sequences have similar potentials for multimeric G4 formation. X-ray crystallography and NMR methods are challenged by these longer sequences (>30 nt), so other tools are needed to explore the possible multimeric G4 landscape. We have implemented an integrated structural biology approach to address this challenge. This approach integrates experimental biophysical results with atomic-level molecular modeling and molecular dynamics simulations that provide quantitatively testable model structures. In every long sequence we have studied so far, we found that multimeric G4 structures readily form, with a surprising diversity of structures dependent on the exact native sequence used. In some cases, stable hairpin duplexes form along with G4 units to provide an even richer landscape. This Account provides an overview of our approach and recent progress and provides a new perspective on the G-quadruplex folding landscape.

Drug discovery of small molecules targeting the higher-order hTERT promoter G-quadruplex.

DNA G-quadruplexes (G4s) are now widely accepted as viable targets in the pursuit of anticancer therapeutics. To date, few small molecules have been identified that exhibit selectivity for G4s over alternative forms of DNA, such as the ubiquitous duplex. We posit that the lack of current ligand specificity arises for multiple reasons: G4 atomic models are often small, monomeric, single quadruplex structures with few or no druggable pockets; targeting G-tetrad faces frequently results in the enrichment of extended electron-deficient polyaromatic end-pasting scaffolds; and virtual drug discovery efforts often under-sample chemical search space. We show that by addressing these issues we can enrich for non-standard molecular templates that exhibit high selectivity towards G4s over other forms of DNA. We performed an extensive virtual screen against the higher-order hTERT core promoter G4 that we have previously characterized, targeting 12 of its unique loop and groove pockets using libraries containing 40 million drug-like compounds for each screen. Using our drug discovery funnel approach, which utilizes high-throughput fluorescence thermal shift assay (FTSA) screens, microscale thermophoresis (MST), and orthogonal biophysical methods, we have identified multiple unique G4 binding scaffolds. We subsequently used two rounds of catalogue-based SAR to increase the affinity of a disubstituted 2-aminoethyl-quinazoline that stabilizes the higher-order hTERT G-quadruplex by binding across its G4 junctional sites. We show selectivity of its binding affinity towards hTERT is virtually unaffected in the presence of near-physiological levels of duplex DNA, and that this molecule downregulates hTERT transcription in breast cancer cells.

Long promoter sequences form higher-order G-quadruplexes: an integrative structural biology study of c-Myc, k-Ras and c-Kit promoter sequences.

We report on higher-order G-quadruplex structures adopted by long promoter sequences obtained by an iterative integrated structural biology approach. Our approach uses quantitative biophysical tools (analytical ultracentrifugation, small-angle X-ray scattering, and circular dichroism spectroscopy) combined with modeling and molecular dynamics simulations, to derive self-consistent structural models. The formal resolution of our approach is 18 angstroms, but in some cases structural features of only a few nucleotides can be discerned. We report here five structures of long (34-70 nt) wild-type sequences selected from three cancer-related promoters: c-Myc, c-Kit and k-Ras. Each sequence studied has a unique structure. Three sequences form structures with two contiguous, stacked, G-quadruplex units. One longer sequence from c-Myc forms a structure with three contiguous stacked quadruplexes. A longer c-Kit sequence forms a quadruplex-hairpin structure. Each structure exhibits interfacial regions between stacked quadruplexes or novel loop geometries that are possible druggable targets. We also report methodological advances in our integrated structural biology approach, which now includes quantitative CD for counting stacked G-tetrads, DNaseI cleavage for hairpin detection and SAXS model refinement. Our results suggest that higher-order quadruplex assemblies may be a common feature within the genome, rather than simple single quadruplex structures.

Identification and characterization of potent, selective, and efficacious inhibitors of human arylamine N-acetyltransferase 1.

Arylamine N-acetyltransferase 1 (NAT1) plays a pivotal role in the metabolism of carcinogens and is a drug target for cancer prevention and/or treatment. A protein-ligand virtual screening of 2 million chemicals was ranked for predicted binding affinity towards the inhibition of human NAT1. Sixty of the five hundred top-ranked compounds were tested experimentally for inhibition of recombinant human NAT1 and N-acetyltransferase 2 (NAT2). The most promising compound 9,10-dihydro-9,10-dioxo-1,2-anthracenediyl diethyl ester (compound 10) was found to be a potent and selective NAT1 inhibitor with an in vitro IC50 of 0.75 µM. Two structural analogs of this compound were selective but less potent for inhibition of NAT1 whereas a third structural analog 1,2-dihydroxyanthraquinone (a compound 10 hydrolysis product also known as Alizarin) showed comparable potency and efficacy for human NAT1 inhibition. Compound 10 inhibited N-acetylation of the arylamine carcinogen 4-aminobiphenyl (ABP) both in vitro and in DNA repair-deficient Chinese hamster ovary (CHO) cells in situ stably expressing human NAT1 and CYP1A1. Compound 10 and Alizarin effectively inhibited NAT1 in cryopreserved human hepatocytes whereas inhibition of NAT2 was not observed. Compound 10 caused concentration-dependent reductions in DNA adduct formation and DNA double-strand breaks following metabolism of aromatic amine carcinogens beta-naphthylamine and/or ABP in CHO cells. Compound 10 inhibited proliferation and invasion in human breast cancer cells and showed selectivity towards tumorigenic versus non-tumorigenic cells. In conclusion, our study identifies potent, selective, and efficacious inhibitors of human NAT1. Alizarin's ability to inhibit NAT1 could reduce breast cancer metastasis particularly to bone.

Nuclear Pyruvate Kinase M2 (PKM2) Contributes to Phosphoserine Aminotransferase 1 (PSAT1)-Mediated Cell Migration in EGFR-Activated Lung Cancer Cells.

An elevated expression of phosphoserine aminotransferase 1 (PSAT1) has been observed in multiple tumor types and is associated with poorer clinical outcomes. Although PSAT1 is postulated to promote tumor growth through its enzymatic function within the serine synthesis pathway (SSP), its role in cancer progression has not been fully characterized. Here, we explore a putative non-canonical function of PSAT1 that contributes to lung tumor progression. Biochemical studies found that PSAT1 selectively interacts with pyruvate kinase M2 (PKM2). Amino acid mutations within a PKM2-unique region significantly reduced this interaction. While PSAT1 loss had no effect on cellular pyruvate kinase activity and PKM2 expression in non-small-cell lung cancer (NSCLC) cells, fractionation studies demonstrated that the silencing of PSAT1 in epidermal growth factor receptor (EGFR)-mutant PC9 or EGF-stimulated A549 cells decreased PKM2 nuclear translocation. Further, PSAT1 suppression abrogated cell migration in these two cell types whereas PSAT1 restoration or overexpression induced cell migration along with an elevated nuclear PKM2 expression. Lastly, the nuclear re-expression of the acetyl-mimetic mutant of PKM2 (K433Q), but not the wild-type, partially restored cell migration in PSAT1-silenced cells. Therefore, we conclude that, in response to EGFR activation, PSAT1 contributes to lung cancer cell migration, in part, by promoting nuclear PKM2 translocation.

Novel small molecule inhibitor of Kpnß1 induces cell cycle arrest and apoptosis in cancer cells.

Karyopherin beta 1 (Kpnß1) is a major nuclear import receptor that mediates the import of cellular cargoes into the nucleus. Recently it has been shown that Kpnß1 is highly expressed in several cancers, and its inhibition by siRNA induces apoptotic cancer cell death, while having little effect on non-cancer cells. This study investigated the effect of a novel small molecule, Inhibitor of Nuclear Import-60 (INI-60), on cancer cell biology, as well as nuclear import activities associated with Kpnß1, and cancer progression in vivo using cervical and oesophageal cancer cell lines. INI-60 treatment resulted in the inhibition of cancer cell proliferation, colony formation, migration and invasion, and induced a G1/S cell cycle arrest, followed by cancer cell death via apoptosis. Non-cancer cells were minimally affected by INI-60 at concentrations that inhibited cancer cells. INI-60 treatment altered the localisation of Kpnß1 and its cargoes, NFκB/p65, NFAT and AP-1, and the overexpression of Kpnß1 reduced INI-60 cytotoxicity. INI-60 also inhibited KYSE 30 oesophageal cancer cell line growth in vivo. Taken together, these results show that INI-60 inhibits the nuclear import of Kpnß1 cargoes and interferes with cancer cell biology. INI-60 presents as a potential therapeutic approach for cancers of different tissue origins and warrants further investigation as a novel anti-cancer agent.

POT1 stability and binding measured by fluorescence thermal shift assays.

The protein POT1 (Protection of Telomeres 1) is an integral part of the shelterin complex that protects the ends of human chromosomes from degradation or end fusions. It is the only component of shelterin that binds single-stranded DNA. We describe here the application of two separate fluorescent thermal shift assays (FTSA) that provide quantitative biophysical characterization of POT1 stability and its interactions. The first assay uses Sypro Orange™ and monitors the thermal stability of POT1 and its binding under a variety of conditions. This assay is useful for the quality control of POT1 preparations, for biophysical characterization of its DNA binding and, potentially, as an efficient screening tool for binding of small molecule drug candidates. The second assay uses a FRET-labeled human telomeric G-quadruplex structure that reveals the effects of POT1 binding on thermal stability from the DNA frame of reference. These complementary assays provide efficient biophysical approaches for the quantitative characterization of multiple aspects of POT1 structure and function. The results from these assays provide thermodynamics details of POT1 folding, the sequence selectivity of its DNA binding and the thermodynamic profile for its binding to its preferred DNA binding sequence. Most significantly, results from these assays elucidate two mechanisms for the inhibition of POT1 -DNA interactions. The first is by competitive inhibition at the POT1 DNA binding site. The second is indirect and is by stabilization of G-quadruplex formation within the normal POT1 single-stranded DNA sequence to prevent POT1 binding.

The solution structures of higher-order human telomere G-quadruplex multimers.

Human telomeres contain the repeat DNA sequence 5'-d(TTAGGG), with duplex regions that are several kilobases long terminating in a 3' single-stranded overhang. The structure of the single-stranded overhang is not known with certainty, with disparate models proposed in the literature. We report here the results of an integrated structural biology approach that combines small-angle X-ray scattering, circular dichroism (CD), analytical ultracentrifugation, size-exclusion column chromatography and molecular dynamics simulations that provide the most detailed characterization to date of the structure of the telomeric overhang. We find that the single-stranded sequences 5'-d(TTAGGG)n, with n = 8, 12 and 16, fold into multimeric structures containing the maximal number (2, 3 and 4, respectively) of contiguous G4 units with no long gaps between units. The G4 units are a mixture of hybrid-1 and hybrid-2 conformers. In the multimeric structures, G4 units interact, at least transiently, at the interfaces between units to produce distinctive CD signatures. Global fitting of our hydrodynamic and scattering data to a worm-like chain (WLC) model indicates that these multimeric G4 structures are semi-flexible, with a persistence length of ∼34 Å. Investigations of its flexibility using MD simulations reveal stacking, unstacking, and coiling movements, which yield unique sites for drug targeting.

Identification of Small-Molecule Inhibitors Targeting Porphyromonas gingivalis Interspecies Adherence and Determination of Their In Vitro and In Vivo Efficacies.

Porphyromonas gingivalis is one of the primary causative agents of periodontal disease and initially colonizes the oral cavity by adhering to commensal streptococci. Adherence requires the interaction of a minor fimbrial protein (Mfa1) of P. gingivalis with streptococcal antigen I/II (AgI/II). Our previous work identified an AgI/II peptide that potently inhibited adherence and significantly reduced P. gingivalis virulence in vivo, suggesting that this interaction represents a potential target for drug discovery. To develop targeted small-molecule inhibitors of this protein-protein interaction, we performed a virtual screen of the ZINC databases to identify compounds that exhibit structural similarity with the two functional motifs (NITVK and VQDLL) of the AgI/II peptide. Thirty three compounds were tested for in vitro inhibition of P. gingivalis adherence and the three most potent compounds, namely, N7, N17, and V8, were selected for further analysis. The in vivo efficacy of these compounds was evaluated in a murine model of periodontitis. Treatment of mice with each of the compounds significantly reduced maxillary alveolar bone resorption in infected animals. Finally, a series of cytotoxicity tests were performed against human and murine cell lines. Compounds N17 and V8 exhibited no significant cytotoxic activity toward any of the cell lines, whereas compound N7 was cytotoxic at the highest concentrations that were tested (20 and 40 µM). These results identify compounds N17 and V8 as potential lead compounds that will facilitate the design of more potent therapeutic agents that may function to limit or prevent P. gingivalis colonization of the oral cavity.

Identification of a Heme Activation Site on the MD-2/TLR4 Complex.

Myeloid differentiation factor-2 (MD-2) binds lipopolysaccharide (LPS) and initiates toll-like receptor-4 (TLR4) pro-inflammatory signaling. Heme also activates TLR4 signaling, but it is unknown if heme interacts with MD-2. Therefore, we examined MD-2 for a potential heme activation site. Heme-agarose and biotin-heme/streptavidin-agarose pulled down recombinant MD-2, which was inhibited by excess free heme. UV/visible spectroscopy confirmed MD-2-heme binding. To determine whether MD-2 was required for heme-mediated TLR4 signaling, HEK293 cells were transfected with MD-2, TLR4, CD14, and an NF-κB luciferase reporter, and then stimulated with heme or LPS. Heme or LPS treatment elicited robust reporter activity. Absence of MD-2, TLR4 or CD14 plasmid abolished NF-κB reporter responses to heme or LPS. In silico analysis identified two potential heme docking sites on MD-2 near conserved amino acids W23/S33/Y34 and Y36/C37/I44. Heme-induced NF-κB activity was reduced by 39 and 78% in HEK293 cells transfected with MD-2 mutants W23A and Y34A, respectively, compared to WT-MD-2. NF-κB activation by LPS was not affected by the same mutants. Biotinyl-heme/streptavidin-agarose pulled down 68% less W23A and 80% less W23A/S33A/Y34A mutant MD-2 than WT-MD-2. In contrast, at the Y36/C37/I44 MD-2 site, heme-induced NF-κB activity was significantly increased by mutants Y36A (191% of WT-MD-2) and unchanged by mutants C37A and I44A (95 and 92%, respectively, of WT-MD-2). In conclusion, these data suggest that heme binds and activates TLR4 signaling at amino acids W23 and Y34 on MD-2.

Human POT1 unfolds G-quadruplexes by conformational selection.

The reaction mechanism by which the shelterin protein POT1 (Protection of Telomeres 1) unfolds human telomeric G-quadruplex structures is not fully understood. We report here kinetic, thermodynamic, hydrodynamic and computational studies that show that a conformational selection mechanism, in which POT1 binding is coupled to an obligatory unfolding reaction, is the most plausible mechanism. Stopped-flow kinetic and spectroscopic titration studies, along with isothermal calorimetry, were used to show that binding of the single-strand oligonucleotide d[TTAGGGTTAG] to POT1 is both fast (80 ms) and strong (-10.1 ± 0.3 kcal mol-1). In sharp contrast, kinetic studies showed the binding of POT1 to an initially folded 24 nt G-quadruplex structure is four orders of magnitude slower. Fluorescence, circular dichroism and analytical ultracentrifugation studies showed that POT1 binding is coupled to quadruplex unfolding, with a final complex with a stoichiometry of 2 POT1 per 24 nt DNA. The binding isotherm for the POT1-quadruplex interaction was sigmoidal, indicative of a complex reaction. A conformational selection model that includes equilibrium constants for both G-quadruplex unfolding and POT1 binding to the resultant single-strand provided an excellent quantitative fit to the experimental binding data. POT1 unfolded and bound to any conformational form of human telomeric G-quadruplex (antiparallel, hybrid, parallel monomers or a 48 nt sequence with two contiguous quadruplexes), but did not avidly interact with duplex DNA or with other G-quadruplex structures. Finally, molecular dynamics simulations provided a detailed structural model of a 2:1 POT1:DNA complex that is fully consistent with experimental biophysical results.

The hTERT core promoter forms three parallel G-quadruplexes.

The structure of the 68 nt sequence with G-quadruplex forming potential within the hTERT promoter is disputed. One model features a structure with three stacked parallel G-quadruplex units, while another features an unusual duplex hairpin structure adjoined to two stacked parallel and antiparallel quadruplexes. We report here the results of an integrated structural biology study designed to distinguish between these possibilities. As part of our study, we designed a sequence with an optimized hairpin structure and show that its biophysical and biochemical properties are inconsistent with the structure formed by the hTERT wild-type sequence. By using circular dichroism, thermal denaturation, nuclear magnetic resonance spectroscopy, analytical ultracentrifugation, small-angle X-ray scattering, molecular dynamics simulations and a DNase I cleavage assay we found that the wild type hTERT core promoter folds into a stacked, three-parallel G-quadruplex structure. The hairpin structure is inconsistent with all of our experimental data obtained with the wild-type sequence. All-atom models for both structures were constructed using molecular dynamics simulations. These models accurately predicted the experimental hydrodynamic properties measured for each structure. We found with certainty that the wild-type hTERT promoter sequence does not form a hairpin structure in solution, but rather folds into a compact stacked three-G-quadruplex conformation.

Identification of functional domains of the minor fimbrial antigen involved in the interaction of Porphyromonas gingivalis with oral streptococci.

Porphyromonas gingivalis is associated with chronic periodontitis and may initially colonize the oral cavity by adhering to streptococci. Adhesion to streptococci is driven by interaction of the minor fimbrial antigen (Mfa1) with streptococcal antigen I/II. We identified the region of antigen I/II required for this interaction and developed small molecule mimetics that inhibited P. gingivalis adherence. However, the functional motifs of Mfa1 involved in the interaction with antigen I/II remain uncharacterized. A series of N- and C-terminal peptide fragments of Mfa1 were expressed and tested for inhibition of P. gingivalis adherence to S. gordonii. This approach identified residues 225-400 of Mfa1 as essential for P. gingivalis adherence. Using the three-dimensional structure of Mfa1, a putative binding cleft was identified using SiteMap and five small molecule mimetics could dock in this site. Site-specific mutation of residues in the predicted cleft, including R240A, W275A, D321A and A357P inhibited the interaction of Mfa1 with streptococci, whereas mutation of residues not in the predicted cleft (V238A, I252F and ΔK253) had no effect. Complementation of an Mfa1-deficient P. gingivalis strain with wild-type mfa1 restored adherence to streptococci, whereas complementation with full-length mfa1 containing the R240A or A357P mutations did not restore adherence. The mutations did not affect polymerization of Mfa1, suggesting that the complemented strains produced intact minor fimbriae. These results identified specific residues and structural motifs required for the Mfa1-antigen I/II interaction and will facilitate the design of small molecule therapeutics to prevent P. gingivalis colonization of the oral cavity.

Solution structure and functional investigation of human guanylate kinase reveals allosteric networking and a crucial role for the enzyme in cancer.

Human guanylate kinase (hGMPK) is the only known enzyme responsible for cellular GDP production, making it essential for cellular viability and proliferation. Moreover, hGMPK has been assigned a critical role in metabolic activation of antiviral and antineoplastic nucleoside-analog prodrugs. Given that hGMPK is indispensable for producing the nucleotide building blocks of DNA, RNA, and cGMP and that cancer cells possess elevated GTP levels, it is surprising that a detailed structural and functional characterization of hGMPK is lacking. Here, we present the first high-resolution structure of hGMPK in the apo form, determined with NMR spectroscopy. The structure revealed that hGMPK consists of three distinct regions designated as the LID, GMP-binding (GMP-BD), and CORE domains and is in an open configuration that is nucleotide binding-competent. We also demonstrate that nonsynonymous single-nucleotide variants (nsSNVs) of the hGMPK CORE domain distant from the nucleotide-binding site of this domain modulate enzymatic activity without significantly affecting hGMPK's structure. Finally, we show that knocking down the hGMPK gene in lung adenocarcinoma cell lines decreases cellular viability, proliferation, and clonogenic potential while not altering the proliferation of immortalized, noncancerous human peripheral airway cells. Taken together, our results provide an important step toward establishing hGMPK as a potential biomolecular target, from both an orthosteric (ligand-binding sites) and allosteric (location of CORE domain-located nsSNVs) standpoint.

Folding Landscape of a Parallel G-Quadruplex.

Circular dichroism and stopped-flow UV spectroscopies were used to investigate the thermodynamic stability and the folding pathway of d[TGAG3TG3TAG3TG3TA2] at 25 °C in solutions containing 25 mM KCl. Under these conditions the oligonucleotide adopts a thermally stable, all-parallel G-quadruplex topography containing three stacked quartets. K+-induced folding shows three resolved relaxation times, each with distinctive spectral changes. Folding is complete within 200 s. These data indicate a folding pathway that involves at least two populated intermediates, one of which seems to be an antiparallel structure that rearranges to the final all-parallel conformation. Molecular dynamics reveals a stereochemically plausible folding pathway that does not involve complete unfolding of the intermediate. The rate of unfolding was determined using complementary DNA to trap transiently unfolded states to form a stable duplex. As assessed by 1D-1H NMR and fluorescence spectroscopy, unfolding is extremely slow with only one observable rate-limiting relaxation time.

G-quadruplex virtual drug screening: A review.

Over the past two decades biologists and bioinformaticians have unearthed substantial evidence supporting a role for G-quadruplexes as important mediators of biological processes. This includes telomere damage signaling, transcriptional activity, and splicing. Both their structural heterogeneity and their abundance in oncogene promoters makes them ideal targets for drug discovery. Currently, there are hundreds of deposited DNA and RNA quadruplex atomic structures which have allowed researchers to begin using in silico drug screening approaches to develop novel stabilizing ligands. Here we provide a review of the past decade of G-quadruplex virtual drug discovery approaches and campaigns. With this we introduce relevant virtual screening platforms followed by a discussion of best practices to assist future G4 VS campaigns.

Identification of G-quadruplex forming sequences in three manatee papillomaviruses.

The Florida manatee (Trichechus manatus latirotris) is a threatened aquatic mammal in United States coastal waters. Over the past decade, the appearance of papillomavirus-induced lesions and viral papillomatosis in manatees has been a concern for those involved in the management and rehabilitation of this species. To date, three manatee papillomaviruses (TmPVs) have been identified in Florida manatees, one forming cutaneous lesions (TmPV1) and two forming genital lesions (TmPV3 and TmPV4). We identified DNA sequences with the potential to form G-quadruplex structures (G4) across the three genomes. G4 were located on both DNA strands and across coding and non-coding regions on all TmPVs, offering multiple targets for viral control. Although G4 have been identified in several viral genomes, including human PVs, most research has focused on canonical structures comprised of three G-tetrads. In contrast, the vast majority of sequences we identified would allow the formation of non-canonical structures with only two G-tetrads. Our biophysical analysis confirmed the formation of G4 with parallel topology in three such sequences from the E2 region. Two of the structures appear comprised of multiple stacked two G-tetrad structures, perhaps serving to increase structural stability. Computational analysis demonstrated enrichment of G4 sequences on all TmPVs on the reverse strand in the E2/E4 region and on both strands in the L2 region. Several G4 sequences occurred at similar regional locations on all PVs, most notably on the reverse strand in the E2 region. In other cases, G4 were identified at similar regional locations only on PVs forming genital lesions. On all TmPVs, G4 sequences were located in the non-coding region near putative E2 binding sites. Together, these findings suggest that G4 are possible regulatory elements in TmPVs.

A rapid fluorescent indicator displacement assay and principal component/cluster data analysis for determination of ligand-nucleic acid structural selectivity.

We describe a rapid fluorescence indicator displacement assay (R-FID) to evaluate the affinity and the selectivity of compounds binding to different DNA structures. We validated the assay using a library of 30 well-known nucleic acid binders containing a variety chemical scaffolds. We used a combination of principal component analysis and hierarchical clustering analysis to interpret the results obtained. This analysis classified compounds based on selectivity for AT-rich, GC-rich and G4 structures. We used the FID assay as a secondary screen to test the binding selectivity of an additional 20 compounds selected from the NCI Diversity Set III library that were identified as G4 binders using a thermal shift assay. The results showed G4 binding selectivity for only a few of the 20 compounds. Overall, we show that this R-FID assay, coupled with PCA and HCA, provides a useful tool for the discovery of ligands selective for particular nucleic acid structures.

G-Quadruplex Secondary Structure Obtained from Circular Dichroism Spectroscopy.

A curated library of circular dichroism spectra of 23 G-quadruplexes of known structure was built and analyzed. The goal of this study was to use this reference library to develop an algorithm to derive quantitative estimates of the secondary structure content of quadruplexes from their experimental CD spectra. Principal component analysis and singular value decomposition were used to characterize the reference spectral library. CD spectra were successfully fit to obtain estimates of the amounts of base steps in anti-anti, syn-anti or anti-syn conformations, in diagonal or lateral loops, or in other conformations. The results show that CD spectra of nucleic acids can be analyzed to obtain quantitative structural information about secondary structure content in an analogous way to methods used to analyze protein CD spectra.