Increasing the conformational stability by replacement of heme axial ligand in c-type cytochrome

To investigate the role of the heme axial ligand in the conformational stability of c-type cytochrome, we constructed M58C and M58H mutants of the red alga Porphyra yezoensis cytochrome c(6) in which the sixth heme iron ligand (Met58) was replaced with Cys and His residues, respectively. The Gibbs free energy change for unfolding of the M58H mutant in water (DeltaG degrees (unf)=1.48 kcal/mol) was lower than that of the wild-type (2.43 kcal/mol), possibly due to the steric effects of the mutation on the apoprotein structure. On the other hand, the M58C mutant exhibited a DeltaG degrees (unf) of 5.45 kcal/mol, a significant increase by 3.02 kcal/mol compared with that of wild-type. This increase was possibly responsible for the sixth heme axial bond of M58C mutant being more stable than that of wild-type according to the heme-bound denaturation curve. Based on these observations, we propose that the sixth heme axial ligand is an important key to determine the conformational stability of c-type cytochromes, and the sixth Cys heme ligand will give stabilizing effects.     

Lysines in the amino-terminal alpha-helix are important to the stability of Rhodobacter capsulatus cytochrome c2.

The protein stabilities of wild type and four site-directed mutants of Rhodobacter capsulatus cytochrome c2 have been characterized. The integrity of the cytochrome c2 iron-sulfur environment was ascertained by titration of the 696-nm absorbance band with alkali, and the conformational stability was determined by titration of the 220-nm circular dichroism signal with Gdn-HCl. Analysis of the alkaline transition pK value of K12D (lysine-12 substituted by aspartate) indicated that the K12D iron-sulfur environment was destabilized by 0.6 kcal/mol relative to the wild-type cytochrome c2 at low ionic strength. In contrast, the alkaline transition pK values of K14E (lysine-14 substituted by glutamate), K32E (lysine-32 substituted by glutamate), and K14E/K32E (lysines-14 and -32 substituted by glutamates) were indistinguishable from the wild type, indicating that these substitutions have no effect on the stability of the iron-sulfur environment. Gdn-HCl denaturation of K12D and K14E indicated that both these mutations decreased conformational stability by 1.3 kcal/mol. In contrast, mutant K32E exhibited a small stabilizing effect of 0.2 kcal/mol. Gdn-HCl denaturation of K14E/K32E indicated that this mutation decreased conformational stability by 1.3 kcal/mol, which is consistent with the additive effects of the single charge mutations at positions 14 and 32. The conformational instability of mutants possessing negative charges at position 12 or 14 is best explained by their positioning at the carboxy-terminal region of the amino-terminal alpha-helix of R. capsulatus cytochrome c2. Accordingly, introduction of negatively charged groups into this region appears to destabilize cytochrome c2 through energetically unfavorable interactions with the dipole of the amino-terminal helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

Increasing Protein Stability by Polar Surface Residues: Domain-Wide Consequences of Interactions Within a Loop

We have examined the influence of surface hydrogen bonds on the stability of proteins by studying the effects of mutations of human immunoglobulin light chain variable domain (V_L). In addition to the variants Y27dD, N28F, and T94H of protein κIV Len that were previously described, we characterized mutants M4L, L27cN, L27cQ, and K39T, double mutant M4L/Y27dD, and triple mutant M4L/Y27dD/T94H. The triple mutant had an enhanced thermodynamic stability of 4.2 kcal/mol. We determined the structure of the triple mutant by x-ray diffraction and correlated the changes in stability due to the mutations with changes in the three-dimensional structure. Y27dD mutant had increased stability of Len by 2.7 kcal/mol, a large value for a single mutation. Asp27d present in CDR1 formed hydrogen bonds with the side-chain and main-chain atoms within the loop. In the case of the K39T mutant, which reduces stability by 2 kcal/mol, Lys39 in addition to forming a hydrogen bond with a carbonyl oxygen of a neighboring loop may also favorably influence the surface electrostatics of the molecule. We showed that hydrogen bonds between residues in surface loops can add to the overall stability of the V_L domains. The contribution to stability is further increased if the surface residue makes more than one hydrogen bond or if it forms a hydrogen bond between neighboring turns or loops separated from each other in the amino acid sequence. Based on our experiments we suggest that stabilization of proteins might be systematically accomplished by introducing additional hydrogen bonds on the surface. These substitutions are more straightforward to predict than core-packing interactions and can be selected to avoid affecting the protein’s function.     

Effect of varying polyglutamate chain length on the structure and stability of ferricytochrome c

The effect of varying polyglutamate chain length on local and global stability of horse heart ferricytochrome c was studied using scanning calorimetry and spectroscopy methods. Spectral data indicate that polyglutamate chain lengths equal or greater than eight monomer units significantly change the apparent pK(a) for the alkaline transition of cytochrome c. The change in pK(a) is comparable to the value when cytochrome c is complexed with cytochrome bc(1). Glutamate and diglutamate do not significantly alter the temperature transition for cleavage of the Met(80)-heme iron bond of cytochrome c. At low ionic strength, polyglutamates consisting of eight or more glutamate monomers increase midpoint of the temperature transition from 57.3+/-0.2 to 66.9+/-0.2 degrees C. On the other hand, the denaturation temperature of cytochrome c decreases from 85.2+/-0.2 to 68.8+/-0.2 degrees C in the presence of polyglutamates with number of glutamate monomers n >or approximately equal 8. The rate constant for cyanide binding to the heme iron of cytochrome c of cytochrome c-polyglutamate complex also decreases by approximately 42.5% with n>or approximately equal 8. The binding constant for the binding of octaglutamate (m.w. approximately 1000) to cyt c was found to be 1.15 x 10(5) M(-1) at pH 8.0 and low ionic strength. The results indicate that the polyglutamate (n>or approximately equal 8) is able to increase the stability of the methionine sulfur-heme iron bond of cytochrome c in spite of structural differences that weaken the overall stability of the cyt c at neutral and slightly alkaline pH.     

Folding-unfolding of goat alpha-lactalbumin studied by stopped-flow circular dichroism and molecular dynamics simulations

Folding reaction of goat alpha-lactalbumin has been studied by stopped-flow circular dichroism and molecular dynamics simulations. The effects of four single mutations and a double mutation on the stability of the protein under a native condition were studied. The mutations were introduced into residues located at a hydrophobic core in the alpha-domain of the molecule. Here we show that an amino acid substitution (T29I) increases the native-state stability of goat alpha-lactalbumin against the guanidine hydrochloride-induced unfolding by 3.5 kcal/mol. Kinetic refolding and unfolding of wild-type and mutant goat alpha-lactalbumin measured by stopped-flow circular dichroism showed that the local structure around the Thr29 side chain was not constructed in the transition state of the folding reaction. To characterize the local structural change around the Thr29 side chain to an atomic level of resolution, we performed high-temperature (at 400 K and 600 K) molecular dynamics simulations and studied the structural change at an initial stage of unfolding observed in the simulation trajectories. The Thr29 portion of the molecule experienced structural disruption accompanied with the loss of inter-residue contacts and with the water molecule penetration in the 400-K simulation as well as in four of the six 600-K simulations. Disruption of the N-terminal portion was also observed and was consistent with the results of kinetic refolding/unfolding experiments shown in our previous report.     

Cis proline mutants of ribonuclease A. I. Thermal stability.

A chemically synthesized gene for ribonuclease A has been expressed in Escherichia coli using a T7 expression system (Studier, F.W., Rosenberg, A.H., Dunn, J.J., & Dubendorff, J.W., 1990, Methods Enzymol. 185, 60-89). The expressed protein, which contains an additional N-terminal methionine residue, has physical and catalytic properties close to those of bovine ribonuclease A. The expressed protein accumulates in inclusion bodies and has scrambled disulfide bonds; the native disulfide bonds are regenerated during purification. Site-directed mutations have been made at each of the two cis proline residues, 93 and 114, and a double mutant has been made. In contrast to results reported for replacement of trans proline residues, replacement of either cis proline is strongly destabilizing. Thermal unfolding experiments on four single mutants give delta Tm approximately equal to 10 degrees C and delta delta G0 (apparent) = 2-3 kcal/mol. The reason is that either the substituted amino acid goes in cis, and cis<==>trans isomerization after unfolding pulls the unfolding equilibrium toward the unfolded state, or else there is a conformational change, which by itself is destabilizing relative to the wild-type conformation, that allows the substituted amino acid to form a trans peptide bond.     

Contributions of a hydrogen bond/salt bridge network to the stability of secondary and tertiary structure in lambda repressor.

In the N-terminal domain of lambda repressor, the Asp 14 side chain forms an intrahelical, hydrogen bond/salt bridge with the Arg 17 side chain and a tertiary hydrogen bond with the Ser 77 side chain. By measuring the stabilities to urea denaturation of the wild-type N-terminal domain and variants containing single, double, and triple alanine substitutions at positions 14, 17, and 77, the side-chain interaction energies, the coupling energy between interactions, and the intrinsic effects of each wild-type side chain on protein stability have been estimated. These studies indicate that the Asp 14-Arg 17 and Asp 14-Ser 77 interactions are stabilizing by roughly 0.8 and 1.5 kcal/mol, respectively, but that Asp 14, by itself, is destabilizing by roughly 0.9 kcal/mol. We also show that a peptide model of alpha-helix 1, which contains Asp 14 and Arg 17, forms a reasonably stable, monomeric helix in solution and responds to alanine mutations at positions 14 and 17 in the fashion expected from the intact protein studies. These studies suggest that it is possible to view the stability effects of mutations in intact proteins in a hierarchical fashion, with the stability of units of secondary structure being distinguishable from the stability of tertiary structure.     

A nuclear mutation prevents processing of a mitochondrially encoded membrane protein in Saccharomyces cerevisiae.

Subunit II of cytochrome oxidase is encoded by the mitochondrial OXI1 gene in Saccharomyces cerevisiae. The temperature-sensitive nuclear pet mutant ts2858 has an apparent higher mol. wt. subunit II when analyzed on lithium dodecylsulfate (LiDS) polyacrylamide gels. However, on LiDS-6M urea gels the apparent mol. wt. of the wild-type protein exceeds that of the mutant. Partial revertants of mutant ts2858 that produce both the wild-type and mutant form of subunit II were isolated. The two forms of subunit II differ at the N-terminal part of the molecule as shown by constructing and analyzing nuclear ts2858 and mitochondrial chain termination double mutants. The presence of the primary translation product in the mutant and of the processed form in the wild-type lacking 15 amino-terminal residues was demonstrated by radiolabel protein sequencing. Comparison of the known DNA sequence with the partial protein sequence obtained reveals that six of the 15 residues are hydrophilic and, unlike most signal sequences, this transient sequence does not contain extended hydrophobic parts. The nuclear mutation ts2858 preventing post-translational processing of cytochrome oxidase subunit II lies either in the gene for a protease or an enzyme regulating a protease.     

High-resolution structural and thermodynamic analysis of extreme stabilization of human procarboxypeptidase by computational protein design

Recent efforts to design de novo or redesign the sequence and structure of proteins using computational techniques have met with significant success. Most, if not all, of these computational methodologies attempt to model atomic-level interactions, and hence high-resolution structural characterization of the designed proteins is critical for evaluating the atomic-level accuracy of the underlying design force-fields. We previously used our computational protein design protocol RosettaDesign to completely redesign the sequence of the activation domain of human procarboxypeptidase A2. With 68% of the wild-type sequence changed, the designed protein, AYEdesign, is over 10 kcal/mol more stable than the wild-type protein. Here, we describe the high-resolution crystal structure and solution NMR structure of AYEdesign, which show that the experimentally determined backbone and side-chains conformations are effectively superimposable with the computational model at atomic resolution. To isolate the origins of the remarkable stabilization, we have designed and characterized a new series of procarboxypeptidase mutants that gain significant thermodynamic stability with a minimal number of mutations; one mutant gains more than 5 kcal/mol of stability over the wild-type protein with only four amino acid changes. We explore the relationship between force-field smoothing and conformational sampling by comparing the experimentally determined free energies of the overall design and these focused subsets of mutations to those predicted using modified force-fields, and both fixed and flexible backbone sampling protocols.     

Contribution of hydrogen bonds to the conformational stability of human lysozyme: calorimetry and X-ray analysis of six Ser --> Ala mutants.

To further examine the contribution of hydrogen bonds to the conformational stability of the human lysozyme, six Ser to Ala mutants were constructed. The thermodynamic parameters for denaturation of these six Ser mutant proteins were investigated by differential scanning calorimetry (DSC), and the crystal structures were determined by X-ray analysis. The denaturation Gibbs energy (DeltaG) of the Ser mutant proteins was changed from 2.0 to -5.7 kJ/mol, compared to that of the wild-type protein. With an analysis in which some factors that affected the stability due to mutation were considered, the contribution of hydrogen bonds to the stability (Delta DeltaGHB) was extracted on the basis of the structures of the mutant proteins. The results showed that hydrogen bonds between protein atoms and between a protein atom and a water bound with the protein molecule favorably contribute to the protein stability. The net contribution of one intramolecular hydrogen bond to protein stability (DeltaGHB) was 8.9 +/- 2.6 kJ/mol on average. However, the contribution to the protein stability of hydrogen bonds between a protein atom and a bound water molecule was smaller than that for a bond between protein atoms.     

Effect of the tyrosine 96 hydrogen bond on the inactivation of cytochrome P-450cam induced by hydrostatic pressure

The effects of removal of the tyrosine 96 hydrogen bond on the stability and conformational events of cytochrome P-450cam are presented in this communication. Hydrostatic pressure has been used as a tool to perturbe the structure leading to the formation of cytochrome P-420, an inactivated but soluble and undenatured form of the enzyme. We show that the spin transition of cytochrome P-450cam, which is known to be influenced by hydrostatic pressure, is affected by this single mutation. The free energy of stabilisation of native substrate-free cytochrome P-450cam is not affected by the removal of the tyrosine 96 hydrogen bond via mutagenesis to phenylalanine, whereas the substrate-bound protein shows a difference of 21 kJ/mol. These results, as well as an observed 110 ml/mol difference for the volume of the inactivation reaction between substrate-bound native and mutant proteins, have been interpreted in terms of a more hydrated heme pocket for the site-directed mutant at position 96 compared to the wild-type protein where camphor is tightly bound via the tyrosine 96 hydrogen bond and water excluded from the active site.     

Increasing the thermostability of staphylococcal nuclease: implications for the origin of protein thermostability

Seven hyper-stable multiple mutants have been constructed in staphylococcal nuclease by various combinations of eight different stabilizing single mutants. The stabilities of these multiple mutants determined by guanidine hydrochloride denaturation were 3.4 to 5.6 kcal/mol higher than that of the wild-type. Their thermal denaturation midpoint temperatures were 12.6 to 22.9 deg. C higher than that of the wild-type. These are among the greatest increases in protein stability and thermal denaturation midpoint temperature relative to the wild-type yet attained. There has been great interest in understanding how proteins found in thermophilic organisms are stabilized. One frequently cited theory is that the packing of hydrophobic side-chains is improved in the cores of proteins isolated from thermophiles when compared to proteins from mesophiles. The crystal structures of four single and five multiple stabilizing mutants of staphylococcal nuclease were solved to high resolution. No large overall structural change was found, with most changes localized around the sites of mutation. Rearrangements were observed in the packing of side-chains in the major hydrophobic core, although none of the mutations was in the core. It is surprising that detailed structural analysis showed that packing had improved, with the volume of the mutant protein's hydrophobic cores decreasing as protein stability increased. Further, the number of van der Waals interactions in the entire protein showed an experimentally significant increase correlated with increasing stability. These results indicate that optimization of packing follows as a natural consequence of increased protein thermostability and that good packing is not necessarily the proximate cause of high stability. Another popular theory is that thermostable proteins have more electrostatic and hydrogen bonding interactions and these are responsible for the high stabilities. The mutants here show that increased numbers of electrostatic and hydrogen bonding interactions are not obligatory for large increases in protein stability.     

Conformational stability and activity of ribonuclease T1 and mutants. Gln25----Lys, Glu58----Ala, and the double mutant

Ribonuclease T1 (RNase T1) and mutants Gln25----Lys, Glu58----Ala, and the double mutant were prepared from a chemically synthesized gene, cloned and expressed in Escherichia coli. The wild-type RNase T1 prepared from the cloned gene was identical in every functional and physical property examined to RNase T1 prepared from Aspergillus oryzae. Urea and thermal unfolding experiments show that Gln25----Lys is 0.9 kcal/mol more stable and Glu58----Ala is 0.8 kcal/mol less stable than wild-type RNase T1. In the double mutant, these contributions cancel and the stability does not differ significantly from that of wild-type RNase T1. For the double mutant, the dependence of delta G on urea concentration is significantly greater than for wild-type RNase T1 or the single mutants. This suggests that the double mutant unfolds more completely in urea than the other proteins. The activity of Gln25----Lys is identical with that of wild-type RNase T1. The activities of Glu58----Ala and the double mutant are 7% of wild-type when GpC hydrolysis is measured (due to a 35-fold decrease in kcat), and 37% of wild-type when RNA hydrolysis is measured. Thus, Glu58 is important, but not essential to the activity of RNase T1.     

Role of hydrogen bond networks and dynamics in positive and negative cooperative stabilization of a protein

Cooperativity mediated through hydrogen bond networks in yeast iso-1-cytochrome c was studied using a thermodynamic triple mutant cycle. Three known stabilizing mutations, Asn 26 to His, Asn 52 to Ile, and Tyr 67 to Phe, were used to construct the triple mutant cycle. The side chain of His 26, a wild-type residue, forms two hydrogen bonds that bridge two substructures of the wild-type protein, and Tyr 67 and Asn 52 are part of an extensive buried hydrogen bond network. The stabilities of all variants in the triple mutant cycle were determined by guanidine hydrochloride denaturation methods and used to determine the pairwise, Delta(2)G(int), and triple interaction energies. His 26 and Ile 52 interact cooperatively (Delta(2)G(int) is 1-2 kcal/mol), whereas the two other pairs of mutations interact anticooperatively (Delta(2)G(int) is -0.5 to -1.5 kcal/mol). Previously reported structural data for iso-1-cytochrome c variants containing these mutations show that changes in the strength of the His 26 to Glu 44 hydrogen bond, apparently caused by changes in main chain dynamics, provide a mechanism for the long distance (His 26 to Phe 67 and His 26 to Ile 52) propagation of pairwise interaction energies. Opposing changes in the strength of the His 26 to Glu 44 hydrogen bond caused by the N52I and Y67F mutations generate a negative triple interaction energy (-0.9 +/-0.7 kcal/mol) that combined with cancellation of cooperative and anticooperative pairwise interactions produce apparent additivity for the stabilizing effects of the single mutations in the triple mutant variant.     

Structure and stability effects of the mutation of glycine 34 to serine in Rhodobacter capsulatus cytochrome c(2)

Gly 34 and the adjacent Pro 35 of Rhodobacter capsulatus cytochrome c(2) (or Gly 29 and Pro 30 in vertebrate cytochrome c) are highly conserved side chains among the class I c-type cytochromes. The mutation of Gly 34 to Ser in Rb. capsulatus cytochrome c(2) has been characterized in terms of physicochemical properties and NMR in both redox states. A comparison of the wild-type cytochrome c(2), the G34S mutation, and the P35A mutation is presented in the context of differences in chemical shifts, the differences in NOE patterns, and structural changes resulting from oxidation of the reduced cytochrome. G34S is substantially destabilized relative to wild-type (2.2 kcal/mol in the oxidized state) but similarly destabilized relative to P35A. Nevertheless, differences in terms of the impact of the mutations on specific structural regions are found when comparing G34S and P35A. Although available data indicates that the overall secondary structure of G34S and wild-type cytochrome c(2) are similar, a number of both perturbations of hydrogen bond networks and interactions with internal waters are found. Thus, the impact of the mutation at position 35 is propagated throughout the cytochrome but with alterations at defined sites within the molecule. Interestingly, we find that the substitution of serine at position 34 results in a perturbation of the heme beta meso and the methyl-5 protons. This suggests that the hydroxyl and beta carbon are positioned away from the solvent and toward the heme. This has the consequence of preferentially stabilizing the oxidized state in G34S, thus, altering hydrogen bond networks which involve the heme propionate, internal waters, and key amino acid side chains. The results presented provide important new insights into the stability and solution structure of the cytochrome c(2).     

A single disulfide bond restores thermodynamic and proteolytic stability to an extensively mutated protein

The potential for engineering stable proteins with multiple amino acid substitutions was explored. Eleven lysine, five methionine, two tryptophan, one glycine, and three threonine substitutions were simultaneously made in barley chymotrypsin inhibitor-2 (CI-2) to substantially improve the essential amino acid content of the protein. These substitutions were chosen based on the three-dimensional structure of CI-2 and an alignment of homologous sequences. The initial engineered protein folded into a wild-type-like structure, but had a free energy of unfolding of only 2.2 kcal/mol, considerably less than the wild-type value of 7.5 kcal/mol. Restoration of the lysine mutation at position 67 to the wild-type arginine increased the free energy of unfolding to 3.1 kcal/mol. Subsequent cysteine substitutions at positions 22 and 82 resulted in disulfide bond formation and a protein with nearly wild-type thermodynamic stability (7.0 kcal/mol). None of the engineered proteins retained inhibitory activity against chymotrypsin or elastase, and all had substantially reduced inhibitory activity against subtilisin. The proteolytic stabilities of the proteins correlated with their thermodynamic stabilities. Reduction of the disulfide bond resulted in substantial loss of both thermodynamic and proteolytic stabilities, confirming that the disulfide bond, and not merely the cysteine substitutions, was responsible for the increased stability. We conclude that it is possible to replace over a third of the residues in CI-2 with minimal disruption of stability and structural integrity.     

Evaluation of cooperative interactions between substructures of iso-1-cytochrome c using double mutant cycles

A double mutant cycle has been used to evaluate interaction energies between the global stabilizer mutation asparagine 52 --> isoleucine (N52I) in iso-1-cytochrome c and mutations producing single surface histidines at positions 26, 33, 39, 54, 73, 89, and 100. These histidine mutation sites are distributed through the four cooperative folding units of cytochrome c. The double mutant cycle starts with the iso-1-cytochrome c variant AcTM, a variant with no surface histidines and with asparagine at position 52. Isoleucine is added singly at position 52, AcTMI52 variant, as are the surface histidines, AcHX variants, where X indicates the histidine sequence position. The double mutant variants, AcHXI52, provide the remaining corner of the double mutant cycle. The stabilities of all variants were determined by guanidine hydrochloride denaturation and interaction energies were calculated between position 52 and each histidine site. Six of the seven double mutants show additive (AcH33I52, AcH39I52, AcH54I52, AcH89I52, and AcH100I52) stability effects or weak interaction energies (AcH73I52) of the histidine mutations and the N52I mutation, consistent with cooperative effects on protein folding and stability being sparsely distributed through the protein structure. The AcH26I52 variant shows a strong favorable interaction energy, 2.0 +/- 0.5 kcal/mol, between the N52I mutation in one substructure and the addition of His 26 to an adjacent substructure. The data are consistent with an entropic stabilization of the intersubstructure hydrogen bond between His 26 and Glu 44 by the Ile 52 mutation.     

Filling small, empty protein cavities: structural and energetic consequences

Most proteins contain small cavities that can be filled by replacing cavity-lining residues by larger ones. Since shortening mutations in hydrophobic cores tend to destabilize proteins, it is expected that cavity-filling mutations may conversely increase protein stability. We have filled three small cavities in apoflavodoxin and determined by NMR and equilibrium unfolding analysis their impact in protein structure and stability. The smallest cavity (14 A3) has been filled, at two different positions, with a variety of residues and, in all cases, the mutant proteins are locally unfolded, their structure and energetics resembling those of an equilibrium intermediate of the thermal unfolding of the wild-type protein. In contrast, two slightly larger cavities of 20 A3 and 21 A3 have been filled with Val to Ile or Val to Leu mutations and the mutants preserve both the native fold and the equilibrium unfolding mechanism. From the known relationship, observed in shortening mutations, between stability changes and the differential hydrophobicity of the exchanged residues and the volume of the cavities, the filling of these apoflavodoxin cavities is expected to stabilize the protein by approximately 1.5 kcal mol(-1). However, both urea and thermal denaturation analysis reveal much more modest stabilizations, ranging from 0.0 kcal mol(-1) to 0.6 kcal mol(-1), which reflects that the accommodation of single extra methyl groups in small cavities requires some rearrangement, necessarily destabilizing, that lowers the expected theoretical stabilization. As the size of these cavities is representative of that of the typical small, empty cavities found in most proteins, it seems unlikely that filling this type of cavities will give rise to large stabilizations.     

Flexible-geometry conformational energy maps for the amino acid residue preceding a proline.

Previously calculated conformational energy maps suggest that the alpha-helical conformation for the residue preceding a proline is disfavored relative to the extended conformation by more than 7 kcal/mol. In known protein structures this conformation is observed, however, to occur for about 9% of all prolines. In addition, introduction or removal of prolines at theoretically unfavorable positions in proteins and peptides can have modest effects on stability and structure. To investigate the discrepancy between calculation and experiment, we have determined how the conformation of the proline affects the calculated energy. We have also explored the effect of bond length and bond angle relaxation on the conformational energy map. The conformational energy of the preceding residue is found to be unaffected by the conformation of the proline, but the effect of allowing covalent bond relaxation is dramatic. If bond lengths and angles, and dihedral angles within the pyrrolidine ring, are allowed to relax, a calculated energy difference between the alpha and beta conformations of 1.1 kcal/mol is obtained, in reasonable agreement with experiment. The detailed shape of the calculated energy surface is also in excellent agreement with the observed conformational distributions in known protein structures.     

Tolerance of point substitution of methionine for isoleucine in hen egg white lysozyme

X-ray structure determination of proteins by using the multiple-wavelength anomalous dispersion method targeting selenomethionine is now widely employed. Isoleucine was examined for the second choice of the substitution of methionine next to leucine. We performed a systematic mutational study of the substitutions of methionine for isoleucine. All mutated lysozymes were less stable than the wild-type by about 1 kcal/mol and it is suggested that this instability was caused by the change in residual hydrophobicity from isoleucine to methionine. The X-ray structures of all mutant lysozymes were very similar to that of the wild-type. In addition, both the accessible surface areas and the conformation of the side chain of methionine in all mutant lysozymes were similar to those of the side chain at the respective isoleucine in the wild-type. Therefore, it is suggested that the mutation from isoleucine to methionine in a protein can be considered as a "safe" substitution.