Journal of Computer-Aided Molecular Design最新文献

{"title":"Molecular dynamics simulations as a guide for modulating small molecule aggregation","authors":"Azam Nesabi,&nbsp;Jas Kalayan,&nbsp;Sara Al-Rawashdeh,&nbsp;Mohammad A. Ghattas,&nbsp;Richard A. Bryce","doi":"10.1007/s10822-024-00557-1","DOIUrl":"10.1007/s10822-024-00557-1","url":null,"abstract":"<div><p>Small colloidally aggregating molecules (SCAMs) can be problematic for biological assays in drug discovery campaigns. However, the self-associating properties of SCAMs have potential applications in drug delivery and analytical biochemistry. Consequently, the ability to predict the aggregation propensity of a small organic molecule is of considerable interest. Chemoinformatics-based filters such as ChemAGG and Aggregator Advisor offer rapid assessment but are limited by the assay quality and structural diversity of their training set data. Complementary to these tools, we explore here the ability of molecular dynamics (MD) simulations as a physics-based method capable of predicting the aggregation propensity of diverse chemical structures. For a set of 32 molecules, using simulations of 100 ns in explicit solvent, we find a success rate of 97% (one molecule misclassified) as opposed to 75% by Aggregator Advisor and 72% by ChemAGG. These short timescale MD simulations are representative of longer microsecond trajectories and yield an informative spectrum of aggregation propensities across the set of solutes, capturing the dynamic behaviour of weakly aggregating compounds. Implicit solvent simulations using the generalized Born model were less successful in predicting aggregation propensity. MD simulations were also performed to explore structure-aggregation relationships for selected molecules, identifying chemical modifications that reversed the predicted behaviour of a given aggregator/non-aggregator compound. While lower throughput than rapid cheminformatics-based SCAM filters, MD-based prediction of aggregation has potential to be deployed on the scale of focused subsets of moderate size, and, depending on the target application, provide guidance on removing or optimizing a compound’s aggregation propensity.</p><h3>Graphical Abstract</h3><div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-03-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10933209/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140108733","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Molecule auto-correction to facilitate molecular design","authors":"Alan Kerstjens,&nbsp;Hans De Winter","doi":"10.1007/s10822-024-00549-1","DOIUrl":"10.1007/s10822-024-00549-1","url":null,"abstract":"<div><p>Ensuring that computationally designed molecules are chemically reasonable is at best cumbersome. We present a molecule correction algorithm that morphs invalid molecular graphs into structurally related valid analogs. The algorithm is implemented as a tree search, guided by a set of policies to minimize its cost. We showcase how the algorithm can be applied to molecular design, either as a post-processing step or as an integral part of molecule generators.</p><h3>Graphical abstract</h3>\u0000<div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC10873457/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139740104","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Rethinking the applicability domain analysis in QSAR models","authors":"Jose R. Mora,&nbsp;Edgar A. Marquez,&nbsp;Noel Pérez-Pérez,&nbsp;Ernesto Contreras-Torres,&nbsp;Yunierkis Perez-Castillo,&nbsp;Guillermin Agüero-Chapin,&nbsp;Felix Martinez-Rios,&nbsp;Yovani Marrero-Ponce,&nbsp;Stephen J. Barigye","doi":"10.1007/s10822-024-00550-8","DOIUrl":"10.1007/s10822-024-00550-8","url":null,"abstract":"<div><p>Notwithstanding the wide adoption of the OECD principles (or best practices) for QSAR modeling, disparities between <i>in silico</i> predictions and experimental results are frequent, suggesting that model predictions are often too optimistic. Of these OECD principles, the applicability domain (AD) estimation has been recognized in several reports in the literature to be one of the most challenging, implying that the actual reliability measures of model predictions are often unreliable. Applying tree-based error analysis workflows on 5 QSAR models reported in the literature and available in the QsarDB repository, i.e., androgen receptor bioactivity (agonists, antagonists, and binders, respectively) and membrane permeability (highest membrane permeability and the intrinsic permeability), we demonstrate that predictions erroneously tagged as reliable (AD prediction errors) overwhelmingly correspond to instances in subspaces (cohorts) with the highest prediction error rates, highlighting the inhomogeneity of the AD space. In this sense, we call for more stringent AD analysis guidelines which require the incorporation of model error analysis schemes, to provide critical insight on the reliability of underlying AD algorithms. Additionally, any selected AD method should be rigorously validated to demonstrate its suitability for the model space over which it is applied. These steps will ultimately contribute to more accurate estimations of the reliability of model predictions. Finally, error analysis may also be useful in “rational” model refinement in that data expansion efforts and model retraining are focused on cohorts with the highest error rates.</p></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139728695","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey","authors":"Md Fulbabu Sk,&nbsp;Sunanda Samanta,&nbsp;Sayan Poddar,&nbsp;Parimal Kar","doi":"10.1007/s10822-023-00548-8","DOIUrl":"10.1007/s10822-023-00548-8","url":null,"abstract":"<div><p>The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop’s flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG<i>-in</i> conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-<i>out</i> conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-<i>out</i> inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI’s binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.</p></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139696630","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A deep neural network: mechanistic hybrid model to predict pharmacokinetics in rat","authors":"Florian Führer,&nbsp;Andrea Gruber,&nbsp;Holger Diedam,&nbsp;Andreas H. Göller,&nbsp;Stephan Menz,&nbsp;Sebastian Schneckener","doi":"10.1007/s10822-023-00547-9","DOIUrl":"10.1007/s10822-023-00547-9","url":null,"abstract":"<div><p>An important aspect in the development of small molecules as drugs or agrochemicals is their systemic availability after intravenous and oral administration. The prediction of the systemic availability from the chemical structure of a potential candidate is highly desirable, as it allows to focus the drug or agrochemical development on compounds with a favorable kinetic profile. However, such predictions are challenging as the availability is the result of the complex interplay between molecular properties, biology and physiology and training data is rare. In this work we improve the hybrid model developed earlier (Schneckener in J Chem Inf Model 59:4893–4905, 2019). We reduce the median fold change error for the total oral exposure from 2.85 to 2.35 and for intravenous administration from 1.95 to 1.62. This is achieved by training on a larger data set, improving the neural network architecture as well as the parametrization of mechanistic model. Further, we extend our approach to predict additional endpoints and to handle different covariates, like sex and dosage form. In contrast to a pure machine learning model, our model is able to predict new end points on which it has not been trained. We demonstrate this feature by predicting the exposure over the first 24 h, while the model has only been trained on the total exposure.</p></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139641396","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Identification of a druggable site on GRP78 at the GRP78-SARS-CoV-2 interface and virtual screening of compounds to disrupt that interface","authors":"Maria Lazou,&nbsp;Jonathan R. Hutton,&nbsp;Arijit Chakravarty,&nbsp;Diane Joseph-McCarthy","doi":"10.1007/s10822-023-00546-w","DOIUrl":"10.1007/s10822-023-00546-w","url":null,"abstract":"<div><p>SARS-CoV-2, the virus that causes COVID-19, led to a global health emergency that claimed the lives of millions. Despite the widespread availability of vaccines, the virus continues to exist in the population in an endemic state which allows for the continued emergence of new variants. Most of the current vaccines target the spike glycoprotein interface of SARS-CoV-2, creating a selection pressure favoring viral immune evasion. Antivirals targeting other molecular interactions of SARS-CoV-2 can help slow viral evolution by providing orthogonal selection pressures on the virus. GRP78 is a host auxiliary factor that mediates binding of the SARS-CoV-2 spike protein to human cellular ACE2, the primary pathway of cell infection. As GRP78 forms a ternary complex with SARS-CoV-2 spike protein and ACE2, disrupting the formation of this complex is expected to hinder viral entry into host cells. Here, we developed a model of the GRP78-Spike RBD-ACE2 complex. We then used that model together with hot spot mapping of the GRP78 structure to identify the putative binding site for spike protein on GRP78. Next, we performed structure-based virtual screening of known drug/candidate drug libraries to identify binders to GRP78 that are expected to disrupt spike protein binding to the GRP78, and thereby preventing viral entry to the host cell. A subset of these compounds has previously been shown to have some activity against SARS-CoV-2. The identified hits are starting points for the further development of novel SARS-CoV-2 therapeutics, potentially serving as proof-of-concept for GRP78 as a potential drug target for other viruses.</p></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.0,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139540966","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Molecular dynamics study on micelle-small molecule interactions: developing a strategy for an extensive comparison","authors":"Aleksei Kabedev,&nbsp;Christel A. S. Bergström,&nbsp;Per Larsson","doi":"10.1007/s10822-023-00541-1","DOIUrl":"10.1007/s10822-023-00541-1","url":null,"abstract":"<div><p>Theoretical predictions of the solubilizing capacity of micelles and vesicles present in intestinal fluid are important for the development of new delivery techniques and bioavailability improvement. A balance between accuracy and computational cost is a key factor for an extensive study of numerous compounds in diverse environments. In this study, we aimed to determine an optimal molecular dynamics (MD) protocol to evaluate small-molecule interactions with micelles composed of bile salts and phospholipids. MD simulations were used to produce free energy profiles for three drug molecules (danazol, probucol, and prednisolone) and one surfactant molecule (sodium caprate) as a function of the distance from the colloid center of mass. To address the challenges associated with such tasks, we compared different simulation setups, including freely assembled colloids versus pre-organized spherical micelles, full free energy profiles versus only a few points of interest, and a coarse-grained model versus an all-atom model. Our findings demonstrate that combining these techniques is advantageous for achieving optimal performance and accuracy when evaluating the solubilization capacity of micelles.</p><h3>Graphical abstract</h3><p>All-atom (AA) and coarse-grained (CG) umbrella sampling (US) simulations and point-wise free energy (FE) calculations were compared to their efficiency to computationally analyze the solubilization of active pharmaceutical ingredients in intestinal fluid colloids.</p>\u0000<div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.5,"publicationDate":"2023-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10822-023-00541-1.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138689379","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"QM assisted ML for 19F NMR chemical shift prediction","authors":"Patrick Penner,&nbsp;Anna Vulpetti","doi":"10.1007/s10822-023-00542-0","DOIUrl":"10.1007/s10822-023-00542-0","url":null,"abstract":"<div><div><h3>Background</h3><p>Ligand-observed 19F NMR detection is an efficient method for screening libraries of fluorinated molecules in fragment-based drug design campaigns. Screening fluorinated molecules in large mixtures makes 19F NMR a high-throughput method. Typically, these mixtures are generated from pools of well-characterized fragments. By predicting 19F NMR chemical shift, mixtures could be generated for arbitrary fluorinated molecules facilitating for example focused screens.</p><h3>Methods</h3><p>In a previous publication, we introduced a method to predict 19F NMR chemical shift using rooted fluorine fingerprints and machine learning (ML) methods. Having observed that the quality of the prediction depends on similarity to the training set, we here propose to assist the prediction with quantum mechanics (QM) based methods in cases where compounds are not well covered by a training set.</p><h3>Results</h3><p>Beyond similarity, the performance of ML methods could be associated with individual features in compounds. A combination of both could be used as a procedure to split input data sets into those that could be predicted by ML and those that required QM processing. We could show on a proprietary fluorinated fragment library, known as LEF (Local Environment of Fluorine), and a public Enamine data set of 19F NMR chemical shifts that ML and QM methods could synergize to outperform either method individually. Models built on Enamine data, as well as model building and QM workflow tools, can be found at https://github.com/PatrickPenner/lefshift and https://github.com/PatrickPenner/lefqm.</p></div></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.5,"publicationDate":"2023-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138570832","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Open-ComBind: harnessing unlabeled data for improved binding pose prediction","authors":"Andrew T. McNutt,&nbsp;David Ryan Koes","doi":"10.1007/s10822-023-00544-y","DOIUrl":"10.1007/s10822-023-00544-y","url":null,"abstract":"<div><p>Determination of the bound pose of a ligand is a critical first step in many in silico drug discovery tasks. Molecular docking is the main tool for the prediction of non-covalent binding of a protein and ligand system. Molecular docking pipelines often only utilize the information of one ligand binding to the protein despite the commonly held hypothesis that different ligands share binding interactions when bound to the same receptor. Here we describe Open-ComBind, an easy-to-use, open-source version of the ComBind molecular docking pipeline that leverages information from multiple ligands without known bound structures to enhance pose selection. We first create distributions of feature similarities between ligand pose pairs, comparing near-native poses with all sampled docked poses. These distributions capture the likelihood of observing similar features, such as hydrogen bonds or hydrophobic contacts, in different pose configurations. These similarity distributions are then combined with a per-ligand docking score to enhance overall pose selection by 5% and 4.5% for high-affinity and congeneric series helper ligands, respectively. Open-ComBind reduces the average RMSD of ligands in our benchmark dataset by 9.0%. We provide Open-ComBind as an easy-to-use command line and Python API to increase pose prediction performance at www.github.com/drewnutt/open_combind.</p></div>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.5,"publicationDate":"2023-12-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10822-023-00544-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138552568","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Correction: Exploring DrugCentral: from molecular structures to clinical effects","authors":"Liliana Halip,&nbsp;Sorin Avram,&nbsp;Ramona Curpan,&nbsp;Ana Borota,&nbsp;Alina Bora,&nbsp;Cristian Bologa,&nbsp;Tudor I. Oprea","doi":"10.1007/s10822-023-00545-x","DOIUrl":"10.1007/s10822-023-00545-x","url":null,"abstract":"","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"38 1","pages":""},"PeriodicalIF":3.5,"publicationDate":"2023-12-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10822-023-00545-x.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138469687","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}