Search

The inhibition of the hERG potassium channel is closely related to the prolonged QT interval, and thus assessing this risk could greatly facilitate the development of therapeutic compounds and the withdrawal of hazardous marketed drugs. The recent increase in SAR information about hERG inhibitors in public databases has led to many successful applications of machine learning techniques to predict hERG inhibition. However, most of these reports constructed their prediction models based on only one SAR database because the differences in the data format and ontology hindered the integration of the databases. In this study, we curated the hERG-related data in ChEMBL, PubChem, GOSTAR, and hERGCentral, and integrated them into the largest database about hERG inhibition by small molecules. Assessment of structural diversity using Murcko frameworks revealed that the integrated database contains more than twice as many chemical scaffolds for hERG inhibitors than any of the individual databases, and covers 18.2% of the Murcko framework-based chemical space occupied by the compounds in ChEMBL. The database provides the most comprehensive information about hERG inhibitors and will be useful to design safer compounds for drug discovery. The database is freely available at . [ABSTRACT FROM AUTHOR]

We previously reported the structure-based design of a highly potent hematopoietic cell kinase (HCK) inhibitor, a pyrrolo-pyrimidine compound designated RK-20449, for treatment of recurrent leukemia. Herein we report the synthesis and structure–activity relationships of some amino acid derivatives of 7-substituted pyrrolo-pyrimidine. Although these derivatives had the same predicted binding conformation as RK-20449, their IC 50 values were 100–1000 times larger than that of the parent compound. We assumed that the basicity of the amine nitrogen, which formed an ionic bond with Asp348 of HCK, markedly affected inhibitory activity against HCK. The pKa values of the nitrogen were predicted by means of an ab initio quantum mechanical method, and complexes of the derivatives with HCK were analyzed by X-ray crystallography. We observed a significant correlation between the predicted pKa and IC 50 values, and the crystal structures of the less potent derivatives showed various types of defects around the ionic bond. [ABSTRACT FROM AUTHOR]

Abstract: In drug discovery process, improvement of ADME/Tox properties of lead compounds including metabolic stability is critically important. Cytochrome P450 (CYP) is one of the major metabolizing enzymes and the prediction of sites of metabolism (SOM) on the given lead compounds is key information to modify the compounds to be more stable against metabolism. There are two factors essentially important in SOM prediction. First is accessibility of each substrate atom to the oxygenated Fe atom of heme in a CYP protein, and the other is the oxidative reactivity of each substrate atom. To predict accessibility of substrate atoms to the heme iron, conventional protein-rigid docking simulations have been applied. However, the docking simulations without consideration of protein flexibility often lead to incorrect answers in the case of very flexible proteins such as CYP3A4. In this study, we demonstrated an approach utilizing molecular dynamics (MD) simulation for SOM prediction in which multiple MD runs were executed using different initial structures. We applied this strategy to CYP3A4 and carbamazepine (CBZ) complex. Through 10ns MD simulations started from five different CYP3A4-CBZ complex models, our approach correctly predicted SOM observed in experiments. The experimentally known epoxidized sites of CBZ by CYP3A4 were successfully predicted as the most accessible sites to the heme iron that was judged from a numerical analysis of calculated ΔG binding and the frequency of appearance. In contrast, the predictions using protein-rigid docking methods hardly provided the correct SOM due to protein flexibility or inaccuracy of the scoring functions. Our strategy using MD simulation with multiple initial structures will be one of the reliable methods for SOM prediction. [Copyright &y& Elsevier]

Abstract: To ensure a continuing pipeline in pharmaceutical research, lead candidates must possess appropriate metabolic stability in the drug discovery process. In vitro ADMET (absorption, distribution, metabolism, elimination, and toxicity) screening provides us with useful information regarding the metabolic stability of compounds. However, before the synthesis stage, an efficient process is required in order to deal with the vast quantity of data from large compound libraries and high-throughput screening. Here we have derived a relationship between the chemical structure and its metabolic stability for a data set of in-house compounds by means of various in silico machine learning such as random forest, support vector machine (SVM), logistic regression, and recursive partitioning. For model building, 1952 proprietary compounds comprising two classes (stable/unstable) were used with 193 descriptors calculated by Molecular Operating Environment. The results using test compounds have demonstrated that all classifiers yielded satisfactory results (accuracy>0.8, sensitivity>0.9, specificity>0.6, and precision>0.8). Above all, classification by random forest as well as SVM yielded kappa values of approximately 0.7 in an independent validation set, slightly higher than other classification tools. These results suggest that nonlinear/ensemble-based classification methods might prove useful in the area of in silico ADME modeling. [Copyright &y& Elsevier]

No explicit π-π interaction term has been incorporated in the conventional molecular dynamics (MD) simulation programs in spite of its significant role in the folding of biomolecules and the clustering of organic chemicals. In this article, we propose a technique to emphasize the effect of π-π interactions using a function of energy and implement it into an MD simulation program. Several trial calculations show that the π-π incorporated program gives improved results consistent with experimental data on atom geometry and has no unfavorable interference with the conventional computational framework. This indicates an importance of the explicit consideration of π-π interactions in MD simulation. © 2007 Wiley Periodicals, Inc.J Comput Chem 2007 [ABSTRACT FROM AUTHOR]

In neurological and neuropsychiatric diseases, different brain regions are affected, and differences in gene expression patterns could potentially explain this mechanism. However, limited studies have precisely explored gene expression in different regions of the human brain. In this study, we performed long-read RNA sequencing on three different brain regions of the same individuals: the cerebellum, hypothalamus, and temporal cortex. Despite stringent filtering criteria excluding isoforms predicted to be artifacts, over half of the isoforms expressed in multiple samples across multiple regions were found to be unregistered in the GENCODE reference. We then especially focused on genes with different major isoforms in each brain region, even with similar overall expression levels, and identified that many of such genes including GAS7 might have distinct roles in dendritic spine and neuronal formation in each region. We also found that DNA methylation might, in part, drive different isoform expressions in different regions. These findings highlight the significance of analyzing isoforms expressed in disease-relevant sites. [ABSTRACT FROM AUTHOR]

Cytochrome P450 (CYP) is deeply involved in the metabolism of chemicals including pharmaceuticals. Therefore, polymorphisms of this enzyme have been widely studied to avoid unfavorable side effects of drugs in chemotherapy. In this work, we performed computational analysis of the mechanism of the decrease in enzymatic activity for three typical polymorphisms in CYP 2C9 species: *2, *3, and *5. Based on the equilibrated structure obtained by molecular dynamics simulation, the volume of the binding pocket and the fluctuation of amino residues responsible for substrate holding were compared between the wild type and the three variants. Further docking simulation was carried out to evaluate the appropriateness of the binding pocket to accommodate substrate chemicals. Every polymorphic variant was suggested to be inferior to the wild type in enzymatic ability from the structural viewpoint. F-G helices were obviously displaced outward in CYP2C9*2. Expansion of the binding pocket, especially the space near F′ helix, was remarkable in CYP2C9*3. Disappearance of the hydrogen bond between K helix and β4 loop was observed in CYP2C9*5. The reduction of catalytic activity of those variants can be explained from the deformation of the binding pocket and the consequent change in binding mode of substrate chemicals. The computational approach is effective for predicting the enzymatic activity of polymorphic variants of CYP. This prediction will be helpful for advanced drug design because calculations forecast unexpected change in drug efficacy for individuals. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010 [ABSTRACT FROM AUTHOR]

Chagas disease is caused by infection with the protozoan parasite Trypanosoma cruzi (T. cruzi). It was originally a Latin American endemic health problem, but now is expanding worldwide as a result of increasing migration. The currently available drugs for Chagas disease, benznidazole and nifurtimox, provoke severe adverse effects, and thus the development of new drugs is urgently required. Ubiquinone (UQ) is essential for respiratory chain and redox balance in trypanosomatid protozoans, therefore we aimed to provide evidence that inhibitors of the UQ biosynthesis have trypanocidal activities. In this study, inhibitors of the human COQ7, a key enzyme of the UQ synthesis, were tested for their trypanocidal activities because they were expected to cross-react and inhibit trypanosomal COQ7 due to their genetic homology. We show the trypanocidal activity of a newly found human COQ7 inhibitor, an oxazinoquinoline derivative. The structurally similar compounds were selected from the commercially available compounds by 2D and 3D ligand-based similarity searches. Among 38 compounds selected, 12 compounds with the oxazinoquinoline structure inhibited significantly the growth of epimastigotes of T. cruzi. The most effective 3 compounds also showed the significant antitrypanosomal activity against the mammalian stage of T. cruzi at lower concentrations than benznidazole, a commonly used drug today. We found that epimastigotes treated with the inhibitor contained reduced levels of UQ9. Further, the growth of epimastigotes treated with the inhibitors was partially rescued by UQ10 supplementation to the culture medium. These results suggest that the antitrypanosomal mechanism of the oxazinoquinoline derivatives results from inhibition of the trypanosomal UQ synthesis leading to a shortage of the UQ pool. Our data indicate that the UQ synthesis pathway of T. cruzi is a promising drug target for Chagas disease. [ABSTRACT FROM AUTHOR]

Over the past few decades, obesity has become a public health issue of global concern. Even though disparities exist between human populations, e.g., the higher liver fat content of the Japanese despite a lower body mass index (BMI), studies on the genetics of obesity still largely focus on populations of European descent, leading to a dearth of genetic data on non-European populations. In this context, this study aimed to establish a broad picture of the genetic attributes of the Japanese population, by examining a representative sample of 18,889 individuals participating in the Tohoku Medical Megabank Project cohort. We applied linear mixed model methods to 17 traits related to obesity and associated diseases to estimate the heritabilities explained by common genetic variants and the genetic correlations between each pair of traits. These analyses allowed us to quantify the SNP heritability of health indicators such as BMI (0.248 6 0.032) and HDL cholesterol (0.324 6 0.031), and to provide one of the few estimates of the SNP heritability of cystatin C in unrelated individuals (0.260 6 0.025). We discuss potential differences between the Japanese and people of European ancestry with respect to the genetic correlations between urinary biomarkers and adiposity traits, for which large estimates were obtained. For instance, the genetic correlations between urine potassium level and the values for weight, BMI, waist circumference, and waist-to-height ratio ranged from 0.290 to 0.559, much higher than the corresponding estimates in the UK Biobank. [ABSTRACT FROM AUTHOR]

Assessing the hERG liability in the early stages of drug discovery programs is important. The recent increase of hERG-related information in public databases enabled various successful applications of machine learning techniques to predict hERG inhibition. However, most of these researches constructed the datasets from only one database, limiting the predictability and scope of the models. In this study, a hERG classification model was constructed using the largest dataset for hERG inhibition built by integrating multiple databases. The integrated dataset consisted of more than 291,000 structurally diverse compounds derived from ChEMBL, GOSTAR, PubChem, and hERGCentral. The prediction model was built by support vector machine (SVM) with descriptor selection based on Non-dominated Sorting Genetic Algorithm-II (NSGA-II) to optimize the descriptor set for maximum prediction performance with the minimal number of descriptors. The SVM classification model using 72 selected descriptors and ECFP_4 structural fingerprints recorded kappa statistics of 0.733 and accuracy of 0.984 for the test set, substantially outperforming the prediction performance of the current commercial applications for hERG prediction. Finally, the applicability domain of the prediction model was assessed based on the molecular similarity between the training set and test set compounds. [ABSTRACT FROM AUTHOR]

A series of novel pyrrolo[2,3- d ]pyrimidines were synthesized by introducing 15 different amino acids to 7-cyclohexyl-5-(4-phenoxyphenyl)-7 H -pyrrolo[2,3- d ]pyrimidine-4-amine. Compounds with potent activities against HCK and FLT3-ITD were evaluated in viability studies with acute myeloid leukemia cell line MV4-11. Our structure activity relationship analyses lead to the identification of compound 31 , which exhibited potent HCK and FLT3-ITD inhibition and activity against the MV4-11 cell line. [ABSTRACT FROM AUTHOR]

[Display omitted] The identification, structure–activity relationships (SARs), and biological effects of new antimalarials consisting of a 2,3,4,9-tetrahydro-1 H -β-carboline core, a coumarin ring, and an oxyalkanoyl linker are described. A cell-based phenotypic approach was employed in this search for novel antimalarial drugs with unique modes of action. Our screening campaign of the RIKEN compound library succeeded in the identification of the known tetrahydro-β-carboline derivative (4e) as a hit compound showing significant in vitro activity. SAR studies on this chemical series led to the discovery of compound 4h having a (R)-methyl group on the oxyacetyl linker with potent inhibition of parasite growth (IC 50 = 2.0 nM). Compound 4h was also found to exhibit significant in vivo antimalarial effects in mouse models. Furthermore, molecular modeling studies on 4e , 4h , and its diastereomer (4j) suggested that the (R)-methyl group of 4h forces the preferential adoption of a specific conformer which is considered to be an active conformer. [ABSTRACT FROM AUTHOR]

Abstract: EGFR is a target protein for the treatment of non small cell lung cancer (NSCLC). The mutations associated with the activation of EGFR kinase activity, such as L858R and G719S, destabilize the inactive conformation of EGFR and are closely linked with the development of NSCLC. The additional T790M mutation reportedly causes drug resistance against the commercially available EGFR inhibitors, gefitinib and erlotinib. In this study, we searched for novel G719S/T790M EGFR inhibitors by a new in silico screening strategy, using two datasets. The results of in silico screening using protein-ligand docking are affected by the selection of 3D structure of the target protein. As the first strategy, we chose the 3D structures for in silico screening by test dockings using the G719S/T790M crystal structure, its molecular dynamics snapshots, and known inhibitors of the drug-resistant EGFR. In the second strategy, we selected the 3D structures by test dockings using all of the EGFR structures, regardless of the mutations, and all of the known EGFR inhibitors. Using each of the 3D structures selected by the strategies, 1000 compounds were chosen from the 71,588 compounds. Kinase assays identified 15 G719S/T790M EGFR inhibitors, including two compounds with novel scaffolds. Analyses of their structure–activity relationships revealed that interactions with the mutated Met790 residue specifically increase the inhibitory activity against G719S/T790M EGFR. [Copyright &y& Elsevier]

Abstract: A new screening method using fluorescent correlation spectroscopy was developed to select kinase inhibitors that competitively inhibit the binding of a fluorescently labeled substrate peptide. Using the method, among approximately 700 candidate compounds selected by virtual screening, we identified a novel Pim-1 kinase inhibitor targeting its peptide binding residues. X-ray crystal analysis of the complex structure of Pim-1 with the inhibitor indicated that the inhibitor actually binds to the ATP-binding site and also forms direct interactions with residues (Asp128 and Glu171) that bind the substrate peptide. These interactions, which cause small side-chain movements, seem to affect the binding ability of the fluorescently labeled substrate. The compound inhibited Pim-1 kinase in vitro, with an IC50 value of 150 nM. Treatment of cultured leukemia cells with the compound reduced the amount of p21 and increased the amount of p27, due to Pim-1 inhibition, and then triggered apoptosis after cell-cycle arrest at the G1/S phase. This screening method may be widely applicable for the identification of various new Pim-1 kinase inhibitors targeting the residues that bind the substrate peptide. [Copyright &y& Elsevier]

Abstract: Human CYP3A4 catalyzes the 10,11-epoxidation of carbamazepine (CBZ). However, the epoxide is less stable in terms of potential energy than hydroxides of the six-membered aromatic ring. To clarify the reason why CYP3A4 produces such an energetically unfavorable compound, the mechanism of epoxidation of CBZ by CYP3A4 was investigated by theoretical calculations. The reaction consisted of two elementary processes in which two C–O bonds were generated stepwise. The rate-determining step was the first one and the activation energy was 21.3kcal/mol at the DFT (B3LYP/6-31G∗∗) level. The activation energy level of the first step of the 10,11-epoxidation was lower than that of the hydroxylation of the aromatic ring. For this reason, 10,11-epoxidation is more probable than hydroxylation of the aromatic ring, and only 10,11-epoxide is formed. [Copyright &y& Elsevier]