PoPMuSiC,rationally designing point mutations in protein structures

PoPMuSiC is an efficient tool for rational computer-aided design of single-site mutations in proteins and peptides. Two types of queries can be submitted. The first option allows to estimate the changes in folding free energy for specific point mutations given by the user. In the second option, all possible point mutations in a given protein or protein region are performed and the most stabilizing or destabilizing mutations, or the neutral mutations with respect to thermodynamic stability, are selected. For each sequence position or secondary structure the deviation from the most stable sequence is moreover evaluated, which helps to identify the most suitable sites for the introduction of mutations.     

Hydrogen-exchange stabilities of RNase T1 and variants with buried and solvent-exposed Ala --> Gly mutations in the helix

Hydrogen-exchange rates were measured for RNase T1 and three variants with Ala --> Gly substitutions at a solvent-exposed (residue 21) and a buried (residue 23) position in the helix: A21G, G23A, and A21G + G23A. These results were used to measure the stabilities of the proteins. The hydrogen-exchange stabilities (DeltaG(HX)) for the most stable residues in each variant agree with the equilibrium conformational stability measured by urea denaturation (DeltaG(U)), if the effects of D(2)O and proline isomerization are included [Huyghues-Despointes, B. M. P., Scholtz, J. M., and Pace, C. N. (1999) Nat. Struct. Biol. 6, 210-212]. These residues also show similar changes in DeltaG(HX) upon Ala --> Gly mutations (DeltaDeltaG(HX)) as compared to equilibrium measurements (DeltaDeltaG(U)), indicating that the most stable residues are exchanging from the globally unfolded ensemble. Alanine is stabilizing compared to glycine by 1 kcal/mol at a solvent-exposed site 21 as seen by other methods for the RNase T1 protein and peptide helix [Myers, J. K., Pace, C. N., and Scholtz, J. M. (1997) Proc. Natl. Acad. Sci. U.S.A. 94, 3833-2837], while it is destabilizing at the buried site 23 by the same amount. For the A21G variant, only local NMR chemical shift perturbations are observed compared to RNase T1. For the G23A variant, large chemical shift changes are seen throughout the sequence, although X-ray crystal structures of the variant and RNase T1 are nearly superimposable. Ala --> Gly mutations in the helix of RNase T1 at both helical positions alter the native-state hydrogen-exchange stabilities of residues throughout the sequence.     

Genetic analysis of the polyomavirus DnaJ domain

Polyomavirus T antigens share a common N-terminal sequence that comprises a DnaJ domain. DnaJ domains activate DnaK molecular chaperones. The functions of J domains have primarily been tested by mutation of their conserved HPD residues. Here, we report detailed mutagenesis of the polyomavirus J domain in both large T (63 mutants) and middle T (51 mutants) backgrounds. As expected, some J mutants were defective in binding DnaK (Hsc70); other mutants retained the ability to bind Hsc70 but were defective in stimulating its ATPase activity. Moreover, the J domain behaves differently in large T and middle T. A given mutation was twice as likely to render large T unstable as it was to affect middle T stability. This apparently arose from middle T's ability to bind stabilizing proteins such as protein phosphatase 2A (PP2A), since introduction of a second mutation preventing PP2A binding rendered some middle T J-domain mutants unstable. In large T, the HPD residues are critical for Rb-dependent effects on the host cell. Residues Q32, A33, Y34, H49, M52, and N56 within helix 2 and helix 3 of the large T J domain were also found to be required for Rb-dependent transactivation. Cyclin A promoter assays showed that J domain function also contributes to large T transactivation that is independent of Rb. Single point mutations in middle T were generally without effect. However, residue Q37 is critical for middle T's ability to form active signaling complexes. The Q37A middle T mutant was defective in association with pp60(c-src) and in transformation.     

Helix mutations stabilize a late productive intermediate on the folding pathway of ubiquitin

We have investigated the relative placement of rate-limiting energy barriers and the role of productive or obstructive intermediates on the folding pathway of yeast wild-type ubiquitin ( wt-Ub) containing the F45W mutation. To manipulate the folding barriers, we have designed a family of mutants in which stabilizing substitutions have been introduced incrementally on the solvent-exposed surface of the main alpha-helix (residues 23-34), which has a low intrinsic helical propensity in the native sequence. Although the U --> I and I --> N transitions are not clearly delineated in the kinetics of wt-Ub, we show that an intermediate becomes highly populated and more clearly resolved as the predicted stability of the helix increases. The observed acceleration in the rate of folding correlates with helix stability and is consistent with the I-state representing a productive rather than misfolded state. A Leffler analysis of the effects on kinetics of changes in stability within the family of helix mutants results in a biphasic correlation in both the refolding and unfolding rates that suggest a shift from a nucleation-condensation mechanism (weakly stabilized helix) toward a diffusion-collision model (highly stabilized helix). Through the introduction of helix-stabilizing mutations, we are able to engineer a well-resolved I-state on the folding pathway of ubiquitin which is likely to be structurally distinct from that which is only weakly populated on the folding pathway of wild-type ubiquitin.     

Optimization of designed armadillo repeat proteins by molecular dynamics simulations and NMR spectroscopy

A multidisciplinary approach based on molecular dynamics (MD) simulations using homology models, NMR spectroscopy, and a variety of biophysical techniques was used to efficiently improve the thermodynamic stability of armadillo repeat proteins (ArmRPs). ArmRPs can form the basis of modular peptide recognition and the ArmRP version on which synthetic libraries are based must be as stable as possible. The 42-residue internal Arm repeats had been designed previously using a sequence-consensus method. Heteronuclear NMR revealed unfavorable interactions present at neutral but absent at high pH. Two lysines per repeat were involved in repulsive interactions, and stability was increased by mutating both to glutamine. Five point mutations in the capping repeats were suggested by the analysis of positional fluctuations and configurational entropy along multiple MD simulations. The most stabilizing single C-cap mutation Q240L was inferred from explicit solvent MD simulations, in which water penetrated the ArmRP. All mutants were characterized by temperature- and denaturant-unfolding studies and the improved mutants were established as monomeric species with cooperative folding and increased stability against heat and denaturant. Importantly, the mutations tested resulted in a cumulative decrease of flexibility of the folded state in silico and a cumulative increase of thermodynamic stability in vitro. The final construct has a melting temperature of about 85°C, 14.5° higher than the starting sequence. This work indicates that in silico studies in combination with heteronuclear NMR and other biophysical tools may provide a basis for successfully selecting mutations that rapidly improve biophysical properties of the target proteins.     

De novo design of transmembrane helix–helix interactions and measurement of stability in a biological membrane

Membrane proteins regulate a large number of cellular functions, and have great potential as tools for manipulation of biological systems. Developing these tools requires a robust and quantitative understanding of membrane protein folding and interactions within the bilayer. With this in mind, we have designed a series of proteins to probe the net thermodynamic contribution of well-known sequence motifs to transmembrane helix-helix association in a biological membrane. The proteins were designed from first principles (de novo) using current knowledge about membrane insertion and stabilizing interaction motifs. A simple poly-Leu “scaffold” was decorated with individual helix interaction motifs (G-XXX-G, polar residues, heptad repeat) to create transmembrane helix–helix interactions of increasing strength. The GALLEX assay, an in vivo assay for measurement of transmembrane helix self-association, was combined with computational methods to characterize the relative strength and mode of interaction for each sequence. In addition, the apparent free energy contribution (ΔΔG^app) of each motif to transmembrane helix self-association was measured in a biological membrane, results that are the first of their kind for these de novo designed sequences, and suggest that the free energy barrier to overcoming weak association is quite small (< 1.4 kcal mol^− 1) in a natural membrane. By quantifying and rationalizing the contribution of key motifs to transmembrane helix association, our work offers a route to direct the design of novel sequences for use in biotechnology or synthetic biology (e.g. molecular switches) and to predict the effects of sequence modification in known transmembrane domains (for control of cellular processes).     

NMR structure of the N-terminal J domain of murine polyomavirus T antigens. Implications for DnaJ-like domains and for mutations of T antigens

The NMR structure of the N-terminal, DnaJ-like domain of murine polyomavirus tumor antigens (PyJ) has been determined to high precision, with root mean square deviations to the mean structure of 0.38 A for backbone atoms and 0.94 A for all heavy atoms of ordered residues 5-41 and 50-69. PyJ possesses a three-helix fold, in which anti-parallel helices II and III are bridged by helix I, similar to the four-helix fold of the J domains of DnaJ and human DnaJ-1. PyJ differs significantly in the lengths of N terminus, helix I, and helix III. The universally conserved HPD motif appears to form a His-Pro C-cap of helix II. Helix I features a stabilizing Schellman C-cap that is probably conserved universally among J domains. On the helix II surface where positive charges of other J domains have been implicated in binding of hsp70s, PyJ contains glutamine residues. Nonetheless, chimeras that replace the J domain of DnaJ with PyJ function like wild-type DnaJ in promoting growth of Escherichia coli. This activity can be modulated by mutations of at least one of these glutamines. T antigen mutations reported to impair cellular transformation by the virus, presumably via interactions with PP2A, cluster in the hydrophobic folding core and at the extreme N terminus, remote from the HPD loop.     

Side-chain contributions to membrane protein structure and stability

The molecular forces that stabilize membrane protein structure are poorly understood. To investigate these forces we introduced alanine substitutions at 24 positions in the B helix of bacteriorhodopsin and examined their effects on structure and stability. Although most of the results can be rationalized in terms of the folded structure, there are a number of surprises. (1) We find a remarkably high frequency of stabilizing mutations (17%), indicating that membrane proteins are not highly optimized for stability. (2) Helix B is kinked, with the kink centered around Pro50. The P50A mutation has no effect on stability, however, and a crystal structure reveals that the helix remains bent, indicating that tertiary contacts dominate in the distortion of this helix. (3) We find that the protein is stabilized by about 1kcal/mol for every 38A(2) of surface area buried, which is quite similar to soluble proteins in spite of their dramatically different environments. (4) We find little energetic difference, on average, in the burial of apolar surface or polar surface area, implying that van der Waals packing is the dominant force that drives membrane protein folding.     

Thermal Stabilization of Dihydrofolate Reductase Using Monte Carlo Unfolding Simulations and Its Functional Consequences

Design of proteins with desired thermal properties is important for scientific and biotechnological applications. Here we developed a theoretical approach to predict the effect of mutations on protein stability from non-equilibrium unfolding simulations. We establish a relative measure based on apparent simulated melting temperatures that is independent of simulation length and, under certain assumptions, proportional to equilibrium stability, and we justify this theoretical development with extensive simulations and experimental data. Using our new method based on all-atom Monte-Carlo unfolding simulations, we carried out a saturating mutagenesis of Dihydrofolate Reductase (DHFR), a key target of antibiotics and chemotherapeutic drugs. The method predicted more than 500 stabilizing mutations, several of which were selected for detailed computational and experimental analysis. We find a highly significant correlation of r = 0.65–0.68 between predicted and experimentally determined melting temperatures and unfolding denaturant concentrations for WT DHFR and 42 mutants. The correlation between energy of the native state and experimental denaturation temperature was much weaker, indicating the important role of entropy in protein stability. The most stabilizing point mutation was D27F, which is located in the active site of the protein, rendering it inactive. However for the rest of mutations outside of the active site we observed a weak yet statistically significant positive correlation between thermal stability and catalytic activity indicating the lack of a stability-activity tradeoff for DHFR. By combining stabilizing mutations predicted by our method, we created a highly stable catalytically active E. coli DHFR mutant with measured denaturation temperature 7.2°C higher than WT. Prediction results for DHFR and several other proteins indicate that computational approaches based on unfolding simulations are useful as a general technique to discover stabilizing mutations.     

How Ala-->Gly mutations in different helices affect the stability of the apomyoglobin molten globule

The apomyoglobin molten globule has a complex, partly folded structure with a folded A[B]GH subdomain; the factors determining its stability are not yet known in detail. Ala-->Gly mutations, made at solvent-exposed positions, are used to probe the role of helix propensity of individual helices in stabilizing the molten globule. Molten globule stability is measured by reversible urea unfolding, monitored both by circular dichroism and by tryptophan fluorescence. Two-state unfolding is tested by superposition of these two unfolding curves, and stability data are reported only for variants which satisfy the superposition test. Results for sites Q8 in the A helix and E109 in the G helix confirm that the helix propensities of the A and G helices both strongly affect molten globule stability, in contrast to results for the G65A/G73A double mutant which show that changing the helix propensity of the E-helix sequence has no significant stabilizing effect. Changing the helix propensity of the B-helix sequence with the G23A/G25A double mutant affects molten globule stability to an intermediate extent, confirming an earlier report that this mutant has increased stability. These results are consistent with the bipartite structure for the molten globule in which the A, G, and H helices are stably folded, while the long E helix is unfolded and the B helix has intermediate stability. Some differences are found in the shapes of the unfolding curves of different mutants even though they satisfy the superposition test for two-state unfolding, and possible explanations are discussed.     

Point mutations in transmembrane helices 2 and 3 of ExbB and TolQ affect their activities in Escherichia coli K-12

Replacement of glutamate 176, the only charged amino acid in the third transmembrane helix of ExbB, with alanine (E176A) abolished ExbB activity in all determined ExbB-dependent functions of Escherichia coli. Combination of the mutations T148A in the second transmembrane helix and T181A in the third transmembrane helix, proposed to form part of a proton pathway through ExbB, also resulted in inactive ExbB. E176 and T148 are strictly conserved in ExbB and TolQ proteins, and T181 is almost strictly conserved in ExbB, TolQ, and MotA.     

Relation between protein stability, evolution and structure, as probed by carboxylic acid mutations

Native proteins are marginally stable. Low thermodynamic stability may actually be advantageous, although the accumulation of neutral, destabilizing mutations may have also contributed to it. In any case, once marginal stability has been reached, it appears plausible that mutations at non-constrained positions become fixed in the course of evolution (due to random drift) with frequencies that roughly reflect the mutation effects on stability ("pseudo-equilibrium hypothesis"). We have found that all glutamate-->aspartate mutations in wild-type Escherichia coli thioredoxin are destabilizing, as well as most of the aspartate-->glutamate mutations. Furthermore, the effect of these mutations on thioredoxin thermodynamic stability shows a robust correlation with the frequencies of occurrence of the involved residues in several-hundred sequence alignments derived from a BLAST search. These results provide direct and quantitative experimental evidence for the pseudo-equilibrium hypothesis and should have general consequences for the interpretation of mutation effects on protein stability, as they suggest that residue environments in proteins may be optimized for stabilizing interactions to a remarkable degree of specificity. We also provide evidence that such stabilizing interactions may be detected in sequence alignments, and briefly discuss the implications of this possibility for the derivation of structural information (on native and denatured states) from comparative sequence analyses.     

Engineering and characterization of a stabilized alpha1/alpha2 module of the class I major histocompatibility complex product Ld

The major histocompatibility complex (MHC) is the most polymorphic locus known, with thousands of allelic variants. There is considerable interest in understanding the diversity of structures and peptide-binding features represented by this class of proteins. Although many MHC proteins have been crystallized, others have not been amenable to structural or biochemical studies due to problems with expression or stability. In the present study, yeast display was used to engineer stabilizing mutations into the class I MHC molecule, Ld. The approach was based on previous studies that showed surface levels of yeast-displayed fusion proteins are directly correlated with protein stability. To engineer a more stable Ld, we selected Ld mutants with increased surface expression from randomly mutated yeast display libraries using anti-Ld antibodies or high affinity, soluble T-cell receptors (TCRs). The most stable Ld mutant, Ld-m31, consisted of a single-chain MHC module containing only the alpha1 and alpha2 domains. The enhanced stability was in part due to a single mutation (Trp-97 --> Arg), shown previously to be present in the allele Lq. Mutant Ld-m31 could bind to Ld peptides, and the specific peptide.Ld-m31 complex (QL9.Ld-m31) was recognized by alloreactive TCR 2C. A soluble form of the Ld-m31 protein was expressed in Escherichia coli and refolded from inclusion bodies at high yields. Surface plasmon resonance showed that TCRs bound to peptide.Ld-m31 complexes with affinities similar to those of native full-length Ld. The TCR and QL9.Ld-m31 formed complexes that could be resolved by native gel electrophoresis, suggesting that stabilized alpha1/alpha2 class I platforms may enable various structural studies.     

Dual Role of an N-terminal Amyloidogenic Mutation in Apolipoprotein A-I: DESTABILIZATION OF HELIX BUNDLE AND ENHANCEMENT OF FIBRIL FORMATION*

A number of naturally occurring mutations of apolipoprotein (apo) A-I, the major protein of HDL, are known to be associated with hereditary amyloidosis and atherosclerosis. Here, we examined the effects of the G26R point mutation in apoA-I (apoA-I_Iowa) on the structure, stability, and aggregation propensity to form amyloid fibril of full-length apoA-I and the N-terminal fragment of apoA-I. Circular dichroism and fluorescence measurements demonstrated that the G26R mutation destabilizes the N-terminal helix bundle domain of full-length protein, leading to increased hydrophobic surface exposure, whereas it has no effect on the initial structure of the N-terminal 1–83 fragment, which is predominantly a random coil structure. Upon incubation for extended periods at neutral pH, the N-terminal 1–83 variants undergo a conformational change to β-sheet-rich structure with a great increase in thioflavin T fluorescence, whereas no structural change is observed in full-length proteins. Comparison of fibril-forming propensity among substituted mutants at Gly-26 position of 1–83 fragments demonstrated that the G26R mutation enhances the nucleation step of fibril formation, whereas G26K and G26E mutations have small or inhibiting effects on the formation of fibrils. These fibrils of the 1–83 variants have long and straight morphology as revealed by atomic force microscopy and exhibited significant toxicity with HEK293 cells. Our results indicate dual critical roles of the arginine residue at position 26 in apoA-I_Iowa: destabilization of the N-terminal helix bundle structure in full-length protein and enhancement of amyloid fibril formation by the N-terminal 1–83 fragment.     

Stabilizing proteins from sequence statistics: the interplay of conservation and correlation in triosephosphate isomerase stability

Understanding the determinants of protein stability remains one of protein science's greatest challenges. There are still no computational solutions that calculate the stability effects of even point mutations with sufficient reliability for practical use. Amino acid substitutions rarely increase the stability of native proteins; hence, large libraries and high-throughput screens or selections are needed to stabilize proteins using directed evolution. Consensus mutations have proven effective for increasing stability, but these mutations are successful only about half the time. We set out to understand why some consensus mutations fail to stabilize, and what criteria might be useful to predict stabilization more accurately. Overall, consensus mutations at more conserved positions were more likely to be stabilizing in our model, triosephosphate isomerase (TIM) from Saccharomyces cerevisiae. However, positions coupled to other sites were more likely not to stabilize upon mutation. Destabilizing mutations could be removed both by removing sites with high statistical correlations to other positions and by removing nearly invariant positions at which "hidden correlations" can occur. Application of these rules resulted in identification of stabilizing mutations in 9 out of 10 positions, and amalgamation of all predicted stabilizing positions resulted in the most stable yeast TIM variant we produced (+8 °C). In contrast, a multimutant with 14 mutations each found to stabilize TIM independently was destabilized by 2 °C. Our results are a practical extension to the consensus concept of protein stabilization, and they further suggest the importance of positional independence in the mechanism of consensus stabilization.     

Crystal structure of the human prion protein reveals a mechanism for oligomerization

The pathogenesis of transmissible encephalopathies is associated with the conversion of the cellular prion protein, PrP(C), into a conformationally altered oligomeric form, PrP(Sc). Here we report the crystal structure of the human prion protein in dimer form at 2 A resolution. The dimer results from the three-dimensional swapping of the C-terminal helix 3 and rearrangement of the disulfide bond. An interchain two-stranded antiparallel beta-sheet is formed at the dimer interface by residues that are located in helix 2 in the monomeric NMR structures. Familial prion disease mutations map to the regions directly involved in helix swapping. This crystal structure suggests that oligomerization through 3D domain-swapping may constitute an important step on the pathway of the PrP(C) --> PrP(Sc) conversion.     

The folding state of the lumenal loop determines the thermal stability of light-harvesting chlorophyll a/b protein

The major light-harvesting protein of photosystem II (LHCIIb) is the most abundant chlorophyll-binding protein in the thylakoid membrane. It contains three membrane-spanning alpha helices; the first and third one closely interact with each other to form a super helix, and all three helices bind most of the pigment cofactors. The protein loop domains connecting the alpha helices also play an important role in stabilizing the LHCIIb structure. Single amino acid exchanges in either loop were found to be sufficient to significantly destabilize the complex assembled in vitro [Heinemann, B., and Paulsen, H. (1999) Biochemistry 38, 14088-14093. Mick, V., Eggert, K., Heinemann, B., Geister, S., and Paulsen, H (2004) Biochemistry 43, 5467-5473]. This work presents an analysis of such point mutations in the lumenal loop with regard to the extent and nature of their effect on LHCIIb stability to obtain detailed information on the contribution of this loop to stabilizing the complex. Most of the mutant proteins yielded pigment-protein complexes if their reconstitution and/or isolation was performed under mild conditions; however, the yields were significantly different. Several mutations in the vicinity of W97 in the N-proximal section of the loop gave low reconstitution yields even under very mild conditions. This confirms our earlier notion that W97 may be of particular relevance in stabilizing LHCIIb. The same amino acid exchanges accelerated thermal complex dissociation in the absence of lithium dodecyl sulfate (LDS) and raised the accessibility of the lumenal loop to protease; both effects were well correlated with the reduction in reconstitution yields. We conclude that a detachment of the lumenal loop is a possible first step in the dissociation of LHCIIb. Dramatically reduced complex yields in the presence but not in the absence of LDS were observed for some but not all mutants, particularly those near the C-proximal end of the loop. We conclude that complex stabilities in the absence and in the presence of LDS do not correlate and most likely are determined by different structural characteristics, at least in LHCIIb but maybe also in other membrane proteins.     

How protein stability and new functions trade off

Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (DeltaDeltaG), thus suggesting a tradeoff between the acquisition of new enzymatic functions and stability. However, since most mutations are destabilizing (DeltaDeltaG>0), one should ask how destabilizing mutations that confer new or altered enzymatic functions relative to all other mutations are. We applied DeltaDeltaG computations by FoldX to analyze the effects of 548 mutations that arose from the directed evolution of 22 different enzymes. The stability effects, location, and type of function-altering mutations were compared to DeltaDeltaG changes arising from all possible point mutations in the same enzymes. We found that mutations that modulate enzymatic functions are mostly destabilizing (average DeltaDeltaG = +0.9 kcal/mol), and are almost as destabilizing as the "average" mutation in these enzymes (+1.3 kcal/mol). Although their stability effects are not as dramatic as in key catalytic residues, mutations that modify the substrate binding pockets, and thus mediate new enzymatic specificities, place a larger stability burden than surface mutations that underline neutral, non-adaptive evolutionary changes. How are the destabilizing effects of functional mutations balanced to enable adaptation? Our analysis also indicated that many mutations that appear in directed evolution variants with no obvious role in the new function exert stabilizing effects that may compensate for the destabilizing effects of the crucial function-altering mutations. Thus, the evolution of new enzymatic activities, both in nature and in the laboratory, is dependent on the compensatory, stabilizing effect of apparently "silent" mutations in regions of the protein that are irrelevant to its function.