New capabilities in determining the binding parameters for ligand-receptor interaction

A new method for determining the binding parameters of ligand-receptor interaction is suggested. The method is based on the application of the so-called coordinate of dilution, suggested by us earlier. We demonstrated that it is possible to determine the binding characteristics of ligand-receptor interaction using either the measurement of the concentration of the ligand-receptor complex at a state of equilibrium or the concentration of free receptors at different dilutions of the studying ligand-receptor mixture. The method also allows the determination of the concentration of the ligand in a pre-existing ligand-receptor mixture without preliminary separation of the interacting counterparts. For this reason the suggested method could be especially useful when the studying very labile receptors for which purification from the corresponding ligand is very difficult or impossible.     

Determining the parameters for receptor-ligand interaction by serial dilution method for the case when the ligand and receptor are in a pre-existing mixture

New methods of determining the binding parameters for ligand-receptor interaction are considered. The considered approaches are based on the earlier suggested method of serial dilution and application of so-called coordinates of dilution. It was shown that the suggested methods allow to evaluate affinity constant and ligand concentration even for the case, when the receptor and corresponding ligand of unknown concentration are in a mixture and their separation from each other is impossible. In this connection the suggested methods are especially useful for studying the ligand-receptor interaction if the receptor is very liable and its purification from the ligand would cause drastic changes of its binding properties.     

Consequences of ligand bivalency in interactions involving particulate receptors: equilibrium and kinetic studies with Sephadex-concanavalin A, butylagarose-phosphorylase b, and Fc receptor-IgG dimer interactions as model systems

Theoretical consideration is given to the interaction of a bivalent ligand with particulate receptor sites, not only from the viewpoint of quantitatively describing the binding behavior but also from that of the kinetics of ligand release upon infinite dilution of a receptor-ligand mixture. In the latter regard, a general expression is derived that describes the time dependence of the amount of ligand bound as a function of two rate constants for the stepwise dissociation of cross-linked ligand-receptor complex and a thermodynamic parameter expressing the initial ratio of singly linked to doubly linked ligand-receptor complexes. An experimental study of the interaction between Sephadex and concanavalin A is then used to illustrate application of this recommended theoretical approach for characterizing the binding behavior and dissociation kinetics of a bivalent ligand for a system in which all ligand-receptor interactions may be described by a single intrinsic association constant. Published results on the interaction of phosphorylase b with butylagarose are also shown to comply with this simplest model of the bivalent ligand hypothesis; but those for the interaction between immunoglobulin G (IgG) dimers and Fc receptors require modification of the model by incorporation of different intrinsic association constants for the successive binding of receptor sites to the bivalent ligand. These results emphasize the need to consider ligand bivalency as a potential phenomenon in studies of interactions between protein ligands and particulate receptors and illustrate procedures by which the effects of ligand bivalency may be identified and characterized.     

Physico-chemical study of the heterogeneous population of opioid receptors

A physico-chemical analysis of the heterogenous population of opioid receptors was carried out. A new difference approach to the analysis of heterogenous receptor systems was developed. This procedure permits to estimate even in high stringency conditions of parameters of three or more types of receptors from the binding isotherms or competitive replacement curves as well as to determine the total receptor concentration in the system. A DELTA computer program based on this difference approach has been developed. The experimental results obtained through the use of the difference technique led to a kinetic model of interaction between morphine and D-Ala2, D-Leu5-enkephalin with rat brain receptors. This model based on the interaction of each ligand with three types of specific binding sites can be used for the determination of kinetic and thermodynamic parameters.     

Stoichiometry of interaction between interferon gamma and its receptor.

The biological response of interferon gamma is mediated by binding to a specific cell-surface receptor. We investigated the stoichiometry of this binding using soluble receptors produced in prokaryotic and eukaryotic expression systems comprising the extracellular ligand-binding domain of the native protein. The ligand-receptor complexes were analyzed by cross-linking, chromatography, analytical ultracentrifugation and laser-light scattering. Cross-linking and chromatography showed that the stoichiometry of the interaction between ligand and receptor depends on the molar ratios of the two components mixed. All approaches confirmed that mixtures of ligand-receptor complexes are formed with one interferon-gamma dimer bound by one or two receptors. The soluble receptor produced in Escherichia coli mainly showed a ligand/receptor stoichiometry of 1:1, while the receptors produced in eukaryotic cells showed a stoichiometry of binding of 1:2. This apparent discrepancy is most likely due to the conformational heterogeneity of the Escherichia-coli-derived protein.     

Reply to Claudio Cobelli and Gianna Toffolo,identifiability from parameter bounds. Structural and numerical aspects. Math. Biosci. 71:237–243 (1984)

We present an analysis of receptor mediated endocytosis which includes the following elements: ligand binding to receptors, interaction of the ligand-receptor complex with coated pits, internalization of coated pit contents, recycling of receptors, and degradation of ligand. The model accounts quantitatively for epidermal growth factor binding and clustering in coated pits at 4°C, for its internalization and degradation at 37°C, and for EGF receptor down-regulation. Steady state analysis of the model indicates that the slope and intercept of a Scatchard plot are functions of the kinetic parameters of the endocytic loop and do not necessarily reflect the affinity and number of receptors in metabolically active cells. Moreover, the model predicts that for homogeneous receptors, a Scatchard plot can be either linear or nonlinear, depending on the concentration of proteins in coated pits which interact with ligand-receptor complexes. A slight generalization of the model in which phorbol ester-receptor complexes compete with EGF-receptor complexes for the same coated pit proteins provides a quantitative explanation for the loss of the high affinity portion of the EGF Scatchard plot subsequent to preincubation with phorbol esters. This explanation leads to the prediction of a local homology between a portion of the phorbol ester receptor sequence and a portion of the EGF receptor sequence.     

A Tunable Coarse-Grained Model for Ligand-Receptor Interaction

Cell-surface receptors are the most common target for therapeutic drugs. The design and optimization of next generation synthetic drugs require a detailed understanding of the interaction with their corresponding receptors. Mathematical approximations to study ligand-receptor systems based on reaction kinetics strongly simplify the spatial constraints of the interaction, while full atomistic ligand-receptor models do not allow for a statistical many-particle analysis, due to their high computational requirements. Here we present a generic coarse-grained model for ligand-receptor systems that accounts for the essential spatial characteristics of the interaction, while allowing statistical analysis. The model captures the main features of ligand-receptor kinetics, such as diffusion dependence of affinity and dissociation rates. Our model is used to characterize chimeric compounds, designed to take advantage of the receptor over-expression phenotype of certain diseases to selectively target unhealthy cells. Molecular dynamics simulations of chimeric ligands are used to study how selectivity can be optimized based on receptor abundance, ligand-receptor affinity and length of the linker between both ligand subunits. Overall, this coarse-grained model is a useful approximation in the study of systems with complex ligand-receptor interactions or spatial constraints.     

Structure analysis of bone morphogenetic protein-2 type I receptor complexes reveals a mechanism of receptor inactivation in juvenile polyposis syndrome

Bone morphogenetic proteins regulate many developmental processes during embryogenesis as well as tissue homeostasis in the adult. Signaling of bone morphogenetic proteins (BMPs) is accomplished by binding to two types of serine/threonine kinase transmembrane receptors termed type I and type II. Because a large number of ligands signal through a limited number of receptors, ligand-receptor interaction in the BMP superfamily is highly promiscuous, with a ligand binding to various receptors and a receptor binding many different BMP ligands. In this study we investigate the interaction of BMP-2 with its two high affinity type I receptors, BMP receptors IA (BMPR-IA) and BMPR-IB. Interestingly, 50% of the residues in the BMP-2 binding epitope of the BMPR-IA receptor are exchanged in BMPR-IB without a decrease in binding affinity or specificity for BMP-2. Our structural and functional analyses show that promiscuous binding of BMP-2 to both type I receptors is achieved by inherent backbone and side-chain flexibility as well as by variable hydration of the ligand-receptor interface enabling the BMP-2 surface to adapt to different receptor geometries. Despite the high degree of amino acid variability found in BMPR-IA and BMPR-IB binding equally to BMP-2, three single point missense mutations in the ectodomain of BMPR-IA cannot be tolerated. In juvenile polyposis syndrome these mutations have been shown to inactivate BMPR-IA. On the basis of our biochemical and biophysical analyses, we can show that the mutations, which are located outside the ligand binding epitope, alter the local or global fold of the receptor, thereby inactivating BMPR-IA and causing a loss of the BMP-2 tumor suppressor function in colon epithelial cells.     

Stacking interaction and its role in kynurenic acid binding to glutamate ionotropic receptors

Stacking interaction is known to play an important role in protein folding, enzyme-substrate and ligand-receptor complex formation. It has been shown to make a contribution into the aromatic antagonists binding with glutamate ionotropic receptors (iGluRs), in particular, the complex of NMDA receptor NR1 subunit with the kynurenic acid (KYNA) derivatives. The specificity of KYNA binding to the glutamate receptors subtypes might partially result from the differences in stacking interaction. We have calculated the optimal geometry and binding energy of KYNA dimers with the four types of aromatic amino acid residues in Rattus and Drosophila ionotropic iGluR subunits. All ab initio quantum chemical calculations were performed taking into account electron correlations at MP2 and MP4 perturbation theory levels. We have also investigated the potential energy surfaces (PES) of stacking and hydrogen bonds (HBs) within the receptor binding site and calculated the free energy of the ligand-receptor complex formation. The energy of stacking interaction depends both on the size of aromatic moieties and the electrostatic effects. The distribution of charges was shown to determine the geometry of polar aromatic ring dimers. Presumably, stacking interaction is important at the first stage of ligand binding when HBs are weak. The freedom of ligand movements and rotation within receptor site provides the precise tuning of the HBs pattern, while the incorrect stacking binding prohibits the ligand-receptor complex formation.     

Reevaluation of ANS binding to human and bovine serum albumins: key role of equilibrium microdialysis in ligand - receptor binding characterization

In this work we return to the problem of the determination of ligand-receptor binding stoichiometry and binding constants. In many cases the ligand is a fluorescent dye which has low fluorescence quantum yield in free state but forms highly fluorescent complex with target receptor. That is why many researchers use dye fluorescence for determination of its binding parameters with receptor, but they leave out of account that fluorescence intensity is proportional to the part of the light absorbed by the solution rather than to the concentration of bound dye. We showed how ligand-receptor binding parameters can be determined by spectrophotometry of the solutions prepared by equilibrium microdialysis. We determined the binding parameters of ANS - human serum albumin (HSA) and ANS - bovine serum albumin (BSA) interaction, absorption spectra, concentration and molar extinction coefficient, as well as fluorescence quantum yield of the bound dye. It was found that HSA and BSA have two binding modes with significantly different affinity to ANS. Correct determination of the binding parameters of ligand-receptor interaction is important for fundamental investigations and practical aspects of molecule medicine and pharmaceutics. The data obtained for albumins are important in connection with their role as drugs transporters.     

Steric effects on multivalent ligand-receptor binding: exclusion of ligand sites by bound cell surface receptors.

Steric effects can influence the binding of a cell surface receptor to a multivalent ligand. To account for steric effects arising from the size of a receptor and from the spacing of binding sites on a ligand, we extend a standard mathematical model for ligand-receptor interactions by introducing a steric hindrance factor. This factor gives the fraction of unbound ligand sites that are accessible to receptors, and thus available for binding, as a function of ligand site occupancy. We derive expressions for the steric hindrance factor for various cases in which the receptor covers a compact region on the ligand surface and the ligand expresses sites that are distributed regularly or randomly in one or two dimensions. These expressions are relevant for ligands such as linear polymers, proteins, and viruses. We also present numerical algorithms that can be used to calculate steric hindrance factors for other cases. These theoretical results allow us to quantify the effects of steric hindrance on ligand-receptor kinetics and equilibria.     

The N-terminal globular domain of Eph receptors is sufficient for ligand binding and receptor signaling.

The Eph family of receptor protein-tyrosine kinases (RTKs) have recently been implicated in patterning and wiring events in the developing nervous system. Eph receptors are unique among other RTKs in that they fall into two large subclasses that show distinct ligand specificities and for the fact that they themselves might function as ''ligands'', thereby activating bidirectional signaling. To gain insight into the mechanisms of ligand-receptor interaction, we have mapped the ligand binding domain in Eph receptors. By using a series of deletion and domain substitution mutants, we now report that an N-terminal globular domain of the Nuk/Cek5 receptor is the ligand binding domain of the transmembrane ligand Lerk2. Using focus formation assays, we show that the Cek5 globular domain is sufficient to confer Lerk2-dependent transforming activity on the Cek9 orphan receptor. Extending our binding studies to other members of both subclasses of receptors, it became apparent that the same domain is used for binding of both transmembrane and glycosylphosphatidyl-anchored ligands. Our studies have determined the first structural elements involved in ligand-receptor interaction and will allow more fine-tuned genetic experiments to elucidate the mechanism of action of these important guidance molecules.     

Effect of monensin on receptor recycling during continuous endocytosis of asialoorosomucoid

The binding of asialoglycoproteins to their liver cell receptor results in internalization of the ligand-receptor complex. These complexes rapidly appear in intracellular compartments termed endosomes whose acidification results in ligand-receptor dissociation. Ligand and receptor subsequently segregate: ligand is transported to lysosomes and is degraded while receptor recycles to the cell surface. The proton ionophore monensin prevents acidification of endosomes and reversibly inhibits this acid-dependent dissociation of ligand from receptor. The present study determined the effect of monensin treatment of short-term cultured rat hepatocytes on cell-surface-receptor content, determined both by their binding activity and immunologically, following continuous endocytosis of asialoorosomucoid. Inclusion of 5 microM monensin in the incubation medium reduced the number of immunologically detectable cell-surface receptors by 20% in the absence of ligand. During continuous endocytosis of asialoorosomucoid, inclusion of monensin resulted in a 30-40% reduction of cell-surface receptor detectable either by ligand binding or immunologically. These results suggest that the reduced liver-cell-surface content of receptor in monensin is due to intracellular trapping of ligand-receptor complexes. The reduction of surface receptor during monensin incubation in the absence of ligand suggests that "constitutive recycling" of plasma membrane components also requires intracellular acidification.     

Analysis ligand-receptor interactions from the molecular to the organism levels

A system of quantitative analysis is proposed for evaluation of ligand-receptor interaction on models of different levels of complexity. For two discrete receptor pools, binding of radio-labelled ligands to specific receptors and the magnitude of physiological response for an effector system with two discrete pools of receptors with different affinities are described with the aid of developed respective equations. The equations' parameters characterise properties of the effector system under study: the number of receptor pools differing in their affinity to a ligand; the number of active receptors of the maximum response magnitude; and the number of ligand's molecules bound to the receptor. The derived parameter's efficiency provides a general characteristic of affinity for the effector system under study. The described method of analysis of the ligand-receptor interactions is applicable to studies of any biological responses yielding quantitative results.     

Assembly of TbetaRI:TbetaRII:TGFbeta ternary complex in vitro with receptor extracellular domains is cooperative and isoform-dependent

Transforming growth factor-beta (TGFbeta) isoforms initiate signaling by assembling a heterotetrameric complex of paired type I (TbetaRI) and type II (TbetaRII) receptors on the cell surface. Because two of the ligand isoforms (TGFbetas 1, 3) must first bind TbetaRII to recruit TbetaRI into the complex, and a third (TGFbeta2) requires a co-receptor, assembly is known to be sequential, cooperative and isoform-dependent. However the source of the cooperativity leading to recruitment of TbetaRI and the universality of the assembly mechanism with respect to isoforms remain unclear. Here, we show that the extracellular domain of TbetaRI (TbetaRI-ED) binds in vitro with high affinity to complexes of the extracellular domain of TbetaRII (TbetaRII-ED) and TGFbetas 1 or 3, but not to either ligand or receptor alone. Thus, recruitment of TbetaRI requires combined interactions with TbetaRII-ED and ligand, but not membrane attachment of the receptors. Cell-based assays show that TbetaRI-ED, like TbetaRII-ED, acts as an antagonist of TGFbeta signaling, indicating that receptor-receptor interaction is sufficient to compete against endogenous, membrane-localized receptors. On the other hand, neither TbetaRII-ED, nor TbetaRII-ED and TbetaRI-ED combined, form a complex with TGFbeta2, showing that receptor-receptor interaction is insufficient to compensate for weak ligand-receptor interaction. However, TbetaRII-ED does bind with high affinity to TGFbeta2-TM, a TGFbeta2 variant substituted at three positions to mimic TGFbetas 1 and 3 at the TbetaRII binding interface. This proves both necessary and sufficient for recruitment of TbetaRI-ED, suggesting that the three different TGFbeta isoforms induce assembly of the heterotetrameric receptor complex in the same general manner.     

The capping of receptors for the epidermal growth factor and the participation of the cytoskeleton in this process in A-431 cells

Epidermal growth factor (EGF) receptor capping results from the interaction between the receptors and polyvalent ligands in A-431 cells examined in suspension at 22 degrees C. Colocalization of actin and spectrin with the ligand-receptor complexes during the redistribution was shown using double immunofluorescence. The obtained data show that the cortical microfilaments are involved in capping. EGF receptors become associated with the Triton-insoluble cytoskeleton as a consequence of ligand binding. EGF-receptor capping is not sensitive to the action of cytochalasin B. Capping in A-431 cells is discussed as a new model for studying the redistribution of the ligand-receptor complex.     

Multiple receptor states are required to describe both kinetic binding and activation of neutrophils via N-formyl peptide receptor ligands

It is well-established that the binding of N-formyl peptides to the N-formyl peptide receptor on neutrophils can be described by a kinetic scheme that involves two ligand-bound receptor states, both a low affinity ligand-receptor complex and a high affinity ligand-receptor complex, and that the rate constants describing ligand-receptor binding and receptor affinity state interconversion are ligand-specific. Here we examine whether differences due to these rate constants, i.e. differences in the numbers and lifetimes of particular receptor states, are correlated with neutrophil responses, namely actin polymerization and oxidant production. We find that an additional receptor state, one not discerned from kinetic binding assays, is required to account for these responses. This receptor state is interpreted as the number of low affinity bound receptors that are capable of activating G proteins; in other words, the accumulation of these active receptors correlates with the extent of both responses. Furthermore, this analysis allows for the quantification of a parameter that measures the relative strength of a ligand to bias the receptor into the active conformation. A model with this additional receptor state is sufficient to describe response data when two ligands (agonist/agonist or agonist/antagonist pairs) are added simultaneously, suggesting that cells respond to the accumulation of active receptors regardless of the identity of the ligand(s).