In silico exploration of aryl sulfonamide analogs as voltage-gated sodium channel 1.7 inhibitors by using 3D-QSAR, molecular docking study, and molecular dynamics simulations

Design and evaluation of novel triazole derivatives as potential anti-gout inhibitors: a comprehensive molecular modeling study

Introduction

Gout is the most common inflammatory arthritis, characterized by hyperuricemia, tophus formation, joint disease, and kidney stones. Uric acid, the final byproduct of purine catabolism, is eliminated via the kidneys and digestive system. Xanthine oxidase (XO) catalyzes the conversion of hypoxanthine and xanthine into uric acid, making XO inhibitors crucial for treating hyperuricemia and gout. Currently, three XO inhibitors are clinically used, showing significant efficacy. A molecular modeling study on triazole derivatives aims to identify novel XO inhibitors using 3D-QSAR, molecular docking, MD simulations, ADMET analysis, and DFT calculations. These computational approaches facilitate drug discovery while reducing research costs.

Methods

Our work focuses on a series of synthesized anti-xanthine oxidase inhibitors, aiming to develop new inhibitors. A computational study was carried out to identify the xanthine oxidase inhibitory structural features of a series of triazole inhibitors using computational method.

Results

A model based on CoMFA and CoMSIA/SEA has been built to predict new triazole derivatives.

Discussion

The optimal model established from CoMFA and CoMSIA/SEA was successfully evaluated for its predictive capability. Visualization of the contour maps of both models showed that modifying the substituents plays a key role in enhancing the biological activity of anti-gout inhibitors. Molecular docking results for complexes N°8-3NVY and N°22-3NVY showed scores of -7.22 kcal/mol and -8.36 kcal/mol, respectively, indicating substantial affinity for the enzyme. Complex N°8-3NVY forms two hydrogen bonds with SER 69 and ASN 71, three alkyl bonds with ALA 70, LEU 74, and ALA 75, and one Pi-Pi T-shaped bond with PHE 68. Complex N°22-3NVY forms three hydrogen bonds with HIS 99, ARG 29, and ILE 91, and one halogen bond with LEU 128 at 3.60 Å. A MD study revealed that the N°22-3NVY complex remained highly stable throughout the simulation. Therefore, we proposed six new molecules, their anti-gout inhibitory activities were predicted using two models, and they were evaluated for Lipinski's rule, and ADMET properties. The results show that both Pred 4 and Pred 5 have better pharmacokinetic properties than the height potent molecule in the studied series, making these two compounds valuable candidates for new anti-gout drugs. Subsequently, using DFT study to evaluate the chemical reactivity properties of these two proposed compounds, the energy gap results revealed that both molecules exhibit moderate chemical stability and reactivity.

Carbenoid-involved reactions integrated with scaffold-based screening generates a Nav1.7 inhibitor

The discovery of selective Nav1.7 inhibitors is a promising approach for developing anti-nociceptive drugs. In this study, we present a novel oxindole-based readily accessible library (OREAL), which is characterized by readily accessibility, unique chemical space, ideal drug-like properties, and structural diversity. We used a scaffold-based approach to screen the OREAL and discovered compound C4 as a potent Nav1.7 inhibitor. The bioactivity characterization of C4 reveals that it is a selective Nav1.7 inhibitor and effectively reverses Paclitaxel-induced neuropathic pain (PINP) in rodent models. Preliminary toxicology study shows C4 is negative to hERG. The consistent results of molecular docking and molecular simulations further support the reasonability of the in-silico screening and show the insight of the binding mode of C4. Our discovery of C4 paves the way for pushing the Nav1.7-based anti-nociceptive drugs forward to the clinic.

Design of potential anti-cancer agents as COX-2 inhibitors, using 3D-QSAR modeling, molecular docking, oral bioavailability proprieties, and molecular dynamics simulation

Modeling the structural properties of novel morpholine-bearing 1, 5-diaryl-diazole derivatives as potent COX-2 inhibitor, two proposed models based on CoMFA and CoMSIA were evaluated by external and internal validation methods. Partial least squares analysis produced statistically significant models with Q 2 values of 0.668 and 0.652 for CoMFA and CoMSIA, respectively, and also a significant non-validated correlation coefficient R² with values of 0.882 and 0.878 for CoMFA and CoMSIA, respectively. Both models met the requirements of Golbraikh and Tropsha, which means that both models are consistent with all validation techniques. Analysis of the CoMFA and CoMSIA contribution maps and molecular docking revealed that the R1 substituent has a very significant effect on their biological activity. The most active molecules were evaluated for their thermodynamic stability by performing MD simulations for 100 ns; it was revealed that the designed macromolecular ligand complex with 3LN1 protein exhibits a high degree of structural and conformational stability. Based on these results, we predicted newly designed compounds, which have acceptable oral bioavailability properties and would have high synthetic accessibility.

A biochemistry-oriented drug design: synthesis, anticancer activity, enzymes inhibition, molecular docking studies of novel 1,2,4-triazole derivatives

In this study, we researched the reactions of 5-(5-bromofuran-2-yl)-4-methyl-1,2,4-triazole-3-thiol and 5-thiophene-(3-ylmethyl)-4R-1,2,4-triazole-3-thiols with some halogen-containing compounds, a number of new compounds were synthesized (1.1-1.5 and 2.1-2.8). These compounds showed excellent to good inhibitory activities on acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) enzymes. For obtaining the effects of these compounds on AChE and BChE enzymes were determined spectrophotometrically according to Ellman. IC50 values of these enzymes were ranging between 1.63 and 17.68 nM for AChE and 8.71 and 84.02 nM for BChE. After, prostate cancer is the second leading cause of cancer-related mortality for men over the age of 65 in developed countries. Current treatment options remain limited in the treatment of advanced-stage prostate cancer leading to biochemical recurrence in almost 40% of the patients. Therefore, there is an urgent need for development of novel therapeutic tools for treatment of prostate cancer patients. In this study, we aimed at analyzing the potential of all compounds against prostate cancer cells. We found that, of the tested compounds, 2.1, 2.2 and 2.3 showed significant cytotoxic activities against PC3 prostate cancer cells, although their effect on the viability of normal prostate cells was limited. These findings suggest their selective targeting potential for prostate cancer cells and offer them as candidate therapeutic agents against prostate cancer. The inhibitory activities of some chemical compounds, such as (1.1-1.5 and 2.1-2.8) were assessed by performing the molecular docking study in the presence of AChE, BChE and prostate cancer protein. MM/GBSA methods are calculated binding free energy. Finally, ADME/T analysis was performed to examine the drug properties of the 13 studied molecules.Communicated by Ramaswamy H. Sarma.

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.

In-silico investigation of 4-nitro-N-1H-pyrazol-3-ylbenzamide towards its potential use against SARS-CoV-2: a DFT, molecular docking and molecular dynamics study

In the present research work, we report the synthesis and characterization of novel pyrazole derivative obtained by the condensation reaction of 4-nitro benzaldehyde group with one equivalent of the 2-amino pyrazole yielding 4-nitro-N-1H-pyrazol-3-ylbenzamide with high yield. The two symmetry-independent molecules (molecule A and molecule B) differ about the central C-N bond, with the dihedral angles between the pyrazole ring system and the nitrobenzene ring being 13.90° and 18.64°, respectively. By optimizing the symmetry-independent dimer molecules, the rotational barrier between the conformers is found to be within the 2.5-5.5 kcal/mol range. QTAIM and RDG based NCI isosurface revealed the presence of strong N-H…N and C-H…O hydrogen bonds which stabilize the two independent centrosymmetric inversion-related dimers. Further, weak and short directional interactions such as C-H…N, H…H and C-H…π were also analyzed systematically using various topological parameters. The compound is found to adhere to the Lipinski's rule of five and exhibit good pharmacokinetic properties. The results of molecular docking studies performed against SARS-CoV-2 virus main protease (PDB IDs: 6LU7, 6W9C and 6WQF) revealed that the compound showed better docking scores. Molecular docking studies verified the inhibition activity of the synthesized novel compound. Finally, the binding free energy and contributed energies were calculated using MM-GBSA method. The 6LU7-ligand complex showed highest binding free energy and among all other interactions, the contributions of the covalent binding and van der Waals energy are found to be significant.Communicated by Ramaswamy H. Sarma.

Novel Dual-Target Kinase Inhibitors of EGFR and ALK Were Designed, Synthesized, and Induced Cell Apoptosis in Non-Small Cell Lung Cancer

ALK-positive NSCLC coexisting with EGFR mutations is a frequently occurring clinical phenomenon. Targeting ALK and EGFR simultaneously may be an effective way to treat these cancer patients. In this study, we designed and synthesized ten new dual-target EGFR/ALK inhibitors. Among them, the optimal compound 9j exhibited good activity with IC₅₀ values of 0.07829 ± 0.03 μM and 0.08183 ± 0.02 μM against H1975 (EGFR ^T790M/L858R) and H2228 (EML4-ALK) cells, respectively. Immunofluorescence assays indicated that the compound could simultaneously inhibit the expression of phosphorylated EGFR and ALK proteins. A kinase assay demonstrated that compound 9j could inhibit both EGFR and ALK kinases; thus, exerting an antitumor effect. Additionally, compound 9j induced apoptosis in a dose-dependent manner and inhibited the invasion and migration of tumor cells. All of these results indicate that 9j is worthy of further study.

Design of novel anti-cancer drugs targeting TRKs inhibitors based 3D QSAR, molecular docking and molecular dynamics simulation

Tropomyosin receptor kinase (TRK) enzymes are responsible for different types of tumors caused by neurotrophic tyrosine receptor kinase gene fusion and have been identified as an effective target for anticancer therapy. The study of the mechanism between polo-like kinase (PLKs) and pyrazol inhibitors was performed using 3D-QSAR modeling, molecular docking, and MD simulations in order to design high-activity inhibitors. The HQSAR (Q2 = 0.793, R2 = 0.917, R2ext = 0.961), CoMFA (Q2 = 0.582, R2 = 0.722, R2ext = 0.951), CoMSIA/SE (Q2 = 0.603, R2 = 0.801, R2ext = 0.849), and Topomer CoMFA (Q2 = 0.726, R2 = 0.992, R2ext = 0.717) showed good reliability and predictability. All models have been successfully tested by external validation, so all five established models are reliable. The analysis of the different contour maps of different models gives structural information to improve the inhibitory function. Molecular docking results show that the amino acids Met 592, GLU 590, LEU 657, VAL 524, and PHE 589 are the active sites of the tropomyosin receptor TRKs. The results obtained by MD showed that compound 19i could form a more stable complex protein (PDB id: 5KVT). Based on these results, we developed new compounds and their expected inhibitory activities. The results of physicochemical and ADME-Tox properties showed that the four proposed molecules are orally bioavailable, and they are not toxic in the Ames test. Thus, these results would provide modeling information that could help experimental researchers find TRK type I inhibitors more efficiently.Communicated by Ramaswamy H. Sarma.

Design, synthesis, and biological evaluation of low-toxic lappaconitine derivatives as potential analgesics

The C₁₈-diterpenoid alkaloid lappaconitine (LA) is a non-addictive analgesic used in China. The toxicity (LD₅₀ = 11.7 mg/kg) limits its application. Two series of LA derivatives, including amides and sulfonamides (1-93), were designed and synthesized by modification on their C4 acetamidobenzoate side chains in this work. In vivo analgesic activity and toxicity of all derivatives were evaluated, and the structure-activity relationship was summarized. Six lead compounds (35, 36, 39, 49, 70, and 89) exhibited approximate analgesic activity to LA but with significantly reduced toxicity. The therapeutic index of these compounds is 14-30 times that of LA. In vivo metabolism study of the lead compounds 39, 49, 70, and 89 were conducted by UPLC-MS^E, indicating the reason for the low toxicity of the potential derivatives might be they are difficult to metabolize to toxic metabolite N-deacetyllappaconitine compared to LA. The effects of lead compounds on sodium channels and hERG channels were also studied by ion channel reader (ICR) which further revealed their analgesic and toxicity-attenuating mechanisms. Sodium channel assay revealed that the analgesic mechanism of these lead compounds was inhibiting the Na_v 1.7 channels. Taken together, compound 39 was provided as a new analgesic lead compound with significantly low toxicity and comparable activity to LA.

Design, Synthesis, and Evaluation of New Mesenchymal–Epithelial Transition Factor (c-Met) Kinase Inhibitors with Dual Chiral Centers

A series of tepotinib derivatives with two chiral centers was designed, synthesized, and evaluated as anticancer agents. The optimal compound (R, S)-12a strongly exhibited antiproliferative activity against MHCC97H cell lines with an IC50 value of 0.002 μM, compared to tepotinib (IC50 = 0.013 μM). Mechanistic studies revealed that compound (R, S)-12a significantly inhibited c-Met activation, as well as the downstream AKT signaling pathway, and suppressed wound closure. Moreover, compound (R, S)-12a induced cellular apoptosis and cell cycle arrest at the G1 phase in a dose-dependent fashion.

Molecular insights into Mmpl3 leads to the development of novel indole-2-carboxamides as antitubercular agents

Tuberculosis (TB) is an air-borne infectious disease and is the leading cause of death among all infectious diseases globally. The current treatment regimen for TB is overtly long and patient non-compliance often leads to drug resistant TB resulting in a need to develop new drugs that will act via novel mechanisms. In this research work, we selected Mycobacterium membrane protein large (MmpL3) as the drug target and indole-2-carboximide as our molecule of interest for further designing new molecules. A homology model was prepared for the Mycobacterium tuberculosis MmpL3 from the crystal structure of Mycobacterium smegmatis MmpL3. A series of indoles which are known to be MmpL3 inhibitors were docked in the prepared protein and the binding site properties were identified. Based on that, 10 molecules were designed and synthesized and their antitubercular activities evaluated. We identified four hits among which the highest potency candidate possessed a minimum inhibitory concentration (MIC) of 1.56 μM at 2-weeks. Finally, molecular dynamics simulation studies were done with 3b and a previously reported MmpL3 inhibitor to understand the intricacies of their binding in real time and to correlate the experimental findings with the simulation data.

Genomic diversity and molecular dynamics interaction on mutational variances among RB domains of SARS-CoV-2 interplay drug inactivation

The scientific community has been releasing whole genomic sequences of SARS-CoV-2 to facilitate the investigation of molecular features and evolutionary history. We retrieved 36 genomes of 18 prevalent countries of Asia, Europe and America for genomic diversity and mutational analysis. Besides, we studied mutations in the RBD regions of Spike (S) proteins to analyze the drug efficiency against these mutations. In this research, phylogenenetic analysis, evolutionary modeling, substitution pattern analysis, molecular docking, dynamics simulation, etc. were performed. The genomic sequences showed >99% similarity with the reference sequence of China.TN93 + G was predicted as a best nucleotide substitution model. It was revealed that effective transition from the co-existing SARS genome to the SARS-CoV-2 and a noticeable positive selection in the SARS-CoV-2 genomes occurred. Moreover, three mutations in RBD domain, Val/ Phe367, Val/ Leu 382 and Ala/ Val522, were discovered in the genomes from Netherland, Bangladesh and the USA, respectively. Molecular docking and dynamics study showed RBD with mutation Val/Leu382 had the lowest binding affinity with remdesivir. In conclusion, the SARS-CoV-2 genomes are similar, but multiple degrees of transitions and transversions occurred. The mutations cause a significant conformational change, which are needed to be investigated during drug and vaccine development.

CK-3, A Novel Methsulfonyl Pyridine Derivative, Suppresses Hepatocellular Carcinoma Proliferation and Invasion by Blocking the PI3K/AKT/mTOR and MAPK/ERK Pathways

Hepatocellular carcinoma (HCC) is an aggressive tumor with a poor prognosis that highly expresses phosphatidylinositol 3-kinase (PI3K) and mitogen-activated protein kinase (ERK). The PI3K/AKT/mTOR and MAPK/ERK signaling pathways play a crucial role in HCC tumor formation, cell cycle, apoptosis and survival. However, no effective targeted therapies against these pathways is available, mainly due to the extensive and complex negative feedback loops between them. Here we used CK-3, a dual blocker of the PI3K/AKT/mTOR and MAPK/ERK pathways, against HCC cell lines to verify its anti-tumor activity in vitro. CK-3 exhibited cytotoxic activity against HCC, as demonstrated with MTT and colony formation assays. The anti-metastatic potential of CK-3 was demonstrated with wound healing and cell invasion assays. The ability of CK-3 to block both the PI3K/AKT/mTOR and MAPK/ERK pathways was also confirmed. CK-3 induced the apoptosis of Hep3B cells, while Bel7402 cells died via mitotic catastrophe (MC). Oral administration of CK-3 also inhibited the subcutaneous growth of BEL7402 cells in nude mice. Simultaneous PI3K/AKT/mTOR and MAPK/ERK pathway inhibition with CK-3 may be superior to single pathway monotherapies by inhibiting their feedback-regulation, and represents a potential treatment for HCC.

A Multi-target Drug Designing for BTK, MMP9, Proteasome and TAK1 for the Clinical Treatment of Mantle Cell Lymphoma

BACKGROUND

Mantle cell lymphoma (MCL) is a type of non-Hodgkin lymphoma characterized by the mutation and overexpression of the cyclin D1 protein by the reciprocal chromosomal translocation t(11;14)(q13:q32).

AIM

The present study aims to identify potential inhibition of MMP9, Proteasome, BTK, and TAK1 and determine the most suitable and effective protein target for the MCL.

METHODOLOGY

Nine known inhibitors for MMP9, 24 for proteasome, 15 for BTK and 14 for TAK1 were screened. SB-3CT (PubChem ID: 9883002), oprozomib (PubChem ID: 25067547), zanubrutinib (PubChem ID: 135565884) and TAK1 inhibitor (PubChem ID: 66760355) were recognized as drugs with high binding capacity with their respective protein receptors. 41, 72, 102 and 3 virtual screened compounds were obtained after the similarity search with compound (PubChem ID:102173753), PubChem compound SCHEMBL15569297 (PubChem ID:72374403), PubChem compound SCHEMBL17075298 (PubChem ID:136970120) and compound CID: 71814473 with best virtual screened compounds.

RESULT

MMP9 inhibitors show commendable affinity and good interaction profile of compound holding PubChem ID:102173753 over the most effective established inhibitor SB-3CT. The pharmacophore study of the best virtual screened compound reveals its high efficacy based on various interactions. The virtual screened compound's better affinity with the target MMP9 protein was deduced using toxicity and integration profile studies.

CONCLUSION

Based on the ADMET profile, the compound (PubChem ID: 102173753) could be a potent drug for MCL treatment. Similar to the established SB-3CT, the compound was non-toxic with LD50 values for both the compounds lying in the same range.

In silico development of potential therapeutic for the pain treatment by inhibiting voltage-gated sodium channel 1.7

The voltage-gated sodium channel Nav1.7 can be considered as a promising target for the treatment of pain. This research presents conformational-independent and 3D field-based QSAR modeling for a series of aryl sulfonamide acting as Nav1.7 inhibitors. As descriptors used for building conformation-independent QSAR models, SMILES notation and local invariants of the molecular graph were used with the Monte Carlo optimization method as a model developer. Different statistical methods, including the index of ideality of correlation, were used to test the quality of the developed models, robustness and predictability and obtained results were good. Obtained results indicate that there is a very good correlation between 3D QSAR and conformation-independent models. Molecular fragments that account for the increase/decrease of a studied activity were defined and used for the computer-aided design of new compounds as potential analgesics. The final evaluation of the developed QSAR models and designed inhibitors were carried out using molecular docking studies, bringing to light an excellent correlation with the QSAR modeling results.

Identification and Pharmacokinetics of Quinone Reductase 2 Inhibitors after Oral Administration of Garcinia mangostana L. Extract in Rat by LC-MS/MS

Garcinia mangostana L. (mangosteen) is a famous tropical fruit that contains a large number of xanthones. Regular consumption of mangosteen may confer health benefits and prevent some diseases, such as malaria. Quinone reductase 2 (QR-2) is a cytosolic enzyme found in human red blood cells, and it is becoming a target for chemoprevention because it is involved in the mechanisms of several diseases, including malaria. To understand whether the xanthones present in mangosteen might inhibit the activity of QR-2, blood samples were collected from rat following the oral administration of mangosteen extract and then incubated with QR-2 followed by UF-HPLC-QTOF/MS analysis to rapidly screen for and identify the QR-2-inhibiting xanthones. A total of 16 xanthones were identified, and six of these (α-mangostin, γ-mangostin, 8-deoxyartanin, 1,3,7-trihydroxy-2,8-di(3-methylbut-2-enyl)xanthone, garcinone E, and 9-hydroxycalabaxanthone) were subjected to QR-2 inhibition assay. γ-Mangostin exhibited the strongest inhibition, achieving an IC₅₀ value of 3.82 ± 0.51 μM. Its interaction with QR-2 was found to involve hydrogen bond and arene-arene interaction as revealed by molecular docking. The present study could provide new insight into the potential application of mangosteen as functional food ingredients for inhibiting the activity of QR-2. However, the extent of daily intake of mangosteen required and the exact contribution of mangosteen to the prevention and treatment of malaria remain subjects of further study.

Synthesis, extracorporeal nephrotoxicity, and 3D-QSAR of andrographolide derivatives

Andrographolide is a traditional Chinese medicine monomer with many pharmacological activities, but has potential nephrotoxicity. Here, we aim to investigate the relationship between modification of andrographolide structure and its nephrotoxicity. Twenty-three andrographolide derivatives were synthesized, and their structures were confirmed by ¹ H-NMR and HRMS. Nephrotoxicity of these compounds against human renal tubular epithelial (HK-2) cells was evaluated using the MTT assay. The results indicated that most of them had significantly reduced nephrotoxicity, especially compounds III, V, and IX_c , with IC₅₀ values of 1,985, 1,300, and 806.9 μmol/L, respectively, which were obviously superior to andrographolide (IC₅₀ 30.60 μmol/L). However, compounds I_a -I_f (IC₅₀ values < 30 μmol/L), the 14-OH derivatives of andrographolide, showed higher nephrotoxicity than that of andrographolide. Three-dimensional quantitative structure-activity relationship (3D-QSAR) models of COMFA and COMSIA were established (COMFA: q²  = 0.639, r²  = 0.951; COMSIA: q²  = 0.569, r²  = 0.857). This model allowed proposing five new compounds with lower theoretical nephrotoxicity, which would be worthwhile to synthesize and evaluate. We believe that predicted models will help us to understand the structural modification requirements of andrographolide to reduce the nephrotoxicity, and further investigations will be needed to determine the mechanism involved in the effect of less nephrotoxic compounds.

Based on the Virtual Screening of Multiple Pharmacophores, Docking and Molecular Dynamics Simulation Approaches toward the Discovery of Novel HPPD Inhibitors

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is an iron-dependent non-heme oxygenase involved in the catabolic pathway of tyrosine, which is an important enzyme in the transformation of 4-hydroxyphenylpyruvic acid to homogentisic acid, and thus being considered as herbicide target. Within this study, a set of multiple structure-based pharmacophore models for HPPD inhibitors were developed. The ZINC and natural product database were virtually screened, and 29 compounds were obtained. The binding mode of HPPD and its inhibitors obtained through molecular docking study showed that the residues of Phe424, Phe381, His308, His226, Gln307 and Glu394 were crucial for activity. Molecular-mechanics-generalized born surface area (MM/GBSA) results showed that the coulomb force, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. These efforts will greatly contribute to design novel and effective HPPD inhibitory herbicides.

Molecular docking, 3D-QSAR, and molecular dynamics simulations of thieno[3,2-b]pyrrole derivatives against anticancer targets of KDM1A/LSD1

Lysine-specific demethylase 1 (LSD1) is a histone-modifying enzyme, which has been proposed as a promising target for anticancer drug development. Extensive research on LSD1 inhibitors has been performed since its discovery. In order to get more information for lead identification and optimization, we carried out a molecular modeling study on a set of 43 thieno[3,2-b]pyrrole competitive inhibitors of LSD1 using three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations. Based on the co-crystallized conformer-based alignment (CCBA) method, 3D-QSAR model of thieno[3,2-b]pyrrole derivatives as LSD1 inhibitors was established. The significant statistics (q² = 0.595, r² = 0.959, r²_pred = 0.846) of the 3D-QSAR indicated the good predictive power and statistical reliability of this model. Based on the corresponding contour maps six LSD1 inhibitors were designed and their activities were predicted by 3D-QSAR model. Meanwhile, molecular docking was performed to simulate the probable binding modes between ligands and LSD1 protein. The molecular interactions mainly contributions to the binding affinity for LSD1 inhibitions were further supplemented by 100 ns MD simulations and binding free energy calculation.

Gaining deeper insights into the surface binding of bedaquiline analogues with the ATP synthase subunit C of Mycobacterium tuberculosis using molecular docking, molecular dynamics simulation and 3D-QSAR techniques

This computational study exclusively illustrates the key molecular features of bedaquiline and its analogues required for binding to mycobacterial ATP synthase.

Identification of potential platelet-derived growth factor receptor α inhibitors by computational screening and binding simulations

Platelet-derived growth factor receptor α (PDGFRα) is considered as a promising target for the treatment of fibrotic diseases. In this study, two types of pharmacophore model, which generated by ligand-based and receptor-based method, were put forward to identify novel chemical entities as PDGFRα inhibitors. It was found that some pharmacophore characteristics established by the two approaches overlap each other. In order to elucidate detailed interactions, representative molecules were selected to predict the conformations and binding modes of the molecules by molecular docking method. The calculation results of binding free energy illustrated that the van der Waals energy and nonpolar solvation were the most prominent contribution to the interactions between the inhibitors and PDGFRα. To further verify the accuracy of the docking results and the stability of the complexes system, the binding modes of two potent PDGFRα inhibitors were examined by 100 ns molecular dynamics simulations. Herein, we reported the basic structural requirements of PDGFRα inhibitors for the first time through molecular docking and molecular dynamics simulations. Subsequently, the two pharmacophore models were used for virtual screening to query potential active molecules from Food and Drug Administration approved database. The hit molecules here might provide additional scaffolds for further optimization of PDGFRα inhibitors.

Computational investigation reveals Picrasidine C as selective PPARα lead: binding pattern, selectivity mechanism and ADME/tox profile

Natural products and their derivatives have been recognized as an important source of therapeutic agents for many years. Previously we isolated a dimeric β-carboline-type alkaloid Picrasidine C from the root of Picrasma quassioides as subtype-selective peroxisome proliferator-activated receptor α (PPARα) agonist. In order to modify this natural product for better affinity and druggability, we investigated a series of properties exhibited by Picrasidine C, such as its binding mode with PPARα, the selectivity mechanism over PPARγ, as well as ADME/Tox profile through computational methods including sequence alignment, molecular docking, pharmacophore modeling and molecular dynamics simulations. The detailed information of binding pattern and affinity for Picrasidine C elucidated here will be valuable for chemical modification. Besides, the steric hindrance of residue Phe363 in PPARγ pocket was speculated as the main isoform selectivity mechanism for Picrasidine C, which would be helpful for the design of selective derivatives. ADME/Tox prediction was conducted to avoid potential undesirable pharmacokinetic properties for reducing the risk of failure. Finally, novel skeletons were derived from lead compound by core hopping method, validated through molecular dynamic simulations and MM-GBSA calculation. In short, the information obtained from computational strategy would be valuable for us to find more potent, safe and selective PPARα agonists during structural optimization.

In silico studies on p21-activated kinase 4 inhibitors: comprehensive application of 3D-QSAR analysis, molecular docking, molecular dynamics simulations, and MM-GBSA calculation

P21-activated kinase 4 (PAK4) is a serine/threonine protein kinase, which is associated with many cancer diseases, and thus being considered as a potential drug target. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations were performed to explore the structure-activity relationship of a series of pyrropyrazole PAK4 inhibitors. The statistical parameters of comparative molecular field analysis (CoMFA, Q ² = 0.837, R ² = 0.990, and R ² _pred = 0.967) and comparative molecular similarity indices analysis (CoMSIA, Q ² = 0.720, R ² = 0.972, and R ² _pred = 0.946) were obtained from 3D-QSAR model, which exhibited good predictive ability and significant statistical reliability. The binding mode of PAK4 with its inhibitors was obtained through molecular docking study, which indicated that the residues of GLU396, LEU398, LYS350, and ASP458 were important for activity. Molecular mechanics generalized born surface area (MM-GBSA) method was performed to calculate the binding free energy, which indicated that the coulomb, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. Furthermore, through 100 ns MD simulations, we obtained the key amino acid residues and the types of interactions they participated in. Based on the constructed 3D-QSAR model, some novel pyrropyrazole derivatives targeting PAK4 were designed with improved predicted activities. Pharmacokinetic and toxicity predictions of the designed PAK4 inhibitors were obtained by the pkCSM, indicating these compounds had better absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. Above research provided a valuable insight for developing novel and effective pyrropyrazole compounds targeting PAK4.

Molecular determinants for ligand binding at Nav1.4 and Nav1.7 channels: Experimental affinity results analyzed by molecular modeling

Here, we focused on exploring the selectivity mechanism against Nav1.7 over Nav1.4 due to different binding modes of two selected inhibitors. By the superposition of Nav1.7 and Nav1.4 proteins, we selected the most homologous chain of Nav1.7 with Nav1.4, defining the active site of Nav1.4-VSD4 based on the aryl sulfonamide binding site of Nav1.7-VSD4. Comparison of the conformations exhibited by Tyr1386 (Nav1.4) and Tyr1537 (Nav1.7) suggested that the steric hindrance caused by Tyr1386 owned primary influence on inhibition selectivity, which was further verified through molecular docking and MD simulation of two representative inhibitors. Our finding would be helpful for discovery of selective Nav1.7 inhibitors over Nav1.4.

Design, Synthesis and Molecular Docking Analysis of Flavonoid Derivatives as Potential Telomerase Inhibitors

Based on the structural scaffolds of natural products, two series of flavonoid derivatives, for a total of twelve compounds, were designed and synthesized as potential human telomerase inhibitors. Using a modified TRAP-PCR assay, compound 5c exhibited the most potent inhibitory activity against human telomerase with an IC₅₀ value of less than 50 μM. In vitro, the results demonstrated that compound 5c had potent anticancer activity against five classes of tumor cell lines. The molecular docking and molecular dynamics analyses binding to the human telomerase holoenzyme were performed to elucidate the binding mode of active compound 5c. This finding helps the rational design of more potent telomerase inhibitors based on the structural scaffolds of natural products.

Computational investigations of gram-negative bacteria phosphopantetheine adenylyltransferase inhibitors using 3D-QSAR, molecular docking and molecular dynamic simulations

Phosphopantetheine adenylyltransferase (PPAT) has been recognized as a promising target to develop novel antimicrobial agents, which is a hexameric enzyme that catalyzes the penultimate step in coenzyme A biosynthesis. In this work, molecular modeling study was performed with a series of PPAT inhibitors using molecular docking, three-dimensional qualitative structure-activity relationship (3D-QSAR) and molecular dynamic (MD) simulations to reveal the structural determinants for their bioactivities. Molecular docking study was applied to understand the binding mode of PPAT with its inhibitors. Subsequently, 3D-QSAR model was constructed to find the features required for different substituents on the scaffolds. For the best comparative molecular field analysis (CoMFA) model, the Q² and R² values of which were calculated as 0.702 and 0.989, while they were calculated as 0.767 and 0.983 for the best comparative molecular similarity index analysis model. The statistical data verified the significance and accuracy of our 3D-QSAR models. Furthermore, MD simulations were carried out to evaluate the stability of the receptor-ligand contacts in physiological conditions, and the results were consistent with molecular docking studies and 3D-QSAR contour map analysis. Binding free energy was calculated with molecular mechanics generalized born surface area approach, the result of which coincided well with bioactivities and demonstrated that van der Waals accounted for the largest portion. Overall, our study provided a valuable insight for further research work on the recognition of potent PPAT inhibitors.Communicated by Ramaswamy H. Sarma.

Towards further Understanding the Structural Requirements of Combretastatin- like Chalcones as Inhibitors of Microtubule Polymerization

BACKGROUND

Microtubules are dynamic filamentous cytoskeletal structures which play several key roles in cell proliferation and trafficking. They are supposed to contribute in the development of important therapeutic targeting tumor cells. Chalcones are important group of natural compounds abundantly found in fruits & vegetables that are known to possess anticancer activity. We have used QSAR and docking studies to understand the structural requirement of chalcones for understanding the mechanism of microtubule polymerization inhibition.

METHODS

Three dimensional (3D) QSAR (CoMFA and CoMSIA), pharmacophore mapping and molecular docking studies were performed for the generation of structure activity relationship of combretastatin-like chalcones through statistical models and contour maps.

RESULTS

Structure activity relationship revealed that substitution of electrostatic, steric and donor groups may enhance the biological activity of compounds as inhibitors of microtubule polymerization. From the docking study, it was clear that compounds bind at the active site of tubulin protein.

CONCLUSION

The given strategies of modelling could be an encouraging way for designing more potent compounds as well as for the elucidation of protein-ligand interaction.