Relationship between local structure and stability in hen egg white lysozyme mutant with alanine substituted for glycine

An amino-domino model described by a cross-peptide-bond Ramachandran plot defines amino acid pairs as local structural units

Protein structure, both at the global and local level, dictates function. Proteins fold from chains of amino acids, forming secondary structures, α-helices and β-strands, that, at least for globular proteins, subsequently fold into a three-dimensional structure. Here, we show that a Ramachandran-type plot focusing on the two dihedral angles separated by the peptide bond, and entirely contained within an amino acid pair, defines a local structural unit. We further demonstrate the usefulness of this cross-peptide-bond Ramachandran plot by showing that it captures β-turn conformations in coil regions, that traditional Ramachandran plot outliers fall into occupied regions of our plot, and that thermophilic proteins prefer specific amino acid pair conformations. Further, we demonstrate experimentally that the effect of a point mutation on backbone conformation and protein stability depends on the amino acid pair context, i.e., the identity of the adjacent amino acid, in a manner predictable by our method.

Glycine Substitution of Residues with Unfavored Dihedral Angles Improves Protein Thermostability

Single mutations that can substantially enhance stability are highly desirable for protein engineering. However, it is generally rare for this kind of mutant to emerge from directed evolution experiments. This study used computational approaches to identify hotspots in a diacylglycerol-specific lipase for mutagenesis with functional hotspot and sequence consensus strategies, followed by ∆∆G calculations for all possible mutations using the Rosetta ddg_monomer protocol. Single mutants with significant ∆∆G changes (≤−2.5 kcal/mol) were selected for expression and characterization. Three out of seven tested mutants showed a significantly enhanced thermostability, with Q282W and A292G in the catalytic pocket and D245G located on the opposite surface of the protein. Remarkably, A292G increased the T5015 (the temperature at which 50% of the enzyme activity was lost after a 15 min of incubation) by ~7 °C, concomitant with a twofold increase in enzymatic activity at the optimal reaction temperature. Structural analysis showed that both A292 and D245 adopted unfavored dihedral angles in the wild-type (WT) enzyme. Substitution of them by glycine might release a steric strain to increase the stability. In sum, substitution by glycine might be a promising strategy to improve protein thermostability.

Structural Analysis of Hen Egg Lysozyme Refolded after Denaturation at Acidic pH

Protein structures fluctuate in solution; therefore, proteins have multiple stable structures that are slightly different from each other. In this study, we determined the crystal structure of hen egg lysozyme refolded after denaturation at acidic pH (rHEL) and found a structure different from native HEL (nHEL). The different local conformations of the peptide bond between Asp101 and Gly102 found in the crystal structure was supported by the NMR results for nHEL and rHEL. The NMR experiments also showed shifts in the heteronuclear single quantum coherence signals derived from Thr43 and Asp52. The chemical shift change of Asp52 could be explained by the crystal structure of rHEL, showing the conformational change of Tyr53, whose phenol ring directly lies on the main chain of Asp52. The catalytic activity of rHEL was similar to that of nHEL, indicating that the conformational change had little effect on activity. In contrast, conformational changes could be detected by the binding of monoclonal antibodies against HEL. Using multiple methods, we successfully detected the unusual structure of HEL, which might be another stable structure of HEL in solution.

Integrated proteomic, phosphoproteomic and N-glycoproteomic analyses of chicken eggshell matrix

Eggshell matrix (EM) proteins play an important biological role in eggshell mineralization and embryo development. Many studies have demonstrated that some matrix proteins undergo posttranslational modifications, including phosphorylation and glycosylation, which have important regulatory effects on the functional properties of the proteins. Systematic analysis of the proteome, the phosphorylated modified proteome and the glycosylated modified proteome of the chicken EM was performed using a proteomics strategy. A total of 112 phosphorylation sites from 69 phosphoproteins and 297 N-glycosylation sites from 182 N-glycoproteins were identified in the chicken EM. Among all these identified modified proteins, 129 were not identified in the proteome (547 proteins). Therefore, a total of 676 EM proteins were identified in this study. Gene ontology (GO) enrichment analysis indicated that EM proteins and phosphoproteins were mainly enriched in regulation of enzyme activity, while EM N-glycoproteins were enriched in immune response regulation.

Protein crystal structure from non-oriented, single-axis sparse X-ray data

X-ray free-electron lasers (XFELs) have inspired the development of serial femtosecond crystallography (SFX) as a method to solve the structure of proteins. SFX datasets are collected from a sequence of protein microcrystals injected across ultrashort X-ray pulses. The idea behind SFX is that diffraction from the intense, ultrashort X-ray pulses leaves the crystal before the crystal is obliterated by the effects of the X-ray pulse. The success of SFX at XFELs has catalyzed interest in analogous experiments at synchrotron-radiation (SR) sources, where data are collected from many small crystals and the ultrashort pulses are replaced by exposure times that are kept short enough to avoid significant crystal damage. The diffraction signal from each short exposure is so 'sparse' in recorded photons that the process of recording the crystal intensity is itself a reconstruction problem. Using the EMC algorithm, a successful reconstruction is demonstrated here in a sparsity regime where there are no Bragg peaks that conventionally would serve to determine the orientation of the crystal in each exposure. In this proof-of-principle experiment, a hen egg-white lysozyme (HEWL) crystal rotating about a single axis was illuminated by an X-ray beam from an X-ray generator to simulate the diffraction patterns of microcrystals from synchrotron radiation. Millions of these sparse frames, typically containing only ∼200 photons per frame, were recorded using a fast-framing detector. It is shown that reconstruction of three-dimensional diffraction intensity is possible using the EMC algorithm, even with these extremely sparse frames and without knowledge of the rotation angle. Further, the reconstructed intensity can be phased and refined to solve the protein structure using traditional crystallographic software. This suggests that synchrotron-based serial crystallography of micrometre-sized crystals can be practical with the aid of the EMC algorithm even in cases where the data are sparse.

Optimal definition of inter-residual contact in globular proteins based on pairwise interaction energy calculations, its robustness, and applications

Although a contact is an essential measurement for the topology as well as strength of non-covalent interactions in biomolecules and their complexes, there is no general agreement in the definition of this feature. Most of the definitions work with simple geometric criteria which do not fully reflect the energy content or ability of the biomolecular building blocks to arrange their environment. We offer a reasonable solution to this problem by distinguishing between "productive" and "non-productive" contacts based on their interaction energy strength and properties. We have proposed a method which converts the protein topology into a contact map that represents interactions with statistically significant high interaction energies. We do not prove that these contacts are exclusively stabilizing, but they represent a gateway to thermodynamically important rather than geometry-based contacts. The process is based on protein fragmentation and calculation of interaction energies using the OPLS force field and relies on pairwise additivity of amino acid interactions. Our approach integrates the treatment of different types of interactions, avoiding the problems resulting from different contributions to the overall stability and the different effect of the environment. The first applications on a set of homologous proteins have shown the usefulness of this classification for a sound estimate of protein stability.

Increasing protein stability by improving beta-turns

Our goal was to gain a better understanding of how protein stability can be increased by improving beta-turns. We studied 22 beta-turns in nine proteins with 66-370 residues by replacing other residues with proline and glycine and measuring the stability. These two residues are statistically preferred in some beta-turn positions. We studied: Cold shock protein B (CspB), Histidine-containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase alpha-subunit, and Maltose binding protein. Of the 15 single proline mutations, 11 increased stability (Average = 0.8 +/- 0.3; Range = 0.3-1.5 kcal/mol), and the stabilizing effect of double proline mutants was additive. On the basis of this and our previous work, we conclude that proteins can generally be stabilized by replacing nonproline residues with proline residues at the i + 1 position of Type I and II beta-turns and at the i position in Type II beta-turns. Other turn positions can sometimes be used if the phi angle is near -60 degrees for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in beta-turns in related proteins give similar results. Proline substitutions increase stability mainly by decreasing the entropy of the denatured state. In contrast, the large, diverse group of proteins considered here had almost no residues in beta-turns that could be replaced by Gly to increase protein stability. Improving beta-turns by substituting Pro residues is a generally useful way of increasing protein stability.

Crystal structures of K33 mutant hen lysozymes with enhanced activities

Using random mutagenesis, we previously obtained K33N mutant lysozyme that showed a large lytic halo on the plate coating Micrococcus luteus. In order to examine the effects of mutation of K33N on enzyme activity, we prepared K33N and K33A mutant lysozymes from yeast. It was found that the activities of both the mutant lysozymes were higher than those of the wild-type lysozyme based on the results of the activity measurements against M. luteus (lytic activity) and glycol chitin. Moreover, 3D structures of K33N and K33A mutant lysozyme were solved by X-ray crystallographic analyses. The side chain of K33 in the wild-type lysozyme hydrogen bonded with N37 involved in the substrate-binding region, and the orientation of the side chain of N37 in K33 mutant lysozymes were different in the wild-type lysozyme. These results suggest that the enhancement of activity in K33N mutant lysozyme was due to an alteration in the orientation of the side chain of N37. On the other hand, K33N lysozyme was less stable than the wild-type lysozyme. Lysozyme may sacrifice its enzyme activity to acquire the conformational stability at position 33.

Mutational analysis of conserved glycine residues 142, 143 and 146 reveals Gly(142) is critical for tetramerization of CTP synthase from Escherichia coli

CTPS (cytidine 5'-triphosphate synthase) catalyses the ATP-dependent formation of CTP from UTP using either ammonia or L-glutamine as the nitrogen source. Binding of the substrates ATP and UTP, or the product CTP, promotes oligomerization of CTPS from inactive dimers to active tetramers. In the present study, site-directed mutagenesis was used to replace the fully conserved glycine residues 142 and 143 within the UTP-binding site and 146 within the CTP-binding site of Escherchia coli CTPS. CD spectral analyses of wild-type CTPS and the glycine mutants showed a slight reduction of approximately 15% in alpha-helical content for G142A and G143A relative to G146A and wild-type CTPS, suggesting some local alterations in structure. Relative to wild-type CTPS, the values of k(cat)/K(m) for ammonia-dependent and glutamine-dependent CTP formation catalysed by G143A were reduced 22- and 16-fold respectively, whereas the corresponding values for G146A were reduced only 1.4- and 1.8-fold respectively. The glutaminase activity (k(cat)) of G146A was similar to that exhibited by the wild-type enzyme, whereas that of G143A was reduced 7.5-fold. G146A exhibited substrate inhibition at high concentrations of ammonia and a partial uncoupling of glutamine hydrolysis from CTP production. Although the apparent affinity (1/[S](0.5)) of G143A and G146A for UTP was reduced approximately 4-fold, G146A exhibited increased co-operativity with respect to UTP. Thus mutations in the CTP-binding site can affect UTP-dependent activity. Surprisingly, G142A was inactive with both ammonia and glutamine as substrates. Gel-filtration HPLC experiments revealed that both G143A and G146A were able to form active tetramers in the presence of ATP and UTP; however, nucleotide-dependent tetramerization of G142A was significantly impaired. Our observations highlight the sensitivity of the structure of CTPS to mutations in the UTP- and CTP-binding sites, with Gly(142) being critical for nucleotide-dependent oligomerization of CTPS to active tetramers. This 'structural sensitivity' may limit the number and/or types of mutations that could be selected for during the development of resistance to cytotoxic pyrimidine nucleotide analogues.

Increasing protein conformational stability by optimizing beta-turn sequence

Protein conformational stability is an important concern in many fields. Here, we describe a strategy for significantly increasing conformational stability by optimizing beta-turn sequence. Proline and glycine residues are statistically preferred at several beta-turn positions, presumably because their unique side-chains contribute favorably to conformational stability in certain beta-turn positions. However, beta-turn sequences often deviate from preferred proline or preferred glycine. Therefore, our strategy involves replacing non-proline and non-glycine beta-turn residues with preferred proline or preferred glycine residues. Here, we develop guidelines for selecting appropriate mutations, and present results for five mutations (S31P, S42G, S48P, T76P, and Q77G) that significantly increase the conformational stability of RNase Sa. The increases in stability ranged from 0.7 kcal/mol to 1.3 kcal/mol. The strategy was successful in overlapping or isolated beta-turns, at buried (up to 50%) or completely exposed sites, and at relatively flexible or inflexible sites. Considering the significant number of beta-turn residues in every globular protein and the frequent deviation of beta-turn sequences from preferred proline and preferred glycine residues, this simple, efficient strategy will be useful for increasing the conformational stability of proteins.

Average assignment method for predicting the stability of protein mutants

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and it will be helpful for designing stable mutants. In this work, we have analyzed the stability of protein mutants using three different data sets of 1791, 1396, and 2204 mutants, respectively, for thermal stability (DeltaTm), free energy change due to thermal (DeltaDeltaG), and denaturant denaturations (DeltaDeltaGH2O), obtained from the ProTherm database. We have classified the mutants into 380 possible substitutions and assigned the stability of each mutant using the information obtained with similar type of mutations. We observed that this assignment could distinguish the stabilizing and destabilizing mutants to an accuracy of 70-80% at different measures of stability. Further, we have classified the mutants based on secondary structure and solvent accessibility (ASA) and observed that the classification significantly improved the accuracy of prediction. The classification of mutants based on helix, strand, and coil distinguished the stabilizing/destabilizing mutants at an average accuracy of 82% and the correlation is 0.56; information about the location of residues at the interior, partially buried, and surface regions of a protein correctly identified the stabilizing/destabilizing residues at an average accuracy of 81% and the correlation is 0.59. The nine subclassifications based on three secondary structures and solvent accessibilities improved the accuracy of assigning stabilizing/destabilizing mutants to an accuracy of 84-89% for the three data sets. Further, the present method is able to predict the free energy change (DeltaDeltaG) upon mutations within a deviation of 0.64 kcal/mol. We suggest that this method could be used for predicting the stability of protein mutants.

Calorimetric Analysis of Thermodynamic Stability and Aggregation for Apo and Holo Amyotrophic Lateral Sclerosis-associated Gly-93 Mutants of Superoxide Dismutase

Differential scanning calorimetry was used to measure changes in thermodynamic stability and aggregation for glycine 93 mutants of human copper, zinc-superoxide dismutase (SOD). Glycine 93 is a conserved residue at position i + 3 of a tight turn and has been found to be a mutational hot spot in familial amyotrophic lateral sclerosis (fALS). The fALS-associated mutations, G93A, G93S, G93R, G93D, and G93V, were made in a pseudo wild-type background containing no free cysteines, which prevented the formation of aberrant disulfide bonds upon thermal unfolding, and enabled quantitative thermodynamic analysis of the effects of the mutations. Thermal unfolding was highly reversible for all the SODs in both the fully metallated (holo) and metal-free (apo) forms. The data for all the holo-SODs and for the apo-pseudo-wild-type SOD were well fit by a 2-state unfolding model for native dimer (N2) to two unfolded monomers (2U), N2 <--> 2U. The holo- and apo-forms of the mutants are significantly destabilized (by 1.5-3.5 kcal mol(-1) monomer) relative to the corresponding forms of pseudo wild-type, with the relative stabilities being correlated with statistical preferences for amino acids in this structural context. Although van't Hoff (DeltaHvH) to calorimetric (DeltaHcal) enthalpy ratios are close to unity for all the holo-SODs and for apo-pseudo-wild-type, consistent with a 2-state transition, DeltaHvH is considerably larger than DeltaHcal for all the apo-mutants. This suggests that the mutations cause apo-SOD to have an increased propensity to misfold or aggregate, which may be linked to increased toxic mutant SOD aggregation in fALS.

Identification of a key structural element for protein folding within beta-hairpin turns

Specific residues in a polypeptide may be key contributors to the stability and foldability of the unique native structure. Identification and prediction of such residues is, therefore, an important area of investigation in solving the protein folding problem. Atypical main-chain conformations can help identify strains within a folded protein, and by inference, positions where unique amino acids may have a naturally high frequency of occurrence due to favorable contributions to stability and folding. Non-Gly residues located near the left-handed alpha-helical region (L-alpha) of the Ramachandran plot are a potential indicator of structural strain. Although many investigators have studied mutations at such positions, no consistent energetic or kinetic contributions to stability or folding have been elucidated. Here we report a study of the effects of Gly, Ala and Asn substitutions found within the L-alpha region at a characteristic position in defined beta-hairpin turns within human acidic fibroblast growth factor, and demonstrate consistent effects upon stability and folding kinetics. The thermodynamic and kinetic data are compared to available data for similar mutations in other proteins, with excellent agreement. The results have identified that Gly at the i+3 position within a subset of beta-hairpin turns is a key contributor towards increasing the rate of folding to the native state of the polypeptide while leaving the rate of unfolding largely unchanged.

Designed to be stable: crystal structure of a consensus ankyrin repeat protein

Ankyrin repeat (AR) proteins mediate innumerable protein-protein interactions in virtually all phyla. This finding suggested the use of AR proteins as designed binding molecules. Based on sequence and structural analyses, we designed a consensus AR with fixed framework and randomized interacting residues. We generated several combinatorial libraries of AR proteins consisting of defined numbers of this repeat. Randomly chosen library members are expressed in soluble form in the cytoplasm of Escherichia coli constituting up to 30% of total cellular protein and show high thermodynamic stability. We determined the crystal structure of one of those library members to 2.0-A resolution, providing insight into the consensus AR fold. Besides the highly complementary hydrophobic repeat-repeat interfaces and the absence of structural irregularities in the consensus AR protein, the regular and extended hydrogen bond networks in the beta-turn and loop regions are noteworthy. Furthermore, all residues found in the turn region of the Ramachandran plot are glycines. Many of these features also occur in natural AR proteins, but not in this rigorous and standardized fashion. We conclude that the AR domain fold is an intrinsically very stable and well-expressed scaffold, able to display randomized interacting residues. This scaffold represents an excellent basis for the design of novel binding molecules.

Exploring the role of a glycine cluster in cold adaptation of an alkaline phosphatase

In an effort to explore the role of glycine clusters on the cold adaptation of enzymes, we designed point mutations aiming to alter the distribution of glycine residues close to the active site of the psychrophilic alkaline phosphatase from the Antarctic strain TAB5. The mutagenesis targets were residues Gly261 and Gly262. The replacement of Gly262 by Ala resulted in an inactive enzyme. Substitution of Gly261 by Ala resulted to an enzyme with lower stability and increased energy of activation. The double mutant G261A/Y269A designed on the basis of side-chain packing criteria from a modelled structure of the enzyme resulted in restoration of the energy of activation to the levels of the native enzyme and in an increased stability compared to the mutant G261A. It seems therefore, that the Gly cluster in combination with its structural environment plays a significant role in the cold adaptation of the enzyme.

Role of amino acid residues in left-handed helical conformation for the conformational stability of a protein

Our previous study of six non-Gly to Gly/Ala mutant human lysozymes in a left-handed helical region showed that only one non-Gly residue at a rigid site had unfavorable strain energy as compared with Gly at the same position (Takano et al., Proteins 2001; 44:233-243). To further examine the role of left-handed residues in the conformational stability of a protein, we constructed ten Gly to Ala mutant human lysozymes. Most Gly residues in human lysozyme are located in the left-handed helix region. The thermodynamic parameters for denaturation and crystal structures were determined by differential scanning calorimetry and X-ray analysis, respectively. The difference in denaturation Gibbs energy (DeltaDeltaG) for the ten Gly to Ala mutants ranged from + 1.9 to -7.5 kJ/mol, indicating that the effect of the mutation depends on the environment of the residue. We confirm that Gly in a left-handed region is more favorable at rigid sites than non-Gly, but there is little difference in energetic cost between Gly and non-Gly at flexible sites. The present results indicate that dihedral angles in the backbone conformation and also the flexibility at the position should be considered for analyses of protein stability, and protein structural determination, prediction, and design.

Hydration of the peptide backbone largely defines the thermodynamic propensity scale of residues at the C' position of the C-capping box of alpha-helices

The C' position of the C-capping box is the second residue outside of the helix. Statistical analysis of residue distribution at the C' position in the alpha-helices' C-capping box showed that different amino acid residues occur with different probabilities, with the strongest preference being for glycine. To understand the physico-chemical basis for this preference, we studied the effects that 17 amino acid substitutions at the C' position in an alpha-helix of ubiquitin have on the stability of this protein. We determined the following rank order of amino acid residues at the C' position with respect to their effect on the stability: Gly>His>Asn>Arg>Lys>Gln>Ala>Phe>Met>Ser>Asp>Glu>Trp>Thr>Pro>Ile>Val. The effect of the amino acid substitutions on the structure also was evaluated by comparing the (1)H-(15)N heteronuclear sequential quantum correlation spectra and showed no significant changes in the structures of the most stable (Gly) and the least stable (Val) variants. The obtained changes in stability highly correlate (r = 0.85) with the statistical distribution of the residues at the C' position indicating that the measured thermodynamic propensities are unbiased by secondary interactions. We also found that the measured thermodynamic propensities correlate well with the amide hydrogen exchange data on short model peptides (r = 0.85) and the calculated hydration of the peptide backbone (r = 0.88). These results combined with the changes in enthalpy and entropy of unfolding of ubiquitin variants suggest that dehydration of the peptide backbone plays a significant role in defining the thermodynamic propensity scale at the C' position of the C-capping box in alpha-helices. This propensity scale is useful for protein secondary structure predictions and protein design.