PCMRESP: A Method for Polarizable Force Field Parameter Development and Transferability of the Polarizable Gaussian Multipole Models Across Multiple Solvents.

The transferability of force field parameters is a crucial aspect of high-quality force fields. Previous investigations have affirmed the transferability of electrostatic parameters derived from polarizable Gaussian multipole models (pGMs) when applied to water oligomer clusters, polypeptides across various conformations, and different sequences. In this study, we introduce PCMRESP, a novel method for electrostatic parametrization in solution, intended for the development of polarizable force fields. We utilized this method to assess the transferability of three models: a fixed charge model and two variants of pGM models. Our analysis involved testing these models on 377 small molecules and 100 tetra-peptides in five representative dielectric environments: gas, diethyl ether, dichloroethane, acetone, and water. Our findings reveal that the inclusion of atomic polarization significantly enhances transferability and the incorporation of permanent atomic dipoles, in the form of covalent bond dipoles, leads to further improvements. Moreover, our tests on dual-solvent strategies demonstrate consistent transferability for all three models, underscoring the robustness of the dual-solvent approach. In contrast, an evaluation of the traditional HF/6-31G* method indicates poor transferability for the pGM-ind and pGM-perm models, suggesting the limitations of this conventional approach.

Assessment of Amino Acid Electrostatic Parametrizations of the Polarizable Gaussian Multipole Model.

Accurate parametrization of amino acids is pivotal for the development of reliable force fields for molecular modeling of biomolecules such as proteins. This study aims to assess amino acid electrostatic parametrizations with the polarizable Gaussian Multipole (pGM) model by evaluating the performance of the pGM-perm (with atomic permanent dipoles) and pGM-ind (without atomic permanent dipoles) variants compared to the traditional RESP model. The 100-conf-combterm fitting strategy on tetrapeptides was adopted, in which (1) all peptide bond atoms (-CO-NH-) share identical set of parameters and (2) the total charges of the two terminal N-acetyl (ACE) and N-methylamide (NME) groups were set to neutral. The accuracy and transferability of electrostatic parameters across peptides with varying lengths and real-world examples were examined. The results demonstrate the enhanced performance of the pGM-perm model in accurately representing the electrostatic properties of amino acids. This insight underscores the potential of the pGM-perm model and the 100-conf-combterm strategy for the future development of the pGM force field.

Dual function of Zika virus NS2B-NS3 protease.

Zika virus (ZIKV) serine protease, indispensable for viral polyprotein processing and replication, is composed of the membrane-anchored NS2B polypeptide and the N-terminal domain of the NS3 polypeptide (NS3pro). The C-terminal domain of the NS3 polypeptide (NS3hel) is necessary for helicase activity and contains an ATP-binding site. We discovered that ZIKV NS2B-NS3pro binds single-stranded RNA with a Kd of ~0.3 µM, suggesting a novel function. We tested various structural modifications of NS2B-NS3pro and observed that constructs stabilized in the recently discovered "super-open" conformation do not bind RNA. Likewise, stabilizing NS2B-NS3pro in the "closed" (proteolytically active) conformation using substrate inhibitors abolished RNA binding. We posit that RNA binding occurs when ZIKV NS2B-NS3pro adopts the "open" conformation, which we modeled using highly homologous dengue NS2B-NS3pro crystallized in the open conformation. We identified two positively charged fork-like structures present only in the open conformation of NS3pro. These forks are conserved across Flaviviridae family and could be aligned with the positively charged grove on NS3hel, providing a contiguous binding surface for the negative RNA strand exiting helicase. We propose a "reverse inchworm" model for a tightly intertwined NS2B-NS3 helicase-protease machinery, which suggests that NS2B-NS3pro cycles between open and super-open conformations to bind and release RNA enabling long-range NS3hel processivity. The transition to the closed conformation, likely induced by the substrate, enables the classical protease activity of NS2B-NS3pro.

Transferability of the Electrostatic Parameters of the Polarizable Gaussian Multipole Model.

Accuracy and transferability are the two highly desirable properties of molecular mechanical force fields. Compared with the extensively used point-charge additive force fields that apply fixed atom-centered point partial charges to model electrostatic interactions, polarizable force fields are thought to have the advantage of modeling the atomic polarization effects. Previous works have demonstrated the accuracy of the recently developed polarizable Gaussian multipole (pGM) models. In this work, we assessed the transferability of the electrostatic parameters of the pGM models with (pGM-perm) and without (pGM-ind) atomic permanent dipoles in terms of reproducing the electrostatic potentials surrounding molecules/oligomers absent from electrostatic parameterizations. Encouragingly, both the pGM-perm and pGM-ind models show significantly improved transferability than the additive model in the tests (1) from water monomer to water oligomer clusters; (2) across different conformations of amino acid dipeptides and tetrapeptides; (3) from amino acid tetrapeptides to longer polypeptides; and (4) from nucleobase monomers to Watson-Crick base pair dimers and tetramers. Furthermore, we demonstrated that the double-conformation fittings using amino acid tetrapeptides in the αR and ß conformations can result in good transferability not only across different tetrapeptide conformations but also from tetrapeptides to polypeptides with lengths ranging from 1 to 20 repetitive residues for both the pGM-ind and pGM-perm models. In addition, the observation that the pGM-ind model has significantly better accuracy and transferability than the point-charge additive model, even though they have an identical number of parameters, strongly suggest the importance of intramolecular polarization effects. In summary, this and previous works together show that the pGM models possess both accuracy and transferability, which are expected to serve as foundations for the development of next-generation polarizable force fields for modeling various polarization-sensitive biological systems and processes.

Accurate Reproduction of Quantum Mechanical Many-Body Interactions in Peptide Main-Chain Hydrogen-Bonding Oligomers by the Polarizable Gaussian Multipole Model.

A key advantage of polarizable force fields is their ability to model the atomic polarization effects that play key roles in the atomic many-body interactions. In this work, we assessed the accuracy of the recently developed polarizable Gaussian Multipole (pGM) models in reproducing quantum mechanical (QM) interaction energies, many-body interaction energies, as well as the nonadditive and additive contributions to the many-body interactions for peptide main-chain hydrogen-bonding conformers, using glycine dipeptide oligomers as the model systems. Two types of pGM models were considered, including that with (pGM-perm) and without (pGM-ind) permanent atomic dipoles. The performances of the pGM models were compared with several widely used force fields, including two polarizable (Amoeba13 and ff12pol) and three additive (ff19SB, ff15ipq, and ff03) force fields. Encouragingly, the pGM models outperform all other force fields in terms of reproducing QM interaction energies, many-body interaction energies, as well as the nonadditive and additive contributions to the many-body interactions, as measured by the root-mean-square errors (RMSEs) and mean absolute errors (MAEs). Furthermore, we tested the robustness of the pGM models against polarizability parameterization errors by employing alternative polarizabilities that are either scaled or obtained from other force fields. The results show that the pGM models with alternative polarizabilities exhibit improved accuracy in reproducing QM many-body interaction energies as well as the nonadditive and additive contributions compared with other polarizable force fields, suggesting that the pGM models are robust against the errors in polarizability parameterizations. This work shows that the pGM models are capable of accurately modeling polarization effects and have the potential to serve as templates for developing next-generation polarizable force fields for modeling various biological systems.

S-nitrosylated TDP-43 triggers aggregation, cell-to-cell spread, and neurotoxicity in hiPSCs and in vivo models of ALS/FTD.

Rare genetic mutations result in aggregation and spreading of cognate proteins in neurodegenerative disorders; however, in the absence of mutation (i.e., in the vast majority of "sporadic" cases), mechanisms for protein misfolding/aggregation remain largely unknown. Here, we show environmentally induced nitrosative stress triggers protein aggregation and cell-to-cell spread. In patient brains with amyotrophic lateral sclerosis (ALS)/frontotemporal dementia (FTD), aggregation of the RNA-binding protein TDP-43 constitutes a major component of aberrant cytoplasmic inclusions. We identify a pathological signaling cascade whereby reactive nitrogen species cause S-nitrosylation of TDP-43 (forming SNO-TDP-43) to facilitate disulfide linkage and consequent TDP-43 aggregation. Similar pathological SNO-TDP-43 levels occur in postmortem human FTD/ALS brains and in cell-based models, including human-induced pluripotent stem cell (hiPSC)-derived neurons. Aggregated TDP-43 triggers additional nitrosative stress, representing positive feed forward leading to further SNO-TDP-43 formation and disulfide-linked oligomerization/aggregation. Critically, we show that these redox reactions facilitate cell spreading in vivo and interfere with the TDP-43 RNA-binding activity, affecting SNMT1 and phospho-(p)CREB levels, thus contributing to neuronal damage in ALS/FTD disorders.

Interaction of the cryptic fragment of myelin basic protein with mitochondrial voltage-dependent anion-selective channel-1 affects cell energy metabolism.

In demyelinating nervous system disorders, myelin basic protein (MBP), a major component of the myelin sheath, is proteolyzed and its fragments are released in the neural environment. Here, we demonstrated that, in contrast with MBP, the cellular uptake of the cryptic 84-104 epitope (MBP84-104) did not involve the low-density lipoprotein receptor-related protein-1, a scavenger receptor. Our pull-down assay, mass spectrometry and molecular modeling studies suggested that, similar with many other unfolded and aberrant proteins and peptides, the internalized MBP84-104 was capable of binding to the voltage-dependent anion-selective channel-1 (VDAC-1), a mitochondrial porin. Molecular modeling suggested that MBP84-104 directly binds to the N-terminal α-helix located midway inside the 19 ß-blade barrel of VDAC-1. These interactions may have affected the mitochondrial functions and energy metabolism in multiple cell types. Notably, MBP84-104 caused neither cell apoptosis nor affected the total cellular ATP levels, but repressed the aerobic glycolysis (lactic acid fermentation) and decreased the l-lactate/d-glucose ratio (also termed as the Warburg effect) in normal and cancer cells. Overall, our findings implied that because of its interactions with VDAC-1, the cryptic MBP84-104 peptide invoked reprogramming of the cellular energy metabolism that favored enhanced cellular activity, rather than apoptotic cell death. We concluded that the released MBP84-104 peptide, internalized by the cells, contributes to the reprogramming of the energy-generating pathways in multiple cell types.

Resolving the Ligand-Binding Specificity in c-MYC G-Quadruplex DNA: Absolute Binding Free Energy Calculations and SPR Experiment.

We report the absolute binding free energy calculation and surface plasmon resonance (SPR) experiment for ligand binding with the c-MYC G-quadruplex DNA. The unimolecular parallel DNA G-quadruplex formed in nuclease hypersensitivity element III1 of the c-MYC gene promoter regulates the c-MYC transcription and is recognized as an emerging drug target for cancer therapy. Quindoline derivatives have been shown to stabilize the G-quadruplex and inhibit the c-MYC expression in cancer cells. NMR revealed two binding sites located at the 5' and 3' termini of the G-quadruplex. Questions about which site is more favored and the basis for the ligand-induced binding site formation remain unresolved. Here, we employ two absolute binding free energy methods, the double decoupling and the potential of mean force methods, to dissect the ligand-binding specificity in the c-MYC G-quadruplex. The calculated absolute binding free energies are in general agreement with the SPR result and suggest that quindoline has a slight preference for the 5' site. The flanking residues around the two sites undergo significant reorganization as the ligand unbinds, which provides evidence for ligand-induced binding pocket formation. The results help interpret experimental data and inform rational design of small molecules targeting the c-MYC G-quadruplex.

Characterization of the Zika virus two-component NS2B-NS3 protease and structure-assisted identification of allosteric small-molecule antagonists.

The recent re-emergence of Zika virus (ZIKV)1, a member of the Flaviviridae family, has become a global emergency. Currently, there are no effective methods of preventing or treating ZIKV infection, which causes severe neuroimmunopathology and is particularly harmful to the developing fetuses of infected pregnant women. However, the pathology induced by ZIKV is unique among flaviviruses, and knowledge of the biology of other family members cannot easily be extrapolated to ZIKV. Thus, structure-function studies of ZIKV proteins are urgently needed to facilitate the development of effective preventative and therapeutic agents. Like other flaviviruses, ZIKV expresses an NS2B-NS3 protease, which consists of the NS2B cofactor and the NS3 protease domain and is essential for cleavage of the ZIKV polyprotein precursor and generation of fully functional viral proteins. Here, we report the enzymatic characterization of ZIKV protease, and we identify structural scaffolds for allosteric small-molecule inhibitors of this protease. Molecular modeling of the protease-inhibitor complexes suggests that these compounds bind to the druggable cavity in the NS2B-NS3 protease interface and affect productive interactions of the protease domain with its cofactor. The most potent compound demonstrated efficient inhibition of ZIKV propagation in vitro in human fetal neural progenitor cells and in vivo in SJL mice. The inhibitory scaffolds could be further developed into valuable research reagents and, ultimately, provide a roadmap for the selection of efficient inhibitors of ZIKV infection.

Selective function-blocking monoclonal human antibody highlights the important role of membrane type-1 matrix metalloproteinase (MT1-MMP) in metastasis.

The invasion-promoting MT1-MMP is a cell surface-associated collagenase with a plethora of critical cellular functions. There is a consensus that MT1-MMP is a key protease in aberrant pericellular proteolysis in migrating cancer cells and, accordingly, a promising drug target. Because of high homology in the MMP family and a limited success in the design of selective small-molecule inhibitors, it became evident that the inhibitor specificity is required for selective and successful MT1-MMP therapies. Using the human Fab antibody library (over 1.25×109 individual variants) that exhibited the extended, 23-27 residue long, VH CDR-H3 segments, we isolated a panel of the inhibitory antibody fragments, from which the 3A2 Fab outperformed others as a specific and potent, low nanomolar range, inhibitor of MT1-MMP. Here, we report the in-depth characterization of the 3A2 antibody. Our multiple in vitro and cell-based tests and assays, and extensive structural modeling of the antibody/protease interactions suggest that the antibody epitope involves the residues proximal to the protease catalytic site and that, in contrast with tissue inhibitor-2 of MMPs (TIMP-2), the 3A2 Fab inactivates the protease functionality by binding to the catalytic domain outside the active site cavity. In agreement with the studies in metastasis by others, our animal studies in acute pulmonary melanoma metastasis support a key role of MT1-MMP in metastatic process. Conversely, the selective anti-MT1-MMP monotherapy significantly alleviated melanoma metastatic burden. It is likely that further affinity maturation of the 3A2 Fab will result in the lead inhibitor and a proof-of-concept for MT1-MMP targeting in metastatic cancers.

High-Throughput Multiplexed Peptide-Centric Profiling Illustrates Both Substrate Cleavage Redundancy and Specificity in the MMP Family.

Matrix metalloproteinases (MMPs) play incompletely understood roles in health and disease. Knowing the MMP cleavage preferences is essential for a better understanding of the MMP functions and design of selective inhibitors. To elucidate the cleavage preferences of MMPs, we employed a high-throughput multiplexed peptide-centric profiling technology involving the cleavage of 18,583 peptides by 18 proteinases from the main sub-groups of the MMP family. Our results enabled comparison of the MMP substrates on a global scale, leading to the most efficient and selective substrates. The data validated the accuracy of our cleavage prediction software. This software allows us and others to locate, with nearly 100% accuracy, the MMP cleavage sites in the peptide sequences. In addition to increasing our understanding of both the selectivity and the redundancy of the MMP family, our study generated a roadmap for the subsequent MMP structural-functional studies and efficient substrate and inhibitor design.

Caspase cleavage sites in the human proteome: CaspDB, a database of predicted substrates.

Caspases are enzymes belonging to a conserved family of cysteine-dependent aspartic-specific proteases that are involved in vital cellular processes and play a prominent role in apoptosis and inflammation. Determining all relevant protein substrates of caspases remains a challenging task. Over 1500 caspase substrates have been discovered in the human proteome according to published data and new substrates are discovered on a daily basis. To aid the discovery process we developed a caspase cleavage prediction method using the recently published curated MerCASBA database of experimentally determined caspase substrates and a Random Forest classification method. On both internal and external test sets, the ranking of predicted cleavage positions is superior to all previously developed prediction methods. The in silico predicted caspase cleavage positions in human proteins are available from a relational database: CaspDB. Our database provides information about potential cleavage sites in a verified set of all human proteins collected in Uniprot and their orthologs, allowing for tracing of cleavage motif conservation. It also provides information about the positions of disease-annotated single nucleotide polymorphisms, and posttranslational modifications that may modulate the caspase cleaving efficiency.

Basis for substrate recognition and distinction by matrix metalloproteinases.

Genomic sequencing and structural genomics produced a vast amount of sequence and structural data, creating an opportunity for structure-function analysis in silico [Radivojac P, et al. (2013) Nat Methods 10(3):221-227]. Unfortunately, only a few large experimental datasets exist to serve as benchmarks for function-related predictions. Furthermore, currently there are no reliable means to predict the extent of functional similarity among proteins. Here, we quantify structure-function relationships among three phylogenetic branches of the matrix metalloproteinase (MMP) family by comparing their cleavage efficiencies toward an extended set of phage peptide substrates that were selected from ∼ 64 million peptide sequences (i.e., a large unbiased representation of substrate space). The observed second-order rate constants [k(obs)] across the substrate space provide a distance measure of functional similarity among the MMPs. These functional distances directly correlate with MMP phylogenetic distance. There is also a remarkable and near-perfect correlation between the MMP substrate preference and sequence identity of 50-57 discontinuous residues surrounding the catalytic groove. We conclude that these residues represent the specificity-determining positions (SDPs) that allowed for the expansion of MMP proteolytic function during evolution. A transmutation of only a few selected SDPs proximal to the bound substrate peptide, and contributing the most to selectivity among the MMPs, is sufficient to enact a global change in the substrate preference of one MMP to that of another, indicating the potential for the rational and focused redesign of cleavage specificity in MMPs.

S-nitrosylation-mediated redox transcriptional switch modulates neurogenesis and neuronal cell death.

Redox-mediated posttranslational modifications represent a molecular switch that controls major mechanisms of cell function. Nitric oxide (NO) can mediate redox reactions via S-nitrosylation, representing transfer of an NO group to a critical protein thiol. NO is known to modulate neurogenesis and neuronal survival in various brain regions in disparate neurodegenerative conditions. However, a unifying molecular mechanism linking these phenomena remains unknown. Here, we report that S-nitrosylation of myocyte enhancer factor 2 (MEF2) transcription factors acts as a redox switch to inhibit both neurogenesis and neuronal survival. Structure-based analysis reveals that MEF2 dimerization creates a pocket, facilitating S-nitrosylation at an evolutionally conserved cysteine residue in the DNA binding domain. S-Nitrosylation disrupts MEF2-DNA binding and transcriptional activity, leading to impaired neurogenesis and survival in vitro and in vivo. Our data define a molecular switch whereby redox-mediated posttranslational modification controls both neurogenesis and neurodegeneration via a single transcriptional signaling cascade.

Sequence-derived structural features driving proteolytic processing.

Proteolytic signaling, or regulated proteolysis, is an essential part of many important pathways such as Notch, Wnt, and Hedgehog. How the structure of the cleaved substrate regions influences the efficacy of proteolytic processing remains underexplored. Here, we analyzed the relative importance in proteolysis of various structural features derived from substrate sequences using a dataset of more than 5000 experimentally verified proteolytic events captured in CutDB. Accessibility to the solvent was recognized as an essential property of a proteolytically processed polypeptide chain. Proteolytic events were found nearly uniformly distributed among three types of secondary structure, although with some enrichment in loops. Cleavages in α-helices were found to be relatively abundant in regions apparently prone to unfolding, while cleavages in ß-structures tended to be located at the periphery of ß-sheets. Application of the same statistical procedures to proteolytic events divided into separate sets according to the catalytic classes of proteases proved consistency of the results and confirmed that the structural mechanisms of proteolysis are universal. The estimated prediction power of sequence-derived structural features, which turned out to be sufficiently high, presents a rationale for their use in bioinformatic prediction of proteolytic events.

Substrate cleavage profiling suggests a distinct function of Bacteroides fragilis metalloproteinases (fragilysin and metalloproteinase II) at the microbiome-inflammation-cancer interface.

Enterotoxigenic anaerobic Bacteroides fragilis is a significant source of inflammatory diarrheal disease and a risk factor for colorectal cancer. Two distinct metalloproteinase types (the homologous 1, 2, and 3 isoforms of fragilysin (FRA1, FRA2, and FRA3, respectively) and metalloproteinase II (MPII)) are encoded by the B. fragilis pathogenicity island. FRA was demonstrated to be important to pathogenesis, whereas MPII, also a potential virulence protein, remained completely uncharacterized. Here, we, for the first time, extensively characterized MPII in comparison with FRA3, a representative of the FRA isoforms. We employed a series of multiplexed peptide cleavage assays to determine substrate specificity and proteolytic characteristics of MPII and FRA. These results enabled implementation of an efficient assay of MPII activity using a fluorescence-quenched peptide and contributed to structural evidence for the distinct substrate cleavage preferences of MPII and FRA. Our data imply that MPII specificity mimics the dibasic Arg↓Arg cleavage motif of furin-like proprotein convertases, whereas the cleavage motif of FRA (Pro-X-X-Leu-(Arg/Ala/Leu)↓) resembles that of human matrix metalloproteinases. To the best of our knowledge, MPII is the first zinc metalloproteinase with the dibasic cleavage preferences, suggesting a high level of versatility of metalloproteinase proteolysis. Based on these data, we now suggest that the combined (rather than individual) activity of MPII and FRA is required for the overall B. fragilis virulence in vivo.

New details of HCV NS3/4A proteinase functionality revealed by a high-throughput cleavage assay.

BACKGROUND: The hepatitis C virus (HCV) genome encodes a long polyprotein, which is processed by host cell and viral proteases to the individual structural and non-structural (NS) proteins. HCV NS3/4A serine proteinase (NS3/4A) is a non-covalent heterodimer of the N-terminal, ∼180-residue portion of the 631-residue NS3 protein with the NS4A co-factor. NS3/4A cleaves the polyprotein sequence at four specific regions. NS3/4A is essential for viral replication and has been considered an attractive drug target. METHODOLOGY/PRINCIPAL FINDINGS: Using a novel multiplex cleavage assay and over 2,660 peptide sequences derived from the polyprotein and from introducing mutations into the known NS3/4A cleavage sites, we obtained the first detailed fingerprint of NS3/4A cleavage preferences. Our data identified structural requirements illuminating the importance of both the short-range (P1-P1') and long-range (P6-P5) interactions in defining the NS3/4A substrate cleavage specificity. A newly observed feature of NS3/4A was a high frequency of either Asp or Glu at both P5 and P6 positions in a subset of the most efficient NS3/4A substrates. In turn, aberrations of this negatively charged sequence such as an insertion of a positively charged or hydrophobic residue between the negatively charged residues resulted in inefficient substrates. Because NS5B misincorporates bases at a high rate, HCV constantly mutates as it replicates. Our analysis revealed that mutations do not interfere with polyprotein processing in over 5,000 HCV isolates indicating a pivotal role of NS3/4A proteolysis in the virus life cycle. CONCLUSIONS/SIGNIFICANCE: Our multiplex assay technology in light of the growing appreciation of the role of proteolytic processes in human health and disease will likely have widespread applications in the proteolysis research field and provide new therapeutic opportunities.

Synthesis, Physicochemical Studies, Molecular Dynamics Simulations, and Metal-Ion-Dependent Antiproliferative and Antiangiogenic Properties of Cone ICL670-Substituted Calix[4]arenes.

Iron chelators, through their capacity to modulate the iron concentration in cells, are promising molecules for cancer chemotherapy. Chelators with high lipophilicity easily enter into cells and deplete the iron intracellular pool. Consequently, iron-dependent enzymes, such as ribonucleotide reductase, which is over-expressed in cancer cells, become nonfunctional. A series of calix[4]arene derivatives substituted at the lower rim by ICL670, a strong FeIII chelator, have been synthesized. Physicochemical properties and antiproliferative, angiogenesis, and tumorigenesis effects of two calix[4]arenes mono- (5a) or disubstituted (5b) with ICL670 have been studied. These compounds form metal complexes in a ratio of one to two ligands per FeIII atom as shown by combined analyses of the protometric titration curves and ESIMS spectra. The grafting of an ICL670 group on a calix[4]arene core does not significantly alter the acid-base properties, but improves the iron-chelating and lipophilicity properties. The best antiproliferative and anti-angiogenic results were obtained with calix[4]arene ligand 5a, which possesses the highest corresponding properties. Analyses of molecular dynamics simulations performed on the two calix[4]arenes provide three-dimensional structures of the complexes and proved 5a to be the most stable upon complexation.

Matrix metalloproteinase proteolysis of the mycobacterial HSP65 protein as a potential source of immunogenic peptides in human tuberculosis.

Mycobacterium tuberculosis is the causative agent of human tuberculosis (TB). Mycobacterial secretory protein ESAT-6 induces matrix metalloproteinase (MMP)-9 in epithelial cells neighboring infected macrophages. MMP-9 then enhances recruitment of uninfected macrophages, which contribute to nascent granuloma maturation and bacterial growth. Disruption of MMP-9 function attenuates granuloma formation and bacterial growth. The abundant mycobacterial 65 kDa heat shock protein (HSP65) chaperone is the major target for the immune response and a critical component in M. tuberculosis adhesion to macrophages. We hypothesized that HSP65 is susceptible to MMP-9 proteolysis and that the resulting HSP65 immunogenic peptides affect host adaptive immunity. To identify MMPs that cleave HSP65, we used MMP-2 and MMP-9 gelatinases, the simple hemopexin domain MMP-8, membrane-associated MMP-14, MMP-15, MMP-16 and MMP-24, and glycosylphosphatidylinositol-linked MMP-17 and MMP-25. We determined both the relative cleavage efficiency of MMPs against the HSP65 substrate and the peptide sequence of the cleavage sites. Cleavage of the unstructured PAGHG474L C-terminal region initiates the degradation of HSP65 by MMPs. This initial cleavage destroys the substrate-binding capacity of the HSP65 chaperone. Multiple additional cleavages of the unfolded HSP65 then follow. MMP-2, MMP-8, MMP-14, MMP-15 and MMP-16, in addition to MMP-9, generate the known highly immunogenic N-terminal peptide of HSP65. Based on our biochemical data, we now suspect that MMP proteolysis of HSP65 in vivo, including MMP-9 proteolysis, also results in the abundant generation of the N-terminal immunogenic peptide and that this peptide, in addition to intact HSP65, contributes to the complex immunomodulatory interplay in the course of TB infection.

The Wnt/planar cell polarity protein-tyrosine kinase-7 (PTK7) is a highly efficient proteolytic target of membrane type-1 matrix metalloproteinase: implications in cancer and embryogenesis.

PTK7 is an essential component of the Wnt/planar cell polarity (PCP) pathway. We provide evidence that the Wnt/PCP pathway converges with pericellular proteolysis in both normal development and cancer. Here, we demonstrate that membrane type-1 matrix metalloproteinase (MT1-MMP), a key proinvasive proteinase, functions as a principal sheddase of PTK7. MT1-MMP directly cleaves the exposed PKP(621)↓LI sequence of the seventh Ig-like domain of the full-length membrane PTK7 and generates, as a result, an N-terminal, soluble PTK7 fragment (sPTK7). The enforced expression of membrane PTK7 in cancer cells leads to the actin cytoskeleton reorganization and the inhibition of cell invasion. MT1-MMP silencing and the analysis of the uncleavable L622D PTK7 mutant confirm the significance of MT1-MMP proteolysis of PTK7 in cell functions. Our data also demonstrate that a fine balance between the metalloproteinase activity and PTK7 levels is required for normal development of zebrafish (Danio rerio). Aberration of this balance by the proteinase inhibition or PTK7 silencing results in the PCP-dependent convergent extension defects in the zebrafish. Overall, our data suggest that the MT1-MMP-PTK7 axis plays an important role in both cancer cell invasion and normal embryogenesis in vertebrates. Further insight into these novel mechanisms may promote understanding of directional cell motility and lead to the identification of therapeutics to treat PCP-related developmental disorders and malignancy.