Conformational polymorphism in telomeric structures: Loop orientation and interloop pairing in d(G4TnG4)

Sequence repeats constituting the telomeric regions of chromosomes are known to adopt a variety of unusual structures, consisting of a G tetraplex stem and short stretches of thymines or thymines and adenines forming loops over the stem. Detailed model building and molecular mechanics studies have been carried out for these telomeric sequences to elucidate different types of loop orientations and possible conformations of thymines in the loop. The model building studies indicate that a minimum of two thymines have to be interspersed between guanine stretches to form folded-back structures with loops across adjacent strands in a G tetraplex (both over the small as well as large groove), while the minimum number of thymines required to build a loop across the diagonal strands in a G tetraplex is three. For two repeat sequences, these hairpins, resulting from different types of folding, can dimerize in three distinct ways—i.e., with loops across adjacent strands and on same side, with loops across adjacent strands and on opposite sides, and with loops across diagonal strands and on opposite sides—to form hairpin dimer structures. Energy minimization studies indicate that all possible hairpin dimers have very similar total energy values, though different structures are stabilized by different types of interactions. When the two loops are on the same side, in the hairpin dimer structures of d(G₄T_nG₄), the thymines form favorably stacked tetrads in the loop region and there is interloop hydrogen bonding involving two hydrogen bonds for each thymine–thymine pair. Our molecular mechanics calculations on various folded-back as well as parallel tetraplex structures of these telomeric sequences provide a theoretical rationale for the experimentally observed feature that the presence of intervening thymine stretches stabilizes folded-back structures, while isolated stretches of guanines adopt a parallel tetraplex structure. © 1994 John Wiley & Sons, Inc.     

Structural Dynamics of Human Telomeric G-Quadruplex Loops Studied by Molecular Dynamics Simulations

Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.     

Hairpin and parallel quartet structures for telomeric sequences.

The role of thymine residues in the formation of G-quartet structures for telomeric sequences has been investigated using model oligonucleotides of the type d(G4TnG4), with n = 1-4. Sequences d(G4T3G4) and d(G4T4G4) adopt a G-quartet structure formed by hairpin dimerization in 70 mM NaCl as judged by a characteristic circular dichroism signature with a 295 nm positive and 265 nm negative bands while d(G4TG4) adopts a parallel G-quartet structure like d(G12) which exhibits a strong positive band at 260 nm and a negative band at 240 nm. The sequence d(G4T2G4) exhibits a mixture of both conformations. The stability of hairpin G-quartet structures decreases with decrease in the number of intervening thymine residues. Potassium permanganate, a single strand specific probe has been used to establish the presence of loops composed of T residues in the hairpin G quartet structures formed by the oligonucleotides d(G4TnG4) with n = 2-4 in 70 mM NaCl. The formation of hairpin G quartet structure for the above sequences is further supported by the enhanced electrophoretic mobility observed on non-denaturing polyacrylamide gels. Human telomeric sequence d(TTAGGG)4 which showed enhanced electrophoretic mobility like Tetrahymena telomeric sequence d(T2G4)4 also exhibited a characteristic CD spectrum for a folded-back G-quartet structure. A detailed model for G-quartet structure involving hairpin dimer with alternating syn-anti-syn-anti conformation for the guanine residues both along the chain as well as around the G tetrad with at least two thymine residues in the loop is proposed. Intermolecular association of short telomeric sequences reported here provides a possible model for chromosomal pairing.     

Molecular dynamics and principal components analysis of human telomeric quadruplex multimers

Guanine-rich DNA repeat sequences located at the terminal ends of chromosomal DNA can fold in a sequence-dependent manner into G-quadruplex structures, notably the terminal 150-200 nucleotides at the 3′ end, which occur as a single-stranded DNA overhang. The crystal structures of quadruplexes with two and four human telomeric repeats show an all-parallel-stranded topology that is readily capable of forming extended stacks of such quadruplex structures, with external TTA loops positioned to potentially interact with other macromolecules. This study reports on possible arrangements for these quadruplex dimers and tetramers, which can be formed from 8 or 16 telomeric DNA repeats, and on a methodology for modeling their interactions with small molecules. A series of computational methods including molecular dynamics, free energy calculations, and principal components analysis have been used to characterize the properties of these higher-order G-quadruplex dimers and tetramers with parallel-stranded topology. The results confirm the stability of the central G-tetrads, the individual quadruplexes, and the resulting multimers. Principal components analysis has been carried out to highlight the dominant motions in these G-quadruplex dimer and multimer structures. The TTA loop is the most flexible part of the model and the overall multimer quadruplex becoming more stable with the addition of further G-tetrads. The addition of a ligand to the model confirms the hypothesis that flat planar chromophores stabilize G-quadruplex structures by making them less flexible.     

Structure of the intramolecular human telomeric G-quadruplex in potassium solution: a novel adenine triple formation

We report the NMR solution structure of the intramolecular G-quadruplex formed in human telomeric DNA in K⁺. The hybrid-type telomeric G-quadruplex consists of three G-tetrads linked with mixed parallel–antiparallel G-strands, with the bottom two G-tetrads having the same G-arrangement (anti:anti:syn:anti) and the top G-tetrad having the reversed G-arrangement (syn:syn:anti:syn). The three TTA loop segments adopt different conformations, with the first TTA assuming a double-chain-reversal loop conformation, and the second and third TTA assuming lateral loop conformations. The NMR structure is very well defined, including the three TTA loops and the two flanking sequences at 5′- and 3′-ends. Our study indicates that the three loop regions interact with the core G-tetrads in a specific way that defines and stabilizes the unique human telomeric G-quadruplex structure in K⁺. Significantly, a novel adenine triple platform is formed with three naturally occurring adenine residues, A21, A3 and A9, capping the top tetrad of the hybrid-type telomeric G-quadruplex. This adenine triple is likely to play an important role in the formation of a stable human telomeric G-quadruplex structure in K⁺. The unique human telomeric G-quadruplex structure formed in K⁺ suggests that it can be specifically targeted for anticancer drug design.     

Conformational dynamics of the human propeller telomeric DNA quadruplex on a microsecond time scale

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. Although G-quartet stems have been well characterized, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 μs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.     

In vitro properties of the conserved mammalian protein hnRNP D suggest a role in telomere maintenance

Mammalian chromosomes terminate with a 3' tail which consists of reiterations of the G-rich repeat, d(TTAGGG). The telomeric tail is the primer for replication by telomerase, and it may also invade telomeric duplex DNA to form terminal lariat structures, or T loops. Here we show that the ubiquitous and highly conserved mammalian protein hnRNP D interacts specifically with the G-rich strand of the telomeric repeat. A single gene encodes multiple isoforms of hnRNP D. All isoforms bind comparably to the G-rich strand, and certain isoforms can also bind tightly and specifically to the C-rich telomeric strand. G-rich telomeric sequences readily form structures stabilized by G-G pairing, which can interfere with telomere replication by telomerase. We show that hnRNP D binding to the G-rich strand destabilizes intrastrand G-G pairing and that hnRNP D interacts specifically with telomerase in human cell extracts. This biochemical analysis suggest that hnRNP D could function in vivo to destabilize structures formed by telomeric G-rich tails and facilitate their extension by telomerase.     

Loop flexibility in human telomeric quadruplex small-molecule complexes

Quadruplex nucleic acids can be formed at the ends of eukaryotic chromosomes. Their formation and stabilisation by appropriate small molecules can be used as a means of inhibiting the telomere maintenance functions of telomerase in human cancer cells. The crystal structures have been determined for a number of complexes between these small molecules and human telomeric DNA and RNA quadruplexes. The detailed structural characteristics of these complexes have been surveyed here and the variations in conformation for the TTA and UUA loops have been explored. Loop conformations have been classified in terms of a number of discrete types and their distribution among the crystal structures. Sugar conformation and backbone angles have also been examined and trends highlighted. One particular loop class has been found to be most prevalent. Implications for in particular, rational drug design, are discussed.     

8-Oxo-7,8-dihydrodeoxyadenosine: the first example of a native DNA lesion that stabilizes human telomeric G-quadruplex DNA

Native DNA lesions in general destabilize DNA secondary structures such as duplex and G-quadruplex because they disrupt optimized interactions in DNA defined by nature. In this paper, we report the first example of a native DNA lesion (8-oxo-7,8-dihydrodeoxyadenosine, OxodA) that stabilizes human telomeric G-quadruplex DNA. CD thermal denaturation studies explicitly displayed increased melting temperatures of telomeric G-quadruplex DNAs that contain OxodA(s) in different DNA loops, suggesting enhanced thermal stability. Conformation studies of G-quadruplex DNAs containing OxodA(s) in the loops using CD and native PAGE revealed that they adopt a similar antiparallel conformation in Na(+) but have much more versatile conformations in K(+). According to computational calculations, the observed stabilization may result from the tight binding of K(+) into the pocket formed by the O8 of OxodA and its loop. The study reported here may provide better understanding of the effect of DNA lesions on G-quadruplex stability and conformation.     

A comparison of four different conformations adopted by human telomeric G‐quadruplex using computer simulations

The telomeric G‐quadruplexes for their unique structural features are considered as potential anticancer drug targets. These, however, exhibit structural polymorphism as different topology types for the intra‐molecular G‐quadruplexes from human telomeric G‐rich sequences have been reported based on NMR spectroscopy and X‐ray crystallography. These techniques provide detailed atomic‐level information about the molecule but relative conformational stability of the different topologies remains unsolved. Therefore, to understand the conformational preference, we have carried out quantum chemical calculations on G‐quartets; used all‐atom molecular dynamics (MD) simulations and steered molecular dynamics (SMD) simulations to characterize the four human telomeric G‐quadruplex topologies based on its G‐tetrad core‐types, viz., parallel, anti‐parallel, mixed‐(3 + 1)‐form1 and mixed‐(3 + 1)‐form2. We have also studied a non‐telomeric sequence along with these telomeric forms giving a comparison between the two G‐rich forms. The structural properties such as base pairing, stacking geometry and backbone conformations have been analyzed. The quantum calculations indicate that presence of a sodium ion inside the G‐tetrad plane or two potassium ions on both sides of the plane give it an overall planarity which is much needed for good stacking to form a helix. MD simulations indicate that capping of the G‐tetrad core by the TTA loops keep the terminal guanine bases away from water. The SMD simulations along with equilibrium MD studies indicate that the parallel and non‐telomeric forms are comparatively less stable. We could come to the conclusion that the anti‐parallel form and also the mixed‐(3 + 1)‐form1 topology are most likely to represent the major conformation., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 83–99, 2016     

RTEL1 dismantles T loops and counteracts telomeric G4-DNA to maintain telomere integrity

T loops and telomeric G-quadruplex (G4) DNA structures pose a potential threat to genome stability and must be dismantled to permit efficient telomere replication. Here we implicate the helicase RTEL1 in the removal of telomeric DNA secondary structures, which is essential for preventing telomere fragility and loss. In the absence of RTEL1, T loops are inappropriately resolved by the SLX4 nuclease complex, resulting in loss of the telomere as a circle. Depleting SLX4 or blocking DNA replication abolished telomere circles (TCs) and rescued telomere loss in RTEL1(-/-) cells but failed to suppress telomere fragility. Conversely, stabilization of telomeric G4-DNA or loss of BLM dramatically enhanced telomere fragility in RTEL1-deficient cells but had no impact on TC formation or telomere loss. We propose that RTEL1 performs two distinct functions at telomeres: it disassembles T loops and also counteracts telomeric G4-DNA structures, which together ensure the dynamics and stability of the telomere.     

Distinct functional contributions of three potential secondary structures in the phage G4 origin of complementary DNA strand synthesis

Three potential secondary structures, stem-loops I, II, and III, are contained in the phage G4 origin of complementary DNA strand synthesis, G4oric, and are believed to be involved in its recognition by dnaG-encoded primase and the synthesis of primer RNA. In a previous publication [Sakai et al., Gene 71 (1988) 323-330], we suggested that base pairing between the loops of stem-loops I, and II, and/or II and III, might play a role in G4oric function. To test this hypothesis, site-directed mutagenesis was used to construct mutants which carried base substitutions in loops I, II and III that destroyed possible interloop base pairing. These mutations, however, did not seriously affect G4oric activity. This indicates that base pairing between the loops is not essential for G4oric functional activity, and also that base substitutions which do not affect the secondary structure of stem-loops I, II and III, do not affect G4oric activity. To complete an analysis of the effects of altering the structure of the G4oric stem-loops, insertions were made into stem-loop III. In contrast to stem-loops I and II, all insertions into stem-loop III destroyed in vivo G4oric activity.     

Structure of a C-rich strand fragment of the human centromeric satellite III: a pH-dependent intercalation topology

Repetitive DNA sequences may adopt unusual pairing arrangements. At acid to neutral pH, cytidine-rich DNA oligodeoxynucleotides can form the i-motif structure in which two parallel-stranded duplexes with C.C(+) pairs are intercalated head-to-tail. The i-motif may be formed by multimeric associations or by intra-molecular folding, depending on the number of cytidine tracts, the nucleotide sequences between them, and the experimental conditions.We have found that a natural fragment of the human centromeric satellite III, d(CCATTCCATTCCTTTCC), can form two monomeric i-motif structures that differ in their intercalation topology and that are favored at pH values higher (the eta-form) and lower (the lambda-form) than 4.6. The change in intercalation may be related to adenine protonation in the loops.We studied the uridine derivative methylated on the first cytidine base, d(5mCCATTCCAUTCCUTTCC), whose proton spectrum is better resolved. The intercalation topologies are (C7.C17)/(5mC1.C11)/(C6.C16)/(C2.C12) for form lambda and (5mC1.C11)/(C7.C17)/(C2.C12)/(C6.C16) for form eta. We have solved the structure of the eta-form, and we present a model for the lambda-form. The switch from eta to lambda involves disruption of the i-motif. In both forms, the central AUT linker crosses the wide groove, and the first and the third linkers loop across the minor grooves. The i-motif core is extended in the eta-form by the inter-loop reverse Watson-Crick A3.U13 pair, whose dissociation constant is around 10(-2) at 0 degrees C, and in the lambda-form by the interloop T5.T15 pair.In contrast, d(5mCCATTCCTTACCTTTCC) folds into a pH-independent structure that has the same intercalation topology as the lambda-form. The i-motif core is extended below by the interloop T5.T15 pair and closed on top by the T8.A10 pair.Thus, the C-rich strand of the human satellite III tandem repeats, like the G-rich strand, can fold into various compact structures. The relevance of these features to centromeric function remains unknown.     

Structural basis for binding of porphyrin to human telomeres

Maintenance of telomere integrity is a hallmark of human cancer, and the single-stranded 3' ends of telomeric DNA are targets for small-molecule anticancer therapies. We report here the crystal structure of a bimolecular human telomeric quadruplex, of the sequence d(TAGGGTTAGGG), in a complex with the quadruplex-binding ligand 5,10,15,20-tetrakis(N-methyl-4-pyridyl)porphyrin (TMPyP4) to a resolution of 2.09 A. The DNA quadruplex topology is parallel-stranded with external double-chain-reversal propeller loops, consistent with previous structural determinations. The porphyrin molecules bind by stacking onto the TTA nucleotides, either as part of the external loop structure or at the 5' region of the stacked quadruplex. This involves stacked on hydrogen-bonded base pairs, formed from those nucleotides not involved in the formation of G-tetrads, and there are thus no direct ligand interactions with G-tetrads. This is in accord with the relative nonselectivity by TMPyP4 for quadruplex DNAs compared to duplex DNA. Porphyrin binding is achieved by remodeling of loops compared to the ligand-free structures. Implications for the design of quadruplex-binding ligands are discussed, together with a model for the formation of anaphase bridges, which are observed following cellular treatment with TMPyP4.     

Distinct second extracellular loop structures of the brain cannabinoid CB(1) receptor: implication in ligand binding and receptor function

The G-protein-coupled receptor (GPCR) second extracellular loop (E2) is known to play an important role in receptor structure and function. The brain cannabinoid (CB(1)) receptor is unique in that it lacks the interloop E2 disulfide linkage to the transmembrane (TM) helical bundle, a characteristic of many GPCRs. Recent mutation studies of the CB(1) receptor, however, suggest the presence of an alternative intraloop disulfide bond between two E2 Cys residues. Considering the oxidation state of these Cys residues, we determine the molecular structures of the 17-residue E2 in the dithiol form (E2(dithiol)) and in the disulfide form (E2(disulfide)) of the CB(1) receptor in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer, using a combination of simulated annealing and molecular dynamics simulation approaches. We characterize the CB(1) receptor models with these two E2 forms, CB(1)(E2(dithiol)) and CB(1)(E2(disulfide)), by analyzing interaction energy, contact number, core crevice, and cross correlation. The results show that the distinct E2 structures interact differently with the TM helical bundle and uniquely modify the TM helical topology, suggesting that E2 of the CB(1) receptor plays a critical role in stabilizing receptor structure, regulating ligand binding, and ultimately modulating receptor activation. Further studies on the role of E2 of the CB(1) receptor are warranted, particularly comparisons of the ligand-bound form with the present ligand-free form.     

Conformation and Stability of Intramolecular Telomeric G-Quadruplexes: Sequence Effects in the Loops

Telomeres are guanine-rich sequences that protect the ends of chromosomes. These regions can fold into G-quadruplex structures and their stabilization by G-quadruplex ligands has been employed as an anticancer strategy. Genetic analysis in human telomeres revealed extensive allelic variation restricted to loop bases, indicating that the variant telomeric sequences maintain the ability to fold into G-quadruplex. To assess the effect of mutations in loop bases on G-quadruplex folding and stability, we performed a comprehensive analysis of mutant telomeric sequences by spectroscopic techniques, molecular dynamics simulations and gel electrophoresis. We found that when the first position in the loop was mutated from T to C or A the resulting structure adopted a less stable antiparallel topology; when the second position was mutated to C or A, lower thermal stability and no evident conformational change were observed; in contrast, substitution of the third position from A to C induced a more stable and original hybrid conformation, while mutation to T did not significantly affect G-quadruplex topology and stability. Our results indicate that allelic variations generate G-quadruplex telomeric structures with variable conformation and stability. This aspect needs to be taken into account when designing new potential anticancer molecules.     

Heterochromatin regions of the chromosomes of the hen and Japanese quail in mitosis and at the lampbrush stage

The mitotic and lampbrush chromosomes of the domestic fowl and Japanese quail were analysed by fluorochrome staining technique. The lampbrush chromosomes of both the subjects displayed a typical "loop-chromomere" structure. Three distinct kinds of loops were distinguished in Gallus g. domesticus--normal, telomeric bows, and lumps. The former are distributed along the whole chromosome length. The latter and the bows were observed in subtelomeric and telomeric regions. By DNA/RNA specific acridine orange staining it was shown that each loop (especially, "lumpy" loops) contained a rich RNP matrix. A comparative analysis of the chromomycin A3/distamycin A banding pattern of mitotic and lampbrush chromosomes shows that the telomeric "bows" and "lumps" are special loops developed in telomeric heterochromatic bands. In Coturnix c. japonica, the CMA/DA-positive bands were not observed in telomeres of mitotic macrochromosomes, except a smallest band in the 2p-arm telomere. The absence of telomeric heterochromatic bands which can be visualized in the quail mitotic chromosomes coincides with the absence of "bow"-like loops. Only small lump-like structures were seen in some telomeres of macroautosomes. The biological significance of loop formation and RNA synthesis in heterochromatic band loops in growing oocytes is briefly discussed.     

Antisense RNA regulation in prokaryotes: rapid RNA/RNA interaction facilitated by a general U-turn loop structure

Efficient gene control by antisense RNA requires rapid bi-molecular interaction with a cognate target RNA. A comparative analysis revealed that a YUNR motif (Y=pyrimidine, R=purine) is ubiquitous in RNA recognition loops in antisense RNA-regulated gene systems. The (Y)UNR sequence motif specifies two intraloop hydrogen bonds forming U-turn structures in many anticodon-loops and all T-loops of tRNAs, the hammerhead ribozyme and in other conserved RNA loops. This structure creates a sharp bend in the RNA phosphate-backbone and presents the following three to four bases in a solvent-exposed, stacked configuration providing a scaffold for rapid interaction with complementary RNA. Sok antisense RNA from plasmid R1 inhibits translation of the hok mRNA by preventing ribosome entry at the mok Shine & Dalgarno element. The 5' single-stranded region of Sok-RNA recognizes a loop in the hok mRNA. We show here, that the initial pairing between Sok antisense RNA and its target in hok mRNA occurs with an observed second-order rate-constant of 2 x 10(6) M(-1) s(-1). Mutations that eliminate the YUNR motif in the target loop of hok mRNA resulted in reduced antisense RNA pairing kinetics, whereas mutations maintaining the YUNR motif were silent. In addition, RNA phosphate-backbone accessibility probing by ethylnitrosourea was consistent with a U-turn structure formation promoted by the YUNR motif. Since the YUNR U-turn motif is present in the recognition units of many antisense/target pairs, the motif is likely to be a generally employed enhancer of RNA pairing rates. This suggestion is consistent with the re-interpretation of the mutational analyses of several antisense control systems including RNAI/RNAII of ColE1, CopA/CopT of R1 and RNA-IN/RNA-OUT of IS10.     

Molecular dynamics simulations of small DNA plasmids: Effects of sequence and supercoiling on intramolecular motions

Small (600 base pair) DNA plasmids were modeled with a simplified representation (3DNA) and the intramolecular motions were studied using molecular mechanics and molecular dynamics techniques. The model is detailed enough to incorporate sequence effects. At the same time, it is simple enough to allow long molecular dynamics simulations. The simulations revealed that large-scale slithering occurs in a homogeneous sequence. In a heterogeneous sequence, containing numerous small intrinsic curves, the centers of the curves are preferentially positioned at the tips of loops. With more curves than loop tips (two in unbranched supercoiled DNA), the heterogeneous sequence plasmid slithers short distances to reposition other curves into the loop tips. However, the DNA is immobilized most of the time, with the loop tips positioned over a few favored curve centers. Branching or looping also appears in the heterogeneous sequence as a new method of repositioning the loop tips. Instead of a smooth progression of increasing writhing with increasing linking difference, theoretical studies have predicted that there is a threshold between unwrithed and writhed DNA at a linking difference between one and two. This has previously been observed in simulations of static structures and is demonstrated here for dynamic homogeneous closed DNA. Such an abrupt transition is not found in the heterogeneous sequence in both the static and dynamic cases. © 1996 John Wiley & Sons, Inc.