Photochemical and electrophysical production of radicals on millisecond timescales to probe the structure, dynamics and interactions of proteins.

The reaction of hydroxyl and other oxygen-based radicals with the side chains of proteins on millisecond timescales has been used to probe the structure of proteins, their dynamics in solution and interactions with other macromolecules. Radicals are generated in high flux within microseconds from synchrotron radiation and discharge sources and react with proteins on timescales that are less than those often attributed to structural reorganisation and folding. The oxygen-based radicals generated in aqueous solution react with proteins to effect limited oxidation at specific amino acids throughout the sequence of the protein. The extent of oxidation at these residue markers is highly influenced by the accessibility of the reaction site to the bulk solvent. The extent of oxidation allows protection levels to be measured based on the degree to which a reaction occurs. A map of a protein's three-dimensional structure is subsequently assembled as in a footprinting experiment. Protein solutions that contain various concentrations of substrates that either promote or disrupt structural transitions can be investigated to facilitate site-specific equilibrium and time-resolved studies of protein folding. The radical-based strategies can also be employed in the study of protein-protein interactions to provide a new avenue for investigating protein complexes and assemblies with high structural resolution. The urea-induced unfolding of apomyoglobin, and the binding domains within the ribonuclease S and calmodulin-melittin protein-peptide complexes are presented to illustrate the approach.     

Unfolding of apomyoglobin helices by synchrotron radiolysis and mass spectrometry.

The synchrotron X-ray protein radiolysis technique is based on a quantitative determination of the extent and the site of millisecond radiolytic oxidation of amino-acid side chains by mass spectrometry. The amino acids most susceptible to radiolytic oxidation are cysteine, methionine, phenylalanine, tyrosine, tryptophan, proline, histidine, and leucine. These residues serve as reactive markers within a protein structure that can be used to monitor changes in solvent accessibility during folding or as part of macromolecular interactions. To monitor the unfolding, the extent of radiolytic products of side chains of reactive amino acids is quantitatively measured by mass spectrometry as a function of the denaturant concentration following proteolysis. This approach provides site-specific unfolding isotherms for various segments of a protein without the use of mutation or labeling techniques. Application of this technique to the equilibrium urea unfolding of apomyoglobin at pH 7.8 has demonstrated the cooperative unfolding of helices A to C consistent with midpoints, DeltaG, and m values derived from fluorescence data. The G helix, in contrast, showed a local unfolding behavior. The similarity of the thermodynamic data derived by this synchrotron-based method for helix A (containing two oxidizable tryptophan residues) to that of the fluorescence data indicates that the limited oxidation of proteins by exposure to X-rays on millisecond timescales does not alter the structure of apomyglobin. This supports the viability of the method for the study of protein folding and the mapping of protein interaction sites.     

Reversible unfolding and refolding behavior of a monomeric aldolase from Staphylococcus aureus.

Thermal and GdmCl-induced unfolding transitions of aldolase from Staphylococcus aureus are reversible under a variety of solvent conditions. Analysis of the transitions reveals that no partially folded intermediates can be detected under equilibrium conditions. The stability of the enzyme is very low with a delta G0 value of -9 +/- 2 kJ/mol at 20 degrees C. The kinetics of unfolding and refolding of aldolase are complex and comprise at least one fast and two slow reactions. This complexity arises from prolyl isomerization reactions in the unfolded chain, which are kinetically coupled to the actual folding reaction. Comparison with model calculations shows that at least two prolyl peptide bonds give rise to the observed slow folding reactions of aldolase and that all of the involved bonds are presumably in the trans conformation in the native state. The rate constant of the actual folding reaction is fast with a relaxation time of about 15 s at the midpoint of the folding transition at 15 degrees C. The data presented on the folding and stability of aldolase are comparable to the properties of much smaller proteins. This might be connected with the simple and highly repetitive tertiary structure pattern of the enzyme, which belongs to the group of alpha/beta barrel proteins.     

Kinetic pathways of beta-hairpin (un)folding in explicit solvent

We examine the dynamical (un)folding pathways of the C-terminal beta-hairpin of protein G-B1 at room temperature in explicit solvent, by employing transition path sampling algorithms. The path ensembles contain information on the folding kinetics, including solvent motion. We determine the transition state ensembles for the two main transitions: 1), the hydrophobic collapse; and 2), the backbone hydrogen bond formation. In both cases the transition state ensembles are characterized by a layer (1) or a strip (2) of water molecules in between the two hairpin strands, supporting the hypothesis of the solvent as lubricant in the folding process. The transition state ensembles do not correspond with saddle points in the equilibrium free-energy landscapes. The kinetic pathways are thus not completely determined by the free-energy landscape. This phenomenon can occur if the order parameters obey different timescales. Using the transition interface sampling technique, we calculate the rate constants for (un)folding and find them in reasonable agreement with experiments, thus supporting the validation of using all-atom force fields to study protein folding.     

Rapid unfolding of a domain populates an aggregation-prone intermediate that can be recognized by GroEL

Some amino acid substitutions in phage P22 coat protein cause a temperature-sensitive folding (tsf) phenotype. In vivo, these tsf amino acid substitutions cause coat protein to aggregate and form intracellular inclusion bodies when folded at high temperatures, but at low temperatures the proteins fold properly. Here the effects of tsf amino acid substitutions on folding and unfolding kinetics and the stability of coat protein in vitro have been investigated to determine how the substitutions change the ability of coat protein to fold properly. The equilibrium unfolding transitions of the tsf variants were best fit to a three-state model, N if I if U, where all species concerned were monomeric, a result confirmed by velocity sedimentation analytical ultracentrifugation. The primary effect of the tsf amino acid substitutions on the equilibrium unfolding pathway was to decrease the stability (DeltaG) and the solvent accessibility (m-value) of the N if I transition. The kinetics of folding and unfolding of the tsf coat proteins were investigated using tryptophan fluorescence and circular dichroism (CD) at 222 nm. The tsf amino acid substitutions increased the rate of unfolding by 8-14-fold, with little effect on the rate of folding, when monitored by tryptophan fluorescence. In contrast, when folding or unfolding reactions were monitored by CD, the reactions were too fast to be observed. The tsf coat proteins are natural substrates for the molecular chaperones, GroEL/S. When native tsf coat protein monomers were incubated with GroEL, they bound efficiently, indicating that a folding intermediate was significantly populated even without denaturant. Thus, the tsf coat proteins aggregate in vivo because of an increased propensity to populate this unfolding intermediate.     

Thermodynamic basis of site-specific cooperativity

Cooperative phenomena in biological macromolecules arise from the interaction of many distinct subsystems, such as structural domains or binding sites. Cooperative properties of the system as a whole, like protein folding or allosteric transitions, are subject to the restrictions imposed by thermodynamic stability. These restrictions, however, do not apply in the case of individual subsystems open to interactions with the rest of the macromolecule. The site-specific properties of such subsystems can be understood in general thermodynamic terms from those of a multicomponent system under particular conditions. The analogy provides a thermodynamic basis for site-specific Cooperativity. © 1994 John Wiley & Sons, Inc.     

Synchrotron protein footprinting: a technique to investigate protein-protein interactions.

Traditional approaches for macromolecular structure elucidation, including NMR, crystallography and cryo-EM have made significant progress in defining the structures of protein-protein complexes. A substantial number of macromolecular structures, however, have not been examined with atomic detail due to sample size and heterogeneity, or resolution limitations of the technique; therefore, the general applicability of each method is greatly reduced. Synchrotron footprinting attempts to bridge the gap in these methods by monitoring changes in accessible surface areas of discrete macromolecular moieties. As evidenced by our previous studies on RNA folding and DNA-protein interactions, the three-dimensional structure is probed by examining the reactions of these moieties with hydroxyl radicals generated by synchrotron X-rays. Here we report the application of synchrotron footprinting to the investigation of protein- protein interactions, as the novel technique has been utilized to successfully map the contact sites of gelsolin segment-1 in the gelsolin segment 1/actin complex. Footprinting results demonstrate that phenylalanine 104, located on the actin binding helix of gelsolin segment 1, is protected from hydroxyl radical modification in the presence of actin. This change in reactivity results from the specific protection of gelsolin segment-1, consistent with the substantial decrease in solvent accessibility of F104 upon actin binding, as calculated from the crystal structural of the gelsolin segment 1/actin complex. The results presented here establish synchrotron footprinting as a broadly applicable method to probe structural features of macromolecular complexes that are not amenable to conventional approaches.     

Thermodynamics of a diffusional protein folding reaction

The folding reactions of several proteins are well described as diffusional barrier crossing processes, which suggests that they should be analyzed by Kramers' rate theory rather than by transition state theory. For the cold shock protein Bc-Csp from Bacillus caldolyticus, we measured stability and folding kinetics, as well as solvent viscosity as a function of temperature and denaturant concentration. Our analysis indicates that diffusional folding reactions can be treated by transition state theory, provided that the temperature and denaturant dependence of the solvent viscosity is properly accounted for, either at the level of the measured rate constants or of the calculated activation parameters. After viscosity correction the activation barriers for folding become less enthalpic and more entropic. The transition from an enthalpic to an entropic folding barrier with increasing temperature is, however, apparent in the data before and after this correction. It is a consequence of the negative activation heat capacity of refolding, which is independent of solvent viscosity. Bc-Csp and its mesophilic homolog Bs-CspB from Bacillus subtilis differ strongly in stability but show identical enthalpic and entropic barriers to refolding. The increased stability of Bc-Csp originates from additional enthalpic interactions that are established after passage through the activated state. As a consequence, the activation enthalpy of unfolding is increased relative to Bs-CspB.     

Protein folding observed by time-resolved synchrotron x-ray scattering. A feasibility study.

A solution to the "protein folding" problem, the successful prediction of tertiary and quaternary protein structure from amino acid or gene sequence, would be a major advance in biology and biotechnology. Knowledge of any intermediate structure between fully unwound and folded would aid folding calculations. The use of high intensity synchrotron x-rays from the SUNY X21 beamline at National Synchrotron Light Source has been investigated as a probe of structural changes during protein folding and unfolding in solution. A temperature jump apparatus was used to study thermally-induced folding and unfolding. Scattering of solutions of myoglobin in the angular range 20 = 1-50 mrad. was measured during temperature jumps between 26 and 76 degrees C. There are clear signs of time/temperature-dependent structural changes, in the small angle region, consistent with those from other equilibrium techniques. Analysis indicates that this experimental technique can be extended to the higher angle region where theoretical calculations indicate more detailed structural information, for example when alpha-helix formation, is present.     

Local kinetic measures of macromolecular structure reveal partitioning among multiple parallel pathways from the earliest steps in the folding of a large RNA molecule

At the heart of the RNA folding problem is the number, structures, and relationships among the intermediates that populate the folding pathways of most large RNA molecules. Unique insight into the structural dynamics of these intermediates can be gleaned from the time-dependent changes in local probes of macromolecular conformation (e.g. reports on individual nucleotide solvent accessibility offered by hydroxyl radical (()OH) footprinting). Local measures distributed around a macromolecule individually illuminate the ensemble of separate changes that constitute a folding reaction. Folding pathway reconstruction from a multitude of these individual measures is daunting due to the combinatorial explosion of possible kinetic models as the number of independent local measures increases. Fortunately, clustering of time progress curves sufficiently reduces the dimensionality of the data so as to make reconstruction computationally tractable. The most likely folding topology and intermediates can then be identified by exhaustively enumerating all possible kinetic models on a super-computer grid. The folding pathways and measures of the relative flux through them were determined for Mg(2+) and Na(+)-mediated folding of the Tetrahymena thermophila group I intron using this combined experimental and computational approach. The flux during Mg(2+)-mediated folding is divided among numerous parallel pathways. In contrast, the flux during the Na(+)-mediated reaction is predominantly restricted through three pathways, one of which is without detectable passage through intermediates. Under both conditions, the folding reaction is highly parallel with no single pathway accounting for more than 50% of the molecular flux. This suggests that RNA folding is non-sequential under a variety of different experimental conditions even at the earliest stages of folding. This study provides a template for the systematic analysis of the time-evolution of RNA structure from ensembles of local measures that will illuminate the chemical and physical characteristics of each step in the process. The applicability of this analysis approach to other macromolecules is discussed.     

Synchrotron Protein Footprinting: A Technique to Investigate Protein-Protein Interactions

Abstract

Traditional approaches for macromolecular structure elucidation, including NMR, crystallography and cryo-EM have made significant progress in defining the structures of protein-protein complexes. A substantial number of macromolecular structures, however, have not been examined with atomic detail due to sample size and heterogeneity, or resolution limitations of the technique; therefore, the general applicability of each method is greatly reduced. Synchrotron footprinting attempts to bridge the gap in these methods by monitoring changes in accessible surface areas of discrete macromolecular moieties. As evidenced by our previous studies on RNA folding and DNA-protein interactions, the three-dimensional structure is probed by examining the reactions of these moieties with hydroxyl radicals generated by synchrotron X-rays. Here we report the application of synchrotron foot- printing to the investigation of protein-protein interactions, as the novel technique has been utilized to successfully map the contact sites of gelsolin segment-1 in the gelsolin segment 1/actin complex. Footprinting results demonstrate that phenylalanine 104, located on the actin binding helix of gelsolin segment 1, is protected from hydroxyl radical modification in the presence of actin. This change in reactivity results from the specific protection of gel- solin segment-1, consistent with the substantial decrease in solvent accessibility of F104 upon actin binding, as calculated from the crystal structural of the gelsolin segment 1/actin complex. The results presented here establish synchrotron footprinting as a broadly applicable method to probe structural features of macromolecular complexes that are not amenable to conventional approaches.     

Small-angle scattering: a view on the properties, structures and structural changes of biological macromolecules in solution

A self-contained presentation of the main concepts and methods for interpretation of X-ray and neutron-scattering patterns of biological macromolecules in solution, including a reminder of the basics of X-ray and neutron scattering and a brief overview of relevant aspects of modern instrumentation, is given. For monodisperse solutions the experimental data yield the scattering intensity of the macromolecules, which depends on the contrast between the solvent and the particles as well as on their shape and internal scattering density fluctuations, and the structure factor, which is related to the interactions between macromolecules. After a brief analysis of the information content of the scattering intensity, the two main approaches for modelling the shape and/or structure of macromolecules and the global minimization schemes used in the calculations are presented. The first approach is based, in its more advanced version, on the spherical harmonics approximation and relies on few parameters, whereas the second one uses bead models with thousands of parameters. Extensions of bead modelling can be used to model domain structure and missing parts in high-resolution structures. Methods for computing the scattering patterns from atomic models including the contribution of the hydration shell are discussed and examples are given, which also illustrate that significant differences sometimes exist between crystal and solution structures. These differences are in some cases explainable in terms of rigid-body motions of parts of the structures. Results of two extensive studies--on ribosomes and on the allosteric protein aspartate transcarbamoylase--illustrate the application of the various methods. The unique bridge between equilibrium structures and thermodynamic or kinetic aspects provided by scattering techniques is illustrated by modelling of intermolecular interactions, including crystallization, based on an analysis of the structure factor and recent time-resolved work on assembly and protein folding.     

Determination of thermodynamics and kinetics of RNA reactions by force

Single-molecule methods have made it possible to apply force to an individual RNA molecule. Two beads are attached to the RNA; one is on a micropipette, the other is in a laser trap. The force on the RNA and the distance between the beads are measured. Force can change the equilibrium and the rate of any reaction in which the product has a different extension from the reactant. This review describes use of laser tweezers to measure thermodynamics and kinetics of unfolding/refolding RNA. For a reversible reaction the work directly provides the free energy; for irreversible reactions the free energy is obtained from the distribution of work values. The rate constants for the folding and unfolding reactions can be measured by several methods. The effect of pulling rate on the distribution of force-unfolding values leads to rate constants for unfolding. Hopping of the RNA between folded and unfolded states at constant force provides both unfolding and folding rates. Force-jumps and force-drops, similar to the temperature jump method, provide direct measurement of reaction rates over a wide range of forces. The advantages of applying force and using single-molecule methods are discussed. These methods, for example, allow reactions to be studied in non-denaturing solvents at physiological temperatures; they also simplify analysis of kinetic mechanisms because only one intermediate at a time is present. Unfolding of RNA in biological cells by helicases, or ribosomes, has similarities to unfolding by force.     

Characterization of a quaternary-structured folding intermediate of an antibody Fab-fragment.

Antibody folding is a complex process comprising folding and association reactions. Although it is usually difficult to characterize kinetic folding intermediates, in the case of the antibody Fab fragment, domain-domain interactions lead to a rate-limiting step of folding, thus accumulating folding intermediates at a late step of folding. Here, we analyzed a late folding intermediate of the Fab fragment of the monoclonal antibody MAK 33 from mouse (kappa/IgG1). As a strategy for accumulation of this intermediate we used partial denaturation of the native Fab by guanidinium chloride. This denaturation intermediate, which can be populated to about 90%, is indistinguishable from a late-folding intermediate with respect to denaturation and renaturation kinetics. The spectroscopic analysis reveals a native-like secondary structure of this intermediate with aromatic side chains only slightly more solvent exposed than in the native state. The respective partner domains are weekly associated. From these data we conclude that the intramolecular association of the two chains during folding, with all domains in a native-like structure, follows a two-step mechanism. In this mechanism, presumably hydrophobic interactions are followed by rearrangements leading to the exact complementarity of the contact sites of the respective domains.     

Microsecond folding of the cold shock protein measured by a pressure-jump technique

A pressure-jump apparatus was employed in investigating the kinetics of protein unfolding and refolding. In the reaction cell, the pressure can be increased or decreased by 100-160 bar within 50-100 microseconds and then held constant. Thus, unfolding and refolding reactions in the time range from 70 microseconds to 70 s can be followed with this technique. Measurements are possible in the transition regions of thermally or denaturant-induced folding in a wide range of temperatures and solvent conditions. We used this pressure-jump method to determine the temperature dependence of the rate constants of unfolding and refolding of the cold shock protein of Bacillus subtilis and of three variants thereof with Phe --> Ala substitutions in the central beta-sheet region. For all variants, the change in heat capacity occurred in refolding between the unfolded and activated states, suggesting that the overall native-like character of the activated state of folding was not changed by the deletion of individual Phe side chains. The Phe27Ala mutation affected the rate of unfolding only; the Phe15Ala and Phe17Ala mutations changed the kinetics of both unfolding and refolding. Although the activated state of folding of the cold shock protein is overall native-like, individual side chains are still in a non-native environment.     

Kinetic Folding Mechanism of Erythropoietin

This report describes what to our knowledge is the first kinetic folding studies of erythropoietin, a glycosylated four-helical bundle cytokine responsible for the regulation of red blood cell production. Kinetic responses for folding and unfolding reactions initiated by manual mixing were monitored by far-ultraviolet circular dichroism and fluorescence spectroscopy, and folding reactions initiated by stopped-flow mixing were monitored by fluorescence. The urea concentration dependence of the observed kinetics were best described by a three-state model with a transiently populated intermediate species that is on-pathway and obligatory. This folding scheme was further supported by the excellent agreement between the free energy of unfolding and m-value calculated from the microscopic rate constants derived from this model and these parameters determined from separate equilibrium unfolding experiments. Compared to the kinetics of other members of the four-helical bundle cytokine family, erythropoietin folding and unfolding reactions were slower and less susceptible to aggregation. We tentatively attribute these slower rates and protection from association events to the large amount of carbohydrate attached to erythropoietin at four sites.     

The Folding Free-Energy Surface of HIV-1 Protease: Insights into the Thermodynamic Basis for Resistance to Inhibitors

Spontaneous mutations at numerous sites distant from the active site of human immunodeficiency virus type 1 protease enable resistance to inhibitors while retaining enzymatic activity. As a benchmark for probing the effects of these mutations on the conformational adaptability of this dimeric β-barrel protein, the folding free-energy surface of a pseudo-wild-type variant, HIV-PR^?, was determined by a combination of equilibrium and kinetic experiments on the urea-induced unfolding/refolding reactions. The equilibrium unfolding reaction was well described by a two-state model involving only the native dimeric form and the unfolded monomer. The global analysis of the kinetic folding mechanism reveals the presence of a fully folded monomeric intermediate that associates to form the native dimeric structure. Independent analysis of a stable monomeric version of the protease demonstrated that a small-amplitude fluorescence phase in refolding and unfolding, not included in the global analysis of the dimeric protein, reflects the presence of a transient intermediate in the monomer folding reaction. The partially folded and fully folded monomers are only marginally stable with respect to the unfolded state, and the dimerization reaction provides a modest driving force at micromolar concentrations of protein. The thermodynamic properties of this system are such that mutations can readily shift the equilibrium from the dimeric native state towards weakly folded states that have a lower affinity for inhibitors but that could be induced to bind to their target proteolytic sites. Presumably, subsequent secondary mutations increase the stability of the native dimeric state in these variants and, thereby, optimize the catalytic properties of the resistant human immunodeficiency virus type 1 protease.     

Experimental approaches for solution X-ray scattering and fiber diffraction

X-ray scattering and diffraction from non-crystalline systems have gained renewed interest in recent years, as focus shifts from the structural chemistry information gained by high-resolution studies to the context of structural physiology at larger length scales. Such techniques permit the study of isolated macromolecules as well as highly organized macromolecular assemblies as a whole under near-physiological conditions. Time-resolved approaches, made possible by advanced synchrotron instrumentation, add a crucial dimension to many of these investigations. This article reviews experimental approaches in non-crystalline X-ray scattering and diffraction that may be used to illuminate important scientific questions such as protein/nucleic acid folding and structure-function relationships in large macromolecular assemblies.     

Complex folding pathways in a simple beta-hairpin

The determination of the folding mechanisms of proteins is critical to understand the topological change that can propagate Alzheimer and Creutzfeld-Jakobs diseases, among others. The computational community has paid considerable attention to this problem; however, the associated time scale, typically on the order of milliseconds or more, represents a formidable challenge. Ab initio protein folding from long molecular dynamics simulations or ensemble dynamics is not feasible with ordinary computing facilities and new techniques must be introduced. Here we present a detailed study of the folding of a 16-residue beta-hairpin, described by a generic energy model and using the activation-relaxation technique. From a total of 90 trajectories at 300 K, three folding pathways emerge. All involve a simultaneous optimization of the complete hydrophobic and hydrogen bonding interactions. The first two pathways follow closely those observed by previous theoretical studies (folding starting at the turn or by interactions between the termini). The third pathway, never observed by previous all-atom folding, unfolding, and equilibrium simulations, can be described as a reptation move of one strand of the beta-sheet with respect to the other. This reptation move indicates that non-native interactions can play a dominant role in the folding of secondary structures. Furthermore, such a mechanism mediated by non-native hydrogen bonds is not available for study by unfolding and Gō model simulations. The exact folding path followed by a given beta-hairpin is likely to be influenced by its sequence and the solvent conditions. Taken together, these results point to a more complex folding picture than expected for a simple beta-hairpin.     

Environment and exposure to solvent of protein atoms. Lysozyme and insulin

A computer program is described for calculating the environment and the exposure to solvent of atoms of a protein. The computation is based on the atomic co-ordinates of the protein and on assumptions like those of Lee &; Richards (1971). Results for lysozyme and insulin are presented. Changes in exposure to solvent and in the nature of contacts that develop through folding, association reactions and crystallization are described numerically. The computations suggest several generalizations. (a) Lattice contacts within the protein crystal are characterized by a significantly smaller involvement of non-polar side chains and a proportionately greater involvement of ionizable side chains than is found for protein folding or for protein association reactions important for biological function, (b) In helical regions the carbonyl oxygen of the first residue in the helix has high probability of being shielded from solvent, (c) Glycine is among the residues having exposure least affected by folding; this accords with the expectation that it lies at bends of the peptide chain on the surface of the molecule.