In silico druggability assessment of the NUDIX hydrolase protein family as a workflow for target prioritization

Computational chemistry has now been widely accepted as a useful tool for shortening lead times in early drug discovery. When selecting new potential drug targets, it is important to assess the likelihood of finding suitable starting points for lead generation before pursuing costly high-throughput screening campaigns. By exploiting available high-resolution crystal structures, an in silico druggability assessment can facilitate the decision of whether, and in cases where several protein family members exist, which of these to pursue experimentally. Many of the algorithms and software suites commonly applied for in silico druggability assessment are complex, technically challenging and not always user-friendly. Here we applied the intuitive open access servers of DoGSite, FTMap and CryptoSite to comprehensively predict ligand binding pockets, druggability scores and conformationally active regions of the NUDIX protein family. In parallel we analyzed potential ligand binding sites, their druggability and pocket parameter using Schrödinger's SiteMap. Then an in silico docking cascade of a subset of the ZINC FragNow library using the Glide docking program was performed to assess identified pockets for large-scale small-molecule binding. Subsequently, this initial dual ranking of druggable sites within the NUDIX protein family was benchmarked against experimental hit rates obtained both in-house and by others from traditional biochemical and fragment screening campaigns. The observed correlation suggests that the presented user-friendly workflow of a dual parallel in silico druggability assessment is applicable as a standalone method for decision on target prioritization and exclusion in future screening campaigns

AXL and CAV-1 play a role for MTH1 inhibitor TH1579 sensitivity in cutaneous malignant melanoma

Cutaneous malignant melanoma (CMM) is the deadliest form of skin cancer and clinically challenging due to its propensity to develop therapy resistance. Reactive oxygen species (ROS) can induce DNA damage and play a significant role in CMM. MTH1 protein protects from ROS damage and is often overexpressed in different cancer types including CMM. Herein, we report that MTH1 inhibitor TH1579 induced ROS levels, increased DNA damage responses, caused mitotic arrest and suppressed CMM proliferation leading to cell death both in vitro and in an in vivo xenograft CMM zebrafish disease model. TH1579 was more potent in abrogating cell proliferation and inducing cell death in a heterogeneous co-culture setting when compared with CMM standard treatments, vemurafenib or trametinib, showing its broad anticancer activity. Silencing MTH1 alone exhibited similar cytotoxic effects with concomitant induction of mitotic arrest and ROS induction culminating in cell death in most CMM cell lines tested, further emphasizing the importance of MTH1 in CMM cells. Furthermore, overexpression of receptor tyrosine kinase AXL, previously demonstrated to contribute to BRAF inhibitor resistance, sensitized BRAF mutant and BRAF/NRAS wildtype CMM cells to TH1579. AXL overexpression culminated in increased ROS levels in CMM cells. Moreover, silencing of a protein that has shown opposing effects on cell proliferation, CAV-1, decreased sensitivity to TH1579 in a BRAF inhibitor resistant cell line. AXL-MTH1 and CAV-1-MTH1 mRNA expressions were correlated as seen in CMM clinical samples. Finally, TH1579 in combination with BRAF inhibitor exhibited a more potent cell killing effect in BRAF mutant cells both in vitro and in vivo. In summary, we show that TH1579-mediated efficacy is independent of BRAF/NRAS mutational status but dependent on the expression of AXL and CAV-1

Recent advances in electronic structure theory and their influence on the accuracy of ab initio potential energy surfaces

Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F+H2 yields HF+H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces

MutT homologue 1 (MTH1) catalyzes the hydrolysis of mutagenic O6-methyl-dGTP.

Nucleotides in the free pool are more susceptible to nonenzymatic methylation than those protected in the DNA double helix. Methylated nucleotides like O6-methyl-dGTP can be mutagenic and toxic if incorporated into DNA. Removal of methylated nucleotides from the nucleotide pool may therefore be important to maintain genome integrity. We show that MutT homologue 1 (MTH1) efficiently catalyzes the hydrolysis of O6-methyl-dGTP with a catalytic efficiency similar to that for 8-oxo-dGTP. O6-methyl-dGTP activity is exclusive to MTH1 among human NUDIX proteins and conserved through evolution but not found in bacterial MutT. We present a high resolution crystal structure of human and zebrafish MTH1 in complex with O6-methyl-dGMP. By microinjecting fertilized zebrafish eggs with O6-methyl-dGTP and inhibiting MTH1 we demonstrate that survival is dependent on active MTH1 in vivo. O6-methyl-dG levels are higher in DNA extracted from zebrafish embryos microinjected with O6-methyl-dGTP and inhibition of O6-methylguanine-DNA methyl transferase (MGMT) increases the toxicity of O6-methyl-dGTP demonstrating that O6-methyl-dGTP is incorporated into DNA. MTH1 deficiency sensitizes human cells to the alkylating agent Temozolomide, a sensitization that is more pronounced upon MGMT inhibition. These results expand the cellular MTH1 function and suggests MTH1 also is important for removal of methylated nucleotides from the nucleotide pool

A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry

We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical derivation of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications. Moreover, much material that is found scattered in the literature has been put together here to facilitate comprehension and to serve as a handy reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and subeqn package

Transancestral mapping and genetic load in systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is an autoimmune disease with marked gender and ethnic disparities. We report a large transancestral association study of SLE using Immunochip genotype data from 27,574 individuals of European (EA), African (AA) and Hispanic Amerindian (HA) ancestry. We identify 58 distinct non-HLA regions in EA, 9 in AA and 16 in HA (∼50% of these regions have multiple independent associations); these include 24 novel SLE regions (P<5 × 10-8), refined association signals in established regions, extended associations to additional ancestries, and a disentangled complex HLA multigenic effect. The risk allele count (genetic load) exhibits an accelerating pattern of SLE risk, leading us to posit a cumulative hit hypothesis for autoimmune disease. Comparing results across the three ancestries identifies both ancestry-dependent and ancestry-independent contributions to SLE risk. Our results are consistent with the unique and complex histories of the populations sampled, and collectively help clarify the genetic architecture and ethnic disparities in SLE

Calculation of molecular thermochemical data and their availability in databases

Thermodynamic properties of molecules can be obtained by experiment, by statistical mechanics in conjunction with electronic structure theory and by empirical rules like group additivity. The latter two methods are briefly re-viewed in this chapter. The overview of electronic structure methods is intended for readers less experienced in electronic structure theory and focuses on concepts without going into mathematical details. This is followed by a brief description of group additivity schemes; finally, an overview of databases listing reliable thermochemical data is given

Transcriptome and genome sequencing uncovers functional variation in humans

Summary Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of mRNA and miRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project – the first uniformly processed RNA-seq data from multiple human populations with high-quality genome sequences. We discovered extremely widespread genetic variation affecting regulation of the majority of genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on cellular mechanisms of regulatory and loss-of-function variation, and allowed us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome