Exploring additivity effects of double mutations on the binding affinity of protein‐protein complexes

Additivity in binding affinity of protein‐protein complexes refers to the change in free energy of binding (ΔΔG_bind) for double (or multiple) mutations which is approximately equal to the sum of their corresponding single mutation ΔΔG_bind values. In this study, we have explored the additivity effect of double mutants, which shows a linear relationship between the binding affinity of double and sum of single mutants with a correlation of 0.90. However, the comparison of ΔΔG_bind values showed a mean absolute deviation of 0.86 kcal/mol, and 25.6% of the double mutants show a deviation of more than 1 kcal/mol, which are identified as non‐additive. The additivity effects have been analyzed based on the influence of structural features such as accessible surface area, long range order, binding propensity change, surrounding hydrophobicity, flexibility, atomic contacts between the mutations and distance between the 2 mutations. We found that non‐additive mutations tend to be closer to each other and have more contacts. We have also used machine learning methods to discriminate additive and non‐additive mutations using structure‐based features, which showed the accuracies in the range of 0.77–0.92 for protein‐protein complexes belonging to different functions. Further, we have compared the additivity effects of protein stability along with binding affinity and explored the similarities and differences between them. The results obtained in this study provide insights into the effects of various structural features on binding affinity of double mutants, and will aid the development of accurate methods to predict the binding affinity of double mutants.     

Investigation of phosphotyrosine recognition by the SH2 domain of the Src kinase.

The binding of tyrosine phosphorylated targets by SH2 domains is required for propagation of many cellular signals in higher eukaryotes; however, the determinants of phosphotyrosine (pTyr) recognition by SH2 domains are not well understood. In order to identify the attributes of pTyr required for high affinity interaction with SH2 domains, the binding of the SH2 domain of the Src kinase (Src SH2 domain) to a dephosphorylated peptide, a phosphoserine-containing peptide, and the amino acid pTyr was studied using titration calorimetry and compared with the binding of a high affinity tyrosyl phosphopeptide. The dephosphorylated peptide and the phosphoserine containing peptide both bind extremely weakly to the Src SH2 domain (DeltaGo (dephosphorylated)=-3.6 kcal/mol, DeltaGo (phosphoserine) >-3.7 kcal/mol); however, the DeltaGo value of pTyr binding is more favorable (-4.7 kcal/mol, or 50 % of the entire binding free energy of a high affinity tyrosyl phosphopeptide). These results indicate that both the phosphate and the tyrosine ring of the pTyr are critical determinants of high affinity binding. Alanine mutagenesis was also used to evaluate the energetic contribution to binding of ten residues located in the pTyr-binding site. Mutation of the strictly conserved Arg betaB5 resulted in a large increase in DeltaGo (DeltaDeltaGo=3.2 kcal/mol) while elimination of the other examined residues each resulted in a significantly smaller (DeltaDeltaGo<1.4 kcal/mol) reduction in affinity, indicating that Arg betaB5 is the single most important determinant of pTyr recognition. However, mutation of Cys betaC3, a residue unique to the Src SH2 domain, surprisingly increased affinity by eightfold (DeltaDeltaGo=-1.1 kcal/mol). Using a double mutant cycle analysis, it was revealed that residues of the pTyr-binding pocket are not coupled to the peptide residues C-terminal to the pTyr. In addition, comparison of each residue's DeltaDeltaGo value upon mutation with that residue's sequence conservation among SH2 domains revealed only a modest correlation between a residue's energetic contribution to pTyr recognition and its conservation throughout evolution. The results of this investigation highlight the importance of a single critical interaction, the buried ionic bond between the phosphate of the pTyr and Arg betaB5 of the SH2 domain, driving the binding of SH2 domains to tyrosine phosphorylated targets.     

Relative Binding Affinity Prediction of Charge-Changing Sequence Mutations with FEP in Protein–Protein Interfaces

Building on the substantial progress that has been made in using free energy perturbation (FEP) methods to predict the relative binding affinities of small molecule ligands to proteins, we have previously shown that results of similar quality can be obtained in predicting the effect of mutations on the binding affinity of protein–protein complexes. However, these results were restricted to mutations which did not change the net charge of the side chains due to known difficulties with modeling perturbations involving a change in charge in FEP. Various methods have been proposed to address this problem. Here we apply the co-alchemical water approach to study the efficacy of FEP calculations of charge changing mutations at the protein–protein interface for the antibody–gp120 system investigated previously and three additional complexes. We achieve an overall root mean square error of 1.2 kcal/mol on a set of 106 cases involving a change in net charge selected by a simple suitability filter using side-chain predictions and solvent accessible surface area to be relevant to a biologic optimization project. Reasonable, although less precise, results are also obtained for the 44 more challenging mutations that involve buried residues, which may in some cases require substantial reorganization of the local protein structure, which can extend beyond the scope of a typical FEP simulation. We believe that the proposed prediction protocol will be of sufficient efficiency and accuracy to guide protein engineering projects for which optimization and/or maintenance of a high degree of binding affinity is a key objective.     

Calorimetric examination of high-affinity Src SH2 domain-tyrosyl phosphopeptide binding: dissection of the phosphopeptide sequence specificity and coupling energetics

SH2 domains are protein modules which interact with specific tyrosine phosphorylated sequences in target proteins. The SH2 domain of the Src kinase binds with high affinity to a tyrosine phosphorylated peptide containing the amino acids Glu, Glu, and Ile (EEI) at the positions +1, +2, and +3 C-terminal to the phosphotyrosine, respectively. To investigate the degree of selectivity of the Src SH2 domain for each amino acid of the EEI motif, the binding thermodynamics of a panel of substitutions at the +1 (Gln, Asp, Ala, Gly), +2 (Gln, Asp, Ala, Gly), and +3 (Leu, Val, Ala, Gly) positions were examined using titration microcalorimetry. It was revealed that the Src SH2 domain is insensitive (DeltaDeltaG degrees 相似文献   

Identification of key amino acid residues in the assembly of enzymes into the pyruvate dehydrogenase complex of Bacillus stearothermophilus: a kinetic and thermodynamic analysis

Structural studies have shown that electrostatic interactions play a major part in the binding of dihydrolipoyl dehydrogenase (E3) to the peripheral subunit-binding domain (PSBD) of the dihydrolipoyl acyltransferase (E2) in the assembly of the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus. The binding is characterized by a small, unfavorable enthalpy change (deltaH degrees = +2.2 kcal/mol) and a large, positive entropy change (TdeltaS degrees = +14.8 kcal/mol). The contributions of individual surface residues of the PSBD of E2 to its interaction with E3 have been assessed by alanine-scanning mutagenesis, surface plasmon resonance detection, and isothermal titration calorimetry. The mutation R135A in the PSBD gave rise to a significant decrease (120-fold) in the binding affinity; two other mutations (R139A and R156A) were associated with smaller effects. The binding of the R135A mutant to E3 was accompanied by a favorable enthalpy (deltaH degrees = -2.6 kcal/mol) and a less positive entropy change (TdeltaS degrees = +7.2 kcal/mol). The midpoint melting temperature (T(m)) of E3-PSBD complexes was determined by differential scanning calorimetry. The R135A mutation caused a significant decrease (5 degrees C) in the T(m), compared with the wild-type complex. The results reveal the importance of Arg135 of the PSBD as a key residue in the molecular recognition of E3 by E2, and as a major participant in the overall entropy-driven binding process. Further, the effects of mutagenesis on the deltaCp of subunit association illustrate the difficulties in attributing changes in heat capacity to specific classes of interactions.     

Empirical calculation of the relative free energies of peptide binding to the molecular chaperone DnaK

We describe a methodology to calculate the relative free energies of protein-peptide complex formation. The interaction energy was decomposed into nonpolar, electrostatic and entropic contributions. A free energy-surface area relationship served to calculate the nonpolar free energy term. The electrostatic free energy was calculated with the finite difference Poisson-Boltzmann method and the entropic contribution was estimated from the loss in the conformational entropy of the peptide side chains. We applied this methodology to a series of DnaK*peptide complexes. On the basis of the single known crystal structure of the peptide-binding domain of DnaK with a bound heptapeptide, we modeled ten other DnaK*heptapeptide complexes with experimentally measured K(d) values from 0.06 microM to 11 microM, using molecular dynamics to refine the structures of the complexes. Molecular dynamic trajectories, after equilibration, were used for calculating the energies with greater accuracy. The calculated relative binding free energies were compared with the experimentally determined free energies. Linear scaling of the calculated terms was applied to fit them to the experimental values. The calculated binding free energies were between -7.1 kcal/mol and - 9.4 kcal/mol with a correlation coefficient of 0.86. The calculated nonpolar contributions are mainly due to the central hydrophobic binding pocket of DnaK for three amino acid residues. Negative electrostatic fields generated by the protein increase the binding affinity for basic residues flanking the hydrophobic core of the peptide ligand. Analysis of the individual energy contributions indicated that the nonpolar contributions are predominant compared to the other energy terms even for peptides with low affinity and that inclusion of the change in conformational entropy of the peptide side chains does not improve the discriminative power of the calculation. The method seems to be useful for predicting relative binding energies of peptide ligands of DnaK and might be applicable to other protein-peptide systems, particularly if only the structure of one protein-ligand complex is available.     

Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.     

Protein-protein recognition: an experimental and computational study of the R89K mutation in Raf and its effect on Ras binding

Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras-Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [gamma-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9+/-1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein-protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 A away, namely residues 37-41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2-3 specific water molecules. This analysis will provide insights into the large experimental database of Ras-Raf mutations.     

Thermodynamic basis of resistance to HIV-1 protease inhibition: calorimetric analysis of the V82F/I84V active site resistant mutant

One of the most serious side effects associated with the therapy of HIV-1 infection is the appearance of viral strains that exhibit resistance to protease inhibitors. The active site mutant V82F/I84V has been shown to lower the binding affinity of protease inhibitors in clinical use. To identify the origin of this effect, we have investigated the binding thermodynamics of the protease inhibitors indinavir, ritonavir, saquinavir, and nelfinavir to the wild-type HIV-1 protease and to the V82F/I84V resistant mutant. The main driving force for the binding of all four inhibitors is a large positive entropy change originating from the burial of a significant hydrophobic surface upon binding. At 25 degrees C, the binding enthalpy is unfavorable for all inhibitors except ritonavir, for which it is slightly favorable (-2.3 kcal/mol). Since the inhibitors are preshaped to the geometry of the binding site, their conformational entropy loss upon binding is small, a property that contributes to their high binding affinity. The V82F/I84V active site mutation lowers the affinity of the inhibitors by making the binding enthalpy more positive and making the entropy change slightly less favorable. The effect on the enthalpy change is, however, the major one. The predominantly enthalpic effect of the V82F/I84V mutation is consistent with the idea that the introduction of the bulkier Phe side chain at position 82 and the Val side chain at position 84 distort the binding site and weaken van der Waals and other favorable interactions with inhibitors preshaped to the wild-type binding site. Another contribution of the V82F/I84V to binding affinity originates from an increase in the energy penalty associated with the conformational change of the protease upon binding. The V82F/I84V mutant is structurally more stable than the wild-type protease by about 1.4 kcal/mol. This effect, however, affects equally the binding affinity of substrate and inhibitors.     

Free energy simulation to investigate the effect of amino acid sequence environment on the severity of osteogenesis imperfecta by glycine mutations in collagen

Molecular dynamics simulations were carried out to calculate the free energy change difference of two collagen-like peptide models for Gly --> Ser mutations causing two different osteogenesis imperfecta phenotypes. These simulations were performed to investigate the impact of local amino acid sequence environment adjacent to a mutation site on the stability of the collagen. The average free energy differences for a Gly --> Ser mutant relative to a wild type are 3.4 kcal/mol and 8.2 kcal/mol for a nonlethal site and a lethal site, respectively. The free energy change differences of mutant containing two Ser residues relative to the wild type at the nonlethal and lethal mutation sites are 4.6 and 9.8 kcal/mol, respectively. Although electrostatic interactions stabilize mutants containing one or two Ser residues at both mutation sites, van der Waals interactions are of sufficient magnitude to cause a net destabilization. The presence of Gln and Arg near the mutation site, which contain large and polar side chains, provide more destabilization than amino acids containing small and nonpolar side chains.     

iSEE: Interface structure,evolution, and energy-based machine learning predictor of binding affinity changes upon mutations

Quantitative evaluation of binding affinity changes upon mutations is crucial for protein engineering and drug design. Machine learning-based methods are gaining increasing momentum in this field. Due to the limited number of experimental data, using a small number of sensitive predictive features is vital to the generalization and robustness of such machine learning methods. Here we introduce a fast and reliable predictor of binding affinity changes upon single point mutation, based on a random forest approach. Our method, iSEE, uses a limited number of interface Structure, Evolution, and Energy-based features for the prediction. iSEE achieves, using only 31 features, a high prediction performance with a Pearson correlation coefficient (PCC) of 0.80 and a root mean square error of 1.41 kcal/mol on a diverse training dataset consisting of 1102 mutations in 57 protein-protein complexes. It competes with existing state-of-the-art methods on two blind test datasets. Predictions for a new dataset of 487 mutations in 56 protein complexes from the recently published SKEMPI 2.0 database reveals that none of the current methods perform well (PCC < 0.42), although their combination does improve the predictions. Feature analysis for iSEE underlines the significance of evolutionary conservations for quantitative prediction of mutation effects. As an application example, we perform a full mutation scanning of the interface residues in the MDM2–p53 complex.     

In Silico analysis of voltage-gated sodium channel in relation to DDT resistance in vector mosquitoes

The voltage-gated sodium channel (VGSC) is the target site for insecticides such as DDT and synthetic pyrethroids. A single base (A-T) change in the knock-down resistance (kdr) allele leads to an amino acid substitution at position 267 that confers the target-mediated resistance to DDT and synthetic pyrethroids in Anopheles gambiae. A theoretical model of the VGSC domain II that contains the site of mutation was constructed using the K;+ channel protein of Aeropyrum pernix as a template. The validated model with 88.6% residues in the favored region was subjected to the CASTp program that predicted 30 pockets in the modeled domain II for ligand interaction. In the model, at position 267, leucine was manually replaced with phenylalanine. When this altered model was subjected to the CASTp program, the search results showed the same number of pockets. The docking results indicate that DDT interacts with the modeled VGSC domain II at position 275 in the presence of leucine or in the presence of phenylalanine (binding energy =-5.32 kcal/mol, -6.21 kcal/mol). It appears from the results that the mutation at position 267 has no direct influence on the interaction of DDT with the target protein. Therefore, to understand the interaction affinity of DDT with the target and influence of the mutation on the existence of active sites/pockets in relation to ligand binding, a whole VGSC model is necessary.     

Thermodynamic effects of proline introduction on protein stability

The amino acid Pro is more rigid than other naturally occurring amino acids and, in proteins, lacks an amide hydrogen. To understand the structural and thermodynamic effects of Pro substitutions, it was introduced at 13 different positions in four different proteins, leucine-isoleucine-valine binding protein, maltose binding protein, ribose binding protein, and thioredoxin. Three of the maltose binding protein mutants were characterized by X-ray crystallography to confirm that no structural changes had occurred upon mutation. In the remaining cases, fluorescence and CD spectroscopy were used to show the absence of structural change. Stabilities of wild type and mutant proteins were characterized by chemical denaturation at neutral pH and by differential scanning calorimetry as a function of pH. The mutants did not show enhanced stability with respect to chemical denaturation at room temperature. However, 6 of the 13 single mutants showed a small but significant increase in the free energy of thermal unfolding in the range of 0.3-2.4 kcal/mol, 2 mutants showed no change, and 5 were destabilized. In five of the six cases, the stabilization was because of reduced entropy of unfolding. However, the magnitude of the reduction in entropy of unfolding was typically several fold larger than the theoretical estimate of -4 cal K(-1) mol(-1) derived from the relative areas in the Ramachandran map accessible to Pro and Ala residues, respectively. Two double mutants were constructed. In both cases, the effects of the single mutations on the free energy of thermal unfolding were nonadditive.     

Sequence and structural features of binding site residues in protein-protein complexes: comparison with protein-nucleic acid complexes

Background

Protein-protein interactions are important for several cellular processes. Understanding the mechanism of protein-protein recognition and predicting the binding sites in protein-protein complexes are long standing goals in molecular and computational biology.

Methods

We have developed an energy based approach for identifying the binding site residues in protein–protein complexes. The binding site residues have been analyzed with sequence and structure based parameters such as binding propensity, neighboring residues in the vicinity of binding sites, conservation score and conformational switching.

Results

We observed that the binding propensities of amino acid residues are specific for protein-protein complexes. Further, typical dipeptides and tripeptides showed high preference for binding, which is unique to protein-protein complexes. Most of the binding site residues are highly conserved among homologous sequences. Our analysis showed that 7% of residues changed their conformations upon protein-protein complex formation and it is 9.2% and 6.6% in the binding and non-binding sites, respectively. Specifically, the residues Glu, Lys, Leu and Ser changed their conformation from coil to helix/strand and from helix to coil/strand. Leu, Ser, Thr and Val prefer to change their conformation from strand to coil/helix.

Conclusions

The results obtained in this study will be helpful for understanding and predicting the binding sites in protein-protein complexes.

     

Structural and computational characterization of the SHV-1 beta-lactamase-beta-lactamase inhibitor protein interface

Beta-lactamase inhibitor protein (BLIP) binds a variety of class A beta-lactamases with affinities ranging from micromolar to picomolar. Whereas the TEM-1 and SHV-1 beta-lactamases are almost structurally identical, BLIP binds TEM-1 approximately 1000-fold tighter than SHV-1. Determining the underlying source of this affinity difference is important for understanding the molecular basis of beta-lactamase inhibition and mechanisms of protein-protein interface specificity and affinity. Here we present the 1.6A resolution crystal structure of SHV-1.BLIP. In addition, a point mutation was identified, SHV D104E, that increases SHV.BLIP binding affinity from micromolar to nanomolar. Comparison of the SHV-1.BLIP structure with the published TEM-1.BLIP structure suggests that the increased volume of Glu-104 stabilizes a key binding loop in the interface. Solution of the 1.8A SHV D104K.BLIP crystal structure identifies a novel conformation in which this binding loop is removed from the interface. Using these structural data, we evaluated the ability of EGAD, a program developed for computational protein design, to calculate changes in the stability of mutant beta-lactamase.BLIP complexes. Changes in binding affinity were calculated within an error of 1.6 kcal/mol of the experimental values for 112 mutations at the TEM-1.BLIP interface and within an error of 2.2 kcal/mol for 24 mutations at the SHV-1.BLIP interface. The reasonable success of EGAD in predicting changes in interface stability is a promising step toward understanding the stability of the beta-lactamase.BLIP complexes and computationally assisted design of tight binding BLIP variants.     

Thermodynamic stability of the bacteriorhodopsin lattice as measured by lipid dilution

To determine the strength of noncovalent interactions that stabilize a membrane protein complex, we have developed an in vitro method for quantifying the dissociation of the bacteriorhodopsin (BR) lattice, a naturally occurring two-dimensional crystal. A lattice suspension was titrated with a short- and long-chain phosphatidylcholine mixture to dilute BR within the lipid bilayer. The fraction of BR in the lattice form as a function of added lipid was determined by visible circular dichroism spectroscopy and fit with a cooperative self-assembly model to obtain a critical concentration for lattice assembly. Critical concentration values of wild-type and mutant proteins were used to calculate the change in lattice stability upon mutation (DeltaDeltaG). By using this method, a series of mutant proteins was examined in which residues at the BR-BR interface were replaced with smaller amino acids, either Ala or Gly. Most of the mutant lattices were destabilized, with DeltaDeltaG values of 0.2-1.1 kcal/mol at 30 degrees C, consistent with favorable packing of apolar residues in the membrane. One mutant, I45A, was stabilized by approximately 1.0 kcal/mol, possibly due to increased lipid entropy. The DeltaDeltaG values agreed well with previous in vivo measurements, except in the case of I45A. The ability to measure the change in stability of mutant protein complexes in a lipid bilayer may provide a means of determining the contributions of specific protein-protein and protein-lipid interactions to membrane protein structure.     

Dissection of binding interactions in the complex between the anti-lysozyme antibody HyHEL-63 and its antigen

Alanine-scanning mutagenesis, X-ray crystallography, and double mutant cycles were used to characterize the interface between the anti-hen egg white lysozyme (HEL) antibody HyHEL-63 and HEL. Eleven HEL residues in contact with HyHEL-63 in the crystal structure of the antigen-antibody complex, and 10 HyHEL-63 residues in contact with HEL, were individually truncated to alanine in order to determine their relative contributions to complex stabilization. The residues of HEL (Tyr20, Lys96, and Lys97) most important for binding HyHEL-63 (Delta G(mutant) - Delta G(wild type) > 3.0 kcal/mol) form a contiguous patch at the center of the surface contacted by the antibody. Hot spot residues of the antibody (Delta Delta G > 2.0 kcal/mol) are organized in two clusters that juxtapose hot spot residues of HEL, resulting in energetic complementarity across the interface. All energetically critical residues are centrally located, shielded from solvent by peripheral residues that contribute significantly less to the binding free energy. Although HEL hot spot residues Lys96 and Lys97 make similar interactions with antibody in the HyHEL-63/HEL complex, alanine substitution of Lys96 results in a nearly 100-fold greater reduction in affinity than the corresponding mutation in Lys97. To understand the basis for this marked difference, we determined the crystal structures of the HyHEL-63/HEL Lys96Ala and HyHEL-63/HEL Lys97Ala complexes to 1.80 and 1.85 A resolution, respectively. Whereas conformational changes in the proteins and differences in the solvent networks at the mutation sites appear too small to explain the observed affinity difference, superposition of free HEL in different crystal forms onto bound HEL in the wild type and mutant HyHEL-63/HEL complexes reveals that the side-chain conformation of Lys96 is very similar in the various structures, but that the Lys97 side chain displays considerable flexibility. Accordingly, a greater entropic penalty may be associated with quenching the mobility of the Lys97 than the Lys96 side chain upon complex formation, reducing binding. To further dissect the energetics of specific interactions in the HyHEL-63/HEL interface, double mutant cycles were constructed to measure the coupling of 13 amino acid pairs, 11 of which are in direct contact in the crystal structure. A large coupling energy, 3.0 kcal/mol, was found between HEL residue Lys97 and HyHEL-63 residue V(H)Asp32, which form a buried salt bridge surrounded by polar residues of the antigen. Thus, in contrast to protein folding where buried salt bridges are generally destabilizing, salt bridges in protein-protein interfaces, whose residual composition is more hydrophilic than that of protein interiors, may contribute significantly to complex stabilization.     

Analysis of the binding energies of testosterone, 5alpha-dihydrotestosterone,androstenedione and dehydroepiandrosterone sulfate with an antitestosterone antibody

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the binding of testosterone (TES), 5alpha-dihydrotestosterone (5ADHT), androstenedione (AND), and dehydroepiandrosterone sulfate (DHEAS) to the monoclonal antitestosterone antibody 3-C(4)F(5). The relative binding free energy of TES and AND was also calculated with free energy perturbation (FEP) simulations. The antibody 3-C(4)F(5) has a relatively high affinity (3 x 10(8) M(-1)) and on overall good binding profile for testosterone but its cross-reactivity with DHEAS has been the main reason for the failure to use this antibody in clinical immunoassays. The relative binding free energies obtained with the MM-PBSA method were 1.5 kcal/mol for 5ADHT, 3.8 kcal/mol for AND, and 4.3 kcal/mol for DHEAS, as compared to TES. When a water molecule of the ligand binding site, observed in the antibody-TES crystal structure, was explicitly included in MM-PBSA calculations, the relative binding energies were 3.4, 4.9, and 5.4 kcal/mol for 5ADHT, AND, and DHEAS, respectively. The calculated numbers are in correct order but larger than the corresponding experimental energies of 1.3, 1.5, and 2.6 kcal/mol, respectively. The fact that the MM-PBSA method reproduced the relative binding free energies of DHEAS, a steroid having a negatively charged sulfate group, and the neutrally charged TES, 5ADHT, and AND in satisfactory agreement with experiment shows the robustness of the method in predicting relative binding affinities. The 800-ps FEP simulations predicted that the antibody 3-C(4)F(5) binds TES 1.3 kcal/mol tighter than AND. Computational mutagenesis of selected amino acid residues of the ligand binding site revealed that the lower affinities of AND and DHEAS as compared to TES are due to a combined effect of several residues, each contributing a small fraction to the tighter binding of TES. An exception to this is Tyr99H, whose mutation to Ala lowered the binding of DHEAS 0.7 kcal/mol more than the binding of TES. This is probably due to the hydrogen bonding interaction formed between the OH group of Tyr99H and the sulfate group of DHEAS. Computational mutagensis data also showed that the affinity of the steroids to the antitestosterone antibody 3-C(4)F(5) would be enhanced if Trp47H were repositioned so that it would make more extensive contacts with the bound ligands. In addition, the binding of steroids to antitestosterone, antiprogesterone, and antiestradiol antibodies is discussed.