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

Analyses of displacements resulting from a point mutation in proteins

The effects of a single residue substitution on the protein backbone are frequently quite small and there are many other potential sources of structural variation for protein. We present here a methodology considering different sources of distortions in order to isolate the very effect of the mutation. To validate our methodology, we consider a well-studied family with many single mutants: the human lysozyme. Most of the perturbations are expected to be at the very localisation of the mutation, but in many cases the effects are propagated at long range. We show that the distances between the mutated residue and the 5% most disturbed residues exponentially decreases. One third of the affected residues are in direct contact with the mutated position; the remaining two thirds are potential allosteric effects. We confirm the reliability of the residues identified as significantly perturbed by comparing our results to experimental studies. We confirm with the present method all the previously identified perturbations. This study shows that mutations have long-range impact on protein backbone that can be detected, although the displacement of the affected atoms is small.

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.

Effect of the disease-causing mutations identified in human ribonuclease (RNase) H2 on the activities and stabilities of yeast RNase H2 and archaeal RNase HII

Eukaryotic ribonuclease (RNase) H2 consists of one catalytic and two accessory subunits. Several single mutations in any one of these subunits of human RNase H2 cause Aicardi-Goutières syndrome. To examine whether these mutations affect the complex stability and activity of RNase H2, three mutant proteins of His-tagged Saccharomyces cerevisiae RNase H2 (Sc-RNase H2*) were constructed. Sc-G42S*, Sc-L52R*, and Sc-K46W* contain single mutations in Sc-Rnh2Ap*, Sc-Rnh2Bp*, and Sc-Rnh2Cp*, respectively. The genes encoding the three subunits were coexpressed in Escherichia coli, and Sc-RNase H2* and its derivatives were purified in a heterotrimeric form. All of these mutant proteins exhibited enzymatic activity. However, only the enzymatic activity of Sc-G42S* was greatly reduced compared to that of the wild-type protein. Gly42 is conserved as Gly10 in Thermococcus kodakareansis RNase HII. To analyze the role of this residue, four mutant proteins, Tk-G10S, Tk-G10A, Tk-G10L, and Tk-G10P, were constructed. All mutant proteins were less stable than the wild-type protein by 2.9-7.6 degrees C in T(m). A comparison of their enzymatic activities, substrate binding affinities, and CD spectra suggests that the introduction of a bulky side chain into this position induces a local conformational change, which is unfavorable for both activity and substrate binding. These results indicate that Gly10 is required to make the protein fully active and stable.

Structural assessment of glycyl mutations in invariantly conserved motifs

Motifs that are evolutionarily conserved in proteins are crucial to their structure and function. In one of our earlier studies, we demonstrated that the conserved motifs occurring invariantly across several organisms could act as structural determinants of the proteins. We observed the abundance of glycyl residues in these invariantly conserved motifs. The role of glycyl residues in highly conserved motifs has not been studied extensively. Thus, it would be interesting to examine the structural perturbations induced by mutation in these conserved glycyl sites. In this work, we selected a representative set of invariant signature (IS) peptides for which both the PDB structure and mutation information was available. We thoroughly analyzed the conformational features of the glycyl sites and their local interactions with the surrounding residues. Using Ramachandran angles, we showed that the glycyl residues occurring in these IS peptides, which have undergone mutation, occurred more often in the L-disallowed as compared with the L-allowed region of the Ramachandran plot. Short range contacts around the mutation site were analyzed to study the steric effects. With the results obtained from our analysis, we hypothesize that any change of activity arising because of such mutations must be attributed to the long-range interaction(s) of the new residue if the glycyl residue in the IS peptide occurred in the L-allowed region of the Ramachandran plot. However, the mutation of those conserved glycyl residues that occurred in the L-disallowed region of the Ramachandran plot might lead to an altered activity of the protein as a result of an altered conformation of the backbone in the immediate vicinity of the glycyl residue, in addition to long range effects arising from the long side chains of the new residue. Thus, the loss of activity because of mutation in the conserved glycyl site might either relate to long range interactions or to local perturbations around the site depending upon the conformational preference of the glycyl residue.

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.

Ionic Peptide Aggregation: Exploration of Conformational Dynamics in Aqueous Solution by Computational Techniques

The effects of end groups on KEK peptide conformational characteristics and self-assembling properties in water solution are investigated by using long lasting all-atom molecular dynamics simulations. The analysis of the structural macroscopic and microscopic properties and the examination of intra- and intermolecular interactions suggest, in agreement with experimental observations, the role played by side chains and terminal regions in determining the characteristic features of the assemblages. Competition between intra- and interchain interactions greatly affects the diffusivity of peptide molecules and the conformational space that they can sample, ultimately controlling the shape, size, and distribution of the aggregate configurations. Different peptide end groups influence peptide flexibility and seem to play a crucial role in determining the aggregates' supramolecular architectures.

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.

Design of a Hyperstable Protein by Rational Consideration of Unfolded State Interactions

Stabilization of proteins is a long-sought objective. Targeting the unfolded state interactions of a protein is not a method used for this purpose, although many proteins are known to contain such interactions. The N-terminal domain of ribosomal protein L9 (NTL9) has a lysine residue at position 12, which makes strong non-native interactions in the unfolded state. Substitution of a d-alanine for G34 in NTL9 is known to stabilize the protein by reducing the entropy of the unfolded state. Here we combine these two mutations to design a hyperstable protein. The structure of the variant is the same as that of wild-type as judged by 2D NMR. The variant is hyperstable as judged by denaturation experiments, where complete thermal unfolding of the protein does not occur in native buffer.

Imperfect DNA mirror repeats in E. coli TnsA and other protein-coding DNA

DNA imperfect mirror repeats (DNA-IMRs) are ubiquitous in protein-coding DNA. However, they overlap and often have different centers of symmetry, making it difficult to evaluate their relationship to each other and to specific DNA and protein motifs and structures. This paper describes a systematic method of determining a hierarchy for DNA-IMRs and evaluates their relationship to protein structural elements (PSEs)--helices, turns and beta-sheets. DNA-IMRs are identifed by two different methods--DNA-IMRs terminated by reverse dinucleotides (rd-IMRs) and DNA-IMRs terminated by a single (mono) matching nucleotide (m-IMRs). Both rd-IMRs and m-IMRs are evaluated in 17 proteins, and illustrated in detail for TnsA. For each of the proteins, Fisher's exact test (FET) is used to measure the coincidence between the terminal dinucleotides of rd-IMRs and the terminal amino acids of individual PSEs. A significant correlation over a span of about 3 nt was found for each protein. The correlation is robust and for most genes, all rd-IMRs<or=13 nt can be removed without the loss of statistical significance. In TnsA, the protein intervals translated by rd-IMRs>16 nt contain approximately 88% of the potential functional motifs. The protein translation of the longest rd- and m-IMRs span sequences important to the protein's structure and function. In all 17 proteins studied, the population of rd-IMRs is substantially less than the expected number and the population of m-IMRs greater than the expected number, indicating strong selective pressures. The association of rd-IMRs with PSEs restricts their spatial distribution, and therefore, their number. The greater than predicted number of m-IMRs indicates that DNA symmetry exists throughout the entire protein-coding region and may stabilize the sequence.

Structural features of thermozymes

Enzymes synthesized by thermophiles and hyperthermophiles are known as thermozymes. These enzymes are typically thermostable, or resistant to irreversible inactivation at high temperatures, and thermophilic, i.e. optimally active at elevated temperatures between 60 and 125 degrees C. Enzyme thermostability encompasses thermodynamic stability and kinetic stability. Thermodynamic stability is defined by the enzyme's free energy of stabilization (deltaG(stab)) and by its melting temperature (Tm). An enzyme's kinetic stability is often expressed as its halflife (t1/2) at defined temperature. DeltaG(stab) of thermophilic proteins is 5-20 kcal/mol higher than that of mesophilic proteins. The thermostability mechanisms for thermozymes are varied and depend on the enzyme; nevertheless, some common features can be identified as contributing to stability. These features include more interactions (i.e. hydrogen bonds, electrostatic interactions, hydrophobic interactions, disulfide bonds, metal binding) than in less stable enzymes and superior conformational structure (i.e. more rigid, higher packing efficiency, reduced entropy of unfolding, conformational strain release and stability of alpha-helix). Understanding of the stabilizing features will greatly facilitate reengineering of some of the mesozymes to more stable thermozymes.

Exploiting the Right Side of the Ramachandran Plot: Substitution of Glycines by d-Alanine Can Significantly Increase Protein Stability

A major goal of protein engineering is the enhancement of protein stability. Here we demonstrate a rational method for enhancing the stability of globular proteins by targeting glycine residues which adopt conformations with Phi > 0. Replacement of such a glycine by d-alanine can lead to a significant increase in stability. The approach is tested at three sites in two model proteins. NMR and CD indicated that the substitutions do not alter the structure. Replacement of glycine-24 of the N-terminal domain of L9 (NTL9) with d-Ala results in an increase in stability of 1.3 kcal mol-1, while replacement of glycine-34 of NTL9 leads to an increase of 1.9 kcal mol-1. Replacement of glycine-331 of the UBA domain with d-Ala leads to an increase in stability of 0.6 kcal mol-1.

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.