An approach for the assessment of the order of disruption of the elements of protein structure upon protein unfolding: A study of carbonic anhydrase B

An experimental approach named μ-analysis has been developed in order to elucidate the sequence of the loss of ordered structure by elements of a protein during the denaturation of the molecule. This approach is applicable for the analysis of proteins that fold (unfold) in a multistep process that involve the formation (destruction) of a range of intermediate states. The concept of the approach consists in systematic analysis of mutagenized forms of the protein with point substitutions of hydrophobic amino-acid residues and additional cysteine bridges. Importantly, the substitutions of the amino-acid residues must be localized to the same structural elements of the protein. Point substitutions of hydrophobic amino-acid residues mainly provide information on the structural elements of the protein that are disrupted at the final stages of protein denaturation. The addition of cysteine bridges to the surface of the protein molecule allows investigation of structural elements of the protein that are the first to unfold upon protein denaturation. Calorimetric studies of non-equilibrium melting of bovine carbonic anhydrase B yielded information on the rate constants of the unfolding of ten mutant forms of the protein. The analysis of the effects of mutations on the rates of different stages of protein unfolding allowed for elucidation of the order of disruption of structural elements of carbonic anhydrase B upon thermal denaturation.     

Evidence for a Shared Mechanism in the Formation of Urea-Induced Kinetic and Equilibrium Intermediates of Horse Apomyoglobin from Ultrarapid Mixing Experiments

In this study, the equivalence of the kinetic mechanisms of the formation of urea-induced kinetic folding intermediates and non-native equilibrium states was investigated in apomyoglobin. Despite having similar structural properties, equilibrium and kinetic intermediates accumulate under different conditions and via different mechanisms, and it remains unknown whether their formation involves shared or distinct kinetic mechanisms. To investigate the potential mechanisms of formation, the refolding and unfolding kinetics of horse apomyoglobin were measured by continuous- and stopped-flow fluorescence over a time range from approximately 100 μs to 10 s, along with equilibrium unfolding transitions, as a function of urea concentration at pH 6.0 and 8°C. The formation of a kinetic intermediate was observed over a wider range of urea concentrations (0–2.2 M) than the formation of the native state (0–1.6 M). Additionally, the kinetic intermediate remained populated as the predominant equilibrium state under conditions where the native and unfolded states were unstable (at ~0.7–2 M urea). A continuous shift from the kinetic to the equilibrium intermediate was observed as urea concentrations increased from 0 M to ~2 M, which indicates that these states share a common kinetic folding mechanism. This finding supports the conclusion that these intermediates are equivalent. Our results in turn suggest that the regions of the protein that resist denaturant perturbations form during the earlier stages of folding, which further supports the structural equivalence of transient and equilibrium intermediates. An additional folding intermediate accumulated within ~140 μs of refolding and an unfolding intermediate accumulated in <1 ms of unfolding. Finally, by using quantitative modeling, we showed that a five-state sequential scheme appropriately describes the folding mechanism of horse apomyoglobin.     

Decoding Structural Properties of a Partially Unfolded Protein Substrate: En Route to Chaperone Binding

Many proteins comprising of complex topologies require molecular chaperones to achieve their unique three-dimensional folded structure. The E.coli chaperone, GroEL binds with a large number of unfolded and partially folded proteins, to facilitate proper folding and prevent misfolding and aggregation. Although the major structural components of GroEL are well defined, scaffolds of the non-native substrates that determine chaperone-mediated folding have been difficult to recognize. Here we performed all-atomistic and replica-exchange molecular dynamics simulations to dissect non-native ensemble of an obligate GroEL folder, DapA. Thermodynamics analyses of unfolding simulations revealed populated intermediates with distinct structural characteristics. We found that surface exposed hydrophobic patches are significantly increased, primarily contributed from native and non-native β-sheet elements. We validate the structural properties of these conformers using experimental data, including circular dichroism (CD), 1-anilinonaphthalene-8-sulfonic acid (ANS) binding measurements and previously reported hydrogen-deutrium exchange coupled to mass spectrometry (HDX-MS). Further, we constructed network graphs to elucidate long-range intra-protein connectivity of native and intermediate topologies, demonstrating regions that serve as central “hubs”. Overall, our results implicate that genomic variations (or mutations) in the distinct regions of protein structures might disrupt these topological signatures disabling chaperone-mediated folding, leading to formation of aggregates.     

Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

Model of β-Sheet of Muscle Fatty Acid Binding Protein of Locusta migratoria Displays Characteristic Topology

The β-sheet of muscle fatty acid binding protein of Locusta migratoria (Lm-FABP) was modeled by employing 2-D NMR data and the Rigid Body Assembly method. The model shows the β-sheet to comprise ten β-strands arranged anti-parallel to each other. There is a β-bulge between Ser 13 and Gln 14 which is a difference from the published structure of β-sheet of bovine heart Fatty Acid Binding Protein. Also, a hydrophobic patch consisting of Ile 45, Phe 51, Phe 64 and Phe 66 is present on the surface which is characteristic of most Fatty Acid Binding Proteins. A “gap” is present between βD and βE that provides evidence for the presence of a portal or opening between the polypeptide chains which allows ligand fatty acids to enter the protein cavity and bind to the protein.     

Reversible thermal denaturation of a 60-kDa genetically engineered beta-sheet polypeptide

A de novo 687-amino-acid residue polypeptide with a regular 32-amino-acid repeat sequence, (GA)₃GY(GA)₃GE(GA)₃GH(GA)₃GK, forms large β-sheet assemblages that exhibit remarkable folding properties and, as well, form fibrillar structures. This construct is an excellent tool to explore the details of β-sheet formation yielding intimate folding information that is otherwise difficult to obtain and may inform folding studies of naturally occurring materials. The polypeptide assumes a fully folded antiparallel β-sheet/turn structure at room temperature, and yet is completely and reversibly denatured at 125°C, adopting a predominant polyproline II conformation. Deep ultraviolet Raman spectroscopy indicated that melting/refolding occurred without any spectroscopically distinct intermediates, yet the relaxation kinetics depend on the initial polypeptide state, as would be indicative of a non-two-state process. Thermal denaturation and refolding on cooling appeared to be monoexponential with characteristic times of ~1 and ~60 min, respectively, indicating no detectable formation of hairpin-type nuclei in the millisecond timescale that could be attributed to nonlocal “nonnative” interactions. The polypeptide folding dynamics agree with a general property of β-sheet proteins, i.e., initial collapse precedes secondary structure formation. The observed folding is much faster than expected for a protein of this size and could be attributed to a less frustrated free-energy landscape funnel for folding. The polypeptide sequence suggests an important balance between the absence of strong nonnative contacts (salt bridges or hydrophobic collapse) and limited repulsion of charged side chains.     

Aromatic cluster mutations produce focal modulations of β-sheet structure

Site-directed mutagenesis is a powerful tool for altering the structure and function of proteins in a focused manner. Here, we examined how a model β-sheet protein could be tuned by mutation of numerous surface-exposed residues to aromatic amino acids. We designed these aromatic side chain “clusters” at highly solvent-exposed positions in the flat, single-layer β-sheet of Borrelia outer surface protein A (OspA). This unusual β-sheet scaffold allows us to interrogate the effects of these mutations in the context of well-defined structure but in the absence of the strong scaffolding effects of globular protein architecture. We anticipated that the introduction of a cluster of aromatic amino acid residues on the β-sheet surface would result in large conformational changes and/or stabilization and thereby provide new means of controlling the properties of β-sheets. Surprisingly, X-ray crystal structures revealed that the introduction of aromatic clusters produced only subtle conformational changes in the OspA β-sheet. Additionally, despite burying a large degree of hydrophobic surface area, the aromatic cluster mutants were slightly less stable than the wild-type scaffold. These results thereby demonstrate that the introduction of aromatic cluster mutations can serve as a means for subtly modulating β-sheet conformation in protein design.     

Human Lutropin (hLH) and Choriogonadotropin (CG) Are Assembled by Different Pathways: A MODEL OF hLH ASSEMBLY

The glycoprotein hormones are all structurally related heterodimers consisting of an α-subunit and a ligand-specific β-subunit that confers their unique biological activity. Crystal structures showed how the β-subunit surrounds a part of the α-subunit, and we showed the existence of the two mechanisms responsible for that assembly. In human choriogonadotropin, the β-subunit is folded before the subunits dock, and the α-subunit becomes incorporated into the dimer by a mechanism we termed “threading,” passing between parts of the preassembled β-subunit. Here, we show that the human lutropin β-subunit is not folded completely prior to its interaction with the α-subunit and show that docking of the subunits enables the α-subunit to serve as a chaperone to the β-subunit. Based on data described here, we propose that the α-subunit facilitates formation of the human lutropin β-subunit by two mechanisms. First, the cystine knot of the α-subunit potentiates formation of the β-subunit cystine knot, and second, contacts between α-subunit loop 2 and a hydrophobic tail in the β-subunit facilitate formation of the seatbelt latch disulfide, which stabilizes the heterodimer. The primary influence of the α-subunit was seen when the hydrophobic tail was present or absent, but the secondary mechanism was required only when the hydrophobic tail of the β-subunit was present. During the evolution of human choriogonadotropin, neither of these α-subunit roles was necessary for folding of the β-subunit. The complex mechanism for lutropin assembly may be required to provide an additional control on its positive feedback function in vertebrate reproduction.     

Predicting the protein folding nucleus from sequences [correction of a sequence]

Understanding the mechanism of protein folding would allow prediction of the three-dimensional structure from sequence data alone. It has been shown that small proteins fold in a small number of kinetic steps and that significantly populated intermediate states exist for some of them. Studies of these intermediates have demonstrated the existence of specific interactions established during the initial stages of folding. Comparison of the amino acids participating in these specific and essential interactions and constituting the folding nucleus with conserved hydrophobic positions of a given fold shows a striking correspondence. This finding opens the perspective of predicting the folding nucleus knowing only a set of divergent sequences of a protein family.     

Structural Basis for Conformational Dynamics of GTP-bound Ras Protein

Ras family small GTPases assume two interconverting conformations, “inactive” state 1 and “active” state 2, in their GTP-bound forms. Here, to clarify the mechanism of state transition, we have carried out x-ray crystal structure analyses of a series of mutant H-Ras and M-Ras in complex with guanosine 5′-(β,γ-imido)triphosphate (GppNHp), representing various intermediate states of the transition. Crystallization of H-RasT35S-GppNHp enables us to solve the first complete tertiary structure of H-Ras state 1 possessing two surface pockets unseen in the state 2 or H-Ras-GDP structure. Moreover, determination of the two distinct crystal structures of H-RasT35S-GppNHp, showing prominent polysterism in the switch I and switch II regions, reveals a pivotal role of the guanine nucleotide-mediated interaction between the two switch regions and its rearrangement by a nucleotide positional change in the state 2 to state 1 transition. Furthermore, the ³¹P NMR spectra and crystal structures of the GppNHp-bound forms of M-Ras mutants, carrying various H-Ras-type amino acid substitutions, also reveal the existence of a surface pocket in state 1 and support a similar mechanism based on the nucleotide-mediated interaction and its rearrangement in the state 1 to state 2 transition. Intriguingly, the conformational changes accompanying the state transition mimic those that occurred upon GDP/GTP exchange, indicating a common mechanistic basis inherent in the high flexibility of the switch regions. Collectively, these results clarify the structural features distinguishing the two states and provide new insights into the molecular basis for the state transition of Ras protein.     

The Roles of Helix I and Strand 5A in the Folding,Function and Misfolding of α1-Antitrypsin

α₁-Antitrypsin, the archetypal member of the serpin superfamily, is a metastable protein prone to polymerization when exposed to stressors such as elevated temperature, low denaturant concentrations or through the presence of deleterious mutations which, in a physiological context, are often associated with disease. Experimental evidence suggests that α₁-Antitrypsin can polymerize via several alternative mechanisms in vitro. In these polymerization mechanisms different parts of the molecule are proposed to undergo conformational change. Both strand 5 and helix I are proposed to adopt different conformations when forming the various polymers, and possess a number of highly conserved residues however their role in the folding and misfolding of α₁-Antitrypsin has never been examined. We have therefore created a range of α₁Antitypsin variants in order to explore the role of these conserved residues in serpin folding, misfolding, stability and function. Our data suggest that key residues in helix I mediate efficient folding from the folding intermediate and residues in strand 5A ensure native state stability in order to prevent misfolding. Additionally, our data indicate that helix I is involved in the inhibitory process and that both structural elements undergo differing conformational rearrangements during unfolding and misfolding. These findings suggest that the ability of α₁-Antitrypsin to adopt different types of polymers under different denaturing conditions may be due to subtle conformational differences in the transiently populated structures adopted prior to the I and M* states.     

Fine structure analysis of a protein folding transition state; distinguishing between hydrophobic stabilization and specific packing

Developing a detailed understanding of the structure and energetics of protein folding transition states is a key step in describing the folding process. The phi-value analysis approach allows the energetic contribution of side-chains to be mapped out by comparing wild-type with individual mutants where conservative changes are introduced. Studies where multiple substitutions are made at individual sites are much rarer but are potentially very useful for understanding the contribution of each element of a side-chain to transition state formation, and for distinguishing the relative importance of specific packing versus hydrophobic interactions. We have made a series of conservative mutations at multiple buried sites in the N-terminal domain of L9 in order to assess the relative importance of specific side-chain packing versus less specific hydrophobic stabilization of the transition state. A total of 28 variants were prepared using both naturally occurring and non-naturally occurring amino acids at six sites. Analysis of the mutants by NMR and CD showed no perturbation of the structure. There is no correlation between changes in hydrophobicity and changes in stability. In contrast, there is excellent linear correlation between the hydrophobicity of a side-chain and the log of the folding rate, ln(k(f)). The correlation between ln(k(f)) and the change in hydrophobicity holds even for substitutions that change the shape and/or size of a side-chain significantly. For most sites, the correlation with the logarithm of the unfolding rate, ln(k(u)), is much worse. Mutants with more hydrophobic amino acid substitutions fold faster, and those with less hydrophobic amino acid substitutions fold slower. The results show that hydrophobic interactions amongst core residues are an important driving force for forming the transition state, and are more important than specific tight packing interactions. Finally, a number of substitutions lead to negative phi-values and the origin of these effects are described.     

Deciphering the Hidden Informational Content of Protein Sequences: FOLDABILITY OF PROINSULIN HINGES ON A FLEXIBLE ARM THAT IS DISPENSABLE IN THE MATURE HORMONE*

Protein sequences encode both structure and foldability. Whereas the interrelationship of sequence and structure has been extensively investigated, the origins of folding efficiency are enigmatic. We demonstrate that the folding of proinsulin requires a flexible N-terminal hydrophobic residue that is dispensable for the structure, activity, and stability of the mature hormone. This residue (Phe^B1 in placental mammals) is variably positioned within crystal structures and exhibits ¹H NMR motional narrowing in solution. Despite such flexibility, its deletion impaired insulin chain combination and led in cell culture to formation of non-native disulfide isomers with impaired secretion of the variant proinsulin. Cellular folding and secretion were maintained by hydrophobic substitutions at B1 but markedly perturbed by polar or charged side chains. We propose that, during folding, a hydrophobic side chain at B1 anchors transient long-range interactions by a flexible N-terminal arm (residues B1–B8) to mediate kinetic or thermodynamic partitioning among disulfide intermediates. Evidence for the overall contribution of the arm to folding was obtained by alanine scanning mutagenesis. Together, our findings demonstrate that efficient folding of proinsulin requires N-terminal sequences that are dispensable in the native state. Such arm-dependent folding can be abrogated by mutations associated with β-cell dysfunction and neonatal diabetes mellitus.     

Slow and Bimolecular Folding of a De Novo Designed Monomeric Protein DS119

De novo protein design offers a unique means to test and advance our understanding of how proteins fold. However, most current design methods are native structure eccentric and folding kinetics has rarely been considered in the design process. Here, we show that a de novo designed mini-protein DS119, which folds into a βαβ structure, exhibits unusually slow and concentration-dependent folding kinetics. For example, the folding time for 50 μM of DS119 was estimated to be ∼2 s. Stopped-flow fluorescence resonance energy transfer experiments further suggested that its folding was likely facilitated by a transient dimerization process. Taken together, these results highlight the need for consideration of the entire folding energy landscape in de novo protein design and provide evidence suggesting nonnative interactions can play a key role in protein folding.     

Reassessing the folding of the KIX domain: Evidence for a two‐state mechanism

The debate about the presence and role of intermediates in the folding of proteins has been a critical issue, especially for fast folders. One of the classical methodologies to identify such metastable species is the “burst-phase analysis,” whereby the observed signal amplitude from stopped-flow traces is determined as a function of denaturant concentration. However, a complication may arise when folding is sufficiently fast to jeopardize the reliability of the stopped-flow technique. In this study, we reassessed the folding of the KIX domain from cAMP Response Element-Binding (CREB)-binding protein, which has been proposed to involve the formation of an intermediate that accumulates in the dead time of the stopped flow. By using an in-house-built capillary continuous flow with a 50-μs dead time, we demonstrate that this intermediate is not present; the problem arose because of the instrumental limitation of the standard stopped flow to assess very fast refolding rate constants (e.g., ≥500 s⁻¹).     

Easy to Use Plastic Optical Fiber-Based Biosensor for Detection of Butanal

The final goal of this work is to achieve a selective detection of butanal by the realization of a simple, small-size and low cost experimental approach. To this end, a porcine odorant-binding protein was used in connection with surface plasmon resonance transduction in a plastic optical fiber tool for the selective detection of butanal by a competitive assay. This allows to reduce the cost and the size of the sensing device and it offers the possibility to design a “Lab-on-a-chip” platform. The obtained results showed that this system approach is able to selectively detect the presence of butanal in the concentration range from 20 μM to 1000 μM.     

Titration Properties and Thermodynamics of the Transition State for Folding: Comparison of Two-state and Multi-state Folding Pathways

CI2 folds and unfolds as a single cooperative unit by simple two-state kinetics, which enables the properties of the transition state to be measured from both the forward and backward rate constants. We have examined how the free energy of the transition state for the folding of chymotrypsin inhibitor 2 (CI2) changes with pH and temperature. In addition to the standard thermodynamic quantities, we have measured the overall acid-titration properties of the transition state and its heat capacity relative to both the denatured and native states. We were able to determine the latter by a method analogous to a well-established procedure for measuring the change in heat capacity for equilibrium unfolding: the enthalpy of activation of unfolding at different values of acid pH were plotted against the average temperature of each determination. Our results show that the transition state of CI2 has lost most of the electrostatic and van der Waals' interactions that are found in the native state, but it remains compact and this prevents water molecules from entering some parts of the hydrophobic core. The properties of the transition state of CI2 are then compared with the major folding transition state of the larger protein barnase, which folds by a multi-state mechanism, with the accumulation of a partly structured intermediate (D^physor I). CI2 folds from a largely unstructured denatured state under physiological conditionsviaa transition state which is compact but relatively uniformly unstructured, with tertiary and secondary structure being formed in parallel. We term this an expanded pathway. Conversely, barnase folds from a largely structured denatured state in which elements of structure are well formed through a transition state that has islands of folded elements of structure. We term this a compact pathway. These two pathways may correspond to the two extreme ends of a continuous spectrum of protein folding mechanisms. Although the properties of the two transition states are very different, the activation barrier for folding (D^phys→3 ) is very similar for both proteins.     

BCL::Score—Knowledge Based Energy Potentials for Ranking Protein Models Represented by Idealized Secondary Structure Elements

The topology of most experimentally determined protein domains is defined by the relative arrangement of secondary structure elements, i.e. α-helices and β-strands, which make up 50–70% of the sequence. Pairing of β-strands defines the topology of β-sheets. The packing of side chains between α-helices and β-sheets defines the majority of the protein core. Often, limited experimental datasets restrain the position of secondary structure elements while lacking detail with respect to loop or side chain conformation. At the same time the regular structure and reduced flexibility of secondary structure elements make these interactions more predictable when compared to flexible loops and side chains. To determine the topology of the protein in such settings, we introduce a tailored knowledge-based energy function that evaluates arrangement of secondary structure elements only. Based on the amino acid C_β atom coordinates within secondary structure elements, potentials for amino acid pair distance, amino acid environment, secondary structure element packing, β-strand pairing, loop length, radius of gyration, contact order and secondary structure prediction agreement are defined. Separate penalty functions exclude conformations with clashes between amino acids or secondary structure elements and loops that cannot be closed. Each individual term discriminates for native-like protein structures. The composite potential significantly enriches for native-like models in three different databases of 10,000–12,000 protein models in 80–94% of the cases. The corresponding application, “BCL::ScoreProtein,” is available at www.meilerlab.org.     

Hydrophobic core packing in the SH3 domain folding transition state.

How tightly packed is the hydrophobic core of a folding transition state structure? We have addressed this question by characterizing the effects on folding kinetics of > 40 substitutions of both large and small amino acids in the hydrophobic core of the Fyn SH3 domain. Our results show that residues at three positions, which we designate as the 'core folding nucleus', are tightly packed in the transition state, and substitutions at these positions cause the largest changes in the folding rate. The other six positions examined appear to be loosely packed; thus, substitutions at these positions with larger hydrophobic residues generally accelerate folding, presumably by increasing the rate of nonspecific hydrophobic collapse. Surprisingly, the folding rate can be greatly accelerated by residues that also significantly destabilize the native state structure. Furthermore, mutants with identical thermodynamic stability can differ by up to 55-fold in their folding rates. These results highlight the importance of hydrophobic core composition, as opposed to only topology, in determining the folding rate of a protein. They also provide a new explanation for the 'abnormal' phi-values observed in many protein folding kinetics studies.     

Role of Protein Stabilizers on the Conformation of the Unfolded State of Cytochrome c and Its Early Folding Kinetics: INVESTIGATION AT SINGLE MOLECULAR RESOLUTION*

An insight into the conformation and dynamics of unfolded and early intermediate states of a protein is essential to understand the mechanism of its aggregation and to design potent inhibitor molecules. Fluorescence correlation spectroscopy has been used to study the effects of several model protein stabilizers on the conformation of the unfolded state and early folding dynamics of tetramethyl rhodamine-labeled cytochrome c from Saccharomyces cerevisiae at single molecular resolution. Special attention has been given to arginine, which is a widely used stabilizer for improving refolding yield of different proteins. The value of the hydrodynamic radius (r_H) obtained by analyzing the intensity fluctuations of the diffusing molecules has been found to increase in a two-state manner as the protein is unfolded by urea. The results further show that the presence of arginine and other protein stabilizers favors a relatively structured conformation of the unfolded states (r_H of 29 Å) over an extended one (r_H of 40 Å), which forms in their absence. Also, the time constant of a kinetic component (τ_R) of about 30 μs has been observed by analyzing the correlation functions, which represents formation of a collapsed state. This time constant varies with urea concentration representing an inverted Chevron plot that shows a roll-over and behavior in the absence of arginine. To the best of our knowledge, this is one of the first applications of fluorescence correlation spectroscopy to study direct folding kinetics of a protein.