A comparison of the folding kinetics of a small,artificially selected DNA aptamer with those of equivalently simple naturally occurring proteins

The folding of larger proteins generally differs from the folding of similarly large nucleic acids in the number and stability of the intermediates involved. To date, however, no similar comparison has been made between the folding of smaller proteins, which typically fold without well-populated intermediates, and the folding of small, simple nucleic acids. In response, in this study, we compare the folding of a 38-base DNA aptamer with the folding of a set of equivalently simple proteins. We find that, as is true for the large majority of simple, single domain proteins, the aptamer folds through a concerted, millisecond-scale process lacking well-populated intermediates. Perhaps surprisingly, the observed folding rate falls within error of a previously described relationship between the folding kinetics of single-domain proteins and their native state topology. Likewise, similarly to single-domain proteins, the aptamer exhibits a relatively low urea-derived Tanford β, suggesting that its folding transition state is modestly ordered. In contrast to this, however, and in contrast to the behavior of proteins, ϕ-value analysis suggests that the aptamer''s folding transition state is highly ordered, a discrepancy that presumably reflects the markedly more important role that secondary structure formation plays in the folding of nucleic acids. This difference notwithstanding, the similarities that we observe between the two-state folding of single-domain proteins and the two-state folding of this similarly simple DNA presumably reflect properties that are universal in the folding of all sufficiently cooperative heteropolymers irrespective of their chemical details.     

Temperature-gradient gel electrophoresis. Thermodynamic analysis of nucleic acids and proteins in purified form and in cellular extracts

A temperature-gradient gel electrophoresis technique and its application to the study of structural transitions of nucleic acids and protein-nucleic acid complexes are described. The temperature gradient is established in a slab gel by means of a simple ancillary device for a commercial horizontal gel apparatus. The gradient may be freely selected between 10 and 80 degrees C, and is highly reproducible and linear. In a normal application the biopolymers migrate perpendicular to the temperature gradient so that every individual molecule is at constant temperature throughout electrophoresis. The structural transition of a biopolymer is seen as a continuous band which is retarded or speeded up in the temperature range of the transition. Dissociation processes are mostly irreversible under the conditions of electrophoresis and, therefore, show up as discontinuous transitions from a slow-moving to fast-moving band. As examples the conformational transitions of viroids, double-stranded RNA from reovirus, double-stranded satellite RNA from cucumber mosaic virus and repressor-operator complexes have been studied. It could be shown that by this method dsRNA molecules may be differentiated which differ only in one base-pair, or proteins differing in one amino acid only. As a particular advantage, temperature-gradient gel electrophoresis allows the study of conformational transitions of biopolymers which have not been purified. The biopolymer may either be identified by silver staining as a specific band among many others or, if the study is carried out on nucleic acids, these may be recorded by hybridization with a radioactive probe.     

Influence of lipid membranes on the conformational transitions of nucleic acids

The conformational transitions of nucleic acids which were enclosed in reverse phase evaporation vesicles (REV) were studied by thermal denaturation with optical recording. Cloned fragments of double-stranded DNA containing 179 base pairs and 187 base pairs, respectively, and polyA.polyU were enclosed in REV with a yield up to every vesicle containing 50 nucleic acid molecules. With the 179 base pairs DNA enclosed in the vesicle from egg lecithin two well resolved helix-coil transitions could be measured; one is very similar in the midpoint-temperature Tm and halfwidth delta T1/2 to the transition of the free nucleic acid, and the other transition occurs stabilized at a 3.5 degrees C higher Tm-value and with a broader delta T1/2, 2.7 degrees C instead of 0.6 degree C. Both transitions are from nucleic acids inside the vesicles. Varying the surface charge of the lipid membrane by adding the negatively charged phosphatidylserine or phosphatidylglycerol, an optimum in the yield of enclosure and a maximum in the increase in Tm (4.5 degrees C) and delta T1/2 (5.5 degrees C instead of 1.0 degrees C) was obtained at 20% phosphatidylserine or phosphatidylglycerol. In vesicles from pure negatively charged lipids no second population of nucleic acids was observed. Qualitatively, similar effects were observed with polyA.polyU. Stabilization and broadening of the second transition is higher for nucleic acids inside vesicles from lipids with unsaturated fatty acids, as dioleoyl-phosphatidylcholine, than with saturated fatty acids, dipalmitoyl-phosphatidylcholine. Stabilization and broadening decrease with increasing ionic strength, whereas the relative contributions of both transitions to the total hypochromicity remain unchanged; the second transition coincides with the first at 90 mM Na+. From the experimental results it was concluded that the interaction of nucleic acids and lipid membranes is mainly of electrostatic nature. The nucleic acids exist inside the vesicles in two populations, one behaving like nucleic acid free in solution and one influenced by the contact with the membrane. All results are in accordance with a model in which the interaction between the nucleic acid and the membrane is in competition with the dipole-dipole interaction inside the membrane surface.     

Why do proteins divide into domains? Insights from lattice model simulations

It is known that larger globular proteins are built from domains, relatively independent structural units. A domain size seems to be limited, and a single domain consists of from few tens to a couple of hundred amino acids. Based on Monte Carlo simulations of a reduced protein model restricted to the face centered simple cubic lattice, with a minimal set of short-range and long-range interactions, we have shown that some model sequences upon the folding transition spontaneously divide into separate domains. The observed domain sizes closely correspond to the sizes of real protein domains. Short chains with a proper sequence pattern of the hydrophobic and polar residues undergo a two-state folding transition to the structurally ordered globular state, while similar longer sequences follow a multistate transition. Homopolymeric (uniformly hydrophobic) chains and random heteropolymers undergo a continuous collapse transition into a single globule, and the globular state is much less ordered. Thus, the factors responsible for the multidomain structure of proteins are sufficiently long polypeptide chain and characteristic, protein-like, sequence patterns. These findings provide some hints for the analysis of real sequences aimed at prediction of the domain structure of large proteins.     

Pressure-temperature phase diagrams of biomolecules

The pressure-temperature phase diagram of various biomolecules is reviewed. Special attention is focused on the elliptic phase diagram of proteins. The phenomenological thermodynamic theory describing this diagram explains the heat, cold and pressure denaturations in a unified picture. The limitations and possible developments of this theory are discussed as well. It is pointed out that a more complex diagram can be obtained when the intermolecular interactions are also taken into account. In this case metastable states appear on the pressure-temperature (p-T) diagram due to intermolecular interactions. Pressure-temperature phase diagrams of other biopolymers are also discussed. While the p-T diagrams of helix-coil transition of nucleic acids and of gel-liquid crystal transition of lipid bilayers are non-elliptical, those of gelatinization of starch and of phase separation of some synthetic polymers show an elliptic profile, similar to that of proteins. Finally, the p-T diagram of bacterial inactivation is shown to be elliptic. From the point of view of basic science, this fact shows that the key factor of inactivation should be the protein type, and from the viewpoint of practical applications, it serves as the theoretical basis of pressure treatment of biosystems.     

Mechanical unfolding of human telomere G-quadruplex DNA probed by integrated fluorescence and magnetic tweezers spectroscopy

Single-molecule techniques facilitate analysis of mechanical transitions within nucleic acids and proteins. Here, we describe an integrated fluorescence and magnetic tweezers instrument that permits detection of nanometer-scale DNA structural rearrangements together with the application of a wide range of stretching forces to individual DNA molecules. We have analyzed the force-dependent equilibrium and rate constants for telomere DNA G-quadruplex (GQ) folding and unfolding, and have determined the location of the transition state barrier along the well-defined DNA-stretching reaction coordinate. Our results reveal the mechanical unfolding pathway of the telomere DNA GQ is characterized by a short distance (<1 nm) to the transition state for the unfolding reaction. This mechanical unfolding response reflects a critical contribution of long-range interactions to the global stability of the GQ fold, and suggests that telomere-associated proteins need only disrupt a few base pairs to destabilize GQ structures. Comparison of the GQ unfolded state with a single-stranded polyT DNA revealed the unfolded GQ exhibits a compacted non-native conformation reminiscent of the protein molten globule. We expect the capacity to interrogate macromolecular structural transitions with high spatial resolution under conditions of low forces will have broad application in analyses of nucleic acid and protein folding.     

raf RBD and ubiquitin proteins share similar folds, folding rates and mechanisms despite having unrelated amino acid sequences

Recent experimental and theoretical studies in protein folding suggest that the rates and underlying mechanisms by which proteins attain the native state are largely determined by the topological complexity of a specific fold rather than by the fine details of the amino acid sequences. However, such arguments are based upon the examination of a limited number of protein folds. To test this view, we sought to investigate whether proteins belonging to the ubiquitin superfamily display similar folding behavior. To do so, we compared the folding-unfolding transitions of mammalian ubiquitin (mUbi) with those of its close yeast homologue (yUbi), and to those of the structurally related Ras binding domain (RBD) of the serine/threonine kinase raf that displays no apparent sequence homology with the ubiquitin family members. As demonstrated for mUbi [Krantz, B. A., and Sosnick, T. R. (2000) Biochemistry 39, 11696-11701], we show that a two-state transition model with no burst phase intermediate can describe folding of both yUbi and raf RBD. We further demonstrate that (1) all three proteins refold at rates that are within 1 order of magnitude (1800, 1100, and 370 s(-1) for mUbi, raf RBD, and yUbi, respectively), (2) both mUbi and raf RBD display similar refolding heterogeneity, and (3) the folding free energy barriers of both mUbi and raf RBD display a similar temperature dependence and sensitivity to a stabilizing agent or to mutations of a structurally equivalent central core residue. These findings are consistent with the view that rates and mechanisms for protein folding depend mostly on the complexity of the native structure topology rather than on the fine details of the amino acid sequence.     

Fatty acid binding sites of serum albumin probed by non-linear spin-label EPR

A novel form of non-linear EPR spectroscopy, viz. the first harmonic absorption spectrum recorded in phase quadrature with respect to the Zeeman field modulation, is used here to investigate spin-lattice relaxation enhancements of nitroxide spin labels bound to serum albumin that are induced by spin-spin interactions with aqueous paramagnetic ions. The advantage of this EPR method is that it is directly sensitive to spin-lattice relaxation and affected relatively little by other spectral parameters (Livshits et al., J. Magn. Reson. 133 (1998) 79-91). Relaxation enhancements by ferricyanide of bound fatty acids (n-SASL) spin-labelled at different positions, n, in the chain are compared with those of different maleimide spin label derivatives attached at the single free -SH group, as well as with those of the spin labels free in solution. It was found that: (1) the encounter frequency of ferricyanide with 5-SASL and 12-SASL bound to serum albumin is more than two times less than that with 16-SASL; (2) the accessibility of ferricyanide to 16-SASL is comparable to that of the more immobilised covalently bound spin labels; and (3) the absolute values of the encounter frequencies for the bound spin-labelled fatty acids are approximately a factor of ten smaller than for the corresponding free spin labels, but the latter show a dependence on position of labelling that is similar to the bound labels. A kinetic scheme that is consistent with these relative differences involves rapid reversible transitions between an 'open' and 'closed' state, in which interaction with aqueous paramagnetic agents is possible only in the 'open' state. The equilibrium strongly favours the 'closed' state, which is further enhanced at low temperatures.     

Kinetics of a finite one-dimensional spin system as a model for protein folding

Simple spin models are used to analyze the kinetic nature of lowest energy state formation of the spin systems as models of protein folding kinetics. The models employed in the present paper are based on the spin systems as models of biopolymers previously proposed by the author for the analysis of the equilibrium nature of transitions [T. Kikuchi, Biophys. Chem. 65 (1997) 109]. In particular, the effect of frustrations on the kinetics is investigated with the Monte Carlo simulations in this study. The results show that the kinetics of the present systems are characterized by the ratio of foldables (pathways on the energy landscape that to lead to the lowest energy state) and the temperature dependence of the mean first passage time of foldables. We also discuss the free energy profile of the present models and the relation of the present results to the kinetics of actual proteins.     

Role of a solvent-exposed aromatic cluster in the folding of Escherichia coli CspA

Escherichia coli CspA is a member of the cold shock protein family. All cold shock proteins studied to date fold rapidly by an apparent two-state mechanism. CspA contains an unusual cluster of aromatic amino acids on its surface that is necessary for nucleic acid binding and also provides stability to CspA (Hillier et al., 1998). To elucidate the role this aromatic cluster plays in the determining the folding rate and pathway of CspA, we have studied the folding kinetics of mutants containing either leucine or serine substituted for Phe 18, Phe20, and/or Phe31. The leucine substitutions are found to accelerate folding and the serine substitutions to decelerate folding. Because these residues exert effects on the free energy of the folding transition state, they may be necessary for nucleating folding. They are not responsible, however, for the very compact, native-like transition state ensemble seen in the cold shock proteins, as the refolding rates of the mutants all show a similar, weak dependence of unfolding rate on denaturant concentration. Using mutant cycle analysis, we show that there is energetic coupling among the three residues between the unfolded and transition states, suggesting that the cooperative nature of these interactions helps to determine the unfolding rate. Overall, our results suggest that separate evolutionary pressures can act simultaneously on the same group of residues to maintain function, stability, and folding rate.     

Using chimeric immunity proteins to explore the energy landscape for alpha-helical protein folding

To address the role of sequence in the folding of homologous proteins, the folding and unfolding kinetics of the all-helical bacterial immunity proteins Im2 and Im9 were characterised, together with six chimeric derivatives of these proteins. We show that both Im2 and Im9 fold rapidly (k(UN)(H(2)O)) approximately 2000 s(-1) at pH 7.0, 25 degrees C) in apparent two-state transitions, through rate-limiting transition states that are highly compact (beta(TS)0.93 and 0.96, respectively). Whilst the folding and unfolding properties of three of the chimeras (Im2 (1-44)(Im9), Im2 (1-64)(Im9 )and Im2 (25-44)(Im9)) are similar to their parental counterparts, in other chimeric proteins the introduced sequence variation results in altered kinetic behaviour. At low urea concentrations, Im2 (1-29)(Im9) and Im2 (56-64)(Im9) fold in two-state transitions via transition states that are significantly less compact (beta(TS) approximately 0.7) than those characterised for the other immunity proteins presented here. At higher urea concentrations, however, the rate-limiting transition state for these two chimeras switches or moves to a more compact species (beta(TS) approximately 0.9). Surprisingly, Im2 (30-64)(Im9) populates a highly collapsed species (beta(I)=0.87) in the dead-time (2.5 ms) of stopped flow measurements. These data indicate that whilst topology may place significant constraints on the folding process, specific inter-residue interactions, revealed here through multiple sequence changes, can modulate the ruggedness of the folding energy landscape.     

Hydration and packing along the folding pathway of SH3 domains by pressure-dependent NMR

The volumetric properties associated with protein folding transitions reflect changes in protein packing and hydration of the states that participate in the folding reaction. Here, NMR spin relaxation techniques are employed to probe the folding-unfolding kinetics of two SH3 domains as a function of pressure so that the changes in partial molar volumes along the folding pathway can be measured. The two domains fold with rates that differ by approximately 3 orders of magnitude, so their folding dynamics must be probed using different NMR relaxation experiments. In the case of the drkN SH3 domain that folds via a two-state mechanism on a time scale of seconds, nitrogen magnetization exchange spectroscopy is employed, while for the G48M mutant of the Fyn SH3 domain where the folding occurs on the millisecond time scale (three-step reaction), relaxation dispersion experiments are utilized. The NMR methodology is extremely sensitive to even small changes in equilibrium and rate constants, so reliable estimates of partial molar volumes can be obtained using low pressures (1-120 bar), thus minimizing perturbations to any of the states along the folding reaction coordinate. The volumetric data that were obtained are consistent with a similar folding mechanism for both SH3 domains, involving early chain compaction to states that are at least partially hydrated. This work emphasizes the role of NMR spin relaxation in studying dynamic processes over a wide range of time scales.     

Structural characterization of yeast acidic ribosomal P proteins forming the P1A-P2B heterocomplex

Acidic ribosomal P proteins form a distinct lateral protuberance on the 60S ribosomal subunit. In yeast, this structure is composed of two heterocomplexes (P1A-P2B and P1B-P2A) attached to the ribosome with the aid of the P0 protein. In solution, the isolated P proteins P1A and P2B have a flexible structure with some characteristics of a molten globule [Zurdo, J., et al. (2000) Biochemistry 39, 8935-8943]. In this report, the structure of P1A-P2B heterocomplex from Saccharomyces cerevisiae is investigated by means of size-exclusion chromatography, chemical cross-linking, circular dichroism, light scattering, and fluorescence spectroscopy. The circular dichroism experiment shows that the complex could be ranked in the tertiary class of all-alpha proteins, with an average alpha-helical content of approximately 65%. Heat and urea denaturation experiments reveal that the P1A-P2B complex, unlike the isolated proteins, has a full cooperative transition which can be fitted into a two-state folding-unfolding model. The behavior of the complex in the presence of 2,2,2-trifluoroethanol also resembles a two-state folding-unfolding transition, further supporting the idea that the heterocomplex contains well-packed side chains. In conclusion, the P1A-P2B heterocomplex, unlike the isolated proteins, has a well-defined hydrophobic core. Consequently, the complex can put up its structure without additional ribosomal components, so the heterodimeric complex reflects the intrinsic properties of the two analyzed proteins, indicating thus that this is the only possible configuration of the P1A and P2B proteins on the ribosomal stalk structure.     

How large are the volume changes accompanying protein transitions and binding?

We present a simple model to describe volume changes accompanying protein folding and binding events. The model enables one to resolve the changes in volume accompanying conformational transitions of proteins as well as association of proteins with other molecules in terms of the intrinsic, thermal and interaction (hydration) contributions. The thermal contribution to protein volume results from thermally activated mutual vibrational motions of contacting solute and solvent molecules. Our calculations suggest that near zero volume changes accompanying protein folding and binding events reflect compensation between significant changes in the intrinsic, thermal and interaction terms. We have quantitatively estimated these terms as a function of the protein's molecular weight and degree of its unfolding. Results described in this work lay foundation for more reliable and physically justified interpretations of volumetric data on protein folding and binding events. We also discuss potential ways of extending applications of our model to analyzing other macromolecular systems and events, including drug-DNA and protein-DNA interactions and helix-to-helix and helix-to-coil transitions of nucleic acids.     

Photoactivatable Reagents Based on Aryl(trifluoromethyl)diazirines: Synthesis and Application for Studying Nucleic Acid–Protein Interactions

Recent data on the synthesis of photoactivatable derivatives of nucleic acids and proteins on the basis of aryl(trifluoromethyl)diazirines—analogs of nucleosides, nucleotides, oligonucleotides, as well as amino acids and peptides—are reviewed. The synthesis of bi- and polyfunctional photoactivatable reagents, including those containing a cleavable function, designed for postsynthetic modification of biopolymers is described. Data are given on the use of the photoactivatable derivatives for studying nucleic acid–protein interactions by the method of photoaffinity labeling. Special consideration is paid to the results obtained by the authors" team in cooperation with other researchers as well as graduate students of the Chemistry of Natural Products Chair, Chemical Faculty, Moscow State University and destined to solve various scientific tasks in the domain of nucleic acid–protein recognition with the use of photoaffinity crosslinking.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.     

Correspondence between anomalous m- and DeltaCp-values in protein folding

Proteins folding according to a classical two-state system characteristically show V-shaped chevron plots. We have previously interpreted the symmetrically curved chevron plot of the protein U1A as denaturant-dependent movements in the position of the transition state ensemble (TSE). S6, a structural analog of U1A, shows a classical V-shaped chevron plot indicative of straightforward two-state kinetics, but the mutant LA30 has a curved unfolding limb, which is most consistent with TSE mobility. The kinetic m-values (derivatives of the rate constants with respect to denaturant concentration) in themselves depend on denaturant concentration. To obtain complementary information about putative mobile TSEs, we have carried out a thermodynamic analysis of the three proteins, based on data for refolding and unfolding over the range 10 degrees C to 70 degrees C. The data at all temperatures can be fitted to two-state model systems. Importantly, for all three proteins the activation heat capacities are, within error, identical to the heat capacities measured in independent experiments under equilibrium conditions. Although the equilibrium heat capacities are essentially invariant with regard to denaturant concentration, the activation heat capacities, similar to the structurally equivalent kinetic m-values, show marked denaturant dependence. Furthermore, the values of beta++ at different denaturant concentrations measured by m-values and by heat capacity values are very similar. These observations are consistent with significant transition state movements within the framework of two-state folding. The basis for TSE movement appears to be enthalpic rather than entropic, suggesting that the binding energy of denaturant-protein interactions is a major determinant of the response of energy landscape contours to changing environments.     

Is it always possible to distinguish two- and three-state systems by evaluating the van't Hoff enthalpy?

Many small globular proteins are traditionally classified as thermodynamical two-state systems, i.e., the protein is either in the native, active state (folded) or in the denatured state (unfolded). We challenge this view and show that there may exist (protein) systems for which a van't Hoff analysis of experimental data cannot determine whether the system corresponds to two or three thermodynamical states when only temperatures in a narrow temperature region around the transition are considered. We generalize a widely employed two-state protein folding model to include a third, transition state. For this three-state system we systematically study the deviation of the calorimetric enthalpy (heat of transition) from the van't Hoff enthalpy, a measure of the two-stateness of a transition. We show that under certain conditions the heat capacity of the three-state system can be almost indistinguishable from the heat capacity for the two-state system over a broad temperature interval. The consequence may be that some three-state (or even more than three-states) systems have been misinterpreted as two-state systems when the conclusion is drawn solely upon the van't Hoff enthalpy. These findings are important not only for proteins, but also for the interpretation of thermodynamical systems in general.