Molecular docking, validation, dynamics simulations, and pharmacokinetic prediction of natural compounds against the SARS-CoV-2 main-protease

Investigation of the molecular structure of CHBP, biological activities and SARS-CoV-2 protein binding interaction by molecular and biomolecular spectroscopy approaches

The objective of this investigation is to learn more about the structural, electrical, spectroscopic, and physiochemical characteristics of biologically active cyano-4'-hydroxybiphenyl (CHBP). The title molecule's optimized conformational analysis was computed using the DFT/B3LYP/6-311++G (d, p) level of theory. The observed wavenumbers were compared with theoretical FT-IR and FT-Raman spectra. 1H and 13C NMR experimental spectra in CDCl3 solution (solvent phase) were recorded and the chemical shift was calculated. NBO analysis was used to examine the transfer of charge as well as the intermolecular and intramolecular bonding of orbitals. The TD-DFT (time-dependent DFT) approach was used to estimate theoretical values for both the gas and solvent (ethanol) in the corresponding transitional research, which was conducted using UV-Vis's spectra. Energy gap (Eg = 0.26764 eV) implies that the strong potential for charge transfer, and the stability of the CHBP compound. CHBP compound's has bioactive nature, its drug-likeness and biological properties were evaluated. The predicted topological polar surface area of 44.02 \AA2 for the molecule falls within the permissible range of < 140 \AA2. Based on the docking results, the most stable docking score value is -6.84 kcal/mol. In that interaction, MET 165 affects both phenyl rings in a pi-sulphur fashion and a single bond hydrogen with protein moieties GLN 192. This suggests that the pi-alkyl in PRO 168 is a hydroxyl substitutional ring. Our findings demonstrate the CHBP compound is a good inhibitor against the SAR COVID-19 viral protein.

MP-13, a novel chimeric peptide of morphiceptin and pepcan-9, produces potent antinociception with limited side effects

Pharmacological investigations have substantiated the potential of bifunctional opioid/cannabinoid agonists in delivering potent analgesia while minimizing adverse reactions. Peptide modulators of cannabinoid receptors, known as pepcans, have been investigated before. In this study, we designed a series of chimeric peptides based on pepcans and morphiceptin (YPFP-NH2). Here, we combined injections of pepcans and morphiceptin to investigate the combination treatment of opioids and cannabis and compared the analgesic effect with chimeric compounds. Subsequently, we employed computational docking to screen the compounds against opioid and cannabinoid receptors, along with an acute pain model, to identify the most promising peptide. Among these peptides, MP-13, a morphiceptin and pepcan-9 (PVNFKLLSH) construct, exhibited superior supraspinal analgesic efficacy in the tail-flick test, with an ED50 value at 1.43 nmol/mouse, outperforming its parent peptides and other chimeric analogs. Additionally, MP-13 displayed potent analgesic activity mediated by mu-opioid receptor (MOR), delta-opioid receptor (DOR), and cannabinoid type 1 (CB1) receptor pathways. Furthermore, MP-13 did not induce psychological dependence and gastrointestinal motility inhibition at the effective analgesic doses, and it maintained non-tolerance-forming antinociception throughout a 7-day treatment regimen, with an unaltered count of microglial cells in the periaqueductal gray region, supporting this observation. Moreover, intracerebroventricular administration of MP-13 demonstrated dose-dependent antinociception in murine models of neuropathic, inflammatory, and visceral pain. Our findings provide promising insights for the development of opioid/cannabinoid peptide agonists, addressing a crucial gap in the field and holding significant potential for future research and development. PERSPECTIVE: This article offers insights into the combination treatment of pepcans with morphiceptin. Among the chimeric peptides, MP-13 exhibited potent analgesic effects in a series of preclinical pain models with a favorable side-effect profile.

In silico analysis of Diosmetin as an effective chemopreventive agent against prostate cancer: molecular docking, validation, dynamic simulation and pharmacokinetic prediction-based studies

Prostate cancer is the second most dangerous cancer type worldwide. While various treatment options are present i.e. agonists and antagonists, their utilization leads to adverse effects and due to this resistance developing, ultimately the outcome is remission. So, to overcome this issue, we have undertaken an in-silico investigation to identify promising and unique flavonoid candidates for combating prostate cancer. Using GOLD software, the study assessed the effectiveness of 560 natural secondary polyphenols against CDKN2. Protein Data Bank was used to retrieve the 3D crystal structure of CDKN2 (PDB Id: 4EK3) and we retrieved the structure of selected secondary polyphenols from the PubChem database. The compound Diosmetin shows the highest GOLD score with the selected Protein i.e. CDKN2 which is 58.72. To better understand the 2-dimensional and 3-dimensional interactions, the interacting amino acid residues were visualised using Discovery Studio 3.5 and Maestro 13.5. Using Schrodinger-Glide, the Diosmetin and CDKN2 were re-docked, and decoy ligands were docked to CDKN2, which was used to further ascertain the study. The ligands with the highest Gold score were forecasted for pharmacokinetics characteristics, and the results were tabulated and analysed. Utilising the Gromacs software and Desmond packages, 100 ns of Diosmetin molecular dynamics simulations were run to evaluate the structural persistence and variations of protein-ligand complexes. Additionally, our investigation revealed that Diosmetin had a better binding affinity with CDKN2 measuring 58.72, and it also showed remarkable stability across a 100-ns simulation. Thus, following in-vitro and in-vivo clinical studies, diosmetin might lead to the Prostate regimen.Communicated by Ramaswamy H. Sarma.

Exploring α, β-unsaturated carbonyl compounds against bacterial efflux pumps via computational approach

Antibiotic resistance has become a pressing global health crisis, with bacterial infections increasingly difficult to treat due to the emergence of multidrug resistance. This study aims to identify potential chalcone molecules that interact with two key multidrug efflux pumps, AcrB and EmrD, of Escherichia coli, using advanced computational tools. In silico ADMET (absorption, distribution, metabolism, excretion, and toxicity), drug-likeness prediction, molecular docking, and molecular dynamics simulation analyses were conducted on a ligand library comprising 100 chalcone compounds against AcrB (PDB: 4DX5) and EmrD (PDB: 2GFP). The results demonstrated that Elastichalcone A (PubChem CID 102103730) exhibited a remarkable binding affinity of -9.9 kcal/mol against AcrB, while 4'-methoxy-4-hydroxychalcone (PubChem CID 5927890) displayed a binding affinity of -9.8 kcal/mol against EmrD. Both ligands satisfied drug-likeness rules and possessed favorable pharmacokinetic profiles. Molecular dynamics simulation of the AcrB-Elastichalcone A complex remained stable over 100 ns, with minimal fluctuations in root-mean-square deviation and root-mean-square fluctuation. The screened ligand library demonstrated good drug-likeness and pharmacokinetic properties. Moreover, the MM/PB(GB)SA calculation indicated the tight binding and thermodynamic stability of the simulated protein-ligand complexes. Overall, this study highlights the potential of chalcones as promising candidates for targeting multidrug efflux pumps, offering a potential strategy to overcome antibiotic resistance. Further exploration and optimization of these compounds may lead to the development of effective therapeutics against multidrug-resistant bacterial infections.Communicated by Ramaswamy H. Sarma.

Pyridine Derivatives as Insecticides: Part 6. Design, Synthesis, Molecular Docking, and Insecticidal Activity of 3-(Substituted)methylthio-5,6,7,8-tetrahydroisoquinoline-4-carbonitriles Toward Aphis gossypii (Glover, 1887)

Three new series of 3-(substituted)methylthio-4-cyano-5,6,7,8-tetrahydroisoquinolines were designed and synthesized starting from readily available materials, 7-acetyl-4-cyano-1,6-dimethyl-6-hydroxy-8-(4-pyridyl, 3-pyridyl, phenyl, 4-methoxyphenyl, or 4-chlorophenyl)-5,6,7,8-tetrahydrosoquinoline-3(2H)-thiones 2a-e in high yields and very pure states. Thus, compounds 2a-e were reacted with some chloro reagents, namely, N-aryl-2-chloroacetamides 3a-f and N-(naphthalen-2-yl)-2-chloroacetamide (3g) under mild basic conditions to give the first two series of the target compounds, 3-(N-aryl)carbamoylmethylthio-5,6,7,8-tetrahydroisoquinoline-4-carbonitriles 4a-l and 5a-e, respectively. Reaction of compounds 2d,e with ethyl chloroacetate under the same conditions gave the other series, 3-ethoxycarbonyl-methylthio-5,6,7,8-tetrahydroisoquinoline-4-carbonitriles 6d,e. Structural formulas of all of the new compounds were elucidated and confirmed by elemental and spectral analyses. The insecticidal activity of all synthesized 5,6,7,8-tetrahydrosoquinolines toward the nymphs and adults of Aphis gossypii were screened. The results revealed the promising insecticidal activity of some tested compounds. Moreover, the structure-activity relationships as well as molecular docking of some representative compounds were evaluated.

Exploring Substituted Tetrazoloquinazoline: Biological Activities, Molecular Docking Analysis, and Anti-Breast Cancer MCF7/HER2 Effects

Breast cancer is a condition where breast tissue cells grow uncontrollably. Various natural and synthesized compounds, such as quinazoline, have been studied for their potential as anticancer agents. Quinazoline derivatives have shown diverse bioactivities, including antimalarial, antifungal, antimicrobial, and anticancer properties. This research aims to synthesize substituted tetrazoloquinazoline and evaluate its potential as an anticancer agent using molecular docking studies with the Molecular Operating Environment (MOE) software. Furthermore, molecular dynamic was also performed to analyze the binding stability of this protein-ligand complex. Additionally, the physicochemical and pharmacokinetic properties of quinazoline compounds were assessed using the website https://www.swissadme.ch. The cytotoxic activity of the compounds was evaluated using the MTT assay. The docking results revealed that substituted tetrazoloquinazoline exhibited a significantly different range of binding free energy compared to the positive control. Moreover, the substituted tetrazoloquinazoline compounds comply with Lipinski's Rule of Five (Ro5), indicating that they are easily absorbable and have good permeability. The cytotoxic activity of the compounds was found to have an IC50 value of >1000 ppm, classifying them as noncytotoxic. It therefore paved the way for the discovery of promising next-generation drugs against breast cancer.

In-silico screening of missense nsSNPs in Delta-opioid receptor protein and their restoring tendency on MCRT interaction; focusing on dynamic nature

Delta-opioid receptor protein (OPRD1) is one of the potential targets for treating pain. The presently available opioid agonists are known to cause unnecessary side effects. To discover a novel opioid agonist, our research group has synthesized a chimeric peptide MCRT and proved its potential activity through in vivo analysis. Non-synonymous SNPs (nsSNPs) missense mutations affect the functionality and stability of proteins leading to diseases. The current research was focused on understanding the role of MCRT in restoring the binding tendency of OPRD1 nsSNPs missense mutations on dynamic nature in comparison with Deltorphin-II and morphiceptin. The deleterious effects of nsSNPs were analyzed using various bioinformatics tools for predicting structural, functional, and oncogenic influence. The shortlisted nine nsSNPs were predicted for allergic reactions, domain changes, post-translation modification, multiple sequence alignment, secondary structure, molecular dynamic simulation (MDS), and peptide docking influence. Further, the docked complex of three shortlisted deleterious nsSNPs was analyzed using an MDS study, and the highly deleterious shortlisted nsSNP A149T was further analyzed for higher trajectory analysis. MCRT restored the binding tendency influence caused by nsSNPs on the dynamics of stability, functionality, binding affinity, secondary structure, residues connection, motion, and folding of OPRD1 protein.

Molecular interaction analysis of ferulic acid (4-hydroxy-3-methoxycinnamic acid) as main bioactive compound from palm oil waste against MCF-7 receptors: An in silico study

Ferulic acid (4-hydroxy-3-methoxycinnamic acid) is a phytochemical compound that is commonly found in conjugated forms within mono-, di-, polysaccharides and other organic compounds in cell walls of grain, fruits, and vegetables. This compound is highly abundant in the palm oil waste. The aim of the study was to predict the anticancer activity of ferulic acid against the breast cancer cell lines (MCF-7) receptors through a computational analysis. MCF-7 receptors with PDB IDs of 1R5K, 2IOG, 4IV2, 4IW6, 5DUE, 5T92, and 5U2B were selected based on the Simplified Molecular Input Line Entry System (SMILES) similarity of the native ligand. Thereafter, the protein was prepared on Chimera 1.16 and docked with ferulic acid on Autodock Vina 1.2.5. The ligand-protein complex interaction was validated by computing the root mean square fluctuation (RMSF) and radius of gyration (Rg) through molecular dynamic simulation. In addition, an absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction was performed on ferulic acid using the pkCSM platform. The molecular docking revealed that the ferulic acid could interact with all receptors as indicated by the affinity energy <-5 kcal/mol. The compound had the most optimum interaction with receptor 2IOG (affinity energy=-6.96 kcal/mol), involving hydrophobic interaction (n=12) and polar hydrogen interaction (n=4). The molecular dynamic simulation revealed that the complex had an RMSF of 1.713 Å with a fluctuation of Rg value around 1.000 Å. The ADMET properties of ferulic acid suggested that the compound is an ideal drug candidate. In conclusion, this study suggested that ferulic acid, which can be isolated from palm oil waste, has the potential to interact with MCF-7 receptors.

Computational Investigations to Identify Potent Natural Flavonoid Inhibitors of the Nonstructural Protein (NSP) 16/10 Complex Against Coronavirus

INTRODUCTION

Globally, the viral pandemic has spread rapidly, resulting in widespread infections. The coronavirus family (CoVs) is one of the various viral families capable of infecting mammals, causing diseases related to the gastrointestinal, neurological, and respiratory systems. Flavonoid compounds have been identified as potentially effective antiviral agents, specifically targeting the virus's nonstructural protein (NSP) 16/10. Flavonoids have also been shown to inhibit virus replication and viral attachment to host cells, making them a promising candidate for antiviral treatment. Further research is needed to understand the full potential of flavonoids as antiviral agents.

METHODOLOGY

This study investigated natural compounds derived from medicinal plants using in silico screening. In addition to assessing drug-likeness, pharmacokinetics, docking, molecular dynamics simulation, bioavailability assessment, and exploration of molecular targets, the screening process entailed analyses of molecular targets and bioavailability. The molecular properties and potential antiviral efficacy of these phytochemical candidates were determined by analyzing them as drug candidates. The results of the study showed that these compounds had potential antiviral activity and could be developed as therapeutic agents. Furthermore, the study showed that the compounds had good bioavailability, suggesting that they are suitable for use as therapeutic agents.

RESULT

An in silico method was used to identify flavonoid compounds for potent antiviral drug molecules against the coronavirus protein complex NSP16/10 protein. The NSP16/10 complex protein binding energy values were -6.14 for isoquercetin, -6.902 for narirutin, -6.052 for myricetin, -7.10 for hesperidin, -4.392 for silibinin, -3.997 for baicalein, -3.712 for taxifolin, and -3.321 for petunidin. Molecular dynamics simulations showed that isoquercetin, hesperidin, and narirutin flavonoids interacted with the COVID-19 virus protein complex NSP16/10 protease up to 100 nanoseconds.

Hybrid Virtual Screening Approach to Predict Novel Natural Compounds against HIV-1 CCR5

BACKGROUND

Human immunodeficiency virus (HIV) infection continues to pose a major global health challenge. HIV entry into host cells via membrane fusion mediated by the viral envelope glycoprotein gp120/gp41 is a key step in the HIV life cycle. CCR5, expressed on CD4+ T cells and macrophages, acts as a coreceptor facilitating HIV-1 entry. The CCR5 antagonist maraviroc is used to treat HIV infection. However, it can cause adverse effects and has limitations such as only inhibiting CCR5-tropic viruses. There remains a need to develop alternative CCR5 inhibitors with improved safety profiles.

PROBLEM STATEMENT

Natural products may offer advantages over synthetic inhibitors including higher bioavailability, binding affinity, effectiveness, lower toxicity, and molecular diversity. However, screening the vast chemical space of natural compounds to identify novel CCR5 inhibitors presents challenges. This study aimed to address this gap through a hybrid ligand-based pharmacophore modeling and molecular docking approach to virtually screen large natural product databases.

METHODS

A reliable pharmacophore model was developed based on 311 known CCR5 antagonists and validated against an external data set. Five natural product databases containing over 306,000 compounds were filtered based on drug-likeness rules. The validated pharmacophore model screened the databases to identify 611 hits. Key residues of the CCR5 receptor crystal structure were identified for docking. The top hits were docked, and interactions were analyzed. Molecular dynamics simulations were conducted to examine complex stability. Computational prediction evaluated pharmacokinetic properties.

RESULTS

Three compounds exhibited similar interactions and binding energies to maraviroc. MD simulations demonstrated complex stability comparable to maraviroc. One compound showed optimal predicted absorption, minimal metabolism, and a lower likelihood of interactions than maraviroc.

CONCLUSION

This computational screening workflow identified three natural compounds with promising CCR5 inhibition and favorable pharmacokinetic profiles. One compound emerged as a lead based on bioavailability potential and minimal interaction risk. These findings present opportunities for developing alternative CCR5 antagonists and warrant further experimental investigation. Overall, the hybrid virtual screening approach proved effective for mining large natural product spaces to discover novel molecular entities with drug-like properties.

Identification of potential biogenic chalcones against antibiotic resistant efflux pump (AcrB) via computational study

The cases of bacterial multidrug resistance are increasing every year and becoming a serious concern for human health. Multidrug efflux pumps are key players in the formation of antibiotic resistance, which transfer out a broad spectrum of drugs from the cell and convey resistance to the host. Efflux pumps have significantly reduced the efficacy of the previously available antibiotic armory, thereby increasing the frequency of therapeutic failures. In gram-negative bacteria, the AcrAB-TolC efflux pump is the principal transporter of the substrate and plays a major role in the formation of antibiotic resistance. In the current work, advanced computer-aided drug discovery approaches were utilized to find hit molecules from the library of biogenic chalcones against the bacterial AcrB efflux pump. The results of the performed computational studies via molecular docking, drug-likeness prediction, pharmacokinetic profiling, pharmacophore mapping, density functional theory, and molecular dynamics simulation study provided ZINC000004695648, ZINC000014762506, ZINC000014762510, ZINC000095099506, and ZINC000085510993 as stable hit molecules against the AcrB efflux pumps. Identified hits could successfully act against AcrB efflux pumps after optimization as lead molecules.Communicated by Ramaswamy H. Sarma.

Targeting key players in Alzheimer's disease: bioactive compounds from Moringa oleifera, Desmodium gangeticum, and Centella asiatica as potential therapeutics

Alzheimer's Disease (AD) is one of the critical reasons for dementia around the world, with a huge number of cases being reported every year. The breakdown of Amyloid Precursor Protein (APP) plays a crucial role in AD development. The Beta-site APP Cleaving Enzyme 1 (BACE1) is a highly significant proteolytic enzyme found to be critically involved in the APP breakdown process and generates beta-amyloid plaques in the extracellular neuronal membrane. In this study, we have used natural compounds with cognitive and neuroprotective activities from three plants, Centella asiatica, Moringa oleifera, and Desmodium gangeticum to inhibit the activity of BACE1. We have identified nine compounds out of 73 compounds filtered out from the three plants showing high affinity with the catalytic dyad region of BACE1 through molecular docking studies. Interestingly, the 200 ns molecular dynamics simulation study further confirmed the stability of the complexes formed between 9 compounds and the BACE1 protein. Furthermore, the free energy calculations also revealed these complexes possess favorable energies. Astilbin, Delphinidin 3-glucoside, and kaempferol 7-O-glucoside showed good binding affinity and structural stability when compared to other compounds and the control CNP520. Following a preliminary screening, the Astilbin compound was chosen based on the grounds of binding affinity, ADMET Properties, Hbond formation, Molecular Dynamic simulation, and MM-PBSA studies. A subsequent 1microsecond molecular dynamics simulation was conducted for the Astilbin complex. Through microsecond simulation, it was found that Astilbin alters BACE1's behavior and induces conformational rearrangements. Thus, this study opens a gateway to inhibit the activity of BACE1 protein through Astilbin thereby disclosing the possibility of managing Alzheimer's Disease.Communicated by Ramaswamy H. Sarma.

Structure-based virtual screening methods for the identification of novel phytochemical inhibitors targeting furin protease for the management of COVID-19

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, is a highly contagious respiratory disease with widespread societal impact. The symptoms range from cough, fever, and pneumonia to complications affecting various organs, including the heart, kidneys, and nervous system. Despite various ongoing efforts, no effective drug has been developed to stop the spread of the virus. Although various types of medications used to treat bacterial and viral diseases have previously been employed to treat COVID-19 patients, their side effects have also been observed. The way SARS-CoV-2 infects the human body is very specific, as its spike protein plays an important role. The S subunit of virus spike protein cleaved by human proteases, such as furin protein, is an initial and important step for its internalization into a human host. Keeping this context, we attempted to inhibit the furin using phytochemicals that could produce minimal side effects. For this, we screened 408 natural phytochemicals from various plants having antiviral properties, against furin protein, and molecular docking and dynamics simulations were performed. Based on the binding score, the top three compounds (robustaflavone, withanolide, and amentoflavone) were selected for further validation. MM/GBSA energy calculations revealed that withanolide has the lowest binding energy of -57.2 kcal/mol followed by robustaflavone and amentoflavone with a binding energy of -45.2 kcal/mol and -39.68 kcal/mol, respectively. Additionally, ADME analysis showed drug-like properties for these three lead compounds. Hence, these natural compounds robustaflavone, withanolide, and amentoflavone, may have therapeutic potential for the management of SARS-CoV-2 by targeting furin.

Exploring gene network and protein interaction analysis of neurotrophin signaling pathway in ameloblastoma

Ameloblastoma is a non-cancerous but aggressive oral tumor emerging from odontogenic epithelial tissue involved during odontogenesis. Since there is lack in unravelling the complete molecular pathogenesis of ameloblastoma, chemotherapy is less attempted and a lot of disagreement over the optimal treatment option. Hence, till date, wide surgical resection is considered to be the reliable treatment for ameloblastoma. The Neurotrophin Signaling pathway plays an important role in neuron signaling and it is closely related with the MAPK pathway, which on the other hand regulated cell differentiation, apoptosis, proliferation, plasticity and survival. Protein- Protein Interaction analysis was analysed with STRING tool using WNL value, identified that CTNNB1, HRAS, NGFR, NGFR, and SORT1 having high interacting with BDNF, NT4, p75NTR, NGF, and NT3. The results of ontology analysis revealed that Neurotrophin signaling pathway is associated with Cell surface receptor signaling pathway, regulation of cell differentiation, regulation of development process, EGFR tyrosine kinase inhibitor resistance, MAPK signaling pathway, PI3K-Akt signaling pathway and Ras signaling pathway leading to pathogenesis involving genes. Further, clustering coefficient values of proteins BDNF, NT4, p75NTR, NGF & NT3 were identified as 0.627, 0.708, 0.367, 0.644 & 0.415. The results of molecular docking studies revealed among the selected ligands Methyl-ɣ-oresellinate, N-(4-Hydroxy-phenyl)-2-phenyl-N-phenylacetyl-acetamide, Atranorin and Oresellinate exhibited high binding affinity with selected protein. The key genes involved in Neurotrophin signaling pathway leading to ameloblastoma pathogenesis is revealed, which are closely associated with cell differentiation, cell proliferation, pro-apoptosis, and pro-survival regulations. Further it can be concluded that Neurotrophin signaling pathway could be one of the promising pathway to tailor the targeted drug therapy for Ameloblastoma treatment.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00223-2.

Berberine modulates cardiovascular diseases as a multitarget-mediated alkaloid with insights into its downstream signals using in silico prospective screening approaches

Atherosclerosis is potentially correlated with several cardiac disorders that are greatly associated with cellular oxidative stress generation, inflammation, endothelial cells dysfunction, and many cardiovascular complications. Berberine is a natural isoquinoline alkaloid compound that widely modulates pathogenesis of atherosclerosis through its different curative potentials. This in silico screening study was designed to confirm the potent restorative properties of berberine chloride as a multitarget-mediated alkaloid against the CVDs and their complications through screening, identifying, visualizing, and evaluating its binding models, affinities, and interactions toward several CVDs-related targets as direct and/or indirect-mediated signals via inhibiting cellular ER stress and apoptotic signals and activating autophagy pathway. The drug-likeness properties of berberine were predicted using the computational QSAR/ADMET and Lipinski's RO5 analyses as well as in silico molecular docking simulations. The potent berberine-binding modes, residues-interaction patterns, and free energies of binding scores towards several CVDs-related targets were estimated using molecular docking tools. Furthermore, the pharmacokinetic properties and toxicological features of berberine were clearly determined. According to this in silico virtual screening study, berberine chloride could restore cardiac function and improve pathogenic features of atherosclerotic CVDs through alleviating ER stress and apoptotic signals, activating autophagy, improving insulin sensitivity, decreasing hyperglycemia and dyslipidemia, increasing intracellular RCT signaling, attenuating oxidative stress and vascular inflammation, and upregulating cellular antioxidant defenses in many cardiovascular tissues. In this in silico study, berberine chloride greatly modulated several potent CVDs-related targets, including SIGMAR1, GRP78, CASP3, BECN1, PIK3C3, SQSTM1/p62, LC3B, GLUT3, INSR, LDLR, LXRα, PPARγ, IL1β, IFNγ, iNOS, COX-2, MCP-1, IL10, GPx1, and SOD3.

6-Octadecenoic and Oleic Acid in Liquid Smoke Rice Husk Showed COVID-19 Inhibitor Properties

In recent years, liquid smoke rice husk (LSRH) has shown its therapeutic potency to diabetes, wound healing, stomatitis, and periodontitis. The phenol, 6-octadecenoic acid, oleic acid, and 9-octadecanoic acid were responsible for their therapeutic effect. The LSRH also demonstrated their potential for infectious diseases such as coronavirus disease (COVID-19). Therefore, the molecular dynamics (MDs) simulation and pharmacophore analysis was performed to analyse the binding stability of 6-octadecenoic and oleic acid. Based on MD simulation, 6-octadecenoic and oleic acids seemed to retain their interactions with Ser144 and Thr24, respectively, with hydrogen bond distance less than 2.9 Å. This interaction was stable during the simulation and has hydrophobic and hydrogen bonds/acceptors. The 6-octadecenoic acid and oleic acid were confirmed to have great potency as inhibitors for COVID-19. These compounds also showed that the existence of hydrophobic and hydrogen bonds/acceptors could increase biological activity.

Preparation, Agricultural Bioactivity Evaluation, Structure-Activity Relationships Estimation, and Molecular Docking of Some Quinazoline Compounds

Quinazoline compounds have gained significant attention in the fields of agriculture and chemistry due to their diverse activities. In this study, we focused on a series of quinazoline derivatives (4a-l). The objectives involved multiple aspects, including preparation, evaluation of their agricultural bioactivity against the maize aphid (Rhopalosiphum maidis), estimation of the structure-activity relationships (SAR), and conducting molecular docking analysis. The results of the agricultural bioactivities revealed that compound (4b) possesses the highest insecticidal activity, and the other compounds have good potential as insecticidal agents. We conducted the SARs and also molecular docking investigation to elucidate the binding modes and interactions of these compounds with target proteins relevant to the agricultural bioactivity. The docking results provided valuable information on the binding affinities and molecular interactions, aiding in the rationalization of the observed bioactivity trends. The enzyme, acetylcholinesterase (AChE), was docked with the 12 synthetic compounds (4a-l). Among these compounds, (4b), (4i), and (4e)exhibited the highest binding affinity, with docking scores (S) of -7.96, -7.83, and -7.73 kcal/mol, respectively. They were followed by compounds (4d) (S = -7.57 kcal/mol), (4c) (S = -7.53 kcal/mol), (4g) (S = -7.34 kcal/mol), (4f) (S = -7.23 kcal/mol), (4h) (S = -7.14 kcal/mol), (4k) (S = -6.61 kcal/mol), (4j) (S = -6.57 kcal/mol), (4a) (S = -6.28 kcal/mol), and finally (4l) (S = -6.01 kcal/mol). These compounds were shown to have a variety of binding interactions within the 2ACE active site, as evidenced by protein-ligand docking configurations. This study gives evidence that those compounds have AChE-inhibitory capabilities and, hence, may be used for AChE-targeting development. Also, the findings in this study highlight the potential of these compounds as agricultural agents and provide valuable insights for the design and development of some quinazoline derivatives with enhanced bioactivity for crop protection.

Synthesis, Characterization, In Silico DFT, Molecular Docking, and Dynamics Simulation Studies of Phenylhydrazono Phenoxyquinolones for Their Hypoglycemic Efficacy

A series of novel 24 phenylhydrazono phenoxyquinoline derivatives were synthesized with moderate to excellent yield and screened for their efficacy against the α-amylase enzyme through in silico studies. The structures were characterized using spectroscopic techniques such as 1HNMR, 13CNMR, and HREI-MS. Comprehensive computational studies including, drug-likeness and ADMET profiling, quantum chemical calculations, molecular docking, and molecular dynamics (MD) simulation studies, were performed. A density functional theory study of the synthesized compounds indicated a favorable reactivity profile. The synthesized novel analogues were docked against α-amylase (PDB 6OCN) enzymes to investigate the binding interactions. Based on the docking studies, one of the compounds was found to be the hit with the highest negative binding affinity for α-amylase. A MD simulation study indicated stable binding throughout the simulation.

Identification of potential JNK3 inhibitors through virtual screening, molecular docking and molecular dynamics simulation as therapeutics for Alzheimer's disease

Alzheimer's disease (AD) is a complex neurological disorder and no effective drug is available for its treatment. Numerous pathological conditions are believed to be responsible for the initiation and development of AD including c-Jun N-terminal kinases (JNKs). The JNKs are one of the enzymes from the mitogen-activated protein kinase (MAPK) family that controls the phosphorylation of various transcription factors on serine and threonine residues, and hold significant responsibilities in tasks like gene expression, cell proliferation, differentiation, and apoptosis. Since, JNK3 is primarily expressed in the brain hence its increased levels in the brain are associated with the AD pathology promoting neurofibrillary tangles, senile plaques, neuroinflammation, and nerve cell apoptosis. The current research work is focused on the development of novel JNK inhibitors as therapeutics for AD employing a structure-based virtual screening (SBVS) approach. The ZINC database (14634052 compounds) was investigated after employing pan assay interference (PAINs), drug-likeness, and diversity picking filter to distinguish molecules interacting with JNK3 by following three docking precision criteria: High Throughput Virtual Screening (HTVS), Standard Precision (SP), and Extra Precision (XP) & MMGBSA. Five lead molecules showed a better docking score in the range of -13.091 to -14.051 kcal/mol better than the reference compound (- 11.828 kcal/mol). The lead compounds displayed acceptable pharmacokinetic properties and were subjected to molecular dynamic simulations of 100 ns and binding free energy calculations. All the lead molecules showed stable RMSD and hydrogen bond interactions throughout the trajectory. The ∆GMM/PBSA_total score for the lead compounds ZINC220382956, ZINC147071339, ZINC207081127, ZINC205151456, ZINC1228819126, and CC-930 was calculated and found to be - 31.39, - 42.8, - 37.04, - 39.01, - 36.5, - 34.16 kcal/mol, respectively. Thus, it was concluded that the lead molecules identified in these studies have the potential to be explored as potent JNK3 inhibitors.

Flavonostilbenes natural hybrids from Rhamnoneuron balansae as potential antitumors targeting ALDH1A1: molecular docking, ADMET, MM-GBSA calculations and molecular dynamics studies

Several studies have linked Cancer stem cells (CSCs) to cancer resistance development to chemotherapy and radiotherapy. ALDH1A1 is a key enzyme that regulates the gene expression of CSCs and creates an immunosuppressive tumor microenvironment. It was reported that quercetin and resveratrol were among the inhibitors of ALDH1A1. In early 2022, it was reported that new 11 flavonostilbenes (rhamnoneuronal D-N) were isolated from Rhamnoneuron balansae as potential antiaging natural products. Rhamnoneuronal H (5) could be envisioned as a natural hybrid of quercetin and resveratrol. It was therefore hypothesized that 5 and its analogous isolates rhamnoneuronal D-G (1-4) and rhamnoneuronal I-N (6-11) would have potential ALDH1A1 inhibitory activity. To this end, all isolates were subjected to molecular docking, MM-GBSA, ADMET, and molecular dynamics simulations studies to assess their potential as new leads for cancer treatment targeting ALDH1A1. In silico findings revealed that natural hybrid 5 has a similar binding affinity, judged by MM-GBSA, to the ALDH1A1 active site when compared to the co-crystalized ligand (-64.71 kcal/mole and -64.12 kcal/mole, respectively). Despite having lesser affinity than that of the co-crystalized ligand, the rest of the flavonostilbenes, except 2-4, displayed better binding affinities (-37.55 kcal/mole to -58.6 kcal/mole) in comparison to either resveratrol (-34.44 kcal/mole) or quercetin (-36.48 kcal/mole). Molecular dynamic simulations showed that the natural hybrids 1, 5-11 are of satisfactory stability up to 100 ns. ADMET outcomes indicate that these hybrids displayed acceptable properties and hence could represent an ideal starting point for the development of potent ALDH1A1 inhibitors for cancer treatment.Communicated by Ramaswamy H. Sarma.

Integration of ligand and structure-based pharmacophore screening for the identification of novel natural leads against Euchromatic histone lysine methyltransferase 2 (EHMT2/G9a)

Herein, we report a blended ligand and structure-based pharmacophore screening approach to identify new natural leads against the Protein Lysine Methyltransferase 2 (EHMT2/G9a). The EHMT2/G9a has been associated with Cancer, Alzheimer's, and aging and is considered an emerging drug target having no clinically passed inhibitor. Purposefully, we developed the ligand-based pharmacophore (Pharmacophore-L) based on the common features of known inhibitors and the structure-based pharmacophore (Pharmacophore-S) based on the interaction profile of available crystal structures. The Pharmacophore-L and Pharmacophore-S were subjected to multiple tiers of validations and utilized in combination for the screening of total 741543 compounds coming from multiple databases. Additional layers of stringency were applied in the screening process to test drug-likeness (using Lipinski's rule, Veber's rule, SMARTS and ADMET filtration), to rule out any toxicity (TOPKAT analysis). The interaction profiles, stabilities, and comparative analysis against the reference were carried out by flexible docking, MD simulation, and MM-GBSA analysis, which finally led to three leads as potential inhibitors of G9a.Communicated by Ramaswamy H. Sarma.

Computer-aided drug discovery of c-Abl kinase inhibitors from plant compounds against chronic myeloid leukemia

Chronic myeloid leukemia (CML) is a hematological malignancy characterized by the neoplastic transformation of hematopoietic stem cells, driven by the Philadelphia (Ph) chromosome resulting from a translocation between chromosomes 9 and 22. This Ph chromosome harbors the breakpoint cluster region (BCR) and the Abelson (ABL) oncogene (BCR-ABL1) which have a constitutive tyrosine kinase activity. However, the tyrosine kinase activity of BCR-ABL1 have been identified as a key player in CML initiation and maintenance through c-Abl kinase. Despite advancements in tyrosine kinase inhibitors, challenges such as efficacy, safety concerns, and recurring drug resistance persist. This study aims to discover potential c-Abl kinase inhibitors from plant compounds with anti-leukemic properties, employing drug-likeness assessment, molecular docking, in silico pharmacokinetics (ADMET) screening, density function theory (DFT), and molecular dynamics simulations (MDS). Out of 58 screened compounds for drug-likeness, 44 were docked against c-Abl kinase. The top hit compound (isovitexin) and nilotinib (control drug) were subjected to rigorous analyses, including ADMET profiling, DFT evaluation, and MDS for 100 ns. Isovitexin demonstrated a notable binding affinity (-15.492 kcal/mol), closely comparable to nilotinib (-16.826 kcal/mol), showcasing a similar binding pose and superior structural stability and reactivity. While these findings suggest isovitexin as a potential c-Abl kinase inhibitor, further validation through urgent in vitro and in vivo experiments is imperative. This research holds promise for providing an alternative avenue to address existing CML treatment and management challenges.Communicated by Ramaswamy H. Sarma.

An efficient eco-friendly, simple, and green synthesis of some new spiro-N-(4-sulfamoyl-phenyl)-1,3,4-thiadiazole-2-carboxamide derivatives as potential inhibitors of SARS-CoV-2 proteases: drug-likeness, pharmacophore, molecular docking, and DFT exploration

INTRODUCTION

The coronavirus disease 2019 (COVID-19) pandemic has caused a global health crisis. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly contagious virus that can cause severe respiratory illness. There is no specific treatment for COVID-19, and the development of new drugs is urgently needed.

PROBLEM STATEMENT

The SARS-CoV-2 main protease (Mpro) enzyme is a critical viral enzyme that plays a vital role in viral replication. The inhibition of Mpro enzyme can be an effective strategy for developing new COVID-19 drugs.

METHODOLOGY

An efficient operationally simple and convenient green synthesis method had been done towards a series of novel spiro-N-(4-sulfamoylphenyl)-2-carboxamide derivatives, in ethanol at room temperature in green conditions, up to 90% yield. The molecular structures of the synthesized compounds were verified using spectroscopic methods.The title compounds were subjected to in silico analysis, including Lipinski's rule and ADMET prediction, in addition to pharmacophore modeling and molecular docking against the active site of SARS-CoV-2 target main protease (Mpro) enzyme (6LU7). Furthermore, both of the top-ranked compounds (5 and 6) and the standard Nirmatrelvir were subjected to DFT analysis.

FINDINGS

The synthesized compounds exhibited good binding affinity to SARS-CoV-2 Mpro enzyme, with binding energy scores ranging from - 7.33 kcal/mol (compound 6) and - 7.22kcal/mol (compound 5) to - 6.54 kcal/mol (compounds 8 and 9). The top-ranked compounds (5 and 6) had lower HOMO-LUMO energy difference (ΔE) than the standard drug Nirmatrelvir. This highlights the potential and relevance of charge transfer at the molecular level.

RECOMMENDATION

These findings suggest that the synthesized spiro-N-(4-sulfamoylphenyl)-2-carboxamide derivatives could be potential candidates for COVID-19 drug development. To confirm these drugs' antiviral efficacy in vivo, more research is required. With very little possibility of failure, this proven method could aid in the search for the SARS-CoV-2 pandemic's desperately needed medications.

Campesterol and dithymoquinone as a potent inhibitors of SARS cov-2 main proteases-promising drug candidates for targeting its novel variants

The sudden outbreak of the COVID-19 pandemic has currently taken approximately 2.4 million lives, with no specific medication and fast-tracked tested vaccines for prevention. These vaccines have their own adverse effects, which have severely affected the global healthcare system. The discovery of the main protease structure of coronavirus (Mpro/Clpro) has resulted in the identification of compounds having antiviral potential, especially from the herbal system. In this study, the computer-associated drug design tools were utilised to analyze the reported phytoconstituents of Nigella sativa for their antiviral activity against the main protease. Fifty-eight compounds were subjected to pharmacological parameter analysis to determine their lead likeness in comparison to the standard drugs (chloroquine and nirmatrelvir) used in the treatment of SARS-CoV-2. Nearly 31 compounds were docked against five different SARS-CoV-2 main proteases, and all compounds showed better binding affinity and inhibition constant against the proteases. However, dithymoquinone and campesterol displayed the best binding scores and hence were further subjected to dynamics and MMPBSA study for 100 ns. The stability analysis shows that dithymoquinone and campesterol show less variation in fluctuation in residues compared to standard complexes. Moreover, dithymoquinone exhibited higher binding affinity and favorable interaction followed by campesterol as compared to the standard drug. The in silico computational analysis provides a promising hit for regulating the main proteases activity.Communicated by Ramaswamy H. Sarma.

High-throughput computational screening for identification of potential hits against bacterial Acriflavine resistance protein B (AcrB) efflux pump

Antibiotic resistance is a pressing global health challenge, driven in part by the remarkable efflux capabilities of efflux pump in AcrB (Acriflavine Resistance Protein B) protein in Gram-negative bacteria. In this study, a multi-approached computational screening strategy encompassing molecular docking, In silico absorption, distribution, metabolism, excretion and toxicity (ADMET) analysis, druglikeness assessment, molecular dynamics simulations and density functional theory studies was employed to identify novel hits capable of acting against AcrB-mediated antibiotic resistance. Ligand library was acquired from the COCONUT database. Performed computational analyses unveiled four promising hit molecules (CNP0298667, CNP0399927, CNP0321542 and CNP0269513). Notably, CNP0298667 exhibited the highest negative binding affinity of -11.5 kcal/mol, indicating a possibility of strong potential to disrupt AcrB function. Importantly, all four hits met stringent druglikeness criteria and demonstrated favorable in silico ADMET profiles, underscoring their potential for further development. MD simulations over 100 ns revealed that the CNP0321542-4DX5 and CNP0269513-4DX5 complexes formed robust and stable interactions with the AcrB efflux pump. The identified hits represent a promising starting point for the design and optimization of novel therapeutics aimed at combating AcrB-mediated antibiotic resistance in Gram-negative bacteria.Communicated by Ramaswamy H. Sarma.

Revealing Propolis Potential Activity on Inhibiting Estrogen Receptor and Heat Shock Protein 90 Overexpressed in Breast Cancer by Bioinformatics Approaches

Breast cancer is the most commonly diagnosed cancer globally, with the highest incidence of breast cancer occurring in Asian countries including Indonesia. Among the types of breast cancer, the estrogen receptor (ER)-positive subtype which is prominent with estrogen receptor alpha (ERα) and heat shock protein 90 (HSP90) overexpression genes becomes the most prevalent than the others, approximately 75% of all breast cancer cases. ERα and HSP90 play a role in breast cancer activities including breast tumor growth, invasion, and metastasis mechanism. Propolis, a natural bee product, has been explored for its anticancer activity. However, there is lack of studies that evaluated the potential inhibitor from propolis compounds to the ERα and HSP90 proteins. Therefore, this article focuses on examining the correlation between ERα and HSP90's role in breast cancer and investigating the potential of 93 unique propolis compositions in inhibiting these genes in breast cancer using in silico approaches. This study revealed the positive correlation between ERα and HSP90 genes in breast cancer disease development. Furthermore, we also found novel potential bioactive compounds of propolis against breast cancer through binding with ERα and HSP90; they were 3',4',7-trihydroxyisoflavone and baicalein-7-O-β-D glucopyranoside, respectively. Further research on these compounds is needed to elucidate deeper mechanisms and activity in the real biological system to develop new breast cancer drug treatments.

Binding order of substrate and cofactor in sulfonamide monooxygenase during sulfa drug degradation: in silico studies

For decades, sulfonamide antibiotics have been used across industries such as agriculture and animal husbandry. However, the use and inadvertent misuse of these antibiotics have resulted in the advent of sulfonamide-drug-resistant strains due to antibiotic pollution. Enzymatic bioremediation of antibiotics remains a potential emerging solution to combat antibiotic pollution. Here, we propose an enzymatic model for the degradation of sulfonamides by Microbacterium sp. We have employed a multi-pronged computational strategy involving - protein structure modelling, ligand docking and molecular dynamics simulations to decipher a plausible binding order for the enzymatic degradation of sulfonamides by the bacterial sulfonamide monooxygenase, SulX. Our results enable us to predict that this degradation is achieved through the sequential binding of the antibiotic sulfonamide followed by the reduced flavin cofactor FMNH2, thereby laying the computational foundation for further advancements in enzyme-mediated degradation of the antibiotic. We also provide a list of experiments which may be performed to verify and follow-up on our in-silico studies.Communicated by Ramaswamy H. Sarma.

Computational studies and structural insights for discovery of potential natural aromatase modulators for hormone-dependent breast cancer

Introduction

The aromatase enzyme plays an important role in the progress of hormone-dependent breast cancer, especially in estrogen receptor-positive (ER+) breast cancers. In case of postmenopausal women, the aromatization of androstenedione to estrone in adipose tissue is the most important source of estrogen. Generally 60%-75% of pre- and post-menopausal women suffer from estrogen-dependent breast cancer, and thus suppressing estrogen has been recognized to be a successful therapy. Hence, to limit the stimulation of estrogen, aromatase inhibitors (AIs) are used in the second-line treatment of breast cancer.

Methods

The present computational study employed an in silico approach in the identification of natural actives targeting the aromatase enzyme from a structurally diverse set of natural products. Molecular docking, QSAR studies and pharmacophore modeling were carried out using the VLife Molecular Design Suite (version 4.6). The stability of the compounds was confirmed by molecular dynamics.

Results

From molecular docking and analysis of interactions with the amino acid residues of the binding cavity, it was found that the amino acid residues interacting with the non-steroidal inhibitors exhibited π-stacking interactions with PHE134, PHE 221, and TRP 224, while the steroidal drug exemestane lacked π-stacking interactions. QSAR studies were performed using the flavonoid compounds, in order to identify the structural functionalities needed to improve the anti-breast cancer activity. Molecular dynamics of the screened hits confirmed the stability of compounds with the target in the binding cavity. Moreover, pharmacophore modelling presented the pharmacophoric features of the selected scaffolds for aromatase inhibitory activity.

Conclusion

The results presented 23 hit compounds that can be developed as anti-breast cancer modulating agents in the near future. Additionally, anthraquinone compounds with minor structural modification can also serve to be potential aromatase inhibitors. The in silico protocol utilised can be useful in the drug discovery process for development of new leads from structurally diverse set of natural products that are comparable to the drugs used clinically in breast cancer therapy.

In silico and in vitro studies: investigating the chemical composition, DFT, molecular docking, and dynamic simulation of Satureja candidissima (Munby) Briq essential oil as a potential antibacterial agent

This study aimed to investigate the chemical composition and antibacterial properties of the essential oil (EO) derived from the aerial parts of Satureja candidissima (Munby) Briq (SC), as well as the mechanisms of interaction between SCEO chemical components and target proteins related to antibacterial activity mechanisms using a molecular docking approach, and for more accuracy molecular dynamic simulation and DFT calculations were carried out. The GC-MS technique was used to analyze the chemical composition of SCEO. The results showed that SCEO contained various chemical compounds, with pulegone being identified as the major component (53.26%). The results also indicated the presence of (+)-menthone (11.02%), borneol (4.43%), 2-cyclohexen-1-one, 3-methyl-6-(1-methylethylidene) (2.50%), and 3-octanol (2.09%). The study revealed that the SCEO displayed antibacterial activity against all tested gram-positive bacteria. To further understand the mechanism behind its antibacterial activity, in silico molecular docking studies were performed. The results indicated that the antibacterial effect of SCEO compounds could be due to the combination with enoyl-[acyl-carrier-protein] reductase [NADPH] FabI (PDB ID: 4ALL) in a variety of ways. The molecular dynamics simulation analysis yielded favorable outcomes for the docked complex involving 1H-cycloprop[e]azulen-7-ol, decahydro-1,1,7-trimethyl-4-methylene, and 1,4,7-tetramethyldecahydro-1H-cyclopropa[e]azulen-4-ol with enoyl-[acyl-carrier-protein] reductase [NADPH]. Geometry optimization, coupled with Density Functional Theory (DFT), can be employed to assess the importance of quantum chemical descriptors in elucidating potential antibacterial activity. Quantum descriptors were computed based on EHOMO and ELUMO. The results of this study provide important insights into the potential use of Satureja candidissima (Munby) Briq EO as antibacterial agent.Communicated by Ramaswamy H. Sarma.

Karanjin, A Promising Bioactive Compound Possessing Anti-cancer Activity against Experimental Model of Non-small Cell Lung Cancer Cells

AIMS

The aim of this study is to isolate the Millettia pinnata (Karanj) leaf extract for pure compound with anticancer properties and to study the molecular target of the isolates in non-small cell lung cancer cell lines.

BACKGROUND

In our earlier research Millettia pinnata leaf extract has demonstrated potential anticancer activities. Thus, in pursuit of the bioactive compounds, the most potential active extract from our previous study was purified. Furthermore, the anticancer properties of the isolated compound karanjin was studied and aimed for apoptosis and restraining growth.

METHODS

A novel method was developed through column chromatography for isolation and purification of the compound karanjin from leaf chloroform extract. The purified component was then characterised using FTIR, mass spectrometry, and NMR. An MTT-based cytotoxicity assay was used to analyse cell cytotoxicity, whereas fluorescence staining was used for apoptosis and reactive oxygen species inhibition quantification. Furthermore, the real-time PCR assay was used to determine the molecular mechanism of action in cells causing cytotoxicity induced by karanjin dosing.

RESULTS

The anticancer activity of karanjin in A549 cell line exhibited prominent activity revealing IC50 value of 4.85 μM. Conferring the predicted molecular pathway study, karanjin restrains the proliferation of cancer cells through apoptosis, which is controlled by extrinsic pathway proteins FAS/FADD/Caspases 8/3/9. Downregulation of KRAS and dependent gene expression also stopped cell proliferation.

CONCLUSION

Karanjin has been identified as a compound with potential effect in non-small cell lung cancer cells. Molecular mechanism for apoptosis and inhibition of reactive oxygen species induced through H2O2 were observed, concluding karanjin have medicinal and antioxidant properties.

Structural Dynamics of P-Rex1 Complexed with Natural Leads Establishes the Protein as an Attractive Target for Therapeutics to Suppress Cancer Metastasis

Phosphatidylinositol 3,4,5-trisphosphate- (PIP3-) dependent Rac exchanger 1 (P-Rex1) functions as Rho guanine nucleotide exchange factor and is activated by synergistic activity of Gβγ and PIP3 of the heterotrimeric G protein. P-Rex1 activates Rac GTPases for regulating cell invasion and migration and promotes metastasis in several human cancers including breast, prostate, and skin cancer. The protein is a promising therapeutic target because of its multifunction roles in human cancers. Herein, the present study attempts to identify selective P-Rex1 natural inhibitors by targeting PIP3-binding pocket using large-size multiple natural molecule libraries. Each library was filtered subsequently in FAF-Drugs4 based on Lipinski's rule of five (RO5), toxicity, and filter pan assay interference compounds (PAINS). The output hits were virtually screened at the PIP3-binding pocket through PyRx AutoDock Vina and cross-checked by GOLD. The best binders at the PIP3-binding pocket were prioritized using a comparative analysis of the docking scores. Top-ranked two compounds with high GOLD fitness score (>80) and lowest AutoDock binding energy (< -12.7 kcal/mol) were complexed and deciphered for molecular dynamics along with control-P-Rex1 complex to validate compound binding conformation and disclosed binding interaction pattern. Both the systems were seen in good equilibrium, and along the simulation time, the compounds are in strong contact with the P-Rex1 PIP3-binding site. Hydrogen bonding analysis towards simulation end identified the formation of 16 and 22 short- and long-distance hydrogen bonds with different percent of occupancy to the PIP3 residues for compound I and compound 2, respectively. Radial distribution function (RDF) analysis of the key hydrogen bonds between the compound and the PIP3 residues demonstrated a strong affinity of the compounds to the mentioned PIP3 pocket. Additionally, MMGB/PBSA energies were performed that confirmed the dominance of Van der Waals energy in complex formation along with favorable contribution from hydrogen bonding. These findings were also cross-validated by a more robust WaterSwap binding energy predictor, and the results are in good agreement with a strong binding affinity of the compounds for the protein. Lastly, the key contribution of residues in interaction with the compounds was understood by binding free energy decomposition and alanine scanning methods. In short, the results of this study suggest that P-Rex1 is a good druggable target to suppress cancer metastasis; therefore, the screened druglike molecules of this study need in vitro and in vivo anti-P-Rex1 validation and may serve as potent leads to fight cancer.

Co-production of 2,5-dihydroxymethylfuran and furfuralcohol from sugarcane bagasse via chemobiocatalytic approach in a sustainable system

2,5-Dihydroxymethylfuran and furfuryl alcohol serve as versatile building-blocks in pharmaceuticals, polymers, and value-added intermediates. To develop an efficient and sustainable method for their production from biomass, a combined approach using deep eutectic solvent Citric acid:Betaine (CTA:BT) for bagasse catalysis and recombinant E. coli SCFD23 for bioreduction of bagasse-derived 5-hydroxymethylfurfural and furfural was devised. Bagasse was effectively transformed into 5-hydroxymethylfurfural (48 mM) and furfural (14 mM) in CTA:BT (8 wt%)-water at 170 °C for 30 min. Bioreduction of 5-hydroxymethylfurfural and furfural by SCFD23 cell co-expressing formate dehydrogenase and NAD(P)H-dependent aldehyde reductase (SsCR) yielded 2,5-dihydroxymethylfuran (90.0 % yield) and furfuryl alcohol (99.0 % yield) in 6 h, using biomass-derived formic acid, xylose and glucose as co-substrates. Molecular docking confirmed the stable binding and reductase activity of SsCR with the biomass-derived 5-hydroxymethylfurfural and furfural. An efficient and eco-friendly chemobiological approach was applied for co-production of 2,5-dihydroxymethylfuran and furfuryl alcohol from biomass in one-pot two-step reaction.

Novel pseudonucleosides and sulfamoyl-oxazolidinone β-D-glucosamine derivative as anti-COVID-19: design, synthesis, and in silico study

New pseudonucleosides containing cyclic sulfamide moiety and sulfamoyl β-D-glucosamine derivative are described. These pseudonucleosides are synthesized in good yields starting from chlorosulfonyl isocyanate and β-D-glucosamine hydrochloride in five steps; (protection, acetylation, removal of the Boc group, sulfamoylation, and cyclization). Further, novel glycosylated sulfamoyloxazolidin-2-one is prepared in three steps; carbamoylation, sulfamoylation, and intramolecular cyclization. The structures of the synthesized compounds were confirmed by usual spectroscopic and spectrometric methods NMR, IR, MS, and EA. Interesting molecular docking of the prepared pseudonucleosides and (Beclabuvir, Remdesivir) drugs with SARS-CoV-2/Mpro (PDB:5R80) was conducted using the same parameters for a fair comparison. A low binding affinity of the synthesized compounds compared to the Beclabuvir and other analysis showed that pseudonucleosides have the ability to inhibit SARS-CoV-2. After the motivating results of molecular docking study, the complex between the SARS-CoV-2 Mpro and compound 7 was subjected to 100 ns molecular dynamics (MD) simulation using Desmond module of Schrodinger suite, during which the receptor-ligand complex showed substantial stability after 10 ns of MD simulation. Also, we studied the prediction of absorption, distribution, properties of metabolism, excretion, and toxicity (ADMET) of the synthesized compounds.Communicated by Ramaswamy H. Sarma.

Rhamnetin, a nutraceutical flavonoid arrests cell cycle progression of human ovarian cancer (SKOV3) cells by inhibiting the histone deacetylase 2 protein

Overexpression of HDAC 2 promotes cell proliferation in ovarian cancer. HDAC 2 is involved in chromatin remodeling, transcriptional repression, and the formation of condensed chromatin structures. Targeting HDAC 2 presents a promising therapeutic approach for correcting cancer-associated epigenetic abnormalities. Consequently, HDAC 2 inhibitors have evolved as an attractive class of anti-cancer agents. This work intended to investigate the anti-cancer abilities and underlying molecular mechanisms of Rhamnetin in human epithelial ovarian carcinoma cells (SKOV3), which remain largely unexplored. We employed various in vitro methods, including MTT, apoptosis study, cell cycle analysis, fluorescence microscopy imaging, and in vitro enzymatic HDAC 2 protein inhibition, to examine the chemotherapeutic sensitivity of Rhamnetin in SKOV3 cells. Additionally, we conducted in silico studies using molecular docking, MD simulation, MM-GBSA, DFT, and pharmacokinetic analysis to investigate the binding interaction mechanism within Rhamnetin and HDAC 2, alongside the compound's prospective as a lead candidate. The in vitro assay confirmed the cytotoxic effects of Rhamnetin on SKOV3 cells, through its inhibition of HDAC 2 activity. Rhamnetin, a nutraceutical flavonoid, halted at the G1 phase of the cell cycle and triggered apoptosis in SKOV3 cells. Furthermore, computational studies provided additional evidence of its stable binding to the HDAC 2 protein's binding site cavity. Based on our findings, we conclude that Rhamnetin effectively promotes apoptosis and mitigates the proliferation of SKOV3 cells through HDAC 2 inhibition. These results highlight Rhamnetin as a potential lead compound, opening a new therapeutic strategy for human epithelial ovarian cancer.Communicated by Ramaswamy H. Sarma.

Structure-based design, optimization of lead, synthesis, and biological evaluation of compounds active against Trypanosoma cruzi

Chagas' disease affects approximately eight million people throughout the world, especially the poorest individuals. The protozoan that causes this disease-Trypanosoma cruzi-has the enzyme cruzipain, which is the main therapeutic target. As no available medications have satisfactory effectiveness and safety, it is of fundamental importance to design and synthesize novel analogues that are more active and selective. In the present study, molecular docking and the in silico prediction of ADMET properties were used as strategies to optimize the trypanocidal activity of the pyrimidine compound ZN3F based on interactions with the target site in cruzipain. From the computational results, eight 4-amino-5-carbonitrile-pyrimidine analogues were proposed, synthesized (5a-f and 7g-h) and, tested in vitro on the trypomastigote form of the Tulahuen strain of T. cruzi. The in silico study showed that the designed analogues bond favorably to important amino acid residues of the active site in cruzipain. An in vitro evaluation of cytotoxicity was performed on L929 mammal cell lines. All derivatives inhibited the Tulahuen strain of T. cruzi and also exhibited lower toxicity to L929 cells. The 5e product, in particular, proved to be a potent, selective (IC50  = 2.79 ± 0.00 μM, selectivity index = 31.3) inhibitor of T. cruzi. The present results indicated the effectiveness of drugs based on the structure of the receptor, revealing the potential trypanocidal of pyrimidines. This study also provides information on molecular aspects for the inhibition of cruzipain.

Combined In Silico and In Vitro Study to Reveal the Structural Insights and Nucleotide-Binding Ability of the Transcriptional Regulator PehR from the Phytopathogen Ralstonia solanacearum

The transcriptional regulator PehR regulates the synthesis of the extracellular plant cell wall-degrading enzyme polygalacturonase, which is essential in the bacterial wilt of plants caused by one of the most devastating plant phytopathogens, Ralstonia solanacearum. The bacterium has a wide global distribution infecting many different plant species, resulting in massive agricultural and economic losses. Because the PehR molecular structure has not yet been determined and the structural consequences of PehR on ligand binding have not been thoroughly investigated, we have used an in silico approach combined with in vitro experiments for the first time to characterize the PehR regulator from a local isolate (Tezpur, Assam, India) of the phytopathogenic bacterium R. solanacearum F1C1. In this study, an in silico approach was employed to model the 3D structure of the PehR regulator, followed by the binding analysis of different ligands against this regulatory protein. Molecular docking studies suggest that ATP has the highest binding affinity for the PehR regulator. By using molecular dynamics (MD) simulation analysis, involving root-mean-square deviation, root-mean-square fluctuations, hydrogen bonding, radius of gyration, solvent-accessible surface area, and principal component analysis, it was possible to confirm the sudden conformational changes of the PehR regulator caused by the presence of ATP. We used an in vitro approach to further validate the formation of the PehR-ATP complex. In this approach, recombinant DNA technology was used to clone, express, and purify the gene encoding the PehR regulator from R. solanacearum F1C1. Purified PehR was used in ATP-binding experiments using fluorescence spectroscopy and Fourier transform infrared spectroscopy, the outcomes of which showed a potent binding to ATP. The putative PehR-ATP-binding analysis revealed the importance of the amino acids Lys190, Glu191, Arg192, Arg375, and Asp378 for the ATP-binding process, but further study is required to confirm this. It will be simpler to comprehend the catalytic mechanisms of a crucial PehR regulator process in R. solanacearum with the aid of the ATP-binding process hints provided by these structural biology applications.

Identification of hits as anti-obesity agents against human pancreatic lipase via docking, drug-likeness, in-silico ADME(T), pharmacophore, DFT, molecular dynamics, and MM/PB(GB)SA analysis

Obesity, characterized by excessive fat accumulation, is a major health concern. Inhibition of human pancreatic lipase, an enzyme involved in fat digestion, offers a potential strategy for weight loss and obesity treatment. This study aimed to identify polyphenols capable of forming stable complexes with human pancreatic lipase to block its activity. Molecular docking, density functional theory (DFT), molecular dynamics (MD) simulations, and MMPBGBSA calculations were employed to evaluate ligand binding, stability, and energy profiles. Pharmacophore modeling was also performed to identify key structural features for effective inhibition. Virtual screening identified ZINC000015120539, ZINC000000899200, ZINC000001531702, and ZINC000013340267 as potential candidates, exhibiting favorable binding and stable interactions over 100 ns MD simulations. These findings provide insights into the inhibitory potential of selected polyphenols on human pancreatic lipase and support further experimental investigations for obesity treatment.Communicated by Ramaswamy H. Sarma.

Expediting the drug discovery for ideal leads against SARS-CoV-2 via molecular docking of repurposed drugs

SARS-CoV-2, the novel coronavirus spreading worldwide urges the need to repurpose drugs that can quickly enter clinical trials to combat the on-going global pandemic. A cluster of proteins are encoded for by the viral genome, each assuming a critical role in pathogen endurance inside the host. To handle the adverse circumstances, robust virtual strategies such as repurposing are coming to the fore due to being economical, efficient and rapid. Five FDA approved repurposed drugs proposed to act as inhibitors by targeting SARS-CoV-2 were used for initial evaluation via molecular docking. Moreover, a comparative analysis of the selected SARS-CoV-2 proteins against five ligands (Clemizole hydrochloride, Exemestane, Nafamostat, Pregnenolone and Umifenovir) was designed. In this regard, non-structural proteins (nsp3, nsp5, nsp10, nsp12 and nsp15), structural proteins (Spike, Nucleocapsid protein) and accessory proteins (ORF 3a, ORF 7a and ORF 9 b) were selected. Here, we aim to expedite the search for a potential drug from the five FDA approved repurposing drugs already in use for treatment of multiple diseases. Based on docking analysis, Umifenovir and Pregnenolone are suggested to show potential inhibitory effects against most of the SARS-CoV-2 proteins. These drugs are noteworthy since they exhibit high binding towards target proteins and should be used as lead compounds towards in vitro and in vivo studies.Communicated by Ramaswamy H. Sarma.

Potency of Xanthone Derivatives from Garcinia mangostana L. for COVID-19 Treatment through Angiotensin-Converting Enzyme 2 and Main Protease Blockade: A Computational Study

ACE2 and Mpro in the pathology of SARS-CoV-2 show great potential in developing COVID-19 drugs as therapeutic targets, due to their roles as the "gate" of viral entry and viral reproduction. Of the many potential compounds for ACE2 and Mpro inhibition, α-mangostin is a promising candidate. Unfortunately, the potential of α-mangostin as a secondary metabolite with the anti-SARS-CoV-2 activity is hindered due to its low solubility in water. Other xanthone isolates, which also possess the xanthone core structure like α-mangostin, are predicted to be potential alternatives to α-mangostin in COVID-19 treatment, addressing the low drug-likeness of α-mangostin. This study aims to assess the potential of xanthone derivative compounds in the pericarp of mangosteen (Garcinia mangostana L.) through computational study. The study was conducted through screening activity using molecular docking study, drug-likeness prediction using Lipinski's rule of five filtration, pharmacokinetic and toxicity prediction to evaluate the safety profile, and molecular dynamic study to evaluate the stability of formed interactions. The research results showed that there were 11 compounds with high potential to inhibit ACE2 and 12 compounds to inhibit Mpro. However, only garcinone B, in addition to being indicated as active, also possesses a drug-likeness, pharmacokinetic, and toxicity profile that was suitable. The molecular dynamic study exhibited proper stability interaction between garcinone B with ACE2 and Mpro. Therefore, garcinone B, as a xanthone derivative isolate compound, has promising potential for further study as a COVID-19 treatment as an ACE2 and Mpro inhibitor.

In Silico Design of New Dual Inhibitors of SARS-CoV-2 MPRO through Ligand- and Structure-Based Methods

The viral main protease is one of the most attractive targets among all key enzymes involved in the life cycle of SARS-CoV-2. Considering its mechanism of action, both the catalytic and dimerization regions could represent crucial sites for modulating its activity. Dual-binding the SARS-CoV-2 main protease inhibitors could arrest the replication process of the virus by simultaneously preventing dimerization and proteolytic activity. To this aim, in the present work, we identified two series' of small molecules with a significant affinity for SARS-CoV-2 M^PRO, by a hybrid virtual screening protocol, combining ligand- and structure-based approaches with multivariate statistical analysis. The Biotarget Predictor Tool was used to filter a large in-house structural database and select a set of benzo[b]thiophene and benzo[b]furan derivatives. ADME properties were investigated, and induced fit docking studies were performed to confirm the DRUDIT prediction. Principal component analysis and docking protocol at the SARS-CoV-2 M^PRO dimerization site enable the identification of compounds 1b,c,i,l and 2i,l as promising drug molecules, showing favorable dual binding site affinity on SARS-CoV-2 M^PRO.

Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies

The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.

Tenofovir antiviral drug solubility enhancement with β-cyclodextrin inclusion complex and in silico study of potential inhibitor against SARS-CoV-2 main protease (Mpro)

Tenofovir (TFR) is an antiviral drug commonly used to fight against viral diseases infection due to its good potency and high genetic barrier to drug resistance. In physiological conditions, TFR is less water soluble, more unstable, and less permeable, limiting its effective therapeutic applications. In addition to their use in treating the Coronavirus disease 2019 (COVID-19), cyclodextrins (CDs) are also being used as a molecule to develop therapies for other diseases due to its enhance solubility and stability. This study is designed to synthesize and characterization of β-CD:TFR inclusion complex and its interaction against SARS-CoV-2 (M^Pro) protein (PDB ID;7cam). Several techniques were used to characterize the prepared β-CD:TFR inclusion complex, including UV-Visible, FT-IR, XRD, SEM, TGA, and DSC, which provided appropriate evidence to confirm the formation. A 1:1 stoichiometry was determined for β-CD:TFR inclusion complex in aqueous medium from UV-Visible absorption spectra by using the Benesi-Hildebrand method. Phase solubility studies proposed that β-CD enhanced the excellent solubility of TFR and the stability constant was obtained at 863 ± 32 M^-1. Moreover, the molecular docking confirmed the experimental results demonstrated the most desirable mode of TFR encapsulated into the β-CD nanocavity via hydrophobic interactions and possible hydrogen bonds. Moreover, TFR was validated in the β-CD:TFR inclusion complex as potential inhibitors against SARS-CoV-2 main protease (M^pro) receptors by using in silico methods. The enhanced solubility, stability, and antiviral activity against SARS-CoV-2 (M^Pro) suggest that β-CD:TFR inclusion complexes can be further used as feasible water-insoluble antiviral drug carriers in viral disease infection.

Network pharmacology of apigeniflavan: a novel bioactive compound of Trema orientalis Linn. in the treatment of pancreatic cancer through bioinformatics approaches

Pancreatic cancer is the seventh most prevalent cause of mortality globally. Since time immemorial, plant-derived products have been in use as therapeutic agents due to the existence of biologically active molecules called secondary metabolites. Flavonoids obtained from plants participate in cell cycle arrest, induce autophagy and apoptosis, and decrease oxidative stress in pancreatic cancer. The present study involves network pharmacology-based study of the methanolic leaf extract of Trema orientalis (MLETO) Linn. From the high-resolution mass spectrometry (HRMS) analysis, 21 nucleated flavonoids were screened out, of which only apigeniflavan was selected for further studies because it followed Lipinski's rule and showed no toxicity. The pharmacokinetics and physiochemical characteristics of apigeniflavan were performed using the online web servers pkCSM, Swiss ADME, and ProTox-II. This is the first in silico study to report the efficiency of apigeniflavan in pancreatic cancer treatment. The targets of apigeniflavan were fetched from SwissTargetPrediction database. The targets of pancreatic cancer were retrieved from DisGeNET and GeneCards. The protein-protein interaction of the common genes using Cytoscape yielded the top five hub genes: KDR, VEGFA, AKT1, SRC, and ESR1. Upon molecular docking, the lowest binding energies corresponded to best docking score which indicated the highest protein-ligand affinity. Kyoto Encyclopaedia of Genes and Genomes (KEGG) database was employed to see the involvement of hub genes in pathways related to pancreatic cancer. The following, pancreatic cancer pathway, MAPK, VEGF, PI3K-Akt, and ErbB signaling pathways, were found to be significant. Our results indicate the involvement of the hub genes in tumor growth, invasion and proliferation in the above-mentioned pathways, and therefore necessitating their downregulation. Moreover, apigeniflavan can flourish as a promising drug for the treatment of pancreatic cancer in future.

Application of molecular dynamics simulation for exploring the roles of plant biomolecules in promoting environmental health

Understanding the dynamic changes of plant biomolecules is vital for exploring their mechanisms in the environment. Molecular dynamics (MD) simulation has been widely used to study structural evolution and corresponding properties of plant biomolecules at the microscopic scale. Here, this review (i) outlines structural properties of plant biomolecules, and the crucial role of MD simulation in advancing studies of the biomolecules; (ii) describes the development of MD simulation in plant biomolecules, determinants of simulation, and analysis parameters; (iii) introduces the applications of MD simulation in plant biomolecules, including the response of the biomolecules to multiple stresses, their roles in corrosive environments, and their contributions in improving environmental health; (iv) reviews techniques integrated with MD simulation, such as molecular biology, quantum mechanics, molecular docking, and machine learning modeling, which bridge gaps in MD simulation. Finally, we make suggestions on determination of force field types, investigation of plant biomolecule mechanisms, and use of MD simulation in combination with other techniques. This review provides comprehensive summaries of the mechanisms of plant biomolecules in the environment revealed by MD simulation and validates it as an applicable tool for bridging gaps between macroscopic and microscopic behavior, providing insights into the wide application of MD simulation in plant biomolecules.

Structure-Activity Relationship of Methyl 4-Aminobenzoate Derivatives as Being Drug Candidate Targeting Glutathione Related Enzymes: in Vitro and in Silico Approaches

A thiol compound, glutathione, is essential for healthy cell defence against xenobiotics and oxidative stress. Glutathione reductase (GR) and glutathione S-transferase (GST) are two glutathione-related enzymes that function in the antioxidant and the detoxification systems. In this study, potential inhibitory effects of methyl 4-aminobenzoate derivatives on GR and GST were examined in vitro. GR and GST were isolated from human erythrocytes with 7.63 EU/mg protein and 5.66 EU/mg protein specific activity, respectively. It was found that compound 1 (methyl 4-amino-3-bromo-5-fluorobenzoate with K_i value of 0.325±0.012 μM) and compound 5 (methyl 4-amino-2-nitrobenzoate with K_i value of 92.41±22.26 μM) inhibited GR and GST stronger than other derivatives. Furthermore, a computer-aided method was used to predict the binding affinities of derivatives, ADME characteristics, and toxicities. Derivatives 4 (methyl 4-amino-2-bromobenzoate) and 6 (methyl 4-amino-2-chlorobenzoate) were estimated to have the lowest binding energies into GR and GST receptors, respectively according to results of in silico studies.

Exploring biogenic chalcones as DprE1 inhibitors for antitubercular activity via in silico approach

Cases of drug-resistant tuberculosis (TB) have increased worldwide in the last few years, and it is a major threat to global TB control strategies and the human population. Mycobacterium tuberculosis is a common causative agent responsible for increasing cases of TB and as reported by WHO, approximately, 1.5 million death occurred from TB in 2020. Identification of new therapies against drug-resistant TB is an urgent need to be considered primarily. The current investigation aims to find the potential biogenic chalcone against the potential targets of drug-resistant TB via in silico approach. The ligand library of biogenic chalcones was screened against DprE1. Results of molecular docking and in silico ADMET prediction revealed that ZINC000005158606 has lead-like properties against the targeted protein. Pharmacophore modeling was done to identify the pharmacophoric features and their geometric distance present in ZINC000005158606. The binding stability study performed using molecular dynamics (MD) simulation of the DprE1-ZINC000005158606 complex revealed the conformational stability of the complex system over 100 ns with minimum deviation. Further, the in silico anti-TB sensitivity of ZINC000005158606 was found to be higher as compared to the standards against Mycobacterium tuberculosis. The overall in silico investigation indicated the potential of identified hit to act as a lead molecule against Mycobacterium tuberculosis.

Mechanism exploration and prognosis study of Astragali Radix-Spreading hedyotis herb for the treatment of lung adenocarcinoma based on bioinformatics approaches and molecular dynamics simulation

Background: Herb pair of Astragali Radix (AR) and Spreading Hedyotis Herb (SH) has been frequently prescribed in clinical for the treatment of lung cancer owing to its favorable efficacy. Yet, the mechanism under the therapeutic effects remained unveiled, which has limited its clinical applications, and new drug development for lung cancer.

Methods: The bioactive ingredients of AR and SH were retrieved from the Traditional Chinese Medicine System Pharmacology Database, with the targets of obtained components predicted by Swiss Target Prediction. Genes related to lung adenocarcinoma (LUAD) were acquired from GeneCards, OMIM and CTD databases, with the hub genes of LUAD screened by CTD database. The intersected targets of LUAD and AR-SH were obtained by Venn, with David Database employed to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Survival analysis of the hub genes of LUAD was carried out using TCGA-LUAD dataset. Molecular docking of core proteins and active ingredients was performed by Auto-Dock Vina software, followed by molecular dynamics simulations of protein-ligand complexes with well-docked conformations.

Results: 29 active ingredients were screened out with 422 corresponding targets predicted. It is revealed that AR-SH can act on various targets such as EGFR, MAPK1, and KARS by ursolic acid (UA), Astragaloside IV(ASIV), and Isomucronulatol 7,2′-di-O-glucoside (IDOG) to alleviate the symptoms of LUAD. Biological processes involved are protein phosphorylation, negative regulation of apoptotic process, and pathways involved are endocrine resistance, EGFR tyrosine kinase inhibitor resistance, PI3K-Akt, and HIF-1 pathway. Molecular docking analysis indicated that the binding energy of most of the screened active ingredients to proteins encoded by core genes was less than −5.6 kcal/mol, with some active ingredients showing even lower binding energy to EGFR than Gefitinib. Three ligand-receptor complexes including EGFR-UA, MAPK1-ASIV, and KRAS-IDOG were found to bind relatively stable by molecular dynamics simulation, which was consistent with the results of molecule docking.

Conclusion: We suggested that the herb pair of AR-SH can act on targets like EGFR, MAPK1 and KRAS by UA, ASIV and IDOG, to play a vital role in the treatment and the enhancement of prognosis of LUAD.

Identification of novel human nicotinamide N-methyltransferase inhibitors: a structure-based pharmacophore modeling and molecular dynamics approach

Human nicotinamide N-methyltransferase (hNNMT) is a cytosolic enzyme associated in the phase-II metabolism, belonging to the S-adenosyl-L-methionine (SAM)-dependent methyltransferases family. Overexpression of hNNMT was observed in diseases such as metabolic disorders and different types of cancers, which suggest NNMT as a prospective therapeutic target. In this study we propose a structure-based pharmacophore model to understand the structural features responsible for the pharmacological activity. The generated model was validated using the ROC curve (AUC), goodness of hit score (GH), specificity, sensitivity and enrichment factor (EF). The pharmacophore was employed to retrieve active molecules from the ZINC database, followed by virtual-screening and molecular docking. Six molecules with the best pharmfit score, binding energy and ADMET properties were identified in this study. A 150 ns molecular dynamics simulation was performed on the selected molecules complexed with hNNMT protein to validate the results. The molecules ZINC35464499, ZINC13311192, ZINC31159282, ZINC14650833, ZINC14819515 and ZINC00303881 were identified, which could be act as the potential hNNMT inhibitors and can also be used as direct hits for developing novel hNNMT antagonists.Communicated by Ramaswamy H. Sarma.