A detailed TD-DFT and intermolecular interaction study of vitamin K in soluble, poorly soluble and insoluble solvents, as well as an ADME and molecular docking study

Vitamin K is one of the most important fat-soluble vitamins and while there are two main types of vitamin K in nature, known as K1 (phylloquinone) and K2 (menaquinones), there is also a synthetic type of vitamin K known as K3 (menadione). Recent studies have shown that it is crucial to know the non-covalent interactions, ADME and molecular docking of molecules in different solvent media. Therefore, we have performed some quantum chemical calculations, ADME and intra-and intermolecular interaction calculations of a number of K1, K2 and K3 such as K1-water (K1 + W), K1-methanol (K1 + M), K1-triacetin (K1 + T), K2-water (K2 + W), K2-methanol (K2 + M), K2-triacetin (K2 + T), K3-water (K3 + W), K3-methanol (K3 + M), K3-triacetin (K3 + T) performed by Density Functional Theory (DFT) and Multiwfn: A multifunctional wavefunction analyzer. Molecular structures, HOMO-LUMO energies, MEP and electronic properties have been calculated and described using DFT at the level of B3LYP/6-311G (d,p) level. The nature of the molecular interactions between vitamin K and solvents such as water, methanol and triacetin were also investigated using topological analyses such as atoms in molecule (AIM), non-covalent interaction index (NCI), reduced density gradient (RDG), Localized orbital locator (LOL) and electron localization function (ELF). In addition, FMO for electronic transitions, MEP for electrophilic and nucleophilic attack, ADME to investigate how a chemical is processed by a living organism, and Fukui functions to determine electron density are explained. Finally, molecular docking was used to determine the biological activity of the vitamin K.

Application of in-silico drug discovery techniques to discover a novel hit for target-specific inhibition of SARS-CoV-2 Mpro's revealed allosteric binding with MAO-B receptor: A theoretical study to find a cure for post-covid neurological disorder

Several studies have revealed that SARS-CoV-2 damages brain function and produces significant neurological disability. The SARS-CoV-2 coronavirus, which causes COVID-19, may infect the heart, kidneys, and brain. Recent research suggests that monoamine oxidase B (MAO-B) may be involved in metabolomics variations in delirium-prone individuals and severe SARS-CoV-2 infection. In light of this situation, we have employed a variety of computational to develop suitable QSAR model using PyDescriptor and genetic algorithm-multilinear regression (GA-MLR) models (R2 = 0.800-793, Q2LOO = 0.734-0.727, and so on) on the data set of 106 molecules whose anti-SARS-CoV-2 activity was empirically determined. QSAR models generated follow OECD standards and are predictive. QSAR model descriptors were also observed in x-ray-resolved structures. After developing a QSAR model, we did a QSAR-based virtual screening on an in-house database of 200 compounds and found a potential hit molecule. The new hit's docking score (-8.208 kcal/mol) and PIC50 (7.85 M) demonstrated a significant affinity for SARS-CoV-2's main protease. Based on post-covid neurodegenerative episodes in Alzheimer's and Parkinson's-like disorders and MAO-B's role in neurodegeneration, the initially disclosed hit for the SARS-CoV-2 main protease was repurposed against the MAO-B receptor using receptor-based molecular docking, which yielded a docking score of -12.0 kcal/mol. This shows that the compound that inhibits SARS-CoV-2's primary protease may bind allosterically to the MAO-B receptor. We then did molecular dynamic simulations and MMGBSA tests to confirm molecular docking analyses and quantify binding free energy. The drug-receptor complex was stable during the 150-ns MD simulation. The first computational effort to show in-silico inhibition of SARS-CoV-2 Mpro and allosteric interaction of novel inhibitors with MAO-B in post-covid neurodegenerative symptoms and other disorders. The current study seeks a novel compound that inhibits SAR's COV-2 Mpro and perhaps binds MAO-B allosterically. Thus, this study will enable scientists design a new SARS-CoV-2 Mpro that inhibits the MAO-B receptor to treat post-covid neurological illness.

Development of novel ligands against SARS-CoV-2 Mpro enzyme: an in silico and in vitro Study

BACKGROUND

Despite tremendous efforts made by scientific community during the outbreak of COVID-19 pandemic, this disease still remains as a public health concern. Although different types of vaccines were globally used to reduce the mortality, emergence of new variants of SARS-CoV-2 is a challenging issue in COVID-19 pharmacotherapy. In this context, target therapy of SARS-CoV-2 by small ligands is a promising strategy.

METHODS

In this investigation, we applied ligand-based virtual screening for finding novel molecules based on nirmatrelvir structure. Various criteria including drug-likeness, ADME, and toxicity properties were applied for filtering the compounds. The selected candidate molecules were subjected to molecular docking and dynamics simulation for predicting the binding mode and binding free energy, respectively. Then the molecules were experimentally evaluated in terms of antiviral activity against SARS-CoV-2 and toxicity assessment.

RESULTS

The results demonstrated that the identified compounds showed inhibitory activity towards SARS-CoV-2 Mpro .

CONCLUSION

In summary, the introduced compounds may provide novel scaffold for further structural modification and optimization with improved anti SARS-CoV-2 Mpro activity.

Active site determination of novel plant versatile peroxidase extracted from Citrus sinensis and bioconversion of β-naphthol

A ligninolytic peroxidase called versatile peroxidase, VP, (EC 1.11.1.16) is an iron-containing metalloenzyme. The most distinctive feature of this enzyme is its composite molecular framework, which combines lignin peroxidase's capacity to oxidize compounds with high-redox potential with manganese peroxidase's capacity to oxidize Mn2+ to Mn3+. In this study, we have extracted amino acid sequences from the Citrus sinensis source and subjected them to various computation tools to visualize the insight secondary and 3D structure, physicochemical properties, and validation of the structure which have not been studied so far to further investigate the catalytic efficiency and effectiveness of VP. The binding energies of HEME and HEME C (HEC) ligands with produced PDB (6rqf.1. A) have been also assessed, analyzed, and confirmed utilizing AutoDock. Binding energies were calculated using the AutoDock and validated by MD simulation using SCHRODINGER DESMOND. Most stable confirmation was achieved through a protein-ligand interaction study. Bio-technological use of VP in the biotransformation of β-naphthol has also been studied. The findings in the current study will have a substantial impact on proteomics, biochemistry, biotechnology, and possible uses of versatile peroxidase in the bio-remediation of different toxic organic compounds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03758-x.

The novel anti-cancer feature of Brazzein through activating of hTLR5 by integration of biological evaluation: molecular docking and molecular dynamics simulation

Many of plant proteins exhibit the properties similar to the antitumor proteins although the anticancer activity of Brazzein on modulating the autophagy signaling pathway has not been determined so far. The present study aimed to develop a simplified system to enable the rational design of the activating extracellular domain of human Toll-like receptor 5 (hTLR5). To identify the anticancer effect of Brazzein, HADDOCK program and molecular dynamics (MD) simulation were applied to examine the binding of the wild type (WT) and p.A19K mutant of Brazzein to the TLR5. The expression of MAP1S and TNF-α genes was estimated based on real-time PCR. The results clearly confirmed that the WT of Brazzein activated hTLR5 in the MCF-7 cell line since the genes were more and significantly less expressed in the cells treated with the WT and p.A19K mutant than the control, respectively. The snapshots of MD simulation exhibit the consistent close interactions of hTLR5 with the two helices of Brazzein on its lateral side. The results of per residue-free energy decomposition analysis substantiate those of intermolecular contact analysis perfectly one. We propose that the WT of Brazzein can act as an antitumor drug candidate.

Elucidation of the inhibitory activity of plant-derived SARS-CoV inhibitors and their potential as SARS-CoV-2 inhibitors

Several drugs are now being tested as possible therapies due to the necessity of treating SARS-CoV-2 infection. Although approved vaccines bring much hope, a vaccination program covering the entire global population will take a very long time, making the development of effective antiviral drugs an effective solution for the immediate treatment of this dangerous infection. Previous studies found that three natural compounds, namely, tannic acid, 3-isotheaflavin-3-gallate and theaflavin-3,3-digallate, are effective proteinase (3CL^pro) inhibitors of SARS-CoV (IC₅₀ <10 µM). Based on this information and due to the high sequence identity between SARS-CoV and SARS-CoV-2 3CL^pro, these three compounds could be candidate inhibitors of SARS-CoV-2 3CL^pro. This paper explores the structural and energetic features that guided the molecular recognition of these three compounds for dimeric SARS-CoV-2 and SARS-CoV 3CL^pro, the functional state of this enzyme, using docking and MD simulations with the molecular mechanics-generalized-born surface area (MMGBSA) approach. Energetic analysis demonstrated that the three compounds reached good affinities in both systems in the following order: tannic acid > 3-isotheaflavin-3-gallate > theaflavin-3,3-digallate. This tendency is in line with that experimentally reported between these ligands and SARS-CoV 3CL^pro. Therefore, tannic acid may have clinical usefulness against COVID-19 by acting as a potent inhibitor of SARS-CoV-2 3CL^pro.Communicated by Ramaswamy H. Sarma.

New styrylquinoxaline: synthesis, structural, biological evaluation, ADMET prediction and molecular docking investigations

The organic compound (E)-3-(4-methylstyryl)quinoxalin-2(1H)-one (SQO) with molecular formula C₁₇H₁₄N₂O was synthesized and analyzed using single crystal X-ray diffraction, ¹H, ¹³C NMR and FTIR spectroscopic techniques. The geometric parameters of the molecule was optimized by density-functional theory (DFT) choosing B3LYP with 6-31++G(d,p) basis set. For compatibility, the theoretical structure and experimental structure were overlapped with each other. Frontier molecular orbitals of the title compound were made, and energy gap between HOMO and LUMO was calculated. Molecular electrostatic potential map was generated finding electrophilic and nucleophilic attack centers using DFT method. Hirshfeld surface analysis (HSA) confirms active regions at the circumference of N1 atoms and O1 atoms that form intermolecular N1-H1···O1 hydrogen bond. The acute oral toxicity study was carried out according to OECD guideline, which approve that the compound SQO was non-toxic. In addition, this quinoxaline derivative was evaluated for its in vitro antidiabetic activity against α-glucosidase and α-amylase enzymes and for antioxidant activity by utilizing several tests as 1,1-diphenyl-2-picryl hydrazyl, (2,2'-azino-bis(3-ethyl benzthiazoline-6-sulfonicacid), reducing power test (FRAP) and hydrogen peroxide activity H₂O₂. The molecular docking studies were performed to investigate the antidiabetic activity of SQO and compared with the experimental results. SQO is a potent antidiabetic from both the experimental and molecular docking results. Finally, the physicochemical, pharmacokinetic and toxicological properties of SQO have been evaluated by using in silico absorption, distribution, metabolism, excretion and toxicity analysis prediction.

Integration of NMR studies, computational predictions, and in vitro assays in the search of marine diterpenes with antitumor activity

Among the compounds of natural origin, diterpenes have proved useful as drugs for the treatment of cancer. Marine organisms, such as soft corals and algae, are a promising source of diterpenes, being a rich and unexplored source of cytotoxic agents. This study evaluated a library of 32 natural and semisynthetic marine diterpenes, including briarane, cembrane, and dolabellane nuclei, with the aim of determining their cytotoxicity against three human cancer cell lines (A549, MCF7, and PC3). The three most active compounds were submitted to a flow cytometry analysis in order to determine induction of apoptosis against the A549 cell line. An NMR analysis was conducted to determine and evaluate the interactions between active diterpenes and tubulin. These interactions were characterized by a computational study using molecular docking and MD simulations. With these results, two cembrane and one chlorinated briarane diterpenes were active against the three human cancer cell lines, induced apoptosis in the A549 cell line, and showed interactions with tubulin preferably at the taxane-binding site. This study is a starting point for the identification and optimization of the marine diterpenes selected for better antitumor activities. It also highlights the power of integrating NMR studies, computational predictions, and in vitro assays in the search for compounds with antitumor activity.

Evaluation of binding of potential ADMET/tox screened saquinavir analogues for inhibition of HIV-protease via molecular dynamics and binding free energy calculations

Developing novel drug molecules against HIV is a scientific quest necessitated by development of drug resistance against used drugs. We report comparative results of molecular dynamics simulation studies on 11 structural analogues of Saquinavir (SQV) - against HIV-protease that were earlier examined for pharmacodynamic and pharmacokinetic properties. We reported analogues S1, S5 and S8 to qualify the ADMET criterion and may be considered as potential lead molecules. In this study the designed molecules were successively docked with native HIV-protease at AutoDock. Docking scores established relative goodness of the 11 analogues against the benchmark for Saquinavir. The docked complexes were subjected to molecular dynamics simulation studies using GROMACS 4.6.2. Four parameters viz. H-bonding, RMSD, Binding energy and Protein-Ligand Distance were used for comparative analyses of the analogues relative to Saquinavir. The comparison and analysis of the results are indicative that analogues S8, S9 and S1 are promising candidates among all the analogues studied. From our earlier work and present study it is evident that among the three S8 and S1 qualify the ADMET criterion and between S1 and S8, the analogue S8 shows more target efficacy and specificity over S1 and have better molecular dynamics simulation results. Thus, of the 11 de novo Saquinavir analogues, the S8 appears to be the most promising candidate as lead molecule for HIV-protease inhibitor and is best suited for testing under biological system. Further validation of the proposed lead molecules through wet lab studies involving antiviral assays however is required.Communicated by Ramaswamy H. Sarma.

Synthesis, Biological Evaluation, and Molecular Docking of Combretastatin and Colchicine Derivatives and their hCE1-Activated Prodrugs as Antiviral Agents

Recent studies indicate that tubulin can be a host factor for vector-borne flaviviruses like dengue (DENV) and Zika (ZIKV), and inhibitors of tubulin polymerization such as colchicine have been demonstrated to decrease virus replication. However, toxicity limits the application of these compounds. Herein we report prodrugs based on combretastatin and colchicine derivatives that contain an ester cleavage site for human carboxylesterase, a highly abundant enzyme in monocytes and hepatocytes targeted by DENV. Relative to their parent compounds, the cytotoxicity of these prodrugs was reduced by several orders of magnitude. All synthesized prodrugs containing a leucine ester were hydrolyzed by the esterase in vitro. In contrast to previous reports, the phenylglycine esters were not cleaved by human carboxylesterase. The antiviral activity of combretastatin, colchicine, and selected prodrugs against DENV and ZIKV in cell culture was observed at low micromolar and sub-micromolar concentrations. In addition, docking studies were performed to understand the binding mode of the studied compounds to tubulin.

Further Insight into the Interactions of the Cytotoxic Macrolides Laulimalide and Peloruside A with Their Common Binding Site

The binding site of the macrolides laulimalide and peloruside A, which is different from that of the clinically useful drugs paclitaxel/taxol and ixabepilone (tax site), is known to be between two adjacent β-tubulin units (ext site). Here, we report our study of the binding of these molecules to an α1β1/α2β2-tubulin "tetramer" model. AutoDock 4.2.6//AutoDock Vina dockings predicted that the affinities of laulimalide and peloruside A for the tax site are quite similar to those for the ext site. However, molecular dynamics (MD) simulations indicated that only when these two ligands are located at the ext site, there are contacts that help stabilize the system, favoring the β1/β2 interactions. The binding affinity of laulimalide for this site is stronger than that of peloruside A, but this is compensated for by additional β1/β2 contacts that are induced by peloruside A. MD studies also suggested that epothilones at the tax site and either laulimalide or peloruside A at the ext site cause similar stabilizing effects (mainly linking the M-loop of β1 and loop H1-B2 of β2). In a "hexamer" model (3 units of αβ-tubulin), the effects are confirmed. Metadynamics simulations of laulimalide and peloruside A, which are reported for the first time, suggest that peloruside A produces a stronger change in the M-loop, which explains the stabilization of the β1/β2 interaction.

Syntheses, spectroscopic and crystallographic characterizations of cis- and trans-dispirocyclic ferrocenylphosphazenes: molecular dockings, cytotoxic and antimicrobial activities

The crystallographic and molecular docking investigations and the determinations of the cytotoxic activities of the dispirocyclic ferrocenylphosphazenes were performed.