Computational determination of hERG-related cardiotoxicity of drug candidates

MultiCBlo: Enhancing predictions of compound-induced inhibition of cardiac ion channels with advanced multimodal learning

Predicting compound-induced inhibition of cardiac ion channels is crucial and challenging, significantly impacting cardiac drug efficacy and safety assessments. Despite the development of various computational methods for compound-induced inhibition prediction in cardiac ion channels, their performance remains limited. Most methods struggle to fuse multi-source data, relying solely on specific dataset training, leading to poor accuracy and generalization. We introduce MultiCBlo, a model that fuses multimodal information through a progressive learning approach, designed to predict compound-induced inhibition of cardiac ion channels with high accuracy. MultiCBlo employs progressive multimodal information fusion technology to integrate the compound's SMILES sequence, graph structure, and fingerprint, enhancing its representation. This is the first application of progressive multimodal learning for predicting compound-induced inhibition of cardiac ion channels, to our knowledge. The objective of this study was to predict the compound-induced inhibition of three major cardiac ion channels: hERG, Cav1.2, and Nav1.5. The results indicate that MultiCBlo significantly outperforms current models in predicting compound-induced inhibition of cardiac ion channels. We hope that MultiCBlo will facilitate cardiac drug development and reduce compound toxicity risks. Code and data are accessible at: https://github.com/taowang11/MultiCBlo. The online prediction platform is freely accessible at: https://huggingface.co/spaces/wtttt/PCICB.

CardioDPi: An explainable deep-learning model for identifying cardiotoxic chemicals targeting hERG, Cav1.2, and Nav1.5 channels

The cardiotoxic effects of various pollutants have been a growing concern in environmental and material science. These effects encompass arrhythmias, myocardial injury, cardiac insufficiency, and pericardial inflammation. Compounds such as organic solvents and air pollutants disrupt the potassium, sodium, and calcium ion channels cardiac cell membranes, leading to the dysregulation of cardiac function. However, current cardiotoxicity models have disadvantages of incomplete data, ion channels, interpretability issues, and inability of toxic structure visualization. Herein, an interpretable deep-learning model known as CardioDPi was developed, which is capable of discriminating cardiotoxicity induced by the human Ether-à-go-go-related gene (hERG) channel, sodium channel (Na_v1.5), and calcium channel (Ca_v1.5) blockade. External validation yielded promising area under the ROC curve (AUC) values of 0.89, 0.89, and 0.94 for the hERG, Na_v1.5, and Ca_v1.5 channels, respectively. The CardioDPi can be freely accessed on the web server CardioDPipredictor (http://cardiodpi.sapredictor.cn/). Furthermore, the structural characteristics of cardiotoxic compounds were analyzed and structural alerts (SAs) can be extracted using the user-friendly CardioDPi-SAdetector web service (http://cardiosa.sapredictor.cn/). CardioDPi is a valuable tool for identifying cardiotoxic chemicals that are environmental and health risks. Moreover, the SA system provides essential insights for mode-of-action studies concerning cardiotoxic compounds.

Two-Dimensional Deep Learning Frameworks for Drug-Induced Cardiotoxicity Detection

The identification of drug-induced cardiotoxicity remains a pressing challenge with far-reaching clinical and economic ramifications, often leading to patient harm and resource-intensive drug recalls. Current methodologies, including in vivo and in vitro models, have severe limitations in accurate identification of cardiotoxic substances. Pioneering a paradigm shift from these conventional techniques, our study presents two deep learning-based frameworks, STFT-CNN and SST-CNN, to assess cardiotoxicity with markedly improved accuracy and reliability. Leveraging the power of induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) as a more human-relevant cell model, we record mechanical beating signals through impedance measurements. These temporal signals were converted into enriched two-dimensional representations through advanced transformation techniques, specifically short-time Fourier transform (STFT) and synchro-squeezing transform (SST). These transformed data are fed into the proposed frameworks for comprehensive analysis, including drug type classification, concentration classification, cardiotoxicity classification, and new drug identification. Compared to traditional models like recurrent neural network (RNN) and 1-dimensional convolutional neural network (1D-CNN), SST-CNN delivered an impressive test accuracy of 98.55% in drug type classification and 99% in distinguishing cardiotoxic and noncardiotoxic drugs. Its feasibility is further highlighted with a stellar 98.5% average accuracy for classification of various concentrations, and the superiority of our proposed frameworks underscores their promise in revolutionizing drug safety assessments. With a potential for scalability, they represent a major leap in drug safety assessments, offering a pathway to more robust, efficient, and human-relevant cardiotoxicity evaluations.

Enhancing hERG Risk Assessment with Interpretable Classificatory and Regression Models

The human Ether-à-go-go-Related Gene (hERG) is a transmembrane protein that regulates cardiac action potential, and its inhibition can induce a potentially deadly cardiac syndrome. In vitro tests help identify hERG blockers at early stages; however, the high cost motivates searching for alternative, cost-effective methods. The primary goal of this study was to enhance the Pred-hERG tool for predicting hERG blockage. To achieve this, we developed new QSAR models that incorporated additional data, updated existing classificatory and multiclassificatory models, and introduced new regression models. Notably, we integrated SHAP (SHapley Additive exPlanations) values to offer a visual interpretation of these models. Utilizing the latest data from ChEMBL v30, encompassing over 14,364 compounds with hERG data, our binary and multiclassification models outperformed both the previous iteration of Pred-hERG and all publicly available models. Notably, the new version of our tool introduces a regression model for predicting hERG activity (pIC50). The optimal model demonstrated an R2 of 0.61 and an RMSE of 0.48, surpassing the only available regression model in the literature. Pred-hERG 5.0 now offers users a swift, reliable, and user-friendly platform for the early assessment of chemically induced cardiotoxicity through hERG blockage. The tool provides versatile outcomes, including (i) classificatory predictions of hERG blockage with prediction reliability, (ii) multiclassificatory predictions of hERG blockage with reliability, (iii) regression predictions with estimated pIC50 values, and (iv) probability maps illustrating the contribution of chemical fragments for each prediction. Furthermore, we implemented explainable AI analysis (XAI) to visualize SHAP values, providing insights into the contribution of each feature to binary classification predictions. A consensus prediction calculated based on the predictions of the three developed models is also present to assist the user's decision-making process. Pred-hERG 5.0 has been designed to be user-friendly, making it accessible to users without computational or programming expertise. The tool is freely available at http://predherg.labmol.com.br.

Role of Machine Learning and Artificial Intelligence in Arrhythmias and Electrophysiology

Machine learning (ML), a subset of artificial intelligence (AI) centered on machines learning from extensive datasets, stands at the forefront of a technological revolution shaping various facets of society. Cardiovascular medicine has emerged as a key domain for ML applications, with considerable efforts to integrate these innovations into routine clinical practice. Within cardiac electrophysiology, ML applications, especially in the automated interpretation of electrocardiograms, have garnered substantial attention in existing literature. However, less recognized are the diverse applications of ML in cardiac electrophysiology and arrhythmias, spanning basic science research on arrhythmia mechanisms, both experimental and computational, as well as contributions to enhanced techniques for mapping cardiac electrical function and translational research related to arrhythmia management. This comprehensive review delves into various ML applications within the scope of this journal, organized into 3 parts. The first section provides a fundamental understanding of general ML principles and methodologies, serving as a foundational resource for readers interested in exploring ML applications in arrhythmia research. The second part offers an in-depth review of studies in arrhythmia and electrophysiology that leverage ML methodologies, showcasing the broad potential of ML approaches. Each subject is thoroughly outlined, accompanied by a review of notable ML research advancements. Finally, the review delves into the primary challenges and future perspectives surrounding ML-driven cardiac electrophysiology and arrhythmias research.

Expanding Predictive Capacities in Toxicology: Insights from Hackathon-Enhanced Data and Model Aggregation

In the realm of predictive toxicology for small molecules, the applicability domain of QSAR models is often limited by the coverage of the chemical space in the training set. Consequently, classical models fail to provide reliable predictions for wide classes of molecules. However, the emergence of innovative data collection methods such as intensive hackathons have promise to quickly expand the available chemical space for model construction. Combined with algorithmic refinement methods, these tools can address the challenges of toxicity prediction, enhancing both the robustness and applicability of the corresponding models. This study aimed to investigate the roles of gradient boosting and strategic data aggregation in enhancing the predictivity ability of models for the toxicity of small organic molecules. We focused on evaluating the impact of incorporating fragment features and expanding the chemical space, facilitated by a comprehensive dataset procured in an open hackathon. We used gradient boosting techniques, accounting for critical features such as the structural fragments or functional groups often associated with manifestations of toxicity.

Benchmarking of Small Molecule Feature Representations for hERG, Nav1.5, and Cav1.2 Cardiotoxicity Prediction

In the field of drug discovery, there is a substantial challenge in seeking out chemical structures that possess desirable pharmacological, toxicological, and pharmacokinetic properties. Complications arise when drugs interfere with the functioning of cardiac ion channels, leading to serious cardiovascular consequences. The discontinuation and removal of numerous approved drugs from the market or at late development stages in the pipeline due to such inhibitory effects further highlight the urgency of addressing this issue. Consequently, the early prediction of potential blockers targeting cardiac ion channels during the drug discovery process is of paramount importance. This study introduces a deep learning framework that computationally determines the cardiotoxicity associated with the voltage-gated potassium channel (hERG), the voltage-gated calcium channel (Cav1.2), and the voltage-gated sodium channel (Nav1.5) for drug candidates. The predictive capabilities of three feature representations─molecular fingerprints, descriptors, and graph-based numerical representations─are rigorously benchmarked. Additionally, a novel training and evaluation data set framework is presented, enabling predictive model training of drug off-target cardiotoxicity using a comprehensive and large curated data set covering these three cardiac ion channels. To facilitate these predictions, a robust and comprehensive small molecule cardiotoxicity prediction tool named CToxPred has been developed. It is made available as open source under the permissive MIT license at https://github.com/issararab/CToxPred.

Repurposing FDA-approved drugs as NLRP3 inhibitors against inflammatory diseases: machine learning and molecular simulation approaches

Activation of NLRP3 (NOD-like receptor family, pyrin domain-containing protein 3) has been associated with multiple chronic pathologies, including diabetes, atherosclerosis, and rheumatoid arthritis. Moreover, histone deacetylases (HDACs), specifically HDAC6 is required for the NLRP3 inflammasome to assemble and activate. Thus, NLRP3 serves as an attractive target for the development of novel therapeutic approaches. Several companies are now attempting to develop specific modulators of the NLRP3 inflammasome, but only a handful of small molecules of NLRP3 inflammasome inhibitors, such as MCC950 and Tranilast, are currently available for clinical use. However, their use is limited due to severe side effects and short half-lives. Thus, the repurposing of FDA-approved drugs with NLRP3 inhibitory activity is needed. The present study was aimed at repurposing preexisting drugs that might act as safe and effective NLRP3 inhibitors. A library of 2,697 FDA-approved drugs was screened for binding with NLRP3 (PDB: 7ALV) using Glide (Schrödinger). The top seven FDA-approved drugs with potential binding affinities were selected based on docking scores and subjected to ADMET profiling using pkCSM and SwissADME. The binding of the ADMET-favorable FDA-approved drugs to NLRP3 was validated using MMGBSA (Prime) and Molecular Dynamics (Desmond) in the Schrödinger suite. ADMET profiling revealed that of the seven best docking drugs, empagliflozin and citicoline had good drug-likeness properties. Moreover, MMGBSA analysis and molecular dynamics demonstrated that empagliflozin and citicoline exhibited stable ligand-NLRP3 interactions in the presence of solvents. This study sheds light on the ability of various FDA-approved drugs to act as NLRP3 inhibitors.Communicated by Ramaswamy H. Sarma.

Automation of Drug Discovery through Cutting-edge In-silico Research in Pharmaceuticals: Challenges and Future Scope

The rapidity and high-throughput nature of in silico technologies make them advantageous for predicting the properties of a large array of substances. In silico approaches can be used for compounds intended for synthesis at the beginning of drug development when there is either no or very little compound available. In silico approaches can be used for impurities or degradation products. Quantifying drugs and related substances (RS) with pharmaceutical drug analysis (PDA) can also improve drug discovery (DD) by providing additional avenues to pursue. Potential future applications of PDA include combining it with other methods to make insilico predictions about drugs and RS. One possible outcome of this is a determination of the drug potential of nontoxic RS. ADME estimation, QSAR research, molecular docking, bioactivity prediction, and toxicity testing all involve impurity profiling. Before committing to DD, RS with minimal toxicity can be utilised in silico. The efficacy of molecular docking in getting a medication to market is still debated despite its refinement and improvement. Biomedical labs and pharmaceutical companies were hesitant to adopt molecular docking algorithms for drug screening despite their decades of development and improvement. Despite the widespread use of "force fields" to represent the energy exerted within and between molecules, it has been impossible to reliably predict or compute the binding affinities between proteins and potential binding medications.

Integrated Molecular Modeling and Machine Learning for Drug Design

Modern therapeutic development often involves several stages that are interconnected, and multiple iterations are usually required to bring a new drug to the market. Computational approaches have increasingly become an indispensable part of helping reduce the time and cost of the research and development of new drugs. In this Perspective, we summarize our recent efforts on integrating molecular modeling and machine learning to develop computational tools for modulator design, including a pocket-guided rational design approach based on AlphaSpace to target protein-protein interactions, delta machine learning scoring functions for protein-ligand docking as well as virtual screening, and state-of-the-art deep learning models to predict calculated and experimental molecular properties based on molecular mechanics optimized geometries. Meanwhile, we discuss remaining challenges and promising directions for further development and use a retrospective example of FDA approved kinase inhibitor Erlotinib to demonstrate the use of these newly developed computational tools.

DICTrank: The largest reference list of 1318 human drugs ranked by risk of drug-induced cardiotoxicity using FDA labeling

Drug-induced cardiotoxicity (DICT) is a leading cause of drug trial failure and discontinuation. Current drug annotations for cardiotoxicity largely focus on individual outcomes or mechanisms. Considering the broad spectrum of adverse cardiac events, we developed Drug-Induced Cardiotoxicity Rank (DICTrank) using FDA labeling and comprehensively classified 1318 human drugs into four categories: Most-DICT-Concern (n = 341), Less-DICT-Concern (n = 528), No-DICT-Concern (n = 343), and Ambiguous-DICT-Concern (n = 106). Notably, DICTrank covers diverse therapeutic categories, of which several were enriched with Most-DICT-Concern drugs, such as antineoplastic agents, sex hormones, anti-inflammatory drugs, beta-blockers, and cardiac therapy. DICTrank currently presents the largest drug list of DICT annotation, and it could contribute to the development of new approach methods, including AI models for early identification of DICT risk during drug development and beyond.

Preclinical data on morpholine (3,5-di-tertbutyl-4-hydroxyphenyl) methanone induced anxiolysis

Trimetozine is used to be indicated for the treatment of mental illnesses, particularly anxiety. The present study provides data on the pharmacological profile of trimetozine derivative morpholine (3,5-di-tert-butyl-4-hydroxyphenyl) methanone (LQFM289) which was designed from molecular hybridization of trimetozine lead compound and 2,6-di-tert-butyl-hydroxytoluene to develop new anxiolytic drugs. Here, we conduct molecular dynamics simulations, docking studies, receptor binding assays, and in silico ADMET profiling of LQFM289 before its behavioral and biochemical assessment in mice within the dose range of 5-20 mg/kg. The docking of LQFM289 showed strong interactions with the benzodiazepine binding sites and matched well with receptor binding data. With the ADMET profile of this trimetozine derivative that predicts a high intestinal absorption and permeability to blood-brain barrier without being inhibited by the permeability glycoprotein, the oral administration of LQFM289 10 mg/kg consistently induced anxiolytic-like behavior of the mice exposed to the open field and light-dark box apparatus without eliciting motor incoordination in the wire, rotarod, and chimney tests. A decrease in the wire and rotarod´s fall latency coupled with an increase in the chimney test´s climbing time and a decrease in the number of crossings in the open field apparatus at the dose of 20 mg/kg of this trimetozine derivative suggest sedative or motor coordination impairment at this highest dose. The attenuation of the anxiolytic-like effects of LQFM289 (10 mg/kg) by flumazenil pretreatment implicates the participation of benzodiazepine binding sites. The lowering of corticosterone and tumor necrosis factor alpha (cytokine) in LQFM289-treated mice at a single oral (acute) dose of 10 mg/kg suggests that the anxiolytic-like effect of this compound also involves the recruitment of non-benzodiazepine binding sites/GABAergic molecular machinery.

Review of machine learning and deep learning models for toxicity prediction

The ever-increasing number of chemicals has raised public concerns due to their adverse effects on human health and the environment. To protect public health and the environment, it is critical to assess the toxicity of these chemicals. Traditional in vitro and in vivo toxicity assays are complicated, costly, and time-consuming and may face ethical issues. These constraints raise the need for alternative methods for assessing the toxicity of chemicals. Recently, due to the advancement of machine learning algorithms and the increase in computational power, many toxicity prediction models have been developed using various machine learning and deep learning algorithms such as support vector machine, random forest, k-nearest neighbors, ensemble learning, and deep neural network. This review summarizes the machine learning- and deep learning-based toxicity prediction models developed in recent years. Support vector machine and random forest are the most popular machine learning algorithms, and hepatotoxicity, cardiotoxicity, and carcinogenicity are the frequently modeled toxicity endpoints in predictive toxicology. It is known that datasets impact model performance. The quality of datasets used in the development of toxicity prediction models using machine learning and deep learning is vital to the performance of the developed models. The different toxicity assignments for the same chemicals among different datasets of the same type of toxicity have been observed, indicating benchmarking datasets is needed for developing reliable toxicity prediction models using machine learning and deep learning algorithms. This review provides insights into current machine learning models in predictive toxicology, which are expected to promote the development and application of toxicity prediction models in the future.

Derivation of Highly Predictive 3D-QSAR Models for hERG Channel Blockers Based on the Quantum Artificial Neural Network Algorithm

The hERG potassium channel serves as an annexed target for drug discovery because the associated off-target inhibitory activity may cause serious cardiotoxicity. Quantitative structure-activity relationship (QSAR) models were developed to predict inhibitory activities against the hERG potassium channel, utilizing the three-dimensional (3D) distribution of quantum mechanical electrostatic potential (ESP) as the molecular descriptor. To prepare the optimal atomic coordinates of dataset molecules, pairwise 3D structural alignments were carried out in order for the quantum mechanical cross correlation between the template and other molecules to be maximized. This alignment method stands out from the common atom-by-atom matching technique, as it can handle structurally diverse molecules as effectively as chemical derivatives that share an identical scaffold. The alignment problem prevalent in 3D-QSAR methods was ameliorated substantially by dividing the dataset molecules into seven subsets, each of which contained molecules with similar molecular weights. Using an artificial neural network algorithm to find the functional relationship between the quantum mechanical ESP descriptors and the experimental hERG inhibitory activities, highly predictive 3D-QSAR models were derived for all seven molecular subsets to the extent that the squared correlation coefficients exceeded 0.79. Given their simplicity in model development and strong predictability, the 3D-QSAR models developed in this study are expected to function as an effective virtual screening tool for assessing the potential cardiotoxicity of drug candidate molecules.

A computational investigation of thymidylate synthase inhibitors through a combined approach of 3D-QSAR and pharmacophore modelling

Thymidylate synthase (TS) is a crucial target of cancer drug discovery and is mainly involved in the De novo synthesis of the DNA precursor thymine. In the present study, to generate reliable models and identify a few promising molecules, we combined QSAR modelling with the pharmacophore hypothesis-generating technique. Input molecules were clustered on their similarity, and a cluster of 74 molecules with a pyrimidine moiety was chosen as the set for 3D-QSAR and pharmacophore modelling. Atom-based and field-based 3D-QSAR models were generated and statistically validated with R2 > 0.90 and Q2 > 0.75. The common pharmacophore hypothesis(CPH) generation identified the best six-point model ADHRRR. Using these best models, a library of FDA-approved drugs was screened for activity and filtered via molecular docking, ADME profiling, and molecular dynamics simulations. The top ten promising TS-inhibiting candidates were identified, and their chemical features profitable for TS inhibitors were explored.Communicated by Ramaswamy H. Sarma.

Synthesis, Computational, and Anticancer In Vitro Investigations of Aminobenzylnaphthols Derived from 2-Naphtol, Benzaldehydes, and α-Aminoacids via the Betti Reaction

Multicomponent reactions have emerged as an important approach for the synthesis of diverse and complicated chemical compounds. They have various advantages over two-component reactions, including the convenience of one-pot procedures and the ability to modify the structure of agents. Here, we employed in vitro and in silico studies to explore the anticancer potential of novel aminobenzylnaphthols derived from the Betti reaction (MMZ compounds). MTT assay was used to explore the cytotoxic activity of the compounds in pancreatic (BxPC-3 cells) and colorectal (HT-29) cancer cell lines or normal human lung fibroblasts (WI-38 cells). Proapoptotic properties of two derivatives MMZ-45AA and MMZ-140C were explored using AO/EB and annexin V-FITC/PI staining. In silico studies including ADMET profiling, molecular target prediction, docking, and dynamics were employed. The compounds exhibited cytotoxic properties and showed proapoptotic properties in respective IC50 concentrations. As indicated by in silico investigations, anticancer activity of MMZs can be attributed to the inhibition of ADORA1, CDK2, and TRIM24. Furthermore, compounds exhibited favorable ADMET properties. MMZs constitute an interesting scaffold for the potential development of new anticancer agents.

Predictive Models Based on Molecular Images and Molecular Descriptors for Drug Screening

Various toxicity and pharmacokinetic evaluations as screening experiments are needed at the drug discovery stage. Currently, to reduce the use of animal experiments and developmental expenses, the development of high-performance predictive models based on quantitative structure-activity relationship analysis is desired. From these evaluation targets, we selected 50% lethal dose (LD50), blood-brain barrier penetration (BBBP), and the clearance (CL) pathway for this investigation and constructed predictive models for each target using 636-11,886 compounds. First, we constructed predictive models using the DeepSnap-deep learning (DL) method and images of compounds as features. The calculated area under the curve (AUC) and balanced accuracy (BAC) were, respectively, 0.887 and 0.818 for LD50, 0.893 and 0.824 for BBBP, and 0.883 and 0.763 for the CL pathway. Next, molecular descriptors (MDs) of compounds were calculated using Molecular Operating Environment, alvaDesc, and ADMET Predictor to construct predictive models using the MD-based method. Using these MDs, we constructed predictive models using DataRobot. The calculated AUC and BAC were, respectively, 0.931 and 0.805 for LD50, 0.919 and 0.849 for BBBP, and 0.900 and 0.807 for the CL pathway. In this investigation, we constructed predictive models combining the DeepSnap-DL and MD-based methods. In ensemble models using the mean predictive probability of the DeepSnap-DL and MD-based methods, the calculated AUC and BAC were, respectively, 0.942 and 0.842 for LD50, 0.936 and 0.853 for BBBP, and 0.908 and 0.832 for the CL pathway, with improved predictive performance observed for all variables compared with either single method alone. Moreover, in consensus models that adopted only compounds for which the results of the two methods agreed, the calculated BAC for LD50, BBBP, and the CL pathway were 0.916, 0.918, and 0.847, respectively, indicating higher predictive performance than the ensemble models for all three variables. The predictive models combining the DeepSnap-DL and MD-based methods displayed high predictive performance for LD50, BBBP, and the CL pathway. Therefore, the application of this approach to prediction targets in various drug discovery screenings is expected to accelerate drug discovery.

Design, synthesis, antiviral evaluation, and In silico studies of acrylamides targeting nsP2 from Chikungunya virus

The Togaviridae family comprises several New- and Old-World Alphaviruses that have been responsible for thousands of human illnesses, including the RNA arbovirus Chikungunya virus (CHIKV). Firstly, it was reported in Tanzania in 1952 but rapidly it spread to several countries from Europe, Asia, and the Americas. Since then, CHIKV has been circulating in diverse countries around the world, leading to increased morbidity rates. Currently, there are no FDA-approved drugs or licensed vaccines to specifically treat CHIKV infections. Thus, there is a lack of alternatives to fight against this viral disease, making it an unmet need. Structurally, CHIKV is composed of five structural proteins (E3, E2, E1, C, and 6k) and four non-structural proteins (nsP1-4), in which nsP2 represents an attractive antiviral target for designing novel inhibitors since it has an essential role in the virus replication and transcription. Herein, we used a rational drug design strategy to select some acrylamide derivatives to be synthesized and evaluated against CHIKV nsP2 and also screened on CHIKV-infected cells. Thus, two regions of modifications were considered for these types of inhibitors, based on a previous study of our group, generating 1560 possible inhibitors. Then, the 24 most promising ones were synthesized and screened by using a FRET-based enzymatic assay protocol targeting CHIKV nsP2, identifying LQM330, 333, 336, and 338 as the most potent inhibitors, with Ki values of 48.6 ± 2.8, 92.3 ± 1.4, 2.3 ± 1.5, and 181.8 ± 2.5 μM, respectively. Still, their Km and Vmax kinetic parameters were also determined, along with their competitive binding modes of CHIKV nsP2 inhibition. Then, ITC analyses revealed KD values of 127, 159, 198, and 218 μM for LQM330, 333, 336, and 338, respectively. Also, their ΔH, ΔS, and ΔG physicochemical parameters were determined. MD simulations demonstrated that these inhibitors present a stable binding mode with nsP2, interacting with important residues of this protease, according to docking analyzes. Moreover, MM/PBSA calculations displayed that van der Waals interactions are mainly responsible for stabilizing the inhibitor-nsP2 complex, and their binding energies corroborated with their Ki values, having -198.7 ± 15.68, -124.8 ± 17.27, -247.4 ± 23.78, and -100.6 ± 19.21 kcal/mol for LQM330, 333, 336, and 338, respectively. Since Sindbis (SINV) nsP2 is similar to CHIKV nsP2, these best inhibitors were screened against SINV-infected cells, and it was verified that LQM330 presented the best result, with an EC50 value of 0.95 ± 0.09 μM. Even at 50 μM concentration, LQM338 was found to be cytotoxic on Vero cells after 48 h. Then, LQM330, 333, and 336 were evaluated against CHIKV-infected cells in antiviral assays, in which LQM330 was found to be the most promising antiviral candidate in this study, exhibiting an EC50 value of 5.2 ± 0.52 μM and SI of 31.78. The intracellular flow cytometry demonstrated that LQM330 is able to reduce the CHIKV cytopathogenic effect on cells, and also reduce the percentage of CHIKV-positive cells from 66.1% ± 7.05 to 35.8% ± 5.78 at 50 μM concentration. Finally, qPCR studies demonstrated that LQM330 was capable of reducing the number of viral RNA copies/μL, suggesting that CHIKV nsP2 is targeted by this inhibitor as its mechanism of action.

Preclinical Testing Techniques: Paving the Way for New Oncology Screening Approaches

Prior to clinical trials, preclinical testing of oncology drug candidates is performed by evaluating drug candidates with in vitro and in vivo platforms. For in vivo testing, animal models are used to evaluate the toxicity and efficacy of drug candidates. However, animal models often display poor translational results as many drugs that pass preclinical testing fail when tested with humans, with oncology drugs exhibiting especially poor acceptance rates. The FDA Modernization Act 2.0 promotes alternative preclinical testing techniques, presenting the opportunity to use higher complexity in vitro models as an alternative to in vivo testing, including three-dimensional (3D) cell culture models. Three-dimensional tissue cultures address many of the shortcomings of 2D cultures by more closely replicating the tumour microenvironment through a combination of physiologically relevant drug diffusion, paracrine signalling, cellular phenotype, and vascularization that can better mimic native human tissue. This review will discuss the common forms of 3D cell culture, including cell spheroids, organoids, organs-on-a-chip, and 3D bioprinted tissues. Their advantages and limitations will be presented, aiming to discuss the use of these 3D models to accurately represent human tissue and as an alternative to animal testing. The use of 3D culture platforms for preclinical drug development is expected to accelerate as these platforms continue to improve in complexity, reliability, and translational predictivity.

In Vitro and In Silico Investigation of BCI Anticancer Properties and Its Potential for Chemotherapy-Combined Treatments

BACKGROUND

DUSP6 phosphatase serves as a negative regulator of MAPK kinases involved in numerous cellular processes. BCI has been identified as a potential allosteric inhibitor with anticancer activity. Our study was designed to test the anticancer properties of BCI in colon cancer cells, to characterize the effect of this compound on chemotherapeutics such as irinotecan and oxaliplatin activity, and to identify potential molecular targets for this inhibitor.

METHODS

BCI cytotoxicity, proapoptotic activity, and cell cycle distribution were investigated in vitro on three colon cancer cell lines (DLD1, HT-29, and Caco-2). In silico investigation was prepared to assess BCI drug-likeness and identify potential molecular targets.

RESULTS

The exposure of colorectal cancer cells with BCI resulted in antitumor effects associated with cell cycle arrest and induction of apoptosis. BCI exhibited strong cytotoxicity on DLD1, HT-29, and Caco-2 cells. BCI showed no significant interaction with irinotecan, but strongly attenuated the anticancer activity of oxaliplatin when administered together. Analysis of synergy potential further confirmed the antagonistic interaction between these two compounds. In silico investigation indicated CDK5 as a potential new target of BCI.

CONCLUSIONS

Our studies point to the anticancer potential of BCI but note the need for a precise mechanism of action.

hERG agonists pose challenges to web-based machine learning methods for prediction of drug-hERG channel interaction

Pharmacological blockade of the IKr channel (hERG) by diverse drugs in clinical use is associated with the Long QT Syndrome that can lead to life threatening arrhythmia. Various computational tools including machine learning models (MLM) for the prediction of hERG inhibition have been developed to facilitate the throughput screening of drugs in development and optimise thus the prediction of hERG liabilities. The use of MLM relies on large libraries of training compounds for the quantitative structure-activity relationship (QSAR) modelling of hERG inhibition. The focus on inhibition omits potential effects of hERG channel agonist molecules and their associated QT shortening risk. It is instructive, therefore, to consider how known hERG agonists are handled by MLM. Here, two highly developed online computational tools for the prediction of hERG liability, Pred-hERG and HergSPred were probed for their ability to detect hERG activator drug molecules as hERG interactors. In total, 73 hERG blockers were tested with both computational tools giving overall good predictions for hERG blockers with reported IC50s below Pred-hERG and HergSPred cut-off threshold for hERG inhibition. However, for compounds with reported IC50s above this threshold such as disopyramide or sotalol discrepancies were observed. HergSPred identified all 20 hERG agonists selected as interacting with the hERG channel. Further studies are warranted to improve online MLM prediction of hERG related cardiotoxicity, by explicitly taking into account channel agonism as well as inhibition.

Computational modeling and synthesis of pyridine variants of benzoyl-phenoxy-acetamide with high glioblastoma cytotoxicity and brain tumor penetration

Glioblastomas are highly aggressive brain tumors for which therapeutic options are very limited. In a quest for new anti-glioblastoma drugs, we focused on specific structural modifications to the benzoyl-phenoxy-acetamide (BPA) structure present in a common lipid-lowering drug, fenofibrate, and in our first prototype glioblastoma drug, PP1. Here, we propose extensive computational analyses to improve the selection of the most effective glioblastoma drug candidates. Initially, over 100 structural BPA variations were analyzed and their physicochemical properties, such as water solubility (- logS), calculated partition coefficient (ClogP), probability for BBB crossing (BBB_SCORE), probability for CNS penetration (CNS-MPO) and calculated cardiotoxicity (hERG), were evaluated. This integrated approach allowed us to select pyridine variants of BPA that show improved BBB penetration, water solubility, and low cardiotoxicity. Herein the top 24 compounds were synthesized and analyzed in cell culture. Six of them demonstrated glioblastoma toxicity with IC50 ranging from 0.59 to 3.24 µM. Importantly, one of the compounds, HR68, accumulated in the brain tumor tissue at 3.7 ± 0.5 µM, which exceeds its glioblastoma IC50 (1.17 µM) by over threefold.

PredPS: Attention-based graph neural network for predicting stability of compounds in human plasma

Stability of compounds in the human plasma is crucial for maintaining sufficient systemic drug exposure and considered an essential factor in the early stages of drug discovery and development. The rapid degradation of compounds in the plasma can result in poor in vivo efficacy. Currently, there are no open-source software programs for predicting human plasma stability. In this study, we developed an attention-based graph neural network, PredPS to predict the plasma stability of compounds in human plasma using in-house and open-source datasets. The PredPS outperformed the two machine learning and two deep learning algorithms that were used for comparison indicating its stability-predicting efficiency. PredPS achieved an area under the receiver operating characteristic curve of 90.1%, accuracy of 83.5%, sensitivity of 82.3%, and specificity of 84.6% when evaluated using 5-fold cross-validation. In the early stages of drug discovery, PredPS could be a helpful method for predicting the human plasma stability of compounds. Saving time and money can be accomplished by adopting an in silico-based plasma stability prediction model at the high-throughput screening stage. The source code for PredPS is available at https://bitbucket.org/krict-ai/predps and the PredPS web server is available at https://predps.netlify.app.

Artificial Intelligence in Drug Toxicity Prediction: Recent Advances, Challenges, and Future Perspectives

Toxicity prediction is a critical step in the drug discovery process that helps identify and prioritize compounds with the greatest potential for safe and effective use in humans, while also reducing the risk of costly late-stage failures. It is estimated that over 30% of drug candidates are discarded owing to toxicity. Recently, artificial intelligence (AI) has been used to improve drug toxicity prediction as it provides more accurate and efficient methods for identifying the potentially toxic effects of new compounds before they are tested in human clinical trials, thus saving time and money. In this review, we present an overview of recent advances in AI-based drug toxicity prediction, including the use of various machine learning algorithms and deep learning architectures, of six major toxicity properties and Tox21 assay end points. Additionally, we provide a list of public data sources and useful toxicity prediction tools for the research community and highlight the challenges that must be addressed to enhance model performance. Finally, we discuss future perspectives for AI-based drug toxicity prediction. This review can aid researchers in understanding toxicity prediction and pave the way for new methods of drug discovery.

A biosensing system employing nonlinear dynamic analysis-assisted neural network for drug-induced cardiotoxicity assessment

Preclinical investigation of drug-induced cardiotoxicity is of importance for drug development. To evaluate such cardiotoxicity, in vitro high-throughput interdigitated electrode-based recording of cardiomyocytes mechanical beating is widely used. To automatically analyze the features from the beating signals for drug-induced cardiotoxicity assessment, artificial neural network analysis is conventionally employed and signals are segmented into cycles and feature points are located in the cycles. However, signal segmentation and location of feature points for different signal shapes require design of specific algorithms. Consequently, this may lower the efficiency of research and the applications of such algorithms in signals with different morphologies are limited. Here, we present a biosensing system that employs nonlinear dynamic analysis-assisted neural network (NDANN) to avoid the signal segmentation process and directly extract features from beating signal time series. By processing beating time series with fixed time duration to avoid the signal segmentation process, this NDANN-based biosensing system can identify drug-induced cardiotoxicity with accuracy over 0.99. The individual drugs were classified with high accuracies over 0.94 and drug-induced cardiotoxicity levels were accurately predicted. We also evaluated the generalization performance of the NDANN-based biosensing system in assessing drug-induced cardiotoxicity through an independent dataset. This system achieved accuracy of 0.85-0.95 for different drug concentrations in identification of drug-induced cardiotoxicity. This result demonstrates that our NDANN-based biosensing system has the capacity of screening newly developed drugs, which is crucial in practical applications. This NDANN-based biosensing system can work as a new screening platform for drug-induced cardiotoxicity and improve the efficiency of bio-signal processing.

Screening of phytochemicals as potential anti-breast cancer agents targeting HER2: an in-silico approach

Breast cancer is the most common cancer among women around the world. Human Epidermal growth factor Receptor-2 (HER2) is a membrane tyrosine kinase overexpressed in 30% of human breast cancers; thus, it serves as an important drug target. Currently available HER2 inhibitor lapatinib targets the ATP binding site of the cytoplasmic kinase domain, blocking autophosphorylation and activation of HER-2. However, it causes side effects like diarrhea, nausea, rash and possible liver toxicity. As phytochemicals have fewer side effects and are relatively affordable, they offer an effective alternative. Hence, we aimed to identify potential phytochemicals that could act as HER2 inhibitors employing computational methods such as molecular docking, molecular dynamic simulation, and ADMET prediction. Out of 1500 phytochemicals docked to the ATP binding site of the HER2 kinase domain, luxenchalcone, rhinacanthin Q, subtrifloralacton D, and 7,7″-dimethyllanaraflavone exhibited higher binding affinity than the reference inhibitor and satisfied the Lipinski's rule of five. Analysis of molecular dynamics simulation trajectory showed that Rhinacanthin Q, subtrifloralacton D, and 7,7″-dimethyllanaraflavone formed a stable and compact complex without vast conformational fluctuations. MM/PBSA binding free energy analysis revealed that Rhinacanthin Q, subtrifloralacton D, and 7,7″-dimethyllanaraflavone have high binding affinity to HER2. Therefore, Rhinacanthin Q, subtrifloralacton D, and 7,7″-dimethyllanaraflavone could be potential bioactive molecules to act as inhibitor of HER2 protein. Eventually, experimental studies are needed to evaluate the potentials of these phytochemicals further. The development of drug for HER2 positive breast cancer could be accelerated with the findings of our research. Communicated by Ramaswamy H. Sarma.

Investigating cardiotoxicity related with hERG channel blockers using molecular fingerprints and graph attention mechanism

Human ether-a-go-go-related gene (hERG) channel blockade by small molecules is a big concern during drug development in the pharmaceutical industry. Failure or inhibition of hERG channel activity caused by drug molecules can lead to prolonging QT interval, which will result in serious cardiotoxicity. Thus, evaluating the hERG blocking activity of all these small molecular compounds is technically challenging, and the relevant procedures are expensive and time-consuming. In this study, we develop a novel deep learning predictive model named DMFGAM for predicting hERG blockers. In order to characterize the molecule more comprehensively, we first consider the fusion of multiple molecular fingerprint features to characterize its final molecular fingerprint features. Then, we use the multi-head attention mechanism to extract the molecular graph features. Both molecular fingerprint features and molecular graph features are fused as the final features of the compounds to make the feature expression of compounds more comprehensive. Finally, the molecules are classified into hERG blockers or hERG non-blockers through the fully connected neural network. We conduct 5-fold cross-validation experiment to evaluate the performance of DMFGAM, and verify the robustness of DMFGAM on external validation datasets. We believe DMFGAM can serve as a powerful tool to predict hERG channel blockers in the early stages of drug discovery and development.

In Vivo and In Silico Analgesic Activity of Ficus populifolia Extract Containing 2-O-β-D-(3',4',6'-Tri-acetyl)-glucopyranosyl-3-methyl Pentanoic Acid

Natural product-based structural templates have immensely shaped small molecule drug discovery, and new biogenic natural products have randomly provided the leads and molecular targets in anti-analgesic activity spheres. Pain relief achieved through opiates and non-steroidal anti-inflammatory drugs (NSAIDs) has been under constant scrutiny owing to their tolerance, dependency, and other organs toxicities and tissue damage, including harm to the gastrointestinal tract (GIT) and renal tissues. A new, 3',4',6'-triacetylated-glucoside, 2-O-β-D-(3',4',6'-tri-acetyl)-glucopyranosyl-3-methyl pentanoic acid was obtained from Ficus populifolia, and characterized through a detailed NMR spectroscopic analysis, i.e., ¹H-NMR, ¹³C-DEPT-135, and the 2D nuclear magnetic resonance (NMR) correlations. The product was in silico investigated for its analgesic prowess, COX-2 binding feasibility and scores, drug likeliness, ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties, possible biosystem's toxicity using the Discovery Studio^®, and other molecular studies computational software programs. The glycosidic product showed strong potential as an analgesic agent. However, an in vivo evaluation, though at strong levels of pain-relieving action, was estimated on the compound's extract owing to the quantity and yield issues of the glycosidic product. Nonetheless, the F. populifolia extract showed the analgesic potency in eight-week-old male mice on day seven of the administration of the extract's dose in acetic acid-induced writhing and hot-plate methods. Acetic acid-induced abdominal writhing for all the treated groups decreased significantly (p < 0.0001), as compared to the control group (n = 6) by 62.9%, 67.9%, and 70.9% of a dose of 100 mg/kg (n = 6), 200 mg/kg (n = 6), and 400 mg/kg (n = 6), respectively. Similarly, using the analgesia meter, the reaction time to pain sensation increased significantly (p < 0.0001), as compared to the control (n = 6). The findings indicated peripheral and central-nervous-system-mediated analgesic action of the product obtained from the corresponding extract.

Ensemble of structure and ligand-based classification models for hERG liability profiling

Drug-induced cardiotoxicity represents one of the most critical safety concerns in the early stages of drug development. The blockade of the human ether-à-go-go-related potassium channel (hERG) is the most frequent cause of cardiotoxicity, as it is associated to long QT syndrome which can lead to fatal arrhythmias. Therefore, assessing hERG liability of new drugs candidates is crucial to avoid undesired cardiotoxic effects. In this scenario, computational approaches have emerged as useful tools for the development of predictive models able to identify potential hERG blockers. In the last years, several efforts have been addressed to generate ligand-based (LB) models due to the lack of experimental structural information about hERG channel. However, these methods rely on the structural features of the molecules used to generate the model and often fail in correctly predicting new chemical scaffolds. Recently, the 3D structure of hERG channel has been experimentally solved enabling the use of structure-based (SB) strategies which may overcome the limitations of the LB approaches. In this study, we compared the performances achieved by both LB and SB classifiers for hERG-related cardiotoxicity developed by using Random Forest algorithm and employing a training set containing 12789 hERG binders. The SB models were trained on a set of scoring functions computed by docking and rescoring calculations, while the LB classifiers were built on a set of physicochemical descriptors and fingerprints. Furthermore, models combining the LB and SB features were developed as well. All the generated models were internally validated by ten-fold cross-validation on the TS and further verified on an external test set. The former revealed that the best performance was achieved by the LB model, while the model combining the LB and the SB attributes displayed the best results when applied on the external test set highlighting the usefulness of the integration of LB and SB features in correctly predicting unseen molecules. Overall, our predictive models showed satisfactory performances providing new useful tools to filter out potential cardiotoxic drug candidates in the early phase of drug discovery.

A deep learning method for predicting molecular properties and compound-protein interactions

Predicting molecular properties and compound-protein interactions (CPIs) are two important areas of drug design and discovery. They are also an essential way to discover lead compounds in virtual screening. Recently, in silico methods based on deep learning have demonstrated excellent performance in various challenges. It is imperative to develop efficient computational methods to predict accurately both molecular properties and CPIs in drug research using deep learning techniques. In this paper, we propose a deep learning method applicable to both molecular property prediction and CPI prediction based on the idea that both are generally influenced by chemical structure and sequence information of compounds and proteins. Molecular properties are inferred by integrating the molecular structure and sequence information of compounds, and CPIs are predicted by integrating protein sequence and compound structure. The method combines topological structure and sequence fingerprint information of molecules, extracts adequately raw data features, and generates highly representative features for prediction. Molecular property prediction experiments were conducted on BACE, P53 and hERG datasets, and CPI prediction experiments were conducted on Human, C. elegans and KIBA datasets. MG-S achieves outperformance in molecular property prediction on P53, the differences in AUC, Precision and MCC are 0.030, 0.050 and 0.100, respectively, over the suboptimal baseline model, and provides consistently good results on BACE and hERG.The model also achieves impressive performance in CPI prediction, the differences in AUC, Precision and MCC on KIBA are 0.141, 0.138, 0.090 and 0.082, respectively, compared with the state-of-the-art models. The comprehensive results show that the MG-S model has higher performance, better classification ability, and faster convergence. MG-S will serve as a useful method to predict compound properties and CPIs in the early stages of drug design and discovery.Our code and datasets are available at: https://github.com/happay-ending/cpi_cpp.

Prediction of drug-induced liver injury and cardiotoxicity using chemical structure and in vitro assay data

Drug-induced liver injury (DILI) and cardiotoxicity (DICT) are major adverse effects triggered by many clinically important drugs. To provide an alternative to in vivo toxicity testing, the U.S. Tox21 consortium has screened a collection of ∼10K compounds, including drugs in clinical use, against >70 cell-based assays in a quantitative high-throughput screening (qHTS) format. In this study, we compiled reference compound lists for DILI and DICT and compared the potential of Tox21 assay data with chemical structure information in building prediction models for human in vivo hepatotoxicity and cardiotoxicity. Models were built with four different machine learning algorithms (e.g., Random Forest, Naïve Bayes, eXtreme Gradient Boosting, and Support Vector Machine) and model performance was evaluated by calculating the area under the receiver operating characteristic curve (AUC-ROC). Chemical structure-based models showed reasonable predictive power for DILI (best AUC-ROC = 0.75 ± 0.03) and DICT (best AUC-ROC = 0.83 ± 0.03), while Tox21 assay data alone only showed better than random performance. DILI and DICT prediction models built using a combination of assay data and chemical structure information did not have a positive impact on model performance. The suboptimal predictive performance of the assay data is likely due to insufficient coverage of an adequately predictive number of toxicity mechanisms. The Tox21 consortium is currently expanding coverage of biological response space with additional assays that probe toxicologically important targets and under-represented pathways that may improve the prediction of in vivo toxicity such as DILI and DICT.

Intersection of stem cell biology and engineering towards next generation in vitro models of human fibrosis

Human fibrotic diseases constitute a major health problem worldwide. Fibrosis involves significant etiological heterogeneity and encompasses a wide spectrum of diseases affecting various organs. To date, many fibrosis targeted therapeutic agents failed due to inadequate efficacy and poor prognosis. In order to dissect disease mechanisms and develop therapeutic solutions for fibrosis patients, in vitro disease models have gone a long way in terms of platform development. The introduction of engineered organ-on-a-chip platforms has brought a revolutionary dimension to the current fibrosis studies and discovery of anti-fibrotic therapeutics. Advances in human induced pluripotent stem cells and tissue engineering technologies are enabling significant progress in this field. Some of the most recent breakthroughs and emerging challenges are discussed, with an emphasis on engineering strategies for platform design, development, and application of machine learning on these models for anti-fibrotic drug discovery. In this review, we discuss engineered designs to model fibrosis and how biosensor and machine learning technologies combine to facilitate mechanistic studies of fibrosis and pre-clinical drug testing.

cardioToxCSM: A Web Server for Predicting Cardiotoxicity of Small Molecules

The design of novel, safe, and effective drugs to treat human diseases is a challenging venture, with toxicity being one of the main sources of attrition at later stages of development. Failure due to toxicity incurs a significant increase in costs and time to market, with multiple drugs being withdrawn from the market due to their adverse effects. Cardiotoxicity, for instance, was responsible for the failure of drugs such as fenspiride, propoxyphene, and valdecoxib. While significant effort has been dedicated to mitigate this issue by developing computational approaches that aim to identify molecules likely to be toxic, including quantitative structure-activity relationship models and machine learning methods, current approaches present limited performance and interpretability. To overcome these, we propose a new web-based computational method, cardioToxCSM, which can predict six types of cardiac toxicity outcomes, including arrhythmia, cardiac failure, heart block, hERG toxicity, hypertension, and myocardial infarction, efficiently and accurately. cardioToxCSM was developed using the concept of graph-based signatures, molecular descriptors, toxicophore matchings, and molecular fingerprints, leveraging explainable machine learning, and was validated internally via different cross validation schemes and externally via low-redundancy blind sets. The models presented robust performances with areas under ROC curves of up to 0.898 on 5-fold cross-validation, consistent with metrics on blind tests. Additionally, our models provide interpretation of the predictions by identifying whether substructures that are commonly enriched in toxic compounds were present. We believe cardioToxCSM will provide valuable insight into the potential cardiotoxicity of small molecules early on drug screening efforts. The method is made freely available as a web server at https://biosig.lab.uq.edu.au/cardiotoxcsm.

Pyridine-N-Oxide Alkaloids from Allium stipitatum and Their Synthetic Disulfide Analogs as Potential Drug Candidates against Mycobacterium tuberculosis: A Molecular Docking, QSBAR, and ADMET Prediction Approach

In this study, we consider pyridine-N-oxide alkaloids from Allium stipitatum and their synthetic disulfide analogs (PDAs) as candidates for next-generational antimycobacterial agents, in light of growing resistance to existing conventional therapies. In silico studies involving molecular docking simulations of 12 PDAs were carried out against 7 Mycobacterium tuberculosis target proteins (MTs) to determine their theoretical binding affinities. Compounds A3, A6, and B9 demonstrated stronger binding affinities on similar MTs. Molecular descriptors (MDs) describing structural and physicochemical properties of the compounds were also calculated using ChemDes, explored using Pearson's correlation analysis, and principal component analysis (PCA) in comparison with MDs from conventional antitubercular medicines. The PDAs possessed similar scores as isoniazid and pyrazinamide. The MDs were also used to conduct a quantitative structure-binding affinity relationship (QSBAR) study by building good fit and significant models through principal component regression (PCR) and partial least squares regression (PLSR). Leave-one-out cross-validation was adopted in the PLSR, resulting in good predictive models on all MTs (range of R² = 0.7541‐0.8992; range of Q² = 0.6183‐0.8162). Both PCR and PLSR models predicted the significant effects of ndonr, Hy, Mol wt, nhev, nring, ndb, Log P, W, Pol, ISIZ, TIAC, Getov, and UI on the binding of ligands to the MTs. In silico prediction of PDAs' ADMET profiles was conducted with QikProp utility. The ADMET profiles of the compounds were favorable. The outcome of the current study strengthens the significance of these compounds as promising lead candidates for the treatment of multidrug-resistant tuberculosis.

Identification of hydantoin based Decaprenylphosphoryl-β-D-Ribose Oxidase (DprE1) inhibitors as antimycobacterial agents using computational tools

Tuberculosis (TB) is one of the emerging infectious diseases in the world. DprE1 (Decaprenylphosphoryl-β-D-ribose 2'-epimerase), an enzyme accountable for mycobacterial cell wall synthesis was the first drug gable target based on discoveries of inhibitors via HTS (high throughput screening). Since then, many literature reports have been published so far enlightening varieties of chemical scaffolds acting as inhibitors of DprE1. Herein, in our present study, we have developed statistically robust GA-MLR (genetic algorithm multiple linear regression), atom-based as well as field based-3D-QSAR models. Both atom-based as well as field based-3D-QSAR models (internally as well as externally validated) were obtained with robust Training set, R² > 0.69 and Test set, Q² > 0.50. We have also developed top ranked 5 point hypothesis AAAHR_1 among 14 CPHs (common pharmacophore hypotheses). We found that our dataset molecule had more docking score (XP mode = - 9.068 kcal/mol) than the standards isoniazid and ethambutol; when docked into binding pockets of enzyme 4P8C with Glide module. We further queried our best docked dataset molecule 151 for ligand based virtual screening using "SwissSimilarity" platform. Among 9 identified hits, we found ZINC12196803 had best binding energies and docking score (docking score = - 9.437 kcal/mol, MMGBSA dgBind = - 70.508 kcal/mol). Finally, our molecular dynamics studies for 1.2-100 ns depicts that these complexes are stable. We have also carried out in-silico ADMET predictions, Cardiac toxicity, 'SwissTargetPredictions' and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) binding energy calculations for further explorations of dataset as well as hit molecules. Our current studies showed that the hit molecule ZINC12196803 may enlighten the path for future developments of DprE1 inhibitors.

In Silico Pharmacokinetic Profiling of the Identified Bioactive Metabolites of Pergularia tomentosa L. Latex Extract and In Vitro Cytotoxic Activity via the Induction of Caspase-Dependent Apoptosis with S-Phase Arrest

The in vitro cytotoxic efficacy of plant latex from Pergularia tomentosa L. was studied using five human cancer cell lines: HeLa cells (cervical carcinoma cells), A-549 (lung carcinoma), Panc-1 (pancreatic carcinoma cells), MDA-MB-231 (metastatic mammary adenocarcinoma), and MRC-5 (lung fibroblast cell line) cells. The phytonutrient content of plant latex was identified using the liquid chromatography/mass spectra-quadrupole time of flight (LC/MS-QTOF) technique. In silico studies of polyphenols were carried out to clarify the potential mode of action of the plant latex’s constituents. The treatment of different tumor cell lines with different concentrations of plant latex revealed a potent efficacy on the human lung carcinoma cell line (A-549) (IC₅₀ = 3.89 µg/mL) compared with that with vinblastine as a positive control (IC₅₀ = 7.12 µg/mL). The effect of the potent concentration of plant latex on the A-549 cell line induced cell arrest, upregulated the expression of pre-apoptotic markers, and downregulated the expression of antiapoptotic markers. Seven identified polyphenols were selected for the in silico study. A docking assessment using the epidermal growth factor receptor kinase (EGFRk) and eltronib as a positive control showed a higher affinity for the enzyme receptor of the selected polyphenols, except for methyl orsellinate and ginkgotoxin. The ADMET assessment demonstrated the inhibitory effect of the polyphenols on CYP450, except for ouabagenin and xanthyletine. The selected polyphenols obey Lipinski’s drug-likeness with no significant toxicity effect. In conclusion, the plant latex of P. tomentosa L. showed cytotoxic activity on the A-549 cell line, and the selected polyphenols showed a promising prodrug agent with a low profile of toxicity in the study.

Ligand-based prediction of hERG-mediated cardiotoxicity based on the integration of different machine learning techniques

Drug-induced cardiotoxicity is a common side effect of drugs in clinical use or under postmarket surveillance and is commonly due to off-target interactions with the cardiac human-ether-a-go-go-related (hERG) potassium channel. Therefore, prioritizing drug candidates based on their hERG blocking potential is a mandatory step in the early preclinical stage of a drug discovery program. Herein, we trained and properly validated 30 ligand-based classifiers of hERG-related cardiotoxicity based on 7,963 curated compounds extracted by the freely accessible repository ChEMBL (version 25). Different machine learning algorithms were tested, namely, random forest, K-nearest neighbors, gradient boosting, extreme gradient boosting, multilayer perceptron, and support vector machine. The application of 1) the best practices for data curation, 2) the feature selection method VSURF, and 3) the synthetic minority oversampling technique (SMOTE) to properly handle the unbalanced data, allowed for the development of highly predictive models (BAMAX = 0.91, AUCMAX = 0.95). Remarkably, the undertaken temporal validation approach not only supported the predictivity of the herein presented classifiers but also suggested their ability to outperform those models commonly used in the literature. From a more methodological point of view, the study put forward a new computational workflow, freely available in the GitHub repository (https://github.com/PDelre93/hERG-QSAR), as valuable for building highly predictive models of hERG-mediated cardiotoxicity.

hERG Blockade Prediction by Combining Site Identification by Ligand Competitive Saturation and Physicochemical Properties

Human ether-a-go-go-related gene (hERG) potassium channel is well-known contributor to drug-induced cardiotoxicity and therefore an extremely important target when performing safety assessments of drug candidates. Ligand-based approaches in connection with quantitative structure active relationships (QSAR) analyses have been developed to predict hERG toxicity. Availability of the recent published cryogenic electron microscopy (cryo-EM) structure for the hERG channel opened the prospect for using structure-based simulation and docking approaches for hERG drug liability predictions. In recent time, the idea of combining structure- and ligand-based approaches for modeling hERG drug liability has gained momentum offering improvements in predictability when compared to ligand-based QSAR practices alone. The present article demonstrates uniting the structure-based SILCS (site-identification by ligand competitive saturation) approach in conjunction with physicochemical properties to develop predictive models for hERG blockade. This combination leads to improved model predictability based on Pearson's R and percent correct (represents rank-ordering of ligands) metric for different validation sets of hERG blockers involving diverse chemical scaffold and wide range of pIC50 values. The inclusion of the SILCS structure-based approach allows determination of the hERG region to which compounds bind and the contribution of different chemical moieties in the compounds to blockade, thereby facilitating the rational ligand design to minimize hERG liability.

CDI Exerts Anti-Tumor Effects by Blocking the FoxM1-DNA Interaction

The Forkhead box protein M1 (FoxM1) is an appealing target for anti-cancer therapeutics as this cell proliferation-associated transcription factor is overexpressed in most human cancers. FoxM1 is involved in tumor invasion, angiogenesis, and metastasis. To discover novel inhibitors that disrupt the FoxM1-DNA interaction, we identified CDI, a small molecule that inhibits the FoxM1–DNA interaction. CDI was identified through an assay based on the time-resolved fluorescence energy transfer response of a labeled consensus oligonucleotide that was bound to a recombinant FoxM1-dsDNA binding domain (FoxM1-DBD) protein and exhibited potent inhibitory activity against FoxM1-DNA interaction. CDI suppressed cell proliferation and induced apoptosis in MDA-MB-231 cells obtained from a breast cancer patient. Furthermore, it decreased not only the mRNA and protein expression of FoxM1 but also that of downstream targets such as CDC25b. Additionally, global transcript profiling of MDA-MB-231 cells by RNA-Seq showed that CDI decreases the expression of FoxM1-regulated genes. The docking and MD simulation results indicated that CDI likely binds to the DNA interaction site of FoxM1-DBD and inhibits the function of FoxM1-DBD. These results of CDI being a possible effective inhibitor of FoxM1-DNA interaction will encourage its usage in pharmaceutical applications.

Mutation-Specific Differences in Kv7.1 (KCNQ1) and Kv11.1 (KCNH2) Channel Dysfunction and Long QT Syndrome Phenotypes

The electrocardiogram (ECG) empowered clinician scientists to measure the electrical activity of the heart noninvasively to identify arrhythmias and heart disease. Shortly after the standardization of the 12-lead ECG for the diagnosis of heart disease, several families with autosomal recessive (Jervell and Lange-Nielsen Syndrome) and dominant (Romano-Ward Syndrome) forms of long QT syndrome (LQTS) were identified. An abnormally long heart rate-corrected QT-interval was established as a biomarker for the risk of sudden cardiac death. Since then, the International LQTS Registry was established; a phenotypic scoring system to identify LQTS patients was developed; the major genes that associate with typical forms of LQTS were identified; and guidelines for the successful management of patients advanced. In this review, we discuss the molecular and cellular mechanisms for LQTS associated with missense variants in KCNQ1 (LQT1) and KCNH2 (LQT2). We move beyond the "benign" to a "pathogenic" binary classification scheme for different KCNQ1 and KCNH2 missense variants and discuss gene- and mutation-specific differences in K+ channel dysfunction, which can predispose people to distinct clinical phenotypes (e.g., concealed, pleiotropic, severe, etc.). We conclude by discussing the emerging computational structural modeling strategies that will distinguish between dysfunctional subtypes of KCNQ1 and KCNH2 variants, with the goal of realizing a layered precision medicine approach focused on individuals.

BayeshERG: a robust, reliable and interpretable deep learning model for predicting hERG channel blockers

Unintended inhibition of the human ether-à-go-go-related gene (hERG) ion channel by small molecules leads to severe cardiotoxicity. Thus, hERG channel blockage is a significant concern in the development of new drugs. Several computational models have been developed to predict hERG channel blockage, including deep learning models; however, they lack robustness, reliability and interpretability. Here, we developed a graph-based Bayesian deep learning model for hERG channel blocker prediction, named BayeshERG, which has robust predictive power, high reliability and high resolution of interpretability. First, we applied transfer learning with 300 000 large data in initial pre-training to increase the predictive performance. Second, we implemented a Bayesian neural network with Monte Carlo dropout to calibrate the uncertainty of the prediction. Third, we utilized global multihead attentive pooling to augment the high resolution of structural interpretability for the hERG channel blockers and nonblockers. We conducted both internal and external validations for stringent evaluation; in particular, we benchmarked most of the publicly available hERG channel blocker prediction models. We showed that our proposed model outperformed predictive performance and uncertainty calibration performance. Furthermore, we found that our model learned to focus on the essential substructures of hERG channel blockers via an attention mechanism. Finally, we validated the prediction results of our model by conducting in vitro experiments and confirmed its high validity. In summary, BayeshERG could serve as a versatile tool for discovering hERG channel blockers and helping maximize the possibility of successful drug discovery. The data and source code are available at our GitHub repository (https://github.com/GIST-CSBL/BayeshERG).

Novel QSAR Approach for a Regression Model of Clearance That Combines DeepSnap-Deep Learning and Conventional Machine Learning

The toxicity, absorption, distribution, metabolism, and excretion properties of some targets are difficult to predict by quantitative structure-activity relationship analysis. Therefore, there is a need for a new prediction method that performs well for these targets. The aim of this study was to develop a new regression model of rat clearance (CL). We constructed a regression model using 1545 in-house compounds for which we had rat CL data. Molecular descriptors were calculated using molecular operating environment, alvaDesc, and ADMET Predictor software. The classification model of DeepSnap and Deep Learning (DeepSnap-DL) with images of the three-dimensional chemical structures of compounds as features was constructed, and the prediction probabilities for each compound were calculated. For molecular descriptor-based methods that use molecular descriptors and conventional machine learning algorithms selected by DataRobot, the correlation coefficient (R ²) and root mean square error (RMSE) were 0.625-0.669 and 0.295-0.318, respectively. We combined molecular descriptors and prediction probability of DeepSnap-DL as features and developed a novel regression method we called the combination model. In the combination model with these two types of features and conventional algorithms selected by DataRobot, R ² and RMSE were 0.710-0.769 and 0.247-0.278, respectively. This finding shows that the combination model performed better than molecular descriptor-based methods. Our combination model will contribute to the design of more rational compounds for drug discovery. This method may be applicable not only to rat CL but also to other pharmacokinetic and pharmacological activity and toxicity parameters; therefore, applying it to other parameters may help to accelerate drug discovery.

A comprehensive exploration of pharmacological properties, bioactivities and inhibitory potentiality of luteolin from Tridax procumbens as anticancer drug by in-silico approach

Tridax procumbens is a flowering plant of the Asteraceae family with a wide range of medicinal uses like anti-inflammatory, anti-diabetic, anti-microbial, immunomodulatory, etc. This study aimed to investigate the anti-cancerous activity of human lung cancer for targeting luteolin, a phytochemical of Tridax procumbens. The computational study has been done for studying the structural properties of luteolin. The drug-likeness of the molecule has been predicted by virtual screening of ADMET properties. The molecular docking technique of the in-silico method is performed to check the complex formation between protein and ligand. The reactivity and stability of the molecule are investigated with the help of molecular dynamics (MD) simulations. In the present work, we have tried to establish a strong candidature of any of the phytochemical of Tridax Procumbens as an inhibitor against human lung cancer.

PredMS: a random forest model for predicting metabolic stability of drug candidates in human liver microsomes

MOTIVATION

Poor metabolic stability leads to drug development failure. Therefore, it is essential to evaluate the metabolic stability of small compounds for successful drug discovery and development. However, evaluating metabolic stability in vitro and in vivo is expensive, time-consuming and laborious. In addition, only a few free software programs are available for metabolic stability data and prediction. Therefore, in this study, we aimed to develop a prediction model that predicts the metabolic stability of small compounds.

RESULTS

We developed a computational model, PredMS, which predicts the metabolic stability of small compounds as stable or unstable in human liver microsomes. PredMS is based on a random forest model using an in-house database of metabolic stability data of 1917 compounds. To validate the prediction performance of PredMS, we generated external test data of 61 compounds. PredMS achieved an accuracy of 0.74, Matthew's correlation coefficient of 0.48, sensitivity of 0.70, specificity of 0.86, positive predictive value of 0.94 and negative predictive value of 0.46 on the external test dataset. PredMS will be a useful tool to predict the metabolic stability of small compounds in the early stages of drug discovery and development.

AVAILABILITY AND IMPLEMENTATION

The source code for PredMS is available at https://bitbucket.org/krictai/predms, and the PredMS web server is available at https://predms.netlify.app.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In Silico Models for Predicting Acute Systemic Toxicity

In this chapter, we give a brief overview of the regulatory requirements for acute systemic toxicity information in the European Union, and we review structure-based computational models that are available and potentially useful in the assessment of acute systemic toxicity. Emphasis is placed on quantitative structure-activity relationship (QSAR) models implemented by means of a range of software tools. The most recently published literature models for acute systemic toxicity are also discussed, and perspectives for future developments in this field are offered.