Structural and physical basis for the elasticity of elastin.

Elastin function is to endow vertebrate tissues with elasticity so that they can adapt to local mechanical constraints. The hydrophobicity and insolubility of the mature elastin polymer have hampered studies of its molecular organisation and structure-elasticity relationships. Nevertheless, a growing number of studies from a broad range of disciplines have provided invaluable insights, and several structural models of elastin have been proposed. However, many questions remain regarding how the primary sequence of elastin (and the soluble precursor tropoelastin) governs the molecular structure, its organisation into a polymeric network, and the mechanical properties of the resulting material. The elasticity of elastin is known to be largely entropic in origin, a property that is understood to arise from both its disordered molecular structure and its hydrophobic character. Despite a high degree of hydrophobicity, elastin does not form compact, water-excluding domains and remains highly disordered. However, elastin contains both stable and labile secondary structure elements. Current models of elastin structure and function are drawn from data collected on tropoelastin and on elastin-like peptides (ELPs) but at the tissue level, elasticity is only achieved after polymerisation of the mature elastin. In tissues, the reticulation of tropoelastin chains in water defines the polymer elastin that bears elasticity. Similarly, ELPs require polymerisation to become elastic. There is considerable interest in elastin especially in the biomaterials and cosmetic fields where ELPs are widely used. This review aims to provide an up-to-date survey of/perspective on current knowledge about the interplay between elastin structure, solvation, and entropic elasticity.

Highlighting the hygroscopic capacities of apiogalacturonans.

To meet the needs of dehydrated skin, molecules with a high hygroscopic potential are necessary to hydrate it effectively and durably. In this context, we were interested in pectins, and more precisely in apiogalacturonans (AGA), a singular one that is currently only found in a few species of aquatic plants. As key structures in water regulation of these aquatic plants and thanks to their molecular composition and conformations, we hypothesized that they could have beneficial role for skin hydration. Spirodela polyrhiza is a duckweed known to be naturally rich in AGA. The aim of this study was to investigate the hygroscopic potential of AGA. Firstly, AGA models were built based on structural information obtained from previous experimental studies. Molecular dynamics (MD) simulations were performed, and the hygroscopic potential was predicted in silico by analyzing the frequency of interaction of water molecules with each AGA residue. Quantification of interactions identified the presence of 23 water molecules on average in contact with each residue of AGA. Secondly, the hygroscopic properties were investigated directly in vivo. Indeed, the water capture in the skin was measured in vivo by Raman microspectroscopy thanks to the deuterated water (D20) tracking. Investigations revealed that AGA significantly capture and retain more water in the epidermis and deeper than a placebo control. Not only do these original natural molecules interact with water molecules, but they capture and retain them efficiently in the skin.

Modelling and Simulations of Extracellular Glycoproteins.

While the knowledge of protein structure and function has seen vast advances in previous decades, the understanding of how their posttranslational modifications, such as glycosylations, influence their structure and function remains poor. However, advances in in silico methodologies to study glycosylations in recent past have enabled us to study this and understand the role of glycosylations in protein structure and function in ways that would not be possible by conventional experimental methods. In this chapter, we will demonstrate how to leverage these methodologies to study glycoproteins and their structural and dynamic properties using molecular modelling techniques.

Differential MMP-14 targeting by biglycan, decorin, fibromodulin, and lumican unraveled by in silico approach.

Small leucine-rich proteoglycans (SLRPs) are major regulators of extracellular matrix assembly and cell signaling. Lumican, a member of the SLRPs family, and its derived peptides were shown to possess antitumor activity by interacting directly with the catalytic domain of MMP-14 leading to the inhibition of its activity. The aim of the present report was to characterize by in silico three-dimensional (3D) modeling the structure and the dynamics of four SLRPs including their core protein and their specific polysaccharide chains to assess their capacity to bind to MMP-14 and to regulate its activity. Molecular docking experiments were performed to identify the specific amino acids of MMP-14 interacting with each of the four SLRPs. The inhibition of each SLRP (100 nM) on MMP-14 activity was measured and the constants of inhibition (Ki) were evaluated. The impact of the number of glycan chains, structures, and dynamics of lumican on the interaction with MMP-14 was assessed by molecular dynamics simulations. Molecular docking analysis showed that all SLRPs bind to MMP-14 through their concave face, but in different regions of the catalytic domain of MMP-14. Each SLRPs inhibited significantly the MMP-14 activity. Finally, molecular dynamics showed the role of glycan chains in interaction with MMP-14 and shielding effect of SLRPs. Altogether, the results demonstrated that each SLRP exhibited inhibition of MMP-14 activity. However, the differential targeting of MMP-14 by the SLRPs was shown to be related not only to the core protein conformation but also to the glycan chain structures and dynamics.

Multiscale modelling of the extracellular matrix.

The extracellular matrix is a complex three-dimensional network of molecules that provides cells with a complex microenvironment. The major constituents of the extracellular matrix such as collagen, elastin and associated proteins form supramolecular assemblies contributing to its physicochemical properties and organization. The structure of proteins and their supramolecular assemblies such as fibrils have been studied at the atomic level (e.g., by X-ray crystallography, Nuclear Magnetic Resonance and cryo-Electron Microscopy) or at the microscopic scale. However, many protein complexes are too large to be studied at the atomic level and too small to be studied by microscopy. Most extracellular matrix components fall into this intermediate scale, so-called the mesoscopic scale, preventing their detailed characterization. Simulation and modelling are some of the few powerful and promising approaches that can deepen our understanding of mesoscale systems. We have developed a set of modelling tools to study the self-organization of the extracellular matrix and large motion of macromolecules at the mesoscale level by taking advantage of the dynamics of articulated rigid bodies as a mean to study a larger range of motions at the cost of atomic resolution.

The Glypican-1/HGF/C-Met and Glypican-1/VEGF/VEGFR2 Ternary Complexes Regulate Hair Follicle Angiogenesis.

The hair renewal involves changes in the morphology of the hair follicle and its micro-vascularization. In alopecia, the hair cycle is accelerated, resulting in the formation of thinner and shorter hair. In addition, alopecia is associated with a decrease in the micro-vascularization of the hair follicles. In this study, the role of glypicans (GPCs) was analyzed in the regulation of the angiogenesis of human dermal microvascular endothelial cells (HDMEC). The analysis of glypican gene expression showed that GPC1 is the major glypican expressed by human keratinocytes of outer root sheath (KORS), human hair follicle dermal papilla cells (HHFDPC) and HDMEC. KORS were demonstrated to secrete VEGF and HGF. The HDMEC pseudotube formation was induced by KORS conditioned media (KORSCM). It was totally abrogated after GPC1 siRNA transfection of HDMEC. Moreover, when cleaved by phospholipase C (PLC), GPC1 promotes the proliferation of HDMEC. Finally, GPC1 was shown to interact directly with VEGFR2 or c-Met to regulate angiogenesis induced by the activation of these receptors. Altogether, these results showed that GPC1 is a key regulator of microvascular endothelial cell angiogenesis induced by VEGF and HGF secreted by KORS. Thus, GPC1 might constitute an interesting target to tackle alopecia in dermatology research.

Differential MMP-14 Targeting by Lumican-Derived Peptides Unraveled by In Silico Approach.

Lumican, a small leucine-rich proteoglycan (SLRP) of the extracellular matrix (ECM), displays anti-tumor properties through its direct interaction with MMP-14. Lumican-derived peptides, such as lumcorin (17 amino acids) or L9M (10 amino acids), are able to inhibit the proteolytic activity of MMP-14 and melanoma progression. This work aimed to visualize the interactions of lumican-derived peptides and MMP-14. Molecular modeling was used to characterize the interactions between lumican-derived peptides, such as lumcorin, L9M, and cyclic L9M (L9Mc, 12 amino acids), and MMP-14. The interaction of L9Mc with MMP-14 was preferential with the MT-Loop domain while lumcorin interacted more with the catalytic site. Key residues in the MMP-14 amino acid sequence were highlighted for the interaction between the inhibitory SLRP-derived peptides and MMP-14. In order to validate the in silico data, MMP-14 activity and migration assays were performed using murine B16F1 and human HT-144 melanoma cells. In contrast to the HT-144 melanoma cell line, L9Mc significantly inhibited the migration of B16F1 cells and the activity of MMP-14 but with less efficacy than lumican and lumcorin. L9Mc significantly inhibited the proliferation of B16F1 but not of HT-144 cells in vitro and primary melanoma tumor growth in vivo. Thus, the site of interaction between the domains of MMP-14 and lumcorin or L9Mc were different, which might explain the differences in the inhibitory effect of MMP-14 activity. Altogether, the biological assays validated the prediction of the in silico study. Possible and feasible improvements include molecular dynamics results.

LIMONADA: A database dedicated to the simulation of biological membranes.

Cellular membranes are composed of a wide diversity of lipid species in varying proportions and these compositions are representative of the organism, cellular type and organelle to which they belong. Because models of these molecular systems simulated by MD steadily gain in size and complexity, they are increasingly representative of specific compositions and behaviors of biological membranes. Due to the number of lipid species involved, of force fields and topologies and because of the complexity of membrane objects that have been simulated, LIMONADA has been developed as an open database allowing to handle the various aspects of lipid membrane simulation. LIMONADA presents published membrane patches with their simulation files and the cellular membrane it models. Their compositions are then detailed based on the lipid identification from LIPID MAPS database plus the lipid topologies and the force field used. LIMONADA is freely accessible on the web at https://limonada.univ-reims.fr/.

Correction to "A 3D Model of Human Lysyl Oxidase, A Cross-Linking Enzyme".

[This corrects the article DOI: 10.1021/acsomega.9b00317.].

Structural Analysis of Nonapeptides Derived from Elastin.

Elastin-derived peptides are released from the extracellular matrix remodeling by numerous proteases and seem to regulate many biological processes, notably cancer progression. The canonical elastin peptide is VGVAPG, which harbors the XGXXPG consensus pattern, allowing interaction with the elastin receptor complex located at the surface of cells. Besides these elastokines, another class of peptides has been identified. This group of bioactive elastin peptides presents the XGXPGXGXG consensus sequence, but the reason for their bioactivity remains unexplained. To better understand their nature and structure-function relationships, herein we searched the current databases for this nonapeptide motif and observed that the XGXPGXGXG elastin peptides define a specific group of tandemly repeated patterns. Further, we focused on four tandemly repeated human elastin nonapeptides, i.e., AGIPGLGVG, VGVPGLGVG, AGVPGLGVG, and AGVPGFGAG. These peptides were analyzed by means of optical spectroscopies and molecular dynamics. Ultraviolet-circular dichroism and Raman spectra are consistent with a mixture of ß-turn, ß-strand, and random-chain secondary elements in aqueous media. Quantitative analysis of their conformations suggested that turns corresponded to half of the total population of structural elements, whereas the remaining half were equally distributed between ß-strand and unordered chains. These distributions were confirmed by molecular dynamics simulations. Altogether, our data suggest that these highly dynamic peptides harbor a type II ß-turn located in their central part. We hypothesize that this structural element could explain their specific bioactivity.

Umbrella Visualization: A method of analysis dedicated to glycan flexibility with UnityMol.

N-glycosylation is a post-translational modification heavily impacting protein functions. Some alterations of glycosylation, such as sialic acid hydrolysis, are related to protein dysfunction. Because of their high flexibility and the many reactive groups of the glycan chains, studying glycans with in vitro methods is a challenging task. Molecular dynamics is a useful tool and probably the only one in biology able to overcome this problem and gives access to conformational information through exhaustive sampling. To better decipher the impact of N-glycans, the analysis and visualization of their influence over time on protein structure is a prerequisite. We developed the Umbrella Visualization, a graphical method that assigns the glycan intrinsic flexibility during a molecular dynamics trajectory. The density plot generated by this method brought relevant informations regarding glycans dynamics and flexibility, but needs further development in order to integrate an accurate description of the protein topology and its interactions. We propose here to transform this analysis method into a visualization mode in UnityMol. UnityMol is a molecular editor, viewer and prototyping platform, coded in C#. The new representation of glycan chains presented in this study takes into account both the main positions adopted by each antenna of a glycan and their statistical relevance. By displaying the collected data on the protein surface, one is then able to investigate the protein/glycan interactions.

A Three-Dimensional Model of Human Lysyl Oxidase, a Cross-Linking Enzyme.

Lysyl oxidase (LOX) is a cross-linking enzyme identified 50 years ago, but its 3D structure is still unknown. We have thus built a 3D model of human LOX by homology modeling using the X-ray structure of human lysyl oxidase-like 2 as a template. This model is the first one to recapitulate all known biochemical features of LOX, namely, the coordination of the copper ion and the formation of the lysine tyrosylquinone cofactor and the disulfide bridges. Furthermore, this model is stable during a 1 µs molecular dynamics simulation. The catalytic site is located in a groove surrounded by two loops. The distance between these loops fluctuated during the simulations, which suggests that the groove forms a hinge with a variable opening, which is able to accommodate the various sizes of LOX substrates. This 3D model is a pre-requisite to perform docking experiments with LOX substrates and other partners to identify binding sites and to design new LOX inhibitors specific for therapeutic purpose.

Small leucine-rich proteoglycans and matrix metalloproteinase-14: Key partners?

Small leucine-rich proteoglycans (SLRPs) are important regulators of extracellular matrix assembly and cell signaling. They are a family of proteoglycans that are present in extracellular matrix and that share in common multiple repeats of a leucine-rich structural motif. SLRPs have been identified as inhibitors of cancer progression by affecting MMPs, especially MMP-14 activity. Lumican, a member of the SLRPs family, and its derived peptides were shown to possess anti-tumor activity. Interestingly, it was demonstrated recently that lumican interacts directly with the catalytic domain of MMP-14 and inhibits its activity. The aim of this review was to summarize the interactions between SLRPs and MMPs with a special interest to lumican.

Mesoscopic Rigid Body Modelling of the Extracellular Matrix Self-Assembly.

The extracellular matrix (ECM) plays an important role in supporting tissues and organs. It even has a functional role in morphogenesis and differentiation by acting as a source of active molecules (matrikines). Many diseases are linked to dysfunction of ECM components and fragments or changes in their structures. As such it is a prime target for drugs. Because of technological limitations for observations at mesoscopic scales, the precise structural organisation of the ECM is not well-known, with sparse or fuzzy experimental observables. Based on the Unity3D game and physics engines, along with rigid body dynamics, we propose a virtual sandbox to model large biological molecules as dynamic chains of rigid bodies interacting together to gain insight into ECM components behaviour in the mesoscopic range. We have preliminary results showing how parameters such as fibre flexibility or the nature and number of interactions between molecules can induce different structures in the basement membrane. Using the Unity3D game engine and virtual reality headset coupled with haptic controllers, we immerse the user inside the corresponding simulation. Untrained users are able to navigate a complex virtual sandbox crowded with large biomolecules models in a matter of seconds.

Intrinsic dynamics study identifies two amino acids of TIMP-1 critical for its LRP-1-mediated endocytosis in neurons.

The tissue inhibitor of metalloproteinases-1 (TIMP-1) exerts inhibitory activity against matrix metalloproteinases and cytokine-like effects. We previously showed that TIMP-1 reduces neurite outgrowth in mouse cortical neurons and that this cytokine-like effect depends on TIMP-1 endocytosis mediated by the low-density lipoprotein receptor-related protein-1 (LRP-1). To gain insight into the interaction between TIMP-1 and LRP-1, we considered conformational changes that occur when a ligand binds to its receptor. TIMP-1 conformational changes have been studied using biomolecular simulations, and our results provide evidence for a hinge region that is critical for the protein movement between the N- and C-terminal TIMP-1 domains. In silico mutants have been proposed on residues F12 and K47, which are located in the hinge region. Biological analyses of these mutants show that F12A or K47A mutation does not alter MMP inhibitory activity but impairs the effect of TIMP-1 on neurite outgrowth. Interestingly, these mutants bind to LRP-1 but are not endocytosed. We conclude that the intrinsic dynamics of TIMP-1 are not involved in its binding to LRP-1 but rather in the initiation of endocytosis and associated biological effects.

Impact of sialic acids on the molecular dynamic of bi-antennary and tri-antennary glycans.

Sialic acids (SA) are monosaccharides that can be located at the terminal position of glycan chains on a wide range of proteins. The post-translational modifications, such as N-glycan chains, are fundamental to protein functions. Indeed, the hydrolysis of SA by specific enzymes such as neuraminidases can lead to drastic modifications of protein behavior. However, the relationship between desialylation of N-glycan chains and possible alterations of receptor function remains unexplored. Thus, the aim of the present study is to establish the impact of SA removal from N-glycan chains on their conformational behavior. We therefore undertook an in silico investigation using molecular dynamics to predict the structure of an isolated glycan chain. We performed, for the first time, 3 independent 500 ns simulations on bi-antennary and tri-antennary glycan chains displaying or lacking SA. We show that desialylation alters both the preferential conformation and the flexibility of the glycan chain. This study suggests that the behavior of glycan chains induced by presence or absence of SA may explain the changes in the protein function.

Identification of TAX2 peptide as a new unpredicted anti-cancer agent.

The multi-modular glycoprotein thrombospondin-1 (TSP-1) is considered as a key actor within the tumor microenvironment. Besides, TSP-1 binding to CD47 is widely reported to regulate cardiovascular function as it promotes vasoconstriction and angiogenesis limitation. Therefore, many studies focused on targeting TSP-1:CD47 interaction, aiming for up-regulation of physiological angiogenesis to enhance post-ischemia recovery or to facilitate engraftment. Thus, we sought to identify an innovative selective antagonist for TSP-1:CD47 interaction. Protein-protein docking and molecular dynamics simulations were conducted to design a novel CD47-derived peptide, called TAX2. TAX2 binds TSP-1 to prevent TSP-1:CD47 interaction, as revealed by ELISA and co-immunoprecipitation experiments. Unexpectedly, TAX2 inhibits in vitro and ex vivo angiogenesis features in a TSP-1-dependent manner. Consistently, our data highlighted that TAX2 promotes TSP-1 binding to CD36-containing complexes, leading to disruption of VEGFR2 activation and downstream NO signaling. Such unpredicted results prompted us to investigate TAX2 potential in tumor pathology. A multimodal imaging approach was conducted combining histopathological staining, MVD, MRI analysis and µCT monitoring for tumor angiography longitudinal follow-up and 3D quantification. TAX2 in vivo administrations highly disturb syngeneic melanoma tumor vascularization inducing extensive tumor necrosis and strongly inhibit growth rate and vascularization of human pancreatic carcinoma xenografts in nude mice.

Tuning self-assembly in elastin-derived peptides.

Elastin-derived peptides are gaining increasing interest as potential biomaterials. Previous studies have demonstrated that short elastin-derived peptides are able to self-assemble into fibrils as the entire elastin protein. The motif responsible for that is the XGGZG motif at least three-fold repeated. In this work we have synthesized and studied, at molecular and supramolecular levels, four pentadecapeptides obtained by switching the X and Z residue with leucine and/or valine. We found that the four peptides formed different supramolecular structures corresponding to specific molecular conformations. Our results show that not only the residue type but also the exact position occupied by the residue in the motif is crucial in driving the self-aggregation. The aim of this work is to provide the basis for designing elastin-derived peptides with tunable supramolecular architecture.

A crucial role for Lyn in TIMP-1 erythroid cell survival signalling pathway.

TIMP-1, a well-known MMP inhibitor, displays other biological activities such as cell survival, proliferation and differentiation in hematopoietic cells. In this report, we investigated the role of the Src-related kinase Lyn in TIMP-1 induced UT-7 erythroleukemic cell survival. We showed that (i) tyrosine 507 of Lyn was dephosphorylated and Lyn kinase activity enhanced by TIMP-1, (ii) Lyn silencing suppressed TIMP-1 anti-apoptotic activity and (iii) Lyn was activated upstream the JAK2/PI 3-kinase/Akt pathway. Our data suggest a novel role for Lyn in erythroid cell survival.

Large-scale proteomic analysis of the grapevine leaf apoplastic fluid reveals mainly stress-related proteins and cell wall modifying enzymes.

BACKGROUND: The extracellular space or apoplast forms a path through the whole plant and acts as an interface with the environment. The apoplast is composed of plant cell wall and space within which apoplastic fluid provides a means of delivering molecules and facilitates intercellular communications. However, the apoplastic fluid extraction from in planta systems remains challenging and this is particularly true for grapevine (Vitis vinifera L.), a worldwide-cultivated fruit plant. Large-scale proteomic analysis reveals the protein content of the grapevine leaf apoplastic fluid and the free interactive proteome map considerably facilitates the study of the grapevine proteome. RESULTS: To obtain a snapshot of the grapevine apoplastic fluid proteome, a vacuum-infiltration-centrifugation method was optimized to collect the apoplastic fluid from non-challenged grapevine leaves. Soluble apoplastic protein patterns were then compared to whole leaf soluble protein profiles by 2D-PAGE analyses. Subsequent MALDI-TOF/TOF mass spectrometry of tryptically digested protein spots was used to identify proteins. This large-scale proteomic analysis established a well-defined proteomic map of whole leaf and leaf apoplastic soluble proteins, with 223 and 177 analyzed spots, respectively. All data arising from proteomic, MS and MS/MS analyses were deposited in the public database world-2DPAGE. Prediction tools revealed a high proportion of (i) classical secreted proteins but also of non-classical secreted proteins namely Leaderless Secreted Proteins (LSPs) in the apoplastic protein content and (ii) proteins potentially involved in stress reactions and/or in cell wall metabolism. CONCLUSIONS: This approach provides free online interactive reference maps annotating a large number of soluble proteins of the whole leaf and the apoplastic fluid of grapevine leaf. To our knowledge, this is the first detailed proteome study of grapevine apoplastic fluid providing a comprehensive overview of the most abundant proteins present in the apoplast of grapevine leaf that could be further characterized in order to elucidate their physiological function.