Effects of receptor clustering on ligand dissociation kinetics: theory and simulations

Sugar-Coated: Can Multivalent Glycoconjugates Improve upon Nature's Design?

Multivalent interactions between receptors and glycans play an important role in many different biological processes, including pathogen infection, self-recognition, and the immune response. The growth in the number of tools and techniques toward the assembly of multivalent glycoconjugates means it is possible to create synthetic systems that more and more closely resemble the diversity and complexity we observe in nature. In this Perspective we present the background to the recognition and binding enabled by multivalent interactions in nature, and discuss the strategies used to construct synthetic glycoconjugate equivalents. We highlight key discoveries and the current state of the art in their applications to glycan arrays, vaccines, and other therapeutic and diagnostic tools, with an outlook toward some areas we believe are of most interest for future work in this area.

Diversity in Notch ligand-receptor signaling interactions

The Notch signaling pathway uses families of ligands and receptors to transmit signals to nearby cells. These components are expressed in diverse combinations in different cell types, interact in a many-to-many fashion, both within the same cell (in cis) and between cells (in trans), and their interactions are modulated by Fringe glycosyltransferases. A fundamental question is how the strength of Notch signaling depends on which pathway components are expressed, at what levels, and in which cells. Here, we used a quantitative, bottom-up, cell-based approach to systematically characterize trans-activation, cis-inhibition, and cis-activation signaling efficiencies across a range of ligand and Fringe expression levels in two mammalian cell types. Each ligand (Dll1, Dll4, Jag1, and Jag2) and receptor variant (Notch1 and Notch2) analyzed here exhibited a unique profile of interactions, Fringe-dependence, and signaling outcomes. All four ligands were able to bind receptors in cis and in trans, and all ligands trans-activated both receptors, although Jag1-Notch1 signaling was substantially weaker than other ligand-receptor combinations. Cis-interactions were predominantly inhibitory, with the exception of the Dll1- and Dll4-Notch2 pairs, which exhibited cis-activation stronger than trans-activation. Lfng strengthened Delta-mediated trans-activation and weakened Jagged-mediated trans-activation for both receptors. Finally, cis-ligands showed diverse cis-inhibition strengths, which depended on the identity of the trans-ligand as well as the receptor. The map of receptor-ligand-Fringe interaction outcomes revealed here should help guide rational perturbation and control of the Notch pathway.

Mechanistic computational modeling of monospecific and bispecific antibodies targeting interleukin-6/8 receptors

The spread of cancer from organ to organ (metastasis) is responsible for the vast majority of cancer deaths; however, most current anti-cancer drugs are designed to arrest or reverse tumor growth without directly addressing disease spread. It was recently discovered that tumor cell-secreted interleukin-6 (IL-6) and interleukin-8 (IL-8) synergize to enhance cancer metastasis in a cell-density dependent manner, and blockade of the IL-6 and IL-8 receptors (IL-6R and IL-8R) with a novel bispecific antibody, BS1, significantly reduced metastatic burden in multiple preclinical mouse models of cancer. Bispecific antibodies (BsAbs), which combine two different antigen-binding sites into one molecule, are a promising modality for drug development due to their enhanced avidity and dual targeting effects. However, while BsAbs have tremendous therapeutic potential, elucidating the mechanisms underlying their binding and inhibition will be critical for maximizing the efficacy of new BsAb treatments. Here, we describe a quantitative, computational model of the BS1 BsAb, exhibiting how modeling multivalent binding provides key insights into antibody affinity and avidity effects and can guide therapeutic design. We present detailed simulations of the monovalent and bivalent binding interactions between different antibody constructs and the IL-6 and IL-8 receptors to establish how antibody properties and system conditions impact the formation of binary (antibody-receptor) and ternary (receptor-antibody-receptor) complexes. Model results demonstrate how the balance of these complex types drives receptor inhibition, providing important and generalizable predictions for effective therapeutic design.

Mechanistic computational modeling of monospecific and bispecific antibodies targeting interleukin-6/8 receptors

The spread of cancer from organ to organ (metastasis) is responsible for the vast majority of cancer deaths; however, most current anti-cancer drugs are designed to arrest or reverse tumor growth without directly addressing disease spread. It was recently discovered that tumor cell-secreted interleukin-6 (IL-6) and interleukin-8 (IL-8) synergize to enhance cancer metastasis in a cell-density dependent manner, and blockade of the IL-6 and IL-8 receptors (IL-6R and IL-8R) with a novel bispecific antibody, BS1, significantly reduced metastatic burden in multiple preclinical mouse models of cancer. Bispecific antibodies (BsAbs), which combine two different antigen-binding sites into one molecule, are a promising modality for drug development due to their enhanced avidity and dual targeting effects. However, while BsAbs have tremendous therapeutic potential, elucidating the mechanisms underlying their binding and inhibition will be critical for maximizing the efficacy of new BsAb treatments. Here, we describe a quantitative, computational model of the BS1 BsAb, exhibiting how modeling multivalent binding provides key insights into antibody affinity and avidity effects and can guide therapeutic design. We present detailed simulations of the monovalent and bivalent binding interactions between different antibody constructs and the IL-6 and IL-8 receptors to establish how antibody properties and system conditions impact the formation of binary (antibody-receptor) and ternary (receptor-antibody-receptor) complexes. Model results demonstrate how the balance of these complex types drives receptor inhibition, providing important and generalizable predictions for effective therapeutic design.

SARS-CoV-2 Infection-Of Music and Mechanics of Its Spikes! A Perspective

The COVID-19 pandemic has been inflicted upon humanity by the SARS-CoV-2 virus, the latest insidious incarnation of the coronaviruses group. While in its wake intense scientific research has produced breakthrough vaccines and cures, there still exists an immediate need to further understand the origin, mechanobiology and biochemistry, and destiny of this virus so that future pandemics arising from similar coronaviruses may be contained more effectively. In this Perspective, we discuss the various evidential findings of virus propagation and connect them to respective underpinning cellular biomechanical states leading to corresponding manifestations of the viral activity. We further propose avenues to tackle the virus, including from a "musical" vantage point, and contain its relentless strides that are currently afflicting the global populace.

Spatially resolved in silico modeling of NKG2D signaling kinetics suggests a key role of NKG2D and Vav1 Co-clustering in generating natural killer cell activation

Natural Killer (NK) cells provide key resistance against viral infections and tumors. A diverse set of activating and inhibitory NK cell receptors (NKRs) interact with cognate ligands presented by target host cells, where integration of dueling signals initiated by the ligand-NKR interactions determines NK cell activation or tolerance. Imaging experiments over decades have shown micron and sub-micron scale spatial clustering of activating and inhibitory NKRs. The mechanistic roles of these clusters in affecting downstream signaling and activation are often unclear. To this end, we developed a predictive in silico framework by combining spatially resolved mechanistic agent based modeling, published TIRF imaging data, and parameter estimation to determine mechanisms by which formation and spatial movements of activating NKG2D microclusters affect early time NKG2D signaling kinetics in a human cell line NKL. We show co-clustering of NKG2D and the guanosine nucleotide exchange factor Vav1 in NKG2D microclusters plays a dominant role over ligand (ULBP3) rebinding in increasing production of phospho-Vav1(pVav1), an activation marker of early NKG2D signaling. The in silico model successfully predicts several scenarios of inhibition of NKG2D signaling and time course of NKG2D spatial clustering over a short (~3 min) interval. Modeling shows the presence of a spatial positive feedback relating formation and centripetal movements of NKG2D microclusters, and pVav1 production offers flexibility towards suppression of activating signals by inhibitory KIR ligands organized in inhomogeneous spatial patterns (e.g., a ring). Our in silico framework marks a major improvement in developing spatiotemporal signaling models with quantitatively estimated model parameters using imaging data.

Natural Killer cells are lymphocytes of our innate immunity and provide important resistance against viral infections and tumors. NK cells scan the local environment with diverse activating and inhibitory NK cell receptors (NKRs) and remain tolerized or lyse target cells expressing cognate ligands to NKRs. NKRs have been found to form micron sized clusters (or microclusters) as they interact with cognate ligands, and mechanisms regarding how the formation and movements of these microclusters influence NK cell signaling and activation, specifically related to activating NKRs, are often unclear. To this end, we develop a predictive spatially resolved early-time NK cell signaling model to study the interplay between membrane-proximal biochemical signaling events and the kinetics of microclusters of activating NKG2D and inhibitory KIR2DL2 receptors. We used published TIRF imaging data to validate our in silico models and estimate model parameters. Predictions from multiple in silico models are tested against a variety of data obtained from published imaging experiments and immunoassays. Our analysis suggests co-clustering of NKG2D and the guanosine nucleotide exchange factor Vav1 in the microclusters plays a major role in enhancing downstream activating signals. The developed framework can be extended to describe spatiotemporal signaling for other activating NKRs including CD16.

Significance of Receptor Mobility in Multivalent Binding on Lipid Membranes

Numerous key biological processes rely on the concept of multivalency, where ligands achieve stable binding only upon engaging multiple receptors. These processes, like viral entry or immune synapse formation, occur on the diffusive cellular membrane. One crucial, yet underexplored aspect of multivalent binding is the mobility of coupled receptors. Here, we discuss the consequences of mobility in multivalent processes from four perspectives: (I) The facilitation of receptor recruitment by the multivalent ligand due to their diffusivity prior to binding. (II) The effects of receptor preassembly, which allows their local accumulation. (III) The consequences of changes in mobility upon the formation of receptor/ligand complex. (IV) The changes in the diffusivity of lipid environment surrounding engaged receptors. We demonstrate how understanding mobility is essential for fully unravelling the principles of multivalent membrane processes, leading to further development in studies on receptor interactions, and guide the design of new generations of multivalent ligands.

The Future of the COVID-19 Pandemic: How Good (or Bad) Can the SARS-CoV2 Spike Protein Get?

Severe acute respiratory syndrome virus 2 (SARS-CoV2) has infected an estimated 400 million people world-wide, causing approximately 6 million deaths from severe coronavirus disease 2019 (COVID-19). The SARS-CoV2 Spike protein plays a critical role in viral attachment and entry into host cells. The recent emergence of highly transmissible variants of SARS-CoV2 has been linked to mutations in Spike. This review provides an overview of the structure and function of Spike and describes the factors that impact Spike’s ability to mediate viral infection as well as the potential limits to how good (or bad) Spike protein can become. Proposed here is a framework that considers the processes of Spike-mediated SARS-CoV2 attachment, dissociation, and cell entry where the role of Spike, from the standpoint of the virus, is to maximize cell entry with each viral-cell collision. Key parameters are identified that will be needed to develop models to identify mechanisms that new Spike variants might exploit to enhance viral transmission. In particular, the importance of considering secondary co-receptors for Spike, such as heparan sulfate proteoglycans is discussed. Accurate models of Spike-cell interactions could contribute to the development of new therapies in advance of the emergence of new highly transmissible SARS-CoV2 variants.

Mechanical disruption of E-cadherin complexes with epidermal growth factor receptor actuates growth factor-dependent signaling

Increased intercellular tension is associated with enhanced cell proliferation and tissue growth. Here, we present evidence for a force-transduction mechanism that links mechanical perturbations of epithelial (E)-cadherin (CDH1) receptors to the force-dependent activation of epidermal growth factor receptor (EGFR, ERBB1)-a key regulator of cell proliferation. Here, coimmunoprecipitation studies first show that E-cadherin and EGFR form complexes at the plasma membrane that are disrupted by either epidermal growth factor (EGF) or increased tension on homophilic E-cadherin bonds. Although force on E-cadherin bonds disrupts the complex in the absence of EGF, soluble EGF is required to mechanically activate EGFR at cadherin adhesions. Fully quantified spectral imaging fluorescence resonance energy transfer further revealed that E-cadherin and EGFR directly associate to form a heterotrimeric complex of two cadherins and one EGFR protein. Together, these results support a model in which the tugging forces on homophilic E-cadherin bonds trigger force-activated signaling by releasing EGFR monomers to dimerize, bind EGF ligand, and signal. These findings reveal the initial steps in E-cadherin-mediated force transduction that directly link intercellular force fluctuations to the activation of growth regulatory signaling cascades.

The Role of Ligand Rebinding and Facilitated Dissociation on the Characterization of Dissociation Rates by Surface Plasmon Resonance (SPR) and Benchmarking Performance Metrics

Surface plasmon resonance (SPR) is a real-time kinetic measurement principle that can probe the kinetic interactions between ligands and their binding sites, and lies at the backbone of pharmaceutical, biosensing, and biomolecular research. The extraction of dissociation rates from SPR-response signals often relies on several commonly adopted assumptions, one of which is the exponential decay of the dissociation part of the response signal. However, certain conditions, such as high density of binding sites or high concentration fluctuations near the surface as compared to the bulk, can lead to non-exponential decays via ligand rebinding or facilitated dissociation. Consequently, fitting the data with an exponential function can underestimate or overestimate the measured dissociation rates. Here, we describe a set of alternative fit functions that can take such effects into consideration along with plasmonic sensor design principles with key performance metrics, thereby suggesting methods for error-free high-precision extraction of the dissociation rates.

Distinct binding kinetics of E-, P- and L-selectins to CD44

Molecular-level selectin-cluster of differentiation 44 (CD44) interactions are far from clear because of the complexity and diversity of CD44 glycosylation and isoforms expressed on various types of cells. By combining experimental measurements and simulation predictions, the binding kinetics of three selectin members to the recombinant CD44 were quantified and the corresponding microstructural mechanisms were explored, respectively. Experimental results showed that the E-selectin-CD44 interactions mainly mediated the firm adhesion of microbeads under shear flow with the strongest rupture force. P- and L-selectins had similar interaction strength but different association and dissociation rates by mediating stable rolling and transient adhesions of microbeads, respectively. Molecular docking and molecular dynamics (MD) simulations predicted that the binding epitopes of CD44 to selectins are all located at the side face of each selectin, although the interfaces denoted as the hinge region are between lectin and epidermal growth factor domains of E-selectin, Lectin domain side of P-selectin and epidermal growth factor domain side of L-selectin, respectively. The lowest binding free energy, the largest rupture force and the longest lifetime for E-selectin, as well as the comparable values for P- and L-selectins, demonstrated in both equilibration and steered MD simulations, supported the above experimental results. These results offer basic data for understanding the functional differences of selectin-CD44 interactions.

Single-molecule imaging and microfluidic platform reveal molecular mechanisms of leukemic cell rolling

Hematopoietic stem/progenitor cell (HSPC) and leukemic cell homing is an important biological phenomenon that occurs through key interactions between adhesion molecules. Tethering and rolling of the cells on endothelium, the crucial initial step of the adhesion cascade, is mediated by interactions between selectins expressed on endothelium to their ligands expressed on HSPCs/leukemic cells in flow. Although multiple factors that affect the rolling behavior of the cells have been identified, molecular mechanisms that enable the essential slow and stable cell rolling remain elusive. Here, using a microfluidics-based single-molecule live cell fluorescence imaging, we reveal that unique spatiotemporal dynamics of selectin ligands on the membrane tethers and slings, which are distinct from that on the cell body, play an essential role in the rolling of the cell. Our results suggest that the spatial confinement of the selectin ligands to the tethers and slings together with the rapid scanning of a large area by the selectin ligands, increases the efficiency of selectin-ligand interactions during cell rolling, resulting in slow and stable rolling of the cell on the selectins. Our findings provide novel insights and contribute significantly to the molecular-level understanding of the initial and essential step of the homing process.

Mechanobiology predicts raft formations triggered by ligand-receptor activity across the cell membrane

Clustering of ligand-binding receptors of different types on thickened isles of the cell membrane, namely lipid rafts, is an experimentally observed phenomenon. Although its influence on cell's response is deeply investigated, the role of the coupling between mechanical processes and multiphysics involving the active receptors and the surrounding lipid membrane during ligand-binding has not yet been understood. Specifically, the focus of this work is on G-protein-coupled receptors (GPCRs), the widest group of transmembrane proteins in animals, which regulate specific cell processes through chemical signalling pathways involving a synergistic balance between the cyclic Adenosine Monophosphate (cAMP) produced by active GPCRs in the intracellular environment and its efflux, mediated by the Multidrug Resistance Proteins (MRPs) transporters. This paper develops a multiphysics approach based on the interplay among energetics, multiscale geometrical changes and mass balance of species, i.e. active GPCRs and MRPs, including diffusion and kinetics of binding and unbinding. Because the obtained energy depends upon both the kinematics and the changes of species densities, balance of mass and of linear momentum are coupled and govern the space-time evolution of the cell membrane. The mechanobiology involving remodelling and change of lipid ordering of the cell membrane allows to predict dynamics of transporters and active receptors -in full agreement with experimentally observed cAMP levels- and how the latter trigger rafts formation and cluster on such sites. Within the current scientific debate on Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) and on the basis of the ascertained fact that lipid rafts often serve as an entry port for viruses, it is felt that approaches accounting for strong coupling among mechanobiological aspects could even turn helpful in better understanding membrane-mediated phenomena such as COVID-19 virus-cell interaction.

Development of a Simple Kinetic Mathematical Model of Aggregation of Particles or Clustering of Receptors

The process of clustering of plasma membrane receptors in response to their agonist is the first step in signal transduction. The rate of the clustering process and the size of the clusters determine further cell responses. Here we aim to demonstrate that a simple 2-differential equation mathematical model is capable of quantitative description of the kinetics of 2D or 3D cluster formation in various processes. Three mathematical models based on mass action kinetics were considered and compared with each other by their ability to describe experimental data on GPVI or CR3 receptor clustering (2D) and albumin or platelet aggregation (3D) in response to activation. The models were able to successfully describe experimental data without losing accuracy after switching between complex and simple models. However, additional restrictions on parameter values are required to match a single set of parameters for the given experimental data. The extended clustering model captured several properties of the kinetics of cluster formation, such as the existence of only three typical steady states for this system: unclustered receptors, receptor dimers, and clusters. Therefore, a simple kinetic mass-action-law-based model could be utilized to adequately describe clustering in response to activation both in 2D and in 3D.

Evidence for an effect of receptor density on ligand occupancy and agonist EC50

Drug-receptor interaction theory predicts that proportional receptor occupancy is a function of ligand concentration as defined by a ligand-receptor affinity constant, and is independent of receptor density. However, we previously observed that the EC₅₀ of 5-HT reduced as the density of 5-HT₃ receptors increased, suggesting an effect of receptor density on occupancy. The current study was designed to maximise variability in experimentally observed currents and confirm this apparent contradiction prospectively. Xenopus oocytes were injected with RNA encoding 5-HT₃A receptors under conditions designed to achieve varying receptor expression levels and 5-HT-evoked currents measured using two electrode voltage clamp. Results from 99 oocytes showed that as the maximal peak current increased from 0.05 µA to 12.1 µA there was a 3.7-fold reduction in EC₅₀. Since occupancy and conductance are directly related in this system, this indicates that for a given concentration of 5-HT, proportional occupancy increases with increased receptor density. We conclude that normalising data masks this correlation, and can result in reduced accuracy of pharmacological measurements. We propose a mechanistic explanation for our observations.

The alpha-1 subunit of the Na+,K+-ATPase (ATP1A1) is required for macropinocytic entry of respiratory syncytial virus (RSV) in human respiratory epithelial cells

Human respiratory syncytial virus (RSV) is the leading viral cause of acute pediatric lower respiratory tract infections worldwide, with no available vaccine or effective antiviral drug. To gain insight into virus-host interactions, we performed a genome-wide siRNA screen. The expression of over 20,000 cellular genes was individually knocked down in human airway epithelial A549 cells, followed by infection with RSV expressing green fluorescent protein (GFP). Knockdown of expression of the cellular ATP1A1 protein, which is the major subunit of the Na⁺,K⁺-ATPase of the plasma membrane, had one of the strongest inhibitory effects on GFP expression and viral titer. Inhibition was not observed for vesicular stomatitis virus, indicating that it was RSV-specific rather than a general effect. ATP1A1 formed clusters in the plasma membrane very early following RSV infection, which was independent of replication but dependent on the attachment glycoprotein G. RSV also triggered activation of ATP1A1, resulting in signaling by c-Src-kinase activity that transactivated epidermal growth factor receptor (EGFR) by Tyr845 phosphorylation. ATP1A1 signaling and activation of both c-Src and EGFR were found to be required for efficient RSV uptake. Signaling events downstream of EGFR culminated in the formation of macropinosomes. There was extensive uptake of RSV virions into macropinosomes at the beginning of infection, suggesting that this is a major route of RSV uptake, with fusion presumably occurring in the macropinosomes rather than at the plasma membrane. Important findings were validated in primary human small airway epithelial cells (HSAEC). In A549 cells and HSAEC, RSV uptake could be inhibited by the cardiotonic steroid ouabain and the digitoxigenin derivative PST2238 (rostafuroxin) that bind specifically to the ATP1A1 extracellular domain and block RSV-triggered EGFR Tyr845 phosphorylation. In conclusion, we identified ATP1A1 as a host protein essential for macropinocytic entry of RSV into respiratory epithelial cells, and identified PST2238 as a potential anti-RSV drug.

RSV continues to be the most important viral cause of severe bronchiolitis and pneumonia in infants and young children, and also has a substantial impact in the elderly. It is estimated to claim the lives of ~118,000 children under five years of age annually. No vaccine or antiviral drug suitable for general use is available. The involvement of host factors in RSV infection and replication is not well understood, but this knowledge might lead to intervention strategies to prevent infection. Using a genome-wide siRNA screen to knock down the expression of over 20,000 individual cellular genes, we identified ATP1A1, the major subunit of the Na⁺,K⁺-ATPase, as an important host protein for RSV entry. We showed that ATP1A1 activation by RSV resulted in transactivation of EGFR by Src-kinase activity, resulting in the uptake of RSV particles into the host cell through macropinocytosis. We also showed that the cardiotonic steroid ouabain and the synthetic digitoxigenin derivative PST2238, which bind specifically to the extracellular domain of ATP1A1, significantly reduced RSV entry. Taken together, we describe a novel ATP1A1-enabled mechanism used by RSV to enter the host cell, and describe candidate antiviral drugs that block this entry.

Receptor-Ligand Rebinding Kinetics in Confinement

Rebinding kinetics of molecular ligands plays a key role in the operation of biomachinery, from regulatory networks to protein transcription, and is also a key factor in design of drugs and high-precision biosensors. In this study, we investigate initial release and rebinding of ligands to their binding sites grafted on a planar surface, a situation commonly observed in single-molecule experiments and that occurs in vivo, e.g., during exocytosis. Via scaling arguments and molecular dynamic simulations, we analyze the dependence of nonequilibrium rebinding kinetics on two intrinsic length scales: the average separation distance between the binding sites and the total diffusible volume (i.e., height of the experimental reservoir in which diffusion takes place or average distance between receptor-bearing surfaces). We obtain time-dependent scaling laws for on rates and for the cumulative number of rebinding events. For diffusion-limited binding, the (rebinding) on rate decreases with time via multiple power-law regimes before the terminal steady-state (constant on-rate) regime. At intermediate times, when particle density has not yet become uniform throughout the diffusible volume, the cumulative number of rebindings exhibits a novel, to our knowledge, plateau behavior because of the three-dimensional escape process of ligands from binding sites. The duration of the plateau regime depends on the average separation distance between binding sites. After the three-dimensional diffusive escape process, a one-dimensional diffusive regime describes on rates. In the reaction-limited scenario, ligands with higher affinity to their binding sites (e.g., longer residence times) delay entry to the power-law regimes. Our results will be useful for extracting hidden timescales in experiments such as kinetic rate measurements for ligand-receptor interactions in microchannels, as well as for cell signaling via diffusing molecules.

Advances in Receptor-Mediated, Tumor-Targeted Drug Delivery

Receptor-mediated drug delivery presents an opportunity to enhance therapeutic efficiency by accumulating drug within the tissue of interest and reducing undesired, off-target effects. In cancer, receptor overexpression is a platform for binding and inhibiting pathways that shape biodistribution, toxicity, cell binding and uptake, and therapeutic function. This review will identify tumor-targeted drug delivery vehicles and receptors that show promise for clinical translation based on quantitative in vitro and in vivo data. The authors describe the rationale to engineer a targeted drug delivery vehicle based on the ligand, chemical conjugation method, and type of drug delivery vehicle. Recent advances in multivalent targeting and ligand organization on tumor accumulation are discussed. Revolutionizing receptor-mediated drug delivery may be leveraged in the therapeutic delivery of chemotherapy, gene editing tools, and epigenetic drugs.

Link between a high k on for drug binding and a fast clinical action: to be or not to be?

Review articles on binding kinetics essentially focus on drugs that dissociate slowly from their target since this is required for the successful treatment of many pathophysiological conditions. Recently, the therapeutic benefit of a high k on (i.e. the second order association rate constant) has also been linked to fast association and to a fast clinical action. Other studies, however, called this assertion into question since additional factors, like the dosing paradigm and the binding mechanism, are important as well. The still ongoing reticence about integrating binding kinetics in lead optimization programs motivated us to critically review the link between the drug's kinetic rate constants and their in vitro and in vivo target occupancy profile, with special focus on k on. The presented simulations tally with a positive link between a drug's effective/observed association rate (which is quite easy to determine in vitro) and the swiftness of its clinical action. On the other hand, the simulations show that the k on-concept should not be confounded with the effective association process since increasing this parameter only enhances the drug's in vitro and in vivo association under certain conditions: the binding mechanism should be suitable, rebinding (and thus the factors within the target's micro-environment that favour this mechanism) should not be too prominent and the dosage should not be kept in par with the drug's affinity. Otherwise, increasing k on could be ineffective or even be counter-productive.

Ligand-receptor binding kinetics in surface plasmon resonance cells: a Monte Carlo analysis

Surface plasmon resonance (SPR) chips are widely used to measure association and dissociation rates for the binding kinetics between two species of chemicals, e.g., cell receptors and ligands. It is commonly assumed that ligands are spatially well mixed in the SPR region, and hence a mean-field rate equation description is appropriate. This approximation however ignores the spatial fluctuations as well as temporal correlations induced by multiple local rebinding events, which become prominent for slow diffusion rates and high binding affinities. We report detailed Monte Carlo simulations of ligand binding kinetics in an SPR cell subject to laminar flow. We extract the binding and dissociation rates by means of the techniques frequently employed in experimental analysis that are motivated by the mean-field approximation. We find major discrepancies in a wide parameter regime between the thus extracted rates and the known input simulation values. These results underscore the crucial quantitative importance of spatio-temporal correlations in binary reaction kinetics in SPR cell geometries, and demonstrate the failure of a mean-field analysis of SPR cells in the regime of high Damköhler number [Formula: see text], where the spatio-temporal correlations due to diffusive transport and ligand-receptor rebinding events dominate the dynamics of SPR systems.

Effects of target binding kinetics on in vivo drug efficacy: koff , kon and rebinding

BACKGROUND AND PURPOSE

Optimal drug therapy often requires continuing high levels of target occupancy. Besides the traditional pharmacokinetic contribution, target binding kinetics is increasingly considered to play an important role as well. While most attention has been focused on the dissociation rate of the complex, recent reports expressed doubt about the unreserved translatability of this pharmacodynamic property into clinical efficacy. 'Micro'-pharmacokinetic mechanisms like drug rebinding and partitioning into the cell membrane may constitute a potential fix.

EXPERIMENTAL APPROACH

Simulations were based on solving differential equations.

KEY RESULTS

Based on a selected range of association and dissociation rate constants, kon and koff , and rebinding potencies of the drugs as variables, their effects on the temporal in vivo occupancy profile of their targets, after one or multiple repetitive dosings, have here been simulated.

CONCLUSIONS AND IMPLICATIONS

Most strikingly, the simulations show that, when rebinding is also taken into account, increasing kon may produce closely the same outcome as decreasing koff when dosing is performed in accordance with the therapeutically most relevant constant [Lmax ]/KD ratio paradigm. Also, under certain conditions, rebinding may produce closely the same outcome as invoking slow diffusion of the drug between the plasma compartment and a target-containing 'effect' compartment. Although the present simulations should only be regarded as a 'proof of principle', these findings may help pharmacologists and medicinal chemists to devise ex vivo and in vitro binding kinetic assays that are more relevant and translatable to in vivo settings.

Affinity Enhancement by Ligand Clustering Effect Inspired by Peptide Dendrimers-Shank PDZ Proteins Interactions

High-affinity binders are desirable tools to probe the function that specific protein−protein interactions play in cell. In the process of seeking a general strategy to design high-affinity binders, we found a clue from the βPIX (p21-activated kinase interacting exchange factor)−Shank PDZ interaction in synaptic assembly: three PDZ-binding sites are clustered by a parallel coiled-coil trimer but bind to Shank PDZ protein with 1:1 stoichiometry (1 trimer/1 PDZ). Inspired by this architecture, we proposed that peptide dendrimer, mimicking the ligand clustering in βPIX, will also show enhanced binding affinity, yet with 1:1 stoichiometry. This postulation has been proven here, as we synthesized a set of monomeric, dimeric and trimeric peptides and measured their binding affinity and stoichiometry with Shank PDZ domains by isothermal titration calorimetry, native mass spectrometry and surface plasmon resonance. This affinity enhancement, best explained by proximity effect, will be useful to guide the design of high-affinity blockers for protein−protein interactions.

An induced rebinding model of antigen discrimination

•

We propose that pMHC binding to TCR can increase (induce) pMHC rebinding.

•

Published studies are consistent with an induced rebinding model.

•

Induced rebinding improves the ability of T cells to discriminate antigens.

•

Induced rebinding relates 3D to 2D TCR–pMHC binding parameters.

T cells have to detect rare high-affinity ‘foreign’ peptide MHC (pMHC) ligands among abundant low-affinity ‘self’-peptide MHC ligands. It remains unclear how this remarkable discrimination is achieved. Kinetic proofreading mechanisms can provide the required specificity but only at the expense of much reduced sensitivity. A number of recent observations suggest that pMHC engagement of T cell receptors (TCRs) induces changes such as clustering and/or conformational alterations that enhance subsequent rebinding. We show that inclusion of induced rebinding to the same pMHC in kinetic proofreading models enhances the sensitivity of TCR recognition while retaining specificity. Moreover, induced rebinding is able to reproduce the striking, and hitherto unexplained, 2D membrane-binding properties recently reported for the TCR.

Heparan sulfate-protein binding specificity

Heparan sulfate (HS) represents a large class of linear polysaccharides that are required for the function of all mammalian physiological systems. HS is characterized by a repeating disaccharide backbone that is subject to a wide range of modifications, making this class of macromolecules arguably the most information dense in all of biology. The majority of HS functions are associated with the ability to bind and regulate a wide range of proteins. Indeed, recent years have seen an explosion in the discovery of new activities for HS where it is now recognized that this class of glycans functions as co-receptors for growth factors and cytokines, modulates cellular uptake of lipoproteins, regulates protease activity, is critical to amyloid plaque formation, is used by opportunistic pathogens to enter cells, and may even participate in epigenetic regulation. This review will discuss the current state of understanding regarding the specificity of HS-protein binding and will describe the concept that protein binding to HS depends on the overall organization of domains within HS rather than fine structure.

Receptor clustering affects signal transduction at the membrane level in the reaction-limited regime

Many types of membrane receptors are found to be organized as clusters on the cell surface. We investigate the potential effect of such receptor clustering on the intracellular signal transduction stage. We consider a canonical pathway with a membrane receptor (R) activating a membrane-bound intracellular relay protein (G). We use Monte Carlo simulations to recreate biochemical reactions using different receptor spatial distributions and explore the dynamics of the signal transduction. Results show that activation of G by R is severely impaired by R clustering, leading to an apparent blunted biological effect compared to control. Paradoxically, this clustering decreases the half maximal effective dose (ED50) of the transduction stage, increasing the apparent affinity. We study an example of inter-receptor interaction in order to account for possible compensatory effects of clustering and observe the parameter range in which such interactions slightly counterbalance the loss of activation of G. The membrane receptors' spatial distribution affects the internal stages of signal amplification, suggesting a functional role for membrane domains and receptor clustering independently of proximity-induced receptor-receptor interactions.

Multivalent proteoglycan modulation of FGF mitogenic responses in perivascular cells

Sprouting of angiogenic perivascular cells is thought to be highly dependent upon autocrine and paracrine growth factor stimulation. Accordingly, we report that corneal angiogenesis induced by ectopic FGF implantation is strongly impaired in NG2/CSPG4 proteoglycan (PG) null mice known to harbour a putative deficit in pericyte proliferation/mobilization. Conversely, no significant differences were seen between wild type and knockout corneas when VEGF was used as an angiocrine factor. Perturbed responsiveness of NG2-deficient pericytes to paracrine and autocrine stimulation by several FGFs could be confirmed in cells isolated from NG2 null mice, while proliferation induced by other growth factors was equivalent in wild type and knockout cells. Identical results were obtained after siRNA-mediated knock-down of NG2 in human smooth muscle-like cell lines, as also demonstrated by the decreased levels of FGF receptor phosphorylation detected in these NG2 deprived cells. Binding assays with recombinant proteins and molecular interactions examined on live cells asserted that FGF-2 bound to NG2 in a glycosaminoglycan-independent, core protein-mediated manner and that the PG was alone capable of retaining FGF-2 on the cell membrane for subsequent receptor presentation. The use of dominant-negative mutant cells, engineered by combined transduction of NG2 deletion constructs and siRNA knock-down of the endogenous PG, allowed us to establish that the FGF co-receptor activity of NG2 is entirely mediated by its extracellular portion. In fact, forced overexpression of the NG2 ectodomain in human smooth muscle-like cells increased their FGF-2-induced mitosis and compensated for low levels of FGF receptor surface expression, in a manner equivalent to that produced by overexpression of the full-length NG2. Upon FGF binding, the cytoplasmic domain of NG2 is phosphorylated, but there is no evidence that this event elicits signal transductions that could bypass the FGFR-mediated ones. Pull-down experiments, protein-protein binding assays and flow cytometry FRET coherently revealed an elective ligand-independent association of NG2 with FGFR1 and FGFR3. The NG2 cooperation with these receptors was also corroborated functionally by the outcome of FGF-2 treatments of cells engineered to express diverse NG2/FGFR combinations. Comprehensively, the findings suggest that perivascular NG2 may serve as a dual modulator of the availability/accessibility of FGF at the cell membrane, as well as the resulting FGFR transducing activity.

Regulation from within: the cytoskeleton in transmembrane signaling

There is mounting evidence that the plasma membrane is highly dynamic and organized in a complex manner. The cortical cytoskeleton is proving to be a particularly important regulator of plasmalemmal organization, modulating the mobility of proteins and lipids in the membrane, facilitating their segregation, and influencing their clustering. This organization plays a critical role in receptor-mediated signaling, especially in the case of immunoreceptors, which require lateral clustering for their activation. Based on recent developments, we discuss the structures and mechanisms whereby the cortical cytoskeleton regulates membrane dynamics and organization, and how the nonuniform distribution of immunoreceptors and their self-association may affect activation and signaling.

Spatial simulations in systems biology: from molecules to cells

Cells are highly organized objects containing millions of molecules. Each biomolecule has a specific shape in order to interact with others in the complex machinery. Spatial dynamics emerge in this system on length and time scales which can not yet be modeled with full atomic detail. This review gives an overview of methods which can be used to simulate the complete cell at least with molecular detail, especially Brownian dynamics simulations. Such simulations require correct implementation of the diffusion-controlled reaction scheme occurring on this level. Implementations and applications of spatial simulations are presented, and finally it is discussed how the atomic level can be included for instance in multi-scale simulation methods.

Perspective: Dynamics of receptor tyrosine kinase signaling complexes

Textbook descriptions of signal transduction complexes provide a static snapshot view of highly dynamic events. Despite enormous strides in identifying the key components of signaling complexes and the underlying mechanisms of signal transduction, our understanding of the dynamic behavior of these complexes has lagged behind. Using the example of receptor tyrosine kinases, this perspective takes a fresh look at the dynamics of the system and their potential impact on signal processing.

Biomechanics of leukocyte rolling

Leukocyte rolling on endothelial cells and other P-selectin substrates is mediated by P-selectin binding to P-selectin glycoprotein ligand-1 expressed on the tips of leukocyte microvilli. Leukocyte rolling is a result of rapid, yet balanced formation and dissociation of selectin-ligand bonds in the presence of hydrodynamic shear forces. The hydrodynamic forces acting on the bonds may either increase (catch bonds) or decrease (slip bonds) their lifetimes. The force-dependent 'catch-slip' bond kinetics are explained using the 'two pathway model' for bond dissociation. Both the 'sliding-rebinding' and the 'allosteric' mechanisms attribute 'catch-slip' bond behavior to the force-induced conformational changes in the lectin-EGF domain hinge of selectins. Below a threshold shear stress, selectins cannot mediate rolling. This 'shear-threshold' phenomenon is a consequence of shear-enhanced tethering and catch bond-enhanced rolling. Quantitative dynamic footprinting microscopy has revealed that leukocytes rolling at venular shear stresses (>0.6 Pa) undergo cellular deformation (large footprint) and form long tethers. The hydrodynamic shear force and torque acting on the rolling cell are thought to be synergistically balanced by the forces acting on tethers and stressed microvilli, however, their relative contribution remains to be determined. Thus, improvement beyond the current understanding requires in silico models that can predict both cellular and microvillus deformation and experiments that allow measurement of forces acting on individual microvilli and tethers.

Impact of receptor clustering on ligand binding

Background

Cellular response to changes in the concentration of different chemical species in the extracellular medium is induced by ligand binding to dedicated transmembrane receptors. Receptor density, distribution, and clustering may be key spatial features that influence effective and proper physical and biochemical cellular responses to many regulatory signals. Classical equations describing this kind of binding kinetics assume the distributions of interacting species to be homogeneous, neglecting by doing so the impact of clustering. As there is experimental evidence that receptors tend to group in clusters inside membrane domains, we investigated the effects of receptor clustering on cellular receptor ligand binding.

Results

We implemented a model of receptor binding using a Monte-Carlo algorithm to simulate ligand diffusion and binding. In some simple cases, analytic solutions for binding equilibrium of ligand on clusters of receptors are provided, and supported by simulation results. Our simulations show that the so-called "apparent" affinity of the ligand for the receptor decreases with clustering although the microscopic affinity remains constant.

Conclusions

Changing membrane receptors clustering could be a simple mechanism that allows cells to change and adapt its affinity/sensitivity toward a given stimulus.

Monte Carlo simulations of plasma membrane corral-induced EGFR clustering

Experimental evidence suggests that the cell membrane is a highly organized structure that is compartmentalized by the underlying membrane cytoskeleton (MSK). The interaction between the cell membrane and the cytoskeleton led to the "picket-fence" model, which was proposed to explain certain aspects of membrane compartmentalization. This model assumes that the MSK hinders and confines the motion of receptors and lipids to compartments in the membrane. However, the impact of the MSK on receptor clustering, aggregation, and downstream signaling remains unclear. For example, some evidence suggests that the MSK enhances dimerization, while other evidence suggests decreased dimerization and signaling. Herein, we use computational Monte Carlo simulations to examine the effects of MSK density and receptor concentration on receptor dimerization and clustering. Preliminary results suggest that the MSK may have the potential to induce receptor clustering, which is a function of both picket-fence density and receptor concentration.

REEP2 enhances sweet receptor function by recruitment to lipid rafts

Heterologously expressed sensory receptors generally do not achieve the ligand sensitivity observed in vivo, and may require specific accessory proteins to ensure optimal function. We searched for taste cell-expressed receptor transporting protein (RTP) and receptor expression enhancing protein (REEP) family members that might serve as accessory molecules to enhance gustatory receptor function. We determined that REEP2 is an integral membrane protein expressed in taste cells, physically associates with both subunits of the type 1 taste receptor 2 and type 1 taste receptor 3 sweet receptor and specifically enhances responses to tastants of heterologously expressed sweet and bitter taste receptors. Downregulation of endogenously expressed REEP2 in the chemosensory enteroendocrine GLUTag cell line dramatically reduced sensitivity of endogenous sweet receptors. In contrast to the observation that RTP1, RTP2, and REEP1 enhance function of olfactory receptors by promoting their transit to the cell surface, we found that REEP2 does not increase cell surface expression of sweet receptors but instead alters their spatial organization. REEP2 recruits sweet receptors into lipid raft microdomains localized near the taste cell's apical region, thereby improving G-protein-coupled receptor signaling and promoting receptor access to tastants arriving through the apical taste pore.

Long-lasting target binding and rebinding as mechanisms to prolong in vivo drug action

An increasing number of examples in the literature suggest that the in vivo duration of drug action not only depends on macroscopic pharmacokinetic properties like plasma half-life and the time needed to equilibrate between the plasma and the effect compartments, but is also influenced by long-lasting target binding and rebinding. The present review combines information from different research areas and simulations to explore the nature of these mechanisms and the conditions in which they are most prevalent. Simulations reveal that these latter phenomena become especially influential when there is no longer sufficient free drug around to maintain high levels of receptor occupancy. There is not always a direct link between slow dissociation and long-lasting in vivo target protection, as the rate of free drug elimination from the effect compartment is also a key influencing factor. Local phenomena that hinder the diffusion of free drug molecules away from their target may allow them to consecutively bind to the same target and/or targets nearby (denoted as 'rebinding') even when their concentration in the bulk phase has already dropped to insignificant levels. The micro-anatomic properties of many effect compartments are likely to intensify this phenomenon. By mimicking the complexity of tissues, intact cells offer the opportunity to investigate both mechanisms under the same, physiologically relevant conditions.

Quantifying the proteolytic release of extracellular matrix-sequestered VEGF with a computational model

Background

VEGF proteolysis by plasmin or matrix metalloproteinases (MMPs) is believed to play an important role in regulating vascular patterning in vivo by releasing VEGF from the extracellular matrix (ECM). However, a quantitative understanding of the kinetics of VEGF cleavage and the efficiency of cell-mediated VEGF release is currently lacking. To address these uncertainties, we develop a molecular-detailed quantitative model of VEGF proteolysis, used here in the context of an endothelial sprout.

Methodology and Findings

To study a cell's ability to cleave VEGF, the model captures MMP secretion, VEGF-ECM binding, VEGF proteolysis from VEGF₁₆₅ to VEGF₁₁₄ (the expected MMP cleavage product of VEGF₁₆₅) and VEGF receptor-mediated recapture. Using experimental data, we estimated the effective bimolecular rate constant of VEGF₁₆₅ cleavage by plasmin to be 328 M⁻¹s⁻¹ at 25°C, which is relatively slow compared to typical MMP-ECM proteolysis reactions. While previous studies have implicated cellular proteolysis in growth factor processing, we show that single cells do not individually have the capacity to cleave VEGF to any appreciable extent (less than 0.1% conversion). In addition, we find that a tip cell's receptor system will not efficiently recapture the cleaved VEGF due to an inability of cleaved VEGF to associate with Neuropilin-1.

Conclusions

Overall, VEGF₁₆₅ cleavage in vivo is likely to be mediated by the combined effect of numerous cells, instead of behaving in a single-cell-directed, autocrine manner. We show that heparan sulfate proteoglycans (HSPGs) potentiate VEGF cleavage by increasing the VEGF clearance time in tissues. In addition, we find that the VEGF-HSPG complex is more sensitive to proteases than is soluble VEGF, which may imply its potential relevance in receptor signaling. Finally, according to our calculations, experimentally measured soluble protease levels are approximately two orders of magnitude lower than that needed to reconcile levels of VEGF cleavage seen in pathological situations.

Localization and socialization: experimental insights into the functional architecture of IP3 receptors

Inositol 1,4,5-trisphosphate (IP(3))-evoked Ca(2+) signals display great spatiotemporal malleability. This malleability depends on diversity in both the cellular organization and in situ functionality of IP(3) receptors (IP(3)Rs) that regulate Ca(2+) release from the endoplasmic reticulum (ER). Recent experimental data imply that these considerations are not independent, such that-as with other ion channels-the local organization of IP(3)Rs impacts their functionality, and reciprocally IP(3)R activity impacts their organization within native ER membranes. Here, we (i) review experimental data that lead to our understanding of the "functional architecture" of IP(3)Rs within the ER, (ii) propose an updated terminology to span the organizational hierarchy of IP(3)Rs observed in intact cells, and (iii) speculate on the physiological significance of IP(3)R socialization in Ca(2+) dynamics, and consequently the emerging need for modeling studies to move beyond gridded, planar, and static simulations of IP(3)R clustering even over short experimental timescales.

Molecular mechanisms for the persistent bronchodilatory effect of the beta 2-adrenoceptor agonist salmeterol

BACKGROUND

Beta(2)-adrenoceptor agonists are effective bronchodilators. In vitro studies demonstrated long-lasting airway smooth muscle relaxation by salmeterol after washout, the quick disappearance of this effect in presence of antagonists and its recovery after antagonist removal. Current explanations invoke salmeterol accumulation in the membrane ('diffusion microkinetic' model) or the existence of salmeterol-binding 'exosites'. An alternative model based on 'rebinding' of a dissociated ligand to the receptor molecules also produces an apparent decrease in the ligand's dissociation rate in the absence of competing ligands. PURPOSE AND APPROACH: Computer-assisted simulations were performed to follow the receptor-occupation by a salmeterol-like ligand and a competing ligand as a function of time. The aptness of the models to describe the above in vitro findings was evaluated.

KEY RESULTS

The 'diffusion microkinetic' model is sufficient to explain a long-lasting beta(2)-adrenoceptor stimulation and reassertion as long as the membrane harbors a high concentration of the agonist. At lower concentration, 'rebinding' and, in second place, 'exosite' binding are likely to become operational.

CONCLUSIONS AND IMPLICATIONS

The 'rebinding' and 'exosite' binding mechanisms take place at a sub-cellular/molecular scale. Pending their demonstration by experiments on appropriate, simple models such as intact cells or membranes thereof, these mechanisms remain hypothetical in the case of salmeterol. Airway smooth muscle contraction could also be governed by additional mechanisms that are particular to this macroscopic approach.

Unravelling the complex dissociation of [(3)H]-rimonabant from plated CB(1) cannabinoid receptor-expressing cells

The dissociation profile of the antagonist [(3)H]-rimonabant from recombinant CB(1) cannabinoid receptors expressed in plated HEK293 cells followed a complex pattern when measured in medium only. After a rapid decline, the specific binding levelled off at about 20% below the initial value. To unravel the responsible mechanism(s), we examined the relative contribution of binding to cells and walls of the culture wells respectively. Washout was also performed in the presence of an excess of unlabelled ligand and/or bovine serum albumin (BSA). The findings suggest that dissociated [(3)H]-rimonabant molecules not only undergo rebinding to the same or neighbouring receptors but also partition in the cell membranes and fix to the walls. As these non-receptor associations still occur in presence of unlabelled ligand, they can be erroneously regarded to represent 'specific binding'. While the unlabelled ligand was most effective in preventing receptor rebinding, BSA was most effective in preventing non-receptor associations. To measure receptor-dissociation only, washout is best performed in presence of unlabelled ligand and BSA or any other protein that can pick-up free radioligand molecules. Yet, washout in medium only could hint at mechanisms that affect the in vivo residence time of the drug in question.

A catalytic role of heparin within the extracellular matrix

We investigated the mechanism by which heparin enhances the binding of vascular endothelial growth factor (VEGF) to the extracellular matrix protein fibronectin. In contrast to other systems, where heparin acts as a protein scaffold, we found that heparin functions catalytically to modulate VEGF binding site availability on fibronectin. By measuring the binding of VEGF and heparin to surface-immobilized fibronectin, we show that substoichiometric amounts of heparin exposed cryptic VEGF binding sites within fibronectin that remain available after heparin removal. Measurement of association and dissociation kinetics for heparin binding to fibronectin indicated that the interaction is rapid and transient. We localized the heparin-responsive element to the C-terminal 40-kDa Hep2 domain of fibronectin. A mathematical model of this catalytic process was constructed that supports a mechanism whereby the heparin-induced conformational change in fibronectin is accompanied by release of heparin. Experiments with endothelial extracellular matrix suggest that this process may also occur within biological matrices. These results indicate a novel mechanism whereby heparin catalyzes the conversion of fibronectin to an open conformation by transiently interacting with fibronectin and progressively hopping from molecule to molecule. Catalytic activation of the extracellular matrix might be an important mechanism for heparin to regulate function during normal and disease states.

Allovalency revisited: an analysis of multisite phosphorylation and substrate rebinding

The utilization of multiple phosphorylation sites in regulating a biological response is ubiquitous in cell signaling. If each site contributes an additional, equivalent binding site, then one consequence of an increase in the number of phosphorylations may be to increase the probability that, upon dissociation, a ligand immediately rebinds to its receptor. How such effects may influence cell signaling systems is not well understood. Here, a self-consistent integral equation formalism for ligand rebinding, in conjunction with Monte Carlo simulations, is employed to further investigate the effects of multiple, equivalent binding sites on shaping biological responses. Multiple regimes that characterize qualitatively different physics due to the differential prevalence of rebinding effects are predicted. Calculations suggest that when ligand rebinding contributes significantly to the dose response, a purely allovalent model can influence the binding curves nonlinearly. The model also predicts that ligand rebinding in itself appears insufficient to generate a highly cooperative biological response.

What is the biological relevance of the specific bond properties revealed by single-molecule studies?

During the last decade, many authors took advantage of new methodologies based on atomic force microscopy (AFM), biomembrane force probes (BFPs), laminar flow chambers or optical traps to study at the single-molecule level the formation and dissociation of bonds between receptors and ligands attached to surfaces. Experiments provided a wealth of data revealing the complexity of bond response to mechanical forces and the dependence of bond rupture on bond history. These results supported the existence of multiple binding states and/or reaction pathways. Also, single bond studies allowed us to monitor attachments mediated by a few bonds. The aim of this review is to discuss the impact of this new information on our understanding of biological molecules and phenomena. The following points are discussed: (i) which parameters do we need to know in order to predict the behaviour of an encounter between receptors and ligands, (ii) which information is actually yielded by single-molecule studies and (iii) is it possible to relate this information to molecular structure?

Antagonist-radioligand binding to D2L-receptors in intact cells

D(2)-dopamine receptors mediate most of the physiological actions of dopamine and are important recognition sites for antipsychotic drugs. Earlier binding studies were predominantly done with broken cell preparations with the tritiated D(2)-receptor antagonists [(3)H]-raclopride, a hydrophilic benzamide, and [(3)H]-spiperone, a highly hydrophobic butyrophenone. Here we compared [(3)H]-raclopride and [(3)H]-spiperone binding properties in intact Chinese Hamster Ovary cells stably expressing recombinant human D(2L)-receptors. Specific binding of both radioligands occurred to a comparable number of sites. In contrast to the rapid dissociation of [(3)H]-raclopride in both medium only and in the presence of an excess of unlabelled ligand [(3)H]-spiperone dissociation was only observed in the latter condition, and it was still slower than in broken cell preparations. However, this could not explain the pronounced difference in the potency of some unlabelled ligands to compete with both radioligands. To integrate these new findings, a model is proposed in which raclopride approaches the receptor from the aqueous phase, while spiperone approaches the receptor by lateral diffusion within the membrane.

Self-consistent theory of reversible ligand binding to a spherical cell

Ligand binding to receptors is the initial event in many signaling processes, and a quantitative understanding of this interaction is important for modeling cell behavior. In this paper, we study the kinetics of reversible ligand binding to receptors on a spherical cell surface using a self-consistent stochastic theory. Binding, dissociation, diffusion and rebinding of ligands are incorporated into the theory in a systematic manner. We derive explicitly the time evolution of the ligand-bound receptor fraction p(t) in various regimes. Contrary to the commonly accepted view, we find that the well-known Berg-Purcell scaling for the association rate is modified as a function of time. Specifically, the effective on-rate changes non-monotonically as a function of time and equals the intrinsic rate at very early as well as late times, while being approximately equal to the Berg-Purcell value at intermediate times. The effective dissociation rate, as it appears in the binding curve or measured in a dissociation experiment, is strongly modified by rebinding events and assumes the Berg-Purcell value except at very late times, where the decay is algebraic and not exponential. In equilibrium, the ligand concentration everywhere in the solution is the same and equals its spatial mean, thus ensuring that there is no depletion in the vicinity of the cell. Implications of our results for binding experiments and numerical simulations of ligand-receptor systems are also discussed.