The Determinants of Stability and Folding in Evolutionarily Diverged Cytochromes c

Cytochrome c has served as a paradigm for the study of protein stability, folding, and molecular evolution, but it remains unclear how these aspects of the protein are related. For example, while the bovine and equine cytochromes c are known to have different stabilities, and possibly different folding mechanisms, it is not known how these differences arise from just three amino acid substitutions introduced during divergence. Using site-selectively incorporated carbon-deuterium bonds, we show that like the equine protein, bovine cytochrome c is induced to unfold by guanidine hydrochloride via a stepwise mechanism, but it does not populate an intermediate as is observed with the equine protein. The increased stability also results in more similar free energies of unfolding observed at different sites within the protein, giving the appearance of a more concerted mechanism. Furthermore, we show that the differences in stability and folding appear to result from a single amino acid substitution that stabilizes a helix by allowing for increased solvation of its N-terminus.     

A logical OR redundancy within the Asx-Pro-Asx-Gly type I beta-turn motif

Turn secondary structure is essential to the formation of globular protein architecture. Turn structures are, however, much more complex than either α-helix or β-sheet, and the thermodynamics and folding kinetics are poorly understood. Type I β-turns are the most common type of reverse turn, and they exhibit a statistical consensus sequence of Asx-Pro-Asx-Gly (where Asx is Asp or Asn). A comprehensive series of individual and combined Asx mutations has been constructed within three separate type I 3:5 G1 bulge β-turns in human fibroblast growth factor-1, and their effects on structure, stability, and folding have been determined. The results show a fundamental logical OR relationship between the Asx residues in the motif, involving H-bond interactions with main-chain amides within the turn. These interactions can be modulated by additional interactions with residues adjacent to the turn at positions i + 4 and i + 6. The results show that the Asx residues in the turn motif make a substantial contribution to the overall stability of the protein, and the Asx logical OR relationship defines a redundant system that can compensate for deleterious point mutations. The results also show that the stability of the turn is unlikely to be the prime determinant of formation of turn structure in the folding transition state.     

Equilibrium folding of pro-HlyA from Escherichia coli reveals a stable calcium ion dependent folding intermediate

HlyA from Escherichia coli is a member of the repeats in toxin (RTX) protein family, produced by a wide range of Gram-negative bacteria and secreted by a dedicated Type 1 Secretion System (T1SS). RTX proteins are thought to be secreted in an unfolded conformation and to fold upon secretion by Ca^2 + binding. However, the exact mechanism of secretion, ion binding and folding to the correct native state remains largely unknown. In this study we provide an easy protocol for high-level pro-HlyA purification from E. coli. Equilibrium folding studies, using intrinsic tryptophan fluorescence, revealed the well-known fact that Ca^2 + is essential for stability as well as correct folding of the whole protein. In the absence of Ca^2 +, pro-HlyA adopts a non-native conformation. Such molecules could however be rescued by Ca^2 + addition, indicating that these are not dead-end species and that Ca^2 + drives pro-HlyA folding. More importantly, pro-HlyA unfolded via a two-state mechanism, whereas folding was a three-state process. The latter is indicative of the presence of a stable folding intermediate. Analysis of deletion and Trp mutants revealed that the first folding transition, at 6–7 M urea, relates to Ca^2 + dependent structural changes at the extreme C-terminus of pro-HlyA, sensed exclusively by Trp914. Since all Trp residues of HlyA are located outside the RTX domain, our results demonstrate that Ca^2 + induced folding is not restricted to the RTX domain. Taken together, Ca^2 + binding to the pro-HlyA RTX domain is required to drive the folding of the entire protein to its native conformation.     

Interplay between Chaperones and Protein Disorder Promotes the Evolution of Protein Networks

Evolution is driven by mutations, which lead to new protein functions but come at a cost to protein stability. Non-conservative substitutions are of interest in this regard because they may most profoundly affect both function and stability. Accordingly, organisms must balance the benefit of accepting advantageous substitutions with the possible cost of deleterious effects on protein folding and stability. We here examine factors that systematically promote non-conservative mutations at the proteome level. Intrinsically disordered regions in proteins play pivotal roles in protein interactions, but many questions regarding their evolution remain unanswered. Similarly, whether and how molecular chaperones, which have been shown to buffer destabilizing mutations in individual proteins, generally provide robustness during proteome evolution remains unclear. To this end, we introduce an evolutionary parameter λ that directly estimates the rate of non-conservative substitutions. Our analysis of λ in Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens sequences reveals how co- and post-translationally acting chaperones differentially promote non-conservative substitutions in their substrates, likely through buffering of their destabilizing effects. We further find that λ serves well to quantify the evolution of intrinsically disordered proteins even though the unstructured, thus generally variable regions in proteins are often flanked by very conserved sequences. Crucially, we show that both intrinsically disordered proteins and highly re-wired proteins in protein interaction networks, which have evolved new interactions and functions, exhibit a higher λ at the expense of enhanced chaperone assistance. Our findings thus highlight an intricate interplay of molecular chaperones and protein disorder in the evolvability of protein networks. Our results illuminate the role of chaperones in enabling protein evolution, and underline the importance of the cellular context and integrated approaches for understanding proteome evolution. We feel that the development of λ may be a valuable addition to the toolbox applied to understand the molecular basis of evolution.     

Structural Bases for Stability-Function Tradeoffs in Antibiotic Resistance

Preorganization of enzyme active sites for substrate recognition typically comes at a cost to the stability of the folded form of the protein; consequently, enzymes can be dramatically stabilized by substitutions that attenuate the size and preorganization “strain” of the active site. How this stability-activity tradeoff constrains enzyme evolution has remained less certain, and it is unclear whether one should expect major stability insults as enzymes mutate towards new activities or how these new activities manifest structurally. These questions are both germane and easy to study in β-lactamases, which are evolving on the timescale of years to confer resistance to an ever-broader spectrum of β-lactam antibiotics. To explore whether stability is a substantial constraint on this antibiotic resistance evolution, we investigated extended-spectrum mutants of class C β-lactamases, which had evolved new activity versus third-generation cephalosporins. Five mutant enzymes had between 100-fold and 200-fold increased activity against the antibiotic cefotaxime in enzyme assays, and the mutant enzymes all lost thermodynamic stability (from 1.7 kcal mol^− 1 to 4.1 kcal mol^− 1), consistent with the stability-function hypothesis. Intriguingly, several of the substitutions were 10-20 Å from the catalytic serine; the question of how they conferred extended-spectrum activity arose. Eight structures, including complexes with inhibitors and extended-spectrum antibiotics, were determined by X-ray crystallography. Distinct mechanisms of action, including changes in the flexibility and ground-state structures of the enzyme, are revealed for each mutant. These results explain the structural bases for the antibiotic resistance conferred by these substitutions and their corresponding decrease in protein stability, which will constrain the evolution of new antibiotic resistance.     

Thermodynamic and kinetic characterization of apoHmpH, a fast-folding bacterial globin

Despite the widespread presence of the globin fold in most living organisms, only eukaryotic globins have been employed as model proteins in folding/stability studies so far. This work introduces the first thermodynamic and kinetic characterization of a prokaryotic globin, that is, the apo form of the heme-binding domain of flavohemoglobin (apoHmpH) from Escherichia coli. This bacterial globin has a widely different sequence but nearly identical structure to its eukaryotic analogues. We show that apoHmpH is a well-folded monomeric protein with moderate stability at room temperature [apparent ΔG°_UN(w) = − 3.1 ± 0.3 kcal mol^− 1; m_UN = − 1.7 kcal mol^− 1 M^− 1] and predominant α-helical structure. Remarkably, apoHmpH is the fastest-folding globin known to date, as it refolds about 4- to 16-fold more rapidly than its eukaryotic analogues (e.g., sperm whale apomyoglobin and soybean apoleghemoglobin), populating a compact kinetic intermediate (β_I = 0.9 ± 0.2) with significant helical content. Additionally, the single Trp120 (located in the native H helix) becomes locked into a fully native-like environment within 6 ms, suggesting that this residue and its closest spatial neighbors complete their folding at ultrafast (submillisecond) speed. In summary, apoHmpH is a bacterial globin that shares the general folding scheme (i.e., a rapid burst phase followed by slower rate-determining phases) of its eukaryotic analogues but displays an overall faster folding and a kinetic intermediate with some fully native-like traits. This study supports the view that the general folding features of bacterial and eukaryotic globins are preserved through evolution while kinetic details differ.     

Molecular genetics of retinal degeneration: A Drosophila perspective

Inherited retinal degeneration in Drosophila has been explored for insights into similar processes in humans. Based on the mechanisms, I divide these mutations in Drosophila into three classes. The first consists of genes that control the specialization of photoreceptor cells including the morphogenesis of visual organelles (rhabdomeres) that house the visual signaling proteins. The second class contains genes that regulate the activity or level of the major rhodopsin, Rh1, which is the light sensor and also provides a structural role for the maintenance of rhabdomeres. Some mutations in Rh1 (NinaE) are dominant due to constitutive activity or folding defects, like autosomal dominant retinitis pigmentosa (ADRP) in humans. The third class consists of genes that control the Ca²⁺ influx directly or indirectly by promoting the turnover of the second messenger and regeneration of PIP₂, or mediate the Ca²⁺-dependent regulation of the visual response. These gene products are critical for the increase in cytosolic Ca²⁺ following light stimulation to initiate negative regulatory events. Here I will focus on the signaling mechanisms underlying the degeneration in norpA, and in ADRP-type NinaE mutants that produce misfolded Rh1. Accumulation of misfolded Rh1 in the ER triggers the unfolded protein response (UPR), while endosomal accumulation of activated Rh1 may initiate autophagy in norpA. Both autophagy and the UPR are beneficial for relieving defective endosomal trafficking and the ER stress, respectively. However, when photoreceptors fail to cope with the persistence of these stresses, a cell death program is activated leading to retinal degeneration.     

Multiple tryptophan probes reveal that ubiquitin folds via a late misfolded intermediate

Much of our understanding of protein folding mechanisms is derived from experiments using intrinsic fluorescence of natural or genetically inserted tryptophan (Trp) residues to monitor protein refolding and site-directed mutagenesis to determine the energetic role of amino acids in the native (N), intermediate (I) or transition (T) states. However, this strategy has limited use to study complex folding reactions because a single fluorescence probe may not detect all low-energy folding intermediates. To overcome this limitation, we suggest that protein refolding should be monitored with different solvent-exposed Trp probes. Here, we demonstrate the utility of this approach by investigating the controversial folding mechanism of ubiquitin (Ub) using Trp probes located at residue positions 1, 28, 45, 57, and 66. We first show that these Trp are structurally sensitive and minimally perturbing fluorescent probes for monitoring folding/unfolding of the protein. Using a conventional stopped-flow instrument, we show that ANS and Trp fluorescence detect two distinct transitions during the refolding of all five Trp mutants at low concentrations of denaturant: T₁, a denaturant-dependent transition and T₂, a slower transition, largely denaturant-independent. Surprisingly, some Trp mutants (Ub^M1W, Ub^S57W) display Trp fluorescence changes during T₁ that are distinct from the expected U → N transition suggesting that the denaturant-dependent refolding transition of Ub is not a U → N transition but represents the formation of a structurally distinct I-state (U → I). Alternatively, this U → I transition could be also clearly distinguished by using a combination of two Trp mutations Ub^F45W-T66W for which the two Trp probes that display fluorescence changes of opposite sign during T₁ and T₂ (Ub^F45W-T66W). Global fitting of the folding/unfolding kinetic parameters and additional folding-unfolding double-jump experiments performed on Ub^M1W, a mutant with enhanced fluorescence in the I-state, demonstrate that the I-state is stable, compact, misfolded, and on-pathway. These results illustrate how transient low-energy I-states can be characterized efficiently in complex refolding reactions using multiple Trp probes.     

Overcoming hysteresis to attain reversible equilibrium folding for outer membrane phospholipase A in phospholipid bilayers

The free energy of unfolding of a membrane protein from lipids into water (ΔG^o_w,l) describes its equilibrium thermodynamic stability. Knowing this parameter gives insight into a membrane protein's sequence-structure-energy relationships. However, there are few measures of membrane protein stability because of the technical difficulties associated with unfolded and partially folded states. Here, we describe the experimental process that allowed us to measure the ΔG^o_w,l of the outer membrane phospholipase A into large unilamellar vesicles (LUVs) of 1,2-dilauroyl-sn-glycero-3-phosphocholine. To arrive at this reversible folding condition, we screened a large number of experimental variables: temperature, incubation time, salt concentration, pH, lipid composition and liposome morphology. The principal challenge we encountered under most conditions was hysteresis between folding and unfolding titrations. A second factor that compromised reversible folding was the observation that a fraction of the protein population tended to aggregate. We found that hysteresis could be completely eliminated on a feasible timescale by conducting experiments at acidic pH, by the slow dilution of the protein in the initial titration setup and by utilizing a low concentration of a detergent as a temporary “holdase” to solubilize the protein upon its initial dilution into folding conditions. We confirmed that the detergent did not disrupt the LUVs using fluorescence emission of lipid-sensitive dyes and light scattering. The results of our parameter search should be generally useful for efforts to measure ΔG^o_w,l for other membrane proteins.     

Folding and unfolding mechanism of highly stable full-consensus ankyrin repeat proteins

Full-consensus designed ankyrin repeat proteins were designed with one to six identical repeats flanked by capping repeats. These proteins express well in Escherichia coli as soluble monomers. Compared to our previously described designed ankyrin repeat protein library, randomized positions have now been fixed according to sequence statistics and structural considerations. Their stability increases with length and is even higher than that of library members, and those with more than three internal repeats are resistant to denaturation by boiling or guanidine hydrochloride. Full denaturation requires their heating in 5 M guanidine hydrochloride. The folding and unfolding kinetics of the proteins with up to three internal repeats were analyzed, as the other proteins could not be denatured. Folding is monophasic, with a rate that is nearly identical for all proteins (∼ 400-800 s^− 1), indicating that essentially the same transition state must be crossed, possibly the folding of a single repeat. In contrast, the unfolding rate decreases by a factor of about 10⁴ with increasing repeat number, directly reflecting thermodynamic stability in these extraordinarily slow denaturation rates. The number of unfolding phases also increases with repeat number. We analyzed the folding thermodynamics and kinetics both by classical two-state and three-state cooperative models and by an Ising-like model, where repeats are considered as two-state folding units that can be stabilized by interacting with their folded nearest neighbors. This Ising model globally describes both equilibrium and kinetic data very well and allows for a detailed explanation of the ankyrin repeat protein folding mechanism.     

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.     

Structural Determinants in Prion Protein Folding and Stability

Prions are responsible for a heterogeneous group of fatal neurodegenerative diseases, involving post-translational modifications of the cellular prion protein. Epidemiological studies on Creutzfeldt-Jakob disease, a prototype prion disorder, show a majority of cases being sporadic, while the remaining occurrences are either genetic or iatrogenic. The molecular mechanisms by which PrP^C is converted into its pathological isoform have not yet been established. While point mutations and seeds trigger the protein to cross the energy barriers, thus causing genetic and infectious transmissible spongiform encephalopathies, respectively, the mechanism responsible for sporadic forms remains unclear. Since prion diseases are protein-misfolding disorders, we investigated prion protein folding and stability as functions of different milieus. Using spectroscopic techniques and atomistic simulations, we dissected the contribution of major structural determinants, also defining the energy landscape of prion protein. In particular, we elucidated (i) the essential role of the octapeptide region in prion protein folding and stability, (ii) the presence of a very enthalpically stable intermediate in prion-susceptible species, and (iii) the role of the disulfide bridge in prion protein folding.     

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.     

Elimination of a cis-Proline-Containing Loop and Turn Optimization Stabilizes a Protein and Accelerates Its Folding

In the N2 domain of the gene-3-protein of phage fd, two consecutive β-strands are connected by a mobile loop of seven residues (157-163). The stability of this loop is low, and the Asp160-Pro161 bond at its tip shows conformational heterogeneity with 90% being in the cis and 10% in the trans form. The refolding kinetics of N2 are complex because the molecules with cis or trans isomers at Pro161 both fold to native-like conformations, albeit with different rates. We employed consensus design to shorten the seven-residue irregular loop around Pro161 to a four-residue type I′ turn without a proline. This increased the conformational stability of N2 by almost 10 kJ mol^− 1 and abolished the complexity of the folding kinetics. Turn sequences obtained from in vitro selections for increased stability strongly resembled those derived from the consensus design. Two other type I′ turns of N2 could also be stabilized by consensus design. For all three turns, the gain in stability originates from an increase in the rate of refolding. The turns form native-like structures early during refolding and thus stabilize the folding transition state. The crystal structure of the variant with all three stabilized turns confirms that the 157-163 loop was in fact shortened to a type I′ turn and that the other turns maintained their type I′ conformation after sequence optimization.     

Thermostable Bacillus subtilis lipases: in vitro evolution and structural insight

In vitro evolution methods are now being routinely used to identify protein variants with novel and enhanced properties that are difficult to achieve using rational design. However, one of the limitations is in screening for beneficial mutants through several generations due to the occurrence of neutral/negative mutations occurring in the background of positive ones. While evolving a lipase in vitro from mesophilic Bacillus subtilis to generate thermostable variants, we have designed protocols that combine stringent three-tier testing, sequencing and stability assessments on the protein at the end of each generation. This strategy resulted in a total of six stabilizing mutations in just two generations with three mutations per generation. Each of the six mutants when evaluated individually contributed additively to thermostability. A combination of all of them resulted in the best variant that shows a remarkable 15 °C shift in melting temperature and a millionfold decrease in the thermal inactivation rate with only a marginal increase of 3 kcal mol⁻¹ in free energy of stabilization. Notably, in addition to the dramatic shift in optimum temperature by 20 °C, the activity has increased two- to fivefold in the temperature range 25-65 °C. High-resolution crystal structures of three of the mutants, each with 5° increments in melting temperature, reveal the structural basis of these mutations in attaining higher thermostability. The structures highlight the importance of water-mediated ionic networks on the protein surface in imparting thermostability. Saturation mutagenesis at each of the six positions did not result in enhanced thermostability in almost all the cases, confirming the crucial role played by each mutation as revealed through the structural study. Overall, our study presents an efficient strategy that can be employed in directed evolution approaches employed for obtaining improved properties of proteins.     

The folding process of acylphosphatase from Escherichia coli is remarkably accelerated by the presence of a disulfide bond

The acylphosphatase from Escherichia coli (EcoAcP) is the first AcP so far studied with a disulfide bond. A mutational variant of the enzyme lacking the disulfide bond has been produced by substituting the two cysteine residues with alanine (EcoAcP mutational variant C5A/C49A, mutEcoAcP). The native states of the two protein variants are similar, as shown by far-UV and near-UV circular dichroism and dynamic light-scattering measurements. From unfolding experiments at equilibrium using intrinsic fluorescence and far-UV circular dichroism as probes, EcoAcP shows an increased conformational stability as compared with mutEcoAcP. The wild-type protein folds according to a two-state model with a very fast rate constant (k_F^H2O = 72,600 s^− 1), while mutEcoAcP folds ca 1500-fold slower, via the accumulation of a partially folded species. The correlation between the hydrophobicity of the polypeptide chain and the folding rate, found previously in the AcP-like structural family, is maintained only when considering the mutant but not the wild-type protein, which folds much faster than expected from this correlation. Similarly, the correlation between the relative contact order and the folding rate holds only for mutEcoAcP. The correlation also holds for EcoAcP, provided the relative contact order value is recalculated by considering the disulfide bridge as an alternate path for the backbone to determine the shortest sequence separation between contacting residues. These results indicate that the presence of a disulfide bond in a protein is an important determinant of the folding rate and allows its contribution to be determined in quantitative terms.