Folding and stability of a fibronectin type III domain of human tenascin

Design and Synthesis of Monobody Variants with Low Immunogenicity

Mirror-image proteins (d-proteins) are promising scaffolds for drug discovery because of their high proteolytic stability and low immunogenic properties. Facile and reproducible processes for the preparation of functional d-proteins are required for their application in therapeutic biologics. In this study, we designed and synthesized a novel monobody variant with two cysteine substitutions that facilitate the synthetic process via sequential native chemical ligations and improve protein stability by disulfide bond formation. The synthetic anti-GFP monobody in this model study exhibited good binding affinity to the target enhanced green fluorescent protein. In vivo administration of the synthetic anti-GFP monobody (l-monobody) to mice induced antidrug antibody (ADA) production, whereas no ADA production was observed following immunization with the mirror-image anti-GFP monobody (d-monobody). These results suggest that the synthetic d-monobody is a non-antibody protein scaffold with low immunogenic properties.

De Novo Design of a β-Helix Tau Protein Scaffold: An Oligomer-Selective Vaccine Immunogen Candidate for Alzheimer's Disease

Tau pathology is associated with many neurodegenerative disorders, including Alzheimer's disease (AD), where the spatio-temporal pattern of tau neurofibrillary tangles strongly correlates with disease progression, which motivates therapeutics selective for misfolded tau. Here, we introduce a new avidity-enhanced, multi-epitope approach for protein-misfolding immunogen design, which is predicted to mimic the conformational state of an exposed epitope in toxic tau oligomers. A predicted oligomer-selective tau epitope 343KLDFK347 was scaffolded by designing a β-helix structure that incorporated multiple instances of the 16-residue tau fragment 339VKSEKLDFKDRVQSKI354. Large-scale conformational ensemble analyses involving Jensen-Shannon Divergence and the embedding depth D showed that the multi-epitope scaffolding approach, employed in designing the β-helix scaffold, was predicted to better discriminate toxic tau oligomers than other "monovalent" strategies utilizing a single instance of an epitope for vaccine immunogen design. Using Rosetta, 10,000 sequences were designed and screened for the linker portions of the β-helix scaffold, along with a C-terminal stabilizing α-helix that interacts with the linkers, to optimize the folded structure and stability of the scaffold. Structures were ranked by energy, and the lowest 1% (82 unique sequences) were verified using AlphaFold. Several selection criteria involving AlphaFold are implemented to obtain a lead-designed sequence. The structure was further predicted to have free energetic stability by using Hamiltonian replica exchange molecular dynamics (MD) simulations. The synthesized β-helix scaffold showed direct binding in surface plasmon resonance (SPR) experiments to several antibodies that were raised to the structured epitope using a designed cyclic peptide. Moreover, the strength of binding of these antibodies to in vitro tau oligomers correlated with the strength of binding to the β-helix construct, suggesting that the construct presents an oligomer-like conformation and may thus constitute an effective oligomer-selective immunogen.

Dry Molten Globule-Like Intermediates in Protein Folding, Function, and Disease

The performance of a protein depends on its correct folding to the final functional native form. Hence, understanding the process of protein folding has remained an important field of research for the scientific community for the past five decades. Two important intermediate states, namely, wet molten globule (WMG) and dry molten globule (DMG), have emerged as critical milestones during protein folding-unfolding reactions. While much has been discussed about WMGs as a common unfolding intermediate, the evidence for DMGs has remained elusive owing to their near-native features, which makes them difficult to probe using global structural probes. This Review puts together the available literature and new evidence on DMGs to give a broader perspective on the universality of DMGs and discuss their significance in protein folding, function, and disease.

Surface salt bridges contribute to the extreme thermal stability of an FN3-like domain from a thermophilic bacterium

This study uses differential scanning calorimetry, X-ray crystallography, and molecular dynamics simulations to investigate the structural basis for the high thermal stability (melting temperature 97.5°C) of a FN3-like protein domain from thermophilic bacteria Thermoanaerobacter tengcongensis (FN3tt). FN3tt adopts a typical FN3 fold with a three-stranded beta sheet packing against a four-stranded beta sheet. We identified three solvent exposed arginine residues (R23, R25, and R72), which stabilize the protein through salt bridge interactions with glutamic acid residues on adjacent strands. Alanine mutation of the three arginine residues reduced melting temperature by up to 22°C. Crystal structures of the wild type (WT) and a thermally destabilized (∆Tm -19.7°C) triple mutant (R23L/R25T/R72I) were found to be nearly identical, suggesting that the destabilization is due to interactions of the arginine residues. Molecular dynamics simulations showed that the salt bridge interactions in the WT were stable and provided a dynamical explanation for the cooperativity observed between R23 and R25 based on calorimetry measurements. In addition, folding free energy changes computed using free energy perturbation molecular dynamics simulations showed high correlation with melting temperature changes. This work is another example of surface salt bridges contributing to the enhanced thermal stability of thermophilic proteins. The molecular dynamics simulation methods employed in this study may be broadly useful for in silico surface charge engineering of proteins.

Engineering Calreticulin-Targeting Monobodies to Detect Immunogenic Cell Death in Cancer Chemotherapy

Surface-exposed calreticulin (ecto-CRT) plays a crucial role in the phagocytic removal of apoptotic cells during immunotherapy. Ecto-CRT is an immunogenic signal induced in response to treatment with chemotherapeutic agents such as doxorubicin (DOX) and mitoxantrone (MTX), and two peptides (KLGFFKR (Integrin-α) and GQPMYGQPMY (CRT binding peptide 1, Hep-I)) are known to specifically bind CRT. To engineer CRT-specific monobodies as agents to detect immunogenic cell death (ICD), we fused these peptide sequences at the binding loops (BC and FG) of human fibronectin domain III (FN3). CRT-specific monobodies were purified from E. coli by affinity chromatography. Using these monobodies, ecto-CRT was evaluated in vitro, in cultured cancer cell lines (CT-26, MC-38, HeLa, and MDA-MB-231), or in mice after anticancer drug treatment. Monobodies with both peptide sequences (CRT3 and CRT4) showed higher binding to ecto-CRT than those with a single peptide sequence. The binding affinity of the Rluc8 fusion protein-engineered monobodies (CRT3-Rluc8 and CRT4-Rluc8) to CRT was about 8 nM, and the half-life in serum and tumor tissue was about 12 h. By flow cytometry and confocal immunofluorescence of cancer cell lines, and by in vivo optical bioluminescence imaging of tumor-bearing mice, CRT3-Rluc8 and CRT4-Rluc8 bound specifically to ecto-CRT and effectively detected pre-apoptotic cells after treatment with ICD-inducing agents (DOX and MTX) but not a non-ICD-inducing agent (gemcitabine). Using CRT-specific monobodies, it is possible to detect ecto-CRT induction in cancer cells in response to drug exposure. This technique may be used to predict the therapeutic efficiency of chemo- and immuno-therapeutics early during anticancer treatment.

Type III secretion system effector proteins are mechanically labile

Multiple gram-negative bacteria encode type III secretion systems (T3SS) that allow them to inject effector proteins directly into host cells to facilitate colonization. To be secreted, effector proteins must be at least partially unfolded to pass through the narrow needle-like channel (diameter <2 nm) of the T3SS. Fusion of effector proteins to tightly packed proteins-such as GFP, ubiquitin, or dihydrofolate reductase (DHFR)-impairs secretion and results in obstruction of the T3SS. Prior observation that unfolding can become rate-limiting for secretion has led to the model that T3SS effector proteins have low thermodynamic stability, facilitating their secretion. Here, we first show that the unfolding free energy ([Formula: see text]) of two Salmonella effector proteins, SptP and SopE2, are 6.9 and 6.0 kcal/mol, respectively, typical for globular proteins and similar to published [Formula: see text] for GFP, ubiquitin, and DHFR. Next, we mechanically unfolded individual SptP and SopE2 molecules by atomic force microscopy (AFM)-based force spectroscopy. SptP and SopE2 unfolded at low force (F _unfold ≤ 17 pN at 100 nm/s), making them among the most mechanically labile proteins studied to date by AFM. Moreover, their mechanical compliance is large, as measured by the distance to the transition state (Δx ^‡ = 1.6 and 1.5 nm for SptP and SopE2, respectively). In contrast, prior measurements of GFP, ubiquitin, and DHFR show them to be mechanically robust (F _unfold > 80 pN) and brittle (Δx ^‡ < 0.4 nm). These results suggest that effector protein unfolding by T3SS is a mechanical process and that mechanical lability facilitates efficient effector protein secretion.

Conserved tryptophan mutation disrupts structure and function of immunoglobulin domain revealing unusual tyrosine fluorescence

Immunoglobulin (Ig) domains are the most prevalent protein domain structure and share a highly conserved folding pattern; however, this structural family of proteins is also the most diverse in terms of biological roles and tissue expression. Ig domains vary significantly in amino acid sequence but share a highly conserved tryptophan in the hydrophobic core of this beta-stranded protein. Palladin is an actin binding and bundling protein that has five Ig domains and plays an important role in normal cell adhesion and motility. Mutation of the core tryptophan in one Ig domain of palladin has been identified in a pancreatic cancer cell line, suggesting a crucial role for this sole tryptophan in palladin Ig domain structure, stability, and function. We found that actin binding and bundling was not completely abolished with removal of this tryptophan despite a partially unfolded structure and significantly reduced stability of the mutant Ig domain as shown by circular dichroism investigations. In addition, this mutant palladin domain displays a tryptophan-like fluorescence attributed to an anomalous tyrosine emission at 341 nm. Our results indicate that this emission originates from a tyrosinate that may be formed in the excited ground state by proton transfer to a nearby aspartic acid residue. Furthermore, this study emphasizes the importance of tryptophan in protein structural stability and illustrates how tyrosinate emission contributions may be overlooked during the interpretation of the fluorescence properties of proteins.

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.

Structure and unfolding of the third type III domain from human fibronectin

Fibronectin is a modular extracellular matrix protein that is essential for vertebrate development. The third type III domain (3FN3) in fibronectin interacts with other parts of fibronectin and with anastellin, a protein fragment that causes fibronectin aggregation. 3FN3 opens readily both as an isolated domain in solution and when part of fibronectin in stretched fibrils, and it was proposed that this opening is important for anastellin binding. We determined the structure of 3FN3 using nuclear magnetic resonance spectroscopy, and we investigated its stability, folding, and unfolding. Similar to most other FN3 domains, 3FN3 contains two antiparallel β-sheets that are composed of three (A, B, and E) and four (C, D, F, and G) β-strands, respectively, and are held together by a conserved hydrophobic interface. cis-trans isomerization of P847 at the end of β-strand C leads to observable conformational heterogeneity in 3FN3, with a cis peptide bond present in almost one-quarter of the molecules. The chemical stability of 3FN3 is relatively low, but the folding rate constant in the absence of denaturant is in the same range as those of other, more stable FN3 domains. Interestingly, the unfolding rate constant in the absence of denaturant is several orders of magnitude higher than the unfolding rate constants of other FN3 domains investigated to date. This unusually fast rate is comparable to the rate of binding of 3FN3 to anastellin at saturating anastellin concentrations, consistent with the model in which 3FN3 has to unfold to interact with anastellin.

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.

Backtracking due to residual structure in the unfolded state changes the folding of the third fibronectin type III domain from tenascin-C

Residual structure in the unfolded state of a protein may play a crucial role in folding and stability. In the present study, using an all (heavy)-atom structure based model and replica exchange molecular dynamics simulations, we explored the folding landscape of the third fibronectin type III domain from tenascin-C (TNfn3). Specifically, both the wild type (WT) and a variant with two additional amino acids, Gly-Leu (GL), at the C-terminus (WT(+GL)) were studied. We found that, although both domains of TNfn3 are topologically frustrated, the early formation of the native contacts from the C-terminal end of WT(+GL) causes more "backtracking" than in the WT. As a result, the WT exhibits a two-state folding behavior with a broad transition-state ensemble, whereas the WT(+GL) folds through a metastable intermediate state. Furthermore, our study confirmed that the core of both proteins is conformationally heterogeneous and noncompact, and folds late mainly due to backtracking of the part of the core. Finally, in agreement with the previous experimental studies, our results clearly demonstrated distinct thermodynamic behavior of the two proteins with WT(+GL) being more stable.

Design of novel FN3 domains with high stability by a consensus sequence approach

The use of consensus design to produce stable proteins has been applied to numerous structures and classes of proteins. Here, we describe the engineering of novel FN3 domains from two different proteins, namely human fibronectin and human tenascin-C, as potential alternative scaffold biotherapeutics. The resulting FN3 domains were found to be robustly expressed in Escherichia coli, soluble and highly stable, with melting temperatures of 89 and 78°C, respectively. X-ray crystallography was used to confirm that the consensus approach led to a structure consistent with the FN3 design despite having only low-sequence identity to natural FN3 domains. The ability of the Tenascin consensus domain to withstand mutations in the loop regions connecting the β-strands was investigated using alanine scanning mutagenesis demonstrating the potential for randomization in these regions. Finally, rational design was used to produce point mutations that significantly increase the stability of one of the consensus domains. Together our data suggest that consensus FN3 domains have potential utility as alternative scaffold therapeutics.

Improvement of structure-based potentials for protein folding by native and nonnative hydrogen bonds

Pure Gō models (where every native interaction equally stabilizes the folded state) have widely proved their convenience in the computational investigation of protein folding. However, a chemistry-based description of the real interactions also provides a desirable tune in the analysis of the folding process, and thus some hybrid Gō potentials that combine both aspects have been proposed. Among all the noncovalent interactions that contribute to protein folding, hydrogen bonds are the only ones with a partial covalent character. This feature makes them directional and, thus, more difficult to model as part of the coarse-grained descriptions that are typically employed in Gō models. Thanks to a simplified but rigorous representation of backbone hydrogen bonds that we have recently proposed, we present in this article a combined potential (Gō + backbone hydrogen bond) to study the thermodynamics of protein folding in the frame of very simple simulation models. We show that the explicit inclusion of hydrogen bonds leads to a systematic improvement in the description of protein folding. We discuss a representative set of examples (from two-state folders to downhill proteins, with different types of native structures) that reveal a relevant agreement with experimental data.

Kinetic partitioning mechanism governs the folding of the third FnIII domain of tenascin-C: evidence at the single-molecule level

Statistical mechanics and molecular dynamics simulations proposed that the folding of proteins can follow multiple parallel pathways on a rugged energy landscape from unfolded state en route to their folded native states. Kinetic partitioning mechanism is one of the possible mechanisms underlying such complex folding dynamics. Here, we use single-molecule atomic force microscopy technique to directly probe the multiplicity of the folding pathways of the third fibronectin type III domain from the extracellular matrix protein tenascin-C (TNfn3). By stretching individual (TNfn3)(8) molecules, we forced TNfn3 domains to undergo mechanical unfolding and refolding cycles, allowing us to directly observe the folding pathways of TNfn3. We found that, after being mechanically unraveled and then relaxed to zero force, TNfn3 follows multiple parallel pathways to fold into their native states. The majority of TNfn3 fold into the native state in a simple two-state fashion, while a small percentage of TNfn3 were found to be trapped into kinetically stable folding intermediate states with well-defined three-dimensional structures. Furthermore, the folding of TNfn3 was also influenced by its neighboring TNfn3 domains. Complex misfolded states of TNfn3 were observed, possibly due to the formation of domain-swapped dimeric structures. Our studies revealed the ruggedness of the folding energy landscape of TNfn3 and provided direct experimental evidence that the folding dynamics of TNfn3 are governed by the kinetic partitioning mechanism. Our results demonstrated the unique capability of single-molecule AFM to probe the folding dynamics of proteins at the single-molecule level.

Fibronectin forms the most extensible biological fibers displaying switchable force-exposed cryptic binding sites

Rather than maximizing toughness, as needed for silk and muscle titin fibers to withstand external impact, the much softer extracellular matrix fibers made from fibronectin (Fn) can be stretched by cell generated forces and display extraordinary extensibility. We show that Fn fibers can be extended more than 8-fold (>700% strain) before 50% of the fibers break. The Young's modulus of single fibers, given by the highly nonlinear slope of the stress-strain curve, changes orders of magnitude, up to MPa. Although many other materials plastically deform before they rupture, evidence is provided that the reversible breakage of force-bearing backbone hydrogen bonds enables the large strain. When tension is released, the nano-sized Fn domains first contract in the crowded environment of fibers within seconds into random coil conformations (molten globule states), before the force-bearing hydrogen bond networks that stabilize the domain's secondary structures are reestablished within minutes (double exponential). The exposure of cryptic binding sites on Fn type III modules increases steeply upon stretching. Thus fiber extension steadily up-regulates fiber rigidity and cryptic epitope exposure, both of which are known to differentially alter cell behavior. Finally, since stress-strain relationships cannot directly be measured in native extracellular matrix (ECM), the stress-strain curves were correlated with stretch-induced alterations of intramolecular fluorescence resonance energy transfer (FRET) obtained from trace amounts of Fn probes (mechanical strain sensors) that can be incorporated into native ECM. Physiological implications of the extraordinary extensibility of Fn fibers and contraction kinetics are discussed.

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.

Manipulating the stability of fibronectin type III domains by protein engineering

We have previously shown that a 'weak' fibronectin type III (fnIII) domain can be engineered to have enhanced mechanical strength by replacing the hydrophobic core with the core of a homologous 'strong' fnIII domain. Here we show that engineering the core is a robust method for manipulating the mechanical strength of this class of proteins. We performed an experiment that is the reverse of one described earlier. The hydrophobic core of a 'weak' domain (FNfn10) was grafted into a 'strong' fnIII domain, TNfn3. This newly engineered protein, TNoFNc, is indeed much less mechanically resistant than TNfn3. Interestingly, TNoFNc is very unstable, approximately 10 kcal mol(-1) less stable than FNfn10, yet its mechanical stability is very similar-a clear reflection of the fact that thermodynamic and mechanical stability are unrelated properties, even where they are both assumed to reflect properties of the hydrophobic core.

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.

Plasticity within the obligatory folding nucleus of an immunoglobulin-like domain

A number of β-sandwich immunoglobulin-like domains have been shown to fold using a set of structurally equivalent residues that form a folding nucleus deep within the core of the protein. Formation of this nucleus is sufficient to establish the complex Greek key topology of the native state. These nucleating residues are highly conserved within the immunoglobulin superfamily, but are less well conserved in the fibronectin type III (fnIII) superfamily, where the requirement is simply to have four interacting hydrophobic residues. However, there are rare examples where this nucleation pattern is absent. In this study, we have investigated the folding of a novel member of the fnIII superfamily whose nucleus appears to lack one of the four buried hydrophobic residues. We show that the folding mechanism is unaltered, but the folding nucleus has moved within the hydrophobic core.

Effect of finite size on cooperativity and rates of protein folding

We analyze the dependence of cooperativity of the thermal denaturation transition and folding rates of globular proteins on the number of amino acid residues, N, using lattice models with side chains, off-lattice Go models, and the available experimental data. A dimensionless measure of cooperativity, Omega(c) (0 < Omega(c) < infinity), scales as Omega(c) approximately N(zeta). The results of simulations and the analysis of experimental data further confirm the earlier prediction that zeta is universal with zeta = 1 + gamma, where exponent gamma characterizes the susceptibility of a self-avoiding walk. This finding suggests that the structural characteristics in the denaturated state are manifested in the folding cooperativity at the transition temperature. The folding rates k(F) for the Go models and a dataset of 69 proteins can be fit using k(F) = k(F)0 exp(-cN(beta)). Both beta = 1/2 and 2/3 provide a good fit of the data. We find that k(F) = k(F)0 exp(-cN(1/2)), with the average (over the dataset of proteins) k(F)0 approximately (0.2 micros)(-1) and c approximately 1.1, can be used to estimate folding rates to within an order of magnitude in most cases. The minimal models give identical N dependence with c approximately 1. The prefactor for off-lattice Go models is nearly 4 orders of magnitude larger than the experimental value.

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.

Unfolding of beta-sheet proteins in SDS

Beta-sheet proteins are particularly resistant to denaturation by sodium dodecyl sulfate (SDS). Here we compare unfolding of two beta-sandwich proteins TNfn3 and TII27 in SDS. The two proteins show different surface electrostatic potential. Correspondingly, TII27 unfolds below the critical micelle concentration via the formation of hemimicelles on the protein surface, whereas TNfn3 only unfolds around the critical micelle concentration. Isothermal titration calorimetry confirms that unfolding of TII27 sets in at lower SDS concentrations, although the total number of bound SDS molecules is similar at the end of unfolding. In mixed micelles with the nonionic detergent dodecyl maltoside, where the concentration of monomeric SDS is insignificant, the behavior of the two proteins converges. TII27 unfolds more slowly than TNfn3 in SDS and follows a two-mode behavior. Additionally TNfn3 shows inhibition of SDS unfolding at intermediate SDS concentrations. Mutagenic analysis suggests that the overall unfolding mechanism is similar to that observed in denaturant for both proteins. Our data confirm the kinetic robustness of beta-sheet proteins toward SDS. We suggest this is related to the inability of SDS to induce significant amounts of alpha-helix structure in these proteins as part of the denaturation process, forcing the protein to denature by global rather than local unfolding.

Influence of the Internal Disulfide Bridge on the Folding Pathway of the CL Antibody Domain

Disulfide bridges are one of the most important factors stabilizing the native structure of a protein. Whereas the basis for their stabilizing effect is well understood, their role in a protein folding reaction still seems to require further attention. We used the constant domain of the antibody light chain (C(L)), a representative of the ubiquitous immunoglobulin (Ig)-superfamily, to delineate the kinetic role of its single buried disulfide bridge. Independent of its redox state, the monomeric C(L) domain adopts a typical Ig-fold under native conditions and does not retain significant structural elements when unfolded. Interestingly, its folding pathway is strongly influenced by the disulfide bridge. The more stable oxidized protein folds via a highly structured on-pathway intermediate, whereas the destabilized reduced protein populates a misfolded off-pathway species on its way to the native state. In both cases, the formation of the intermediate species is shown to be independent of the isomerization state of the Tyr(141)-Pro(142) bond. Our results demonstrate that the internal disulfide bridge in an antibody domain restricts the folding pathway by bringing residues of the folding nucleus into proximity thus facilitating the way to the native state.

BPPred: a Web-based computational tool for predicting biophysical parameters of proteins

We exploit the availability of recent experimental data on a variety of proteins to develop a Web-based prediction algorithm (BPPred) to calculate several biophysical parameters commonly used to describe the folding process. These parameters include the equilibrium m-values, the length of proteins, and the changes upon unfolding in the solvent-accessible surface area, in the heat capacity, and in the radius of gyration. We also show that the knowledge of any one of these quantities allows an estimate of the others to be obtained, and describe the confidence limits with which these estimations can be made. Furthermore, we discuss how the kinetic m-values, or the Beta Tanford values, may provide an estimate of the solvent-accessible surface area and the radius of gyration of the transition state for protein folding. Taken together, these results suggest that BPPred should represent a valuable tool for interpreting experimental measurements, as well as the results of molecular dynamics simulations.

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.

Folding rate prediction using n-order contact distance for proteins with two- and three-state folding kinetics

It is a challenging task to understand the relationship between sequences and folding rates of proteins. Previous studies are found that one of contact order (CO), long-range order (LRO), total contact distance (TCD), chain topology parameter (CTP), and effective length (Leff) has a significant correlation with folding rate of proteins. In this paper, we introduce a new parameter called n-order contact distance (nOCD) and use it to predict folding rate of proteins with two- and three-state folding kinetics. A good linear correlation between the folding rate logarithm lnkf and nOCD with n=1.2, alpha=0.6 is found for two-state folders (correlation coefficient is -0.809, P-value<0.0001) and n=2.8, alpha=1.5 for three-state folders (correlation coefficient is -0.816, P-value<0.0001). However, this correlation is completely absent for three-state folders with n=1.2, alpha=0.6 (correlation coefficient is 0.0943, P-value=0.661) and for two-state folders with n=2.8, alpha=1.5 (correlation coefficient is -0.235, P-value=0.2116). We also find that the average number of contacts per residue Pm in the interval of m for two-state folders is smaller than that for three-state folders. The probability distribution P(gamma) of residue having gamma pairs of contacts fits a Gaussian distribution for both two- and three-state folders. We observe that the correlations between square radius of gyration S2 and number of residues for two- and three-state folders are both good, and the correlation coefficient is 0.908 and 0.901, and the slope of the fitting line is 1.202 and 0.795, respectively. Maybe three-state folders are more compact than two-state folders. Comparisons with nTCD and nCTP are also made, and it is found that nOCD is the best one in folding rate prediction.

Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.

Protein stability and surface electrostatics: a charged relationship

Engineering proteins to withstand a broad range of conditions continues to be a coveted objective, holding the potential to advance biomedicine, industry, and our understanding of disease. One way of achieving this goal lies in elucidating the underlying interactions that define protein stability. It has been shown that the hydrophobic effect, hydrogen bonding, and packing interactions between residues in the protein interior are dominant factors that define protein stability. The role of surface residues in protein stability has received much less attention. It has been believed that surface residues are not important for protein stability particularly because their interactions with the solvent should be similar in the native and unfolded states. In the case of surface charged residues, it was sometimes argued that solvent exposure meant that the high dielectric of the solvent will further decrease the strength of the charge-charge interactions. In this paper, we challenge the notion that the surface charged residues are not important for protein stability. We computationally redesigned sequences of five different proteins to optimize the surface charge-charge interactions. All redesigned proteins exhibited a significant increase in stability relative to their parent proteins, as experimentally determined by circular dichroism spectroscopy and differential scanning calorimetry. These results suggest that surface charge-charge interactions are important for protein stability and that rational optimization of charge-charge interactions on the protein surface can be a viable strategy for enhancing protein stability.

Entropy capacity determines protein folding

Search and study of the general principles that govern kinetics and thermodynamics of protein folding generate a new insight into the factors controlling this process. Here, based on the known experimental data and using theoretical modeling of protein folding, we demonstrate that there exists an optimal relationship between the average conformational entropy and the average energy of contacts per residue-that is, an entropy capacity-for fast protein folding. Statistical analysis of conformational entropy and number of contacts per residue for 5829 protein structures from four general structural classes (all-alpha, all-beta, alpha/beta, alpha+beta) demonstrates that each class of proteins has its own class-specific average number of contacts (class alpha/beta has the largest number of contacts) and average conformational entropy per residue (class all-alpha has the largest number of rotatable angles phi, psi, and chi per residue). These class-specific features determine the folding rates: alpha proteins are the fastest folding proteins, then follow beta and alpha+beta proteins, and finally alpha/beta proteins are the slowest ones. Our result is in agreement with the experimental folding rates for 60 proteins. This suggests that structural and sequence properties are important determinants of protein folding rates.

Domain Unfolding Plays a Role in Superfibronectin Formation

Superfibronectin (sFN) is a fibronectin (FN) aggregate that is formed by mixing FN with anastellin, a fragment of the first type III domain of FN. However, the mechanism of this aggregation has not been clear. In this study, we found that anastellin co-precipitated with FN in a ratio of approximately 4:1, anastellin:FN monomer. The primary binding site for anastellin was in the segment (III)1-3, which bound three molecules of anastellin and was able to form a precipitate without the rest of the FN molecule. Anastellin binding to (III)3 caused a conformational change in that domain that exposed a cryptic thermolysin-sensitive site. An additional anastellin binds to (III)11, where it enhances thermolysin digestion of (III)11. An engineered disulfide bond in (III)3 inhibited both aggregation and protease digestion, suggesting that the stability of (III)3 is a key factor in sFN formation. We propose a three-step model for sFN formation: 1) FN-III domains spontaneously unfold and refold; 2) anastellin binds to an unfolded domain, preventing its refolding and leaving it with exposed hydrophobic surfaces and beta-sheet edges; and 3) these exposed elements bind to similar exposed elements on other molecules, leading to aggregation. The model is consistent with our observation that the kinetics of aggregation are first order, with a reaction time of 500-700 s. Similar mechanisms may contribute to the assembly of the native FN matrix.

Influence of Trp Mutation on Native, Intermediate, and Transition States of Goat α-Lactalbumin: An Equilibrium and Kinetic Study

Equilibrium circular dichroism and kinetic stopped-flow fluorescence studies on the stability and the folding kinetics of a set of Trp to Phe mutants of goat alpha-lactalbumin (GLA) were used to characterize the native, intermediate, and transition states of these constructs. GLA contains four tryptophan residues, three of which, Trp26, Trp104, and Trp118, are located in the alpha-domain, while the fourth, Trp60, is located in the beta-domain. Trp26, Trp60, and Trp104 are part of a hydrophobic cluster, whereas Trp118 is situated in a more flexible region near the C-terminal end of the protein. In each case, the mutation leads to a reduction in the overall stability, but only for W26F and W60F is an equilibrium intermediate observed in guanidine hydrochloride-induced unfolding experiments. In kinetic refolding experiments, however, for all samples a burst phase is observed, the amplitude of which depends on the specific mutation. Refolding and unfolding kinetics can adequately be described by a sequential three-state mechanism. phi value analysis showed that the local structure around Trp26, Trp60, and Trp104 is formed in the intermediate and in the transition state of the folding reaction, while around Trp118 no persistent native contacts are observed. From these findings, we conclude that, although hydrophobicity is a major driving force for folding, minor steric changes induced by point mutation can considerably influence the overall stability and the folding process of the protein.

Fibronectin fibrillogenesis, a cell-mediated matrix assembly process

The extracellular matrix provides a framework for cell adhesion, supports cell movement, and serves to compartmentalize tissues into functional units. Fibronectin is a core component of many extracellular matrices where it regulates a variety of cell activities through direct interactions with cell surface integrin receptors. Fibronectin is synthesized by many adherent cells which then assemble it into a fibrillar network. The assembly process is integrin-dependent and fibronectin-integrin interactions initiate a step-wise process involving conformational activation of fibronectin outside and organization of the actin cytoskeleton inside. During assembly, fibronectin undergoes conformational changes that expose fibronectin-binding sites and promote intermolecular interactions needed for fibril formation. In this review, the main steps of fibronectin assembly are described and recent studies on fibronectin conformational changes are discussed.

A Statistical Model for Predicting Protein Folding Rates from Amino Acid Sequence with Structural Class Information

Prediction of protein folding rates from amino acid sequences is one of the most important challenges in molecular biology. In this work, I have related the protein folding rates with physical-chemical, energetic and conformational properties of amino acid residues. I found that the classification of proteins into different structural classes shows an excellent correlation between amino acid properties and folding rates of two- and three-state proteins, indicating the importance of native state topology in determining the protein folding rates. I have formulated a simple linear regression model for predicting the protein folding rates from amino acid sequences along with structural class information and obtained an excellent agreement between predicted and experimentally observed folding rates of proteins; the correlation coefficients are 0.99, 0.96 and 0.95, respectively, for all-alpha, all-beta and mixed class proteins. This is the first available method, which is capable of predicting the protein folding rates just from the amino acid sequence with the aid of generic amino acid properties and structural class information.

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.

Topological determinants of protein unfolding rates

For proteins that fold by two-state kinetics, the folding and unfolding processes are believed to be closely related to their native structures. In particular, folding and unfolding rates are influenced by the native structures of proteins. Thus, we focus on finding important topological quantities from a protein structure that determine its unfolding rate. After constructing graphs from protein native structures, we investigate the relationships between unfolding rates and various topological quantities of the graphs. First, we find that the correlation between the unfolding rate and the contact order is not as prominent as in the case of the folding rate and the contact order. Next, we investigate the correlation between the unfolding rate and the clustering coefficient of the graph of a protein native structure, and observe no correlation between them. Finally, we find that a newly introduced quantity, the impact of edge removal per residue, has a good overall correlation with protein unfolding rates. The impact of edge removal is defined as the ratio of the change of the average path length to the edge removal probability. From these facts, we conclude that the protein unfolding process is closely related to the protein native structure.

A flexible approach for understanding protein stability

A distance constraint model (DCM) is presented that identifies flexible regions within protein structure consistent with specified thermodynamic condition. The DCM is based on a rigorous free energy decomposition scheme representing structure as fluctuating constraint topologies. Entropy non-additivity is problematic for naive decompositions, limiting the success of heat capacity predictions. The DCM resolves non-additivity by summing over independent entropic components determined by an efficient network-rigidity algorithm. A minimal 3-parameter DCM is demonstrated to accurately reproduce experimental heat capacity curves. Free energy landscapes and quantitative stability-flexibility relationships are obtained in terms of global flexibility. Several connections to experiment are made.

Time-dependent conformational changes in fibrinogen measured by atomic force microscopy

Tapping-mode atomic force microscopy was used to study the time-dependent changes in the structure of fibrinogen under aqueous conditions following adsorption on two model surfaces: hydrophobic graphite and hydrophilic mica. Fibrinogen was observed in the characteristic trinodular form, and the dimensions of the adsorbed molecules were consistent with previously reported values for these surfaces. On the basis of the differences in the relative heights of the D and the E domains, four orientation states were observed for fibrinogen adsorbed on both the surfaces. On graphite, the initial asymmetric orientation states disappeared with spreading over time. Some small lateral movements of the adsorbed proteins were observed on mica during repeated scanning, whereas no such movement was observed on graphite, indicating strong adhesion of fibrinogen to a hydrophobic surface. Spreading kinetics of fibrinogen on the two surfaces was determined by measuring the heights of the D and E domains over a time period of approximately 2 h. On graphite, the heights of both the D and E domains decreased with time to a lower plateau value of 1.0 nm. On mica, the heights of both the D and E domains showed an increase, rising to an upper plateau value of approximately 2.1 nm. The spreading of the D and E domains on graphite was analyzed using an 'exponential-decay-of-height' model. A spreading rate constant of approximately 4.7 x 10(-4) s(-1) was observed for the whole fibrinogen molecule adsorbed on graphite, corresponding to a free energy of unfolding of approximately 37 kT. Extrapolation of the exponential curve in the model to t = 0 yielded values of 2.3 and 2.2 nm for the heights of the D and the E domains at the time of contact with the hydrophobic graphite substrate, significantly less than their free solution diameters. A two-step spreading model is proposed to explain this observation.

Unification of the Folding Mechanisms of Non-two-state and Two-state Proteins

We have collected the kinetic folding data for non-two-state and two-state globular proteins reported in the literature, and investigated the relationships between the folding kinetics and the native three-dimensional structure of these proteins. The rate constants of formation of both the intermediate and the native state of non-two-state folders were found to be significantly correlated with protein chain length and native backbone topology, which is represented by the absolute contact order and sequence-distant native pairs. The folding rate of two-state folders, which is known to be correlated with the native backbone topology, apparently does not correlate significantly with protein chain length. On the basis of a comparison of the folding rates of the non-two-state and two-state folders, it was found that they are similarly dependent on the parameters that reflect the native backbone topology. This suggests that the mechanisms behind non-two-state and two-state folding are essentially identical. The present results lead us to propose a unified mechanism of protein folding, in which folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the intermediate. Apparently, two-state folding is merely a simplified version of hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.

Outlining folding nuclei in globular proteins

Our theoretical approach for prediction of folding/unfolding nuclei in three-dimensional protein structures is based on a search for free energy saddle points on networks of protein unfolding pathways. Under some approximations, this search is performed rapidly by dynamic programming and results in prediction of Phi values, which can be compared with those found experimentally. In this study, we optimize some details of the model (specifically, hydrogen atoms are taken into account in addition to heavy atoms), and compare the theoretically obtained and experimental Phi values (which characterize involvement of residues in folding nuclei) for all 17 proteins, where Phi values are now known for many residues. We show that the model provides good Phi value predictions for proteins whose structures have been determined by X-ray analysis (the average correlation coefficient is 0.65), with a more limited success for proteins whose structures have been determined by NMR techniques only (the average correlation coefficient is 0.34), and that the transition state free energies computed from the same model are in a good anticorrelation with logarithms of experimentally measured folding rates at mid-transition (the correlation coefficient is -0.73).

Contact order revisited: influence of protein size on the folding rate

Guided by the recent success of empirical model predicting the folding rates of small two-state folding proteins from the relative contact order (CO) of their native structures, by a theoretical model of protein folding that predicts that logarithm of the folding rate decreases with the protein chain length L as L(2/3), and by the finding that the folding rates of multistate folding proteins strongly correlate with their sizes and have very bad correlation with CO, we reexamined the dependence of folding rate on CO and L in attempt to find a structural parameter that determines folding rates for the totality of proteins. We show that the Abs_CO = CO x L, is able to predict rather accurately folding rates for both two-state and multistate folding proteins, as well as short peptides, and that this Abs_CO scales with the protein chain length as L(0.70 +/- 0.07) for the totality of studied single-domain proteins and peptides.