Search

PURPOSE: Patients with relapsed or refractory T-cell acute lymphoblastic leukemia (T-ALL) or lymphoblastic lymphoma (T-LBL) have limited therapeutic options. Clinical use of genomic profiling provides an opportunity to identify targetable alterations to inform therapy. EXPERIMENTAL DESIGN: We describe a cohort of 14 pediatric patients with relapsed or refractory T-ALL enrolled on the Leukemia Precision-based Therapy (LEAP) Consortium trial (NCT02670525) and a patient with T-LBL, discovering alterations in platelet-derived growth factor receptor-α (PDGFRA) in 3 of these patients. We identified a novel mutation in PDGFRA, p.D842N, and used an integrated structural modeling and molecular biology approach to characterize mutations at D842 to guide therapeutic targeting. We conducted a preclinical study of avapritinib in a mouse patient-derived xenograft (PDX) model of FIP1L1-PDGFRA and PDGFRA p.D842N leukemia. RESULTS: Two patients with T-ALL in the LEAP cohort (14%) had targetable genomic alterations affecting PDGFRA, a FIP1-like 1 protein/PDGFRA (FIP1L1-PDGFRA) fusion and a novel mutation in PDGFRA, p.D842N. The D842N mutation resulted in PDGFRA activation and sensitivity to tested PDGFRA inhibitors. In a T-ALL PDX model, avapritinib treatment led to decreased leukemia burden, significantly prolonged survival, and even cured a subset of mice. Avapritinib treatment was well tolerated and yielded clinical benefit in a patient with refractory T-ALL. CONCLUSIONS: Refractory T-ALL has not been fully characterized. Alterations in PDGFRA or other targetable kinases may inform therapy for patients with refractory T-ALL who otherwise have limited treatment options. Clinical genomic profiling, in real time, is needed for fully informed therapeutic decision making.

Preventing the biogenesis of disease-relevant proteins is an attractive therapeutic strategy, but attempts to target essential protein biogenesis factors have been hampered by excessive toxicity. Here we describe KZR-8445, a cyclic depsipeptide that targets the Sec61 translocon and selectively disrupts secretory and membrane protein biogenesis in a signal peptide-dependent manner. KZR-8445 potently inhibits the secretion of pro-inflammatory cytokines in primary immune cells and is highly efficacious in a mouse model of rheumatoid arthritis. A cryogenic electron microscopy structure reveals that KZR-8445 occupies the fully opened Se61 lateral gate and blocks access to the lumenal plug domain. KZR-8445 binding stabilizes the lateral gate helices in a manner that traps select signal peptides in the Sec61 channel and prevents their movement into the lipid bilayer. Our results establish a framework for the structure-guided discovery of novel therapeutics that selectively modulate Sec61-mediated protein biogenesis.

The ability to sense and respond to changes in cellular oxygen levels is critical for aerobic organisms and requires a molecular oxygen sensor. The prototypical sensor is the oxygen-dependent enzyme PHD: hypoxia inhibits its ability to hydroxylate the transcription factor HIF, causing HIF to accumulate and trigger the classic HIF-dependent hypoxia response. A small handful of other oxygen sensors are known, all of which are oxygen-dependent enzymes. However, hundreds of oxygen-dependent enzymes exist among aerobic organisms, raising the possibility that additional sensors remain to be discovered. This review summarizes known and potential hypoxia sensors among human O2-dependent enzymes and highlights their possible roles in hypoxia-related adaptation and diseases.

BackgroundAutosomal dominant mutations in α-synuclein, TDP-43 and tau are thought to predispose to neurodegeneration by enhancing protein aggregation. While a subset of α-synuclein, TDP-43 and tau mutations has been shown to increase the structural propensity of these proteins toward self-association, rates of aggregation are also highly dependent on protein steady state concentrations, which are in large part regulated by their rates of lysosomal degradation. Previous studies have shown that lysosomal proteases operate precisely and not indiscriminately, cleaving their substrates at very specific linear amino acid sequences. With this knowledge, we hypothesized that certain coding mutations in α-synuclein, TDP-43 and tau may lead to increased protein steady state concentrations and eventual aggregation by an alternative mechanism, that is, through disrupting lysosomal protease cleavage recognition motifs and subsequently conferring protease resistance to these proteins.ResultsTo test this possibility, we first generated comprehensive proteolysis maps containing all of the potential lysosomal protease cleavage sites for α-synuclein, TDP-43 and tau. In silico analyses of these maps indicated that certain mutations would diminish cathepsin cleavage, a prediction we confirmed utilizing in vitro protease assays. We then validated these findings in cell models and induced neurons, demonstrating that mutant forms of α-synuclein, TDP-43 and tau are degraded less efficiently than wild type despite being imported into lysosomes at similar rates.ConclusionsTogether, this study provides evidence that pathogenic mutations in the N-terminal domain of α-synuclein (G51D, A53T), low complexity domain of TDP-43 (A315T, Q331K, M337V) and R1 and R2 domains of tau (K257T, N279K, S305N) directly impair their own lysosomal degradation, altering protein homeostasis and increasing cellular protein concentrations by extending the degradation half-lives of these proteins. These results also point to novel, shared, alternative mechanism by which different forms of neurodegeneration, including synucleinopathies, TDP-43 proteinopathies and tauopathies, may arise. Importantly, they also provide a roadmap for how the upregulation of particular lysosomal proteases could be targeted as potential therapeutics for human neurodegenerative disease.

Malaria is among the tropical diseases that cause the most deaths in Africa. Around 500,000 malaria deaths are reported yearly among African children under the age of five. Chloroquine (CQ) is a low-cost antimalarial used worldwide for the treatment of Plasmodium vivax malaria. Due to resistance mechanisms, CQ is no longer effective against most malaria cases caused by P. falciparum. The World Health Organization recommends artemisinin combination therapies for P. falciparum malaria, but resistance is emerging in Southeast Asia and some parts of Africa. Therefore, new medicines for treating malaria are urgently needed. Previously, our group identified the 4-aminoquinoline DAQ, a CQ analog containing an acetylenic bond in its side chain, which overcomes CQ resistance in K1 P. falciparum strains. In this work, the antiplasmodial profile, drug-like properties, and pharmacokinetics of DAQ were further investigated. DAQ showed no cross-resistance against standard CQ-resistant strains (e.g., Dd2, IPC 4912, RF12) nor against P. falciparum and P. vivax isolates from patients in the Brazilian Amazon. Using drug pressure assays, DAQ showed a low propensity to generate resistance. DAQ showed considerable solubility but low metabolic stability. The main metabolite was identified as a mono N-deethylated derivative (DAQM), which also showed significant inhibitory activity against CQ-resistant P. falciparum strains. Our findings indicated that the presence of a triple bond in CQ-analogues may represent a low-cost opportunity to overcome known mechanisms of resistance in the malaria parasite.

Cysteine side chains can exist in distinct oxidation states depending on the pH and redox potential of the environment, and cysteine oxidation plays important yet complex regulatory roles. Compared with the effects of post-translational modifications such as phosphorylation, the effects of oxidation of cysteine to sulfenic, sulfinic, and sulfonic acid on protein structure and function remain relatively poorly characterized. We present an analysis of the role of cysteine reactivity as a regulatory factor in proteins, emphasizing the interplay between electrostatics and redox potential as key determinants of the resulting oxidation state. A review of current computational approaches suggests underdeveloped areas of research for studying cysteine reactivity through molecular simulations.

KIT is a type 3 receptor tyrosine kinase that plays a crucial role in cellular growth and proliferation. Mutations in KIT can dysregulate its active-inactive equilibrium. Activating mutations drive cancer growth, while deactivating mutations result in the loss of skin and hair pigmentation in a disease known as piebaldism. Here, we propose a method based on molecular dynamics and free energy calculations to predict the functional effect of KIT mutations. Our calculations may have important clinical implications by defining the functional significance of previously uncharacterized KIT mutations and guiding targeted therapy.

Because nonopioid analgesics are much sought after, we computationally docked more than 301 million virtual molecules against a validated pain target, the α2A-adrenergic receptor (α2AAR), seeking new α2AAR agonists chemotypes that lack the sedation conferred by known α2AAR drugs, such as dexmedetomidine. We identified 17 ligands with potencies as low as 12 nanomolar, many with partial agonism and preferential Gi and Go signaling. Experimental structures of α2AAR complexed with two of these agonists confirmed the docking predictions and templated further optimization. Several compounds, including the initial docking hit '9087 [mean effective concentration (EC50) of 52 nanomolar] and two analogs, '7075 and PS75 (EC50 4.1 and 4.8 nanomolar), exerted on-target analgesic activity in multiple in vivo pain models without sedation. These newly discovered agonists are interesting as therapeutic leads that lack the liabilities of opioids and the sedation of dexmedetomidine.

Cancer persister cells are able to survive otherwise lethal doses of drugs through nongenetic mechanisms, which can lead to cancer regrowth and drug resistance. The broad spectrum of molecular differences observed between persisters and their treatment-naïve counterparts makes it challenging to identify causal mechanisms underlying persistence. Here, we modulate environmental signals to identify cellular mechanisms that promote the emergence of persisters and to pinpoint actionable vulnerabilities that eliminate them. We found that interferon-γ (IFNγ) can induce a pro-persistence signal that can be specifically eliminated by inhibition of type I protein arginine methyltransferase (PRMT) (PRMTi). Mechanistic investigation revealed that signal transducer and activator of transcription 1 (STAT1) is a key component connecting IFNγ's pro-persistence and PRMTi's antipersistence effects, suggesting a previously unknown application of PRMTi to target persisters in settings with high STAT1 expression. Modulating environmental signals can accelerate the identification of mechanisms that promote and eliminate cancer persistence.

Human platelet 12-(S)-Lipoxygenase (12-LOX) is a fatty acid metabolizing oxygenase that plays an important role in platelet activation and cardiometabolic disease. ML355 is a specific 12-LOX inhibitor that has been shown to decrease thrombosis without prolonging hemostasis and protect human pancreatic islets from inflammatory injury. It has an amenable drug-like scaffold with nM potency and encouraging ADME and PK profiles, but its binding mode to the active site of 12-LOX remains unclear. In the current work, we combined computational modeling and experimental mutagenesis to propose a model in which ML355 conforms to the "U" shape of the 12-LOX active site, with the phenyl linker region wrapping around L407. The benzothiazole of ML355 extends into the bottom of the active site cavity, pointing towards residues A417 and V418. However, reducing the active site depth alone did not affect ML355 potency. In order to lower the potency of ML355, the cavity needed to be reduced in both length and width. In addition, H596 appears to position ML355 in the active site through an interaction with the 2-methoxy phenol moiety of ML355. Combined, this binding model suggested that the benzothiazole of ML355 could be enlarged. Therefore, a naphthyl-benzothiazole derivative of ML355, Lox12Slug001, was synthesized and shown to have 7.2-fold greater potency than ML355. This greater potency is proposed to be due to additional van der Waals interactions and pi-pi stacking with F414 and F352. Lox12Slug001 was also shown to be highly selective against 12-LOX relative to the other LOX isozymes and more importantly, it showed activity in rescuing human islets exposed to inflammatory cytokines with comparable potency to ML355. Further studies are currently being pursued to derivatize ML355 in order to optimize the additional space in the active site, while maintaining acceptable drug-like properties.

Human epithelial 15-lipoxygenase-2 (h15-LOX-2, ALOX15B) is expressed in many tissues and has been implicated in atherosclerosis, cystic fibrosis and ferroptosis. However, there are few reported potent/selective inhibitors that are active ex vivo. In the current work, we report newly discovered molecules that are more potent and structurally distinct from our previous inhibitors, MLS000545091 and MLS000536924 (Jameson et al, PLoS One, 2014, 9, e104094), in that they contain a central imidazole ring, which is substituted at the 1-position with a phenyl moiety and with a benzylthio moiety at the 2-position. The initial three molecules were mixed-type, non-reductive inhibitors, with IC50 values of 0.34 ± 0.05 μM for MLS000327069, 0.53 ± 0.04 μM for MLS000327186 and 0.87 ± 0.06 μM for MLS000327206 and greater than 50-fold selectivity versus h5-LOX, h12-LOX, h15-LOX-1, COX-1 and COX-2. A small set of focused analogs was synthesized to demonstrate the validity of the hits. In addition, a binding model was developed for the three imidazole inhibitors based on computational docking and a co-structure of h15-LOX-2 with MLS000536924. Hydrogen/deuterium exchange (HDX) results indicate a similar binding mode between MLS000536924 and MLS000327069, however, the latter restricts protein motion of helix-α2 more, consistent with its greater potency. Given these results, we designed, docked, and synthesized novel inhibitors of the imidazole scaffold and confirmed our binding mode hypothesis. Importantly, four of the five inhibitors mentioned above are active in an h15-LOX-2/HEK293 cell assay and thus they could be important tool compounds in gaining a better understanding of h15-LOX-2's role in human biology. As such, a suite of similar pharmacophores that target h15-LOX-2 both in vitro and ex vivo are presented in the hope of developing them as therapeutic agents.

The nuclear receptor-related 1 protein, Nurr1, is a transcription factor critical for the development and maintenance of dopamine-producing neurons in the substantia nigra pars compacta, a cell population that progressively loses the ability to make dopamine and degenerates in Parkinson's disease. Recently, we demonstrated that Nurr1 binds directly to and is regulated by the endogenous dopamine metabolite 5,6-dihydroxyindole (DHI). Unfortunately, DHI is an unstable compound, and thus a poor tool for studying Nurr1 function. Here, we report that 5-chloroindole, an unreactive analog of DHI, binds directly to the Nurr1 ligand binding domain with micromolar affinity and stimulates the activity of Nurr1, including the transcription of genes governing the synthesis and packaging of dopamine.

In this paper, human platelet 12-lipoxygenase [h12-LOX (ALOX12)], human reticulocyte 15-lipoxygenase-1 [h15-LOX-1 (ALOX15)], and human epithelial 15-lipoxygenase-2 [h15-LOX-2 (ALOX15B)] were observed to react with docosahexaenoic acid (DHA) and produce 17S-hydroperoxy-4Z,7Z,10Z,13Z,15E,19Z-docosahexaenoic acid (17S-HpDHA). The kcat/KM values with DHA for h12-LOX, h15-LOX-1, and h15-LOX-2 were 12, 0.35, and 0.43 s-1 μM-1, respectively, which demonstrate h12-LOX as the most efficient of the three. These values are comparable to their counterpart kcat/KM values with arachidonic acid (AA), 14, 0.98, and 0.24 s-1 μM-1, respectively. Comparison of their product profiles with DHA demonstrates that the three LOX isozymes produce 11S-HpDHA, 14S-HpDHA, and 17S-HpDHA, to varying degrees, with 17S-HpDHA being the majority product only for the 15-LOX isozymes. The effective kcat/KM values (kcat/KM × percent product formation) for 17S-HpDHA of the three isozymes indicate that the in vitro value of h12-LOX was 2.8-fold greater than that of h15-LOX-1 and 1.3-fold greater than that of h15-LOX-2. 17S-HpDHA was an effective substrate for h12-LOX and h15-LOX-1, with four products being observed under reducing conditions: protectin DX (PDX), 16S,17S-epoxy-4Z,7Z,10Z,12E,14E,19Z-docosahexaenoic acid (16S,17S-epoxyDHA), the key intermediate in neuroprotection D1 biosynthesis [NPD1, also known as protectin D1 (PD1)], 11,17S-diHDHA, and 16,17S-diHDHA. However, h15-LOX-2 did not react with 17-HpDHA. With respect to their effective kcat/KM values, h12-LOX was markedly less effective than h15-LOX-1 in reacting with 17S-HpDHA, with a 55-fold lower effective kcat/KM in producing 16S,17S-epoxyDHA and a 27-fold lower effective kcat/KM in generating PDX. This is the first direct demonstration of h15-LOX-1 catalyzing this reaction and reveals an in vitro pathway for PDX and NPD1 intermediate biosynthesis. In addition, epoxide formation from 17S-HpDHA and h15-LOX-1 was negatively affected via allosteric regulation by 17S-HpDHA (Kd = 5.9 μM), 12S-hydroxy-5Z,8Z,10E,14Z-eicosatetraenoic acid (12S-HETE) (Kd = 2.5 μM), and 17S-hydroxy-13Z,15E,19Z-docosatrienoic acid (17S-HDTA) (Kd = 1.4 μM), suggesting a possible regulatory pathway in reducing epoxide formation. Finally, 17S-HpDHA and PDX inhibited platelet aggregation, with EC50 values of approximately 1 and 3 μM, respectively. The in vitro results presented here may help advise in vivo PDX and NPD1 intermediate (i.e., 16S,17S-epoxyDHA) biosynthetic investigations and support the benefits of DHA rich diets.

Despite a remarkable increase in the genomic profiling of cancer, integration of genomic discoveries into clinical care has lagged behind. We report the feasibility of rapid identification of targetable mutations in 153 pediatric patients with relapsed/refractory or high-risk leukemias enrolled on a prospective clinical trial conducted by the LEAP Consortium. Eighteen percent of patients had a high confidence Tier 1 or 2 recommendation. We describe clinical responses in the 14% of patients with relapsed/refractory leukemia who received the matched targeted therapy. Further, in order to inform future targeted therapy for patients, we validated variants of uncertain significance, performed ex vivo drug-sensitivity testing in patient leukemia samples, and identified new combinations of targeted therapies in cell lines and patient-derived xenograft models. These data and our collaborative approach should inform the design of future precision medicine trials. SIGNIFICANCE: Patients with relapsed/refractory leukemias face limited treatment options. Systematic integration of precision medicine efforts can inform therapy. We report the feasibility of identifying targetable mutations in children with leukemia and describe correlative biology studies validating therapeutic hypotheses and novel mutations.See related commentary by Bornhauser and Bourquin, p. 1322.This article is highlighted in the In This Issue feature, p. 1307.

As organelles of the innate immune system, inflammasomes activate caspase-1 and other inflammatory caspases that cleave gasdermin D (GSDMD). Caspase-1 also cleaves inactive precursors of the interleukin (IL)-1 family to generate mature cytokines such as IL-1β and IL-18. Cleaved GSDMD forms transmembrane pores to enable the release of IL-1 and to drive cell lysis through pyroptosis1-9. Here we report cryo-electron microscopy structures of the pore and the prepore of GSDMD. These structures reveal the different conformations of the two states, as well as extensive membrane-binding elements including a hydrophobic anchor and three positively charged patches. The GSDMD pore conduit is predominantly negatively charged. By contrast, IL-1 precursors have an acidic domain that is proteolytically removed by caspase-110. When permeabilized by GSDMD pores, unlysed liposomes release positively charged and neutral cargoes faster than negatively charged cargoes of similar sizes, and the pores favour the passage of IL-1β and IL-18 over that of their precursors. Consistent with these findings, living-but not pyroptotic-macrophages preferentially release mature IL-1β upon perforation by GSDMD. Mutation of the acidic residues of GSDMD compromises this preference, hindering intracellular retention of the precursor and secretion of the mature cytokine. The GSDMD pore therefore mediates IL-1 release by electrostatic filtering, which suggests the importance of charge in addition to size in the transport of cargoes across this large channel.

Noninvasive methods to study glucocorticoid receptor (GR) signaling are urgently needed to elaborate the complexity of GR signaling in normal physiology and human disorders, as well as to identify selective GR modulators to treat diseases. Here, we report evidence supporting translational studies with (±)-11C-5-(4-fluorobenzyl)-10-methoxy-2,2,4-trimethyl-2,5-dihydro-1H-chromeno[3,4-f]-quinoline ((±)-11C-YJH08), a radioligand for PET that engages the ligand binding domain on GR. Methods: (±)-11C-YJH08 was synthesized by reacting the phenol precursor with 11C-methyl iodide. The biodistribution was studied in vivo. Specific binding was tested in vivo with adrenalectomy and ligand competition. A library of analogs was synthesized and studied in vitro and in vivo to understand the (±)-11C-YJH08 structure-activity relationship. Rodent dosimetry studies were performed to estimate the human-equivalent doses of (±)-11C-YJH08. Results: (±)-11C-YJH08 was synthesized by reaction of the phenolic precursor with 11C-methyl iodide, giving a radiochemical yield of 51.7% ± 4.7% (decay-corrected to starting 11C-methyl iodide). Specific binding was observed in many tissues, including the brain. An analysis of the (±)-YJH08 structure-activity relationship showed that (R)- and (S)-enantiomers are equally avid for GR by occupying discrete binding modes. A focused chemical screen revealed that the aryl fluoride motif on YJH08 is essential for high-affinity GR binding in vitro, high tissue uptake in vivo, and efficient passage across the blood-brain barrier. Lastly, we performed dosimetry studies on rodents, from which we estimated the human-equivalent doses of (±)-11C-YJH08 to be commensurate with the widely used 11C and 18F tracers. Conclusion: These studies reveal the molecular determinants of a high-affinity and high-selectivity ligand-receptor interaction and support the use of (±)-11C-YJH08 PET to make the first measurements of GR expression in human subjects.

The unfolded protein response (UPR) maintains protein folding homeostasis in the endoplasmic reticulum (ER). In metazoan cells, the Ire1 branch of the UPR initiates two functional outputs-non-conventional mRNA splicing and selective mRNA decay (RIDD). By contrast, Ire1 orthologs from Saccharomyces cerevisiae and Schizosaccharomyces pombe are specialized for only splicing or RIDD, respectively. Previously, we showed that the functional specialization lies in Ire1's RNase activity, which is either stringently splice-site specific or promiscuous (Li et al., 2018). Here, we developed an assay that reports on Ire1's RNase promiscuity. We found that conversion of two amino acids within the RNase domain of S. cerevisiae Ire1 to their S. pombe counterparts rendered it promiscuous. Using biochemical assays and computational modeling, we show that the mutations rewired a pair of salt bridges at Ire1 RNase domain's dimer interface, changing its protomer alignment. Thus, Ire1 protomer alignment affects its substrates specificity.

Small molecule colloidal aggregates adsorb and partially denature proteins, inhibiting them artifactually. Oddly, this inhibition is typically time-dependent. Two mechanisms might explain this: low concentrations of the colloid and enzyme might mean low encounter rates, or colloid-based protein denaturation might impose a kinetic barrier. These two mechanisms should have different concentration dependencies. Perplexingly, when enzyme concentration was increased, incubation times actually lengthened, inconsistent with both models and with classical chemical kinetics of solution species. We therefore considered molecular crowding, where colloids with lower protein surface density demand a shorter incubation time than more crowded colloids. To test this, we grew and shrank colloid surface area. As the surface area shrank, the incubation time lengthened, while as it increased, the converse was true. These observations support a crowding effect on protein binding to colloidal aggregates. Implications for drug delivery and for detecting aggregation-based inhibition will be discussed.

It was previously shown that human platelet 12S-lipoxygenase (h12-LOX) exists as a dimer; however, the specific structure is unknown. In this study, we create a model of the dimer through a combination of computational methods, experimental mutagenesis, and hydrogen-deuterium exchange (HDX) investigations. Initially, Leu183 and Leu187 were replaced by negatively charged glutamate residues and neighboring aromatic residues were replaced with alanine residues (F174A/W176A/L183E/L187E/Y191A). This quintuple mutant disrupted both the hydrophobic and π-π interactions, generating an h12-LOX monomer. To refine the determinants for dimer formation further, the L183E/L187E mutant was generated and the equilibrium shifted mostly toward the monomer. We then submitted the predicted monomeric structure to protein-protein docking to create a model of the dimeric complex. A total of nine of the top 10 most energetically favorable docking conformations predict a TOP-to-TOP dimeric arrangement of h12-LOX, with the α-helices containing a Leu-rich region (L172, L183, L187, and L194), corroborating our experimental results showing the importance of these hydrophobic interactions for dimerization. This model was supported by HDX investigations that demonstrated the stabilization of four, non-overlapping peptides within helix α2 of the TOP subdomain for wt-h12-LOX, consistent with the dimer interface. Most importantly, our data reveal that the dimer and monomer of h12-LOX have distinct biochemical properties, suggesting that the structural changes due to dimerization have allosteric effects on active site catalysis and inhibitor binding.

Glycosaminoglycans (GAGs) such as heparan sulfate and chondroitin sulfate decorate all mammalian cell surfaces. These mucopolysaccharides act as coreceptors for extracellular ligands, regulating cell signaling, growth, proliferation, and adhesion. In glioblastoma, the most common type of primary malignant brain tumor, dysregulated GAG biosynthesis results in altered chain length, sulfation patterns, and the ratio of contributing monosaccharides. These events contribute to the loss of normal cellular function, initiating and sustaining malignant growth. Disruption of the aberrant cell surface GAGs with small molecule inhibitors of GAG biosynthetic enzymes is a potential therapeutic approach to blocking the rogue signaling and proliferation in glioma, including glioblastoma. Previously, 4-azido-xylose-α-UDP sugar inhibited both xylosyltransferase (XYLT-1) and β-1,4-galactosyltransferase-7 (β-GALT-7)-the first and second enzymes of GAG biosynthesis-when microinjected into a cell. In another study, 4-deoxy-4-fluoro-β-xylosides inhibited β-GALT-7 at 1 mM concentration in vitro. In this work, we seek to solve the enduring problem of drug delivery to human glioma cells at low concentrations. We developed a library of hydrophobic, presumed prodrugs 4-deoxy-4-fluoro-2,3-dibenzoyl-(α- or β-) xylosides and their corresponding hydrophilic inhibitors of XYLT-1 and β-GALT-7 enzymes. The prodrugs were designed to be activatable by carboxylesterase enzymes overexpressed in glioblastoma. Using a colorimetric MTT assay in human glioblastoma cell lines, we identified a prodrug-drug pair (4-nitrophenyl-α-xylosides) as lead drug candidates. The candidates arrest U251 cell growth at an IC50 = 380 nM (prodrug), 122 μM (drug), and U87 cells at IC50 = 10.57 μM (prodrug). Molecular docking studies were consistent with preferred binding of the α- versus β-nitro xyloside conformer to XYLT-1 and β-GALT-7 enzymes.

The oxylipins, 5S,12S-dihydroxy-6E,8Z,10E,14Z-eicosatetraenoic acid (5S,12S-diHETE) and 5S,15S-dihydroxy-6E,8Z,11Z,13E-eicosatetraenoic acid (5S,15S-diHETE), have been identified in cell exudates and have chemotactic activity toward eosinophils and neutrophils. Their biosynthesis has been proposed to occur by sequential oxidations of arachidonic acid (AA) by lipoxygenase enzymes, specifically through oxidation of AA by h5-LOX followed by h12-LOX, h15-LOX-1, or h15-LOX-2. In this work, h15-LOX-1 demonstrates altered positional specificity when reacting with 5S-HETE, producing 90% 5S,12S-diHETE, instead of 5S,15S-diHETE, with kinetics 5-fold greater than that of h12-LOX. This is consistent with previous work in which h15-LOX-1 reacts with 7S-HDHA, producing the noncanonical, DHA-derived, specialized pro-resolving mediator, 7S,14S-diHDHA. It is also determined that oxygenation of 5S-HETE by h15-LOX-2 produces 5S,15S-diHETE and its biosynthetic kcat/KM flux is 2-fold greater than that of h15-LOX-1, suggesting that h15-LOX-2 may have a greater role in lipoxin biosynthesis than previously thought. In addition, it is shown that oxygenation of 12S-HETE and 15S-HETE by h5-LOX is kinetically slow, suggesting that the first step in the in vitro biosynthesis of both 5S,12S-diHETE and 5S,15S-diHETE is the production of 5S-HETE.

Natural products serve as chemical blueprints for most antibiotics in clinical use. The evolutionary process by which these molecules arise is inherently accompanied by the co-evolution of resistance mechanisms that shorten the clinical lifetime of any given class of antibiotics1. Virginiamycin acetyltransferase (Vat) enzymes are resistance proteins that provide protection against streptogramins2, potent antibiotics against Gram-positive bacteria that inhibit the bacterial ribosome3. Owing to the challenge of selectively modifying the chemically complex, 23-membered macrocyclic scaffold of group A streptogramins, analogues that overcome the resistance conferred by Vat enzymes have not been previously developed2. Here we report the design, synthesis, and antibacterial evaluation of group A streptogramin antibiotics with extensive structural variability. Using cryo-electron microscopy and forcefield-based refinement, we characterize the binding of eight analogues to the bacterial ribosome at high resolution, revealing binding interactions that extend into the peptidyl tRNA-binding site and towards synergistic binders that occupy the nascent peptide exit tunnel. One of these analogues has excellent activity against several streptogramin-resistant strains of Staphylococcus aureus, exhibits decreased rates of acetylation in vitro, and is effective at lowering bacterial load in a mouse model of infection. Our results demonstrate that the combination of rational design and modular chemical synthesis can revitalize classes of antibiotics that are limited by naturally arising resistance mechanisms.

The two oxylipins 7S,14S-dihydroxydocosahexaenoic acid (diHDHA) and 7S,17S-diHDHA [resolvin D5 (RvD5)] have been found in macrophages and infectious inflammatory exudates and are believed to function as specialized pro-resolving mediators (SPMs). Their biosynthesis is thought to proceed through sequential oxidations of DHA by lipoxygenase (LOX) enzymes, specifically, by human 5-LOX (h5-LOX) first to 7(S)-hydroxy-4Z,8E,10Z,13Z,16Z,19Z-DHA (7S-HDHA), followed by human platelet 12-LOX (h12-LOX) to form 7(S),14(S)-dihydroxy-4Z,8E,10Z,12E,16Z,19Z-DHA (7S,14S-diHDHA) or human reticulocyte 15-LOX-1 (h15-LOX-1) to form RvD5. In this work, we determined that oxidation of 7(S)-hydroperoxy-4Z,8E,10Z,13Z,16Z,19Z-DHA to 7S,14S-diHDHA is performed with similar kinetics by either h12-LOX or h15-LOX-1. The oxidation at C14 of DHA by h12-LOX was expected, but the noncanonical reaction of h15-LOX-1 to make over 80% 7S,14S-diHDHA was larger than expected. Results of computer modeling suggested that the alcohol on C7 of 7S-HDHA hydrogen bonds with the backbone carbonyl of Ile399, forcing the hydrogen abstraction from C12 to oxygenate on C14 but not C17. This result raised questions regarding the synthesis of RvD5. Strikingly, we found that h15-LOX-2 oxygenates 7S-HDHA almost exclusively at C17, forming RvD5 with faster kinetics than does h15-LOX-1. The presence of h15-LOX-2 in neutrophils and macrophages suggests that it may have a greater role in biosynthesizing SPMs than previously thought. We also determined that the reactions of h5-LOX with 14(S)-hydroperoxy-4Z,7Z,10Z,12E,16Z,19Z-DHA and 17(S)-hydroperoxy-4Z,7Z,10Z,13Z,15E,19Z-DHA are kinetically slow compared with DHA, suggesting that these reactions may be minor biosynthetic routes in vivo. Additionally, we show that 7S,14S-diHDHA and RvD5 have anti-aggregation properties with platelets at low micromolar potencies, which could directly regulate clot resolution.

Human reticulocyte 15-lipoxygenase-1 (h15-LOX-1 or ALOX15) and platelet 12-lipoxygenase (h12-LOX or ALOX12) catalysis of docosahexaenoic acid (DHA) and the maresin precursor, 14S-hydroperoxy-4Z,7Z,10Z,12E,16Z,19Z-docosahexaenoic acid (14S-HpDHA), were investigated to determine their product profiles and relative rates in the biosynthesis of the key maresin intermediate, 13S,14S-epoxy-4Z,7Z,9E,11E,16Z,19Z-docosahexaenoic acid (13S,14S-epoxy-DHA). Both enzymes converted DHA to 14S-HpDHA, with h12-LOX having a 39-fold greater kcat/KM value (14.0 ± 0.8 s-1 μM-1) than that of h15-LOX-1 (0.36 ± 0.08 s-1 μM-1) and a 1.8-fold greater 14S-HpDHA product selectivity, 81 and 46%, respectively. However, h12-LOX was markedly less effective at producing 13S,14S-epoxy-DHA from 14S-HpDHA than h15-LOX-1, with a 4.6-fold smaller kcat/KM value, 0.0024 ± 0.0002 and 0.11 ± 0.006 s-1 μM-1, respectively. This is the first evidence of h15-LOX-1 to catalyze this reaction and reveals a novel in vitro pathway for maresin biosynthesis. In addition, epoxidation of 14S-HpDHA is negatively regulated through allosteric oxylipin binding to h15-LOX-1 and h12-LOX. For h15-LOX-1, 14S-HpDHA (Kd = 6.0 μM), 12S-hydroxy-5Z,8Z,10E,14Z-eicosatetraenoic acid (12S-HETE) (Kd = 3.5 μM), and 14S-hydroxy-7Z,10Z,12E,16Z,19Z-docosapentaenoic acid (14S-HDPAω-3) (Kd = 4.0 μM) were shown to decrease 13S,14S-epoxy-DHA production. h12-LOX was also shown to be allosterically regulated by 14S-HpDHA (Kd = 3.5 μM) and 14S-HDPAω-3 (Kd = 4.0 μM); however, 12S-HETE showed no effect, indicating for the first time an allosteric response by h12-LOX. Finally, 14S-HpDHA inhibited platelet aggregation at a submicrololar concentration, which may have implications in the benefits of diets rich in DHA. These in vitro biosynthetic pathways may help guide in vivo maresin biosynthetic investigations and possibly direct therapeutic interventions.

KIT is a type-3 receptor tyrosine kinase that is frequently mutated at exon 11 or 17 in a variety of cancers. First-generation KIT tyrosine kinase inhibitors (TKI) are ineffective against KIT exon 17 mutations, which favor an active conformation that prevents these TKIs from binding. The ATP-competitive inhibitors, midostaurin and avapritinib, which target the active kinase conformation, were developed to inhibit exon 17-mutant KIT. Because secondary kinase domain mutations are a common mechanism of TKI resistance and guide ensuing TKI design, we sought to define problematic KIT kinase domain mutations for these emerging therapeutics. Midostaurin and avapritinib displayed different vulnerabilities to secondary kinase domain substitutions, with the T670I gatekeeper mutation being selectively problematic for avapritinib. Although gatekeeper mutations often directly disrupt inhibitor binding, we provide evidence that T670I confers avapritinib resistance indirectly by inducing distant conformational changes in the phosphate-binding loop. These findings suggest combining midostaurin and avapritinib may forestall acquired resistance mediated by secondary kinase domain mutations. SIGNIFICANCE: This study identifies potential problematic kinase domain mutations for next-generation KIT inhibitors midostaurin and avapritinib.

The progressive failure of protein homeostasis is a hallmark of aging and a common feature in neurodegenerative disease. As the enzymes executing the final stages of autophagy, lysosomal proteases are key contributors to the maintenance of protein homeostasis with age. We previously reported that expression of granulin peptides, the cleavage products of the neurodegenerative disease protein progranulin, enhance the accumulation and toxicity of TAR DNA binding protein 43 (TDP-43) in Caenorhabditis elegans (C. elegans). In this study we show that C. elegans granulins are produced in an age- and stress-dependent manner. Granulins localize to the endolysosomal compartment where they impair lysosomal protease expression and activity. Consequently, protein homeostasis is disrupted, promoting the nuclear translocation of the lysosomal transcription factor HLH-30/TFEB, and prompting cells to activate a compensatory transcriptional program. The three C. elegans granulin peptides exhibited distinct but overlapping functional effects in our assays, which may be due to amino acid composition that results in distinct electrostatic and hydrophobicity profiles. Our results support a model in which granulin production modulates a critical transition between the normal, physiological regulation of protease activity and the impairment of lysosomal function that can occur with age and disease.

Progranulin (PGRN) is an evolutionarily conserved glycoprotein associated with several disease states, including neurodegeneration, cancer, and autoimmune disorders. This protein has recently been implicated in the regulation of lysosome function, whereby PGRN may bind to and promote the maturation and activity of the aspartyl protease cathepsin D (proCTSD, inactive precursor; matCTSD, mature, enzymatically active form). As the full-length PGRN protein can be cleaved into smaller peptides, called granulins, we assessed the function of these granulin peptides in binding to proCTSD and stimulating matCTSD enzyme activity in vitro. Here, we report that full-length PGRN and multi-granulin domain peptides bound to proCTSD with low to submicromolar binding affinities. This binding promoted proCTSD destabilization, the magnitude of which was greater for multi-granulin domain peptides than for full-length PGRN. Such destabilization correlated with enhanced matCTSD activity at acidic pH. The presence and function of multi-granulin domain peptides have typically been overlooked in previous studies. This work provides the first in vitro quantification of their binding and activity on proCTSD. Our study highlights the significance of multi-granulin domain peptides in the regulation of proCTSD maturation and enzymatic activity and suggests that attention to PGRN processing will be essential for the future understanding of the molecular mechanisms leading to neurodegenerative disease states with loss-of-function mutations in PGRN.

Post-translational modifications of the RNA polymerase II C-terminal domain (CTD) coordinate the transcription cycle. Crosstalk between different modifications is poorly understood. Here, we show how acetylation of lysine residues at position 7 of characteristic heptad repeats (K7ac)-only found in higher eukaryotes-regulates phosphorylation of serines at position 5 (S5p), a conserved mark of polymerases initiating transcription. We identified the regulator of pre-mRNA-domain-containing (RPRD) proteins as reader proteins of K7ac. K7ac enhanced CTD peptide binding to the CTD-interacting domain (CID) of RPRD1A and RPRD1B proteins in isothermal calorimetry and molecular modeling experiments. Deacetylase inhibitors increased K7ac- and decreased S5-phosphorylated polymerases, consistent with acetylation-dependent S5 dephosphorylation by an RPRD-associated S5 phosphatase. Consistent with this model, RPRD1B knockdown increased S5p but enhanced K7ac, indicating that RPRD proteins recruit K7 deacetylases, including HDAC1. We also report autoregulatory crosstalk between K7ac and S5p via RPRD proteins and their interactions with acetyl- and phospho-eraser proteins.

Nurr1, a nuclear receptor essential for the development, maintenance, and survival of midbrain dopaminergic neurons, is a potential therapeutic target for Parkinson's disease, a neurological disorder characterized by the degeneration of these same neurons. Efforts to identify Nurr1 agonists have been hampered by the recognition that it lacks several classic regulatory elements of nuclear receptor function, including the canonical ligand-binding pocket. Here we report that the dopamine metabolite 5,6-dihydroxyindole (DHI) binds directly to and modulates the activity of Nurr1. Using biophysical assays and X-ray crystallography, we show that DHI binds to the ligand-binding domain within a non-canonical pocket, forming a covalent adduct with Cys566. In cultured cells and zebrafish, DHI stimulates Nurr1 activity, including the transcription of target genes underlying dopamine homeostasis. These findings suggest avenues for developing synthetic Nurr1 ligands to ameliorate the symptoms and progression of Parkinson's disease.

Tau, a member of the MAP2/tau family of microtubule-associated proteins, stabilizes and organizes axonal microtubules in healthy neurons. In neurodegenerative tauopathies, tau dissociates from microtubules and forms neurotoxic extracellular aggregates. MAP2/tau family proteins are characterized by three to five conserved, intrinsically disordered repeat regions that mediate electrostatic interactions with the microtubule surface. Here, we used molecular dynamics, microtubule-binding experiments, and live-cell microscopy, revealing that highly-conserved histidine residues near the C terminus of each microtubule-binding repeat are pH sensors that can modulate tau-microtubule interaction strength within the physiological intracellular pH range. We observed that at low pH (7.5), tau deprotonation decreased binding to microtubules both in vitro and in cells. Electrostatic and hydrophobic characteristics of histidine were both required for tau-microtubule binding, as substitutions with constitutively and positively charged nonaromatic lysine or uncharged alanine greatly reduced or abolished tau-microtubule binding. Consistent with these findings, tau-microtubule binding was reduced in a cancer cell model with increased intracellular pH but was rapidly restored by decreasing the pH to normal levels. These results add detailed insights into the intracellular regulation of tau activity that may be relevant in both normal and pathological conditions.

Human immunodeficiency virus-1 viral infectivity factor (Vif) is an intrinsically disordered protein responsible for the ubiquitination of the APOBEC3 (A3) antiviral proteins. Vif folds when it binds Cullin-RING E3 ligase 5 and the transcription cofactor CBF-β. A five-protein complex containing the substrate receptor (Vif, CBF-β, Elongin-B, Elongin-C (VCBC)) and Cullin5 (CUL5) has a published crystal structure, but dynamics of this VCBC-CUL5 complex have not been characterized. Here, we use molecular dynamics (MD) simulations and NMR to characterize the dynamics of the VCBC complex with and without CUL5 and an A3 protein bound. Our simulations show that the VCBC complex undergoes global dynamics involving twisting and clamshell opening of the complex, whereas VCBC-CUL5 maintains a more static conformation, similar to the crystal structure. This observation from MD is supported by methyl-transverse relaxation-optimized spectroscopy NMR data, which indicates that the VCBC complex without CUL5 is dynamic on the μs-ms timescale. Our NMR data also show that the VCBC complex is more conformationally restricted when bound to the antiviral APOBEC3F (one of the A3 proteins), consistent with our MD simulations. Vif contains a flexible linker region located at the hinge of the VCBC complex, which changes conformation in conjunction with the global dynamics of the complex. Like other substrate receptors, VCBC can exist alone or in complex with CUL5 and other proteins in cells. Accordingly, the VCBC complex could be a good target for therapeutics that would inhibit full assembly of the ubiquitination complex by stabilizing an alternate VCBC conformation.

Single-copy loss-of-function mutations in the progranulin gene (PGRN) underlie the neurodegenerative disease frontotemporal lobar degeneration, while homozygous loss-of-function of PGRN results in the lysosomal storage disorder neuronal ceroid lipofuscinosis. Despite evidence that normal PGRN levels are critical for neuronal health, the function of this protein is not yet understood. Here, we show that PGRN stimulates the in vitro maturation of the lysosomal aspartyl protease cathepsin D (CTSD). CTSD is delivered to the endolysosomal system as an inactive precursor (proCTSD) and requires sequential cleavage steps via intermediate forms to achieve the mature state (matCTSD). In co-immunoprecipitation experiments, PGRN interacts predominantly with immature pro- and intermediate forms of CTSD. PGRN enhances in vitro conversion of proCTSD to matCTSD in a concentration-dependent manner. Differential scanning fluorimetry shows a destabilizing effect induced by PGRN on proCTSD folding (∆Tm = -1.7 °C at a 3:1 molar ratio). We propose a mechanism whereby PGRN binds to proCTSD, destabilizing the propeptide from the enzyme catalytic core and favoring conversion to mature forms of the enzyme. Further understanding of the role of PGRN in CTSD maturation will assist in the development of targeted therapies for neurodegenerative disease.

Human platelet ALOX12 (hALOX12 or h12-LOX) has been implicated in a variety of human diseases. The present study investigates the active site of hALOX12 to more thoroughly understand how it positions the substrate and achieves nearly perfect regio- and stereospecificities (i.e., 100 ± 5% of the 12(S)-hydroperoxide product), utilizing site-directed mutagenesis. Specifically, we have determined that Arg402 is not as important in substrate binding as previously seen for hALOX15 but that His596 may play a role in anchoring the carboxy terminal of the arachidonic acid during catalysis. In addition, Phe414 creates a π-stacking interaction with a double bond of arachidonic acid (Δ11), and Ala417/Val418 define the bottom of the cavity. However, the influence of Ala417/Val418 on the profile is markedly less for hALOX12 than that seen in hALOX15. Mutating these two residues to larger amino acids (Ala417Ile/Val418Met) only increased the generation of 15-HpETE by 24 ± 2%, but conversely, smaller residues at these positions converted hALOX15 to almost 100% hALOX12 reactivity [Gan et al. (1996) J. Biol. Chem. 271, 25412-25418]. However, we were able to increase 15-HpETE to 46 ± 3% by restricting the width of the active site with the Ala417Ile/Val418Met/Ser594Thr mutation, indicating both depth and width of the active site are important. Finally, residue Leu407 is shown to play a critical role in positioning the substrate correctly, as seen by the increase of 15-HpETE to 21 ± 1% for the single Leu407Gly mutant. These results outline critical differences between the active site requirements of hALOX12 relative to hALOX15 and explain both their product specificity and inhibitory differences.

Sodium-dependent glucose transporters (SGLTs) exploit sodium gradients to transport sugars across the plasma membrane. Due to their role in renal sugar reabsorption, SGLTs are targets for the treatment of type 2 diabetes. Current therapeutics are phlorizin derivatives that contain a sugar moiety bound to an aromatic aglycon tail. Here, we develop structural models of human SGLT1/2 in complex with inhibitors by combining computational and functional studies. Inhibitors bind with the sugar moiety in the sugar pocket and the aglycon tail in the extracellular vestibule. The binding poses corroborate mutagenesis studies and suggest a partial closure of the outer gate upon binding. The models also reveal a putative Na+ binding site in hSGLT1 whose disruption reduces the transport stoichiometry to the value observed in hSGLT2 and increases inhibition by aglycon tails. Our work demonstrates that subtype selectivity arises from Na+-regulated outer gate closure and a variable region in extracellular loop EL5.

Cancer cells rely on the chaperone heat shock protein 70 (Hsp70) for survival and proliferation. Recently, benzothiazole rhodacyanines have been shown to bind an allosteric site on Hsp70, interrupting its binding to nucleotide-exchange factors (NEFs) and promoting cell death in breast cancer cell lines. However, proof-of-concept molecules, such as JG-98, have relatively modest potency (EC50 ≈ 0.7-0.4 μM) and are rapidly metabolized in animals. Here, we explored this chemical series through structure- and property-based design of ∼300 analogs, showing that the most potent had >10-fold improved EC50 values (∼0.05 to 0.03 μM) against two breast cancer cells. Biomarkers and whole genome CRISPRi screens confirmed members of the Hsp70 family as cellular targets. On the basis of these results, JG-231 was found to reduce tumor burden in an MDA-MB-231 xenograft model (4 mg/kg, ip). Together, these studies support the hypothesis that Hsp70 may be a promising target for anticancer therapeutics.

Protein-protein interactions (PPIs) are an important category of putative drug targets. Improvements in high-throughput screening (HTS) have significantly accelerated the discovery of inhibitors for some categories of PPIs. However, methods suitable for screening multiprotein complexes (e.g. those composed of three or more different components) have been slower to emerge. Here, we explored an approach that uses reconstituted multiprotein complexes (RMPCs). As a model system, we chose heat shock protein 70 (Hsp70), which is an ATP-dependent molecular chaperone that interacts with co-chaperones, including DnaJA2 and BAG2. The PPIs between Hsp70 and its co-chaperones stimulate nucleotide cycling. Thus, to re-create this ternary protein system, we combined purified human Hsp70 with DnaJA2 and BAG2 and then screened 100,000 diverse compounds for those that inhibited co-chaperone-stimulated ATPase activity. This HTS campaign yielded two compounds with promising inhibitory activity. Interestingly, one inhibited the PPI between Hsp70 and DnaJA2, whereas the other seemed to inhibit the Hsp70-BAG2 complex. Using secondary assays, we found that both compounds inhibited the PPIs through binding to allosteric sites on Hsp70, but neither affected Hsp70's intrinsic ATPase activity. Our RMPC approach expands the toolbox of biochemical HTS methods available for studying difficult-to-target PPIs in multiprotein complexes. The results may also provide a starting point for new chemical probes of the Hsp70 system.

Chordoid glioma is a rare brain tumor thought to arise from specialized glial cells of the lamina terminalis along the anterior wall of the third ventricle. Despite being histologically low-grade, chordoid gliomas are often associated with poor outcome, as their stereotypic location in the third ventricle makes resection challenging and efficacious adjuvant therapies have not been developed. Here we performed genomic profiling on 13 chordoid gliomas and identified a recurrent D463H missense mutation in PRKCA in all tumors, which localizes in the kinase domain of the encoded protein kinase C alpha (PKCα). Expression of mutant PRKCA in immortalized human astrocytes led to increased phospho-ERK and anchorage-independent growth that could be blocked by MEK inhibition. These studies define PRKCA as a recurrently mutated oncogene in human cancer and identify a potential therapeutic vulnerability in this uncommon brain tumor.

The functions of most proteins are yet to be determined. The function of an enzyme is often defined by its interacting partners, including its substrate and product, and its role in larger metabolic networks. Here, we describe a computational method that predicts the functions of orphan enzymes by organizing them into a linear metabolic pathway. Given candidate enzyme and metabolite pathway members, this aim is achieved by finding those pathways that satisfy structural and network restraints implied by varied input information, including that from virtual screening, chemoinformatics, genomic context analysis, and ligand -binding experiments. We demonstrate this integrative pathway mapping method by predicting the L-gulonate catabolic pathway in Haemophilus influenzae Rd KW20. The prediction was subsequently validated experimentally by enzymology, crystallography, and metabolomics. Integrative pathway mapping by satisfaction of structural and network restraints is extensible to molecular networks in general and thus formally bridges the gap between structural biology and systems biology.

The functions of most proteins are yet to be determined. The function of an enzyme is often defined by its interacting partners, including its substrate and product, and its role in larger metabolic networks. Here, we describe a computational method that predicts the functions of orphan enzymes by organizing them into a linear metabolic pathway. Given candidate enzyme and metabolite pathway members, this aim is achieved by finding those pathways that satisfy structural and network restraints implied by varied input information, including that from virtual screening, chemoinformatics, genomic context analysis, and ligand -binding experiments. We demonstrate this integrative pathway mapping method by predicting the L-gulonate catabolic pathway in Haemophilus influenzae Rd KW20. The prediction was subsequently validated experimentally by enzymology, crystallography, and metabolomics. Integrative pathway mapping by satisfaction of structural and network restraints is extensible to molecular networks in general and thus formally bridges the gap between structural biology and systems biology.

Central nervous system (CNS) chemical protection depends upon discrete control of small-molecule access by the blood-brain barrier (BBB). Curiously, some drugs cause CNS side-effects despite negligible transit past the BBB. To investigate this phenomenon, we asked whether the highly BBB-enriched drug efflux transporter MDR1 has dual functions in controlling drug and endogenous molecule CNS homeostasis. If this is true, then brain-impermeable drugs could induce behavioral changes by affecting brain levels of endogenous molecules. Using computational, genetic, and pharmacologic approaches across diverse organisms, we demonstrate that BBB-localized efflux transporters are critical for regulating brain levels of endogenous steroids and steroid-regulated behaviors (sleep in Drosophila and anxiety in mice). Furthermore, we show that MDR1-interacting drugs are associated with anxiety-related behaviors in humans. We propose a general mechanism for common behavioral side effects of prescription drugs: pharmacologically challenging BBB efflux transporters disrupts brain levels of endogenous substrates and implicates the BBB in behavioral regulation.

The intracellular pH (pHi) of most cancers is constitutively higher than that of normal cells and enhances proliferation and cell survival. We found that increased pHi enabled the tumorigenic behaviors caused by somatic arginine-to-histidine mutations, which are frequent in cancer and confer pH sensing not seen with wild-type proteins. Experimentally raising the pHi increased the activity of R776H mutant epidermal growth factor receptor (EGFR-R776H), thereby increasing proliferation and causing transformation in fibroblasts. An Arg-to-Gly mutation did not confer these effects. Molecular dynamics simulations of EGFR suggested that decreased protonation of His776 at high pH causes conformational changes in the αC helix that may stabilize the active form of the kinase. An Arg-to-His, but not Arg-to-Lys, mutation in the transcription factor p53 (p53-R273H) decreased its transcriptional activity and attenuated the DNA damage response in fibroblasts and breast cancer cells with high pHi. Lowering pHi attenuated the tumorigenic effects of both EGFR-R776H and p53-R273H. Our data suggest that some somatic mutations may confer a fitness advantage to the higher pHi of cancer cells.

Nicotinic acid adenine dinucleotide phosphate (NAADP) and cyclic ADP-ribose (cADPR) are Ca2+-mobilizing messengers important for modulating cardiac excitation-contraction coupling and pathophysiology. CD38, which belongs to the ADP-ribosyl cyclase family, catalyzes synthesis of both NAADP and cADPR in vitro However, it remains unclear whether this is the main enzyme for their production under physiological conditions. Here we show that membrane fractions from WT but not CD38-/- mouse hearts supported NAADP and cADPR synthesis. Membrane permeabilization of cardiac myocytes with saponin and/or Triton X-100 increased NAADP synthesis, indicating that intracellular CD38 contributes to NAADP production. The permeabilization also permitted immunostaining of CD38, with a striated pattern in WT myocytes, whereas CD38-/- myocytes and nonpermeabilized WT myocytes showed little or no staining, without striation. A component of β-adrenoreceptor signaling in the heart involves NAADP and lysosomes. Accordingly, in the presence of isoproterenol, Ca2+ transients and contraction amplitudes were smaller in CD38-/- myocytes than in the WT. In addition, suppressing lysosomal function with bafilomycin A1 reduced the isoproterenol-induced increase in Ca2+ transients in cardiac myocytes from WT but not CD38-/- mice. Whole hearts isolated from CD38-/- mice and exposed to isoproterenol showed reduced arrhythmias. SAN4825, an ADP-ribosyl cyclase inhibitor that reduces cADPR and NAADP synthesis in mouse membrane fractions, was shown to bind to CD38 in docking simulations and reduced the isoproterenol-induced arrhythmias in WT hearts. These observations support generation of NAADP and cADPR by intracellular CD38, which contributes to effects of β-adrenoreceptor stimulation to increase both Ca2+ transients and the tendency to disturb heart rhythm.

Macrocycles are appealing drug candidates due to their high affinity, specificity, and favorable pharmacological properties. In this study, we explored the effects of chemical modifications to a natural product macrocycle upon its activity, 3D geometry, and conformational entropy. We chose thiocillin as a model system, a thiopeptide in the ribosomally encoded family of natural products that exhibits potent antimicrobial effects against Gram-positive bacteria. Since thiocillin is derived from a genetically encoded peptide scaffold, site-directed mutagenesis allows for rapid generation of analogues. To understand thiocillin's structure-activity relationship, we generated a site-saturation mutagenesis library covering each position along thiocillin's macrocyclic ring. We report the identification of eight unique compounds more potent than wild-type thiocillin, the best having an 8-fold improvement in potency. Computational modeling of thiocillin's macrocyclic structure revealed a striking requirement for a low-entropy macrocycle for activity. The populated ensembles of the active mutants showed a rigid structure with few adoptable conformations while inactive mutants showed a more flexible macrocycle which is unfavorable for binding. This finding highlights the importance of macrocyclization in combination with rigidifying post-translational modifications to achieve high-potency binding.