Two proteins with the same structure respond very differently to mutation: the role of plasticity in protein stability

Design for Solubility May Reveal Induction of Amide Hydrogen/Deuterium Exchange by Protein Self-Association

Structural heterogeneity often constrains the characterization of aggregating proteins to indirect or low-resolution methods, obscuring mechanistic details of association. Here, we report progress in understanding the aggregation of Adnectins, engineered binding proteins with an immunoglobulin-like fold. We rationally design Adnectin solubility and measure amide hydrogen/deuterium exchange (HDX) under conditions that permit transient protein self-association. Protein-protein binding commonly slows rates of HDX; in contrast, we find that Adnectin association may induce faster HDX for certain amides, particularly in the C-terminal β-strand. In aggregation-prone proteins, we identify a pattern of very different rates of amide HDX for residues linked by reciprocal hydrogen bonds in the native structure. These results may be explained by local loss of native structure and formation of an inter-protein interface. Amide HDX induced by self-association, detected here by deliberate modulation of propensity for such interactions, may be a general phenomenon with the potential to expose mechanisms of aggregation by diverse proteins.

A study on non-synonymous mutational patterns in structural proteins of SARS-CoV-2

SARS-CoV-2 is mutating and creating divergent variants across the world. An in-depth investigation of the amino acid substitutions in the genomic signature of SARS-CoV-2 proteins is highly essential for understanding its host adaptation and infection biology. A total of 9587 SARS-CoV-2 structural protein sequences collected from 49 different countries are used to characterize protein-wise variants, substitution patterns (type and location), and major substitution changes. The majority of the substitutions are distinct, mostly in a particular location, and lead to a change in an amino acid's biochemical properties. In terms of mutational changes, envelope (E) and membrane (M) proteins are relatively more stable than nucleocapsid (N) and spike (S) proteins. Several co-occurrence substitutions are observed, particularly in S and N proteins. Substitution specific to active sub-domains reveals that heptapeptide repeat, fusion peptides, transmembrane in S protein, and N-terminal and C-terminal domains in the N protein are remarkably mutated. We also observe a few deleterious mutations in the above domains. The overall study on non-synonymous mutation in structural proteins of SARS-CoV-2 at the start of the pandemic indicates a diversity amongst virus sequences.

Spectroscopic Evidence of the Salt-Induced Conformational Change around the Localized Electric Charges on the Protein Surface of Fibronectin Type III

The effect of salt on the electrostatic interaction of a protein is an important issue, because addition of salt affects protein stability and association/aggregation. Although adding salt is a generally recognized strategy to improve protein stability, this improvement does not necessarily occur. The lack of an effect upon the addition of salt was previously confirmed for the tenth fibronectin type III domain from human fibronectin (FN3) by thermal stability analysis. However, the detailed molecular mechanism is unknown. In the present study, by employing the negatively charged carboxyl triad on the surface of FN3 as a case study, the molecular mechanism of the inefficient NaCl effect on protein stability was experimentally addressed using spectroscopic methods. Complementary analysis using Raman spectroscopy and 8-anilino-1-naphthalenesulfonic acid fluorescence revealed the three-phase behavior of the salt-protein interaction between NaCl and FN3 over a wide salt concentration range from 100 mM to 4.0 M, suggesting that the Na⁺-specific binding to the negatively charged carboxyl triad causes a local conformational change around the binding site with an accompanying structural change in the overall protein, which contributes to the protein's structural destabilization. This spectroscopic evidence clarifies the molecular understanding of the inefficiency of salt to improve protein stability. The findings will inform the optimization of formulation conditions.

Toxin Neutralization Using Alternative Binding Proteins

Animal toxins present a major threat to human health worldwide, predominantly through snakebite envenomings, which are responsible for over 100,000 deaths each year. To date, the only available treatment against snakebite envenoming is plasma-derived antivenom. However, despite being key to limiting morbidity and mortality among snakebite victims, current antivenoms suffer from several drawbacks, such as immunogenicity and high cost of production. Consequently, avenues for improving envenoming therapy, such as the discovery of toxin-sequestering monoclonal antibodies against medically important target toxins through phage display selection, are being explored. However, alternative binding protein scaffolds that exhibit certain advantages compared to the well-known immunoglobulin G scaffold, including high stability under harsh conditions and low cost of production, may pose as possible low-cost alternatives to antibody-based therapeutics. There is now a plethora of alternative binding protein scaffolds, ranging from antibody derivatives (e.g., nanobodies), through rationally designed derivatives of other human proteins (e.g., DARPins), to derivatives of non-human proteins (e.g., affibodies), all exhibiting different biochemical and pharmacokinetic profiles. Undeniably, the high level of engineerability and potentially low cost of production, associated with many alternative protein scaffolds, present an exciting possibility for the future of snakebite therapeutics and merit thorough investigation. In this review, a comprehensive overview of the different types of binding protein scaffolds is provided together with a discussion on their relevance as potential modalities for use as next-generation antivenoms.

Protein folding, misfolding and aggregation: The importance of two-electron stabilizing interactions

Proteins associated with neurodegenerative diseases are highly pleiomorphic and may adopt an all-α-helical fold in one environment, assemble into all-β-sheet or collapse into a coil in another, and rapidly polymerize in yet another one via divergent aggregation pathways that yield broad diversity of aggregates’ morphology. A thorough understanding of this behaviour may be necessary to develop a treatment for Alzheimer’s and related disorders. Unfortunately, our present comprehension of folding and misfolding is limited for want of a physicochemical theory of protein secondary and tertiary structure. Here we demonstrate that electronic configuration and hyperconjugation of the peptide amide bonds ought to be taken into account to advance such a theory. To capture the effect of polarization of peptide linkages on conformational and H-bonding propensity of the polypeptide backbone, we introduce a function of shielding tensors of the C^α atoms. Carrying no information about side chain-side chain interactions, this function nonetheless identifies basic features of the secondary and tertiary structure, establishes sequence correlates of the metamorphic and pH-driven equilibria, relates binding affinities and folding rate constants to secondary structure preferences, and manifests common patterns of backbone density distribution in amyloidogenic regions of Alzheimer’s amyloid β and tau, Parkinson’s α-synuclein and prions. Based on those findings, a split-intein like mechanism of molecular recognition is proposed to underlie dimerization of Aβ, tau, αS and PrP^C, and divergent pathways for subsequent association of dimers are outlined; a related mechanism is proposed to underlie formation of PrP^Sc fibrils. The model does account for: (i) structural features of paranuclei, off-pathway oligomers, non-fibrillar aggregates and fibrils; (ii) effects of incubation conditions, point mutations, isoform lengths, small-molecule assembly modulators and chirality of solid-liquid interface on the rate and morphology of aggregation; (iii) fibril-surface catalysis of secondary nucleation; and (iv) self-propagation of infectious strains of mammalian prions.

Predicting the impact of Lynch syndrome-causing missense mutations from structural calculations

Accurate methods to assess the pathogenicity of mutations are needed to fully leverage the possibilities of genome sequencing in diagnosis. Current data-driven and bioinformatics approaches are, however, limited by the large number of new variations found in each newly sequenced genome, and often do not provide direct mechanistic insight. Here we demonstrate, for the first time, that saturation mutagenesis, biophysical modeling and co-variation analysis, performed in silico, can predict the abundance, metabolic stability, and function of proteins inside living cells. As a model system, we selected the human mismatch repair protein, MSH2, where missense variants are known to cause the hereditary cancer predisposition disease, known as Lynch syndrome. We show that the majority of disease-causing MSH2 mutations give rise to folding defects and proteasome-dependent degradation rather than inherent loss of function, and accordingly our in silico modeling data accurately identifies disease-causing mutations and outperforms the traditionally used genetic disease predictors. Thus, in conclusion, in silico biophysical modeling should be considered for making genotype-phenotype predictions and for diagnosis of Lynch syndrome, and perhaps other hereditary diseases.

The protein quality control system targets misfolded proteins for degradation. So far it has not been possible from sequence or structural data to predict the biological stability of a misfolded protein, or the effect of mutations on intracellular protein levels. Here we demonstrate that in silico saturation mutagenesis and biophysical calculations of the structural stability of the human mismatch repair protein MSH2 correlate with cellular protein levels, turnover and function. Of 24 different MSH2 variants, some of which are linked to Lynch syndrome, a destabilization of as little as 3 kcal/mol is sufficient to cause rapid degradation via the ubiquitin-proteasome pathway. Thus, biophysical modeling can, to a large extent, predict the metabolic stability of proteins. We also show that the same biophysical calculations can be used to distinguish with high accuracy neutral sequence variation from pathogenic variants, and that the calculations outperform several traditionally used disease predictors. We therefore suggest the method to be of potential value for patient stratification in Lynch syndrome, and perhaps other hereditary diseases.

Spontaneous Unfolding-Refolding of Fibronectin Type III Domains Assayed by Thiol Exchange: THERMODYNAMIC STABILITY CORRELATES WITH RATES OF UNFOLDING RATHER THAN FOLDING

Globular proteins are not permanently folded but spontaneously unfold and refold on time scales that can span orders of magnitude for different proteins. A longstanding debate in the protein-folding field is whether unfolding rates or folding rates correlate to the stability of a protein. In the present study, we have determined the unfolding and folding kinetics of 10 FNIII domains. FNIII domains are one of the most common protein folds and are present in 2% of animal proteins. FNIII domains are ideal for this study because they have an identical seven-strand β-sandwich structure, but they vary widely in sequence and thermodynamic stability. We assayed thermodynamic stability of each domain by equilibrium denaturation in urea. We then assayed the kinetics of domain opening and closing by a technique known as thiol exchange. For this we introduced a buried Cys at the identical location in each FNIII domain and measured the kinetics of labeling with DTNB over a range of urea concentrations. A global fit of the kinetics data gave the kinetics of spontaneous unfolding and refolding in zero urea. We found that the folding rates were relatively similar, ∼0.1-1 s^-1, for the different domains. The unfolding rates varied widely and correlated with thermodynamic stability. Our study is the first to address this question using a set of domains that are structurally homologous but evolved with widely varying sequence identity and thermodynamic stability. These data add new evidence that thermodynamic stability correlates primarily with unfolding rate rather than folding rate. The study also has implications for the question of whether opening of FNIII domains contributes to the stretching of fibronectin matrix fibrils.

A many-body term improves the accuracy of effective potentials based on protein coevolutionary data

The study of correlated mutations in alignments of homologous proteins proved to be successful not only in the prediction of their native conformation but also in the development of a two-body effective potential between pairs of amino acids. In the present work, we extend the effective potential, introducing a many-body term based on the same theoretical framework, making use of a principle of maximum entropy. The extended potential performs better than the two-body one in predicting the energetic effect of 308 mutations in 14 proteins (including membrane proteins). The average value of the parameters of the many-body term correlates with the degree of hydrophobicity of the corresponding residues, suggesting that this term partly reflects the effect of the solvent.

Predicting protein thermal stability changes upon point mutations using statistical potentials: Introducing HoTMuSiC

The accurate prediction of the impact of an amino acid substitution on the thermal stability of a protein is a central issue in protein science, and is of key relevance for the rational optimization of various bioprocesses that use enzymes in unusual conditions. Here we present one of the first computational tools to predict the change in melting temperature ΔT_m upon point mutations, given the protein structure and, when available, the melting temperature T_m of the wild-type protein. The key ingredients of our model structure are standard and temperature-dependent statistical potentials, which are combined with the help of an artificial neural network. The model structure was chosen on the basis of a detailed thermodynamic analysis of the system. The parameters of the model were identified on a set of more than 1,600 mutations with experimentally measured ΔT_m. The performance of our method was tested using a strict 5-fold cross-validation procedure, and was found to be significantly superior to that of competing methods. We obtained a root mean square deviation between predicted and experimental ΔT_m values of 4.2 °C that reduces to 2.9 °C when ten percent outliers are removed. A webserver-based tool is freely available for non-commercial use at soft.dezyme.com.

Structural and dynamic properties that govern the stability of an engineered fibronectin type III domain

Consensus protein design is a rapid and reliable technique for the improvement of protein stability, which relies on the use of homologous protein sequences. To enhance the stability of a fibronectin type III (FN3) domain, consensus design was employed using an alignment of 2123 sequences. The resulting FN3 domain, FN3con, has unprecedented stability, with a melting temperature >100°C, a ΔG_D−N of 15.5 kcal mol⁻¹ and a greatly reduced unfolding rate compared with wild-type. To determine the underlying molecular basis for stability, an X-ray crystal structure of FN3con was determined to 2.0 Å and compared with other FN3 domains of varying stabilities. The structure of FN3con reveals significantly increased salt bridge interactions that are cooperatively networked, and a highly optimized hydrophobic core. Molecular dynamics simulations of FN3con and comparison structures show the cooperative power of electrostatic and hydrophobic networks in improving FN3con stability. Taken together, our data reveal that FN3con stability does not result from a single mechanism, but rather the combination of several features and the removal of non-conserved, unfavorable interactions. The large number of sequences employed in this study has most likely enhanced the robustness of the consensus design, which is now possible due to the increased sequence availability in the post-genomic era. These studies increase our knowledge of the molecular mechanisms that govern stability and demonstrate the rising potential for enhancing stability via the consensus method.

High-Throughput Ligand Discovery Reveals a Sitewise Gradient of Diversity in Broadly Evolved Hydrophilic Fibronectin Domains

Discovering new binding function via a combinatorial library in small protein scaffolds requires balance between appropriate mutations to introduce favorable intermolecular interactions while maintaining intramolecular integrity. Sitewise constraints exist in a non-spatial gradient from diverse to conserved in evolved antibody repertoires; yet non-antibody scaffolds generally do not implement this strategy in combinatorial libraries. Despite the fact that biased amino acid distributions, typically elevated in tyrosine, serine, and glycine, have gained wider use in synthetic scaffolds, these distributions are still predominantly applied uniformly to diversified sites. While select sites in fibronectin domains and DARPins have shown benefit from sitewise designs, they have not been deeply evaluated. Inspired by this disparity between diversity distributions in natural libraries and synthetic scaffold libraries, we hypothesized that binders resulting from discovery and evolution would exhibit a non-spatial, sitewise gradient of amino acid diversity. To identify sitewise diversities consistent with efficient evolution in the context of a hydrophilic fibronectin domain, >105 binders to six targets were evolved and sequenced. Evolutionarily favorable amino acid distributions at 25 sites reveal Shannon entropies (range: 0.3-3.9; median: 2.1; standard deviation: 1.1) supporting the diversity gradient hypothesis. Sitewise constraints in evolved sequences are consistent with complementarity, stability, and consensus biases. Implementation of sitewise constrained diversity enables direct selection of nanomolar affinity binders validating an efficient strategy to balance inter- and intra-molecular interaction demands at each site.

Characterization of monobody scaffold interactions with ligand via force spectroscopy and steered molecular dynamics

Monobodies are antibody alternatives derived from fibronectin that are thermodynamically stable, small in size, and can be produced in bacterial systems. Monobodies have been engineered to bind a wide variety of target proteins with high affinity and specificity. Using alanine-scanning mutagenesis simulations, we identified two scaffold residues that are critical to the binding interaction between the monobody YS1 and its ligand, maltose-binding protein (MBP). Steered molecular dynamics (SMD) simulations predicted that the E47A and R33A mutations in the YS1 scaffold substantially destabilize the YS1-MBP interface by reducing the bond rupture force and the lifetime of single hydrogen bonds. SMD simulations further indicated that the R33A mutation weakens the hydrogen binding between all scaffold residues and MBP and not just between R33 and MBP. We validated the simulation data and characterized the effects of mutations on YS1-MBP binding by using single-molecule force spectroscopy and surface plasmon resonance. We propose that interfacial stability resulting from R33 of YS1 stacking with R344 of MBP synergistically stabilizes both its own bond and the interacting scaffold residues of YS1. Our integrated approach improves our understanding of the monobody scaffold interactions with a target, thus providing guidance for the improved engineering of monobodies.

Stabilization of the third fibronectin type III domain of human tenascin-C through minimal mutation and rational design

Non-antibody scaffolds are increasingly used to generate novel binding proteins for both research and therapeutic applications. Our group has developed the tenth fibronectin type III domain of human tenascin-C (TNfn3) as one such scaffold. As a scaffold, TNfn3 must tolerate extensive mutation to introduce novel binding sites. However, TNfn3's marginal stability (T(m) ∼ 59°C, ΔG(unfolding) = 5.7 kcal/mol) stands as a potential obstacle to this process. To address this issue, we sought to engineer highly stable TNfn3 variants. We used two parallel strategies. Using insights gained from structural analysis of other FN3 family members, we (1) rationally designed stabilizing point mutations or (2) introduced novel stabilizing disulfide bonds. Both strategies yielded highly stable TNfn3 variants with T(m) values as high as 83°C and ΔG(unfolding) values as high as 9.4 kcal/mol. Notably, only three or four mutations were required to achieve this level of stability with either approach. These results validate our rational design strategies and illustrate that substantial stability increases can be achieved with minimal mutation. One TNfn3 variant reported here has now been successfully used as a scaffold to develop two promising therapeutic molecules. We anticipate that other variants described will exhibit similar utility.

Energetic and structural details of the trigger-loop closing transition in RNA polymerase II

An evolutionarily conserved element in RNA polymerase II, the trigger loop (TL), has been suggested to play an important role in the elongation rate, fidelity of selection of the matched nucleoside triphosphate (NTP), catalysis of transcription elongation, and translocation in both eukaryotes and prokaryotes. In response to NTP binding, the TL undergoes large conformational changes to switch between distinct open and closed states to tighten the active site and avail catalysis. A computational strategy for characterizing the conformational transition pathway is presented to bridge the open and closed states of the TL. Information from a large number of independent all-atom molecular dynamics trajectories from Hamiltonian replica exchange and targeted molecular dynamics simulations is gathered together to assemble a connectivity map of the conformational transition. The results show that with a cognate NTP, TL closing should be a spontaneous process. One major intermediate state is identified along the conformational transition pathway, and the key structural features are characterized. The complete pathway from the open TL to the closed TL provides a clear picture of the TL closing.

SLLISWD sequence in the 10FNIII domain initiates fibronectin fibrillogenesis

Fibronectin (FN) assembly into extracellular matrix is tightly regulated and essential to embryogenesis and wound healing. FN fibrillogenesis is initiated by cytoskeleton-derived tensional forces transmitted across transmembrane integrins onto RGD binding sequences within the tenth FN type III (10FNIII) domains. These forces unfold 10FNIII to expose cryptic FN assembly sites; however, a specific sequence has not been identified in 10FNIII. Our past steered molecular dynamics simulations modeling 10FNIII unfolding by force at its RGD loop predicted a mechanical intermediate with a solvent-exposed N terminus spanning the A and B β-strands. Here, we experimentally confirm that the predicted 23-residue cryptic peptide 1 (CP1) initiates FN multimerization, which is mediated by interactions with 10FNIII that expose hydrophobic surfaces that support 8-anilino-1-napthalenesulfonic acid binding. Localization of multimerization activity to the C terminus led to the discovery of a minimal 7-amino acid "multimerization sequence" (SLLISWD), which induces polymerization of FN and the clotting protein fibrinogen in addition to enhancing FN fibrillogenesis in fibroblasts. A point mutation at Trp-6 that reduces exposure of hydrophobic sites for 8-anilino-1-napthalenesulfonic acid binding and β-structure formation inhibits FN multimerization and prevents physiological cell-based FN assembly in culture. We propose a model for cell-mediated fibrillogenesis whereby cell traction force initiates a cascade of intermolecular exchange starting with the unfolding of 10FNIII to expose the multimerization sequence, which interacts with strand B of another 10FNIII domain via a Trp-mediated β-strand exchange to stabilize a partially unfolded intermediate that propagates FN self-assembly.

Atomic-level description of ubiquitin folding

Equilibrium molecular dynamics simulations, in which proteins spontaneously and repeatedly fold and unfold, have recently been used to help elucidate the mechanistic principles that underlie the folding of fast-folding proteins. The extent to which the conclusions drawn from the analysis of such proteins, which fold on the microsecond timescale, apply to the millisecond or slower folding of naturally occurring proteins is, however, unclear. As a first attempt to address this outstanding issue, we examine here the folding of ubiquitin, a 76-residue-long protein found in all eukaryotes that is known experimentally to fold on a millisecond timescale. Ubiquitin folding has been the subject of many experimental studies, but its slow folding rate has made it difficult to observe and characterize the folding process through all-atom molecular dynamics simulations. Here we determine the mechanism, thermodynamics, and kinetics of ubiquitin folding through equilibrium atomistic simulations. The picture emerging from the simulations is in agreement with a view of ubiquitin folding suggested from previous experiments. Our findings related to the folding of ubiquitin are also consistent, for the most part, with the folding principles derived from the simulation of fast-folding proteins, suggesting that these principles may be applicable to a wider range of proteins.

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6±0.3 kcal/mol per -CH(2)- group), than to the stability of a large protein, VlsE (1.6±0.3 kcal/mol per -CH(2)- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH(2)- group on folding contributes, on average, 1.1±0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG(tr) values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60±4% and hydrogen bonds contribute 40±4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.

Adnectins: engineered target-binding protein therapeutics

Adnectins™ are a new family of therapeutic proteins based on the 10th fibronectin type III domain, and designed to bind with high affinity and specificity to therapeutically relevant targets. Adnectins share with antibody variable domains a beta-sheet sandwich fold with diversified loops, but differ from antibodies in primary sequence and have a simpler, single-domain structure without disulfide bonds. As a consequence, Adnectins bind targets with affinity and specificity as high as those of antibodies, but are easier to manipulate genetically and compatible with bacterial expression systems. Adnectins that bind macromolecular targets with nanomolar and picomolar affinity have been selected using in vitro evolution methods, including mRNA display, phage display and yeast display. CT-322, a PEGylated, anti-angiogenic Adnectin that binds vascular endothelial growth factor (VEGF) receptor 2 and blocks its interaction with VEGF A, C and D, is being evaluated in Phase II clinical trials for efficacy in several oncology indications.

Water's role in the force-induced unfolding of ubiquitin

In atomic force spectroscopic studies of the elastomeric protein ubiquitin, the β-strands 1-5 serve as the force clamp. Simulations show how the rupture force in the force-induced unfolding depends on the kinetics of water molecule insertion into positions where they can eventually form hydrogen bonding bridges with the backbone hydrogen bonds in the force-clamp region. The intrusion of water into this region is slowed down by the hydrophobic shielding effect of carbonaceous groups on the surface residues of β-strands 1-5, which thereby regulates water insertion prior to hydrogen bond breakage. The experiments show that the unfolding of the mechanically stressed protein is nonexponential due to static disorder. Our simulations show that different numbers and/or locations of bridging water molecules give rise to a long-lived distribution of transition states and static disorder. We find that slowing down the translational (not rotational) motions of the water molecules by increasing the mass of their oxygen atoms, which leaves the force field and thereby the equilibrium structure of the solvent unchanged, increases the average rupture force; however, the early stages of the force versus time behavior are very similar for our "normal" and fictitious "heavy" water models. Finally, we construct six mutant systems to regulate the hydrophobic shielding effect of the surface residues in the force-clamp region. The mutations in the two termini of β-sheets 1-5 are found to determine a preference for different unfolding pathways and change mutant's average rupture force.

Underexposed polar residues and protein stabilization

Increasing protein stability is interesting for practical reasons and because it tests our understanding of protein energetics. We explore here the feasibility of stabilizing proteins by replacing underexposed polar residues by apolar ones of similar size and shape. We have compared the stability of wild-type apoflavodoxin with that of a few carefully selected mutants carrying Y → F, Q → L, T → V or K → M replacements. Although a clear inverse correlation between native solvent exposures of replaced polar residues and stability of mutants is observed, most mutations fail to stabilize the protein. The promising exceptions are the two Q → L mutations tested, which characteristically combine the greatest reduction in polar burial with the greatest increase in apolar burial relative to wild type. Analysis of published stability data corresponding to a variety of mutant proteins confirms that, unlike Y → F or T → V replacements, Q → L mutations tend to be stabilizing, and it suggests that N → L mutations might be stabilizing as well. On the other hand, we show that the stability changes associated to the apoflavodoxin mutations can be rationalized in terms of differential polar and apolar burials upon folding plus a generic destabilizing penalty term. Simple equations combining these contributions predict stability changes in a large data set of 113 mutants (Y → F, Q → L or T → V) similarly well as more complex algorithms available on the Internet.

The HD-exchange motions of ribosomal protein S6 are insensitive to reversal of the protein-folding pathway

An increasing number of protein structures are found to encompass multiple folding nuclei, allowing their structures to be formed by several competing pathways. A typical example is the ribosomal protein S6, which comprises two folding nuclei (sigma1 and sigma2) defining two competing pathways in the folding energy landscape: sigma1 --> sigma2 and sigma2 --> sigma1. The balance between the two pathways, and thus the order of folding events, is easily controlled by circular permutation. In this study, we make use of this ability to manipulate the folding pathway to demonstrate that the dynamic motions of the S6 structure are independent of how the protein folds. The HD-exchange protection factors remain the same upon complete reversal of the folding order. The phenomenon arises because the HD-exchange motions and the high-energy excitations controlling the folding pathway occur at separated free-energy levels: the Boltzmann distribution of unproductive unfolding attempts samples all unfolding channels in parallel, even those that end up in excessively high barriers. Accordingly, the findings provide a simple rationale for how to interpret native-state dynamics without the need to invoke fluctuations off the normal unfolding reaction coordinate.

Transient opening of fibronectin type III (FNIII) domains: the interaction of the third FNIII domain of FN with anastellin

We previously reported that the fibronectin (FN) type III domains of FN may unfold to interact with anastellin and form FN aggregates. In the present study, we have focused on the interaction between anastellin and the third FN type III domain (III3), which is a key anastellin binding site on FN. Anastellin binding to III3 was monitored by 8-anilino-1-naphthalene sulfonate (ANS) fluorescence. ANS binding to anastellin dramatically increased its emission intensity, but this was reduced to half by the addition of III3, suggesting that ANS and III3 share a common hydrophobic binding site on anastellin. An engineered mutant of III3 that was stabilized by an intrachain disulfide bond did not interact with anastellin, as seen by its failure to interfere with ANS binding to anastellin. We also mutated hydrophobic core residues to destabilize III3 and found that these mutants were still capable of interacting with anastellin. Anastellin binding to III3 was also monitored using an intramolecular green fluorescent protein (GFP)-based fluorescence resonance energy transfer (FRET) construct, in which III3 was flanked by two GFP variants (III3-FRET). Anastellin bound to III3-FRET and caused an increase in the FRET signal. The dissociation constant was estimated to be approximately 210 nM. The binding kinetics of anastellin to III3-FRET fit a first-order reaction with a half-time of approximately 30 s; the kinetics with destabilized III3 mutants were even faster. Matrix-assisted laser desorption ionization-time-of-flight (MALDI-TOF) mass spectrometry suggested that the middle part of III3 became destabilized and protease sensitive upon anastellin binding. Thus, the stability of III3 seems to be a key factor in anastellin binding.

The Kir channel immunoglobulin domain is essential for Kir1.1 (ROMK) thermodynamic stability, trafficking and gating

The renal inward rectifying potassium channel Kir1.1 plays key roles in regulating electrolyte homeostasis and blood pressure. Loss-of-function mutations in the channel cause a life-threatening salt and water balance disorder in infants called antenatal Bartter syndrome (ABS). Of more than 30 ABS mutations identified, approximately half are located in the intracellular domain of the channel. The mechanisms underlying channel dysfunction for most of these mutations are unknown. By mapping intracellular mutations onto an atomic model of Kir1.1, we found that several of these are localized to a phylogenetically ancient immunoglobulin (Ig)-like domain (IgLD) that has not been characterized previously, prompting us to examine this structure in detail. The IgLD is assembled from two beta-pleated sheets packed face-to-face, creating a beta-sheet interface or core, populated by highly conserved side chains. Thermodynamic calculations on computationally mutated channels suggest that IgLD core residues are among the most important residues for determining cytoplasmic domain stability. Consistent with this notion, we show that two ABS mutations (A198T and Y314C) located within the IgLD core impair channel biosynthesis and trafficking in mammalian cells. A fraction of core mutant channels reach the cell surface, but are electrically silent due to closure of the helix-bundle gate. Compensatory mutation-induced rescue of channel function revealed that IgLD core mutants fail to rectify. Our study sheds new light on the pathogenesis of ABS and establishes the IgLD as an essential structure within the Kir channel family.

High-throughput analysis of the protein sequence-stability landscape using a quantitative yeast surface two-hybrid system and fragment reconstitution

Stability evaluation of many mutants can lead to a better understanding of the sequence determinants of a structural motif and of factors governing protein stability and protein evolution. The traditional biophysical analysis of protein stability is low throughput, limiting our ability to widely explore sequence space in a quantitative manner. In this study, we have developed a high-throughput library screening method for quantifying stability changes, which is based on protein fragment reconstitution and yeast surface display. Our method exploits the thermodynamic linkage between protein stability and fragment reconstitution and the ability of the yeast surface display technique to quantitatively evaluate protein-protein interactions. The method was applied to a fibronectin type III (FN3) domain. Characterization of fragment reconstitution was facilitated by the co-expression of two FN3 fragments, thus establishing a yeast surface two-hybrid method. Importantly, our method does not rely on competition between clones and thus eliminates a common limitation of high-throughput selection methods in which the most stable variants are recovered predominantly. Thus, it allows for the isolation of sequences that exhibit a desired level of stability. We identified more than 100 unique sequences for a beta-bulge motif, which was significantly more informative than natural sequences of the FN3 family in revealing the sequence determinants for the beta-bulge. Our method provides a powerful means for the rapid assessment of the stability of many variants, for the systematic assessment of the contribution of different factors to protein stability, and for enhancement of the protein stability.

Role of a tryptophan anchor in human topoisomerase I structure, function and inhibition

Human Top1 (topoisomerase I) relaxes supercoiled DNA during cell division and transcription. Top1 is composed of 765 amino acids and contains an unstructured N-terminal domain of 200 amino acids, and a structured functional domain of 565 amino acids that binds and relaxes supercoiled DNA. In the present study we examined the region spanning the junction of the N-terminal domain and functional domain (junction region). Analysis of several published Top1 structures revealed that three tryptophan residues formed a network of aromatic stacking interactions and electrostatic interactions that anchored the N-terminus of the functional domain to sub-domains containing the nose cone and active site. Mutation of the three tryptophan residues (Trp(203)/Trp(205)/Trp(206)) to an alanine residue, either individually or together, in silico revealed that the individual tryptophan residue's contribution to the tryptophan 'anchor' was additive. When the three tryptophan residues were mutated to alanine in vitro, the resulting mutant Top1 differed from wild-type Top1 in that it lacked processivity, exhibited resistance to camptothecin and was inactivated by urea. The results indicated that the tryptophan anchor stabilized the N-terminus of the functional domain and prevented the loss of Top1 structure and function.

Crosstalk between the protein surface and hydrophobic core in a core-swapped fibronectin type III domain

Two homologous fibronectin type III (fnIII) domains, FNfn10 (the 10th fnIII domain of human fibronectin) and TNfn3 (the third fnIII domain of human tenascin), have essentially the same backbone structure, although they share only ∼ 24% sequence identity. While they share a similar folding mechanism with a common core of key residues in the folding transition state, they differ in many other physical properties. We use a chimeric protein, FNoTNc, to investigate the molecular basis for these differences. FNoTNc is a core-swapped protein, containing the “outside” (surface and loops) of FNfn10 and the hydrophobic core of TNfn3. Remarkably, FNoTNc retains the structure of the parent proteins despite the extent of redesign, allowing us to gain insight into which components of each parent protein are responsible for different aspects of its behaviour. Naively, one would expect properties that appear to depend principally on the core to be similar to TNfn3, for example, the response to mutations, folding kinetics and side-chain dynamics, while properties apparently determined by differences in the surface and loops, such as backbone dynamics, would be more like FNfn10. While this is broadly true, it is clear that there are also unexpected crosstalk effects between the core and the surface. For example, the anomalous response of FNfn10 to mutation is not solely a property of the core as we had previously suggested.

An experimental study of GFP-based FRET, with application to intrinsically unstructured proteins

We have experimentally studied the fluorescence resonance energy transfer (FRET) between green fluorescent protein (GFP) molecules by inserting folded or intrinsically unstructured proteins between CyPet and Ypet. We discovered that most of the enhanced FRET signal previously reported for this pair was due to enhanced dimerization, so we engineered a monomerizing mutation into each. An insert containing a single fibronectin type III domain (3.7 nm end-to-end) gave a moderate FRET signal while a two-domain insert (7.0 nm) gave no FRET. We then tested unstructured proteins of various lengths, including the charged-plus-PQ domain of ZipA, the tail domain of alpha-adducin, and the C-terminal tail domain of FtsZ. The structures of these FRET constructs were also studied by electron microscopy and sedimentation. A 12 amino acid linker and the N-terminal 33 amino acids of the charged domain of the ZipA gave strong FRET signals. The C-terminal 33 amino acids of the PQ domain of the ZipA and several unstructured proteins with 66-68 amino acids gave moderate FRET signals. The 150 amino acid charged-plus-PQ construct gave a barely detectable FRET signal. FRET efficiency was calculated from the decreased donor emission to estimate the distance between donor and acceptor. The donor-acceptor distance varied for unstructured inserts of the same length, suggesting that they had variable stiffness (persistence length). We conclude that GFP-based FRET can be useful for studying intrinsically unstructured proteins, and we present a range of calibrated protein inserts to experimentally determine the distances that can be studied.

Experiments suggest that simulations may overestimate electrostatic contributions to the mechanical stability of a fibronectin type III domain

Steered molecular dynamics simulations have previously been used to investigate the mechanical properties of the extracellular matrix protein fibronectin. The simulations suggest that the mechanical stability of the tenth type III domain from fibronectin (FNfn10) is largely determined by a number of critical hydrogen bonds in the peripheral strands. Interestingly, the simulations predict that lowering the pH from 7 to ∼4.7 will increase the mechanical stability of FNfn10 significantly (by ∼33 %) due to the protonation of a few key acidic residues in the A and B strands. To test this simulation prediction, we used single-molecule atomic force microscopy (AFM) to investigate the mechanical stability of FNfn10 at neutral pH and at lower pH where these key residues have been shown to be protonated. Our AFM experimental results show no difference in the mechanical stability of FNfn10 at these different pH values. These results suggest that some simulations may overestimate the role played by electrostatic interactions in determining the mechanical stability of proteins.

Designing an extracellular matrix protein with enhanced mechanical stability

The extracellular matrix proteins tenascin and fibronectin experience significant mechanical forces in vivo. Both contain a number of tandem repeating homologous fibronectin type III (fnIII) domains, and atomic force microscopy experiments have demonstrated that the mechanical strength of these domains can vary significantly. Previous work has shown that mutations in the core of an fnIII domain from human tenascin (TNfn3) reduce the unfolding force of that domain significantly: The composition of the core is apparently crucial to the mechanical stability of these proteins. Based on these results, we have used rational redesign to increase the mechanical stability of the 10th fnIII domain of human fibronectin, FNfn10, which is directly involved in integrin binding. The hydrophobic core of FNfn10 was replaced with that of the homologous, mechanically stronger TNfn3 domain. Despite the extensive substitution, FNoTNc retains both the three-dimensional structure and the cell adhesion activity of FNfn10. Atomic force microscopy experiments reveal that the unfolding forces of the engineered protein FNoTNc increase by approximately 20% to match those of TNfn3. Thus, we have specifically designed a protein with increased mechanical stability. Our results demonstrate that core engineering can be used to change the mechanical strength of proteins while retaining functional surface interactions.

Conformational dynamics in loop swap mutants of homologous fibronectin type III domains

Fibronectin type III (FN-III) domains are autonomously folded modules found in a variety of multidomain proteins. The 10th FN-III domain from fibronectin (fnFN10) and the 3rd FN-III domain from tenascin-C (tnFN3) have 27% sequence identity and the same overall fold; however, the CC' loop has a different pattern of backbone hydrogen bonds and the FG loop is longer in fnFN10 compared to tnFN3. To examine the influence of length, sequence, and context in determining dynamical properties of loops, CC' and FG loops were swapped between fnFN10 and tnFN3 to generate four mutant proteins and backbone conformational dynamics on ps-ns and mus-ms timescales were characterized by solution (15)N-NMR spin relaxation spectroscopy. The grafted loops do not strongly perturb the properties of the protein scaffold; however, specific effects of the mutations are observed for amino acids that are proximal in space to the sites of mutation. The amino acid sequence primarily dictates conformational dynamics when the wild-type and grafted loop have the same length, but both sequence and context contribute to conformational dynamics when the loop lengths differ. The results suggest that changes in conformational dynamics of mutant proteins must be considered in both theoretical studies and protein design efforts.

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.

Design, expression, and stability of a diverse protein library based on the human fibronectin type III domain

Protein libraries based on natural scaffolds enable the generation of novel molecular tools and potential therapeutics by directed evolution. Here, we report the design and construction of a high complexity library (30 x 10(13) sequences) based on the 10th fibronectin type III domain of human fibronectin (10FnIII). We examined the bacterial expression characteristics and stability of this library using a green fluorescent protein (GFP)-reporter screen, SDS-PAGE analysis, and chemical denaturation, respectively. The high throughput GFP reporter screen demonstrates that a large fraction of our library expresses significant levels of soluble protein in bacteria. However, SDS-PAGE analysis of expression cultures indicates the ratio of soluble to insoluble protein expressed varies greatly for randomly chosen library members. We also tested the stabilities of several representative variants by guanidinium chloride denaturation. All variants tested displayed cooperative unfolding transitions similar to wild-type, and two exhibited free energies of unfolding equal to wild-type 10FnIII. This work demonstrates the utility of GFP-based screening as a tool for analysis of high-complexity protein libraries. Our results indicate that a vast amount of protein sequence space surrounding the 10FnIII scaffold is accessible for the generation of novel functions by directed as well as natural evolution.

Using Model Proteins to Quantify the Effects of Pathogenic Mutations in Ig-like Proteins

It has proved impossible to purify some proteins implicated in disease in sufficient quantities to allow a biophysical characterization of the effect of pathogenic mutations. To overcome this problem we have analyzed 37 different disease-causing mutations located in the L1 and IL2Rgamma proteins in well characterized related model proteins in which mutations that are identical or equivalent to pathogenic mutations were introduced. We show that data from these models are consistent and that changes in stability observed can be correlated to severity of disease, to correct trafficking within the cell and to in vitro ligand binding studies. Interestingly, we find that any mutations that cause a loss of stability of more than 2 kcal/mol are severely debilitating, even though some model proteins with these mutations can be easily expressed and analyzed. Furthermore we show that the severity of mutation can be predicted by a DeltaDeltaG(evolution) scale, a measure of conservation. Our results demonstrate that model proteins can be used to analyze disease-causing mutations when wild-type proteins are not stable enough to carry mutations for biophysical analysis.

High-affinity fragment complementation of a fibronectin type III domain and its application to stability enhancement

The tenth fibronectin type III (FN3) domain of human fibronectin (FNfn10), a prototype of the ubiquitous FN3 domain, is a small, monomeric beta-sandwich protein. In this study, we have bisected FNfn10 in each loop to generate a total of six fragment pairs. We found that fragment pairs bisected at multiple loops of FNfn10 show complementation in vivo as tested with a yeast two-hybrid system. The dissociation constant of these fragment pairs determined in vitro were as low as 3 nM, resulting in one of the tightest fragment complementation systems reported so far. Furthermore, we show that the affinity of fragment complementation is correlated with the stability of the uncut parent protein. Exploring this correlation, we screened a yeast two-hybrid library of one fragment and identified mutations that suppress the effect of a destabilizing mutation in the other fragment. One of the identified mutations significantly increased the stability of the uncut wild-type protein, proving that fragment complementation can be used as a novel strategy for the selection of proteins with enhanced stability.

Antibody mimics based on human fibronectin type three domain engineered for thermostability and high-affinity binding to vascular endothelial growth factor receptor two

The tenth human fibronectin type three domain ((10)Fn3) is a small (10 kDa), extremely stable and soluble protein with an immunoglobulin-like fold, but without cysteine residues. Selections from (10)Fn3-based libraries of proteins with randomized loops have yielded high-affinity, target-specific antibody mimics. However, little is known about the biophysical properties of such antibody mimics, which will determine their suitability for in vitro and medical applications. We characterized target binding and biophysical properties of two related (10)Fn3-based antibody mimics that bind vascular endothelial growth factor receptor two (VEGF-R2). The first antibody mimic, which has a dissociation constant (K(d)) of 13 nM, is highly stable [melting temperature (T(m))=62 degrees C] and soluble, whereas the second, which binds VEGF-R2 with 40 x higher affinity, is less stable (T(m) < 40 degrees C) and relatively insoluble. We used our understanding of these two (10)Fn3 derivatives and of wild-type (10)Fn3 structure to engineer the next generation of antibody mimics, which have an improved combination of high affinity (K(d)=0.59 nM), stability (T(m)=53 degrees C) and solubility. Our findings illustrate that (10)Fn3-based antibody mimics can be engineered for favorable biophysical properties even when 20% of the wild-type (10)Fn3 sequence is mutated in order to satisfy target-binding requirements.

Mechanical Unfolding of TNfn3: The Unfolding Pathway of a fnIII Domain Probed by Protein Engineering, AFM and MD Simulation

Protein engineering Phi-value analysis combined with single molecule atomic force microscopy (AFM) was used to probe the molecular basis for the mechanical stability of TNfn3, the third fibronectin type III domain from human tenascin. This approach has been adopted previously to solve the forced unfolding pathway of a titin immunoglobulin domain, TI I27. TNfn3 and TI I27 are members of different protein superfamilies and have no sequence identity but they have the same beta-sandwich structure consisting of two antiparallel beta-sheets. TNfn3, however, unfolds at significantly lower forces than TI I27. We compare the response of these proteins to mechanical force. Mutational analysis shows that, as is the case with TI I27, TNfn3 unfolds via a force-stabilised intermediate. The key event in forced unfolding in TI I27 is largely the breaking of hydrogen bonds and hydrophobic interactions between the A' and G-strands. The mechanical Phi-value analysis and molecular dynamics simulations reported here reveal that significantly more of the TNfn3 molecule contributes to its resistance to force. Both AFM experimental data and molecular dynamics simulations suggest that the rate-limiting step of TNfn3 forced unfolding reflects a transition from the extended early intermediate to an aligned intermediate state. As well as losses of interactions of the A and G-strands and associated loops there are rearrangements throughout the core. As was the case for TI I27, the forced unfolding pathway of TNfn3 is different from that observed in denaturant studies in the absence of force.

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.

Comparison of the transition states for folding of two Ig-like proteins from different superfamilies

In the "fold approach" proteins with a similar fold but different sequences are compared in order to investigate the relationship between native state structure and folding behaviour. Here we compare the properties of the transition states for folding of TI I27, the 27th immunoglobulin domain from human cardiac titin, and that of TNfn3, the third fibronectin type III domain from human tenascin. Experimental phi-values were used as restraints in molecular dynamics simulations to determine the structures that make up the transition state ensembles (TSEs) for folding of the two proteins. The restrained simulations that we present allow a detailed structural comparison of the two TSEs to be made. Further calculations show explicitly that for both proteins the formation of the interactions involving the residues in the folding nucleus is sufficient for the establishment of the topology of the Ig-like fold. We found that, although the folding nuclei of the two proteins are similar, the packing of the folding nucleus of TI I27 is much tighter than that of TNfn3, reflecting the higher experimental phi-values and beta(T) (Tanford Beta) of TI I27. These results suggest that the folding nucleus can be significantly deformed to accommodate extensive sequence variation while conserving the same folding mechanism.

Biophysical investigations of engineered polyproteins: implications for force data

Dynamic force spectroscopy is rapidly becoming a standard biophysical technique. Significant advances in the methods of analysis of force data have resulted in ever more complex systems being studied. The use of cloning systems to produce homologous tandem repeats rather than the use of endogenous multidomain proteins has facilitated these developments. What is poorly addressed are the physical properties of these constructed polyproteins. Are the properties of the individual domains in the construct independent of one another or attenuated by adjacent domains? We present data for a construct of eight fibronectin type III domains from the human form of tenascin that exhibits approximately 1 kcal mol(-1) increase in stability compared to the monomer. This effect is salt and pH dependent, suggesting that the stabilization results from electrostatic interactions, possibly involving charged residues at the interfaces of the domains. Kinetic analysis shows that this stabilization reflects a slower unfolding rate. Clearly, if domain-domain interactions affect the unfolding force, this will have implications for the comparison of absolute forces between types of domains. Mutants of the tenascin 8-mer construct exhibit the same change in stability as that observed for the corresponding mutation in the monomer. And when Phi-values are calculated for the 8-mer construct, the pattern is similar to that observed for the monomer. Therefore, mutational analyses to resolve mechanical unfolding pathways appear valid. Importantly, we show that interactions between the domains may be masked by changes in experimental conditions.

Long-range dynamic effects of point mutations propagate through side chains in the serine protease inhibitor eglin c

Long-range interactions are fundamental to protein behaviors such as cooperativity and allostery. In an attempt to understand the role protein flexibility plays in such interactions, the distribution of local fluctuations in a globular protein was monitored in response to localized, nonelectrostatic perturbations. Two valine-to-alanine mutations were introduced into the small serine protease inhibitor eglin c, and the (15)N and (2)H NMR spin relaxation properties of these variants were analyzed in terms of the Lipari-Szabo dynamics formalism and compared to those of the wild type. Significant changes in picosecond to nanosecond dynamics were observed in side chains located as much as 13 A from the point of mutation. Additionally, those residues experiencing altered dynamics appear to form contiguous surfaces within the protein. In the case of V54A, the large-to-small mutation results in a rigidification of connected residues, even though this mutation decreases the global stability. These findings suggest that dynamic perturbations arising from single mutations may propagate away from the perturbed site through networks of interacting side chains. That this is observed in eglin c, a classically nonallosteric protein, suggests that such behavior will be observed in many, if not all, globular proteins. Differences in behavior between the two mutants suggest that dynamic responses will be context-dependent.

Hydrophobic core fluidity of homologous protein domains: relation of side-chain dynamics to core composition and packing

The side-chain dynamics of methyl groups in two structurally related proteins from the fibronectin type III (fnIII) superfamily, the third fnIII domain from human tenascin (TNfn3) and the tenth fnIII domain from human fibronectin (FNfn10), have been studied by NMR spectroscopy. Side-chain order parameters reveal that the hydrophobic cores of the two proteins have substantially different mobilities. The core of TNfn3 is very dynamic, with exceptionally low order parameters for the most deeply buried residues, while that of FNfn10 is more like those of other proteins which have been studied with this technique, having a relatively rigid core with uniformly distributed dynamics. The unusually dynamic core of TNfn3 appears to be related to its amino acid composition, which makes it more fluid-like. A further explanation for the mobility of the TNfn3 core may be found in the negative correlation between the order parameter and excess packing volume, which shows that the core of TNfn3 is less densely packed and consequently has lower methyl order parameters for its buried residues. Rotameric transitions, presumably facilitated by the lower packing density, appear to make an important contribution to lowering the order parameters, and have been probed by measuring three-bond scalar couplings. Overall, although backbone dynamics is generally similar for proteins with the same topology on a fast time scale (picoseconds to nanoseconds), this study shows that a single fold can accommodate a wide variation in the dynamic properties of its core.

Comparison of the transition state ensembles for folding of Im7 and Im9 determined using all-atom molecular dynamics simulations with ϕ value restraints

Delineation of the structural properties of transition states is key to deriving models for protein folding. Here we describe the structures of the transition states of the bacterial immunity proteins Im7 and Im9 obtained by all-atom molecular dynamics simulations with phi value restraints derived from protein engineering experiments. This pair of proteins is of special interest because, at pH 7 and 10 degrees C, Im7 folds via an intermediate while Im9 folds with a two-state transition. The structures of the transition states for Im7 and Im9, together with their radii of gyration and distances from the native state, are similar. The typical distance between any two members of the transition state ensemble of both proteins is large, with that of Im9 nearly twice that of Im7. Thus, a broad range of structures make up the transition state ensembles of these proteins. The ensembles satisfy the set of rather low phi values and yet are consistent with high beta(T) values (> 0.85 for both proteins). For both Im7 and Im9 the inter-helical angles are highly variable in the transition state ensembles, although the native contacts between helices I and IV are well conserved. By measuring the distribution of the accessible surface area for each residue we show that the hydrophobic residues that are buried in the native state remain buried in the transition state, corresponding to a hydrophobic collapse to a relatively ordered globule. The data provide new insights into the structural properties of the transition states of these proteins at an atomic level of detail and show that molecular dynamics simulations with phi value restraints can significantly enhance the knowledge of the transition state ensembles (TSE) provided by the experimental phi values alone.

The relationship between conservation, thermodynamic stability, and function in the SH3 domain hydrophobic core

To investigate the relationships between sequence conservation, protein stability, and protein function, we have measured the thermodynamic stability, folding kinetics, and in vitro peptide-binding activity of a large number of single-site substitutions in the hydrophobic core of the Fyn SH3 domain. Comparison of these data to that derived from an analysis of a large alignment of SH3 domain sequences revealed a very good correlation between the distinct pattern of conservation observed at each core position and the thermodynamic stability of mutants. Conservation was also found to correlate well with the unfolding rates of mutants, but not to the folding rates, suggesting that evolution selects more strongly for optimal native state packing interactions than for maximal folding rates. Structural analysis suggests that residue-residue core packing interactions are very similar in all SH3 domains, which provides an explanation for the correlation between conservation and mutant stability effects studied in a single SH3 domain. We also demonstrate a correlation between stability and the in vivo activity of mutants, and between conservation and activity. However, the relationship between conservation and activity was very strong only for the three most conserved hydrophobic core positions. The weaker correlation between activity and conservation seen at the other seven core positions indicates that maintenance of protein stability is the dominant selective pressure at these positions. In general, the pattern of conservation at hydrophobic core positions appears to arise from conserved packing constraints, and can be effectively utilized to predict the destabilizing effects of amino acid substitutions.

Self-consistent determination of the transition state for protein folding: application to a fibronectin type III domain

We present a general approach in which theory and experiments are combined in an iterative manner to provide a detailed description of the transition state ensemble (TSE) for folding. The method is illustrated by applying it to TNfn3, a fibronectin type III domain protein. In the first iteration, a coarse-grained determination of the TSE is carried out by using a limited set of experimental phi values as constraints in a molecular dynamics sampling simulation. The resulting model of the TSE is used to determine the additional residues whose phi value measurement would provide the most information for refining the TSE. Successive iterations with an increasing number of phi value measurements are carried out until no further changes in the properties of the TSE are detected or there are no additional residues whose phi values can be measured. In the study of TNfn3 three iterations were necessary to achieve self-consistency. A retrospective application of the method can be used to determine the accuracy of the TSE results and to find "key residues" for folding, i.e., those that are most important for the formation of the TSE. The approach reported here is an efficient method for finding the structures that make up the TSEs for protein folding. Its use will improve future efforts for their experimental determination and refinement.

Exploring the potential of the monobody scaffold: effects of loop elongation on the stability of a fibronectin type III domain

The tenth fibronectin type III domain of human fibronectin (FNfn10) is a small, monomeric beta-sandwich protein, similar to the immunoglobulins. We have developed small antibody mimics, 'monobodies', using FNfn10 as a scaffold. We initially altered two loops of FNfn10 that are structurally equivalent to two of the hypervariable loops of the immunoglobulin domain. In order to assess the possibility of utilizing other loops in FNfn10 for target binding, we determined the effects of the elongation of each loop on the conformational stability of FNfn10. We found that all six loops of FNfn10 allowed the introduction of four glycine residues while retaining the global fold. Insertions in the AB and FG loops exhibited very small degrees of destabilization, comparable to or less than predicted entropic penalties due to the elongation, suggesting the absence of stabilizing interactions in these loops in wild-type FNfn10. Insertions in the BC, CD and DE loops, respectively, resulted in modest destabilization. In contrast, the EF loop elongation was highly destabilizing, consistent with previous studies showing the presence of stabilizing interactions in this loop. These results suggest that all loops, except for the EF loop, can be used for engineering a binding site, thus demonstrating excellent properties of the monobody scaffold.

Validity of Gō models: comparison with a solvent-shielded empirical energy decomposition

Do Gō-type model potentials provide a valid approach for studying protein folding? They have been widely used for this purpose because of their simplicity and the speed of simulations based on their use. The essential assumption in such models is that only contact interactions existing in the native state determine the energy surface of a polypeptide chain, even for non-native configurations sampled along folding trajectories. Here we use an all-atom molecular mechanics energy function to investigate the adequacy of Gō-type potentials. We show that, although the contact approximation is accurate, non-native contributions to the energy can be significant. The assumed relation between residue-residue interaction energies and the number of contacts between them is found to be only approximate. By contrast, individual residue energies correlate very well with the number of contacts. The results demonstrate that models based on the latter should give meaningful results (e.g., as used to interpret phi values), whereas those that depend on the former are only qualitative, at best.

Determination of a transition state at atomic resolution from protein engineering data

We present a method for determining the structure of the transition state ensemble (TSE) of a protein by using phi values derived from protein engineering experiments as restraints in molecular dynamics simulations employing a realistic all-atom molecular mechanics energy function. The method uses a biasing potential to select an ensemble of structures having phi values in agreement with the experimental data set. An application to acylphosphatase (AcP), a protein for which phi values have been measured for 24 out of 98 residues, illustrates the approach. The properties of the TSE determined in this way are compared with those of a coarse-grained model obtained using a Monte Carlo (MC) sampling method based on a C(alpha) representation of the structure. The two TSEs determined at different structural resolution are consistent and complementary. While the C(alpha) model allows better sampling of the conformation space occupied by the transition state, the all-atom model offers a more detailed description of the structural and energetic properties of the conformations included in the TSE. The combination of low-resolution C(alpha) results with all-atom molecular dynamics simulations provides a powerful and general method for determining the nature of TSEs from protein engineering data.