The stability and dynamics of computationally designed proteins.

Protein stability, dynamics and function are intricately linked. Accordingly, protein designers leverage dynamics in their designs and gain insight to their successes and failures by analyzing their proteins' dynamics. Molecular dynamics (MD) simulations are a powerful computational tool for quantifying both local and global protein dynamics. This review highlights studies where MD simulations were applied to characterize the stability and dynamics of designed proteins and where dynamics were incorporated into computational protein design. First, we discuss the structural basis underlying the extreme stability and thermostability frequently observed in computationally designed proteins. Next, we discuss examples of designed proteins, where dynamics were not explicitly accounted for in the design process, whose coordinated motions or active site dynamics, as observed by MD simulation, enhanced or detracted from their function. Many protein functions depend on sizeable or subtle conformational changes, so we finally discuss the computational design of proteins to perform a specific function that requires consideration of motion by multi-state design.

Thermostabilization mechanisms in thermophilic versus mesophilic three-helix bundle proteins.

The engineered three-helix bundle, UVF, is thermostabilized entropically due to heightened, native-state dynamics. However, it is unclear whether this thermostabilization strategy is observed in natural proteins from thermophiles. We performed all-atom, explicit solvent molecular dynamics simulations of two three-helix bundles from thermophilic H. butylicus (2lvsN and 2lvsC) and compared their dynamics to a mesophilic three-helix bundle, the Engrailed homeodomain (EnHD). Like UVF, 2lvsC had heightened native dynamics, which it maintained without unfolding at 100°C. Shortening and rigidification of loops in 2lvsN and 2lvsC and increased surface hydrogen bonds in 2lvsN were observed, as is common in thermophilic proteins. A buried disulfide and salt bridge in 2lvsN and 2lvsC, respectively, provided some stabilization, and addition of a homologous disulfide bond in EnHD slowed unfolding. The transferability and commonality of stabilization strategies among members of the three-helix bundle fold suggest that these strategies may be general and deployable in designing thermostable proteins.

The Role of Electrostatics and Folding Kinetics on the Thermostability of Homologous Cold Shock Proteins.

Understanding which aspects contribute to the thermostability of proteins is a challenge that has persisted for decades, and it is of great relevance for protein engineering. Several types of interactions can influence the thermostability of a protein. Among them, the electrostatic interactions have been a target of particular attention. Aiming to explore how this type of interaction can affect protein thermostability, this paper investigated four homologous cold shock proteins from psychrophilic, mesophilic, thermophilic, and hyperthermophilic organisms using a set of theoretical methodologies. It is well-known that electrostatics as well as hydrophobicity are key-elements for the stabilization of these proteins. Therefore, both interactions were initially analyzed in the native structure of each protein. Electrostatic interactions present in the native structures were calculated with the Tanford-Kirkwood model with solvent accessibility, and the amount of hydrophobic surface area buried upon folding was estimated by measuring both folded and extended structures. On the basis of Energy Landscape Theory, the local frustration and the simplified alpha-carbon structure-based model were modeled with a Debye-Hückel potential to take into account the electrostatics and the effects of an implicit solvent. Thermodynamic data for the structure-based model simulations were collected and analyzed using the Weighted Histogram Analysis and Stochastic Diffusion methods. Kinetic quantities including folding times, transition path times, folding routes, and Φ values were also obtained. As a result, we found that the methods are able to qualitatively infer that electrostatic interactions play an important role on the stabilization of the most stable thermophilic cold shock proteins, showing agreement with the experimental data.

A Dynamic Hydrophobic Core and Surface Salt Bridges Thermostabilize a Designed Three-Helix Bundle.

Thermostable proteins are advantageous in industrial applications, as pharmaceuticals or biosensors, and as templates for directed evolution. As protein-design methodologies improve, bioengineers are able to design proteins to perform a desired function. Although many rationally designed proteins end up being thermostable, how to intentionally design de novo, thermostable proteins is less clear. UVF is a de novo-designed protein based on the backbone structure of the Engrailed homeodomain (EnHD) and is highly thermostable (Tm > 99°C vs. 52°C for EnHD). Although most proteins generally have polar amino acids on their surfaces and hydrophobic amino acids buried in their cores, protein engineers followed this rule exactly when designing UVF. To investigate the contributions of the fully hydrophobic core versus the fully polar surface to UVF's thermostability, we built two hybrid, chimeric proteins combining the sets of buried and surface residues from UVF and EnHD. Here, we determined a structural, dynamic, and thermodynamic explanation for UVF's thermostability by performing 4 µs of all-atom, explicit-solvent molecular dynamics simulations at 25 and 100°C, Tanford-Kirkwood solvent accessibility Monte Carlo electrostatic calculations, and a thermodynamic analysis of 40 temperature runs by the weighted-histogram analysis method of heavy-atom, structure-based models of UVF, EnHD, and both chimeric proteins. Our models showed that UVF was highly dynamic because of its fully hydrophobic core, leading to a smaller loss of entropy upon folding. The charged residues on its surface made favorable electrostatic interactions that contributed enthalpically to its thermostability. In the chimeric proteins, both the hydrophobic core and charged surface independently imparted thermostability.

Molecular dynamics simulations suggest stabilizing mutations in a de novo designed α/ß protein.

Designing functional proteins that can withstand extreme heat is beneficial for industrial and protein therapeutic applications. Thus, elucidating the atomic-level determinants of thermostability is a major interest for rational protein design. To that end, we compared the structure and dynamics of a set of previously designed, thermostable proteins based on the activation domain of human procarboxypeptidase A2 (AYEwt). The mutations in these designed proteins were intended to increase hydrophobic core packing and inter-secondary-structure interactions. To evaluate whether these design strategies were successfully deployed, we performed all-atom, explicit-solvent molecular dynamics (MD) simulations of AYEwt and three designed variants at both 25 and 100°C. Our MD simulations agreed with the relative experimental stabilities of the designs based on their secondary structure content, Cα root-mean-square deviation/fluctuation, and buried-residue solvent accessible surface area. Using a contact analysis, we found that the designs stabilize inter-secondary structure interactions and buried hydrophobic surface area, as intended. Based on our analysis, we designed three additional variants to test the role of helix stabilization, core packing, and a Phe â†’ Met mutation on thermostability. We performed the additional MD simulations and analysis on these variants, and these data supported our predictions.

A Carboxylate to Amide Substitution That Switches Protein Folds.

Metamorphic proteins are biomolecules prone to adopting alternative conformations. Because of this feature, they represent ideal systems to investigate the general rules allowing primary structure to dictate protein topology. A comparative molecular dynamics study was performed on the denatured states of two proteins, sharing nearly identical amino-acid sequences (88 %) but different topologies, namely an all-α-helical bundle protein named GA 88 and an α+ß-protein named GB 88. The analysis allowed successful design of and experimental validation of a site-directed mutant that promotes, at least in part, the switch in folding from GB 88 to GA 88. The mutated position, in which a glutamic acid was replaced by a glutamine, does not make any intramolecular interactions in the native state of GA 88, such that its stabilization can be explained by considering the effects on the denatured state. The results represent a direct demonstration of the role of the denatured state in sculpting native structure.

Structural and functional consequences of cardiac troponin C L57Q and I61Q Ca(2+)-desensitizing variants.

Two cTnC variants, L57Q and I61Q, both of which are located on helix C within the N domain of cTnC, were originally reported in the skeletal muscle system [Tikunova, Davis, J. Biol. Chem. 279 (2004) 35341-35352], as the analogous L58Q and I62Q sTnC, and demonstrated a decreased Ca(2+) binding affinity. Here, we provide detailed characterization of structure-function relationships for these two cTnC variants, to determine if they behave differently in the cardiac system and as a framework for determining similarities and differences with other cTnC mutations that have been associated with DCM. We have used an integrative approach to study the structure and function of these cTnC variants both in solution and in silico, to understand how the L57Q and I61Q mutations influence Ca(2+) binding at site II, the subsequent effects on the interaction with cTnI, and the structural changes which are associated with these changes. Steady-state and stopped flow fluorescence spectroscopy confirmed that a decrease in Ca(2+) affinity for recombinant cTnC and cTn complexes containing the L57Q or I61Q variants. The L57Q variant was intermediate between WT and I61Q cTnC and also did not significantly alter cTnC-cTnI interaction in the absence of Ca(2+), but did decrease the interaction in the presence of Ca(2+). In contrast, I61Q decreased the cTnC-cTnI interaction in both the absence and presence of Ca(2+). This difference in the absence of Ca(2+) suggests a greater structural change in cNTnC may occur with the I61Q mutation than the L57Q mutation. MD simulations revealed that the decreased Ca(2+) binding induced by I61Q may result from destabilization of the Ca(2+) binding site through interruption of intra-molecular interactions when residue 61 forms new hydrogen bonds with G70 on the Ca(2+) binding loop. The experimentally observed interruption of the cTnC-cTnI interaction caused by L57Q or I61Q is due to the disruption of key hydrophobic interactions between helices B and C in cNTnC. This study provides a molecular basis of how single mutations in the C helix of cTnC can reduce Ca(2+) binding affinity and cTnC-cTnI interaction, which may provide useful insights for a better understanding of cardiomyopathies and future gene-based therapies.

Promiscuous contacts and heightened dynamics increase thermostability in an engineered variant of the engrailed homeodomain.

A thermostabilized variant (UVF) of the engrailed homeodomain (EnHD) was previously engineered by Mayo and co-workers. The melting temperature of the non-natural, designed protein is 50°C higher than the natural wild-type protein (>99 vs. 52°C), and the two proteins share 22% sequence identity. We have performed extensive (1 µs) all-atom, explicit solvent molecular dynamics simulations of the wild-type and engineered proteins to investigate their structural and dynamic properties at room temperature and at 100°C. Our simulations are in good agreement with nuclear magnetic resonance data available for the two proteins [nuclear Overhauser effect crosspeaks (NOEs), J-coupling constants and order parameters for EnHD; and NOEs for UVF], showing that we reproduce the backbone dynamics and side chain packing in the native state of both proteins. UVF was more dynamic at room temperature than EnHD, with respect to both its backbone and side chain motion. When the temperature was raised, the thermostable protein maintained this mobility while retaining its native conformation. EnHD, on the other hand, was unable to maintain its more rigid native structure at higher temperature and began to unfold. Heightened protein dynamics leading to promiscuous and dynamically interchangeable amino acid contacts makes UVF more tolerant to increasing temperature, providing a molecular explanation for heightened thermostability of this protein.

Multimolecule test-tube simulations of protein unfolding and aggregation.

Molecular dynamics simulations of protein folding or unfolding, unlike most in vitro experimental methods, are performed on a single molecule. The effects of neighboring molecules on the unfolding/folding pathway are largely ignored experimentally and simply not modeled computationally. Here, we present two all-atom, explicit solvent molecular dynamics simulations of 32 copies of the Engrailed homeodomain (EnHD), an ultrafast-folding and -unfolding protein for which the folding/unfolding pathway is well-characterized. These multimolecule simulations, in comparison with single-molecule simulations and experimental data, show that intermolecular interactions have little effect on the folding/unfolding pathway. EnHD unfolded by the same mechanism whether it was simulated in only water or also in the presence of other EnHD molecules. It populated the same native state, transition state, and folding intermediate in both simulation systems, and was in good agreement with experimental data available for each of the three states. Unfolding was slowed slightly by interactions with neighboring proteins, which were mostly hydrophobic in nature and ultimately caused the proteins to aggregate. Protein-water hydrogen bonds were also replaced with protein-protein hydrogen bonds, additionally contributing to aggregation. Despite the increase in protein-protein interactions, the protein aggregates formed in simulation did not do so at the total exclusion of water. These simulations support the use of single-molecule techniques to study protein unfolding and also provide insight into the types of interactions that occur as proteins aggregate at high temperature at an atomic level.

Structural and functional consequences of the cardiac troponin C L48Q Ca(2+)-sensitizing mutation.

Calcium binding to the regulatory domain of cardiac troponin C (cNTnC) causes a conformational change that exposes a hydrophobic surface to which troponin I (cTnI) binds, prompting a series of protein-protein interactions that culminate in muscle contraction. A number of cTnC variants that alter the Ca(2+) sensitivity of the thin filament have been linked to disease. Tikunova and Davis engineered a series of cNTnC mutations that altered Ca(2+) binding properties and studied the effects on the Ca(2+) sensitivity of the thin filament and contraction [Tikunova, S. B., and Davis, J. P. (2004) J. Biol. Chem. 279, 35341-35352]. One of the mutations they engineered, the L48Q variant, resulted in a pronounced increase in the cNTnC Ca(2+) binding affinity and Ca(2+) sensitivity of cardiac muscle force development. In this work, we sought structural and mechanistic explanations for the increased Ca(2+) sensitivity of contraction for the L48Q cNTnC variant, using an array of biophysical techniques. We found that the L48Q mutation enhanced binding of both Ca(2+) and cTnI to cTnC. Nuclear magnetic resonance chemical shift and relaxation data provided evidence that the cNTnC hydrophobic core is more exposed with the L48Q variant. Molecular dynamics simulations suggest that the mutation disrupts a network of crucial hydrophobic interactions so that the closed form of cNTnC is destabilized. The findings emphasize the importance of cNTnC's conformation in the regulation of contraction and suggest that mutations in cNTnC that alter myofilament Ca(2+) sensitivity can do so by modulating Ca(2+) and cTnI binding.

The denatured state dictates the topology of two proteins with almost identical sequence but different native structure and function.

The protein folding problem is often studied by comparing the mechanisms of proteins sharing the same structure but different sequence. The recent design of the two proteins G(A)88 and G(B)88, displaying different structures and functions while sharing 88% sequence identity (49 out of 56 amino acids), allows the unique opportunity for a complementary approach. At which stage of its folding pathway does a protein commit to a given topology? Which residues are crucial in directing folding mechanisms to a given structure? By using a combination of biophysical and computational techniques, we have characterized the folding of both G(A)88 and G(B)88. We show that, contrary to expectation, G(B)88, characterized by a native α+ß fold, displays in the denatured state a content of native-like helical structure greater than G(A)88, which is all-α in its native state. Both experiments and simulations indicate that such residual structure may be tuned by changing pH. Thus, despite the high sequence identity, the folding pathways for these two proteins appear to diverge as early as in the denatured state. Our results suggest a mechanism whereby protein topology is committed very early along the folding pathway, being imprinted in the residual structure of the denatured state.

Refolding the engrailed homeodomain: structural basis for the accumulation of a folding intermediate.

The ultrafast folding pathway of the engrailed homeodomain has been exceptionally well characterized by experiment and simulation. Helices II and III of the three-helix bundle protein form the native helix-turn-helix motif as an on-pathway intermediate within a few microseconds. The slow step is then the proper docking of the helices in approximately 15 mus. However, there is still the unexplained puzzle of why helix docking is relatively slow, which is part of the more general question as to why rearrangements of intermediates occur slowly. To address this problem, we performed 46 all-atom molecular dynamics refolding simulations in explicit water, for a total of 15 micros of simulation time. The simulations started from an intermediate state structure that was generated in an unfolding simulation at 498 K and was then quenched to folding-permissive temperatures. The protein refolded successfully in only one of the 46 simulations, and in that case the refolding pathway mirrored the unfolding pathway at high temperature. In the 45 simulations in which the protein did not fully fold, nonnative salt bridges trapped the protein, which explains why the protein folds relatively slowly from the intermediate state.

Microscopic reversibility of protein folding in molecular dynamics simulations of the engrailed homeodomain.

The principle of microscopic reversibility states that at equilibrium the number of molecules entering a state by a given path must equal those exiting the state via the same path under identical conditions or, in structural terms, that the conformations along the two pathways are the same. There has been some indirect evidence indicating that protein folding is such a process, but there have been few conclusive findings. In this study, we performed molecular dynamics simulations of an ultrafast unfolding and folding protein at its melting temperature to observe, on an atom-by-atom basis, the pathways the protein followed as it unfolded and folded within a continuous trajectory. In a total of 0.67 micros of simulation in water, we found six transient denaturing events near the melting temperature (323 and 330 K) and an additional refolding event following a previously identified unfolding event at a high temperature (373 K). In each case, unfolding and refolding transition state ensembles were identified, and they agreed well with experiment on the basis of a comparison of S and Phi values. On the basis of several structural properties, these 13 transition state ensembles agreed very well with each other and with four previously identified transition states from high-temperature denaturing simulations. Thus, not only were the unfolding and refolding transition states part of the same ensemble, but in five of the seven cases, the pathway the protein took as it unfolded was nearly identical to the subsequent refolding pathway. These events provide compelling evidence that protein folding is a microscopically reversible process. In the other two cases, the folding and unfolding transition states were remarkably similar to each other but the paths deviated.

Structure/function analysis of the interaction of phosphatidylinositol 4,5-bisphosphate with actin-capping protein: implications for how capping protein binds the actin filament.

The heterodimeric actin-capping protein (CP) can be inhibited by polyphosphoinositides, which may be important for actin polymerization at membranes in cells. Here, we have identified a conserved set of basic residues on the surface of CP that are important for the interaction with phosphatidylinositol 4,5-bisphosphate (PIP(2)). Computational docking studies predicted the identity of residues involved in this interaction, and functional and physical assays with site-directed mutants of CP confirmed the prediction. The PIP(2) binding site overlaps with the more important of the two known actin-binding sites of CP. Correspondingly, we observed that loss of PIP(2) binding correlated with loss of actin binding among the mutants. Using TIRF (total internal reflection fluorescence) microscopy, we observed that PIP(2) rapidly converted capped actin filaments to a growing state, consistent with uncapping. Together, these results extend our understanding of how CP binds to the barbed end of the actin filament, and they support the idea that CP can "wobble" when bound to the barbed end solely by the C-terminal "tentacle" of its beta-subunit.