Allosteric effects of monoclonal antibodies on human growth hormone

We have previously shown that a monoclonal antibody (MAb) recognizing the human growth hormone (hGH) antigenic domain left exposed after binding to lactogenic receptors enhanced hGH binding probably through allosteric effects on the hormone binding site. Since receptors displaying different specificities would not recognize exactly the same hGH region, we explored whether some of our MAb could affect hGH binding to somatogenic receptors from rabbit liver and to human liver hGH-specific receptors.The effect of MAbAE5, AC8 and F11 on hGH binding was measured by determining the formation of¹²⁵I-MAb:hGH:receptor complexes using two different experimental approaches. Results from procedure A, which involved the previous binding of the hormone to microsomes before adding¹²⁵I-MAb, indicated that the hGH domain defined by epitopes AE5, AC8 and F11 is uncovered in the various hormone:receptor complexes.Procedure B was devised to reveal any alteration in the hGH molecule induced by the MAb. In this case preformed¹²⁵I-MAb:hGH complexes were added to microsomes. Data showed that¹²⁵I-MAb AE5:hGH complexes bound better to the various receptors than¹²⁵I-MAb AE5 to hGH:receptor complexes. On the contrary, hGH previously bound to¹²⁵I-MAb AC8 or¹²⁵I-MAb F11 was less recognized by the receptors than the free hormone. Furthermore, binding of MAb AE5 or MAb F11 to hGH 20 K (a natural hGH variant lacking residues 32–46) also enhanced its affinity to the various receptors whereas MAb AC8 did not inhibit hGH 20 K binding.Results indicated that MAb recognizing the hGH antigenic area that remains unmasked after binding to different membrane-bound receptors are able to affect hormone binding site. MAb would induce either positive or negative allosteric changes in the hormone region involved in its binding to lactogenic, somatogenic and hGH-specific receptors.     

Positive cooperative effects between receptors induced by an anti-human growth hormone allosteric monoclonal antibody

Monoclonal antibodies (MAb) anti-human growth hormone (hGH) termed MAb AE5, AC8 and F11 recognize a cluster of epitopes left exposed after hormone binding to receptors. Since these MAb were able to produce either positive (MAb AE5) or negative (MAb AC8 and F11) allosteric effects on hGH binding, the purpose of this work was to further characterize MAb behavior. Results indicated a straight correlation between MAb allosteric effects and affinity constant values for binding of different hGH:MAb complexes to lactogenic receptors from rat liver. Affinity of hGH:MAb AE5 as well as hGH:Fab AE5 complexes enhanced proportionally to the fraction of occupied receptors and Hill coefficients higher than 1 were obtained, suggesting the induction of positive cooperative effects between membrane-bound receptors. On the other hand, hGH:MAb AC8 and hGH:MAb F11 complexes binding affinity to lactogenic sites could not be related to receptor occupancy degree. It is proposed that binding of hGH:MAb AE5 complexes to receptors would elicit a conformational change on adjacent receptor molecules leading to an increase of their affinity to bind subsequent hGH:MAb AE5 complexes.     

Multi-resolution contour-based fitting of macromolecular structures

A novel contour-based matching criterion is presented for the quantitative docking of high-resolution structures of components into low-resolution maps of macromolecular complexes. The proposed Laplacian filter is combined with a six-dimensional search using fast Fourier transforms to rapidly scan the rigid-body degrees of freedom of a probe molecule relative to a fixed target density map. A comparison of the docking performance with the standard cross-correlation criterion demonstrates that contour matching with the Laplacian filter significantly extends the viable resolution range of correlation-based fitting to resolutions as low as 30 A. The gain in docking precision at medium to low resolution (15-30 A) is critical for image reconstructions from electron microscopy (EM). The new algorithm enables for the first time the reliable docking of smaller molecular components into EM densities of large biomolecular assemblies at such low resolutions. As an example of the practical effectiveness of contour-based fitting, a new pseudo-atomic model of a microtubule was constructed from a 20 A resolution EM map and from atomic structures of alpha and beta tubulin subunits.     

Alternative views of functional protein binding epitopes obtained by combinatorial shotgun scanning mutagenesis

Combinatorial shotgun scanning mutagenesis was used to analyze two large, related protein binding sites to assess the specificity and importance of individual side chain contributions to binding affinity. The strategy allowed for cost-effective generation of a plethora of functional data. The ease of the technology promoted comprehensive investigations, in which the classic alanine-scanning approach was expanded with two additional strategies, serine- and homolog-scanning. Binding of human growth hormone (hGH) to the hGH receptor served as the model system. The entire high affinity receptor-binding sites (site 1) of wild-type hGH (hGH_wt) and of an affinity-improved variant (hGH_v) were investigated and the results were compared. The contributions that 35 residue positions make to binding were assessed on each hormone molecule by both serine- and homolog-scanning. The hormone molecules were displayed on the surfaces of bacteriophage, and the 35 positions were randomized simultaneously to allow equal starting frequencies of the wild-type residue and either serine or a homologous mutation in separate libraries. Functional selections for binding to the hGH receptor shifted the relative wild-type/mutant frequencies at each position to an extent characteristic of the functional importance of the side chain. Functional epitope maps were created and compared to previous maps obtained by alanine-scanning. Comparisons between the different scans provide insights into the affinity maturation process that produced hGH_v. The serine and homolog-scanning results expand upon and complement the alanine-scanning results and provide additional data on the robustness of the high affinity receptor-binding site of hGH.     

Design, docking, and evaluation of multiple libraries against multiple targets

We present a general approach to the design, docking, and virtual screening of multiple combinatorial libraries against a family of proteins. The method consists of three main stages: docking the scaffold, selecting the best substituents at each site of diversity, and comparing the resultant molecules within and between the libraries. The core "divide-and-conquer" algorithm for side-chain selection, developed from an earlier version (Sun et al., J Comp Aided Mol Design 1998;12:597-604), provides a way to explore large lists of substituents with linear rather than combinatorial time dependence. We have applied our method to three combinatorial libraries and three serine proteases: trypsin, chymotrypsin, and elastase. We show that the scaffold docking procedure, in conjunction with a novel vector-based orientation filter, reproduces crystallographic binding modes. In addition, the free-energy-based scoring procedure (Zou et al., J Am Chem Soc 1999;121:8033-8043) is able to reproduce experimental binding data for P1 mutants of macromolecular protease inhibitors. Finally, we show that our method discriminates between a peptide library and virtual libraries built on benzodiazepine and tetrahydroisoquinolinone scaffolds. Implications of the docking results for library design are explored.     

Conformational and molecular interaction studies of glucagon-like peptide-2 with its N-terminal extracellular receptor domain

Glucagon-like peptide-2 (GLP-2) is a therapeutic target used in the treatment of short bowel syndrome. In this paper, we present the three dimensional solution structure of GLP-2 peptide determined using nuclear magnetic resonance (NMR) and molecular modelling. The GLP-2 adopts an α-helical conformation similar to that of secretin family of hormones. In order to understand the molecular details governing the ligand binding and receptor activation, macromolecular docking studies were performed between the N-terminal extracellular domain of GLP-2 receptor and the GLP-2 hormone using a data driven docking program.     

Flexible docking using tabu search and an empirical estimate of binding affinity

This article describes the implementation of a new docking approach. The method uses a Tabu search methodology to dock flexibly ligand molecules into rigid receptor structures. It uses an empirical objective function with a small number of physically based terms derived from fitting experimental binding affinities for crystallographic complexes. This means that docking energies produced by the searching algorithm provide direct estimates of the binding affinities of the ligands. The method has been tested on 50 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. All water molecules are removed from the structures and ligand molecules are minimized in vacuo before docking. The lowest energy geometry produced by the docking protocol is within 1.5 Å root-mean square of the experimental binding mode for 86% of the complexes. The lowest energies produced by the docking are in fair agreement with the known free energies of binding for the ligands. Proteins 33:367–382, 1998. © 1998 Wiley-Liss, Inc.     

Increased site 1 affinity improves biopotency of porcine growth hormone. Evidence against diffusion dependent receptor dimerization

Based on phage display optimization studies with human growth hormone (GH), it is thought that the biopotency of GH cannot be increased. This is proposed to be a result of the affinity of the first receptor for hormone far exceeding that which is required to trap the hormone long enough to allow diffusion of the second receptor to form the ternary complex, which initiates signaling. We report here that despite similar site 1 kinetics to the hGH/hGH receptor interaction, the potency of porcine GH for its receptor can be increased up to 5-fold by substituting hGH residues involved in site 1 binding into pGH. Based on extensive mutations and BIAcore studies, we show that the higher potency and site 1 affinity of hGH for the pGHR is primarily a result of a decreased off-rate associated with residues in the extended loop between helices 1 and 2 that interact with the two key tryptophans Trp104 and Trp169 in the receptor binding hot spot. Our mutagenic analysis has also identified a second determinant (Lys165), which in addition to His169, restricts the ability of non-primate hormones to activate hGH receptor. The increased biopotency of GH that we observe can be explained by a model for GH receptor activation where subunit alignment is critical for effective signaling.     

A flexible docking scheme to explore the binding selectivity of PDZ domains

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTA LIGAND , we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density‐95/Dlg/ZO‐1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately.     

Comparison of the intermediate complexes of human growth hormone bound to the human growth hormone and prolactin receptors.

The crystal structures of complexes of human growth hormone (hGH) with the growth hormone and prolactin receptors (hGHR and hPRLR, respectively), together with the mutational data available for these systems, suggest that an extraordinary combination of conformational adaptability, together with finely tuned specificity, governs the molecular recognition processes operative in these systems. On the one hand, in the active 1:2 ligand-receptor complexes, 2 copies of the same receptor use the identical set of binding determinants to recognize topographically different surfaces on the hormone. On the other hand, comparing the 1:1 hGH-hGHR and hGH-hPRLR complexes, 2 distinct receptors use this same set of binding determinants to interact with the identical binding site on the ligand, even though few residues among the binding determinants are conserved. The structural evidence demonstrates that this versatility is accomplished by local conformational flexibility of the binding loops, allowing adaptation to different binding environments, together with rigid-body movements of the receptor domains, necessary for the creation of specific interactions with the same binding site.     

Binding mode of ecdysone agonists to the receptor: comparative modeling and docking studies

Three-dimensional structure models of the ligand-binding domain of the ecdysone receptor of Heliothis virescens were built by the homology modeling technique from the crystal structures of nuclear receptors. Two models were created based both on known ligand-binding domain structures of the receptors with the highest sequence identity to the ecdysone receptor, and on those of steroid hormone receptors. The latter model, which was found to have better stereochemical quality and be in good agreement with the binding of the steroidal framework of the endogenous agonist 20-hydroxyecdysone, was used for docking studies. The docking of 20-hydroxyecdysone to the receptor model revealed that the ligand molecule can interact with the receptor in a similar manner to other steroid hormone-receptor complexes. The docking of a dibenzoylhydrazine agonist, chromafenozide, was performed based on the correspondences between the molecule and 20-dydroxyecdysone expected by molecular comparison. The interactions of the ligands with the receptor in the complexes modeled were investigated and found to be consistent with known structure-activity relationships.     

Identification of residues outside the two binding sites that are critical for activation of the lactogenic activity of human growth hormone

Human growth hormone (hGH) binds lactogenic or somatotrophic receptors, creating active heterotrimeric complexes. Comparison of hGH structures, either free or bound to a single lactogenic or somatotrophic receptor, shows binding is associated with structural changes. Changes in hGH structure are unique when binding either lactogenic or somatotrophic receptors and they influence the spatial arrangement of residues constituting the second receptor-binding site. Using site-directed mutagenesis, we identified a contiguous set of largely hydrophobic residues that forms a motif communicating between the two receptor-binding sites of hGH. The residues are external to the receptor-binding epitopes and were identified when their mutation reduced site 2 function without changing site 1 function. The motif includes Phe44, Leu93, Tyr160, Leu163, and Tyr164, located in two hydrophobic clusters between the receptor-binding sites. Their mutation to Glu disrupts hydrophobic interactions and reduces lactogenic activity between 4.7- and 85-fold with little effect on somatotrophic activity or spectroscopic properties. These differential effects indicate that loss of lactogenic activity is not a result of global mis-folding. We propose the loss of lactogenic activity results from disruption of specific hydrophobic clusters that disables the site 1 binding-induced structuring of the second receptor-binding site.     

Receptor binding properties of four-helix-bundle growth factors deduced from electrostatic analysis

Hormones of the hematopoietin class mediate signal transduction by binding to specific transmembrane receptors. Structural data show that the human growth hormone (hGH) forms a complex with a homodimeric receptor and that hGH is a member of a class of hematopoietins possessing an antiparallel 4-α-helix bundle fold. Mutagenesis experiments suggest that electrostatic interactions may have an important influence on hormone-receptor recognition. In order to examine the specificity of hormone-receptor complexation, an analysis was made of the electrostatic potentials of hGH, interleukin-2 (IL-2), interleukin-4 (IL-4), granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage colony-stimulating factor (GM-CSF), and the hGH and IL-4 receptors. The binding surfaces of hGH and its receptor, and of IL-4 and its receptor, show complementary electrostatic potentials. The potentials of the hGH and its receptor display approximately 2-fold rotational symmetry because the receptor subunits are identical. In contrast, the potentials of GM-CSF and IL-2 lack such symmetry, consistent with their known high affinity for hetero-oligomeric receptors. Analysis of the electrostatic potentials supports a recently proposed hetero-oligomeric model for a high-affinity IL-4 receptor and suggests a possible new receptor binding mode for G-CSF; it also provides valuable information for guiding structural and mutagenesis studies of signal-transducing proteins and their receptors.     

Distance-restrained docking of rifampicin and rifamycin SV to RNA polymerase using systematic FRET measurements: developing benchmarks of model quality and reliability

We are developing distance-restrained docking strategies for modeling macromolecular complexes that combine available high-resolution structures of the components and intercomponent distance restraints derived from systematic fluorescence resonance energy transfer (FRET) measurements. In this article, we consider the problem of docking small-molecule ligands within macromolecular complexes. Using simulated FRET data, we have generated a series of benchmarks that permit estimation of model accuracy based on the quantity and quality of FRET-derived distance restraints, including the number, random error, systematic error, distance distribution, and radial distribution of FRET-derived distance restraints. We find that expected model accuracy is 10 A or better for models based on: i), > or =20 restraints with up to 15% random error and no systematic error, or ii), > or =20 restraints with up to 15% random error, up to 10% systematic error, and a symmetric radial distribution of restraints. Model accuracies can be improved to 5 A or better by increasing the number of restraints to > or =40 and/or by optimizing the distance distribution of restraints. Using experimental FRET data, we have defined the positions of the binding sites within bacterial RNA polymerase of the small-molecule inhibitors rifampicin (Rif) and rifamycin SV (Rif SV). The inferred binding sites for Rif and Rif SV were located with accuracies of, respectively, 7 and 10 A relative to the crystallographically defined binding site for Rif. These accuracies agree with expectations from the benchmark simulations and suffice to indicate that the binding sites for Rif and Rif SV are located within the RNA polymerase active-center cleft, overlapping the binding site for the RNA-DNA hybrid.     

Homology modeling of rabbit prolactin hormone complexed with its receptor

The pituitary hormone prolactin (prl) is implicated in a number of biological functions, especially lactation, which is mediated through specific lactogenic receptors (PrlR). Human growth hormone (hGH) is also a pituitary hormone responsible for linear growth. While the growth hormone receptor (hGHR) binds only hGH, hPrlR can interact with both hGH and hPrl. Using structural information from the human growth hormone (hGH)/receptor (hGHR) complex, we modeled by homology a complex between rabbit prolactin hormone (rbPrl) and its receptor (rbPrlR). While the somatogenic hormone/somatogenic receptor (hGH/hGHR) and somatogenic hormone/lactogenic receptor (hGH/hPrlR) interactions are now known and well studied, here we propose a model for the interaction of the lactogenic hormone with its receptor (rbPrl/rbPrlR), and compare these three kinds of ligand/receptor interaction. We identified residues contributing to the active site and tested the potential dimerization of the receptor. Biochemical studies and information deduced from the modeled complex do not exclude a homodimeric form but point to a functional heterodimeric complex. Proteins 27: 459–468, 1997. © 1997 Wiley-Liss, Inc.     

Selecting high-affinity binding proteins by monovalent phage display

Variants of human growth hormone (hGH) with increased affinity and specificity for the hGH receptor were isolated using an improved phage display system. Nearly one million random mutants of hGH were generated at 12 sites previously shown to modulate binding to the hGH receptor or human prolactin (hPRL) receptor. The mutant hormones were displayed in a monovalent fashion from filamentous phage particles as fusions to the gene III product of M13 packaged within each particle. After three to six cycles of enrichment for hGH-phage particles that bound to hGH receptor beads, we isolated hGH mutants that exhibited consensus binding sequences for the hGH receptor. Residues previously identified as important for hGH receptor binding by alanine-scanning mutagenesis were more highly conserved by this selection method. However, other residues nearby were not optimal, and by mutating them, hormone variants having greater affinity and selectivity for the hGH receptor were isolated. This approach should be useful for those who wish to modify and understand the energetics of protein-ligand interfaces.     

Model of three-dimensional structure of vitamin D receptor and its binding mechanism with 1alpha,25-dihydroxyvitamin D(3).

Comparative modeling of the vitamin D receptor three-dimensional structure and computational docking of 1alpha,25-dihydroxyvitamin D(3) into the putative binding pocket of the two deletion mutant receptors: (207-423) and (120-422, Delta [164-207]) are reported and evaluated in the context of extensive mutagenic analysis and crystal structure of holo hVDR deletion protein published recently. The obtained molecular model agrees well with the experimentally determined structure. Six different conformers of 1alpha,25-dihydroxyvitamin D(3) were used to study flexible docking to the receptor. On the basis of values of conformational energy of various complexes and their consistency with functional activity, it appears that 1alpha,25-dihydroxyvitamin D(3) binds the receptor in its 6-s-trans form. The two lowest energy complexes obtained from docking the hormone into the deletion protein (207-423) differ in conformation of ring A and orientation of the ligand molecule in the VDR pocket. 1alpha,25-Dihydroxyvitamin D(3) possessing the A-ring conformation with axially oriented 1alpha-hydroxy group binds receptor with its 25-hydroxy substituent oriented toward the center of the receptor cavity, whereas ligand possessing equatorial conformation of 1alpha-hydroxy enters the pocket with A ring directed inward. The latter conformation and orientation of the ligand is consistent with the crystal structure of hVDR deletion mutant (118-425, Delta [165-215]). The lattice model of rVDR (120-422, Delta [164-207]) shows excellent agreement with the crystal structure of the hVDR mutant. The complex obtained from docking the hormone into the receptor has lower energy than complexes for which homology modeling was used. Thus, a simple model of vitamin D receptor with the first two helices deleted can be potentially useful for designing a general structure of ligand, whereas the advanced lattice model is suitable for examining binding sites in the pocket.     

Two EGF molecules contribute additively to stabilization of the EGFR dimer.

Receptor dimerization is generally considered to be the primary signaling event upon binding of a growth factor to its receptor at the cell surface. Little, however, is known about the precise molecular details of ligand-induced receptor dimerization, except for studies of the human growth hormone (hGH) receptor. We have analyzed the binding of epidermal growth factor (EGF) to the extracellular domain of its receptor (sEGFR) using titration calorimetry, and the resulting dimerization of sEGFR using small-angle X-ray scattering. EGF induces the quantitative formation of sEGFR dimers that contain two EGF molecules. The data obtained from the two approaches suggest a model in which one EGF monomer binds to one sEGFR monomer, and that receptor dimerization involves subsequent association of two monomeric (1:1) EGF-sEGFR complexes. Dimerization may result from bivalent binding of both EGF molecules in the dimer and/or receptor-receptor interactions. The requirement for two (possibly bivalent) EGF monomers distinguishes EGF-induced sEGFR dimerization from the hGH and interferon-gamma receptors, where multivalent binding of a single ligand species (either monomeric or dimeric) drives receptor oligomerization. The proposed model of EGF-induced sEGFR dimerization suggests possible mechanisms for both ligand-induced homo- and heterodimerization of the EGFR (or erbB) family of receptors.     

A three-dimensional model of the U1 small nuclear ribonucleoprotein particle

Most of the pre-mRNAs in the eukaryotic cell are comprised of protein-coding exons and non-protein-coding introns. The introns are removed and the exons are ligated together, or spliced, by a large, macromolecular complex known as the spliceosome. This RNA-protein assembly is made up of five uridine-rich small nuclear RNAs (U1-, U2-, U4-, U5- and U6-snRNA) as well over 300 proteins, which form small nuclear ribonucleoprotein particles (snRNPs). Initial recognition of the 5′ exon/intron splice site is mediated by the U1 snRNP, which is composed of the U1 snRNA as well as at least ten proteins. By combining structural informatics tools with the available biochemical and crystallographic data, we attempted to simulate a complete, three dimensional U1 snRNP from the silk moth, Bombyx mori. Comparison of our model with empirically derived crystal structures and electron micrographs pinpoints both the strengths and weaknesses in the in silico determination of macromolecular complexes. One of the most striking differences between our model and experimentally generated structures is in the positioning of the U1 snRNA stem-loops. This highlights the continuing difficulties in generating reliable, complex RNA structures; however, three-dimensional modeling of individual protein subunits by threading provided models of biological significance and the use of both automated and manual docking strategies generated a complex that closely reflects the assembly found in nature. Yet, without utilizing experimentally-derived contacts to select the most likely docking scenario, ab initio docking would fall short of providing a reliable model. Our work shows that the combination of experimental data with structural informatics tools can result in generation of near-native macromolecular complexes.     

Effects of 22K or 20K human growth hormone on lipolysis, leptin production in adipocytes in the presence and absence of human growth hormone binding protein.

OBJECTIVE AND METHOD: We studied the effects of human growth hormone (hGH) on leptin production and lipolysis stimulation in the presence or absence of human growth hormone binding protein (hGHBP) using 3T3- L1-hGHR adipocytes which efficiently express human growth hormone receptor. RESULTS AND CONCLUSION: It was clarified that (1) hGH decreases leptin secretion after hGH-induced lipolysis stimulation, and (2) the reduction of leptin production and lipolysis stimulation by 22K hGH was attenuated with hGHBP, whereas that by 20K hGH, which is a naturally occurring isoform of 22K hGH, was not affected with hGHBP.