Increasing protein stability using a rational approach combining sequence homology and structural alignment: Stabilizing the WW domain

Mega-scale experimental analysis of protein folding stability in biology and design

Advances in DNA sequencing and machine learning are providing insights into protein sequences and structures on an enormous scale1. However, the energetics driving folding are invisible in these structures and remain largely unknown2. The hidden thermodynamics of folding can drive disease3,4, shape protein evolution5-7 and guide protein engineering8-10, and new approaches are needed to reveal these thermodynamics for every sequence and structure. Here we present cDNA display proteolysis, a method for measuring thermodynamic folding stability for up to 900,000 protein domains in a one-week experiment. From 1.8 million measurements in total, we curated a set of around 776,000 high-quality folding stabilities covering all single amino acid variants and selected double mutants of 331 natural and 148 de novo designed protein domains 40-72 amino acids in length. Using this extensive dataset, we quantified (1) environmental factors influencing amino acid fitness, (2) thermodynamic couplings (including unexpected interactions) between protein sites, and (3) the global divergence between evolutionary amino acid usage and protein folding stability. We also examined how our approach could identify stability determinants in designed proteins and evaluate design methods. The cDNA display proteolysis method is fast, accurate and uniquely scalable, and promises to reveal the quantitative rules for how amino acid sequences encode folding stability.

CovET: A covariation-evolutionary trace method that identifies protein structure-function modules

Measuring the relative effect that any two sequence positions have on each other may improve protein design or help better interpret coding variants. Current approaches use statistics and machine learning but rarely consider phylogenetic divergences which, as shown by Evolutionary Trace studies, provide insight into the functional impact of sequence perturbations. Here, we reframe covariation analyses in the Evolutionary Trace framework to measure the relative tolerance to perturbation of each residue pair during evolution. This approach (CovET) systematically accounts for phylogenetic divergences: at each divergence event, we penalize covariation patterns that belie evolutionary coupling. We find that while CovET approximates the performance of existing methods to predict individual structural contacts, it performs significantly better at finding structural clusters of coupled residues and ligand binding sites. For example, CovET found more functionally critical residues when we examined the RNA recognition motif and WW domains. It correlates better with large-scale epistasis screen data. In the dopamine D2 receptor, top CovET residue pairs recovered accurately the allosteric activation pathway characterized for Class A G protein-coupled receptors. These data suggest that CovET ranks highest the sequence position pairs that play critical functional roles through epistatic and allosteric interactions in evolutionarily relevant structure-function motifs. CovET complements current methods and may shed light on fundamental molecular mechanisms of protein structure and function.

Computational and synthetic biology approaches for the biosynthesis of antiviral and anticancer terpenoids from Bacillus subtilis

Recent advancements in medicinal research have identified several antiviral and anticancer terpenoids that are usually deployed as a source of flavor, fragrances and pharmaceuticals. Under the current COVID-19 pandemic conditions, natural therapeutics with least side effects are the need of the hour to save the patients, especially, which are pre-affected with other medical complications. Although, plants are the major sources of terpenoids; however, for the environmental concerns, the global interest has shifted to the biocatalytic production of molecules from microbial sources. The gram-positive bacterium Bacillus subtilis is a suitable host in this regard due to its GRAS (generally regarded as safe) status, ease in genetic manipulations and wide industrial acceptability. The B. subtilis synthesizes its terpenoid molecules from 1-deoxy-d-xylulose-5-phosphate (DXP) pathway, a common route in almost all microbial strains. Here, we summarize the computational and synthetic biology approaches to improve the production of terpenoid-based therapeutics from B. subtilis by utilizing DXP pathway. We focus on the in-silico approaches for screening the functionally improved enzyme-variants of the two crucial enzymes namely, the DXP synthase (DXS) and farnesyl pyrophosphate synthase (FPPS). The approaches for engineering the active sites are subsequently explained. It will be helpful to construct the functionally improved enzymes for the high-yield production of terpenoid-based anticancer and antiviral metabolites, which would help to reduce the cost and improve the availability of such therapeutics for the humankind.

Consensus mutagenesis approach improves the thermal stability of system xc - transporter, xCT, and enables cryo-EM analyses

System x_c⁻ is an amino acid antiporter that imports L‐cystine into cells and exports intracellular L‐glutamate, at a 1:1 ratio. As L‐cystine is an essential precursor for glutathione synthesis, system x_c⁻ supports tumor cell growth through glutathione‐based oxidative stress resistance and is considered as a potential therapeutic target for cancer treatment. System x_c⁻ consists of two subunits, the light chain subunit SLC7A11 (xCT) and the heavy chain subunit SLC3A2 (also known as CD98hc or 4F2hc), which are linked by a conserved disulfide bridge. Although the recent structures of another SLC7 member, L‐type amino acid transporter 1 (LAT1) in complex with CD98hc, have provided the structural basis toward understanding the amino acid transport mechanism, the detailed molecular mechanism of xCT remains unknown. To revealthe molecular mechanism, we performed single‐particle analyses of the xCT‐CD98hc complex. As wild‐type xCT‐CD98hc displayed poor stability and could not be purified to homogeneity, we applied a consensus mutagenesis approach to xCT. The consensus mutated construct exhibited increased stability as compared to the wild‐type, and enabled the cryoelectron microscopy (cryo‐EM) map to be obtained at 6.2 Å resolution by single‐particle analysis. The cryo‐EM map revealed sufficient electron density to assign secondary structures. In the xCT structure, the hash and arm domains are well resolved, whereas the bundle domain shows some flexibility. CD98hc is positioned next to the xCT transmembrane domain. This study provides the structural basis of xCT, and our consensus‐based strategy could represent a good choice toward solving unstable protein structures.

Convolution Neural Network-Based Prediction of Protein Thermostability

Most natural proteins exhibit poor thermostability, which limits their industrial application. Computer-aided rational design is an efficient purpose-oriented method that can improve protein thermostability. Numerous machine-learning-based methods have been designed to predict the changes in protein thermostability induced by mutations. However, all of these methods have certain limitations due to existing mutation coding methods that overlook protein sequence features. Here we propose a method to predict protein thermostability using convolutional neural networks based on an in-depth study of thermostability-related protein properties. This method comprises a three-dimensional coding algorithm, including protein mutation information and a strategy to extract neighboring features at protein mutation sites based on multiscale convolution. The accuracies on the S1615 and S388 data sets, which are widely used for protein thermostability predictions, reached 86.4 and 87%, respectively. The Matthews correlation coefficient was nearly double those produced using other methods. Furthermore, a model was constructed to predict the thermostability of Rhizomucor miehei lipase mutants based on the S3661 data set, a single amino acid mutation data set screened from the ProTherm protein thermodynamics database. Compared with the RIF strategy, which consists of three algorithms, i.e., Rosetta ddg monomer, I Mutant 3.0, and FoldX, the accuracy of the proposed method was higher (75.0 vs 66.7%), and the negative sample resolution was simultaneously enhanced. These results indicate that our prediction method more effectively assessed the protein thermostability and distinguished its features, making it a powerful tool to devise mutations that enhance the thermostability of proteins, particularly enzymes.

Global analysis of protein folding using massively parallel design, synthesis, and testing

Proteins fold into unique native structures stabilized by thousands of weak interactions that collectively overcome the entropic cost of folding. Although these forces are "encoded" in the thousands of known protein structures, "decoding" them is challenging because of the complexity of natural proteins that have evolved for function, not stability. We combined computational protein design, next-generation gene synthesis, and a high-throughput protease susceptibility assay to measure folding and stability for more than 15,000 de novo designed miniproteins, 1000 natural proteins, 10,000 point mutants, and 30,000 negative control sequences. This analysis identified more than 2500 stable designed proteins in four basic folds-a number sufficient to enable us to systematically examine how sequence determines folding and stability in uncharted protein space. Iteration between design and experiment increased the design success rate from 6% to 47%, produced stable proteins unlike those found in nature for topologies where design was initially unsuccessful, and revealed subtle contributions to stability as designs became increasingly optimized. Our approach achieves the long-standing goal of a tight feedback cycle between computation and experiment and has the potential to transform computational protein design into a data-driven science.

Electrostatic effects on the folding stability of FKBP12

The roles of electrostatic interactions in protein folding stability have been a matter of debate, largely due to the complexity in the theoretical treatment of these interactions. We have developed computational methods for calculating electrostatic effects on protein folding stability. To rigorously test and further refine these methods, here we carried out experimental studies into electrostatic effects on the folding stability of the human 12-kD FK506 binding protein (FKBP12). This protein has a close homologue, FKBP12.6, with amino acid substitutions in only 18 of their 107 residues. Of the 18 substitutions, 8 involve charged residues. Upon mutating FKBP12 residues at these 8 positions individually into the counterparts in FKBP12.6, the unfolding free energy (ΔGu) of FKBP12 changed by -0.3 to 0.7 kcal/mol. Accumulating stabilizing substitutions resulted in a mutant with a 0.9 kcal/mol increase in stability. Additional charge mutations were grafted from a thermophilic homologue, MtFKBP17, which aligns to FKBP12 with 31% sequence identity over 89 positions. Eleven such charge mutations were studied, with ΔΔGu varying from -2.9 to 0.1 kcal/mol. The predicted electrostatic effects by our computational methods with refinements herein had a root-mean-square deviation of 0.9 kcal/mol from the experimental ΔΔGu values on 16 single mutations of FKBP12. The difference in ΔΔGu between mutations grafted from FKBP12.6 and those from MtFKBP17 suggests that more distant homologues are less able to provide guidance for enhancing folding stability.

Consensus protein design

A popular and successful strategy in semi-rational design of protein stability is the use of evolutionary information encapsulated in homologous protein sequences. Consensus design is based on the hypothesis that at a given position, the respective consensus amino acid contributes more than average to the stability of the protein than non-conserved amino acids. Here, we review the consensus design approach, its theoretical underpinnings, successes, limitations and challenges, as well as providing a detailed guide to its application in protein engineering.

Autoinhibitory structure of the WW domain of HYPB/SETD2 regulates its interaction with the proline-rich region of huntingtin

Huntington's disease (HD) is an autosomally dominant neurodegenerative disorder caused by expansion of polyglutamine (polyQ) in the huntingtin (Htt) protein. Htt yeast two-hybrid protein B (HYPB/SETD2), a histone methyltransferase, directly interacts with Htt and is involved in HD pathology. Using NMR techniques, we characterized a polyproline (polyP) stretch at the C terminus of HYPB, which directly interacts with the following WW domain and leads this domain predominantly to be in a closed conformational state. The solution structure shows that the polyP stretch extends from the back and binds to the WW core domain in a typical binding mode. This autoinhibitory structure regulates interaction between the WW domain of HYPB and the proline-rich region (PRR) of Htt, as evidenced by NMR and immunofluorescence techniques. This work provides structural and mechanistic insights into the intramolecular regulation of the WW domain in Htt-interacting partners and will be helpful for understanding the pathology of HD.

Engineering and kinetic stabilization of the therapeutic enzyme Anabeana variabilis phenylalanine ammonia lyase

Anabeana variabilis phenylalanine ammonia lyase has just recently been discovered and introduced in clinical trials of phenylketonuria enzyme replacement therapy for its outstanding kinetic properties. In the present study, kinetic stabilization of this therapeutically important enzyme has been explored by introduction of a disulfide bond into the structure. Site-directed mutagenesis was performed with quick-change PCR method. Recombinant wild-type and mutated enzymes were expressed in Escherichia coli, and his-tagged proteins were affinity purified. Formation of disulfide bond was confirmed by Ellman's method, and then chemical unfolding, kinetic behavior, and thermal inactivation of mutated enzyme were compared with the wild type. Based on our results, the Q292C mutation resulted in a significant improvement in kinetic stability and resistance against chemical unfolding of the enzyme while kinetic parameters and pH profile of enzyme activity were remained unaffected. The results of the present study provided an insight towards designing phenylalanine ammonia lyases with higher stability.

Selection of a high-affinity WW domain against the extracellular region of VEGF receptor isoform-2 from a combinatorial library using CIS display

WW domains are small β-sheet motifs that are involved in intracellular signalling through the recognition of proline-rich or phosphorylated linear peptide sequences. Here, we describe modification of this motif to provide a framework for engineering the side chains exposed on its concave surface. This non-natural scaffold incorporates an additional tryptophan, has a shorter loop 1 and supports modification of 25% of the natural protein to form a novel affinity reagent. We demonstrate the utility of this structure by selecting a high-affinity binder to the extracellular region of human vascular endothelial growth factor receptor isoform 2 (VEGFR-2) from a library of modifications, using a cell-free molecular display platform, CIS display. The isolate has low nanomolar affinity to VEGFR-2 and inhibits binding of human VEGF to its receptor with nanomolar activity. The structure is amenable to cyclisation to improve its proteolytic stability and has advantages over larger protein scaffolds in that it can be synthesised chemically to high yields offering potential for therapeutic and non-therapeutic applications.

A fundamental protein property, thermodynamic stability, revealed solely from large-scale measurements of protein function

The ability of a protein to carry out a given function results from fundamental physicochemical properties that include the protein's structure, mechanism of action, and thermodynamic stability. Traditional approaches to study these properties have typically required the direct measurement of the property of interest, oftentimes a laborious undertaking. Although protein properties can be probed by mutagenesis, this approach has been limited by its low throughput. Recent technological developments have enabled the rapid quantification of a protein's function, such as binding to a ligand, for numerous variants of that protein. Here, we measure the ability of 47,000 variants of a WW domain to bind to a peptide ligand and use these functional measurements to identify stabilizing mutations without directly assaying stability. Our approach is rooted in the well-established concept that protein function is closely related to stability. Protein function is generally reduced by destabilizing mutations, but this decrease can be rescued by stabilizing mutations. Based on this observation, we introduce partner potentiation, a metric that uses this rescue ability to identify stabilizing mutations, and identify 15 candidate stabilizing mutations in the WW domain. We tested six candidates by thermal denaturation and found two highly stabilizing mutations, one more stabilizing than any previously known mutation. Thus, physicochemical properties such as stability are latent within these large-scale protein functional data and can be revealed by systematic analysis. This approach should allow other protein properties to be discovered.

Redesign of a WW domain peptide for selective recognition of single-stranded DNA

A β-sheet miniprotein based on the FBP11 WW1 domain sequence has been redesigned for the molecular recognition of ssDNA. A previous report showed that a β-hairpin peptide dimer, (WKWK)(2), binds ssDNA with low micromolar affinity but with little selectivity over duplex DNA. This report extends those studies to a three-stranded β-sheet miniprotein designed to mimic the OB-fold. The new peptide binds ssDNA with low micromolar affinity and shows about 10-fold selectivity for ssDNA over duplex DNA. The redesigned peptide no longer binds its native ligand, the polyproline helix, confirming that the peptide has been redesigned for the function of binding ssDNA. Structural studies provide evidence that this peptide consists of a well-structured β-hairpin made of strands 2 and 3 with a less structured first strand that provides affinity for ssDNA but does not improve the stability of the full peptide. These studies provide insight into protein-DNA interactions as well as a novel example of protein redesign.

High-resolution mapping of protein sequence-function relationships

We present a large-scale approach to investigate the functional consequences of sequence variation in a protein. The approach entails the display of hundreds of thousands of protein variants, moderate selection for activity and high-throughput DNA sequencing to quantify the performance of each variant. Using this strategy, we tracked the performance of >600,000 variants of a human WW domain after three and six rounds of selection by phage display for binding to its peptide ligand. Binding properties of these variants defined a high-resolution map of mutational preference across the WW domain; each position had unique features that could not be captured by a few representative mutations. Our approach could be applied to many in vitro or in vivo protein assays, providing a general means for understanding how protein function relates to sequence.

Sequence determinants of thermodynamic stability in a WW domain--an all-beta-sheet protein

The stabilities of 66 sequence variants of the human Pin1 WW domain have been determined by equilibrium thermal denaturation experiments. All 34 residues composing the hPin1 WW three-stranded beta-sheet structure could be replaced one at a time with at least one different natural or non-natural amino acid residue without leading to an unfolded protein. Alanine substitutions at only four positions within the hPin1 WW domain lead to a partially or completely unfolded protein-in the absence of a physiological ligand. The side chains of these four residues form a conserved, partially solvent-inaccessible, continuous hydrophobic minicore comprising the N- and C-termini. Ala mutations at five other residues, three of which constitute the ligand binding patch on the concave side of the beta-sheet, significantly destabilize the hPin1 WW domain without leading to an unfolded protein. The remaining mutations affect protein stability only slightly, suggesting that only a small subset of side chain interactions within the hPin1 WW domain are mandatory for acquiring and maintaining a stable, cooperatively folded beta-sheet structure.

Evaluating beta-turn mimics as beta-sheet folding nucleators

Beta-turns are common conformations that enable proteins to adopt globular structures, and their formation is often rate limiting for folding. Beta-turn mimics, molecules that replace the i + 1 and i + 2 amino acid residues of a beta-turn, are envisioned to act as folding nucleators by preorganizing the pendant polypeptide chains, thereby lowering the activation barrier for beta-sheet formation. However, the crucial kinetic experiments to demonstrate that beta-turn mimics can act as strong nucleators in the context of a cooperatively folding protein have not been reported. We have incorporated 6 beta-turn mimics simulating varied beta-turn types in place of 2 residues in an engineered beta-turn 1 or beta-bulge turn 1 of the Pin 1 WW domain, a three-stranded beta-sheet protein. We present 2 lines of kinetic evidence that the inclusion of beta-turn mimics alters beta-sheet folding rates, enabling us to classify beta-turn mimics into 3 categories: those that are weak nucleators but permit Pin WW folding, native-like nucleators, and strong nucleators. Strong nucleators accelerate folding relative to WW domains incorporating all alpha-amino acid sequences. A solution NMR structure reveals that the native Pin WW beta-sheet structure is retained upon incorporating a strong E-olefin nucleator. These beta-turn mimics can now be used to interrogate protein folding transition state structures and the 2 kinetic analyses presented can be used to assess the nucleation capacity of other beta-turn mimics.

Compensatory evolution of a WW domain variant lacking the strictly conserved Trp residue

Replacement of conserved amino acid residues during evolution of proteins can lead to divergence and the formation of new families with novel functions, but is often deleterious to both protein structure and function. Using the WW domain, we experimentally examined whether and to what degree second-site mutations can compensate for the reduction of function and loss of structure that accompany substitution of a strictly conserved amino acid residue. The W17F mutant of the WW domain, with substitution of the most strictly conserved Trp residue, is known to lack a specific three-dimensional structure and shows reduced binding affinity in comparison to the wild type. To obtain second-site revertants, we performed a selection experiment based on the proline-rich peptide (PY ligand) binding affinity using the W17F mutant as the initial sequence. After selection by ribosome display, we were able to select revertants that exhibited a maximum ninefold higher affinity to the PY ligand than the W17F mutant and showed an even better affinity than the wild type. In addition, we found that the functional restoration resulted in increased binding specificity in selected revertants, and the structures were more compact, with increased amounts of secondary structure, in comparison to the W17F mutant. Our results suggest that the defective structure and function of the proteins caused by mutations in highly conserved residues occurring through divergent evolution not only can be restored but can be further improved by compensatory mutations.

The role of the turn in beta-hairpin formation during WW domain folding

The folding of WW domains is rate limited by formation of a beta-hairpin comprising residues from strands 1 and 2. Residues in the turn of this hairpin have reported Phi-values for folding close to 1 and have been proposed to nucleate folding. High Phi-values do not necessarily imply that the energetics of formation are a driving force for initiating folding. We demonstrate by NMR studies and molecular dynamics simulations that the first turn of the hYAP, FBP28, and PIN1 WW domains is structurally dynamic and solvent exposed in the native and folding transition states. It is, therefore, unlikely that the formation of the beta-turn per se provides the energetic driving force for hairpin folding. It is more likely that the turn acts as an easily formed hinge that facilitates the formation of the hairpin; it is a nucleus as defined by the nucleation-condensation mechanism whereby a diffuse nucleus is stabilized by associated interactions.

The creation of a novel fluorescent protein by guided consensus engineering

Consensus engineering has been used to increase the stability of a number of different proteins, either by creating consensus proteins from scratch or by modifying existing proteins so that their sequences more closely match a consensus sequence. In this paper we describe the first application of consensus engineering to the ab initio creation of a novel fluorescent protein. This was based on the alignment of 31 fluorescent proteins with >62% homology to monomeric Azami green (mAG) protein, and used the sequence of mAG to guide amino acid selection at positions of ambiguity. This consensus green protein is extremely well expressed, monomeric and fluorescent with red shifted absorption and emission characteristics compared to mAG. Although slightly less stable than mAG, it is better expressed and brighter under the excitation conditions typically used in single molecule fluorescence spectroscopy or confocal microscopy. This study illustrates the power of consensus engineering to create stable proteins using the subtle information embedded in the alignment of similar proteins and shows that the benefits of this approach may extend beyond stability.

Design and NMR conformational study of a beta-sheet peptide based on Betanova and WW domains

A good approach to test our current knowledge on formation of protein beta-sheets is de novo protein design. To obtain a three-stranded beta-sheet mini-protein, we have built a series of chimeric peptides by taking as a template a previously designed beta-sheet peptide, Betanova-LLM, and incorporating N- and/or C-terminal extensions taken from WW domains, the smallest natural beta-sheet domain that is stable in absence of disulfide bridges. Some Betanova-LLM strand residues were also substituted by those of a prototype WW domain. The designed peptides were cloned and expressed in Escherichia coli. The ability of the purified peptides to adopt beta-sheet structures was examined by circular dichroism (CD). Then, the peptide showing the highest beta-sheet population according to the CD spectra, named 3SBWW-2, was further investigated by 1H and 13C NMR. Based on NOE and chemical shift data, peptide 3SBWW-2 adopts a well defined three-stranded antiparallel beta-sheet structure with a disordered C-terminal tail. To discern between the contributions to beta-sheet stability of strand residues and the C-terminal extension, the structural behavior of a control peptide with the same strand residues as 3SBWW-2 but lacking the C-terminal extension, named Betanova-LYYL, was also investigated. beta-Sheet stability in these two peptides, in the parent Betanova-LLM and in WW-P, a prototype WW domain, decreased in the order WW-P > 3SBWW-2 > Betanova-LYYL > Betanova-LLM. Conclusions about the contributions to beta-sheet stability were drawn by comparing structural properties of these four peptides.

Φ-Analysis at the Experimental Limits: Mechanism of β-Hairpin Formation

The 37-residue Formin-binding protein, FBP28, is a canonical three-stranded beta-sheet WW domain. Because of its small size, it is so insensitive to chemical denaturation that it is barely possible to determine accurately a denaturation curve, as the transition spans 0-7 M guanidinium hydrochloride (GdmCl). It is also only marginally stable, with a free energy of denaturation of just 2.3 kcal/mol at 10 degrees Celsius so only small changes in energy upon mutation can be tolerated. But these properties and relaxation times for folding of 25 micros-400 micros conspire to allow the rapid acquisition of accurate and reproducible kinetic data for Phi-analysis using classical temperature-jump methods. The transition state for folding is highly polarized with some regions having Phi-values of 0 and others 1, as readily seen in chevron plots, with Phi-values of 0 having the refolding arms overlaying and those of 1 the unfolding arms superimposable. Good agreement is seen with transition state structures identified from independent molecular dynamics (MD) simulations at 60, 75, and 100 degrees Celsius, which allows us to explore further the details of the folding and unfolding pathway of FBP28. The first beta-turn is near native-like in the transition state for folding (experimental) and unfolding (MD and experiment). The simulations show that there are transient contacts between the aromatic side-chains of the beta-strands in the denatured state and that these interactions provide the driving force for folding of the first beta-hairpin of this three-stranded sheet. Only after the backbone hydrogen bonds are formed between beta1 and beta2 does a hydrogen bond form to stabilize the intervening turn, or the first beta-turn.

Relationship between local structural entropy and protein thermostability

We developed a technique to compute structural entropy directly from protein sequences. We explored the possibility of using structural entropy to identify residues involved in thermal stabilization of various protein families. Examples include methanococcal adenylate kinase, Ribonuclease HI and holocytochrome c(551). Our results show that the positions of the largest structural entropy differences between wild type and mutant usually coincide with the residues relevant to thermostability. We also observed a good linear relationship between the average structural entropy and the melting temperatures for adenylate kinase and its chimeric constructs. To validate this linear relationship, we compiled a large dataset comprised of 1153 sequences and found that most protein families still display similar linear relationships. Our results suggest that the multitude of interactions involved in thermal stabilization may be generalized into the tendency of proteins to maintain local structural conservation. The linear relationship between structural entropy and protein thermostability should be useful in the study of protein thermal stabilization.

Probing the Ligand-Binding Specificity and Analyzing the Folding State of SPOT-Synthesized FBP28 WW Domain Variants

The WW domains are known as the smallest naturally occurring, monomeric, triple-stranded, antiparallel beta-sheet domains. Hence, we chose the FBP28 WW domain as a model to investigate the stability of the beta-sheet structure at the amino acid level in the context of its function (ligand binding). The structure-function relationship was investigated through a complete substitution analysis of the FBP28 WW domain, with variants synthesized as a cellulose-bound peptide array. The functionality of the FBP28 WW domain variants was examined by probing the peptide array for ligand binding. In addition, selected FBP28 WW domain variants were investigated by CD measurements to determine the stability of the antiparallel beta-sheet structure. We discuss the correlation between structure stability and functionality for the FBP28 WW domain, as well as the effect of ligand-induced structure stabilization.

Conservation and amyloid formation: a study of the gelsolin-like family

The mechanism through which globular proteins transform into amyloid fibrils is still not understood. Here we analyze the structure and sequence conservation to assess the differential stability of segments from two structurally related protein families: the amyloidogenic gelsolin-like and its structurally related cofilin-like. The two families belong to the actin depolymerizing proteins, with a central beta-sheet stacked between 2 and 4 alpha-helices. Although sequentially remote, the two families share regions of high and low conservation and stability. Our results show a highly conserved hydrophobic and aromatic cluster, located at a central buried beta-hairpin. The geometry of the aromatic residues with respect to each other is strictly conserved, suggesting involvement in strand registering and beta-sheet stabilization. Consistent with experiment, we find a region of weak conservation and stability at one of the exposed beta-strands (strand B in the gelsolin-like family). This region was recently found to be affected by a point mutation-mediated destabilization of the human gelsolin domain 2, which facilitates the first proteolytic event in the formation of the amyloidogenic fragment. Thus, both experimental and computational conservation analyses suggest that this unstable region may constitute a first step in amyloid formation. Our analysis uses a recently developed multiple-structure comparison algorithm in which molecules are aligned simultaneously.

Design and characterization of a hyperstable p16INK4a that restores Cdk4 binding activity when combined with oncogenic mutations

Cyclin-dependent kinase inhibitor p16(INK4a) is the founding member of the INK4 family of tumor suppressors capable of arresting mammalian cell division. Missense mutations in the p16(INK4a) gene (INK4a/CDKN2A/MTS1) are strongly linked to several types of human cancer. These mutations are evenly distributed throughout this small, ankyrin repeat protein and the majority of them disrupt the native secondary and/or tertiary structure, leading to protein unfolding, aggregation and loss of function. We report here the use of multiple stabilizing substitutions to increase the stability of p16(INK4a) and furthermore, to restore Cdk4 binding activity of several defective, cancer-related mutant proteins. Stabilizing substitutions were predicted using four different techniques. The three most effective substitutions were combined to create a hyperstable p16(INK4a) variant that is 1.4 kcal/mol more stable than wild-type. This engineered construct is monomeric in solution with wild-type-like secondary and tertiary structure and cyclin-dependent kinase 4 binding activity. Interestingly, these hyperstable substitutions, when combined with oncogenic mutations R24P, P81L or V126D, can significantly restore Cdk4 binding activity, despite the divergent features of each destabilizing mutation. Extensive biophysical studies indicate that the hyperstable substitutions enhance the binding activity of mutant p16 through several different mechanisms, including an increased amount of secondary structure and thermostability, reduction in exposed hydrophobic surface(s) and/or a reduced tendency to aggregate. This apparent global suppressor effect suggests that increasing the thermodynamic stability of p16 can be used as a general strategy to restore the biological activity to defective mutants of this important tumor suppressor protein.

Definition of an inhibitory juxtamembrane WW-like domain in the platelet-derived growth factor beta receptor

A variety of tumors contain activating mutations in the cytoplasmic juxtamembrane domain of the type III family of receptor-tyrosine kinases, and some constructed mutations in this domain induce ligand-independent receptor activation. To explore the role of this domain in regulation of receptor activity, we subjected the juxtamembrane domain of the murine platelet-derived growth factor (PDGF) beta receptor to alanine-scanning mutagenesis. The mutant receptors were expressed in Ba/F3 cells and tested for constitutive tyrosine phosphorylation, association with phosphatidylinositol 3'-kinase, and their ability to induce cell survival and proliferation in the absence of interleukin-3. The mutant receptors accumulated to similar levels and appeared to undergo a normal PDGF-induced increase in tyrosine phosphorylation. Alanine substitutions at numerous positions located throughout the juxtamembrane domain caused constitutive receptor activation, as did an alanine insertion in the membrane-proximal segment of the juxtamembrane domain and a six-amino acid deletion in the center of the domain. It is possible to model the PDGF receptor juxtamembrane domain as a short alpha-helix followed by a three-stranded beta-sheet very similar to the known structures of WW domains. Strikingly, the activating mutations clustered in the central portions of the first and second beta strands and along one face of the beta-sheet, whereas the loops connecting the strands were largely devoid of mutationally sensitive positions. These findings provide strong support for the model that the activating mutations in the juxtamembrane region stimulate receptor activity by disrupting an inhibitory WW-like domain.

The carboxyl segment of the mumps virus V protein associates with Stat proteins in vitro via a tryptophan-rich motif

Viruses of the Paramyxovirinae, similar to other viruses, have evolved specific proteins that interdict IFN action as part of a general strategy to counteract host innate immunity. In many (but not all) cases, this interdiction is accompanied by a lowering of the intracellular levels of the STAT proteins. Among rubulaviruses, there is a notable variation in how they interfere with IFN action. Whereas SV41, SV5, and MuV all act by lowering Stat1, hPIV2 acts by lowering Stat2. Here, we show that the mumps and hPIV2 V proteins both form a complex with several Stat proteins in a mixed-extract assay. This suggests that the specific degradation of these Stat proteins is not determined by complex formation, but presumably at some later stage of the degradation pathway. V/Stat complex formation requires a specific carboxyl segment of V. However, a previously unrecognized trp-rich motif, rather than the Zn(++)-binding cys-cluster of this segment, appears to be required for V/Stat interaction. The C protein of Sendai (respiro-) virus, another P gene encoded protein, also forms a complex with Stat1, and prebinding of MuV V to Stat1 prevents the subsequent binding of SeV C. Our results suggest that rubulavirus V proteins may be related to both the C and the V proteins of respiroviruses.

BetaCore, a designed water soluble four-stranded antiparallel beta-sheet protein

BetaCore is a designed approximately 50-residue protein in which two BPTI-derived core modules, CM I and CM II, are connected by a 22-atom cross-link. At low temperature and pH 3, homo- and heteronuclear NMR data report a dominant folded ('f') conformation with well-dispersed chemical shifts, i, i+1 periodicity, numerous long-range NOEs, and slowed amide hydrogen isotope exchange patterns that is a four-stranded antiparallel beta-sheet with nonsymmetrical and specific association of CM I and CM II. BetaCore 'f' conformations undergo reversible, global, moderately cooperative, non-two-state thermal transitions to an equilibrium ensemble of unfolded 'u' conformations. There is a significant energy barrier between 'f' and 'u' conformations. This is the first designed four-stranded antiparallel beta-sheet that folds in water.

The consensus concept for thermostability engineering of proteins: further proof of concept

Previously, we calculated a consensus amino acid sequence from 13 homologous fungal phytases. A synthetic gene was constructed and recombinantly expressed. Surprisingly, consensus phytase-1 was 15-26 degrees C more thermostable than all parent phytases used in its design [Lehmann et al. (2000)Protein Eng., 13, 49-57]. In the present study, inclusion of six further phytase sequences in the amino acid sequence alignment resulted in the replacement of 38 amino acid residues in either one or both of the new consensus phytases-10 and -11. Since consensus phytase-10, again, was 7.4 degrees C more thermostable than consensus phytase-1, the thermostability effects of most of the 38 amino acid substitutions were tested by site-directed mutagenesis. Both stabilizing and destabilizing mutations were identified, but all affected the stability of the enzyme by <3 degrees C. The combination of all stabilizing amino acid exchanges in a multiple mutant of consensus phytase-1 increased the unfolding temperature from 78.0 to 88.5 degrees C. Likewise, back-mutation of four destabilizing amino acids and introduction of an additional stabilizing amino acid in consensus phytase-10 further increased the unfolding temperature from 85.4 to 90.4 degrees C. The thermostabilization achieved is the result of a combination of slight improvements from multiple amino acid exchanges rather than being the effect of a single or of just a few dominating mutations that have been introduced by chance. The present findings support the general validity of the consensus concept for thermostability engineering of proteins.

WW and SH3 domains, two different scaffolds to recognize proline-rich ligands

WW domains are small protein modules composed of approximately 40 amino acids. These domains fold as a stable, triple stranded beta-sheet and recognize proline-containing ligands. WW domains are found in many different signaling and structural proteins, often localized in the cytoplasm as well as in the cell nucleus. Based on analyses of seven structures of WW domains, we discuss their diverse binding preferences and sequence conservation patterns. While modeling WW domains for which structures have not been determined we uncovered a case of potential molecular and functional convergence between WW and SH3 domains. The binding surface of the modeled WW domain of Npw38 protein shows a remarkable similarity to the SH3 domain of Sem5 protein, confirming biochemical data on similar binding predilections of both domains.

NMR solution structure of the isolated Apo Pin1 WW domain: comparison to the x-ray crystal structures of Pin1

The NMR solution structure of the isolated Apo Pin1 WW domain (6-39) reveals that it adopts a twisted three-stranded antiparallel beta-sheet conformation, very similar to the structure exhibited by the crystal of this domain in the context of the two domain Pin1 protein. While the B factors in the apo x-ray crystal structure indicate that loop 1 and loop 2 are conformationally well defined, the solution NMR data suggest that loop 1 is quite flexible, at least in the absence of the ligand. The NMR chemical shift and nuclear Overhauser effect pattern exhibited by the 6-39 Pin1 WW domain has proven to be diagnostic for demonstrating that single site variants of this domain adopt a normally folded structure. Knowledge of this type is critical before embarking on time-consuming kinetic and thermodynamic studies required for a detailed understanding of beta-sheet folding.

Structure and stability effects of mutations designed to increase the primary sequence symmetry within the core region of a beta-trefoil

Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil hyperfamily and exhibits a characteristic threefold symmetry of the tertiary structure. However, evidence of this symmetry is not readily apparent at the level of the primary sequence. This suggests that while selective pressures may exist to retain (or converge upon) a symmetric tertiary structure, other selective pressures have resulted in divergence of the primary sequence during evolution. Using intra-chain and homologue sequence comparisons for 19 members of this family of proteins, we have designed mutants of FGF-1 that constrain a subset of core-packing residues to threefold symmetry at the level of the primary sequence. The consequences of these mutations regarding structure and stability were evaluated using a combination of X-ray crystallography and differential scanning calorimetry. The mutational effects on structure and stability can be rationalized through the characterization of "microcavities" within the core detected using a 1.0A probe radius. The results show that the symmetric constraint within the primary sequence is compatible with a well-packed core and near wild-type stability. However, despite the general maintenance of overall thermal stability, a noticeable increase in non-two-state denaturation follows the increase in primary sequence symmetry. Therefore, properties of folding, rather than stability, may contribute to the selective pressure for asymmetric primary core sequences within symmetric protein architectures.