Mechanically unfolding the small, topologically simple protein L

beta-sheet proteins are generally more able to resist mechanical deformation than alpha-helical proteins. Experiments measuring the mechanical resistance of beta-sheet proteins extended by their termini led to the hypothesis that parallel, directly hydrogen-bonded terminal beta-strands provide the greatest mechanical strength. Here we test this hypothesis by measuring the mechanical properties of protein L, a domain with a topology predicted to be mechanically strong, but with no known mechanical function. A pentamer of this small, topologically simple protein is resistant to mechanical deformation over a wide range of extension rates. Molecular dynamics simulations show the energy landscape for protein L is highly restricted for mechanical unfolding and that this protein unfolds by the shearing apart of two structural units in a mechanism similar to that proposed for ubiquitin, which belongs to the same structural class as protein L, but unfolds at a significantly higher force. These data suggest that the mechanism of mechanical unfolding is conserved in proteins within the same fold family and demonstrate that although the topology and presence of a hydrogen-bonded clamp are of central importance in determining mechanical strength, hydrophobic interactions also play an important role in modulating the mechanical resistance of these similar proteins.     

Mechanically unfolding proteins: the effect of unfolding history and the supramolecular scaffold

The mechanical resistance of a folded domain in a polyprotein of five mutant I27 domains (C47S, C63S I27)(5)is shown to depend on the unfolding history of the protein. This observation can be understood on the basis of competition between two effects, that of the changing number of domains attempting to unfold, and the progressive increase in the compliance of the polyprotein as domains unfold. We present Monte Carlo simulations that show the effect and experimental data that verify these observations. The results are confirmed using an analytical model based on transition state theory. The model and simulations also predict that the mechanical resistance of a domain depends on the stiffness of the surrounding scaffold that holds the domain in vivo, and on the length of the unfolded domain. Together, these additional factors that influence the mechanical resistance of proteins have important consequences for our understanding of natural proteins that have evolved to withstand force.     

Light-induced unfolding of photoactive yellow protein mutant M100L

Light-activation of the PAS domain protein photoactive yellow protein (PYP) is believed to trigger a negative phototactic response in the phototropic bacterium Halorhodospira halophila. To investigate transient conformational changes of the PYP photocycle, we utilized the PYP mutant M100L that displays an increased lifetime of the putative signaling-state photointermediate PYP(M) by 3 orders of magnitude, as previously reported for the M100A mutant [Devanathan, S., Genick, U. K., Canestrelli, I. L., Meyer, T. E., Cusanovich, M. A., Getzoff, E. D., and Tollin, G. Biochemistry (1998) 37, 11563-11568]. The FTIR difference spectrum of PYP(M) and the ground state of M100L demonstrated extensive peptide-backbone structural changes as observed in the FTIR difference spectrum of the wild-type protein and PYP(M). The conformational change investigated by CD spectroscopy in the far-UV region showed reduction of the alpha-helical content by approximately 40%, indicating a considerable amount of changes in the secondary structure. The optical activity of the p-coumaric acid chromophore completely vanished upon PYP(M) in contrast to the dark state, indicating deformation of the binding pocket structure in PYP(M). The tertiary structural changes were further monitored by small-angle X-ray scattering measurements, which demonstrated a significant increase of the radius of gyration of the molecule by approximately 5% in PYP(M). These structural changes were reversed concomitantly with the chromophore anionization upon the dark state recovery. The observed changes of the quantities provided a more vivid view of the structural changes of the mutant PYP in going from PYP(M) to PYP(dark), which can be regarded as a process of folding of the secondary and the tertiary structures of the "PAS" domain structure, coupled with the p-coumaric acid chromophore deprotonation and isomerization.     

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq-UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq-UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq-UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq-UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq-UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq-UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq-UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq-UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq-UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq-UIM, appears to behave very close to an equilibrium two-state.     

Kinetic study of protein unfolding and refolding using urea gradient electrophoresis

When total cytoplasmic RNA from mouse Friend cells is fractionated using oligo(dT)-cellulose or poly(U)-Sepharose chromatography, approximately 20% of the messenger RNA activity (as measured in the reticulocyte lysate cell-free system) remains in the unbound fraction, even though this contains < 0.5% of the poly(A) (as measured by titration with poly(U)). This RNA, operationally defined as poly(A)^?, is found almost entirely in polysome structures in vivo. Its major translation products, as shown by one-dimensional sodium dodecyl sulphate-containing gels, are the histones and actin. Two-dimensional gels (isoelectric focusing: sodium dodecyl sulphate/gel electrophoresis) show that, with the exception of the mRNAs coding for histones, poly(A)^? mRNA encodes similar proteins to poly(A)⁺ mRNA, though in very different abundances. This is directly confirmed by the arrest of the translation of the abundant poly(A)^? mRNAs after hybridization with a complementary DNA transcribed from poly(A)⁺ RNA.RNA sequences which are rare in the poly(A)⁺ RNA are also found in poly(A)^? RNA, as shown by hybridizing a cDNA transcribed from poly(A)⁺ RNA to total and poly(A)^? polysomal RNA. That this does not simply represent a flow-through of poly(A)⁺ RNA is indicated by (i) the lack of poly(A) by hybridizing to poly(U) in this fraction, (ii) the fact that further passage through poly(U)-Sepharose does not remove the hybridizing sequences, (iii) the very different quantitative distribution of proteins encoded by poly(A)⁺ and poly(A)^? RNAs. We also think that it does not result from removal of poly(A) from polyadenylated RNAs during extraction because RNAs prepared using the minimum of manipulations give similar results. The distribution of both total mRNA and α and β globin mRNAs between poly(A)⁺ and poly(A)^? RNA does not change significantly during the dimethyl sulphoxide-induced differentiation of Friend cells.     

Mechanical unfolding of a beta-hairpin using molecular dynamics

Single-molecule mechanical unfolding experiments have the potential to provide insights into the details of protein folding pathways. To investigate the relationship between force-extension unfolding curves and microscopic events, we performed molecular dynamics simulations of the mechanical unfolding of the C-terminal hairpin of protein G. We have studied the dependence of the unfolding pathway on pulling speed, cantilever stiffness, and attachment points. Under conditions that generate low forces, the unfolding trajectory mimics the untethered, thermally accessible pathway previously proposed based on high-temperature studies. In this stepwise pathway, complete breakdown of backbone hydrogen bonds precedes dissociation of the hydrophobic cluster. Under more extreme conditions, the cluster and hydrogen bonds break simultaneously. Transitions between folding intermediates can be identified in our simulations as features of the calculated force-extension curves.     

Equilibrium folding and unfolding pathways for a model protein

The folding–unfolding process of reduced bovine pancreatic trypsin inhibitor was investigated with an idealized model employing approximate free energies. The protein is regarded to consist of only C^α and C^β atoms. The backbone dihedral angles are the only conformational variables and are permitted to take discrete values at every 10°. Intraresidue energies consist of two terms: an empirical part taken from the observed frequency distributions of (?,ψ) and an additional favorable energy assigned to the native conformation of each residue. Interresidue interactions are simplified by assuming that there is an attractive energy operative only between residue pairs in close contact in the native structure. A total of 230,000 molecular conformations, with no atomic overlaps, ranging from the native state to the denatured state, are randomly generated by changing the sampling bias. Each conformation is classified according to its conformational energy, F; a conformational entropy, S(F) is estimated for each value of F from the number of samples. The dependence of S(F) on energy reveals that the folding–unfolding transition for this idealized model is an “all-or-none” type; this is attributable to the specific long-range interactions. Interresidue contact probabilities, averaged over samples representing various stages of folding, serve to characterize folding intermediates. Most probable equilibrium pathways for the folding–unfolding transition are constructed by connecting conformationally similar intermediates. The specific details obtained for bovine pancreatic trypsin inhibitor are as follows: (1) Folding begins with the appearance of nativelike medium-range contacts at a β-turn and at the α-helix. (2) These grow to include the native pair of interacting β-strands. This state includes intact regular secondary conformations, as well as the interstrand sheet contacts, and corresponds to an activated state with the highest free energy on the pathway. (3) Additional native long-range contacts are completely formed either toward the amino terminus or toward the carboxyl terminus. (4) In a final step, the missing contacts appear. Although these folding pathways for this model are not consistent with experimental reports, it does indicate multiple folding pathways. The method is general and can be applied to any set of calculated conformational energies and furthermore permits investigation of gross folding features.     

The unfolding action of GroEL on a protein substrate

A molecular dynamics simulation of the active unfolding of denatured rhodanese by the chaperone GroEL is presented. The compact denatured protein is bound initially to the cis cavity and forms stable contacts with several of the subunits. As the cis ring apical domains of GroEL undergo the transition from the closed to the more open (ATP-bound) state, they exert a force on rhodanese that leads to the increased unfolding of certain loops. The contacts between GroEL and rhodanese are analyzed and their variation during the GroEL transition is shown. The major contacts, which give rise to the stretching force, are found to be similar to those observed in crystal structures of peptides bound to the apical domains. The results of the simulation show that multidomain interactions play an essential role, in accord with experiments. Implications of the results for mutation experiments and for the action of GroEL are discussed.     

Native-state interconversion of a metamorphic protein requires global unfolding

Lymphotactin (Ltn) is a unique chemokine that under physiological solution conditions displays large-scale structural heterogeneity, defining a new category of "metamorphic proteins". Previous Ltn studies have indicated that each form is required for proper function, but the mechanism of interconversion remains unknown. Here we have investigated the temperature dependence of kinetic rates associated with interconversion and unfolding by stopped-flow fluorescence to determine transition-state free energies. Comparisons of derived thermodynamic parameters revealed striking similarities between interconversion and protein unfolding. We conclude that Ltn native-state rearrangement proceeds by way of a large-scale unfolding process rather than a unique intermediate structure.     

Diffusional barrier in the unfolding of a small protein

To determine how the dynamics of the polypeptide chain in a protein molecule are coupled to the bulk solvent viscosity, the unfolding by urea of the small protein barstar was studied in the presence of two viscogens, xylose and glycerol. Thermodynamic studies of unfolding show that both viscogens stabilize barstar by a preferential hydration mechanism, and that viscogen and urea act independently on protein stability. Kinetic studies of unfolding show that while the rate-limiting conformational change during unfolding is dependent on the bulk solvent viscosity, eta, its rate does not show an inverse dependence on eta, as expected by Kramers' theory. Instead, the rate is found to be inversely proportional to an effective viscosity, eta + xi, where xi is an adjustable parameter which needs to be included in the rate equation. xi is found to have a value of -0.7 cP in xylose and -0.5 cP in glycerol, in the case of unfolding, at constant urea concentration as well as under isostability conditions. Hence, the unfolding protein chain does not experience the bulk solvent viscosity, but instead an effective solvent viscosity, which is lower than the bulk solvent viscosity by either 0.7 cP or 0.5 cP. A second important result is the validation of the isostability assumption, commonly used in protein folding studies but hitherto untested, according to which if a certain concentration of urea can nullify the effect of a certain concentration of viscogen on stability, then the same concentrations of urea and viscogen will also not perturb the free energy of activation of the unfolding of the protein.     

Increasing temperature accelerates protein unfolding without changing the pathway of unfolding

We have traditionally relied on extremely elevated temperatures (498K, 225 degrees C) to investigate the unfolding process of proteins within the timescale available to molecular dynamics simulations with explicit solvent. However, recent advances in computer hardware have allowed us to extend our thermal denaturation studies to much lower temperatures. Here we describe the results of simulations of chymotrypsin inhibitor 2 at seven temperatures, ranging from 298K to 498K. The simulation lengths vary from 94ns to 20ns, for a total simulation time of 344ns, or 0.34 micros. At 298K, the protein is very stable over the full 50ns simulation. At 348K, corresponding to the experimentally observed melting temperature of CI2, the protein unfolds over the first 25ns, explores partially unfolded conformations for 20ns, and then refolds over the last 35ns. Above its melting temperature, complete thermal denaturation occurs in an activated process. Early unfolding is characterized by sliding or breathing motions in the protein core, leading to an unfolding transition state with a weakened core and some loss of secondary structure. After the unfolding transition, the core contacts are rapidly lost as the protein passes on to the fully denatured ensemble. While the overall character and order of events in the unfolding process are well conserved across temperatures, there are substantial differences in the timescales over which these events take place. We conclude that 498K simulations are suitable for elucidating the details of protein unfolding at a minimum of computational expense.     

A measure of conformational entropy change during thermal protein unfolding using neutron spectroscopy

Thermal unfolding of proteins at high temperatures is caused by a strong increase of the entropy change which lowers Gibbs free energy change of the unfolding transition (DeltaG(unf) = DeltaH - TDeltaS). The main contributions to entropy are the conformational entropy of the polypeptide chain itself and ordering of water molecules around hydrophobic side chains of the protein. To elucidate the role of conformational entropy upon thermal unfolding in more detail, conformational dynamics in the time regime of picoseconds was investigated with neutron spectroscopy. Confined internal structural fluctuations were analyzed for alpha-amylase in the folded and the unfolded state as a function of temperature. A strong difference in structural fluctuations between the folded and the unfolded state was observed at 30 degrees C, which increased even more with rising temperatures. A simple analytical model was used to quantify the differences of the conformational space explored by the observed protein dynamics for the folded and unfolded state. Conformational entropy changes, calculated on the basis of the applied model, show a significant increase upon heating. In contrast to indirect estimates, which proposed a temperature independent conformational entropy change, the measurements presented here, demonstrated that the conformational entropy change increases with rising temperature and therefore contributes to thermal unfolding.     

NMR unfolding studies on a liver bile acid binding protein reveal a global two-state unfolding and localized singular behaviors

The folding properties of a bile acid binding protein, belonging to a subfamily of the fatty acid binding proteins, have been here investigated both by hydrogen exchange measurements, using the SOFAST NMR approach, and urea denaturation experiments. The urea unfolding profiles of individual residues, acting as single probes, were simultaneously analyzed through a global fit, according to a two-state unfolding model. The resulting conformational stability ΔG_U(H₂O) = 7.2 ± 0.25 kcal mol⁻¹ is in good agreement with hydrogen exchange stability ΔG_op. While the majority of protein residues satisfy this model, few amino-acids display a singular behavior, not directly amenable to the presence of a folding intermediate, as reported for other fatty acid binding proteins. These residues are part of a protein patch characterized by enhanced plasticity. To explain this singular behavior a tentative model has been proposed which takes into account the interplay between the dynamic features and the formation of transient aggregates. A functional role for this plasticity, related to translocation across the nuclear membrane, is discussed.     

Topological determinants of protein unfolding rates

For proteins that fold by two-state kinetics, the folding and unfolding processes are believed to be closely related to their native structures. In particular, folding and unfolding rates are influenced by the native structures of proteins. Thus, we focus on finding important topological quantities from a protein structure that determine its unfolding rate. After constructing graphs from protein native structures, we investigate the relationships between unfolding rates and various topological quantities of the graphs. First, we find that the correlation between the unfolding rate and the contact order is not as prominent as in the case of the folding rate and the contact order. Next, we investigate the correlation between the unfolding rate and the clustering coefficient of the graph of a protein native structure, and observe no correlation between them. Finally, we find that a newly introduced quantity, the impact of edge removal per residue, has a good overall correlation with protein unfolding rates. The impact of edge removal is defined as the ratio of the change of the average path length to the edge removal probability. From these facts, we conclude that the protein unfolding process is closely related to the protein native structure.     

Monitoring protein unfolding transitions by NMR-spectroscopy

NMR-spectroscopy has certain unique advantages for recording unfolding transitions of proteins compared e.g. to optical methods. It enables per-residue monitoring and separate detection of the folded and unfolded state as well as possible equilibrium intermediates. This allows a detailed view on the state and cooperativity of folding of the protein of interest and the correct interpretation of subsequent experiments. Here we summarize in detail practical and theoretical aspects of such experiments. Certain pitfalls can be avoided, and meaningful simplification can be made during the analysis. Especially a good understanding of the NMR exchange regime and relaxation properties of the system of interest is beneficial. We show by a global analysis of signals of the folded and unfolded state of GB1 how accurate values of unfolding can be extracted and what limits different NMR detection and unfolding methods. E.g. commonly used exchangeable amides can lead to a systematic under determination of the thermodynamic protein stability. We give several perspectives of how to deal with more complex proteins and how the knowledge about protein stability at residue resolution helps to understand protein properties under crowding conditions, during phase separation and under high pressure.

     

Machinery of protein folding and unfolding

During the past two years, a large amount of biochemical, biophysical and low- to high-resolution structural data have provided mechanistic insights into the machinery of protein folding and unfolding. It has emerged that dual functionality in terms of folding and unfolding might exist for some systems. The majority of folding/unfolding machines adopt oligomeric ring structures in a cooperative fashion and utilise the conformational changes induced by ATP binding/hydrolysis for their specific functions.     

LDL protein nitration: implication for LDL protein unfolding

Oxidatively- or enzymatically-modified low-density lipoprotein (LDL) is intimately involved in the initiation and progression of atherosclerosis. The in vivo modified LDL is electro-negative (LDL⁻) and consists of peroxidized lipid and unfolded apoB-100 protein. This study was aimed at establishing specific protein modifications and conformational changes in LDL⁻ assessed by liquid chromatography/tandem mass spectrometry (LC/MS/MS) and circular dichroism analyses, respectively. The functional significance of these chemical modifications and structural changes were validated with binding and uptake experiments to- and by bovine aortic endothelial cells (BAEC).The plasma LDL⁻ fraction showed increased nitrotyrosine and lipid peroxide content as well as a greater cysteine oxidation as compared with native- and total-LDL. LC/MS/MS analyses of LDL⁻ revealed specific modifications in the apoB-100 moiety, largely involving nitration of tyrosines in the α-helical structures and β₂ sheet as well as cysteine oxidation to cysteic acid in β₁ sheet. Circular dichroism analyses showed that the α-helical content of LDL⁻ was substantially lower (∼25%) than that of native LDL (∼90%); conversely, LDL⁻ showed greater content of β-sheet and random coil structure, in agreement with unfolding of the protein. These results were mimicked by treatment of LDL subfractions with peroxynitrite (ONOO⁻) or SIN-1: similar amino acid modifications as well as conformational changes (loss of α-helical structure and gain in β-sheet structure) were observed. Both LDL⁻ and ONOO⁻-treated LDL showed a statistically significant increase in binding and uptake to- and by BAEC compared to native LDL. We further found that most binding and uptake in control-LDL was through LDL-R with minimal oxLDL-R-dependent uptake. ONOO⁻-treated LDL was significantly bound and endocytosed by LOX-1, CD36, and SR-A with minimal contribution from LDL-R.It is suggested that lipid peroxidation and protein nitration may account for the mechanisms leading to apoB-100 protein unfolding and consequential increase in modified LDL binding and uptake to and by endothelial cells that is dependent on oxLDL scavenger receptors.     

Viscoelastic study of the mechanical unfolding of a protein by AFM

We have applied a dynamic force modulation technique to the mechanical unfolding of a homopolymer of immunoglobulin (Ig) domains from titin, (C47S C63S I27)5, [(I27)5] to determine the viscoelastic response of single protein molecules as a function of extension. Both the stiffness and the friction of the homopolymer system show a sudden decrease when a protein domain unfolds. The decrease in measured friction suggests that the system is dominated by the internal friction of the (I27)5 molecule and not solvent friction. In the stiffness-extension spectrum we detected an abrupt feature before each unfolding event, the amplitude of which decreased with each consecutive unfolding event. We propose that these features are a clear indication of the formation of the known unfolding intermediate of I27, which has been observed previously in constant velocity unfolding experiments. This simple force modulation AFM technique promises to be a very useful addition to constant velocity experiments providing detailed viscoelastic characterization of single molecules under extension.     

Probing membrane protein unfolding with pulse proteolysis

Technical challenges have greatly impeded the investigation of membrane protein folding and unfolding. To develop a new tool that facilitates the study of membrane proteins, we tested pulse proteolysis as a probe for membrane protein unfolding. Pulse proteolysis is a method to monitor protein folding and unfolding, which exploits the significant difference in proteolytic susceptibility between folded and unfolded proteins. This method requires only a small amount of protein and, in many cases, may be used with unpurified proteins in cell lysates. To evaluate the effectiveness of pulse proteolysis as a probe for membrane protein unfolding, we chose Halobacterium halobium bacteriorhodopsin (bR) as a model system. The denaturation of bR in SDS has been investigated extensively by monitoring the change in the absorbance at 560 nm (A₅₆₀). In this work, we demonstrate that denaturation of bR by SDS results in a significant increase in its susceptibility to proteolysis by subtilisin. When pulse proteolysis was applied to bR incubated in varying concentrations of SDS, the remaining intact protein determined by electrophoresis shows a cooperative transition. The midpoint of the cooperative transition (C_m) shows excellent agreement with that determined by A₅₆₀. The C_m values determined by pulse proteolysis for M56A and Y57A bRs are also consistent with the measurements made by A₅₆₀. Our results suggest that pulse proteolysis is a quantitative tool to probe membrane protein unfolding. Combining pulse proteolysis with Western blotting may allow the investigation of membrane protein unfolding in situ without overexpression or purification.