Statistical Mechanics Metrics in Pairing and Parsing In Silico and Phenotypic Data of a Novel Genetic NFκB1 (c.T638A) Variant

(1) Background: Mutations in NFκB1, a transcriptional regulator of immunomodulating proteins, are a known cause of inborn errors of immunity. Our proband is a 22-year-old male with a diagnosis of common variable immunodeficiency (CVID), cytopenias with massive splenomegaly, and nodular regenerative hyperplasia of the liver. Genetic studies identified a novel, single-point mutation variant in NFκB1, c. T638A p. V213E. (2) Methods: Next-generation panel sequencing of the patient uncovered a novel single-point mutation in the NFκB1 gene that was modeled using the I-TASSER homology-modeling software, and molecular dynamics were assessed using the YASARA2 software (version 20.14.24). (3) Results: This variant replaces valine with glutamic acid at position 213 in the NFκB1 sequence. Molecular modeling and molecular dynamic studies showed altered dynamics in and around the rel homology domain, ankyrin regions, and death domain of the protein. We postulate that these changes alter overall protein function. (4) Conclusions: This case suggests the pathogenicity of a novel variant using protein-modeling techniques and molecular dynamic simulations.

Protein structure-based in-silico approaches to drug discovery: Guide to COVID-19 therapeutics

With more than 5 million fatalities and close to 300 million reported cases, COVID-19 is the first documented pandemic due to a coronavirus that continues to be a major health challenge. Despite being rapid, uncontrollable, and highly infectious in its spread, it also created incentives for technology development and redefined public health needs and research agendas to fast-track innovations to be translated. Breakthroughs in computational biology peaked during the pandemic with renewed attention to making all cutting-edge technology deliver agents to combat the disease. The demand to develop effective treatments yielded surprising collaborations from previously segregated fields of science and technology. The long-standing pharmaceutical industry's aversion to repurposing existing drugs due to a lack of exponential financial gain was overrun by the health crisis and pressures created by front-line researchers and providers. Effective vaccine development even at an unprecedented pace took more than a year to develop and commence trials. Now the emergence of variants and waning protections during the booster shots is resulting in breakthrough infections that continue to strain health care systems. As of now, every protein of SARS-CoV-2 has been structurally characterized and related host pathways have been extensively mapped out. The research community has addressed the druggability of a multitude of possible targets. This has been made possible due to existing technology for virtual computer-assisted drug development as well as new tools and technologies such as artificial intelligence to deliver new leads. Here in this article, we are discussing advances in the drug discovery field related to target-based drug discovery and exploring the implications of known target-specific agents on COVID-19 therapeutic management. The current scenario calls for more personalized medicine efforts and stratifying patient populations early on for their need for different combinations of prognosis-specific therapeutics. We intend to highlight target hotspots and their potential agents, with the ultimate goal of using rational design of new therapeutics to not only end this pandemic but also uncover a generalizable platform for use in future pandemics.

Unveiling an indole alkaloid diketopiperazine biosynthetic pathway that features a unique stereoisomerase and multifunctional methyltransferase

The 2,5-diketopiperazines are a prominent class of bioactive molecules. The nocardioazines are actinomycete natural products that feature a pyrroloindoline diketopiperazine scaffold composed of two D-tryptophan residues functionalized by N- and C-methylation, prenylation, and diannulation. Here we identify and characterize the nocardioazine B biosynthetic pathway from marine Nocardiopsis sp. CMB-M0232 by using heterologous biotransformations, in vitro biochemical assays, and macromolecular modeling. Assembly of the cyclo-L-Trp-L-Trp diketopiperazine precursor is catalyzed by a cyclodipeptide synthase. A separate genomic locus encodes tailoring of this precursor and includes an aspartate/glutamate racemase homolog as an unusual D/L isomerase acting upon diketopiperazine substrates, a phytoene synthase-like prenyltransferase as the catalyst of indole alkaloid diketopiperazine prenylation, and a rare dual function methyltransferase as the catalyst of both N- and C-methylation as the final steps of nocardioazine B biosynthesis. The biosynthetic paradigms revealed herein showcase Nature's molecular ingenuity and lay the foundation for diketopiperazine diversification via biocatalytic approaches.

Mycolactone: A Broad Spectrum Multitarget Antiviral Active in the Picomolar Range for COVID-19 Prevention and Cure

We have previously shown computationally that Mycolactone (MLN), a toxin produced by Mycobacterium ulcerans, strongly binds to Munc18b and other proteins, presumably blocking degranulation and exocytosis of blood platelets and mast cells. We investigated the effect of MLN on endocytosis using similar approaches, and it bound strongly to the N-terminal of the clathrin protein and a novel SARS-CoV-2 fusion protein. Experimentally, we found 100% inhibition up to 60 nM and 84% average inhibition at 30 nM in SARS-CoV-2 live viral assays. MLN was also 10× more potent than remdesivir and molnupiravir. MLN's toxicity against human alveolar cell line A549, immortalized human fetal renal cell line HEK293, and human hepatoma cell line Huh7.1 were 17.12%, 40.30%, and 36.25%, respectively. The cytotoxicity IC₅₀ breakpoint ratio versus anti-SARS-CoV-2 activity was more than 65-fold. The IC₅₀ values against the alpha, delta, and Omicron variants were all below 0.020 µM, and 134.6 nM of MLN had 100% inhibition in an entry and spread assays. MLN is eclectic in its actions through its binding to Sec61, AT2R, and the novel fusion protein, making it a good drug candidate for treating and preventing COVID-19 and other similarly transmitted enveloped viruses and pathogens.

Utilizing genetic code expansion to modify N-TIMP2 specificity towards MMP-2, MMP-9, and MMP-14

Matrix metalloproteinases (MMPs) regulate the degradation of extracellular matrix (ECM) components in biological processes. MMP activity is controlled by natural tissue inhibitors of metalloproteinases (TIMPs) that non-selectively inhibit the function of multiple MMPs via interaction with the MMPs' Zn2+-containing catalytic pocket. Recent studies suggest that TIMPs engineered to confer MMP specificity could be exploited for therapeutic purposes, but obtaining specific TIMP-2 inhibitors has proved to be challenging. Here, in an effort to improve MMP specificity, we incorporated the metal-binding non-canonical amino acids (NCAAs), 3,4-dihydroxyphenylalanine (L-DOPA) and (8-hydroxyquinolin-3-yl)alanine (HqAla), into the MMP-inhibitory N-terminal domain of TIMP2 (N-TIMP2) at selected positions that interact with the catalytic Zn2+ ion (S2, S69, A70, L100) or with a structural Ca2+ ion (Y36). Evaluation of the inhibitory potency of the NCAA-containing variants towards MMP-2, MMP-9 and MMP-14 in vitro revealed that most showed a significant loss of inhibitory activity towards MMP-14, but not towards MMP-2 and MMP-9, resulting in increased specificity towards the latter proteases. Substitutions at S69 conferred the best improvement in selectivity for both L-DOPA and HqAla variants. Molecular modeling provided an indication of how MMP-2 and MMP-9 are better able to accommodate the bulky NCAA substituents at the intermolecular interface with N-TIMP2. The models also showed that, rather than coordinating to Zn2+, the NCAA side chains formed stabilizing polar interactions at the intermolecular interface with MMP-2 and MMP-9. Our findings illustrate how incorporation of NCAAs can be used to probe-and possibly exploit-differential tolerance for substitution within closely related protein-protein complexes as a means to improve specificity.

Utilizing genetic code expansion to modify N-TIMP2 specificity towards MMP-2, MMP-9, and MMP-14

Matrix metalloproteinases (MMPs) regulate the degradation of extracellular matrix (ECM) components in biological processes. MMP activity is controlled by natural tissue inhibitors of metalloproteinases (TIMPs) that non-selectively inhibit the function of multiple MMPs via interaction with the MMPs' Zn ²⁺ -containing catalytic pocket. Recent studies suggest that TIMPs engineered to confer MMP specificity could be exploited for therapeutic purposes, but obtaining specific TIMP-2 inhibitors has proved to be challenging. Here, in an effort to improve MMP specificity, we incorporated the metal-binding non-canonical amino acids (NCAAs), 3,4-dihydroxyphenylalanine (L-DOPA) and (8-hydroxyquinolin-3-yl)alanine (HqAla), into the MMP-inhibitory N-terminal domain of TIMP2 (N-TIMP2) at selected positions that interact with the catalytic Zn ²⁺ ion (S2, S69, A70, L100) or with a structural Ca ²⁺ ion (Y36). Evaluation of the inhibitory potency of the NCAA-containing variants towards MMP-2, MMP-9 and MMP-14 in vitro revealed that most showed a significant loss of inhibitory activity towards MMP-14, but not towards MMP-2 and MMP-9, resulting in increased specificity towards the latter proteases. Substitutions at S69 conferred the best improvement in selectivity for both L-DOPA and HqAla variants. Molecular modeling revealed how MMP-2 and MMP-9 are better able to accommodate the bulky NCAA substituents at the intermolecular interface with N-TIMP2. The models also showed that, rather than coordinating to Zn ²⁺ , the NCAA side chains formed stabilizing polar interactions at the intermolecular interface with MMP-2 and MMP-9. The findings illustrate how incorporation of NCAAs can be used to probe and exploit differential tolerance for substitution within closely related protein-protein complexes to achieve improved specificity.

Substitution of PINK1 Gly411 modulates substrate receptivity and turnover

The ubiquitin (Ub) kinase-ligase pair PINK1-PRKN mediates the degradation of damaged mitochondria by macroautophagy/autophagy (mitophagy). PINK1 surveils mitochondria and upon stress accumulates on the mitochondrial surface where it phosphorylates serine 65 of Ub to activate PRKN and to drive mitochondrial turnover. While loss of either PINK1 or PRKN is genetically linked to Parkinson disease (PD) and activating the pathway seems to have great therapeutic potential, there is no formal proof that stimulation of mitophagy is always beneficial. Here we used biochemical and cell biological methods to study single nucleotide variants in the activation loop of PINK1 to modulate the enzymatic function of this kinase. Structural modeling and in vitro kinase assays were used to investigate the molecular mechanism of the PINK1 variants. In contrast to the PD-linked PINK1^G411S mutation that diminishes Ub kinase activity, we found that the PINK1^G411A variant significantly boosted Ub phosphorylation beyond levels of PINK1 wild type. This resulted in augmented PRKN activation, mitophagy rates and increased viability after mitochondrial stress in midbrain-derived, gene-edited neurons. Mechanistically, the G411A variant stabilizes the kinase fold of PINK1 and transforms Ub to adopt the preferred, C-terminally retracted conformation for improved substrate turnover. In summary, we identify a critical role of residue 411 for substrate receptivity that may now be exploited for drug discovery to increase the enzymatic function of PINK1. The genetic substitution of Gly411 to Ala increases mitophagy and may be useful to confirm neuroprotection in vivo and might serve as a critical positive control during therapeutic development.Abbreviations: ATP: adenosine triphosphate; CCCP: carbonyl cyanide m-chlorophenyl hydrazone; Ub-CR: ubiquitin with C-terminally retracted tail; CTD: C-terminal domain (of PINK1); ELISA: enzyme-linked immunosorbent assay; HCI: high-content imaging; IB: immunoblot; IF: immunofluorescence; NPC: neuronal precursor cells; MDS: molecular dynamics simulation; PD: Parkinson disease; p-S65-Ub: ubiquitin phosphorylated at Ser65; RMSF: root mean scare fluctuation; TOMM: translocase of outer mitochondrial membrane; TVLN: ubiquitin with T66V and L67N mutation, mimics Ub-CR; Ub: ubiquitin; WT: wild-type.

Phenotypic screening models for rapid diagnosis of genetic variants and discovery of personalized therapeutics

Precision medicine strives for highly individualized treatments for disease under the notion that each individual's unique genetic makeup and environmental exposures imprints upon them not only a disposition to illness, but also an optimal therapeutic approach. In the realm of rare disorders, genetic predisposition is often the predominant mechanism driving disease presentation. For such, mostly, monogenic disorders, a causal gene to phenotype association is likely. As a result, it becomes important to query the patient's genome for the presence of pathogenic variations that are likely to cause the disease. Determining whether a variant is pathogenic or not is critical to these analyses and can be challenging, as many disease-causing variants are novel and, ergo, have no available functional data to help categorize them. This problem is exacerbated by the need for rapid evaluation of pathogenicity, since many genetic diseases present in young children who will experience increased morbidity and mortality without rapid diagnosis and therapeutics. Here, we discuss the utility of animal models, with a focus mainly on C. elegans, as a contrast to tissue culture and in silico approaches, with emphasis on how these systems are used in determining pathogenicity of variants with uncertain significance and then used to screen for novel therapeutics.

Attacking COVID-19 Progression Using Multi-Drug Therapy for Synergetic Target Engagement

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.

Endotheliitis, endothelin, and endothelin receptor blockers in COVID-19

We summarize the role of endothelin as a potent vasoconstrictor, pro-inflammatory, pro-oxidative agent in the pathophysiologic effects and end-organ dysfunction of coronavirus disease 2019 (COVID-19). Endotheliitis is an under-recognized pathophysiologic process that causes various types of dysfunction in end organs, including heart, lung, kidney, and brain. Endothelin receptor blockers, such as bosentan and sitaxentan, can pave a path ahead in the realm of COVID-19 therapies. These agents have a potential role against COVID-19 and should be studied in research trials to determine their efficacy in treatment of this severe disease.

Clinical description & molecular modeling of novel MAX pathogenic variant causing pheochromocytoma in family, supports paternal parent-of-origin effect

The titular member of the MAX network of proteins, MYC-associated factor X (MAX), serves an important regulatory function in transcription of E-box genes associated with cell proliferation, differentiation, and apoptosis. Wild type MAX dimerizes with both MYC and MAD, both of which are members of the MAX network, and can promote or repress cell functions as needed. However, pathogenic variants in MAX are known to upset this balance, leading to uncontrolled oncogenic activity and disease phenotypes such as paragangliomas and pheochromocytomas. We report a 58-year-old male and his 32-year-old daughter, both of which have a history of pheochromocytoma and the unique nonsense MAX variant c.271C>T (p.Q91X). These individuals were diagnosed with pheochromocytomas in their early twenties that were later removed through corrective surgery. The father now presents with recurrent symptoms of hypertension, hyperhidrosis, and headaches, which accompany new pheochromocytomas of his remaining adrenal gland. Pathogenicity of this MAX variant is proven through molecular modeling. The case of this father-daughter pair supports both heritability of pheochromocytoma and the paternal parent-of-origin effect for MAX pathogenic variants.

A Virtual Screening Platform Identifies Chloroethylagelastatin A as a Potential Ribosomal Inhibitor

Chloroethylagelastatin A (CEAA) is an analogue of agelastatin A (AA), a natural alkaloid derived from a marine sponge. It is under development for therapeutic use against brain tumors as it has excellent central nervous system (CNS) penetration and pre-clinical therapeutic activity against brain tumors. Recently, AA was shown to inhibit protein synthesis by binding to the ribosomal A-site. In this study, we developed a novel virtual screening platform to perform a comprehensive screening of various AA analogues showing that AA analogues with proven therapeutic activity including CEAA have significant ribosomal binding capacity whereas therapeutically inactive analogues show poor ribosomal binding and revealing structural fingerprint features essential for drug-ribosome interactions. In particular, CEAA was found to have greater ribosomal binding capacity than AA. Biological tests showed that CEAA binds the ribosome and contributes to protein synthesis inhibition. Our findings suggest that CEAA may possess ribosomal inhibitor activity and that our virtual screening platform may be a useful tool in discovery and development of novel ribosomal inhibitors.

An activating germline IDH1 variant associated with a tumor entity characterized by unilateral and bilateral chondrosarcoma of the mastoid

Chondrogenic tumors involving the temporal bone are rare and typically arise spontaneously with unilateral presentation. Somatic IDH mutations are common in these tumors, but germline inheritance has not been documented to our knowledge. We describe familial chondrosarcoma, grade 1, of the mastoid with unilateral presentation in the mother and bilateral presentation in each of her two children. Each individual presented with headaches, facial paresis, and conductive hearing loss between the ages of 9–12. Exome sequencing of all three affected family members identified a shared germline heterozygous c.299G>A (p.Arg100Gln) missense variant in IDH1. The p.Arg100Gln variant has only rarely been observed as a somatic mutation in glial tumors, and previous in vitro experiments have shown that p.Arg100Gln produces small amounts of the oncometabolite D-2-hydroxyglutarate (D2HG). Biochemical testing in all three affected family members on urine and plasma was unable to detect increases in D2HG in these sample types. Due to insufficient tumor for methylation studies, we performed genome-wide methylation analysis of an IDH1 p.Arg100Gln mutant brain tumor from an unrelated individual to functionally evaluate this variant. These studies demonstrated a global hypermethylation phenotype consistent with other known isocitrate dehydrogenase (IDH) mutant brain tumors, suggesting that this variant has neomorphic activity despite low-level production of D2HG. The bones of the facial skeleton are formed by membranous ossification and we hypothesize that abnormal embryonic cartilage that rests within the suture lines may be involved in this tumor entity. Testing of additional individuals with similar presentations is needed to confirm this finding and clarify the associated phenotypes.

Genomics combined with a protein informatics platform to assess a novel pathogenic variant c.1024 A>G (p.K342E) in OPA1 in a patient with autosomal dominant optic atrophy

BACKGROUND

Autosomal Dominant Optic Atrophy (ADOA) is caused by mutations in the Optic Atrophy 1 Gene which disrupts the OPA1 protein. This disruption affects the normal function of the protein; impairs fusion of the mitochondrial inner membrane; and prevents normal OPA1 protein degradation. These events cause damage in retinal ganglion cells that could affect the patients with symptoms ranging from none to legally blind.

MATERIALS AND METHODS

Our study identifies a missense variant mutation, c.1024 A > G (p.K342E), in OPA1 gene causing ADOA. Diagnosed clinically in three family members and the presence of this mutation was confirmed in two members by genetic testing. Pathogenic variants in OPA1 impact the secondary protein structure and function by causing non-conservative amino acid substitutions. We also modeled this mutation and compared it to the wild type using statistical mechanics.

RESULTS AND CONCLUSIONS

The proband's pathogenic variant, c.1024 A > G (p.K342E), is located in the GTPase domain of OPA1 and causes changes in the protein structure by affecting the oligomerization pattern thus resulting in ADOA. Identifying the pathogenic potential of the missense mutations in the OPA1 gene using neoteric protein modeling techniques would help in the early detection of ADOA in patients who have family history of blindness. This action would help in providing early follow up, possible treatment in the future, and genetic counseling. Abbreviations: ADOA: Autosomal Dominant Optic Atrophy; CYCS: Caspase Activator Cytochrome C; OPA1: Optic Atrophy Gene 1; RGC: Retinal Ganglion Cells; VUS: Variant of Uncertain Significance.

Structural Models for the Dynamic Effects of Loss-of-Function Variants in the Human SIM1 Protein Transcriptional Activation Domain

Single-minded homologue 1 (SIM1) is a transcription factor with numerous different physiological and developmental functions. SIM1 is a member of the class I basic helix-loop-helix-PER-ARNT-SIM (bHLH-PAS) transcription factor family, that includes several other conserved proteins, including the hypoxia-inducible factors, aryl hydrocarbon receptor, neuronal PAS proteins, and the CLOCK circadian regulator. Recent studies of HIF-a-ARNT and CLOCK-BMAL1 protein complexes have revealed the organization of their bHLH, PASA, and PASB domains and provided insight into how these heterodimeric protein complexes form; however, experimental structures for SIM1 have been lacking. Here, we describe the first full-length atomic structural model for human SIM1 with its binding partner ARNT in a heterodimeric complex and analyze several pathogenic variants utilizing state-of-the-art simulations and algorithms. Using local and global positional deviation metrics, deductions to the structural basis for the individual mutants are addressed in terms of the deleterious structural reorganizations that could alter protein function. We propose new experiments to probe these hypotheses and examine an interesting SIM1 dynamic behavior. The conformational dynamics demonstrates conformational changes on local and global regions that represent a mechanism for dysfunction in variants presented. In addition, we used our ab initio hybrid model for further prediction of variant hotspots that can be engineered to test for counter variant (restoration of wild-type function) or basic research probe.

Structural And Computational Perspectives of Selectively Targeting Mutant Proteins

Diseases are often caused by mutant proteins. Many drugs have limited effectiveness and/or toxic side effects because of a failure to selectively target the disease-causing mutant variant, rather than the functional wild type protein. Otherwise, the drugs may even target different proteins with similar structural features. Designing drugs that successfully target mutant proteins selectively represents a major challenge. Decades of cancer research have led to an abundance of potential therapeutic targets, often touted to be "master regulators". For many of these proteins, there are no FDA-approved drugs available; for others, off-target effects result in dose-limiting toxicity. Cancer-related proteins are an excellent medium to carry the story of mutant-specific targeting, as the disease is both initiated and sustained by mutant proteins; furthermore, current chemotherapies generally fail at adequate selective distinction. This review discusses some of the challenges associated with selective targeting from a structural biology perspective, as well as some of the developments in algorithm approach and computational workflow that can be applied to address those issues. One of the most widely researched proteins in cancer biology is p53, a tumor suppressor. Here, p53 is discussed as a specific example of a challenging target, with contemporary drugs and methodologies used as examples of burgeoning successes. The oncogene KRAS, which has been described as "undruggable", is another extensively investigated protein in cancer biology. This review also examines KRAS to exemplify progress made towards selective targeting of diseasecausing mutant proteins. Finally, possible future directions relevant to the topic are discussed.

Functional Analysis of the SIM1 Variant p.G715V in 2 Patients With Obesity

CONTEXT

Single-minded homologue 1 (SIM1) is a transcription factor with several physiological and developmental functions. Haploinsufficiency of SIM1 is associated with early-onset obesity with or without Prader-Willi-like (PWL) features and may exhibit incomplete penetrance.

CASE DESCRIPTION

Next-generation sequencing was performed for 2 male patients with obesity, including 1 man presenting with intellectual disability (ID), body mass index (BMI) of 47.4, and impulse-control disorder, and the other man with early obesity (BMI of 36); sequencing revealed a missense variant in SIM1 (c.2144G>T; p.G715V) in both individuals. Previous studies have identified several disease-associated variants that fall near the p.G715V variant within the C-terminal domain of SIM1. We examined p.G715V variant stability and activity in a doxycycline-inducible stable cell line transfected with an artificial reporter construct and either ARNT or ARNT2 as a partner protein.

CONCLUSIONS

Functional testing of the p.G715V variant revealed a significant reduction in SIM1-mediated transcriptional activity. We also generated the first ab initio hybrid protein model for full-length SIM1 to show the predicted spatial relationship between p.G715V and other previously described variants in this region and identified a putative mutation hotspot within the C-terminus. Significant clinical heterogeneity has been observed in patients with SIM1 variants, particularly with regards to the PWL phenotype. In the patient with ID, a second variant of uncertain significance in CHD2 was identified that may contribute to his ID and behavioral disturbances, emphasizing the role of additional genetic modifiers.

Molecular Dissection of FUS Points at Synergistic Effect of Low-Complexity Domains in Toxicity

RNA-binding protein aggregation is a pathological hallmark of several neurodegenerative disorders, including amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). To gain better insight into the molecular interactions underlying this process, we investigated FUS, which is mutated and aggregated in both ALS and FTLD. We generated a Drosophila model of FUS toxicity and identified a previously unrecognized synergistic effect between the N-terminal prion-like domain and the C-terminal arginine-rich domain to mediate toxicity. Although the prion-like domain is generally considered to mediate aggregation of FUS, we find that arginine residues in the C-terminal low-complexity domain are also required for maturation of FUS in cellular stress granules. These data highlight an important role for arginine-rich domains in the pathology of RNA-binding proteins.

•

Both QGSY and RGG2 domains are necessary for FUS-induced neurodegeneration in flies

•

Arginine-rich domains interact with QGSY hydrogels and liquid droplets

•

RGG2 arginines promote phase separation of FUS in vitro and in cells

•

FUS phase separation behavior in vitro correlates with neurodegeneration in vivo

Development of multi-drug loaded PEGylated nanodiamonds to inhibit tumor growth and metastasis in genetically engineered mouse models of pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease. Nanomedicine, however, offers new opportunities to facilitate drug delivery in PDAC. Our previous work has shown that poly(ethylene glycol)-functionalized nanodiamond (ND) mediated drug delivery offered a considerable improvement over free drug in PDAC. Inspired by this result and guided by molecular simulations, we opted for simultaneous loading of irinotecan and curcumin in ultra-small PEGylated NDs (ND-IRT + CUR). We observed that ND-IRT + CUR was more efficacious in killing AsPC-1 and PANC-1 cells than NDs with single drugs. Using NDs functionalized with a near-infrared (NIR) dye, we demonstrated the preferential localization of the NDs in tumors and metastatic lesions. We further demonstrate that ND-IRT + CUR is capable of producing pronounced anti-tumor effects in two different clinically relevant, immune-competent genetic models of PDAC. Cytokine profiling indicated that NDs with or without drugs downregulated the expression of IL-10, a key modulator of the tumor microenvironment. Thus, using a combination of in silico, in vitro, and in vivo approaches, we show for the first time the remarkable anti-tumor efficacy of PEGylated NDs carrying a dual payload of irinotecan plus curcumin. These results highlight the potential use of such nano-carriers in the treatment of patients with pancreatic cancer.

Design and Evaluation of PEGylated Liposomal Formulation of a Novel Multikinase Inhibitor for Enhanced Chemosensitivity and Inhibition of Metastatic Pancreatic Ductal Adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) has one of the highest mortality rates among cancers. Chemotherapy is the standard first-line treatment, but only modest survival benefits are observed. With the advent of targeted therapies, epidermal growth factor receptor (EGFR) has been acknowledged as a prospective target in PDAC since it is overexpressed in up to 60% of cases. Similarly, the tyrosine-protein kinase Met (cMET) is also overexpressed in PDAC (27-60%) and is a prognostic marker for poor survival. Interestingly, EGFR and cMET share some common signaling pathways including PI3K/Akt and MAPK pathways. Small molecule inhibitors or bispecific antibodies that can target both EGFR and cMET are therefore emerging as novel options for cancer therapy. We previously developed a dual EGFR and cMET inhibitor (N19) that was able to inhibit tumor growth in nonsmall cell lung cancer models resistant to EGFR tyrosine kinase inhibitors (TKI). Here, we report the development of a novel liposomal formulation of N19 (LN19) and showed significant growth inhibition and increased sensitivity toward gemcitabine in the pancreatic adenocarcinoma orthotopic xenograft model. Taken together, our results suggest that LN19 can be valued as an effective combination therapy with conventional chemotherapy such as gemcitabine for PDAC patients.

Molecular Inhibitor of QSOX1 Suppresses Tumor Growth In Vivo

Quiescin sulfhydryl oxidase 1 (QSOX1) is an enzyme overexpressed by many different tumor types. QSOX1 catalyzes the formation of disulfide bonds in proteins. Because short hairpin knockdowns (KD) of QSOX1 have been shown to suppress tumor growth and invasion in vitro and in vivo, we hypothesized that chemical compounds inhibiting QSOX1 enzymatic activity would also suppress tumor growth, invasion, and metastasis. High throughput screening using a QSOX1-based enzymatic assay revealed multiple potential QSOX1 inhibitors. One of the inhibitors, known as "SBI-183," suppresses tumor cell growth in a Matrigel-based spheroid assay and inhibits invasion in a modified Boyden chamber, but does not affect viability of nonmalignant cells. Oral administration of SBI-183 inhibits tumor growth in 2 independent human xenograft mouse models of renal cell carcinoma. We conclude that SBI-183 warrants further exploration as a useful tool for understanding QSOX1 biology and as a potential novel anticancer agent in tumors that overexpress QSOX1.

Integrative data fusion for comprehensive assessment of a novel CHEK2 variant using combined genomics, imaging, and functional-structural assessments via protein informatics

The CHEK2 gene and its encoded protein Chk2 have a well-known role in cancers, especially those related to breast cancer mediated through the BRCA1 gene. Additionally Chk2 has a crucial role in DNA repair, apoptosis and the cell cycle, which is why classification of variants of uncertain significance (VUS) is an area highly sought for a better elucidation of the "genomic effect" that results. Because it can often take years before enough clinical data is accumulated, and the costly and expensive functional analysis for individual variants presents a significant hurdle, it is important to identify other tools to help aid in clarifying the impact of specific variants on a protein's function and eventually the patient's health outcome. Here we describe a newly identified CHEK2 variant and analyze with an integrated approach combining genomics (whole exome analysis), clinical study, radiographic imaging, and protein informatics to identify and predict the functional impact of the VUS on the protein's behavior and predicted impact on the related pathways. The observed and analyzed defects in the protein were consistent with the expected clinical effect. Here, we support the use of personalized protein modeling and informatics and further our goal of developing a large-scale protein deposition archive for all protein-level VUS.

Disulfide engineering of human Kunitz-type serine protease inhibitors enhances proteolytic stability and target affinity toward mesotrypsin

Serine protease inhibitors of the Kunitz-bovine pancreatic trypsin inhibitor (BPTI) family are ubiquitous biological regulators of proteolysis. These small proteins are resistant to proteolysis, but can be slowly cleaved within the protease-binding loop by target proteases, thereby compromising their activity. For the human protease mesotrypsin, this cleavage is especially rapid. Here, we aimed to stabilize the Kunitz domain structure against proteolysis through disulfide engineering. Substitution within the Kunitz inhibitor domain of the amyloid precursor protein (APPI) that incorporated a new disulfide bond between residues 17 and 34 reduced proteolysis by mesotrypsin 74-fold. Similar disulfide engineering of tissue factor pathway inhibitor-1 Kunitz domain 1 (^KD1TFPI1) and bikunin Kunitz domain 2 (^KD2bikunin) likewise stabilized these inhibitors against mesotrypsin proteolysis 17- and 6.6-fold, respectively. Crystal structures of disulfide-engineered APPI and ^KD1TFPI1 variants in a complex with mesotrypsin at 1.5 and 2.0 Å resolution, respectively, confirmed the formation of well-ordered disulfide bonds positioned to stabilize the binding loop. Long all-atom molecular dynamics simulations of disulfide-engineered Kunitz domains and their complexes with mesotrypsin revealed conformational stabilization of the primed side of the inhibitor-binding loop by the engineered disulfide, along with global suppression of conformational dynamics in the Kunitz domain. Our findings suggest that the Cys-17-Cys-34 disulfide slows proteolysis by dampening conformational fluctuations in the binding loop and minimizing motion at the enzyme-inhibitor interface. The generalizable approach developed here for the stabilization against proteolysis of Kunitz domains, which can serve as important scaffolds for therapeutics, may thus find applications in drug development.

Protein informatics combined with multiple data sources enriches the clinical characterization of novel TRPV4 variant causing an intermediate skeletal dysplasia

Background

Transient receptor potential cation channel subfamily V member 4 (TRPV4) is an ion channel permeable to Ca²⁺ that is sensitive to physical, hormonal, and chemical stimuli. This protein is expressed in many cell types, including osteoclasts, chondrocytes, and sensory neurons. As such, pathogenic variants of this gene are associated with skeletal dysplasias and neuromuscular disorders. Pathogenesis of these phenotypes is not yet completely understood, but it is known that genotype–phenotype correlations for TRPV4 pathogenic variants often are not present.

Methods

Newly characterized, suspected pathogenic variant in TRPV4 was analyzed using protein informatics and personalized protein‐level molecular studies, genomic exome analysis, and clinical study.

Results

This statement is demonstrated in the family of our proband, a 47‐year‐old female having the novel c.2401A>G (p.K801E) variant of TRPV4. We discuss the common symptoms between the proband, her father, and her daughter, and compare her phenotype to known TRPV4‐associated skeletal dysplasias.

Conclusions

Protein informatics and molecular modeling are used to confirm the pathogenicity of the unique TRPV4 variant found in this family. Multiple data were combined in a comprehensive manner to give complete overall perspective on the patient disease and prognosis.

Protein modeling and clinical description of a novel in-frame GLB1 deletion causing GM1 gangliosidosis type II

Background

Beta‐galactosidase‐1 (GLB1) is a lysosomal hydrolase that is responsible for breaking down specific glycoconjugates, particularly GM1 (monosialotetrahexosylganglioside). Pathogenic variants in GLB1 cause two different lysosomal storage disorders: GM1 gangliosidosis and mucopolysaccharidosis type IVB. In GM1 gangliosidosis, decreased β‐galactosidase‐1 enzymatic activity leads to the accumulation of GM1 gangliosides, predominantly within the CNS. We present a 22‐month‐old proband with GM1 gangliosidosis type II (late‐infantile form) in whom a novel homozygous in‐frame deletion (c.1468_1470delAAC, p.Asn490del) in GLB1 was detected.

Methods

We used an experimental protein structure of β‐galactosidase‐1 to generate a model of the p.Asn490del mutant and performed molecular dynamic simulations to determine whether this mutation leads to altered ligand positioning compared to the wild‐type protein. In addition, residual mutant enzyme activity in patient leukocytes was evaluated using a fluorometric assay.

Results

Molecular dynamics simulations showed the deletion to alter the catalytic site leading to misalignment of the catalytic residues and loss of collective motion within the model. We predict this misalignment will lead to impaired catalysis of β‐galactosidase‐1 substrates. Enzyme assays confirmed diminished GLB1 enzymatic activity (~3% of normal activity) in the proband.

Conclusions

We have described a novel, pathogenic in‐frame deletion of GLB1 in a patient with GM1 gangliosidosis type II.

Protein molecular modeling techniques investigating novel TAB2 variant R347X causing cardiomyopathy and congenital heart defects in multigenerational family

BACKGROUND

Haploinsufficiency of TAB2 is known to cause congenital heart defects and cardiomyopathy due to its important roles in cardiovascular tissue, both during development and through adult life. We report a sibling pair displaying adult-onset cardiomyopathy, hypermobility, and mild myopia. Our proband, a 39-year-old male, presents only with the above symptoms, while his 36-year-old sister was also notable for a ventricular septal defect in her infancy.

METHODS

Whole-exome sequencing was utilized to identify the molecular basis of the phenotype found in two siblings. A molecular modeling technique that takes advantage of conformational sampling advances (Maxwell's demon molecular dynamics and Monte Carlo) were used to make a model of the mutant variant for comparative analytics to the wild-type.

RESULTS

Exome sequencing revealed a novel, heterogeneous pathogenic variant in TAB2, c.1039 C>T (p.R347X), that was present in both individuals. This pathogenic variant removes just over half the residues from the TAB2 protein and severely impacts its functional ability, which we describe in detail.

CONCLUSIONS

Analysis of the proband's family showed a history of cardiomyopathy, but no congenital heart defects or connective tissue disease. We highlight the heterogeneity in phenotype of TAB2 pathogenic variants and confirm the pathogenicity of this new variant through neoteric protein modeling techniques.

Protein molecular modeling shows residue T599 is critical to wild-type function of POLG and description of a novel variant associated with the SANDO phenotype

Sensory ataxic neuropathy with dysarthria and ophthalmoparesis (SANDO) is a rare phenotype resulting from pathogenic variants of mitochondrial DNA polymerase gamma (POLG). We modeled a novel POLG variant, T599P, that causes the SANDO phenotype and another variant at the same residue, p.T599E, to observe their effect on protein function and confirm the pathogenicity of T599P. Through neoteric molecular modeling techniques, we show that changes at the T599 residue position introduce extra rigidity into the surrounding helix–loop–helix, which places steric pressure on nearby nucleotides. We also provide a clinical description of the T599P variant, which was found in a 42-year-old female proband. The proband presented a 1-year history of progressive gait instability, dysarthria and foot numbness. Her neurologic examination revealed ataxic dysarthria, restricted eye movements, head and palatal tremors, reduced lower limb reflexes, distal multimodal sensory loss and a wide, unsteady ataxic gait. Electromyography studies indicated a sensory neuropathy. Whole-exome sequencing was pursued after tests for infectious, inflammatory and paraneoplastic causes were negative.

Multifaceted peptide assisted one-pot synthesis of gold nanoparticles for plectin-1 targeted gemcitabine delivery in pancreatic cancer

An astute modification of the plectin-1-targeting peptide KTLLPTP by introducing a C-terminal cysteine preceded by a tyrosine residue imparted a reducing property to the peptide. This novel property is then exploited to fabricate gold nanoparticles (GNP) via an in situ reduction of gold(iii) chloride in a one-pot, green synthesis. The modified peptide KTLLPTPYC also acts as a template to generate highly monodispersed, spherical GNPs with a narrow size distribution and improved stability. Plectin-1 is known to be aberrantly expressed in the surface of pancreatic ductal adenocarcinoma (PDAC) cells while showing cytoplasmic expression in normal cells. The synthesized GNPs are thus in situ surface modified with the peptides via the cysteine residue leaving the N-terminal KTLLPTP sequence free for targeting plectin-1. The visual molecular dynamics based simulations support the experimental observations like particle size, gemcitabine conjugation and architecture of the peptide-grafted nanoassembly. Additionally, GNPs conjugated to gemcitabine demonstrate significantly higher cytotoxicity in vitro in two established PDAC cell lines (AsPC-1 and PANC-1) and an admirable in vivo antitumor efficacy in a PANC-1 orthotopic xenograft model through selective uptake in PDAC tumor tissues. Altogether, this strategy represents a unique method for the fabrication of a GNP based targeted drug delivery platform using a multifaceted peptide that acts as reducing agent, template for GNP synthesis and targeting agent to display remarkable selectivity towards PDAC.

Accelerated bottom-up drug design platform enables the discovery of novel stearoyl-CoA desaturase 1 inhibitors for cancer therapy

Here we present an innovative computational-based drug discovery strategy, coupled with machine-based learning and functional assessment, for the rational design of novel small molecule inhibitors of the lipogenic enzyme stearoyl-CoA desaturase 1 (SCD1). Our methods resulted in the discovery of several unique molecules, of which our lead compound SSI-4 demonstrates potent anti-tumor activity, with an excellent pharmacokinetic and toxicology profile. We improve upon key characteristics, including chemoinformatics and absorption/distribution/metabolism/excretion (ADME) toxicity, while driving the IC50 to 0.6 nM in some instances. This approach to drug design can be executed in smaller research settings, applied to a wealth of other targets, and paves a path forward for bringing small-batch based drug programs into the Clinic.

A novel splice site variant in CYP11A1 in trans with the p.E314K variant in a male patient with congenital adrenal insufficiency

BACKGROUND

The CYP11A1 gene encodes the cytochrome P450 side-chain cleavage enzyme, which is essential for steroid formation. Recessive variants in this gene can lead to impairment of sexual differentiation caused by a complete or partial loss of steroid hormone production. The phenotypic spectrum in affected 46XY males may vary from surgically repairable defects including cryptorchidism and hypospadias to complete feminization of external gonads, accompanied by symptoms of adrenal dysfunction.

METHODS

Whole-exome sequencing (WES) of a 12-year-old male proband and his parents was performed after a protracted diagnostic odyssey failed to uncover the cause of his primary adrenal insufficiency. Of note, the proband had early symptomatology and corrective surgery for hypospadias, raising suspicion for a disorder of steroidogenesis.

RESULTS

WES identified compound heterozygous variants in CYP11A1 including a novel canonical splice site variant (c.425+1G>A) and a previously reported p.E314K variant, which were consistent with a diagnosis of congenital adrenal insufficiency with partial 46XY sex reversal.

CONCLUSION

Congenital adrenal insufficiency with 46XY sex reversal is a rare disorder that is characterized by dysregulation of steroid hormone synthesis, leading to adrenal and gonadal dysfunction. In this report, we describe a patient with adrenal insufficiency, hypospadias, and skin hyperpigmentation who was found to have a novel c.425+1G>A variant in trans with the p.E314K variant in CYP11A1. We performed structural analyses to examine the effect of the p.E314K variant on protein function and show that it falls in the core of the protein may disrupt cholesterol binding in the active site.

Small molecule inhibitors of mesotrypsin from a structure-based docking screen

PRSS3/mesotrypsin is an atypical isoform of trypsin, the upregulation of which has been implicated in promoting tumor progression. To date there are no mesotrypsin-selective pharmacological inhibitors which could serve as tools for deciphering the pathological role of this enzyme, and could potentially form the basis for novel therapeutic strategies targeting mesotrypsin. A virtual screen of the Natural Product Database (NPD) and Food and Drug Administration (FDA) approved Drug Database was conducted by high-throughput molecular docking utilizing crystal structures of mesotrypsin. Twelve high-scoring compounds were selected for testing based on lowest free energy docking scores, interaction with key mesotrypsin active site residues, and commercial availability. Diminazene (CID22956468), along with two similar compounds presenting the bis-benzamidine substructure, was validated as a competitive inhibitor of mesotrypsin and other human trypsin isoforms. Diminazene is the most potent small molecule inhibitor of mesotrypsin reported to date with an inhibitory constant (Ki) of 3.6±0.3 μM. Diminazene was subsequently co-crystalized with mesotrypsin and the crystal structure was solved and refined to 1.25 Å resolution. This high resolution crystal structure can now offer a foundation for structure-guided efforts to develop novel and potentially more selective mesotrypsin inhibitors based on similar molecular substructures.

The PINK1 p.I368N mutation affects protein stability and ubiquitin kinase activity

Background

Mutations in PINK1 and PARKIN are the most common causes of recessive early-onset Parkinson’s disease (EOPD). Together, the mitochondrial ubiquitin (Ub) kinase PINK1 and the cytosolic E3 Ub ligase PARKIN direct a complex regulated, sequential mitochondrial quality control. Thereby, damaged mitochondria are identified and targeted to degradation in order to prevent their accumulation and eventually cell death. Homozygous or compound heterozygous loss of either gene function disrupts this protective pathway, though at different steps and by distinct mechanisms. While structure and function of PARKIN variants have been well studied, PINK1 mutations remain poorly characterized, in particular under endogenous conditions. A better understanding of the exact molecular pathogenic mechanisms underlying the pathogenicity is crucial for rational drug design in the future.

Methods

Here, we characterized the pathogenicity of the PINK1 p.I368N mutation on the clinical and genetic as well as on the structural and functional level in patients’ fibroblasts and in cell-based, biochemical assays.

Results

Under endogenous conditions, PINK1 p.I368N is expressed, imported, and N-terminally processed in healthy mitochondria similar to PINK1 wild type (WT). Upon mitochondrial damage, however, full-length PINK1 p.I368N is not sufficiently stabilized on the outer mitochondrial membrane (OMM) resulting in loss of mitochondrial quality control. We found that binding of PINK1 p.I368N to the co-chaperone complex HSP90/CDC37 is reduced and stress-induced interaction with TOM40 of the mitochondrial protein import machinery is abolished. Analysis of a structural PINK1 p.I368N model additionally suggested impairments of Ub kinase activity as the ATP-binding pocket was found deformed and the substrate Ub was slightly misaligned within the active site of the kinase. Functional assays confirmed the lack of Ub kinase activity.

Conclusions

Here we demonstrated that mutant PINK1 p.I368N can not be stabilized on the OMM upon mitochondrial stress and due to conformational changes in the active site does not exert kinase activity towards Ub. In patients’ fibroblasts, biochemical assays and by structural analyses, we unraveled two pathomechanisms that lead to loss of function upon mutation of p.I368N and highlight potential strategies for future drug development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13024-017-0174-z) contains supplementary material, which is available to authorized users.

PINK1, Parkin, and Mitochondrial Quality Control: What can we Learn about Parkinson's Disease Pathobiology?

The first clinical description of Parkinson's disease (PD) will embrace its two century anniversary in 2017. For the past 30 years, mitochondrial dysfunction has been hypothesized to play a central role in the pathobiology of this devastating neurodegenerative disease. The identifications of mutations in genes encoding PINK1 (PTEN-induced kinase 1) and Parkin (E3 ubiquitin ligase) in familial PD and their functional association with mitochondrial quality control provided further support to this hypothesis. Recent research focused mainly on their key involvement in the clearance of damaged mitochondria, a process known as mitophagy. It has become evident that there are many other aspects of this complex regulated, multifaceted pathway that provides neuroprotection. As such, numerous additional factors that impact PINK1/Parkin have already been identified including genes involved in other forms of PD. A great pathogenic overlap amongst different forms of familial, environmental and even sporadic disease is emerging that potentially converges at the level of mitochondrial quality control. Tremendous efforts now seek to further detail the roles and exploit PINK1 and Parkin, their upstream regulators and downstream signaling pathways for future translation. This review summarizes the latest findings on PINK1/Parkin-directed mitochondrial quality control, its integration and cross-talk with other disease factors and pathways as well as the implications for idiopathic PD. In addition, we highlight novel avenues for the development of biomarkers and disease-modifying therapies that are based on a detailed understanding of the PINK1/Parkin pathway.

Coinhibition of the deubiquitinating enzymes, USP14 and UCHL5, with VLX1570 is lethal to ibrutinib- or bortezomib-resistant Waldenstrom macroglobulinemia tumor cells

The survival of Waldenstrom macroglobulinemia (WM) tumor cells hinges on aberrant B-cell receptor (BCR) and MYD88 signaling. WM cells upregulate the proteasome function to sustain the BCR-driven growth while maintaining homeostasis. Clinically, two treatment strategies are used to disrupt these complementary yet mutually exclusive WM survival pathways via ibrutinib (targets BTK/MYD88 node) and bortezomib (targets 20 S proteasome). Despite the success of both agents, WM patients eventually become refractory to treatment, highlighting the adaptive plasticity of WM cells and underscoring the need for development of new therapeutics. Here we provide a comprehensive preclinical report on the anti-WM activity of VLX1570, a novel small-molecule inhibitor of the deubiquitinating enzymes (DUBs), ubiquitin-specific protease 14 (USP14) and ubiquitin carboxyl-terminal hydrolase isozyme L5 (UCHL5). Both DUBs reside in the 19 S proteasome cap and their inhibition by VLX1570 results in rapid and tumor-specific apoptosis in bortezomib- or ibrutinib-resistant WM cells. Notably, treatment of WM cells with VLX1570 downregulated BCR-associated elements BTK, MYD88, NFATC, NF-κB and CXCR4, the latter whose dysregulated function is linked to ibrutinib resistance. VLX1570 administered to WM-xenografted mice resulted in decreased tumor burden and prolonged survival (P=0.0008) compared with vehicle-treated mice. Overall, our report demonstrates significant value in targeting USP14/UCHL5 with VLX1570 in drug-resistant WM and carries a high potential for clinical translation.

An Acrobatic Substrate Metamorphosis Reveals a Requirement for Substrate Conformational Dynamics in Trypsin Proteolysis

The molecular basis of enzyme catalytic power and specificity derives from dynamic interactions between enzyme and substrate during catalysis. Although considerable effort has been devoted to understanding how conformational dynamics within enzymes affect catalysis, the role of conformational dynamics within protein substrates has not been addressed. Here, we examine the importance of substrate dynamics in the cleavage of Kunitz-bovine pancreatic trypsin inhibitor protease inhibitors by mesotrypsin, finding that the varied conformational dynamics of structurally similar substrates can profoundly impact the rate of catalysis. A 1.4-Å crystal structure of a mesotrypsin-product complex formed with a rapidly cleaved substrate reveals a dramatic conformational change in the substrate upon proteolysis. By using long all-atom molecular dynamics simulations of acyl-enzyme intermediates with proteolysis rates spanning 3 orders of magnitude, we identify global and local dynamic features of substrates on the nanosecond-microsecond time scale that correlate with enzymatic rates and explain differential susceptibility to proteolysis. By integrating multiple enhanced sampling methods for molecular dynamics, we model a viable conformational pathway between substrate-like and product-like states, linking substrate dynamics on the nanosecond-microsecond time scale with large collective substrate motions on the much slower time scale of catalysis. Our findings implicate substrate flexibility as a critical determinant of catalysis.

Heterozygous PINK1 p.G411S increases risk of Parkinson's disease via a dominant-negative mechanism

See Gandhi and Plun-Favreau (doi:10.1093/aww320) for a scientific commentary on this article.

Heterozygous mutations in recessive Parkinson’s disease genes have been postulated to increase disease risk. Puschmann et al. report a genetic association between heterozygous PINK1 p.G411S and Parkinson’s disease. They provide structural and functional explanations for a partial dominant-negative effect of the mutant protein, which impairs wild-type PINK1 activity through hetero-dimerization.

See Gandhi and Plun-Favreau (doi:10.1093/aww320) for a scientific commentary on this article.

It has been postulated that heterozygous mutations in recessive Parkinson’s genes may increase the risk of developing the disease. In particular, the PTEN-induced putative kinase 1 (PINK1) p.G411S (c.1231G>A, rs45478900) mutation has been reported in families with dominant inheritance patterns of Parkinson’s disease, suggesting that it might confer a sizeable disease risk when present on only one allele. We examined families with PINK1 p.G411S and conducted a genetic association study with 2560 patients with Parkinson’s disease and 2145 control subjects. Heterozygous PINK1 p.G411S mutations markedly increased Parkinson’s disease risk (odds ratio = 2.92, P = 0.032); significance remained when supplementing with results from previous studies on 4437 additional subjects (odds ratio = 2.89, P = 0.027). We analysed primary human skin fibroblasts and induced neurons from heterozygous PINK1 p.G411S carriers compared to PINK1 p.Q456X heterozygotes and PINK1 wild-type controls under endogenous conditions. While cells from PINK1 p.Q456X heterozygotes showed reduced levels of PINK1 protein and decreased initial kinase activity upon mitochondrial damage, stress-response was largely unaffected over time, as expected for a recessive loss-of-function mutation. By contrast, PINK1 p.G411S heterozygotes showed no decrease of PINK1 protein levels but a sustained, significant reduction in kinase activity. Molecular modelling and dynamics simulations as well as multiple functional assays revealed that the p.G411S mutation interferes with ubiquitin phosphorylation by wild-type PINK1 in a heterodimeric complex. This impairs the protective functions of the PINK1/parkin-mediated mitochondrial quality control. Based on genetic and clinical evaluation as well as functional and structural characterization, we established p.G411S as a rare genetic risk factor with a relatively large effect size conferred by a partial dominant-negative function phenotype.

Selective Targeting of Extracellular Insulin-Degrading Enzyme by Quasi-Irreversible Thiol-Modifying Inhibitors

Many therapeutically important enzymes are present in multiple cellular compartments, where they can carry out markedly different functions; thus, there is a need for pharmacological strategies to selectively manipulate distinct pools of target enzymes. Insulin-degrading enzyme (IDE) is a thiol-sensitive zinc-metallopeptidase that hydrolyzes diverse peptide substrates in both the cytosol and the extracellular space, but current genetic and pharmacological approaches are incapable of selectively inhibiting the protease in specific subcellular compartments. Here, we describe the discovery, characterization, and kinetics-based optimization of potent benzoisothiazolone-based inhibitors that, by virtue of a unique quasi-irreversible mode of inhibition, exclusively inhibit extracellular IDE. The mechanism of inhibition involves nucleophilic attack by a specific active-site thiol of the enzyme on the inhibitors, which bear an isothiazolone ring that undergoes irreversible ring opening with the formation of a disulfide bond. Notably, binding of the inhibitors is reversible under reducing conditions, thus restricting inhibition to IDE present in the extracellular space. The identified inhibitors are highly potent (IC50(app) = 63 nM), nontoxic at concentrations up to 100 μM, and appear to preferentially target a specific cysteine residue within IDE. These novel inhibitors represent powerful new tools for clarifying the physiological and pathophysiological roles of this poorly understood protease, and their unusual mechanism of action should be applicable to other therapeutic targets.

Abstract 4459: Accelerated drug discovery platform yields synthesis of novel stearoyl-CoA desaturase 1 inhibitors that demonstrate anti-tumor efficacy in several models of aggressive cancer

Background

Intra-tumor heterogeneity is the driving force behind disease progression and metastasis, drug resistance, and relapse in cancer patients. Defining personalized treatment regimens and increasing the arsenal of targeted therapies available will provide the best chance of prolonging survival and finding curative treatments. Recent work by our group and others has shed light on aberrant de novo lipogenesis as a feature of many cancers including kidney, breast, and lung cancer. Of the molecules involved in fatty acid metabolism, stearoyl-CoA desaturase 1 (SCD1) is frequently over-expressed. Targeted inhibition of SCD1 demonstrates a loss of tumor cell viability in a variety of cancer models. While there currently are specific small molecule inhibitors for SCD1 none are under investigation as anti-cancer therapeutics. In order to evaluate the efficacy of SCD1 inhibitors against ccRCC in both a preclinical and eventually a clinical setting we sought to generate unique compounds and test their biological efficacy.

Methods

Utilizing an innovative in silico approach, we designed new inhibitors for SCD1 via a scaffold hopping approach while searching for different “core” scaffolds. Critical interaction moieties from the chemical R-groups were held fixed, while generated cores were rapidly scanned for best-fit criteria (shape, docking, pharmacophore fit (QSAR)) into our Z-matrix. Following de novo compound generation and synthesis, we implemented iterative analog generation for the top hits. Using this method and diversity of chemical space, we pursued three divergent species of compounds into the nanomolar regime.

Biological testing of the novel compounds included high-throughput proliferative-based screening, oleic acid rescue, and evaluation of the endoplasmic stress response. Enzymatic inhibition of SCD1 was evaluated by LC-MS. Preclinical pharmacokinetics was established using a range of doses administered intravenously and orally. Anti-tumor response was evaluated in an ectopic model of renal carcinoma, both as a single agent and in combination with SOC.

Results

A cohort of potential unique small molecule inhibitors against SCD1 were designed and synthesized. Of these, a select group of compounds demonstrated potent and selective enzyme inhibition and reduced tumor growth in in vitro and in vivo models of ccRCC. These novel inhibitors reliably reproduced activation of ER stress markers, authenticating this biological response as a harbinger for successful abrogation of SCD1 activity.

Conclusions

We have established an effective, expedited, and economical format for in silico modeling of novel agents based on a scaffold-hopping methodology. Our findings provide compelling evidence supporting the therapeutic benefit of SCD1 inhibitors for clinical use as broad-spectrum anti-tumor agents alone or in combination with SOC.

Citation Format: Christina A. Von Roemeling, Thomas R. Caulfield, Derek C. Radisky, Ilah Bok, Laura A. Marlow, James Miller, Mojda Sidiqi, Anthony B. Pinkerton, Winston W. Tan, Amy L. Lane, Han W. Tun, John A. Copland. Accelerated drug discovery platform yields synthesis of novel stearoyl-CoA desaturase 1 inhibitors that demonstrate anti-tumor efficacy in several models of aggressive cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4459. doi:10.1158/1538-7445.AM2015-4459

(Patho-)physiological relevance of PINK1-dependent ubiquitin phosphorylation

Mutations in PINK1 and PARKIN cause recessive, early-onset Parkinson's disease (PD). Together, these two proteins orchestrate a protective mitophagic response that ensures the safe disposal of damaged mitochondria. The kinase PINK1 phosphorylates ubiquitin (Ub) at the conserved residue S65, in addition to modifying the E3 ubiquitin ligase Parkin. The structural and functional consequences of Ub phosphorylation (pS65-Ub) have already been suggested from in vitro experiments, but its (patho-)physiological significance remains unknown. We have generated novel antibodies and assessed pS65-Ub signals in vitro and in cells, including primary neurons, under endogenous conditions. pS65-Ub is dependent on PINK1 kinase activity as confirmed in patient fibroblasts and postmortem brain samples harboring pathogenic mutations. We show that pS65-Ub is reversible and barely detectable under basal conditions, but rapidly induced upon mitochondrial stress in cells and amplified in the presence of functional Parkin. pS65-Ub accumulates in human brain during aging and disease in the form of cytoplasmic granules that partially overlap with mitochondrial, lysosomal, and total Ub markers. Additional studies are now warranted to further elucidate pS65-Ub functions and fully explore its potential for biomarker or therapeutic development.

Structural and Functional Impact of Parkinson Disease-Associated Mutations in the E3 Ubiquitin Ligase Parkin

Mutations in the PARKIN/PARK2 gene that result in loss-of-function of the encoded, neuroprotective E3 ubiquitin ligase Parkin cause recessive, familial early-onset Parkinson disease. As an increasing number of rare Parkin sequence variants with unclear pathogenicity are identified, structure-function analyses will be critical to determine their disease relevance. Depending on the specific amino acids affected, several distinct pathomechanisms can result in loss of Parkin function. These include disruption of overall Parkin folding, decreased solubility, and protein aggregation. However pathogenic effects can also result from misregulation of Parkin autoinhibition and of its enzymatic functions. In addition, interference of binding to coenzymes, substrates, and adaptor proteins can affect its catalytic activity too. Herein, we have performed a comprehensive structural and functional analysis of 21 PARK2 missense mutations distributed across the individual protein domains. Using this combined approach, we were able to pinpoint some of the pathogenic mechanisms of individual sequence variants. Similar analyses will be critical in gaining a complete understanding of the complex regulations and enzymatic functions of Parkin. These studies will not only highlight the important residues, but will also help to develop novel therapeutics aimed at activating and preserving an active, neuroprotective form of Parkin.

Synthesis and evaluation of derivatives of the proteasome deubiquitinase inhibitor b-AP15

The ubiquitin-proteasome system (UPS) is increasingly recognized as a therapeutic target for the development of anticancer therapies. The success of the 20S proteasome core particle (20S CP) inhibitor bortezomib in the clinical management of multiple myeloma has raised the possibility of identifying other UPS components for therapeutic intervention. We previously identified the small molecule b-AP15 as an inhibitor of 19S proteasome deubiquitinase (DUB) activity. Building upon our previous data, we performed a structure-activity relationship (SAR) study on b-AP15 and identified VLX1570 as an analog with promising properties, including enhanced potency and improved solubility in aqueous solution. In silico modeling was consistent with interaction of VLX1570 with key cysteine residues located at the active sites of the proteasome DUBs USP14 and UCHL5. VLX1570 was found to inhibit proteasome deubiquitinase activity in vitro in a manner consistent with competitive inhibition. Furthermore, using active-site-directed probes, VLX1570 also inhibited proteasome DUB activity in exposed cells. Importantly, VLX1570 did not show inhibitory activity on a panel of recombinant non-proteasome DUBs, on recombinant kinases, or on caspase-3 activity, suggesting that VLX1570 is not an overtly reactive general enzyme inhibitor. Taken together, our data shows the chemical and biological properties of VLX1570 as an optimized proteasome DUB inhibitor.

Targeted inhibition of the deubiquitinating enzymes, USP14 and UCHL5, induces proteotoxic stress and apoptosis in Waldenström macroglobulinaemia tumour cells

Deubiquitinase enzymes (DUBs) of the proteasomal 19S regulatory particle are emerging as important therapeutic targets in several malignancies. Here we demonstrate that inhibition of two proteasome-associated DUBs (USP14 and UCHL5) with the small molecule DUB inhibitor b-AP15, results in apoptosis of human Waldenström macroglobulinaemia (WM) cell lines and primary patient-derived WM tumour cells. Importantly, b-AP15 produced proteotoxic stress and apoptosis in WM cells that have acquired resistance to the proteasome inhibitor bortezomib. In silico modelling identified protein residues that were critical for the binding of b-AP15 with USP14 or UCHL5 and proteasome enzyme activity assays confirmed that b-AP15 does not affect the proteolytic capabilities of the 20S proteasome β-subunits. In vitro toxicity from b-AP15 appeared to result from a build-up of ubiquitinated proteins and activation of the endoplasmic reticulum stress response in WM cells, an effect that also disrupted the mitochondria. Focused transcriptome profiling of b-AP15-treated WM cells revealed modulation of several genes regulating cell stress and NF-κB signalling, the latter whose protein translocation and downstream target activation was reduced by b-AP15 in vitro. This is the first report to define the effects and underlying mechanisms associated with inhibition of USP14 and UCHL5 DUB activity in WM tumour cells.

Nimbolide targets BCL2 and induces apoptosis in preclinical models of Waldenströms macroglobulinemia

Neem leaf extract (NLE) has medicinal properties, which have been attributed to its limonoid content. We identified the NLE tetranorterpenoid, nimbolide, as being the key limonoid responsible for the cytotoxicity of NLE in various preclinical models of human B-lymphocyte cancer. Of the models tested, Waldenströms macroglobulinemia (WM) cells were most sensitive to nimbolide, undergoing significant mitochondrial mediated apoptosis. Notably, nimbolide toxicity was also observed in drug-resistant (bortezomib or ibrutinib) WM cells. To identify putative targets of nimbolide, relevant in WM, we used chemoinformatics-based approaches comprised of virtual in silico screening, molecular modeling and target–ligand reverse docking. In silico analysis revealed the antiapoptotic protein BCL2 was the preferential binding partner of nimbolide. The significance of this finding was further tested in vitro in RS4;11 (BCL2-dependent) tumor cells, in which nimbolide induced significantly more apoptosis compared with BCL2 mutated (Jurkat BCL2^Ser70-Ala) cells. Lastly, intraperitoneal administration of nimbolide in WM tumor xenografted mice, significantly reduced tumor growth and IgM secretion in vivo, while modulating the expression of several proteins as seen on immunohistochemistry. Overall, our data demonstrate that nimbolide is highly active in WM cells, as well as other B-cell cancers, and engages BCL2 to exert its cytotoxic activity.

Phosphorylation by PINK1 releases the UBL domain and initializes the conformational opening of the E3 ubiquitin ligase Parkin

Loss-of-function mutations in PINK1 or PARKIN are the most common causes of autosomal recessive Parkinson's disease. Both gene products, the Ser/Thr kinase PINK1 and the E3 Ubiquitin ligase Parkin, functionally cooperate in a mitochondrial quality control pathway. Upon stress, PINK1 activates Parkin and enables its translocation to and ubiquitination of damaged mitochondria to facilitate their clearance from the cell. Though PINK1-dependent phosphorylation of Ser65 is an important initial step, the molecular mechanisms underlying the activation of Parkin's enzymatic functions remain unclear. Using molecular modeling, we generated a complete structural model of human Parkin at all atom resolution. At steady state, the Ub ligase is maintained inactive in a closed, auto-inhibited conformation that results from intra-molecular interactions. Evidently, Parkin has to undergo major structural rearrangements in order to unleash its catalytic activity. As a spark, we have modeled PINK1-dependent Ser65 phosphorylation in silico and provide the first molecular dynamics simulation of Parkin conformations along a sequential unfolding pathway that could release its intertwined domains and enable its catalytic activity. We combined free (unbiased) molecular dynamics simulation, Monte Carlo algorithms, and minimal-biasing methods with cell-based high content imaging and biochemical assays. Phosphorylation of Ser65 results in widening of a newly defined cleft and dissociation of the regulatory N-terminal UBL domain. This motion propagates through further opening conformations that allow binding of an Ub-loaded E2 co-enzyme. Subsequent spatial reorientation of the catalytic centers of both enzymes might facilitate the transfer of the Ub moiety to charge Parkin. Our structure-function study provides the basis to elucidate regulatory mechanisms and activity of the neuroprotective Parkin. This may open up new avenues for the development of small molecule Parkin activators through targeted drug design.

Parkinson's disease (PD) is a devastating neurological condition caused by the selective and progressive degeneration of dopaminergic neurons in the brain. Loss-of-function mutations in the PINK1 or PARKIN genes are the most common causes of recessively inherited PD. Together the encoded proteins coordinate a protective cellular quality control pathway that allows elimination of impaired mitochondria in order to prevent further cellular damage and ultimately death. Although it is known that the kinase PINK1 operates upstream and activates the E3 Ubiquitin ligase Parkin, the molecular mechanisms remain elusive. Here, we combined state-of-the art computational and functional biological methods to demonstrate that Parkin is sequentially activated through PINK1-dependent phosphorylation and subsequent structural rearrangement. The induced motions result in release of Parkin's closed, auto-inhibited conformation to liberate its enzymatic functions. We provide for the first time a complete protein structure of Parkin at an all atom resolution and a comprehensive molecular dynamics simulation of its activation and opening conformations. The generated models will allow uncovering the exact mechanisms of regulation and enzymatic activity of Parkin and potentially the development of novel therapeutics through a structure-function-based drug design.

Targeted manipulation of the sortilin-progranulin axis rescues progranulin haploinsufficiency

Progranulin (GRN) mutations causing haploinsufficiency are a major cause of frontotemporal lobar degeneration (FTLD-TDP). Recent discoveries demonstrating sortilin (SORT1) is a neuronal receptor for PGRN endocytosis and a determinant of plasma PGRN levels portend the development of enhancers targeting the SORT1–PGRN axis. We demonstrate the preclinical efficacy of several approaches through which impairing PGRN's interaction with SORT1 restores extracellular PGRN levels. Our report is the first to demonstrate the efficacy of enhancing PGRN levels in iPSC neurons derived from frontotemporal dementia (FTD) patients with PGRN deficiency. We validate a small molecule preferentially increases extracellular PGRN by reducing SORT1 levels in various mammalian cell lines and patient-derived iPSC neurons and lymphocytes. We further demonstrate that SORT1 antagonists and a small-molecule binder of PGRN_588–593, residues critical for PGRN–SORT1 binding, inhibit SORT1-mediated PGRN endocytosis. Collectively, our data demonstrate that the SORT1–PGRN axis is a viable target for PGRN-based therapy, particularly in FTD-GRN patients.