Valine 44 and valine 45 of human glutathione synthetase are key for subunit stability and negative cooperativity

It was hypothesized that residues Val44 and Val45 serve as important residues for human glutathione synthetase (hGS) function and stability given their location at the dimer interface of this enzyme. Computational studies suggest that mutation at Val45 has more impact on the structure and stability of hGS than does mutation at Val44. Experimentally, enzymes with mutations at the 44 and or 45 positions of hGS were prepared, purified and assayed for initial activity. Val45 position mutations (either to alanine or tryptophan) have a greater impact on enzyme activity than do mutations at Val44. Differential scanning calorimetry experiments reveal a loss of stability in all mutant enzymes, with V45 mutations being less stable than the corresponding Val44 mutations. The γ-GluABA substrate affinity remains unaltered in V44A and V45A mutant enzymes, but increases when tryptophan is introduced at either of these positions. Hill coefficients trend towards less negative cooperativity with the exception of V45W mutant hGS. These results imply that residues V44 and V45 are located along the allosteric pathway of this negatively cooperative dimeric enzyme, that their mutation impacts the allosteric pathway more than it does the active site of hGS, and that these residues (and by extension the dimer interface in which they are located) are integral to the stability of human glutathione synthetase.     

Conservative mutation Met8 --> Leu affects the folding process and structural stability of squash trypsin inhibitor CMTI-I

Protein molecules can accommodate a large number of mutations without noticeable effects on their stability and folding kinetics. On the other hand, some mutations can have quite strong effects on protein conformational properties. Such mutations either destabilize secondary structures, e.g., alpha-helices, are incompatible with close packing of protein hydrophobic cores, or lead to disruption of some specific interactions such as disulfide cross links, salt bridges, hydrogen bonds, or aromatic-aromatic contacts. The Met8 --> Leu mutation in CMTI-I results in significant destabilization of the protein structure. This effect could hardly be expected since the mutation is highly conservative, and the side chain of residue 8 is situated on the protein surface. We show that the protein destabilization is caused by rearrangement of a hydrophobic cluster formed by side chains of residues 8, Ile6, and Leu17 that leads to partial breaking of a hydrogen bond formed by the amide group of Leu17 with water and to a reduction of a hydrophobic surface buried within the cluster. The mutation perturbs also the protein folding. In aerobic conditions the reduced wild-type protein folds effectively into its native structure, whereas more then 75% of the mutant molecules are trapped in various misfolded species. The main conclusion of this work is that conservative mutations of hydrophobic residues can destabilize a protein structure even if these residues are situated on the protein surface and partially accessible to water. Structural rearrangement of small hydrophobic clusters formed by such residues can lead to local changes in protein hydration, and consequently, can affect considerably protein stability and folding process.     

Impact of intra-subunit domain-domain interactions on creatine kinase activity and stability

Creatine kinase (CK) is a key enzyme in vertebrate excitable tissues. In this research, five conserved residues located on the intra-subunit domain-domain interface were mutated to explore their role in the activity and structural stability of CK. The mutations of Val72 and Gly73 decreased both the activity and stability of CK. The mutations of Cys74 and Val75, which had no significant effect on CK activity and structure, gradually decreased the stability and reactivation of CK. Our results suggested that the mutations might modify the correct positioning of the loop contributing to domain-domain interactions, and result in decreased stability against denaturation.     

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.     

A New Method to Characterize Hydrophobic Organization of Proteins: Application to Rational Protein Engineering of Barnase

Abstract

We present a new algorithm for characterization of protein spatial structure basing on the molecular hydrophobicity potential approach. The method is illustrated by the analysis of three-dimensional structure of barnase and barnase-barstar complex. Current approach enables identification of amino acid residues situated in unfavorable environment (these residues may be “active” for binding), and to map quantitatively hydrophobic, hydrophilic and unfavorable hydrophobic-hydrophilic intra-and inter-molecular contacts involving backbone and side-chain segments of amino acid residues. Calculation of individual contributions of amino acid residues to such contacts permits identification of structurally-important residues. The contact plots obtained with molecular hydrophobicity potential calculations, provide easy rules to choose sites for mutations, which can increase a strength of intra- or inter-molecular hydrophobic interactions. The unfavorable hydrophobic-hydrophilic contact can be mutated to favorable hydrophobic, and already existing weak hydrophobic contact can be strengthen by increasing hydrophobicity of residues in contact. Basing on the analysis of the contact plots, we suggest several mutations of barnase which are supposed to increase intramolecular hydrophobic interactions, and thus might lead to increased stability of the protein. Part of these mutations was studied previously experimentally, and indeed stabilized barnase. The other of predicted mutations were not studied experimentally yet. Several new mutations of barnase and barstar are also proposed to enhance the hydrophobic interactions on their binding interface.     

Simulations of SIN mutations and histone variants in human nucleosomes reveal altered protein-DNA and core histone interactions

Alterations in the stability of a nucleosome exert predominant influence on chromatin structure and eukaryotic gene expression. In an attempt to investigate the mononucleosome stability using computational approaches, we have simulated the structure of a human mononucleosome and have compared their energies under the influence of core mutations, tail substitutions, variant histones, and orthologs. We observe that mutant nucleosomes carrying SIN (SWI Independent) mutations do not alter the overall nucleosomal structure but cause local structural changes leading to significant changes in energy and hence the stability. We observe that the nucleosome stability is altered by the substitution of only certain critical lysine residues on the H3 tails. Interestingly, the incorporation of variants H2A.Z and H3.3 lower nucleosome stability as evidenced by small energy changes. However, the substitution of histone orthologs did not alter structural stability. Our simulations to determine the nucleosome stability using energy trends emphasize the role of mutations, variants, and orthologs as determinants of chromatin structure at the nucleosome core particle level. The destabilization we observe on the human nucleosome with core mutations show similar trends of instability as validated experimentally in yeast.     

Stabilization of proteins by enhancement of inter-residue hydrophobic contacts: lessons of T4 lysozyme and barnase

Although the hydrophobic interactions are considered as the main contributors to the protein stability, not much examples of protein stabilization by rational increasing of this type of interactions still can be found in literature. This is partly due to the lack of proper theoretical "measure" of hydrophobic interactions and their changes upon mutations. In the present paper the molecular hydrophobicity potential approach is used to assess how the changes in type and the strength of inter-residue contacts upon single amino acid mutations are correlated with the changes in thermodynamic stability of T4 lysozyme and barnase mutants, and which factors affect these correlations. Mutations changing unfavorable hydrophilic-to-hydrophobic contacts into favorable hydrophobic were found to enhance the thermodynamic stability in more than 81 % of cases, if these mutations do not create steric bumps and do not involve proline residues and hydrogen-bonded side-chains. Mutations increasing hydrophobic contributions (according to molecular hydrophobicity potential formalism) lead to increase of thermodynamic stability in more than 94% of cases for certain type of mutations (i.e., mutations not involving charged residues, Pro and residues with side-chain hydrogen bonds, when these mutations do not introduce steric bumps and do not involve strongly exposed residues and residues situated at helix N- and C-cap positions). For this type of mutations the correlation was found between the change in hydrophobic contributions of mutated residues deltaCphob and thermodynamic parameters deltaTm (change in melting temperature) and deltadeltaG (change in free energy of unfolding). Although the correlation coefficients were larger if the experimental structures of mutants were used for the calculations (correlation coefficients r(exp) deltaC,deltaT = .85 and r(exp) deltaC,deltadeltaG = 0.87) than if the modeled structures were used instead (r(mod) deltaC,deltaT = 0.74 and r(mod)deltaC,deltadeltaG = 0.76), the modelled structures of mutants in the vast majority of cases can be used for qualitative predictition of the protein stabilization. Basing on the analysis of mutations increasing hydrophobic contributions in T4 lysozyme the substitution matrix was derived, which can be used to decide which new residue should be put instead the old one to increase the stability of protein. The estimation shows that the number of potential mutation sites for enhancement of hydrophobic interactions in T4 lysozyme is quite large, and only approximately 10 per cent of them were studied thus far. Basing on the current analysis of T4 lysozyme and barnase mutations the algorithm for increasing of protein stability via increasing of hydrophobic interactions for the proteins with known spatial structure is proposed.     

Simulations of SIN Mutations and Histone Variants in Human Nucleosomes Reveal Altered Protein-DNA and Core Histone Interactions

Abstract Alterations in the stability of a nucleosome exert predominant influence on chromatin structure and eukaryotic gene expression. In an attempt to investigate the mononucleosome stability using computational approaches, we have simulated the structure of a human mononucleosome and have compared their energies under the influence of core mutations, tail substitutions, variant histones, and orthologs. We observe that mutant nucleosomes carrying SIN (SWI Independent) mutations do not alter the overall nucleosomal structure but cause local structural changes leading to significant changes in energy and hence the stability. We observe that the nucleosome stability is altered by the substitution of only certain critical lysine residues on the H3 tails. Interestingly, the incorporation of variants H2A.Z and H3.3 lower nucleosome stability as evidenced by small energy changes. However, the substitution of histone orthologs did not alter structural stability. Our simulations to determine the nucleosome stability using energy trends emphasize the role of mutations, variants, and orthologs as determinants of chromatin structure at the nucleosome core particle level. The de-stabilization we observe on the human nucleosome with core mutations show similar trends of instability as validated experimentally in yeast.     

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

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

Impact of intra-subunit interactions on the dimeric arginine kinase activity and structural stability

Arginine kinase (AK) catalyzes the reversible phosphorylation of arginine by ATP, yielding the phosphoarginine. In this research, six conserved residues located on the intra-subunit domain-domain interfaces were mutated to explore their roles in the activity and structural stability of dimer AK. The mutations D69A, E70A, E71A and F80A led to pronounced loss of AK activity and structural stability. Although the mutations V75A and F76A had little effect on AK activity and structure, they caused gradually decreased the stability and reactivation of dimer AK. Our results suggested that the mutations might affect the correct positioning of the N-loop and C-loop thus disrupted the efficient recognition and interactions between the N-terminal domain and C-terminal domain which may influence the compact dimer structure, and result in decreased activity and structural stability.     

Toward Fast Determination of Protein Stability Maps: Experimental and Theoretical Analysis of Mutants of a Nocardiopsis prasina Serine Protease

The stability of serine proteases is of major importance for their application in industrial processes. Here we study the determinants of the stability of a Nocardiopsis prasina serine protease using fast residual activity assays, a feature classification algorithm, and structure-based energy calculation algorithms for 121 micropurified mutant enzyme clones containing multiple point mutations. Using a multivariate regression analysis, we deconvolute the data for the mutant clones and find that mutations of residues Asn47 and Pro124 are deleterious to the stability of the enzyme. Both of these residues are situated in loops that are known to be important for the stability of the highly homologous α-lytic protease. Structure-based energy calculations with PEATSA give a good general agreement with the trend of experimentally measured values but also identify a number of clones that the algorithm fails to predict correctly. We discuss the significance of the results in relation to the structure and function of closely related proteases, comment on the optimal experimental design when performing high-throughput experiments for characterizing the determinants of protein stability, and discuss the performance of structure-based energy calculations with complex data sets such as the one presented here.     

Cysteine‐free rop: A four‐helix bundle core mutant has wild‐type stability and structure but dramatically different unfolding kinetics

Cysteine residues can complicate the folding and storage of proteins due to improper formation of disulfide bonds or oxidation of residues that are natively reduced. Wild‐type Rop is a homodimeric four‐helix bundle protein and an important model for protein design in the understanding of protein stability, structure and folding kinetics. In the native state, Rop has two buried, reduced cysteine residues in its core, but these are prone to oxidation in destabilized variants, particularly upon extended storage. To circumvent this problem, we designed and characterized a Cys‐free variant of Rop, including solving the 2.3 Å X‐ray crystal structure. We show that the C38A C52V variant has similar structure, stability and in vivo activity to wild‐type Rop, but that it has dramatically faster unfolding kinetics like virtually every other mutant of Rop that has been characterized. This cysteine‐free Rop has already proven useful for studies on solution topology and on the relationship of core mutations to stability. It also suggests a general strategy for removal of pairs of Cys residues in proteins, both to make them more experimentally tractable and to improve their storage properties for therapeutic or industrial purposes.     

Folding of bovine pancreatic trypsin inhibitor (BPTI) variants in which almost half the residues are alanine

Recent studies indicate that a fraction of the information contained in an amino acid sequence may be sufficient for specifying a native protein structure. An earlier alanine-scanning experiment conducted on bovine pancreatic trypsin inhibitor (BPTI; 58 residues) suggested that if cumulative mutations have additive effects on protein stability, a native protein structure could be built from BPTI sequences that contained many alanine residues distributed throughout the protein. To test this hypothesis, we designed and produced six BPTI mutants containing from 21 to 29 alanine residues. We found that the melting temperature of mutants containing up to 27 alanine residues (48 % of the total number of residues) could be predicted quite well by the sum of the change in melting temperature for the single mutations. Additionally, these same mutants folded into a native-like structure, as judged by their cooperative thermal denaturation curves and heteronuclear multiple quantum correlation (HMQC) NMR spectra. A BPTI mutant containing 22 alanine residues was further shown by 2D and 3D-NMR to fold into a structure very similar to that of native BPTI, and to be a functional trypsin inhibitor. These results provide insight into the extent to which native protein structure and function can be achieved with a highly simplified amino acid sequence.     

Key sites residues involved in interacting with chemicals of pheromone‐binding proteins from Lymantria dispar

Pheromone‐binding proteins (PBPs) are distributed widely on the antennae of insects, and they are believed to be involved in the process of chemical signal transduction, but their interaction with chemicals is largely unknown. Here, we present our findings on the key amino acid residues of PBPs in the gypsy moth, Lymantria dispar. Potential key residues were screened with the Calculate Mutation Energy program and molecular docking methods. Mutated proteins were obtained by mutating residues to alanine via site‐directed mutagenesis. Circular dichroism (CD) spectroscopy showed that the mutated proteins formed α‐helix, and the stability of protein structure was influenced due to mutations. Fluorescence binding assays were further conducted with the mutated proteins, sex pheromones and analogues. Results showed that to PBP 1, tyrosine at position 41 and phenylalanine at position 76 could be the key amino acid residues influencing the stability of structure; in addition, phenylalanine at 36 and lysine at position 94 could be key amino acid residues interacting with chemicals. To PBP 2, glycine at position 49, phenylalanine at position 76 and lysine at position 121 could be the key amino acid residues in the structural stability. These results shed light on the relationship between the specific amino acids and functions of PBPs in transmitting the chemical signals.     

Importance of surrounding residues for protein stability of partially buried mutations

For understanding the factors influencing protein stability, we have analyzed the relationship between changes in protein stability caused by partially buried mutations and changes in 48 physico-chemical, energetic and conformational properties of amino acid residues. Multiple regression equations were derived to predict the stability of protein mutants and the efficiency of the method has been verified with both back-check and jack-knife tests. We observed a good agreement between experimental and computed stabilities. Further, we have analyzed the effect of sequence window length from 1 to 12 residues on each side of the mutated residue to include the sequence information for predicting protein stability and we found that the preferred window length for obtaining the highest correlation is different for each secondary structure; the preferred window length for helical, strand and coil mutations are, respectively, 0, 9 and 4 residues on both sides of the mutant residues. However, all the secondary structures have significant correlation for a window length of one residue on each side of the mutant position, implying the role of short-range interactions. Extraction of surrounding residue information for various distances (3 to 20A) around the mutant position showed the highest correlation at 8A, 6A and 7A, respectively, for mutations in helical, strand and coil segments. Overall, the information about the surrounding residues within the sphere of 7 to 8A, may explain better the stability in all subsets of partially buried mutations implying that this distance is sufficient to accommodate the residues influenced by major intramolecular interactions for the stability of protein structures.     

Mutations can cause light chains to be too stable or too unstable to form amyloid fibrils

Light chain (AL) amyloidosis is an incurable human disease, where the amyloid precursor is a misfolding‐prone immunoglobulin light‐chain. Here, we identify the role of somatic mutations in the structure, stability and in vitro fibril formation for an amyloidogenic AL‐12 protein by restoring four nonconservative mutations to their germline (wild‐type) sequence. The single restorative mutations do not affect significantly the native structure, the unfolding pathway, and the reversibility of the protein. However, certain mutations either decrease (H32Y and H70D) or increase (R65S and Q96Y) the protein thermal stability. Interestingly, the most and the least stable mutants, Q96Y and H32Y, do not form amyloid fibrils under physiological conditions. Thus, Q96 and H32 are key residues for AL‐12 stability and fibril formation and restoring them to the wild‐type residues preclude amyloid formation. The mutants whose equilibrium is shifted to either the native or unfolded states barely sample transient partially folded states, and therefore do not form fibrils. These results agree with previous observations by our laboratory and others that amyloid formation occurs because of the sampling of partially folded states found within the unfolding transition (Blancas‐Mejia and Ramirez‐Alvarado, Ann Rev Biochem 2013;82:745–774). Here we provide a new insight on the AL amyloidosis mechanism by demonstrating that AL‐12 does not follow the established thermodynamic hypothesis of amyloid formation. In this hypothesis, thermodynamically unstable proteins are more prone to amyloid formation. Here we show that within a thermal stability range, the most stable protein in this study is the most amyloidogenic protein.     

Cataract-causing mutation R233H affects the stabilities of βB1- and βA3/βB1-crystallins with different pH-dependence

Disease-causing mutations can be stabilizing or destabilizing. Missense mutations of structural residues are generally destabilizing, while stabilizing mutations are usually linked to alterations in protein functions. Stabilizing mutations are rarely identified in mutations linked to congenital cataract, a disease caused by the opacification of the lens. In this research, we found that R233H mutation had little impact on βB1-crystallin structure, solubility and thermal stability under neutral solution pH conditions. The mutation increased βB1 stability against guanidine hydrochloride-induced denaturation, suggesting that Arg233 might be a functional residue. Further analysis indicated that the R233H mutation did not affect the formation of βA3/βB1 heteromer, but significantly reduced heteromer stability against heat- and guanidine hydrochloride-induced denaturation. The R233H mutation negatively affected the thermal stabilities and aggregatory propensities of βB1 and βA3/βB1 with different pH-dependence, implying that the protonation of His side chains during acidification played a regulatory role in crystallin stability and aggregation. Molecular dynamic simulations indicated that Arg233 is one of the residues forming an inter-subunit ion-pairing network with intrinsically dynamic nature. Based on these observations, we proposed that the highly dynamic ion-pairing network contributed to the tradeoff among βB1 solubility, stability, aggregatory propensity and function of protecting βA3.     

Cover Image,Volume 116, Number 2, February 2019

A key point of protein stability engineering is to identify specific target residues whose mutations can stabilize the protein structure without negatively affecting the function or activity of the protein. Here, we propose a method called RiSLnet (Rapid i dentification of Smart mutant Library using residue network) to identify such residues by combining network analysis for protein residue interactions, identification of conserved residues, and evaluation of relative solvent accessibility. To validate its performance, the method was applied to four proteins, that is, T4 lysozyme, ribonuclease H, barnase, and cold shock protein B. Our method predicted beneficial mutations in thermal stability with ~62% average accuracy when the thermal stability of the mutants was compared with the ones in the Protherm database. It was further applied to lysine decarboxylase (CadA) to experimentally confirm its accuracy and effectiveness. RiSLnet identified mutations increasing the thermal stability of CadA with the accuracy of ~60% and significantly reduced the number of candidate residues (~99%) for mutation. Finally, combinatorial mutations designed by RiSLnet and in silico saturation mutagenesis yielded a thermally stable triple mutant with the half-life (T _1/2) of 114.9 min at 58°C, which is approximately twofold higher than that of the wild-type.