Insight into G-DNA structural polymorphism and folding from sequence and loop connectivity through free energy analysis

The lengths of G-tracts and their connecting loop sequences determine G-quadruplex folding and stability. Complete understanding of the sequence-structure relationships remains elusive. Here, single-loop G-quadruplexes were investigated using explicit solvent molecular dynamics (MD) simulations to characterize the effect of loop length, loop sequence, and G-tract length on the folding topologies and stability of G-quadruplexes. Eight loop types, including different variants of lateral, diagonal, and propeller loops, and six different loop sequences [d0 (i.e., no intervening residues in the loop), dT, dT(2), dT(3), dTTA, and dT(4)] were considered through MD simulation and free energy analysis. In most cases the free energetic estimates agree well with the experimental observations. The work also provides new insight into G-quadruplex folding and stability. This includes reporting the observed instability of the left propeller loop, which extends the rules for G-quadruplex folding. We also suggest a plausible explanation why human telomere sequences predominantly form hybrid-I and hybrid-II type structures in K(+) solution. Overall, our calculation results indicate that short loops generally are less stable than longer loops, and we hypothesize that the extreme stability of sequences with very short loops could possibly derive from the formation of parallel multimers. The results suggest that free energy differences, estimated from MD and free energy analysis with current force fields and simulation protocols, are able to complement experiment and to help dissect and explain loop sequence, loop length, and G-tract length and orientation influences on G-quadruplex structure.     

G-quadruplex DNA assemblies: loop length, cation identity, and multimer formation

G-rich DNA sequences are able to fold into structures called G-quadruplexes. To obtain general trends in the influence of loop length on the structure and stability of G-quadruplex structures, we studied oligodeoxynucleotides with random bases in the loops. Sequences studied are dGGGW(i)GGGW(j)GGGW(k)GGG, with W = thymine or adenine with equal probability, and i, j, and k comprised between 1 and 4. All were studied by circular dichroism, native gel electrophoresis, UV-monitored thermal denaturation, and electrospray mass spectrometry, in the presence of 150 mM potassium, sodium, or ammonium cations. Parallel conformations are favored by sequences with short loops, but we also found that sequences with short loops form very stable multimeric quadruplexes, even at low strand concentration. Mass spectrometry reveals the formation of dimers and trimers. When the loop length increases, preferred quadruplex conformations tend to be more intramolecular and antiparallel. The nature of the cation also has an influence on the adopted structures, with K(+) inducing more parallel multimers than NH4(+) and Na(+). Structural possibilities are discussed for the new quadruplex higher-order assemblies.     

Discrete state model and accurate estimation of loop entropy of RNA secondary structures

Conformational entropy makes important contribution to the stability and folding of RNA molecule, but it is challenging to either measure or compute conformational entropy associated with long loops. We develop optimized discrete k-state models of RNA backbone based on known RNA structures for computing entropy of loops, which are modeled as self-avoiding walks. To estimate entropy of hairpin, bulge, internal loop, and multibranch loop of long length (up to 50), we develop an efficient sampling method based on the sequential Monte Carlo principle. Our method considers excluded volume effect. It is general and can be applied to calculating entropy of loops with longer length and arbitrary complexity. For loops of short length, our results are in good agreement with a recent theoretical model and experimental measurement. For long loops, our estimated entropy of hairpin loops is in excellent agreement with the Jacobson-Stockmayer extrapolation model. However, for bulge loops and more complex secondary structures such as internal and multibranch loops, we find that the Jacobson-Stockmayer extrapolation model has large errors. Based on estimated entropy, we have developed empirical formulae for accurate calculation of entropy of long loops in different secondary structures. Our study on the effect of asymmetric size of loops suggest that loop entropy of internal loops is largely determined by the total loop length, and is only marginally affected by the asymmetric size of the two loops. Our finding suggests that the significant asymmetric effects of loop length in internal loops measured by experiments are likely to be partially enthalpic. Our method can be applied to develop improved energy parameters important for studying RNA stability and folding, and for predicting RNA secondary and tertiary structures. The discrete model and the program used to calculate loop entropy can be downloaded at http://gila.bioengr.uic.edu/resources/RNA.html.     

Structural dynamics and cation interactions of DNA quadruplex molecules containing mixed guanine/cytosine quartets revealed by large-scale MD simulations.

Large-scale molecular dynamics (MD) simulations have been utilized to study G-DNA quadruplex molecules containing mixed GCGC and all-guanine GGGG quartet layers. Incorporation of mixed GCGC quartets into G-DNA stems substantially enhances their sequence variability. The mixed quadruplexes form rigid assemblies that require integral monovalent cations for their stabilization. The interaction of cations with the all-guanine quartets is the leading contribution for the stability of the four-stranded assemblies, while the mixed quartets are rather tolerated within the structure. The simulations predict that two cations are preferred to stabilize a four-layer quadruplex stem composed of two GCGC and two all-guanine quartets. The distribution of cations in the structure is influenced by the position of the GCGC quartets within the quadruplex, the presence and arrangement of thymidine loops connecting the guanine/cytosine stretches forming the stems, and the cation type present (Na(+) or K(+)). The simulations identify multiple nanosecond-scale stable arrangements of the thymidine loops present in the molecules investigated. In these thymidine loops, several structured pockets are identified capable of temporarily coordinating cations. However, no stable association of cations to a loop has been observed. The simulations reveal several paths through the thymidine loop regions that can be followed by the cations when exchanging between the central ion channel in the quadruplex stem and the surrounding solvent. We have carried out 20 independent simulations while the length of simulations reaches a total of 90 ns, rendering this study one of the most extensive MD investigations carried out on nucleic acids so far. The trajectories provide a largely converged characterization of the structural dynamics of these four-stranded G-DNA molecules.     

The stability of intramolecular DNA quadruplexes with extended loops forming inter- and intra-loop duplexes

We have examined the formation of intramolecular quadruplex DNA structures in which the loops have been extended so as to generate short DNA duplexes. Fluorescence melting and DNase I cleavage studies show that duplexes can be formed within each loop, but that duplexes between the loops are not stable.     

Time-resolved infrared spectroscopy of thiopeptide isomerization and hydrogen-bond breaking

The photoisomerization of the protected tetrathioxopeptide Boc-Ala-Gly(=S)-Ala-Aib-OMe was followed using time-resolved infrared spectroscopy in the amide I region in combination with isotope labeling. In acetonitrile at room temperature, approximately half of the molecules are found in a loop conformation, restrained by an intramolecular hydrogen bond, while the other half adopts more extended conformations. UV-excitation of the thioxopeptide unit immediately weakens the intramolecular hydrogen bond. After the molecules have relaxed to the electronic ground state with a 130 ps time-constant, a delayed re-formation of the intramolecular hydrogen bond is observed for molecules returning to the initial trans conformation of the thioamide bond, while the loop structure is permanently broken when the molecules isomerize to the cis conformation.     

Visualisation of inter loop tertiary base pairing and stacking interactions in an unusual slipped loop DNA structure

The conformation of an unusual slipped loop DNA structure exhibited by the sequence d(GAATTCCCGAATTC)₂ is determined using a combination of geometrical and molecular mechanics methods. This sequence is known to form a B-DNA-like duplex with the central non-complementary cytosines extruded into single stranded loop regions. The unusual feature is that the interior guanine does not pair with the cytosine across, instead, it pairs with the cytosine upstream by skipping two cytosines, leading to a slipped loop DNA structure with the loops staggered by two base pairs. The two loops, despite being very small, can fold across minor or major groove symmetrically or asymmetrically disposed, with one of the loop bases partially blocking the major or minor groove. Most interestingly, for certain conformations, the loop bases approach one another at close proximity so as to engage even in base pairing as well as base stacking interactions across the major groove. While such pairing and stacking are common in the tertiary folds of RNA, this is the first time that such an interaction is visualized in a DNA. This observation demonstrates that a W-C pair can readily be accomplished in a typical slipped loop structure postulated for DNA. Such tertiary loop interaction may prevent access to regulatory proteins across the major groove of the duplex DNA, thus providing a structure-function relation for the occurrence of slipped loop structure in DNA. Contribution no. 839 from this department     

A proteome-wide analysis of domain architectures of prokaryotic single-spanning transmembrane proteins

We performed a proteome-wide survey of the domain architectures in single-spanning transmembrane (TM) proteins (single-spannings) from 87 sequenced prokaryotic (Bacterial and Archaean) genomes by assigning Pfam domains to their N-tail and C-tail loops. Out of 14,625 single-spannings, 3,516 sequences have at least one domain assigned, and no domains were assigned to 7,850, with the remaining 3,259 with less reliable assignment. In the domain-assigned sequences, 3116 sequences are with at most two domains, and the other 400 sequences with more than two. The assigned domains distribute over 651 Pfam families, which account for 11.4% of the total Pfam-A families. Among the 651 families are mostly soluble-protein-originated ones, but only 21 families are unique to TM proteins. The occurrence frequency of the individual domain families follows a power-law, that is, 264 families occur only once, 106 just twice, and the families appeared more than 30 times are counted by only 39. It is found that the great majority of the sequences having one or two domains are of the type II topology with the C-tail loop containing domains on it. On the contrary, the N-tail loop of the same type topology seldom carries domains. Importantly, the assigned domains are always found on the tail loops longer than 60 residues, even for the small domains with less than 30 residues. There are still as many as 5,800 sequences without assigned domains in spite of having at least one long tail, on which no less than 1,000 novel domain families are expected most likely to lie concealed unknown yet. We also investigated the domain arrangement preference and the domain family combination patterns in 'singlets' (single-spannings with one assigned domain) and 'doublets' (with two domains).     

Expanding the Topological Landscape by a G-Column Flip of a Parallel G-Quadruplex

Canonical G-quadruplexes can adopt a variety of different topologies depending on the arrangement of propeller, lateral, or diagonal loops connecting the four G-columns. A novel intramolecular G-quadruplex structure is derived through inversion of the last G-tract of a three-layered parallel fold, associated with the transition of a single propeller into a lateral loop. The resulting (3+1) hybrid fold features three syn⋅anti⋅anti⋅anti G-tetrads with a 3’-terminal all-syn G-column. Although the ability of forming a duplex stem-loop between G-tracts seems beneficial for a propeller-to-lateral loop rearrangement, unmodified G-rich sequences resist folding into the new (3+1) topology. However, refolding can be driven by the incorporation of syn-favoring guanosine analogues into positions of the fourth G-stretch. The presented hybrid-type G-quadruplex structure as determined by NMR spectroscopy may provide for an additional scaffold in quadruplex-based technologies.     

End-to-end vs interior loop formation kinetics in unfolded polypeptide chains

The conformational search for favorable intramolecular interactions during protein folding is limited by intrachain diffusion processes. Recent studies on the dynamics of loop formation in unfolded polypeptide chains have focused on loops involving residues near the chain ends. During protein folding, however, most contacts are formed between residues in the interior of the chain. We compared the kinetics of end-to-end loop formation (type I loops) to the formation of end-to-interior (type II loops) and interior-to-interior loops (type III loops) using triplet-triplet energy transfer from xanthone to naphthylalanine. The results show that formation of type II and type III loops is slower compared to type I loops of the same size and amino acid sequence. The rate constant for type II loop formation decreases with increasing overall chain dimensions up to a limiting value, at which loop formation is about 2.5-fold slower for type II loops compared to type I loops. Comparing type II loops of different loop size and amino acid sequence shows that the ratio of loop dimension over total chain dimension determines the rate constant for loop formation. Formation of type III loops is 1.7-fold slower than formation of type II loops, indicating that local chain motions are strongly coupled to motions of other chain segments which leads to faster dynamics toward the chain ends. Our results show that differences in the kinetics of formation of type I, type II, and type III loops are mainly caused by differences in internal flexibility at the different positions in the polypeptide chain. Interactions of the polypeptide chain with the solvent contribute to the kinetics of loop formation, which are strongly viscosity-dependent. However, the observed differences in the kinetics of formation of type I, type II, and type III loops are not due to the increased number of peptide-solvent interactions in type II and type III loops compared to type I loops as indicated by identical viscosity dependencies for the kinetics of formation of the different types of loops.     

Topology variation and loop structural homology in crystal and simulated structures of a bimolecular DNA quadruplex

The topology of DNA quadruplexes depends on the nature and number of the nucleotides linking G-quartet motifs. To assess the effects of a three-nucleotide TTT linker, the crystal structure of the DNA sequence d(G(4)T(3)G(4)) has been determined at 1.5 A resolution, together with that of the brominated analogue d(G(4)(Br)UTTG(4)) at 2.4 A resolution. Both sequences form bimolecular intermolecular G-quadruplexes with lateral loops. d(G(4)(Br)UTTG(4)) crystallized in the monoclinic space group P2(1) with three quadruplex molecules in the asymmetric unit, two associating together as a head-to-head stacked dimer, and the third as a single head-to-tail dimer. The head-to-head dimers have two lateral loops on the same G-quadruplex face and form an eight-G-quartet stack, with a linear array of seven K(+) ions between the quartets. d(G(4)T(3)G(4)) crystallized in the orthorhombic space group C222 and has a structure very similar to the head-to-tail dimer in the P2(1) unit cell. The sequence studied here is able to form several different folds; however, all four quadruplexes in the two structures have lateral loops, in contrast to the diagonal loops reported for the analogous quadruplex with T(4) loops. A total of seven independent T(3) loops were observed in the two structures. These can be classified into two discrete conformational classes, suggesting that these represent preferred loop conformations that are independent of crystal-packing forces.     

Response analysis in biochemical chain reactions with negative feedforward and feedbackward loops

Finite difference equations can be used to study the responses of biochemical chain reactions at any step of the chain to an external stimulus. In this study, we developed mathematical models for two hypothetical chain reactions involving loops to study the responses in the chain as the length of the chain gets longer, so called transient and steady state responses. The first model is for a chain with a negative feedforward loop, and the second one is for a chain that has a negative feedback loop. Although both of the models have the same steady state equations and values, we showed that the chain with negative feedforward and negative feedback loops can produce significantly different behaviors. The former can bring the chain into oscillations with various periods and eventually chaos when the feedback is strong enough as the length of the reaction chain increases, whereas the latter is not capable of producing oscillations and more complicated dynamics.     

Structure of an unprecedented G-quadruplex scaffold in the human c-kit promoter

The c-kit oncogene is an important target in the treatment of gastrointestinal tumors. A potential approach to inhibition of the expression of this gene involves selective stabilization of G-quadruplex structures that may be induced to form in the c-kit promoter region. Here we report on the structure of an unprecedented intramolecular G-quadruplex formed by a G-rich sequence in the c-kit promoter in K+ solution. The structure represents a new folding topology with several unique features. Most strikingly, an isolated guanine is involved in G-tetrad core formation, despite the presence of four three-guanine tracts. There are four loops: two single-residue double-chain-reversal loops, a two-residue loop, and a five-residue stem-loop, which contain base-pairing alignments. This unique structural scaffold provides a highly specific platform for the future design of ligands specifically targeted to the promoter DNA of c-kit.     

Monomer–Dimer Equilibrium for the 5′–5′ Stacking of Propeller‐Type Parallel‐Stranded G‐Quadruplexes: NMR Structural Study

Guanine‐rich sequence motifs, which contain tracts of three consecutive guanines connected by single non‐guanine nucleotides, are abundant in the human genome and can form a robust G‐quadruplex structure with high stability. Herein, by using NMR spectroscopy, we investigate the equilibrium between monomeric and 5′–5′ stacked dimeric propeller‐type G‐quadruplexes that are formed by DNA sequences containing GGGT motifs. We show that the monomer–dimer equilibrium depends on a number of parameters, including the DNA concentration, DNA flanking sequences, the concentration and type of cations, and the temperature. We report on the high‐definition structure of a simple monomeric G‐quadruplex containing three single‐residue loops, which could serve as a reference for propeller‐type G‐quadruplex structures in solution.     

Structural basis of telomeric RNA quadruplex--acridine ligand recognition

Human telomeric DNA is now known to be transcribed into noncoding RNA sequences, termed TERRA. These sequences, which are believed to play roles in the regulation of telomere function, can form higher-order quadruplex structures and may themselves be the target of therapeutic intervention. The crystal structure of a TERRA quadruplex-acridine small-molecule complex at a resolution of 2.60 ?, is reported here and contrasts remarkably with the structure of the analogous DNA quadruplex complex. The bimolecular RNA complex has a parallel-stranded topology with propeller-like UUA loops. These loops are held in particular conformations by multiple hydrogen bonds involving the O2' hydroxyl groups of the ribonucleotide sugars and play an active role in binding the acridine molecules to the RNA quadruplex. By contrast, the analogous DNA quadruplex complex has simpler 1:1 acridine binding, with no loop involvement. There are significant loop conformational changes in the RNA quadruplex compared to the native TERRA quadruplex (Collie, G. W.; Haider, S. M.; Neidle, S.; Parkinson, G. N. Nucleic Acids Res. 2010, 38, 5569 - 5580), which have implications for the future design of small molecules targeting TERRA quadruplexes, and RNA quadruplexes more generally.     

Control of the Conductance of Engineered Protein Nanopores through Concerted Loop Motions

Protein nanopores have attracted much interest for nucleic acid sequencing, chemical sensing, and protein folding at the single molecule level. The outer membrane protein OmpG from E. coli stands out because it forms a nanopore from a single polypeptide chain. This property allows the separate engineering of each of the seven extracellular loops that control access to the pore. The longest of these loops, loop 6, has been recognized as the main gating loop that closes the pore at low pH values and opens it at high pH values. A method was devised to pin each of the loops to the embedding membrane and measure the single‐pore conductances of the resulting constructs. The electrophysiological and complementary NMR measurements show that the pinning of individual loops alters the structure and dynamics of neighboring and distant loops in a correlated fashion. Pinning loop 6 generates a constitutively open pore and patterns of concerted loop motions control access to the OmpG nanopore.