Hydrophobic core substitutions in calbindin D9k: effects on stability and structure

Engineering a calcium-dependent conformational change in Calbindin D9k by secondary elements replacement

EF-hand is a common motif in Ca²⁺-binding proteins, some of which present a conformational change upon Ca²⁺-binding, a relevant property for signal transduction. In the present work, we investigated the behavior of Calbindin D_9k, a modulator protein with a high affinity for Ca²⁺ but structurally insensitive to its presence. Its non-canoncal N-terminal EF-hand was replaced by chimeric motifs, containing increasing structural elements from the sensor troponin C SCIII motif. We demonstrated that the loop and helix II were the necessary elements for a conformational change promoted by calcium in chimeric Calbindin D_9k. Fusion of the isolated chimeric motifs to an activity reporter gene showed the loop as the minimal element to promote a conformational change. The discrepancy between these results is discussed in the light of inter-motif interactions and helix I participation in modulating the Ca²⁺ affinity and restricting motif conformation.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Functional and structural analysis of the conserved EFhd2 protein

EFhd2 is a novel protein conserved from C. elegans to H. sapiens. This novel protein was originally identified in cells of the immune and central nervous systems. However, it is most abundant in the central nervous system, where it has been found associated with pathological forms of the microtubule-associated protein tau. The physiological or pathological roles of EFhd2 are poorly understood. In this study, a functional and structural analysis was carried to characterize the molecular requirements for EFhd2's calcium binding activity. The results showed that mutations of a conserved aspartate on either EF-hand motif disrupted the calcium binding activity, indicating that these motifs work in pair as a functional calcium binding domain. Furthermore, characterization of an identified single-nucleotide polymorphisms (SNP) that introduced a missense mutation indicates the importance of a conserved phenylalanine on EFhd2 calcium binding activity. Structural analysis revealed that EFhd2 is predominantly composed of alpha helix and random coil structures and that this novel protein is thermostable. EFhd2's thermo stability depends on its N-terminus. In the absence of the N-terminus, calcium binding restored EFhd2's thermal stability. Overall, these studies contribute to our understanding on EFhd2 functional and structural properties, and introduce it into the family of canonical EF-hand domain containing proteins.

On the mechanism of protein fold-switching by a molecular sensor

Alternate frame folding (AFF) is a mechanism by which conformational change can be engineered into a protein. The protein structure switches from the wild-type fold (N) to a circularly-permuted fold (N'), or vice versa, in response to a signaling event such as ligand binding. Despite the fact that the two native states have similar structures, their interconversion involves folding and unfolding of large parts of the molecule. This rearrangement is reported by fluorescent groups whose relative proximities change as a result of the order-disorder transition. The nature of the conformational change is expected to be similar from protein to protein; thus, it may be possible to employ AFF as a general method to create optical biosensors. Toward that goal, we test basic aspects of the AFF mechanism using the AFF variant of calbindin D(9k). A simple three-state model for fold switching holds that N and N' interconvert through the unfolded state. This model predicts that the fundamental properties of the switch--calcium binding affinity, signal response (i.e., fluorescence change upon binding), and switching rate--can be controlled by altering the relative stabilities of N and N'. We find that selectively destabilizing N or N' changes the equilibrium properties of the switch (binding affinity and signal response) in accordance with the model. However, kinetic data indicate that the switching pathway does not require whole-molecule unfolding. The rate is instead limited by unfolding of a portion of the protein, possibly in concert with folding of a corresponding region.

Protein stability, flexibility and function

Proteins rely on flexibility to respond to environmental changes, ligand binding and chemical modifications. Potentially, a perturbation that changes the flexibility of a protein may interfere with its function. Millions of mutations have been performed on thousands of proteins in quests for a delineation of the molecular details of their function. Several of these mutations interfered with the binding of a specific ligand with a concomitant effect on the stability of the protein scaffold. It has been ambiguous and not straightforward to recognize if any relationships exist between the stability of a protein and the affinity for its ligand. In this review, we present examples of proteins where changes in stability results in changes in affinity and of proteins where stability and affinity are uncorrelated. We discuss the possibility for a relationship between stability and binding. From the data presented is it clear that there are specific sites (flexibility hotspots) in proteins that are important for both binding and stability. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

FH8--a small EF-hand protein from Fasciola hepatica

Vaccine and drug development for fasciolasis rely on a thorough understanding of the mechanisms involved in parasite-host interactions. FH8 is an 8 kDa protein secreted by the parasite Fasciola hepatica in the early stages of infection. Sequence analysis revealed that FH8 has two EF-hand Ca(2+)-binding motifs, and our experimental data show that the protein binds Ca(2+) and that this induces conformational alterations, thus causing it to behave like a sensor protein. Moreover, FH8 displays low affinity for Ca(2+) (K(obs) = 10(4) m(-1)) and is highly stable in its apo and Ca(2+)-loaded states. Homology models were built for FH8 in both states. It has only one globular domain, with two binding sites and appropriate groups in the positions for coordination of the metal ions. However, an unusually high content of positively charged amino acids in one of the binding sites, when compared with the prototypical sensor proteins, potentially affects the protein's affinity for Ca(2+). The only Cys present in FH8, conserved in the homologous proteins of other helminth parasites, is located on the surface, allowing the formation of dimers, detected on SDS gels. These findings reflect specificities of FH8, which are most probably related to its roles both in the parasite and in the host.

A Ca2+-sensing molecular switch based on alternate frame protein folding

Existing strategies for creating biosensors mainly rely on large conformational changes to transduce a binding event to an output signal. Most molecules, however, do not exhibit large-scale structural changes upon substrate binding. Here, we present a general approach (alternate frame folding, or AFF) for engineering allosteric control into ligand binding proteins. AFF can in principle be applied to any protein to establish a binding-induced conformational change, even if none exists in the natural molecule. The AFF design duplicates a portion of the amino acid sequence, creating an additional "frame" of folding. One frame corresponds to the wild-type sequence, and folding produces the normal structure. Folding in the second frame yields a circularly permuted protein. Because the two native structures compete for a shared sequence, they fold in a mutually exclusive fashion. Binding energy is used to drive the conformational change from one fold to the other. We demonstrate the approach by converting the protein calbindin D(9k) into a molecular switch that senses Ca2+. The structures of Ca2+-free and Ca2+-bound calbindin are nearly identical. Nevertheless, the AFF mechanism engineers a robust conformational change that we detect using two covalently attached fluorescent groups. Biological fluorophores can also be employed to create a genetically encoded sensor. AFF should be broadly applicable to create sensors for a variety of small molecules.

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.

A Carboxyl-terminal Hydrophobic Interface Is Critical to Sodium Channel Function

Perturbation of sodium channel inactivation, a finely tuned process that critically regulates the flow of sodium ions into excitable cells, is a common functional consequence of inherited mutations associated with epilepsy, skeletal muscle disease, autism, and cardiac arrhythmias. Understanding the structural basis of inactivation is key to understanding these disorders. Here we identify a novel role for a structural motif in the COOH terminus of the heart NaV1.5 sodium channel in determining channel inactivation. Structural modeling predicts an interhelical hydrophobic interface between paired EF hands in the proximal region of the NaV1.5 COOH terminus. The predicted interface is conserved among almost all EF hand-containing proteins and is the locus of a number of disease-associated mutations. Using the structural model as a guide, we provide biochemical and biophysical evidence that the structural integrity of this interface is necessary for proper Na+ channel inactivation gating. We thus demonstrate a novel role of the sodium channel COOH terminus structure in the control of channel inactivation and in pathologies caused by inherited mutations that disrupt it.

Detecting cryptic epitopes created by nanoparticles

As potential applications of nanotechnology and nanoparticles increase, so too does the likelihood of human exposure to nanoparticles. Because of their small size, nanoparticles are easily taken up into cells (by receptor-mediated endocytosis), whereupon they have essentially free access to all cellular compartments. Similarly to macroscopic biomaterial surfaces (that is, implants), nanoparticles become coated with a layer of adsorbed proteins immediately upon contact with physiological solutions (unless special efforts are taken to prevent this). The process of adsorption often results in conformational changes of the adsorbed protein, which may be affected by the larger curvature of nanoparticles compared with implant surfaces. Protein adsorption may result in the exposure at the surface of amino acid residues that are normally buried in the core of the native protein, which are recognized by the cells as "cryptic epitopes." These cryptic epitopes may trigger inappropriate cellular signaling events (as opposed to being rejected by the cells as foreign bodies). However, identification of such surface-exposed epitopes is nontrivial, and the molecular nature of the adsorbed proteins should be investigated using biological and physical science methods in parallel with systems biology studies of the induced alterations in cell signaling.

Reconstitution of calmodulin from domains and subdomains: Influence of target peptide

Reconstitution studies of a protein from domain fragments can furnish important insights into the distinctive role of particular domain interactions and how they affect biophysical properties important for function. Using isothermal titration calorimetry (ITC) and a number of spectroscopic and chromatographic tools, including CD, fluorescence and NMR spectroscopy, size-exclusion chromatography and non-denaturing agarose gel electrophoresis, we have investigated the reconstitution of the ubiquitous Ca2+-sensor protein calmodulin (CaM) and its globular domains from fragments comprising one or two EF-hands. The studies were carried out with and without the target peptide from smooth muscle myosin light chain kinase (smMLCKp). The CaM-target complex can be reconstituted from the three components consisting of the target peptide and the globular domains TR1C and TR2C. In the absence of peptide, there is no evidence for association of the globular domains. The globular domains can further be reconstituted from their corresponding native subdomains. The dissociation constant, K(D), in 2 mM Tris-HCl (pH 7.5), for the subdomain complexes, EF1:EF2 and EF3:EF4, was determined with ITC to 9.3 x 10(-7) M and 5.9 x 10(-8) M, respectively. Thus, the affinity between the two C-terminal subdomains, located within TR2C, is stronger by a factor of 16 than that between the corresponding subdomains within TR1C. These observations are corroborated by the spectroscopic and chromatographic investigations.

Proteins QSAR with Markov average electrostatic potentials

Classic physicochemical and topological indices have been largely used in small molecules QSAR but less in proteins QSAR. In this study, a Markov model is used to calculate, for the first time, average electrostatic potentials xik for an indirect interaction between aminoacids placed at topologic distances k within a given protein backbone. The short-term average stochastic potential xi1 for 53 Arc repressor mutants was used to model the effect of Alanine scanning on thermal stability. The Arc repressor is a model protein of relevance for biochemical studies on bioorganics and medicinal chemistry. A linear discriminant analysis model developed correctly classified 43 out of 53, 81.1% of proteins according to their thermal stability. More specifically, the model classified 20/28, 71.4% of proteins with near wild-type stability and 23/25, 92.0% of proteins with reduced stability. Moreover, predictability in cross-validation procedures was of 81.0%. Expansion of the electrostatic potential in the series xi0, xi1, xi2, and xi3, justified the use of the abrupt truncation approach, being the overall accuracy >70.0% for xi0 but equal for xi1, xi2, and xi3. The xi1 model compared favorably with respect to others based on D-Fire potential, surface area, volume, partition coefficient, and molar refractivity, with less than 77.0% of accuracy [Ramos de Armas, R.; González-Díaz, H.; Molina, R.; Uriarte, E. Protein Struct. Func. Bioinf.2004, 56, 715]. The xi1 model also has more tractable interpretation than others based on Markovian negentropies and stochastic moments. Finally, the model is notably simpler than the two models based on quadratic and linear indices. Both models, reported by Marrero-Ponce et al., use four-to-five time more descriptors. Introduction of average stochastic potentials may be useful for QSAR applications; having xik amenable physical interpretation and being very effective.

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

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

Predicting stability of Arc repressor mutants with protein stochastic moments

As more and more protein structures are determined and applied to drug manufacture, there is increasing interest in studying their stability. In this study, the stochastic moments ((SR)pi(k)) of 53 Arc repressor mutants were introduced as molecular descriptors modeling protein stability. The Linear Discriminant Analysis model developed correctly classified 43 out of 53, 81.13% of proteins according to their thermal stability. More specifically, the model classified 20/28 (71.4%) proteins with near wild-type stability and 23/25 (92%) proteins with reduced stability. Moreover, validation of the model was carried out by re-substitution procedures (81.0%). In addition, the stochastic moments based model compared favorably with respect to others based on physicochemical and geometric parameters such as D-Fire potential, surface area, volume, partition coefficient, and molar refractivity, which presented less than 77% of accuracy. This result illustrates the possibilities of the stochastic moments' method for the study of bioorganic and medicinal chemistry relevant proteins.

Designing sequence to control protein function in an EF-hand protein

The extent of conformational change that calcium binding induces in EF-hand proteins is a key biochemical property specifying Ca(2+) sensor versus signal modulator function. To understand how differences in amino acid sequence lead to differences in the response to Ca(2+) binding, comparative analyses of sequence and structures, combined with model building, were used to develop hypotheses about which amino acid residues control Ca(2+)-induced conformational changes. These results were used to generate a first design of calbindomodulin (CBM-1), a calbindin D(9k) re-engineered with 15 mutations to respond to Ca(2+) binding with a conformational change similar to that of calmodulin. The gene for CBM-1 was synthesized, and the protein was expressed and purified. Remarkably, this protein did not exhibit any non-native-like molten globule properties despite the large number of mutations and the nonconservative nature of some of them. Ca(2+)-induced changes in CD intensity and in the binding of the hydrophobic probe, ANS, implied that CBM-1 does undergo Ca(2+) sensorlike conformational changes. The X-ray crystal structure of Ca(2+)-CBM-1 determined at 1.44 A resolution reveals the anticipated increase in hydrophobic surface area relative to the wild-type protein. A nascent calmodulin-like hydrophobic docking surface was also found, though it is occluded by the inter-EF-hand loop. The results from this first calbindomodulin design are discussed in terms of progress toward understanding the relationships between amino acid sequence, protein structure, and protein function for EF-hand CaBPs, as well as the additional mutations for the next CBM design.

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

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

Stability scale and atomic solvation parameters extracted from 1023 mutation experiments

The stability scale of 20 amino acid residues is derived from a database of 1023 mutation experiments on 35 proteins. The resulting scale of hydrophobic residues has an excellent correlation with the octanol-to-water transfer free energy corrected with an additional Flory-Huggins molar-volume term (correlation coefficient r = 0.95, slope = 1.05, and a near zero intercept). Thus, hydrophobic contribution to folding stability is characterized remarkably well by transfer experiments. However, no corresponding correlation is found for hydrophilic residues. Both the hydrophilic portion and the entire scale, however, correlate strongly with average burial accessible surface (r = 0.76 and 0.97, respectively). Such a strong correlation leads to a near uniform value of the atomic solvation parameters for atoms C, S, O/N, O(-0.5), and N(+0.5,1). All are in the range of 12-28 cal x mol(-1) A(-2), close to the original estimate of hydrophobic contribution of 25-30 cal x mol(-1) A(-2) to folding stability. Without any adjustable parameters, the new stability scale and new atomic solvation parameters yielded an accurate prediction of protein-protein binding free energy for a separate database of 21 protein-protein complexes (r = 0.80 and slope = 1.06, and r = 0.83 and slope = 0.93, respectively).

The EF-hand domain: a globally cooperative structural unit

EF-hand Ca(2+)-binding proteins participate in both modulation of Ca(2+) signals and direct transduction of the ionic signal into downstream biochemical events. The range of biochemical functions of these proteins is correlated with differences in the way in which they respond to the binding of Ca(2+). The EF-hand domains of calbindin D(9k) and calmodulin are homologous, yet they respond to the binding of calcium ions in a drastically different manner. A series of comparative analyses of their structures enabled the development of hypotheses about which residues in these proteins control the calcium-induced changes in conformation. To test our understanding of the relationship between protein sequence and structure, we specifically designed the F36G mutation of the EF-hand protein calbindin D(9k) to alter the packing of helices I and II in the apoprotein. The three-dimensional structure of apo F36G was determined in solution by nuclear magnetic resonance spectroscopy and showed that the design was successful. Surprisingly, significant structural perturbations also were found to extend far from the site of mutation. The observation of such long-range effects provides clear evidence that four-helix EF-hand domains should be treated as a single globally cooperative unit. A hypothetical mechanism for how the long-range effects are transmitted is described. Our results support the concept of energetic and structural coupling of the key residues that are crucial for a protein's fold and function.

Structure determination of human and murine beta-defensins reveals structural conservation in the absence of significant sequence similarity

Defensins are cationic and cysteine-rich peptides that play a crucial role in the host defense against microorganisms of many organisms by their capability to permeabilize bacterial membranes. The low sequence similarity among the members of the large mammalian beta-defensin family suggests that their antimicrobial activity is largely independent of their primary structure. To investigate to what extent these defensins share a similar fold, the structures of the two human beta-defensins, hBD-1 and hBD-2, as well as those of two novel murine defensins, termed mBD-7 and mBD-8, were determined by nuclear magnetic resonance spectroscopy. All four defensins investigated share a striking similarity on the level of secondary and tertiary structure including the lack of a distinct hydrophobic core, suggesting that the fold is mainly stabilized by the presence of three disulfide bonds. In addition to the overall shape of the molecules, the ratio of solvent-exposed polar and hydrophobic side chains is also very similar among the four defensins investigated. It is significant that beta-defensins do not exhibit a common pattern of charged and hydrophobic residues on the protein surface and that the beta-defensin-specific fold appears to accommodate a wide range of different amino acids at most sequence positions. In addition to the implications for the mode of biological defensin actions, these findings are of particular interest because beta-defensins have been suggested as lead compounds for the development of novel peptide antibiotics for the therapy of infectious diseases.

Structure of rat parvalbumin with deleted AB domain: implications for the evolution of EF hand calcium-binding proteins and possible physiological relevance

Among the EF-hand Ca(2+)-binding proteins, parvalbumin (PV) and calbindin D9k (CaB) have the function of Ca(2+) buffers. They evolved from an ancestor protein through two phylogenetic pathways, keeping one pair of EF-hands. They differ by the extra helix-loop-helix (AB domain) found in PV and by the linker between the binding sites. To investigate whether the deletion of AB in PV restores a CaB-like structure, we prepared and solved the structure of the truncated rat PV (PVratDelta37) by X-ray and NMR. PVratDelta37 keeps the PV fold, but is more compact, having a well-structured linker, which differs remarkably from CaB. PvratDelta37 has no stable apo-form, has lower affinity for Ca(2+) than full-length PV, and does not bind Mg(2+), in contrast to CaB. Structural differences of the hydrophobic core are partially responsible for lowering the calcium-binding affinity of the truncated protein. It can be concluded that the AB domain, like the linker of CaB, plays a role in structural stabilization. The AB domain of PV protects the hydrophobic core, and is required to maintain high affinity for divalent cation binding. Therefore, the AB domain possibly modulates PV buffer function.

Four-body potentials reveal protein-specific correlations to stability changes caused by hydrophobic core mutations

Mutational experiments show how changes in the hydrophobic cores of proteins affect their stabilities. Here, we estimate these effects computationally, using four-body likelihood potentials obtained by simplicial neighborhood analysis of protein packing (SNAPP). In this procedure, the volume of a known protein structure is tiled with tetrahedra having the center of mass of one amino acid side-chain at each vertex. Log-likelihoods are computed for the 8855 possible tetrahedra with equivalent compositions from structural databases and amino acid frequencies. The sum of these four-body potentials for tetrahedra present in a given protein yields the SNAPP score. Mutations change this sum by changing the compositions of tetrahedra containing the mutated residue and their related potentials. Linear correlation coefficients between experimental mutational stability changes, Delta(DeltaG(unfold)), and those based on SNAPP scoring range from 0.70 to 0.94 for hydrophobic core mutations in five different proteins. Accurate predictions for the effects of hydrophobic core mutations can therefore be obtained by virtual mutagenesis, based on changes to the total SNAPP likelihood potential. Significantly, slopes of the relation between Delta(DeltaG(unfold)) and DeltaSNAPP for different proteins are statistically distinct, and we show that these protein-specific effects can be estimated using the average SNAPP score per residue, which is readily derived from the analysis itself. This result enhances the predictive value of statistical potentials and supports previous suggestions that "comparable" mutations in different proteins may lead to different Delta(DeltaG(unfold)) values because of differences in their flexibility and/or conformational entropy.

The dimerization interface of the metastasis-associated protein S100A4 (Mts1): in vivo and in vitro studies

The S100 calcium-binding proteins are implicated in signal transduction, motility, and cytoskeletal dynamics. The three-dimensional structure of several S100 proteins revealed that the proteins form non-covalent dimers. However, the mechanism of the S100 dimerization is still obscure. In this study we characterized the dimerization of S100A4 (also named Mts1) in vitro and in vivo. Analytical ultracentrifugation revealed that apoS100A4 was present in solution as a mixture of monomers and dimers in a rapidly reversible equilibrium (K(d) = 4 +/- 2 microm). The binding of calcium promoted dimerization. Replacement of Tyr-75 by Phe resulted in the stabilization of the dimer. Helix IV is known to form the major part of the dimerization interface in homologous S100 proteins. By using the yeast two-hybrid system we showed that only a few residues of helix IV, namely Phe-72, Tyr-75, Phe-78, and Leu-79, are essential for dimerization in vivo. A homology model demonstrated that these residues form a hydrophobic cluster on helix IV. Their role is to stabilize the structure of individual subunits rather than provide specific interactions across the dimerization surface. Our mutation data showed that the specificity at the dimerization surface is not particularly stringent, which is consistent with recent data indicating that S100 proteins can form heterodimers.

Solvation energetics and conformational change in EF-hand proteins

Calmodulin and other members of the EF-hand protein family are known to undergo major changes in conformation upon binding Ca(2+). However, some EF-hand proteins, such as calbindin D9k, bind Ca(2+) without a significant change in conformation. Here, we show the importance of a precise balance of solvation energetics to conformational change, using mutational analysis of partially buried polar groups in the N-terminal domain of calmodulin (N-cam). Several variants were characterized using fluorescence, circular dichroism, and NMR spectroscopy. Strikingly, the replacement of polar side chains glutamine and lysine at positions 41 and 75 with nonpolar side chains leads to dramatic enhancement of the stability of the Ca(2+)-free state, a corresponding decrease in Ca(2+)-binding affinity, and an apparent loss of ability to change conformation to the open form. The results suggest a paradigm for conformational change in which energetic strain is accumulated in one state in order to modulate the energetics of change to the alternative state.

Residue depth: a novel parameter for the analysis of protein structure and stability

BACKGROUND

Accessible surface area is a parameter that is widely used in analyses of protein structure and stability. Accessible surface area does not, however, distinguish between atoms just below the protein surface and those in the core of the protein. In order to differentiate between such buried residues we describe a computational procedure for calculating the depth of a residue from the protein surface.

RESULTS

Residue depth correlates significantly better than accessibility with effects of mutations on protein stability and on protein-protein interactions. The deepest residues in the native state invariably undergo hydrogen exchange by global unfolding of the protein and are often significantly protected in the corresponding molten-globule states.

CONCLUSIONS

Depth is often a more useful gage of residue burial than accessibility. This is probably related to the fact that the protein interior and surrounding solvent differ significantly in polarity and packing density. Hence, the strengths of van der Waals and electrostatic interactions between residues in a protein might be expected to depend on the distance of the residue(s) from the protein surface.