Molecular dynamics study of water channels in natural and synthetic amyloid-β fibrils

We compared all-atom explicit solvent molecular dynamics simulations of three types of Aβ(1-40) fibrils: brain-seeded fibrils (2M4J, with a threefold axial symmetry) and the other two, all-synthetic fibril polymorphs (2LMN and 2LMP, made under different fibrillization conditions). Fibril models were constructed using either a finite or an infinite number of layers made using periodic images. These studies yielded four conclusions. First, finite fibrils tend to unravel in a manner reminiscent of fibril dissolution, while infinite fibrils were more stable during simulations. Second, salt bridges in these fibrils remained stable in those fibrils that contained them initially, and those without salt bridges did not develop them over the time course of the simulations. Third, all fibrils tended to develop a "stagger" or register shift of β-strands along the fibril axis. Fourth and most importantly, the brain-seeded, 2M4J, infinite fibrils allowed bidirectional transport of water in and out of the central longitudinal core of the fibril by rapidly developing gaps at the fibril vertices. 2LMP fibrils also showed this behavior, although to a lesser extent. The diffusion of water molecules in the fibril core region involved two dynamical states: a localized state and directed diffusion in the presence of obstacles. These observations provided support for the hypothesis that Aβ fibrils could act as nanotubes. At least some Aβ oligomers resembled fibrils structurally in having parallel, in-register β-sheets and a sheet-turn-sheet motif. Thus, our findings could have implications for Aβ cytotoxicity, which may occur through the ability of oligomers to form abnormal water and ion channels in cell membranes.

Ionic strength independence of charge distributions in solvation of biomolecules

Electrostatic forces enormously impact the structure, interactions, and function of biomolecules. We perform all-atom molecular dynamics simulations for 5 proteins and 5 RNAs to determine the dependence on ionic strength of the ion and water charge distributions surrounding the biomolecules, as well as the contributions of ions to the electrostatic free energy of interaction between the biomolecule and the surrounding salt solution (for a total of 40 different biomolecule/solvent combinations). Although water provides the dominant contribution to the charge density distribution and to the electrostatic potential even in 1M NaCl solutions, the contributions of water molecules and of ions to the total electrostatic interaction free energy with the solvated biomolecule are comparable. The electrostatic biomolecule/solvent interaction energies and the total charge distribution exhibit a remarkable insensitivity to salt concentrations over a huge range of salt concentrations (20 mM to 1M NaCl). The electrostatic potentials near the biomolecule's surface obtained from the MD simulations differ markedly, as expected, from the potentials predicted by continuum dielectric models, even though the total electrostatic interaction free energies are within 11% of each other.

Characterization of tertiary folding of RNA by circular dichroism and urea

CD spectroscopy can be used to monitor RNA tertiary folding transitions that may not be observable by absorbance spectroscopy. With the use of computer-controlled titrators, data can be acquired rapidly, and accurate thermodynamic properties can be obtained over a wide variety of conditions. Thus, CD spectroscopy provides a useful complement to site-resolved or chemical modification methods.

Reduced contact order and RNA folding rates

We investigated the relationship between RNA structure and folding rates accounting for hierarchical structural formation. Folding rates of two-state folding proteins correlate well with relative contact order, a quantitative measure of the number and sequence distance between tertiary contacts. These proteins do not form stable structures prior to the rate-limiting step. In contrast, most secondary structures are stably formed prior to the rate-limiting step in RNA folding. Accordingly, we introduce "reduced contact order", a metric that reflects only the number of residues available to participate in the conformational search after the formation of secondary structure. Plotting the folding rates and the reduced contact order from ten different RNAs suggests that RNA folding can be divided into two classes. To examine this division, folding rates of circularly permutated isomers are compared for two RNAs, one from each class. Folding rates vary by tenfold for circularly permuted Bacillus subtilis RNase P RNA isomers, whereas folding rates vary by only 1.2-fold for circularly permuted catalytic domains. This difference is likely related to the dissimilar natures of their rate-limiting steps.

Stepwise conversion of a mesophilic to a thermophilic ribozyme

A fundamental question in RNA folding is the mechanism of thermodynamic stability. We investigated the equilibrium folding of a series of sequence variants in which one to three motifs of a 255-nucleotide mesophilic ribozyme were substituted with the corresponding motifs from its thermophilic homologue. Substitution of three crucial motifs individually or in groups results in a continual increase in the stability and folding cooperativity in a stepwise fashion. We find an unexpected relationship between stability and folding cooperativity. Without changing the folding cooperativity, RNAs having a similar native structure can only achieve moderate change in stability and likewise, without changing stability, RNAs having a similar native structure can only achieve moderate change in folding cooperativity. This intricate relationship must be included in the predictions of tertiary RNA stability.

The rate-limiting step in the folding of a large ribozyme without kinetic traps

A fundamental question in RNA folding is the nature of the rate-limiting step. Folding of large RNAs often is trapped by the need to undo misfolded structures, which precludes the study of the other, potentially more interesting aspects in the rate-limiting step, such as conformational search, metal ion binding, and the role of productive intermediates. The catalytic domain of the Bacillus subtilis RNase P RNA folds without a kinetic trap, thereby providing an ideal system to elucidate these steps. We analyzed the folding kinetics by using fluorescence and absorbance spectroscopies, catalytic activity, and synchrotron small-angle x-ray scattering. Folding begins with the rapid formation of early intermediates wherein the majority of conformational search occurs, followed by the slower formation of subsequent intermediates. Before the rate-limiting step, more than 98% of the total structure has formed. The rate-limiting step is a small-scale structural rearrangement involving prebound metal ions.

Engineered metal binding sites map the heterogeneous folding landscape of a coiled coil

To address whether proteins fold along multiple pathways, i,i+4 bi-histidine metal binding sites are introduced into dimeric and crosslinked versions of the leucine zipper region of the growth control transcription factor, GCN4. Divalent metal ion binding enhances both the equilibrium and folding activation free energies for GCN4. The enhancement of folding rates quantifies the fraction of molecules that have the binding site in a helical geometry in the transition state. Hence, this new method, termed Psi-analysis, identifies the degree of pathway heterogeneity for a protein that folds in a two-state manner, a capability that is generally unavailable even with single molecule methods. Adjusting metal ion concentration continuously varies the stability of the bi-histidine region without additional structural perturbation to the protein. For dimeric and crosslinked versions, the accompanying changes in kinetic barrier heights at each metal ion concentration maps the folding landscape as well as establishes the importance of connectivity in pathway selection. Furthermore, this method can be generalized to other biophysical studies, where the ability to continuously tune the stability of a particular region with no extraneous structural perturbation is advantageous.

Modular construction of a tertiary RNA structure: the specificity domain of the Bacillus subtilis RNase P RNA

The structure of the specificity domain (S-domain) of the Bacillus subtilis RNase P RNA has been proposed to be composed of a core and a buttress module, analogous to the bipartite structure of the P4-P6 domain of the Tetrahymena group I ribozyme. The core module is the functional unit of the S-domain and contains the binding site for the T stem-loop of a tRNA. The buttress module provides structural stability to the core module and consists of a GA3 tetraloop and its receptor. To explicitly test the hypothesis that modular construction can describe the structure of the S-domain and is a useful RNA design strategy, we analyzed the equilibrium folding and substrate binding of three classes of S-domain mutants. Addition or deletion of a base pair in the helical linker region between the modules only modestly destabilizes the tertiary structure. tRNA binding selectivity is affected in one but not in two other mutants of this class. Elimination of the GA3 tetraloop-receptor interactions significantly destabilizes the core module and results in the loss of tRNA binding selectivity. Replacing the buttress module with that of a homologous RNase P RNA maintains the tRNA binding selectivity. Overall, we have observed that the linker regions between the two modules can tolerate moderate structural changes and that the buttress modules can be shuffled between homologous S-domains. These results suggest that it is feasible to design an RNA using a buttress module to stabilize a functional module.

PAS domain receptor photoactive yellow protein is converted to a molten globule state upon activation

Biological signaling generally involves the activation of a receptor protein by an external stimulus followed by protein-protein interactions between the activated receptor and its downstream signal transducer. The current paradigm for the relay of signals along a signal transduction chain is that it occurs by highly specific interactions between fully folded proteins. However, recent results indicate that many regulatory proteins are intrinsically unstructured, providing a serious challenge to this paradigm and to the nature of structure-function relationships in signaling. Here we study the structural changes that occur upon activation of the blue light receptor photoactive yellow protein (PYP). Activation greatly reduces the tertiary structure of PYP but leaves the level secondary structure largely unperturbed. In addition, activated PYP exposes previously buried hydrophobic patches and allows significant solvent penetration into the core of the protein. These traits are the distinguishing hallmarks of molten globule states, which have been intensively studied for their role in protein folding. Our results show that receptor activation by light converts PYP to a molten globule and indicate stimulus-induced unfolding to a partially unstructured molten globule as a novel theme in signaling.

The thermodynamic origin of the stability of a thermophilic ribozyme

Understanding the mechanism of thermodynamic stability of an RNA structure has significant implications for the function and design of RNA. We investigated the equilibrium folding of a thermophilic ribozyme and its mesophilic homologue by using hydroxyl radical protection, small-angle x-ray scattering, and circular dichroism. Both RNAs require Mg(2+) to fold to their native structures that are very similar. The stability is measured as a function of Mg(2+) and urea concentrations at different temperatures. The enhanced stability of the thermophilic ribozyme primarily is derived from a tremendous increase in the amount of structure formed in the ultimate folding transition. This increase in structure formation and cooperativity arises because the penultimate and the ultimate folding transitions in the mesophilic ribozyme become linked into a single transition in the folding of the thermophilic ribozyme. Therefore, the starting point, or reference state, for the transition to the native, functional thermophilic ribozyme is significantly less structured. The shift in the reference state, and the resulting increase in folding cooperativity, is likely due to the stabilization of selected native interactions that only form in the ultimate transition. This mechanism of using a less structured intermediate and increased cooperativity to achieve higher functional stability for tertiary RNAs is fundamentally different from that commonly proposed to explain the increased stability of thermophilic proteins.

Altering the intermediate in the equilibrium folding of unmodified yeast tRNAPhe with monovalent and divalent cations

The isothermal equilibrium folding of the unmodified yeast tRNA(Phe) is studied as a function of Na(+), Mg(2+), and urea concentration with hydroxyl radical protection, circular dichroism, and diethyl pyrocarbonate (DEPC) modification. These assays indicate that this tRNA folds in Na(+) alone. Similar to folding in Mg(2+), folding in Na(+) can be described by two transitions, unfolded-to-intermediate-to-native. The I-to-N transition has a Na(+) midpoint of approximately 0.5 M and a Hill constant of approximately 4. Unexpectedly, the urea m-value, the dependence of free energy on urea concentration, for the I-to-N transition is significantly smaller in Na(+) than in Mg(2+), 0.4 versus 1.7 kcal mol(-1) M(-1), indicating that more structure is formed in the Mg(2+)-induced transition. DEPC modification indicates that the I state in Na(+)-induced folding contains all four helices of tRNA and the I-to-N transition primarily corresponds to the formation of the tertiary structure. In contrast, the intermediate in Mg(2+)-induced folding contains only three helices, and the I-to-N transition corresponds to the formation of the acceptor stem plus tertiary structure. The cation dependence of the intermediates arises from the differences in the stability of the acceptor stem and the tertiary structure. The acceptor stem is stable at a lower Na(+) concentration than required for the tertiary structure formation. The relative stability is reversed in Mg(2+) so that the acceptor stem and the tertiary structure form simultaneously in the I-to-N transition. These results demonstrate that formation of the RNA secondary structure can be independent or coupled to the formation of the tertiary structure depending on their relative stability in monovalent and divalent ions.

The Bacillus subtilis RNase P holoenzyme contains two RNase P RNA and two RNase P protein subunits

Ribonuclease P (RNase P) catalyzes the 5' maturation of precursor tRNA transcripts and, in bacteria, is composed of a catalytic RNA and a protein. We investigated the oligomerization state and the shape of the RNA alone and the holoenzyme of Bacillus subtilis RNase P in the absence of substrate by synchrotron small-angle X-ray scattering and affinity retention. The B. subtilis RNase P RNA alone is a monomer; however, the scattering profile changes upon the addition of monovalent ions, possibly suggesting different interdomain angles. To our surprise, the X-ray scattering data combined with the affinity retention results indicate that the holoenzyme contains two RNase P RNA and two RNase P protein molecules. We propose a structural model of the holoenzyme with a symmetrical arrangement of the two RNA subunits, consistent with the X-ray scattering results. This (P RNA)2(P protein)2 complex likely binds substrate differently than the conventional (P RNA)1(P protein)1 complex; therefore, the function of the B. subtilis RNase P holoenzyme may be more diverse than previously thought. These revisions to our knowledge of the RNase P holoenzyme suggest a more versatile role for proteins in ribonucleoprotein complexes.

Distinguishing between two-state and three-state models for ubiquitin folding

Conflicting results exist regarding whether the folding of mammalian ubiquitin at 25 degrees C is a simple, two-state kinetic process or a more complex, three-state process with a defined kinetic intermediate. We have measured folding rate constants up to about 1000 s(-1) using conventional rapid mixing methods in single-jump, double-jump, and continuous-flow modes. The linear dependence of folding rates on denaturant concentration and the lack of an unaccounted "burst-phase" change for the fluorescence signal indicate that a two-state folding model is adequate to describe the folding pathway. This behavior also is seen for folding in the presence of the stabilizing additives 0.23 M sodium sulfate and 1 M sodium chloride. These results stress the need for caution in interpreting deviations from ideal two-state "chevron" behavior when folding is heterogeneous or folding rate constants are near the detection limit.

Mg2+-dependent compaction and folding of yeast tRNAPhe and the catalytic domain of the B. subtilis RNase P RNA determined by small-angle X-ray scattering

We apply synchrotron-based small-angle X-ray scattering to investigate the relationship between compaction, metal binding, and structure formation of two RNAs at 37 degrees C: the 76 nucleotide yeast tRNA(Phe) and the 255 nucleotide catalytic domain of the Bacillus subtilis RNase P RNA. For both RNAs, this method provides direct evidence for the population of a distinct folding intermediate. The relative compaction between the intermediate and the native state does not correlate with the size of the RNA but does correlate well with the amount of surface burial as quantified previously by the urea-dependent m-value. The total compaction process can be described in two major stages. Starting from a completely unfolded state (4-8 M urea, no Mg(2+)), the major amount of compaction occurs upon the dilution of the denaturant and the addition of micromolar amounts of Mg(2+) to form the intermediate. The native state forms in a single transition from the intermediate state upon cooperative binding of three to four Mg(2+) ions. The characterization of this intermediate by small-angle X-ray scattering lends strong support for the cooperative Mg(2+)-binding model to describe the stability of a tertiary RNA.

D/H amide kinetic isotope effects reveal when hydrogen bonds form during protein folding

We have exploited a procedure to identify when hydrogen bonds (H-bonds) form under two-state folding conditions using equilibrium and kinetic deuterium/hydrogen amide isotope effects. Deuteration decreases the stability of equine cytochrome c and the dimeric and crosslinked versions of the GCN4-p1 coiled coil by approximately 0. 5 kcal mol-1. For all three systems, the decrease in equilibrium stability is reflected by a decrease in refolding rates and a near equivalent increase in unfolding rates. This apportionment indicates that approximately 50% of the native H-bonds are formed in the transition state of these helical proteins. In contrast, an alpha/beta protein, mammalian ubiquitin, exhibits a small isotope effect only on unfolding rates, suggesting its folding pathway may be different. These four proteins recapitulate the general trend that approximately 50% of the surface buried in the native state is buried in the transition state, leading to the hypothesis that H-bond formation in the transition state is cooperative, with alpha-helical proteins forming a number of H-bonds proportional to the amount of surface buried in the transition state.

Applicability of urea in the thermodynamic analysis of secondary and tertiary RNA folding

The equilibrium folding of a series of self-complementary RNA duplexes and the unmodified yeast tRNA(Phe) is studied as a function of urea and Mg(2+) concentration with optical spectroscopies and chemical modification under isothermal conditions. Via application of standard methodologies from protein folding, the folding free energy and its dependence on urea concentration, the m value, are determined. The free energies of the RNA duplexes obtained from the urea titrations are in good agreement with those calculated from thermal melting studies [Freier, S. I., et al. (1986) Proc. Natl. Acad. Sci. U.S.A. 83, 9373]. The m value correlates with the length of the RNA duplex and is not sensitive to ionic conditions and temperature. The folding of the unmodified yeast tRNA(Phe) can be described by two Mg(2+)-dependent transitions, the second of which corresponds to the formation of the native tertiary structure as confirmed by hydroxyl radical protection and partial nuclease digestion. Both transitions are sensitive to urea and have m values of 0.94 and 1.70 kcal mol(-)(1) M(-)(1), respectively. Although the precise chemical basis of urea denaturation of RNA is uncertain, the m values for the duplexes and tRNA(Phe) are proportional to the amount of the surface area buried in the folding transition. This proportionality, 0.099 cal mol(-)(1) M(-)(1) A(-)(2), is very similar to that observed for proteins, 0.11 cal mol(-)(1) M(-)(1) A(-)(2) [Myers, J., Pace, N., and Scholtz, M. (1995) Protein Sci. 4, 2138]. These results indicate that urea titration can be used to measure both the free energy and the magnitude of an RNA folding transition.

A thermodynamic framework and cooperativity in the tertiary folding of a Mg2+-dependent ribozyme

The folding thermodynamics of the catalytic domain from the Bacillus subtilis RNase P RNA is analyzed using circular dichroism and fluorescence spectroscopies, hydroxyl radical protection, and catalytic activity. Folding of this 255-nucleotide ribozyme can be described with three populated species: unfolded (U), intermediate (I), and native (N) states. The U-to-I transition primarily involves secondary structure formation, whereas the I-to-N transition is dominated by tertiary structure formation. The I-to-N transition is highly cooperative as indicated by the coincidence of the four probes applied here. Two isothermal methods are used to determine the stability of the N state relative to the I state at 10 and 37 degrees C. The first method measures the extent of Mg(2+)-induced folding without urea or at constant urea concentrations. The second method measures the extent of urea-induced unfolding at constant Mg(2+) concentrations. Via application of a cooperative binding analysis, the Mg(2+) transition midpoint (K(Mg)), the Hill constant (n), and the urea-dependent surface burial parameter (m value) determined by both methods are identical, indicating that they report the same, reversible folding event. Three conclusions can be drawn from these results. (i) The folding free energy of a Mg(2+)-dependent tertiary RNA structure can be described by the K(Mg) and n parameters according to a cooperative Mg(2+) binding model. (ii) The Hill constant for this tertiary RNA structure probably represents the differential number of Mg(2+) ions bound in the I-to-N transition. (iii) Under physiological conditions, the stability of this large ribozyme is similar to that of small globular proteins.

Mg2+-dependent folding of a large ribozyme without kinetic traps

The folding kinetics of the catalytic domain of Bacillus subtilis ribonuclease P is analyzed here by fluorescence and catalytic activity. The folding pathway is apparently free of kinetic traps, as indicated by a decrease in folding rates upon the addition of urea. We apply Mg2+ and urea chevron analysis to fully describe the folding and unfolding kinetics of this ribozyme. A folding scheme containing two kinetic intermediates completely accounts for the free energy, the Mg2+ Hill coefficient and the surface buried in the equilibrium transition. At saturating Mg 2+concentrations, folding is limited by a barrier that is independent of Mg2+ and urea. These results describe the first trap-free folding pathway of a large ribozyme and indicate that kinetic traps are not an obligate feature of RNA folding.

Transition state heterogeneity in GCN4 coiled coil folding studied by using multisite mutations and crosslinking

We have investigated the folding behavior of dimeric and covalently crosslinked versions of the 33-residue alpha-helical GCN4-p1 coiled coil derived from the leucine zipper region of the transcriptional activator GCN4. The effects of multisite substitutions indicate that folding occurs along multiple routes with nucleation sites located throughout the protein. The similarity in activation energies of the different routes together with an analysis of intrinsic helical propensities indicate that minimal helix is present before a productive collision of the two chains. However, approximately one-third to one-half of the total helical structure is formed in the postcollision transition state ensemble. For the crosslinked, monomeric version, folding occurs along a single robust pathway. Here, the region nearest the crosslink, with the least helical propensity, is structured in the transition state whereas the region farthest from the tether, with the most propensity, is completely unstructured. Hence, the existence of transition state heterogeneity and the selection of folding routes critically depend on chain topology.

Folding of a large ribozyme during transcription and the effect of the elongation factor NusA

We compared in vitro transcription-initiated folding of the ribozyme from Bacillus subtilis RNase P to refolding from the full-length, denatured state by monitoring the appearance of its catalytic activity. At 37 degrees C, Mg(2+)-initiated refolding of the wild type and a circularly permutate ribozyme takes minutes and is limited by a kinetic trap. Transcription by T7 RNA polymerase alters the folding pathway of both RNAs and introduces new kinetic traps. Transcription by the core Escherichia coli RNA polymerase yields the same result, in spite of its 4-fold-slower elongation rate. However, the presence of its elongation factor NusA accelerates more than 10-fold the transcription-initiated folding of the circularly, permutated ribozyme by E. coli RNA polymerase. The effect of NusA likely is caused by its enhancement of transcriptional pausing because NusA did not accelerate transcription-initiated folding using a mutant RNA polymerase that failed to pause or respond to NusA during ribozyme synthesis. We conclude that both transcription and specific pausing therein can alter RNA-folding pathways.

Viscosity dependence of the folding kinetics of a dimeric and monomeric coiled coil

We measured whether solvent viscosity, and hence chain diffusion, plays a role in the rate-limiting step of the folding reactions of GCN4-p2', a simple alpha-helical coiled coil derived from the leucine zipper region of bZIP transcriptional activator GCN4. To deconvolute the dual effects of viscosogenic solvents on both viscosity, eta, and stability, earlier attempts assumed that the cosolvent and denaturant interact to the same degree in the transition state. Applying this analysis to GCN4-p2' yielded a nearly 1/eta dependence between folding rates and viscosity for both the dimeric and the cross-linked, monomeric versions of the coiled coil, but it revealed no such coherent relationship for cytochrome c. We also developed a method to determine the relative viscosity dependence of the dimeric and monomeric forms of the coiled coil independent of the assumption concerning the transition state's relative interaction with cosolvents and denaturants. Application of this method indicated that the effect of viscosity on both the folding and the unfolding rates was the same for the dimeric and monomeric versions, further supporting the view that the folding of the dimeric version is folding-limited rather than encounter-limited. The finding that GCN4-p2' folding appears to exhibit a 1/eta viscosity dependence implies that the rate-limiting step in folding is opposed predominantly by solvent-derived rather than internal frictional forces. These results are interpreted in relation to various models for protein folding.

Trifluoroethanol promotes helix formation by destabilizing backbone exposure: desolvation rather than native hydrogen bonding defines the kinetic pathway of dimeric coiled coil folding

We measure the effects of low concentrations of helix-stabilizing cosolvents, including 2,2,2-trifluoroethanol (TFE), on the thermodynamics and kinetics of folding of the dimeric alpha-helical coiled coil derived from the leucine zipper region of bZIP transcriptional activator GCN4. The change in kinetic behavior upon addition of 5% (v/v) TFE indicates that it stabilizes the transition state to the same degree as the fully helical native state. However, folding rates are largely insensitive to alanine to glycine mutagenesis, indicating that the majority of helical structure is formed after the transition state. Equilibrium hydrogen isotope partitioning measurements indicate that intramolecular hydrogen bonds are not strengthened by TFE and that amide hydrogen bonds in the transition state are nearly the same strength as those in the unfolded state. Thus, the mechanism by which TFE exerts its helix-stabilizing effects can be divorced from helix formation and does not depend on the strengthening of intrahelical hydrogen bonds. Rather, TFE increases the structure of the binary alcohol/water solvent, thereby increasing the energetic cost associated with solvation of the polypeptide backbone. At low concentrations, TFE destabilizes the unfolded species and thereby indirectly enhances the kinetics and thermodynamics of folding of the coiled coil. A high degree of polypeptide backbone desolvation, and not the formation of regular helical structure and native strength hydrogen bonds, is the critical feature of the transition state for folding of this small dimeric protein.

The burst phase in ribonuclease A folding and solvent dependence of the unfolded state

Submillisecond burst phase signals measured in kinetic protein folding experiments have been widely interpreted in terms of the fast formation of productive folding intermediates. Experimental comparisons with non-folding polypeptide chains show that, for ribonuclease A and cytochrome c, these signals in fact reflect a shift from one biased ensemble of the unfolded state to another as a function of change in denaturant concentration.

Intermediates and kinetic traps in the folding of a large ribozyme revealed by circular dichroism and UV absorbance spectroscopies and catalytic activity

The folding thermodynamics and kinetics for the ribozyme from Bacillus subtilis RNase P are analyzed using circular dichroism and UV absorbance spectroscopies and catalytic activity. At 37 degrees C, the addition of Mg2+ (Kd approximately 50 microM) to the unfolded state produces an intermediate state within 1 ms which contains a comparable amount of secondary structure as the native ribozyme. The subsequent transition to the native state (Kd[Mg] approximately 0.8 mM, Hill coefficient approximately 3.5) has a half-life of hundreds of seconds as measured by circular dichroism at 278 nm and by a ribozyme activity assay. Surprisingly, the formation of the native structure is accelerated strongly by the addition of a denaturant; approximately 30-fold at 4.5 M urea. Thus, the rate-limiting step entails the disruption of a considerable number of interactions. The folding of this, and presumably other large RNAs, is slow due to the structural rearrangement of kinetically trapped species. Taken together with previous submillisecond relaxation kinetics of tRNA tertiary structure, we suggest that error-free RNA folding can be on the order of milliseconds.

Ultrafast signals in protein folding and the polypeptide contracted state

To test the significance of ultrafast protein folding signals (<<1 msec), we studied cytochrome c (Cyt c) and two Cyt c fragments with major C-terminal segments deleted. The fragments remain unfolded under all conditions and so could be used to define the unfolded baselines for protein fluorescence and circular dichroism (CD) as a function of denaturant concentration. When diluted from high to low denaturant in kinetic folding experiments, the fragments readjust to their new baseline values in a "burst phase" within the mixing dead time. The fragment burst phase reflects a contraction of the polypeptide from a more extended unfolded condition at high denaturant to a more contracted unfolded condition in the poorer, low denaturant solvent. Holo Cyt c exhibits fluorescence and CD burst phase signals that are essentially identical to the fragment signals over the whole range of final denaturant concentrations, evidently reflecting the same solvent-dependent, relatively nonspecific contraction and not the formation of a specific folding intermediate. The significance of fast folding signals in Cyt c and other proteins is discussed in relation to the hypothesis of an initial rate-limiting search-nucleation-collapse step in protein folding [Sosnick, T. R., Mayne, L. & Englander, S. W. (1996) Proteins Struct. Funct. Genet. 24, 413-426].

Thymosin beta 4 binds actin in an extended conformation and contacts both the barbed and pointed ends

The beta-thymosins are a family of highly polar peptides which serve in vivo to maintain a reservoir of unpolymerized actin monomers. In vitro, beta-thymosins form 1:1 complexes with actin monomers and inhibit both polymerization and exchange of the bound nucleotide. Circular dichroism data indicate that free thymosin beta 4 is predominantly unstructured, containing at most six residues of alpha-helix, and that up to six additional residues may adopt an alpha-helical conformation upon binding actin. NMR data indicate that many parts of thymosin beta 4 are not in tight contact with actin. Contacts between specific residues in actin and thymosin beta 4 were identified by zero-length cross-linking followed by isolation and sequencing of cross-linked peptides. After carbodiimide-mediated cross-linking, Lys-3 of thymosin beta 4 was cross-linked to Glu-167 of actin, and Lys-18 of thymosin beta 4 was cross-linked to one of the the N-terminal acidic residues of actin (Asp-1-Glu-4); the cross-linked actin residues lie within subdomains 3 and 1, respectively. These two contacts flank the alpha-helical region of thymosin beta 4 and place it on the barbed end; thymosin beta 4 can thus block actin polymerization sterically. After transglutaminase-mediated cross-linking, Lys-38 of thymosin beta 4 was cross-linked to Gln-41 of actin, placing the C-terminal region of thymosin beta 4 in contact with subdomain 2 on the pointed end; thymosin beta 4 may sterically block actin polymerization at the pointed end as well as the barbed end of the monomer. The distance between the pointed-end and barbed-end contacts requires that the C-terminal half of thymosin beta 4 be in a predominantly extended conformation.

Hydrogen exchange: the modern legacy of Linderstrøm-Lang

This discussion, prepared for the Protein Society's symposium honoring the 100th anniversary of Kaj Linderstrøm-Lang, shows how hydrogen exchange approaches initially conceived and implemented by Lang and his colleagues some 50 years ago are contributing to current progress in structural biology. Examples are chosen from the active protein folding field. Hydrogen exchange methods now make it possible to define the structure of protein folding intermediates in various contexts: as tenuous molten globule forms at equilibrium under destabilizing conditions, in kinetic intermediates that exist for less than one second, and as infinitesimally populated excited state forms under native conditions. More generally, similar methods now find broad application in studies of protein structure, energetics, and interactions. This article considers the rise of these capabilities from their inception at the Carlsberg Labs to their contemporary role as a significant tool of modern structural biology.

Molecular collapse: the rate-limiting step in two-state cytochrome c folding

Experiments with cytochrome c (cyt c) show that an initial folding event, molecular collapse, is not an energetically downhill continuum as commonly presumed but represents a large-scale, time-consuming, cooperative barrier-crossing process. In the absence of later misfold-reorganization barriers, the early collapse barrier limits cyt c folding to a time scale of milliseconds. The collapse process itself appears to be limited by an uphill search for some coarsely determined transition state structure that can nucleate subsequent energetically downhill folding events. An earlier "burst phase" event at strongly native conditions appears to be a non-specific response of the unfolded chain to reduced denaturant concentration. The molecular collapse process may or may not require the co-formation of the amino- and carboxyl-terminal helices, which are present in an initial metastable intermediate directly following the rate-limiting collapse. After the collapse-nucleation event, folding can proceed rapidly in an apparent two-state manner, probably by way of a predetermined sequence of metastable intermediates that leads to the native protein structure (Bai et al., Science 269:192-197, 1995).

The role of helix formation in the folding of a fully alpha-helical coiled coil

To determine when secondary structure forms as two chains coalesce to form an alpha-helical dimer, the folding rates of variants of the coiled coil region of GCN4 were compared. Residues at non-perturbing positions along the exterior length of the helices were substituted one at a time with alanine and glycine to vary helix propensity and therefore dimer stability. For all variants, the bimolecular folding rate remains largely unchanged; the unfolding rate changes to largely account for the change in stability. Thus, contrary to most folding models, widespread helix is not yet formed at the rate-limiting step in the folding pathway. The high-energy transition state is a collapsed form that contains little if any secondary structure, as suggested for the globular protein cytochrome c (Sosnick et al., Proteins 24: 413-426, 1996).

Mechanisms and uses of hydrogen exchange

Recent work has largely completed our understanding of the hydrogen-exchange chemistry of unstructured proteins and nucleic acids. Some of the high-energy structural fluctuations that determine the hydrogen-exchange behavior of native macromolecules have been explained; others remain elusive. A growing number of applications are exploiting hydrogen-exchange behavior to study difficult molecular systems and elicit otherwise inaccessible information on protein structure, dynamics and energetics.

Protein folding intermediates: native-state hydrogen exchange

The hydrogen exchange behavior of native cytochrome c in low concentrations of denaturant reveals a sequence of metastable, partially unfolded forms that occupy free energy levels reaching up to the fully unfolded state. The step from one form to another is accomplished by the unfolding of one or more cooperative units of structure. The cooperative units are entire omega loops or mutually stabilizing pairs of whole helices and loops. The partially unfolded forms detected by hydrogen exchange appear to represent the major intermediates in the reversible, dynamic unfolding reactions that occur even at native conditions and thus may define the major pathway for cytochrome c folding.

Structures of fd gene 5 protein.nucleic acid complexes: a combined solution scattering and electron microscopy study

Small-angle scattering and electron microscopy studies of fd gene 5 protein (g5p) and reconstituted g5p.nucleic acid complexes have been used to test models for the complexes and evaluate their uniqueness. In addition, we have obtained new information on the dependence of nucleotide type and protein/nucleotide (P/N) ratio on the structure of the complexes. Reconstituted complexes were made with single-stranded fd viral DNA (fd ssDNA), poly[d(A)] and poly[r(A)]. All complexes form similar left-handed, flexible superhelices having approximately the same diameter, but the pitch differs among these complexes. The g5p protein is a dimer in solution and the dimers associate to form a superhelical framework to which the polynucleotide is attached. The combined X-ray and neutron scattering data confirm the nucleic acid is inside the protein superhelix. A Monte Carlo integration modeling procedure applied to the scattering data was used to systematically test large numbers of possible models for each complex, and previously proposed models based on parameters obtained from electron microscopy were found to be essentially correct and unique. The data on the complexes with different P/N ratios showed that mass per unit length values decreased while the rise per dimer and pitch of the superhelix increased for g5p.fd-ssDNA complexes with decreasing P/N ratios.

The barriers in protein folding

Elimination of an interaction which forms in denatured cytochrome c enables the majority of the molecules to fold to the native state on a 15 ms time scale, without populating observable intermediates. These results are contrary to the current view that particular steps in protein folding, including the supposedly rate-limiting molten globule to native transition, are intrinsically slow. Instead it appears that intermediates characterized so far may be kinetically trapped by barriers that are optional rather than integral to the folding process. Major barriers may result from misorganization of the chain in the initial condensation step.

Solution structure of the cAMP-dependent protein kinase catalytic subunit and its contraction upon binding the protein kinase inhibitor peptide

Small-angle X-ray scattering and Fourier transform infrared (FTIR) spectroscopy experiments have been completed on the catalytic subunit of the cAMP-dependent protein kinase. Measurements were made both with and without the protein kinase inhibitor peptide, PKI alpha(5-22)amide. Binding of the peptide results in an overall contraction of the structure that is characterized by a decrease of 9% in radius of gyration and about 16% in the maximum linear dimension. Both the secondary structure content of the protein/peptide complex, as determined by FTIR, and the solution structure of this binary complex, as determined by X-ray scattering, agree well with the structural characteristics of this complex as elucidated by the crystal structure [Knighton, D.R., Zheng, J., Ten Eyck, L. F., Ashford, V.A., Xuong, N.H., Taylor, S.S., & Sowadsi, J. M. (1991a) Science 253, 407-414]. Further, the contraction of the structure observed by X-ray scattering upon inhibitor peptide binding is not accompanied by any detectable change in secondary structure content of the kinase. We have modeled the contraction of the kinase upon inhibitor peptide binding as a simple rotation of the large and small lobes seen in the crystal structure such that the cleft between them is closed. For a substrate these changes would then allow catalysis to ensue. The hinge for this movement occurs around a glycine that is one of the protein kinase family consensus amino acids.

Denatured states of ribonuclease A have compact dimensions and residual secondary structure

Using small-angle X-ray scattering and Fourier transform infrared spectroscopy, we have determined that the thermally denatured state of native ribonuclease A is on average a compact structure having residual secondary structure. Under strongly reducing conditions, the protein further unfolds into a looser structure with larger dimensions but still retains a comparable amount of secondary structure. The dimensions of the thermally and chemically denatured states of the reduced protein are different but both are more compact than is predicted for a random coil of the same length. These results demonstrate that thermal denaturation in ribonuclease A is not a simple two-state transition from a native to a completely disordered random coil state.

Distances between the antigen-binding sites of three murine antibody subclasses measured using neutron and X-ray scattering

For three different murine immunoglobulins (IgG subclasses 1, 2a, and 2b), the distances between their antigen-binding sites have been measured using neutron scattering from deuterated antigens complexed with proteated IgG. Neutron-scattering data were measured for each antibody-antigen complex in a 41% D2O solvent. Unlike the proteated antibody molecule, the deuterated antigens are strongly contrasted against the 41% D2O solvent and give rise to a scattering profile that contains an interference term related to the distance between the deuterated antigens. For all three subclasses, the damping of this interference term, which gives information on the relative flexibility of the antigen-binding sites, indicates that a single distance is inadequate to describe the observed scattering and a distribution of distances is needed. The scattering profile has been modeled for each subclass to give the mean distance between the antigens and the variance of this distance. For all three IgG subclasses, the mean distance is between 117 and 134 A, and the variance is large (approximately 40 A), indicating a high degree of flexibility of the Fab arms. Small-angle X-ray scattering measurements on the same samples are consistent with the neutron-scattering results.