Dissecting cooperative and additive binding energetics in the affinity maturation pathway of a protein-protein interface

When two proteins associate they form a molecular interface that is a structural and energetic mosaic. Within such interfaces, individual amino acid residues contribute distinct binding energies to the complex. In combination, these energies are not necessarily additive, and significant positive or negative cooperative effects often exist. The basis of reliable algorithms to predict the specificities and energies of protein-protein interactions depends critically on a quantitative understanding of this cooperativity. We have used a model protein-protein system defined by an affinity maturation pathway, comprising variants of a T cell receptor Vbeta domain that exhibit an overall affinity range of approximately 1500-fold for binding to the superantigen staphylococcal enterotoxin C3, in order to dissect the cooperative and additive energetic contributions of residues within an interface. This molecular interaction has been well characterized previously both structurally, by x-ray crystallographic analysis, and energetically, by scanning alanine mutagenesis. Through analysis of group and individual maturation and reversion mutations using surface plasmon resonance spectroscopy, we have identified energetically important interfacial residues, determined their cooperative and additive energetic properties, and elucidated the kinetic and thermodynamic bases for molecular evolution in this system. The summation of the binding free energy changes associated with the individual mutations that define this affinity maturation pathway is greater than that of the fully matured variant, even though the affinity gap between the end point variants is relatively large. Two mutations in particular, both located in the complementarity determining region 2 loop of the Vbeta domain, exhibit negative cooperativity.     

Molecular cooperativity of drebrin(1-300) binding and structural remodeling of f-actin

Drebrin A, an actin-binding protein, is a key regulatory element in synaptic plasticity of neuronal dendrites. Understanding how drebrin binds and remodels F-actin is important for a functional analysis of their interactions. Conventionally, molecular models for protein-protein interactions use binding parameters derived from bulk solution measurements with limited spatial resolution, and the inherent assumption of homogeneous binding sites. In the case of actin filaments, their structural and dynamic states—as well as local changes in those states—may influence their binding parameters and interaction cooperativity. Here, we probed the structural remodeling of single actin filaments and the binding cooperativity of DrebrinA_1-300 –F–actin using AFM imaging. We show direct evidence of DrebrinA_1-300-induced cooperative changes in the helical structure of F-actin and observe the binding cooperativity of drebrin to F-actin with nanometer resolution. The data confirm at the in vitro molecular level that variations in the F-actin helical structure can be modulated by cooperative binding of actin-binding proteins.     

Structural,energetic, and functional analysis of a protein-protein interface at distinct stages of affinity maturation

Due to a paucity of studies that synthesize structural, energetic, and functional analyses of a series of protein complexes representing distinct stages in an affinity maturation pathway, the biophysical basis for the molecular evolution of protein-protein interactions is poorly understood. Here, we combine crystal structures and binding-free energies of a series of variant superantigen (SAG)-major histocompatibility complex (MHC) class II complexes exhibiting increasingly higher affinity to reveal that this affinity maturation pathway is controlled largely by two biophysical factors: shape complementarity and buried hydrophobic surface. These factors, however, do not contribute equivalently to the affinity maturation of the interface, as the former dominates the early steps of the maturation process while the latter is responsible for improved binding in later steps. Functional assays reveal how affinity maturation of the SAG-MHC interface corresponds to T cell activation by SAGs.     

Mutational analysis of the affinity maturation of antibody 48G7

The affinity maturation of antibody 48G7 from its germline predecessor 48G7g has been studied at a molecular level through a combination of structural and biochemical means. Each of the nine somatic mutations accumulated during affinity maturation has been assessed for gain or loss of function in both the germline and affinity-matured antibodies. Individual somatic mutations were found to be either positive or neutral in their effects on affinity for hapten JWJ1, with a marked context-dependence for some sites of mutation. In a number of cases significant cooperativity was found between pairs of somatically mutated residues. Interpretation of the structural changes introduced by many of the point mutations has been possible due to the availability of high-resolution crystal structures of 48G7g and 48G7, and mechanisms by which these structural changes may result in enhanced affinity for hapten have been identified. Precise dissection of structure-function relationships in this system provides additional insights into the role of cooperativity in the evolution of antibody affinity. Comparison of 48G7 with previously characterized systems provides a varied view of the structure-function mechanisms by which the humoral immune system produces large increases in affinity.     

Reverse Engineering the Cooperative Machinery of Human Hemoglobin

Hemoglobin transports molecular oxygen from the lungs to all human tissues for cellular respiration. Its α₂β₂ tetrameric assembly undergoes cooperative binding and releasing of oxygen for superior efficiency and responsiveness. Over past decades, hundreds of hemoglobin structures were determined under a wide range of conditions for investigation of molecular mechanism of cooperativity. Based on a joint analysis of hemoglobin structures in the Protein Data Bank (Ren, companion article), here I present a reverse engineering approach to elucidate how two subunits within each dimer reciprocate identical motions that achieves intradimer cooperativity, how ligand-induced structural signals from two subunits are integrated to drive quaternary rotation, and how the structural environment at the oxygen binding sites alter their binding affinity. This mechanical model reveals the intricate design that achieves the cooperative mechanism and has previously been masked by inconsistent structural fluctuations. A number of competing theories on hemoglobin cooperativity and broader protein allostery are reconciled and unified.     

Probing mechanisms of resistance to the tuberculosis drug isoniazid: Conformational changes caused by inhibition of InhA, the enoyl reductase from Mycobacterium tuberculosis

The frontline tuberculosis drug isoniazid (INH) inhibits InhA, the NADH-dependent fatty acid biosynthesis (FAS-II) enoyl reductase from Mycobacterium tuberculosis (MTB), via formation of a covalent adduct with NAD(+) (the INH-NAD adduct). Resistance to INH can be correlated with many mutations in MTB, some of which are localized in the InhA cofactor binding site. While the InhA mutations cause a substantial decrease in the affinity of InhA for NADH, surprisingly the same mutations result in only a small impact on binding of the INH-NAD adduct. Based on the knowledge that InhA interacts in vivo with other components of the FAS-II pathway, we have initiated experiments to determine whether enzyme inhibition results in structural changes that could affect protein-protein interactions involving InhA and how these ligand-induced conformational changes are modulated in the InhA mutants. Significantly, while NADH binding to wild-type InhA is hyperbolic, the InhA mutants bind the cofactor with positive cooperativity, suggesting that the mutations permit access to a second conformational state of the protein. While cross-linking studies indicate that enzyme inhibition causes dissociation of the InhA tetramer into dimers, analytical ultracentrifugation and size exclusion chromatography reveal that ligand binding causes a conformational change in the protein that prevents cross-linking across one of the dimer-dimer interfaces in the InhA tetramer. Interestingly, a similar ligand-induced conformational change is also observed for the InhA mutants, indicating that the mutations modulate communication between the subunits without affecting the two conformational states of the protein that are present.     

X-ray snapshots of the maturation of an antibody response to a protein antigen

The process whereby the immune system generates antibodies of higher affinities during a response to antigen (affinity maturation) is a prototypical example of molecular evolution. Earlier studies have been confined to antibodies specific for small molecules (haptens) rather than for proteins. We compare the structures of four antibodies bound to the same site on hen egg white lysozyme (HEL) at different stages of affinity maturation. These X-ray snapshots reveal that binding is enhanced, not through the formation of additional hydrogen bonds or van der Waals contacts or by an increase in total buried surface, but by burial of increasing amounts of apolar surface at the expense of polar surface, accompanied by improved shape complementarity. The increase in hydrophobic interactions results from highly correlated rearrangements in antibody residues at the interface periphery, adjacent to the central energetic hot spot. This first visualization of the maturation of antibodies to protein provides insights into the evolution of high affinity in other protein-protein interfaces.     

Independent tyrosyl contributions to the CD of Ff gene 5 protein and the distinctive effects of Y41H and Y41F mutants on protein-protein cooperative interactions

The gene 5 protein (g5p) of the Ff virus contains five Tyr, individual mutants of which have now all been characterized by CD spectroscopy. The protein has a dominant tyrosyl 229-nm L(a) CD band that is shown to be approximately the sum of the five individual Tyr contributions. Tyr41 is particularly important in contributing to the high cooperativity with which the g5p binds to ssDNA, and Y41F and Y41H mutants are known to differ in dimer-dimer packing interactions in crystal structures. We compared the solution structures and binding properties of the Y41F and Y41H mutants using CD spectroscopy. Secondary structures of the mutants were similar by CD analyses and close to those derived from the crystal structures. However, there were significant differences in the binding properties of the two mutant proteins. The Y41H protein had an especially low binding affinity and perturbed the spectrum of poly[d(A)] in 2 mM Na(+) much less than did Y41F and the wild-type gene 5 proteins. Moreover, a change in the Tyr 229 nm band, assigned to the perturbation of Tyr34 at the dimer-dimer interface, was absent in titrations with the Y41H mutant under low salt conditions. In contrast, titrations with the Y41H mutant in 50 mM Na(+) exhibited typical CD changes of both the nucleic acid and the Tyr 229-nm band. Thus, protein-protein and g5p-ssDNA interactions appeared to be mutually influenced by ionic strength, indicative of correlated changes in the ssDNA binding and cooperativity loops of the protein or of indirect structural constraints.     

Switch-like Responses of Two Cholesterol Sensors Do Not Require Protein Oligomerization in Membranes

Many cellular processes are sensitive to levels of cholesterol in specific membranes and show a strongly sigmoidal dependence on membrane composition. The sigmoidal responses of the cholesterol sensors involved in these processes could arise from several mechanisms, including positive cooperativity (protein effects) and limited cholesterol accessibility (membrane effects). Here, we describe a sigmoidal response that arises primarily from membrane effects due to sharp changes in the chemical activity of cholesterol. Our models for eukaryotic membrane-bound cholesterol sensors are soluble bacterial toxins that show an identical switch-like specificity for endoplasmic reticulum membrane cholesterol. We show that truncated versions of these toxins fail to form oligomers but still show sigmoidal binding to cholesterol-containing membranes. The nonlinear response emerges because interactions between bilayer lipids control cholesterol accessibility to toxins in a threshold-like fashion. Around these thresholds, the affinity of toxins for membrane cholesterol varies by >100-fold, generating highly cooperative lipid-dependent responses independently of protein-protein interactions. Such lipid-driven cooperativity may control the sensitivity of many cholesterol-dependent processes.     

Affinity mass spectrometry-based approaches for the analysis of protein-protein interaction and complex mixtures of peptide-ligands.

Combined applications of affinity purification procedures and mass-spectrometric analyses (affinity mass spectrometry or affinity-directed mass spectrometry) have gained broad interest in various fields of biological sciences. We have extended these techniques to the purification and analysis of closely related peptides from complex mixtures and to the characterization of binding motifs and relative affinities in protein-protein interactions. The posttranslational modifications in the carboxy-terminal region of porcine brain tubulin are used as an example for the applicability of affinity mass spectrometry in the characterization of complex patterns of related peptides. We also show that affinity mass spectrometry allows the mapping of sequential binding motifs of two interacting proteins. Using the ActA/Mena protein-protein complex as a model system, we show that we can selectively purify Mena-binding peptides from a tryptic digest of ActA. The results from this assay are compared to data sets obtained earlier by classical methods using synthetic peptides and molecular genetic experiments. As a further expansion of affinity mass spectrometry, we have established an internally standardized system that allows comparison of the affinities of related ligands for a given protein. Here the affinities of two peptide ligands for the monoclonal tubulin-specific antibody YL1/2 are determined in terms of half-maximal competition.     

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.     

A likelihood-based index of protein protein binding affinities with application to influenza HA escape from antibodies

In many biological systems, proteins interact with other organic molecules to produce indispensable functions, in which molecular recognition phenomena are essential. Proteins have kept or gained their functions during molecular evolution. Their functions seem to be flexible, and a few amino acid substitutions sometimes cause drastic changes in function. In order to monitor and predict such drastic changes in the early stages in target populations, we need to identify patterns of structural changes during molecular evolution causing decreases or increases in the binding affinity of protein complexes. In previous work, we developed a likelihood-based index to quantify the degree to which a sequence fits a given structure. This index was named the sequence-structure fitness (SSF) and is calculated empirically based on amino acid preferences and pairwise interactions in the structural environment present in template structures. In the present work, we used the SSF to develop an index to measure the binding affinity of protein-protein complexes defined as the log likelihood ratio, contrasting the fitness of the sequences to the structure of the complex and that of the uncomplexed proteins. We applied the developed index to the complexes formed between influenza A hemagglutinin (HA) and four antibodies. The antibody-antigen binding region of HA is under strong selection pressure by the host immune system. Hence, examination of the long-term adaptation of HA to the four antibodies could reveal the strategy of the molecular evolution of HA. Two antibodies cover the HA receptor-binding region, while the other two bind away from the receptor-binding region. By focusing on branches with a significant decline in binding ability, we could detect key amino acid replacements and investigate the mechanism via conditional probabilities. The contrast between the adaptations to the two types of antibodies suggests that the virus adapts to the immune system at the cost of structural change.     

Molecular Basis of eRF3 Recognition by the MLLE Domain of Poly(A)-Binding Protein

PABPC1 (cytosolic poly(A)-binding protein 1) is an RNA-binding protein that binds to the poly(A) tail of mRNAs to promote translation and mRNA turnover. In addition to RNA-binding domains, PABPC1 contains a unique protein-protein interaction domain, MLLE (also known as PABC) that binds regulatory proteins and translation factors that contain a conserved 12 amino acid peptide motif termed PAM2. Eukaryotic Release Factor 3 (eRF3/GSPT1) contains two overlapping PAM2 sequences, which are required for its activity. Here, we determined the crystal structures of the MLLE domain from PABPC1 in complex with the two PAM2 regions of eRF3. The structures reveal a mechanism of cooperativity between the two PAM2 sites that increases the binding affinity but prevents the binding of more than one molecule of eRF3 to PABPC1. Relative to previous structures, the high-resolution crystal structures force a re-evaluation of the PAM2 motif and improve our understanding of the molecular basis of MLLE peptide recognition.     

Computational Assessment of Protein–protein Binding Affinity by Reversely Engineering the Energetics in Protein Complexes

The cellular functions of proteins are maintained by forming diverse complexes. The stability of these complexes is quantified by the measurement of binding affinity, and mutations that alter the binding affinity can cause various diseases such as cancer and diabetes. As a result, accurate estimation of the binding stability and the effects of mutations on changes of binding affinity is a crucial step to understanding the biological functions of proteins and their dysfunctional consequences. It has been hypothesized that the stability of a protein complex is dependent not only on the residues at its binding interface by pairwise interactions but also on all other remaining residues that do not appear at the binding interface. Here, we computationally reconstruct the binding affinity by decomposing it into the contributions of interfacial residues and other non-interfacial residues in a protein complex. We further assume that the contributions of both interfacial and non-interfacial residues to the binding affinity depend on their local structural environments such as solvent-accessible surfaces and secondary structural types. The weights of all corresponding parameters are optimized by Monte-Carlo simulations. After cross-validation against a large-scale dataset, we show that the model not only shows a strong correlation between the absolute values of the experimental and calculated binding affinities, but can also be an effective approach to predict the relative changes of binding affinity from mutations. Moreover, we have found that the optimized weights of many parameters can capture the first-principle chemical and physical features of molecular recognition, therefore reversely engineering the energetics of protein complexes. These results suggest that our method can serve as a useful addition to current computational approaches for predicting binding affinity and understanding the molecular mechanism of protein–protein interactions.     

Structural basis of DNA recognition by p53 tetramers

The tumor-suppressor protein p53 is among the most effective of the cell's natural defenses against cancer. In response to cellular stress, p53 binds as a tetramer to diverse DNA targets containing two decameric half-sites, thereby activating the expression of genes involved in cell-cycle arrest or apoptosis. Here we present high-resolution crystal structures of sequence-specific complexes between the core domain of human p53 and different DNA half-sites. In all structures, four p53 molecules self-assemble on two DNA half-sites to form a tetramer that is a dimer of dimers, stabilized by protein-protein and base-stacking interactions. The protein-DNA interface varies as a function of the specific base sequence in correlation with the measured binding affinities of the complexes. The new data establish a structural framework for understanding the mechanisms of specificity, affinity, and cooperativity of DNA binding by p53 and suggest a model for its regulation by regions outside the sequence-specific DNA binding domain.     

Structural adaptation of extreme halophilic proteins through decrease of conserved hydrophobic contact surface

Background

Disrupting protein-protein interactions by small organic molecules is nowadays a promising strategy employed to block protein targets involved in different pathologies. However, structural changes occurring at the binding interfaces make difficult drug discovery processes using structure-based drug design/virtual screening approaches. Here we focused on two homologous calcium binding proteins, calmodulin and human centrin 2, involved in different cellular functions via protein-protein interactions, and known to undergo important conformational changes upon ligand binding.

Results

In order to find suitable protein conformations of calmodulin and centrin for further structure-based drug design/virtual screening, we performed in silico structural/energetic analysis and molecular docking of terphenyl (a mimicking alpha-helical molecule known to inhibit protein-protein interactions of calmodulin) into X-ray and NMR ensembles of calmodulin and centrin. We employed several scoring methods in order to find the best protein conformations. Our results show that docking on NMR structures of calmodulin and centrin can be very helpful to take into account conformational changes occurring at protein-protein interfaces.

Conclusions

NMR structures of protein-protein complexes nowadays available could efficiently be exploited for further structure-based drug design/virtual screening processes employed to design small molecule inhibitors of protein-protein interactions.     

Computational Assessment of Protein–protein Binding Affinity by Reversely Engineering the Energetics in Protein Complexes

The cellular functions of proteins are maintained by forming diverse complexes. The stability of these complexes is quantified by the measurement of binding affinity, and mutations that alter the binding affinity can cause various diseases such as cancer and diabetes. As a result, accurate estimation of the binding stability and the effects of mutations on changes of binding affinity is a crucial step to understanding the biological functions of proteins and their dysfunctional consequences. It has been hypothesized that the stability of a protein complex is dependent not only on the residues at its binding interface by pairwise interactions but also on all other remaining residues that do not appear at the binding interface. Here, we computationally reconstruct the binding affinity by decomposing it into the contributions of interfacial residues and other non-interfacial residues in a protein complex. We further assume that the contributions of both interfacial and non-interfacial residues to the binding affinity depend on their local structural environments such as solvent-accessible surfaces and secondary structural types. The weights of all corresponding parameters are optimized by Monte-Carlo simulations. After cross-validation against a large-scale dataset, we show that the model not only shows a strong correlation between the absolute values of the experimental and calculated binding affinities, but can also be an effective approach to predict the relative changes of binding affinity from mutations. Moreover, we have found that the optimized weights of many parameters can capture the first-principle chemical and physical features of molecular recognition, therefore reversely engineering the energetics of protein complexes. These results suggest that our method can serve as a useful addition to current computational approaches for predicting binding affinity and understanding the molecular mechanism of protein–protein interactions.     

Enthalpy-entropy compensation and cooperativity as thermodynamic epiphenomena of structural flexibility in ligand-receptor interactions

Ligand binding is a thermodynamically cooperative process in many biochemical systems characterized by the conformational flexibility of the reactants. However, the contribution of conformational entropy to cooperativity of ligation needs to be elucidated. Here, we perform kinetic and thermodynamic analyses on a panel of cycle-mutated peptides, derived from influenza H3 HA(306-319), interacting with wild type and a mutant HLA-DR. We observe that, within a certain range of peptide affinity, this system shows isothermal entropy-enthalpy compensation (iEEC). The incremental increases in conformational entropy measured as disruptive mutations are added in the ligand or receptor are more than sufficient in magnitude to account for the experimentally observed lack of free-energy decrease cooperativity. Beyond this affinity range, compensation is not observed, and therefore, the ability of the residual interactions to form a stable complex decreases in an exponential fashion. Taken together, our results indicate that cooperativity and iEEC constitute the thermodynamic epiphenomena of the structural fluctuation that accompanies ligand-receptor complex formation in flexible systems. Therefore, ligand binding affinity prediction needs to consider how each source of binding energy contributes synergistically to the folding and kinetic stability of the complex in a process based on the trade-off between structural tightening and restraint of conformational mobility.