Interaction Between SARS-CoV-2 Spike Protein S1 Subunit and Oyster Heat Shock Protein 70

There is growing evidence that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) contaminates the marine environment and is bioaccumulated in filter-feeding shellfish. Previous study shows the Pacific oyster tissues can bioaccumulate the SARS-CoV-2, and the oyster heat shock protein 70 (oHSP70) may play as the primary attachment receptor to bind SARS-CoV-2's recombinant spike protein S1 subunit (rS1). However, detailed information about the interaction between rS1 and oHSP70 is still unknown. In this study, we confirmed that the affinity of recombinant oHSP70 (roHSP70) for rS1 (KD = 20.4 nM) is comparable to the receptor-binding affinity of rACE2 for rS1 (KD = 16.7 nM) by surface plasmon resonance (SPR)-based Biacore and further validated by enzyme-linked immunosorbent assay (ELISA). Three truncated proteins (roHSP70-N/C/M) and five mutated proteins (p.I229del, p.D457del, p.V491_K495del, p.K556I, and p.ΣroHSP70) were constructed according to the molecular docking results. All three truncated proteins have significantly lower affinity for rS1 than the full-length roHSP70, indicating that all three segments of roHSP70 are involved in binding to rS1. Further, the results of SPR and ELISA showed that all five mutant proteins had significantly lower affinity for rS1 than roHSP70, suggesting that amino acids at these sites are involved in binding to rS1. This study provides a preliminary theoretical basis for the bioaccumulation of SARS-CoV-2 in oyster tissues or using roHSP70 as the capture unit to selectively enrich virus particles for detection.

Reconstitution of cytolinker-mediated crosstalk between actin and vimentin

Cell shape and motility are determined by the cytoskeleton, an interpenetrating network of actin filaments, microtubules, and intermediate filaments. The biophysical properties of each filament type individually have been studied extensively by cell-free reconstitution. By contrast, the interactions between the three cytoskeletal networks are relatively unexplored. They are coupled via crosslinkers of the plakin family such as plectin. These are challenging proteins for reconstitution because of their giant size and multidomain structure. Here we engineer a recombinant actin-vimentin crosslinker protein called 'ACTIF' that provides a minimal model system for plectin, recapitulating its modular design with actin-binding and intermediate filament-binding domains separated by a coiled-coil linker for dimerisation. We show by fluorescence and electron microscopy that ACTIF has a high binding affinity for vimentin and actin and creates mixed actin-vimentin bundles. Rheology measurements show that ACTIF-mediated crosslinking strongly stiffens actin-vimentin composites. Finally, we demonstrate the modularity of this approach by creating an ACTIF variant with the intermediate filament binding domain of Adenomatous Polyposis Coli. Our protein engineering approach provides a new cell-free system for the biophysical characterization of intermediate filament-binding crosslinkers and for understanding the mechanical synergy between actin and vimentin in mesenchymal cells.

Protein-Protein Binding Kinetics by Biolayer Interferometry

The specific kinetics and thermodynamics of protein-protein interactions underlie the molecular mechanisms of cellular functions; hence the characterization of these interaction parameters is central to the quantitative understanding of physiological and pathological processes. Many methods have been developed to study protein-protein interactions, which differ in various features including the interaction detection principle, the sensitivity, whether the method operates in vivo, in vitro, or in silico, the temperature control, the use of labels, immobilization, the amount of sample required, the number of measurements that can be accomplished simultaneously, or the cost. Bio-Layer Interferometry (BLI) is a label-free biophysical method to measure the kinetics of protein-protein interactions. Label-free interaction assays are a broad family of methods that do not require protein modifications (other than immobilization) or labels such as fusions with fluorescent proteins or transactivating domains or chemical modifications like biotinylation or reaction with radionuclides. Besides BLI, other label-free techniques that are widely used for determining protein-protein interactions include surface plasmon resonance (SPR), thermophoresis, and isothermal titration calorimetry (ITC), among others.

SARS-CoV-2S-Protein-Ace2 Binding Analysis Using Surface Plasmon Resonance

Surface plasmon resonance (SPR) allows for the label-free determination of the binding affinity and rate constants of bimolecular interactions. Here, we describe the method used for the analysis of the Ace2-SARS-CoV2 S-protein interaction using indirect capture of the S-protein onto the SPR surface, and flowing monomeric Ace2. This method will allow for the determination of the rate constants for affinity, with additional analysis that is achievable using S-protein capture levels in conjunction with the sensorgram response for relative activity benchmarking.

Real-time tracking of drug binding to influenza A M2 reveals a high energy barrier

The drug Rimantadine binds to two different sites in the M2 protein from influenza A, a peripheral site and a pore site that is the primary site of efficacy. It remained enigmatic that pore binding did not occur in certain detergent micelles, and in particular incomplete binding was observed in a mixture of lipids selected to match the viral membrane. Here we show that two effects are responsible, namely changes in the protein upon pore binding that prevented detergent solubilization, and slow binding kinetics in the lipid samples. Using 55-100 kHz magic-angle spinning NMR, we characterize kinetics of drug binding in three different lipid environments: DPhPC, DPhPC with cholesterol and viral mimetic membrane lipid bilayers. Slow pharmacological binding kinetics allowed the characterization of spectral changes associated with non-specific binding to the protein periphery in the kinetically trapped pore-apo state. Resonance assignments were determined from a set of proton-detected 3D spectra. Chemical shift changes associated with functional binding in the pore of M2 were tracked in real time in order to estimate the activation energy. The binding kinetics are affected by pH and the lipid environment and in particular cholesterol. We found that the imidazole-imidazole hydrogen bond at residue histidine 37 is a stable feature of the protein across several lipid compositions. Pore binding breaks the imidazole-imidazole hydrogen bond and limits solubilization in DHPC detergent.

Parameter estimation and identifiability analysis for a bivalent analyte model of monoclonal antibody-antigen binding

Surface plasmon resonance (SPR) is an extensively used technique to characterize antigen-antibody interactions. Affinity measurements by SPR typically involve testing the binding of antigen in solution to monoclonal antibodies (mAbs) immobilized on a chip and fitting the kinetics data using 1:1 Langmuir binding model to derive rate constants. However, when it is necessary to immobilize antigens instead of the mAbs, a bivalent analyte (1:2) binding model is required for kinetics analysis. This model is lacking in data analysis packages associated with high throughput SPR instruments and the packages containing this model do not explore multiple local minima and parameter identifiability issues that are common in non-linear optimization. Therefore, we developed a method to use a system of ordinary differential equations for analyzing 1:2 binding kinetics data. Salient features of this method include a grid search on parameter initialization and a profile likelihood approach to determine parameter identifiability. Using this method we found a non-identifiable parameter in data set collected under the standard experimental design. A simulation-guided improved experimental design led to reliable estimation of all rate constants. The method and approach developed here for analyzing 1:2 binding kinetics data will be valuable for expeditious therapeutic antibody discovery research.

An integrated modelling approach for targeted degradation: insights on optimization, data requirements and PKPD predictions from semi- or fully-mechanistic models and exact steady state solutions

The value of an integrated mathematical modelling approach for protein degraders which combines the benefits of traditional turnover models and fully mechanistic models is presented. Firstly, we show how exact solutions of the mechanistic models of monovalent and bivalent degraders can provide insight on the role of each system parameter in driving the pharmacological response. We show how on/off binding rates and degradation rates are related to potency and maximal effect of monovalent degraders, and how such relationship can be used to suggest a compound optimization strategy. Even convoluted exact steady state solutions for bivalent degraders provide insight on the type of observations required to ensure the predictive capacity of a mechanistic approach. Specifically for PROTACs, the structure of the exact steady state solution suggests that the total remaining target at steady state, which is easily accessible experimentally, is insufficient to reconstruct the state of the whole system at equilibrium and observations on different species (such as binary/ternary complexes) are necessary. Secondly, global sensitivity analysis of fully mechanistic models for PROTACs suggests that both target and ligase baselines (actually, their ratio) are the major sources of variability in the response of non-cooperative systems, which speaks to the importance of characterizing their distribution in the target patient population. Finally, we propose a pragmatic modelling approach which incorporates the insights generated with fully mechanistic models into simpler turnover models to improve their predictive ability, hence enabling acceleration of drug discovery programs and increased probability of success in the clinic.

Conformational antigenic heterogeneity as a cause of the persistent fraction in HIV-1 neutralization

BACKGROUND

Neutralizing antibodies (NAbs) protect against HIV-1 acquisition in animal models and show promise in treatment of infection. They act by binding to the viral envelope glycoprotein (Env), thereby blocking its receptor interactions and fusogenic function. The potency of neutralization is largely determined by affinity. Less well explained is the persistent fraction, the plateau of remaining infectivity at the highest antibody concentrations.

RESULTS

We observed different persistent fractions for neutralization of pseudovirus derived from two Tier-2 isolates of HIV-1, BG505 (Clade A) and B41 (Clade B): it was pronounced for B41 but not BG505 neutralization by NAb PGT151, directed to the interface between the outer and transmembrane subunits of Env, and negligible for either virus by NAb PGT145 to an apical epitope. Autologous neutralization by poly- and monoclonal NAbs from rabbits immunized with soluble native-like B41 trimer also left substantial persistent fractions. These NAbs largely target a cluster of epitopes lining a hole in the dense glycan shield of Env around residue 289. We partially depleted B41-virion populations by incubating them with PGT145- or PGT151-conjugated beads. Each depletion reduced the sensitivity to the depleting NAb and enhanced it to the other. Autologous neutralization by the rabbit NAbs was decreased for PGT145-depleted and enhanced for PGT151-depleted B41 pseudovirus. Those changes in sensitivity encompassed both potency and the persistent fraction. We then compared soluble native-like BG505 and B41 Env trimers affinity-purified by each of three NAbs: 2G12, PGT145, or PGT151. Surface plasmon resonance showed differences among the fractions in antigenicity, including kinetics and stoichiometry, congruently with the differential neutralization. The large persistent fraction after PGT151 neutralization of B41 was attributable to low stoichiometry, which we explained structurally by clashes that the conformational plasticity of B41 Env causes.

CONCLUSION

Distinct antigenic forms even of clonal HIV-1 Env, detectable among soluble native-like trimer molecules, are distributed over virions and may profoundly mold neutralization of certain isolates by certain NAbs. Affinity purifications with some antibodies may yield immunogens that preferentially expose epitopes for broadly active NAbs, shielding less cross-reactive ones. NAbs reactive with multiple conformers will together reduce the persistent fraction after passive and active immunization.

Chloroquine-analogues block anthrax protective antigen channels in steady-state and kinetic studies

The tripartite anthrax toxin from Bacillus anthracis represents the prototype of A-B type of toxins, where the effector A (an enzymatic subunit) is transported with the help of a binding component B into a target cell. Anthrax toxin consists of three different molecules, two effectors, lethal factor (LF) and edema factor (EF) and the binding component also known as protective antigen (PA). PA forms heptamers or octamers following binding to host cell's receptors and mediates the translocation of the effectors into the cytosol via the endosomal pathway. The cation-selective PA₆₃-channel is able to reconstitute in lipid membranes and can be blocked by chloroquine and other heterocyclic compounds. This suggests that the PA₆₃-channel contains a binding site for quinolines. In this study, we investigated the structure-function relationship of different quinolines for the block of the PA₆₃-channel. The affinity of the different chloroquine analogues to the PA₆₃-channel as provided by the equilibrium dissociation constant was measured using titrations. Some quinolines had a much higher affinity to the PA₆₃-channel than chloroquine itself. We also performed ligand-induced current noise measurements using fast Fourier transformation to get insight in the kinetics of the binding of some quinolines to the PA₆₃-channel. The on-rate constants of ligand binding were around 10⁸M^-1·s^-1 at 150mM KCl and were only little dependent on the individual quinoline. The off-rates varied between 4s^-1 and 160s^-1 and depended much more on the structure of the molecules than the on-rate constants. The possible use of the 4-aminoquinolines as a therapy is discussed.

The Importance of Binding Kinetics and Drug-Target Residence Time in Pharmacology

A dominant assumption of pharmacology throughout the 20^th century has been that in vivo target occupancy - and attendant pharmacodynamics - depends on the systemic concentration of drug relative to the equilibrium dissociation constant for the drug-target complex. In turn, the duration of pharmacodynamics is temporally linked to the systemic pharmacokinetics of the drug. Yet, there are many examples of drugs for which pharmacodynamic effect endures long after the systemic concentration of drug has waned to (equilibrium) insignificant levels. To reconcile such data, the drug-target residence time model was formulated, positing that it is the lifetime (or residence time) of the binary drug-target complex, and not its equilibrium affinity per se, that determines the extent and duration of drug pharmacodynamics. Here we review this model, its evolution over time, and its applications to natural ligand-macromolecule biology and synthetic drug-target pharmacology.

Fast Protocols for Characterizing Antibody-Peptide Binding

Microarray assay formats gained popularity in the 1990s, first implemented in DNA-based arrays but later adopted for use with proteins, namely antibodies, peptides, low molecular weight (LMW) molecules, such as lipids, and even tissues. In nucleic acid-based affinity assays and arrays, but not in protein or peptide arrays, the specificity and affinity of complementary strand interactions can be deduced from or adjusted through modifications to the nucleotide sequence. Arrays of LMW molecules are characterized by largely uniform but low binding affinities. Multiplexed protein-based affinity assays, such as microarrays, might present an additional challenge due to heterogeneity of antigen properties and of their binding affinities. The use of peptides instead of proteins reduces physical heterogeneity of these reagents through either the widened peptide selection options or rational sequence engineering. However, rational engineering of binding affinities remains an unmet need, and peptide-binding affinities to the respective antipeptide antibodies could vary by orders of magnitude. Hence, multiplexing of such assays by using a microarray format and data analysis and interpretation requires some knowledge of their binding affinities. Low-throughput binding assays to characterize such peptide-antipeptide antibodies interactions are widely available, but scaling-up of traditional protein- and peptide-binding assays might present practical challenges. Here, we describe fast label-free practical approach especially suitable for estimating peptide-binding affinities. The method in question relies on commercially available biolayer interferometry-based equipment with a protocol which can be easily scaled-up, subject to user needs and equipment availability.

Humoral immune response characterization of heterologous prime-boost vaccination with CoronaVac and BNT162b2

Background

Vaccination against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has proven to be a successful strategy for prevent severe infections. CoronaVac and BNT162b2 are the most used vaccines worldwide, but their use in heterologous vaccination schedules is still subjected to evaluation.

Methods

Fifty healthy individuals who received heterologous prime-boost vaccination with CoronaVac and BNT162b2 were enrolled in a post-vaccination serological follow-up longitudinal prospective study. We evaluated specific serum anti-receptor binding domain (RBD) IgG antibody levels, and their capacity to block RBD-ACE2 interaction with a surrogate neutralization assay. In 20 participants, we assessed antibody binding kinetics by surface plasmon resonance, and Fc-mediated functions by ADCC and ADCP reporter assays.

Results

Our baseline seronegative cohort, displayed seroconversion after two doses of CoronaVac and an important decrease in serum anti-RBD IgG antibodies levels 80 days post-second dose. These levels increased significantly early after the third dose with BNT162b2, but 73 days after the booster we found a new fall. Immunoglobulin functionalities showed a similar behavior.

Conclusions

The heterologous prime-boost vaccination with CoronaVac and BNT162b2 generated an impressive increase in serum anti-RBD specific antibody levels followed by a drop. Nevertheless, these titers remained well above those found in individuals only vaccinated with CoronaVac in the same elapsed time. Serum IgG levels showed high correlation with antibody binding analysis, their capacity to block RBD-ACE2 interaction, and Fc-effectors mechanisms. Our work sheds light on the humoral immune response to heterologous vaccination with CoronaVac and BNT162b2, to define a post-vaccination correlate of protection against SARS-CoV-2 infection and to discuss the scheduling of future vaccine boosters in general population.

Toll-like receptors change morphine-induced antinociception, tolerance and dependence: Studies using male and female TLR and signalling gene KO mice

Toll-like receptors (TLR) have been proposed as a site of action that alters opioid pharmacodynamics. However, a comprehensive assessment of acute opioid antinociception, tolerance and withdrawal behaviours in genetic null mutant strains with altered innate immune signalling has not been performed. Nor has the impact of genetic deletion of TLR2/4 on high-affinity opioid receptor binding. Here we show that diminished TLR signalling potentiates acute morphine antinociception equally in male and female mice. However, only male TIR8 null mutant mice showed reduced morphine analgesia. Analgesic tolerance was prevented in TLR2 and TLR4 null mutants, but not MyD88 animals. Withdrawal behaviours were only protected in TLR2^-/- mice. In silico docking simulations revealed opioid ligands bound preferentially to the LPS binding pocket of MD-2 rather than TLR4. There was no binding of [³H](-)-naloxone or [³H]diprenorphine to TLR4 in the concentrations explored. These data confirm that opioids have high efficacy activity at innate immune pattern recognition binding sites but do not bind to TLR4 and identify critical pathway and sex-specific effects of the complex innate immune signalling contributions to opioid pharmacodynamics. These data further support the behavioural importance of the TLR-opioid interaction but fail to demonstrate direct evidence for high-affinity binding of the TLR4 signalling complex to ligands.

Kinetic profiling and functional characterization of 8-phenylxanthine derivatives as A2B adenosine receptor antagonists

A_2B adenosine receptor (A_2BAR) antagonists have therapeutic potential in inflammation-related diseases such as asthma, chronic obstructive pulmonary disease and cancer. However, no drug is currently clinically approved, creating a demand for research on novel antagonists. Over the last decade, the study of target binding kinetics, along with affinity and potency, has been proven valuable in early drug discovery stages, as it is associated with improved in vivo drug efficacy and safety. In this study, we report the synthesis and biological evaluation of a series of xanthine derivatives as A_2BAR antagonists, including an isothiocyanate derivative designed to bind covalently to the receptor. All 28 final compounds were assessed in radioligand binding experiments, to evaluate their affinity and for those qualifying, kinetic binding parameters. Both structure-affinity and structure-kinetic relationships were derived, providing a clear relationship between affinity and dissociation rate constants. Two structurally similar compounds, 17 and 18, were further evaluated in a label-free assay due to their divergent kinetic profiles. An extended cellular response was associated with long A_2BAR residence times. This link between a ligand's A_2BAR residence time and its functional effect highlights the importance of binding kinetics as a selection parameter in the early stages of drug discovery.

Electrochemically differentiated human MSCs biosensing platform for quantification of nestin and β-III tubulin as whole-cell system

Polydimethylsiloxane (PDMS) on ITO substrate was used to create a well with conducting surface to adhere human mesenchymal stem cells (hMSCs) and provide electrochemical stimulation for inducing their differentiation into neural-like cells. The cells that received electrochemical stimulation did not show any noticeable change in their viability and proliferation. The cell morphology of the differentiated hMSCs adherent on ITO showed outgrowth and elongation in one dimension, resembling neural-like cells. Immunocytochemistry assessment by quantifying the expression of nestin and β-III tubulin also confirmed the differentiation of hMSCs. These differentiated hMSCs adherent on ITO were used as electrochemical biosensing platform for differential pulse voltammetry (DPV) measurement for selectively quantifying cell surface markers expressed by neural stem cells and mature neurons. The variation of nestin antibodies concentrations from 9 μU to 27 μU showed a linear increase in DPV current with a detection sensitivity of ∼28 nA/μU of antibody. Varying concentrations of β-III tubulin antibodies from 30 μU to 210 μU showed a linear increase in DPV current with a detection sensitivity of ∼2.0 nA/μU of antibody. The highest expression level of cell surface marker corresponding to β-III tubulin in total adherent cells on ITO was calculated. It was in the order of 10^-8 U of antibodies/cell, representing the total population of mature neuron cells. This new way of detection may rapidly assess the quantitative expression of cell surface markers/antigens.

Functional organization of RNA polymerase II in nuclear subcompartments

Distinct clusters of RNA polymerase II are responsible for gene transcription inside eukaryotic cell nuclei. Despite the functional implications of such subnuclear organization, the attributes of these clusters and the mechanisms underlying their formation remain only partially understood. Recently, the concept of proteins and RNA phase-separating into liquid-like droplets was proposed to drive the formation of transcriptionally-active subcompartments. Here, we attempt to reconcile previous with more recent findings, and discuss how the different ways of assembling the active RNA polymerase II transcriptional machinery relate to nuclear compartmentalization.

vCOMBAT: a novel tool to create and visualize a computational model of bacterial antibiotic target-binding

Background

As antibiotic resistance creates a significant global health threat, we need not only to accelerate the development of novel antibiotics but also to develop better treatment strategies using existing drugs to improve their efficacy and prevent the selection of further resistance. We require new tools to rationally design dosing regimens from data collected in early phases of antibiotic and dosing development. Mathematical models such as mechanistic pharmacodynamic drug-target binding explain mechanistic details of how the given drug concentration affects its targeted bacteria. However, there are no available tools in the literature that allow non-quantitative scientists to develop computational models to simulate antibiotic-target binding and its effects on bacteria.

Results

In this work, we have devised an extension of a mechanistic binding-kinetic model to incorporate clinical drug concentration data. Based on the extended model, we develop a novel and interactive web-based tool that allows non-quantitative scientists to create and visualize their own computational models of bacterial antibiotic target-binding based on their considered drugs and bacteria. We also demonstrate how Rifampicin affects bacterial populations of Tuberculosis bacteria using our vCOMBAT tool.

Conclusions

The vCOMBAT online tool is publicly available at https://combat-bacteria.org/.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04536-3.

Predicting Residence Time of GPCR Ligands with Machine Learning

Drug-target residence time, the duration of binding at a given protein target, has been shown in some protein families to be more significant for conferring efficacy than binding affinity. To carry out efficient optimization of residence time in drug discovery, machine learning models that can predict that value need to be developed. One of the main challenges with predicting residence time is the paucity of data. This chapter outlines all of the currently available ligand kinetic data, providing a repository that contains the largest publicly available source of GPCR-ligand kinetic data to date. To help decipher the features of kinetic data that might be beneficial to include in computational models for the prediction of residence time, the experimental evidence for properties that influence residence time are summarized. Finally, two different workflows for predicting residence time with machine learning are outlined. The first is a single-target model trained on ligand features; the second is a multi-target model trained on features generated from molecular dynamics simulations.

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.

Charge Sensitive Optical Detection for Measurement of Small-Molecule Binding Kinetics

Charge sensitive optical detection (CSOD) technique is a label-free method for real-time measurement of molecular interactions. Traditional label-free optical detection techniques mostly measure the mass of a molecule, and they are less sensitive to small molecules. In contrast, CSOD detects the charge of a molecule, where the signal does not diminish with the size of the molecule, thus capable for studying small molecules. In addition, CSOD is compatible with the standard microplate platform, making it suitable for high-throughput screening of drug candidates. In CSOD, an optical fiber functionalized with the probe molecule is dipped into a well of a microplate where an alternate perpendicular electrical field is applied to the fiber, which drives the fiber into oscillation because of the presence of surface charge on the fiber. The binding of the target molecules changes the charge of the fiber, and thus the amplitude and phase of the oscillating fiber, which are precisely measured through tracking of the optical images of the fiber tip.

Quantification of Receptor Association, Dissociation, and NA-Dependent Motility of Influenza A Particles by Biolayer Interferometry

We describe a method for real-time analysis and quantification of influenza A virus (IAV)-receptor interactions by biolayer interferometry (BLI). Biotinylated synthetic sialoglycans or sialoglycoproteins (biotinylated or Fc-tagged) were immobilized on the tip of biosensors (coated with streptavidin or protein A) that were subsequently dipped into IAV particle solutions in 96-well plates. Association and/or dissociation of IAV particles was recorded in consecutive steps in buffers of choice. From the association and dissociation curves, parameters can be derived that describe IAV particle-receptor interactions in absence or presence of neuraminidase activity. Overall, the method provides a quantitative description of the hemagglutinin-neuraminidase balance that determines the interaction kinetics of IAV with specific sialoglycan receptors.

Analytical and functional aspects of protein-ligand interactions: Beyond induced fit and conformational selection

Ligand-dependent changes in protein conformation are foundational to biology. Historical mechanistic models for substrate-specific proteins are induced fit (IF) and conformational selection (CS), which invoke a change in protein conformation after ligand binds or before ligand binds, respectively. These mechanisms have important, but rarely discussed, functional relevance because IF vs. CS can differentially affect a protein's substrate specificity or promiscuity, and its regulatory properties. The modern view of proteins as conformational ensembles in both ligand free and bound states, together with the realization that most proteins exhibit some substrate promiscuity, demands a deeper interpretation of the historical models and provides an opportunity to improve mechanistic analyses. Here we describe alternative analytical strategies for distinguishing the historical models, including the more complex expanded versions of IF and CS. Functional implications of the different models are described. We provide an alternative perspective based on protein ensembles interacting with ligand ensembles that clarifies how a single protein can 'apparently' exploit different mechanisms for different ligands. Mechanistic information about protein ensembles can be optimized when they are probed with multiple ligands.

Comparison of methods for quantitative biomolecular interaction analysis

In order to perform good kinetic experiments, not only the experimental conditions have to be optimized, but the evaluation procedure as well. The focus of this work is the in-depth comparison of different approaches and algorithms to determine kinetic rate constants for biomolecular interaction analysis (BIA). The different algorithms are applied not only to flawless simulated data, but also to real-world measurements. We compare five mathematical approaches for the evaluation of binding curves following pseudo-first-order kinetics with different noise levels. In addition, reflectometric interference spectroscopy (RIfS) measurements of two antibodies are evaluated to determine their binding kinetics. The advantages and disadvantages of the individual approach will be investigated and discussed in detail. In summary, we will raise awareness on how to evaluate and judge results from BIA by using different approaches rather than having to rely on “black box” closed (commercial) software packages.

Kinetic Analysis of a Protein-protein Complex to Determine its Dissociation Constant (KD) and the Effective Concentration (EC50) of an Interplaying Effector Molecule Using Bio-layer Interferometry

Biolayer interferometry (BLI) is an emerging analytical tool that allows the study of protein complexes in real time to determine protein complex kinetic parameters. This article describes a protocol to determine the K_D of a protein complex using a 6×His tagged fusion protein as bait immobilized on the NTA sensor chip of the FortéBio^® Octet K2 System (Sartorius). We also describe how to determine the half maximal effective concentration (EC₅₀, also known as IC₅₀ for inhibiting effectors) of a metabolite. The complete protocol allows the determination of protein complex K_D and small molecular effector EC₅₀ within 8 h, measured in triplicates. Graphic abstract: Principle of the Biolayer interferometry measurement. (Middle, top) Exemplary result of the BLI measurement using Octet^® (Raw Data). Wavelength shift (Δλ) against time. (A) Baseline 1. Baseline measurement. When the sensor is equilibrated in the kinetics buffer. The light is reflected with no difference. (B) Load. The his-tagged proteins (ligand) are loaded onto the sensor surface. The light is reflected with a shift of the wavelength. (C) Baseline 2. The loaded sensor is equilibrated in the kinetics buffer. No further wavelength shift appears. (D) Association. The loaded sensor is dipped into the analyte solution. The analyte binds to the immobilized ligand along with an increased wavelength shift. (E) Dissociation. Afterward, the sensor is dipped again into the kinetics buffer without the analyte. Some analyte molecules dissociate. The wavelength shift decreases. (Subfigures A-E) The left side shows the position of the sensor during the measurement seen in the representative BLI measurement, marked with the figure label. The right side shows the light path in the sensor. Black waves represent the light emitted to the sensor surface. The red waves show the light reflected from the sensor surface back to the detector.

Analyze impact of tumor-associated kinetics on antibody delivery in solid tumors with a physiologically based pharmacokinetics/pharmacodynamics model

Monoclonal antibody (mAb)-based drugs are critical anti-cancer therapies. Unfortunately, therapeutic efficacy can be compromised by spatially heterogeneous intratumoral Ab deposition. Binding-site barriers arising from Ab and tumor-associated kinetics often underlie this phenomenon. Quantitative insight into these issues may lead to more efficient drug delivery. Difficulties in addressing this issue include (1) lack of techniques to quantify critical kinetic events, (2) lack of a pharmacokinetic/pharmacodynamic (PK/PD) model to assess important parameters for specific tumor types, and (3) uncertainty or variability of critical kinetic factors even within a single tumor type. This study developed a mechanism-based PK/PD model to profile heterogeneous distribution of Ab within tumors and tested this model using real-life experimental data. Model simulations incorporating several uncertainties were used to determine how mAb and tumor-associated kinetics influence receptor occupancy. Simulations were also used to predict the potential impact of these findings in preclinical tumor models and human tumors. We found significant differences in tumor-associated kinetics between groups in which mAb therapy was effective versus groups in which it was ineffective. These kinetic differences included rates of tumor-associated antigen (TAA) degradation, TAA expression, apparent flow rates of interstitial fluid, and ratios of Ab-TAA complex internalization to TAA degradation. We found less significant differences in mAb kinetics, including rates of clearance or affinity for target antigens. In conclusion, our mechanism-based PK/PD model suggests that TAA-associated kinetic factors participate more significantly than those associated with the Ab in generating barriers to mAb delivery and distribution in tumors.

Binding kinetics of liposome conjugated E-selectin and P-selectin glycoprotein ligand-1 measured with atomic force microscopy

Various ligand-functionalized liposomes have been developed for targeted therapies. Typically, the binding properties of the ligands and targeted proteins are measured with surface plasmon resonance (SPR), where the proteins are immobilized on a rigid surface. However, the difference of protein-ligand binding kinetics between liposome-conjugated protein and rigid surface-conjugated protein is not fully understood. In this work, the binding kinetics of P-selectin glycoprotein ligand-1 (PSGL-1) and E-selectin conjugated on liposome and on rigid surfaces are investigated with Atomic Force Microscopy (AFM). The results suggest that protein orientation and diffusion on liposomal membrane can alter the binding kinetics of the protein-ligand interaction. Specifically, the association and dissociation rate constant of AFM probe-conjugated E-selectin and glass-conjugated PSGL-1 are measured as 9.32 × 10⁴ M^-1s^-1 and 1.54 s^-1, respectively. While for the liposome-conjugated E-selectin and glass-conjugated PSGL-1, the kinetic constants are measured as 5.00 × 10⁷ M^-1s^-1 and 2.76 s^-1, respectively. Thus, there is an order's magnitude increase of binding affinity (from k_d = 16.51 μM to k_d = 0.06 μM) when protein is attached to liposome compared to attached to a rigid surface. The results might provide better understanding and pave the way for the future design of the ligand-targeted liposomes.

KD determination from time-resolved experiments on live cells with LigandTracer and reconciliation with end-point flow cytometry measurements

Design of next-generation therapeutics comes with new challenges and emulates technology and methods to meet them. Characterizing the binding of either natural ligands or therapeutic proteins to cell-surface receptors, for which relevant recombinant versions may not exist, represents one of these challenges. Here we report the characterization of the interaction of five different antibody therapeutics (Trastuzumab, Rituximab, Panitumumab, Pertuzumab, and Cetuximab) with their cognate target receptors using LigandTracer. The method offers the advantage of being performed on live cells, alleviating the need for a recombinant source of the receptor. Furthermore, time-resolved measurements, in addition to allowing the determination of the affinity of the studied drug to its target, give access to the binding kinetics thereby providing a full characterization of the system. In this study, we also compared time-resolved LigandTracer data with end-point K_D determination from flow cytometry experiments and hypothesize that discrepancies between these two approaches, when they exist, generally come from flow cytometry titration curves being acquired prior to full equilibration of the system. Our data, however, show that knowledge of the kinetics of the interaction allows to reconcile the data obtained by flow cytometry and LigandTracer and demonstrate the complementarity of these two methods.

Accurate quantitative determination of affinity and binding kinetics for tight binding inhibition of xanthine oxidase

The accurate quantitative determination of affinity and binding kinetics (BK) for tight binding inhibition is extraordinary important from both the continuous optimization of compounds, particularly in developing structure-activity relationships (SAR), and the prediction of in vivo target occupancy (TO). Due to the unique properties for tight binding inhibition that the inhibitors are characterized by the ultrahigh-affinity, relatively fast association to the target enzyme combined with extremely slow dissociation of the inhibitor-enzyme binary complex, the classical steady state equilibrium methods are no longer valid. Here, we made several recommendations of how to design the optimal experiments and apply special mathematical calculation approaches to quantitatively evaluate the accurate affinity and BK as the examples of two tight binding inhibitors against the xanthine oxidase (XO), as well as compared the differences in the results calculated from the different data analytical methods and analyzed the influence of these differences on the XO engagement in human. Analysis of the results displayed that the accurate apparent dissociation constant (K_i^*,app) was 0.2 ± 0.06 nM for topiroxotstat and was 0.45 ± 0.2 nM for febuxostat; that on-rate (k_on) was (4.3 ± 1.1) × 10⁶ M^-1s^-1 for topiroxotstat and was(133.3 ± 3.5) × 10⁶ M^-1s^-1 for febuxostat, and off-rate (k_off) was (1.0±0.2) × 10^-5 s^-1 for topiroxotstat and was ≤ 0.16 × 10^-5 s^-1for febuxostat. Moreover, there were significant differences in the K_i^*,app and k_off values estimated using the appropriate specialized methods for tight binding inhibition versus classical steady state equilibrium methods, with the substantial differences of 14-fold and 32-fold reduction for topiroxostat, respectively, and of 9.6-fold and ≥ 213-fold reduction for febuxostat, while the k_on values remain the moderate differences for the two inhibitors. The obvious greater AUC of XO engagement time courses and longer durations of above 70% engagement by the appropriate specialized methods for tight binding inhibition were observed that the results display the differences of 70.1% and 88%, respectively for topiroxostat and of 38.1% and 35.0%, respectively for febuxostat in human liver cell than by classical steady state equilibrium methods. Again, our studies provide several valuable recommendations of the optimal experiment protocols and appropriate analytical approaches for accurately quantitatively assessing the affinity and BK parameters as well as demonstrate the ability of our recommended methods to generate reliable data for tight binding inhibitors against XO.

Impact of the N501Y substitution of SARS-CoV-2 Spike on neutralizing monoclonal antibodies targeting diverse epitopes

The emergence and rapid spread of the B.1.1.7 lineage (VOC-202012/01) SARS-CoV-2 variant has aroused global concern. The N501Y substitution is the only mutation in the interface between the RBD of B.1.1.7 and ACE2, raising concerns that its recognition by neutralizing antibodies may be affected. Here, we assessed the neutralizing activity and binding affinity of a panel of 12 monoclonal antibodies against the wild type and N501Y mutant SARS-CoV-2 pseudovirus and RBD protein, respectively. We found that the neutralization activity and binding affinity of most detected antibodies (10 out of 12) were unaffected, although the N501Y substitution decreased the neutralizing and binding activities of CB6 and increased that of BD-23. These findings could be of value in the development of therapeutic antibodies.

Simultaneous Imaging of Single Protein Size, Charge, and Binding Using A Protein Oscillation Approach

Electrophoresis and Western blot are important tools in protein research for detection and identification of proteins. These traditional techniques separate the proteins based on size and charge differences and identify the proteins by antibody binding. Over the past decade, the emergence of single-molecule techniques has shown great potential in improving the resolution of the traditional protein analysis methods to the single-molecule level. However, such single-molecule techniques measure either size or charge, and it is challenging to measure both at the same time. Recently, we have developed a single-molecule approach to address this problem. We tether the single proteins to a surface with a polymer linker and drive them into oscillation with an electric field. By tracking the electromechanical response of the proteins to the field using an optical imaging method, the size and charge can be obtained simultaneously. Binding of antibodies or ions to the tethered protein also changes the size and charge, which allows us to probe the interactions. This protocol includes fabrication of protein oscillators, configuration of the optical detection system, and analysis of the oscillation signal for quantification of protein size and charge. We wish this protocol will enable researchers to perform comprehensive single-protein analysis on a single platform.

Influence of bacterial culture medium on peptidoglycan binding of cell wall lytic enzymes

The bacteriolysin lysostaphin (Lst) and endolysin PlyPH are potent modular lytic enzymes with activity against clinically-relevant Gram-positive Staphylococcus aureus and Bacillus cereus, respectively. Both enzymes possess an N-terminal catalytic domain and C-terminal binding domain, with the latter conferring significant enzyme specificity. Lst and PlyPH show reduced activity in the presence of bacterial growth-supporting conditions, such as complex media. Here, we hypothesize that Lst and PlyPH bind poorly to their targets in growth media, which may influence their use in antimicrobial applications in the food industry, as therapeutics, and for control of microbial communities. To this end, binding of isolated Lst and PlyPH binding domains to target bacteria was quantified in the presence of three increasingly complex media - phosphate buffered saline (PBS), defined growth medium (AAM) and undefined complex medium (TSB) by surface plasmon resonance (SPR) and flow cytometry. Evaluation of binding kinetics by SPR demonstrated that PlyPH binding was particularly sensitive to medium composition, with 8-fold lower association and 3.4-fold lower dissociation rate constants to B. cereus in TSB compared to PBS. Flow cytometry studies indicated a decrease in the binding-dependent fluorescent populations of S. aureus and B. cereus, for lysostaphin binding domain and PlyPH binding domain, respectively, in TSB compared to PBS. Enzyme binding behavior was consistent with the enzymes' catalytic activity in the three media, thereby suggesting that compromised enzyme binding could be responsible for poor activity in more complex growth media.

Protein Ligand Interactions Using Surface Plasmon Resonance

Surface Plasmon Resonance (SPR) is a powerful biophysical method for characterizing small molecule binding to proteins. Owing to its ability to characterize binary inteactions between warheads and E3 ligases or substrates, SPR is a useful tool for the development of targeted protein degraders. SPR is also an effective method for optimizing linkers and characterizing ternary complex interactions that are mediated by heterobifunctional ligands (Roy et al. ACS Chem Biol 14:361-368, 2019). Recent advances in the throughput of modern instruments have improved the ability of SPR to rapidly triage ligands based on binding kinetics and affinity, making this technique invaluable for driving degrader optimization. This chapter describes the characterization of ligands binding to the Thalidomide Binding Domain of mouse Cereblon (mCRBN-TBD) using the Biacore 8K+.

Isothermal Titration Calorimetry

Calorimetry is a classical biophysical method that by definition measures heat. In isothermal titration calorimetry (ITC), the heat is the result of titrating interacting components together and allows direct determination of the thermodynamics for this process. The measured heat reflects the enthalpy change (ΔH), and the prospect of determining this in biological systems where high-resolution structural information is available has led to the possibility of rational thermodynamics-guided design of ligands. Although there are limitations to this approach due to the participation of solvent in the thermodynamics, ITC has become an established technique in many labs providing a valuable tool with which to quantify protein-protein interactions. With careful use, ITC can also provide additional insights into the binding process or be used in increasingly complex systems and where interaction is coupled to other molecular events.

Luciferase Complementation Approaches to Measure GPCR Signaling Kinetics and Bias

An understanding of the kinetic contributions to G protein-coupled receptor pharmacology and signaling is increasingly important in compound profiling. Nonequilibrium conditions are commonly present in vivo, for example, as the drug competes with dynamic changes in hormone or neurotransmitter concentration for the receptor. Under such conditions individual binding kinetic properties of the ligands can influence duration of action, local ligand concentration, and functional properties such as the degree of insurmountable inhibition. Mapping the kinetic patterns of GPCR signaling events elicited by agonists, rather than a peak response at a single timepoint, is often key to predicting their functional impact. This is also a path to a better understanding of the origins of ligand bias, and whether such ligands demonstrate their effects through selection of distinct GPCR conformations, or via their kinetic properties. Recent developments in complementation approaches, based on a small bright shrimp luciferase Nanoluc, provide a new route to kinetic analysis of GPCR signaling in living cells that is amenable to the throughput required for compound profiling. In the NanoBiT luciferase complementation system, GPCRs and effector proteins are tagged with Nanoluc fragments optimized for their low interacting affinity and stability. The interactions brought about by GPCR recruitment of the effector are reproduced by a rapid and reversible increase in NanoBiT luminescence, in the presence of its substrate furimazine. Here we discuss the methods for optimizing and validating the GPCR NanoBiT assays, and protocols for their application to study endpoint and kinetic aspects of agonist and antagonist pharmacology. We also describe how timecourse families of agonist concentration response curves, derived from a single NanoBiT assay experiment, can be used to evaluate the kinetic components in operational model derived parameters of ligand bias.

Fluorescence Anisotropy-Based Assay for Characterization of Ligand Binding Dynamics to GPCRs: The Case of Cy3B-Labeled Ligands Binding to MC4 Receptors in Budded Baculoviruses

During the past decade, fluorescence methods have become valuable tools for characterizing ligand binding to G protein-coupled receptors (GPCRs). However, only a few of the assays enable studying wild-type receptors and monitor the ligand binding in real time. One of the approaches that is inherently suitable for this purpose is the fluorescence anisotropy (FA) assay. In the FA assay, the change of ligand's rotational freedom connected with its binding to the receptor can be monitored with a conventional fluorescence plate reader equipped with suitable optical filters. To achieve the high receptor concentration required for the assay and the low autofluorescence levels essential for reliable results, budded baculoviruses that display GPCRs on their surfaces can be used. The monitoring process generates a substantial amount of kinetic data, which is usually stored as a proprietary file format limiting the flexibility of data analysis. To solve this problem, we propose the use of the data curation software Aparecium ( http://gpcr.ut.ee/aparecium.html ), which integrates experimental data with metadata in a Minimum Information for Data Analysis in Systems Biology (MIDAS) format. Aparecium enables data export to different software packages for fitting to suitable kinetic or equilibrium models. A combination of the FA assay with the novel data analysis strategy is suitable for screening new active compounds, but also for modeling complex systems of ligand binding to GPCRs. We present the proposed approach using different fluorescent probes and assay types to characterize ligand binding to melanocortin 4 (MC₄) receptor.

Understanding Binding-Induced Folding by Temperature Jump

Temperature jump is a powerful technique for the characterization of fast kinetics and can be readily employed to understand both binding and folding reactions. Here we summarize briefly a temperature-jump prototypical experiment between an intrinsically disordered protein and its physiological partner. The model used is the N_TAIL domain from Measles virus Nucleoprotein and its natural ligand, the globular P_XD domain from Measles virus Phosphoprotein. We recapitulate how to set up the experiment and how to analyze data in order to extract the kinetic parameters of the reaction.

Revisit ligand-receptor interaction at the human vasopressin V2 receptor: A kinetic perspective

The vasopressin V₂ receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target for water balance disorders such as polycystic kidney disease. Traditionally, the discovery of novel agents for the vasopressin V₂ receptor has been guided by evaluating their receptor affinity, largely ignoring the binding kinetics. However, the latter is receiving increasing attention in the drug research community and has been proved to be a more complete descriptor of the dynamic process of ligand-receptor interaction. Herein we aim to revisit the molecular basis of ligand-vasopressin V₂ receptor interaction from the less-investigated kinetic perspective. A homogenous time-resolved fluorescence resonance energy transfer (TR-FRET) assay was set up and optimized, which enabled accurate kinetic profiling of unlabeled vasopressin V₂ receptor ligands. Receptor occupancy profiles of two representative antagonists with distinct target residence time were simulated. Their functional effects were further explored in cAMP assays. Our results showed that the antagonist with longer receptor residence time (lixivaptan) displayed sustained target occupancy than the antagonist with shorter receptor residence time (mozavaptan). In accordance, lixivaptan displayed insurmountable antagonism and wash-resistant inhibitory effect on the cellular cAMP level, while not so for mozavaptan. Together, our data provide evidence that binding kinetics, next to their affinity, offers additional information for the dynamic process of ligand-receptor interaction. Hopefully, this study may lead to more kinetics-directed medicinal chemistry efforts and aid the design and discovery of different-in-class of vasopressin V₂ receptor ligands for clinical applications.

Multi-scale modeling of drug binding kinetics to predict drug efficacy

Optimizing drug therapies for any disease requires a solid understanding of pharmacokinetics (the drug concentration at a given time point in different body compartments) and pharmacodynamics (the effect a drug has at a given concentration). Mathematical models are frequently used to infer drug concentrations over time based on infrequent sampling and/or in inaccessible body compartments. Models are also used to translate drug action from in vitro to in vivo conditions or from animal models to human patients. Recently, mathematical models that incorporate drug-target binding and subsequent downstream responses have been shown to advance our understanding and increase predictive power of drug efficacy predictions. We here discuss current approaches of modeling drug binding kinetics that aim at improving model-based drug development in the future. This in turn might aid in reducing the large number of failed clinical trials.

The impact of different human IgG capture molecules on the kinetics analysis of antibody-antigen interaction

Surface plasmon resonance (SPR) is a well-established method to characterize biomolecular interactions and is widely used in drug discovery and development. Here, we demonstrate that capture surfaces profoundly impact the binding kinetics parameters that are measured for antibody-antigen interactions. Six unique antibody-antigen interactions were characterized using eight different anti-human IgG capture surfaces. The antigen binding affinities for six different human monoclonal antibodies (hmAbs) captured using three different goat anti-human Fc (AHC) polyclonal antibody (pAb) surfaces were in reasonable agreement (3-7-fold weaker) with those measured by kinetic exclusion assay (KinExA). In contrast, up to 81, 32, 489, 2826, and 219-fold weaker antigen binding affinities were measured using mouse AHC mAb, Protein G, Protein A, Protein A/G, and Protein L surfaces, respectively. Protein A, Protein A/G and Protein G interacted with the Fab of hmAbs, possibly affecting antigen binding to hmAbs captured over these surfaces. Additional studies revealed that mouse AHC mAb binds hmAbs with a weak affinity (5.5-36.3 nM) and t_½ values of 1.4-3.3min, compared to the sub-nanomolar affinities of the goat AHC pAbs. These results emphasize the value of measuring binding kinetics of the capture molecule before immobilizing them onto the sensor surface to perform capture kinetics assays on label-free biosensors.

Simultaneous evaluation of multiple microarray surface chemistries through real-time interferometric imaging

Surface chemistry is a crucial aspect for microarray modality biosensor development. The immobilization capability of the functionalized surface is indeed a limiting factor for the final yield of the binding reaction. In this work, we were able to simultaneously compare the functionality of protein ligands that were locally immobilized on different polymers, while on the same solid support, therefore demonstrating a new way of multiplexing. Our goal was to investigate, in a single experiment, both the immobilization efficiency of a group of reactive polymers and the resulting affinity of the tethered molecules. This idea was demonstrated by spotting many reactive polymers on a Si/SiO₂ chip and depositing the molecular probes on the spots immediately after. As a proof of concept, we focused on which polymers would better immobilize a model protein (α-Lactalbumin) and a peptide (LAC-1). We successfully showed that this protocol is applicable to proteins and peptides with a good efficiency. By means of real-time binding measurements performed with the interferometric reflectance imaging sensor (IRIS), local functionalization proved to be comparable to the classical flat coating solution. The final outcome highlights the multiplexing power of this method: first, it allows to characterize dozens of polymers at once. Secondly, it removes the limitation, related to coated surfaces, that only molecules with the same functional groups can be tethered to the same solid support. By applying this protocol, many types of molecules can be studied simultaneously and immobilization for each probe can be individually optimized.

Measuring Small Molecule Binding to Escherichia coli AcrB by Surface Plasmon Resonance

Antimicrobial resistance (AMR) is rapidly becoming one of the great healthcare challenges. A common mechanism employed by pathogenic bacteria to avoid the action of certain antibiotics is to overexpress efflux pumps that can extrude these drugs from the cell rendering them ineffective. Small molecule inhibitors that target bacterial efflux pumps provide a route toward reversing AMR. Here, we describe the application of surface plasmon resonance (SPR) technology to characterize protein:small molecule interactions between the inner membrane protein AcrB subunit of the Escherichia coli AcrA-AcrB-TolC efflux pump and its substrates and novel inhibitors. The SPR assay provides quantitative data about the kinetics of binding that can help guide the development of new chemotherapies to combat AMR.

Expression, purification and stabilization of human serotonin transporter from E. coli

The serotonin transporter belongs to the family of sodium-chloride coupled neurotransmitter transporter and is related to depression in humans. It is therefore an important drug target to support treatment of depression. Recently, structures of human serotonin transporter in complex with inhibitor molecules have been published. However, the production of large protein amounts for crystallization experiments remains a bottleneck. Here we present the possibility to obtain purified serotonin transporter from E. coli. Fos-choline 12 solubilized target protein was obtained with a purity of >95% and a yield of 1.2 mg L^-1 culture in autoinduction medium. CD spectroscopic analysis of protein stability allowed identifying CHS and POPX as stabilizing components, which increased hSERT thermostability by 7 °C. The kinetic dissociation constant K_D of 2.8 μM (±0.05) for of the inhibitor Desipramine was determined with a ka of 10,848 M - 1 s-1 (±220) and a kd of 0.03 s-1 (±4.7 × 10-5).

Exopolysaccharide from Streptococcus thermophilus as stabilizer in fermented dairy: Binding kinetics and interactions with casein of milk

Exopolysaccharides (EPSs) from lactic acid bacteria have great effect on the quality of fermented milk products. However, the mechanism for the quality improvement has not been well described. This study aimed to investigate the molecular binding kinetics and interactions between EPS obtained from Streptococcus thermophilus AR333 (EPS333) and casein of milk (CM) in a simulated acidifying process. The results indicated that EPS333 had a significant effect on the stability of casein micelles at acidic pH (6.0-4.5) according to the turbidity, ζ-potential, particle size and distribution analysis. The adsorption-desorption study by bio-layer interferometry identified the direct affinity binding between EPS333 and CM, the interactive moiety of casein was α-casein, rather than β- or κ-casein. Fluorescence quenching analysis revealed that the force types of interaction between EPS333 and CM were dynamically changeable during the acidifying process, mainly from electrostatic interaction at pH 7.0-6.5, to hydrophobic or hydrogen bonding at pH 6.5-5.5, and then transferred to electrostatic interaction again at pH 5.5-5.0. Conclusively, EPS333 could bind with CM directly via different binding forces during acidifying process to stabilize the properties of casein micelles.

Ligand-Binding Calculations with Metadynamics

All-atom molecular dynamics simulations can capture the dynamic degrees of freedom that characterize molecular recognition, the knowledge of which constitutes the cornerstone of rational approaches to drug design and optimization. In particular, enhanced sampling algorithms, such as metadynamics, are powerful tools to dramatically reduce the computational cost required for a mechanistic description of the binding process. Here, we describe the essential details characterizing these simulation strategies, focusing on the critical step of identifying suitable reaction coordinates, as well as on the different analysis algorithms to estimate binding affinity and residence times. We conclude with a survey of published applications that provides explicit examples of successful simulations for several targets.

Biased agonism of clinically approved μ-opioid receptor agonists and TRV130 is not controlled by binding and signaling kinetics

Binding and signaling kinetics have previously proven important in validation of biased agonism at GPCRs. Here we provide a comprehensive kinetic pharmacological comparison of clinically relevant μ-opioid receptor agonists, including the novel biased agonist oliceridine (TRV130) which is in clinical trial for pain management. We demonstrate that the bias profile observed for the selected agonists is not time-dependent and that agonists with dramatic differences in their binding kinetic properties can display the same degree of bias. Binding kinetics analyses demonstrate that buprenorphine has 18-fold higher receptor residence time than oliceridine. This is thus the largest pharmacodynamic difference between the clinically approved drug buprenorphine and the clinical candidate oliceridine, since their bias profiles are similar. Further, we provide the first pharmacological characterization of (S)-TRV130 demonstrating that it has a similar pharmacological profile as the (R)-form, oliceridine, but displays 90-fold lower potency than the (R)-form. This difference is driven by a significantly slower association rate. Finally, we show that the selected agonists are differentially affected by G protein-coupled receptor kinase 2 and 5 (GRK2 and GRK5) expression. GRK2 and GRK5 overexpression greatly increased μ-opioid receptor internalization induced by morphine, but only had modest effects on buprenorphine and oliceridine-induced internalization. Overall, our data reveal that the clinically available drug buprenorphine displays a similar pharmacological bias profile in vitro compared to the clinical candidate drug oliceridine and that this bias is independent of binding kinetics suggesting a mechanism driven by receptor-conformations. This article is part of the Special Issue entitled 'New Vistas in Opioid Pharmacology'.

Binding kinetics of ligands acting at GPCRs

The influence of drug-receptor binding kinetics has often been overlooked during the development of new therapeutics that target G protein-coupled receptors (GPCRs). Over the last decade there has been a growing understanding that an in-depth knowledge of binding kinetics at GPCRs is required to successfully target this class of proteins. Ligand binding to a GPCR is often not a simple single step process with ligand freely diffusing in solution. This review will discuss the experiments and equations that are commonly used to measure binding kinetics and how factors such as allosteric regulation, rebinding and ligand interaction with the plasma membrane may influence these measurements. We will then consider the molecular characteristics of a ligand and if these can be linked to association and dissociation rates.

Fibroblast Growth Factor-2 Binding to Heparan Sulfate Proteoglycans Varies with Shear Stress in Flow-Adapted Cells

Fibroblast growth factor 2 (FGF2), an important regulator of angiogenesis, binds to endothelial cell (EC) surface FGF receptors (FGFRs) and heparan sulfate proteoglycans (HSPGs). FGF2 binding kinetics have been predominantly studied in static culture; however, the endothelium is constantly exposed to flow which may affect FGF2 binding. We therefore used experimental and computational techniques to study how EC FGF2 binding changes in flow. ECs adapted to 24 h of flow demonstrated biphasic FGF2-HSPG binding, with FGF2-HSPG complexes increasing up to 20 dynes/cm2 shear stress and then decreasing at higher shear stresses. To understand how adaptive EC surface remodeling in response to shear stress may affect FGF2 binding to FGFR and HSPG, we implemented a computational model to predict the relative effects of flow-induced surface receptor changes. We then fit the computational model to the experimental data using relationships between HSPG availability and FGF2-HSPG dissociation and flow that were developed from a basement membrane study, as well as including HSPG production. These studies suggest that FGF2 binding kinetics are altered in flow-adapted ECs due to changes in cell surface receptor quantity, availability, and binding kinetics, which may affect cell growth factor response.

A live cell NanoBRET binding assay allows the study of ligand-binding kinetics to the adenosine A3 receptor

There is a growing interest in understanding the binding kinetics of compounds that bind to G protein-coupled receptors prior to progressing a lead compound into clinical trials. The widely expressed adenosine A₃ receptor (A₃AR) has been implicated in a range of diseases including immune conditions, and compounds that aim to selectively target this receptor are currently under development for arthritis. Kinetic studies at the A₃AR have been performed using a radiolabelled antagonist, but due to the kinetics of this probe, they have been carried out at 10 °C in membrane preparations. In this study, we have developed a live cell NanoBRET ligand binding assay using fluorescent A₃AR antagonists to measure kinetic parameters of labelled and unlabelled compounds at the A₃AR at physiological temperatures. The kinetic profiles of four fluorescent antagonists were determined in kinetic association assays, and it was found that XAC-ser-tyr-X-BY630 had the longest residence time (RT = 288 ± 62 min) at the A₃AR. The association and dissociation rate constants of three antagonists PSB-11, compound 5, and LUF7565 were also determined using two fluorescent ligands (XAC-ser-tyr-X-BY630 or AV039, RT = 6.8 ± 0.8 min) as the labelled probe and compared to those obtained using a radiolabelled antagonist ([³H]PSB-11, RT = 44.6 ± 3.9 min). There was close agreement in the kinetic parameters measured with AV039 and [³H]PSB-11 but significant differences to those obtained using XAC-S-ser-S-tyr-X-BY630. These data indicate that selecting a probe with the appropriate kinetics is important to accurately determine the kinetics of unlabelled ligands with markedly different kinetic profiles.

Chemical modification of protein A chromatography ligands with polyethylene glycol. I: Effects on IgG adsorption equilibrium, kinetics, and transport

Chemical modification of Protein A (ProA) chromatography ligands with polyethylene glycol (PEGylation) has been proposed as a strategy to increase the process selectivity and resin robustness by providing the ligand with a steric repulsion barrier against non-specific binding. This article comprises a comprehensive study of IgG adsorption and transport in Repligen CaptivA PriMAB resin with PEGylated ProA ligands that are modified using 5.2 and 21.5 kDa PEG chains. We studied the impact of the molecular weight of the PEG as well as the extent of PEGylation for the 5.2 kDa PEG modification. In all cases, PEGylation of ProA ligands decreases the resin average pore size, particle porosity, and static binding capacity for IgG proportional to the volume of conjugated PEG in the resin. Resin batch uptake experiments conducted in bulk via a stirred-tank system and with individual resin particles under confocal laser scanning microscopy suggests that PEGylation introduces heterogeneity into IgG binding kinetics: a fraction of the IgG binding sites are transformed from typical fast association kinetic behavior to slow kinetic behavior. pH gradient elution experiments of an IgG molecule on the modified resins show an increase in IgG elution pH for all modified resins, implying a decrease in IgG-ProA binding affinity on modification. Despite losses in static binding capacity for all resins with PEGylated ligands, the loss of dynamic binding capacity at 10% breakthrough (DBC_10%) ranged more broadly from almost 0-47% depending on the PEG molecular weight and the extent of PEGylation. Minimal losses in DBC_10% were observed with a low extent of PEGylation with a smaller molecular weight PEG, while higher losses were observed at higher extents of PEGylation and with higher molecular weight PEG due to decreased static binding capacity and increased mass transfer resistance. This work provides insight into the practical implications for resin performance if PEGylation is considered as a strategy for selectivity enhancement in affinity chromatography with macromolecular ligands.

A binding kinetics study of human adenosine A3 receptor agonists

The human adenosine A₃ (hA₃) receptor has been suggested as a viable drug target in inflammatory diseases and in cancer. So far, a number of selective hA₃ receptor agonists (e.g. IB-MECA and 2-Cl-IB-MECA) inducing anti-inflammatory or anticancer effects are under clinical investigation. Drug-target binding kinetics is increasingly recognized as another pharmacological parameter, next to affinity, for compound triage in the early phases of drug discovery. However, such a kinetics-driven analysis has not yet been performed for the hA₃ receptor. In this study, we first validated a competition association assay for adenosine A₃ receptor agonists to determine the target interaction kinetics. Affinities and Kinetic Rate Index (KRI) values of 11 ribofurano and 10 methanocarba nucleosides were determined in radioligand binding assays. Afterwards, 15 analogues were further selected (KRI <0.70 or KRI >1.35) for full kinetics characterization. The structure-kinetics relationships (SKR) were derived and longer residence times were associated with methanocarba and enlarged adenine N⁶ and C2 substitutions. In addition, from a k_on-k_off-K_D kinetic map we divided the agonists into three subgroups. A residence time "cliff" was observed, which might be relevant to (N)-methanocarba derivatives' rigid C2-arylalkynyl substitutions. Our findings provide substantial evidence that, next to affinity, additional knowledge of binding kinetics is useful for developing and selecting new hA₃R agonists in the early phase of the drug discovery process.