Exploring QSTR and toxicophore of hERG K+ channel blockers using GFA and HypoGen techniques

Virtual screening, molecular docking, molecular dynamics, and MM-GBSA approaches identify prospective fructose-1,6-bisphosphatase inhibitors from pineapple for diabetes management

Diabetes affects millions globally and poses treatment challenges. Targeting the enzyme fructose-1,6-bisphosphatase (FBPase) in gluconeogenesis and exploring plant-based therapies offer potential solutions for improving diabetes management while supporting sustainability and medicinal advancements. Utilizing pineapple (Ananas comosus L. Merr.) waste as a source of drug precursors could be valuable for health and environmental care due to its medicinal benefits and abundant yearly biomass production. Therefore, this study conducted a virtual screening to identify potential natural compounds from pineapple that could inhibit FBPase activity. A total of 112 compounds were screened for drug-likeness and ADMET properties, and molecular docking simulations were performed on 20 selected compounds using blind docking. The lead compound, butane-2,3-diyl diacetate, was subjected to 100 ns MD simulations, revealing a binding energy of -5.4 kcal/mol comparable to metformin (-5.6 kcal/mol). The MD simulation also confirmed stable complexes with crucial hydrogen bonds. Glu20, Ala24, Thr27, Gly28, Glu29, Leu30, Val160, Met177, Asp178, and Cys179 were identified as key amino acids that stabilized the human liver FBPase-butane-2,3-diyl diacetate complex, while Tyr215 and Asp218 played a crucial role in the human liver FBPase-Metformin complex. Our study indicates that the lead compound has high intestinal solubility. Therefore, it would show rapid bloodstream distribution and effective action on the target protein, making butane-2,3-diyl diacetate a potential antidiabetic drug candidate. However, further investigations in vitro, preclinical, and clinical trials are required to thoroughly assess its efficacy and safety.Communicated by Ramaswamy H. Sarma.

Applications of the Novel Quantitative Pharmacophore Activity Relationship Method QPhAR in Virtual Screening and Lead-Optimisation

Pharmacophores are an established concept for the modelling of ligand–receptor interactions based on the abstract representations of stereoelectronic molecular features. They became widely popular as filters for the fast virtual screening of large compound libraries. A lot of effort has been put into the development of sophisticated algorithms and strategies to increase the computational efficiency of the screening process. However, hardly any focus has been put on the development of automated procedures that optimise pharmacophores towards higher discriminatory power, which still has to be done manually by a human expert. In the age of machine learning, the researcher has become the decision-maker at the top level, outsourcing analysis tasks and recurrent work to advanced algorithms and automation workflows. Here, we propose an algorithm for the automated selection of features driving pharmacophore model quality using SAR information extracted from validated QPhAR models. By integrating the developed method into an end-to-end workflow, we present a fully automated method that is able to derive best-quality pharmacophores from a given input dataset. Finally, we show how the QPhAR-generated models can be used to guide the researcher with insights regarding (un-)favourable interactions for compounds of interest.

On modeling and utilizing chemical compound information with deep learning technologies: A task-oriented approach

A large number of chemical compounds are available in databases such as PubChem and ZINC. However, currently known compounds, though large, represent only a fraction of possible compounds, which is known as chemical space. Many of these compounds in the databases are annotated with properties and assay data that can be used for drug discovery efforts. For this goal, a number of machine learning algorithms have been developed and recent deep learning technologies can be effectively used to navigate chemical space, especially for unknown chemical compounds, in terms of drug-related tasks. In this article, we survey how deep learning technologies can model and utilize chemical compound information in a task-oriented way by exploiting annotated properties and assay data in the chemical compounds databases. We first compile what kind of tasks are trying to be accomplished by machine learning methods. Then, we survey deep learning technologies to show their modeling power and current applications for accomplishing drug related tasks. Next, we survey deep learning techniques to address the insufficiency issue of annotated data for more effective navigation of chemical space. Chemical compound information alone may not be powerful enough for drug related tasks, thus we survey what kind of information, such as assay and gene expression data, can be used to improve the prediction power of deep learning models. Finally, we conclude this survey with four important newly developed technologies that are yet to be fully incorporated into computational analysis of chemical information.

Discovery of potent apoptosis signal-regulating kinase 1 inhibitors via integrated computational strategy and biological evaluation

Apoptosis signal-regulating Kinase 1 (ASK1) has been confirmed as a potential therapeutic target for the treatment of non-alcoholic steatohepatitis (NASH) disorder and the discovery of ASK1 inhibitors has attracted increasing attention. In this work, a series of in silico methods including pharmacophore screening, docking binding site analysis, protein-ligand interaction fingerprint (PLIF) similarity investigation and molecular docking were applied to find the potential hits from commercial compound databases. Five compounds with potential inhibitory activity were purchased and submitted to biological activity validation. Thus, one hit compound was discovered with micromolar IC₅₀ value (10.59 μM) against ASK1. Results demonstrated that the integration of computation methods and biological test was quite reliable for the discovery of potent ASK1 inhibitors and the strategy could be extended to other similar targets of interest.

Identification of potential leukocyte antigen-related protein (PTP-LAR) inhibitors through 3D QSAR pharmacophore-based virtual screening and molecular dynamics simulation

Owing to its negative regulatory role in insulin signaling, protein tyrosine phosphatase of leukocyte antigen-related protein (PTP-LAR) was widely thought as a potential drug target for diabetes. Now, it was urgent to search for potential LAR inhibitors targeting diabetes. Initially, the pharmacophore models of LAR inhibitors were established with the application of the HypoGen module. The cost analysis, test set validation, as well as Fischer's test was used to verify the efficiency of pharmacophore model. Then, the best pharmacophore model (Hypo-1-LAR) was applied for the virtual screening of the ZINC database. And 30 compounds met the Lipinski's rule of five. Among them, 10 compounds with better binding affinity than the known LAR inhibitor (BDBM50296375) were discovered by docking studies. Finally, molecular dynamics simulations and post-analysis experiments (RMSD, RMSF, PCA, DCCM and RIN) were conducted to explore the effect of ligands (ZINC97018474 and Compound 1) on LAR and preliminary understand why ZINC97018474 had better inhibitory activity than Compound 1 (BDBM50296375). Communicated by Ramaswamy H. Sarma.

Experimentally Validated Pharmacoinformatics Approach to Predict hERG Inhibition Potential of New Chemical Entities

The hERG (human ether-a-go-go-related gene) encoded potassium ion (K⁺) channel plays a major role in cardiac repolarization. Drug-induced blockade of hERG has been a major cause of potentially lethal ventricular tachycardia termed Torsades de Pointes (TdPs). Therefore, we presented a pharmacoinformatics strategy using combined ligand and structure based models for the prediction of hERG inhibition potential (IC₅₀) of new chemical entities (NCEs) during early stages of drug design and development. Integrated GRid-INdependent Descriptor (GRIND) models, and lipophilic efficiency (LipE), ligand efficiency (LE) guided template selection for the structure based pharmacophore models have been used for virtual screening and subsequent hERG activity (pIC₅₀) prediction of identified hits. Finally selected two hits were experimentally evaluated for hERG inhibition potential (pIC₅₀) using whole cell patch clamp assay. Overall, our results demonstrate a difference of less than ±1.6 log unit between experimentally determined and predicted hERG inhibition potential (IC₅₀) of the selected hits. This revealed predictive ability and robustness of our models and could help in correctly rank the potency order (lower μM to higher nM range) against hERG.

3D QSAR Pharmacophore Based Virtual Screening for Identification of Potential Inhibitors for CDC25B

Owing to its fundamental roles in cell cycle phases, the cell division cycle 25B (CDC25B) was broadly considered as potent clinical drug target for cancers. In this study, 3D QSAR pharmacophore models for CDC25B inhibitors were developed by the module of Hypogen. Three methods (cost analysis, test set prediction, and Fisher's test) were applied to validate that the models could be used to predict the biological activities of compounds. Subsequently, 26 compounds satisfied Lipinski's rule of five were obtained by the virtual screening of the Hypo-1-CDC25B against ZINC databases. It was then discovered that 9 identified molecules had better binding affinity than a known CDC25B inhibitors-compound 1 using docking studies. The molecular dynamics simulations showed that the compound had favorable conformations for binding to the CDC25B. Thus, our findings here would be helpful to discover potent lead compounds for the treatment of cancers.

QSAR Models towards Cholinesterase Inhibitors for the Treatment of Alzheimer's Disease

Alzheimer's Disease (AD) is a multifactorial neurological syndrome with the combination of aging, genetic, and environmental factors triggering the pathological decline. Interestingly, the importance of the Acetylcholinesterase (AChE) enzyme has increased due to its involvement in the ß-amyloid peptide fibril formation during AD pathogenesis. In silico technique, QSAR has proven its usefulness in pharmaceutical research for the design/optimization of new chemical entities. Further, QSAR method advanced the scope of rational drug design and the search for the mechanism of drug action. It is a well-established fact that the chemical and pharmaceutical effects of a compound are closely related to its physico-chemical properties, which can be calculated by various methods from the compound structure. This chapter focuses on different Quantitative Structure-Activity Relationship (QSAR) studies carried out for a variety of cholinesterase inhibitors for the treatment of AD. These predictive models will be potentially used for further designing better and safer drugs against AD.

ADMET Evaluation in Drug Discovery. 16. Predicting hERG Blockers by Combining Multiple Pharmacophores and Machine Learning Approaches

Blockade of human ether-à-go-go related gene (hERG) channel by compounds may lead to drug-induced QT prolongation, arrhythmia, and Torsades de Pointes (TdP), and therefore reliable prediction of hERG liability in the early stages of drug design is quite important to reduce the risk of cardiotoxicity-related attritions in the later development stages. In this study, pharmacophore modeling and machine learning approaches were combined to construct classification models to distinguish hERG active from inactive compounds based on a diverse data set. First, an optimal ensemble of pharmacophore hypotheses that had good capability to differentiate hERG active from inactive compounds was identified by the recursive partitioning (RP) approach. Then, the naive Bayesian classification (NBC) and support vector machine (SVM) approaches were employed to construct classification models by integrating multiple important pharmacophore hypotheses. The integrated classification models showed improved predictive capability over any single pharmacophore hypothesis, suggesting that the broad binding polyspecificity of hERG can only be well characterized by multiple pharmacophores. The best SVM model achieved the prediction accuracies of 84.7% for the training set and 82.1% for the external test set. Notably, the accuracies for the hERG blockers and nonblockers in the test set reached 83.6% and 78.2%, respectively. Analysis of significant pharmacophores helps to understand the multimechanisms of action of hERG blockers. We believe that the combination of pharmacophore modeling and SVM is a powerful strategy to develop reliable theoretical models for the prediction of potential hERG liability.

Toxicophore exploration as a screening technology for drug design and discovery: techniques, scope and limitations

Toxicity is a common drawback of newly designed chemotherapeutic agents. With the exception of pharmacophore-induced toxicity (lack of selectivity at higher concentrations of a drug), the toxicity due to chemotherapeutic agents is based on the toxicophore moiety present in the drug. To date, methodologies implemented to determine toxicophores may be broadly classified into biological, bioanalytical and computational approaches. The biological approach involves analysis of bioactivated metabolites, whereas the computational approach involves a QSAR-based method, mapping techniques, an inverse docking technique and a few toxicophore identification/estimation tools. Being one of the major steps in drug discovery process, toxicophore identification has proven to be an essential screening step in drug design and development. The paper is first of its kind, attempting to cover and compare different methodologies employed in predicting and determining toxicophores with an emphasis on their scope and limitations. Such information may prove vital in the appropriate selection of methodology and can be used as screening technology by researchers to discover the toxicophoric potentials of their designed and synthesized moieties. Additionally, it can be utilized in the manipulation of molecules containing toxicophores in such a manner that their toxicities might be eliminated or removed.

Computational investigations of hERG channel blockers: New insights and current predictive models

Identification of potential human Ether-a-go-go Related-Gene (hERG) potassium channel blockers is an essential part of the drug development and drug safety process in pharmaceutical industries or academic drug discovery centers, as they may lead to drug-induced QT prolongation, arrhythmia and Torsade de Pointes. Recent reports also suggest starting to address such issues at the hit selection stage. In order to prioritize molecules during the early drug discovery phase and to reduce the risk of drug attrition due to cardiotoxicity during pre-clinical and clinical stages, computational approaches have been developed to predict the potential hERG blockage of new drug candidates. In this review, we will describe the current in silico methods developed and applied to predict and to understand the mechanism of actions of hERG blockers, including ligand-based and structure-based approaches. We then discuss ongoing research on other ion channels and hERG polymorphism susceptible to be involved in LQTS and how systemic approaches can help in the drug safety decision.

In silico prediction of hERG inhibition

The voltage-gated potassium channel encoded by hERG carries a delayed rectifying potassium current (IKr) underlying repolarization of the cardiac action potential. Pharmacological blockade of the hERG channel results in slowed repolarization and therefore prolongation of action potential duration and an increase in the QT interval as measured on an electrocardiogram. Those are possible to cause sudden death, leading to the withdrawals of many drugs, which is the reason for hERG screening. Computational in silico prediction models provide a rapid, economic way to screen compounds during early drug discovery. In this review, hERG prediction models are classified as 2D and 3D quantitative structure–activity relationship models, pharmacophore models, classification models, and structure based models (using homology models of hERG).

Computational methods in drug discovery

Computer-aided drug discovery/design methods have played a major role in the development of therapeutically important small molecules for over three decades. These methods are broadly classified as either structure-based or ligand-based methods. Structure-based methods are in principle analogous to high-throughput screening in that both target and ligand structure information is imperative. Structure-based approaches include ligand docking, pharmacophore, and ligand design methods. The article discusses theory behind the most important methods and recent successful applications. Ligand-based methods use only ligand information for predicting activity depending on its similarity/dissimilarity to previously known active ligands. We review widely used ligand-based methods such as ligand-based pharmacophores, molecular descriptors, and quantitative structure-activity relationships. In addition, important tools such as target/ligand data bases, homology modeling, ligand fingerprint methods, etc., necessary for successful implementation of various computer-aided drug discovery/design methods in a drug discovery campaign are discussed. Finally, computational methods for toxicity prediction and optimization for favorable physiologic properties are discussed with successful examples from literature.

Chemical structure representations and applications in computational toxicity

Efficient storage and retrieval of chemical structures is one of the most important prerequisite for solving any computational-based problem in life sciences. Several resources including research publications, text books, and articles are available on chemical structure representation. Chemical substances that have same molecular formula but several structural formulae, conformations, and skeleton framework/scaffold/functional groups of the molecule convey various characteristics of the molecule. Today with the aid of sophisticated mathematical models and informatics tools, it is possible to design a molecule of interest with specified characteristics based on their applications in pharmaceuticals, agrochemicals, biotechnology, nanomaterials, petrochemicals, and polymers. This chapter discusses both traditional and current state of art representation of chemical structures and their applications in chemical information management, bioactivity- and toxicity-based predictive studies.

Molecular Docking Guided Comparative GFA, G/PLS, SVM and ANN Models of Structurally Diverse Dual Binding Site Acetylcholinesterase Inhibitors

Recently discovered 42 AChE inhibitors binding at the catalytic and peripheral anionic site were identified on the basis of molecular docking approach, and its comparative quantitative structure-activity relationship (QSAR) models were developed. These structurally diverse inhibitors were obtained by our previously reported high-throughput in vitro screening technique using 384-well plate's assay based on colorimetric method of Ellman. QSAR models were developed using (i) genetic function algorithm, (ii) genetic partial least squares, (iii) support vector machine and (iv) artificial neural network techniques. The QSAR model robustness and significance was critically assessed using different cross-validation techniques on test data set. The generated QSAR models using thermodynamic, electrotopological and electronic descriptors showed that nonlinear methods are more robust than linear methods, and provide insight into the structural features of compounds that are important for AChE inhibition.

Development, validation, and use of quantitative structure-activity relationship models of 5-hydroxytryptamine (2B) receptor ligands to identify novel receptor binders and putative valvulopathic compounds among common drugs

Some antipsychotic drugs are known to cause valvular heart disease by activating serotonin 5-HT(2B) receptors. We have developed and validated binary classification QSAR models capable of predicting potential 5-HT(2B) actives. The classification accuracies of the models built to discriminate 5-HT(2B) actives from the inactives were as high as 80% for the external test set. These models were used to screen in silico 59,000 compounds included in the World Drug Index, and 122 compounds were predicted as actives with high confidence. Ten of them were tested in radioligand binding assays and nine were found active, suggesting a success rate of 90%. All validated actives were then tested in functional assays, and one compound was identified as a true 5-HT(2B) agonist. We suggest that the QSAR models developed in this study could be used as reliable predictors to flag drug candidates that are likely to cause valvulopathy.

hERG-related drug toxicity and models for predicting hERG liability and QT prolongation

BACKGROUND

hERG K(+) channels have been recognized as a primary antitarget in safety pharmacology. Their blockade, caused by several drugs with different therapeutic indications, may lead to QT prolongation and, eventually, to potentially fatal arrhythmia, namely torsade de pointes. Therefore, a number of preclinical models have been developed to predict hERG liability early in the drug development process.

OBJECTIVE

The aim of this review is to outline the present state of the art on drug-induced hERG blockade, providing insights on the predictive value of in vitro and in silico models for hERG liability.

METHODS

On the basis of latest reports, high-throughput preclinical models have been discussed outlining advantages and limitations.

CONCLUSION

Although no single model has an absolute value, an integrated risk assessment is recommended to predict the pro-arrhythmic risk of a given drug. This prediction requires expertise from different areas and should encompass emerging issues such as interference with hERG trafficking and QT shortening.

Ligand-based molecular design of 4-benzylpiperidinealkylureas and amides as CCR3 antagonists

Asthma is an inflammatory disease of the lungs. Clinical studies suggest that eotaxin and chemokine receptor-3 (CCR3) play a primary role in the recruitment of eosinophils in allergic asthma. Development of novel and potent CCR3 antagonists could provide a novel mechanism for inhibition of this recruitment process, thereby preventing asthma. With the intention of designing new ligands with enhanced inhibitor potencies against CCR3, a 3D-QSAR CoMFA study was carried out on 41 4-benzylpiperidinealkylureas and amide derivatives. The best statistics of the developed CoMFA model were r (2) = 0.960, r(2)cv, n = 32 for the training set and r(2)pred, n = 9 for the test set. The generated 3D-QSAR contribution maps shed some light on the effects of the substitution pattern related to CCR3 antagonist activity.

Comparative molecular field analysis and comparative molecular similarity indices analysis of hydroxyethylamine derivatives as selective human BACE-1 inhibitor

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed based on comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), on a series of 43 hydroxyethylamine derivatives, acting as potent inhibitors of beta-site amyloid precursor protein (APP) cleavage enzyme (BACE-1). The crystal structure of the BACE-1 enzyme (PDB ID: 2HM1) with one of the most active compound 28 was available, and we assumed it to be the bioactive conformation of the studied series, for 3D-QSAR analysis. Statistically significant 3D-QSAR model was established on a training set of 34 compounds, which were validated by a test set of 9 compounds. For the best CoMFA model, the statistics are, r2 = 0.998, r2 cv =0.810, n = 34 for the training set and r2 pred = 0.934, n = 9 for the test set. For the best CoMSIA model (combined steric, electrostatic, hydrophobic, and hydrogen bond donor fields), the statistics are r2 = 0.978, r2 cv = 0.754, n = 34 for the training set and r2 pred = 0.750, n = 9 for the test set. The resulting contour maps, produced by the best CoMFA and CoMSIA models, were used to identify the structural features relevant to the biological activity in this series of analogs. The data generated from the present study will further help to design novel, potent, and selective BACE-1 inhibitors.