Positive selection dictates the choice between kinetic and thermodynamic protein folding and stability in subtilases

Subtilisin E (SbtE) is a member of the ubiquitous superfamily of serine proteases called subtilases and serves as a model for understanding propeptide-mediated protein folding mechanisms. Unlike most proteins that adopt thermodynamically stable conformations, the native state of SbtE is trapped into a kinetically stable conformation. While kinetic stability offers distinct functional advantages to the native state, the constraints that dictate the selection between kinetic and thermodynamic folding and stability remain unknown. Using highly conserved subtilases, we demonstrate that adaptive evolution of sequence dictates selection of folding pathways. Intracellular and extracellular serine proteases (ISPs and ESPs, respectively) constitute two subfamilies within the family of subtilases that have highly conserved sequences, structures, and catalytic activities. Our studies on the folding pathways of subtilisin E (SbtE), an ESP, and its homologue intracellular serine protease 1 (ISP1), an ISP, show that although topology, contact order, and hydrophobicity that drive protein folding reactions are conserved, ISP1 and SbtE fold through significantly different pathways and kinetics. While SbtE absolutely requires the propeptide to fold into a kinetically trapped conformer, ISP1 folds to a thermodynamically stable state more than 1 million times faster and independent of a propeptide. Furthermore, kinetics establish that ISP1 and SbtE fold through different intermediate states. An evolutionary analysis of folding constraints in subtilases suggests that observed differences in folding pathways may be mediated through positive selection of specific residues that map mostly onto the protein surface. Together, our results demonstrate that closely related subtilases can fold through distinct pathways and mechanisms, and suggest that fine sequence details can dictate the choice between kinetic and thermodynamic folding and stability.     

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.     

Probing the folding intermediate of Rd-apocyt b562 by protein engineering and infrared T-jump

Small proteins often fold in an apparent two-state manner with the absence of detectable early-folding intermediates. Recently, using native-state hydrogen exchange, intermediates that exist after the rate-limiting transition state have been identified for several proteins. However, little is known about the folding kinetics from these post-transition intermediates to their corresponding native states. Herein, we have used protein engineering and a laser-induced temperature-jump (T-jump) technique to investigate this issue and have applied it to Rd-apocyt b(562) , a four-helix bundle protein. Previously, it has been shown that Rd-apocyt b(562) folds via an on-pathway hidden intermediate, which has only the N-terminal helix unfolded. In the present study, a double mutation (V16G/I17A) in the N-terminal helix of Rd-apocyt b(562) was made to further increase the relative population of this intermediate state at high temperature by selectively destabilizing the native state. In the circular dichroism thermal melting experiment, this mutant showed apparent two-state folding behavior. However, in the T-jump experiment, two kinetic phases were observed. Therefore, these results are in agreement with the idea that a folding intermediate is populated on the folding pathway of Rd-apocyt b(562) . Moreover, it was found that the exponential growth rate of the native state from this intermediate state is roughly (25 microsec)(-1) at 65 degrees C.     

Searching for multiple folding pathways of a nearly symmetrical protein: temperature dependent phi-value analysis of the B domain of protein A

The B domain of protein A (BdpA) is a popular paradigm for simulating protein folding pathways. The discrepancies between so many simulations and subsequent experimental testing may be attributable to the protein being highly symmetrical: changing experimental conditions could perturb the subtle interplay between the effects of symmetry in the native structure and the effects of asymmetry from specific interactions in a given sequence. If the protein folds via multiple pathways, perturbations, such as temperature, denaturant concentration, and mutation, should change the flux of micro pathways, leading to changes in the bulk properties of the transition state. We tested this hypothesis by conducting a Phi-analysis of BdpA as a function of temperature from 25.0 degrees C to 60.0 degrees C. The Phi-values had no significant dependence on temperature and the values at 55.0 degrees C (denaturing conditions) are very similar to those at 25.0 degrees C (folding conditions), indicating the structure of the transition state does not significantly change although the experimental conditions are considerably altered. The results suggest that BdpA folds via a single dominant folding pathway.     

The folding energy landscape of the dimerization domain of Escherichia coli Trp repressor: a joint experimental and theoretical investigation

Enhanced structural insights into the folding energy landscape of the N-terminal dimerization domain of Escherichia coli tryptophan repressor, [2-66]2 TR, were obtained from a combined experimental and theoretical analysis of its equilibrium folding reaction. Previous studies have shown that the three intertwined helices in [2-66]2 TR are sufficient to drive the formation of a stable dimer for the full-length protein, [2-107]2 TR. The monomeric and dimeric folding intermediates that appear during the folding reactions of [2-66]2 TR have counterparts in the folding mechanism of the full-length protein. The equilibrium unfolding energy surface on which the folding and dimerization reactions occur for [2-66]2 TR was examined with a combination of native-state hydrogen exchange analysis, pepsin digestion and matrix-assisted laser/desorption mass spectrometry performed at several concentrations of protein and denaturant. Peptides corresponding to all three helices in [2-66]2 TR show multi-layered protection patterns consistent with the relative stabilities of the dimeric and monomeric folding intermediates. The observation of protection exceeding that offered by the dimeric intermediate in segments from all three helices implies that a segment-swapping mechanism may be operative in the monomeric intermediate. Protection greater than that expected from the global stability for a single amide hydrogen in a peptide from the C-helix possibly and another from the A-helix may reflect non-random structure, possibly a precursor for segment swapping, in the urea-denatured state. Native topology-based model simulations that correspond to a funnel energy landscape capture both the monomeric and dimeric intermediates suggested by the HX MS data and provide a rationale for the progressive acquisition of secondary structure in their conformational ensembles.     

Opposing roles for DNA structure-specific proteins Rad1, Msh2, Msh3, and Sgs1 in yeast gene targeting

Targeted gene replacement (TGR) in yeast and mammalian cells is initiated by the two free ends of the linear targeting molecule, which invade their respective homologous sequences in the chromosome, leading to replacement of the targeted locus with a selectable gene from the targeting DNA. To study the postinvasion steps in recombination, we examined the effects of DNA structure-specific proteins on TGR frequency and heteroduplex DNA formation. In strains deleted of RAD1, MSH2, or MSH3, we find that the frequency of TGR is reduced and the mechanism of TGR is altered while the reverse is true for deletion of SGS1, suggesting that Rad1 and Msh2:Msh3 facilitate TGR while Sgs1 opposes it. The altered mechanism of TGR in the absence of Msh2:Msh3 and Rad1 reveals a separate role for these proteins in suppressing an alternate gene replacement pathway in which incorporation of both homology regions from a single strand of targeting DNA into heteroduplex with the targeted locus creates a mismatch between the selectable gene on the targeting DNA and the targeted gene in the chromosome.     

Morphogenesis of a protein: folding pathways and the energy landscape

Current knowledge on the reaction whereby a protein acquires its native three-dimensional structure was obtained by and large through characterization of the folding mechanism of simple systems. Given the multiplicity of amino acid sequences and unique folds, it is not so easy, however, to draw general rules by comparing folding pathways of different proteins. In fact, quantitative comparison may be jeopardized not only because of the vast repertoire of sequences but also in view of a multiplicity of structures of the native and denatured states. We have tackled the problem of the relationships between the sequence information and the folding pathway of a protein, using a combination of kinetics, protein engineering and computational methods, applied to relatively simple systems. Our strategy has been to investigate the folding mechanism determinants using two complementary approaches, i.e. (i) the study of members of the same family characterized by a common fold, but substantial differences in amino acid sequence, or (ii) heteromorphic pairs characterized by largely identical sequences but with different folds. We discuss some recent data on protein-folding mechanisms by presenting experiments on different members of the PDZ domain family and their circularly permuted variants. Characterization of the energetics and structures of intermediates and TSs (transition states), obtained by Φ-value analysis and restrained MD (molecular dynamics) simulations, provides a glimpse of the malleability of the dynamic states and of the role of the topology of the native states and of the denatured states in dictating folding and misfolding pathways.     

Structural comparison of the two alternative transition states for folding of TI I27

TI I27, a beta-sandwich domain from the human muscle protein titin, has been shown to fold via two alternative pathways, which correspond to a change in the folding mechanism. Under physiological conditions, TI I27 folds by a classical nucleation-condensation mechanism (diffuse transition state), whereas at extreme conditions of temperature and denaturant it switches to having a polarized transition state. We have used experimental Phi-values as restraints in ensemble-averaged molecular dynamics simulations to determine the ensembles of structures representing the two transition states. The comparison of these ensembles indicates that when native interactions are substantially weakened, a protein may still be able to fold if it can access an alternative transition state characterized by a much larger entropic contribution. Analysis of the probability distribution of Phi-values derived from ensemble averaged simulations, enables us to identify residues that form contacts in some members of the ensemble but not in others illustrating that many interactions present in transition states are not strictly required for the successful completion of the folding process.     

Denatured-state energy landscapes of a protein structural database reveal the energetic determinants of a framework model for folding

Position-specific denatured-state thermodynamics were determined for a database of human proteins by use of an ensemble-based model of protein structure. The results of modeling denatured protein in this manner reveal important sequence-dependent thermodynamic properties in the denatured ensembles as well as fundamental differences between the denatured and native ensembles in overall thermodynamic character. The generality and robustness of these results were validated by performing fold-recognition experiments, whereby sequences were matched with their respective folds based on amino acid propensities for the different energetic environments in the protein, as determined through cluster analysis. Correlation analysis between structure and energetic information revealed that sequence segments destined for β-sheet in the final native fold are energetically more predisposed to a broader repertoire of states than are sequence segments destined for α-helix. These results suggest that within the subensemble of mostly unstructured states, the energy landscapes are dominated by states in which parts of helices adopt structure, whereas structure formation for sequences destined for β-strand is far less probable. These results support a framework model of folding, which suggests that, in general, the denatured state has evolutionarily evolved to avoid low-energy conformations in sequences that ultimately adopt β-strand. Instead, the denatured state evolved so that sequence segments that ultimately adopt α-helix and coil will have a high intrinsic structure formation capability, thus serving as potential nucleation sites.     

Correct folding of the beta-barrel of the human membrane protein VDAC requires a lipid bilayer

Spontaneous membrane insertion and folding of beta-barrel membrane proteins from an unfolded state into lipid bilayers has been shown previously only for few outer membrane proteins of Gram-negative bacteria. Here we investigated membrane insertion and folding of a human membrane protein, the isoform 1 of the voltage-dependent anion-selective channel (hVDAC1) of mitochondrial outer membranes. Two classes of transmembrane proteins with either alpha-helical or beta-barrel membrane domains are known from the solved high-resolution structures. VDAC forms a transmembrane beta-barrel with an additional N-terminal alpha-helix. We demonstrate that similar to bacterial OmpA, urea-unfolded hVDAC1 spontaneously inserts and folds into lipid bilayers upon denaturant dilution in the absence of folding assistants or energy sources like ATP. Recordings of the voltage-dependence of the single channel conductance confirmed folding of hVDAC1 to its active form. hVDAC1 developed first beta-sheet secondary structure in aqueous solution, while the alpha-helical structure was formed in the presence of lipid or detergent. In stark contrast to bacterial beta-barrel membrane proteins, hVDAC1 formed different structures in detergent micelles and phospholipid bilayers, with higher content of beta-sheet and lower content of alpha-helix when inserted and folded into lipid bilayers. Experiments with mixtures of lipid and detergent indicated that the content of beta-sheet secondary structure in hVDAC1 decreased at increased detergent content. Unlike bacterial beta-barrel membrane proteins, hVDAC1 was not stable even in mild detergents such as LDAO or dodecylmaltoside. Spontaneous folding of outer membrane proteins into lipid bilayers indicates that in cells, the main purpose of membrane-inserted or associated assembly factors may be to select and target beta-barrel membrane proteins towards the outer membrane instead of actively assembling them under consumption of energy as described for the translocons of cytoplasmic membranes.     

Branching in the sequential folding pathway of cytochrome c

Previous results indicate that the folding pathways of cytochrome c and other proteins progressively build the target native protein in a predetermined stepwise manner by the sequential formation and association of native-like foldon units. The present work used native state hydrogen exchange methods to investigate a structural anomaly in cytochrome c results that suggested the concerted folding of two segments that have little structural relationship in the native protein. The results show that the two segments, an 18-residue omega loop and a 10-residue helix, are able to unfold and refold independently, which allows a branch point in the folding pathway. The pathway that emerges assembles native-like foldon units in a linear sequential manner when prior native-like structure can template a single subsequent foldon, and optional pathway branching is seen when prior structure is able to support the folding of two different foldons.     

Distribution of rare triplets along mRNA and their relation to protein folding

It is believed that pausing during mRNA translation plays some role in ensuring proper folding of newly synthesized sections of a protein chain. Such pausing occurs when rare triplets are encountered in the mRNA, as it takes additional time for the corresponding rare species of tRNA to be delivered. To determine whether pause sites are non-randomly distributed along prokaryotic mRNA (cDNA), we have located clusters of rare triplets in cDNA sequences from 21 different bacteria. From the individual profiles of local codon frequencies calculated with various windows, the positions of the clusters of the rarest codons were taken for generation of the combined histograms of positional preferences of the pause sites. The histograms show that in the prokaryotic sequences, the pause sites are located preferentially at the start positions and at about 155 triplets from the starts. To verify the generality of these observations, the data are grouped in six independent sets about 500 sequences each, all revealing the same features. A less prominent maximum is also seen at the triplet position 75. Judging by the amplitude of the peak at 155 triplets, an optimal cluster size is estimated to equal 18 triplets. The distance 155 closely corresponds to the sizes of typical protein folds and to earlier estimated prokaryotic protein sequence segments. This supports the suggestion of a role for translation pausing in the cotranslational folding of protein domains. The profiles of rare codons in mRNA can serve in the detection or prediction of boundaries between protein domains.     

Protein folding: then and now

Over the past three decades the protein folding field has undergone monumental changes. Originally a purely academic question, how a protein folds has now become vital in understanding diseases and our abilities to rationally manipulate cellular life by engineering protein folding pathways. We review and contrast past and recent developments in the protein folding field. Specifically, we discuss the progress in our understanding of protein folding thermodynamics and kinetics, the properties of evasive intermediates, and unfolded states. We also discuss how some abnormalities in protein folding lead to protein aggregation and human diseases.     

Sequence and structural homology among membrane-associated domains of CFTR and certain transporter proteins

A sequence comparison of the two membrane-associated (MA) domains of the cystic fibrosis transmembrane conductance regulator (CFTR), multidrug resistance transporter (MDR), and -factor pheromone export system (STE6) proteins, each of which are believed to contain a total of 12 transmembrane (TM) segments, reveals significant amino acid homology and length conservation in the loop regions that connect individual TM sequences. Similar structural homology is observed between these proteins, hemolysin B (HLYB) and the major histocompatibility-linked peptide transporter, HAM1, the latter two which contain a single MA domain composed of six TM segments. In addition, there are specific sequences that are conserved within the TM segments of the five different membrane proteins. This observation suggests that the folding topologies of the MA domains of MDR, STE6, and CFTR in the plasma membrane are likely to be very similar. The sequence analysis also reveals that there are three characteristic motifs (a pair of aromatic residues, LTLXXXXXXP and GXXL) that are conserved in MDR, STE6, HLYB, HAM1, but not in CFTR. We propose that although CFTR may be evolutionarily related to these other membrane proteins, it belongs to a separate subclass.     

Clinical, Molecular, and Genetic Characteristics of PAPA Syndrome: A Review

PAPA syndrome (Pyogenic Arthritis, Pyoderma gangrenosum, and Acne) is an autosomal dominant, hereditary auto-inflammatory disease arising from mutations in the PSTPIP1/CD2BP1 gene on chromosome 15q. These mutations produce a hyper-phosphorylated PSTPIP1 protein and alter its participation in activation of the "inflammasome" involved in interleukin-1 (IL-1β) production. Overproduction of IL-1β is a clear molecular feature of PAPA syndrome. Ongoing research is implicating other biochemical pathways that may be relevant to the distinct pyogenic inflammation of the skin and joints characteristic of this disease. This review summarizes the recent and rapidly accumulating knowledge on these molecular aspects of PAPA syndrome and related disorders.     

Intimate view of a kinetic protein folding intermediate: residue-resolved structure, interactions, stability, folding and unfolding rates, homogeneity

A cytochrome c kinetic folding intermediate was studied by hydrogen exchange (HX) pulse labeling. Advances in the technique and analysis made it possible to define the structured and unstructured regions, equilibrium stability, and kinetic opening and closing rates, all at an amino acid-resolved level. The entire N-terminal and C-terminal helices are formed and docked together at their normal native positions. They fray in both directions from the interaction region, due to a progression in both unfolding and refolding rates, leading to the surprising suggestion that helix propagation may proceed very slowly in the condensed milieu. Several native-like beta turns are formed. Some residues in the segment that will form the native 60s helix are protected but others are not, suggesting energy minimization to some locally non-native conformation in the transient intermediate. All other regions are unprotected, presumably dynamically disordered. The intermediate resembles a partially constructed native state. It is early, on-pathway, and all of the refolding molecules pass through it. These and related results consistently point to distinct, homogeneous, native-like intermediates in a stepwise sequential pathway, guided by the same factors that determine the native structure. Previous pulse labeling efforts have always assumed EX2 exchange during the labeling pulse, often leading to the suggestion of heterogeneous intermediates in alternative parallel pathways. The present work reveals a dominant role for EX1 exchange in the high pH labeling pulse, which will mimic heterogeneous behavior when EX2 exchange is assumed. The general problem of homogeneous versus heterogeneous intermediates and pathways is discussed.     

Developmental plasticity as a cohesive evolutionary process between sympatric alternate-year insect cohorts

Many species, particularly insects, pass through a series of distinct phases during their life history, with the developmental timing directed towards appropriate resources. Any factor that creates variation in developmental timing may partition a population into discrete populations-or 'cohorts'. Where there is continued failure to recruit outside the natal cohort then alternate cohorts will have their own internal dynamics, eventually leading to independent demographic and evolutionary trajectories. By contrast, continued variation in development rates within a cohort-cohort splitting-may homogenise otherwise independent demographic units. Using a panel of 14 microsatellite loci, we quantify the genetic signature of apparent demographic isolation between coexisting, but alternate, semivoltine cohorts of the damselfly Coenagrion mercuriale at locations that span its distribution in the UK. We find consistently low levels of genetic divergence between sympatric cohorts of C. mercuriale, indicative of developmental plasticity during the larval stage (unregulated development) whereby some individuals complete their development outside the predominant 2-year (semivoltine) period. Thus, individuals that alter their developmental rate successfully recruit to a different cohort. Despite maintaining contrasting population sizes, gene flow between alternate cohorts broadly is sufficient to place them on a similar evolutionary trajectory and also buffers against loss of genetic diversity. Such flexible larval development permits a response to local conditions and may facilitate response to environmental change.