Contribution of unfolding and intermolecular architecture to fibronectin fiber extensibility

The extracellular matrix contains components with remarkable mechanical properties, including fibronectin (Fn) fibers with extensibilities of >700% strain. We utilized what we consider a novel technique to quantify the extent of molecular unfolding that contributes to Fn fiber extension, and we compared this behavior with stochastic models of Fn fibers with different molecular arrangements. In vitro unfolding as a function of strain was measured by fluorescently labeling cysteines in modules FnIII7 and III15 in artificial Fn fibers. A calibration technique we also consider novel made it possible to demonstrate that 44% of cysteines in these modules were exposed in Fn fibers strained to 421% extension, up from 8% exposure without strain. In silico unfolding was measured by applying a constant strain rate to a fiber represented by a network of wormlike chain springs, each representing an individual Fn molecule. Unfolding rates were calculated with a tension-dependent stochastic model applied to FnIII modules in each molecule. A comparison of these approaches revealed that only a molecular arrangement permitting unequal mechanical loading of Fn molecules recapitulates in vitro unfolding. These data have implications for Fn-dependent mechanotransduction and give insight into how the molecular architecture of natural materials permits such remarkable extensibility.     

Cells in 3D matrices under interstitial flow: Effects of extracellular matrix alignment on cell shear stress and drag forces

Interstitial flow is an important regulator of various cell behaviors both in vitro and in vivo, yet the forces that fluid flow imposes on cells embedded in a 3D extracellular matrix (ECM), and the effects of matrix architecture on those forces, are not well understood. Here, we demonstrate how fiber alignment can affect the shear and pressure forces on the cell and ECM. Using computational fluid dynamics simulations, we show that while the solutions of the Brinkman equation accurately estimate the average fluid shear stress and the drag forces on a cell within a 3D fibrous medium, the distribution of shear stress on the cellular surface as well as the peak shear stresses remain intimately related to the pericellular fiber architecture and cannot be estimated using bulk-averaged properties. We demonstrate that perpendicular fiber alignment of the ECM yields lower shear stress and pressure forces on the cells and higher stresses on the ECM, leading to decreased permeability, while parallel fiber alignment leads to higher stresses on cells and increased permeability, as compared to a cubic lattice arrangement. The Spielman–Goren permeability relationships for fibrous media agreed well with CFD simulations of flow with explicitly considered fibers. These results suggest that the experimentally observed active remodeling of ECM fibers by fibroblasts under interstitial flow to a perpendicular alignment could serve to decrease the shear and drag forces on the cell.     

Integrin α5/fibronectin1 and focal adhesion kinase are required for lens fiber morphogenesis in zebrafish

Lens fiber formation and morphogenesis requires a precise orchestration of cell– extracellular matrix (ECM) and cell–cell adhesive changes in order for a lens epithelial cell to adopt a lens fiber fate, morphology, and migratory ability. The cell–ECM interactions that mediate these processes are largely unknown, and here we demonstrate that fibronectin1 (Fn1), an ECM component, and integrin α5, its cellular binding partner, are required in the zebrafish lens for fiber morphogenesis. Mutations compromising either of these proteins lead to cataracts, characterized by defects in fiber adhesion, elongation, and packing. Loss of integrin α5/Fn1 does not affect the fate or viability of lens epithelial cells, nor does it affect the expression of differentiation markers expressed in lens fibers, although nucleus degradation is compromised. Analysis of the intracellular mediators of integrin α5/Fn1 activity focal adhesion kinase (FAK) and integrin-linked kinase (ILK) reveals that FAK, but not ILK, is also required for lens fiber morphogenesis. These results support a model in which lens fiber cells use integrin α5 to migrate along a Fn-containing substrate on the apical side of the lens epithelium and on the posterior lens capsule, likely activating an intracellular signaling cascade mediated by FAK in order to orchestrate the cytoskeletal changes in lens fibers that facilitate elongation, migration, and compaction.     

Evidence that αC region is origin of low modulus, high extensibility, and strain stiffening in fibrin fibers

Fibrin fibers form the structural scaffold of blood clots and perform the mechanical task of stemming blood flow. Several decades of investigation of fibrin fiber networks using macroscopic techniques have revealed remarkable mechanical properties. More recently, the microscopic origins of fibrin's mechanics have been probed through direct measurements on single fibrin fibers and individual fibrinogen molecules. Using a nanomanipulation system, we investigated the mechanical properties of individual fibrin fibers. The fibers were stretched with the atomic force microscope, and stress-versus-strain data was collected for fibers formed with and without ligation by the activated transglutaminase factor XIII (FXIIIa). We observed that ligation with FXIIIa nearly doubled the stiffness of the fibers. The stress-versus-strain behavior indicates that fibrin fibers exhibit properties similar to other elastomeric biopolymers. We propose a mechanical model that fits our observed force extension data, is consistent with the results of the ligation data, and suggests that the large observed extensibility in fibrin fibers is mediated by the natively unfolded regions of the molecule. Although some models attribute fibrin's force-versus-extension behavior to unfolding of structured regions within the monomer, our analysis argues that these models are inconsistent with the measured extensibility and elastic modulus.     

Walking the Line: A Fibronectin Fiber-Guided Assay to Probe Early Steps of (Lymph)angiogenesis

Angiogenesis and lymphangiogenesis are highly complex morphogenetic processes, central to many physiological and pathological conditions, including development, cancer metastasis, inflammation and wound healing. While it is described that extracellular matrix (ECM) fibers are involved in the spatiotemporal regulation of angiogenesis, current angiogenesis assays are not specifically designed to dissect and quantify the underlying molecular mechanisms of how the fibrillar nature of ECM regulates vessel sprouting. Even less is known about the role of the fibrillar ECM during the early stages of lymphangiogenesis. To address such questions, we introduced here an in vitro (lymph)angiogenesis assay, where we used microbeads coated with endothelial cells as simple sprouting sources and deposited them on single Fn fibers used as substrates to mimic fibrillar ECM. The fibers were deposited on a transparent substrate, suitable for live microscopic observation of the ensuing cell outgrowth events at the single cell level. Our proof-of-concept studies revealed that fibrillar Fn, compared to Fn-coated surfaces, provides far stronger sprouting and guidance cues to endothelial cells, independent of the tested mechanical strains of the Fn fibers. Additionally, we found that VEGF-A, but not VEGF-C, stimulates the collective outgrowth of lymphatic endothelial cells (LEC), while the collective outgrowth of blood vascular endothelial cells (HUVEC) was prominent even in the absence of these angiogenic factors. In addition to the findings presented here, the modularity of our assay allows for the use of different ECM or synthetic fibers as substrates, as well as of other cell types, thus expanding the range of applications in vascular biology and beyond.     

Fibrin-fiber architecture influences cell spreading and differentiation

ABSTRACT

The mechanical and structural properties of the extracellular matrix (ECM) play an important role in regulating cell fate. The natural ECM has a complex fibrillar structure and shows nonlinear mechanical properties, which are both difficult to mimic synthetically. Therefore, systematically testing the influence of ECM properties on cellular behavior is very challenging. In this work we show two different approaches to tune the fibrillar structure and mechanical properties of fibrin hydrogels. Addition of extra thrombin before gelation increases the protein density within the fibrin fibers without significantly altering the mechanical properties of the resulting hydrogel. On the other hand, by forming a composite hydrogel with a synthetic biomimetic polyisocyanide network the protein density within the fibrin fibers decreases, and the mechanics of the composite material can be tuned by the PIC/fibrin mass ratio. The effect of the changes in gel structure and mechanics on cellular behavior are investigated, by studying human mesenchymal stem cell (hMSC) spreading and differentiation on these gels. We find that the trends observed in cell spreading and differentiation cannot be explained by the bulk mechanics of the gels, but correlate to the density of the fibrin fibers the gels are composed of. These findings strongly suggest that the microscopic properties of individual fibers in fibrous networks play an essential role in determining cell behavior.     

Heparin-dependent regulation of fibronectin matrix conformation

Extracellular matrix (ECM) conformation is regulated by a variety of stimuli in vivo, including mechanical forces and allosteric binding partners, and these conformational changes contribute to the regulation of cell behavior. Heparin and heparan sulfate, for example, have been shown to regulate the sequestration and presentation of numerous growth factors, including vascular endothelial growth factor, on the heparin 2 binding domain in fibronectin (Fn). However, mechanical force also alters Fn conformation, indicating that the growth factor binding region may be co-regulated by both heparin and mechanical force. Herein, we describe a simple antibody-based method for evaluating the conformation of the heparin 2 binding domain in Fn, and use it to determine the relative contributions of heparin and mechanical strain to the regulation of Fn conformation. We achieved specificity in quantifying conformational changes in this region of Fn by measuring the ratio of two fluorescent monoclonal antibodies, one that is insensitive to Fn conformational changes and a second whose binding is reduced or enhanced by non-equilibrium conformational changes. Importantly, this technique is shown to work on Fn adsorbed on surfaces, single Fn fibers, and Fn matrix fibers in cell culture. Using our dual antibody approach, we show that heparin and mechanical strain co-regulate Fn conformation in matrix fibrils, which is the first demonstration of heparin-dependent regulation of Fn in its physiologically-relevant fibrillar state. Furthermore, the dual antibody approach utilizes commercially available antibodies and simple immunohistochemistry, thus making it accessible to a wide range of scientists interested in Fn mechanobiology.     

Silk self-assembly mechanisms and control from thermodynamics to kinetics

Silkworms and spiders generate fibers that exhibit high strength and extensibility. The underlying mechanisms involved in processing silk proteins into fiber form remain incompletely understood, resulting in the failure to fully recapitulate the remarkable properties of native fibers in vitro from regenerated silk solutions. In the present study, the extensibility and high strength of regenerated silks were achieved by mimicking the natural spinning process. Conformational transitions inside micelles, followed by aggregation of micelles and their stabilization as they relate to the metastable structure of silk are described. Subsequently, the mechanisms to control the formation of nanofibrous structures were elucidated. The results clarify that the self-assembly of silk in aqueous solution is a thermodynamically driven process where kinetics also play a key role. Four key factors, molecular mobility, charge, hydrophilic interactions, and concentration underlie the process. Adjusting these factors can balance nanostructure and conformational composition, and be used to achieve silk-based materials with properties comparable to native fibers. These mechanisms suggest new directions to design silk-based multifunctional materials.     

Single-molecule force spectroscopy reveals a stepwise unfolding of Caenorhabditis elegans giant protein kinase domains

Myofibril assembly and disassembly are complex processes that regulate overall muscle mass. Titin kinase has been implicated as an initiating catalyst in signaling pathways that ultimately result in myofibril growth. In titin, the kinase domain is in an ideal position to sense mechanical strain that occurs during muscle activity. The enzyme is negatively regulated by intramolecular interactions occurring between the kinase catalytic core and autoinhibitory/regulatory region. Molecular dynamics simulations suggest that human titin kinase acts as a force sensor. However, the precise mechanism(s) resulting in the conformational changes that relieve the kinase of this autoinhibition are unknown. Here we measured the mechanical properties of the kinase domain and flanking Ig/Fn domains of the Caenorhabditis elegans titin-like proteins twitchin and TTN-1 using single-molecule atomic force microscopy. Our results show that these kinase domains have significant mechanical resistance, unfolding at forces similar to those for Ig/Fn β-sandwich domains (30-150 pN). Further, our atomic force microscopy data is consistent with molecular dynamic simulations, which show that these kinases unfold in a stepwise fashion, first an unwinding of the autoinhibitory region, followed by a two-step unfolding of the catalytic core. These data support the hypothesis that titin kinase may function as an effective force sensor.     

ECM Protein Nanofibers and Nanostructures Engineered Using Surface-initiated Assembly

The extracellular matrix (ECM) in tissues is synthesized and assembled by cells to form a 3D fibrillar, protein network with tightly regulated fiber diameter, composition and organization. In addition to providing structural support, the physical and chemical properties of the ECM play an important role in multiple cellular processes including adhesion, differentiation, and apoptosis. In vivo, the ECM is assembled by exposing cryptic self-assembly (fibrillogenesis) sites within proteins. This process varies for different proteins, but fibronectin (FN) fibrillogenesis is well-characterized and serves as a model system for cell-mediated ECM assembly. Specifically, cells use integrin receptors on the cell membrane to bind FN dimers and actomyosin-generated contractile forces to unfold and expose binding sites for assembly into insoluble fibers. This receptor-mediated process enables cells to assemble and organize the ECM from the cellular to tissue scales. Here, we present a method termed surface-initiated assembly (SIA), which recapitulates cell-mediated matrix assembly using protein-surface interactions to unfold ECM proteins and assemble them into insoluble fibers. First, ECM proteins are adsorbed onto a hydrophobic polydimethylsiloxane (PDMS) surface where they partially denature (unfold) and expose cryptic binding domains. The unfolded proteins are then transferred in well-defined micro- and nanopatterns through microcontact printing onto a thermally responsive poly(N-isopropylacrylamide) (PIPAAm) surface. Thermally-triggered dissolution of the PIPAAm leads to final assembly and release of insoluble ECM protein nanofibers and nanostructures with well-defined geometries. Complex architectures are possible by engineering defined patterns on the PDMS stamps used for microcontact printing. In addition to FN, the SIA process can be used with laminin, fibrinogen and collagens type I and IV to create multi-component ECM nanostructures. Thus, SIA can be used to engineer ECM protein-based materials with precise control over the protein composition, fiber geometry and scaffold architecture in order to recapitulate the structure and composition of the ECM in vivo.     

A fiber-based constitutive model predicts changes in amount and organization of matrix proteins with development and disease in the mouse aorta

Decreased elastin in mice (Eln+/?) yields a functioning vascular system with elevated blood pressure and increased arterial stiffness that is morphologically distinct from wild-type mice (WT). Yet, function is retained enough that there is no appreciable effect on life span and some mechanical properties are maintained constant. It is not understood how the mouse modifies the normal developmental process to produce a functioning vascular system despite a deficiency in elastin. To quantify changes in mechanical properties, we have applied a fiber-based constitutive model to mechanical data from the ascending aorta during postnatal development of WT and Eln+/? mice. Results indicate that the fiber-based constitutive model is capable of distinguishing elastin amounts and identifying trends during development. We observe an increase in predicted circumferential stress contribution from elastin with age, which correlates with increased elastin amounts from protein quantification data. The model also predicts changes in the unloaded collagen fiber orientation with age, which must be verified in future work. In Eln+/? mice, elastin amounts are decreased at each age, along with the predicted circumferential stress contribution of elastin. Collagen amounts in Eln+/? aorta are comparable to WT, but the predicted circumferential stress contribution of collagen is increased. This may be due to altered organization or structure of the collagen fibers. Relating quantifiable changes in arterial mechanics with changes in extracellular matrix (ECM) protein amounts will help in understanding developmental remodeling and in producing treatments for human diseases affecting ECM proteins.     

A Three-Dimensional Computational Model of Collagen Network Mechanics

Extracellular matrix (ECM) strongly influences cellular behaviors, including cell proliferation, adhesion, and particularly migration. In cancer, the rigidity of the stromal collagen environment is thought to control tumor aggressiveness, and collagen alignment has been linked to tumor cell invasion. While the mechanical properties of collagen at both the single fiber scale and the bulk gel scale are quite well studied, how the fiber network responds to local stress or deformation, both structurally and mechanically, is poorly understood. This intermediate scale knowledge is important to understanding cell-ECM interactions and is the focus of this study. We have developed a three-dimensional elastic collagen fiber network model (bead-and-spring model) and studied fiber network behaviors for various biophysical conditions: collagen density, crosslinker strength, crosslinker density, and fiber orientation (random vs. prealigned). We found the best-fit crosslinker parameter values using shear simulation tests in a small strain region. Using this calibrated collagen model, we simulated both shear and tensile tests in a large linear strain region for different network geometry conditions. The results suggest that network geometry is a key determinant of the mechanical properties of the fiber network. We further demonstrated how the fiber network structure and mechanics evolves with a local formation, mimicking the effect of pulling by a pseudopod during cell migration. Our computational fiber network model is a step toward a full biomechanical model of cellular behaviors in various ECM conditions.     

Elastic modulus and hydraulic permeability of MDCK monolayers

The critical role of cell mechanics in tissue health has led to the development of many in vitro methods that measure the elasticity of the cytoskeleton and whole cells, yet the connection between these local cell properties and bulk measurements of tissue mechanics remains unclear. To help bridge this gap, we have developed a monolayer indentation technique for measuring multi-cellular mechanics in vitro. Here, we measure the elasticity of cell monolayers and uncover the role of fluid permeability in these multi-cellular systems, finding that the resistance of fluid transport through cells controls their force–response at long times.     

Elucidation of extracellular matrix mechanics from muscle fibers and fiber bundles

The importance of the extracellular matrix (ECM) in muscle is widely recognized, since ECM plays a central role in proper muscle development (Buck and Horwitz, 1987), tissue structural support (Purslow, 2002), and transmission of mechanical signals between fibers and tendon (Huijing, 1999). Since substrate biomechanical properties have been shown to be critical in the biology of tissue development and remodeling (Engler et al., 2006; Gilbert et al., 2010), it is likely that mechanics are critical for ECM to perform its function. Unfortunately, there are almost no data available regarding skeletal muscle ECM viscoelastic properties. This is primarily due to the impossibility of isolating and testing muscle ECM. Therefore, this note presents a new method to quantify viscoelastic ECM modulus by combining tests of single muscle fibers and fiber bundles. Our results demonstrate that ECM is a highly nonlinearly elastic material, while muscle fibers are linearly elastic.     

Thermal Adaptation of Conformational Dynamics in Ribonuclease H

The relationship between inherent internal conformational processes and enzymatic activity or thermodynamic stability of proteins has proven difficult to characterize. The study of homologous proteins with differing thermostabilities offers an especially useful approach for understanding the functional aspects of conformational dynamics. In particular, ribonuclease HI (RNase H), an 18 kD globular protein that hydrolyzes the RNA strand of RNA:DNA hybrid substrates, has been extensively studied by NMR spectroscopy to characterize the differences in dynamics between homologs from the mesophilic organism E. coli and the thermophilic organism T. thermophilus. Herein, molecular dynamics simulations are reported for five homologous RNase H proteins of varying thermostabilities and enzymatic activities from organisms of markedly different preferred growth temperatures. For the E. coli and T. thermophilus proteins, strong agreement is obtained between simulated and experimental values for NMR order parameters and for dynamically averaged chemical shifts, suggesting that these simulations can be a productive platform for predicting the effects of individual amino acid residues on dynamic behavior. Analyses of the simulations reveal that a single residue differentiates between two different and otherwise conserved dynamic processes in a region of the protein known to form part of the substrate-binding interface. Additional key residues within these two categories are identified through the temperature-dependence of these conformational processes.     

The Unfolding MD Simulations of Cyclophilin: Analyzed by Surface Contact Networks and Their Associated Metrics

Currently, considerable interest exists with regard to the dissociation of close packed aminoacids within proteins, in the course of unfolding, which could result in either wet or dry moltenglobules. The progressive disjuncture of residues constituting the hydrophobic core ofcyclophilin from L. donovani (LdCyp) has been studied during the thermal unfolding of the molecule, by molecular dynamics simulations. LdCyp has been represented as a surface contactnetwork (SCN) based on the surface complementarity (S_m) of interacting residues within themolecular interior. The application of S_m to side chain packing within proteins make it a very sensitive indicator of subtle perturbations in packing, in the thermal unfolding of the protein. Network based metrics have been defined to track the sequential changes in the disintegration ofthe SCN spanning the hydrophobic core of LdCyp and these metrics prove to be highly sensitive compared to traditional metrics in indicating the increased conformational (and dynamical) flexibility in the network. These metrics have been applied to suggest criteria distinguishing DMG, WMG and transition state ensembles and to identify key residues involved in crucial conformational/topological events during the unfolding process.     

Rate-Dependent Behavior of the Amorphous Phase of Spider Dragline Silk

The time-dependent stress-strain behavior of spider dragline silk was already observed decades ago, and has been attributed to the disordered sequences in silk proteins, which compose the soft amorphous matrix. However, the actual molecular origin and magnitude of internal friction within the amorphous matrix has remained inaccessible, because experimentally decomposing the mechanical response of the amorphous matrix from the embedded crystalline units is challenging. Here, we used atomistic molecular dynamics simulations to obtain friction forces for the relative sliding of peptide chains of Araneus diadematus spider silk within bundles of these chains as a representative unit of the amorphous matrix in silk fibers. We computed the friction coefficient and coefficient of viscosity of the amorphous phase to be in the order of 10⁻⁶ Ns/m and 10⁴ Ns/m², respectively, by extrapolating our simulation data to the viscous limit. Finally, we used a finite element method for the amorphous phase, solely based on parameters derived from molecular dynamics simulations including the newly determined coefficient of viscosity. With this model the time scales of stress relaxation, creep, and hysteresis were assessed, and found to be in line with the macroscopic time-dependent response of silk fibers. Our results suggest the amorphous phase to be the primary source of viscosity in silk and open up the avenue for finite element method studies of silk fiber mechanics including viscous effects.     

Submillisecond Elastic Recoil Reveals Molecular Origins of Fibrin Fiber Mechanics

Fibrin fibers form the structural scaffold of blood clots. Thus, their mechanical properties are of central importance to understanding hemostasis and thrombotic disease. Recent studies have revealed that fibrin fibers are elastomeric despite their high degree of molecular ordering. These results have inspired a variety of molecular models for fibrin’s elasticity, ranging from reversible protein unfolding to rubber-like elasticity. An important property that has not been explored is the timescale of elastic recoil, a parameter that is critical for fibrin’s mechanical function and places a temporal constraint on molecular models of fiber elasticity. Using high-frame-rate imaging and atomic force microscopy-based nanomanipulation, we measured the recoil dynamics of individual fibrin fibers and found that the recoil was orders of magnitude faster than anticipated from models involving protein refolding. We also performed steered discrete molecular-dynamics simulations to investigate the molecular origins of the observed recoil. Our results point to the unstructured αC regions of the otherwise structured fibrin molecule as being responsible for the elastic recoil of the fibers.     

Rate-Dependent Behavior of the Amorphous Phase of Spider Dragline Silk

The time-dependent stress-strain behavior of spider dragline silk was already observed decades ago, and has been attributed to the disordered sequences in silk proteins, which compose the soft amorphous matrix. However, the actual molecular origin and magnitude of internal friction within the amorphous matrix has remained inaccessible, because experimentally decomposing the mechanical response of the amorphous matrix from the embedded crystalline units is challenging. Here, we used atomistic molecular dynamics simulations to obtain friction forces for the relative sliding of peptide chains of Araneus diadematus spider silk within bundles of these chains as a representative unit of the amorphous matrix in silk fibers. We computed the friction coefficient and coefficient of viscosity of the amorphous phase to be in the order of 10⁻⁶ Ns/m and 10⁴ Ns/m², respectively, by extrapolating our simulation data to the viscous limit. Finally, we used a finite element method for the amorphous phase, solely based on parameters derived from molecular dynamics simulations including the newly determined coefficient of viscosity. With this model the time scales of stress relaxation, creep, and hysteresis were assessed, and found to be in line with the macroscopic time-dependent response of silk fibers. Our results suggest the amorphous phase to be the primary source of viscosity in silk and open up the avenue for finite element method studies of silk fiber mechanics including viscous effects.