Structure-based design and activity modeling of novel epidermal growth factor receptor kinase inhibitors; an in silico approach

Ligand and structure-based virtual screening approaches in drug discovery: minireview

The compilation of ligand and structure-based molecular modeling methods has become an important practice in virtual screening applied to drug discovery. This systematic review addresses and ranks various virtual screening strategies to drive the selection of the optimal method for studies that have as their starting point a multi-ligand investigation and investigation based on the protein structure of a therapeutic target. This study shows examples of applications and an evaluation based on the objective and problematic of a series of virtual screening studies present in the ScienceDirect® database. The results showed that the molecular docking technique is widely used in scientific production, indicating that approaches that use protein structure as a starting point are the most promising strategy for drug discovery that relies on virtual screening-based research.

Molecular modelling studies of substituted indole derivatives as novel influenza a virus inhibitors

The emergence of drug-resistant strains motivate researchers to find new innovative anti-IAV candidates with a different mode of action. In this work, molecular modelling strategies, such as 2D-QSAR, 3D-QSAR, molecular docking, molecular dynamics, FMOs, and ADMET were applied to some substituted indoles as IAV inhibitors. The best-developed 2D-QSAR models, MLR (Q2 = 0.7634, R2train = 0.8666) and ANN[4-3-1] (Q2 = 0.8699, R2train = 0.8705) revealed good statistical validation for the inhibitory response predictions. The 3D-QSAR models, CoMFA (Q2 = 0.504, R2train = 0.805) and CoMSIA/SEDHA (Q2 = 0.619, R2train = 0.813) are selected as the best 3D models following the global thresholds. In addition, the contour maps generated from the CoMFA and CoMSIA models illustrate the relationship between the molecular fields and the inhibitory effects of the studied molecules. The results of the studies led to the design of five new molecules (24a-e) with enhanced anti-IAV activities and binding potentials using the most active molecule (24) as the template scaffold. The conformational stability of the best-designed molecules with the NA protein showed hydrophobic and H-bonds with the key residues from the molecular dynamics simulations of 100 ns. Furthermore, the global reactivity indices from the DFT calculations portrayed the relevance of 24c in view of its smaller band gap as also justified by our QSAR and molecular simulation studies.Communicated by Ramaswamy H. Sarma.

Computer-aided discovery of novel SmDHODH inhibitors for schistosomiasis therapy: Ligand-based drug design, molecular docking, molecular dynamic simulations, drug-likeness, and ADMET studies

Schistosomiasis, also known as bilharzia or snail fever, is a tropical parasitic disease resulting from flatworms of the Schistosoma genus. This often overlooked disease has significant impacts in affected regions, causing enduring morbidity, hindering child development, reducing productivity, and creating economic burdens. Praziquantel (PZQ) is currently the only treatment option for schistosomiasis. Given the potential rise of drug resistance and the limited treatment choices available, there is a need to develop more effective inhibitors for this neglected tropical disease (NTD). In view of this, quantitative structure-activity relationship studies (QSAR), molecular docking, molecular dynamics simulations, drug-likeness, and ADMET predictions were applied to 31 inhibitors of Schistosoma mansoni Dihydroorotate dehydrogenase (SmDHODH). The designed QSAR model demonstrated robust statistical parameters including an R2 of 0.911, R2adj of 0.890, Q2cv of 0.686, R2pred of 0.807, and cR2p of 0.825, confirming its robustness. Compound 26, identified as the most active derivative, emerged as a lead candidate for new potential inhibitors through ligand-based drug design. Subsequently, 12 novel compounds (26A-26L) were designed with enhanced inhibition activity and binding affinity. Molecular docking studies revealed strong and stable interactions, including hydrogen bonding and hydrophobic interactions, between the designed compounds and the target receptor. Molecular dynamics simulations over 100 nanoseconds and MM-PBSA free binding energy (ΔGbind) calculations validated the stability of the two best-designed molecules (26A and 26L). Furthermore, drug-likeness and ADMET prediction analyses affirmed the potential of these designed compounds, suggesting their promise as innovative agents for treating schistosomiasis.

A QSAR study for predicting malformation in zebrafish embryo

BACKGROUND

Developmental toxicity tests are extremely expensive, require a large number of animals, and are time-consuming. It is necessary to develop a new approach to simplify the analysis of developmental endpoints. One of these endpoints is malformation, and one group of ongoing methods for simplifying is in silico models. In this study, we aim to develop a quantitative structure-activity relationship (QSAR) model and identify the best algorithm for predicting malformations, as well as the most important and effective physicochemical properties associated with malformation.

METHODS

The dataset was extracted from a reliable database called COMPTOX. Physicochemical properties (descriptors) were calculated using Mordred and RDKit chemoinformatics software. The data were cleaned, preprocessed, and then split into training and testing sets. Machine learning algorithms, such as gradient boosting model (GBM) and logistic regression (LR), as well as deep learning models, including multilayer perceptron (MLP) and neural networks (NNs) trained with train set data and different sets of descriptors. The models were then validated with test set and various statistical parameters, such as Matthew's correlation coefficient (MCC) and balanced accuracy (BAC) score, were used to compare the models.

RESULTS

A set of descriptors containing with 78% AUC was identified as the best set of descriptors. Gradient boosting was determined to be the best algorithm with 78% predictive power.

CONCLUSIONS

The descriptors that were the most effective for developing models directly impact the mechanism of malformation, and GBM is the best model due to its MCC and BAC.

Unveiling potent inhibitors for schistosomiasis through ligand-based drug design, molecular docking, molecular dynamics simulations and pharmacokinetics predictions

Schistosomiasis is a neglected tropical disease which imposes a considerable and enduring impact on affected regions, leading to persistent morbidity, hindering child development, diminishing productivity, and imposing economic burdens. Due to the emergence of drug resistance and limited management options, there is need to develop additional effective inhibitors for schistosomiasis. In view of this, quantitative structure-activity relationship studies, molecular docking, molecular dynamics simulations, drug-likeness and pharmacokinetics predictions were applied to 39 Schistosoma mansoni Thioredoxin Glutathione Reductase (SmTGR) inhibitors. The chosen QSAR model demonstrated robust statistical parameters, including an R2 of 0.798, R2adj of 0.767, Q2cv of 0.681, LOF of 0.930, R2test of 0.776, and cR2p of 0.746, confirming its reliability. The most active derivative (compound 40) was identified as a lead candidate for the development of new potential non-covalent inhibitors through ligand-based design. Subsequently, 12 novel compounds (40a-40l) were designed with enhanced anti-schistosomiasis activity and binding affinity. Molecular docking studies revealed strong and stable interactions, including hydrogen bonding, between the designed compounds and the target receptor. Molecular dynamics simulations over 100 nanoseconds and MM-PBSA free binding energy (ΔGbind) calculations validated the stability of the two best-designed molecules. Furthermore, drug-likeness and pharmacokinetics prediction analyses affirmed the potential of these designed compounds, suggesting their promise as innovative agents for the treatment of schistosomiasis.

In-silico molecular modelling studies of some camphor imine based compounds as anti-influenza A (H1N1) pdm09 virus agents

The advent of influenza A (H1N1) drug-resistant strains led to the search quest for more potent inhibitors of the influenza A virus, especially in this devastating COVID-19 pandemic era. Hence, the present research utilized some molecular modelling strategies to unveil new camphor imine-based compounds as anti-influenza A (H1N1) pdm09 agents. The 2D-QSAR results revealed GFA-MLR (R2train = 0.9158, Q2=0.8475) and GFA-ANN (R2train = 0.9264, Q2=0.9238) models for the anti-influenza A (H1N1) pdm09 activity prediction which have passed the QSAR model acceptability thresholds. The results from the 3D-QSAR studies also revealed CoMFA (R2train =0.977, Q2=0.509) and CoMSIA_S (R2train =0.976, Q2=0.527) models for activity predictions. Based on the notable information derived from the 2D-QSAR, 3D-QSAR, and docking analysis, ten (10) new camphor imine-based compounds (22a-22j) were designed using the most active compound 22 as the template. Furthermore, the high predicted activity and binding scores of compound 22j were further justified by the high reactive sites shown in the electrostatic potential maps and other quantum chemical calculations. The MD simulation of 22j in the active site of the influenza hemagglutinin (HA) receptor confirmed the dynamic stability of the complex. Moreover, the appraisals of drug-likeness and ADMET properties of the proposed compounds showed zero violation of Lipinski's criteria with good pharmacokinetic profiles. Hence, the outcomes in this work recommend further in-depth in vivo and in-vitro investigations to validate these theoretical findings.Communicated by Ramaswamy H. Sarma.

2D-QSAR, molecular docking, drug-likeness, and ADMET/pharmacokinetic predictions of some non-small cell lung cancer therapeutic agents

Objectives

Non-small cell lung cancer (NSCLC) is the most common type of lung cancer, with nearly 2 million diagnoses and a 17% 5-year survival rate. The aim of this study was to use computer-aided techniques to identify potential therapeutic agents for NSCLC.

Methods

The two dimensional-quantitative structure-activity relationship (2D-QSAR) modeling was employed on some potential NSCLC therapeutic agents to develop a highly predictive model. Molecular docking-based virtual screening were conducted on the same set of compounds to identify potential hit compounds. The pharmacokinetic features of the best hits were evaluated using SWISSADME and pkCSM online web servers, respectively.

Results

The model generated via 2D-QSAR modeling was highly predictive with R2= 0.798, R2adj = 0.754, Q2CV = 0.673, R2 test = 0.531, and cRp2 = 0.627 assessment parameters. Molecular docking-based virtual screening identified compounds 25, 32, 15, 21, and 23 with the highest MolDock scores as the best hits, of which compound 25 had the highest MolDock score of -138.329 kcal/mol. All of the identified hits had higher MolDock scores than the standard drug (osimertinib). The best hit compounds were ascertained to be drug-like in nature following the Lipinski's rule of five. Also, their ADMET features displayed average pharmacokinetic profiles.

Conclusion

After successful preclinical testing, the hit compounds identified in this study may serve as potential NSCLC therapeutic agents due to their safety and efficacy with the exception of compound 23, which was found to be toxic. They can also serve as a template for designing novel NSCLC therapeutic agents.

Computational modelling studies of some 1,3-thiazine derivatives as anti-influenza inhibitors targeting H1N1 neuraminidase via 2D-QSAR, 3D-QSAR, molecular docking, and ADMET predictions

Background

Influenza virus disease remains one of the most contagious diseases that aided the deaths of many patients, especially in this COVID-19 pandemic era. Recent discoveries have shown that the high prevalence of influenza and SARS-CoV-2 coinfection can rapidly increase the death rate of patients. Hence, it became necessary to search for more potent inhibitors for influenza disease therapy. The present study utilized some computational modeling concepts such as 2D-QSAR, 3D-QSAR, molecular docking simulation, and ADMET predictions of some 1,3-thiazine derivatives as inhibitors of influenza neuraminidase (NA).

Results

The 2D-QSAR modeling results showed GFA-MLR ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9192, Q2 = 0.8767, R2adj = 0.8991, RMSE = 0.0959, $$R_{{{\text{test}}}}^{2}$$ R test 2  = 0.8943, $$R_{{{\text{pred}}}}^{2}$$ R pred 2  = 0.7745) and GFA-ANN ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9227, Q2 = 0.9212, RMSE = 0.0940, $$R_{{{\text{test}}}}^{2}$$ R test 2  = 0.8831, $$R_{{{\text{pred}}}}^{2}$$ R pred 2  = 0.7763) models with the computed descriptors as ATS7s, SpMax5_Bhv, nHBint6, and TDB9m for predicting the NA inhibitory activities of compounds which have passed the global criteria of accepting QSAR model. The 3D-QSAR modeling was carried out based on the comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA). The CoMFA_ES ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.9620, Q2 = 0.643) and CoMSIA_SED ($$R_{{\text{train }}}^{2}$$ R train 2  = 0.8770, Q2 = 0.702) models were found to also have good and reliable predicting ability. The compounds were also virtually screened based on their binding scores via molecular docking simulations with the active site of the NA (H1N1) target receptor which also confirms their resilient potency. Four potential lead compounds (4, 7, 14, and 15) with the relatively high inhibitory rate (> 50%) and docking (> − 6.3 kcal/mol) scores were identified as the possible lead candidates for in silico exploration of improved anti-influenza agents.

Conclusion

The drug-likeness and ADMET predictions of the lead compounds revealed non-violation of Lipinski’s rule and good pharmacokinetic profiles as important guidelines for rational drug design. Hence, the outcome of this research set a course for the in silico design and exploration of novel NA inhibitors with improved potency.

In-silico modelling studies of 5-benzyl-4-thiazolinone derivatives as influenza neuraminidase inhibitors via 2D-QSAR, 3D-QSAR, molecular docking, and ADMET predictions

Influenza virus disease is one of the most infectious diseases responsible for many human deaths, and the high mutability of the virus causes drug resistance effects in recent times. As such, it became necessary to explore more inhibitors that could avert future influenza pandemics. The present research utilized some in-silico modelling concepts such as 2D-QSAR, 3D-QSAR, molecular docking simulation, and ADMET predictions on some 5-benzyl-4-thiazolinone derivatives as influenza neuraminidase (NA) inhibitors. The 2D-QSAR modelling results revealed GFA-MLR (R train 2 =0.8414, Q2 = 0.7680) and GFA-ANN (R train ​ 2 =0.8754, Q2 = 0.8753) models with the most relevant descriptors (MATS3i, SpMax5_Bhe, minsOH and VE3_D) for predicting the inhibitory activities of the molecules which has passed the global criteria of accepting QSAR models. The results of the 3D-QSAR modelling results showed that CoMFA_ES (R train ​ 2 =0.9030, Q2 = 0.5390) and CoMSIA_EA (R train 2 =0.880, Q2 = 0.547) models are having good predicting ability among other developed models. The molecules were virtually screened via molecular docking simulation with the active site of NA protein receptor (pH1N1) which confirms their resilient potency when compared with zanamivir standard drug. Molecule 11 as the most potent molecule formed more H-bond interactions with the key residues such as TRP178, ARG152, ARG292, ARG371, and TYR406 that triggered the catalytic reactions for NA inhibition. Furthermore, six (6) molecules (9, 10, 11, 17, 22, and 31) with relatively high inhibitory activities and docking scores were identified as the possible leads for in-silico exploration of novel NA inhibitors. The drug-likeness and ADMET predictions of the lead molecules revealed non-violation of Lipinski's rule and good pharmacokinetic profiles respectively, which are important guidelines for rational drug design. Hence, the outcome of this study overlaid a solid foundation for the in-silico design and exploration of novel NA inhibitors with improved potency.

Computer-aided design of some quinazoline analogues as epidermal growth factor receptor inhibitors

Background

The treatment of epidermal growth factor receptor (EGFR)-muted non-small cell lung cancer (NSCLC) remains among the utmost important unachieved therapeutic need worldwide. Development of EGFR inhibitors to treat NSCLC mutations has been among the difficult tasks faced by researchers in this area. As such, there is a need to discover more EGFR inhibitors. The purpose of this work is to perform computer-aided/structure-based design of novel EGFR inhibitors, elucidate their nature of interactions with their target, and also assess their ADMET properties as well as their drug-likeness, respectively. Compound 17 with a highest binding affinity of −9.5kcal/mol was identified as the template hit compound using molecular docking virtual screening in our previous work. The compound interacted with the active site of the EGFR receptor via hydrogen bond with the following amino acid residues MET793, MET793, THR854, and ASP855 with bond distances of 2.61394 (Å), 2.18464 (Å), 2.57601 (Å), and 2.68794 (Å), respectively. It also interacted with the active site of the EGFR receptor via halogen bond (GLN791), hydrophobic bond (LEU718, CYS797, LYS745, ALA743, ALA743, and VAL726), electrostatic bond (LYS745), and others (MET766), respectively. Furthermore, from our previous study, the following descriptors (ATSC6m, ATSC8e, MATS7m, SpMax3_Bhp, SpMax5_Bhs, and MaxHBint10) contained in the reported model were found to be responsible for the inhibitory activities of the studied compounds. In this research, the template (compound 17) was modified manually by attaching halo-phenyl and halo-phenyl-amino rings on the para position of the flouro-nitro-benzamide moiety of the template compound, respectively.

Results

A computer-aided design/structure-based approach was used to design six new EGFR inhibitors using molecule 17 as the template compound for the design identified in our previously reported work. Molecular docking investigation was performed to elucidate the binding mode of these newly designed EGFR inhibitors with the binding pose of EGFR receptor (pdb code 4ZAU) and found to have better affinities which range from −9.5 to −10.4 kcal/mol than the template compound and gefitinib, the control, respectively. The ADMET property assessment of these newly designed EGFR inhibitors indicated that they were orally bioavailable with good absorption, distribution, metabolism, and excretory properties with no toxicity. And for their drug-likeness, they were seen to have a higher molecular weight which might be as a result of halo-phenyl-amino ring attachments. Based on this finding, halo-phenyl-amino rings might be responsible for the inhibitory activities of these newly designed compounds.

Conclusion

The six newly designed EGFR inhibitors were found to have higher binding affinities toward their target EGFR receptor than the template compound and gefitinib which was used as the control in this research. They were seen to have good ADMET and drug-like properties which indicate that they might be orally bioavailable. Furthermore, according to their synthetic accessibility score, they can be easily synthesized in the laboratory because the values were found to be less than five which fall within the easy portion of the scale. Therefore, this research recommends that these newly designed EGFR inhibitors should be synthesized most especially those with higher binding affinities, good ADMET, and drug-likeness properties than the template compound.

Design of more potent quinazoline derivatives as EGFRWT inhibitors for the treatment of NSCLC: a computational approach

Background

Lung cancer remains the leading and deadly type of cancer worldwide. It was estimated to account for about 25% of the 7 million people that died as a result of cancer-related issues/mortality every year in the world. Non-small cell lung cancer (NSCLC) is the lethal/deadly class of lung cancer with nearly 1.5 million reported cases and less than 20% survival rate. Therefore, it becomes necessary to explore more effective NSCLC drugs.

Result

A computational approach was employed here to design ten new EGFRWT inhibitors using compound 18 as a template for the design identified with the best binding affinity and good pharmacokinetic properties previously reported in our work. The modeled inhibitory activities of these newly designed EGFRWT inhibitors (range from 7.746966 to 11.09261) were better than that of the hit compound with pIC50 of 7.5639 and gefitinib the positive control with pIC50 of 5.879426. The ligand-binding interaction between these newly designed EGFRWT inhibitors and the EGFR tyrosine kinase receptor as shown in Table 3 was investigated and elucidated using molecular docking protocol. Based on the molecular docking results, the binding affinities of these newly designed EGFRWT inhibitors were found to be between − 8.8 and − 9.5 kcal/mol. The designed compound SFD10 has the highest binding affinity of − 9.5 kcal/mol followed by compound SFD8 (with a binding affinity of − 9.3 kcal/mol), then by compound SFD9 and 4 (each with a binding affinity of − 9.3 kcal/mol). None of them was found to have more than one violation of the filtering criterion used in this study thereby showing good ADMET properties.

Conclusion

The modeled inhibitory activities and binding affinities of these newly designed EGFRWT inhibitors were found to be higher than that of the template compound and the control (gefitinib) used in this research. They were also seen to be non-toxic with good pharmacokinetic properties.

Biologically active phthalocyanine metal complexes: Preparation, evaluation of α-glycosidase and anticholinesterase enzyme inhibition activities, and molecular docking studies

In this study, preparation, as well as investigation of α-glycosidase and cholinesterase (ChE) enzyme inhibition activities of furan-2-ylmethoxy-substituted compounds 1-7, are reported. Peripherally, tetra-substituted copper and manganese phthalocyanines (5 and 6) were synthesized for the first time. The substitution of furan-2-ylmethoxy groups provides remarkable solubility to the complex and redshift of the phthalocyanines Q-band. Besides, the inhibitory effects of these compounds on acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and α-glycosidase (α-Gly) enzymes have been investigated. The AChE was inhibited by these compounds (1-7) in low micromolar levels, and K _i values were recorded between 11.17 ± 1.03 and 83.28 ± 11.08 µM. Against the BChE, the compounds demonstrated K _i values from 7.55 ± 0.98 to 81.35 ± 12.80 µM. Also, these compounds (1-7) effectively inhibited α-glycosidase, with K _i values in the range of 744.87 ± 67.33 to 1094.38 ± 88.91 µM. For α-glycosidase, the most effective K _i values of phthalocyanines 3 and 6 were with K _i values of 744.87 ± 67.33 and 880.36 ± 56.77 µM, respectively. Moreover, the studied metal complexes were docked with target proteins PDB ID: 4PQE, 1P0I, and 3WY1. Pharmacokinetic parameters and secondary chemical interactions that play an active role in interaction were predicted with docking simulation results. Overall, furan-2-ylmethoxy-substituted phthalocyanines can be considered as potential agents for the treatment of Alzheimer's diseases and diabetes mellitus.