Structural effects of point mutations in proteins

A structural database of 11 families of chains differing by a single amino acid substitution has been built. Another structural dataset of 5 families with identical sequences has been used for comparison. The RMSD computed after a global superimposition of the mutated protein on each native one is smaller than the RMSD calculated among proteins of identical sequences. The effect of the perturbation is very local, and not necessarily the highest at the position of the mutation. A RMSD between mutated and native proteins is computed over a 3‐residue or a 7‐residue window at each position. To separate the effects of structural fluctuations due to point mutations from other sources, pair RMSD have been translated into P values which themselves are included in a score called P‐RANK. This score allows highlighting small backbone distortions by comparing these RMSD between mutated and native positions to the RMSD at the same positions in the absence of a mutation. It results from the P‐RANK that 38% of all mutations produce a significant effect on the displacement. When compared with a random distribution of RMSD at un‐mutated positions, we show that, even if the RMSD is greater when the mutation is in loops than in regular secondary structure, the relative effect is more important for regular secondary structures and for buried positions. We confirm the absence of correlation between RMSD and the predicted variation of free energy of folding but we found a small correlation between high RMSD and the error in the prediction of ΔΔG.     

The stability of Taq DNA polymerase results from a reduced entropic folding penalty; identification of other thermophilic proteins with similar folding thermodynamics

The thermal stability of Taq DNA polymerase is well known, and is the basis for its use in PCR. A comparative thermodynamic characterization of the large fragment domains of Taq (Klentaq) and E. coli (Klenow) DNA polymerases has been performed by obtaining full Gibbs‐Helmholtz stability curves of the free energy of folding (ΔG) versus temperature. This analysis provides the temperature dependencies of the folding enthalpy and entropy (ΔH and ΔS), and the heat capacity (ΔC_p) of folding. If increased or enhanced non‐covalent bonding in the native state is responsible for enhanced thermal stabilization of a protein, as is often proposed, then an enhanced favourable folding enthalpy should, in general, be observed for thermophilic proteins. However, for the Klenow–Klentaq homologous pair, the folding enthalpy (ΔH_fold) of Klentaq is considerably less favorable than that of Klenow at all temperatures. In contrast, it is found that Klentaq's extreme free energy of folding (ΔG_fold) originates from a significantly reduced entropic penalty of folding (ΔS_fold). Furthermore, the heat capacity changes upon folding are similar for Klenow and Klentaq. Along with this new data, comparable extended analysis of available thermodynamic data for 17 other mesophilic–thermophilic protein pairs (where enough applicable thermodynamic data exists) shows a similar pattern in seven of the 18 total systems. When analyzed with this approach, the more familiar “reduced ΔC_p mechanism” for protein thermal stabilization (observed in a different six of the 18 systems) frequently manifests as a temperature dependent shift from enthalpy driven stabilization to a reduced‐entropic‐penalty model. Proteins 2014; 82:785–793. © 2013 Wiley Periodicals, Inc.     

Iterative annealing mechanism explains the functions of the GroEL and RNA chaperones

Molecular chaperones are ATP‐consuming machines, which facilitate the folding of proteins and RNA molecules that are kinetically trapped in misfolded states. Unassisted folding occurs by the kinetic partitioning mechanism according to which folding to the native state, with low probability as well as misfolding to one of the many metastable states, with high probability, occur rapidly. GroEL is an all‐purpose stochastic machine that assists misfolded substrate proteins to fold. The RNA chaperones such as CYT‐19, which are ATP‐consuming enzymes, help the folding of ribozymes that get trapped in metastable states for long times. GroEL does not interact with the folded proteins but CYT‐19 disrupts both the folded and misfolded ribozymes. The structures of GroEL and RNA chaperones are strikingly different. Despite these differences, the iterative annealing mechanism (IAM) quantitatively explains all the available experimental data for assisted folding of proteins and ribozymes. Driven by ATP binding and hydrolysis and GroES binding, GroEL undergoes a catalytic cycle during which it samples three allosteric states, T (apo), R (ATP bound), and R^″ (ADP bound). Analyses of the experimental data show that the efficiency of the GroEL–GroES machinery and mutants is determined by the resetting rate k _{R ″ → T} , which is largest for the wild‐type (WT) GroEL. Generalized IAM accurately predicts the folding kinetics of Tetrahymena ribozyme and its variants. Chaperones maximize the product of the folding rate and the steady‐state native state fold by driving the substrates out of equilibrium. Neither the absolute yield nor the folding rate is optimized.     

Role of electrostatic interactions for the stability and folding behavior of cold shock protein

The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Learning about protein folding via potential functions

Over the last few years we have developed an empirical potential function that solves the protein structure recognition problem: given the sequence for an n-residue globular protein and a collection of plausible protein conformations, including the native conformation for that sequence, identify the correct, native conformation. Having determined this potential on the basis of only some 6500 native/nonnative pairs of structures for 58 proteins, we find it recognizes the native conformation for essentially all compact, soluble, globular proteins having known native conformations in comparisons with 10⁴ to 10⁶ reasonable alternative conformations apiece. In this sense, the potential encodes nearly all the essential features of globular protein conformational preference. In addition it “knows” about many additional factors in protein folding, such as the stabilization of multimeric proteins, quaternary structure, the role of disulfide bridges and ligands, proproteins vs. processed proteins, and minimal strand lengths in globular proteins. Comparisons are made with other sorts of protein folding problems, and applications in protein conformational determination and prediction are discussed. © 1994 Wiley-Liss, Inc.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. III. Dynamics of long‐range hydrophobic interactions

A 20‐residue peptide, IG(42–61), derived from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptoccocus was studied using circular dichroism, nuclear magnetic resonance (NMR) spectroscopy at various temperatures and by differential scanning calorimetry (DSC). Unlike other related peptides studied so far, this peptide displays two heat capacity peaks in DSC measurements (at a scanning rate of 1.5 deg/min at a peptide concentration of 0.07 mM), which suggests a three‐state folding/unfolding process. The results from DSC and NMR measurements suggest the formation of a dynamic network of hydrophobic interactions stabilizing the structure, which resembles a β‐hairpin shape over a wide range of temperatures (283–313 K). Our results show that IG (42–61) possesses a well‐organized three‐dimensional structure stabilized by long‐range hydrophobic interactions (Tyr50 ··· Phe57 and Trp48 ··· Val59) at T = 283 K and (Trp48 ··· Val59) at 305 and 313 K. The mechanism of β‐hairpin folding and unfolding, as well as the influence of peptide length on its conformational properties, are also discussed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Conformations of folded proteins in restricted spaces

A new method is presented to examine the complete range of folded topologies accessible in the compact state of globular proteins. The procedure is to generate all conformations, with volume exclusion, upon a lattice in a space restricted to the individual protein's known compact conformational space. Using one lattice point per residue, we find 10(2)-10(4) possible compact conformations for the five small globular proteins studied. Subsequently, these conformations are evaluated in terms of residue-specific, pairwise contact energies that favor nonbonded, hydrophobic interactions. Native structures for the five proteins are always found within the best 2% of all conformers generated. This novel method is simple and general and can be used to determine a small group of most favorable overall arrangements for the folding of specific amino acid sequences within a restricted space.     

Unique biochemical and behavioral alterations in drosophilashibirets1 mutants imply a conformational state affecting dynamin subcellular distribution and synaptic vesicle cycling

ABSTRACT Dynamin is a GTPase protein that is essential for clathrin‐mediated endocytosis of synaptic vesicle membranes. The Drosophila dynamin mutation shi^ts1 changes a single residue (G273D) at the boundary of the GTPase domain. In cell fractionation of homogenized fly heads without monovalent cations, all dynamin was in pellet fractions and was minimally susceptible to Triton‐X extraction. Addition of Na⁺ or K⁺ can extract dynamin to the cytosolic (supernatant) fraction. The shi^ts1 mutation reduced the sensitivity of dynamin to salt extraction compared with other temperature‐sensitive alleles or wild type. Sensitivity to salt extraction in shi^ts1 was enhanced by GTP and nonhydrolyzable GTP‐γS. The shi^ts1 mutation may therefore induce a conformational change, involving the GTP binding site, that affects dynamin aggregation. Temperature‐sensitive shibire mutations are known to arrest endocytosis at restrictive temperatures, with concomitant accumulation of presynaptic collared pits. Consistent with an effect upon dynamin aggregation, intact shi^ts1 flies recovered much more slowly from heat‐induced paralysis than did other temperature‐sensitive shibire mutants. Moreover, a genetic mutation that lowers GTP abundance (awd^msf15), which reduces the paralytic temperature threshold of other temperature‐sensitive shibire mutations that lie closer to consensus GTPase motifs, did not reduce the paralytic threshold of shi^ts1. Taken together, the results may link the GTPase domain to conformational shifts that influence aggregation in vitro and endocytosis in vivo, and provide an unexpected point of entry to link the biophysical properties of dynamin to physiological processes at synapses. © 2002 Wiley Periodicals, Inc. J Neurobiol 53: 319–329, 2002     

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.     

The conformation of a nascent polypeptide inside the ribosome tunnel affects protein targeting and protein folding

In this report, we describe insights into the function of the ribosome tunnel that were obtained through an analysis of an unusual 25 residue N‐terminal motif (EspP_1‐25) associated with the signal peptide of the Escherichia coli EspP protein. It was previously shown that EspP_1‐25 inhibits signal peptide recognition by the signal recognition particle, and we now show that fusion of EspP_1‐25 to a cytoplasmic protein causes it to aggregate. We obtained two lines of evidence that both of these effects are attributable to the conformation of EspP_1‐25 inside the ribosome tunnel. First, we found that mutations in EspP_1‐25 that abolished its effects on protein targeting and protein folding altered the cross‐linking of short nascent chains to ribosomal components. Second, we found that a mutation in L22 that distorts the tunnel mimicked the effects of the EspP_1‐25 mutations on protein biogenesis. Our results provide evidence that the conformation of a polypeptide inside the ribosome tunnel can influence protein folding under physiological conditions and suggest that ribosomal mutations might increase the solubility of at least some aggregation‐prone proteins produced in E. coli.     

Co‐evolutionary analysis implies auxiliary functions of HSP110 in Plasmodium falciparum

Plasmodium falciparum encounters frequent environmental challenges during its life cycle which makes productive protein folding immensely challenging for its metastable proteome. To identify the important components of protein folding machinery involved in maintaining P. falciparum proteome, we performed a proteome‐wide phylogenetic profiling across various species. We found that except HSP110, the parasite lost all other cytosolic nucleotide exchange factors essential for regulating HSP70 which is the centrum of the protein folding network. Evolutionary and structural analysis shows that besides its canonical interaction with HSP70, PfHSP110 has acquired sequence insertions for additional dynamic interactions. Molecular co‐evolution profile depicts that the co‐evolving proteins of PfHSP110 belong to distinct pathways like genetic variation, DNA repair, fatty acid biosynthesis, protein modification/trafficking, molecular motions, and apoptosis. These proteins exhibit unique physiochemical properties like large size, high iso‐electric point, low solubility, and antigenicity, hence require PfHSP110 chaperoning to attain functional state. Co‐evolving protein interaction network suggests that PfHSP110 serves as an important hub to coordinate protein quality control, survival, and immune evasion pathways in the parasite. Overall, our findings highlight potential accessory roles of PfHSP110 that may provide survival advantage to the parasite during its lifecycle and febrile conditions. The data also open avenues for experimental validation of auxiliary functions of PfHSP110 and their exploration for design of better antimalarial strategies. Proteins 2015; 83:1513–1525. © 2015 Wiley Periodicals, Inc.     

How the folding rates of two‐ and multistate proteins depend on the amino acid properties

Proteins fold by either two‐state or multistate kinetic mechanism. We observe that amino acids play different roles in different mechanism. Many residues that are easy to form regular secondary structures (α helices, β sheets and turns) can promote the two‐state folding reactions of small proteins. Most of hydrophilic residues can speed up the multistate folding reactions of large proteins. Folding rates of large proteins are equally responsive to the flexibility of partial amino acids. Other properties of amino acids (including volume, polarity, accessible surface, exposure degree, isoelectric point, and phase transfer energy) have contributed little to folding kinetics of the proteins. Cysteine is a special residue, it triggers two‐state folding reaction and but inhibits multistate folding reaction. These findings not only provide a new insight into protein structure prediction, but also could be used to direct the point mutations that can change folding rate. Proteins 2014; 82:2375–2382. © 2014 Wiley Periodicals, Inc.     

Export of a hyperexpressed mammalian globular cytochrome b5 precursor in Escherichia coli is dramatically affected by the nature of the amino acid flanking the secretory signal sequence cleavage bond

A chimeric mammalian globular cytochrome b₅ fused to Escherichia coli alkaline phosphatase signal sequence (SS) was used as a model probe to investigate the influence of substituting each one of the standard 20 amino acids at its N‐terminus on the Sec‐dependent export of the precursor to the periplasmic space of E. coli. Substituting the native Met⁺¹ of the passenger protein flanking the SS with any one of the remaining 19 amino acids introduced significant changes in the export of cytochrome b₅ without jamming the Sec‐dependent translocon. Acidic and hydrophilic residues proved to be the most efficient promoters of export. Small, nonbulky and basic residues yielded intermediate levels of the hemoprotein export. Replacement with a Cys⁺¹ residue generated significant quantities of both monomeric and disulfide‐linked dimeric forms. However, bulky, aromatic and hydrophobic residues caused a significant decline in the rates of secretion. In expectation with their absences in the natural periplasmically secreted proteins, Pro and Ile‐tagged cytochrome b₅ precursors failed to generate any detectable secreted recombinant products. Although Ala, amongst the native E. coli periplasmic proteins, is the preferred X⁺¹ residue with an occurrence of 50% frequency, it proved half as effective in promoting export when inserted proximally to the SS of cytochrome b₅. The mechanisms involved for these export variations are discussed. The findings will prove beneficial for high‐level generation of recombinant proteins by secretory means for pharmaceutical and related biotechnological applications.     

Nucleation‐based prediction of the protein folding rate and its correlation with the folding nucleus size

The problem of protein self‐organization is in the focus of current molecular biology studies. Although the general principles are understood, many details remain unclear. Specifically, protein folding rates are of interest because they dictate the rate of protein aggregation which underlies many human diseases. Here we offer predictions of protein folding rates and their correlation with folding nucleus sizes. We calculated free energies of the transition state and sizes of folding nuclei for 84 proteins and peptides whose other parameters were measured at the point of thermodynamic equilibrium between their unfolded and native states. We used the dynamic programming method where each residue was considered to be either as folded as in its native state or completely disordered. The calculated and measured folding rates showed a good correlation at the temperature mid‐transition point (the correlation coefficient was 0.75). Also, we pioneered in demonstrating a moderate (‐0.57) correlation coefficient between the calculated sizes of folding nuclei and the folding rates. Predictions made by different methods were compared. The established good correlation between the estimated free energy barrier and the experimentally found folding rate of each studied protein/peptide indicates that our model gives reliable results for the considered data set. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Noncollagenous region of the streptococcal collagen‐like protein is a trimerization domain that supports refolding of adjacent homologous and heterologous collagenous domains

Proper folding of the (Gly‐Xaa‐Yaa)_n sequence of animal collagens requires adjacent N‐ or C‐terminal noncollagenous trimerization domains which often contain coiled‐coil or beta sheet structure. Collagen‐like proteins have been found recently in a number of bacteria, but little is known about their folding mechanism. The Scl2 collagen‐like protein from Streptococcus pyogenes has an N‐terminal globular domain, designated V_sp, adjacent to its triple‐helix domain. The V_sp domain is required for proper refolding of the Scl2 protein in vitro. Here, recombinant V_sp domain alone is shown to form trimers with a significant α‐helix content and to have a thermal stability of T_m = 45°C. Examination of a new construct shows that the V_sp domain facilitates efficient in vitro refolding only when it is located N‐terminal to the triple‐helix domain but not when C‐terminal to the triple‐helix domain. Fusion of the V_sp domain N‐terminal to a heterologous (Gly‐Xaa‐Yaa)_n sequence from Clostridium perfringens led to correct folding and refolding of this triple‐helix, which was unable to fold into a triple‐helical, soluble protein on its own. These results suggest that placement of a functional trimerization module adjacent to a heterologous Gly‐Xaa‐Yaa repeating sequence can lead to proper folding in some cases but also shows specificity in the relative location of the trimerization and triple‐helix domains. This information about their modular nature can be used in the production of novel types of bacterial collagen for biomaterial applications.     

Information and redundancy in the burial folding code of globular proteins within a wide range of shapes and sizes

Recent ab initio folding simulations for a limited number of small proteins have corroborated a previous suggestion that atomic burial information obtainable from sequence could be sufficient for tertiary structure determination when combined to sequence‐independent geometrical constraints. Here, we use simulations parameterized by native burials to investigate the required amount of information in a diverse set of globular proteins comprising different structural classes and a wide size range. Burial information is provided by a potential term pushing each atom towards one among a small number L of equiprobable concentric layers. An upper bound for the required information is provided by the minimal number of layers L^min still compatible with correct folding behavior. We obtain L^min between 3 and 5 for seven small to medium proteins with 50 ≤ N_r ≤ 110 residues while for a larger protein with N_r = 141 we find that L ≥ 6 is required to maintain native stability. We additionally estimate the usable redundancy for a given L ≥ L^min from the burial entropy associated to the largest folding‐compatible fraction of “superfluous” atoms, for which the burial term can be turned off or target layers can be chosen randomly. The estimated redundancy for small proteins with L = 4 is close to 0.8. Our results are consistent with the above‐average quality of burial predictions used in previous simulations and indicate that the fraction of approachable proteins could increase significantly with even a mild, plausible, improvement on sequence‐dependent burial prediction or on sequence‐independent constraints that augment the detectable redundancy during simulations. Proteins 2016; 84:515–531. © 2016 Wiley Periodicals, Inc.     

Pairwise contact energy statistical potentials can help to find probability of point mutations

To adopt a particular fold, a protein requires several interactions between its amino acid residues. The energetic contribution of these residue–residue interactions can be approximated by extracting statistical potentials from known high resolution structures. Several methods based on statistical potentials extracted from unrelated proteins are found to make a better prediction of probability of point mutations. We postulate that the statistical potentials extracted from known structures of similar folds with varying sequence identity can be a powerful tool to examine probability of point mutation. By keeping this in mind, we have derived pairwise residue and atomic contact energy potentials for the different functional families that adopt the (α/β)₈ TIM‐Barrel fold. We carried out computational point mutations at various conserved residue positions in yeast Triose phosphate isomerase enzyme for which experimental results are already reported. We have also performed molecular dynamics simulations on a subset of point mutants to make a comparative study. The difference in pairwise residue and atomic contact energy of wildtype and various point mutations reveals probability of mutations at a particular position. Interestingly, we found that our computational prediction agrees with the experimental studies of Silverman et al. (Proc Natl Acad Sci 2001;98:3092–3097) and perform better prediction than i_Mutant and Cologne University Protein Stability Analysis Tool. The present work thus suggests deriving pairwise contact energy potentials and molecular dynamics simulations of functionally important folds could help us to predict probability of point mutations which may ultimately reduce the time and cost of mutation experiments. Proteins 2016; 85:54–64. © 2016 Wiley Periodicals, Inc.     

Specificity and cooperativity at β-lactamase position 104 in TEM-1/BLIP and SHV-1/BLIP interactions

Establishing a quantitative understanding of the determinants of affinity in protein–protein interactions remains challenging. For example, TEM‐1/β‐lactamase inhibitor protein (BLIP) and SHV‐1/BLIP are homologous β‐lactamase/β‐lactamase inhibitor protein complexes with disparate K_d values (3 nM and 2 μM, respectively), and a single substitution, D104E in SHV‐1, results in a 1000‐fold enhancement in binding affinity. In TEM‐1, E104 participates in a salt bridge with BLIP K74, whereas the corresponding SHV‐1 D104 does not in the wild type SHV‐1/BLIP co‐structure. Here, we present a 1.6 Å crystal structure of the SHV‐1 D104E/BLIP complex that demonstrates that this point mutation restores this salt bridge. Additionally, mutation of a neighboring residue, BLIP E73M, results in salt bridge formation between SHV‐1 D104 and BLIP K74 and a 400‐fold increase in binding affinity. To understand how this salt bridge contributes to complex affinity, the cooperativity between the E/K or D/K salt bridge pair and a neighboring hot spot residue (BLIP F142) was investigated using double mutant cycle analyses in the background of the E73M mutation. We find that BLIP F142 cooperatively stabilizes both interactions, illustrating how a single mutation at a hot spot position can drive large perturbations in interface stability and specificity through a cooperative interaction network. Proteins 2011. © 2011 Wiley‐Liss, Inc.