Caulfield, Thomas R.

Matrix metalloproteinases (MMPs) are a family of zinc‐dependent endopeptidases that function to mediate tissue remodeling; however, they are also implicated in pathological conditions such as cancer, inflammatory and degenerative diseases. There are 23 human MMPs and significant sequence and structural homology exists among them and with other metalloproteinases. This represents a substantial challenge in developing therapeutically viable selective inhibitors due to off‐target effects. Tissue inhibitors of metalloproteinases (TIMPs) are a family of four endogenous protein protease inhibitors that regulate MMP activity. Our lab and others have utilized TIMPs as scaffolds, engineering them to selectively target specific MMPs for potential therapeutic applications. We previously employed a yeast surface display platform to identify TIMP‐1 variants possessing enhanced MMP‐3 affinity. In the current study, we used structures of MMP‐3 co‐crystallized with either WT TIMP‐1 or an engineered TIMP‐1 variant as seeds for molecular dynamics simulations. Simulations were also conducted on unbound WT and engineered TIMP‐1 to examine MMP‐independent dynamics that distinguish the engineered variant. The simulation data were analyzed with root‐mean‐square deviation/fluctuation, time‐dependent secondary structure, temporally resolved hydrogen bonding, conformational clustering, and dynamical network algorithms to discern structural and dynamical determinants of the increased affinity. Key observations that demonstrate the increased effectiveness of engineered TIMP‐1 as an MMP‐3 inhibitor were secondary structural changes in the MMP‐binding loops in proximity to mutations, attenuated global and mutation‐adjacent dynamics, decreased conformational dynamics, and increased persistence of MMP‐3‐compatible conformations. There were also indications of greater interdomain cooperativity within the engineered TIMP‐1 variant: increased interdomain hydrogen bonding, coordinated motions of the N‐terminal region with the C‐terminal domain, and increased intercommunity information flow/motion transfer. From these observations, we hypothesize the increased interdomain interactions represent a mechanism of dynamical stabilization such that mainly cooperative motion between domains presides, rather than random independent motions. Our findings suggest propagated motions between domains may influence the MMP affinity of TIMPs. Furthermore, these data suggest future strategies to design next‐generation mutant TIMP libraries by targeting regions dynamically coupled to the MMP binding loops for variation, thus exploiting natural protein dynamics to achieve desired selectivity.

Leishmaniasis remains a serious public health problem in many tropical regions of the world. Among neglected tropical diseases, the mortality rate of leishmaniasis is second only to malaria. All currently approved therapeutics have toxic side effects and face rapidly increasing resistance. To identify existing drugs with antileishmanial activity and predict the mechanism of action, we designed a drug-discovery pipeline utilizing both in-silico and in-vitro methods. First, we screened compounds from the Selleckchem Bio-Active Compound Library containing ~1,622 FDA-approved drugs and narrowed these down to 96 candidates based on data mining for possible anti-parasitic properties. Next, we completed preliminary in-vitro testing of compounds against Leishmania amastigotes and selected the most promising active compounds, Lansoprazole and Posaconazole. We identified possible Leishmania drug targets of Lansoprazole and Posaconazole using several available servers. Our in-silico screen identified likely Lansoprazole targets as the closely related calcium-transporting ATPases (LdBPK_352080.1, LdBPK_040010.1, and LdBPK_170660.1), and the Posaconazole target as lanosterol 14-alpha-demethylase (LdBPK_111100.1). Further validation showed LdBPK_352080.1 to be the most plausible target based on induced-fit docking followed by long (100ns) MD simulations to confirm the stability of the docked complexes. We present a likely ion channel-based mechanism of action of Lansoprazole against Leishmania calcium-transporting ATPases, which are essential for parasite metabolism and infectivity. The LdBPK_111100.1 interaction with Posaconazole is very similar to the known fungal orthologue. Herein, we present two novel anti-leishmanial agents, Posaconazole and Lansoprazole, already approved by the FDA for different indications and propose plausible mechanisms of action for their antileishmanial activity.

A rare APOE variant APOE3 -Jac reduces self-aggregation and promotes healthy brain aging.

PLEXIND1 is upregulated in several cancers, including pancreatic ductal adenocarcinoma (PDAC). It is an established mediator of semaphorin signaling, and neuropilins are its known coreceptors. Herein, we report data to support the proposal that PLEXIND1 acts as a transforming growth factor beta (TGFβ) coreceptor, modulating cell growth through SMAD3 signaling. Our findings demonstrate that PLEXIND1 plays a pro-tumorigenic role in PDAC cells with oncogenic KRAS (KRASmut). We show in KRASmut PDAC cell lines (PANC-1, AsPC-1,4535) PLEXIND1 downregulation results in decreased cell viability (in vitro) and reduced tumor growth (in vivo). Conversely, PLEXIND1 acts as a tumor suppressor in the PDAC cell line (BxPC-3) with wild-type KRAS (KRASwt), as its reduced expression results in higher cell viability (in-vitro) and tumor growth (in vivo). Additionally, we demonstrate that PLEXIND1-mediated interactions can be selectively disrupted using a peptide based on its C-terminal sequence (a PDZ domain-binding motif), an outcome that may possess significant therapeutic implications. To our knowledge, this is the first report showing that (1) PLEXIND1 acts as a TGFβ coreceptor and mediates SMAD3 signaling, and (2) differential roles of PLEXIND1 in PDAC cell lines correlate with KRASmut and KRASwt status.

The mitochondrial ubiquitin (Ub) kinase PINK1 and the cytosolic E3 Ub ligase PRKN (Parkin) together mediate the selective degradation of damaged mitochondria via the autophagy‐lysosome system (mitophagy). Complete loss of either gene function results in failure of this mitochondrial quality control and thus accumulation of damaged organelles. While homozygous or compound heterozygous mutations in PINK1 or PRKN are the most common causes of early‐onset Parkinson's disease (PD), the contribution of heterozygous variants to PD is debated. However, recently we identified a partial dominant‐negative effect of a heterozygous PINK1 missense mutations that is associated with a significant increased risk and an earlier onset compared to late‐onset, sporadic PD. This underscores that not all missense mutations are alike and suggests critical threshold functions of this stress‐activated, neuroprotective pathway. Here, we analyzed the underlying mechanisms in greater detail and found that the pathogenic PINK1 mutant, which can dimerize with wild‐type PINK1 protein, is aberrantly auto‐phosphorylated, impacting its enzymatic activity towards the substrates Ub and PRKN. Phosphorylation of both PINK1 targets however, is critical for PRKN recruitment as well as efficient labeling of mitochondria with phosphorylated poly‐Ub chains that serve as the mitophagy tag. In addition to cells from PD patients, we confirmed the altered phosphorylation pattern of mutant PINK1 under endogenous conditions in isogenic CRISPR/Cas9 genome engineered cells. Excitingly, while studying the pathogenic mechanisms in detail with a structure‐function approach, we identified a novel PINK1 variant that exerts even greater kinase activity than wild type PINK1. Our results further suggest that this hyperactive form of PINK1 likely modulates substrate receptivity thereby resulting in enhanced kinase activity and greater mitophagy rates. Going forward, it will be interesting to study if enhanced PINK1 signaling translates into greater neuroprotection. An in‐depth understanding of the regulation and enzyme functions of PINK1 may provide a rationale for a structure‐function based drug design in the future. Support or Funding Information This work was partially supported by funding from the Mayo Clinic Foundation, the National Institutes of Health (NIH)/National Institute of Neurological Disorders and Stroke (NINDS) [R01#NS085070] and the Michael J. Fox Foundation for Parkinson's Research to W.S. This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal .

Serine protease inhibitors of the Kunitz‐bovine pancreatic trypsin inhibitor (BPTI) family, such as tissue factor pathway inhibitor‐1 Kunitz domain 1 ( KD1 TFPI1) and amyloid precursor protein inhibitor (APPI), are biological regulators of proteolysis and serve as important scaffolds for developing therapeutics. Mesotrypsin is a unique human protease with cancer‐promoting properties which rapidly cleaves and inactivates serine protease inhibitors; because of this, it is a challenging therapeutic target and an excellent model for developing therapeutic inhibitors with improved proteolytic resistance. We hypothesize that stabilization of the Kunitz inhibitor binding loop by disulfide bond anchoring to the scaffold can reduce proteolytic susceptibility and improve binding efficacy. Here we have engineered a new disulfide bond into both KD1 TFPI1 and an APPI variant previously engineered for mesotrypsin affinity. The new disulfide bond anchors the binding loop to the scaffold, resulting in reduced mesotrypsin proteolytic cleavage by 17‐fold and 3‐fold, respectively. To gain insight into the structural basis for improved proteolytic stability, we crystallized mutant KD1 TFPI1 and APPI with mesotrypsin and solved these novel structures. We then performed all‐atom molecular dynamics simulations using these mesotrypsin‐mutant inhibitor complexes compared to models lacking the engineered disulfide bond. Our simulation results support our hypothesis that the cause of decreased proteolysis rates is dampened conformational fluctuations in the binding loop and minimized motion at the enzyme‐inhibitor interface. We also performed simulations using unbound inhibitors to show that stabilization of the binding loop to states similar to the bound configuration correlates with stronger binding affinity. Our disulfide engineering approach is potentially a useful tool in biopharmaceutical development for mesotrypsin, other serine proteases, and possibly for other enzyme targets. Support or Funding Information NIH grant R01CA154387 to E.S.R., European Research Council grant 336041 to N.P., US‐Israel Binational Science Foundation grant BSF2015134 to E.S.R. and N.P. Xray diffraction data were measured at the Berkeley Center for Structural Biology at the Advanced Light Source and at the European Synchrotron Radiation Facility. This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal .

The prohormone convertases PC1/3 and PC2 are eukaryotic serine proteases involved in proteolytic maturation of peptide hormone precursors and implicated in a variety of pathological conditions, including obesity, diabetes, and neurodegenerative diseases. In this work, we screened 45 compounds obtained by derivatization of a 2,5‐dideoxystreptamine scaffold. We found four promising PC1/3 competitive inhibitors and three PC2 inhibitors exhibiting various inhibition mechanisms with sub‐ and low micromolar inhibitory potency against a fluorogenic substrate. Low micromolar concentrations of certain compounds blocked the processing of the physiological substrate proglucagon. Interestingly, we identified compounds that were able to stimulate both 87 kDa PC1/3 and PC2 activity. Molecular modeling studies revealed interactions of the PC1/3 inhibitors with the active site that suggest structural modifications to further enhance potency. In support of experimental data suggesting that PC2 inhibition likely occurs via an allosteric mechanism, we identified several possible allosteric binding sites using computational searches. This work was supported by DA050854 to Iris Lindberg.

There is evidence in the literature that many serine proteases, as typified by thrombin, possess an autoinhibitory loop (AIL) that controls substrate access to the active site. The presence of an AIL has not been documented in trypsins. Mesotrypsin is a human trypsin implicated in tumor cell invasiveness for pancreatic and other cancers. However, it is a challenging therapeutic target because of high structural similarity to other trypsin‐like serine proteases, typically lower affinity for inhibitors compared to related serine proteases, and unique ability to cleave proteinaceous serine protease inhibitors. We have recently solved a high‐resolution crystal structure of unbound mesotrypsin which revealed that residues 215‐220 partially occluded the substrate binding site. Comparison of this new structure with previous inhibitor‐bound crystal structures of mesotrypsin demonstrate the conformational transition of the AIL between closed and open states. These two conformational states were investigated with molecular dynamics simulations to discern the dynamical behavior of the AIL in different conformations. The AIL has increased local dynamics in the closed position, as well as a larger number of conformational clusters; this suggests there is a more well‐defined open state ready to accept substrate, while there are a dynamic variety of closed states that preclude substrate binding to various extents. Additionally, we identified novel binding sites that are exposed by the conformational shift of the AIL to the closed state. We have conducted initial in silico screening for compounds predicted to bind preferentially to the AIL closed (but not open) structure and stabilize the closed conformation. These efforts have resulted in more than 40 compounds as initial hits. Exploiting the AIL as a novel allosteric binding site represents a new mechanism for targeting mesotrypsin and preventing its protumorigenic function.

COVID-19 is a devastating respiratory and inflammatory illness caused by a new coronavirus that is rapidly spreading throughout the human population. Over the past 12 months, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19, has already infected over 160 million (>20% located in United States) and killed more than 3.3 million people around the world (>20% deaths in USA). As we face one of the most challenging times in our recent history, there is an urgent need to identify drug candidates that can attack SARS-CoV-2 on multiple fronts. We have therefore initiated a computational dynamics drug pipeline using molecular modeling, structure simulation, docking and machine learning models to predict the inhibitory activity of several million compounds against two essential SARS-CoV-2 viral proteins and their host protein interactors—S/Ace2, Tmprss2, Cathepsins L and K, and Mpro—to prevent binding, membrane fusion and replication of the virus, respectively. All together, we generated an ensemble of structural conformations that increase high-quality docking outcomes to screen over >6 million compounds including all FDA-approved drugs, drugs under clinical trial (>3000) and an additional >30 million selected chemotypes from fragment libraries. Our results yielded an initial set of 350 high-value compounds from both new and FDA-approved compounds that can now be tested experimentally in appropriate biological model systems. We anticipate that our results will initiate screening campaigns and accelerate the discovery of COVID-19 treatments.