Assembly and self-association of oxytricha telomeric nucleoprotein complexes

Two types of specific telomeric protein-DNA complex are reconstituted upon incubation of purified Oxytricha telomere protein with (T4G4)4, an oligodeoxynucleotide of telomeric sequence. The complexes differ in electrophoretic mobility, in protein-DNA contacts, and in the rate of DNA exchange. The patterns of protein-DNA interaction determined by modification interference suggest a model in which the protein can bind either to the two T4G4 repeats at the 3' end or to two internal repeats; in the latter case, it can make a different set of contacts with the terminal repeat to form the more stable complex. Native telomeric chromatin isolated from Oxytricha contains both types of complexes. The reconstituted monomeric complexes associate to give a high molecular weight form that has an altered chemical footprint. Such interactions may mediate the association of chromosomal telomeres in vivo.     

Monovalent cation-induced structure of telomeric DNA: the G-quartet model

We have investigated the structures formed by oligonucleotides composed of two or four repeats of the telomeric sequences from Oxytricha and Tetrahymena. The Oxytricha four-repeat molecule (d(T4G4)4 = Oxy-4) forms structures with increased electrophoretic mobility in nondenaturing gels containing Na+, K+, or Cs+, but not in gels containing Li+ or no added salt. Formation of the folded structure results in protection of a set of dG's from methylation by dimethyl sulfate. Efficient UV-induced cross-links are observed in Oxy-4 and the related sequence from Tetrahymena (d(T2G4)4 = Tet-4), and join thymidine residues in different repeats. Models proposed to account for these data involve G-quartets, hydrogen-bonded structures formed from four guanosine residues in a square-planar array. We propose that the G-quartet structure must be dealt with in vivo by the telomere replication machinery.     

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.     

Frustrated folding of guanine quadruplexes in telomeric DNA

Human chromosomes terminate in long, single-stranded, DNA overhangs of the repetitive sequence (TTAGGG)n. Sets of four adjacent TTAGGG repeats can fold into guanine quadruplexes (GQ), four-stranded structures that are implicated in telomere maintenance and cell immortalization and are targets in cancer therapy. Isolated GQs have been studied in detail, however much less is known about folding in long repeat sequences. Such chains adopt an enormous number of configurations containing various arrangements of GQs and unfolded gaps, leading to a highly frustrated energy landscape. To better understand this phenomenon, we used mutagenesis, thermal melting, and global analysis to determine stability, kinetic, and cooperativity parameters for GQ folding within chains containing 8–12 TTAGGG repeats. We then used these parameters to simulate the folding of 32-repeat chains, more representative of intact telomeres. We found that a combination of folding frustration and negative cooperativity between adjacent GQs increases TTAGGG unfolding by up to 40-fold, providing an abundance of unfolded gaps that are potential binding sites for telomeric proteins. This effect was most pronounced at the chain termini, which could promote telomere extension by telomerase. We conclude that folding frustration is an important and largely overlooked factor controlling the structure of telomeric DNA.     

Conformational stability and mechanism of folding of ribonuclease T1

Urea and thermal unfolding curves for ribonuclease T1 (RNase T1) were determined by measuring several different physical properties. In all cases, steep, single-step unfolding curves were observed. When these results were analyzed by assuming a two-state folding mechanism, the plots of fraction unfolded protein versus denaturant were coincident. The dependence of the free energy of unfolding, delta G (in kcal/mol), on urea concentration is given by delta G = 5.6 - 1.21 (urea). The parameters characterizing the thermodynamics of unfolding are: midpoint of the thermal unfolding curve, Tm = 48.1 degrees C, enthalpy change at Tm, delta Hm = 97 kcal/mol, and heat capacity change, delta Cp = 1650 cal/mol deg. A single kinetic phase was observed for both the folding and unfolding of RNase T1 in the transition and post-transition regions. However, two slow kinetic phases were observed during folding in the pre-transition region. These two slow phases account for about 90% of the observed amplitude, indicating that a faster kinetic phase is also present. The slow phases probably result from cis-trans isomerization at the 2 proline residues that have a cis configuration in folded RNase T1. These results suggest that RNase T1 folds by a highly cooperative mechanism with no structural intermediates once the proline residues have assumed their correct isomeric configuration. At 25 degrees C, the folded conformation is more stable than the unfolded conformations by 5.6 kcal/mol at pH 7 and by 8.9 kcal/mol at pH 5, which is the pH of maximum stability. At pH 7, the thermodynamic data indicate that the maximum conformational stability of 8.3 kcal/mol will occur at -6 degrees C.     

Multifrequency calorimetry of the folding/unfolding transition of cytochrome c.

The folding-unfolding transition of Fe(III) cytochrome c has been studied with the new technique of multifrequency calorimetry. Multifrequency calorimetry is aimed at measuring directly the dynamics of the energetic events that take place during a thermally induced transition by measuring the frequency dispersion of the heat capacity. This is done by modulating the folding/unfolding equilibrium using a variable frequency, small oscillatory temperature perturbation (approximately 0.05-0.1 degrees C) centered at the equilibrium temperature of the system. Fe(III) cytochrome c at pH 4 undergoes a fully reversible folding/unfolding transition centered at 67.7 degrees C and characterized by an enthalpy change of 81 kcal/mol and heat capacity difference between unfolded and folded states of 0.9 kcal/K*mol. By measuring the temperature dependence of the frequency dispersion of the heat capacity in the frequency range of 0.1-1 Hz it has been possible to examine the time regime of the enthalpic events associated with the transition. The multifrequency calorimetry results indicate that approximately 85% of the excess heat capacity associated with the folding/unfolding transition relaxes with a single relaxation time of 326 +/- 68 ms at the midpoint of the transition region. This is the first time that the time regime in which heat is absorbed and released during protein folding/unfolding has been measured.     

Stability of telomeric G-quadruplexes

In most eukaryotes, telomeric DNA consists of repeats of a short motif that includes consecutive guanines and may hence fold into G-quadruplexes. Budding yeasts have telomeres composed of longer repeats and show variation in the degree of repeat homogeneity. Although telomeric sequences from several organisms have been shown to fold into G-quadruplexes in vitro, surprisingly, no study has been dedicated to the comparison of G-quadruplex folding and stability of known telomeric sequences. Furthermore, to our knowledge, folding of yeast telomeric sequences into intramolecular G-quadruplexes has never been investigated. Using biophysical and biochemical methods, we studied sequences mimicking about four repetitions of telomeric motifs from a variety of organisms, including yeasts, with the aim of comparing the G-quadruplex folding potential of telomeric sequences among eukaryotes. G-quadruplex folding did not appear to be a conserved feature among yeast telomeric sequences. By contrast, all known telomeric sequences from eukaryotes other than yeasts folded into G-quadruplexes. Nevertheless, while G(3)T(1-4)A repeats (found in a variety of organisms) and G(4)T(2,4) repeats (found in ciliates) folded into stable G-quadruplexes, G-quadruplexes formed by repetitions of G(2)T(2)A and G(2)CT(2)A motifs (found in many insects and in nematodes, respectively) appeared to be in equilibrium with non-G-quadruplex structures (likely hairpin-duplexes).     

Equilibrium unfolding mechanism of chicken muscle triose phosphate isomerase

Triose phosphate isomerase (TIM) was prepared and purified from chicken breast muscle. The equilibrium unfolding of TIM by urea was investigated by following the changes of intrinsic fluorescence and circular dichroism spectroscopy, and the equilibrium thermal unfolding by differential scanning calorimetry (DSC). Results show that the unfolding of TIM in urea is highly cooperative and no folding intermediate was detected in the experimental conditions used. The thermodynamic parameters of TIM during its urea induced unfolding were calculated as DeltaG degrees =3.54 kcal.mol(-1), and m(G) = 0.67 kcal.mol(-1)M(-1), which just reflect the unfolding of dissociated folded monomer to fully unfolded monomer transition, while the dissociation energy of folded dimer to folded monomer is probe silence. DSC results indicate that TIM unfolding follows an irreversible two-state step with a slow aggregation process. The cooperative unfolding ratio, DeltaH(cal)/DeltaH(vH), was measured close to 2, indicating that the two subunits of chicken muscle TIM unfold independently. The van't Hoff enthalpy, DeltaH(vH), was estimated as about 200 kcal.mol(-1). These results support the unfolding mechanism with a folded monomer formation before its tertiary structure and secondary structure unfolding.     

pH-induced denaturation of proteins: a single salt bridge contributes 3-5 kcal/mol to the free energy of folding of T4 lysozyme

The energetics of a salt bridge formed between the side chains of aspartic acid 70 (Asp70) and histidine 31 (His31) of T4 lysozyme have been examined by nuclear magnetic resonance techniques. The pKa values of the residues in the native state are perturbed from their values in the unfolded protein such that His31 has a pKa value of 9.1 in the native state and 6.8 in the unfolded state at 10 degrees C in moderate salt. Similarly, the aspartate pKa is shifted to a value of about 0.5 in the native state from its value of 3.5-4.0 in the unfolded state. These shifts in pKa show that the salt bridge is stabilized 3-5 kcal/mol. This implies that the salt bridge stabilizes the native state by 3-5 kcal/mol as compared to the unfolded state. This is reflected in the thermodynamic stability of mutants of the protein in which Asp70, His31, or both are replaced by asparagine. These observations and consideration of the thermodynamic coupling of protonation state to folding of proteins suggest a mechanism of acid denaturation in which the unfolded state is progressively stabilized by protonation of its acid residues as pH is lowered below pH 4. The unfolded state is stabilized only if acidic groups in the folded state have lower pKa values than in the unfolded state. When the pH is sufficiently low, the acid groups of both the native and unfolded states are fully protonated, and the apparent unfolding equilibrium constant becomes pH independent. Similar arguments apply to base-induced unfolding.(ABSTRACT TRUNCATED AT 250 WORDS)     

Thermodynamic analysis of the structural stability of the shiga toxin B-subunit

The conformational stability of Shiga toxin B-subunit (STxB), a pentameric protein from Shigella dysenteriae has been characterized by high sensitivity differential scanning calorimetry and circular dichroism spectroscopy under different solvent conditions. It is shown that the thermal folding/unfolding of STxB is a reversible process involving a highly cooperative transition between folded pentamer and unfolded monomers. The conformational stability of STxB is pH dependent and because of its pentameric nature is also concentration dependent. STxB is maximally stable in the pH range from 5 to 9 (Delta G upon unfolding is close to 13 kcal per mol of monomer at 25 degrees C), and its stability decreases both at lower and at higher pH values. The pH dependence of the Gibbs energy of stabilization between pH 2.5 and 5 is consistent with the change in the ionizable state of an average of four groups per monomer upon unfolding. Structural thermodynamic calculations show that the stabilization of the STxB pentamer is primarily due to the interactions established between monomers rather than intramonomer interactions. The folding of an isolated monomer into the conformation existing in the pentamer is unfavorable and expected to be characterized by a free-energy change upon folding in the order of 2.5 kcal mol(-1) at 25 degrees C. On the average, intersubunit interaction induced upon oligomerization of folded monomers should contribute close to -13.4 kcal per mol of monomer to bring the overall Gibbs energy to the experimentally determined value at this temperature.     

Correlation between mutational destabilization of phage T4 lysozyme and increased unfolding rates

The thermodynamics and kinetics of unfolding of 28 bacteriophage T4 lysozyme variants were compared by using urea gradient gel electrophoresis. The mutations studied cause a variety of sequence changes at different residues throughout the polypeptide chain and result in a wide range of thermodynamic stabilities. A striking relationship was observed between the thermodynamic and kinetic effects of the amino acid replacements: All the substitutions that destabilized the native protein by 2 kcal/mol or more also increased the rate of unfolding. The observed increases in unfolding rate corresponded to a decrease in the activation energy of unfolding (delta Gu) at least 35% as large as the decrease in thermodynamic stability (delta Gu). Thus, the destabilizing lesions bring the free energy of the native state closer to that of both the unfolded state and the transition state for folding and unfolding. Since a large fraction of the mutational destabilization is expressed between the transition state and the native conformation, the changes in folding energetics cannot be accounted for by effects on the unfolded state alone. The results also suggest that interactions throughout much of the folded structure are altered in the formation of the transition state during unfolding.     

The yeast telomere length regulator TEL2 encodes a protein that binds to telomeric DNA.

TEL2 is required for telomere length regulation and viability in Saccharomyces cerevisiae. To investigate the mechanism by which Tel2p regulates telomere length, the majority (65%) of the TEL2 ORF was fused to the 3'-end of the gene for maltose binding protein, expressed in bacteria and the purified protein used in DNA binding studies. Rap1p, the major yeast telomere binding protein, recognizes a 13 bp duplex site 5'-GGTGTGTGGGTGT-3' in yeast telomeric DNA with high affinity. Gel shift experiments revealed that the MBP-Tel2p fusion binds the double-stranded yeast telomeric Rap1p site in a sequence-specific manner. Analysis of mutated sites showed that MBP-Tel2p could bind 5'-GTGTGTGG-3' within this 13 bp site. Methylation interference analysis revealed that Tel2p contacts the 5'-terminal guanine in the major groove. MBP-Tel2p did not bind duplex telomeric DNA repeats from vertebrates, Tetrahymena or Oxytricha. These results suggest that Tel2p is a DNA binding protein that recognizes yeast telomeric DNA.     

A plant gene encoding a Myb-like protein that binds telomeric GGTTTAG repeats in vitro

A gene (AtTRP1) encoding a telomeric repeat-binding protein has been isolated from Arabidopsis thaliana. AtTRP1 is a single copy gene located on chromosome 5 of A. thaliana. The protein AtTRP1 encoded by this gene is not only homologous to the Myb DNA-binding motifs of other telomere-binding proteins but also is similar to several initiator-binding proteins in plants. Gel retardation assay revealed that the 115 residues on the C terminus of this protein, including the Myb motif, are sufficient for binding to the double-stranded plant telomeric sequence. The isolated DNA-binding domain of AtTRP1 recognizes each telomeric repeat centered on the sequence GGTTTAG. The almost full-length protein of AtTRP1 does not form any complex at all with the DNA fragments carrying four or fewer GGTTTAG repeats. However, it forms a complex with the sequence (GGTTTAG)(8) more efficiently than with the sequence (GGTTTAG)(5). These data suggest that the minimum length of a telomeric DNA for AtTRP1 binding consists of five GGTTTAG repeats and that the optimal AtTRP1 binding may require eight or more GGTTTAG repeats. It also implies that this protein AtTRP1 may bind in vivo primarily to the ends of plant chromosomes, which consist of long stretches of telomeric repeats.     

PLP-1 binds nematode double-stranded telomeric DNA

The integrity and proper functioning of telomeres require association of telomeric DNA sequences with specific binding proteins. We have characterized PLP-1, a PURa homolog encoded by F45E4.2, which we previously identified as a candidate double stranded telomere binding protein, by affinity chromatography followed by mass spectrometry. PLP-1 bound double-stranded telomeric DNA in vitro as shown by competition assays. Core binding was provided by the third and fourth nucleotides of the TTAGGC telomeric repeat. This is quite different from the binding sequence of CEH-37, another C. elegans telomere binding protein, suggesting that multiple proteins may bind nematode telomeric DNA simultaneously in vivo.     

Force-induced unfolding of human telomeric G-quadruplex: A steered molecular dynamics simulation study

We study the unfolding of a parallel G-quadruplex from human telomeric DNA by mechanical stretching using steered molecular dynamics (MD) simulation. We find that the force curves and unfolding processes strongly depend on the pulling sites. With pulling sites located on the sugar-phosphate backbone, the force-extension curve shows a single peak and the unfolding proceeds sequentially. Pulling sites located on the terminal nucleobases lead to a force-extension curve with two peaks and the unfolding is more cooperative. Simulations of the refolding of partially unfolded quadruplexes show very different behavior for the two different pulling modalities. In particular, starting from an unfolded state prepared by nucleobase pulling leads to a long-lived intermediate state whose existence is also corroborated by the free energy profile computed with the Jarzynski equation. Based on this observation, we propose a novel folding pathway for parallel G-quadruplexes with the human telomere sequence.     

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10 degrees C by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG=2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding rate of less than 4 min(-1) measured by stopped-flow circular dichroism. The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E.coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding.     

G148-GA3: a streptococcal virulence module with atypical thermodynamics of folding optimally binds human serum albumin at physiological temperatures

The third albumin binding domain of streptococcal protein G strain 148 (G148-GA3) belongs to a novel class of prokaryotic albumin binding modules that is thought to support virulence in several bacterial species. Here, we characterize G148-GA3 folding and albumin binding by using differential scanning calorimetry and isothermal titration calorimetry to obtain the most complete set of thermodynamic state functions for any member of this medically significant module. When buffered at pH 7.0 the 46-amino acid alpha-helical domain melts at 72 degrees C and exhibits marginal stability (15 kJ/mol) at 37 degrees C. G148-GA3 unfolding is characterized by small contributions to entropy from non-hydrophobic forces and a low DeltaCp (1.1 kJ/(deg mol)). Isothermal titration calorimetry reveals that the domain has evolved to optimally bind human serum albumin near 37 degrees C with a binding constant of 1.4 x 10 7 M(-1). Analysis of G148-GA3 thermodynamics suggests that the domain experiences atypically small per residue changes in structural dynamics and heat capacity while transiting between folded and unfolded states.     

An aggregation-prone intermediate species is present in the unfolding pathway of the monomeric portal protein of bacteriophage P22: implications for portal assembly

The head of the P22 bacteriophage is interrupted by a unique dodecameric portal vertex that serves as a conduit for the entrance and exit of the DNA. Here, the in vitro unfolding/refolding processes of the portal protein of P22 were investigated at different temperatures (1, 25, and 37 degrees C) through the use of urea and high hydrostatic pressure (HHP) combined with spectroscopic techniques. We have characterized an intermediate species, IU, which forms at 25 degrees C during unfolding or refolding of the portal protein in 2-4 M urea. IU readily forms amorphous aggregates, rendering the folding process irreversible. On the other hand, at 1 degrees C, a two-state process is observed (DeltaGf = -2.2 kcal/mol). When subjected to HHP at 25 or 37 degrees C, the portal monomer undergoes partial denaturation, also forming an intermediate species, which we call IP. IP also tends to aggregate but, differently from IU, aggregates into a ring-like structure as seen by size-exclusion chromatography and electron microscopy. Again, at 1 degrees C the unfolding induced by HHP proved to be reversible, with DeltaGf = -2.4 kcal/mol and DeltaV = 72 mL/mol. Interestingly, at 25 degrees C, the binding of the hydrophobic probe bis-ANS to the native portal protein destabilizes it and completely blocks its aggregation under HHP. These data are relevant to the process by which the portal protein assembles into dodecamers in vivo, since species such as IP must prevail over IU in order to guarantee the proper ring formation.     

Mutational analysis of the BPTI folding pathway: II. Effects of aromatic-->leucine substitutions on folding kinetics and thermodynamics.

The rates of the individual steps in the disulfide-coupled folding and unfolding of eight BPTI variants, each containing a single aromatic to leucine amino acid replacement, were measured. From this analysis, the contributions of the four phenylalanine and four tyrosine residues to the stabilities of the native protein and the disulfide-bonded folding intermediates were determined. While the substitutions were found to destabilize the native protein by 2 to 7 kcal/mol, they had significantly smaller effects on the intermediates that represent the earlier stages of folding, even when the site of the substitution was located within the ordered regions of the intermediates. These results suggest that stabilizing interactions contribute less to conformational stability in the context of a partially folded intermediate than in a fully folded native protein, perhaps because of decreased cooperativity among the individual interactions. The kinetic analysis also provides new information about the transition states associated with the slowest steps in folding and unfolding, supporting previous suggestions that these transition states are extensively unfolded. Although the substitutions caused large changes in the distribution of folding intermediates and in the rates of some steps in the folding pathway, the kinetically-preferred pathway for all of the variants involved intramolecular disulfide rearrangements, as observed previously for the wild-type protein. These results suggest that the predominance of the rearrangement mechanism reflects conformational constraints present relatively early in the folding pathway.     

Sequence-specific DNA recognition by the Myb-like domain of plant telomeric protein RTBP1

We have identified a rice gene encoding a DNA-binding protein that specifically recognizes the telomeric repeat sequence TTTAGGG found in plants. This gene, which we refer to as RTBP1 (rice telomere-binding protein 1), encodes a polypeptide with a predicted molecular mass of 70 kDa. RTBP1 is ubiquitously expressed in various organs and binds DNA with two or more duplex TTTAGGG repeats. The predicted protein sequence includes a single domain at the C terminus with extensive homology to Myb-like DNA binding motif. The Myb-like domain of RTBP1 is very closely related to that of other telomere-binding proteins, including TRF1, TRF2, Taz1p, and Tbf1p, indicating that DNA-binding domains of telomere-binding proteins are well conserved among evolutionarily distant species. To obtain precise information on the sequence of the DNA binding site recognized by RTBP1, we analyzed the sequence-specific binding properties of the isolated Myb-like domain of RTBP1. The isolated Myb-like domain was capable of sequence-specific DNA binding as a homodimer. Gel retardation analysis with a series of mutated telomere probes revealed that the internal GGGTTT sequence in the two-telomere repeats is critical for binding of Myb-like domain of RTBP1, which is consistent with the model of the TRF1.DNA complex showing that base-specific contacts are made within the sequence GGGTTA. To the best of our knowledge, RTBP1 is the first cloned gene in which the product is able to bind double-stranded telomeric DNA in plants. Because the Myb-like domain appears to be a significant motif for a large class of proteins that bind the duplex telomeric DNA, RTBP1 may play important roles in plant telomere function in vivo.