An investigation of G-quadruplex structural polymorphism in the human telomere using a combined approach of hydrodynamic bead modeling and molecular dynamics simulation

Probing the Effects of Chemical Modifications on Anticoagulant and Antiproliferative Activity of Thrombin Binding Aptamer

Thrombin binding aptamer (TBA) is one of the best-known G-quadruplex (G4)-forming aptamers that efficiently binds to thrombin, resulting in anticoagulant effects. TBA also possesses promising antiproliferative properties. As with most therapeutic oligonucleotides, chemical modifications are critical for therapeutic applications, particularly to improve thermodynamic stability, resistance in biological environment, and target affinity. To evaluate the effects of nucleobase and/or sugar moiety chemical modifications, five TBA analogues have been designed and synthesized considering that the chair-like G4 structure is crucial for biological activity. Their structural and biological properties have been investigated by Circular Dichroism (CD), Nuclear Magnetic Resonance (NMR), native polyacrylamide gel electrophoresis (PAGE) techniques, and PT and MTT assays. The analogue TBAB contains 8-bromo-2'-deoxyguanosine (B) in G-syn glycosidic positions, while TBAL and TBAM contain locked nucleic acid guanosine (L) or 2'-O-methylguanosine (M) in G-anti positions, respectively. Instead, both the two types of modifications have been introduced in TBABL and TBABM with the aim of obtaining synergistic effects. In fact, both derivatives include B in syn positions, exhibiting in turn L and M in the anti ones. The most appealing results have been obtained for TBABM, which revealed an interesting cytotoxic activity against breast and prostate cancer cell lines, while in the case of TBAB, extraordinary thermal stability (Tm approximately 30 °C higher than that of TBA) and an anticoagulant activity higher than original aptamer were observed, as expected. These data indicate TBAB as the best TBA anticoagulant analogue here investigated and TBABM as a promising antiproliferative derivative.

Strategy for modeling higher-order G-quadruplex structures recalcitrant to NMR determination

Guanine-rich nucleic acids can form intramolecularly folded four-stranded structures known as G-quadruplexes (G4s). Traditionally, G4 research has focused on short, highly modified DNA or RNA sequences that form well-defined homogeneous compact structures. However, the existence of longer sequences with multiple G4 repeats, from proto-oncogene promoters to telomeres, suggests the potential for more complex higher-order structures with multiple G4 units that might offer selective drug-targeting sites for therapeutic development. These larger structures present significant challenges for structural characterization by traditional high-resolution methods like multi-dimensional NMR and X-ray crystallography due to their molecular complexity. To address this current challenge, we have developed an integrated structural biology (ISB) platform, combining experimental and computational methods to determine self-consistent molecular models of higher-order G4s (xG4s). Here we outline our ISB method using two recent examples from our lab, an extended c-Myc promoter and long human telomere G4 repeats, that highlights the utility and generality of our approach to characterizing biologically relevant xG4s.

Effect of molecular crowders on ligand binding kinetics with G-quadruplex DNA probed by fluorescence correlation spectroscopy

Guanine-rich single-stranded DNA folds into G-quadruplex DNA (GqDNA) structures, which play crucial roles in various biological processes. These structures are also promising targets for ligands, potentially inducing antitumor effects. While thermodynamic parameters of ligand/DNA interactions are well-studied, the kinetics of ligand interaction with GqDNA, particularly in cell-like crowded environments, remain less explored. In this study, we investigate the impact of molecular crowding agents (glucose, sucrose, and ficoll 70) at physiologically relevant concentrations (20% w/v) on the association and dissociation rates of the benzophenoxazine-core based ligand, cresyl violet (CV), with human telomeric antiparallel-GqDNA. We utilized fluorescence correlation spectroscopy (FCS) along with other techniques. Our findings reveal that crowding agents decrease the binding affinity of CV to GqDNA, with the most significant effect-a nearly three-fold decrease-observed with ficoll 70. FCS measurements indicate that this decrease is primarily due to a viscosity-induced slowdown of ligand association in the crowded environment. Interestingly, dissociation rates remain largely unaffected by smaller crowders, with only small effect observed in presence of ficoll 70 due to direct but weak interaction between the ligand and ficoll. These results along with previously reported data provide valuable insights into ligand/GqDNA interactions in cellular contexts, suggesting a conserved mechanism of saccharide crowder influence, regardless of variations in GqDNA structure and ligand binding mode. This underscores the importance of considering crowding effects in the design and development of GqDNA-targeted drugs for potential cancer treatment.

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'G_xL_yG_xL_yG_xL_yG_x, 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.

The Binding Pocket at the Interface of Multimeric Telomere G-quadruplexes: Myth or Reality?

Human telomeric DNA with hundreds of repeats of the 5'-TTAGGG-3' motif plays a crucial role in several biological processes. It folds into G-quadruplex (G4) structures and features a pocket at the interface of two contiguous G4 blocks. Up to now no structural NMR and crystallographic data are available for ligands interacting with contiguous G4s. Naphthalene diimide monomers and dyads were investigated as ligands of a dimeric G4 of human telomeric DNA comparing the results with those of the model monomeric G4. Time-resolved fluorescence, circular dichroism, isothermal titration calorimetry and molecular modeling were used to elucidate binding features. Ligand fluorescence lifetime and induced circular dichroism unveiled occupancy of the binding site at the interface. Thermodynamic parameters confirmed the hypothesis as they remarkably change for the dyad complexes of the monomeric and dimeric telomeric G4. The bi-functional ligand structure of the dyads is a fundamental requisite for binding at the G4 interface as only the dyads engage in complexes with 1 : 1 stoichiometry, lodging in the pocket at the interface and establishing multiple interactions with the DNA skeleton. In the absence of NMR and crystallographic data, our study affords important proofs of binding at the interface pocket and clues on the role played by the ligand structure.

Supramolecular Polymorphism of (G4C2)n Repeats Associated with ALS and FTD

Guanine-rich DNA sequences self-assemble into highly stable fourfold structures known as DNA-quadruplexes (or G-quadruplexes). G-quadruplexes have furthermore the tendency to associate into one-dimensional supramolecular aggregates termed G-wires. We studied the formation of G-wires in solutions of the sequences d(G₄C₂)_n with n = 1, 2, and 4. The d(G₄C₂)_n repeats, which are associated with some fatal neurological disorders, especially amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), represent a challenging research topic due to their extensive structural polymorphism. We used dynamic light scattering (DLS) to measure translational diffusion coefficients and consequently resolve the length of the larger aggregates formed in solution. We found that all three sequences assemble into longer structures than previously reported. The d(G₄C₂) formed extremely long G-wires with lengths beyond 80 nm. The d(G₄C₂)₂ formed a relatively short stacked dimeric quadruplex, while d(G₄C₂)₄ formed multimers corresponding to seven stacked intramolecular quadruplexes. Profound differences between the multimerization properties of the investigated sequences were also confirmed by the AFM imaging of surface films. We propose that π-π stacking of the basic G-quadruplex units plays a vital role in the multimerization mechanism, which might be relevant for transformation from the regular medium-length to disease-related long d(G₄C₂)_n repeats.

Analyzing Criteria Affecting the Functionality of G-Quadruplex-Based DNA Aptazymes as Colorimetric Biosensors and Development of Quinine-Binding Aptazymes

A DNA aptazyme consists of an aptamer domain and a DNAzyme module, in which the DNAzyme activity can be regulated by the aptamer-target interaction. The complex of G-quadruplex (GQ) and hemin is a peroxidase-mimicking DNAzyme and has become increasingly popular as a reporter system for biosensing applications. The development of GQ-based aptazymes is of high interest as they can be used as label-free biosensors for the real-time detection of pathogens. Herein, we rationally designed ca. 200 GQ-based aptazyme candidates and evaluated the suitability of 14 aptamers targeting quinine, Protein A, Staphylococcus enterotoxin B, and ATP for this detection concept. As a result, six novel aptazymes were developed for the specific detection of quinine based on two quinine-binding aptamers. The rest of designed probes, however, hardly showed significant functionality. To uncover the reasons, we performed enzyme-linked oligonucleotide assays to find how the affinity of aptamers is affected once conjugated to the DNAzyme sequence or upon integration into the aptazyme probe. Furthermore, we investigated the impact of the structure-switching functionality in the parent aptamer and the effect of the reaction matrix on the efficiency of probes.

Structure of two G-quadruplexes in equilibrium in the KRAS promoter

KRAS is one of the most mutated oncogenes and still considered an undruggable target. An alternative strategy would consist in targeting its gene rather than the protein, specifically the formation of G-quadruplexes (G4) in its promoter. G4 are secondary structures implicated in biological processes, which can be formed among G-rich DNA (or RNA) sequences. Here we have studied the major conformations of the commonly known KRAS 32R, or simply 32R, a 32 residue sequence within the KRAS Nuclease Hypersensitive Element (NHE) region. We have determined the structure of the two major stable conformers that 32R can adopt and which display slow equilibrium (>ms) with each other. By using different biophysical methods, we found that the nucleotides G9, G25, G28 and G32 are particularly implicated in the exchange between these two conformations. We also showed that a triad at the 3' end further stabilizes one of the G4 conformations, while the second conformer remains more flexible and less stable.

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.

Polarizable Molecular Dynamics Simulations of Two c-kit Oncogene Promoter G-Quadruplexes: Effect of Primary and Secondary Structure on Loop and Ion Sampling

G-quadruplexes (GQs) are highly ordered nucleic acid structures that play fundamental roles in regulating gene expression and maintaining genomic stability. GQs are topologically diverse and enriched in promoter sequences of growth regulatory genes and proto-oncogenes, suggesting that they may serve as attractive targets for drug design at the level of transcription rather than inhibiting the activity of the protein products of these genes. The c-kit promoter contains three adjacent GQ-forming sequences that have proposed antagonistic effects on gene expression and thus are promising drug targets for diseases such as gastrointestinal stromal tumors, mast cell disease, and leukemia. Because GQ stability is influenced by primary structure, secondary structure, and ion interactions, a greater understanding of GQ structure, dynamics, and ion binding properties is needed to develop novel, GQ-targeting therapeutics. Here, we performed molecular dynamics simulations to systematically study the c-kit2 and c-kit* GQs, evaluating nonpolarizable and polarizable force fields (FFs) and examining the effects of base substitutions and cation type (K⁺, Na⁺, and Li⁺) on the dynamics of their isolated and linked structures. We found that the Drude polarizable FF outperformed the additive CHARMM36 FF in two- and three-tetrad GQs and solutions of KCl, NaCl, and LiCl. Drude simulations with different cations agreed with the known GQ stabilization preference (K⁺ > Na⁺ > Li⁺) and illustrated that tetrad core-ion coordination differs as a function of cation type. Finally, we showed that differences in primary and secondary structure influence loop sampling, ion binding, and core-ion energetics of GQs.

Exploiting σ-hole interaction to design small uncharged ligand molecules to stabilize G-quadruplex-DNA: a computational study

Small neutral seleno molecules were designed to stabilize G-quadruplex-DNA to accelerate the death of cancer cells. A new approach was considered in this study to design new ligands to stabilize G-quadruplex-DNA using the non-covalent σ-hole interaction. The systematic study has been performed with ligands in the absence and presence of σ-hole interaction to stabilize G-tetrad. Fluorine-substituted seleno ligands interact with the bases of G-quadruplex strongly with the σ-holes present in the ligands. The binding of the fluorinated ligands FSeCF₂SeCF₂ SeF2(4) (~75.0 kcal/mol) is stronger than that of BRACO-19 (~70.0 kcal/mol) calculated at the same level of theory with G-tetrad. It is reported that BRACO-19 is employed to inhibit the enzymatic activity of telomerase. The calculated results also reveal the importance of optimum chain length of ligands to achieve a better binding ability with G-tetrad. The MESP calculations and AIM analysis corroborate the trends of binding of these ligand molecules with G-quadruplexes. The small neutral molecules possess considerable advantage to pass through the lipid bilayer of the cell membrane by passive transport and are of choice in biological research and for clinical trials. Graphical AbtractSmall neutral ligands with σ-holes can stabilize G-quadruplex to inhibit enzymatic activity of telomerase.

Tuning the ring strain effect in acridine derivatives on binding affinity with G-quadruplex-DNA: A computational and experimental study

Search for inhibitors to stabilize the telomeric G-quadruplex in order to deter telomerase activity is an active area of research. Inhibitors play an important role to initiate the tumor cell mortalization process. This work reports for the first time of acridine derivative with four membered ammonium rings at the side chain to surpass the binding ability against BRACO-19 with G-quadruplex-DNA. It is known in the literature that acridine based molecule BRACO-19 can effectively bind with G-quadruplex-DNA. The computational study performed in this study revealed that the binding ability of acridine based molecule can be augmented with subtle variation in the molecular structure of the drug like candidates. Steered molecular dynamics (SMD) performed with the acridine derivatives and G-quadruplex DNA showed the importance of ring strain to the side chain of those ligand molecules. The rupture force analysis, hydrogen bonding interactions and the calculated free energies in MM-PBSA method suggest that ligand 3 is superior than that of BRACO-19. The synthesized ligand 3 and BRACO-19 showed the binding constants obtained from ITC measurements are 4 × 10⁶ mol^-1 and 2.6 × 10⁶, which corroborates the computational findings.

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.

Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes

G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]₄ GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]₄ simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na⁺ ions move inside the GQs in a concerted manner, while larger relocations of the K⁺ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K⁺ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.

The effects of monovalent metal ions on the conformation of human telomere DNA using analytical ultracentrifugation

A human telomere DNA segment (HT-DNA) can fold into a G-quadruplex in the presence of some monovalent cations. These cations can interact with the phosphate groups of the DNA segment and/or with the O6 oxygen atom of guanines, which are called non-specific interactions and specific interactions, respectively. However, until now how these two interactions affect the structure of HT-DNA has not been well understood. In this study, a combination of analytical ultracentrifugation (AUC) and circular dichroism (CD) was used to explore the effects of these two interactions on the structure of a 22-mer single-stranded DNA with a sequence of 5'-AGGG(TTAGGG)3-3'. The results showed that the standard sedimentation coefficient (s20,w) of HT-DNA starts to increase when the concentration of potassium ions (CK(+)) is higher than 10.0 µM due to the formation of a G-quadruplex through specific interactions. Whereas, for a control DNA, a higher CK(+) value of 1.0 mM was needed for increasing s20,w due to non-specific interactions. Moreover, potassium ions could promote the formation of the G-quadruplex much more easily than lithium, sodium and cesium ions, presumably due to its appropriate size in the dehydrated state and easier dehydration. The molar mass of DNA at different cation concentrations was nearly a constant and close to the theoretical value of the molar mass of monomeric HT-DNA, indicating that what we observed is the structural change of individual DNA chains.

Role of Alkali Metal Ions in G-Quadruplex Nucleic Acid Structure and Stability

G-quadruplexes are guanine-rich nucleic acids that fold by forming successive quartets of guanines (the G-tetrads), stabilized by intra-quartet hydrogen bonds, inter-quartet stacking, and cation coordination. This specific although highly polymorphic type of secondary structure deviates significantly from the classical B-DNA duplex. G-quadruplexes are detectable in human cells and are strongly suspected to be involved in a number of biological processes at the DNA and RNA levels. The vast structural polymorphism exhibited by G-quadruplexes, together with their putative biological relevance, makes them attractive therapeutic targets compared to canonical duplex DNA. This chapter focuses on the essential and specific coordination of alkali metal cations by G-quadruplex nucleic acids, and most notably on studies highlighting cation-dependent dissimilarities in their stability, structure, formation, and interconversion. Section 1 surveys G-quadruplex structures and their interactions with alkali metal ions while Section 2 presents analytical methods used to study G-quadruplexes. The influence of alkali cations on the stability, structure, and kinetics of formation of G-quadruplex structures of quadruplexes will be discussed in Sections 3 and 4. Section 5 focuses on the cation-induced interconversion of G-quadruplex structures. In Sections 3 to 5, we will particularly emphasize the comparisons between cations, most often K(+) and Na(+) because of their prevalence in the literature and in cells.

Folding dynamics and conformational heterogeneity of human telomeric G-quadruplex structures in Na+ solutions by single molecule FRET microscopy

G-quadruplex structures can occur throughout the genome, including at telomeres. They are involved in cellular regulation and are potential drug targets. Human telomeric G-quadruplex structures can fold into a number of different conformations and show large conformational diversity. To elucidate the different G-quadruplex conformations and their dynamics, we investigated telomeric G-quadruplex folding using single molecule FRET microscopy in conditions where it was previously believed to yield low structural heterogeneity. We observed four FRET states in Na⁺ buffers: an unfolded state and three G-quadruplex related states that can interconvert between each other. Several of these states were almost equally populated at low to medium salt concentrations. These observations appear surprising as previous studies reported primarily one G-quadruplex conformation in Na⁺ buffers. Our results permit, through the analysis of the dynamics of the different observed states, the identification of a more stable G-quadruplex conformation and two transient G-quadruplex states. Importantly these results offer a unique view into G-quadruplex topological heterogeneity and conformational dynamics.

Site specific replacements of a single loop nucleoside with a dibenzyl linker may switch the activity of TBA from anticoagulant to antiproliferative

Many antiproliferative G-quadruplexes (G4s) arise from the folding of GT-rich strands. Among these, the Thrombin Binding Aptamer (TBA), as a rare example, adopts a monomolecular well-defined G4 structure. Nevertheless, the potential anticancer properties of TBA are severely hampered by its anticoagulant action and, consequently, no related studies have appeared so far in the literature. We wish to report here that suitable chemical modifications in the TBA sequence can preserve its antiproliferative over anticoagulant activity. Particularly, we replaced one residue of the TT or TGT loops with a dibenzyl linker to develop seven new quadruplex-forming TBA based sequences (TBA-bs), which were studied for their structural (CD, CD melting, 1D NMR) and biological (fibrinogen, PT and MTT assays) properties. The three-dimensional structures of the TBA-bs modified at T13 (TBA-bs13) or T12 (TBA-bs12), the former endowed with selective antiproliferative activity, and the latter acting as potently as TBA in both coagulation and MTT assays, were further studied by 2D NMR restrained molecular mechanics. The comparative structural analyses indicated that neither the stability, nor the topology of the G4s, but the different localization of the two benzene rings of the linker was responsible for the loss of the antithrombin activity for TBA-bs13.

Extended molecular dynamics of a c-kit promoter quadruplex

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

Hydrodynamic Models of G-Quadruplex Structures

G-quadruplexes are noncannonical four-stranded DNA or RNA structures formed by guanine-rich repeating sequences. Guanine nucleotides can hydrogen bond to form a planar tetrad structure. Such tetrads can stack to form quadruplexes of various molecularities with a variety of types of single-stranded loops joining the tetrads. High-resolution structures may be obtained by X-ray crystallography or NMR spectroscopy for quadruplexes formed by short (≈25 nt) sequences but these methods have yet to succeed in characterizing higher order quadruplex structures formed by longer sequences. An integrated computational and experimental approach was implemented in our laboratory to obtain structural models for higher order quadruplexes that might form in longer telomeric or promoter sequences. In our approach, atomic-level models are built using folding principles gleaned from available high-resolution structures and then optimized by molecular dynamics. The program HYDROPRO is then used to construct bead models of these structures to predict experimentally testable hydrodynamic properties. Models are validated by comparison of these properties with measured experimental values obtained by analytical ultracentrifugation or other biophysical tools. This chapter describes our approach and practical procedures.

Separation of Quadruplex Polymorphism in DNA Sequences by Reversed-Phase Chromatography

This unit describes a method for the separation of a mixture of quadruplex conformations formed from the same parent sequence via reversed-phase chromatography (RPC). Polymorphism is inherent to quadruplex formation and even relatively simple quadruplex-forming sequences can fold into a cornucopia of possible conformations and topologies. Isolation of a specific conformation for study can be problematic. This is especially true for conformations of the human telomere sequence d(GGG(TTAGGG)3). High performance liquid chromatography (HPLC), especially reversed-phase chromatography, has been a mainstay of nucleic acid research and purification for many decades. We have successfully applied this method to the problem of separating individual quadruplex species in the ensemble from the same parent sequence.