Energetics of aliphatic deletions in protein cores

Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Despite advances in artificial intelligence methods, protein folding remains in many ways an enigma to be solved. Accurate computation of protein folding energetics could help drive fields such as protein and drug design and genetic interpretation. However, the challenge of calculating the state functions governing protein folding from first-principles remains unaddressed. We present here a simple approach that allows us to accurately calculate the energetics of protein folding. It is based on computing the energy of the folded and unfolded states at different temperatures using molecular dynamics simulations. From this, two essential quantities (ΔH and ΔCp) are obtained and used to calculate the conformational stability of the protein (ΔG). With this approach, we have successfully calculated the energetics of two- and three-state proteins, representatives of the major structural classes, as well as small stability differences (ΔΔG) due to changes in solution conditions or variations in an amino acid residue.

Protposer: the web server that readily proposes protein stabilizing mutations with high PPV

Protein stability is a requisite for most biotechnological and medical applications of proteins. As natural proteins tend to suffer from a low conformational stability ex vivo, great efforts have been devoted toward increasing their stability through rational design and engineering of appropriate mutations. Unfortunately, even the best currently used predictors fail to compute the stability of protein variants with sufficient accuracy and their usefulness as tools to guide the rational stabilisation of proteins is limited. We present here Protposer, a protein stabilising tool based on a different approach. Instead of quantifying changes in stability, Protposer uses structure- and sequence-based screening modules to nominate candidate mutations for subsequent evaluation by a logistic regression model, carefully trained to avoid overfitting. Thus, Protposer analyses PDB files in search for stabilization opportunities and provides a ranked list of promising mutations with their estimated success rates (eSR), their probabilities of being stabilising by at least 0.5 kcal/mol. The agreement between eSRs and actual positive predictive values (PPV) on external datasets of mutations is excellent. When Protposer is used with its Optimal kappa selection threshold, its PPV is above 0.7. Even with less stringent thresholds, Protposer largely outperforms FoldX, Rosetta and PoPMusiC. Indicating the PDB file of the protein suffices to obtain a ranked list of mutations, their eSRs and hints on the likely source of the stabilization expected. Protposer is a distinct, straightforward and highly successful tool to design protein stabilising mutations, and it is freely available for academic use at http://webapps.bifi.es/the-protposer.

Net Evolutionary Loss of Residue Polarity in Drosophilid Protein Cores Indicates Ongoing Optimization of Amino Acid Composition

Amino acid frequencies in proteins may not be at equilibrium. We consider two possible explanations for the nonzero net residue fluxes in drosophilid proteins. First, protein interiors may have a suboptimal residue composition and be under a selective pressure favoring stability, that is, leading to the loss of polar (and the gain of large) amino acids. One would then expect stronger net fluxes on the protein interior than at the exposed sites. Alternatively, if most of the polarity loss occurs at the exposed sites and the selective constraint on amino acid composition at such sites decreases over time, net loss of polarity may be neutral and caused by disproportionally high occurrence of polar residues at exposed, least constrained sites. We estimated net evolutionary fluxes of residue polarity and volume at sites with different solvent accessibility in conserved protein families from 12 species of Drosophila. Net loss of polarity, miniscule in magnitude, but consistent across all lineages, occurred at all sites except the most exposed ones, where net flux of polarity was close to zero or, in membrane proteins, even positive. At the intermediate solvent accessibility the net fluxes of polarity and volume were similar to neutral predictions, whereas much of the polarity loss not attributable to neutral expectations occurred at the buried sites. These observations are consistent with the hypothesis that residue composition in many proteins is structurally suboptimal and continues to evolve toward lower polarity in the protein interior, in particular in proteins with intracellular localization. The magnitude of polarity and volume changes was independent from the protein’s evolutionary age, indicating that the approach to equilibrium has been slow or that no such single equilibrium exists.

Contribution of buried distal amino acid residues in horse liver alcohol dehydrogenase to structure and catalysis

The dynamics of enzyme catalysis range from the slow time scale (∼ms) for substrate binding and conformational changes to the fast time (∼ps) scale for reorganization of substrates in the chemical step. The contribution of global dynamics to catalysis by alcohol dehydrogenase was tested by substituting five different, conserved amino acid residues that are distal from the active site and located in the hinge region for the conformational change or in hydrophobic clusters. X-ray crystallography shows that the structures for the G173A, V197I, I220 (V, L, or F), V222I, and F322L enzymes complexed with NAD⁺ and an analogue of benzyl alcohol are almost identical, except for small perturbations at the sites of substitution. The enzymes have very similar kinetic constants for the oxidation of benzyl alcohol and reduction of benzaldehyde as compared to the wild-type enzyme, and the rates of conformational changes are not altered. Less conservative substitutions of these amino acid residues, such as G173(V, E, K, or R), V197(G, S, or T), I220(G, S, T, or N), and V222(G, S, or T) produced unstable or poorly expressed proteins, indicating that the residues are critical for global stability. The enzyme scaffold accommodates conservative substitutions of distal residues, and there is no evidence that fast, global dynamics significantly affect the rate constants for hydride transfers. In contrast, other studies show that proximal residues significantly participate in catalysis.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

Effects of cavities at the nicotinamide binding site of liver alcohol dehydrogenase on structure, dynamics and catalysis

A role for protein dynamics in enzymatic catalysis of hydrogen transfer has received substantial scientific support, but the connections between protein structure and catalysis remain to be established. Valine residues 203 and 207 are at the binding site for the nicotinamide ring of the coenzyme in liver alcohol dehydrogenase and have been suggested to facilitate catalysis with “protein-promoting vibrations” (PPV). We find that the V207A substitution has small effects on steady-state kinetic constants and the rate of hydrogen transfer; the introduced cavity is empty and is tolerated with minimal effects on structure (determined at 1.2 Å for the complex with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol). Thus, no evidence is found to support a role for Val-207 in the dynamics of catalysis. The protein structures and ligand geometries (including donor–acceptor distances) in the V203A enzyme complexed with NAD⁺ and 2,3,4,5,6-pentafluorobenzyl alcohol or 2,2,2-trifluoroethanol (determined at 1.1 Å) are very similar to those for the wild-type enzyme, except that the introduced cavity accommodates a new water molecule that contacts the nicotinamide ring. The structures of the V203A enzyme complexes suggest, in contrast to previous studies, that the diminished tunneling and decreased rate of hydride transfer (16-fold, relative to that of the wild-type enzyme) are not due to differences in ground-state ligand geometries. The V203A substitution may alter the PPV and the reorganization energy for hydrogen transfer, but the protein scaffold and equilibrium thermal motions within the Michaelis complex may be more significant for enzyme catalysis.

Structural and energetic basis of carbohydrate-aromatic packing interactions in proteins

Carbohydrate-aromatic interactions mediate many biological processes. However, the structure-energy relationships underpinning direct carbohydrate-aromatic packing interactions in aqueous solution have been difficult to assess experimentally and remain elusive. Here, we determine the structures and folding energetics of chemically synthesized glycoproteins to quantify the contributions of the hydrophobic effect and CH-π interactions to carbohydrate-aromatic packing interactions in proteins. We find that the hydrophobic effect contributes significantly to protein-carbohydrate interactions. Interactions between carbohydrates and aromatic amino acid side chains, however, are supplemented by CH-π interactions. The strengths of experimentally determined carbohydrate CH-π interactions do not correlate with the electrostatic properties of the involved aromatic residues, suggesting that the electrostatic component of CH-π interactions in aqueous solution is small. Thus, tight binding of carbohydrates and aromatic residues is driven by the hydrophobic effect and CH-π interactions featuring a dominating dispersive component.

Interplay between drying and stability of a TIM barrel protein: a combined simulation-experimental study

Recent molecular dynamics simulations have suggested important roles for nanoscale dewetting in the stability, function, and folding dynamics of proteins. Using a synergistic simulation-experimental approach on the αTS TIM barrel protein, we validated this hypothesis by revealing the occurrence of drying inside hydrophobic amino acid clusters and its manifestation in experimental measures of protein stability and structure. Cavities created within three clusters of branched aliphatic amino acids [isoleucine, leucine, and valine (ILV) clusters] were found to experience strong water density fluctuations or intermittent dewetting transitions in simulations. Individually substituting 10 residues in the large ILV cluster at the N-terminus with less hydrophobic alanines showed a weakening or diminishing effect on dewetting that depended on the site of the mutation. Our simulations also demonstrated that replacement of buried leucines with isosteric, polar asparagines enhanced the wetting of the N- and C-terminal clusters. The experimental results on the stability, secondary structure, and compactness of the native and intermediate states for the asparagine variants are consistent with the preferential drying of the large N-terminal cluster in the intermediate. By contrast, the region encompassing the small C-terminal cluster experiences only partial drying in the intermediate, and its structure and stability are unaffected by the asparagine substitution. Surprisingly, the structural distortions required to accommodate the replacement of leucine by asparagine in the N-terminal cluster revealed the existence of alternative stable folds in the native basin. This combined simulation-experimental study demonstrates the critical role of drying within hydrophobic ILV clusters in the folding and stability of the αTS TIM barrel.

Defining the nature of thermal intermediate in 3 state folding proteins: apoflavodoxin, a study case

The early stages of the thermal unfolding of apoflavodoxin have been determined by using atomistic multi microsecond-scale molecular dynamics (MD) simulations complemented with a variety of experimental techniques. Results strongly suggest that the intermediate is reached very early in the thermal unfolding process and that it has the properties of an “activated” form of the native state, where thermal fluctuations in the loops break loop-loop contacts. The unrestrained loops gain then kinetic energy corrupting short secondary structure elements without corrupting the core of the protein. The MD-derived ensembles agree with experimental observables and draw a picture of the intermediate state inconsistent with a well-defined structure and characteristic of a typical partially disordered protein. Our results allow us to speculate that proteins with a well packed core connected by long loops might behave as partially disordered proteins under native conditions, or alternatively behave as three state folders. Small details in the sequence, easily tunable by evolution, can yield to one or the other type of proteins.

A simplistic view of protein structure tends to emphasize the opposition between the native state and the denatured ensemble of unfolded conformations. In addition to these extreme conformations, proteins subjected to a variety of perturbations often populate alternative partly unfolded conformations, some of which are close in energy to the native state and, accordingly, can be populated under native or quasi-native conditions. There is increasing evidence that these “perturbed” conformations participate in protein function or, in some cases, are related to the outcome of folding diseases. We have used the “state of the art” molecular dynamics combined with a variety of experimental techniques to characterize for the first time, to our knowledge, the thermal intermediate of a three-state folding protein (apoflavodoxin). Based on our results we have been able to suggest a general mechanism of thermal unfolding in complex proteins and to determine interesting links between thermal intermediates and partially unfolded 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.

Novel modulation factor quantifies the role of water molecules in protein interactions

Water molecules decrease the potential of mean force of a hydrogen bond (H-bond), as well as modulate (de)solvation forces, but exactly how much has not been easy to determine. Crystallographic water molecules provide snapshots of optimal solutions for the role of solvent in protein interactions, information that is often ignored by implicit solvent models. Motivated by high-resolution crystal structures, we describe a simple quantitative approach to explicitly incorporate the role of molecular water in protein interactions. Applications to protein-DNA interactions show that the accuracy of binding free-energy estimates improves significantly if a distinction is made between H-bonds that are desolvated (or only contact crystal waters), solvated by mobile waters trapped at the binding interface, or partially solvated through connections to bulk water. These different environments are modeled by a unique "water" scaling factor that decreases or increases the strength of hydrogen bonds depending on whether water contacts the acceptor or donor atoms or the bond is fully desolvated, respectively. Our empirical energies are fully consistent with mobile water molecules having a strong polarization effect in direct intermolecular interactions.

Structural analysis of an equilibrium folding intermediate in the apoflavodoxin native ensemble by small-angle X-ray scattering

Intermediate conformations are crucial to our understanding of how proteins fold into their native structures and become functional. Conventional spectroscopic measurements of thermal denaturation transitions allow the detection of equilibrium intermediates but often provide little structural detail; thus, application of more informative techniques is required. Here we used small-angle X-ray scattering (SAXS) to study the thermal denaturation of four variants of Anabaena PCC 7119 flavodoxin, including the wild-type apo and holo forms, and two mutants, E20K/E72K and F98N. Denaturation was monitored from changes in SAXS descriptors. Although the starting and final points of the denaturation were similar for the flavodoxin variants tested, substantial differences in the unfolding pathway were apparent between them. In agreement with calorimetric data, analysis of the SAXS data sets indicated a three-state unfolding equilibrium for wild-type apoflavodoxin, a two-state equilibrium for the F98N mutant, and increased thermostability of the E20K/E72K mutant and holoflavodoxin. Although the apoflavodoxin intermediate consistently appeared mixed with significant amounts of either native or unfolded conformations, its SAXS profile was derived from the deconvolution of the temperature-dependent SAXS data set. The apoflavodoxin thermal intermediate was structurally close to the native state but less compact, thereby indicating incipient unfolding. The residues that foster denaturation were explored by an ensemble of equilibrium ϕ-value restrained molecular dynamics. These simulations pointed to residues located in the cofactor and partner-protein recognition regions as the initial sites of denaturation and suggest a conformational adaptation as the mechanism of action in apoflavodoxin.

Optimization of the Gβ1 Domain by Computational Design and by in Vitro Evolution: Structural and Energetic Basis of Stabilization

Computational design and in vitro evolution are major strategies for stabilizing proteins. For the four critical positions 16, 18, 25, and 29 of the B domain of the streptococcal protein G (Gbeta1), they identified the same optimal residues at positions 16 and 25, but not at 18 and 29. Here we analyzed the energetic contributions of the residues from these two approaches by single and double mutant analyses and determined crystal structures for a variant from the calculation (I16/L18/E25/K29) and from the selection (I16/I18/E25/F29). The structural analysis explains the observed differences in stabilization. Residues 16, 18, and 29 line an invagination, which results from a packing defect between the helix and the beta-sheet of Gbeta1. In all stabilized variants, residues with larger side-chains occur at these positions and packing is improved. In the selected variant, packing is better optimized than in the computed variant. Such differences in side-chain packing strongly affect stability but are difficult to evaluate by computation.

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.

Protein-protein interactions and disease: use of S. cerevisiae as a model system

Disease-causing mutations are increasingly being studied to see if they cause the loss or gain of protein-protein interactions. Because the interaction network of humans is poorly understood and difficult to investigate, here we propose the use of Saccharomyces cerevisiae as a model system for understanding the impact of disease-causing mutations on protein-protein interactions. Alignments of human disease-associated proteins and 379 yeast orthologs showed that 124 of these proteins have >40% sequence identity, with some orthologs having up to 89% identity. A total of 1826 amino acid mutations associated with human disease were found to map to invariant amino acids in yeast. These mutations were proportionately more likely to be non-conservative than non-disease associated polymorphisms for the same proteins (p=0.016). Importantly, 73 of the mutations mapped to protein-protein interaction domains, implying a direct link between mutation and changes in protein interactivity. In the manuscript, all alignment information and tables that map mutations and diseases to yeast orthologs are given. This will help researchers experimentally test the impact of mutations on protein-protein interactions in S. cerevisiae and, by homology, explore the role of such mutations in the genesis of human disease.