To understand structural and thermodynamic features of disulfides within an alpha-helix, a non-redundant dataset comprising of 5025 polypeptide chains containing 2311 disulfides was examined. Thirty-five examples were found of intrahelical disulfides involving a CXXC motif between the N-Cap and third helical positions. GLY and PRO were the most common amino acids at positions 1 and 2, respectively. The N-Cap residue for disulfide bonded CXXC motifs had average (phi,psi) values of (-112 +/- 25.2 degrees , 106 +/- 25.4 degrees ). To further explore conformational requirements for intrahelical disulfides, CYS pairs were introduced at positions N-Cap-3; 1,4; 7,10 in two helices of an Escherichia coli thioredoxin mutant lacking its active site disulfide (nSS Trx). In both helices, disulfides formed spontaneously during purification only at positions N-Cap-3. Mutant stabilities were characterized by chemical denaturation studies (in both oxidized and reduced states) and differential scanning calorimetry (oxidized state only). All oxidized as well as reduced mutants were destabilized relative to nSS Trx. All mutants were redox active, but showed decreased activity relative to wild-type thioredoxin. Such engineered disulfides can be used to probe helix start sites in proteins of unknown structure and to introduce redox activity into proteins. Conversely, a protein with CYS residues at positions N-Cap and 3 of an alpha-helix is likely to have redox activity. Proteins 2010.(c) 2009 Wiley-Liss, Inc.

SecB, a soluble cytosolic chaperone component of the Sec export pathway, binds to newly synthesized precursor proteins and prevents their premature aggregation and folding and subsequently targets them to the translocation machinery on the membrane. PreMBP, the precursor form of maltose binding protein, has a 26-residue signal sequence attached to the N-terminus of MBP and is a physiological substrate of SecB. We examine the effect of macromolecular crowding and SecB on the stability and refolding of denatured preMBP and MBP. PreMBP was less stable than MBP (DeltaT(m )= 7 +/- 0.5 K) in both crowded and uncrowded solutions. Crowding did not cause any substantial changes in the thermal stability of MBP (DeltaT(m )= 1 +/- 0.4 K) or preMBP (DeltaT(m )= 0 +/- 0.6 K), as observed in spectroscopically monitored thermal unfolding experiments. However, both MBP and preMBP were prone to aggregation while refolding under crowded conditions. In contrast to MBP aggregates, which were amorphous, preMBP aggregates form amyloid fibrils. Under uncrowded conditions, a molar excess of SecB was able to completely prevent aggregation and promote disaggregation of preformed aggregates of MBP. When a complex of the denatured protein and SecB was preformed, SecB could completely prevent aggregation and promote folding of MBP and preMBP even in crowded solution. Thus, in addition to maintaining substrates in an unfolded, export-competent conformation, SecB also suppresses the aggregation of its substrates in the crowded intracellular environment. SecB is also able to promote passive disaggregation of macroscopic aggregates of MBP in the absence of an energy source such as ATP or additional cofactors. These experiments also demonstrate that signal peptide can greatly influence protein stability and aggregation propensity.

Understanding the molecular basis for the enhanced stability of proteins from thermophiles has been hindered by a lack of structural data for homologous pairs of proteins from thermophiles and mesophiles. To overcome this difficulty, complete genome sequences from 9 thermophilic and 21 mesophilic bacterial genomes were aligned with protein sequences with known structures from the protein data bank. Sequences with high homology to proteins with known structures were chosen for further analysis. High quality models of these chosen sequences were obtained using homology modeling. The current study is based on a data set of models of 900 mesophilic and 300 thermophilic protein single chains and also includes 178 templates of known structure. Structural comparisons of models of homologous proteins allowed several factors responsible for enhanced thermostability to be identified. Several statistically significant, specific amino acid substitutions that occur going from mesophiles to thermophiles are identified. Most of these are at solvent-exposed sites. Salt bridges occur significantly more often in thermophiles. The additional salt bridges in thermophiles are almost exclusively in solvent-exposed regions, and 35% are in the same element of secondary structure. Helices in thermophiles are stabilized by intrahelical salt bridges and by an increase in negative charge at the N-terminus. There is an approximate decrease of 1% in the overall loop content and a corresponding increase in helical content in thermophiles. Previously overlooked cation-pi interactions, estimated to be twice as strong as ion-pairs, are significantly enriched in thermophiles. At buried sites, statistically significant hydrophobic amino acid substitutions are typically consistent with decreased side chain conformational entropy.

Accurate identification of cavities is important in the study of protein structure, stability, design, and ligand binding. Identification and quantitation of cavities is a nontrivial problem because most cavities are connected to the protein exterior. We describe a computational procedure for quantitating cavity volumes and apply this to derive an estimate of the hydrophobic driving force in protein folding. A grid-based Monte Carlo procedure is used to position water molecules on the surface of a protein. A Voronoi procedure is used to identify and quantitate empty space within the solvated protein. Additional cavities not detected by other existing procedures can be identified. Most of these are close to surface concavities. Residue volumes for both the interior and the surface residues as well as cavity volumes are in good agreement with volumes calculated from fully hydrated protein structures obtained from molecular dynamic simulations. We show that the loss of stability because of cavity-creating mutations correlates better with cavity volumes determined by this procedure than with cavity volumes determined by other methods. Available structural and thermodynamic data for a number of cavity-containing mutants were analyzed to obtain estimates of 26.1 cal x mol(-1) x A(-3) and 18.5 cal x mol(-1) x A(-2) for the relative contributions of cavity formation and the hydrophobic effect to the observed stability changes. The present estimate for the hydrophobic driving force is at the lower end of estimates derived from model compound studies and considerably lower than previous estimates of approximately 50 cal x mol(-1) x A(-2) derived from protein mutational data. In the absence of structural rearrangement, on average, deletion of a single methylene group is expected to result in losses in stability of 0.41 and 0.70 kcal x mol(-1) resulting from decrease in hydrophobicity and packing, respectively.

The amino acid Pro is more rigid than other naturally occurring amino acids and, in proteins, lacks an amide hydrogen. To understand the structural and thermodynamic effects of Pro substitutions, it was introduced at 13 different positions in four different proteins, leucine-isoleucine-valine binding protein, maltose binding protein, ribose binding protein, and thioredoxin. Three of the maltose binding protein mutants were characterized by X-ray crystallography to confirm that no structural changes had occurred upon mutation. In the remaining cases, fluorescence and CD spectroscopy were used to show the absence of structural change. Stabilities of wild type and mutant proteins were characterized by chemical denaturation at neutral pH and by differential scanning calorimetry as a function of pH. The mutants did not show enhanced stability with respect to chemical denaturation at room temperature. However, 6 of the 13 single mutants showed a small but significant increase in the free energy of thermal unfolding in the range of 0.3-2.4 kcal/mol, 2 mutants showed no change, and 5 were destabilized. In five of the six cases, the stabilization was because of reduced entropy of unfolding. However, the magnitude of the reduction in entropy of unfolding was typically several fold larger than the theoretical estimate of -4 cal K(-1) mol(-1) derived from the relative areas in the Ramachandran map accessible to Pro and Ala residues, respectively. Two double mutants were constructed. In both cases, the effects of the single mutations on the free energy of thermal unfolding were nonadditive. Proteins 2007.(c) 2006 Wiley-Liss, Inc.

Calculations predict that cation-&pgr; interactions make an important contribution to protein stability. While there have been some attempts to experimentally measure strengths of cation-&pgr; interactions using peptide model systems, much less experimental data are available for globular proteins. We have attempted to determine the magnitude of cation-&pgr; interactions of Lys with aromatic amino acids in four different proteins (LIVBP, MBP, RBP, and Trx). In each case, Lys was replaced with Gln and Met. In a separate series of experiments, the aromatic amino acid in each cation-&pgr; pair was replaced by Leu. Stabilities of wild-type (WT) and mutant proteins were characterized by both thermal and chemical denaturation. Gln and aromatic --> Leu mutants were consistently less stable than corresponding Met mutants, reflecting the nonisosteric nature of these substitutions. The strength of the cation-&pgr; interaction was assessed by the value of the change in the free energy of unfolding [&Dgr;&Dgr;G degrees =&Dgr;G degrees (Met)-&Dgr;G degrees (WT)]. This ranged from +1.1 to -1.9 kcal/mol (average value -0.4 kcal/mol) at 298 K and +0.7 to -2.6 kcal/mol (average value -1.1 kcal/mol) at the Tm of each WT. It therefore appears that the strength of cation-&pgr; interactions increases with temperature. In addition, the experimentally measured values are appreciably smaller in magnitude than calculated values with an average difference |&Dgr;G degrees expt -&Dgr;G degrees calc|av of 2.9 kcal/mol. At room temperature, the data indicate that cation-&pgr; interactions are at best weakly stabilizing and in some cases are clearly destabilizing. However, at elevated temperatures, close to typical Tm's, cation-&pgr; interactions are generally stabilizing.

gp120 is a subunit of the envelope glycoprotein of HIV-1. The third variable loop region of gp120 (V3 loop) contains multiple immunodominant epitopes and is also functionally important for deciding cell-tropism of the virus. 447-52D is a monoclonal antibody that recognizes the conserved tip of the V3 loop in a beta turn conformation. This antibody has previously been shown to neutralize diverse strains of the virus. In an attempt to generate an immunogen competent to generate 447-52D like antibodies, the known epitope of 447-52D was inserted at three different surface loop locations in the small, stable protein E. coli Thioredoxin (Trx). At one of the three locations (between residues 74 and 75) the insertion was tolerated, the resulting protein was stable and soluble, and bound 447-52D with an affinity similar to that of intact gp120. Upon immunization, the V3 peptide-inserted Trx scaffold was able to generate anti-V3 antibodies that could compete out 447-52D binding to gp120. Epitope mapping studies demonstrated that these anti-V3 antibodies recognized the same epitope as 447-52d. Although the 447-52D type antibodies were estimated to be present at concentrations of 50-400 mug/ ml of serum, these were not able to effect neutralization of strains like JRFL and BAL but could neutralize the sensitive MN strain. The data suggest that because of the low accessibility of the V3 loop on primary isolates such as JRFL, it will be difficult to elicit a V3-specific, 447-52D-like antibody response to effectively neutralize such isolates.

Every residue of the 101-aa Escherichia coli toxin CcdB was substituted with Ala, Asp, Glu, Lys, and Arg by using site-directed mutagenesis. The activity of each mutant in vivo was characterized as a function of Controller of Cell Division or Death B protein (CcdB) transcriptional level. The mutation data suggest that an accessibility value of 5% is an appropriate cutoff for definition of buried residues. At all buried positions, introduction of Asp results in an inactive phenotype at all CcdB transcriptional levels. The average amount of destabilization upon substitution at buried positions decreases in the order Asp>Glu>Lys>Arg>Ala. Asp substitutions at buried sites in two other proteins, maltose-binding protein and thioredoxin, also were shown to be severely destabilizing. Ala and Asp scanning mutagenesis, in combination with dose-dependent expression phenotypes, was shown to yield important information on protein structure and activity. These results also suggest that such scanning mutagenesis data can be used to rank order sequence alignments and their corresponding homology models, as well as to distinguish between correct and incorrect structural alignments. With continuous reductions in oligonucleotide costs and increasingly efficient site-directed mutagenesis procedures, comprehensive scanning mutagenesis experiments for small proteins/domains are quite feasible.

Disulfides cross-link residues in a protein that are separated in primary sequence and stabilize the protein through entropic destabilization of the unfolded state. While the removal of naturally occurring disulfides leads to protein destabilization, introduction of engineered disulfides does not always lead to significant stabilization of a protein. We have analyzed naturally occurring disulfides that span adjacent antiparallel strands of beta sheets (cross-strand disulfides). Cross-strand disulfides have recently been implicated as redox-based conformational switches in proteins such as gp120 and CD4. The propensity of these disulfides to act as conformational switches was postulated on the basis of the hypothesis that this class of disulfide is conformationally strained. In the present analysis, there was no evidence to suggest that cross-strand disulfides are more strained compared to other disulfides as assessed by their torsional energy. It was also observed that these disulfides occur solely at non-hydrogen-bonded (NHB) registered pairs of adjacent antiparallel strands and not at hydrogen-bonded (HB) positions as suggested previously. One of the half-cystines involved in cross-strand disulfide formation often occurs at an edge strand. Experimental confirmation of the stabilizing effects of such disulfides was carried out in Escherichia coli thioredoxin. Four pairs of cross-strand cysteines were introduced, two at HB and two at NHB pairs. Disulfides were formed in all four cases. However, as predicted from our analysis, disulfides at NHB positions resulted in an increase in melting temperature of 7-10 degrees C, while at HB positions there was a corresponding decrease of -7 degrees C. The reduced state of all proteins had similar stability.

When incorporated into a polypeptide chain, proline (Pro) differs from all other naturally occurring amino acid residues in two important respects. The varphi dihedral angle of Pro is constrained to values close to -65 degrees and Pro lacks an amide hydrogen. Consequently, mutations which result in introduction of Pro can significantly affect protein stability. In the present work, we describe a procedure to accurately predict the effect of Pro introduction on protein thermodynamic stability. Seventy-seven of the 97 non-Pro amino acid residues in the model protein, CcdB, were individually mutated to Pro, and the in vivo activity of each mutant was characterized. A decision tree to classify the mutation as perturbing or nonperturbing was created by correlating stereochemical properties of mutants to activity data. The stereochemical properties including main chain dihedral angle varphi and main chain amide H-bonds (hydrogen bonds) were determined from 3D models of the mutant proteins built using MODELLER. We assessed the performance of the decision tree on a large dataset of 163 single-site Pro mutations of T4 lysozyme, 74 nsSNPs, and 52 other Pro substitutions from the literature. The overall accuracy of this algorithm was found to be 81% in the case of CcdB, 77% in the case of lysozyme, 76% in the case of nsSNPs, and 71% in the case of other Pro substitution data. The accuracy of Pro scanning mutagenesis for secondary structure assignment was also assessed and found to be at best 69%. Our prediction procedure will be useful in annotating uncharacterized nsSNPs of disease-associated proteins and for protein engineering and design.

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of DeltaCp of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Little is known about mechanisms responsible for the temperature-sensitive (ts) phenotype, or of the transferability of ts mutants of a specific gene between organisms. Using a structure-based approach, nine ts mutants of Gal4 were generated in yeast by mutating four DNA binding residues. Two of these nine yeast ts mutants were cloned into P element vectors under control of the Elav and GMR promoters and transgenic Drosophila lines were generated. These were crossed to UAS reporter lines and progeny were characterized for reporter gene expression as a function of temperature. Both of these yeast ts mutants show a ts phenotype in Drosophila and result in rapid induction of reporter gene expression upon shifting to the permissive temperature. Exposed, functional residues involved in protein-ligand or protein-protein interactions appear to be attractive candidate sites for generating ts mutants that are transferable between organisms.