From complex data to clear insights: visualizing molecular dynamics trajectories.

Advances in simulations, combined with technological developments in high-performance computing, have made it possible to produce a physically accurate dynamic representation of complex biological systems involving millions to billions of atoms over increasingly long simulation times. The analysis of these computed simulations is crucial, involving the interpretation of structural and dynamic data to gain insights into the underlying biological processes. However, this analysis becomes increasingly challenging due to the complexity of the generated systems with a large number of individual runs, ranging from hundreds to thousands of trajectories. This massive increase in raw simulation data creates additional processing and visualization challenges. Effective visualization techniques play a vital role in facilitating the analysis and interpretation of molecular dynamics simulations. In this paper, we focus mainly on the techniques and tools that can be used for visualization of molecular dynamics simulations, among which we highlight the few approaches used specifically for this purpose, discussing their advantages and limitations, and addressing the future challenges of molecular dynamics visualization.

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.

In silico validation of hyaluronic acid - drug conjugates based targeted drug delivery for the treatment of COVID-19.

The impact of COVID-19 urges scientists to develop targeted drug delivery to manage Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) viral infections with a fast recovery rate. The aim of the study is to develop Hyaluronic Acid (HA) drug conjugates of viral drugs to target two important enzymes (Mpro and PLpro) of SARS-CoV-2. Three antiviral drugs, namely Dexamethasone (DEX), Favipiravir (FAV), and Remdesivir (REM), were chosen for HA conjugation due to their reactive functional groups. Free forms of drugs (DEX, FAV, REM) and HA drug conjugates (HA-DEX, HA-FAV, HA-REM, HA-RHA, HA-RHE) were validated against Mpro (PDB ID 6LU7) and PLpro (PDB 7LLZ), which play an essential role in the replication and reproduction of the SARS-CoV-2 virus. The results of the present study revealed that HA-drug conjugates possess higher binding affinity and the best docking score towards the Mpro and PLpro target proteins of SARS-CoV-2 than their free forms of drugs. ADMET screening resulted that HA-drug conjugates exhibited better pharmacokinetic profiles than their pure forms of drugs. Further, molecular dynamic simulation studies, essential dynamics and free energy landscape analyses show that HA antiviral drug conjugates possess good trajectories and energy status, with the PLpro target protein (PDB 7LLZ) of SARS-CoV-2 through long-distance (500 ns) simulation screening. The research work recorded the best drug candidate for Cell-Targeted Drug Delivery (CTDD) for SARS-CoV-2-infected cells through hyaluronic acid conjugates of antiviral drugs.Communicated by Ramaswamy H. Sarma.

Sialic acids cleavage induced by elastin-derived peptides impairs the interaction between insulin and its receptor in adipocytes 3T3-L1.

The insulin receptor (IR) plays an important role in insulin signal transduction, the defect of which is believed to be the root cause of type 2 diabetes. In 3T3-L1 adipocytes as in other cell types, the mature IR is a heterotetrameric cell surface glycoprotein composed of two α subunits and two ß subunits. Our objective in our study, is to understand how the desialylation of N-glycan chains, induced by elastin-derived peptides, plays a major role in the function of the IR. Using the 3T3-L1 adipocyte line, we show that removal of the sialic acid from N-glycan chains (N893 and N908), induced by the elastin receptor complex (ERC) and elastin derived-peptides (EDPs), leads to a decrease in the autophosphorylation activity of the insulin receptor. We demonstrate by molecular dynamics approaches that the absence of sialic acids on one of these two sites is sufficient to generate local and general modifications of the structure of the IR. Biochemical approaches highlight a decrease in the interaction between insulin and its receptor when ERC sialidase activity is induced by EDPs. Therefore, desialylation by EDPs is synonymous with a decrease of IR sensitivity in adipocytes and could thus be a potential source of insulin resistance associated with diabetic conditions.

Unraveling the inhibitory mechanism of adenylyl cyclase 8E: New insights into regulatory pathways of cAMP signal integration.

Adenylyl Cyclase 8E (AC8E), which lacks part of M1 transmembrane domain, has been previously shown to dimerize with AC3 and down-regulate its activity, but the molecular mechanism of this inhibitory effect has remained elusive. Here, we first show that AC8E also inhibits AC2 and AC6, highlighting the functional importance of this novel regulatory mechanism in the cAMP signaling pathway across AC families. We then completed the partial structure of Bos taurus AC9 using combinations of comparative modeling and fold recognition methods, and used this as a template to build the first full 3D-models of AC8 and AC8E. These models evidenced that the lack of M1 transmembrane domain of AC8E shifts the N-terminal domain, which impacts the orientation of the helical domains, thus affecting the catalytic site. This was confirmed in living cells with cAMP imaging, where we showed that the N-terminal domain is required for reducing cAMP production. Our data also show that AC8E prevents the translocation of other ACs towards the plasma membrane, further reducing the cAMP responsiveness to extracellular signals. This newly discovered dual inhibitory mechanism provides an additional level of regulation of cAMP-dependent signals integration.

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.

How molecular modelling can better broaden the understanding of glycosylations.

Glycosylations are among the most ubiquitous post-translational modifications (PTMs) in proteins, and the effects of their perturbations are seen in various diseases such as cancers, diabetes and arthritis to name a few. Yet they remain one of the most enigmatic aspects of protein structure and function. On the other hand, molecular modelling techniques have been rapidly bridging this knowledge gap since the last decade. In this review, we discuss how these techniques have proven to be indispensable for a better understanding of the role of glycosylations in glycoprotein structure and function.

Early Alterations of Intra-Mural Elastic Lamellae Revealed by Synchrotron X-ray Micro-CT Exploration of Diabetic Aortas.

Diabetes is a major concern of our society as it affects one person out of 11 around the world. Elastic fiber alterations due to diabetes increase the stiffness of large arteries, but the structural effects of these alterations are poorly known. To address this issue, we used synchrotron X-ray microcomputed tomography with in-line phase contrast to image in three dimensions C57Bl6J (control) and db/db (diabetic) mice with a resolution of 650 nm/voxel and a field size of 1.3 mm3. Having previously shown in younger WT and db/db mouse cohorts that elastic lamellae contain an internal supporting lattice, here we show that in older db/db mice the elastic lamellae lose this scaffold. We coupled this label-free method with automated image analysis to demonstrate that the elastic lamellae from the arterial wall are structurally altered and become 11% smoother (286,665 measurements). This alteration suggests a link between the loss of the 3D lattice-like network and the waviness of the elastic lamellae. Therefore, waviness measurement appears to be a measurable elasticity indicator and the 3D lattice-like network appears to be at the origin of the existence of this waviness. Both could be suitable indicators of the overall elasticity of the aorta.

Anti-Toxoplasma gondii effect of lupane-type triterpenes from the bark of black alder (Alnus glutinosa) and identification of a potential target by reverse docking.

Toxoplasmosis is a worldwide parasitosis that is generally benign. The infestation may pose a risk to immunocompromized patients and to fetuses when pregnant women have recently seroconverted. Current treatments have numerous side effects and chemoresistance is emerging, hence the need to find new anti-Toxoplasma gondii substances. This study focuses on the antiparasitic potential of lupane-type pentacyclic triterpenes isolated from the bark of black alder (Alnus glutinosa), as well as the hypothesis of their macromolecular target by an original method of reverse docking. Among the isolated triterpenes, betulone was the most active compound with an IC50 of 2.7 ± 1.2 µM, a CC50 greater than 80 µM, and a selectivity index of over 29.6. An additional study of the anti-T. gondii potential of commercially available compounds (betulonic acid methyl ester and betulonic acid) showed the important role of the C3 ketone function and the C28 oxidation level on the lupane-type triterpene in the antiparasitic activity since their IC50 and CC50 were similar to that of betulone. Finally, the most active compounds were subjected to the AMIDE reverse docking workflow. A dataset of 87 T. gondii proteins from the Protein Data Bank was created. It identified calcium-dependent protein kinase CDPK3 as the most likely target of betulin derivatives.

TITLE: Potentiel anti-Toxoplasma gondii de triterpènes de type lupane de l'écorce de l'aulne glutineux, Alnus glutinosa, et identification d'une cible potentielle par docking inverse. ABSTRACT: La toxoplasmose est une parasitose mondiale, généralement bénigne. Les personnes à risque sont les patients immunodéprimés et les fœtus chez les femmes enceintes nouvellement séroconverties. Les traitements actuels ont de nombreux effets secondaires et des phénomènes de chimiorésistance apparaissent, d'où la nécessité de trouver de nouvelles substances actives contre T. gondii. Cette étude porte sur le potentiel antiparasitaire des triterpènes pentacycliques de type lupane isolés de l'écorce de l'aulne glutineux (Alnus glutinosa) et formule une hypothèse quant à leur cible protéique par l'utilisation d'une méthode originale de docking inverse. Parmi les triterpènes isolés, la bétulone s'est révélée être la plus active avec une CI50 de 2,7 µM ± 1,2 µM, une CC50 supérieure à 80 µM et un indice de sélectivité supérieur à 29,6. L'étude complémentaire du potentiel anti-T. gondii de composés disponibles commercialement et analogues à la bétulone (acide bétulonique et methyl ester de l'acide bétulonique) a montré le rôle important de la fonction cétone en C3 et du degré d'oxydation de la position 28 du squelette triterpénique de type lupane dans l'activité antiparasitaire puisque leurs CI50 et CC50 étaient similaires aux valeurs rencontrées pour la bétulone. Enfin, les composés les plus actifs ont été soumis au flux de travail de docking inverse d'AMIDE. Un ensemble de 87 protéines de T. gondii de la Protein Data Bank a été créé. La protéine kinase calcium dépendante CDPK3 a été identifiée comme la cible la plus probable des dérivés de la bétuline.

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.

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.

X-ray microtomography reveals a lattice-like network within aortic elastic lamellae.

The arterial wall consists of three concentric layers: intima, media, and adventitia. Beyond their resident cells, these layers are characterized by an extracellular matrix (ECM), which provides both biochemical and mechanical support. Elastin, the major component of arterial ECM, is present in the medial layer and organized in concentric elastic lamellae that confer resilience to the wall. We explored the arterial wall structures from C57Bl6 (control), db/db (diabetic), and ApoE-/- (atherogenic) mice aged 3 months using synchrotron X-ray computed microtomography on fixed and unstained tissues with a large image field (8 mm3 ). This approach combined a good resolution (0.83 µm/voxel), large 3D imaging field. and an excellent signal to noise ratio conferred by phase-contrast imaging. We determined from 2D virtual slices that the thickness of intramural ECM structures was comparable between strains but automated image analysis of the 3D arterial volumes revealed a lattice-like network within concentric elastic lamellae. We hypothesize that this network could play a role in arterial mechanics. This work demonstrates that phase-contrast synchrotron X-ray computed microtomography is a powerful technique which to characterize unstained soft tissues.

Effects of changes in glycan composition on glycoprotein dynamics: example of N-glycans on insulin receptor.

Glycosylation is among the most common post-translational modifications in proteins, although it is observed in only about 10% of all the protein structures in protein data bank (PDB). Modifications of sugar composition in glycoproteins profoundly impact the overall physiology of the organism. One such example is the development of insulin resistance, which has been attributed to the removal of sialic acid residues from N-glycans of insulin receptor (IR) from various experimental studies. How such modifications affect the glycan-glycoprotein dynamics, and ultimately their function is not clearly understood to date. In this study, we performed molecular dynamics simulations of glycans in different environments. We studied the effects of removal of sialic acid on the glycan, as well as on the dynamics of leucine-rich repeat L1 domain of the IR ectodomain. We observed perturbations in L1 domain dynamics as a result of the removal of sialic acid. The perturbations include an increase in the flexibility of insulin-binding residues, which may affect insulin binding with IR. These changes are accompanied by perturbations in glycan-protein interactions and perturbation of long-range allosteric dynamics. Our observations will further aid in understanding the role of sugars in maintaining homeostasis and how changes in glycan composition may lead to perturbations in homeostasis, ultimately leading to conditions such as insulin resistance.

AMIDE v2: High-Throughput Screening Based on AutoDock-GPU and Improved Workflow Leading to Better Performance and Reliability.

Molecular docking is widely used in computed drug discovery and biological target identification, but getting fast results can be tedious and often requires supercomputing solutions. AMIDE stands for AutoMated Inverse Docking Engine. It was initially developed in 2014 to perform inverse docking on High Performance Computing. AMIDE version 2 brings substantial speed-up improvement by using AutoDock-GPU and by pulling a total revision of programming workflow, leading to better performances, easier use, bug corrections, parallelization improvements and PC/HPC compatibility. In addition to inverse docking, AMIDE is now an optimized tool capable of high throughput inverse screening. For instance, AMIDE version 2 allows acceleration of the docking up to 12.4 times for 100 runs of AutoDock compared to version 1, without significant changes in docking poses. The reverse docking of a ligand on 87 proteins takes only 23 min on 1 GPU (Graphics Processing Unit), while version 1 required 300 cores to reach the same execution time. Moreover, we have shown an exponential acceleration of the computation time as a function of the number of GPUs used, allowing a significant reduction of the duration of the inverse docking process on large datasets.

Blending Ferulic Acid Derivatives and Polylactic Acid into Biobased and Transparent Elastomeric Materials with Shape Memory Properties.

Thanks to its remarkable properties such as sustainability, compostability, biocompatibility, and transparency, poly-l-lactic acid (PLA) would be a suitable replacement for oil-based polymers should it not suffer from low flexibility and poor toughness, restricting its use to rigid plastic by excluding elastomeric applications. Indeed, there are few fully biobased and biodegradable transparent elastomers-PLA-based or not-currently available. In the last decades, many strategies have been investigated to soften PLA and enhance its toughness and elongation at break by using plasticizers, oligomers, or polymers. This work shows how a ferulic acid-derived biobased additive (BDF) blends with a common rigid and brittle commercial grade of polylactic acid to provide a transparent non-covalently cross-linked elastomeric material with shape memory behavior exhibiting an elongation at break of 434% (vs 6% for pristine PLA). Through a structure-activity relationship analysis conducted with BDF analogues and a modeling study, we propose a mechanism based on π-π stacking to account for the elastomeric properties. Blending ferulic acid derivatives with polylactic acid generates a new family of fully sustainable transparent elastomeric materials with functional properties such as shape memory.

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.