Molecular Mechanics Generalized Born Surface Area Research Articles

Isolation, Virtual Screening, and Evaluation of Hazelnut-Derived Immunoactive Peptides for the Inhibition of SARS-CoV-2 Main Protease.

The aim of this study is to validate the activity of hazelnut (Corylus avellana L.)-derived immunoactive peptides inhibiting the main protease (Mpro) of SARS-CoV-2 and further unveil their interaction mechanism using in vitro assays, molecular dynamics (MD) simulations, and binding free energy calculations. In general, the enzymatic hydrolysis components, especially molecular weight < 3 kDa, possess good immune activity as measured by the proliferation ability of mouse splenic lymphocytes and phagocytic activity of mouse peritoneal macrophages. Over 866 unique peptide sequences were isolated, purified, and then identified by nanohigh-performance liquid chromatography/tandem mass spectrometry (NANO-HPLC-MS/MS) from hazelnut protein hydrolysates, but Trp-Trp-Asn-Leu-Asn (WWNLN) and Trp-Ala-Val-Leu-Lys (WAVLK) in particular are found to increase the cell viability and phagocytic capacity of RAW264.7 macrophages as well as promote the secretion of the cytokines nitric oxide (NO), tumor necrosis factor-α (TNF-α), and interleukin-1β (IL-1β). Fluorescence resonance energy transfer assay elucidated that WWNLN and WAVLK exhibit excellent inhibitory potency against Mpro, with IC50 values of 6.695 and 16.750 μM, respectively. Classical all-atom MD simulations show that hydrogen bonds play a pivotal role in stabilizing the complex conformation and protein-peptide interaction. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculation indicates that WWNLN has a lower binding free energy with Mpro than WAVLK. Furthermore, adsorption, distribution, metabolism, excretion, and toxicity (ADMET) predictions illustrate favorable drug-likeness and pharmacokinetic properties of WWNLN compared to WAVLK. This study provides a new understanding of the immunomodulatory activity of hazelnut hydrolysates and sheds light on peptide inhibitors targeting Mpro.

Cheminformatics approach to identify andrographolide derivatives as dual inhibitors of methyltransferases (nsp14 and nsp16) of SARS-CoV-2

The Covid-19 pandemic outbreak has accelerated tremendous efforts to discover a therapeutic strategy that targets severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to control viral infection. Various viral proteins have been identified as potential drug targets, however, to date, no specific therapeutic cure is available against the SARS-CoV-2. To address this issue, the present work reports a systematic cheminformatic approach to identify the potent andrographolide derivatives that can target methyltransferases of SARS-CoV-2, i.e. nsp14 and nsp16 which are crucial for the replication of the virus and host immune evasion. A consensus of cheminformatics methodologies including virtual screening, molecular docking, ADMET profiling, molecular dynamics simulations, free-energy landscape analysis, molecular mechanics generalized born surface area (MM-GBSA), and density functional theory (DFT) was utilized. Our study reveals two new andrographolide derivatives (PubChem CID: 2734589 and 138968421) as natural bioactive molecules that can form stable complexes with both proteins via hydrophobic interactions, hydrogen bonds and electrostatic interactions. The toxicity analysis predicts class four toxicity for both compounds with LD50 value in the range of 500–700 mg/kg. MD simulation reveals the stable formation of the complex for both the compounds and their average trajectory values were found to be lower than the control inhibitor and protein alone. MMGBSA analysis corroborates the MD simulation result and showed the lowest energy for the compounds 2734589 and 138968421. The DFT and MEP analysis also predicts the better reactivity and stability of both the hit compounds. Overall, both andrographolide derivatives exhibit good potential as potent inhibitors for both nsp14 and nsp16 proteins, however, in-vitro and in vivo assessment would be required to prove their efficacy and safety in clinical settings. Moreover, the drug discovery strategy aiming at the dual target approach might serve as a useful model for inventing novel drug molecules for various other diseases.

Identification of small molecule inhibitors against Lin28/let-7 to suppress tumor progression and its alleviation role in LIN28-dependent glucose metabolism.

The reciprocal suppression of an RNA-binding protein LIN28 (human abnormal cell lineage 28) and miRNA Let-7 (Lethal 7) is considered to have a prime role in hepatocellular carcinoma (HCC). Though targeting this inhibition interaction is effective for therapeutics, it causes other unfavorable effects on glucose metabolism and increased insulin resistance. Hence, this study aims to identify small molecules targeting Lin28/let-7 interaction along with additional potency to improve insulin sensitivity. Of 22,14,996 small molecules screened by high throughput virtual screening, 6 molecules, namely 41354, 1558, 12437, 23837, 15710, and 8319 were able to block the LIN28 interaction with let-7 and increase the insulin sensitivity via interacting with PPARγ (peroxisome proliferator-activated receptors γ). MM-GBSA (Molecular Mechanics-Generalized Born Surface Area) analysis is used to re-score the binding affinity of docked complexes. Upon further analysis, it is also seen that these molecules have superior ADME (Absorption, Distribution, Metabolism, and Excretion) properties and form stable complexes with the targets for a significant period in a biologically simulated environment (Molecular Dynamics simulation) for 100 ns. From our results, we hypothesize that these identified 6 small molecules can be potential candidates for HCC treatment and the glucose metabolic disorder caused by the HCC treatment.

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

Background The coronavirus disease (COVID-19) pandemic continues to be a public health problem worldwide. Several therapeutic targets of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been identified, whereas the main protease (Mpro) is necessary for virus replication. Since SARS- CoV-2 Mpro mutation rates are inherently high, searching for new inhibitors remains challenging. Herein, this work aimed to evaluate 84 natural compounds from Peruvian sources against different mutations on the Mpro target. Methods We applied virtual screening, all-atom molecular dynamics simulations, and binding free energy estimation by Molecular Mechanics/Generalized Born Surface Area (MM/GBSA). Results The virtual screening results helped us identify rutin as the top compound against different Mpro mutations. Likewise, the computational simulations demonstrated the high structural stability of the Mpro-rutin system. Conclusions his research evaluated the antiviral capacity of Peruvian sources against SARS-CoV-2 Mpro and its mutations, which could be important in preventing and treating SARS-CoV-2 infection.

Membrane-assisted tariquidar access and binding mechanisms of human ATP-binding cassette transporter P-glycoprotein.

The human multidrug transporter P-glycoprotein (P-gp) is physiologically essential and of key relevance to biomedicine. Recent structural studies have shed light on the mode of inhibition of the third-generation inhibitors for human P-gp, but the molecular mechanism by which these inhibitors enter the transmembrane sites remains poorly understood. In this study, we utilized all-atom molecular dynamics (MD) simulations to characterize human P-gp dynamics under a potent inhibitor, tariquidar, bound condition, as well as the atomic-level binding pathways in an explicit membrane/water environment. Extensive unbiased simulations show that human P-gp remains relatively stable in tariquidar-free and bound states, while exhibiting a high dynamic binding mode at either the drug-binding pocket or the regulatory site. Free energy estimations by partial nudged elastic band (PNEB) simulations and Molecular Mechanics Generalized Born Surface Area (MM/GBSA) method identify two energetically favorable binding pathways originating from the cytoplasmic gate with an extended tariquidar conformation. Interestingly, free tariquidar in the lipid membrane predominantly adopts extended conformations similar to those observed at the regulatory site. These results suggest that membrane lipids may preconfigure tariquidar into an active ligand conformation for efficient binding to the regulatory site. However, due to its conformational plasticity, tariquidar ultimately moves toward the drug-binding pocket in both pathways, explaining how it acts as a substrate at low concentrations. Our molecular findings propose a membrane-assisted mechanism for the access and binding of the third-generation inhibitors to the binding sites of human P-gp, and offer deeper insights into the molecule design of more potent inhibitors against P-gp-mediated drug resistance.

Discovering potential WRN inhibitors from natural product database through computational methods

Microsatellite instability (MSI) is a relatively common feature associated with multiple cancers, and Werner syndrome (WRN) ATP-dependent helicase has been recognized as a novel target for treating MSI cancers, such as colorectal cancer. A small-molecule inhibitor targeting WRN would be a promising strategy for treating colorectal cancer with high MSI expression. In this study, we employed a computer-assisted drug discovery strategy to screen over 30,000 natural product molecules. By using a combination of docking, ligand efficiency, Molecular Mechanics/Generalized Born Surface Area (MM/GBSA), and thermodynamic integration (TI) calculations, we identified MOL008980, MOL010740, MOL011832, T4743, TN1166, and TNP-002173 as potential WRN inhibitors. Subsequent molecular dynamics simulation revealed that these screened natural products possessed better binding dynamic characteristics than ATP substrate and were capable of inhibiting the dynamic process of WRN, making them potential strong ATP competitive inhibitors. In conclusion, our computational approach revealed potential WRN inhibitors from a natural product database, providing a theoretical basis for future research.

Integrating virtual screening, pharmacoinformatics profiling, and molecular dynamics: identification of promising inhibitors targeting 3CLpro of SARS-CoV-2.

Introduction: The pursuit of effective therapeutic solutions for SARS-CoV-2 infections and COVID-19 necessitates the repurposing of existing compounds. This study focuses on the detailed examination of the central protease, 3-chymotrypsin-like protease (3CLpro), a pivotal player in virus replication. The combined approach of molecular dynamics simulations and virtual screening is employed to identify potential inhibitors targeting 3CLpro. Methods: A comprehensive virtual screening of 7120 compounds sourced from diverse databases was conducted. Four promising inhibitors, namely EN1036, F6548-4084, F6548-1613, and PUBT44123754, were identified. These compounds exhibited notable attributes, including high binding affinity (ranging from -5.003 to -5.772Kcal/mol) and superior Induced Fit Docking scores (ranging from -671.66 to -675.26Kcal/mol) compared to co-crystallized ligands. Results: In-depth analysis revealed that F6548-1613 stood out, demonstrating stable hydrogen bonds with amino acids His41 and Thr62. Notably, F6548-1613 recorded a binding energy of -65.72kcal/mol in Molecular Mechanics Generalized Born Surface Area (MMGBSA) simulations. These findings were supported by Molecular Dynamics simulations, highlighting the compound's efficacy in inhibiting 3CLpro. Discussion: The identified compounds, in compliance with Lipinski's rule of five and exhibiting functional molecular interactions with 3CLpro, present promising therapeutic prospects. The integration of in silico methodologies significantly expedites drug discovery, laying the foundation for subsequent experimental validation and optimization. This approach holds the potential to develop effective therapeutics for SARS-CoV-2.

Design and In Silico and In Vitro Evaluations of a Novel Nicotinamide Derivative as a VEGFR-2 Inhibitor

A new nicotinamide derivative, (E)-N-(4-(1-(2-(4-benzamidobenzoyl)hydrazone)ethyl)phenyl)nicotinamide, was designed as a VEGFR-2 inhibitor. Utilizing the density functional theory (DFT) calculations, the three-dimensional structure of the designed compound was determined, shedding light on its stability and reactivity. Molecular docking revealed its capability to inhibit VEGFR-2, which was further supported by molecular dynamics (MD) simulations confirming its binding to the target protein. In addition, molecular mechanics-generalized born surface area (MM-GBSA), protein-ligand interactions profiler (PLIP), and essential dynamics studies provided further validation of the compound’s precise binding with optimal energy. Then, the “compound 10” was synthesized and subjected to in vitro assays. Compound 10 inhibited VEGFR-2 with an IC50 value of 105.4 ± 0.896 nM, comparing sorafenib’s IC50 value of 61.65 ± 0.934 nM. Besides, it exhibited cytotoxicity against HepG2 and MCF-7 cancer cell lines, with IC50 values of 35.78 ± 0.863 μM and 57.62 μM ± 0.871, comparing sorafenib’s IC50 values of 5.95 ± 0.917 μM and 8.45 ± 0.912 μM. Furthermore, compound 10 demonstrated a lower level of toxicity towards Vero cell lines, with an IC50 value of 127.3 μM. Likewise, compound 10 induced apoptosis in HepG2 cell lines through a flow cytometric analysis in addition to an increase in the levels of caspase-3 and caspase-9. Moreover, compound 10 hindered the migration and healing abilities of HepG2 cells. In conclusion, our study positions compound 10 as a promising candidate for further chemical modifications and biological evaluations.

Proteomics-based vaccine targets annotation and design of multi-epitope vaccine against antibiotic-resistant Streptococcus gallolyticus

Streptococcus gallolyticus is a non-motile, gram-positive bacterium that causes infective endocarditis. S. gallolyticus has developed resistance to existing antibiotics, and no vaccine is currently available. Therefore, it is essential to develop an effective S. gallolyticus vaccine. Core proteomics was used in this study together with subtractive proteomics and reverse vaccinology approach to find antigenic proteins that could be utilized for the design of the S. gallolyticus multi-epitope vaccine. The pipeline identified two antigenic proteins as potential vaccine targets: penicillin-binding protein and the ATP synthase subunit. T and B cell epitopes from the specific proteins were forecasted employing several immunoinformatics and bioinformatics resources. A vaccine (360 amino acids) was created using a combination of seven cytotoxic T cell lymphocyte (CTL), three helper T cell lymphocyte (HTL), and five linear B cell lymphocyte (LBL) epitopes. To increase immune responses, the vaccine was paired with a cholera enterotoxin subunit B (CTB) adjuvant. The developed vaccine was highly antigenic, non-allergenic, and stable for human use. The vaccine's binding affinity and molecular interactions with the human immunological receptor TLR4 were studied using molecular mechanics/generalized Born surface area (MMGBSA), molecular docking, and molecular dynamic (MD) simulation analyses. Escherichia coli (strain K12) plasmid vector pET-28a ( +) was used to examine the ability of the vaccine to be expressed. According to the outcomes of these computer experiments, the vaccine is quite promising in terms of developing a protective immunity against diseases. However, in vitro and animal research are required to validate our findings.

Tepotinib and tivantinib as potential inhibitors for the serine/threonine kinase of the mpox virus: insights from structural bioinformatics analysis

The serine/threonine kinase (STK) plays a central role as the primary kinase in poxviruses, directing phosphoryl transfer reactions. Such reactions are pivotal for the activation of certain proteins during viral replication, assembly, and maturation. Therefore, targeting this key protein is anticipated to impede virus replication. In this work, a structural bioinformatics approach was employed to evaluate the potential of drug-like kinase inhibitors in binding to the ATP-binding pocket on the STK of the Mpox virus. Virtual screening of known kinase inhibitors revealed that the top 10 inhibitors exhibited binding affinities ranging from −8.59 to −12.05 kcal/mol. The rescoring of compounds using the deep-learning default model in GNINA was performed to predict accurate binding poses. Subsequently, the top three inhibitors underwent unbiased molecular dynamics (MD) simulations for 100 ns. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) analysis and Principal Component Analysis (PCA) suggested tepotinib as a competitive inhibitor for Mpox virus STK as evidenced by its binding free energy and the induction of similar conformational behavior of the enzyme. Nevertheless, it is sensible to experimentally test all top 10 compounds, as scoring functions and energy calculations may not consistently align with experimental findings. These insights are poised to provide an attempt to identify an effective inhibitor for the Mpox virus. Communicated by Ramaswamy H. Sarma

Harnessing Brazilian biodiversity database: identification of flavonoids as potential inhibitors of SARS-CoV-2 main protease using computational approaches and all-atom molecular dynamics simulation.

SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) is the etiological agent responsible for the global outbreak of COVID-19 (Coronavirus Disease 2019). The main protease of SARS-CoV-2, Mpro, is a key enzyme that plays a vital role in mediating viral replication and transcription. In this study, a comprehensive computational approach was employed to investigate the binding affinity, selectivity, and stability of natural product candidates as potential new antivirals acting on the viral polyprotein processing mediated by SARS-CoV-2 Mpro. A library of 288 flavonoids extracted from Brazilian biodiversity was screened to select potential Mpro inhibitors. An initial filter based on Lipinski's rule of five was applied, and 204 compounds that did not violate any of the Lipinski rules were selected. The compounds were then docked into the active site of Mpro using the GOLD program, and the poses were subsequently re-scored using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) binding free energy calculations performed by AmberTools23. The top five flavonoids with the best MM-GBSA binding free energy values were selected for analysis of their interactions with the active site residues of the protein. Next, we conducted a toxicity and drug-likeness analysis, and non-toxic compounds were subjected to molecular dynamics simulation and free energy calculation using the MM-PBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) method. It was observed that the five selected flavonoids had lower MM-GBSA binding free energy with Mpro than the co-crystal ligand. Furthermore, these compounds also formed hydrogen bonds with two important residues, Cys145 and Glu166, in the active site of Mpro. Two compounds that passed the drug-likeness filter showed stable conformations during the molecular dynamics simulations. Among these, NuBBE_867 exhibited the best MM-PBSA binding free energy value compared to the crystallographic inhibitor. Therefore, this study suggests that NuBBE_867 could be a potential inhibitor against the main protease of SARS-CoV-2 and may be further examined to confirm our results.

Potential anti-obesity effect of saponin metabolites from adzuki beans: A computational approach.

In contrast to its widespread traditional and popular culinary use to reduce weight, Vigna angularis (adzuki beans) was not subjected to sufficient scientific scrutiny. Particularly, its saponins whose role was never investigated before to unveil the beans' antidiabetic and anti-obesity effects. Four vital pancreatic and intestinal carbohydrate enzymes were selected to assess the potency of the triterpenoidal saponins of V. angularis to bind and activate these proteins through high-precision molecular modeling and dynamics mechanisms with accurate molecular mechanics Generalized Born Surface Area (MMGBSA) energy calculations; thus, recognizing their anti-obesity potential. Our results showed that adzukisaponin VI and adzukisaponin IV were the best compounds in the α-amylase and α-glucosidase enzymatic grooves, respectively. Adzukisaponin VI and angulasaponin C were the best fitting in the N-termini of sucrase-isomaltose (SI) enzyme, and angulasaponin C was the best scoring compound in maltase-glucoamylase C-termini. All of them outperformed the standard drug acarbose. These compounds in their protein complexes were selected to undergo molecular simulations of the drug-bound protein compared to the apo-protein through 100 ns, which confirmed the consistency of binding to the key amino acid residues in the four enzyme pockets with the least propensity of unfolding. Detailed analysis is given of the different polar and hydrophobic binding interactions of docked compounds. While maltase-adzukisaponin VI complex scored the lowest MMGBSA free energy of -67.77 Kcal/mol, α-amylase complex with angulasaponin B revealed the free binding energy of -74.18 Kcal/mol with a dominance of van der Waals energy (ΔEVDW) and the least change from the start to the end of the simulation time. This study will direct researchers to the significance of isolating the pure adzuki saponin components to conduct future invitro and invivo experimental works and even clinical trials.

Network pharmacology, molecular simulation, and binding free energy calculation-based investigation of Neosetophomone B revealed key targets for the treatment of cancer.

In the current study, Neosetophomone B (NSP-B) was investigated for its anti-cancerous potential using network pharmacology, quantum polarized ligand docking, molecular simulation, and binding free energy calculation. Using SwissTarget prediction, and Superpred, the molecular targets for NSP-B were predicted while cancer-associated genes were obtained from DisGeNet. Among the total predicted proteins, only 25 were reported to overlap with the disease-associated genes. A protein-protein interaction network was constructed by using Cytoscape and STRING databases. MCODE was used to detect the densely connected subnetworks which revealed three sub-clusters. Cytohubba predicted four targets, i.e., fibroblast growth factor , FGF20, FGF22, and FGF23 as hub genes. Molecular docking of NSP-B based on a quantum-polarized docking approach with FGF6, FGF20, FGF22, and FGF23 revealed stronger interactions with the key hotspot residues. Moreover, molecular simulation revealed a stable dynamic behavior, good structural packing, and residues' flexibility of each complex. Hydrogen bonding in each complex was also observed to be above the minimum. In addition, the binding free energy was calculated using the MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) and MM/PBSA (Molecular Mechanics/Poisson-Boltzmann Surface Area) approaches. The total binding free energy calculated using the MM/GBSA approach revealed values of -36.85kcal/mol for the FGF6-NSP-B complex, -43.87kcal/mol for the FGF20-NSP-B complex, and -37.42kcal/mol for the FGF22-NSP-B complex, and -41.91kcal/mol for the FGF23-NSP-B complex. The total binding free energy calculated using the MM/PBSA approach showed values of -30.05kcal/mol for the FGF6-NSP-B complex, -39.62kcal/mol for the FGF20-NSP-B complex, -34.89kcal/mol for the FGF22-NSP-B complex, and -37.18kcal/mol for the FGF23-NSP-B complex. These findings underscore the promising potential of NSP-B against FGF6, FGF20, FGF22, and FGF23, which are reported to be essential for cancer signaling. These results significantly bolster the potential of NSP-B as a promising candidate for cancer therapy.

PHARMACOINFORMATICS ANALYSIS OF MORUS MACROURA FOR DRUG DISCOVERY AND DEVELOPMENT

Objective: Pharmacoinformatics is an innovative approach rapidly evolving in pharmaceutical research and drug development. This study focuses on analysing Morus macroura, a plant species with untapped pharmacological potential. This investigation aims to leverage pharmacoinformatics techniques to unveil the hidden potential of Morus macroura in drug discovery and development.&#x0D; Methods: The study includes analyses of protein-protein interactions, deep learning docking, adsorption tests, distribution, metabolism, excretion, molecular dynamics simulations and free energy calculation using Molecular Mechanics Generalized Born Surface Area (MMGBSA).&#x0D; Results: Nine active compounds were identified in Morus macroura, namely Andalasin A, Guangsangon K, Guangsangon L, Guangsangon M, Guangsangon N, Macrourone C, Mulberrofuran G, Mulberrofuran K, and Mulberroside C. These compounds exhibit protein-protein interaction activities against a cytochrome P450 monooxygenase that catalyses the conversion of C19 androgens. These plant compounds influence aromatase excess syndrome, deficiency, and ovarian dysgenesis. Regarding drug-likeness, Mulberroside C and Macrourone C demonstrated good absorption potential by adhering to Lipinski's rule of five. Deep learning docking simulations yielded affinity results of-9.62 kcal/mol for Guangsangon M,-10.44 kcal/mol for Macrourone C, and-10.99 kcal/mol for Guangsangon L. Subsequent molecular dynamics simulations indicated that Guangsangon L and Macrourone C remained stable during a 100 ns simulation.&#x0D; Conclusion: Morus macroura interacts with important proteins, particularly CYP19A1, which might influence health conditions like aromatase excess syndrome and ovarian dysgenesis. These findings provide potential paths for addressing specific health issues and advancing drug development. Molecular dynamics simulations indicated that Guangsangon L and Macrourone C remained stable during simulation.

The Interaction Mechanism of Picolinamide Fungicide Targeting on the Cytochrome bc1 Complex and Its Structural Modification.

Picolinamide fungicides, structurally related to UK-2A and antimycin-A, bind into the Qi-site in the bc1 complex. However, the detailed binding mode of picolinamide fungicides remains unknown. In the present study, antimycin-A and UK-2A were selected to study the binding mode of picolinamide inhibitors with four protonation states in the Qi-site by integrating molecular dynamics simulation, molecular docking, and molecular mechanics Generalized Born surface area (MM/GBSA) calculations. Subsequently, a series of new picolinamide derivatives were designed and synthesized to further understand the effects of substituents on the tail phenyl ring. The computational results indicated that the substituted aromatic rings in antimycin-A and UK-2A were the pharmacophore fragments and made the primary contribution when bound to a protein. Compound 9g-hydrolysis formed H-bonds with Hie201 and Ash228 and showed an IC50 value of 6.05 ± 0.24 μM against the porcine bc1 complex. Compound 9c, with a simpler chemical structure, showed higher control effects than florylpicoxamid against cucumber downy mildew and expanded the fungicidal spectrum of picolinamide fungicides. The structural and mechanistic insights obtained from the present study will provide a valuable clue for the future designing of new promising Qi-site inhibitors.

Bioactivity and molecular docking studies of selected plant compounds.

The shortcomings of synthetic pesticides are becoming increasingly evident. In the pursuit of looking for alternative, environmentally friendly pesticides, the potential use of two naturally occurring sesquiterpenes, α-bisabolol and bisabolol oxide A, as bactericides or acaricides, was investigated against three gram-positive and three gram-negative bacteria, and the eggs and larvae of the cattle tick Hyalomma scupense, respectively. The primary targeted site of synthetic insecticides is acetylcholinesterase (AChE), so the inhibitory activity of both sesquiterpenes was evaluated. Then, molecular docking and molecular mechanics generalized Born surface area (MM-GBSA)studies were performed to gain insight into their binding mode in the AChE active site. Bisabolol oxide A showed good activity against Staphylococcus aureus, Escherichia coli, and Salmonella enteritidis with inhibition zones ranging from 16 to 34 mm and minimal inhibitory concentration (CIM) of 3.125  mg mL-1 . Ovicidal and larvicidal assays revealed promising activity of α-bisabolol and bisabolol oxide A against H. scupense, with bisabolol oxide A being more effective against larval ticks with lethal concentration (LC50 ) value of 0.78%. Both compounds also displayed potent anti-AChE activity with inhibition concentration (IC50 ) values of 37.09 and 28.14 μg mL-1 , respectively. Furthermore, α-bisabolol and bisabolol oxide A exhibited good and comparable docking scores (-7.289 and -7.139 Kcal mol-1 , respectively) and were found to accommodate in the active-site gorge of AChE via hydrogen bonding and hydrophobic interactions. Bisabolol oxide A and, to a lesser degree, α-bisabolol are active against bacteria and ectoparasites and may represent an economical and sustainable alternative to toxic synthetic pesticides to control pathogens. © 2024 Society of Chemical Industry.

The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 − CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 − CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation. Communicated by Ramaswamy H. Sarma

In silico and in vitro screening of pyrrole-based Hydrazide-Hydrazones as novel acetylcholinesterase inhibitors

Virtual screening is emerging as a highly applied technique and gained prominence as widely used method for the search and identification of potential hits, significantly reducing the time needed to discover novel and effective compounds compared to high-throughput screening. Recently, the superiority of simulations with multiple programs compared to a single software docking has been discussed. The aim of this work was to apply consensus docking, molecular mechanics/generalized Born surface area (MM/GBSA) free binding energy recalculations, and in vitro evaluations on an in-house dataset of recently synthesized pyrrole-based hydrazide-hydrazones in the search for novel acetylcholinesterase (AChE) inhibitors. Two licensed softwares – GOLD 5.3 and Glide, were employed for the virtual screenings, and several chemotherapeutic potential hits were identified. Furthermore, MM/GBSA free binding energy recalculations were provided to enhance the robustness of the in silico results. The MM/GBSA scores of the top ten pyrrole-based hydrazide-hydrazones were ranging from -60.44 to -70.93 Kcal/mol. Subsequent, in vitro evaluations of the top ranked compounds revealed that 12d exhibited the highest AChE inhibitory activity, with a 55% inhibition rate at a concentration of 10 μM. Moreover, this prominent pyrrole-based AChE inhibitor formed stable complex with the active site of the enzyme. Interactions with the active amino residues Tyr72 and Tyr286 indicated that 12d was located near the peripheral anionic site of the enzyme. Additionally, in silicoADME investigations using QikProp demonstrated that 12d possesses optimal pharmacokinetic properties. In conclusion, this study identified a novel pyrrole-based AChE inhibitor 12d through a combination of computational and experimental findings.

Binding behavior of receptor binding domain of the SARS-CoV-2 virus and ivermectin

The COVID-19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), sparked an international debate on effective ways to prevent and treat the virus. Specifically, there were many varying opinions on the use of ivermectin (IVM) throughout the world, with minimal research to support either side. IVM is an FDA-approved antiparasitic drug that was discovered in the 1970s and was found to show antiviral activity. The objective of this study is to examine the binding behavior and rates of association and dissociation between SARS-CoV-2 receptor binding domain (RBD), IVM, and their combination using aminopropylsilane (APS) biosensors as surrogates for the hydrophobic interaction between the viral protein and human angiotensin-converting enzyme 2 (ACE2) receptors to determine the potential of IVM as a repurposed drug for SARS-CoV-2 prevention and treatment. The IVM, RBD, and combination binding kinetics were analyzed using biolayer interferometry (BLI) and validated with multiple in silico techniques including protein–ligand docking, molecular dynamics simulation, molecular mechanics-generalized Born surface area (MM-GBSA), and principal component analysis (PCA). Our results suggest that with increasing IVM concentrations the association rate with the hydrophobic biosensor increases with a simultaneous decrease in dissociation. Significant kinetic changes to RBD, when combined with IVM, were found only at a concentration a thousand times the approved dosage with minimal changes found over a 35-min time period. Our study suggests that IVM is not an effective preventative or treatment method at the currently approved dosage.

Evaluation of novel Anti-SARS-CoV-2 compounds by targeting nucleoprotein and envelope protein through homology modeling, docking simulations, ADMET, and molecular dynamic simulations with the MM/GBSA calculation

The current prominent virus that induces severe acute respiratory syndrome is SARS-CoV-2. The incidence of COVID-19 cases is increasing, necessitating the immediate development of effective treatments. Our objective was to employ an in-silico approach to evaluate the effectiveness of conventional compounds against COVID-19's nucleoprotein and envelope protein. A docking simulation was performed on 9 compounds as SARS-coronavirus inhibitors using AMDock software. Anti-covid-19 activities were further evaluated for the compounds. Based on docking results, the binding affinity of "N-(4-carbamoylphenyl)-8-cyclopropyl-7-(naphthalen-1-ylmethyl)-5-oxo-2,3-dihydro-5H-thiazolo[3,2-a]pyridine-3-carboxamide,” also called compound 36 in this research, was found to be −8.8 ​kcal/mol for the modeled envelope protein and −7.3 ​kcal/mol for the template envelope protein, while −10.1 ​kcal/mol for the modeled nucleocapsid proteins (NP) and −8.7 ​kcal/mol for the template nucleocapsid proteins (NP) of SARS-coronavirus, respectively. The ligand and control drug (ritonavir) with high docking scores were subjected to pharmacological screening, molecular dynamic simulations, and Molecular Mechanics-generalized Born Surface Area (MM/GBSA) calculations. Furthermore, the jobs of pharmacokinetics were assessed, and the outcomes acquired show that the proposed compound 36 includes great oral bioavailability and a capacity to diffuse through various organic boundaries. The protein-ligand complexes were subjected to dynamic simulation analyses with a re-enactment time of 10 ns, likewise, their free binding energy was inspected operating the MM/GBSA approach. The docking (MD simulation) results acquired emphasize the pivotal residues answerable for the protein-ligand interaction, giving an understanding of the method of association. The MD simulation analysis verifies the structural stability of the selected complexes during the MD trajectory, with minor changes detected. The MM/GBSA data show that compound 36 has the lowest free energy of −12.498 ​kcal/mol for EP and −57.5185 ​kcal/mol for NP proteins of SARS-coronavirus, confirming the molecular docking result. As a result, the identified chemical can be used to develop a new family of antiviral medications against SARS-coronavirus-2.