Coarse-grained strategy for modeling protein stability in concentrated solutions. III: directional protein interactions

We extend our coarse-grained modeling strategy described in parts I and II of this investigation to account for nonuniform spatial distributions of hydrophobic residues on the solvent-exposed surfaces of native proteins. Within this framework, we explore how patchy surfaces can influence the solvent-mediated protein-protein interactions, and the unfolding and self-assembly behaviors of proteins in solution. In particular, we compare the equilibrium unfolding and self-assembly trends for three model proteins that share the same overall sequence hydrophobicity, but exhibit folded configurations with different solvent-exposed native-state surface morphologies. Our model provides new insights into how directional interactions can affect native-state protein stability in solution. We find that strongly-directional attractions between native molecules with patchy surfaces can help stabilize the folded conformation through the formation of self-assembled clusters. In contrast, native proteins with more uniform surfaces are destabilized by protein-protein attractions involving the denatured state. Finally, we discuss how the simulation results provide insights into the experimental solution behaviors of several proteins that display directional interactions in their native states.     

Coarse-grained strategy for modeling protein stability in concentrated solutions. II: phase behavior

We use highly efficient transition-matrix Monte Carlo simulations to determine equilibrium unfolding curves and fluid phase boundaries for solutions of coarse-grained globular proteins. The model we analyze derives the intrinsic stability of the native state and protein-protein interactions from basic information about protein sequence using heteropolymer collapse theory. It predicts that solutions of low hydrophobicity proteins generally exhibit a single liquid phase near their midpoint temperatures for unfolding, while solutions of proteins with high sequence hydrophobicity display the type of temperature-inverted, liquid-liquid transition associated with aggregation processes of proteins and other amphiphilic molecules. The phase transition occurring in solutions of the most hydrophobic protein we study extends below the unfolding curve, creating an immiscibility gap between a dilute, mostly native phase and a concentrated, mostly denatured phase. The results are qualitatively consistent with the solution behavior of hemoglobin (HbA) and its sickle variant (HbS), and they suggest that a liquid-liquid transition resulting in significant protein denaturation should generally be expected on the phase diagram of high-hydrophobicity protein solutions. The concentration fluctuations associated with this transition could be a driving force for the nonnative aggregation that can occur below the midpoint temperature.     

Protein-DNA interactions: structural, thermodynamic and clustering patterns of conserved residues in DNA-binding proteins

Amino acid residues, which play important roles in protein function, are often conserved. Here, we analyze thermodynamic and structural data of protein-DNA interactions to explore a relationship between free energy, sequence conservation and structural cooperativity. We observe that the most stabilizing residues or putative hotspots are those which occur as clusters of conserved residues. The higher packing density of the clusters and available experimental thermodynamic data of mutations suggest cooperativity between conserved residues in the clusters. Conserved singlets contribute to the stability of protein-DNA complexes to a lesser extent. We also analyze structural features of conserved residues and their clusters and examine their role in identifying DNA-binding sites. We show that about half of the observed conserved residue clusters are in the interface with the DNA, which could be identified from their amino acid composition; whereas the remaining clusters are at the protein-protein or protein-ligand interface, or embedded in the structural scaffolds. In protein-protein interfaces, conserved residues are highly correlated with experimental residue hotspots, contributing dominantly and often cooperatively to the stability of protein-protein complexes. Overall, the conservation patterns of the stabilizing residues in DNA-binding proteins also highlight the significance of clustering as compared to single residue conservation.     

Relation between protein stability, evolution and structure, as probed by carboxylic acid mutations

Native proteins are marginally stable. Low thermodynamic stability may actually be advantageous, although the accumulation of neutral, destabilizing mutations may have also contributed to it. In any case, once marginal stability has been reached, it appears plausible that mutations at non-constrained positions become fixed in the course of evolution (due to random drift) with frequencies that roughly reflect the mutation effects on stability ("pseudo-equilibrium hypothesis"). We have found that all glutamate-->aspartate mutations in wild-type Escherichia coli thioredoxin are destabilizing, as well as most of the aspartate-->glutamate mutations. Furthermore, the effect of these mutations on thioredoxin thermodynamic stability shows a robust correlation with the frequencies of occurrence of the involved residues in several-hundred sequence alignments derived from a BLAST search. These results provide direct and quantitative experimental evidence for the pseudo-equilibrium hypothesis and should have general consequences for the interpretation of mutation effects on protein stability, as they suggest that residue environments in proteins may be optimized for stabilizing interactions to a remarkable degree of specificity. We also provide evidence that such stabilizing interactions may be detected in sequence alignments, and briefly discuss the implications of this possibility for the derivation of structural information (on native and denatured states) from comparative sequence analyses.     

Sequence-structure signals of 3D domain swapping in proteins

Three-dimensional domain swapping occurs when two or more identical proteins exchange identical parts of their structure to generate an oligomeric unit. It affects proteins with diverse sequences and structures, and is expected to play important roles in evolution, functional regulation and even conformational diseases. Here, we search for traces of domain swapping in the protein sequence, by means of algorithms that predict the structure and stability of proteins using database-derived potentials. Regions whose sequences are not optimal with regard to the stability of the native structure, or showing marked intrinsic preferences for non-native conformations in absence of tertiary interactions are detected in most domain-swapping proteins. These regions are often located in areas crucial in the swapping process and are likely to influence it on a kinetic or thermodynamic level. In addition, cation-pi interactions are frequently observed to zip up the edges of the interface between intertwined chains or to involve hinge loop residues, thereby modulating stability. We end by proposing a set of mutations altering the swapping propensities, whose experimental characterization would contribute to refine our in silico derived hypotheses.     

Evolutionary Capacitance and Control of Protein Stability in Protein-Protein Interaction Networks

In addition to their biological function, protein complexes reduce the exposure of the constituent proteins to the risk of undesired oligomerization by reducing the concentration of the free monomeric state. We interpret this reduced risk as a stabilization of the functional state of the protein. We estimate that protein-protein interactions can account for of additional stabilization; a substantial contribution to intrinsic stability. We hypothesize that proteins in the interaction network act as evolutionary capacitors which allows their binding partners to explore regions of the sequence space which correspond to less stable proteins. In the interaction network of baker''s yeast, we find that statistically proteins that receive higher energetic benefits from the interaction network are more likely to misfold. A simplified fitness landscape wherein the fitness of an organism is inversely proportional to the total concentration of unfolded proteins provides an evolutionary justification for the proposed trends. We conclude by outlining clear biophysical experiments to test our predictions.     

Quantitation of protein-protein interactions by thermal stability shift analysis

Thermal stability shift analysis is a powerful method for examining binding interactions in proteins. We demonstrate that under certain circumstances, protein-protein interactions can be quantitated by monitoring shifts in thermal stability using thermodynamic models and data analysis methods presented in this work. This method relies on the determination of protein stabilities from thermal unfolding experiments using fluorescent dyes such as SYPRO Orange that report on protein denaturation. Data collection is rapid and straightforward using readily available real-time polymerase chain reaction instrumentation. We present an approach for the analysis of the unfolding transitions corresponding to each partner to extract the affinity of the interaction between the proteins. This method does not require the construction of a titration series that brackets the dissociation constant. In thermal shift experiments, protein stability data are obtained at different temperatures according to the affinity- and concentration-dependent shifts in unfolding transition midpoints. Treatment of the temperature dependence of affinity is, therefore, intrinsic to this method and is developed in this study. We used the interaction between maltose-binding protein (MBP) and a thermostable synthetic ankyrin repeat protein (Off7) as an experimental test case because their unfolding transitions overlap minimally. We found that MBP is significantly stabilized by Off7. High experimental throughput is enabled by sample parallelization, and the ability to extract quantitative binding information at a single partner concentration. In a single experiment, we were able to quantify the affinities of a series of alanine mutants, covering a wide range of affinities (～ 100 nM to ～ 100 μM).     

Exploring local and non-local interactions for protein stability by structural motif engineering

In order to probe the relative contribution of local and non-local interactions to the thermodynamic stability of proteins, we have devised an experimental approach based on a combination of motif engineering and sequence shuffling. Candidate chain segments in an immunoglobulin V(L) domain were identified whose conformation is proposed to be dominated by non-local interactions. Locally interacting structural motifs of a different conformation were then constructed as replacements, by introducing motif consensus sequences. We find that all nine replacements we constructed systematically reduce the folding cooperativity. By comparing this destabilising effect with the folding transitions of shuffled sequences for three of these motifs, we estimate the contribution of local, native interactions to the free energy of folding. Our results suggest that local and non-local interactions contribute to stability by an approximately equal amount, but that local interactions stabilise by increasing the resistance to denaturation while non-local interactions increase folding cooperativity. The systematic loss of stability by sequence shuffling in these host-guest experiments suggests that the designed interactions indeed are present in the native state, thus consensus sequence engineering may be a useful tool in structure design, but non-local interactions must be taken into account for global stability engineering. Statistical approaches are powerful tools for engineering protein structure and stability, but an analysis based on local sequence propensities alone does not adequately represent the balance of sequence and context in protein structures.     

First-step mutations for adaptation at elevated temperature increase capsid stability in a virus

The relationship between mutation, protein stability and protein function plays a central role in molecular evolution. Mutations tend to be destabilizing, including those that would confer novel functions such as host-switching or antibiotic resistance. Elevated temperature may play an important role in preadapting a protein for such novel functions by selecting for stabilizing mutations. In this study, we test the stability change conferred by single mutations that arise in a G4-like bacteriophage adapting to elevated temperature. The vast majority of these mutations map to interfaces between viral coat proteins, suggesting they affect protein-protein interactions. We assess their effects by estimating thermodynamic stability using molecular dynamic simulations and measuring kinetic stability using experimental decay assays. The results indicate that most, though not all, of the observed mutations are stabilizing.     

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

Analysis of protein self-association under conditions of the thermodynamic non-ideality

Analysis of protein-protein interactions in highly concentrated solutions requires a consideration of the non-ideality in such solutions which is expressed by the virial coefficients. Different equations are presented to estimate effects of the thermodynamic non-ideality on the macromolecular interaction of self-associating proteins in sedimentation equilibrium experiments. Usually the influence of thermodynamic non-ideal behavior are described by concentration power series. The convergence of such power series is limited at high solute concentration. When expressing the thermodynamic non-ideality by an activity power series this disadvantage can be minimized. The developed centrifuge equations are the basis for a global analysis to estimate equilibrium constants and the corresponding thermodynamic activities of the reactants. Based on fit analysis of synthetic concentration profiles it was established that marked deviations from the expected association constants are observed for proteins with strong association forces between solute molecules. Considerable differences were also observed in weakly interacting systems. This was due to the excluded volume of the protein which is similar in magnitude to the binding constant. For interactions with moderate affinities values extremely close to the true binding values were obtained, as confirmed by experimental results with concanavalin A.     

Structural analyses of proteins electroblotted from native polyacrylamide gels onto polyvinyldiene difluoride membranes. A method for determining the stoichiometry of protein-protein interaction in solution.

The usefulness of a native gel electroblotting technique in the study of protein-protein interactions was demonstrated by the determination of the stoichiometry of the interaction between interleukin-2 (IL-2) and the alpha subunit of IL-2 receptor (IL-2R alpha) in solution. Complexes formed between the recombinant forms of the two proteins in solution were separated from the noncomplexed protein molecules by electrophoresis in a native polyacrylamide gel and the protein bands were electroblotted quantitatively onto polyvinyldiene difluoride membranes for further structural analysis. The data obtained from sequence and amino acid analyses of the blotted proteins provided direct evidence that IL-2 binds to IL-2R alpha in a 1:1 ratio. This methodology should be applicable to the study of other structure/function aspects of protein-protein interactions in solution.     

Protein fusions of BPI with CETP retain functions inherent to each

Cholesteryl ester transfer protein (CETP), bactericidal/permeability inducing protein (BPI), and lipopolysaccharide binding protein (LBP) are members of the lipid transfer/lipopolysaccharide binding protein (LT/LBP) family of proteins that share a common secondary/tertiary structure. Despite this commonality of structure, very different patterns of lipid binding and protein-protein interactions are observed among the family members. BPI was previously shown to retain aspects of its own function when part of it was fused with LBP to form a chimeric protein. We have extended those observations to CETP. Some aspects of cholesteryl ester transfer function can be maintained in a chimeric protein even when over 40% of the sequence is from BPI. Further replacement of an additional 60 amino acids resulted in a complete loss of CETP function even though the chimera was able to retain some BPI-like properties. These artificial fusions retain BPI functions such as lipopolysaccharide (LPS) binding and protein-protein interactions that are not observed with native CETP. BPI-CETP chimeras are inhibited by LPS but cannot be inhibited by small molecule CETP inhibitors as effectively as native CETP. These results localize the site of LPS binding in BPI to a region no larger than the amino terminal 155 amino acids. This region can participate in some protein-protein interactions similar to intact BPI. Chimeras containing the amino terminus of CETP and the carboxy terminus of BPI did not retain any observable CETP function. These results further confirm the modular nature of the LT/LBP family of proteins but also highlight the discrete nature of their individual functions.     

Major structural determinants of transmembrane proteins identified by principal component analysis

We identify amino acid characteristics important in determining the secondary structures of transmembrane proteins, and compare them with characteristics important for cytoplasmic proteins. Using information derived from multiple sequence alignments, we perform a principal component analysis (PCA) to identify the directions in the 20-dimensional amino acid frequency space that comprise the most variance within each protein secondary structure. These vectors represent the important position-specific properties of the amino acids for coils, turns, beta sheets, and alpha helices. As expected, the most important axis for most of the datasets was hydrophobicity. Additional axes, distinct from hydrophobicity, are surprising, especially in the case of transmembrane alpha helices, where the effects of aromaticity and beta-branching are the next two most significant characteristics. The axis representing beta-branching also has equal importance in cytoplasmic and transmembrane helices, a finding that contrasts with some experimental results in membrane-like environments. In a further analysis, we examine trends for some of the PCA axes over averaged transmembrane alpha helices, and find interesting results for aromaticity.     

Beta-turn propensities as paradigms for the analysis of structural motifs to engineer protein stability.

The thermodynamic stability of a protein provides an experimental metric for the relationship of protein sequence and native structure. We have investigated an approach based on an analysis of the structural database for stability engineering of an immunoglobulin variable domain. The most frequently occurring residues in specific positions of beta-turn motifs were predicted to increase the folding stability of mutants that were constructed by site-directed mutagenesis. Even in positions in which different residues are conserved in immunoglobulin sequences, the predictions were confirmed. Frequently, mutants with increased beta-turn propensities display increased folding cooperativities, suggesting pronounced effects on the unfolded state independent of the expected effect on conformational entropy. We conclude that structural motifs with predominantly local interactions can serve as templates with which patterns of sequence preferences can be extracted from the database of protein structures. Such preferences can predict the stability effects of mutations for protein engineering and design.     

Molecular dissection,tissue localization and Ca2+ binding of the ryanodine receptor of Caenorhabditis elegans

We investigated the possible role of residues at the Ccap position in an alpha-helix on protein stability. A set of 431 protein alpha-helices containing a C'-Gly from the Protein Data Bank (PDB) was analyzed, and the normalized frequencies for finding particular residues at the Ccap position, the average fraction of buried surface area, and the hydrogen bonding patterns of the Ccap residue side-chain were calculated. We found that on average the Ccap position is 70% buried and noted a significant correlation (R=0.8) between the relative burial of this residue and its hydrophobicity as defined by the Gibbs energy of transfer from octanol or cyclohexane to water. Ccap residues with polar side-chains are commonly involved in hydrogen bonding. The hydrogen bonding pattern is such that, the longer side-chains of Glu, Gln, Arg, Lys, His form hydrogen bonds with residues distal (>+/-4) in sequence, while the shorter side-chains of Asp, Asn, Ser, Thr exhibit hydrogen bonds with residues close in sequence (<+/-4), mainly involving backbone atoms. Experimentally we determined the thermodynamic propensities of residues at the Ccap position using the protein ubiquitin as a model system. We observed a large variation in the stability of the ubiquitin variants depending on the nature of the Ccap residue. Furthermore, the measured changes in stability of the ubiquitin variants correlate with the hydrophobicity of the Ccap residue. The experimental results, together with the statistical analysis of protein structures from the PDB, indicate that the key hydrophobic capping interactions between a helical residue (C3 or C4) and a residue outside the helix (C", C3' or C4') are frequently enhanced by the hydrophobic interactions with Ccap residues.     

Thermodynamic databases for proteins and protein-nucleic acid interactions

Thermodynamic data regarding proteins and their interactions are important for understanding the mechanisms of protein folding, protein stability, and molecular recognition. Although there are several structural databases available for proteins and their complexes with other molecules, databases for experimental thermodynamic data on protein stability and interactions are rather scarce. Thus, we have developed two electronically accessible thermodynamic databases. ProTherm, Thermodynamic Database for Proteins and Mutants, contains numerical data of several thermodynamic parameters of protein stability, experimental methods and conditions, along with structural, functional, and literature information. ProNIT, Thermodynamic Database for Protein-Nucleic Acid Interactions, contains thermodynamic data for protein-nucleic acid binding, experimental conditions, structural information of proteins, nucleic acids and the complex, and literature information. These data have been incorporated into 3DinSight, an integrated database for structure, function, and properties of biomolecules. A WWW interface allows users to search for data based on various conditions, with different display and sorting options, and to visualize molecular structures and their interactions. These thermodynamic databases, together with structural databases, help researchers gain insight into the relationship among structure, function, and thermodynamics of proteins and their interactions, and will become useful resources for studying proteins in the postgenomic era.     

Identifying transient protein-protein interactions in EphB2 signaling by blue native PAGE and mass spectrometry

Receptor tyrosine kinases (RTKs) are proteins that upon ligand stimulation undergo dimerization and autophosphorylation. Eph receptors (EphRs) are RTKs that are found in different cell types, from both tissues that are developing and from mature tissues, and play important roles in the development of the central nervous system and peripheral nervous system. EphRs also play roles in synapse formation, neural crest formation, angiogenesis and in remodeling the vascular system. Interaction of EphRs with their ephrin ligands lead to activation of signal transduction pathways and formation of many transient protein-protein interactions that ultimately leads to cytoskeletal remodeling. However, the sequence of events at the molecular level is not well understood. We used blue native PAGE and MS to analyze the transient protein-protein interactions that resulted from the stimulation of EphB2 receptors by their ephrinB1-Fc ligands. We analyzed the phosphotyrosine-containing protein complexes immunoprecipitated from the cell lysates of both unstimulated (-) and ephrinB1-Fc-stimulated (+) NG108 cells. Our experiments allowed us to identify many signaling proteins, either known to be part of EphB2 signaling or new for this pathway, which are involved in transient protein-protein interactions upon ephrinB1-Fc stimulation. These data led us to investigate the roles of proteins such as FAK, WAVEs and Nischarin in EphB2 signaling.