Molecular Modeling Studies of C-Glycosylfavone Derivatives as GSK-3β Inhibitors Based on QSAR and Docking Analysis

Synthesis, spectroscopic analysis, DFT, docking, MD and antioxidant activity of tetrahydrocurcumin

In recent years, numerous researchers have made local chemical modifications to the structure of curcumin while its basic structure remains unchanged, thus, producing curcumin derivatives. In this article, tetrahydrocurcumin was obtained by hydrogenation of curcumin, DFT calculation and characterization at the theoretical level of B3LYP/6 -311++G(d,p) were carried out. The observed IR and Raman spectra are in good agreement with the theoretical spectra. The FMO and ESP of tetrahydrocurcumin are predicted. The interaction in the system is shown graphically and analyzed by IGMH. Compared with curcumin, tetrahydrocurcumin lacks the unsaturated C = C bond, which makes it more stable and more bioavailable. Molecular docking with antioxidant targets elucidated the ligand-protein interaction and molecular dynamics simulation showed the antioxidant activity of tetrahydrocurcumin. The antioxidant activity of tetrahydrocurcumin was proved by DPPH• and •OH radical scavenging experiments. In essence, these derivatives exhibit enhanced physiological activity in certain aspects compared to the original curcumin. Moreover, the computational pharmacology techniques lay a theoretical groundwork for the development and modification of high-efficiency, low-toxicity drugs that interface with various targets of curcumin in the future.Communicated by Ramaswamy H. Sarma.

Structure-based virtual screening of novel natural products as chalcone derivatives against SARS-CoV-2 Mpro

Coronavirus disease 2019 (COVID-19), which is caused by SARS-CoV-2, has spread quickly around the world, causing a global pandemic. It has infected more than 500 million people as of April 28, 2022. Much research has been reported to stop the virus from spreading, but there are currently no approved medicines to treat COVID-19. In this work, a dataset of 142 natural products collected from various medicinal plants was used to perform structure-based virtual screening (SBVS) through the combined application of molecular docking and molecular dynamics (MD) simulation methods. First, the dataset of compounds was optimized using the density functional theory (DFT) approach. The optimized compounds were then submitted to the first screening, which was done by the pKCM web server to look for drug-likeness and the PyRx to look for binding affinity. Among the 142 natural substances, 10 compounds were selected for docking validation. Compounds that interact with CYS145 and LEU141, the essential catalytic residues, as well as compounds with binding affinities less than -8.0 kcal/mol, are considered promising anti-SARS-CoV-2 drug candidates. The top-ranked compounds were then evaluated by MD simulations and MM-GBSA method. These results could help researchers come up with new natural compounds that could be used to treat SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Insight parameter drug design for human β-tryptase inhibition integrated molecular docking, QSAR, molecular dynamics simulation, and pharmacophore modelling studies of α-keto-[1,2,4]-oxadiazoles

Dengue hemorrhagic fever (DHF) is severe dengue with a hallmark of vascular leakage. β-tryptase has been found to promote vascular leakage in DHF patients, which could be a potential target for DHF treatment. This study aims to develop a theoretical background for designing and selecting human β-tryptase inhibitors through computational studies. Thirty-four α-keto-[1,2,3]-oxadiazoles scaffold-based compounds were used to generate 2D-QSAR models and for molecular docking studies with β-tryptase (PDB Code 4A6L). In addition, molecular dynamics (MD) simulation and molecular mechanics generalised born surface area (MM-GBSA) analysis on the binding of the reported most active compound, compound 11e, towards β-tryptase were performed. Finally, a structure-based pharmacophore model was generated. The selected 2D-QSAR models have statistically proven good models by internal and external validation as well as the y-randomization test. The docking results of compound 11e showed lower CDOCKER energy than the 4A6L co-crystallised ligand and a similar binding pattern as the 4A6L co-crystallised ligand. From molecular dynamics simulation, 4A6L in compound 11e bound state has RMSD below 2 Å throughout the 500 ns simulation, indicating the docked complex is stable. Besides, MM-GBSA analysis suggested the 4A6L-compound 11e docked complex (-66.04 Kcal/mol) is structurally as stable as the 4A6L-native ligand co-crystallized structure (-66.84 Kcal/mol). The best pharmacophore model identified features included hydrogen bond acceptor, ionic interaction, hydrophobic interaction, and aromatic ring, which contribute to the inhibitory potency of a compound. This study supplied insight and knowledge for developing novel chemical compounds with improved inhibition of β-tryptase.Communicated by Ramaswamy H. Sarma.

Homology modeling, virtual screening, molecular docking, molecular dynamic (MD) simulation, and ADMET approaches for identification of natural anti-Parkinson agents targeting MAO-B protein

Monoamine oxidase-B (MAO-B) is a flavin-dependent enzyme involved in various neurodegenerative disorders. Here, a dataset of 142 chalcone derivatives, collected from various natural plants, was screened by combining structure-based virtual screening and ADMET approaches. The goal is to discover novel natural chalcones as potential MAO-B inhibitors. With the help of the Gaussian 09.5 software, the 3D chemical structures of compounds were optimized using the DFT method. The 3D structure of the hMAO-B enzyme was built using the Modeller software. The optimized structures were subjected to virtual screening by Autodock Vina, implicated in PyRx software. Among the 142 natural substances, 43 were selected based on their binding affinity. Then, the pharmacokinetic proprieties and toxicity of these compounds were evaluated using ADMET analysis. Ten compounds were predicted to have ADMET characteristics with no side effects. The binding modes and interactions of the top selected compounds were then evaluated using AutoDock 4.2. Compounds P60 and P81 were found to be potential inhibitors of MAO-B compared to rasagiline, safinamid, and selegiline, the reference drugs. The stability of the selected compounds was confirmed by MD simulation. Based on this finding, compounds P60 and P81 could be considered potential hMAO-B inhibitors.

Development of novel monoamine oxidase B (MAO-B) inhibitors by combined application of docking-based alignment, 3D-QSAR, ADMET prediction, molecular dynamics simulation, and MM_GBSA binding free energy

Unsaturated ketone derivatives are known as inhibitors of monoamine oxidase B (MAO-B), a potential drug target of Parkinson's disease. Here, docking-based alignment, 3 D-QSAR (three-dimensional quantitative structure-activity relationship) studies, ADMET (absorption, distribution, metabolism, excretion, and toxicity) prediction, molecular dynamics (MD) simulation, and MM_GBSA binding free energy were performed on a novel series of MAO-B inhibitors. The objective is to predict new MAO-B inhibitors with high potency activity. The 3 D-QSAR models were created using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). Molecular docking findings indicated that compounds with strong inhibitory efficacy also had a high binding affinity. 3 D-QSAR studies showed the importance of steric, electrostatic, and H-bond acceptor fields on the inhibitory activity of MAO-B. Based on the appropriate 3 D-QSAR model, a new series of MAO-B inhibitors were predicted and their pharmacokinetic characteristics were evaluated using in silico ADMET prediction. All screened compounds show good oral bioavailability without any side effects. Moreover, the dynamic behavior and stability of the most active compounds were evaluated using MD simulations. The results showed that unsaturated ketone derivatives are stable and compact during the 100 ns of MD simulation. Finally, the binding free energy of complexes was determined using the MM_GBSA method; the findings indicated that the T1 compound is more stable (ΔG_binding = -409.506 KJ/mol) than the data set's highest active compound (ΔG_binding = -31.883 KJ/mol).Communicated by Ramaswamy H. Sarma.

Computer-aided drug design applied to a series of pyridinyl imidazole derivatives targeting p38α MAP kinase: 2D-QSAR, docking, MD simulation, and ADMET investigations

The p38 mitogen-activated protein kinase (MAPK) is a crucial target for chronic inflammation.