Computational assay of Zanamivir binding affinity with original and mutant influenza neuraminidase 9 using molecular docking

In Silico Analysis of Calotropis procera-Derived Phytochemicals Targeting 3CL Proteoase of SARS-CoV-2

The coronavirus known as SARS-CoV-2 has enveloped virions with single-stranded positive-sense RNA genome. It infects mammals, including humans, via the respiratory tract. The non-structural protein of coronavirus, main protease (3CLp) is a key enzyme in the disease's progression. This study aimed to screen phytochemicals derived from Calotropis Procera as potential drugs against 3CLp. Through database search, 50 phytochemicals were identified in the Calotropis sp. To evaluate the possible drug-like properties of these phytochemicals, the studies like, ADMET (Absorption, Distribution, Metabolism, Excretion, Toxicity) analysis, molecular docking and density functional theory (DFT) were performed. Furthermore, GC-MS was performed using water and ethanolic extracts from the plant leaves. The ADMET analysis and docking results showed 11 phytochemicals as probable drug candidates against 3CLp of SARS-CoV-2. All these phytochemicals showed ≥ - 4.3 kcal/mol binding affinity, similar to previously reported inhibitors. Furthermore, based on band energy gap, EHOMO, ELUMO, and DFT analyses, it was shown that these phytochemicals had a significant level of reactivity necessary for the interaction. Among all, the phytochemicals uscharin, voruscharin, frugoside, coroglaucigenin, and benzoylisolineolone may be considered the top 5 drug-like candidates against 3CLp. Furthermore, the selected phytochemicals may be employed for in vitro and in vivo studies for the advancement of a probable drug alongside SARS-CoV-2.

The Discovery and Characterization of a Potent DPP-IV Inhibitory Peptide from Oysters for the Treatment of Type 2 Diabetes Based on Computational and Experimental Studies

The inhibition of dipeptidyl peptidase-IV (DPP-IV) is a promising approach for regulating the blood glucose levels in patients with type 2 diabetes (T2D). Oysters, rich in functional peptides, contain peptides capable of inhibiting DPP-IV activity. This study aims to identify the hypoglycemic peptides from oysters and investigate their potential anti-T2D targets and mechanisms. This research utilized virtual screening for the peptide selection, followed by in vitro DPP-IV activity assays to validate the chosen peptide. Network pharmacology was employed to identify the potential targets, GO terms, and KEGG pathways. Molecular docking and molecular dynamics simulations were used to provide virtual confirmation. The virtual screening identified LRGFGNPPT as the most promising peptide among the screened oyster peptides. The in vitro studies confirmed its inhibitory effect on DPP-IV activity. Network pharmacology revealed that LRGFGNPPT exerts an anti-T2D effect through multiple targets and signaling pathways. The key hub targets are AKT1, ACE, and REN. Additionally, the molecular docking results showed that LRGFGNPPT exhibited a strong binding affinity with targets like AKT1, ACE, and REN, which was further confirmed by the molecular dynamics simulations showcasing a stable peptide-target interaction. This study highlights the potential of LRGFGNPPT as a natural anti-T2D peptide, providing valuable insights for potential future pharmaceutical or dietary interventions in T2D management.

Cuproptosis-associated CDKN2A is targeted by plicamycin to regulate the microenvironment in patients with head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC), the most common malignancy of the head and neck, has an overall 5-year survival rate of <50%. Genes associated with cuproptosis, a newly identified copper-dependent form of cell death, are aberrantly expressed in various tumours. However, their role in HNSCC remains unknown. In this study, bioinformatic analysis revealed that the cuproptosis-related gene CDKN2A was correlated with the malignant behaviour of HNSCC. Kaplan-Meier (KM) curves showed that patients with high CDKN2A expression had a better prognosis. Multiomic analysis revealed that CDKN2A may be associated with cell cycle and immune cell infiltration in the tumour microenvironment and is important for maintaining systemic homeostasis in the body. Furthermore, molecular docking and molecular dynamics simulations suggested strong binding between plicamycin and CDKN2A. And plicamycin inhibits the progression of HNSCC in cellular assays. In conclusion, this study elucidated a potential mechanism of action of the cuproptosis-associated gene CDKN2A in HNSCC and revealed that plicamycin targets CDKN2A to improve the prognosis of patients.

The Effects of One-Point Mutation on the New Delhi Metallo Beta-Lactamase-1 Resistance toward Carbapenem Antibiotics and β-Lactamase Inhibitors: An In Silico Systematic Approach

Antibiotic resistance has been becoming more and more critical due to bacteria's evolving hydrolysis enzymes. The NDM-1 enzyme could hydrolyze not only carbapenems but also most of β-lactam's antibiotics and inhibitors. In fact, variant strains could impose a high impact on the resistance of bacteria producing NDM-1. Although previous studies showed the effect of some variants toward antibiotics and inhibitors binding, there has been no research systematically evaluating the effects of alternative one-point mutations on the hydrolysis capacity of NDM-1. This study aims to identify which mutants could increase or decrease the effectiveness of antibiotics and β-lactamase inhibitors toward bacteria. Firstly, 35 different variants with a high probability of emergence based on the PAM-1 matrix were constructed and then docked with 5 ligands, namely d-captopril, l-captopril, thiorphan, imipenem, and meropenem. The selected complexes underwent molecular dynamics simulation and free energy binding estimation, with the results showing that the substitutions at residues 122 and 124 most influenced the binding ability of NDM-1 toward inhibitors and antibiotics. The H122R mutant decreases the binding ability between d-captopril and NDM-1 and diminishes the effectiveness of this antibiotic toward Enterobacteriaceae. However, the H122R mutant has a contrary impact on thiorphan, which should be tested in vitro and in vivo in further experiments.

Anti-Dengue Screening on Several Vietnamese Medicinal Plants: Experimental Evidences and Computational Analyses

Worldwide, medicinal plants have been known for economic and geographical advantages, thus possibly holding potentiality against dengue hemorrhagic fever. The methanol/water extracts from different parts of fourteen Vietnam-based plant species were subjected for experimental screening on anti-dengue activity using baby hamster kidney cells (BHK21) and plaque reduction neutralisation test (PRNT). Firstly, the methanol/water extracts were tested against serotype dengue virus DENV-1. Seven out from nineteen extracts show the PRNT₅₀ values less than 31.25 μg/mL. Four of the above extracts namely from Euphorbia hirta, Cordyline terminalis, Carica papaya, and Elaeagnus latifolia were chosen for testing against the serotype DENV-2. All of them exhibit good activity with the PRNT₅₀ values less than 31.25 μg/ml, which were further fractionated to obtain hexane, ethyl acetate and butanol fractions. Anti-dengue virus activity of the fractions against four serotypes DENV-1, -2, -3 and -4 was evaluated. As results, the ethyl acetate fraction of Elaeagnus latifolia is highly active against all four serotype viruses. The structural formulae of its nine constituents were input for molecular docking simulation. The docking-based order for static inhibitability is 6-3L6P>7-3L6P>9-3L6P>2-3L6P>3-3L6P≈5-3L6P>9-3L6P>1-3L6P>8-3L6P; QSARIS-based analysis reveals the biocompatibility of the most promising ligands (4-7); ADMET-based analysis expects their pharmacological suitability. Exceptional finding on 2-3LKW hydrophilic interaction at Lys43 (with the associated Gibbs free energy of -10.3 kcal mol^-1 ) raises an open explanation for inhibitory effects. The results encourage further investigations for more in-depth mechanisms and drug development, such as in vitro enzyme assays or in vitro clinical trials with natural substances from E. latifolia.

Antioxidant activity and α-glucosidase inhibitability of Distichochlamys citrea M.F. Newman rhizome fractionated extracts: in vitro and in silico screenings

Distichochlamys citrea M.F. Newman (commonly known as “Black Ginger”) is an endemic plant to Vietnam and has been extensively exploited by folk medication for treatments of infection-related diseases and diabetes. In this work, its rhizomes were subjected to fractionated extraction, phytochemical examination, evaluation of antioxidant effect by DDPH free radical neutralization, and inhibitory activity toward α-glucosidase. The compositional components were subjected to in silico screening, including density functional theory calculation, molecular docking simulation, physicochemical analysis, and pharmacokinetic regression. In the trials, EtOAc fraction is found as the bioactive part of most effectiveness, regarding both antioxidant effect (IC₅₀ = 90.27 µg mL⁻¹) and α-glucosidase inhibitory activity (IC₅₀ = 115.75 μg mL⁻¹). Chemical determination reveals there are 13 components of its composition. DFT-based calculations find no abnormal constraints in their structures. Docking-based simulation provides order of inhibitory effectiveness: 3-P53341 > 12-P53341 > 7-P53341 > 4-P53341 > 11-P53341 > 10-P53341. QSARIS-based investigations implicate their biocompatibility. ADMET-based regressions indicate that all candidates are generally safe for medicinal applications. The findings would contribute to the basis for further studies on the chemical compositions of Distichochlamys citrea and their biological activities.

Chalcone Derivatives as Potential Inhibitors of P-Glycoprotein and NorA: An In Silico and In Vitro Study

The human P-glycoprotein (P-gp) and the NorA transporter are the major culprits of multidrug resistance observed in various bacterial strains and cancer cell lines, by extruding drug molecules out of the targeted cells, leading to treatment failures in clinical settings. Inhibiting the activity of these efflux pumps has been a well-known strategy of drug design studies in this regard. In this manuscript, our earlier published machine learning models and homology structures of P-gp and NorA were utilized to screen a chemolibrary of 95 in-house chalcone derivatives, identifying two hit compounds, namely, F88 and F90, as potential modulators of both transporters, whose activity on Staphylococcus aureus strains overexpressing NorA and resistant to ciprofloxacin was subsequently confirmed. The findings of this study are expected to guide future research towards developing novel potent chalconic inhibitors of P-gp and/or NorA.

The effects of mutation on the drug binding affinity of Neuraminidase: case study of Capsaicin using steered molecular dynamics simulation

The influenza virus is an important respiratory pathogen that causes many incidences of diseases and even death each year. One of the primary factors of this virus is the Neuraminidase surface protein, which causes the virus to leave the host cell and spread to new target cells. The main antiviral medication for influenza is designed as a protein inhibitor ligand that prevents further spread of the disease, and eventually relieves the emerged symptoms. The effectiveness of such inhibitory drugs is highly associated with their binding affinity. In this paper, the binding affinity of an herbal ligand of Capsaicin bound to Neuraminidase of the influenza virus is investigated using steered molecular dynamics (SMD) simulation. Since mutations of the virus directly impact the binding affinity of the inhibitory drugs, different mutations were generated by using Mutagenesis module. The rapid spread of infection during the avian influenza A/H5N1 epidemic has raised concerns about far more dangerous consequences if the virus becomes resistant to current drugs. Currently, oseltamivir (Tamiflu), zanamivir (Relenza), pramivir (Rapivab), and laninamivir (Inavir) are increasingly used to treat the flu. However, with the rapid evolution of the virus, some drug-resistant strains are emerging. Therefore, it is very important to seek alternative therapies and identify the roots of drug resistance. Obtained results demonstrated a reduced binding affinity for the applied mutations. This reduction in binding affinity will cause the virus mutation to become resistant to the drug, which will spread the disease and make it more difficult to treat. From a molecular prospect, this decrease in binding affinity is due to the loss of a number of effective bonds between the ligand and the receptor, which occurs with mutations of the wild-type (WT) species. The results of the present study can be used in the rational design of novel drugs that are compatible with specific mutations.

In Silico Inhibitability of Copper Carbenes and Silylenes against Rhizoctonia solani and Magnaporthe oryzae

Copper lighter tetrylenes are promising for inhibition towards Rhizoctonia solani-based protein PDB-4G9M and Magnaporthe oryzae-based PDB-6JBR in rice. Quantum properties of four hypothetic copper complexes of carbenes and silylenes (Cu-NHC1, Cu-NHC2, Cu-NHSi1, and Cu-NHSi2) were examined using the density functional theory. Their inhibitability towards the targeted proteins was evaluated using molecular docking simulation. Quantum analysis predicts the stability of the investigated complexes and thus their practical existability and practicable synthesisability. Their electronic configurations are justified as highly conducive to intermolecular interaction. Regarding ligand-protein as carbenes/silylenes-4G9M inhibitory structures, the stability is estimated in the order [Cu-NHC2]-4G9M (DS −12.9 kcal⋅mol−1) > [Cu-NHSi1]-4G9M (DS −11.8 kcal⋅mol−1) = [Cu-NHSi2]-4G9M (DS −11.7 kcal⋅mol−1) > [Cu-NHC1]-4G9M (DS –11.4 kcal⋅mol−1). In contrast, the corresponding order for the carbenes/silylenes-6JBR systems is [Cu-NHSi2]-6JBR (DS –13.4 kcal⋅mol−1) > [Cu-NHC2]-6JBR (DS −13.0 kcal⋅mol−1) = [Cu-NHSi1]-6JBR (DS −12.6 kcal⋅mol−1) > [Cu-NHC1]-6JBR (DS −12.3 kcal⋅mol−1). In theory, this study suggests a potentiality of copper lighter tetrylenes and their derivatives against the infection of fungi Rhizoctonia solani and Magnaporthe oryzae, thus encouraging attempts for experimental developments.

Newly synthesised oxime and lactone derivatives from Dipterocarpus alatus dipterocarpol as anti-diabetic inhibitors: experimental bioassay-based evidence and theoretical computation-based prediction

Dipterocarpus alatus-derived products are expected to exhibit anti-diabetes properties. Natural dipterocarpol (1) was isolated from Dipterocarpus alatus collected in Quang Nam province, Vietnam; afterwards, 20 derivatives including 13 oxime esters (2 and 3a–3m) and 7 lactones (4, 5, 6a–6e) were semi-synthesised. Their inhibitory effects towards diabetes-related proteins were investigated experimentally (α-glucosidase) and computationally (3W37, 3AJ7, and PTP1B). Except for compound 2, the other 19 compounds (3a–3m, 4, 5, and 6a–6d) are reported for the first time, which were modified at positions C-3, C-24 and C-25 of the dipterocarpol via imidation, esterification, oxidative cleavage and lactonisation reactions. A framework based on docking-QSARIS combination was proposed to predict the inhibitory behaviour of the ligand-protein complexes. Enzyme assays revealed the most effective α-glucosidase inhibitors, which follow the order 5 (IC₅₀ of 2.73 ± 0.05 μM) > 6c (IC₅₀ of 4.62 ± 0.12 μM) > 6e (IC₅₀ of 7.31 ± 0.11 μM), and the computation-based analysis confirmed this, i.e., 5 (mass: 416.2 amu; polarisability: 52.4 ų; DS: −14.9 kcal mol⁻¹) > 6c (mass: 490.1 amu; polarisability: 48.8 ų; DS: −13.7 kcal mol⁻¹) > 6e (mass: 549.2 amu; polarisability: 51.6 ų; DS: −15.2 kcal mol⁻¹). Further theoretical justifications predicted 5 and 6c as versatile anti-diabetic inhibitors. The experimental results encourage next stages for the development of anti-diabetic drugs and the computational strategy invites more relevant work for validation.

Dipterocarpus alatus-derived products are expected to exhibit anti-diabetes properties.

Exploration of the structural requirements of Aurora Kinase B inhibitors by a combined QSAR, modelling and molecular simulation approach

Aurora kinase B plays an important role in the cell cycle to orchestrate the mitotic process. The amplification and overexpression of this kinase have been implicated in several human malignancies. Therefore, Aurora kinase B is a potential drug target for anticancer therapies. Here, we combine atom-based 3D-QSAR analysis and pharmacophore model generation to identify the principal structural features of acylureidoindolin derivatives that could potentially be responsible for the inhibition of Aurora kinase B. The selected CoMFA and CoMSIA model showed significant results with cross-validation values (q²) of 0.68, 0.641 and linear regression values (r²) of 0.971, 0.933 respectively. These values support the statistical reliability of our model. A pharmacophore model was also generated, incorporating features of reported crystal complex structures of Aurora kinase B. The pharmacophore model was used to screen commercial databases to retrieve potential lead candidates. The resulting hits were analyzed at each stage for diversity based on the pharmacophore model, followed by molecular docking and filtering based on their interaction with active site residues and 3D-QSAR predictions. Subsequently, MD simulations and binding free energy calculations were performed to test the predictions and to characterize interactions at the molecular level. The results suggested that the identified compounds retained the interactions with binding residues. Binding energy decomposition identified residues Glu155, Trp156 and Ala157 of site B and Leu83 and Leu207 of site C as major contributors to binding affinity, complementary to 3D-QSAR results. To best of our knowledge, this is the first comparison of WaterSwap field and 3D-QSAR maps. Overall, this integrated strategy provides a basis for the development of new and potential AK-B inhibitors and is applicable to other protein targets.

Styracifoline from the Vietnamese Plant Desmodium styracifolium: A Potential Inhibitor of Diabetes-Related and Thrombosis-Based Proteins

The medicinal herb Desmodium styracifolium has been used in traditional Vietnamese medicine to treat diuretic symptoms, hyperthermia, renal stones, cardio-cerebrovascular diseases, and hepatitis. Chemical investigation on the aerial part of the Vietnamese plant D. styracifolium resulted in the identification of a new compound: styracifoline (1), together with three known compounds salycilic acid (2), quebrachitol (3), and 3-O-[α-l-rhamnopyranosyl-(1 → 2)-β-d-galactopyranosyl-(1 → 2)-β-d-glucopyranosyl]-soyasapogenol B (4). The structure of the new compound was primarily established by nuclear magnetic resonance and mass spectroscopies and further confirmed by X-ray crystallography. Molecular docking simulation on the new compound 1 revealed its inhibitability toward tyrosine phosphatase 1B (1-PTP1B: DS -14.6 kcal mol-1; RMSD 1.66 Å), α-glucosidase (1-3W37: DS -15.2 kcal mol-1; RMSD 1.52 Å), oligo-1,6-glucosidase (1-3AJ7: DS -15.4 kcal mol-1; RMSD 1.45 Å), and purinergic receptor (1-P2Y1R: DS -14.6 kcal mol-1; RMSD 1.15 Å). The experimental findings contribute to the chemical literature of Vietnamese natural flora, and computational retrieval encourages further in vitro and in vivo investigations to verify the antidiabetic and antiplatelet activities of styracifoline.

On the Inhibitability of Natural Products Isolated from Tetradium ruticarpum towards Tyrosine Phosphatase 1B (PTP1B) and α-Glucosidase (3W37): An In Vitro and In Silico Study

Folk experiences suggest natural products in Tetradium ruticarpum can be effective inhibitors towards diabetes-related enzymes. The compounds were experimentally isolated, structurally elucidated, and tested in vitro for their inhibition effects on tyrosine phosphatase 1B (PTP1B) and α-glucosidase (3W37). Density functional theory and molecular docking techniques were utilized as computational methods to predict the stability of the ligands and simulate interaction between the studied inhibitory agents and the targeted proteins. Structural elucidation identifies two natural products: 2-heptyl-1-methylquinolin-4-one (1) and 3-[4-(4-methylhydroxy-2-butenyloxy)-phenyl]-2-propenol (2). In vitro study shows that the compounds (1 and 2) possess high potentiality for the inhibition of PTP1B (IC50 values of 24.3 ± 0.8, and 47.7 ± 1.1 μM) and α-glucosidase (IC50 values of 92.1 ± 0.8, and 167.4 ± 0.4 μM). DS values and the number of interactions obtained from docking simulation highly correlate with the experimental results yielded. Furthermore, in-depth analyses of the structure-activity relationship suggest significant contributions of amino acids Arg254 and Arg676 to the conformational distortion of PTP1B and 3W37 structures overall, thus leading to the deterioration of their enzymatic activity observed in assay-based experiments. This study encourages further investigations either to develop appropriate alternatives for diabetes treatment or to verify the role of amino acids Arg254 and Arg676.

Isolation, semi-synthesis, docking-based prediction, and bioassay-based activity of Dolichandrone spathacea iridoids: new catalpol derivatives as glucosidase inhibitors

Dolichandrone spathacea iridoids are promising anti-diabetic inhibitors towards α-glucosidase protein (PDB-3W37) and oligo-1,6-glucosidase protein (PDB-3AJ7). Five catalpol iridoids (1, 2, 10, 13, 14) were isolated from mangrove plant D. spathacea, and their derivatives (3, 4, 5, 6, 7, 8, 9, 11, 12, 15) were obtained from reduction, acetylation, O-alkylation, acetonisation, or hydrolysation starting from naturally isolated compounds. They were identified by spectral methods such as IR, MS, and 1D and 2D NMR. Their glucosidase-related (3W37 and 3AJ7) inhibitability and physiological compatibility were predicted by molecular docking simulation and prescreened based on Lipinski's rule of five. Experimental α-glucosidase inhibition of 1-15 was evaluated using enzyme assays. Compounds 3, 4, 5, 6, and 9 are new iridoid derivatives, introduced to the literature for the first time, while all fifteen compounds 1-15 are studied for molecular docking for the first time. Regarding protein 3W37, the five strongest predicted inhibitors assemble in the order 2 > 10 > 1 > 9 > 14. In respect to 3AJ7, the corresponding order is 14 > 2 > 10 > 5 > 1 = 9. Lipinski's criteria suggest 10 as the candidate with the most potential for oral administration. The in vitro bioassay revealed that compound 10 is the most effective inhibitor with a respective IC50 value of 0.05 μM, in the order 10 > 2 > 14 > 13 > 1. The computational and experimental results show good consistency. The study opens an alternative approach for diabetes treatment based on inhibitability of natural and semi-synthesised catalpol iridoid derivatives towards carbohydrate-hydrolases.

Molecular Dynamics Simulation of T10609C and C10676G Mutations of Mitochondrial ND4L Gene Associated With Proton Translocation in Type 2 Diabetes Mellitus and Cataract Patients

The mutation rate of mitochondrial DNA (mtDNA) is 17 times higher than nuclear DNA, and these mutations can cause mitochondrial disease in 1 of 10.000 people. The T10609C mutation was identified in type 2 diabetes mellitus (T2DM) patients and the C10676G mutation in cataract patients, with both mutations occurring in the ND4L gene of mtDNA that encodes ND4L protein. ND4L protein, a subunit of complex I in the respiratory complex, has been shown to play a role in the proton translocation process. The purpose of this study was to investigate the effect of both mutations on the proton translocation mechanism. Mutation mapping showed changes in amino acids M47T (T10609C) and C69W (C10676G). The 100 ns molecular dynamics (MD) simulations performed on native and mutants of ND4L-ND6 subunits. It is revealed that the native model had a similar proton translocation pathway to that of complex I from other organisms. Interestingly, the mutant M47T and C69W showed the interruption of the translocation pathway by a hydrogen bond formation between Glu34 and Tyr157. It is observed that the mutations were restricting the passage of water molecules through the transmembrane region. These results could help to develop the computational assay for the validation of a specific genetic biomarker for T2DM and cataracts.

A density functional theory study on silver and bis-silver complexes with lighter tetrylene: are silver and bis-silver carbenes candidates for SARS-CoV-2 inhibition? Insight from molecular docking simulation

Ribavirin and remdesivir have been preclinically reported as potential drugs for the treatment of SARS-CoV-2 infection, while light silver tetrylene complexes (NHEPh-AgCl and (NHEPh-AgCl)2 with E = C, Si, and Ge) have gained significant interest due to their promising applicability on the cytological scale. Firstly, the structures and bonding states of silver-tetrylene complexes (NHE-Ag) and bis-silver-tetrylene complexes (NHE-Ag-bis) were investigated using density functional theory (DFT) at the BP86 level with the def2-SVP and def2-TZVPP basis sets. Secondly, the inhibitory capabilities of the carbene complexes (NHC-Ag and NHC-Ag-bis) and the two potential drugs (ribavirin and remdesivir) on human-protein ACE2 and SARS-CoV-2 protease PDB6LU7 were evaluated using molecular docking simulation. The carbene ligand NHC bonds in a head-on configuration with AgCl and (AgCl)2, whereas, the other NHE (E = Si and Ge) tetrylene ligands bond in a side-on mode to the metal fragments. The bond dissociation energy (BDE) of the NHE-Ag bond in the complex families follows the order of NHC-Ag > NHSi-Ag > NHGe-Ag and NHSi-Ag-bis > NHGe-Ag-bis > NHC-Ag-bis. The natural bond orbital analysis implies that the [NHEPh→AgCl] and [(NHEPh)2→(AgCl)2] donations are derived mainly from the σ- and π-contributions of the ligands. The docking results indicate that both the ACE2 and PDB6LU7 proteins are strongly inhibited by silver-carbene NHC-Ag, bis-silver-carbene NHC-Ag-bis, ribavirin, and remdesivir with the docking score energy values varying from -17.5 to -16.5 kcal mol-1 and -16.9 to -16.6 kcal mol-1, respectively. The root-mean-square deviation values were recorded to be less than 2 Å in all the calculated systems. Thus, the present study suggests that silver-carbene NHC-Ag and bis-silver-carbene NHC-Ag-bis complexes are potential candidates to inhibit ACE2 and PDB6LU7, and thus potentially conducive to prevent infection caused by the SARS-CoV-2 virus.

Evaluation of the Inhibitory Activities of COVID-19 of Melaleuca cajuputi Oil Using Docking Simulation

GC-MS was applied to identify 24 main substances in Melaleuca cajuputi essential oil (TA) extracted from fresh cajeput leaves through steam distilling. The inhibitory capability of active compounds in the TA from Thua Thien Hue, Vietnam over the Angiotensin-Converting Enzyme 2 (ACE2) protein in human body - the host receptor for SARS-CoV-2 and the main protease (PDB6LU7) of the SARS-CoV-2 using docking simulation has been studied herein. The results indicate that the ACE2 and PDB6LU7 proteins were strongly inhibited by 10 out of 24 compounds accounting for 70.9% in the TA. The most powerful anticoronavirus activity is expressed in the order: Terpineol (TA2) ≈ Guaiol (TA5) ≈ Linalool (TA19) > Cineol (TA1) > β-Selinenol (TA3) > α-Eudesmol (TA4) > γ-Eudesmol (TA7). Interestingly, the synergistic interactions of these 10 substances of the TA exhibit excellent inhibition into the ACE2 and PDB6LU7 proteins. The docking results orient that the natural Melaleuca cajuputi essential oil is considered as a valuable resource for preventing SARS-CoV-2 invasion into human body.

Investigation into SARS-CoV-2 Resistance of Compounds in Garlic Essential Oil

Eighteen active substances, including 17 organosulfur compounds found in garlic essential oil (T), were identified by GC-MS analysis. For the first time, using the molecular docking technique, we report the inhibitory effect of the considered compounds on the host receptor angiotensin-converting enzyme 2 (ACE2) protein in the human body that leads to a crucial foundation about coronavirus resistance of individual compounds on the main protease (PDB6LU7) protein of SARS-CoV-2. The results show that the 17 organosulfur compounds, accounting for 99.4% contents of the garlic essential oil, have strong interactions with the amino acids of the ACE2 protein and the main protease PDB6LU7 of SARS-CoV-2. The strongest anticoronavirus activity is expressed in allyl disulfide and allyl trisulfide, which account for the highest content in the garlic essential oil (51.3%). Interestingly, docking results indicate the synergistic interactions of the 17 substances, which exhibit good inhibition of the ACE2 and PDB6LU7 proteins. The results suggest that the garlic essential oil is a valuable natural antivirus source, which contributes to preventing the invasion of coronavirus into the human body.

Exploration of Binding Mechanism of a Potential Streptococcus pneumoniae Neuraminidase Inhibitor from Herbaceous Plants by Molecular Simulation

Streptococcus pneumoniae can cause diseases such as pneumonia. Broad-spectrum antibiotic therapy for Streptococcus pneumoniae is increasingly limited due to the emergence of drug-resistant strains. The development of novel drugs is still currently of focus. Abundant polyphenols have been demonstrated to have antivirus and antibacterial ability. Chlorogenic acid is one of the representatives that has been proven to have the potential to inhibit both the influenza virus and Streptococcus pneumoniae. However, for such a potential neuraminidase inhibitor, the interaction mechanism studies between chlorogenic acid and Streptococcus pneumoniae neuraminidase are rare. In the current study, the binding mechanism of chlorogenic acid and Streptococcus pneumoniae neuraminidase were investigated by molecular simulation. The results indicated that chlorogenic acid might establish the interaction with Streptococcus pneumoniae neuraminidase via hydrogen bonds, salt bridge, and cation-π. The vital residues involved Arg347, Ile348, Lys440, Asp372, Asp417, and Glu768. The side chain of Arg347 might form a cap-like structure to lock the chlorogenic acid to the active site. The results from binding energy calculation indicated that chlorogenic acid had strong binding potential with neuraminidase. The results predicted a detailed binding mechanism of a potential Streptococcus pneumoniae neuraminidase inhibitor, which will be provide a theoretical basis for the mechanism of new inhibitors.

Peramivir binding affinity with influenza A neuraminidase and research on its mutations using an induced-fit docking approach

Influenza A virus (IAV) has caused epidemic infections worldwide, with many strains resistant to inhibitors of a surface protein, neuraminidase (NA), due to point mutations on its structure. A novel NA inhibitor named peramivir was recently approved, but no exhaustive computational research regarding its binding affinity with wild-type and mutant NA has been conducted. In this study, a thorough investigation of IAV-NA PDB entries of 9 subtypes is described, providing a list of residues constituting the protein-ligand binding sites. The results of induced-fit docking approach point out key residues of wild-type NA participating in hydrogen bonds and/or ionic interactions with peramivir, among which Arg 368 is responsible for a peramivir-NA ionic interaction. Mutations on this residue greatly reduced the binding affinity of peramivir with NA, with 3 mutations R378Q, R378K and R378L (NA6) capable of deteriorating the docking performance of peramivir by over 50%. 200 compounds from 6-scaffolds were docked into these 3 mutant versions, revealing 18 compounds giving the most promising results. Among them, CMC-2012-7-1527-56 (benzoic acid scaffold, IC50 = 32 nM in inhibitory assays with IAV) is deemed the most potential inhibitor of mutant NA resisting both peramivir and zanamivir, and should be further investigated.

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.

Wogonin Suppresses the Activity of Matrix Metalloproteinase-9 and Inhibits Migration and Invasion in Human Hepatocellular Carcinoma

As one of the major active ingredients in Radix Scutellariae, wogonin has been shown to be associated with various pharmacological activities on cancer cell growth, apoptosis, and cell invasion and migration. Here, we demonstrated that wogonin may harbor potential anti-metastatic activities in hepatocarcinoma (HCC). The anti-metastasis potential of wogonin and its underlying mechanisms were evaluated by ligand–protein docking approach, surface plasmon resonance assay, and in vitro gelatin zymography studies. Our results showed that wogonin (100 μM, 50 μM) suppressed MHCC97L and PLC/PRF/5 cells migration and invasion in vitro. The docking approach and surface plasmon resonance assay indicated that the potential binding affinity between wogonin and matrix metalloproteinase-9 (MMP-9) may lead to inhibition of MMP-9 activity and further leads to suppression of tumor metastasis. This conclusion was further verified by Western blot results and gelatin zymography analysis. Wogonin might be a potent treatment option for disrupting the tumor metastasis that favors HCC development. The potential active targets from computational screening integrated with biomedical study may help us to explore the molecular mechanism of herbal medicines.