In silico detection of potential inhibitors from vitamins and their derivatives compounds against SARS-CoV-2 main protease by using molecular docking, molecular dynamic simulation and ADMET profiling

Structure-based virtual screening, ADMET analysis, and molecular dynamics simulation of Moroccan natural compounds as candidates for the SARS-CoV-2 inhibitors

The lack of treatments and vaccines effective against SARS-CoV-2 has forced us to explore natural compounds that could potentially inhibit this virus. Additionally, Morocco is renowned for its rich plant diversity and traditional medicinal uses, which inspires us to leverage our cultural heritage and the abundance of natural resources in our country for therapeutic purposes. In this study, an extensive investigation was conducted to gather a collection of phytoconstituents extracted from Moroccan plants, aiming to evaluate their ability to inhibit the proliferation of the SARS-CoV-2 virus. Molecular docking of the studied compounds was performed at the active sites of the main protease (6lu7) and spike (6m0j) proteins to assess their binding affinity to these target proteins. Compounds exhibiting high affinity to the proteins underwent further evaluation based on Lipinski's rule and ADME-Tox analysis to gain insights into their oral bioavailability and safety. The results revealed that the two compounds demonstrated strong binding affinity to the target proteins, making them potential candidates for oral antiviral drugs against SARS-CoV-2. The molecular dynamics results from this computational analysis supported the overall stability of the resulting complex.

Identification of compounds from natural Peruvian sources as potential inhibitors of SARS-CoV-2 Mpro mutations by virtual screening and computational simulations

Background

Although the COVID-19 pandemic has diminished in intensity, the virus continues to circulate globally. The SARS-CoV-2 main protease (Mpro) is a key enzyme in the life cycle of the virus, making it important for the development of treatments against future variants of the virus. In this work, Peruvian natural compounds were evaluated against different mutations of the SARS-CoV-2 Mpro.

Methods

In silico techniques such as virtual screening, all-atom molecular dynamics simulations, and energy estimation analysis were applied.

Results

Of the tested compounds by virtual screening, rutin was identified as the best binding agent against the different proposed Mpro mutations. In addition, computational simulations and energy estimation analysis demonstrated the high structural and energetic stability between the Mpro-rutin systems.

Conclusions

Overall, our study identified rutin as the most promising compound with a strong affinity for various Mpro mutations, potentially playing a key role in the development of new treatments for emerging viral variants.

Computational exploration of compounds in Xylocarpus granatum as a potential inhibitor of Plasmodium berghei using docking, molecular dynamics, and DFT studies

Malaria remains a global health concern, with the emergence of resistance to the antimalarial drug atovaquone through cytochrome b (cyt b) being well-documented. This study was prompted by the presence of this mutation in cyt b to enable new drug candidates capable of overcoming drug resistance. Our objective was to identify potential drug candidates from compounds of Xylocarpus granatum by computationally assessing their interactions with Plasmodium berghei cyt b. Using computational methods, we modeled cyt b (GenBank: AF146076.1), identified the binding cavity, and analyzed the Ramachandran plot against cyt b. Additionally, we conducted drug-likeness and absorption, distribution, metabolism, excretion, and toxicity (ADMET) studies, along with density functional theory (DFT) analysis of the compounds. Molecular docking and molecular dynamics simulation (MDS) were used to evaluate the binding energy and stability of the cyt b-ligand complex. Notably, our investigation highlighted kaempferol as a promising compound due to its high binding energy of 7.67 kcal/mol among all X. granatum compounds, coupled with favorable pharmacological properties (ADMET) and antiprotozoal properties at Pa 0.345 > Pi 0.009 (PASS value). DFT analysis showed that kaempferol has an energy gap of 4.514 eV. MDS indicated that all tested ligands caused changes in bonding and affected the structural conformation of cyt b, as observed before MDS (0 ns) and after MDS (100 ns). The most notable differences were observed in the types of hydrogen bonds between 0 and 100 ns. Nevertheles, MDS results from a 100 ns simulation revealed consistent behavior for kaempferol across various parameters including root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA), molecular mechanics-Poisson Boltzmann surface area (MM-PBSA), and hydrogen bonds. The cyt b-kaempferol complex demonstrated favorable energy stability, as supported by the internal energy distribution values observed in principal component analysis (PCA), which closely resembled those of the atovaquone control. Additionally, trajectory stability analysis indicated structural stability, with a cumulative eigenvalue of 24.7 %. Dynamic cross-correlation matrix (DCCM) analysis revealed a positive correlation among catalytic cytochrome residues within the amino acid residues range 119-268. The results of our research indicate that the structure of kaempferol holds promise as a potential candidate against Plasmodium.

Molecular docking and molecular dynamics simulation of wheat gluten-derived antioxidant peptides acting through the Keap1-Nrf2 pathway

BACKGROUND

In our previous study, we successfully identified five peptides from wheat gluten: Ala-Pro-Ser-Tyr (APSY), Leu-Tyr (LY), Pro-Tyr (PY), Arg-Gly-Gly-Tyr (RGGY) and Tyr-Gln (YQ). Molecular docking and molecular dynamics simulation methods were employed to investigate the interaction between these antioxidant peptides and the Kelch-like ECH-associated protein 1 (Keap1 protein), revealing the molecular mechanism of their non-competitive binding. In addition, the total antioxidant capacity of the five peptides was determined using the 2,2'-azinobis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt (ABTS) method.

RESULTS

The affinities of APSY, LY, PY, RGGY and YQ were -8.9, -8.3, -8.5, -9.1 and - 7.9 kcal mol-1, respectively. The five peptides effectively bound to Keap1 protein through hydrogen, π-σ, π-alkyl and alkyl interactions. Significant roles were observed for the P1 pocket residue ARG-415 and the P3 pocket residue ALA-556 in the interactions of the Keap1-peptide complexes. Molecular dynamics simulations further elucidated the dynamic process of peptide binding to the Keap1 protein. All five peptides formed stable complexes with Keap1 protein, with van der Waals forces playing crucial roles in these complex systems, indicative of the peptides' strong binding ability to Keap1 protein. The van der Waals forces were -178.74, -123.11, -134.36, -132.59, and -121.44 kJ mol-1 for the Keap1-APSY, Keap1-LY, Keap1-PY, Keap1-RGGY and Keap1-YQ complexes, respectively. These peptides exhibited excellent antioxidant effects. Among them, the YQ peptide exhibited the highest total antioxidant capacity, with an activity value of 1.18 ± 0.06 mmol Trolox equivalent (TE) L-1 at a concentration of 0.10 mg mL-1. The RGGY, PY, LY and APSY peptides followed in descending order, with activity values of 0.91 ± 0.05, 0.72 ± 0.06, 0.62 ± 0.04 and 0.60 ± 0.05 mmol TE L-1, respectively.

CONCLUSION

These results unveiled the molecular mechanism by which the five antioxidant peptides act on active pockets through the Keap1-Nrf2 signaling pathway, providing a theoretical basis for the development of antioxidants. © 2024 Society of Chemical Industry.

Molecular insight into binding affinities and blockade effects of selected flavonoid compounds on the PD-1/PD-L1 pathway

This study investigated the binding mechanisms of the flavonoids apigenin (Api), kaempferol (Kmp), and quercetin (Que) to the PD-L1 dimer using a combination of molecular modeling and experimental techniques. The binding free energy results demonstrated that the flavonoids could tightly bind to the PD-L1 dimer, with the binding abilities following the trend Que > Kmp > Api. Key residues Ile54, Tyr56, Met115, Ala121, and Tyr123 were identified as important for binding. The flavonoids primarily bind to the C-, F-, and G-sheet domains. The spontaneous formation of the complex systems was mainly driven by hydrophobic forces. Dynamic cross-correlation matrix and secondary structure analyses further indicated that the studied flavonoids could stably interact with the binding sites. ELISA results showed that the flavonoids could effectively block PD-1/PD-L1 interactions, although the inhibitory activity of Api was weaker. Therefore, flavonols might be more effective inhibitors compared to flavones. The findings of this study are expected to contribute to the development of novel flavonoids targeting the PD-1/PD-L1 pathway.

Unveiling potent inhibitors for schistosomiasis through ligand-based drug design, molecular docking, molecular dynamics simulations and pharmacokinetics predictions

Schistosomiasis is a neglected tropical disease which imposes a considerable and enduring impact on affected regions, leading to persistent morbidity, hindering child development, diminishing productivity, and imposing economic burdens. Due to the emergence of drug resistance and limited management options, there is need to develop additional effective inhibitors for schistosomiasis. In view of this, quantitative structure-activity relationship studies, molecular docking, molecular dynamics simulations, drug-likeness and pharmacokinetics predictions were applied to 39 Schistosoma mansoni Thioredoxin Glutathione Reductase (SmTGR) inhibitors. The chosen QSAR model demonstrated robust statistical parameters, including an R2 of 0.798, R2adj of 0.767, Q2cv of 0.681, LOF of 0.930, R2test of 0.776, and cR2p of 0.746, confirming its reliability. The most active derivative (compound 40) was identified as a lead candidate for the development of new potential non-covalent inhibitors through ligand-based design. Subsequently, 12 novel compounds (40a-40l) were designed with enhanced anti-schistosomiasis activity and binding affinity. Molecular docking studies revealed strong and stable interactions, including hydrogen bonding, between the designed compounds and the target receptor. Molecular dynamics simulations over 100 nanoseconds and MM-PBSA free binding energy (ΔGbind) calculations validated the stability of the two best-designed molecules. Furthermore, drug-likeness and pharmacokinetics prediction analyses affirmed the potential of these designed compounds, suggesting their promise as innovative agents for the treatment of schistosomiasis.

Structure and Ligand Based Virtual Screening and MPO Topological Analysis of Triazolo Thiadiazepine-fused Coumarin Derivatives as Anti-Parkinson Drug Candidates

Parkinson's disease (PD) is a debilitating condition that can cause locomotor problems in affected patients, such as tremors and body rigidity. PD therapy often includes the use of monoamine oxidase B (MAOB) inhibitors, particularly phenylhalogen compounds and coumarin-based semi-synthetic compounds. The objective of this study was to analyze the structural, pharmacokinetic, and pharmacodynamic profile of a series of Triazolo Thiadiazepine-fused Coumarin Derivatives (TDCDs) against MAOB, in comparison with the inhibitor safinamide. To achieve this goal, we utilized structure-based virtual screening techniques, including target prediction and absorption, distribution, metabolism, and excretion (ADME) prediction based on multi-parameter optimization (MPO) topological analysis, as well as ligand-based virtual screening techniques, such as docking and molecular dynamics. The findings indicate that the TDCDs exhibit structural similarity to other bioactive compounds containing coumarin and MAOB-binding azoles, which are present in the ChEMBL database. The topological analyses suggest that TDCD3 has the best ADME profile, particularly due to the alignment between low lipophilicity and high polarity. The coumarin and triazole portions make a strong contribution to this profile, resulting in a permeability with Papp estimated at 2.15 × 10-5 cm/s, indicating high cell viability. The substance is predicted to be metabolically stable. It is important to note that this is an objective evaluation based on the available data. Molecular docking simulations showed that the ligand has an affinity energy of - 8.075 kcal/mol with MAOB and interacts with biological substrate residues such as Pro102 and Phe103. The results suggest that the compound has a safe profile in relation to the MAOB model, making it a promising active ingredient for the treatment of PD.

Exploring the Efficacy of Noncovalent SARS-CoV-2 Main Protease Inhibitors: A Computational Simulation Analysis Study

The SARS-CoV-2 main protease, as a key target for antiviral therapeutics, is instrumental in maintaining virus stability, facilitating translation, and enabling the virus to evade innate immunity. Our research focused on designing non-covalent inhibitors to counteract the action of this protease. Utilizing a 3D-QSAR model and contour map, we successfully engineered eight novel non-covalent inhibitors. Further evaluation and comparison of these novel compounds through methodologies including molecular docking, ADMET analysis, frontier molecular orbital studies, molecular dynamics simulations, and binding free energy revealed that the inhibitors N02 and N03 demonstrated superior research performance (N02 ΔGbind=-206.648 kJ/mol, N03 ΔGbind=-185.602 kJ/mol). These findings offer insightful guidance for the further refinement of molecular structures and the development of more efficacious inhibitors. Consequently, future investigations can draw upon these findings to unearth more potent inhibitors, thereby amplifying their impact in the treatment and prevention of associated diseases.

2,4-disubstituted 6-fluoroquinolines as potent antiplasmodial agents: QSAR, homology modeling, molecular docking and ADMET studies

Objective

This work was designed to study 2,4-disubstituted 6-fluoroquinolines as antiplasmodial agents by using in silico techniques, to aid in the design of novel analogs with high potency against malaria and high inhibition of Plasmodium falciparum translation elongation factor 2 (PfeEF2), a novel drug target.

Methods

Quantitative structure-activity relationships (QSAR) of 2,4-disubstituted 6-fluoroquinolines were studied with the genetic function approximation technique in Material Studio software. The 3D structure of PfeEF2 was modeled in the SWISS-MODEL workspace through homology modeling. A molecular docking study of the modeled PfeEF2 and 2,4-disubstituted 6-fluoroquinolines was conducted with Autodock Vina in Pyrx software. Furthermore, the in silico pharmacokinetic properties of selected compounds were investigated.

Results

A robust, reliable and predictive QSAR model was developed that related the chemical structures of 2,4-disubstituted 6-fluoroquinolines to their antiplasmodium activities. The model had an internal squared correlation coefficient R2 of 0.921, adjusted squared correlation coefficient R2adj of 0.878, leave-one-out cross-validation coefficient Q2cv of 0.801 and predictive squared correlation coefficient R2pred of 0.901. The antiplasmodium activity of 6-fluoroquinolines was found to depend on the n5Ring, GGI9, TDB7u, TDB8u and RDF75i physicochemical properties: n5Ring, TDB8u and RDF75i were positively associated, whereas GGI9 and TDB7u were negatively associated, with the antiplasmodium activity of the compounds. Stable complexes formed between the compounds and modeled PfeEF2, with binding affinity ranging from -8.200 to -10.700 kcal/mol. Compounds 5, 11, 16, 22 and 24 had better binding affinities than quinoline-4-carboxamide (DDD107498), as well as good pharmacokinetic properties, and therefore may be better inhibitors of this novel target.

Conclusion

QSAR and docking studies provided insight into designing novel 2,4-disubstituted 6-fluoroquinolines with high antiplasmodial activity and good structural properties for inhibiting a novel antimalarial drug target.

Evaluation of molecular mechanisms of riboflavin anti-COVID-19 action reveals anti-inflammatory efficacy rather than antiviral activity

BACKGROUND

Riboflavin (vitamin B2) is one of the most important water-soluble vitamins and a coenzyme involved in many biochemical processes. It has previously been shown that adjuvant therapy with flavin mononucleotide (a water-soluble form of riboflavin) correlates with normalization of clinically relevant immune markers in patients with COVID-19, but the mechanism of this effect remains unclear. Here, the antiviral and anti-inflammatory effects of riboflavin were investigated to elucidate the molecular mechanisms underlying the riboflavin-induced effects.

METHODS

Riboflavin was evaluated for recombinant SARS-CoV-2 PLpro inhibition in an enzyme kinetic assay and for direct inhibition of SARS-CoV-2 replication in Vero E6 cells, as well as for anti-inflammatory activity in polysaccharide-induced inflammation models, including endothelial cells in vitro and acute lung inflammation in vivo.

RESULTS

For the first time, the ability of riboflavin at high concentrations (above 50 μM) to inhibit SARS-CoV-2 PLpro protease in vitro was demonstrated; however, no inhibition of viral replication in Vero E6 cells in vitro was found. At the same time, riboflavin exerted a pronounced anti-inflammatory effect in the polysaccharide-induced inflammation model, both in vitro, preventing polysaccharide-induced cell death, and in vivo, reducing inflammatory markers (IL-1β, IL-6, and TNF-α) and normalizing lung histology.

CONCLUSIONS

It is concluded that riboflavin reveals anti-inflammatory rather than antiviral activity for SARS-CoV-2 infection.

GENERAL SIGNIFICANCE

Riboflavin could be suggested as a promising compound for the therapy of inflammatory diseases of broad origin.

Genotoxic effects and mitosis aberrations of chromium (VI) on root cells of Vicia faba and its molecular docking analysis

Like other heavy metals, Cr (VI) is a powerful carcinogen and mutagen agent. Its toxic effects on plants are well considered. In order to elucidate its adverse effects, the present work aims to study the mitosis aberrations of Cr (VI) on the Vicia faba root-cells and its molecular docking analysis to understand the genotoxicity mechanisms. In-vivo, Vicia faba plants were exposed to 50 and 100 μM Cr (VI) for 48 h. In-silico, molecular docking and molecular dynamics simulation were used to study the interactions between dichromate and tubulin tyrosine ligase T2R-TTL (PDBID: 5XIW) with reference to Colchicine (microtubule inhibitor). According to our results, Cr (VI) affects growth and cell division and also induces many mitosis aberrations such as chromosome sticking, anaphase/telophase bridges, lagging chromosomes and fragmentation during all phases of mitosis. On the one hand, Cr (VI) reduces mitotic index and promotes micronuclei induction. The in-silico results showed that dichromate establishes very strong bonds at the binding site of the tubulin tyrosine ligase T2R-TTL, with a binding affinity of -5.17 Kcal/Mol and an inhibition constant of 163.59 μM. These interactions are similar to those of colchicine with this protein, so dichromate could be a very potent inhibitor of this protein's activity. TTL plays a fundamental role in the tyrosination/detyrosination of tubulin, which is crucial to the regulation of the microtubule cytoskeleton. Its inhibition leads to the appearance of many morphogenic abnormalities such as mitosis aberrations. In conclusion, our data confirm the highest genotoxicity effects of Cr (VI) on Vicia faba root-cells.

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.

Computational Estimation of Residues Involving Resistance to the SARS-CoV-2 Main Protease Inhibitor Ensitrelvir Based on Virtual Alanine Scan of the Active Site

Ensitrelvir is a noncovalent inhibitor of the main protease (Mpro) of severe acute respiratory syndrome coronavirus 2. Acquisition of drug resistance in virus-derived proteins is a serious therapeutic concern, and drug resistance occurs due to amino acid mutations. In this study, we computationally constructed 24 mutants, in which one residue around the active site was replaced with alanine and performed molecular dynamics simulations to the complex of Mpro and ensitrelvir to predict the residues involved in drug resistance. We evaluated the changes in the entire protein structure and ligand configuration in each of these mutants and estimated which residues were involved in ensitrelvir recognition. This method is called a virtual alanine scan. In nine mutants (S1A, T26A, H41A, M49A, L141A, H163A, E166A, V186A, and R188A), although the entire protein structure and catalytic dyad (cysteine (Cys)145 and histidine (His)41) were not significantly moved, the ensitrelvir configuration changed. Thus, it is considered that these mutants did not recognize ensitrelvir while maintaining Mpro enzymatic activities, and Ser1, Thr26, His41, Met49, Leu141, His163, Glu166, Val186, and Arg188 may be related to ensitrelvir resistance. The ligand shift noted in M49A was similar to that observed in M49I, which has been shown to be experimentally ensitrelvir resistant. These findings suggest that our research approach can predict mutations that incite drug resistance.

2-Amino Thiazole Derivatives as Prospective Aurora Kinase Inhibitors against Breast Cancer: QSAR, ADMET Prediction, Molecular Docking, and Molecular Dynamic Simulation Studies

The aurora kinase is a key enzyme that is implicated in tumor growth. Research revealed that small molecules that target aurora kinase have beneficial effects as anticancer agents. In the present study, in order to identify potential antibreast cancer agents with aurora kinase inhibitory activity, we employed QSARINS software to perform the quantitative structure-activity relationship (QSAR). The statistical values resulted from the study include R2 = 0.8902, CCCtr = 0.7580, Q2 LOO = 0.7875, Q2LMO = 0.7624, CCCcv = 0.7535, R2ext = 0.8735, and CCCext = 0.8783. Among the four generated models, the two best models encompass five important variables, including PSA, EstateVSA5, MoRSEP3, MATSp5, and RDFC24. The parameters including the atomic volume, atomic charges, and Sanderson's electronegativity played an important role in designing newer lead compounds. Based on the above data, we have designed six series of compounds including 1a-e, 2a-e, 3a-e, 4a-e, 5a-e, and 6a-e. All these compounds were subjected to molecular docking studies by using AutoDock v4.2.6 against the aurora kinase protein (1MQ4). Among the above 30 compounds, the 2-amino thiazole derivatives 1a, 2a, 3e, 4d, 5d, and 6d have excellent binding interactions with the active site of 1MQ4. Compound 1a had the highest docking score (-9.67) and hence was additionally subjected to molecular dynamic simulation investigations for 100 ns. The stable binding of compound 1a with 1MQ4 was verified by RMSD, RMSF, RoG, H-bond, molecular mechanics-generalized Born surface area (MM-GBSA), free binding energy calculations, and solvent-accessible surface area (SASA) analyses. Furthermore, newly designed compound 1a exhibited excellent ADMET properties. Based on the above findings, we propose that the designed compound 1a may be utilized as the best theoretical lead for future experimental research of selective inhibition of aurora kinase, therefore assisting in the creation of new antibreast cancer drugs.

SARS-CoV-2 main protease mutation analysis via a kinematic method

The Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) is the virus responsible for the COVID-19 pandemic. COVID-19 continues to cause millions of deaths globally in part due to immune-evading mutations. SARS-CoV-2 main protease (Mpro) is an important enzyme for viral replication and potentially an effective drug target. Mutations affect the dynamics of enzymes and thereby their activity and ability to bind ligands. Here, we use kinematic flexibility analysis (KFA) to identify how mutations and ligand binding changes the conformational flexibility of Mpro. KFA decomposes macromolecules into regions of different flexibility near-instantly from a static structure, allowing conformational dynamics analysis at scale. Altogether, we analyzed 47 mutation sites across 69 Mpro-ligand complexes resulting in more than 3300 different structures which includes 69 mutated structures with all 47 sites mutated simultaneously and 3243 single residue mutated structures. We found that mutations generally increased the conformational flexibility of the protein. Understanding the impact of mutations on the flexibility of Mpro is essential for identifying potential drug targets in the treatment of SARS-CoV-2. Further studies in this area can offer valuable insights into the mechanisms of molecular recognition.

Biological investigations of Aspergillus ficuum via in vivo, in vitro and in silico analyses

Serious human health impacts have been observed worldwide due to several life-threatening diseases such as cancer, candidiasis, hepatic coma, and gastritis etc. Exploration of nature for the treatment of such fatal diseases is an area of immense interest for the scientific community. Based on this idea, the genus Aspergillus was selected to discover its hidden therapeutic potential. The genus Aspergillus is known to possess several biologically active compounds. The current research aimed to assess the biological and pharmacological potency of the extracts of less-studied Aspergillus ficuum (FCBP-DNA-1266) (A. ficuum) employing experimental and bioinformatics approaches. The disc diffusion method was used for the antifungal investigation, and the MTT assay was performed to assess the anticancer effects. Mice were employed as an in vivo model to evaluate the antispasmodic effects. A standard spectrophotometric technique was applied to gauge the urease inhibitory activity. The antifungal studies indicate that both n-hexane and ethyl acetate extracts were significantly active against Candida albicans (C. albicans) with their zone of inhibitions (ZOI) values reported as 19 ± 1.06 mm and 25 ± 0.55 mm, respectively at a dose of 30 µg.mL-1. In vitro cytotoxicity assay against HeLa, fibroblast 3T3, prostate PC3, and breast MCF-7 cancer cell lines was performed. The ethyl acetate extract of A. ficuum was found to be significantly active against MCF-7 with its IC50 value of 43.88 µg.mL-1. However, no substantial effects on the percent cell death of HeLa cancer cell lines were observed. In addition, the A. ficuum extracts also inhibited the urease enzyme compared to standard thiourea. The antispasmodic activity of A. ficuum extract was assessed by an in vivo model and the results demonstrated promising activity at 150 mg.kg-1. Molecular docking results also supported the antifungal, anticancer, and antiurease potency of A. ficuum extract. Overall, the results display promising aspects of A. ficuum extract as a future pharmacological source.

Surface enhanced Raman scattering investigation of tecovirimat on silver, gold and platinum loaded silica nanocomposites: Theoretical analysis (DFT) and molecular modeling

As of today, there have been 612 million confirmed cases of coronavirus disease (COVID-19) around the world, with over 6 million fatalities. Tecovirimat (TPOXX) is an anti-viral drug, and it was the first drug approved for the treatment of anti-pox virus in the US. However, the effectiveness of this drug against COVID-19 has not yet been explored. Since TPOXX is an anti-viral drug, an attempt has been made to determine its ability to act as a COVID inhibitor. Recent medical advances have resulted in the development of nano cage-based drug delivery. Drug delivery clusters based on nano cages have recently been used in the medical industry. As such, we used DFT coupled to the B3LYP/LANL2DZ basis set to study the adsorption behavior of the anti-viral drug TPOXX on Au/Ag/Pt⋯SiO2loaded silica nanocomposites. In order to identify the active site of the molecule, we have used the frontier molecular orbital (FMO) theory of molecular electrostatic potential (MEP). The compound and its complexes obey Lipinski's rule of five and have good drug-likeness properties based on the bioactivity evaluation. The biological properties of organic molecules and nano metal clusters were compared. The TPOXX with its nanocomposites was also studied in terms of Electron Localization Function (ELF) and Localized Orbital Locator (LOL). Molecular docking was performed for both pure molecule and its silica nanocomposites-doped derivatives with the chosen proteins to discuss the protein-ligand binding properties. These results could be more helpful in designing the drug and exploring its application for the inhibition of SARS-CoV-2.

Research in the Field of Drug Design and Development

The processes used by academic and industrial scientists to discover new drugs have recently experienced a true renaissance, with many new and exciting techniques being developed over the past 5-10 years alone. Drug design and discovery, and the search for new safe and well-tolerated compounds, as well as the ineffectiveness of existing therapies, and society's insufficient knowledge concerning the prophylactics and pharmacotherapy of the most common diseases today, comprise a serious challenge. This can influence not only the quality of human life, but also the health of whole societies, which became evident during the COVID-19 pandemic. In general, the process of drug development consists of three main stages: drug discovery, preclinical development using cell-based and animal models/tests, clinical trials on humans and, finally, forward moving toward the step of obtaining regulatory approval, in order to market the potential drug. In this review, we will attempt to outline the first three most important consecutive phases in drug design and development, based on the experience of three cooperating and complementary academic centers of the Visegrád group; i.e., Medical University of Lublin, Poland, Masaryk University of Brno, Czech Republic, and Comenius University Bratislava, Slovak Republic.

New Insight into the Pharmacological Importance of Atropine as the Potential Inhibitor of AKR1B1 via Detailed Computational Investigations: DFTs, ADMET, Molecular Docking, and Molecular Dynamics Studies

The aim of this research is to investigate the quantum geometric properties and chemical reactivity of atropine, a pharmaceutically active tropane alkaloid. Using density functional theory (DFT) computations with the B3LYP/SVP functional theory basis set, the most stable geometry of atropine was determined. Additionally, a variety of energetic molecular parameters were calculated, such as the optimized energy, atomic charges, dipole moment, frontier molecular orbital energies, HOMO-LUMO energy gap, molecular electrostatic potential, chemical reactivity descriptors, and molecular polarizability. To determine atropine's inhibitory potential, molecular docking was used to analyze ligand interactions within the active pockets of aldo-keto reductase (AKR1B1 and AKR1B10). The results of these studies showed that atropine has greater inhibitory action against AKR1B1 than AKR1B10, which was further validated through molecular dynamic simulations by analyzing root mean square deviation (RMSD) and root mean square fluctuations (RMSF). The results of the molecular docking simulation were supplemented with simulation data, and the ADMET characteristics were also determined to predict the drug likeness of a potential compound. In conclusion, the research suggests that atropine has potential as an inhibitor of AKR1B1 and could be used as a parent compound for the synthesis of more potent leads for the treatment of colon cancer associated with the sudden expression of AKR1B1.

Galactoside-Based Molecule Enhanced Antimicrobial Activity through Acyl Moiety Incorporation: Synthesis and In Silico Exploration for Therapeutic Target

In this study, a series of galactoside-based molecules, compounds of methyl β-d-galactopyranoside (MDGP, 1), were selectively acylated using 2-bromobenzoyl chloride to obtain 6-O-(2-bromobenzoyl) substitution products, which were then transformed into 2,3,4-tri-O-6-(2-bromobenzoyl) compounds (2-7) with various nontraditional acyl substituents. The chemical structures of the synthesized analogs were characterized by spectroscopic methods and physicochemical and elemental data analyses. The antimicrobial activities of the compounds against five human pathogenic bacteria and two phyto-fungi were evaluated in vitro and it was found that the acyl moiety-induced synthesized analogs exhibited varying levels of antibacterial activity against different bacteria, with compounds 3 and 6 exhibiting broad-spectrum activity and compounds 2 and 5 exhibiting activity against specific bacteria. Compounds 3 and 6 were tested for MIC (minimum inhibitory concentration) and MBC (minimum bactericidal concentration) based on their activity. The synthesized analogs were also found to have potential as a source of new antibacterial agents, particularly against gram-positive bacteria. The antifungal results suggested that the synthesized analogs could be a potential source of novel antifungal agents. Moreover, cytotoxicity testing revealed that the compounds are less toxic. A structure-activity relationship (SAR) investigation revealed that the lauroyl chain [CH₃(CH₂)₁₀CO-] and the halo-aromatic chain [3(/4)-Cl.C₆H₄CO-] in combination with sugar, had the most potent activity against bacterial and fungal pathogens. Density functional theory (DFT)-calculated thermodynamic and physicochemical parameters, and molecular docking, showed that the synthesized molecule may block dengue virus 1 NS2B/NS3 protease (3L6P). A 150 ns molecular dynamic simulation indicated stable conformation and binding patterns in a stimulating environment. In silico ADMET calculations suggested that the designed (MDGP, 1) had good drug-likeness values. In summary, the newly synthesized MDGP analogs exhibit potential antiviral activity and could serve as a therapeutic target for dengue virus 1 NS2B/NS3 protease.

Studying the Mechanism of Interaction of Doxofylline with Human Lysozyme: A Biophysical and In Silico Approach

In this study, multiple spectroscopic and computational methods were utilized to investigate the binding mechanism of doxofylline with lysozyme. The in vitro methods were used to obtain the binding kinetics and thermodynamics. UV-vis spectroscopy indicated the formation of complex between doxofylline and lysozyme. The Gibb's free energy and binding constant from UV-vis data was obtained as -7.20 kcal M^-1 and 1.929 × 10⁵ M^-1, respectively. Doxofylline successfully quenched the fluorescence of lysozyme, confirming the formation of complex. The k_q and K_sv values for the quenching of lysozyme's fluorescence by doxofylline were 5.74 × 10¹¹ M^-1 s^-1 and 3.32 × 10³ M^-1, respectively. These values signified a moderate binding affinity between doxofylline and lysozyme. In synchronous spectroscopy, red shifts were observed for indicating the changes in microenvironment of lysozyme following the binding of doxofylline. The secondary structural analysis was determined using circular dichroism (CD) which revealed an increase in % α-helical as a result of doxofylline interaction. The binding affinity and flexibility of lysozyme upon complexation have been revealed via molecular docking and molecular dynamic (MD) simulations, respectively. According to the many parameters of the MD simulation, the lysozyme-doxofylline complex was stable under physiological conditions. All during the simulation time, hydrogen bonds were continuously present. The MM-PBSA binding energy for lysozyme and doxofylline binding was found to be -30.55 kcal mol^-1.

Calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2) inhibitors: a novel approach in small molecule discovery

The calcium/calmodulin dependent protein kinase kinase 2 (CAMKK2) plays a key role in regulation of intracellular calcium levels and signaling pathways. It is involved in activation of downstream signaling pathways that regulate various cellular processes. Dysregulation of CAMKK2 activity has been linked to various diseases including cancer, suggesting that CAMKK2 inhibitors might be beneficial in oncological, metabolic and inflammatory indications. The most pressing issues in small molecule discovery are synthesis feasibility, novel chemical structure and desired biological characteristics. To circumvent this constraint, we employed 'DrugspaceX' for rapid lead identification, followed by repositioning seven FDA-approved drugs for CAMKK2 inhibition. Further, first-level transformation (Set1 analogues) was performed in 'DrugspaceX', followed by virtual screening. The t-SNE visualization revealed that the transformations surrounding Rucaparib, Treprostinil and Canagliflozin are more promising for developing CAMKK2 inhibitors. Second, using the top-ranked Set1 analogues, Set2 analogues were generated, and virtual screening revealed the top-ranked five analogues. Among the top five Set2 analogues, DE273038_5 had the lowest docking score of -11.034 kcal/mol and SA score of 2.59, retaining the essential interactions with Hotspot residues LYS194 and VAL270 across 250 ns simulation period. When compared to the other four compounds, the ligand effectiveness score was 0.409, and the number of rotatable penalties was only three. Further, DE273038_5 after two rounds of transformations was discovered to be novel and had not been previously described in other databases. These data suggest that the new candidate DE273038_5 is likely to have inhibitory activity at the CAMKK2 active site, implying potential therapeutic use.Communicated by Ramaswamy H. Sarma.

Multidimensional virtual screening approaches combined with drug repurposing to identify potential covalent inhibitors of SARS-CoV-2 3CL protease

The outbreak of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused an unprecedented global pandemic, and new cases are still on the rise due to the absence of effective medicines. However, developing new drugs within a short time is extremely difficult. Repurposing the existing drugs provides a fast and effective strategy to identify promising inhibitors. Here we focus on the screening of drugs database for discovering potential covalent inhibitors that target 3-chymotrypsin-like protease (3CLpro), an essential enzyme mediating viral replication and transcription. Firstly, we constructed a receptor-ligand pharmacophore model and verified it through decoy set. The importance of pharmacophore features was evaluated by combining molecular dynamics simulation with interaction analyses. Then, covalent docking was used to perform further screening. According to docking score and Prime/Molecular Mechanics Generalized Born Surface Area (MM-GBSA) score, total ten compounds obtained good scores and successfully established covalent bonds with the catalytic Cys145 residue. They also formed favorable interactions with key residues in active sites and closely integrated with 3CLpro with binding modes similar to known 3CLpro inhibitor. Finally, the top four hits DB08732, DB04653, DB01871 and DB07299 were further subjected to 100 ns molecular dynamics (MD) simulation and MM-GBSA binding free energy calculations. The results suggest that the four candidates show good binding affinities for 3CLpro, which warrants further evaluation for their in-vitro/in-vivo activities. Overall, our research methods provide a valuable reference for discovering promising inhibitors against SARS-CoV-2 and help to fight against the epidemic.Communicated by Ramaswamy H. Sarma.

Cyclohexane-1,3-dione Derivatives as Future Therapeutic Agents for NSCLC: QSAR Modeling, In Silico ADME-Tox Properties, and Structure-Based Drug Designing Approach

The abnormal expression of the c-Met tyrosine kinase has been linked to the proliferation of several human cancer cell lines, including non-small-cell lung cancer (NSCLC). In this context, the identification of new c-Met inhibitors based on heterocyclic small molecules could pave the way for the development of a new cancer therapeutic pathway. Using multiple linear regression (MLR)-quantitative structure-activity relationship (QSAR) and artificial neural network (ANN)-QSAR modeling techniques, we look at the quantitative relationship between the biological inhibitory activity of 40 small molecules derived from cyclohexane-1,3-dione and their topological, physicochemical, and electronic properties against NSCLC cells. In this regard, screening methods based on QSAR modeling with density-functional theory (DFT) computations, in silico pharmacokinetic/pharmacodynamic (ADME-Tox) modeling, and molecular docking with molecular electrostatic potential (MEP) and molecular mechanics-generalized Born surface area (MM-GBSA) computations were used. Using physicochemical (stretch-bend, hydrogen bond acceptor, Connolly molecular area, polar surface area, total connectivity) and electronic (total energy, highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy levels) molecular descriptors, compound 6d is identified as the optimal scaffold for drug design based on in silico screening tests. The computer-aided modeling developed in this study allowed us to design, optimize, and screen a new class of 36 small molecules based on cyclohexane-1,3-dione as potential c-Met inhibitors against NSCLC cell growth. The in silico rational drug design approach used in this study led to the identification of nine lead compounds for NSCLC therapy via c-Met protein targeting. Finally, the findings are validated using a 100 ns series of molecular dynamics simulations in an aqueous environment on c-Met free and complexed with samples of the proposed lead compounds and Foretinib drug.

Identification of new dihydrophenanthrene derivatives as promising anti-SARS-CoV-2 drugs through in silico investigations

To research, evaluate, and invent novel compounds that inhibit SARS-CoV-2 activity, a series of reported 39 substituted 9, 10-dihydrophenanthrene derivatives were subjected to a quantitative structure-activity relationship (QSAR) study. Gaussian 09 and ChemOffice programs were used to calculate the molecular descriptors employed to determine their impact on the studied activity. Then we reduced the number of descriptors by eliminating the redundant information using principal component analysis (PCA). The creation of molecular models was done by using multiple linear regression (MLR) according to the principles established by the Organization for Economic Co-operation and Development (OECD) and the validation by using external and internal validation, Y-randomization tests, and domain of applicability. Moreover, we evaluated the toxicity of developed compounds using ADMET and Molecular docking to determine their optimal position to form a stable complex. As a result, four molecules may be used to develop a novel drug that can inhibit SARS-CoV-2 without causing the side effect.

An in-silico pharmacophore-based molecular docking study to evaluate the inhibitory potentials of novel fungal triterpenoid Astrakurkurone analogues against a hypothetical mutated main protease of SARS-CoV-2 virus

BACKGROUND

The main protease is an important structural protein of SARS-CoV-2, essential for its survivability inside a human host. Considering current vaccines' limitations and the absence of approved therapeutic targets, Mpro may be regarded as the potential candidate drug target. Novel fungal phytocompound Astrakurkurone may be studied as the potential Mpro inhibitor, considering its medicinal properties reported elsewhere.

METHODS

In silico molecular docking was performed with Astrakurkurone and its twenty pharmacophore-based analogues against the native Mpro protein. A hypothetical Mpro was also constructed with seven mutations and targeted by Astrakurkurone and its analogues. Furthermore, multiple parameters such as statistical analysis (Principal Component Analysis), pharmacophore alignment, and drug likeness evaluation were performed to understand the mechanism of protein-ligand molecular interaction. Finally, molecular dynamic simulation was done for the top-ranking ligands to validate the result.

RESULT

We identified twenty Astrakurkurone analogues through pharmacophore screening methodology. Among these twenty compounds, two analogues namely, ZINC89341287 and ZINC12128321 showed the highest inhibitory potentials against native and our hypothetical mutant Mpro, respectively (-7.7 and -7.3 kcal mol-1) when compared with the control drug Telaprevir (-5.9 and -6.0 kcal mol-1). Finally, we observed that functional groups of ligands namely two aromatic and one acceptor groups were responsible for the residual interaction with the target proteins. The molecular dynamic simulation further revealed that these compounds could make a stable complex with their respective protein targets in the near-native physiological condition.

CONCLUSION

To conclude, Astrakurkurone analogues ZINC89341287 and ZINC12128321 can be potential therapeutic agents against the highly infectious SARS-CoV-2 virus.

Repurposing of Chemotherapeutics to Combat COVID-19

Severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) is a novel strain of SARS coronavirus. The COVID-19 disease caused by this virus was declared a pandemic by the World Health Organization (WHO). SARS-CoV-2 mainly spreads through droplets sprayed by coughs or sneezes of the infected to a healthy person within the vicinity of 6 feet. It also spreads through asymptomatic carriers and has negative impact on the global economy, security and lives of people since 2019. Numerous lives have been lost to this viral infection; hence there is an emergency to build up a potent measure to combat SARS-CoV-2. In view of the non-availability of any drugs or vaccines at the time of its eruption, the existing antivirals, antibacterials, antimalarials, mucolytic agents and antipyretic paracetamol were used to treat the COVID-19 patients. Still there are no specific small molecule chemotherapeutics available to combat COVID-19 except for a few vaccines approved for emergency use only. Thus, the repurposing of chemotherapeutics with the potential to treat COVID-19 infected people is being used. The antiviral activity for COVID-19 and biochemical mechanisms of the repurposed drugs are being explored by the biological assay screening and structure-based in silico docking simulations. The present study describes the various US-FDA approved chemotherapeutics repositioned to combat COVID-19 along with their screening for biological activity, pharmacokinetic and pharmacodynamic evaluation.

Rational identification of small molecules derived from 9,10-dihydrophenanthrene as potential inhibitors of 3CLpro enzyme for COVID-19 therapy: a computer-aided drug design approach

Small molecules such as 9,10-dihydrophenanthrene derivatives have remarkable activity toward inhibition of SARS-CoV-2 3CLpro and COVID-19 proliferation, which show a strong correlation between their structures and bioactivity. Therefore, these small compounds could be suitable for clinical pharmaceutical use against COVID-19. The objective of this study was to remodel the structures of 9,10-dihydrophenanthrene derivatives to achieve a powerful biological activity against 3CLpro and favorable pharmacokinetic properties for drug design and discovery. Therefore, by the use of bioinformatics techniques, we developed robust 3D-QSAR models that are capable of describing the structure-activity relationship for 46 molecules based on 9,10-dihydrophenanthrene derivatives using CoMFA/SE (R 2 = 0.97, Q 2 = 0.81, R 2 pred = 0.95, c R 2 p = 0.71) and CoMSIA/SEHDA (R 2 = 0.94, Q 2 = 0.76, R 2 pred = 0.91, c R 2 p = 0.65) techniques. Accordingly, 96 lead compounds were generated based on a template molecule that showed the highest observed activity in vitro (T40, pIC50 = 5.81) and predicted their activities and bioavailability in silico. The rational screening outputs of 3D-QSAR, Molecular docking, ADMET, and MM-GBSA led to the identification of 9 novel modeled molecules as potent noncovalent drugs against SARS-CoV-2-3CLpro. Finally, by molecular dynamics simulations, the stability and structural dynamics of 3CLpro free and complex (PDB code: 6LU7) were discussed in the presence of samples of 9,10-dihydrophenanthrene derivative in an aqueous environment. Overall, the retrosynthesis of the proposed drug compounds in this study and the evaluation of their bioactivity in vitro and in vivo may be interesting for designing and discovering a new drug effective against COVID-19.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11224-022-02004-z.

Riboflavin, a Potent Neuroprotective Vitamin: Focus on Flavivirus and Alphavirus Proteases

Several neurotropic viruses are members of the flavivirus and alphavirus families. Infections caused by these viruses may cause long-term neurological sequelae in humans. The continuous emergence of infections caused by viruses around the world, such as the chikungunya virus (CHIKV) (Alphavirus genus), the zika virus (ZIKV) and the yellow fever virus (YFV) (both of the Flavivirus genus), warrants the development of new strategies to combat them. Our study demonstrates the inhibitory potential of the water-soluble vitamin riboflavin against NS2B/NS3pro of ZIKV and YFV and nsP2pro of CHIKV. Riboflavin presents a competitive inhibition mode with IC50 values in the medium µM range of 79.4 ± 5.0 µM for ZIKV NS2B/NS3pro and 45.7 ± 2.9 μM for YFV NS2B/NS3pro. Against CHIKV nsP2pro, the vitamin showed a very strong effect (93 ± 5.7 nM). The determined dissociation constants (KD) are significantly below the threshold value of 30 µM. The ligand binding increases the thermal stability between 4 °C and 8 °C. Unexpectedly, riboflavin showed inhibiting activity against another viral protein; the molecule was also able to inhibit the viral entry of CHIKV. Molecular dynamics simulations indicated great stability of riboflavin in the protease active site, which validates the repurposing of riboflavin as a promising molecule in drug development against the viruses presented here.

Computational study of quinoline-based thiadiazole compounds as potential antileishmanial inhibitors

Leishmaniasis is a severe disease caused by protozoan parasites of the genus Leishmania and it is accountable for sizable morbidity and mortality worldwide.