Dissection of protein structure and folding by directed mutagenesis

pPerturb: A Server for Predicting Long-Distance Energetic Couplings and Mutation-Induced Stability Changes in Proteins via Perturbations

The strength of intraprotein interactions or contact network is one of the dominant factors determining the thermodynamic stabilities of proteins. The nature and the extent of connectivity of this network also play a role in allosteric signal propagation characteristics upon ligand binding to a protein domain. Here, we develop a server for rapid quantification of the strength of an interaction network by employing an experimentally consistent perturbation approach previously validated against a large data set of 375 mutations in 19 different proteins. The web server can be employed to predict the extent of destabilization of proteins arising from mutations in the protein interior in experimentally relevant units. Moreover, coupling distances-a measure of the extent of percolation on perturbation-and overall perturbation magnitudes are predicted in a residue-specific manner, enabling a first look at the distribution of energetic couplings in a protein or its changes upon ligand binding. We show specific examples of how the server can be employed to probe for the distribution of local stabilities in a protein, to examine changes in side chain orientations or packing before and after ligand binding, and to predict changes in stabilities of proteins upon mutations of buried residues. The web server is freely available at http://pbl.biotech.iitm.ac.in/pPerturb and supports recent versions of all major browsers.

Engineering Order and Cooperativity in a Disordered Protein

Structural disorder in proteins arises from a complex interplay between weak hydrophobicity and unfavorable electrostatic interactions. The extent to which the hydrophobic effect contributes to the unique and compact native state of proteins is, however, confounded by large compensation between multiple entropic and energetic terms. Here we show that protein structural order and cooperativity arise as emergent properties upon hydrophobic substitutions in a disordered system with non-intuitive effects on folding and function. Aided by sequence-structure analysis, equilibrium, and kinetic spectroscopic studies, we engineer two hydrophobic mutations in the disordered DNA-binding domain of CytR that act synergistically, but not in isolation, to promote structure, compactness, and stability. The double mutant, with properties of a fully ordered domain, exhibits weak cooperativity with a complex and rugged conformational landscape. The mutant, however, binds cognate DNA with an affinity only marginally higher than that of the wild type, though nontrivial differences are observed in the binding to noncognate DNA. Our work provides direct experimental evidence of the dominant role of non-additive hydrophobic effects in shaping the molecular evolution of order in disordered proteins and vice versa, which could be generalized to even folded proteins with implications for protein design and functional manipulation.

Nature's Shortcut to Protein Folding

This Feature Article presents a view of the protein folding transition based on the hypothesis that Nature has built features within the sequences that enable a Shortcut to efficient folding. Nature's Shortcut is proposed to be the early establishment of a set of nonlocal weak contacts, constituting protein loops that significantly constrain regions of the collapsed disordered protein into a native-like low-resolution fluctuating topology of major sections of the backbone. Nature's establishment of this scaffold of nonlocal contacts is claimed to bypass what would otherwise be a nearly hopeless unaided search for the final three-dimensional structure in proteins longer than ∼100 amino acids. To support this main contention of the Feature Article, the loop hypothesis (LH) description of early folding events is experimentally tested with time-resolved Förster resonance energy transfer techniques for adenylate kinase, and the data are shown to be consistent with theoretical predictions from the sequential collapse model (SCM). The experimentally based LH and the theoretically founded SCM are argued to provide a unified picture of the role of nonlocal contacts as constituting Nature's Shortcut to protein folding. Importantly, the SCM is shown to reliably predict key nonlocal contacts utilizing only primary sequence information. This view on Nature's Shortcut is open to the protein community for further detailed assessment, including its practical consequences, by suitable application of advanced experimental and computational techniques.

Protein evolution revisited

Antifreeze proteins (AFPs) protect marine fishes from freezing in icy seawater. They evolved relatively recently, most likely in response to the formation of sea ice and Cenozoic glaciations that occurred less than 50 million years ago, following a greenhouse Earth event. Based on their diversity, AFPs have independently evolved on many occasions to serve the same function, with some remarkable examples of convergent evolution at the structural level, and even instances of lateral gene transfer. For some AFPs, the progenitor gene is recognizable. The intense selection pressure exerted by icy seawater, which can rapidly kill unprotected fish, has led to massive AFP gene amplification, as well as some partial gene duplications that have increased the size and activity of the antifreeze. The many protein evolutionary processes described in Gordon H. Dixon's Essays in Biochemistry article will be illustrated here by examples from studies on AFPs. Abbreviations: AFGP: antifreeze glycoproteins; AFP: antifreeze proteins; GHD: Gordon H. Dixon; SAS: sialic acid synthase; TH: thermal hysteresis.

In silico mutational analysis and identification of stability centers in human interleukin-4

Interleukin-4 (IL-4) is a multifunctional cytokine that plays a critical role in apoptosis, differentiation and proliferation. The intensity of IL4 response depends upon binding to its receptor, IL-4R. The therapeutic efficiency of interleukins can be increased by generating structural mutants having greater stability. In the present work, attempts were made to increase the stability of human IL-4 using in-silico site directed mutagenesis. Different orthologous sequences of IL4 from Pan troglodytes, Aotusnigriceps, Macacamulatta, Papiohamadryas, Chlorocebusaethiops, Vicugnapacos, Susscrofa and Homo sapiens were aligned using Clustal Omega that revealed the conserved and non-conserved positions. For each non-conserved position, possible favorable and stabilizing mutations were found using CUPSAT with predicted ΔΔG (kcal/mol). The one with highest ΔΔG (kcal/mol) among all possible mutations, for each non-conserved position was selected and introduced manually in human IL-4 sequence resulting in multiple mutants of IL-4. Mutant proteins were modeled using structure of IL4 (PDB ID: 2B8U) as a template by SWISS MODEL. The mutants A49L and Q106T were identified to have stability centre using SCide. Molecular dynamics and docking analysis also confirmed the mutants stability and binding respectively. Mutants A49L and Q106T had -7.580079 kcal/mol and -39.418124 kcal/mol respectively lesser energy value than the wild type IL4. The result suggested that, the stability of human IL-4 has been increased by mutation.

Dissecting protein function: an efficient protocol for identifying separation-of-function mutations that encode structurally stable proteins

Mutations that confer the loss of a single biochemical property (separation-of-function mutations) can often uncover a previously unknown role for a protein in a particular biological process. However, most mutations are identified based on loss-of-function phenotypes, which cannot differentiate between separation-of-function alleles vs. mutations that encode unstable/unfolded proteins. An alternative approach is to use overexpression dominant-negative (ODN) phenotypes to identify mutant proteins that disrupt function in an otherwise wild-type strain when overexpressed. This is based on the assumption that such mutant proteins retain an overall structure that is comparable to that of the wild-type protein and are able to compete with the endogenous protein (Herskowitz 1987). To test this, the in vivo phenotypes of mutations in the Est3 telomerase subunit from Saccharomyces cerevisiae were compared with the in vitro secondary structure of these mutant proteins as analyzed by circular-dichroism spectroscopy, which demonstrates that ODN is a more sensitive assessment of protein stability than the commonly used method of monitoring protein levels from extracts. Reverse mutagenesis of EST3, which targeted different categories of amino acids, also showed that mutating highly conserved charged residues to the oppositely charged amino acid had an increased likelihood of generating a severely defective est3(-) mutation, which nevertheless encoded a structurally stable protein. These results suggest that charge-swap mutagenesis directed at a limited subset of highly conserved charged residues, combined with ODN screening to eliminate partially unfolded proteins, may provide a widely applicable and efficient strategy for generating separation-of-function mutations.

Thermal precipitation fluorescence assay for protein stability screening

A simple and reliable method of protein stability assessment is desirable for high throughput expression screening of recombinant proteins. Here we described an assay termed thermal precipitation fluorescence (TPF) which can be used to compare thermal stabilities of recombinant protein samples directly from cell lysate supernatants. In this assay, target membrane proteins are expressed as recombinant fusions with a green fluorescence protein tag and solubilized with detergent, and the fluorescence signals are used to report the quantity of the fusion proteins in the soluble fraction of the cell lysate. After applying a heat shock, insoluble protein aggregates are removed by centrifugation. Subsequently, the amount of remaining protein in the supernatant is quantified by in-gel fluorescence analysis and compared to samples without a heat shock treatment. Over 60 recombinant membrane proteins from Escherichia coli were subject to this screening in the presence and absence of a few commonly used detergents, and the results were analyzed. Because no sophisticated protein purification is required, this TPF technique is suitable to high throughput expression screening of recombinant membrane proteins as well as soluble ones and can be used to prioritize target proteins based on their thermal stabilities for subsequent large scale expression and structural studies.

Membrane topology and identification of critical amino acid residues in the Wzx O-antigen translocase from Escherichia coli O157:H4

Wzx belongs to a family of membrane proteins involved in the translocation of isoprenoid lipid-linked glycans, which is loosely related to members of the major facilitator superfamily. Despite Wzx homologs performing a conserved function, it has been difficult to pinpoint specific motifs of functional significance in their amino acid sequences. Here, we elucidate the topology of the Escherichia coli O157 Wzx (Wzx(EcO157)) by a combination of bioinformatics and substituted cysteine scanning mutagenesis, as well as targeted deletion-fusions to green fluorescent protein and alkaline phosphatase. We conclude that Wzx(EcO157) consists of 12 transmembrane (TM) helices and six periplasmic and five cytosolic loops, with N and C termini facing the cytoplasm. Four TM helices (II, IV, X, and XI) contain polar residues (aspartic acid or lysine), and they may form part of a relatively hydrophilic core. Thirty-five amino acid replacements to alanine or serine were targeted to five native cysteines and most of the aspartic acid, arginine, and lysine residues. From these, only replacements of aspartic acid-85, aspartic acid-326, arginine-298, and lysine-419 resulted in a protein unable to support O-antigen production. Aspartic acid-85 and lysine-419 are located in TM helices II and XI, while arginine-298 and aspartic acid-326 are located in periplasmic and cytosolic loops 4, respectively. Further analysis revealed that the charge at these positions is required for Wzx function since conservative substitutions maintaining the same charge polarity resulted in a functional protein, whereas those reversing or eliminating polarity abolished function. We propose that the functional requirement of charged residues at both sides of the membrane and in two TM helices could be important to allow the passage of the Und-PP-linked saccharide substrate across the membrane.

XylS-Pm promoter interactions through two helix-turn-helix motifs: identifying XylS residues important for DNA binding and activation

The XylS protein is the positive transcription regulator of the TOL plasmid meta-cleavage pathway operon Pm. XylS belongs to the AraC family of transcriptional regulators and exhibits an N-terminal domain involved in effector recognition, and a C-terminal domain, made up of seven alpha-helices conforming two helix-turn-helix DNA-binding domains. alpha-Helix 3 and alpha-helix 6 are the recognition helices. In consonance with XylS structural organization, Pm exhibits a bipartite DNA-binding motif consisting of two boxes, called A and B, whose sequences are TGCA and GGNTA, respectively. This bipartite motif is repeated at the Pm promoter so that one of the XylS monomers binds to each of the repeats. An extensive series of genetic epistasis assays combining mutant Pm promoters and XylS single substitution mutant proteins revealed that alpha-helix 3 contacts A box nucleotides, whereas alpha-helix 6 residues contact B box nucleotides. In alpha-helix 3, Asn246 and Arg242 are involved in specific contacts with the TG dinucleotide at box A, whereas Arg296 and Glu299 contact the second G and T nucleotides at box B. On the basis of our results and of the three-dimensional model of the XylS C-terminal domain, we propose that Ser243, Glu249 and Lys250 in alpha-helix 3, and Asn299 and Arg302 in alpha-helix 6 contact the phosphate backbones. Alanine substitutions at the predicted phosphate backbone-contacting residues yielded mutants with low levels of activity, suggesting that XylS-Pm binding specificity not only involves specific amino acid-base interactions, but also relies on secondary DNA structure, which, although at another level, also comes into play. We propose a model in which a XylS dimer binds to the direct repeats in Pm in a head-to-tail conformation that allows the direct interaction of the XylS proximal subunit with the RNA polymerase sigma factor.

SIMPLE estimate of the free energy change due to aliphatic mutations: Superior predictions based on first principles

The bioinformatics revolution of the last decade has been instrumental in the development of empirical potentials to quantitatively estimate protein interactions for modeling and design. Although computationally efficient, these potentials hide most of the relevant thermodynamics in 5-to-40 parameters that are fitted against a large experimental database. Here, we revisit this longstanding problem and show that a careful consideration of the change in hydrophobicity, electrostatics, and configurational entropy between the folded and unfolded state of aliphatic point mutations predicts 20-30% less false positives and yields more accurate predictions than any published empirical energy function. This significant improvement is achieved with essentially no free parameters, validating past theoretical and experimental efforts to understand the thermodynamics of protein folding. Our first principle analysis strongly suggests that both the solute-solute van der Waals interactions in the folded state and the electrostatics free energy change of exposed aliphatic mutations are almost completely compensated by similar interactions operating in the unfolded ensemble. Not surprisingly, the problem of properly accounting for the solvent contribution to the free energy of polar and charged group mutations, as well as of mutations that disrupt the protein backbone remains open.

Energetics of aliphatic deletions in protein cores

Although core residues can sometimes be replaced by shorter ones without introducing significant changes in protein structure, the energetic consequences are typically large and destabilizing. Many efforts have been devoted to understand and predict changes in stability from analysis of the environment of mutated residues, but the relationships proposed for individual proteins have often failed to describe additional data. We report here 17 apoflavodoxin large-to-small mutations that cause overall protein destabilizations of 0.6-3.9 kcal.mol(-1). By comparing two-state urea and three-state thermal unfolding data, the overall destabilizations observed are partitioned into effects on the N-to-I and on the I-to-U equilibria. In all cases, the equilibrium intermediate exerts a "buffering" effect that reduces the impact of the overall destabilization on the N-to-I equilibrium. The performance of several structure-energetics relationships, proposed to explain the energetics of hydrophobic shortening mutations, has been evaluated by using an apoflavodoxin data set consisting of 14 mutations involving branching-conservative aliphatic side-chain shortenings and a larger data set, including similar mutations implemented in seven model proteins. Our analysis shows that the stability changes observed for any of the different types of mutations (LA, IA, IV, and VA) in either data set are best explained by a combination of differential hydrophobicity and of the calculated volume of the modeled cavity (as previously observed for LA and IA mutations in lysozyme T4). In contrast, sequence conservation within the flavodoxin family, which is a good predictor for charge-reversal stabilizing mutations, does not perform so well for aliphatic shortening ones.

Comprehensive alanine-scanning mutagenesis of Escherichia coli CsrA defines two subdomains of critical functional importance

The RNA-binding protein CsrA (carbon storage regulator) of Escherichia coli is a global regulator of gene expression and is representative of the CsrA/RsmA family of bacterial proteins. These proteins act by regulating mRNA translation and stability and are antagonized by binding to small noncoding RNAs. Although the RNA target sequence and structure for CsrA binding have been well defined, little information exists concerning the protein requirements for RNA recognition. The three-dimensional structures of three CsrA/RsmA proteins were recently solved, revealing a novel protein fold consisting of two interdigitated monomers. Here, we performed comprehensive alanine-scanning mutagenesis on csrA of E. coli and tested the 58 resulting mutants for regulation of glycogen accumulation, motility, and biofilm formation. Quantitative effects of these mutations on expression of glgCA'-'lacZ, flhDC'-'lacZ, and pgaA'-'lacZ translational fusions were also examined, and eight of the mutant proteins were purified and tested for RNA binding. These studies identified two regions of the amino acid sequence that were critical for regulation and RNA binding, located within the first (beta1, residues 2-7) and containing the last (beta5, residues 40-47) beta-strands of CsrA. The beta1 and beta5 strands of opposite monomers lie adjacent and parallel to each other in the three-dimensional structure of this protein. Given the symmetry of the CsrA dimer, these findings imply that two distinct RNA binding surfaces or functional subdomains lie on opposite sides of the protein.

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.

Solvation in protein folding analysis: combination of theoretical and experimental approaches

An effort to combine theoretical analyses and protein engineering methods has been made to probe the folding mechanism of SH3 by using Energy Landscape Theory and a phi-value analysis. Particular emphasis was given to core residues and the effect of desolvation during the folding event by replacing the core valines with isosteric threonines. These mutations have the advantage of keeping the core structurally invariant while affecting core stability relative to the unfolded state. Although the valines that form the core appear spatially invariant, the folding kinetics of their threonine mutants varies, indicating their different extent of solvation in the transition-state ensemble. Theoretical studies predicted the distribution of folding kinetics of threonine mutants without previous knowledge of the measured rates. This initial success encourages further investigations of the molecular details behind these macroscopic phenomena and of the role of solvation in the folding mechanism.

Non-standard insulin design: structure-activity relationships at the periphery of the insulin receptor

The design of insulin analogues has emphasized stabilization or destabilization of structural elements according to established principles of protein folding. To this end, solvent-exposed side-chains extrinsic to the receptor-binding surface provide convenient sites of modification. An example is provided by an unfavorable helical C-cap (Thr(A8)) whose substitution by favorable amino acids (His(A8) or Arg(A8)) has yielded analogues of improved stability. Remarkably, these analogues also exhibit enhanced activity, suggesting that activity may correlate with stability. Here, we test this hypothesis by substitution of diaminobutyric acid (Dab(A8)), like threonine an amino acid of low helical propensity. The crystal structure of Dab(A8)-insulin is similar to those of native insulin and the related analogue Lys(A8)-insulin. Although no more stable than native insulin, the non-standard analogue is twice as active. Stability and affinity can therefore be uncoupled. To investigate alternative mechanisms by which A8 substitutions enhance activity, multiple substitutions were introduced. Surprisingly, diverse aliphatic, aromatic and polar side-chains enhance receptor binding and biological activity. Because no relationship is observed between activity and helical propensity, we propose that local interactions between the A8 side-chain and an edge of the hormone-receptor interface modulate affinity. Dab(A8)-insulin illustrates the utility of non-standard amino acids in hypothesis-driven protein design.

Tolerance of Arc repressor to multiple-alanine substitutions

Arc repressor mutants containing from three to 15 multiple-alanine substitutions have spectral properties expected for native Arc proteins, form heterodimers with wild-type Arc, denature cooperatively with Tms equal to or greater than wild type, and, in some cases, fold as much as 30-fold faster and unfold as much as 50-fold slower than wild type. Two of the mutants, containing a total of 14 different substitutions, also footprint operator DNA in vitro. The stability of some of the proteins with multiple-alanine mutations is significantly greater than that predicted from the sum of the single substitutions, suggesting that a subset of the wild-type residues in Arc may interact in an unfavorable fashion. Overall, these results show that almost half of the residues in Arc can be replaced by alanine en masse without compromising the ability of this small, homodimeric protein to fold into a stable, native-like structure.

The native state of apomyoglobin described by proton NMR spectroscopy: the A-B-G-H interface of wild-type sperm whale apomyoglobin

Proton nuclear magnetic resonance spectroscopy was applied to sperm whale apomyoglobin to describe the conformation adopted by the protein under native conditions. The study focused on the A-B-G-H interface, a region known to form a compact subdomain in the apoprotein (Hughson and Baldwin, Biochemistry 28:4415-4422, 1989). Two histidine residues located in this subdomain, His24 and His119, interact and are thought to play a role in the acid denaturation process (Barrick et al., J. Mol. Biol. 237:588-601, 1994). A stable double mutant at these positions (His24Val/His119Phe sperm whale apomyoglobin) was compared with wild-type apomyoglobin. The amino acid replacements result in chemical shift perturbations near the mutations, in particular in the AB interhelical region, and in a deceleration of backbone amide hydrogen exchange in the B helix from position 27 to position 33. The double mutant data were used to expand and confirm the wild-type spectral analysis. Signals from the D helix were identified that demonstrate the formation of holoprotein-like structure. The assigned wild-type nuclear Overhauser effects, although in small number, were sufficient to construct a model of the compact subdomain of the apoprotein. This was achieved by using the structure of the holoprotein and restraining it with the geometrical information on the apoprotein in a simulated annealing procedure. The experimental restraints define a low-resolution model of the A-B-G-H interface in apomyoglobin.

Stability and function: two constraints in the evolution of barstar and other proteins

BACKGROUND

Barstar is the intracellular inhibitor of barnase, an extracellular RNAse of Bacillus amyloliquefaciens. The dissociation constant of the barnase-barstar complex is 10(-14) M with an association rate constant between barnase and barstar of 3.7 x 10(8) s-1 M-1. The rapid association arises in part from the clustering of four acidic residues (Asp35, Asp39, Glu76 and Glu80) on the barnase-binding surface of barstar. The negatively charged barnase-binding surface of barstar effectively 'steers' the inhibitor towards the positively charged active site of barnase.

RESULTS

Mutating any one of the four acidic side chains of barstar to an alanine results in an approximately two-fold decrease in the association rate constant, while the dissociation rate constant increases from five orders of magnitude for Asp39-->Ala, to no significant change for Glu80-->Ala. The stability of barstar is increased by all four mutations, the increase ranging from 0.3 kcal mol-1 for Asp35-->Ala or Asp39-->Ala, to 2.1 kcal mol-1 for Glu80-->Ala.

CONCLUSIONS

The evolutionary pressure on barstar for rapid binding of barnase is so strong that glutamate is preferred over alanine at position 80, even though it does not directly interact with barnase in the complex and significantly destabilizes the inhibitor structure. This, and other examples from the literature, suggest that proteins evolve primarily to optimize their function in vivo, with relatively little evolutionary pressure to increase stability above a certain threshold, thus allowing greater latitude in the evolution of enzyme activity.