Mechano-Redox Control of Mac-1 De-Adhesion by PDI Promotes Directional Movement Under Flow.

BACKGROUND: Neutrophil migration is critical to the initiation and resolution of inflammation. Macrophage-1 antigen (Mac-1; CD11b/CD18, αMß2) is a leukocyte integrin essential for firm adhesion to endothelial ICAM-1 (intercellular adhesion molecule 1) and migration of neutrophils in the shear forces of the circulation. PDI (protein disulfide isomerase) has been reported to influence neutrophil adhesion and migration. We aimed to elucidate the molecular mechanism of PDI control of Mac-1 affinity for ICAM-1 during neutrophil migration under fluid shear. METHODS: Neutrophils isolated from whole blood were perfused over microfluidic chips coated with ICAM-1. Colocalization of Mac-1 and PDI on neutrophils was visualized by fluorescently labeled antibodies and confocal microscopy. The redox state of Mac-1 disulfide bonds was mapped by differential cysteine alkylation and mass spectrometry. Wild-type or disulfide mutant Mac-1 was expressed recombinantly in Baby Hamster Kidney cells to measure ligand affinity. Mac-1 conformations were measured by conformation-specific antibodies and molecular dynamics simulations. Neutrophils crawling on immobilized ICAM-1 were measured in presence of oxidized or reduced PDI, and the effect of PDI inhibition using isoquercetin on neutrophil crawling on inflamed endothelial cells was examined. Migration indices in the X- and Y-direction were determined and the crawling speed was calculated. RESULTS: PDI colocalized with high-affinity Mac-1 at the trailing edge of stimulated neutrophils when crawling on ICAM-1 under fluid shear. PDI cleaved 2 allosteric disulfide bonds, C169-C176 and C224-C264, in the ßI domain of the ß2 subunit, and cleavage of the C224-C264 disulfide bond selectively controls Mac-1 disengagement from ICAM-1 under fluid shear. Molecular dynamics simulations and conformation-specific antibodies reveal that cleavage of the C224-C264 bond induces conformational change and mechanical stress in the ßI domain. This allosterically alters the exposure of an αI domain epitope associated with a shift of Mac-1 to a lower-affinity state. These molecular events promote neutrophil motility in the direction of flow at high shear stress. Inhibition of PDI by isoquercetin reduces neutrophil migration in the direction of flow on endothelial cells during inflammation. CONCLUSIONS: Shear-dependent PDI cleavage of the neutrophil Mac-1 C224-C264 disulfide bond triggers Mac-1 de-adherence from ICAM-1 at the trailing edge of the cell and enables directional movement of neutrophils during inflammation.

O-glycans Expand Lubricin and Attenuate Its Viscosity and Shear Thinning.

Lubricin, an intrinsically disordered glycoprotein, plays a pivotal role in facilitating smooth movement and ensuring the enduring functionality of synovial joints. The central domain of this protein serves as a source of this excellent lubrication and is characterized by its highly glycosylated, negatively charged, and disordered structure. However, the influence of O-glycans on the viscosity of lubricin remains unclear. In this study, we employ molecular dynamics simulations in the absence and presence of shear, along with continuum simulations, to elucidate the intricate interplay between O-glycans and lubricin and the impact of O-glycans on lubricin's conformational properties and viscosity. We found the presence of O-glycans to induce a more extended conformation in fragments of the disordered region of lubricin. These O-glycans contribute to a reduction in solution viscosity but at the same time weaken shear thinning at high shear rates, compared to nonglycosylated systems with the same density. This effect is attributed to the steric and electrostatic repulsion between the fragments, which prevents their conglomeration and structuring. Our computational study yields a mechanistic mechanism underlying previous experimental observations of lubricin and paves the way to a more rational understanding of its function in the synovial fluid.

Mechanical forces control the valency of the malaria adhesin VAR2CSA by exposing cryptic glycan binding sites.

Plasmodium falciparum (Pf) is responsible for the most lethal form of malaria. VAR2CSA is an adhesin protein expressed by this parasite at the membrane of infected erythrocytes for attachment to the placenta, leading to pregnancy-associated malaria. VAR2CSA is a large 355 kDa multidomain protein composed of nine extracellular domains, a transmembrane helix, and an intracellular domain. VAR2CSA binds to Chondroitin Sulphate A (CSA) of the proteoglycan matrix of the placenta. Shear flow, as the one occurring in blood, has been shown to enhance the (VAR2CSA-mediated) adhesion of Pf-infected erythrocytes on the CSA-matrix. However, the underlying molecular mechanism governing this enhancement has remained elusive. Here, we address this question by using equilibrium, force-probe, and docking-based molecular dynamics simulations. We subjected the VAR2CSA protein-CSA sugar complex to a force mimicking the tensile force exerted on this system due to the shear of the flowing blood. We show that upon this force exertion, VAR2CSA undergoes a large opening conformational transition before the CSA sugar chain dissociates from its main binding site. This preferential order of events is caused by the orientation of the molecule during elongation, as well as the strong electrostatic attraction of the sugar to the main protein binding site. Upon opening, two additional cryptic CSA binding sites get exposed and a functional dodecameric CSA molecule can be stably accommodated at these force-exposed positions. Thus, our results suggest that mechanical forces increase the avidity of VAR2CSA by turning it from a monovalent to a multivalent state. We propose this to be the molecular cause of the observed shear-enhanced adherence. Mechanical control of the valency of VAR2CSA is an intriguing hypothesis that can be tested experimentally and which is of relevance for the understanding of the malaria infection and for the development of anti placental-malaria vaccines targeting VAR2CSA.

FERM domains recruit ample PI(4,5)P2s to form extensive protein-membrane attachments.

The four-point-one ezrin-radixin-moesin homology (FERM) protein domain is a multifunctional protein-lipid binding site, constituting an integral part of numerous membrane-associated proteins. Its interaction with the lipid phosphatidylinositol-4,5-bisphosphate (PIP2), located at the inner leaflet of eukaryotic plasma membranes, is important for localization, anchorage, and activation of FERM-containing proteins. FERM-PIP2 complexes structurally determined so far exclusively feature a 1:1 binding stoichiometry of protein and lipid, with a few basic FERM residues neutralizing the -4 charge of the bound PIP2. Whether this picture from static crystal structures also applies to the dynamic interaction of FERM domains on PIP2 membranes is unknown. We here quantified the stoichiometry of FERM-PIP2 binding in a lipid bilayer using atomistic molecular dynamics simulations and experiments on solid supported membranes for the FERM domains of focal adhesion kinase and ezrin. In contrast to the structural data, we find much higher average stoichiometries of FERM-PIP2 binding, amounting to 1:3 or 1:4 ratios, respectively. In simulations, the full set of basic residues at the membrane interface, 7 and 15 residues for focal adhesion kinase and ezrin, respectively, engages in PIP2 interactions. In addition, Na ions enter the FERM-membrane binding interface, compensating negative PIP2 charges in case of high charge surpluses from bound PIP2. We propose the multivalent binding of FERM domains to PIP2 in lipid bilayers to significantly enhance the stability of FERM-membrane binding and to render the FERM-membrane linkage highly adjustable.

The interplay between adsorption and aggregation of von Willebrand factor chains in shear flows.

Von Willebrand factor (VWF) is a giant extracellular glycoprotein that carries out a key adhesive function during primary hemostasis. Upon vascular injury and triggered by the shear of flowing blood, VWF establishes specific interactions with several molecular partners in order to anchor platelets to collagen on the exposed subendothelial surface. VWF also interacts with itself to form aggregates that, adsorbed on the surface, provide more anchor sites for the platelets. However, the interplay between elongation and subsequent exposure of cryptic binding sites, self-association, and adsorption on the surface remained unclear for VWF. In particular, the role of shear flow in these three processes is not well understood. In this study, we address these questions by using Brownian dynamics simulations at a coarse-grained level of resolution. We considered a system consisting of multiple VWF-like self-interacting chains that also interact with a surface under a shear flow. By a systematic analysis, we reveal that chain-chain and chain-surface interactions coexist nontrivially to modulate the spontaneous adsorption of VWF and the posterior immobilization of secondary tethered chains. Accordingly, these interactions tune VWF's extension and its propensity to form shear-assisted functional adsorbed aggregates. Our data highlight the collective behavior VWF self-interacting chains have when bound to the surface, distinct from that of isolated or flowing chains. Furthermore, we show that the extension and the exposure to solvent have a similar dependence on shear flow, at a VWF-monomer level of resolution. Overall, our results highlight the complex interplay that exists between adsorption, cohesion, and shear forces and their relevance for the adhesive hemostatic function of VWF.

The role of flow in the self-assembly of dragline spider silk proteins.

Hydrodynamic flow in the spider duct induces conformational changes in dragline spider silk proteins (spidroins) and drives their assembly, but the underlying physical mechanisms are still elusive. Here we address this challenging multiscale problem with a complementary strategy of atomistic and coarse-grained molecular dynamics simulations with uniform flow. The conformational changes at the molecular level were analyzed for single-tethered spider silk peptides. Uniform flow leads to coiled-to-stretch transitions and pushes alanine residues into ß sheet and poly-proline II conformations. Coarse-grained simulations of the assembly process of multiple semi-flexible block copolymers using multi-particle collision dynamics reveal that the spidroins aggregate faster but into low-order assemblies when they are less extended. At medium-to-large peptide extensions (50%-80%), assembly slows down and becomes reversible with frequent association and dissociation events, whereas spidroin alignment increases and alanine repeats form ordered regions. Our work highlights the role of flow in guiding silk self-assembly into tough fibers by enhancing alignment and kinetic reversibility, a mechanism likely relevant also for other proteins whose function depends on hydrodynamic flow.

Structure and dynamics of the von Willebrand Factor C6 domain.

Von Willebrand disease (VWD) is a bleeding disorder with different levels of severity. VWD-associated mutations are located in the von Willebrand factor (VWF) gene, coding for the large multidomain plasma protein VWF with essential roles in hemostasis and thrombosis. On the one hand, a variety of mutations in the C-domains of VWF are associated with increased bleeding upon vascular injury. On the other hand, VWF gain-of-function (GOF) mutations in the C4 domain have recently been identified, which induce an increased risk of myocardial infarction. Mechanistic insights into how these mutations affect the molecular behavior of VWF are scarce and holistic approaches are challenging due to the multidomain and multimeric character of this large protein. Here, we determine the structure and dynamics of the C6 domain and the single nucleotide polymorphism (SNP) variant G2705R in C6 by combining nuclear magnetic resonance spectroscopy, molecular dynamics simulations and aggregometry. Our findings indicate that this mutation mostly destabilizes VWF by leading to a more pronounced hinging between both subdomains of C6. Hemostatic parameters of variant G2705R are close to normal under static conditions, but the missense mutation results in a gain-of-function under flow conditions, due to decreased VWF stem stability. Together with the fact that two C4 variants also exhibit GOF characteristics, our data underline the importance of the VWF stem region in VWF's hemostatic activity and the risk of mutation-associated prothrombotic properties in VWF C-domain variants due to altered stem dynamics.

Thermodynamic stabilization of von Willebrand factor A1 domain induces protein loss of function.

The von Willebrand disease (vWD) is the most common hereditary bleeding disorder caused by defects of the von Willebrand Factor (vWF), a large extracellular protein in charge of adhering platelets to sites of vascular lesions. vWF performs this essential homeostatic task via specific protein-protein interactions between the vWF A1 domain and the platelet receptor, the glycoprotein Ib alpha (GPIBα). The two naturally occurring vWF A1 domain mutations G1324A and G1324S, near the GPIBα binding site, induce a dramatic decrease in platelet adhesion, resulting in a bleeding disorder classified as type 2M vWD. However, the reason for the drastic phenotypic response induced by these two supposedly minor modifications remains unclear. We addressed this question using a combination of equilibrium-molecular dynamics (MD) and nonequilibrium MD-based free energy simulations. Our data confirms that both mutations maintain the highly stable Rossmann fold of the vWF A1 domain. G1324A and G1324S mutations hardly changed the per-residue flexibility of the A1 domain but induced a global conformational change affecting the region near the binding site to GPIBα. Furthermore, we observed two significant changes in the vWF A1 domain upon mutation, the global redistribution of the internal mechanical stress and the increased thermodynamic stability of the A1 domain. These observations are consistent with previously reported mutations increasing the melting temperature. Overall, our results support the idea of thermodynamic conformational restriction of A1-before the binding to GPIBα-as a crucial factor determining the loss-of-function of the G1324A(S) vWD mutants.

Martini 3 coarse-grained force field for poly(para-phenylene ethynylene)s.

Poly(para-phenylene ethynylene)s, or short PPEs, are a class of conjugated and semi-flexible polymers with a strongly delocalized π electron system and increased chain stiffness. Due to this, PPEs have a wide range of technological applications. Although the material properties of single-chains or mixtures of few PPE chains have been studied in detail, the properties of large assemblies remain to be fully explored. Here, we developed a coarse-grained model for PPEs with the Martini 3 force field to enable computational studies of PPEs in large-scale assembly. We used an optimization geometrical approach to take the shape of the π conjugated backbone into account and also applied an additional angular potential to tune the mechanical bending stiffness of the polymer. Our Martini 3 model reproduces key structural and thermodynamic observables of single PPE chains and mixtures, such as persistence length, density, packing and stacking. We show that chain entanglement increases with the expense of nematic ordering with growing PPE chain length. With the Martini 3 PPE model at hand, we are now able to cover large spatio-temporal scales and thereby to uncover key aspects for the structural organization of PPE bulk systems. The model is also predicted to be of high applicability to investigate out-of-equilibrium behavior of PPEs under mechanical force.

DNA binds to a specific site of the adhesive blood-protein von Willebrand factor guided by electrostatic interactions.

Neutrophils release their intracellular content, DNA included, into the bloodstream to form neutrophil extracellular traps (NETs) that confine and kill circulating pathogens. The mechanosensitive adhesive blood protein, von Willebrand Factor (vWF), interacts with the extracellular DNA of NETs to potentially immobilize them during inflammatory and coagulatory conditions. Here, we elucidate the previously unknown molecular mechanism governing the DNA-vWF interaction by integrating atomistic, coarse-grained, and Brownian dynamics simulations, with thermophoresis, gel electrophoresis, fluorescence correlation spectroscopy (FCS), and microfluidic experiments. We demonstrate that, independently of its nucleotide sequence, double-stranded DNA binds to a specific helix of the vWF A1 domain, via three arginines. This interaction is attenuated by increasing the ionic strength. Our FCS and microfluidic measurements also highlight the key role shear-stress has in enabling this interaction. Our simulations attribute the previously-observed platelet-recruitment reduction and heparin-size modulation, upon establishment of DNA-vWF interactions, to indirect steric hindrance and partial overlap of the binding sites, respectively. Overall, we suggest electrostatics-guiding DNA to a specific protein binding site-as the main driving force defining DNA-vWF recognition. The molecular picture of a key shear-mediated DNA-protein interaction is provided here and it constitutes the basis for understanding NETs-mediated immune and hemostatic responses.

Structure and dynamics of the platelet integrin-binding C4 domain of von Willebrand factor.

Von Willebrand factor (VWF) is a key player in the regulation of hemostasis by promoting recruitment of platelets to sites of vascular injury. An array of 6 C domains forms the dimeric C-terminal VWF stem. Upon shear force activation, the stem adopts an open conformation allowing the adhesion of VWF to platelets and the vessel wall. To understand the underlying molecular mechanism and associated functional perturbations in disease-related variants, knowledge of high-resolution structures and dynamics of C domains is of paramount interest. Here, we present the solution structure of the VWF C4 domain, which binds to the platelet integrin and is therefore crucial for the VWF function. In the structure, we observed 5 intra- and inter-subdomain disulfide bridges, of which 1 is unique in the C4 domain. The structure further revealed an unusually hinged 2-subdomain arrangement. The hinge is confined to a very short segment around V2547 connecting the 2 subdomains. Together with 2 nearby inter-subdomain disulfide bridges, this hinge induces slow conformational changes and positional alternations of both subdomains with respect to each other. Furthermore, the structure demonstrates that a clinical gain-of-function VWF variant (Y2561) is more likely to have an effect on the arrangement of the C4 domain with neighboring domains rather than impairing platelet integrin binding.

Molecular Dynamics Simulations of Molecules in Uniform Flow.

Flow at the molecular level induces shear-induced unfolding of single proteins and can drive their assembly, the mechanisms of which are not completely understood. To be able to analyze the role of flow on molecules, we present uniform-flow molecular dynamics simulations at atomic level. The pull module of the GRoningen MAchine for Chemical Simulations package was extended to be able to force-group atoms within a defined layer of the simulation box. Application of this external enforcement to explicit water molecules, together with the coupling to a thermostat, led to a uniform terminal velocity of the solvent water molecules. We monitored the density of the whole system to establish the conditions under which the simulated flow is well-behaved. A maximal velocity of 1.3 m/s can be generated if a pull slice of 8 nm is used, and high velocities would require larger pull slices to still maintain a stable density. As expected, the target velocity increases linearly with the total external force applied. Finally, we suggest an appropriate setup to stretch a protein by uniform flow, in which protein extensions depend on the flow conditions. Our implementation provides an efficient computational tool to investigate the effect of the flow at the molecular level.

GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations.

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations. Most importantly, it allows the comparison between computed and reference maps (e.g., experimental) through calculation of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quantitatively contrast perturbed and control simulation data sets and to examine how much biomolecular systems resemble both synthetic and experimental density maps. This was especially advantageous for multimolecule systems in which standard comparisons like RMSDs are difficult to compute. In addition, GROmaρs incorporates absolute and relative spatial free-energy estimates to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the density calculation. Furthermore, masking of regions was implemented to speed up calculations and to facilitate the comparison with experimental maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and average charge distribution in biomolecular systems. We employed all these functionalities to inspect the localization of lipid and water molecules in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, we anticipate a high applicability of GROmaρs for the analysis of molecular dynamics simulations and their comparison with experimentally determined densities.

Understanding Conformational Dynamics of Complex Lipid Mixtures Relevant to Biology.

This is a perspective article entitled "Frontiers in computational biophysics: understanding conformational dynamics of complex lipid mixtures relevant to biology" which is following a CECAM meeting with the same name.

von Willebrand factor is dimerized by protein disulfide isomerase.

Multimeric von Willebrand factor (VWF) is essential for primary hemostasis. The biosynthesis of VWF high-molecular-weight multimers requires spatial separation of each step because of varying pH value requirements. VWF is dimerized in the endoplasmic reticulum by formation of disulfide bonds between the C-terminal cysteine knot (CK) domains of 2 monomers. Here, we investigated the basic question of which protein catalyzes the dimerization. We examined the putative interaction of VWF and the protein disulfide isomerase PDIA1, which has previously been used to visualize endoplasmic reticulum localization of VWF. Excitingly, we were able to visualize the PDI-VWF dimer complex by high-resolution stochastic optical reconstruction microscopy and atomic force microscopy. We proved and quantified direct binding of PDIA1 to VWF, using microscale thermophoresis and fluorescence correlation spectroscopy (dissociation constants KD = 236 ± 66 nM and KD = 282 ± 123 nM by microscale thermophoresis and fluorescence correlation spectroscopy, respectively). The similar KD (258 ± 104 nM) measured for PDI interaction with the isolated CK domain and the atomic force microscopy images strongly indicate that PDIA1 binds exclusively to the CK domain, suggesting a key role of PDIA1 in VWF dimerization. On the basis of protein-protein docking and molecular dynamics simulations, combined with fluorescence microscopy studies of VWF CK-domain mutants, we suggest the following mechanism of VWF dimerization: PDI initiates VWF dimerization by forming the first 2 disulfide bonds Cys2771-2773' and Cys2771'-2773. Subsequently, the third bond, Cys2811-2811', is formed, presumably to protect the first 2 bonds from reduction, thereby rendering dimerization irreversible. This study deepens our understanding of the mechanism of VWF dimerization and the pathophysiological consequences of its inhibition.

Mutation G1629E Increases von Willebrand Factor Cleavage via a Cooperative Destabilization Mechanism.

The large multimeric glycoprotein von Willebrand Factor (VWF) plays a pivotal adhesive role during primary hemostasis. VWF is cleaved by the protease ADAMTS13 as a down-regulatory mechanism to prevent excessive VWF-mediated platelet aggregation. For each VWF monomer, the ADAMTS13 cleavage site is located deeply buried inside the VWF A2 domain. External forces in vivo or denaturants in vitro trigger the unfolding of this domain, thereby leaving the cleavage site solvent-exposed and ready for cleavage. Mutations in the VWF A2 domain, facilitating the cleavage process, cause a distinct form of von Willebrand disease (VWD), VWD type 2A. In particular, the VWD type 2A Gly1629Glu mutation drastically accelerates the proteolytic cleavage activity, even in the absence of forces or denaturants. However, the effect of this mutation has not yet been quantified, in terms of kinetics or thermodynamics, nor has the underlying molecular mechanism been revealed. In this study, we addressed these questions by using fluorescence correlation spectroscopy, molecular dynamics simulations, and free energy calculations. The measured enzyme kinetics revealed a 20-fold increase in the cleavage rate for the Gly1629Glu mutant compared with the wild-type VWF. Cleavage was found cooperative with a cooperativity coefficient n = 2.3, suggesting that the mutant VWF gives access to multiple cleavage sites of the VWF multimer at the same time. According to our simulations and free energy calculations, the Gly1629Glu mutation causes structural perturbation in the A2 domain and thereby destabilizes the domain by ∼10 kJ/mol, promoting its unfolding. Taken together, the enhanced proteolytic activity of Gly1629Glu can be readily explained by an increased availability of the ADAMTS13 cleavage site through A2-domain-fold thermodynamic destabilization. Our study puts forward the Gly1629Glu mutant as a very efficient enzyme substrate for ADAMTS13 activity assays.

Mutual A domain interactions in the force sensing protein von Willebrand factor.

The von Willebrand factor (VWF) is a glycoprotein in the blood that plays a central role in hemostasis. Among other functions, VWF is responsible for platelet adhesion at sites of injury via its A1 domain. Its adjacent VWF domain A2 exposes a cleavage site under shear to degrade long VWF fibers in order to prevent thrombosis. Recently, it has been shown that VWF A1/A2 interactions inhibit the binding of platelets to VWF domain A1 in a force-dependent manner prior to A2 cleavage. However, whether and how this interaction also takes place in longer VWF fragments as well as the strength of this interaction in the light of typical elongation forces imposed by the shear flow of blood remained elusive. Here, we addressed these questions by using single molecule force spectroscopy (SMFS), Brownian dynamics (BD), and molecular dynamics (MD) simulations. Our SMFS measurements demonstrate that the A2 domain has the ability to bind not only to single A1 domains but also to VWF A1A2 fragments. SMFS experiments of a mutant [A2] domain, containing a disulfide bond which stabilizes the domain against unfolding, enhanced A1 binding. This observation suggests that the mutant adopts a more stable conformation for binding to A1. We found intermolecular A1/A2 interactions to be preferred over intramolecular A1/A2 interactions. Our data are also consistent with the existence of two cooperatively acting binding sites for A2 in the A1 domain. Our SMFS measurements revealed a slip-bond behavior for the A1/A2 interaction and their lifetimes were estimated for forces acting on VWF multimers at physiological shear rates using BD simulations. Complementary fitting of AFM rupture forces in the MD simulation range adequately reproduced the force response of the A1/A2 complex spanning a wide range of loading rates. In conclusion, we here characterized the auto-inhibitory mechanism of the intramolecular A1/A2 bond as a shear dependent safeguard of VWF, which prevents the interaction of VWF with platelets.

Stress Propagation through Biological Lipid Bilayers in Silico.

Membrane tension plays various critical roles in the cell. We here asked how fast and how far localized pulses of mechanical stress dynamically propagate through biological lipid bilayers. In both coarse-grained and all-atom molecular dynamics simulations of a dipalmitoylphosphatidylcholine lipid bilayer, we observed nanometer-wide stress pulses, propagating very efficiently longitudinally at a velocity of approximately 1.4 ± 0.5 nm/ps (km/s), in close agreement with the expected speed of sound from experiments. Remarkably, the predicted characteristic attenuation time of the pulses was in the order of tens of picoseconds, implying longitudinal stress propagation over length scales up to several tens of nanometers before damping. Furthermore, the computed dispersion relation leading to such damping was consistent with proposed continuum viscoelastic models of propagation. We suggest this mode of stress propagation as a potential ultrafast mechanism of signaling that may quickly couple mechanosensitive elements in crowded biological membranes.

Mechanism of Focal Adhesion Kinase Mechanosensing.

Mechanosensing at focal adhesions regulates vital cellular processes. Here, we present results from molecular dynamics (MD) and mechano-biochemical network simulations that suggest a direct role of Focal Adhesion Kinase (FAK) as a mechano-sensor. Tensile forces, propagating from the membrane through the PIP2 binding site of the FERM domain and from the cytoskeleton-anchored FAT domain, activate FAK by unlocking its central phosphorylation site (Tyr576/577) from the autoinhibitory FERM domain. Varying loading rates, pulling directions, and membrane PIP2 concentrations corroborate the specific opening of the FERM-kinase domain interface, due to its remarkably lower mechanical stability compared to the individual alpha-helical domains and the PIP2-FERM link. Analyzing downstream signaling networks provides further evidence for an intrinsic mechano-signaling role of FAK in broadcasting force signals through Ras to the nucleus. This distinguishes FAK from hitherto identified focal adhesion mechano-responsive molecules, allowing a new interpretation of cell stretching experiments.