Furanone derivatives as new inhibitors of CDC7 kinase: development of structure activity relationship model using 3D QSAR, molecular docking, and in silico ADMET

Leveraging new approach methodologies: ecotoxicological modelling of endocrine disrupting chemicals to Danio rerio through machine learning and toxicity studies

New approach methodologies (NAMs) offer information tailored to the intended application while reducing the use of animals. NAMs aim to develop quantitative structure-activity relationship (QSAR) and quantitive-Read-Across structure-activity relationship (q-RASAR) models to predict and categorize the acute toxicity of known and unknown endocrine-disrupting chemicals (EDCs) against zebrafish. EDCs are a diverse group of toxic substances that disrupt the endocrine system of humans and animals. The q-RASAR model was constructed and verified using validation metrics (R2 = 0.886 and Q2 = 0.814) which found to be more reliable model compare to QSAR model. The substructure fingerprint was well-fitted for the classification model and it was validated using 10-fold average accuracy (Q = 86.88%), specificity (Sp = 88.89%), Matthew's correlation curve (MCC = 0.621) and receiver operating characteristics (ROC = 0.828). The dataset of unknown substances revealed that phenolphthalein (Php) exhibited a significant level of toxicity based on q-RASAR model. The docking and simulation study indicated that the computationally derived important features successfully bound to the target zebrafish sex hormone binding globulin (zfSHBG). The experimental LC50 value of 0.790 mg L-1 was very close to the predicted value of 0.763 mg L-1, which provides high confidence to the developed model.

Structure-based virtual screening, ADMET analysis, and molecular dynamics simulation of Moroccan natural compounds as candidates for the SARS-CoV-2 inhibitors

The lack of treatments and vaccines effective against SARS-CoV-2 has forced us to explore natural compounds that could potentially inhibit this virus. Additionally, Morocco is renowned for its rich plant diversity and traditional medicinal uses, which inspires us to leverage our cultural heritage and the abundance of natural resources in our country for therapeutic purposes. In this study, an extensive investigation was conducted to gather a collection of phytoconstituents extracted from Moroccan plants, aiming to evaluate their ability to inhibit the proliferation of the SARS-CoV-2 virus. Molecular docking of the studied compounds was performed at the active sites of the main protease (6lu7) and spike (6m0j) proteins to assess their binding affinity to these target proteins. Compounds exhibiting high affinity to the proteins underwent further evaluation based on Lipinski's rule and ADME-Tox analysis to gain insights into their oral bioavailability and safety. The results revealed that the two compounds demonstrated strong binding affinity to the target proteins, making them potential candidates for oral antiviral drugs against SARS-CoV-2. The molecular dynamics results from this computational analysis supported the overall stability of the resulting complex.

Py-CoMFA, docking, and molecular dynamics simulations of Leishmania (L.) amazonensis arginase inhibitors

Leishmaniasis is a disease caused by a protozoan of the genus Leishmania, affecting millions of people, mainly in tropical countries, due to poor social conditions and low economic development. First-line chemotherapeutic agents involve highly toxic pentavalent antimonials, while treatment failure is mainly due to the emergence of drug-resistant strains. Leishmania arginase (ARG) enzyme is vital in pathogenicity and contributes to a higher infection rate, thus representing a potential drug target. This study helps in designing ARG inhibitors for the treatment of leishmaniasis. Py-CoMFA (3D-QSAR) models were constructed using 34 inhibitors from different chemical classes against ARG from L. (L.) amazonensis (LaARG). The 3D-QSAR predictions showed an excellent correlation between experimental and calculated pIC50 values. The molecular docking study identified the favorable hydrophobicity contribution of phenyl and cyclohexyl groups as substituents in the enzyme allosteric site. Molecular dynamics simulations of selected protein-ligand complexes were conducted to understand derivatives' interaction modes and affinity in both active and allosteric sites. Two cinnamide compounds, 7g and 7k, were identified, with similar structures to the reference 4h allosteric site inhibitor. These compounds can guide the development of more effective arginase inhibitors as potential antileishmanial drugs.

2,4-disubstituted 6-fluoroquinolines as potent antiplasmodial agents: QSAR, homology modeling, molecular docking and ADMET studies

Objective

This work was designed to study 2,4-disubstituted 6-fluoroquinolines as antiplasmodial agents by using in silico techniques, to aid in the design of novel analogs with high potency against malaria and high inhibition of Plasmodium falciparum translation elongation factor 2 (PfeEF2), a novel drug target.

Methods

Quantitative structure-activity relationships (QSAR) of 2,4-disubstituted 6-fluoroquinolines were studied with the genetic function approximation technique in Material Studio software. The 3D structure of PfeEF2 was modeled in the SWISS-MODEL workspace through homology modeling. A molecular docking study of the modeled PfeEF2 and 2,4-disubstituted 6-fluoroquinolines was conducted with Autodock Vina in Pyrx software. Furthermore, the in silico pharmacokinetic properties of selected compounds were investigated.

Results

A robust, reliable and predictive QSAR model was developed that related the chemical structures of 2,4-disubstituted 6-fluoroquinolines to their antiplasmodium activities. The model had an internal squared correlation coefficient R2 of 0.921, adjusted squared correlation coefficient R2adj of 0.878, leave-one-out cross-validation coefficient Q2cv of 0.801 and predictive squared correlation coefficient R2pred of 0.901. The antiplasmodium activity of 6-fluoroquinolines was found to depend on the n5Ring, GGI9, TDB7u, TDB8u and RDF75i physicochemical properties: n5Ring, TDB8u and RDF75i were positively associated, whereas GGI9 and TDB7u were negatively associated, with the antiplasmodium activity of the compounds. Stable complexes formed between the compounds and modeled PfeEF2, with binding affinity ranging from -8.200 to -10.700 kcal/mol. Compounds 5, 11, 16, 22 and 24 had better binding affinities than quinoline-4-carboxamide (DDD107498), as well as good pharmacokinetic properties, and therefore may be better inhibitors of this novel target.

Conclusion

QSAR and docking studies provided insight into designing novel 2,4-disubstituted 6-fluoroquinolines with high antiplasmodial activity and good structural properties for inhibiting a novel antimalarial drug target.

A Structural and In Silico Investigation of Potential CDC7 Kinase Enzyme Inhibitors

A crucial role in the regulation of DNA replication is played by the highly conserved CDC kinase. The CDC7 kinase could serve as a target for therapeutic intervention in cancer. The primary heterocyclic substance is pyrazole, and its derivatives offer great potential as treatments for cancer cell lines. Here, we synthesized the two pyrazole derivatives: 4-(2-(4-chlorophenyl)hydrazinyl)-5-methyl-2-tosyl-1H-pyrazol-3(2H)-one (PYRA-1) and 4-(2-(2,4-difluorophenyl)hydrazinyl)-5-methyl-2-tosyl-1H-pyrazol-3(2H)-one (PYRA-2). The structural confirmation of both the compounds at the three-dimensional level is characterized using single crystal X-ray diffraction and density functional theory. Furthermore, the in silico chemical biological properties were derived using molecular docking and molecular dynamics (MD) simulations. PYRA-1 and PYRA-2 crystallize in the P-1 (a = 8.184(9), b = 14.251(13), c = 15.601(15), α = 91.57(8), β = 97.48(9), 92.67(9), V = 1801.1(3) 3, and Z = 2) and P21/n (a = 14.8648(8), b = 8.5998(4), c = 15.5586(8), β = 116.47(7), V = 1780.4(19) 3, and Z = 4), space groups, respectively. In both PYRA-1 and PYRA-2 compounds, C-H···O intermolecular connections are common to stabilize the crystal structure. In addition, short intermolecular interactions stabilizes with C-H···π and π-π stacking. Crystal packing analysis was quantified using Hirshfeld surface analysis resulting in C···H, O···H, and H···H contacts in PYRA-1 exhibiting more contribution than in PYRA-2. The conformational stabilities of the molecules are same in the gas and liquid phases (water and DMSO). The docking scores measured for PYRA-1 and PYRA-2 with CDC7 kinase complexes are -5.421 and -5.884 kcal/mol, respectively. The MD simulations show that PYRA-2 is a more potential inhibitor than PYRA-1 against CDC7 kinase.

Combined computational approaches for developing new anti-Alzheimer drug candidates: 3D-QSAR, molecular docking and molecular dynamics studies of liquiritigenin derivatives

Butyrylcholinesterase is an acetylcholine-degrading enzyme involved in the memorization process, which is becoming an interesting target for the symptomatic treatment of Alzheimer's disease. In the present investigation, the structure-activity relationship of a set of Liquiritigenin derivatives recently revealed to be Butyrylcholinesterase inhibitors was studied basing on comparative field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMISA). As a result, performant models with high predictive capability have been developed (CoMFA model: R² = 0.91, Q² = 0.62, R² _pred = 0.85; CoMISA model: R² = 0.92, Q² = 0.59, R² _pred = 0.83) and implemented to design new Liquiritigenin derivatives with improved activity. Besides, the affinity of the designed derivatives towards the active site of Butyrylcholinesterase, was confirmed by molecular docking and molecular dynamics studies. Moreover, they exhibited good pharmacokinetics properties. Accordingly, the outcomes of the present investigations can provide important direction for the development of new anti-Alzheimer's drug candidates.

Computer-Aided Drug Design Applied to Secondary Metabolites as Anticancer Agents

Computer-Aided Drug Design (CADD) techniques have garnered a great deal of attention in academia and industry because of their great versatility, low costs, possibilities of cost reduction in in vitro screening and in the development of synthetic steps; these techniques are compared with highthroughput screening, in particular for candidate drugs. The secondary metabolism of plants and other organisms provide substantial amounts of new chemical structures, many of which have numerous biological and pharmacological properties for virtually every existing disease, including cancer. In oncology, compounds such as vimblastine, vincristine, taxol, podophyllotoxin, captothecin and cytarabine are examples of how important natural products enhance the cancer-fighting therapeutic arsenal. In this context, this review presents an update of Ligand-Based Drug Design and Structure-Based Drug Design techniques applied to flavonoids, alkaloids and coumarins in the search of new compounds or fragments that can be used in oncology. A systematical search using various databases was performed. The search was limited to articles published in the last 10 years. The great diversity of chemical structures (coumarin, flavonoids and alkaloids) with cancer properties, associated with infinite synthetic possibilities for obtaining analogous compounds, creates a huge chemical environment with potential to be explored, and creates a major difficulty, for screening studies to select compounds with more promising activity for a selected target. CADD techniques appear to be the least expensive and most efficient alternatives to perform virtual screening studies, aiming to selected compounds with better activity profiles and better "drugability".

Piperazine-substituted derivatives of favipiravir for Nipah virus inhibition: What do in silico studies unravel?

Favipiravir is found to show excellent in-vitro inhibition activity against Nipah virus. To explore the structure-property relationship of Favipiravir, in silico designing of a series of piperazine substituted Favipiravir derivatives are attempted and computational screening has been done to evaluate its bimolecular interactions with Nipah virus. The geometrical features of all the molecules have been addressed from Density Functional Theory calculations. Chemical reactivity descriptor analysis was carried out to understand various reactivity parameters. The drug-likeness properties were estimated by a detailed ADMET study. The binding ability and the mode of binding of these derivatives into the Nipah virus are obtained from molecular docking studies. Our calculations show greater binding ability for the designed inhibitors compared to that of the experimentally reported molecule. Overall, the present work proves to offers new insights and guidelines for synthetic chemists to develop new drugs using piperazine substituted Favipiravir in the treatment of Nipah virus.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s42452-020-04051-9.

Whether the Validation of the Predictive Potential of Toxicity Models is a Solved Task?

Different kinds of biological activities are defined by complex biochemical interactions, which are termed as a "mathematical function" not only of the molecular structure but also for some additional circumstances, such as physicochemical conditions, interactions via energy and information effects between a substance and organisms, organs, cells. These circumstances lead to the great complexity of prediction for biochemical endpoints, since all "details" of corresponding phenomena are practically unavailable for the accurate registration and analysis. Researchers have not a possibility to carry out and analyse all possible ways of the biochemical interactions, which define toxicological or therapeutically attractive effects via direct experiment. Consequently, a compromise, i.e. the development of predictive models of the above phenomena, becomes necessary. However, the estimation of the predictive potential of these models remains a task that is solved only partially. This mini-review presents a collection of attempts to be used for the above-mentioned task, two special statistical indices are proposed, which may be a measure of the predictive potential of models. These indices are (i) Index of Ideality of Correlation; and (ii) Correlation Contradiction Index.

Virtual identification of novel peroxisome proliferator-activated receptor (PPAR) α/δ dual antagonist by 3D-QSAR, molecule docking, and molecule dynamics simulation

The therapeutic potential of PPARs antagonists extends beyond diabetes. PPARs antagonists represent a new drug class that holds promise as a broadly applicable therapeutic approach for cancer treatment. Thus, there is a strong need to develop a rational design strategy for creating PPARs antagonists. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR) models of PPARα receptor (CoMFA-1, q ² = 0.636, r ² = 0.953; CoMSIA-1, q ² = 0.779, r ² = 0.999) and PPARδ receptor (CoMFA-2, q ² = 0.624, r ² = 0.906; CoMSIA-2, q ² = 0.627, r ² = 0.959) were successfully constructed using 35 triazolone ring derivatives. Contour map analysis revealed that the electrostatic and hydrophobic fields played vital roles in the bioactivity of dual antagonists. Molecular docking studies suggested that the hydrogen bonding, electrostatic and hydrophobic interactions all influenced the binding of receptor-ligand complex. Based on the information obtained above, we designed a series of compounds. The docking results were mutually validated with 3D-QSAR results. Three-dimensional-QSAR and absorption, distribution, metabolism, excretion and toxicity (ADMET) predictions indicated that 19 newly designed compounds possessed excellent biological activity and physicochemical properties. In summary, this research could provide theoretical guidance for the structural optimization of novel PPARα and δ dual antagonists. Communicated by Ramaswamy H. Sarma.

QSAR study and rustic ligand-based virtual screening in a search for aminooxadiazole derivatives as PIM1 inhibitors

Background

Quantitative structure–activity relationship (QSAR) was carried out to study a series of aminooxadiazoles as PIM1 inhibitors having pk_i ranging from 5.59 to 9.62 (k_i in nM). The present study was performed using Genetic Algorithm method of variable selection (GFA), multiple linear regression analysis (MLR) and non-linear multiple regression analysis (MNLR) to build unambiguous QSAR models of 34 substituted aminooxadiazoles toward PIM1 inhibitory activity based on topological descriptors.

Results

Results showed that the MLR and MNLR predict activity in a satisfactory manner. We concluded that both models provide a high agreement between the predicted and observed values of PIM1 inhibitory activity. Also, they exhibit good stability towards data variations for the validation methods. Furthermore, based on the similarity principle we performed a database screening to identify putative PIM1 candidates inhibitors, and predict their inhibitory activities using the proposed MLR model.

Conclusions

This approach can be easily handled by chemists, to distinguish, which ones among the future designed aminooxadiazoles structures could be lead-like and those that couldn’t be, thus, they can be eliminated in the early stages of drug discovery process.