Folding and Hydrodynamics of a DNA i-Motif from the c-MYC Promoter Determined by Fluorescent Cytidine Analogs

The four-stranded i-motif (iM) conformation of cytosine-rich DNA has importance to a wide variety of biochemical systems that range from their use in nanomaterials to potential roles in oncogene regulation. The iM structure is formed at slightly acidic pH, where hemiprotonation of cytosine results in a stable C-C⁺ basepair. Here, we performed fundamental studies to examine iM formation from a C-rich strand from the promoter of the human c-MYC gene. We used a number of biophysical techniques to characterize both the hydrodynamic properties and folding kinetics of a folded iM. Our hydrodynamic studies using fluorescence anisotropy decay and analytical ultracentrifugation show that the iM structure has a compact size in solution and displays the rigidity of a double strand. By studying the rates of circular dichroism spectral changes and quenching of fluorescent cytidine analogs, we also established a mechanism for the folding of a random coil oligo into the iM. In the course of determining this folding pathway, we established that the fluorescent dC analogs tC° and PdC can be used to monitor individual residues of an iM structure and to determine the pK_a of an iM. We established that the C-C⁺ hydrogen bonding of certain bases initiates the folding of the iM structure. We also showed that substitutions in the loop regions of iMs give a distinctly different kinetic signature during folding compared with bases that are intercalated. Our data reveal that the iM passes through a distinct intermediate form between the unfolded and folded forms. Taken together, our results lay the foundation for using fluorescent dC analogs to follow structural changes during iM formation. Our technique may also be useful for examining folding and structural changes in more complex iMs.     

Hysteresis in poly‐2′‐deoxycytidine i‐motif folding is impacted by the method of analysis as well as loop and stem lengths

In DNA, i‐motif (iM) folds occur under slightly acidic conditions when sequences rich in 2′‐deoxycytidine (dC) nucleotides adopt consecutive dC self base pairs. The pH stability of an iM is defined by the midpoint in the pH transition (pH_T) between the folded and unfolded states. Two different experiments to determine pH_T values via circular dichroism (CD) spectroscopy were performed on poly‐dC iMs of length 15, 19, or 23 nucleotides. These experiments demonstrate two points: (1) pH_T values were dependent on the titration experiment performed, and (2) pH‐induced denaturing or annealing processes produced isothermal hysteresis in the pH_T values. These results in tandem with model iMs with judicious mutations of dC to thymidine to favor particular folds found the hysteresis was maximal for the shorter poly‐dC iMs and those with an even number of base pairs, while the hysteresis was minimal for longer poly‐dC iMs and those with an odd number of base pairs. Experiments to follow the iM folding via thermal changes identified thermal hysteresis between the denaturing and annealing cycles. Similar trends were found to those observed in the CD experiments. The results demonstrate that the method of iM analysis can impact the pH_T parameter measured, and hysteresis was observed in the pH_T and T_m values.     

Sensitivity of Superfolder GFP to Ionic Agents

Superfolder variant of the green fluorescent protein (sfGFP) became a favorite probe for examination of the unfolding–refolding processes of fluorescent proteins with beta-barrel structure owing to its reversible unfolding in comparison with other fluorescent proteins. Its benefit is the proper folding even in fusion constructions with poorly folded polypeptides. We noticed that guanidine thiocyanate affects not only the structure of protein but its chromophore directly. Therefore we studied the influence of ionic denaturants and salts including guanidine thiocyanate, guanidine hydrochloride, sodium chloride and sodium thiocyanate on spectral features of sfGFP. It was shown that moderate amounts of the studied agents do not disrupt sfGFP structure but provoke pronounced alteration of its spectral characteristics. Changes in absorption and CD spectra in visible spectral range indicate the specific binding of SCN⁻ and Cl⁻ anions in the sfGFP chromophore vicinity. The anion binding results in the redistribution of sfGFP molecules with neutral and anionic chromophores. This also hinders the proton transfer in the chromophore excited state, considerably decreasing the fluorescence intensity of sfGFP. Our results indicate that when ionic denaturants are used in the studies of fluorescent protein folding their effect on fluorophore charge state should be taken into account.     

Salt-Dependent Folding Energy Landscape of RNA Three-Way Junction

RNAs are highly negatively charged chain molecules. Metal ions play a crucial role in RNA folding stability and conformational changes. In this work, we employ the recently developed tightly bound ion (TBI) model, which accounts for the correlation between ions and the fluctuation of ion distributions, to investigate the ion-dependent free energy landscape for the three-way RNA junction in a 16S rRNA domain. The predicted electrostatic free energy landscape suggests that 1), ion-mediated electrostatic interactions cause an ensemble of unfolded conformations narrowly populated around the maximally extended structure; and 2), Mg²⁺ ion-induced correlation effects help bring the helices to the folded state. Nonelectrostatic interactions, such as noncanonical interactions within the junctions and between junctions and helix stems, might further limit the conformational diversity of the unfolded state, resulting in a more ordered unfolded state than the one predicted from the electrostatic effect. Moreover, the folded state is predominantly stabilized by the coaxial stacking force. The TBI-predicted folding stability agrees well with the experimental measurements for the different Na⁺ and Mg²⁺ ion concentrations. For Mg²⁺ solutions, the TBI model, which accounts for the Mg²⁺ ion correlation effect, gives more improved predictions than the Poisson-Boltzmann theory, which tends to underestimate the role of Mg²⁺ in stabilizing the folded structure. Detailed control tests indicate that the dominant ion correlation effect comes from the charge-charge Coulombic correlation rather than the size (excluded volume) correlation between the ions. Furthermore, the model gives quantitative predictions for the ion size effect in the folding energy landscape and folding cooperativity.     

''Hairpin'' catalytic RNA model: evidence for helices and sequence requirement for substrate RNA.

We have identified the catalytic domain within the sequence of the negative strand of the satellite RNA of tobacco ringspot virus. Minimum energy RNA folding calculations predict a two dimensional model with four major helical regions which are supported by mutagenesis experiments. This model for the catalytic complex consists of a 50 base catalytic RNA and a 14 base substrate RNA folded together in a type of hairpin two dimensional structure. Part of the recognition region between the catalyst and substrate is two helices of 6 bases and 4 bases respectively. Catalytic activity remains when the bases in these two helices are changed but base pairing is maintained. Thus an appropriately engineered 'hairpin' catalyst is capable of cleaving heterologous RNA.     

The chicken-egg scenario of protein folding revisited

What is the first step in protein folding - hydrophobic collapse (compaction) or secondary structure formation? It is still not clear if the major driving force in protein folding is hydrogen bonding or hydrophobic interactions or both. We analyzed data on the conformational characteristics of 41 globular proteins in native and partially folded conformational states. Our analysis shows that a good correlation exists between relative decrease in hydrodynamic volume and increase in secondary structure content. No compact equilibrium intermediates lacking secondary structure, or highly ordered non-compact species, were found. This correlation provides experimental support for the hypothesis that hydrophobic collapse occurs simultaneously with formation of secondary structure in the early stages of the protein folding.     

Single-Molecule TPM Studies on the Conversion of Human Telomeric DNA

Human telomere contains guanine-rich (G-rich) tandem repeats of single-stranded DNA sequences at its 3′ tail. The G-rich sequences can be folded into various secondary structures, termed G-quadruplexes (G4s), by Hoogsteen basepairing in the presence of monovalent cations (such as Na⁺ and K⁺). We developed a single-molecule tethered particle motion (TPM) method to investigate the unfolding process of G4s in the human telomeric sequence AGGG(TTAGGG)₃ in real time. The TPM method monitors the DNA tether length change caused by formation of the G4, thus allowing the unfolding process and structural conversion to be monitored at the single-molecule level. In the presence of its antisense sequence, the folded G4 structure can be disrupted and converted to the unfolded conformation, with apparent unfolding time constants of 82 s and 3152 s. We also observed that the stability of the G4 is greatly affected by different monovalent cations. The folding equilibrium constant of G4 is strongly dependent on the salt concentration, ranging from 1.75 at 5 mM Na⁺ to 3.40 at 15 mM Na⁺. Earlier spectral studies of Na⁺- and K⁺-folded states suggested that the spectral conversion between these two different folded structures may go through a structurally unfolded intermediate state. However, our single-molecule TPM experiments did not detect any totally unfolded intermediate within our experimental resolution when sodium-folded G4 DNA molecules were titrated with high-concentration, excess potassium ions. This observation suggests that a totally unfolding pathway is likely not the major pathway for spectral conversion on the timescale of minutes, and that interconversion among folded states can be achieved by the loop rearrangement. This study also demonstrates that TPM experiments can be used to study conformational changes in single-stranded DNA molecules.     

In vitro synthesis of uniform poly(dG)-poly(dC) by Klenow exo- fragment of polymerase I

In this paper, we describe a production procedure of the one-to-one double helical complex of poly(dG)–poly(dC), characterized by a well-defined length (up to 10 kb) and narrow size distribution of molecules. Direct evidence of strands slippage during poly(dG)–poly(dC) synthesis by Klenow exo⁻ fragment of polymerase I is obtained by fluorescence resonance energy transfer (FRET). We show that the polymer extension results in an increase in the separation distance between fluorescent dyes attached to 5′ ends of the strands in time and, as a result, losing communication between the dyes via FRET. Analysis of the products of the early steps of the synthesis by high-performance liquid chromatography and mass spectroscopy suggest that only one nucleotide is added to each of the strand composing poly(dG)–poly(dC) in the elementary step of the polymer extension. We show that proper pairing of a base at the 3′ end of the primer strand with a base in sequence of the template strand is required for initiation of the synthesis. If the 3′ end nucleotide in either poly(dG) or poly(dC) strand is substituted for A, the polymer does not grow. Introduction of the T-nucleotide into the complementary strand to permit pairing with A-nucleotide results in the restoration of the synthesis. The data reported here correspond with a slippage model of replication, which includes the formation of loops on the 3′ ends of both strands composing poly(dG)–poly(dC) and their migration over long-molecular distances (μm) to 5′ ends of the strands.     

Dynamics and stability of polymorphic human telomeric G-quadruplex under tension

As critical DNA structures capping the human chromosome ends, the stability and structural polymorphism of human telomeric G-quadruplex (G4) have drawn increasing attention in recent years. This work characterizes the equilibrium transitions of single-molecule telomeric G4 at physiological K⁺ concentration. We report three folded states of telomeric G4 with markedly different lifetime and mechanical stability. Our results show that the kinetically favored folding pathway is through a short-lived intermediate state to a longer-lived state. By examining the force dependence of transition rates, the force-dependent transition free energy landscape for this pathway is determined. In addition, an ultra-long-lived form of telomeric G4 structure with a much stronger mechanical stability is identified.     

Folding of the multidomain human immunodeficiency virus type-I integrase.

Protein folding conditions were established for human immunodeficiency virus integrase (IN) obtained from purified bacterial inclusion bodies. IN was denatured by 6 M guanidine.HCl-5 mM dithiothreitol, purified by gel filtration, and precipitated by ammonium sulfate. The reversible solvation of precipitated IN by 6 M guanidine.HCl allowed for wide variation of protein concentration in the folding reaction. A 6-fold dilution of denatured IN by 1 M NaCl buffer followed by dialysis produced enzymatically active IN capable of 3' OH end processing, strand transfer, and disintegration using various human immunodeficiency virus-1 (HIV-1) long terminal repeat DNA substrates. The specific activities of folded IN preparations for these enzymatic reactions were comparable to those of soluble IN purified directly from bacteria. The subunit composition and enzymatic activities of IN were affected by the folding conditions. Standard folding conditions were defined in which monomers and protein aggregates sedimenting as dimers and tetramers wree produced. These protein aggregates were enzymatically active, whereas monomers had reduced strand transfer activity. Temperature modifications of the folding conditions permitted formation of mainly monomers. Upon assaying, these monomers were efficient for strand transfer and disintegration, but the oligomeric state of IN under the conditions of the assay is determinate. Our results suggest that monomers of the multidomain HIV-1 IN are folded correctly for various catalytic activities, but the conditions for specific oligomerization in the absence of catalytic activity are undefined.     

Characterizing intermolecular interactions that initiate native-like protein aggregation

Folded proteins can access aggregation-prone states without the need for transitions that cross the energy barriers for unfolding. In this study we characterized the initial steps of aggregation from a native-like state of the acylphosphatase from Sulfolobus solfataricus (Sso AcP). Using computer simulations restrained by experimental hydrogen/deuterium (H/D) exchange data, we provide direct evidence that under aggregation-promoting conditions Sso AcP populates a conformational ensemble in which native-like structure is retained throughout the sequence in the absence of local unfolding (N¹), although the protein exhibits an increase in hydrodynamic radius and dynamics. This transition leads an edge strand to experience an increased affinity for a specific unfolded segment of the protein. Direct measurements by means of H/D exchange rates, isothermal titration calorimetry, and intermolecular relaxation enhancements show that after formation of N¹, an intermolecular interaction with an antiparallel arrangement is established between the edge strand and the unfolded segment of the protein. However, under conditions that favor the fully native state of Sso AcP, such an interaction is not established. Thus, these results reveal a novel (to our knowledge) self-assembly mechanism for a folded protein that is based on the increased flexibility of highly aggregation-prone segments in the absence of local unfolding.     

NMR structural study on the self-trimerization of d(GTTAGG) into a dynamic trimolecular G-quadruplex assembly preferentially in Na+ solution with a moderate K+ tolerance

Vast G-quadruplexes (GQs) are primarily folded by one, two, or four G-rich oligomers, rarely with an exception. Here, we present the first NMR solution structure of a trimolecular GQ (tri-GQ) that is solely assembled by the self-trimerization of d(GTTAGG), preferentially in Na⁺ solution tolerant to an equal amount of K⁺ cation. Eight guanines from three asymmetrically folded strands of d(GTTAGG) are organized into a two-tetrad core, which features a broken G-column and two width-irregular grooves. Fast strand exchanges on a timescale of second at 17°C spontaneously occur between folded tri-GQ and unfolded single-strand of d(GTTAGG) that both species coexist in dynamic equilibrium. Thus, this tri-GQ is not just simply a static assembly but rather a dynamic assembly. Moreover, another minor tetra-GQ that has putatively tetrameric (2+2) antiparallel topology becomes noticeable only at an extremely high strand concentration above 18 mM. The major tri-GQ and minor tetra-GQ are considered to be mutually related, and their reversible interconversion pathways are proposed accordingly. The sequence d(GTTAGG) could be regarded as either a reading frame shifted single repeat of human telomeric DNA or a 1.5 repeat of Bombyx mori telomeric DNA. Overall, our findings provide new insight into GQs and expect more functional applications.     

Synthesis and antitumor activity evaluation of a novel combi-nitrosourea prodrug: Designed to release a DNA cross-linking agent and an inhibitor of O6-alkylguanine-DNA alkyltransferase

The drug resistance of CENUs induced by O⁶-alkylguanine-DNA alkyltransferase (AGT), which repairs the O⁶-alkylated guanine and subsequently inhibits the formation of dG–dC cross-links, hinders the application of CENU chemotherapies. Therefore, the discovery of CENU analogs with AGT inhibiting activity is a promising approach leading to novel CENU chemotherapies with high therapeutic index. In this study, a new combi-nitrosourea prodrug 3-(3-(((2-amino-9H-purin-6-yl)oxy)methyl)benzyl)-1-(2-chloroethyl)-1-nitrosourea (6), designed to release a DNA cross-linking agent and an inhibitor of AGT, was synthesized and evaluated for its antitumor activity and ability to induce DNA interstrand cross-links (ICLs). The results indicated that 6 exhibited higher cytotoxicity against mer⁺ glioma cells compared with ACNU, BCNU, and their respective combinations with O⁶-benzylguanine (O⁶-BG). Quantifications of dG–dC cross-links induced by 6 were performed using HPLC–ESI-MS/MS. Higher levels of dG–dC cross-link were observed in 6-treated human glioma SF763 cells (mer⁺), whereas lower levels of dG–dC cross-link were observed in 6-treated calf thymus DNA, when compared with the groups treated with BCNU and ACNU. The results suggested that the superiority of 6 might result from the AGT inhibitory moiety, which specifically functions in cells with AGT activity. Molecular docking studies indicated that five hydrogen bonds were formed between the O⁶-BG analogs released from 6 and the five residues in the active pocket of AGT, which provided a reasonable explanation for the higher AGT-inhibitory activity of 6 than O⁶-BG.     

Distinct differences in metal ion specificity of RNA and DNA G-quadruplexes

RNA G-quadruplexes, as their well-studied DNA analogs, require the presence of cations to fold and remain stable. This is the first comprehensive study on the interaction of RNA quadruplexes with metal ions. We investigated the formation and stability of two highly conserved and biologically relevant RNA quadruplex-forming sequences (24nt-TERRA and 18nt-NRAS) in the presence of several monovalent and divalent metal ions, namely Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺, NH₄ ⁺, Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺. Circular dichroism was used to probe the influence of these metal ions on the folded fraction of the parallel G-quadruplexes, and UV thermal melting experiments allowed to assess the relative stability of the structures in each cationic condition. Our results show that the RNA quadruplexes are more stable than their DNA counterparts under the same buffer conditions. We have observed that the addition of mainly Na⁺, K⁺, Rb⁺, NH₄ ⁺, as well as Sr²⁺ and Ba²⁺ in water, shifts the equilibrium to the folded quadruplex form, whereby the NRAS sequence responds stronger than TERRA. However, only K⁺ and Sr²⁺ lead to a significant increase in the stability of the folded structures, which is consistent with their coordination to the O6 atoms from the G-quartet guanosines. Compared to the respective DNA motives, dNRAS and htelo, the RNA sequences are not stabilized by Na⁺ ions. Finally, the difference in response between NRAS and TERRA, as well as to the corresponding DNA sequences with respect to different metal ions, could potentially be exploited for selective targeting purposes.     

A-like guanine-guanine stacking in the aqueous DNA duplex of d(GGGGCCCC)

We have used CD spectroscopy, NMR spectroscopy and unrestrained molecular dynamics to study conformational properties of a DNA duplex formed by the self-complementary octamer d(GGGGCCCC). Its unusual CD spectrum contains features indicating A-like stacking of half of the bases, whereas the other half stack in a B-like fashion. Unrestrained molecular dynamics simulations converged to a stable B-like double-helix of d(GGGGCCCC). However, the double-helix contained a central hole whose size was half of that occurring in structure A. In the canonical structure B, the hole does not exist at all because the base-pairs cross the double-helix centre. The cytosine bases were stacked in the duplex of d(GGGGCCCC) as in structure B, while stacking of the guanine bases displayed features characteristic for structure A. NMR spectroscopy revealed that the A-like guanine-guanine stacking was accompanied by an increased tendency of the deoxyribose rings attached to the guanine bases to be puckered in an A-like fashion. Otherwise, the duplex of d(GGGGCCCC) showed no clash, no bend and no other significant deviation from structure B. The present analysis demonstrates a remarkable propensity of the guanine runs to stack in an A-like fashion even within the B-DNA framework. This property explains why the oligo(dG). oligo(dC) tracts switch into structure A so easily. Secondly, this property may influence replication, because structure A is replicated more faithfully than structure B. Thirdly, the oligo(dG) runs might have played an important role in early evolution, when DNA took on functions that originally evolved on RNA. Fourthly, the present study extends the vocabulary of DNA secondary structures by the heteronomous duplex of d(GGGGCCCC) in which the B-like strand of oligo(dC) is bound to the A-like strand of oligo(dG).     

Modulating protein folding rates in vivo and in vitro by side-chain interactions between the parallel beta strands of green fluorescent protein

We have identified pairs of residues across the two parallel beta strands of green fluorescent protein that facilitate native strand register of the surface-exposed beta barrel. After constructing a suitable host environment around two guest residues, minimizing interactions of the guest residues with surrounding side-chains yet maintaining the wild-type protein structure and the chromophore environment, we introduced a library of cross-strand pairings by cassette mutagenesis. Colonies of Escherichia coli transformed with the library differ in intracellular fluorescence. Most of the fluorescent pairs have predominantly charged and polar guest site residues. The magnitude and the rate of fluorescence acquisition in vivo from transformed E. coli cells varies among the mutants despite comparable levels of protein expression. Spectroscopic measurements of purified mutants show that the native protein structure is maintained. Kinetic studies using purified protein with fully matured chromophores demonstrate that the mutants span a 10-fold range in folding rates with undetectable differences in unfolding rates. Thus, green fluorescent protein provides an ideal system for monitoring determinants of in vivo protein folding. Cross-strand pairings affect both protein stability and folding kinetics by favoring the formation of native strand register preferentially to non-native strand alignments.     

Displacement synthesis and cloning of symmetrical complementary DNA containing duck globin mRNA sequences

A displacement synthesis procedure was used to construct symmetrical recombinant cDNAs. A double-stranded cDNA containing a hairpin loop was extended by the addition of a homopolymer to the 3′ end. This was followed by displacement, or “third strand” synthesis that was primed by an oligonucleotide hybridized to the homopolymer. Ideally, the product should be a DNA containing an inverted repeat with twofold rotational symmetry about nonsymmetrical sequences representing the hairpin loop in the original double-stranded cDNA. Duck globin cDNAs were synthesized by the displacement mode of construction and cloned in pBR322. An α-globin recombinant, pDGPα-2, was isolated and sequenced. This recombinant was found to have two 3′ half regions (mRNA sequence sense) inverted about a nonsymmetrical 5′ half region and an adjacent oligo(dGṡdC) homopolymer. The sequence arrangement indicates that the cDNA folded back on itself, forming a large loop, to prime synthesis of a second strand. We propose that the internal oligo(dGṡ dC) arose through dynamic shifts in cDNA intrastrand structure during the course of synthesis.