Search

VAF concordance

Supplementary Data file including supplementary methods and clinical case vignettes.

Reasons for not using recommended targeted therapy.

Patient sequencing data, by patient number and diagnosis.

This document contains supplementary methods, tables, figures and references for the manuscript "ATP-competitive inhibitors midostaurin and avapritinib have distinct resistance profiles in exon 17-mutant KIT".

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 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.

CREBBP (CBP/KAT3A) and its paralogue EP300 (KAT3B) are lysine acetyltransferases (KATs) that are essential for human development. They each comprise 10 domains through which they interact with >400 proteins, making them important transcriptional co-activators and key nodes in the human protein-protein interactome. The bromodomains of CREBBP and EP300 enable the binding of acetylated lysine residues from histones and a number of other important proteins, including p53, p73, E2F, and GATA1. Here, we report a work to develop a high-affinity, small-molecule ligand for the CREBBP and EP300 bromodomains [(-)-OXFBD05] that shows >100-fold selectivity over a representative member of the BET bromodomains, BRD4(1). Cellular studies using this ligand demonstrate that the inhibition of the CREBBP/EP300 bromodomain in HCT116 colon cancer cells results in lowered levels of c-Myc and a reduction in H3K18 and H3K27 acetylation. In hypoxia (

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

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.

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.

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.

In spite of significant efforts to reduce malaria deaths, this disease still kills around 445,000 people every year. Overcoming drug resistance is one of the main goals of current malaria research programs. This is challenging, since the biology of Plasmodium is not fully understood, requiring the development of advanced models for data analysis in the search for new antimalarials. Areas covered: In this review the authors introduce the importance of computational models to address the challenges of drug discovery, presenting examples of pioneering systems biology approaches in the search for new antimalarial drugs and their role in the future of drug research programs. Other related topics are discussed, e.g. regulation of malaria pathogenesis by epigenetics and the importance of new platforms for malaria network. Expert opinion: The use of a systems biology approach in antimalarial drug discovery emerges in a scenario where the most efficient antimalarial chemotherapies are showing resistance in Southeast Asia. New models for a better understanding of Plasmodium cell function have already proved to be powerful tools for uncovering complex mechanisms of resistance, and have great potential to inform the design of novel small molecules with both high antimalarial activity and transmission-blocking potential to improve the control of malaria.

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., Author summary Progressive decline in maintenance of protein homeostasis clearly contributes to the development of neurodegenerative disorders, yet the molecular basis of this decline is poorly understood. Here, we take advantage of molecular genetic techniques available in the model organism C. elegans to investigate the mechanism underlying neurodegenerative disease due to mutations in the progranulin gene. We find that age, gene mutation and physiological stress lead to the accumulation of lysosomal granulins (the cleavage products of the progranulin protein) thereby disrupting cellular protein homeostasis. Granulin expression impairs animal fitness, resistance to stress and neuronal function, and stimulates a lysosomal stress response in an attempt to up-regulate lysosomal genes and restore normal function. Our findings are particularly important because they suggest a new, rational target—inhibition of progranulin cleavage into granulins—for neurodegenerative disease therapy.

Regulation of drug prices that ensures adequate access to effective treatments and promotes innovation is a global challenge. In the United States, the government does not regulate drug prices when they come onto market. On the other hand, in countries such as France and Brazil, government agencies are responsible for setting up price limits by leveraging the interests of the companies and the countries' population. In Brazil, safety and efficacy of drugs are regulated by the Brazilian Health Regulatory Agency, and drug prices are regulated by the Pharmaceutical Market Regulation Chamber with a participation of Brazilian Health Regulatory Agency. Here, we introduce the current challenges faced by the Brazilian government in the drug price regulation and present proposed initiatives aiming to streamline access to innovative treatments for its citizens.

Protein kinase CK2 has emerged as an attractive cancer therapeutic target. Previous studies have highlighted the challenge of optimizing CK2 ATP-competitive inhibitors that have low druggability due to their polycyclic ring scaffolds. Therefore the development of novel inhibitors with non-polycyclic scaffolds emerges as a promising strategy for drug discovery targeting CK2. In this current study, based on the similar predicted binding poses of the linear 2-propenone scaffold of isoliquiritigenin with that of the polycyclic inhibitor CX-4945, a series of 2-propenone derivatives containing an amine-substituted five-membered heterocycle and a benzoic acid were designed, synthesized and evaluated for their in vitro CK2 inhibition and anti-cancer activity. Compound 8b was found to be the most potent CK2 inhibitor (IC50 = 0.6 μM) with the anti-proliferative activity on HepG2 cancer cells (IC50 = 14 μM), compared to the activity of isoliquiritigenin (IC50 = 17 μM and 51 μM, respectively). Molecular docking was performed to understand the binding modes of the newly designed 2-propenone derivatives with CK2. Compound 8b formed the most favorable network of hydrogen bonds with both the hinge region and positive area. Our results indicate that CK2 derivatives with a linear 2-propenone scaffold are promising candidates for anti-cancer drug discovery.

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.

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.

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.

Protein kinase CK2 is considered as an emerging target in cancer therapy, and recent efforts have been made to develop its ATP-competitive inhibitors, but achieving selectivity with respect to related kinases remains challenging because of the highly conserved ATP-binding pocket of kinases. Non-ATP competitive inhibitors might solve this challenge; one such strategy is to identify compounds that target the CK2α/CK2β interface as CK2 holoenzyme antagonists. Here we improved the binding affinity to CK2α and cell-based anti-cancer activity of a CK2β-derived cyclic peptide (Pc) by combining structure-based computational design with experimental evaluation. By analyzing molecular dynamics simulations of Pc bound to CK2α, a series of Pc-derived peptides was rationally designed and synthesized to evaluate their binding affinity to CK2α, as well as anti-proliferative and pro-apoptotic effects against HepG2 cancer cell line. One amino acid substitutions on Pc, I192F, exhibited over 10-fold improvement in the predicted binding affinity to CK2α when compared to Pc, and a cell-permeable version, I192F-Tat, also demonstrated more potent anti-proliferative and pro-apoptotic effects against HepG2 compared to Pc. A second modification of Pc, H193W, also led to more potent cell-based activity, despite having weaker binding affinity (∼5×) to CK2α. The discovery of the I192F and H193W peptides provides new insights for further optimization of CK2 antagonist candidates as anti-cancer leads.

In spite of recent efforts to eradicate malaria in the world, this parasitic disease is still considered a major public health problem, with a total of 216 million cases of malaria and 445,000 deaths in 2016. Artemisinin-based combination therapies remain effective in most parts of the world, but recent cases of resistance in Southeast Asia have urged for novel approaches to treat malaria caused by Plasmodium falciparum. In this work, we present chloroquine analogs that exhibited high activity against sensitive and chloroquine-resistant P. falciparum blood parasites and were also active against P. berghei infected mice. Among the compounds tested, DAQ, a chloroquine analog with a more linear side chain, was shown to be the most active in vitro and in vivo, with low cytotoxicity, and therefore may serve as the basis for the development of more effective chloroquine analogs to aid malaria eradication., Graphical abstract Synthesis of N-(4-(dimethylamino)but-2-enyl)-7-chloro-quinolin-4-amine (DAQ): DAQ, a chloroquine analog with a more linear side chain, exhibited high activity against sensitive and chloroquine-resistant P. falciparum blood parasites and was also active against P. berghei infected mice.Image 1, Highlights • Resistance to recommended malaria therapies by WHO is emerging. • Chloroquine analogs with modified side chains may overcome P. falciparum resistance. • DAQ, a CQ-analog with a more linear side chain, is the most promising tested compound. • DAQ may serve as a lead compound for the design of more effective antimalarials.

We have studied the cation–p interactions of neutral aromatic ligands with the cationic amino acid residues arginine, histidine and lysine using ab initio calculations, symmetry adapted perturbation theory (SAPT), and a systematic meta-analysis of all available Protein Data Bank (PDB) X-ray structures. Quantum chemical potential energy surfaces for these interactions were obtained at the DLPNO-CCSD(T) level of theory and compared against the empirical distribution of 2012 unique protein–ligand cation–p interactions found in X-ray crystal structures. We created a workflow to extract these structures from the PDB,filtering by interaction type and residue pKa. The gas phase cation–p interaction of lysine is the strongest by more than 10 kcal mol1, but the empirical distribution of 582 X-ray structures lies away from the minimum on the interaction PES. In contrast, 1381 structures involving arginine match the underlying calculated PES with good agreement. SAPT analysis revealed that underlying differences in the balance of electrostatic and dispersion contributions are responsible for this behavior in the context of theprotein environment. The lysine–arene interaction, dominated by electrostatics, is greatly weakened by a surrounding dielectric medium and causes it to become essentially negligible in strength and without a well-defined equilibrium separation. The arginine–arene interaction involves a near equal mix ofdispersion and electrostatic attraction, which is weakened to a much smaller degree by the surrounding medium. Our results account for the paucity of cation–p interactions involving lysine, even though this is a more common residue than arginine. Aromatic ligands are most likely to interact with cationicarginine residues as this interaction is stronger than for lysine in higher polarity surroundings.

Bromodomains are protein modules that bind to acetylated lysine residues. Their interaction with histone proteins suggests that they function as "readers" of histone lysine acetylation, a component of the proposed "histone code". Bromodomain-containing proteins are often found as components of larger protein complexes with roles in fundamental cellular process including transcription. The publication of two potent ligands for the BET bromodomains in 2010 demonstrated that small molecules can inhibit the bromodomain-acetyl-lysine protein-protein interaction. These molecules display strong phenotypic effects in a number of cell lines and affect a range of cancers in vivo. This work stimulated intense interest in developing further ligands for the BET bromodomains and the design of ligands for non-BET bromodomains. Here we review the recent progress in the field with particular attention paid to ligand design, the assays employed in early ligand discovery, and the use of computational approaches to inform ligand design.

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.

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.

The fear of biological warfare agents (BWA) use by terrorists is the major concern of the security agencies and health authorities worldwide today. The non-existence of vaccines or drugs against most BWA and the possibility of genetic modified strains has turned the search for new drugs to a state of urgency. Fast in silico techniques are, therefore, perfect tools for this task once they can quickly provide structures of several new lead compounds for further experimental work. Here we try to present a mini-review on docking and molecular dynamics simulations studies applied to the drug design against the BWA Bacillus anthracis and Yersinia pestis.

CREBBP bromodomains, epigenetic “reader” proteins that recognize acetylated histone lysine residues, are a current target for cancer therapy. We show that experimental CREBBP binding affinities of small-molecules with aromatic or heteroaromatic functional groups are strongly influenced by a cation–π interaction with a positively charged arginine residue. For a series of fifteen 5-isoxazolylbenzimidazole derivatives, the strength of this non-covalent interaction is directly related to improvements in binding to CREBBP. The aromatic substituents’ inductive and resonance effects are not obviously correlated with observed structure and affinity relationships. In contrast, a coulombic electrostatic model can quantitatively predict the interaction strength. We have assessed different Molecular Mechanics (MM) and Quantum Mechanics (QM) descriptions of the protein–ligand interaction. Quantitative models for binding affinity were generated from: (1) Poisson Boltzmann Surface Area (MM-PBSA) and Generalized Born Surface Area (MM-GBSA) scoring functions that incorporated the entire ligand and (2) QM-complexation energies and (3) Electrostatic Potential Surface values (ESPs) that analyzed the varying aromatic group. A linear relationship between QM-computed ESP values is established for the cation–π interaction strength, and gives the best correlation (R2 = 0.84) with experimental binding affinities. This model also ranks ligand affinity most accurately (rs = 0.91) from the models tested. Consideration of the electrostatic potential in response to the local effects of substituents in addition to that of the aromatic ring is necessary to understand and describe the interaction with the cationic guanidinium ion. This leads to an improved understanding and the ability to quantitatively predict the magnitude of non-covalent interactions in the CREBBP active site.

The benzoxazinone and dihydroquinoxalinone fragments were employed as novel acetyl lysine mimics in the development of CREBBP bromodomain ligands. While the benzoxazinone series showed low affinity for the CREBBP bromodomain, expansion of the dihydroquinoxalinone series resulted in the first potent inhibitors of a bromodomain outside the BET family. Structural and computational studies reveal that an internal hydrogen bond stabilizes the protein‐bound conformation of the dihydroquinoxalinone series. The side chain of this series binds in an induced‐fit pocket forming a cation–π interaction with R1173 of CREBBP. The most potent compound inhibits binding of CREBBP to chromatin in U2OS cells.

JMJD2A catalyses the demethylation of di- and trimethylated lysine residues in histone tails and is a target for the development of new anticancer medicines. Mechanistic details of demethylation are yet to be elucidated and are important for the understanding of epigenetic processes. We have evaluated the initial step of histone demethylation by JMJD2A and demonstrate the dramatic effect of the protein environment upon oxygen binding using quantum mechanics/molecular mechanics (QM/MM) calculations. The changes in electronic structure have been studied for possible spin states and different conformations of O2 , using a combination of quantum and classical simulations. O2 binding to this histone demethylase is computed to occur preferentially as an end-on superoxo radical bound to a high-spin ferric centre, yielding an overall quintet ground state. The favourability of binding is strongly influenced by the surrounding protein: we have quantified this effect using an energy decomposition scheme into electrostatic and dispersion contributions. His182 and the methylated lysine assist while Glu184 and the oxoglutarate cofactor are deleterious for O2 binding. Charge separation in the superoxo-intermediate benefits from the electrostatic stabilization provided by the surrounding residues, stabilizing the binding process significantly. This work demonstrates the importance of the extended protein environment in oxygen binding, and the role of energy decomposition in understanding the physical origin of binding/recognition.

Epigenetic pathways are involved in a wide range of diseases, including cancer and neurological disorders. Specifically, histone modifying and reading processes are the most broadly studied and are targeted by several licensed drugs. Although there have been significant advances in understanding the mechanistic aspects underlying epigenetic regulation, the development of selective small-molecule inhibitors remains a challenge. Experimentally, it is generally difficult to elucidate the atomistic basis for substrate recognition, as well as the sequence of events involved in binding and the subsequent chemical processes. In this regard, computational modelling is particularly valuable, since it can provide structural features (including transition state structures along with kinetic and thermodynamic parameters) that enable both qualitative and quantitative evaluation of the mechanistic details involved. Here, we summarize knowledge gained from computational modelling studies elucidating the role of the protein environment in histone-lysine modifying and reading mechanisms. We give a perspective on the importance of calculations to aid and advance the understanding of these processes and for the future development of selective inhibitors for epigenetic regulators.

Compounds similar to lapatinib and gefitinib have been investigated as potential inhibitors of the intracellular receptor tyrosine kinase (RTK) domain of the human epidermal receptor 2 (HER2), which is a promising molecular target to the drug design of new chemotherapies for breast, lung, ovarian and colorectal cancers. In this study, we have searched potential HER2 inhibitors used for treatment of other illnesses such as hepatitis, bacterial infections and sexual impotence screened in the DrugBank. The compounds selected were subjected to virtual screening docking in order to evaluate the main interactions between them and the RTK domain of HER2. The selected compounds were investigated by flexible docking, molecular dynamics studies and ΔG bind calculations. The results suggest that antrafenine, saprisartan, reserpine, irinotecan and udenafil are potential candidates to inhibit the RTK domain of HER2.

The literature has reported that ferriprotoporphyrin IX (hematin) intoxicates the malarial parasite through competition with NADH for the active site of the enzyme lactate dehydrogenase (LDH). In order to avoid this, the parasite polymerizes hematin to hemozoin. The quinoline derivatives are believed to form complexes with dimeric hematin, avoiding the formation of hemozoin and still inhibiting LDH. In order to investigate this hypothesis we calculated the docking energies of NADH and some quinoline derivatives (in the free forms and in complex with dimeric hematin) in the active site of the Plasmodium falciparum LDH (PfLDH). Ours results showed better docking score values to the complexes when compared to the free compounds, pointing them as more efficient inhibitors of Pf_LDH. Further we performed Molecular Dynamics (MD) simulations studies on the best docking conformation of the complex chloroquine-dimeric hematin with PfLDH. Our in silico results corroborate experimental data suggesting a possible action route for the quinoline derivatives in the inhibition of PfLDH.

In spite of recent progress, cancer is still one of the most serious health problems of mankind. Recently, it has been discovered that tumor hypoxia can be exploited for selective anticancer treatment using radiosensitizers that are activated only under hypoxic conditions. The most commonly used radiosensitizers are the 5-nitroimidazole derivatives. The toxicity of bioreductive anticancer drugs, such as radiosensitizers is associated to their interaction with DNA. In this work, we have investigated the interaction between the model radiosensitizers metronizole, nimorazole and secnidazole with salmon DNA in order to get insights on the drug–macromolecule interactions. To this end, we have employed NMR techniques (PFG NMR spectra and spin–lattice relaxation rates) in combination with theoretical tools, such as docking calculations and MD simulations. Initially, results show that the δ values are not the most appropriated NMR parameters to map the interaction topology of drug–macromolecule complexes. Furthermore our data indicate that radiosensitizers, in the inactive form, interact considerably with DNA, significantly increasing its toxicity. In fact, we obtained a good agreement between that technique and docking and MD simulations. This suggests that improvements in the structures of these molecules in order to achieve new and more selective bioreductive anticancer drugs are still necessary.

DNA topoisomerase I from Plasmodium falciparum (PfTopoI), a potential selective target for chemotherapy and drug development against malaria, is used here, together with human Topo I (HssTopoI), for docking, molecular dynamics (MD) studies and experimental assays. Six synthetic isoflavonoid derivatives and the known PfTopoI inhibitors camptothecin and topotecan were evaluated in parallel. Theoretical results suggest that these compounds dock in the binding site of camptothecin and topotecan inside both enzymes and that LQB223 binds selectively in PfTopoI. In vitro tests against P. falciparum blood parasites corroborated the theoretical findings. The selectivity index (SI) of LQB223 ≥ 98 suggests that this molecule is the most promising in the group of compounds tested. In vivo experiments in mice infected with P. berghei showed that LQB223 has an antimalarial activity similar to that of chloroquine.

The suitable computation of accurate atomic charges for the GROMACS topology *.itp files of small molecules, generated in the PRODRG server, has been a tricky task nowadays because it does not calculate atomic charges using an ab initio method. Usually additional steps of structure optimization and charges calculation, followed by a tedious manual replacement of atomic charges in the *.itp file, are needed. In order to assist this task, we report here the ITP Adjuster 1.0, a utility program developed to perform the replacement of the PRODRG charges in the *.itp files of small molecules by ab initio charges.

Previously, we have reported that nicotinic acid adenine dinucleotide phosphate (NAADP) and cyclic ADP-ribose (cADPR) are calcium mobilising messengers that have significant roles in excitation-contraction coupling and disease development of the heart. The enzyme responsible for producing these molecules is a potential therapeutic target, but its identity remains unclear. CD38 is an ADP-ribosyl cyclase (ARC) that catalyses synthesis of NAADP and cADPR in various tissues, with one caveat: the targets of NAADP and cADPR are intracellular while the active site of CD38 is generally thought to be located on the extracellular surface, making the physiological relevance of this ecto-enzyme questionable. Permeabilised cardiac myocytes showed greater ARC activity than intact cells, suggesting the presence of intracellular ARCs. Immunostaining experiments showed that cell permeabilisation was critical for labelling CD38 with a striated staining pattern. We further showed that sarcoplasmic reticulum enriched membrane preparations from sheep hearts exhibited ARC activity with CD38-like properties. Genetic knock-out of CD38 in mice suppressed NAADP and cADPR synthesis. Isolated ventricular myocytes from CD38-/- mice showed reduced peak calcium transients and contractions, both under basal conditions and when stimulated by the β-adrenoceptor agonist, isoprenaline. CD38-/- hearts were more resistant to arrhythmias during over-stimulation of β-adrenoceptor pathway. A novel ARC inhibitor, SAN4825, suppressed the in vitro synthesis of NAADP and cADPR, and in blind docking simulations was shown to target the active site of CD38. The observations are consistent with an important role for CD38 in producing cADPR and NAADP intracellularly under basal conditions; the production of both these substances is increased by β-adrenoceptor stimulation and the elevated levels enhance arrhythmogenicity. In conclusion, we provide evidence that CD38 generates NAADP and cADPR in situ intracellularly and this enzymatic activity plays a role in both physiology and pathophysiology of the heart.

The Plasmodium falciparum lactate dehydrogenase enzyme (PfLDH) has been considered as a potential molecular target for antimalarials due to this parasite's dependence on glycolysis for energy production. Because the LDH enzymes found in P. vivax, P. malariae and P. ovale (pLDH) all exhibit ∼90% identity to PfLDH, it would be desirable to have new anti-pLDH drugs, particularly ones that are effective against P. falciparum, the most virulent species of human malaria. Our present work used docking studies to select potential inhibitors of pLDH, which were then tested for antimalarial activity against P. falciparum in vitro and P. berghei malaria in mice. A virtual screening in DrugBank for analogs of NADH (an essential cofactor to pLDH) and computational studies were undertaken, and the potential binding of the selected compounds to the PfLDH active site was analyzed using Molegro Virtual Docker software. Fifty compounds were selected based on their similarity to NADH. The compounds with the best binding energies (itraconazole, atorvastatin and posaconazole) were tested against P. falciparum chloroquine-resistant blood parasites. All three compounds proved to be active in two immunoenzymatic assays performed in parallel using monoclonals specific to PfLDH or a histidine rich protein (HRP2). The IC(50) values for each drug in both tests were similar, were lowest for posaconazole (

Chloroquine (CQ) is a cost effective antimalarial drug with a relatively good safety profile (or therapeutic index). However, CQ is no longer used alone to treat patients with Plasmodium falciparum due to the emergence and spread of CQ-resistant strains, also reported for P. vivax. Despite CQ resistance, novel drug candidates based on the structure of CQ continue to be considered, as in the present work. One CQ analog was synthesized as monoquinoline (MAQ) and compared with a previously synthesized bisquinoline (BAQ), both tested against P. falciparum in vitro and against P. berghei in mice, then evaluated in vitro for their cytotoxicity and ability to inhibit hemozoin formation. Their interactions with residues present in the NADH binding site of P falciparum lactate dehydrogenase were evaluated using docking analysis software. Both compounds were active in the nanomolar range evaluated through the HRPII and hypoxanthine tests. MAQ and BAQ derivatives were not toxic, and both compounds significantly inhibited hemozoin formation, in a dose-dependent manner. MAQ had a higher selectivity index than BAQ and both compounds were weak PfLDH inhibitors, a result previously reported also for CQ. Taken together, the two CQ analogues represent promising molecules which seem to act in a crucial point for the parasite, inhibiting hemozoin formation.