Molecular dynamics simulation study on the interaction of collagen-like peptides with gelatinase-A (MMP-2)

De novo design of a stapled peptide targeting SARS-CoV-2 spike protein receptor-binding domain

Although effective vaccines have been developed against SARS-CoV-2, many regions in the world still have low rates of vaccination and new variants with mutations in the viral spike protein have reduced the effectiveness of most available vaccines and treatments. There is an urgent need for a drug to cure this disease and prevent infection. The SARS-CoV-2 virus enters the host cell through protein-protein interaction between the virus's spike protein and the host's angiotensin converting enzyme (ACE2). Using protein design software and molecular dynamics simulations, we have designed a 17-residue peptide (pep39), that binds to the spike protein receptor-binding domain (RBD) and blocks interaction of spike protein with ACE2. We have confirmed the binding activity of the designed peptide for the original spike protein and the delta variant spike protein using micro-cantilever and bio-layer interferometry (BLI) based methods. We also confirmed that pep39 strongly inhibits SARS-CoV-2 virus replication in Vero E6 cells. Taken together these data suggest that a newly designed spike protein RBD blocking peptide pep39 has a potential as a SARS-CoV-2 virus inhibitor.

Impact of collagen-like peptide (CLP) heterotrimeric triple helix design on helical thermal stability and hierarchical assembly: a coarse-grained molecular dynamics simulation study

Collagen-like peptides (CLP) are multifunctional materials garnering a lot of recent interest from the biomaterials community due to their hierarchical assembly and tunable physicochemical properties. In this work, we present a computational study that links the design of CLP heterotrimers to the thermal stability of the triple helix and their self-assembly into fibrillar aggregates and percolated networks. Unlike homotrimeric helices, the CLP heterotrimeric triple helices in this study are made of CLP strands of different chain lengths that result in 'sticky' ends with available hydrogen bonding groups. These 'sticky' ends at one end or both ends of the CLP heterotrimer then facilitate inter-helix hydrogen bonding leading to self-assembly into fibrils (clusters) and percolated networks. We consider the cases of three sticky end lengths - two, four, and six repeat units - present entirely on one end or split between two ends of the CLP heterotrimer. We observe in CLP heterotrimer melting curves generated using coarse grained Langevin dynamics simulations at low CLP concentration that increasing sticky end length results in lower melting temperatures for both one and two sticky ended CLP designs. At higher CLP concentrations, we observe non-monotonic trends in cluster sizes with increasing sticky end length with one sticky end but not for two sticky ends with the same number of available hydrogen bonding groups as the one sticky end; this nonmonotonicity stems from the formation of turn structures stabilized by hydrogen bonds at the single, sticky end for sticky end lengths greater than four repeat units. With increasing CLP concentration, heterotrimers also form percolated networks with increasing sticky end length with a minimum sticky end length of four repeat units required to observe percolation. Overall, this work informs the design of thermoresponsive, peptide-based biomaterials with desired morphologies using strand length and dispersity as a handle for tuning thermal stability and formation of supramolecular structures.

JED: a Java Essential Dynamics Program for comparative analysis of protein trajectories

BACKGROUND

Essential Dynamics (ED) is a common application of principal component analysis (PCA) to extract biologically relevant motions from atomic trajectories of proteins. Covariance and correlation based PCA are two common approaches to determine PCA modes (eigenvectors) and their eigenvalues. Protein dynamics can be characterized in terms of Cartesian coordinates or internal distance pairs. In understanding protein dynamics, a comparison of trajectories taken from a set of proteins for similarity assessment provides insight into conserved mechanisms. Comprehensive software is needed to facilitate comparative-analysis with user-friendly features that are rooted in best practices from multivariate statistics.

RESULTS

We developed a Java based Essential Dynamics toolkit called JED to compare the ED from multiple protein trajectories. Trajectories from different simulations and different proteins can be pooled for comparative studies. JED implements Cartesian-based coordinates (cPCA) and internal distance pair coordinates (dpPCA) as options to construct covariance (Q) or correlation (R) matrices. Statistical methods are implemented for treating outliers, benchmarking sampling adequacy, characterizing the precision of Q and R, and reporting partial correlations. JED output results as text files that include transformed coordinates for aligned structures, several metrics that quantify protein mobility, PCA modes with their eigenvalues, and displacement vector (DV) projections onto the top principal modes. Pymol scripts together with PDB files allow movies of individual Q- and R-cPCA modes to be visualized, and the essential dynamics occurring within user-selected time scales. Subspaces defined by the top eigenvectors are compared using several statistical metrics to quantify similarity/overlap of high dimensional vector spaces. Free energy landscapes can be generated for both cPCA and dpPCA.

CONCLUSIONS

JED offers a convenient toolkit that encourages best practices in applying multivariate statistics methods to perform comparative studies of essential dynamics over multiple proteins. For each protein, Cartesian coordinates or internal distance pairs can be employed over the entire structure or user-selected parts to quantify similarity/differences in mobility and correlations in dynamics to develop insight into protein structure/function relationships.

Molecular Dynamics Studies of Matrix Metalloproteases

Matrix metalloproteases are multidomain enzymes with a remarkable proteolytic activity located in the extracellular environment. Their catalytic activity and structural properties have been intensively studied during the last few decades using both experimental and theoretical approaches, but many open questions still remain. Extensive molecular dynamics simulations enable the sampling of the configurational space of a molecular system, thus contributing to the characterization of the structure, dynamics, and ligand binding properties of a particular MMP. Based on previous computational experience, we provide in this chapter technical and methodological guidelines that may be useful to and stimulate other researchers to perform molecular dynamics simulations to help address unresolved questions concerning the molecular mode of action of MMPs.

Molecular underpinnings of corneal angiogenesis: advances over the past decade

The cornea is maintained in an avascular state by maintaining an environment whereby anti-angiogenic factors take the upper hand over factors promoting angiogenesis. Many of the common pathologies affecting the cornea involve the disruption of such equilibrium and the shift towards new vessel formation, leading to corneal opacity and eventually-vision loss. Therefore it is of paramount importance that the molecular underpinnings of corneal neovascularization (CNV) be clearly understood, in order to develop better targeted treatments. This article is a review of the literature on the recent discoveries regarding pro-angiogenic factors of the cornea (such as vascular endothelial growth factors, fibroblast growth factor and matrix metalloproteinases) and anti-angiogenic factors of the cornea (such as endostatins and neostatins). Further, we review the molecular underpinnings of lymphangiogenesis, a process now known to be almost separate from (yet related to) hemangiogenesis.

Importance of the Linker Region in Matrix Metalloproteinase-1 Domain Interactions

Collagenolysis is catalyzed by enzymes from the matrix metalloproteinase (MMP) family, where one of the most studied is MMP-1. The X-ray crystallographic structure of MMP-1 complexed with a collagen-model triple-helical peptide (THP) provided important atomistic information, but few details on the effects of the conformational flexibility on catalysis. In addition, the role of the linker region between the catalytic (CAT) and hemopexin-like (HPX) domains was not defined. In order to reveal the dynamics and correlations of MMP-1 comprehensive atomistic molecular dynamics simulations of an MMP-1•THP complex was performed. To examine the role of the linker region for MMP-1 function simulations with linker regions from MT1-MMP/MMP-14 and MMP-13 replacing the MMP-1 linker region were performed. The MD studies were in good agreement with the experimental observation that in the MMP-1•THP X-ray crystallographic structure MMP-1 is in a "closed" conformation. MD revealed that the interactions of the THP with the both the CAT and HPX domains of MMP-1 are dynamic in nature, and the linker region of MMP-1 influences the interactions and dynamics of both the CAT and HPX domains and collagen binding to MMP-1.

Collagen interactions: Drug design and delivery

Collagen is a major component in a wide range of drug delivery systems and biomaterial applications. Its basic physical and structural properties, together with its low immunogenicity and natural turnover, are keys to its biocompatibility and effectiveness. In addition to its material properties, the collagen triple-helix interacts with a large number of molecules that trigger biological events. Collagen interactions with cell surface receptors regulate many cellular processes, while interactions with other ECM components are critical for matrix structure and remodeling. Collagen also interacts with enzymes involved in its biosynthesis and degradation, including matrix metalloproteinases. Over the past decade, much information has been gained about the nature and specificity of collagen interactions with its partners. These studies have defined collagen sequences responsible for binding and the high-resolution structures of triple-helical peptides bound to its natural binding partners. Strategies to target collagen interactions are already being developed, including the use of monoclonal antibodies to interfere with collagen fibril formation and the use of triple-helical peptides to direct liposomes to melanoma cells. The molecular information about collagen interactions will further serve as a foundation for computational studies to design small molecules that can interfere with specific interactions or target tumor cells. Intelligent control of collagen biological interactions within a material context will expand the effectiveness of collagen-based drug delivery.

A perspective on structural and computational work on collagen

Collagen is the single most abundant protein in the extracellular matrix in the animal kingdom, with remarkable structural and functional diversity and regarded one of the most useful biomaterials.

Proteins in aggregates functionally impact multiple neurodegenerative disease models by forming proteasome-blocking complexes

Age-dependent neurodegenerative diseases progressively form aggregates containing both shared components (e.g., TDP-43, phosphorylated tau) and proteins specific to each disease. We investigated whether diverse neuropathies might have additional aggregation-prone proteins in common, discoverable by proteomics. Caenorhabditis elegans expressing unc-54p/Q40::YFP, a model of polyglutamine array diseases such as Huntington's, accrues aggregates in muscle 2–6 days posthatch. These foci, isolated on antibody-coupled magnetic beads, were characterized by high-resolution mass spectrometry. Three Q40::YFP-associated proteins were inferred to promote aggregation and cytotoxicity, traits reduced or delayed by their RNA interference knockdown. These RNAi treatments also retarded aggregation/cytotoxicity in Alzheimer's disease models, nematodes with muscle or pan-neuronal Aβ_1–42 expression and behavioral phenotypes. The most abundant aggregated proteins are glutamine/asparagine-rich, favoring hydrophobic interactions with other random-coil domains. A particularly potent modulator of aggregation, CRAM-1/HYPK, contributed < 1% of protein aggregate peptides, yet its knockdown reduced Q40::YFP aggregates 72–86% (P < 10⁻⁶). In worms expressing Aβ_1–42, knockdown of cram-1 reduced β-amyloid 60% (P < 0.002) and slowed age-dependent paralysis > 30% (P < 10⁻⁶). In wild-type worms, cram-1 knockdown reduced aggregation and extended lifespan, but impaired early reproduction. Protection against seeded aggregates requires proteasome function, implying that normal CRAM-1 levels promote aggregation by interfering with proteasomal degradation of misfolded proteins. Molecular dynamic modeling predicts spontaneous and stable interactions of CRAM-1 (or human orthologs) with ubiquitin, and we verified that CRAM-1 reduces degradation of a tagged-ubiquitin reporter. We propose that CRAM-1 exemplifies a class of primitive chaperones that are initially protective and highly beneficial for early reproduction, but ultimately impair aggregate clearance and limit longevity.

Matrix metalloproteinase interactions with collagen and elastin

Most abundant in the extracellular matrix are collagens, joined by elastin that confers elastic recoil to the lung, aorta, and skin. These fibrils are highly resistant to proteolysis but can succumb to a minority of the matrix metalloproteinases (MMPs). Considerable inroads to understanding how such MMPs move to the susceptible sites in collagen and then unwind the triple helix of collagen monomers have been gained. The essential role in unwinding of the hemopexin-like domain of interstitial collagenases or the collagen binding domain of gelatinases is highlighted. Elastolysis is also facilitated by the collagen binding domain in the cases of MMP-2 and MMP-9, and remote exosites of the catalytic domain in the case of MMP-12.