Development of a Model Protein Interaction Pair as a Benchmarking Tool for the Quantitative Analysis of 2-Site Protein-Protein Interactions

A significant challenge in the molecular interaction field is to accurately determine the stoichiometry and stepwise binding affinity constants for macromolecules having >1 binding site. The mission of the Molecular Interactions Research Group (MIRG) of the Association of Biomolecular Resource Facilities (ABRF) is to show how biophysical technologies are used to quantitatively characterize molecular interactions, and to educate the ABRF members and scientific community on the utility and limitations of core technologies [such as biosensor, microcalorimetry, or analytic ultracentrifugation (AUC)]. In the present work, the MIRG has developed a robust model protein interaction pair consisting of a bivalent variant of the Bacillus amyloliquefaciens extracellular RNase barnase and a variant of its natural monovalent intracellular inhibitor protein barstar. It is demonstrated that this system can serve as a benchmarking tool for the quantitative analysis of 2-site protein-protein interactions. The protein interaction pair enables determination of precise binding constants for the barstar protein binding to 2 distinct sites on the bivalent barnase binding partner (termed binase), where the 2 binding sites were engineered to possess affinities that differed by 2 orders of magnitude. Multiple MIRG laboratories characterized the interaction using isothermal titration calorimetry (ITC), AUC, and surface plasmon resonance (SPR) methods to evaluate the feasibility of the system as a benchmarking model. Although general agreement was seen for the binding constants measured using solution-based ITC and AUC approaches, weaker affinity was seen for surface-based method SPR, with protein immobilization likely affecting affinity. An analysis of the results from multiple MIRG laboratories suggests that the bivalent barnase-barstar system is a suitable model for benchmarking new approaches for the quantitative characterization of complex biomolecular interactions.     

Thermodynamics of DNA binding and distortion by the hyperthermophile chromatin protein Sac7d

Sac7d is a hyperthermophile chromatin protein which binds non-specifically to the minor groove of duplex DNA and induces a sharp kink of 66 degrees with intercalation of valine and methionine side-chains. We have utilized the thermal stability of Sac7d and the lack of sequence specificity to define the thermodynamics of DNA binding over a wide temperature range. The binding affinity for poly(dGdC) was moderate at 25 degrees C (Ka = 3.5(+/-1.6) x 10(6) M(-1)) and increased by nearly an order of magnitude from 10 degrees C to 80 degrees C. The enthalpy of binding was unfavorable at 25 degrees C, and decreased linearly from 5 degrees C to 60 degrees C. A positive binding heat at 25 degrees C is attributed in part to the energy of distorting DNA, and ensures that the temperature of maximal binding affinity (75.1+/-5.6 degrees C) is near the growth temperature of Sulfolobus acidocaldarius. Truncation of the two intercalating residues to alanine led to a decreased ability to bend and unwind DNA at 25 degrees C with a small decrease in binding affinity. The energy gained from intercalation is slightly greater than the free energy penalty of bending duplex DNA. Surprisingly, reduced distortion from the double alanine substitution did not lead to a significant decrease in the heat of binding at 25 degrees C. In addition, an anomalous positive DeltaCp of binding was observed for the double alanine mutant protein which could not be explained by the change in polar and apolar accessible surface areas. Both the larger than expected binding enthalpy and the positive heat capacity can be explained by a temperature dependent structural transition in the protein-DNA complex with a Tm of 15-20 degrees C and a DeltaH of 15 kcal/mol. Data are discussed which indicate that the endothermic transition in the complex is consistent with DNA distortion.     

A binding free energy hot spot in the ankyrin repeat protein GABPbeta mediated protein-protein interaction

The frequently observed ankyrin repeat motif represents a structural scaffold evolved for mediating protein-protein interactions. As such, these repeats modulate a diverse range of cellular functions. We thermodynamically characterized the heterodimeric GA-binding protein (GABP) alphabeta complex and focused specifically on the interaction mediated by the ankyrin repeat domain of the GABPbeta. Our isothermal titration calorimetric analysis of the interaction between the GABP subunits determined an association constant (K(A)) of 6.0 x 10(8) M(-1) and that the association is favorably driven by a significant change in enthalpy (DeltaH) and a minor change in entropy (-TDeltaS). A total of 16 GABPbeta interface residues were chosen for alanine scanning mutagenesis. The calorimetrically measured differences in the free energy of binding were compared to computationally calculated values resulting in a correlation coefficient r = 0.71. We identified three spatially contiguous hydrophobic and aromatic residues that form a binding free energy hot spot (DeltaDeltaG > 2.0 kcal/mol). One residue provides structural support to the hot spot residues. Three non-hot spot residues are intermediate contributors (DeltaDeltaG approximately 1.0 kcal/mol) and create a canopy-like structure over the hot spot residues to possibly occlude solvent and orientate the subunits. The remaining interface residues are located peripherally and have weak contributions. Finally, our mutational analysis revealed a significant entropy-enthalpy compensation for this interaction.     

Predicting Ligand Binding Affinity with Alchemical Free Energy Methods in a Polar Model Binding Site

We present a combined experimental and modeling study of organic ligand molecules binding to a slightly polar engineered cavity site in T4 lysozyme (L99A/M102Q). For modeling, we computed alchemical absolute binding free energies. These were blind tests performed prospectively on 13 diverse, previously untested candidate ligand molecules. We predicted that eight compounds would bind to the cavity and five would not; 11 of 13 predictions were correct at this level. The RMS error to the measurable absolute binding energies was 1.8 kcal/mol. In addition, we computed “relative” binding free energies for six phenol derivatives starting from two known ligands: phenol and catechol. The average RMS error in the relative free energy prediction was 2.5 kcal/mol (phenol) and 1.1 kcal/mol (catechol). To understand these results at atomic resolution, we obtained x-ray co-complex structures for nine of the diverse ligands and for all six phenol analogs. The average RMSD of the predicted pose to the experiment was 2.0 Å (diverse set), 1.8 Å (phenol-derived predictions), and 1.2 Å (catechol-derived predictions). We found that predicting accurate affinities and rank-orderings required near-native starting orientations of the ligand in the binding site. Unanticipated binding modes, multiple ligand binding, and protein conformational change all proved challenging for the free energy methods. We believe that these results can help guide future improvements in physics-based absolute binding free energy methods.     

Structural and binding studies of the three-metal center in two mycobacterial PPM Ser/Thr protein phosphatases

Phospho-Ser/Thr protein phosphatases (PPs) are dinuclear metalloenzymes classed into two large families, PPP and PPM, on the basis of sequence similarity and metal ion dependence. The archetype of the PPM family is the α isoform of human PP2C (PP2Cα), which folds into an α/β domain similar to those of PPP enzymes. The recent structural studies of three bacterial PPM phosphatases, Mycobacterium tuberculosis MtPstP, Mycobacterium smegmatis MspP, and Streptococcus agalactiae STP, confirmed the conservation of the overall fold and dinuclear metal center in the family, but surprisingly revealed the presence of a third conserved metal-binding site in the active site. To gain insight into the roles of the three-metal center in bacterial enzymes, we report structural and metal-binding studies of MtPstP and MspP. The structure of MtPstP in a new trigonal crystal form revealed a fully active enzyme with the canonical dinuclear metal center but without the third metal ion bound to the catalytic site. The absence of metal correlates with a partially unstructured flap segment, indicating that the third manganese ion contributes to reposition the flap, but is dispensable for catalysis. Studies of metal binding to MspP using isothermal titration calorimetry revealed that the three Mn²⁺-binding sites display distinct affinities, with dissociation constants in the nano- and micromolar range for the two catalytic metal ions and a significantly lower affinity for the third metal-binding site. In agreement, the structure of inactive MspP at acidic pH was determined at atomic resolution and shown to lack the third metal ion in the active site. Structural comparisons of all bacterial phosphatases revealed positional variations in the third metal-binding site that are correlated with the presence of bound substrate and the conformation of the flap segment, supporting a role of this metal ion in assisting enzyme-substrate interactions.     

Predicting absolute ligand binding free energies to a simple model site

A central challenge in structure-based ligand design is the accurate prediction of binding free energies. Here we apply alchemical free energy calculations in explicit solvent to predict ligand binding in a model cavity in T4 lysozyme. Even in this simple site, there are challenges. We made systematic improvements, beginning with single poses from docking, then including multiple poses, additional protein conformational changes, and using an improved charge model. Computed absolute binding free energies had an RMS error of 1.9 kcal/mol relative to previously determined experimental values. In blind prospective tests, the methods correctly discriminated between several true ligands and decoys in a set of putative binders identified by docking. In these prospective tests, the RMS error in predicted binding free energies relative to those subsequently determined experimentally was only 0.6 kcal/mol. X-ray crystal structures of the new ligands bound in the cavity corresponded closely to predictions from the free energy calculations, but sometimes differed from those predicted by docking. Finally, we examined the impact of holding the protein rigid, as in docking, with a view to learning how approximations made in docking affect accuracy and how they may be improved.     

Biophysical characterization of cyclic nucleotide phosphodiesterases

We have compared selected biophysical properties of three phosphodiesterases, from Arabidopsis thaliana, Saccharomyces cerevisiae, and Escherichia coli. All of them belong to a recently identified family of cyclic nucleotide phosphodiesterases. Experiments elucidating folding stability, protein fluorescence, oligomerization behavior, and the effects of substrates were conducted, revealing differences between the plant and the yeast protein. According to CD spectroscopy, the latter protein exhibits an (alpha + beta) fold rather than an (alpha/beta) fold as found with CPDase (A. thaliana). The redox-dependent structural reorganization recently found for the plant protein by X-ray crystallography could not be detected by CD spectroscopy due to its only marginal effect on the total percentage of helical content. However, in the present study a redox-dependent effect was also observed for the yeast CPDase. The enzymatic activity of wild type CPDase (A. thaliana) as well as of four mutants were characterized by isothermal titration calorimetry and the results prove the requirement of all four residues of the previously identified tandem signature motif for the catalytic function. Within the comparison of the three proteins in this study, the PDase Homolog/RNA ligase (E. coli) shares more similarities with the plant than with the yeast protein.     

Arabidopsis cytosolic acyl-CoA-binding proteins ACBP4, ACBP5 and ACBP6 have overlapping but distinct roles in seed development

Eukaryotic cytosolic ACBPs (acyl-CoA-binding proteins) bind acyl-CoA esters and maintain a cytosolic acyl-CoA pool, but the thermodynamics of their protein–lipid interactions and physiological relevance in plants are not well understood. Arabidopsis has three cytosolic ACBPs which have been identified as AtACBP4, AtACBP5 and AtACBP6, and microarray data indicated that all of them are expressed in seeds; AtACBP4 is expressed in early embryogenesis, whereas AtACBP5 is expressed later. ITC (isothermal titration calorimetry) in combination with transgenic Arabidopsis lines were used to investigate the roles of these three ACBPs from Arabidopsis thaliana. The dissociation constants, stoichiometry and enthalpy change of AtACBP interactions with various acyl-CoA esters were determined using ITC. Strong binding of recombinant (r) AtACBP6 with long-chain acyl-CoA (C16- to C18-CoA) esters was observed with dissociation constants in the nanomolar range. However, the affinity of rAtACBP4 and rAtACBP5 to these acyl-CoA esters was much weaker (dissociation constants in the micromolar range), suggesting that they interact with acyl-CoA esters differently from rAtACBP6. When transgenic Arabidopsis expressing AtACBP6pro::GUS was generated, strong GUS (β-glucuronidase) expression in cotyledonary-staged embryos and seedlings prompted us to measure the acyl-CoA contents of the acbp6 mutant. This mutant accumulated higher levels of C18:1-CoA and C18:1- and C18:2-CoAs in cotyledonary-staged embryos and seedlings, respectively, in comparison with the wild type. The acbp4acbp5acbp6 mutant showed the lightest seed weight and highest sensitivity to abscisic acid during germination, suggesting their physiological functions in seeds.