3D-QSAR, molecular docking, DFT and ADMET studies on quinazoline derivatives to explore novel DHFR inhibitors

Amphiphilic amidines as potential plasmic membrane-targeting antifungal agents: synthesis, bio-activities and QSAR

BACKGROUND

Cationic antimicrobial peptides (AMPs) possess broad-spectrum biological activities with less inclination to inducing antibiotic resistance. Herein a battery of amphiphilic amidines were designed by mimicking the characteristics of AMPs. The antifungal activities and the effects to the hyphal morphology and membrane permeability were investigated.

RESULTS

The results indicated the inhibitory rates of ten compounds were over 80% to Botrytis cinerea and ten compounds over 90% to Valsa mali Miyabe et Yamada at 50 mg L-1. The half maximal effective concentration (EC50) values of compound 5g and 6g to V. mali were 1.21 and 1.90 mg L-1 respectively. The protective rate against apple canker of compound 5g reached 93.4% at 100 mg L-1 on twigs, superior to carbendazim (53.3%). When treated with 5g, the cell membrane permeability and leakage of content of V. mali increased, accompanied with the decrease of superoxide dismutase (SOD) and catalase (CAT) level. Concurrently, the mycelial hyphae contracted, wrinkled, and collapsed, providing evidence of membrane perturbation. A three-dimensional quantitative structure-activity relationship (3D-QSAR) between the topic compounds and the EC50 to V. mali was established showing good predictability (r2 = 0.971).

CONCLUSION

Amphiphilic amidines can acquire antifungal activities by acting on the plasmic membrane. Compound 5g could be a promising lead in discovering novel fungicidal candidates. © 2024 Society of Chemical Industry.

Discovery of novel bromodomain-containing protein 4 (BRD4-BD1) inhibitors combined with 3d-QSAR, molecular docking and molecular dynamics in silico

Bromine-containing domain protein 4 (BRD4) plays a crucial role in regulating transcription and genome stability. Selective inhibitors of BRD4-BD1 can specifically target specific bromine domains to affect cell proliferation, apoptosis, and differentiation. In this work, 43 selective benzoazepinone BRD4-BD1 inhibitors were studied using molecular simulations and three-dimensional quantitative conformation relationships (3D-QSAR). A reliable 3D-QSAR model was established based on COMFA (Q2 = 0.532, R2 = 0.981) and COMSIA (S + E + H (Q2 = 0.536, R2 = 0.979) two different analysis methods. Through 3D-QSAR model prediction and quantum chemical analysis, 15 small molecules with stronger inhibitory activity than the template compounds were constructed, and 5 new compounds with higher predictive activity and binding affinity were screened by molecular docking and ADMET methods. According to the molecular dynamics simulation, the key residues that can interact with BRD4-BD1 protein and molecular docking results are consistent, including ASN140, MET132, GLN85, MET105, ASN135 and TYR97. From the MD trajectory, we calculated and analyzed RMSD, RMSF, free binding energy, FECM, DCCM and PCA, the loop region formed by amino acids VAL45∼PRO62 showed α-helix, β-folding and clustering towards the active center with the extension of simulation time. Further optimization of the structure of active candidate compounds A6, A11, A14, and A15 will provide the necessary theoretical basis for the synthesis and activity evaluation of novel BRD4-BD1 inhibitors.Communicated by Ramaswamy H. Sarma.

Design of new molecules against cervical cancer using DFT, theoretical spectroscopy, 2D/3D-QSAR, molecular docking, pharmacophore and ADMET investigations

Cervical cancer is a major health problem of women. Hormone therapy, via aromatase inhibition, has been proposed as a promising way of blocking estrogen production as well as treating the progression of estrogen-dependent cancer. To overcome the challenging complexities of costly drug design, in-silico strategy, integrating Structure-Based Drug Design (SBDD) and Ligand-Based Drug Design (LBDD), was applied to large representative databases of 39 quinazoline and thioquinazolinone compound derivatives. Quantum chemical and physicochemical descriptors have been investigated using density functional theory (DFT) and MM2 force fields, respectively, to develop 2D-QSAR models, while CoMSIA and CoMFA descriptors were used to build 3D-QSAR models. The robustness and predictive power of the reliable models were verified, via several validation methods, leading to the design of 6 new drug-candidates. Afterwards, 2 ligands were carefully selected using virtual screening methods, taking into account the applicability domain, synthetic accessibility, and Lipinski's criteria. Molecular docking and pharmacophore modelling studies were performed to examine potential interactions with aromatase (PDB ID: 3EQM). Finally, the ADMET properties were investigated in order to select potential drug-candidates against cervical cancer for experimental in vitro and in vivo testing.

Regulating pri/pre-microRNA up/down expressed in cancer proliferation, angiogenesis and metastasis using selected potent triterpenoids

MicroRNAs (miRNAs) play a crucial role in cancer progression by selectively inducing translational degradation of messenger RNA (mRNA) via sequence-specific interactions with the 3'-untranslated region (3'-UTR). The potential targeting of miRNA has been recognized as a significant avenue for investigating the biological progression of diverse cancer types. Consequently, targeting of pri-miRNA and pre-miRNA by phytochemicals emerges as a viable strategy in the realm of anticancer therapies. Among phytochemicals, triterpenoids have garnered significant recognition for their chemotherapeutic and chemopreventive capabilities in combating multiple cancers. To date, there is a dearth of literature about the molecular interactions between triterpenoids and miRNAs. The primary objective of this investigation is to discern the potential triterpenoids that can function as modulators for specific miRNAs, namely pri-miRNA-19b-2, pre-miR21, microRNA 20b, pri-miRNA-208a, pri-miRNA-378a, pri-miRNA-320b-2, and pri-miRNA-300, achieved through the use of in silico investigations. The study primarily focused on performing drug-likeness, computer-aided toxicity, and pharmacokinetic prediction studies for triterpenoids. Furthermore, molecular docking and simulation techniques were employed to investigate these compounds. The triterpenoids studied were shown to have drug-likeness characteristics, although asiatic acid, lupeol, and pristimerin were able to pass all toxicity tests. Among the triterpenoids that underwent docking, pristimerin had a significant binding energy of -10.9 kcal/mol during its interaction with pri-miR-378a. The stable interaction between the pristimerin and miRNA complex was demonstrated by molecular dynamics simulation. As a result, pristimerin has the potential to act as a modulator of carcinogenic miRNAs, making it a promising candidate for cancer prevention and treatment due to its tailored modulation of miRNA activity.

Design of new small molecules derived from indolin-2-one as potent TRKs inhibitors using a computer-aided drug design approach

Tropomyosin receptor kinase (TRKs) enzymes are responsible for cancers associated with the neurotrophic tyrosine kinase receptor gene fusion and are identified as effective targets for anticancer drug discovery. A series of small-molecule indolin-2-one derivatives showed remarkable biological activity against TRKs enzymatic activity. These small molecules could have an excellent profile for pharmaceutical application in the treatment of cancers caused by TRKs activity. The aim of this study is to modify the structure of these molecules to obtain new molecules with improved TRK inhibitory activity and pharmacokinetic properties favorable to the design of new drugs. Based on these series, we carried out a 3D-QSAR study. As a result, robust and reliable CoMFA and CoMSIA models are developed and applied to the design of 11 new molecules. These new molecules have a biological activity superior to the most active molecule in the starting series. The eleven designed molecules are screened using drug-likeness, ADMET proprieties, molecular docking, and MM-GBSA filters. The results of this screening identified the T1, T3, and T4 molecules as the best candidates for strong inhibition of TRKs enzymatic activity. In addition, molecular dynamics simulations are performed for TRK free and complexed with ligands T1, T3, and T4 to evaluate the stability of ligand-protein complexes over the simulation time. On the other hand, we proposed experimental synthesis routes for these newly designed molecules. Finally, the designed molecules T1, T2, and T3 have great potential to become reliable candidates for the conception of new drug inhibitors of TRKs.Communicated by Ramaswamy H. Sarma.

Discovery of anti-colon cancer agents targeting wild-type and mutant p53 using computer-aided drug design

Mutations in the p53 gene are common and occur in over 50% of all cancers, as it is involved in DNA damage repair, cell cycle regulation and apoptosis. Moreover, the p53 gene is mutated in 70% of colon cancers. Therefore, the development of drugs to combat this mutation requires urgent attention. With this in mind, in silico drug design approaches were applied on quinoline derivatives with anticancer activity. In 3D-QSAR study, steric, electrostatic, hydrophobic and H-bond acceptor fields (SEHA) play an important role in prediction and design of new colon cancer compounds. Indeed, the two best CoMSIA/SEHA models with (Q2 = 0.737, R2 = 0.914, R pred 2 = 0.720) and (Q2 = 0.738, R2 = 0.919, R pred 2 = 0.739) show good prediction of human colon carcinoma HCT 116 (p53+/+) and (p53-/-) activities, respectively. Furthermore, the predictive ability and robustness of these models were tested by several validation methods. Molecular docking analyses reveal crucial interactions with the active sites of the p53 protein in both wild type and mutant. Based on these theoretical studies, we designed 10 new compounds with good anticancer activity potential, which were evaluated using ADMET properties. Molecular dynamics simulations were performed to confirm the detailed binding mode of the docking results. Finally, the MM-GBSA based on molecular dynamics simulation confirmed that the designed compounds were able to form stable hydrogen bonding interactions with the crucial residues, which are essential to overcome the p53 mutation in colon cancer.Communicated by Ramaswamy H. Sarma.

Computational insights into pediatric adenovirus inhibitors: in silico strategies for drug repurposing

Human adenovirus (HADV) infection can pose a serious threat to children, leading to a variety of respiratory illnesses and other complications. Particularly, children with weak immune systems are vulnerable to severe adenovirus infections with high mortality. The main focus of this study is to propose new antiviral agents as lead HADV inhibitors for children. So, several antiviral agents used in children were subjected to finding new HADV inhibitors using important computational methods of molecular docking, molecular dynamics (MD) simulation, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding free energy calculations, density functional theory (DFT), and pharmacokinetic analysis. Molecular docking of standard cidofovir along with other ligands, suggested that sofosbuvir has the highest binding energy (-10.8 kcal/mol), followed by baloxavir marboxil (-10.36 kcal/mol). Further, the analysis of molecular interactions using MD simulation (100 ns) and MM-PBSA indicated that baloxavir marboxil has formed the most stable protein-ligand complex with HADV, followed by sofosbuvir. The binding free energies of baloxavir marboxil and sofosbuvir were found to be -61.724 kJ/mol and -48.123 kJ/mol, respectively. The DFT and drug-likeness properties of these compounds were also investigated. Overall, two antiviral agents, such as baloxavir marboxil, and sofosbuvir are suggested as lead repurposed candidates against HADV.Communicated by Ramaswamy H. Sarma.

Antioxidant, Volatile Compounds; Antimicrobial, Anti-Inflammatory, and Dermatoprotective Properties of Cedrus atlantica (Endl.) Manetti Ex Carriere Essential Oil: In Vitro and In Silico Investigations

Cedrus atlantica (Endl.) Manetti ex Carriere is an endemic tree possessing valuable health benefits which has been widely used since time immemorial in international traditional pharmacopoeia. The aim of this exploratory investigation is to determine the volatile compounds of C. atlantica essential oils (CAEOs) and to examine their in vitro antimicrobial, antioxidant, anti-inflammatory, and dermatoprotective properties. In silico simulations, including molecular docking and pharmacokinetics absorption, distribution, metabolism, excretion, and toxicity (ADMET), and drug-likeness prediction were used to reveal the processes underlying in vitro biological properties. Gas chromatography-mass spectrophotometry (GC-MS) was used for the chemical screening of CAEO. The antioxidant activity of CAEO was investigated using four in vitro complementary techniques, including ABTS and DPPH radicals scavenging activity, ferric reductive power, and inhibition of lipid peroxidation (β-carotene test). Lipoxygenase (5-LOX) inhibition and tyrosinase inhibitory assays were used for testing the anti-inflammatory and dermatoprotective properties. GC-MS analysis indicated that the main components of CAEO are β-himachalene (28.99%), α-himachalene (14.43%), and longifolene (12.2%). An in vitro antimicrobial activity of CAEO was examined against eleven strains of Gram-positive bacteria (three strains), Gram-negative bacteria (four strains), and fungi (four strains). The results demonstrated high antibacterial and antifungal activity against ten of them (>15 mm zone of inhibition) using the disc-diffusion assay. The microdilution test showed that the lowest values of MIC and MBC were recorded with the Gram-positive bacteria in particular, which ranged from 0.0625 to 0.25 % v/v for MIC and from 0.5 to 0.125 % v/v for MBC. The MIC and MFC of the fungal strains ranged from 0.5 to 4.0% (MIC) and 0.5 to 8.0% v/v (MFC). According to the MBC/MIC and MFC/MIC ratios, CAEO has bactericidal and fungicidal activity. The results of the in vitro antioxidant assays revealed that CAEO possesses remarkable antioxidant activity. The inhibitory effects on 5-LOX and tyrosinase enzymes was also significant (p < 0.05). ADMET investigation suggests that the main compounds of CAEO possess favorable pharmacokinetic properties. These findings provide scientific validation of the traditional uses of this plant and suggest its potential application as natural drugs.

Structure-based drug design, molecular dynamics simulation, ADMET, and quantum chemical studies of some thiazolinones targeting influenza neuraminidase

The genetic mutability of the influenza virus leads to the existence of drug-resistant strains which is dangerous, particularly with the lingering coronavirus disease (COVID-19). This necessitated the need for the search and discovery of more potential anti-influenza agents to avert future outbreaks. In furtherance of our previous in-silico studies on 5-benzyl-4-thiazolinones as anti-influenza neuraminidase (NA) inhibitors, molecule 11 was selected as the template scaffold for the structure-based drug design due to its good binding, pharmacokinetic profiling, and better NA inhibitory activity. As such, eighteen (18) new molecules (11a-r) were designed with better MolDock scores as compared with the template scaffold and the zanamivir reference drug. However, the dynamic stability of molecule 11a in the binding cavity of the NA target (3TI5) showed water-mediated hydrogen and hydrophobic bondings with the active residues such as Arg118, Ile149, Arg152, Ile222, Trp403, and Ile427 after the MD simulation for 100 ns. The drug-likeness and ADMET assessment of all designed molecules predicted non-violation of the stipulated thresholds of Lipinski's rule and good pharmacokinetic properties respectively. In addition, the quantum chemical calculations also suggested the significant chemical reactivity of molecules with their smaller band energy gap, high electrophilicity, high softness, and low hardness. The results obtained in this study proposed a reliable in-silico viewpoint for anti-influenza drug discovery and development.Communicated by Ramaswamy H. Sarma.

Binding Studies and Lead Generation of Pteridin-7(8H)-one Derivatives Targeting FLT3

Ligand modification by substituting chemical groups within the binding pocket is a popular strategy for kinase drug development. In this study, a series of pteridin-7(8H)-one derivatives targeting wild-type FMS-like tyrosine kinase-3 (FLT3) and its D835Y mutant (FL3D835Y) were studied using a combination of molecular modeling techniques, such as docking, molecular dynamics (MD), binding energy calculation, and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies. We determined the protein–ligand binding affinity by employing molecular mechanics Poisson–Boltzmann/generalized Born surface area (MM-PB/GBSA), fast pulling ligand (FPL) simulation, linear interaction energy (LIE), umbrella sampling (US), and free energy perturbation (FEP) scoring functions. The structure–activity relationship (SAR) study was conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), and the results were emphasized as a SAR scheme. In both the CoMFA and CoMSIA models, satisfactory correlation statistics were obtained between the observed and predicted inhibitory activity. The MD and SAR models were co-utilized to design several new compounds, and their inhibitory activities were anticipated using the CoMSIA model. The designed compounds with higher predicted pIC50 values than the most active compound were carried out for binding free energy evaluation to wild-type and mutant receptors using MM-PB/GBSA, LIE, and FEP methods.