Mapping intermolecular interactions and active site conformations: from human MMP-1 crystal structure to molecular dynamics free energy calculations

Molecular Dynamics Simulations of Matrix Metalloproteinase 13 and the Analysis of the Specificity Loop and the S1'-Site

The specificity loop of Matrix Metalloproteinases (MMPs) is known to regulate recognition of their substrates, and the S1'-site surrounded by the loop is a unique place to address the selectivity of ligands toward each MMP. Molecular dynamics (MD) simulations of apo-MMP-13 and its complex forms with various ligands were conducted to identify the role of the specificity loop for the ligand binding to MMP-13. The MD simulations showed the dual role of T247 as a hydrogen bond donor to the ligand, as well as a contributor to the formation of the van der Waal surface area, with T245 and K249 on the S1'-site. The hydrophobic surface area mediated by T247 blocks the access of water molecules to the S1'-site of MMP-13 and stabilizes the ligand in the site. The F252 residue is flexible in order to search for the optimum location in the S1'-site of the apo-MMP-13, but once a ligand binds to the S1'-site, it can form offset π-π or edge-to-π stacking interactions with the ligand. Lastly, H222 and Y244 provide the offset π-π and π-CH(Cβ) interactions on each side of the phenyl ring of the ligand, and this sandwiched interaction could be critical for the ligand binding to MMP-13.

Bioengineered liver crosslinked with nano-graphene oxide enables efficient liver regeneration via MMP suppression and immunomodulation

Decellularized extracellular matrix scaffold, widely utilized for organ engineering, often undergoes matrix decomposition after transplantation and produces byproducts that cause inflammation, leading to clinical failure. Here we propose a strategy using nano-graphene oxide to modify the biophysical properties of decellularized liver scaffolds. Notably, we demonstrate that scaffolds crosslinked with nano-graphene oxide show high resistance to enzymatic degradation via direct inhibition of matrix metalloproteinase activity and increased mechanical rigidity. We find that M2-like macrophage polarization is promoted within the crosslinked scaffolds, which reduces graft-elicited inflammation. Moreover, we show that low activities of matrix metalloproteinases, attributed to both nano-graphene oxide and tissue inhibitors of metalloproteinases expressed by M2c, can protect the crosslinked scaffolds against in vivo degradation. Lastly, we demonstrate that bioengineered livers fabricated with the crosslinked scaffolds remain functional, thereby effectively regenerating damaged livers after transplantation into liver failure mouse models. Overall, nano-graphene oxide crosslinking prolongs allograft survival and ultimately improves therapeutic effects of bioengineered livers, which offer an alternative for donor organs.

Variation in Leishmania chemokine suppression driven by diversification of the GP63 virulence factor

Leishmaniasis is a neglected tropical disease with diverse outcomes ranging from self-healing lesions, to progressive non-healing lesions, to metastatic spread and destruction of mucous membranes. Although resolution of cutaneous leishmaniasis is a classic example of type-1 immunity leading to self-healing lesions, an excess of type-1 related inflammation can contribute to immunopathology and metastatic spread. Leishmania genetic diversity can contribute to variation in polarization and robustness of the immune response through differences in both pathogen sensing by the host and immune evasion by the parasite. In this study, we observed a difference in parasite chemokine suppression between the Leishmania (L.) subgenus and the Viannia (V.) subgenus, which is associated with severe immune-mediated pathology such as mucocutaneous leishmaniasis. While Leishmania (L.) subgenus parasites utilize the virulence factor and metalloprotease glycoprotein-63 (gp63) to suppress the type-1 associated host chemokine CXCL10, L. (V.) panamensis did not suppress CXCL10. To understand the molecular basis for the inter-species variation in chemokine suppression, we used in silico modeling to identify a putative CXCL10-binding site on GP63. The putative CXCL10 binding site is in a region of gp63 under significant positive selection, and it varies from the L. major wild-type sequence in all gp63 alleles identified in the L. (V.) panamensis reference genome. Mutating wild-type L. (L.) major gp63 to the L. (V.) panamensis sequence at the putative binding site impaired cleavage of CXCL10 but not a non-specific protease substrate. Notably, Viannia clinical isolates confirmed that L. (V.) panamensis primarily encodes non-CXCL10-cleaving gp63 alleles. In contrast, L. (V.) braziliensis has an intermediate level of activity, consistent with this species having more equal proportions of both alleles. Our results demonstrate how parasite genetic diversity can contribute to variation in immune responses to Leishmania spp. infection that may play critical roles in the outcome of infection.

Leishmaniasis is a neglected tropical disease caused by Leishmania parasites and spread by the bites of infected sand flies. Most cases of leishmaniasis present as self-healing sores that are resolved by a balanced immune response. Other cases of leishmaniasis involve spread to sites distant from the original bite, including damage of the inner surfaces of the mouth and nose. These cases of leishmaniasis involve an excessive immune response. Leishmania parasites produce virulence factor proteins, such as GP63, to trick the immune system into mounting a weaker response. GP63 specifically degrades signaling proteins that attract and activate certain immune cells. Here, we demonstrate that Leishmania parasite species have evolved to differ in their ability to degrade signaling proteins. In Leishmania species known to cause more immune-mediated tissue damage, the GP63 virulence factor has evolved to not degrade specific immune signaling proteins, thus attracting, and activating more immune cells. Our results demonstrate how diversity among Leishmania parasite species can contribute to variation in immune responses that may play critical roles in the outcome of infection.

Allosteric Communications between Domains Modulate the Activity of Matrix Metalloprotease-1

An understanding of the structure-dynamics relationship is essential for understanding how a protein works. Prior research has shown that the activity of a protein correlates with intradomain dynamics occurring at picosecond to millisecond timescales. However, the correlation between interdomain dynamics and the function of a protein is poorly understood. Here, we show that communications between the catalytic and hemopexin domains of matrix metalloprotease-1 (MMP1) on type 1 collagen fibrils correlate with its activity. Using single-molecule Förster resonance energy transfer, we identified functionally relevant open conformations in which the two MMP1 domains are well separated, which were significantly absent for catalytically inactive point mutant (E219Q) of MMP1 and could be modulated by an inhibitor or an enhancer of activity. The observed relevance of open conformations resolves the debate about the roles of open and closed MMP1 structures in function. We fitted the histograms of single-molecule Förster resonance energy transfer values to a sum of two Gaussians and the autocorrelations to an exponential and power law. We used a two-state Poisson process to describe the dynamics and calculate the kinetic rates from the fit parameters. All-atom and coarse-grained simulations reproduced some of the experimental features and revealed substrate-dependent MMP1 dynamics. Our results suggest that an interdomain separation facilitates opening up the catalytic pocket so that the collagen chains come closer to the MMP1 active site. Coordination of functional conformations at different parts of MMP1 occurs via allosteric communications that can take place via interactions mediated by collagen even if the linker between the domains is absent. Modeling dynamics as a Poisson process enables connecting the picosecond timescales of molecular dynamics simulations with the millisecond timescales of single-molecule measurements. Water-soluble MMP1 interacting with water-insoluble collagen fibrils poses challenges for biochemical studies that the single-molecule tracking can overcome for other insoluble substrates. Interdomain communications are likely important for multidomain proteins.

Velvet antler polypeptide partially rescue facet joint osteoarthritis-like phenotype in adult β-catenin conditional activation mice

BACKGROUND

Wnt/β-catenin signaling pathway is closely related to osteoarthritis. In our preliminary study, β-catenin conditional activation (cAct) mice that specifically over-express β-catenin gene in cartilage chondrocyte exhibits osteoarthritis-like phenotype in the lumbar disc and knee joint. Therefore, we used the mice to model FJ-OA and test the potential curative effect of Velvet Antler Polypeptide (VAP) on this mice model.

METHODS

We tested the effect of VAP on β-catenin conditional activation mice, and used Cre negative littermates as controls. Micro-CT, histology and histomorphometry analysis were performed to evaluate the curative effect of VAP on mice facet joint-like phenotype. Expression of β-catenin and collagen II was detected by immunohistochemistry (IHC) and western-blot., MMP13, ADAMTS4 and ADAMTS5 was detected by immunofluorescence (IF). RT-PCR analysis was preformed to detect mRNA expression of cartilage degrading enzymes, such as MMP13, ADAMTS4 and ADAMTS5.

RESULTS

Results of micro-CT (μCT) analysis showed that VAP could partially reverse lumbar disc osteophyte formation observed in β-catenin(ex3)^Col2ER mice. Histology data revealed VAP partially improved facet joint cartilage tissue invades. Histomorphometry analysis showed an increase in total cartilage area after VAP treatment. IHC show that VAP reduced β-catenin protein levels and moderately up-regulated collagen II protein levels. RT-PCR and IF data showed that VAP down-regulated the expression of extracellular matrix synthesis (ECM) degradation enzymes MMP13, ADAMTS4 and ADAMTS5.

CONCLUSION

Taken together, VAP may modulate ECM by inhibits MMP13, ADAMTS4 and ADAMTS5 via Wnt /β-catenin signaling pathway. Velvet Antler Polypeptide may be a potential medicine for FJ-OA.

Delineating the active site architecture of G9a lysine methyltransferase through substrate and inhibitor binding mode analysis: a molecular dynamics study

Mono- and di-methylation of the H3K9 residue in the histone tail by G9a lysine methyltransferase is associated with transcriptional suppression of genes. Here, we use molecular dynamics simulation and free energy calculations of five different modified/mutated G9a substrate peptides to elucidate the rationale behind the substrate binding to G9a. We also investigated the binding energy contribution based architecture of the active site of G9a to understand substrate and inhibitor binding. Wild-type peptide (H3K9) shows better binding affinity than mono- and di-methylated lysine (K9) and other modified peptides (K9A and R8A). Arg8 of the substrate peptide is crucial for determining the degree of conformational freedom/stability of the wild-type substrate peptide, as well as binding to G9a. Our results also suggest that the G9a active site is segregated into energy rich and low regions, and the energy rich region alone is used by the inhibitors for binding. These insights into the active site architecture should be taken into consideration in virtual screening experiments designed to discover novel inhibitors for G9a. In particular, compounds that could interact with the six residues of G9a - Asp1074, Asp1083, Leu1086, Asp1088, Tyr1154 and Phe1158 - should be preferentially tested in G9a inhibition biological assays. Communicated by Ramaswamy H. Sarma.

Hierarchical virtual screening of the dual MMP-2/HDAC-6 inhibitors from natural products based on pharmacophore models and molecular docking

The dual-target inhibitors tend to improve the response rate in treating tumors, comparing with the single-target inhibitors. Matrix metalloproteinase-2 (MMP-2) and histone deacetylase-6 (HDAC-6) are attractive targets for cancer therapy. In this study, the hierarchical virtual screening of dual MMP-2/HDAC-6 inhibitors from natural products is investigated. The pharmacophore model of MMP-2 inhibitors is built based on ligands, but the pharmacophore model of HDAC-6 inhibitors is built based on the experimental crystal structures of multiple receptor-ligand complexes. The reliability of these two pharmacophore models is validated subsequently. The hierarchical virtual screening, combining these two different pharmacophore models of MMP-2 and HDAC-6 inhibitors with molecular docking, is carried out to identify the dual MMP-2/HDAC-6 inhibitors from a database of natural products. The four potential dual MMP-2/HDAC-6 inhibitors of natural products, STOCK1 N-46177, STOCK1 N-52245, STOCK1 N-55477, and STOCK1 N-69706, are found. The studies of binding modes show that the screened four natural products can simultaneously well bind with the MMP-2 and HDAC-6 active sites by different kinds of interactions, to inhibit the MMP-2 and HDAC-6 activities. In addition, the ADMET properties of screened four natural products are assessed. These found dual MMP-2/HDAC-6 inhibitors of natural products could serve as the lead compounds for designing the new dual MMP-2/HDAC-6 inhibitors having higher biological activities by carrying out structural modifications and optimizations in the future studies.

ForceGen: atomic covalent bond value derivation for Gromacs

A large number of crystallographic protein structures include ligands, small molecules and post-translational modifications. Atomic bond force values for computational atomistic models of post-translational or non-standard amino acids, metal binding active sites, small molecules and drug molecules are not readily available in most simulation software packages. We present ForceGen, a Java tool that extracts the bond stretch and bond angle force values and equilibrium values from the Hessian of a Gaussian vibrational frequency analysis. The parameters are compatible with force fields derived using the second order tensor of the Hessian. The output is formatted with the Gromacs topology in mind. This study further demonstrates the use of ForceGen over the quantum mechanically derived structures of a small organic solvent, a naturally occurring protein crosslink derived from two amino acids following post-translational modification and the amino acid ligands of a zinc ion. We then derive Laplacian bond orders to understand how the resulting force values relate back to the quantum mechanical model. The parameterisation of the organic solvent, toluene, was verified using Molecular Mechanics simulations. The structural data from the simulation compared well with the quantum mechanical structure and the system density compared well with experimental values.