Bioprotective Role of Phytocompounds Against the Pathogenesis of Non-alcoholic Fatty Liver Disease to Non-alcoholic Steatohepatitis: Unravelling Underlying Molecular Mechanisms.

Non-alcoholic fatty liver disease (NAFLD), with a global prevalence of 25%, continues to escalate, creating noteworthy concerns towards the global health burden. NAFLD causes triglycerides and free fatty acids to build up in the liver. The excessive fat build-up causes inflammation and damages the healthy hepatocytes, leading to non-alcoholic steatohepatitis (NASH). Dietary habits, obesity, insulin resistance, type 2 diabetes, and dyslipidemia influence NAFLD progression. The disease burden is complicated due to the paucity of therapeutic interventions. Obeticholic acid is the only approved therapeutic agent for NAFLD. With more scientific enterprise being directed towards the understanding of the underlying mechanisms of NAFLD, novel targets like lipid synthase, farnesoid X receptor signalling, peroxisome proliferator-activated receptors associated with inflammatory signalling, and hepatocellular injury have played a crucial role in the progression of NAFLD to NASH. Phytocompounds have shown promising results in modulating hepatic lipid metabolism and de novo lipogenesis, suggesting their possible role in managing NAFLD. This review discusses the ameliorative role of different classes of phytochemicals with molecular mechanisms in different cell lines and established animal models. These compounds may lead to the development of novel therapeutic strategies for NAFLD progression to NASH. This review also deliberates on phytomolecules undergoing clinical trials for effective management of NAFLD.

Shaping the Future of Obesity Treatment: In Silico Multi-Modeling of IP6K1 Inhibitors for Obesity and Metabolic Dysfunction.

Recent research has uncovered a promising approach to addressing the growing global health concern of obesity and related disorders. The inhibition of inositol hexakisphosphate kinase 1 (IP6K1) has emerged as a potential therapeutic strategy. This study employs multiple ligand-based in silico modeling techniques to investigate the structural requirements for benzisoxazole derivatives as IP6K1 inhibitors. Firstly, we developed linear 2D Quantitative Structure-Activity Relationship (2D-QSAR) models to ensure both their mechanistic interpretability and predictive accuracy. Then, ligand-based pharmacophore modeling was performed to identify the essential features responsible for the compounds' high activity. To gain insights into the 3D requirements for enhanced potency against the IP6K1 enzyme, we employed multiple alignment techniques to set up 3D-QSAR models. Given the absence of an available X-ray crystal structure for IP6K1, a reliable homology model for the enzyme was developed and structurally validated in order to perform structure-based analyses on the selected dataset compounds. Finally, molecular dynamic simulations, using the docked poses of these compounds, provided further insights. Our findings consistently supported the mechanistic interpretations derived from both ligand-based and structure-based analyses. This study offers valuable guidance on the design of novel IP6K1 inhibitors. Importantly, our work exclusively relies on non-commercial software packages, ensuring accessibility for reproducing the reported models.

Structural insights into the interactions of repositioning and known drugs for Alzheimer's disease with hen egg white lysozyme by MM-GBSA.

Six drugs (dapsone, diltiazem, timolol, rosiglitazone, mesalazine, and milnacipran) that were predicted by network-based polypharmacology approaches as potential anti-Alzheimer's drugs, have been subjected in this study for in silico and in vitro evaluation to check their potential against protein fibrillation, which is a causative factor for multiple diseases such as Alzheimer's disease, Parkinson's disease, Huntington disease, cardiac myopathy, type-II diabetes mellitus and many others. Molecular docking and thereafter molecular dynamics (MD) simulations revealed that diltiazem, rosiglitazone, and milnacipran interact with the binding residues such as Asp52, Glu35, Trp62, and Asp101, which lie within the fibrillating region of HEWL. The MM-GBSA analysis revealed -7.86, -5.05, and -10.29 kcal/mol as the binding energy of diltiazem, rosiglitazone, and milnacipran. The RMSD and RMSF calculations revealed significant stabilities of these ligands within the binding pocket of HEWL. While compared with two reported ligands inhibiting HEWL fibrillation, milnacipran depicted almost similar binding potential with one of the known ligands (Ligand binding affinity -10.66 kcal/mol) of HEWL. Furthermore, secondary structure analyses revealed notable inhibition of the secondary structural changes with our candidate ligand; especially regarding retention of the 3/10 α-helix both by DSSP simulation, Circular dichroism, and FESEM-based microscopic image analyses. Taking further into experimental validation, all three ligands inhibited fibrillation in HEWL in simulated conditions as revealed by blue shift in Congo red assay and later expressing percentage inhibition in ThioflavinT assay as well. However, dose-dependent kinetics revealed that the antifibrillatory effects of drugs are more pronounced at low protein concentrations.Communicated by Ramaswamy H. Sarma.

Finding structural requirements of structurally diverse α-glucosidase and α-amylase inhibitors through validated and predictive 2D-QSAR and 3D-QSAR analyses.

Diabetes mellitus (DM) is a chronic metabolic disorder characterized by hyperglycemic state. The α-glucosidase and α-amylase are considered two major targets for the management of Type 2 DM due to their ability of metabolizing carbohydrates into simpler sugars. In the current study, cheminformatics analyses were performed to develop validated and predictive models with a dataset of 187 α-glucosidase and α-amylase dual inhibitors. Separate linear, interpretable and statistically robust 2D-QSAR models were constructed with datasets containing the activities of α-glucosidase and α-amylase inhibitors with an aim to explain the crucial structural and physicochemical attributes responsible for higher activity towards these targets. Consequently, some descriptors of the models pointed out the importance of specific structural moieties responsible for the higher activities for these targets and on the other hand, properties such as ionization potential and mass of the compounds as well as number of hydrogen bond donors in molecules were found to be crucial in determining the binding potentials of the dataset compounds. Statistically significant 3D-QSAR models were developed with both α-glucosidase and α-amylase inhibition datapoints to estimate the importance of 3D electrostatic and steric fields for improved potentials towards these two targets. Molecular docking performed with selected compounds with homology model of α-glucosidase and X-ray crystal structure of α-amylase largely supported the interpretations obtained from the cheminformatic analyses. The current investigation should serve as important guidelines for the design of future α-glucosidase and α-amylase inhibitors. Besides, the current investigation is entirely performed by using non-commercial open-access tools to ensure easy accessibility and reproducibility of the investigation which may help researchers throughout the world to work more on drug design and discovery.

In Silico Modeling and Structural Analysis of Soluble Epoxide Hydrolase Inhibitors for Enhanced Therapeutic Design.

Human soluble epoxide hydrolase (sEH), a dual-functioning homodimeric enzyme with hydrolase and phosphatase activities, is known for its pivotal role in the hydrolysis of epoxyeicosatrienoic acids. Inhibitors targeting sEH have shown promising potential in the treatment of various life-threatening diseases. In this study, we employed a range of in silico modeling approaches to investigate a diverse dataset of structurally distinct sEH inhibitors. Our primary aim was to develop predictive and validated models while gaining insights into the structural requirements necessary for achieving higher inhibitory potential. To accomplish this, we initially calculated molecular descriptors using nine different descriptor-calculating tools, coupled with stochastic and non-stochastic feature selection strategies, to identify the most statistically significant linear 2D-QSAR model. The resulting model highlighted the critical roles played by topological characteristics, 2D pharmacophore features, and specific physicochemical properties in enhancing inhibitory potential. In addition to conventional 2D-QSAR modeling, we implemented the Transformer-CNN methodology to develop QSAR models, enabling us to obtain structural interpretations based on the Layer-wise Relevance Propagation (LRP) algorithm. Moreover, a comprehensive 3D-QSAR analysis provided additional insights into the structural requirements of these compounds as potent sEH inhibitors. To validate the findings from the QSAR modeling studies, we performed molecular dynamics (MD) simulations using selected compounds from the dataset. The simulation results offered crucial insights into receptor-ligand interactions, supporting the predictions obtained from the QSAR models. Collectively, our work serves as an essential guideline for the rational design of novel sEH inhibitors with enhanced therapeutic potential. Importantly, all the in silico studies were performed using open-access tools to ensure reproducibility and accessibility.

Designing multi-target drugs for the treatment of major depressive disorder.

INTRODUCTION: Major depressive disorders (MDD) pose major health burdens globally. Currently available medications have their limitations due to serious adverse effects, long latency periods as well as resistance. Considering the highly complicated pathological nature of this disorder, it has been suggested that multitarget drugs or multi-target-directed ligands (MTDLs) may provide long-term therapeutic solutions for the treatment of MDD. AREAS COVERED: In the current review, recent lead design and lead modification strategies have been covered. Important investigations reported in the last ten years (2013-2022) for the preclinical development of MTDLs (through synthetic medicinal chemistry and biological evaluation) for the treatment of MDD were discussed as case studies to focus on the recent design strategies. The discussions are categorized on the basis of pharmacological targets. Based on these important case studies, the challenges involved in different design strategies were discussed in detail. EXPERT OPINION: Even though large variations were observed in the selection of pharmacological targets, some potential biological targets (NMDA, melatonin receptors) are required to be explored extensively for the design of MTDLs. Similarly, apart from structure activity relationship (SAR), in silico techniques such as multitasking cheminformatic modeling, molecular dynamics simulation and virtual screening should be exploited to a greater extent.

Predicting the ecotoxicity of endocrine disruptive chemicals: Multitasking in silico approaches towards global models.

Manufactured substances known as endocrine disrupting chemicals (EDCs) released in the environment, through the use of cosmetic products or pesticides, can cause severe eco and cytotoxicity that may induce trans-generational as well as long-term deleterious effects on several biological species at relatively low doses, unlike other classical toxins. As the need for effective, affordable and fast EDCs environmental risk assessment has become increasingly pressing, the present work introduces the first moving average-based multitasking quantitative structure-toxicity relationship (MA-mtk QSTR) modeling specifically developed for predicting the ecotoxicity of EDCs against 170 biological species belonging to six groups. Based on 2,301 data-points with high structural and experimental diversity, as well as on the usage of various advanced machine learning methods, the novel most predictive QSTR models display overall accuracies > 87% in both training and prediction sets. However, maximum external predictivity was achieved when a new multitasking consensus modeling approach was applied to these models. Additionally, the developed linear model provided means to investigate the determining factors for eliciting higher ecotoxicity by the EDCs towards different biological species, identifying several factors such as solvation, molecular mass and surface area as well as the number of specific molecular fragments (e.g.: aromatic hydroxy and aliphatic aldehyde). The resource to non-commercial open-access tools to develop the models is a useful step towards library screening to speed up regulatory decision on discovery of safe alternatives to reduce the hazards of EDCs.

Multi-model in silico characterization of 3-benzamidobenzoic acid derivatives as partial agonists of Farnesoid X receptor in the management of NAFLD.

Non-alcoholic fatty liver disease (NAFLD) is a pathological condition which is strongly correlated with fat accumulation in the liver that has become a major health hazard globally. So far, limited treatment options are available for the management of NAFLD and partial agonism of Farnesoid X receptor (FXR) has proven to be one of the most promising strategies for treatment of NAFLD. In present work, a range of validated predictive cheminformatics and molecular modeling studies were performed with a series of 3-benzamidobenzoic acid derivatives in order to recognize their structural requirements for possessing higher potency towards FXR. 2D-QSAR models were able to extract the most significant structural attributes determining the higher activity towards the receptor. Ligand-based pharmacophore model was created with a novel and less-explored open access tool named QPhAR to acquire information regarding important 3D-pharmacophoric features that lead to higher agonistic potential towards the FXR. The alignment of the dataset compounds based on pharmacophore mapping led to 3D-QSAR models that pointed out the most crucial steric and electrostatic influence. Molecular dynamics (MD) simulation performed with the most potent and the least potent derivatives of the current dataset helped us to understand how to link the structural interpretations obtained from 2D-QSAR, 3D-QSAR and pharmacophore models with the involvement of specific amino acid residues in the FXR protein. The current study revealed that hydrogen bond interactions with carboxylate group of the ligands play an important role in the ligand receptor binding but higher stabilization of different helices close to the binding site of FXR (e.g., H5, H6 and H8) through aromatic scaffolds of the ligands should lead to higher activity for these ligands. The present work affords important guidelines towards designing novel FXR partial agonists for new therapeutic options in the management of NAFLD. Moreover, we relied mainly on open-access tools to develop the in-silico models in order to ensure their reproducibility as well as utilization.

In silico characterization of aryl benzoyl hydrazide derivatives as potential inhibitors of RdRp enzyme of H5N1 influenza virus.

RNA-dependent RNA polymerase (RdRp) is a potential therapeutic target for the discovery of novel antiviral agents for the treatment of life-threatening infections caused by newly emerged strains of the influenza virus. Being one of the most conserved enzymes among RNA viruses, RdRp and its inhibitors require further investigations to design novel antiviral agents. In this work, we systematically investigated the structural requirements for antiviral properties of some recently reported aryl benzoyl hydrazide derivatives through a range of in silico tools such as 2D-quantitative structure-activity relationship (2D-QSAR), 3D-QSAR, structure-based pharmacophore modeling, molecular docking and molecular dynamics simulations. The 2D-QSAR models developed in the current work achieved high statistical reliability and simultaneously afforded in-depth mechanistic interpretability towards structural requirements. The structure-based pharmacophore model developed with the docked conformation of one of the most potent compounds with the RdRp protein of H5N1 influenza strain was utilized for developing a 3D-QSAR model with satisfactory statistical quality validating both the docking and the pharmacophore modeling methodologies performed in this work. However, it is the atom-based alignment of the compounds that afforded the most statistically reliable 3D-QSAR model, the results of which provided mechanistic interpretations consistent with the 2D-QSAR results. Additionally, molecular dynamics simulations performed with the apoprotein as well as the docked complex of RdRp revealed the dynamic stability of the ligand at the proposed binding site of the receptor. At the same time, it also supported the mechanistic interpretations drawn from 2D-, 3D-QSAR and pharmacophore modeling. The present study, performed mostly with open-source tools and webservers, returns important guidelines for research aimed at the future design and development of novel anti-viral agents against various RNA viruses like influenza virus, human immunodeficiency virus-1, hepatitis C virus, corona virus, and so forth.

Predicting the Surface Tension of Deep Eutectic Solvents: A Step Forward in the Use of Greener Solvents.

Deep eutectic solvents (DES) are an important class of green solvents that have been developed as an alternative to toxic solvents. However, the large-scale industrial application of DESs requires fine-tuning their physicochemical properties. Among others, surface tension is one of such properties that have to be considered while designing novel DESs. In this work, we present the results of a detailed evaluation of Quantitative Structure-Property Relationships (QSPR) modeling efforts designed to predict the surface tension of DESs, following the Organization for Economic Co-operation and Development (OECD) guidelines. The data set used comprises a large number of structurally diverse binary DESs and the models were built systematically through rigorous validation methods, including 'mixtures-out'- and 'compounds-out'-based data splitting. The most predictive individual QSPR model found is shown to be statistically robust, besides providing valuable information about the structural and physicochemical features responsible for the surface tension of DESs. Furthermore, the intelligent consensus prediction strategy applied to multiple predictive models led to consensus models with similar statistical robustness to the individual QSPR model. The benefits of the present work stand out also from its reproducibility since it relies on fully specified computational procedures and on publicly available tools. Finally, our results not only guide the future design and screening of novel DESs with a desirable surface tension but also lays out strategies for efficiently setting up silico-based models for binary mixtures.

Moving Average-Based Multitasking In Silico Classification Modeling: Where Do We Stand and What Is Next?

Conventional in silico modeling is often viewed as 'one-target' or 'single-task' computer-aided modeling since it mainly relies on forecasting an endpoint of interest from similar input data. Multitasking or multitarget in silico modeling, in contrast, embraces a set of computational techniques that efficiently integrate multiple types of input data for setting up unique in silico models able to predict the outcome(s) relating to various experimental and/or theoretical conditions. The latter, specifically, based upon the Box-Jenkins moving average approach, has been applied in the last decade to several research fields including drug and materials design, environmental sciences, and nanotechnology. The present review discusses the current status of multitasking computer-aided modeling efforts, meanwhile describing both the existing challenges and future opportunities of its underlying techniques. Some important applications are also discussed to exemplify the ability of multitasking modeling in deriving holistic and reliable in silico classification-based models as well as in designing new chemical entities, either through fragment-based design or virtual screening. Focus will also be given to some software recently developed to automate and accelerate such types of modeling. Overall, this review may serve as a guideline for researchers to grasp the scope of multitasking computer-aided modeling as a promising in silico tool.

Synthesis, anticancer activity, SAR and binding mode of interaction studies of substituted pentanoic acids: part II.

Aim: Our previous results suggest that phenyl/naphthylacetyl pentanoic acid derivatives may exhibit dual MMP-2 and HDAC8 inhibitory activities and show effective cytotoxic properties. Methodology: Here, 13 new compounds (C1-C13) were synthesized and characterized. Along with these new compounds, 16 previously reported phenyl/napthylacetyl pentanoic acid derivatives (C14-C29) were biologically evaluated. Results: Compounds C6 and C27 showed good cytotoxicity against leukemia cell line Jurkat E6.1. The mechanisms of cytotoxicity of these compounds were confirmed by DNA deformation assay and reactive oxygen species assay. MMP-2 and HDAC8 expression assays suggested the dual inhibiting property of these two compounds. These findings were supported by results of molecular docking studies. In silico pharmacokinetic properties showed compounds C6 and C27 have high gastrointestinal absorption. Conclusion: This study highlights the action of phenyl/naphthylacetyl pentanoic acid derivatives as anticancer agents.

First multi-target QSAR model for predicting the cytotoxicity of acrylic acid-based dental monomers.

OBJECTIVE: Acrylic acid derivatives are frequently used as dental monomers and their cytotoxicity towards various cell lines is well documented. This study aims to probe the structural and physicochemical attributes responsible for higher toxicity of dental monomers, using quantitative structure-activity relationships (QSAR) modeling approaches. METHODS: A regression-based linear single-target QSAR (st-QSAR) model was developed with a comparatively small dataset containing 39 compounds, the cytotoxicity of which has been assessed over the Hela S3 cell line. By contrast, a classification-based multi-target QSAR model was developed with 138 compounds, the cytotoxicity of which has been reported against 18 different cell lines. Both models were set up following rigorous validation protocols confirming their statistical significance and robustness. RESULTS: The performance of the linear mt-QSAR model, developed with various feature selection and post-selection similarity searching-based schemes, superseded that of all non-linear models produced with six machine learning methods by hyperparameter optimization. The final derived st-QSAR and mt-QSAR linear models are shown to be highly predictive, as well as revealing the crucial structural and physicochemical factors responsible for higher cytotoxicity of the dental monomers. SIGNIFICANCE: This study is the first attempt on unveiling the cytotoxicity of dental monomers over several cell lines by means of a single multi-target QSAR model. Further, such a model is ready to get widespread applicability in the screening of new monomers, judging from its almost accurate predictions over diverse experimental assay conditions.

Multi-Target In Silico Prediction of Inhibitors for Mitogen-Activated Protein Kinase-Interacting Kinases.

The inhibitors of two isoforms of mitogen-activated protein kinase-interacting kinases (i.e., MNK-1 and MNK-2) are implicated in the treatment of a number of diseases including cancer. This work reports, for the first time, a multi-target (or multi-tasking) in silico modeling approach (mt-QSAR) for probing the inhibitory potential of these isoforms against MNKs. Linear and non-linear mt-QSAR classification models were set up from a large dataset of 1892 chemicals tested under a variety of assay conditions, based on the Box-Jenkins moving average approach, along with a range of feature selection algorithms and machine learning tools, out of which the most predictive one (>90% overall accuracy) was used for mechanistic interpretation of the likely inhibition of MNK-1 and MNK-2. Considering that the latter model is suitable for virtual screening of chemical libraries-i.e., commercial, non-commercial and in-house sets, it was made publicly accessible as a ready-to-use FLASK-based application. Additionally, this work employed a focused kinase library for virtual screening using an mt-QSAR model. The virtual hits identified in this process were further filtered by using a similarity search, in silico prediction of drug-likeness, and ADME profiles as well as synthetic accessibility tools. Finally, molecular dynamic simulations were carried out to identify and select the most promising virtual hits. The information gathered from this work can supply important guidelines for the discovery of novel MNK-1/2 inhibitors as potential therapeutic agents.

Density of Deep Eutectic Solvents: The Path Forward Cheminformatics-Driven Reliable Predictions for Mixtures.

Deep eutectic solvents (DES) are often regarded as greener sustainable alternative solvents and are currently employed in many industrial applications on a large scale. Bearing in mind the industrial importance of DES-and because the vast majority of DES has yet to be synthesized-the development of cheminformatic models and tools efficiently profiling their density becomes essential. In this work, after rigorous validation, quantitative structure-property relationship (QSPR) models were proposed for use in estimating the density of a wide variety of DES. These models were based on a modelling dataset previously employed for constructing thermodynamic models for the same endpoint. The best QSPR models were robust and sound, performing well on an external validation set (set up with recently reported experimental density data of DES). Furthermore, the results revealed structural features that could play crucial roles in ruling DES density. Then, intelligent consensus prediction was employed to develop a consensus model with improved predictive accuracy. All models were derived using publicly available tools to facilitate easy reproducibility of the proposed methodology. Future work may involve setting up reliable, interpretable cheminformatic models for other thermodynamic properties of DES and guiding the design of these solvents for applications.

AKT Inhibitors: The Road Ahead to Computational Modeling-Guided Discovery.

AKT, is a serine/threonine protein kinase comprising three isoforms-namely: AKT1, AKT2 and AKT3, whose inhibitors have been recognized as promising therapeutic targets for various human disorders, especially cancer. In this work, we report a systematic evaluation of multi-target Quantitative Structure-Activity Relationship (mt-QSAR) models to probe AKT' inhibitory activity, based on different feature selection algorithms and machine learning tools. The best predictive linear and non-linear mt-QSAR models were found by the genetic algorithm-based linear discriminant analysis (GA-LDA) and gradient boosting (Xgboost) techniques, respectively, using a dataset containing 5523 inhibitors of the AKT isoforms assayed under various experimental conditions. The linear model highlighted the key structural attributes responsible for higher inhibitory activity whereas the non-linear model displayed an overall accuracy higher than 90%. Both these predictive models, generated through internal and external validation methods, were then used for screening the Asinex kinase inhibitor library to identify the most potential virtual hits as pan-AKT inhibitors. The virtual hits identified were then filtered by stepwise analyses based on reverse pharmacophore-mapping based prediction. Finally, results of molecular dynamics simulations were used to estimate the theoretical binding affinity of the selected virtual hits towards the three isoforms of enzyme AKT. Our computational findings thus provide important guidelines to facilitate the discovery of novel AKT inhibitors.

QSAR-Co-X: an open source toolkit for multitarget QSAR modelling.

Quantitative structure activity relationships (QSAR) modelling is a well-known computational tool, often used in a wide variety of applications. Yet one of the major drawbacks of conventional QSAR modelling is that models are set up based on a limited number of experimental and/or theoretical conditions. To overcome this, the so-called multitasking or multitarget QSAR (mt-QSAR) approaches have emerged as new computational tools able to integrate diverse chemical and biological data into a single model equation, thus extending and improving the reliability of this type of modelling. We have developed QSAR-Co-X, an open source python-based toolkit (available to download at https://github.com/ncordeirfcup/QSAR-Co-X ) for supporting mt-QSAR modelling following the Box-Jenkins moving average approach. The new toolkit embodies several functionalities for dataset selection and curation plus computation of descriptors, for setting up linear and non-linear models, as well as for a comprehensive results analysis. The workflow within this toolkit is guided by a cohort of multiple statistical parameters and graphical outputs onwards assessing both the predictivity and the robustness of the derived mt-QSAR models. To monitor and demonstrate the functionalities of the designed toolkit, four case-studies pertaining to previously reported datasets are examined here. We believe that this new toolkit, along with our previously launched QSAR-Co code, will significantly contribute to make mt-QSAR modelling widely and routinely applicable.

Advanced in Silico Methods for the Development of Anti- Leishmaniasis and Anti-Trypanosomiasis Agents.

Leishmaniasis and trypanosomiasis occur primarily in undeveloped countries and account for millions of deaths and disability-adjusted life years. Limited therapeutic options, high toxicity of chemotherapeutic drugs and the emergence of drug resistance associated with these diseases demand urgent development of novel therapeutic agents for the treatment of these dreadful diseases. In the last decades, different in silico methods have been successfully implemented for supporting the lengthy and expensive drug discovery process. In the current review, we discuss recent advances pertaining to in silico analyses towards lead identification, lead modification and target identification of antileishmaniasis and anti-trypanosomiasis agents. We describe recent applications of some important in silico approaches, such as 2D-QSAR, 3D-QSAR, pharmacophore mapping, molecular docking, and so forth, with the aim of understanding the utility of these techniques for the design of novel therapeutic anti-parasitic agents. This review focuses on: (a) advanced computational drug design options; (b) diverse methodologies - e.g.: use of machine learning tools, software solutions, and web-platforms; (c) recent applications and advances in the last five years; (d) experimental validations of in silico predictions; (e) virtual screening tools; and (f) rationale or justification for the selection of these in silico methods.

A unified in silico model based on perturbation theory for assessing the genotoxicity of metal oxide nanoparticles.

Nanomaterials (NMs) are an ever-increasing field of interest, due to their wide range of applications in science and technology. However, despite providing solutions to many societal problems and challenges, NMs are associated with adverse effects with potential severe damages towards biological species and their ecosystems. Particularly, it has been confirmed that NMs may induce serious genotoxic effects on various biological targets. Given the difficulties of experimental assays for estimating the genotoxicity of many NMs on diverse biological targets, development of alternative methodologies is crucial to establish their level of safety. In silico modelling approaches, such as Quantitative Structure-Toxicity Relationships (QSTR), are now considered a promising solution for such purpose. In this work, a perturbation theory machine learning (PTML) based QSTR approach is proposed for predicting the genotoxicity of metal oxide NMs under various experimental assay conditions. The application of such perturbation approach to 6084 NM-NM pair cases, set up from 78 unique NMs, afforded a final PTML-QSTR model with an accuracy better than 96% for both training and test sets. This model was then used to predict the genotoxicity of some NMs not included in the modelling dataset. The results for this independent data set were in excellent agreement with the experimental ones. Overall, that thus suggests that the derived PTML-QSTR model is a reliable in silico tool to rapidly and cost-efficiently assess the genotoxicity of metal oxide NMs. Finally, and most importantly, the model provides important insights regarding the mechanism of the genotoxicity triggered by these NMs.

Multi-Target Chemometric Modelling, Fragment Analysis and Virtual Screening with ERK Inhibitors as Potential Anticancer Agents.

Two isoforms of extracellular regulated kinase (ERK), namely ERK-1 and ERK-2, are associated with several cellular processes, the aberration of which leads to cancer. The ERK-1/2 inhibitors are thus considered as potential agents for cancer therapy. Multitarget quantitative structure-activity relationship (mt-QSAR) models based on the Box-Jenkins approach were developed with a dataset containing 6400 ERK inhibitors assayed under different experimental conditions. The first mt-QSAR linear model was built with linear discriminant analysis (LDA) and provided information regarding the structural requirements for better activity. This linear model was also utilised for a fragment analysis to estimate the contributions of ring fragments towards ERK inhibition. Then, the random forest (RF) technique was employed to produce highly predictive non-linear mt-QSAR models, which were used for screening the Asinex kinase library and identify the most potential virtual hits. The fragment analysis results justified the selection of the hits retrieved through such virtual screening. The latter were subsequently subjected to molecular docking and molecular dynamics simulations to understand their possible interactions with ERK enzymes. The present work, which utilises in-silico techniques such as multitarget chemometric modelling, fragment analysis, virtual screening, molecular docking and dynamics, may provide important guidelines to facilitate the discovery of novel ERK inhibitors.