Prediction of Catch-Slip Bond Transition of Kindlin2/β3 Integrin via Steered Molecular Dynamics Simulation

Kindlin-2 Regulates the Oncogenic Activities of Integrins and TGF-β In Triple Negative Breast Cancer Progression and Metastasis

Background

Kindlin-2, an adaptor protein, is dysregulated in various human cancers, including triple negative breast cancer (TNBC), where it drives tumor progression and metastasis by influencing several cancer hallmarks. One well-established role of Kindlin-2 involves the regulation of integrin signaling, achieved by directly binding to the cytoplasmic tail of the integrin β subunit. In this study, we present novel insights into Kindlin-2's involvement in stabilizing the β1-Integrin:TGF-β type 1 receptor (TβRI) complexes, acting as a physical bridge that links β1-Integrin to TβRI. The loss of Kindlin-2 results in the degradation of this protein complex, leading to the inhibition of downstream oncogenic pathways.

Methods

Our methodology encompassed a diverse range of in vitro assays, including CRISPR/Cas9 gene editing, cell migration, 3D tumorsphere formation and invasion, solid binding, co-immunoprecipitation, cell adhesion and spreading assays, as well as western blot and flow cytometry analyses, utilizing MDA-MB-231 and 4T1 TNBC cell lines. Additionally, preclinical in vivo mouse models of TNBC tumor progression and metastasis were employed to substantiate our findings.

Results

The investigation revealed that the direct interaction between Kindlin-2 and β1-Integrin is mediated through the C-terminal F3 domain of Kindlin-2, while the interaction between Kindlin-2 and TβRI is facilitated through the F2 domain of Kindlin-2. Disruption of this bridge, achieved via CRISPR/Cas9-mediated knockout of Kindlin-2, led to the degradation of β1-Integrin and TβRI, resulting in the inhibition of oncogenic pathways downstream of both proteins, subsequently hindering tumor growth and metastasis. Treatment of Kindlin-2-deficient cells with the proteasome inhibitor MG-132 restored the expression of both β1-Integrin and TβRI. Furthermore, the rescue of Kindlin-2 expression reinstated their oncogenic activities both in vitro and in vivo.

Conclusions

This study identifies a novel function of Kindlin-2 in stabilizing the β1-Integrin:TβR1 complexes and regulating their downstream oncogenic signaling. The translational implications of these findings are substantial, potentially unveiling new therapeutically targeted pathways crucial for the treatment of TNBC tumors.

Molecular mechanisms of catch bonds and their implications for platelet hemostasis

Adhesive molecular bonds between blood cells are essential for thrombosis and hemostasis as they provide means for platelet adhesion, aggregation, and signaling in flowing blood. According to the nowadays conventional definition, a "catch" bond is a type of non-covalent bio-molecular bridge, whose dissociation lifetime counter-intuitively increases with applied tensile force. Following recent experimental findings, such receptor-ligand protein bonds are vital to the blood cells involved in the prevention of bleeding (hemostatic response) and infection (immunity). In this review, we examine the up-to-date experimental discoveries and theoretical insights about catch bonds between the blood cells, their biomechanical principles at the molecular level, and their role in platelet thrombosis and hemostasis.

Multiscale models of integrins and cellular adhesions

Computational models of integrin-based adhesion complexes have revealed important insights into the mechanisms by which cells establish connections with their external environment. However, how changes in conformation and function of individual adhesion proteins regulate the dynamics of whole adhesion complexes remains largely elusive. This is because of the large separation in time and length scales between the dynamics of individual adhesion proteins (nanoseconds and nanometers) and the emergent dynamics of the whole adhesion complex (seconds and micrometers), and the limitations of molecular simulation approaches in extracting accurate free energies, conformational transitions, reaction mechanisms, and kinetic rates, that can inform mechanisms at the larger scales. In this review, we discuss models of integrin-based adhesion complexes and highlight their main findings regarding: (i) the conformational transitions of integrins at the molecular and macromolecular scales and (ii) the molecular clutch mechanism at the mesoscale. Lastly, we present unanswered questions in the field of modeling adhesions and propose new ideas for future exciting modeling opportunities.

A MD Simulation Prediction for Regulation of N-Terminal Modification on Binding of CD47 to CD172a in a Force-Dependent Manner

Cancer cells can evade immune surveillance through binding of its transmembrane receptor CD47 to CD172a on myeloid cells. CD47 is recognized as a promising immune checkpoint for cancer immunotherapy inhibiting macrophage phagocytosis. N-terminal post-translated modification (PTM) via glutaminyl cyclase is a landmark event in CD47 function maturation, but the molecular mechanism underlying the mechano-chemical regulation of the modification on CD47/CD172a remains unclear. Here, we performed so-called "ramp-clamp" steered molecular dynamics (SMD) simulations, and found that the N-terminal PTM enhanced interaction of CD172a with CD47 by inducing a dynamics-driven contraction of the binding pocket of the bound CD172a, an additional constraint on CYS15 on CD47 significantly improved the tensile strength of the complex with or without PTM, and a catch bond phenomenon would occur in complex dissociation under tensile force of 25 pN in a PTM-independent manner too. The residues GLN52 and SER66 on CD172a reinforced the H-bonding with their partners on CD47 in responding to PTM, while ARG69 on CD172 with its partner on CD47 might be crucial in the structural stability of the complex. This work might serve as molecular basis for the PTM-induced function improvement of CD47, should be helpful for deeply understanding CD47-relevant immune response and cancer development, and provides a novel insight in developing of new strategies of immunotherapy targeting this molecule interaction.

Receptor-Ligand Binding: Effect of Mechanical Factors

Gaining insight into the in situ receptor-ligand binding is pivotal for revealing the molecular mechanisms underlying the physiological and pathological processes and will contribute to drug discovery and biomedical application. An important issue involved is how the receptor-ligand binding responds to mechanical stimuli. This review aims to provide an overview of the current understanding of the effect of several representative mechanical factors, such as tension, shear stress, stretch, compression, and substrate stiffness on receptor-ligand binding, wherein the biomedical implications are focused. In addition, we highlight the importance of synergistic development of experimental and computational methods for fully understanding the in situ receptor-ligand binding, and further studies should focus on the coupling effects of these mechanical factors.

Uncovering and engineering the mechanical properties of the adhesion GPCR ADGRG1 GAIN domain

Key cellular functions depend on the transduction of extracellular mechanical signals by specialized membrane receptors including adhesion G-protein coupled receptors (aGPCRs). While recently solved structures support aGPCR activation through shedding of the extracellular GAIN domain, the molecular mechanisms underpinning receptor mechanosensing remain poorly understood. When probed using single-molecule atomic force spectroscopy and molecular simulations, ADGRG1 GAIN dissociated from its tethered agonist at forces significantly higher than other reported signaling mechanoreceptors. Strong mechanical resistance was achieved through specific structural deformations and force propagation pathways under mechanical load. ADGRG1 GAIN variants computationally designed to lock the alpha and beta subdomains and rewire mechanically-induced structural deformations were found to modulate the GPS-Stachel rupture forces. Our study provides unprecedented insights into the molecular underpinnings of GAIN mechanical stability and paves the way for engineering mechanosensors, better understanding aGPCR function, and informing drug-discovery efforts targeting this important receptor class.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Tension Enhances the Binding Affinity of β1 Integrin by Clamping Talin Tightly: An Insight from Steered Molecular Dynamics Simulations

Integrin activation is a predominant step for cell-cell and cell-ECM interactions. Talin and Kindlin are mechanosensitive adaptor proteins that bind to the integrin cytoplasmic tail and mediate integrin activation, cytoskeleton rearrangement, and focal adhesion assembly. However, knowledge about how Talin and Kindlin synergistically assist integrin activation remains unclear. Here, we performed so-called "ramp-clamp" SMD simulations, which modeled the mechanosignaling from Kindlin, to investigate the effect of tension on the interaction of the β1 integrin cytoplasmic tail with the Talin-F3 domain. The present results showed that mild but not excessive stretching enhanced the binding of integrin with Talin. This mechanical regulation on integrin affinity to Talin referred to an event cascade, in which under stretching, the integrin cytoplasmic tail adopted allostery in response to the mechanical stimulus, remodeling of integrin in favor of Talin-association ensued, and finally, a stable, close-knit complex was formed. In the cascade, the torsion angle transition of integrin was the cue for the stable interaction of the complex under tensile force. The present work suggested a model for Talin and Kindlin to synergistically activate integrin. It should help understand integrin activation and its mechanochemical regulation mechanism, integrin-related innate cellular immune responses, cell adhesion, cell-cell interaction, and integrin-related drug development.

Spatiotemporal characteristics of P-selectin-induced β2 integrin activation of human neutrophils under flow

Activation of integrins is crucial for recruitment of flowing leukocytes to inflammatory or injured vascular sites, but their spatiotemporal characteristics are incompletely understood. We discovered that β2-integrin activation over the entire surface of neutrophils on immobilized P-selectin occurred via mitogen-activated protein kinase (MAPK) or non-MAPK signaling with a minute-level timescale in a force-dependent manner. In flow, MAPK signaling required intracellular Ca2+ release to activate integrin within 2 min. Integrin activation via non-MAPK signaling occurred first locally in the vicinity of ligated P-selectin glycoprotein ligand-1 (PSGL-1) within sub-seconds, and then over the entire cell surface within 1 min in an extracellular Ca2+ influx-dependent manner. The transition from a local (but rapid) to global (but slow) activation mode was triggered by ligating the freshly activated integrin. Lipid rafts, moesin, actin, and talin were involved in non-MAPK signaling. Fluid loads had a slight effect on local integrin activation with a second-level timescale, but served as enhancers of global integrin activation.

Force-enhanced biophysical connectivity of platelet β3 integrin signaling through Talin is predicted by steered molecular dynamics simulations

Platelet β3-integrin signaling through Talin is crucial in platelet transmembrane signaling, activation, adhesion, spreading and aggregation, and remains unclear in mechano-microenvironments. In order to examine Talin-β3 integrin biophysical connectivity, a series of “ramp-clamp” steered molecular dynamics (SMD) simulations were performed on complex of F3 domain of Talin and cytoplasmic tail of β3 integrin to imitate different force-loads in platelet. Pull-induced allostery of the hydrophobic pocket in F3 domain might markedly enhance complex rupture-force (> 150pN) and slow down breakage of the complex; the complex should mechano-stable for its conformational conservation under loads (≤ 80pN); increasing force below 60pN would decrease the complex dissociation probability, and force-induced extension of β5 strand on Talin and binding site residues, ASP⁷⁴⁰ and ALA⁷⁴² as well as Asn⁷⁴⁴, on β3-integrin were responsible for the force-enhanced linkage of the Talin-β3 integrin. Force might enhance biophysical connectivity of β3-integrin signaling through Talin by a catch bond mechanism, which be mediated by the force-induced allostery of complex at clamped stage. This work provides a novel insight into the force-regulated transmembrane β3-integrin signaling and its molecular basis for platelet activation, and exhibited a potential power of the present computer strategy in predicting mechanical regulation on ligand-receptor interaction under loads.

Residues R1075, D1090, R1095, and C1130 Are Critical in ADAMTS13 TSP8-Spacer Interaction Predicted by Molecular Dynamics Simulation

ADAMTS13 (A Disintegrin and Metalloprotease with Thrombospondin type 1 repeats, member 13) cleaves von Willebrand Factor (VWF) multimers to limit the prothrombotic function of VWF. The deficiency of ADAMTS13 causes a lethal thrombotic microvascular disease, thrombotic thrombocytopenic purpura (TTP). ADAMTS13 circulates in a “closed” conformation with the distal domain associating the Spacer domain to avoid off-target proteolysis or recognition by auto-antibodies. However, the interactions of the distal TSP8 domain and the Spacer domain remain elusive. Here, we constructed the TSP8-Spacer complex by a combination of homology modelling and flexible docking. Molecular dynamics simulation was applied to map the binding sites on the TSP8 or Spacer domain. The results predicted that R1075, D1090, R1095, and C1130 on the TSP8 domain were key residues that interacted with the Spacer domain. R1075 and R1095 bound exosite-4 tightly, D1090 formed multiple hydrogen bonds and salt bridges with exosite-3, and C1130 interacted with both exosite-3 and exosite-4. Specific mutations of exosite-3 (R568K/F592Y/R660K/Y661F/Y665F) or the four key residues (R1075A/D1090A/R1095A/C1130A) impaired the binding of the TSP8 domain to the Spacer domain. These results shed new light on the understanding of the auto-inhibition of ADAMTS13.

Characterization of the interactions of ADAMTS13 CUB1 domain to WT- and GOF-Spacer domain by molecular dynamics simulation

Metalloprotease ADAMTS13 specifically cleaves VWF (von Willebrand Factor) to prevent excessive platelet aggregation and thrombus formation at the sites of vascular injury. To avoid non-specific cleavage, ADAMTS13 has the auto-inhibition effect in which the Spacer domain in N-terminal interacts with the CUB1 domain in C-terminal, resulting in decreased proteolytic activity. Previous studies reported that exosite-3 in the Spacer domain was a key binding site in the Spacer-CUB1 interaction. When exosite-3 was mutated (R660K/F592Y/R568K/Y661F/Y665F, GOF), the auto-inhibition of ADAMTS13 was disrupted and the enzymatic activity was markedly increased. However, the characteristics of the Spacer-CUB1 interaction is not fully understood. Here, we constructed the model of Spacer-CUB1 complex by homologous modeling and molecular docking to characterize the Spacer-CUB1 binding and predict key amino acid residues via molecular dynamics simulation. Our data showed that G607-S610 was a non-reported potential binding site in the Spacer domain; GOF mutation attenuated the formation of hydrogen bond between exosite-3 and the CUB1 domain; Residues E1231, R1251, L1258, D1259 and T1261 in the CUB1 domain might play an important role in the Spacer-CUB1 interaction. Our study advances the understanding of the structural basis of the auto-inhibition of ADAMTS13 and provides information about the key residues in the binding interface.

Comparison of Linear vs. Cyclic RGD Pentapeptide Interactions with Integrin αvβ3 by Molecular Dynamics Simulations

Simple Summary

The integrin α_vβ₃-RGD motif interaction plays a key role in the progression of malignant tumor. Although two typical cyclic and linear RGD short peptides have been widely used in tumor diagnosis and therapy, little is known about the internal dynamic mechanism for different configurations of RGD peptides with different affinities interacting with the integrin α_vβ₃. Our results showed that the cyclic RGD peptide had a more stable configuration in binding to integrins α_vβ₃, which depended on the higher binding energy and higher static electrical energy, especially in the interaction between Asp^RGD-MIDAS. The steered molecular dynamics simulation showed a stronger interaction for the cyclic RGD-integrin α_vβ₃ system than the linear one, with a larger dissociation force (average peak force) and more time to dissociate. Our findings provide insights into the dynamics of integrin α_vβ₃ interactions with linear and cyclic RGD ligands and offer some new therapeutic approaches for the design and development of novel antitumor drugs.

Abstract

Integrin α_vβ₃ interacting with the short Arg-Gly-Asp (RGD) motif plays a critical role in the progression of several types of tumors. However, the effects of the RGD structure (cyclic or linear) with integrin α_vβ₃ at the atomic level remain poorly understood. Here, we performed association and dissociation dynamic simulations for integrin α_vβ₃ in complex with a linear or cyclic pentapeptide by steered molecular dynamics simulations. Compared with cyclic RGD, the linear RGD peptide triggers instability of the configurational changes, mainly resting with the RGD domain due to its flexibility. The main interaction energy between Mg²⁺ and cyclic RGD is much stronger than that of the linear RGD system by the well shield to lessen attacks by free water molecules. The force-dependent dissociation results show that it is easier for linear RGD peptides to leave the active site and much quicker than the cyclic RGD ligand, whereas it is harder to enter the appropriate active binding site in linear RGD. The Ser¹²³-Asp^RGD bond may play a critical role in the allosteric pathway. Our findings provide insights into the dynamics of α_vβ₃ interactions with linear and cyclic RGD ligands and contribute to the application of RGD-based strategies in preclinical therapy.

MD Simulations on a Well-Built Docking Model Reveal Fine Mechanical Stability and Force-Dependent Dissociation of Mac-1/GPIbα Complex

Interaction of leukocyte integrin macrophage-1 antigen (Mac-1) to platelet glycoprotein Ibα (GPIbα) is critical for platelet-leukocyte crosstalk in hemostasis and inflammatory responses to vessel injuries under hemodynamic environments. The mechano-regulation and its molecular basis for binding of Mac-1 to GPIbα remain unclear, mainly coming from the lack of crystal structure of the Mac-1/GPIbα complex. We herein built a Mac-1/GPIbα complex model through a novel computer strategy, which included a flexible molecular docking and system equilibrium followed by a "force-ramp + snapback" molecular dynamics (MD) simulation. With this model, a series of "ramp-clamp" steered molecular dynamics (SMD) simulations were performed to examine the GPIbα-Mac-1 interaction under various loads. The results demonstrated that the complex was mechano-stable for both the high rupture force (>250 pN) at a pulling velocity of 3 Å/ns and the conformational conservation under various constant tensile forces (≤75 pN); a catch-slip bond transition was predicted through the dissociation probability, examined with single molecular AFM measurements, reflected by the interaction energy and the interface H-bond number, and related to the force-induced allostery of the complex; besides the mutation-identified residues D222 and R218, the residues were also dominant in the binding of Mac-1 to GPIbα. This study recommended a valid computer strategy for building a likely wild-type docking model of a complex, provided a novel insight into the mechanical regulation mechanism and its molecular basis for the interaction of Mac-1 with GPIbα, and would be helpful for understanding the platelet-leukocyte interaction in hemostasis and inflammatory responses under mechano-microenvironments.