Blind Prediction of Charged Ligand Binding Affinities in a Model Binding Site

Predicting absolute protein–ligand binding affinities remains a frontier challenge in ligand discovery and design. This becomes more difficult when ionic interactions are involved because of the large opposing solvation and electrostatic attraction energies. In a blind test, we examined whether alchemical free-energy calculations could predict binding affinities of 14 charged and 5 neutral compounds previously untested as ligands for a cavity binding site in cytochrome c peroxidase. In this simplified site, polar and cationic ligands compete with solvent to interact with a buried aspartate. Predictions were tested by calorimetry, spectroscopy, and crystallography. Of the 15 compounds predicted to bind, 13 were experimentally confirmed, while 4 compounds were false negative predictions. Predictions had a root-mean-square error of 1.95 kcal/mol to the experimental affinities, and predicted poses had an average RMSD of 1.7 Å to the crystallographic poses. This test serves as a benchmark for these thermodynamically rigorous calculations at predicting binding affinities for charged compounds and gives insights into the existing sources of error, which are primarily electrostatic interactions inside proteins. Our experiments also provide a useful set of ionic binding affinities in a simplified system for testing new affinity prediction methods.     

Binding Affinities Controlled by Shifting Conformational Equilibria: Opportunities and Limitations

Conformational selection is an established mechanism in molecular recognition. Despite its power to explain binding events, it is hardly used in protein/ligand design to modulate molecular recognition. Here, we explore the opportunities and limitations of design by conformational selection. Using appropriate thermodynamic cycles, our approach predicts the effects of a conformational shift on binding affinity and also allows one to disentangle the effects induced by a conformational shift from other effects influencing the binding affinity. The method is assessed and applied to explain the contribution of a conformational shift on the binding affinity of six ubiquitin mutants showing different conformational shifts in six different complexes.     

Two Amino Acid Residues Confer Different Binding Affinities of Abelson Family Kinase Src Homology 2 Domains for Phosphorylated Cortactin

The closely related Abl family kinases, Arg and Abl, play important non-redundant roles in the regulation of cell morphogenesis and motility. Despite similar N-terminal sequences, Arg and Abl interact with different substrates and binding partners with varying affinities. This selectivity may be due to slight differences in amino acid sequence leading to differential interactions with target proteins. We report that the Arg Src homology (SH) 2 domain binds two specific phosphotyrosines on cortactin, a known Abl/Arg substrate, with over 10-fold higher affinity than the Abl SH2 domain. We show that this significant affinity difference is due to the substitution of arginine 161 and serine 187 in Abl to leucine 207 and threonine 233 in Arg, respectively. We constructed Abl SH2 domains with R161L and S187T mutations alone and in combination and find that these substitutions are sufficient to convert the low affinity Abl SH2 domain to a higher affinity “Arg-like” SH2 domain in binding to a phospho-cortactin peptide. We crystallized the Arg SH2 domain for structural comparison to existing crystal structures of the Abl SH2 domain. We show that these two residues are important determinants of Arg and Abl SH2 domain binding specificity. Finally, we expressed Arg containing an “Abl-like” low affinity mutant Arg SH2 domain (L207R/T233S) and find that this mutant, although properly localized to the cell periphery, does not support wild type levels of cell edge protrusion. Together, these observations indicate that these two amino acid positions confer different binding affinities and cellular functions on the distinct Abl family kinases.     

An In Silico Analysis of the Binding Modes and Binding Affinities of Small Molecule Modulators of PDZ-Peptide Interactions

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known K_i or K_d values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high K_i or K_d moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.     

Rational Design and Cyclization of MIG6 Peptide to Restore its Binding Affinity for ErbB Family Receptor Tyrosine Kinases

The ErbB receptor tyrosine kinase family consists of four members: EGFR/Her1, Her2, Her3 and Her4; they have been involved in a variety of malignant tumors. The mitogen-inducible gene 6 (MIG6) is a natural tumor-suppressor protein that can inactivate ErbB signaling by directly binding to the catalytic domain of ErbB kinases. Here, a peptide segment s2p was stripped from the MIG6 interaction interface with ErbB, which exhibited a very low or no affinity to the four ErbB kinases. Structural dynamics simulations revealed that the linear peptide is highly flexible in unbound state and would incur a considerable entropy penalty upon binding to kinases. In this respect, the s2p peptide was cyclized by rational design of a disulfide bond across its two termini, resulting in a cyclic peptide s2p-c. Integration of computational analysis and experimental assay found that the cyclization can largely constrain s2p conformation, thus minimizing the entropy penalty and restoring the binding affinity of s2p-c to kinases.     

Stable Expression of Recombinant AMPA Receptor Subunits: Binding Affinities and Effects of Allosteric Modulators

Abstract: Homomeric AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid)-type glutamate receptors (GluRs) were stably expressed in kidney cells from cDNAs encoding GluR1 flop, GluR2 flip, GluR2 flop, and GluR3 flop subunits. The recombinant receptors were of the expected size and showed functional properties in whole-cell recording as previously reported. [³H]AMPA binding to all subunits was increased to a similar extent by the chaotropic ion thiocyanate (SCN^?). Significant differences were found in the Scatchard plots, however, which were linear and of high affinity for GluR1 and -3 receptors (K_D values of 33 and 52 nM, respectively) but showed curvature for GluR2 receptors, indicating the presence of two components with distinct affinities. As with brain AMPA receptors, solubilization of GluR2 receptors reduced the number of lower-affinity sites and correspondingly increased the number of higher-affinity sites. The sulfhydryl reagent p-chloromercuriphenylsulfonic acid, which increases binding to brain receptors, produced only minor changes except in the case of GluR2 flip. These results indicate that GluR2, among the subunits examined here, most closely resembles the native AMPA receptors in brain membranes. [³H]AMPA binding was inhibited in a noncompetitive manner by two drugs that change the desensitization kinetics of the AMPA receptor. In agreement with physiological observations, the apparent affinity of cyclothiazide for GluR2 flip (EC₅₀ = 7 µM) was higher than that for receptors made of flop subunits (49–130 µM). In contrast, BDP-37, a member of the benzamide family of drugs, exhibited a lower potency for GluR2 flip (58 µM) than for any of the flop isoforms (18–40 µM). These results predict that the action of centrally active AMPA-receptor modulators varies across brain regions depending on their flip/flop composition.     

Binding specificities and affinities of egf domains for ErbB receptors

ErbB receptor activation is a complex process and is dependent upon the type and number of receptors expressed on a given cell. Previous studies with defined combinations of ErbB receptors expressed in mammalian cells have helped elucidate specific biological responses for many of the recognized gene products that serve as ligands for these receptors. However, no study has examined the binding of these ligands in a defined experimental system. To address this issue, the relative binding affinities of the egf domains of eleven ErbB ligands were measured on six ErbB receptor combinations using a soluble receptor-ligand binding format. The ErbB2/4 heterodimer was shown to bind all ligands tested with moderate to very high affinity. In contrast, ErbB3 showed much more restrictive ligand binding specificity and measurable binding was observed only with heregulin, neuregulin2beta, epiregulin and the synthetic heregulin/egf chimera, biregulin. These studies also revealed that ErbB2 preferentially enhances ligand binding to ErbB3 or ErbB4 and to a lesser degree to ErbB1.     

Heterogeneity of Estrogen Binding Sites in Mouse Mammary Cancer

Abstract

The recent demonstration in our laboratory of at least two specific estrogen binding sites in the rat uterus prompted us to investigate similar heterogeneity of binding sites in a trans-plantable ovarian dependent mouse mammary tumor (MXT-3590). Saturation analysis of cytoplasmic (protamine sulfate or hydroxylapatite exchange assay) or crude nuclear fractions (protamine sulfate precipitated nuclear exchange assay) revealed two binding components: type I which conforms to the classically described estrogen receptor and type II which has a lower affinity for estradiol but a greater capacity than type I sites. Exposure of cytosol to charcoal partially removes bound ³H-estradiol from type II sites but not from type I sites. Type II sites are specific for estrogens and do not translocate from the cytoplasmic to the nuclear compartment. Although Type II sites undergo dissociation on prelabeled sucrose density gradients, they are readily demonstrable by postlabeling sucrose density gradient fractions and hydroxylapatite adsorption. Since the presence of type II sites interferes with the measurement of the estrogen receptor (type I) which may also undergo dissociation on sucrose gradients, we recommended that the technique of postlabeling be used for the sucrose gradient analysis of type I and II sites. In addition, saturation assays should be performed over a wide range of ³H-es-tradiol concentrations (0.1–120 nM) for proper evaluation of both sites. These considerations may contribute to more accurate predictions about the response of breast cancers to endocrine therapies.     

Conformational Heterogeneity of Cyclosporin A in Cyclophilin 18 Binding

The immunosuppressive drug cyclosporin A (CsA) binds to its receptor protein cyclophilin 18 (Cyp18) in two distinct kinetic phases, while the mechanism remains elusive. Stopped-flow measurements coupled with titration and competition experiments were used to investigate the puzzling two-phase process of CsA and Cyp18 interaction. This study leads to the dissection of different conformational fractions of either direct fast binding or slow binding with rate-limiting conformational inter-conversion and the real-time measurement of k_on value (8.34 ± 0.22 x10⁶ M^-1s^-1) in solution. Furthermore, our study indicates that the structure of CsA during dissociation from the protein possesses a distribution of conformations different from those in solution under equilibrium condition.     

Postmortem Degradation Alters Fluorographic Labeling Patterns and Affinities of Benzodiazepine Binding Proteins

To investigate the effect of endogenous proteolysis on the molecular weights of the benzodiazepine binding proteins, brains of trout, chicken, and rat were removed immediately after death and stored at room temperature for various periods of time before they were frozen. Photoaffinity labeling of membranes with [3H]flunitrazepam, followed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and fluorography, revealed proteolytic fragments of 47K in trout, chicken, and rat. The proteolysis set in rapidly after death. Seemingly in parallel with the degradation observed fluorographically, the affinity for [3H]flunitrazepam increased without systematic changes in receptor density. The degradation pattern was not identical to that of the photolabeled trypsinized benzodiazepine binding proteins. The endogenous proteolytic fragments were deglycosylated in two steps. In conclusion, proteolytic effects must be taken into account when interpreting labeling patterns and binding parameters.     

Molecular Modeling of Cytochrome P450 1A1: Enzyme-Substrate Interactions and Substrate Binding Affinities

Abstract

Human cytochrome P450 1A1, which is present in lungs, plays an important role in the metabolic activation of chemical carcinogens, and in particular, is thought to be linked to lung cancer. The mechanism of carcinogenesis is related to the enzyme's ability to oxidize highly toxic compounds, such as polycyclic aromatic hydrocarbons (PAHs), to their carcinogenic derivatives.

In order to better understand P450 1A1 function, a homology model of this enzyme has been constructed. The model has been based on the structure of P450 2C5, the first mammalian P450 to be crystallized. The coordinates of the model have been calculated using a consensus strategy, and the resulting structure has been evaluated with the ProStat and Profiles-3D programs. P450 1A1 substrates, such as benzo[a]pyrene, ethoxyresorufin and methoxyre- sorufin, were then docked into the active site of the model, and key amino acid residues able to interact with the substrate, have been identified. The analysis of enzyme-substrate interactions indicated that hydrophobic interactions are mainly responsible for binding of these substrates in the active site. Moreover, the non-bond enzyme-substrate interaction energy for ethoxyresorufin was lower than that for methoxyresorufin, which is consistent with higher activity of 1A1 towards the former substrate. Key residue Val-382 may play an important role in these interactions. Additionally, we performed binding free energy calculations for the three substrates. The obtained values were similar to those observed experimentally, which suggests that this approach might be useful for prediction of binding constants.     

Electrophoretic Differences in the Capsid Proteins of Simian Virus 40 Plaque Mutants

The structural proteins of three mutants of simian virus 40 (SV40) which differ in plaque size, temperature sensitivity, oncogenicity, host cell restriction, and immunological properties were studied. The polypeptide components of these SV40 strains could not be distinguished by their polyacrylamide gel electrophoretic patterns. When the dissociated virions of two of the mutants were analyzed by the isoelectric focusing technique in a urea gradient, the capsid protein peaks were found to differ significantly in their isoelectric points. The capsid protein of the small-plaque mutant had an isoelectric point of pH 6.51 as compared with pH 6.28 for the large-plaque strain. Isoelectric focusing of the isolated capsid protein revealed three components, a single major subunit and two minor forms. The coat proteins of two of the mutants, small-plaque and minute-plaque strains, were indistinguishable by this technique. The capsid protein peaks obtained by isoelectric focusing were further analyzed by polyacryalmide gel electrophoresis.     

Affinities or Apparent Affinities of the Transport Adenosine Triphosphatase System

The interactions of potassium ions and ATP on transport ATPase activity are discussed, and the interpretation of these interactions is shown to be often ambiguous. Caldwell''s (1968) Physiological Review model is discussed with particular reference to the observed kinetics of sodium: sodium exchange in red cells. Recent experimental work on the properties of the ouabain-sensitive component of potassium efflux from red cells is described. This component of efflux occurs only if either sodium or potassium are present in the external medium, but the effects of external sodium and potassium are not additive. The relation between ouabain-sensitive potassium efflux and the external concentration of sodium (in a potassium-free medium) or of potassium (in low- and high-sodium media) are described. When starved sodium-poor red cells are poisoned with iodoacetamide, loaded with phosphate, and incubated in high-sodium potassium-free media, the ouabain-sensitive efflux of potassium appears to be accompanied by the reversal of the entire ATPase system. About two to three potassium ions leave by the ouabain-sensitive route for each molecule of ATP synthesized. If potassium is present in the external medium, no ouabain-sensitive synthesis of ATP occurs and the ouabain-sensitive efflux of potassium presumably involves the reversal of only the last part of the ATPase system.     

Lactococcus lactis SpOx Spontaneous Mutants: a Family of Oxidative-Stress-Resistant Dairy Strains

Numerous industrial bacteria generate hydrogen peroxide (H₂O₂), which may inhibit the growth of other bacteria in mixed ecosystems. We isolated spontaneous oxidative-stress-resistant (SpOx) Lactococcus lactis mutants by using a natural selection method with milk-adapted strains on dairy culture medium containing H₂O₂. Three SpOx mutants displayed greater H₂O₂ resistance. One of them, SpOx3, demonstrated better behavior in different oxidative-stress situations: (i) higher long-term survival upon aeration in LM17 and milk and (ii) the ability to grow with H₂O₂-producing Lactobacillus delbrueckii subsp. delbrueckii strains. Furthermore, the transit kinetics of the SpOx3 mutant in the digestive tract of a human flora-associated mouse model was not affected.     

Large Genetic Heterogeneity Within Amoebas of the Species Naegleria gruberi and Evolutionary Affinities to the Other Species of the Genus

The allozyme survey was extended to 7 strains of Naegleria gruberi and N. jadini in order to further characterize the genetic structure of these free-living amoebas. As formerly known for several characters the electrophoretic evidence reveals considerable heterogeneity at the genetic level among N. gruberi strains. Moreover, 2 distinct gene pools, that might likely represent natural taxa, are clearly identified. The single strain of N. jadini appears evolutionarily related to 1 group of N. gruberi which is also related to N. a. australiensis.     

Targeted Protein Engineering Provides Insights into Binding Mechanism and Affinities of Bacterial Collagen Adhesins

The collagen-binding bacterial proteins, Ace and Cna, are well characterized on the biochemical and structural level. Despite overall structural similarity, recombinant forms of the Ace and Cna ligand-binding domains exhibit significantly different affinities and binding kinetics for collagen type I (CI) in vitro. In this study, we sought to understand, in submolecular detail, the bases for these differences. Using a structure-based approach, we engineered Cna and Ace variants by altering specific structural elements within the ligand-binding domains. Surface plasmon resonance-based binding analysis demonstrated that mutations that are predicted to alter the orientation of the Ace and Cna N₁ and N₂ subdomains significantly affect the interaction between the MSCRAMM (microbial surface components recognizing adhesive matrix molecule) and CI in vitro, including affinity, association/dissociation rates and binding ratio. Moreover, we utilized this information to engineer an Ace variant with an 11,000-fold higher CI affinity than the parent protein. Finally, we noted that several engineered proteins that exhibited a weak interaction with CI recognized more sites on CI, suggesting an inverse correlation between affinity and specificity.     

L-Glutamate Receptor Heterogeneity: Labeling of Distinct Receptor Sub-Types Using Radioligand Binding Techniques

Abstract

This study demonstrates (1) that L-[³H]glutamate labels 3 distinct binding sites (types A1, A2 and A4) in isolated rat brain membranes and (2) that only the N-methyl-aspartate (A1) and quisqualate (A2) receptor classes are associated with the postsynaptic density (PSD). L-[³H]glutamate bound to PSDs with K_d 339 nM and B_max 6.1 pmol/mg protein. These sites were resolved into 2 distinct sub-types on the basis of inhibition studies. N-Methyl-aspartate maximally inhibited 57% of PSD-located L-glutamate binding sites (the A1 site) and quisqualate 43% (the A2 site); the effects of both substances were additive. The ligand selectivities of these 2 sites indicated their identity with the N-methyl-D-aspartate and quisqualate receptor classes defined electrophysiologically. The Cl^?-dependent population of L-glutamate binding sites (the A4 site) which predominates in synaptic membranes was absent from PSDs.     

Quantitation of the Effect of ErbB2 on Epidermal Growth Factor Receptor Binding and Dimerization

The epidermal growth factor (EGF) receptor is a member of the ErbB family of receptors that also includes ErbB2, ErbB3, and ErbB4. These receptors form homo- and heterodimers in response to ligand with ErbB2 being the preferred dimerization partner. Here we use (125)I-EGF binding to quantitate the interaction of the EGF receptor with ErbB2. We show that the EGFR/ErbB2 heterodimer binds EGF with a 7-fold higher affinity than the EGFR homodimer. Because it cannot bind a second ligand, the EGFR/ErbB2 heterodimer is not subject to ligand-induced dissociation caused by the negatively cooperative binding of EGF to the second site on the EGFR homodimer. This increases the stability of the heterodimer relative to the homodimer and is associated with enhanced and prolonged EGF receptor autophosphorylation. These effects are independent of the kinase activity of ErbB2 but require back-to-back dimerization of the EGF receptor with ErbB2. Back-to-back dimerization is also required for phosphorylation of ErbB2. These findings provide a molecular explanation for the apparent preference of the EGF receptor for dimerizing with ErbB2 and suggest that the phosphorylation of ErbB2 occurs largely in the context of the EGFR/ErbB2 heterodimer, rather than through lateral phosphorylation of isolated ErbB2 subunits.