Submillisecond folding of the peripheral subunit-binding domain.

Folding and unfolding rates have been measured for the peripheral subunit-binding domain, a small three-helix protein. The protein folds very fast, with rates too rapid to be measured using traditional stopped-flow techniques. Folding and unfolding rates were measured as a function of temperature using dynamic NMR lineshape analysis. At the lowest temperature at which there is sufficient broadening to measure rates, 41 degrees C, the folding rate is 16,050 s(-1). Thus, the halftime required for folding is 43 micros. At the same temperature, the unfolding rate is 2800 s(-1). Identical rates were measured using resolved resonances from Val16 in the loop and Val21 at the end of the 310-helix. Folding rates have been correlated with protein topology, and this correlation is consistent with the rapid folding of the peripheral subunit-binding domain. The results presented here show that the peripheral subunit-binding domain is the third fastest folding protein for which rates have been estimated. The folding rate is the fastest that has been directly measured and provides further support for the importance of chain topology as a major determinant of folding rates.     

Ultra-fast barrier-limited folding in the peripheral subunit-binding domain family

We have determined the solution structures, equilibrium properties and ultra-fast folding kinetics for three bacterial homologues of the peripheral subunit-binding domain (PSBD) family. The mesophilic homologue, BBL, was less stable than the thermophilic and hyper-thermophilic variants (E3BD and POB, respectively). The broad unfolding transitions of each PSBD, when probed by different techniques, were essentially superimposable, consistent with co-operative denaturation. Temperature-jump and continuous-flow fluorescence methods were used to measure the folding kinetics for E3BD, POB and BBL. E3BD folded fairly rapidly at 298K (folding half-time approximately 25 micros) and BBL and POB folded even faster (folding half-times approximately 3-5 micros). The variations in equilibrium and kinetic behaviour observed for the PSBD family resembles that of the homeodomain family, where the folding pattern changes from apparent two-state transitions to multi-state kinetics as the denatured state becomes more structured. The faster folding of POB may be a consequence of its higher propensity to form helical structure in the region corresponding to the folding nucleus of E3BD. The ultra-fast folding of BBL appears to be a consequence of residual structure in the denatured ensemble, as with engrailed homeodomain. We discuss issues concerning "one-state", downhill folding, and find no evidence for, and strong evidence against, it occurring in these PSBDs. The shorter construct used previously for BBL was destabilized significantly and the stability further perturbed by the introduction of fluorescent probes. Thermal titrations for 11 side-chains scattered around the protein, when probed by (13)C-NMR experiments, could be fit globally to a common co-operative transition.     

Conservation of transition state structure in fast folding peripheral subunit-binding domains

Φ-Value analysis was used to characterise the structure of the transition state (TS) for folding of POB L146A Y166W, a peripheral subunit-binding domain that folds in microseconds. Helix 2 was structured in the TS with consolidating interactions from the structured loop that connects the two α-helices. This distribution of Φ-values was very similar to that determined for E3BD F166W, a homologue with high sequence and structural similarity. The extrapolated folding rate constants in water at 298 K were 210,000 s^− 1 for POB and 27,500 s^− 1 for E3BD. A contribution to the faster folding of POB came from its having significantly greater helical propensity in helix 2, the folding nucleus. The folding rate also appeared to be influenced by differences in the sequence and structural properties of the loop connecting the two helices. Unimodal downhill folding has been proposed as a conserved, biologically important property of peripheral subunit-binding domains. POB folds five times faster and E3BD folds slower than a proposed limit of 40,000 s^− 1 for barrier-limited folding. However, experimental evidence strongly suggests that both POB L146A Y166W and E3BD F166W fold in a barrier-limited process through a very similar TS ensemble.     

Effects of varying the local propensity to form secondary structure on the stability and folding kinetics of a rapid folding mixed alpha/beta protein: characterization of a truncation mutant of the N-terminal domain of the ribosomal protein L9.

The N-terminal domain of the ribosomal protein L9 forms a split betaalphabeta structure with a long C-terminal helix. The folding transitions of a 56 residue version of this protein have previously been characterized, here we report the results of a study of a truncation mutant corresponding to residues 1-51. The 51 residue protein adopts the same fold as the 56 residue protein as judged by CD and two-dimensional NMR, but it is less stable as judged by chemical and thermal denaturation experiments. Studies with synthetic peptides demonstrate that the C-terminal helix of the 51 residue version has very little propensity to fold in isolation in contrast to the C-terminal helix of the 56 residue variant. The folding rates of the two proteins, as measured by stopped-flow fluorescence, are essentially identical, indicating that formation of local structure in the C-terminal helix is not involved in the rate-limiting step of folding.     

The transition state for folding of a peripheral subunit-binding domain contains robust and ionic-strength dependent characteristics

The denaturant dependencies of the folding and unfolding kinetics were used to characterize the structure of the transition state for folding of E3BD, a peripheral subunit-binding domain. For the majority of E3BD mutants, the Phi-values calculated at 298 K from the analysis of chevron plots were in good agreement with those previously determined at 325 K using Arrhenius analysis. This agreement further demonstrates the general robustness of Phi-value analyses, since different experiments, methods of denaturation and thermodynamic assumptions were used to determine each set of Phi(F) values. The structure of the transition state for folding was grossly conserved at 298 K and 325 K, with residues in Helix I playing a lesser role in folding than those located in the 3(10) helix, disordered loop and Helix II. However, the energetic contributions of a cluster of basic residues close to the N-terminus and Helix I, which are an integral part of the ligand-binding site, were susceptible to ionic strength effects because of electrostatic strain in native and transition states of E3BD at low ionic strength. We found no evidence of the downhill folding previously proposed for E3BD, even though the conditions employed in this study significantly increased the energetic bias towards the native state.     

Conformational analysis of peptide fragments derived from the peripheral subunit-binding domain from the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus: evidence for nonrandom structure in the unfolded state

There is currently a great deal of interest in the early events in protein folding. Two issues that have generated particular interest are the nature of the unfolded state under native conditions and the role of local interactions in folding. Here, we report the results of a study of a set of peptides derived from a small two-helix protein, the peripheral subunit-binding domain of the pyruvate dehydrogenase multienzyme complex. Five peptides of overlapping sequence were prepared, including sequences corresponding to each of the helices and to the region connecting them. The peptides were characterized by CD and, where possible, nmr. A peptide corresponding to the second helix is between 12 and 17% helical at neutral pH. CD also indicates a lower percentage of helical structure in the peptide corresponding to the first alpha-helix, although the values of the alpha-proton chemical shifts suggest some preference for nonrandom structure. Peptides corresponding to the interhelical loop, which in the full domain contains two overlapping beta-turns and a 5-residue 3(10)-helix, are less structured. There is no significant change in the helicity of any of these peptides with pH. To test for fragment complementation, CD spectra of the two peptides derived from each helix and the long connecting peptide were compared to the spectra of each possible pair, as well as to a mixture containing all three. No increase in structure was observed. We complement our peptide studies by characterizing a point mutant, D34V, which disrupts a critical hydrogen bonding network. This mutant is unable to fold and provides a useful model of the denatured state. The mutant is between 9 and 16% helical as judged by CD. The modest amount of helical structure formed in some of the peptide fragments and in the point mutant suggests that the denatured state of the peripheral subunit binding domain is not completely unstructured. This may contribute to the very rapid folding observed for the intact protein.     

Solution structure of a double mutant of the carboxy-terminal dimerization domain of the HIV-1 capsid protein

As in other retroviruses, the HIV-1 capsid (CA) protein is composed of two domains, the N-terminal domain (NTD) and the C-terminal domain (CTD), joined by a flexible linker. The dimerization of the CTD is thought to be a critical step in the assembly of the immature and mature viral capsids. The precise nature of the functional form of CTD dimerization interface has been a subject of considerable interest. Previously, the CTD dimer was thought to involve a face-to-face dimerization observed in the early crystallographic studies. Recently, the crystallographic structure for a domain-swapped CTD dimer has been determined. This dimer, with an entirely different interface that includes the major homology region (MHR) has been suggested as the functional form during the Gag assembly. The structure determination of the monomeric wt CTD of HIV-1 has not been possible because of the monomer-dimer equilibrium in solution. We report the NMR structure of the [W184A/M185A]-CTD mutant in its monomeric form. These mutations interfere with dimerization without abrogating the assembly activity of Gag and CA. The NMR structure shows some important differences compared to the CTD structure in the face-to-face dimer. Notably, the helix-2 is much shorter, and the kink seen in the crystal structure of the wt CTD in the face-to-face dimer is absent. These NMR studies suggest that dimerization-induced conformational changes may be present in the two crystal structures of the CTD dimers and also suggest a mechanism that can simultaneously accommodate both of the distinctly different dimer models playing functional roles during the Gag assembly of the immature capsids.     

Influence of pH on NMR structure and stability of the human prion protein globular domain

The NMR structure of the globular domain of the human prion protein (hPrP) with residues 121-230 at pH 7.0 shows the same global fold as the previously published structure determined at pH 4.5. It contains three alpha-helices, comprising residues 144-156, 174-194, and 200-228, and a short anti-parallel beta-sheet, comprising residues 128-131 and 161-164. There are slight, strictly localized, conformational changes at neutral pH when compared with acidic solution conditions: helix alpha1 is elongated at the C-terminal end with residues 153-156 forming a 310-helix, and the population of helical structure in the C-terminal two turns of helix alpha 2 is increased. The protonation of His155 and His187 presumably contributes to these structural changes. Thermal unfolding monitored by far UV CD indicates that hPrP-(121-230) is significantly more stable at neutral pH. Measurements of amide proton protection factors map local differences in protein stability within residues 154-157 at the C-terminal end of helix alpha 1 and residues 161-164 of beta-strand 2. These two segments appear to form a separate domain that at acidic pH has a larger tendency to unfold than the overall protein structure. This domain could provide a "starting point" for pH-induced unfolding and thus may be implicated in endosomic PrPC to PrPSc conformational transition resulting in transmissible spongiform encephalopathies.     

Interactions of the peripheral subunit-binding domain of the dihydrolipoyl acetyltransferase component in the assembly of the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus.

The enzymes pyruvate decarboxylase (E1) and dihydrolipoyl dehydrogenase (E3) bind tightly but in a mutually exclusive manner to the peripheral subunit-binding domain (PSBD) of dihydrolipoyl acetyltransferase in the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus. The use of directed mutagenesis, surface plasmon resonance detection and isothermal titration microcalorimetry revealed that several positively charged residues of the PSBD, most notably Arg135, play an important part in the interaction with both E1 and E3, whereas Met131 makes a significant contribution to the binding of E1 only. This indicates that the binding sites for E1 and E3 on the PSBD are overlapping but probably significantly different, and that additional hydrophobic interactions may be involved in binding E1 compared with E3. Arg135 of the PSBD was also replaced with cysteine (R135C), which was then modified chemically by alkylation with increasingly large aliphatic groups (R135C -methyl, -ethyl, -propyl and -butyl). The pattern of changes in the values of DeltaG degrees, DeltaH degrees and TDeltaS degrees that were found to accompany the interaction with the variant PSBDs differed between E1 and E3 despite the similarities in the free energies of their binding to the wild-type. The importance of a positive charge on the side-chain at position 135 for the interaction of the PSBD with E3 and E1 was apparent, although lysine was found to be an imperfect substitute for arginine. The results offer further evidence of entropy-enthalpy compensation ('thermodynamic homeostasis') - a feature of systems involving a multiplicity of weak interactions.     

Crystal structure of yeast mitochondrial peripheral membrane protein Tim44p C-terminal domain

The protein transports from the cell cytosol to the mitochondria matrix are carried out by the translocase of the outer membrane (TOM) complex and the translocase of the inner membrane (TIM) complexes. Tim44p is an essential mitochondrial peripheral membrane protein and a major component of TIM23 translocon. Tim44p can tightly associate with the inner mitochondrial membrane. To investigate the mechanism by which Tim44p functions in the TIM23 translocon to deliver the mitochondrial protein precursors, we have determined the crystal structure of the yeast Tim44p C-terminal domain to 3.2A resolution using the MAD method. The Tim44p C-terminal domain forms a monomer in the crystal structure and contains six alpha-helices and four antiparallel beta-strands. A large hydrophobic pocket was identified on the Tim44p structure surface. The N-terminal helix A1 is positively charged and the helix A1 protrudes out from the Tim44p main body.     

Structural basis of the enhanced stability of a mutant ribozyme domain and a detailed view of RNA--solvent interactions

BACKGROUND: The structure of P4-P6, a 160 nucleotide domain of the self-splicing Tetrahymena thermophila intron, was solved previously. Mutants of the P4-P6 RNA that form a more stable tertiary structure in solution were recently isolated by successive rounds of in vitro selection and amplification. RESULTS: We show that a single-site mutant (Delta C209) possessing greater tertiary stability than wild-type P4-P6 also crystallizes much more rapidly and under a wider variety of conditions. The crystal structure provides a satisfying explanation for the increased stability of the mutant; the deletion of C209 allows the adjacent bulged adenine to enter the P4 helix and form an A-G base pair, presumably attenuating the conformational flexibility of the helix. The structure of another mutant (Delta A210) was also solved and supports this interpretation. The crystals of Delta C209 diffract to a higher resolution limit than those of wild-type RNA (2.25 A versus 2.8 A), allowing assignment of innersphere and outersphere coordination contacts for 27 magnesium ions. Structural analysis reveals an intricate solvent scaffold with a preponderance of ordered water molecules on the inside rather than the surface of the folded RNA domain. CONCLUSIONS: In vitro evolution facilitated the identification of a highly stable, structurally homogeneous mutant RNA that was readily crystallizable. Analysis of the structure suggests that improving RNA secondary structure can stabilize tertiary structure and perhaps promote crystallization. In addition, the higher resolution model provides new details of metal ion-RNA interactions and identifies a core of ordered water molecules that may be integral to RNA tertiary structure formation.     

Increased rigidity of domain structures enhances the stability of a mutant enzyme created by directed evolution

A mutant of kanamycin nucleotidyltransferase (KNT) was previously created by directed evolution. This mutant, HTK, has 19 amino acid substitutions, which increase the thermostability by 20 degrees C. In this study, we have examined to what extent each mutation contributes to the increased stability and analyzed how the mutations affect the structure of KNT at 72 degrees C using molecular dynamics simulations. The effects of some mutations on the stability are simply additive, but those of others are cooperative. Mutations with large effects on the stability are introduced into the N-terminal domain, which appears to be less stable than the C-terminal domain. Results of the molecular dynamics simulations have indicated that the rigidity of the domain structures is increased by the mutations: at 72 degrees C, the intradomain fluctuations of HTK are decreased, and in turn, its interdomain motions are pronounced, whereas the structure of the preevolved KNT fluctuates randomly. Chemical modification experiments of cysteine residues have shown that the cysteine residues of HTK are less accessible to an SH reagent than those of the preevolved KNT. The present results suggest that the 19 mutations of HTK stabilize KNT by affecting the dynamic behavior of the structure of this enzyme without significantly changing its static overall structure.     

Solution structure of a conserved domain of antizyme: a protein regulator of polyamines

Antizyme and its isoforms are members of an unusual yet broadly conserved family of proteins, with roles in regulating polyamine levels within cells. Antizyme has the ability to bind and inhibit the enzyme ornithine decarboxylase (ODC), targeting it for degradation at the proteasome; antizyme is also known to affect the transport of polyamines and interact with the antizyme inhibitor protein (AZI), as well as the cell-cycle protein cyclin D1. In the present work, NMR methods were used to determine the solution structure of a stable, folded domain of mammalian antizyme isoform-1 (AZ-1), consisting of amino acid residues 87-227. The protein was found to contain eight beta strands and two alpha helices, with the strands forming a mixed parallel and antiparallel beta sheet. At the level of primary sequence, antizyme is not similar to any protein of known structure, and results show that antizyme exhibits a novel arrangement of its strands and helices. Interestingly, however, the fold of antizyme is similar to that found in a family of acetyl transferases, as well as translation initiation factor IF3, despite a lack of functional relatedness between these proteins. Structural results, combined with amino acid sequence comparisons, were used to identify conserved features among the various homologues of antizyme and their isoforms. Conserved surface residues, including a cluster of acidic amino acids, were found to be located on a single face of antizyme, suggesting this surface is a possible site of interaction with target proteins such as ODC. This structural model provides an essential framework for an improved future understanding of how the different parts of antizyme play their roles in polyamine regulation.     

Exploring the molecular linkage of protein stability traits for enzyme optimization by iterative truncation and evolution

The stability of proteins is paramount for their therapeutic and industrial use and, thus, is a major task for protein engineering. Several types of chemical and physical stabilities are desired, and discussion revolves around whether each stability trait needs to be addressed separately and how specific and compatible stabilizing mutations act. We demonstrate a stepwise perturbation-compensation strategy, which identifies mutations rescuing the activity of a truncated TEM β-lactamase. Analyses relating structural stress with the external stresses of heat, denaturants, and proteases reveal our second-site suppressors as general stability centers that also improve the full-length enzyme. A library of lactamase variants truncated by 15 N-terminal and three C-terminal residues (Bla-NΔ15CΔ3) was subjected to activity selection and DNA shuffling. The resulting clone with the best in vivo performance harbored eight mutations, surpassed the full-length wild-type protein by 5.3 °C in T(m), displayed significantly higher catalytic activity at elevated temperatures, and showed delayed guanidine-induced denaturation. The crystal structure of this mutant was determined and provided insights into its stability determinants. Stepwise reconstitution of the N- and C-termini increased its thermal, denaturant, and proteolytic resistance successively, leading to a full-length enzyme with a T(m) increased by 15.3 °C and a half-denaturation concentration shifted from 0.53 to 1.75 M guanidinium relative to that of the wild type. These improvements demonstrate that iterative truncation-optimization cycles can exploit stability-trait linkages in proteins and are exceptionally suited for the creation of progressively stabilized variants and/or downsized proteins without the need for detailed structural or mechanistic information.     

Structure-based mutant stability predictions on proteins of unknown structure

The ability to rapidly and accurately predict the effects of mutations on the physicochemical properties of proteins holds tremendous importance in the rational design of modified proteins for various types of industrial, environmental or pharmaceutical applications, as well as in elucidating the genetic background of complex diseases. In many cases, the absence of an experimentally resolved structure represents a major obstacle, since most currently available predictive software crucially depend on it. We investigate here the relevance of combining coarse-grained structure-based stability predictions with a simple comparative modeling procedure. Strikingly, our results show that the use of average to high quality structural models leads to virtually no loss in predictive power compared to the use of experimental structures. Even in the case of low quality models, the decrease in performance is quite limited and this combined approach remains markedly superior to other methods based exclusively on the analysis of sequence features.     

The solution structure of the FATC domain of the protein kinase target of rapamycin suggests a role for redox-dependent structural and cellular stability

The target of rapamycin (TOR) is a highly conserved Ser/Thr kinase that plays a central role in the control of cellular growth. TOR has a characteristic multidomain structure. Only the kinase domain has catalytic function; the other domains are assumed to mediate interactions with TOR substrates and regulators. Except for the rapamycin-binding domain, there are no high-resolution structural data available for TOR. Here, we present a structural, biophysical, and mutagenesis study of the extremely conserved COOH-terminal FATC domain. The importance of this domain for TOR function has been highlighted in several publications. We show that the FATC domain, in its oxidized form, exhibits a novel structural motif consisting of an alpha-helix and a COOH-terminal disulfide-bonded loop between two completely conserved cysteine residues. Upon reduction, the flexibility of the loop region increases dramatically. The structural data, the redox potential of the disulfide bridge, and the biochemical data of a cysteine to serine mutant indicate that the intracellular redox potential can affect the cellular amount of the TOR protein via the FATC domain. Because the amount of TOR mRNA is not changed, the redox state of the FATC disulfide bond is probably influencing the degradation of TOR.     

Prediction of protein mutant stability using classification and regression tool

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and the solving of which would help for designing stable mutants. In this work, we have analyzed the stability of protein mutants using two different datasets of 1396 and 2204 mutants obtained from ProTherm database, respectively for free energy change due to thermal (DeltaDeltaG) and denaturant denaturations (DeltaDeltaG(H(2)O)). We have used a set of 48 physical, chemical energetic and conformational properties of amino acid residues and computed the difference of amino acid properties for each mutant in both sets of data. These differences in amino acid properties have been related to protein stability (DeltaDeltaG and DeltaDeltaG(H(2)O)) and are used to train with classification and regression tool for predicting the stability of protein mutants. Further, we have tested the method with 4 fold, 5 fold and 10 fold cross validation procedures. We found that the physical properties, shape and flexibility are important determinants of protein stability. The classification of mutants based on secondary structure (helix, strand, turn and coil) and solvent accessibility (buried, partially buried, partially exposed and exposed) distinguished the stabilizing/destabilizing mutants at an average accuracy of 81% and 80%, respectively for DeltaDeltaG and DeltaDeltaG(H(2)O). The correlation between the experimental and predicted stability change is 0.61 for DeltaDeltaG and 0.44 for DeltaDeltaG(H(2)O). Further, the free energy change due to the replacement of amino acid residue has been predicted within an average error of 1.08 kcal/mol and 1.37 kcal/mol for thermal and chemical denaturation, respectively. The relative importance of secondary structure and solvent accessibility, and the influence of the dataset on prediction of protein mutant stability have been discussed.     

Electron transfer and stability of the cytochrome b6f complex in a small domain deletion mutant of cytochrome f

The lumen segment of cytochrome f consists of a small and a large domain. The role of the small domain in the biogenesis and stability of the cytochrome b(6)f complex and electron transfer through the cytochrome b(6)f complex was studied with a small domain deletion mutant in Chlamydomonas reinhardtii. The mutant is able to grow photoautotrophically but with a slower rate than the wild type strain. The heme group is covalently attached to the polypeptide, and the visible absorption spectrum of the mutant protein is identical to that of the native protein. The kinetics of electron transfer in the mutant were measured by flash kinetic spectroscopy. Our results show that the rate for the oxidation of cytochrome f was unchanged (t(12) = approximately 100 micros), but the half-time for the reduction of cytochrome f is increased (t(12) = 32 ms; for wild type, t(12) = 2.1 ms). Cytochrome b(6) reduction was slower than that of the wild type by a factor of approximately 2 (t(12) = 8.6 ms; for wild type, t(12) = 4.7 ms); the slow phase of the electrochromic band shift also displayed a slower kinetics (t(12) = 5.5 ms; for wild type, t(12) = 2.7 ms). The stability of the cytochrome b(6)f complex in the mutant was examined by following the kinetics of the degradation of the individual subunits after inhibiting protein synthesis in the chloroplast. The results indicate that the cytochrome b(6)f complex in the small domain deletion mutant is less stable than in the wild type. We conclude that the small domain is not essential for the biogenesis of cytochrome f and the cytochrome b(6)f complex. However, it does have a role in electron transfer through the cytochrome b(6)f complex and contributes to the stability of the complex.     

Estimation of the evolutionary stability of the Influenza A virus: Prediction of variable regions in the domain structure of the M1 protein

Influenza virus is a human pathogen that is responsible for several pandemics with high death rates. Two programs were written for the statistical analysis of the coding regions of genes from extensive samples of Influenza A virus. While inner viral proteins appeared to be evolutionarily stable, surface antigens and the NS1 nonstructural protein are highly variable. Using programs to predict the evolutionarily variable regions inside the M1 protein sequence it was shown that the N-domain is the most conserved domain while the M- and C-domains are the most variable. In addition, the C-domain was shown to be the most unfolded one, according to prediction by the DISOPRED algorithm and, thus, possesses the most structural plasticity.