The two-pathway model for the catch-slip transition in biological adhesion

Some recently studied biological noncovalent bonds have shown increased lifetime when stretched by mechanical force. In each case these counterintuitive "catch-bonds" have transitioned into ordinary "slip-bonds" that become increasingly shorter lived as the tensile force on the bond is further increased. We describe analytically how these results are supported by a physical model whereby the ligand escapes the receptor binding site via two alternative routes, a catch-pathway that is opposed by the applied force and a slip-pathway that is promoted by force. The model predicts under what conditions and at what critical force the catch-to-slip transition would be observed, as well as the degree to which the bond lifetime is enhanced at the critical force. The model is applied to four experimentally studied systems taken from the literature, involving the binding of P- and L-selectins to sialyl Lewis(X) oligosaccharide-containing ligands. Good quantitative fit to the experimental data is obtained, both for experiments with a constant force and for experiments where the force increases linearly with time.     

Altered Dynamics of S. aureus Phosphofructokinase via Bond Restraints at Two Distinct Allosteric Binding Sites

The effect of perturbation at the allosteric site was investigated through several replicas of molecular dynamics (MD) simulations conducted on bacterial phosphofructokinase (SaPFK). In our previous work, an alternative binding site was estimated to be allosteric in addition to the experimentally reported one. To highlight the effect of both allosteric sites on receptor’s dynamics, MD runs were carried out on apo forms with and without perturbation. Perturbation was achieved via incorporating multiple bond restraints for residue pairs located at the allosteric site. Restraints applied to the predicted site caused one dimer to stiffen, whereas an increase in mobility was detected in the same dimer when the experimentally resolved site was restrained. Fluctuations in C_α-C_α distances which is used to disclose residues with high potential of communication indicated a marked increase in signal transmission within each dimer as the receptor switched to a restrained state. Cross-correlation of positional fluctuations indicated an overall decrease in the magnitude of both positive and negative correlations when restraints were employed on the predicted allosteric site whereas an exact opposite effect was observed for the reported site. Finally, mutual correspondence between positional fluctuations noticeably increased with restraints on predicted allosteric site, whereas an opposite effect was observed for restraints applied on experimentally reported one. In view of these findings, it is clear that the perturbation of either one of two allosteric sites effected the dynamics of the receptor with a distinct and contrasting character.     

Model energy landscapes and the force-induced dissociation of ligand-receptor bonds

We discuss models for the force-induced dissociation of a ligand-receptor bond, occurring in the context of cell adhesion or single molecule unbinding force measurements. We consider a bond with a structured energy landscape which is modeled by a network of force dependent transition rates between intermediate states. The behavior of a model with only one intermediate state and a model describing a molecular zipper is studied. We calculate the bond lifetime as a function of an applied force and unbinding forces under an increasing applied load and determine the relationship between both quantities. The dissociation via an intermediate state can lead to distinct functional relations of the bond lifetime on force. One possibility is the occurrence of three force regimes where the lifetime of the bond is determined by different transitions within the energy landscape. This case can be related to recent experimental observations of the force-induced dissociation of single avidin-biotin bonds.     

Benzodiazepine ligands can act as allosteric modulators of the Type 1 cholecystokinin receptor

The cholecystokinin (CCK₁) receptor is a G protein-coupled receptor important for nutrient homeostasis. The molecular basis of CCK-receptor binding has been debated, with one prominent model suggesting occupation of the same region of the intramembranous helical bundle as benzodiazepines. Here, we used a specific assay of allosteric ligand interaction to probe the mode of binding of devazepide, a prototypic benzodiazepine ligand. Devazepide elicited marked slowing of dissociation of pre-bound CCK, only possible through binding to a topographically distinct allosteric site. This effect was disrupted by chemical modification of a cysteine in the benzodiazepine-binding pocket. Application of an allosteric model to the equilibrium interaction between a series of benzodiazepine ligands and CCK yielded quantitative estimates of each modulator’s affinity for the allosteric site, as well as the degree of negative cooperativity for the interaction between occupied orthosteric and allosteric sites. The allosteric nature of benzodiazepine binding to the CCK₁ receptor provides new opportunities for small molecule drug development.     

Mathematical model for the effects of adhesion and mechanics on cell migration speed.

Migration of mammalian blood and tissue cells over adhesive surfaces is apparently mediated by specific reversible reactions between cell membrane adhesion receptors and complementary ligands attached to the substratum. Although in a number of systems these receptors and ligand molecules have been isolated and identified, a theory capable of predicting the effects of their properties on cell migration behavior currently does not exist. We present a simple mathematical model for elucidating the dependence of cell speed on adhesion-receptor/ligand binding and cell mechanical properties. Our model can be applied to propose answers to questions such as: does an optimal adhesiveness exist for cell movement? How might changes in receptor and ligand density and/or affinity affect the rate of migration? Can cell rheological properties influence movement speed? This model incorporates cytoskeletal force generation, cell polarization, and dynamic adhesion as requirements for persistent cell movement. A critical feature is the proposed existence of an asymmetry in some cell adhesion-receptor property, correlated with cell polarity. We consider two major alternative mechanisms underlying this asymmetry: (a) a spatial distribution of adhesion-receptor number due to polarized endocytic trafficking and (b) a spatial variation in adhesion-receptor/ligand bond strength. Applying a viscoelastic-solid model for cell mechanics allows us to represent one-dimensional locomotion with a system of differential equations describing cell deformation and displacement along with adhesion-receptor dynamics. In this paper, we solve these equations under the simplifying assumption that receptor dynamics are at a quasi-steady state relative to cell locomotion. Thus, our results are strictly valid for sufficiently slow cell movement, as typically observed for tissue cells such as fibroblasts. Numerical examples relevant to experimental systems are provided. Our results predict how cell speed might vary with intracellular contractile force, cell rheology, receptor/ligand kinetics, and receptor/ligand number densities. A biphasic dependence is shown to be possible with respect to some of the system parameters, with position of the maxima essentially governed by a balance between transmitted contractile force and adhesiveness. We demonstrate that predictions for the two alternative asymmetry mechanisms can be distinguished and could be experimentally tested using cell populations possessing different adhesion-receptor numbers.     

Triphasic Force Dependence of E-Selectin/Ligand Dissociation Governs Cell Rolling under Flow

During inflammation, flowing leukocytes tether to and roll on vascular surfaces through the association and dissociation of selectin/ligand bonds. The interactions of P- and L- selectins with their respective ligands exhibit catch-slip bonds, such that increasing force initially prolongs and then shortens bond lifetimes. In addition, catch-slip bonds have been shown to govern L-selectin-mediated cell rolling. Using a flow chamber and biomembrane force probe, we show a triphasic force dependence of E-selectin/ligand dissociation that initially behaves as slip bonds, then transitions to catch bonds, and finally transitions again to slip bonds as the force increases. These transitions govern the velocities of neutrophils, HL-60 cells, and Colo-205 cells rolling on E-selectin, as evidenced by the fact that their velocities exhibited a triphasic force dependence that inversely matched the triphasic lifetime-force relationship. At low forces, slip bonds may also precede catch bonds for interactions of P- and L-selectin with their ligands.     

Functional properties of human skeletal muscle acetylcholine receptors expressed by the TE671 cell line

Functional properties of acetylcholine receptors from intact TE671 human medulloblastoma cells were examined using tracer ion flux, ligand competition against 125I-labeled alpha-bungarotoxin binding, and single channel recording measurements. 125I-Labeled alpha-bungarotoxin binds to surface receptors with the forward rate constant 1.8 X 10(5) M-1 s-1 and dissociates with the rate constant 4.6 X 10(-5) s-1, at 21 degrees C; the apparent dissociation constant is 2.6 X 10(-10) M. alpha-Bungarotoxin binds to at least two sites/receptor, but blocks agonist-induced 22Na+ uptake when bound to only one site. The reversible antagonists, dimethyl-d-tubocurarine and gallamine, occupy two sites which exhibit nearly equivalent affinities, but block agonist-induced uptake by occupying only one site. Strong agonists activate rapid sodium uptake with relatively low affinity, but desensitize with a much higher affinity; among agonists, the ratio of low to high affinity dissociation constants ranges from 1600 to 4000. By using the estimated dissociation constants, the allosteric model of Monod, Wyman, and Changeux (MWC) can be fitted to the concentration dependencies of both steady-state agonist occupancy and desensitization. The fitting analysis discloses an allosteric constant of 3 X 10(-5), which is the ratio of activatable to desensitized receptors in the absence of agonist. The rate of recovery from desensitization can exceed the rate of onset of desensitization elicited by low concentrations of agonist, further supporting the general MWC framework. Single channel recordings show that the channel opening probability is greater than 0.7 at high agonist concentrations. Favorable channel opening is shown to only slightly oppose strong desensitization.     

基于弹性网络模型的蛋白质变构路径与关键残基识别研究

目的 变构效应在蛋白质生物学功能执行过程中发挥着重要的调控作用，如何基于蛋白质空间结构，有效识别变构信号的传播路径和关键的残基位点是蛋白质结构-功能关系研究领域的热点科学问题。方法 本研究利用基于弹性网络模型（elastic network model，ENM）的力分布计算方法，通过分析蛋白质对外力的响应过程，来识别体系的变构路径以及变构过程中的关键残基。在该方法中，对蛋白质的关键变构位点施加外力，通过对体系形变以及内力分布情况的分析，有效识别与外力承载区域形变相耦合的关键残基，从而得到力信号在蛋白质结构内的传播路径。结果 利用该方法研究了人类磷酸甘油酸激酶（human phosphoglycerate kinase，hPGK）和蛋白质酪氨酸磷酸酶（protein tyrosine phosphatase，PTP）PDZ2结构域的变构调控路径和关键残基。对于hPGK，识别出从底物结合位点到铰链区的两条变构信号传导路径。对于PTP PDZ2，也成功识别出从配体结合位点传递到蛋白质远端的两条长程变构调控路径。计算结果与实验和分子动力学（molecular dynamics，MD）模拟得到的结果一致。结论 本研究为蛋白质体系关键残基识别及变构路径研究提供了有效的分析方法。     

A stereochemical model of the transpeptidation complex

Molecular models are proposed to describe the relative arrangement of aminoacyl and peptidyl tRNAs when bound to their respective A and P sites on the ribosome. The crystallographically determined structures of tRNAasp and tRNAphe have served as the models for these bound structures, while the imposed steric constraints for the model complexes were based on the results of published experimental data. The constructed models satisfy the stereochemical requirements needed for codon-anticodon interaction and for peptide bond formation. In this paper, the results of the complex containing tRNAphe as the A and P site bound transfer RNAs, is compared to a similarly constructed model which uses tRNAasp as the ribosome-bound transfer RNAs. The models have the following three major features: 1) the aminoacyl and peptidyl transfer RNAs assume an angle of approximately 45 degrees relative to each other; 2) in providing the proper stereochemistry for peptide bond condensation, a significant kink must be present in the messenger RNA between the A site and P site codons; and 3) a comparison of the two model complexes indicates that structural variations between the tRNAs or any allosteric transitions of the transfer RNAs associated with codon-anticodon recognition may be accommodated in the model by way of freedom of rotation about the phosphate backbone bonds in the mRNA between consecutive codons.     

A hydrogen bond in loop A is critical for the binding and function of the 5-HT3 receptor

The binding sites of Cys-loop receptors are formed from at least six loops (A-F). Here we have used mutagenesis, radioligand binding, voltage clamp electrophysiology, and homology modeling to probe the role of two residues in loop A of the 5-HT3 receptor: Asn128 and Glu129. The data show that substitution of Asn128, with a range of alternative natural and unnatural amino acids, changed the EC50 (from approximately 10-fold more potent to approximately 10-fold less potent than that of the wild type), increased the maximal peak current for mCPBG compared to 5-HT (R max) 2-19-fold, and decreased n H, indicating this residue is involved in receptor gating; we propose Asn128 faces away from the binding pocket and plays a role in facilitating transitions between conformational states. Substitutions of Glu129 resulted in functional receptors only when the residue could accept a hydrogen bond, but with both these and other substitutions, no [(3)H]granisetron binding could be detected, indicating a role in ligand binding. We propose that Glu129 faces into the binding pocket, where, through its ability to hydrogen bond, it plays a critical role in ligand binding. Thus, the data support a modified model of the 5-HT3 receptor binding site and show that loop A plays a critical role in both the ligand binding and function of this receptor.     

How the headpiece hinge angle is opened: New insights into the dynamics of integrin activation

How the integrin head transitions to the high-affinity conformation is debated. Although experiments link activation with the opening of the hinge angle between the betaA and hybrid domains in the ligand-binding headpiece, this hinge is closed in the liganded alpha(v)beta3 integrin crystal structure. We replaced the RGD peptide ligand of this structure with the 10th type III fibronectin module (FnIII10) and discovered through molecular dynamics (MD) equilibrations that when the conformational constraints of the leg domains are lifted, the betaA/hybrid hinge opens spontaneously. Together with additional equilibrations on the same nanosecond timescale in which small structural variations impeded hinge-angle opening, these simulations allowed us to identify the allosteric pathway along which ligand-induced strain propagates via elastic distortions of the alpha1 helix to the betaA/hybrid domain hinge. Finally, we show with steered MD how force accelerates hinge-angle opening along the same allosteric pathway. Together with available experimental data, these predictions provide a novel framework for understanding integrin activation.     

Molecular mechanisms of allosteric probe dependence in μ opioid receptor

Allostery is one of the most important features of proteins. It greatly contributes to the complexity of life, since it enables possibility of precise tuning of protein function, as well as performing more than one function per protein. Probe dependence is one of the unique features of allostery. It allows a protein to respond differently to the same allosteric modulator when different drugs or transmitters are bound. Unfortunately, allosteric mechanisms are difficult to investigate experimentally. Instead, they can be reproduced artificially in simulations. We simulated in silico a native-like cell membrane fragment with an active-state human μ opioid receptor (MOR) in order to investigate diverse effects of a receptor’s positive allosteric modulator on various agonists. Particular emphasis on native-likeness of the environment was put. We managed to reproduce the experimentally observed effects, which allowed us to take deeper insight into their underlying mechanisms. We found an allosteric pathway in the receptor, leading from the ligand binding site to the intracellular, effector site. We observed that the modulator affected the pathway, inducing different resultant responses for full and partial agonists.     

Energy landscape of streptavidin-biotin complexes measured by atomic force microscopy

The dissociation of ligand and receptor involves multiple transitions between intermediate states formed during the unbinding process. In this paper, we explored the energy landscape of the streptavidin-biotin interaction by using the atomic force microscope (AFM) to measure the unbinding dynamics of individual ligand-receptor complexes. The rupture force of the streptavidin-biotin bond increased more than 2-fold over a range of loading rates between 100 and 5000 pN/s. Moreover, the force measurements showed two regimes of loading in the streptavidin-biotin force spectrum, revealing the presence of two activation barriers in the unbinding process. Parallel experiments carried out with a streptavidin mutant (W120F) were used to investigate the molecular determinants of the activation barriers. From these experiments, we attributed the outer activation barrier in the energy landscape to the molecular interaction of the '3-4' loop of streptavidin that closes behind biotin.     

Force-induced unfolding simulations of the human Notch1 negative regulatory region: possible roles of the heterodimerization domain in mechanosensing

Notch receptors are core components of the Notch signaling pathway and play a central role in cell fate decisions during development as well as tissue homeostasis. Upon ligand binding, Notch is sequentially cleaved at the S2 site by an ADAM protease and at the S3 site by the γ-secretase complex. Recent X-ray structures of the negative regulatory region (NRR) of the Notch receptor reveal an auto-inhibited fold where three protective Lin12/Notch repeats (LNR) of the NRR shield the S2 cleavage site housed in the heterodimerization (HD) domain. One of the models explaining how ligand binding drives the NRR conformation from a protease-resistant state to a protease-sensitive one invokes a mechanical force exerted on the NRR upon ligand endocytosis. Here, we combined physics-based atomistic simulations and topology-based coarse-grained modeling to investigate the intrinsic and force-induced folding and unfolding mechanisms of the human Notch1 NRR. The simulations support that external force applied to the termini of the NRR disengages the LNR modules from the heterodimerization (HD) domain in a well-defined, largely sequential manner. Importantly, the mechanical force can further drive local unfolding of the HD domain in a functionally relevant fashion that would provide full proteolytic access to the S2 site prior to heterodimer disassociation. We further analyzed local structural features, intrinsic folding free energy surfaces, and correlated motions of the HD domain. The results are consistent with a model in which the HD domain possesses inherent mechanosensing characteristics that could be utilized during Notch activation. This potential role of the HD domain in ligand-dependent Notch activation may have implications for understanding normal and aberrant Notch signaling.     

Regulation of the catalytic function of coagulation factor VIIa by a conformational linkage of surface residue Glu 154 to the active site

Coagulation factor VIIa is an allosterically regulated trypsin-like serine protease that initiates the coagulation pathways upon complex formation with its cellular receptor and cofactor tissue factor (TF). The analysis of a conformation-sensitive monoclonal antibody directed to the macromolecular substrate exosite in the VIIa protease domain demonstrated a conformational link from this exosite to the catalytic cleft that is independent of cofactor-induced allosteric changes. In this study, we identify Glu 154 as a critical surface-exposed exosite residue side chain that undergoes conformational changes upon active site inhibitor binding. The Glu 154 side chain is important for hydrolysis of scissile bond mimicking peptidyl p-nitroanilide substrates, and for inhibition of VIIa's amidolytic function upon antibody binding. This exosite residue is not linked to the catalytic cleft residue Lys 192 which plays an important role in thrombin's allosteric coupling to exosite I. Allosteric linkages between VIIa's active site and the cofactor binding site or between the cofactor binding site and the macromolecular substrate exosite were not influenced by mutation of Glu 154. Glu 154 thus only influences the linkage of the macromolecular substrate binding exosite to the catalytic center. These data provide novel evidence that allosteric regulation of VIIa's catalytic function involves discrete and independent conformational linkages and that allosteric transitions in the VIIa protease domain are not globally coupled.     

Allosteric effect of gallamine on muscarinic cholinergic receptor binding

Gallamine interacts with an allosteric site on muscarinic acetylcholine receptor complexes in rat brain membranes, thereby slowing the dissociation of a radiolabelled ligand ([³H]N-methylscopolamine) from the receptor complex. This effect involves the elimination of the fast component of the biphasic dissociation curve. The allosteric effect of gallamine is equally prominent in membranes containing predominantly Ml (cerebral cortex) and M2 (brainstem) subtypes of muscarinic receptor. Gallamine's action is not affected by a variety of treatments which influence the conformational state of the receptor as reflected by agonist binding affinity, including treatments with heat,N-ethymaleimide and trypsin. A guanine nucleotide (5-guanylylimidodiphosphate), however, moderates the effects of gallamine on muscarinic receptors in brainstem, but not in cortical, membranes.     

Linked-function origins of cooperativity in a symmetrical dimer

The thermodynamic origins of substrate binding cooperativity in a dimeric enzyme that can bind one substrate (A) and one allosteric ligand (X) to each of two identical subunits are discussed. It is assumed that maximal activity is not subject to allosteric modification and that the substrates and allosteric ligands achieve binding equilibrium in the steady state. Each uniquely ligated form is assumed to be capable of exhibiting unique binding properties, and only the principles of thermodynamic linkage are used to constrain the system further. The explicit relationship between the Hill coefficient, the concentration of X, and the magnitudes of the relevant coupling free energies and dissociation constants is derived. In the absence of X only the homotropic coupling between substrate sites contributes to a nonhyperbolic substrate saturation profile. An allosteric ligand, X, can alter the cooperativity in two distinct ways, one mechanism being manifested when X is saturating and the only only when X is present at saturating concentrations. By evaluating the concentration of substrate required to produce half-maximal velocity as a function of [X], as well as the Hill coefficients when X is absent and fully saturating, the dissociation and coupling constants most important for understanding the mechanisms of allosteric action in an enzyme of this type can be determined.     

Study of the fragmentation patterns of the phosphate-arginine noncovalent bond

Our previous work has highlighted the role of certain amino acid residues, mainly two or more adjacent arginine on one peptide and two or more adjacent glutamate, or aspartate, or a phosphorylated residue on the other in the formation of noncovalent complexes (NCX) between peptides. In the present study, we employ ESI-MS to investigate the gas-phase stability and dissociation pathways of the NCX of a basic peptide VLRRRRKRVN, an epitope from the third intracellular loop of the dopamine D(2) receptor, with the phosphopetide SVSTDpTpSAE, an epitope from the cannabinoid CB1 carboxyl terminus. ESI-MS/MS analysis of the NCX between VLRRRRKRVN and SVSTDpTpSAE suggests two dissociation pathways for the NCX. The major pathway is the disruption of the electrostatic interactions between the Arg residues and the phosphate groups, while an alternative pathway is also recorded, in which the complex is dissociated along the covalent bond between the oxygen from either Thr or Ser and HPO(3). To verify the alternative pathway, we have used an ion trap instrument to conduct MS(3) analysis on the product ions of both dissociation pathways.