Molecular docking and toxicity studies of nerve agents against acetylcholinesterase (AChE)

Acetylcholinesterase (AChE) is a cholinergic enzyme that plays an essential role in the autonomic nervous system. This enzyme is often the target of many nerve agents. When this enzyme is inhibited, its function to hydrolyze acetylcholine is stopped, accumulating the acetylcholine in the tissue and causing prolonged stimulation. Some of the significant nerve agents include sarin (GB), soman (GD), tabun (GA), and venomous agent (VX). In order to determine which compound is the most stable and has the best affinity, the nerve agent venomous agent (VX), sarin (GB), soman (GD), and tabun (GA) are docked to the acetylcholinesterase (AChE) enzyme. After that, toxicity tests will be performed on 17 targets for the compound that was chosen. Autodock Vina 1.2.0 is the software used for the docking procedure. should use the Pymol program version 2.5.4 for analysis and the Ligplot software version 2.2 for visualization of the docking findings. The 'Tox Prediction' algorithm from Insilico was used to determine the toxicity of various substances. Based on the results of molecular docking, the free binding energy of Donepezil, sarin (GB), soman (GD), tabun (GA), and venomous agent (VX) in kcal/mol are -12,3, -4.8, -6.0, -5,1, and -6.3 respectively. Finally, four ligands bind strongly to the receptor Donepezil at RMSD 0.327 Å, and the venomous agent (VX) compound binds the most strongly compared to the other test ligands. Furthermore, in the toxicity test of Compound VX, which exhibits neurotoxic activity, no toxic activity was observed on specific organs and targets.

Molecular modeling and simulations of some antiviral drugs, benzylisoquinoline alkaloid, and coumarin molecules to investigate the effects on Mpro main viral protease inhibition

Background

SARS-CoV-2 is a deadly viral disease and uncounted deaths occurs since its first appearance in the year 2019. The antiviral drugs, benzylisoquinoline alkaloids, and coumarin molecules were searched using different online engines for drug repurposing with SARS-CoV-2 and to investigate the effects on main viral protease (Mpro) upon their bindings.

Methods

A database composed of antiviral drugs, benzylisoquinoline alkaloids, and Coumarin molecules was screened through a molecular docking strategy to uncover the interactions of collected molecules with SARS-CoV-2 Mpro. Further, molecular dynamics simulations (MDS) were implemented for 100 ns to calculate the stability of the best complexed molecular scaffold with Mpro. The conformations of the simulated complexes were investigated by using principal component analysis (PCA) and Gibbs energy landscape (FEL) and DSSP together. Next, free binding energy (ΔG_bind) was calculated using the mmpbsa method.

Results

Molecular docking simulations demonstrate 17 molecules exhibited better binding affinity out of 99 molecules present in the database with the viral protease Mpro, followed ADMET properties and were documented. The Coumarin-EM04 molecular scaffold exhibited interactions with catalytical dyad HIS41, CYS145, and neighboring amino acids SER165 and GLN189 in the catalytical site. The crucial factor RMSD was calculated to determine the orientations of Coumarin-EM04. The Coumarin-EM04 complexed with Mpro was found stable in the binding site during MDS. Furthermore, the free energy binding ΔG_bind of Coumarin-EM04 was found to be −187.471 ± 2.230 kJ/mol, and for Remdesivir ΔG_bind was −171.926 ± 2.237 kJ/mol with SARS-CoV-2 Mpro^.

Conclusion

In this study, we identify potent molecules that exhibit interactions with catalytical dyad HIS41 and CYS145 amino acids and unravel Coumarin-EM04 exhibited ΔG_bind higher than Remdesivir against Mpro and thus may serve better antiviral agent against SARS-CoV-2.

Design of a protein-targeted DNA aptamer using atomistic simulation

The concentrations of specific macromolecular species can be quantified using diagnostic tools that rely on molecular recognition by nucleic acid aptamers. One such approach involves the formation of osmium tetroxide 2,2'-bipyridine protein adducts, followed by electrochemical detection of analytes that bind specifically to electrode-tethered aptamers. In conjunction with a 27-mer DNA aptamer that binds specifically to exosite II on human alpha thrombin, this technique permits, in theory, a highly sensitive diagnostic tool for the quantification of serum thrombin levels. However, thrombin's aptamer binding site is lined by two tryptophan residues and the conjugation of bulky osmium groups to these residues weakens aptamer binding by an estimated 4 to 12 kcal/mol, undermining detection sensitivity. Therefore, we have rationally modified this DNA aptamer to strengthen its thrombin binding in the presence of conjugated osmium. Specifically, aptamers carrying long hydrophobic thymine derivatives in place of guanine 21 have binding affinities for osmium-conjugated thrombin that are enhanced by 10 to 15 kcal/mol, suggesting that these modified aptamers may be effective in a highly sensitive electrochemical sensor for the quantification of low concentrations of thrombin. Our approach of using molecular simulation to subtly re-engineer a DNA aptamer may be generally applicable for the optimization of other macromolecular binding interfaces.Communicated by Ramaswamy H. Sarma.

LY3041658/ interleukin-8 complex structure as targets for IL-8 small molecule inhibitors discovery using a combination of in silico methods

Since interleukin-8 (IL-8/CXCL8) and its receptor, CXCR1 and CXCR2, were known in the early 1990s, biological pathways related to these proteins were proven to have high clinical value in cancer and inflammatory/autoimmune conditions treatment. Recently, IL-8 has been identified as biomarker for severe COVID-19 patients and COVID-19 prognosis. Boyles et al. (mAbs 12 (2020), pp. 1831880) have published a high-resolution X-ray crystal structure of the LY3041658 Fab in a complex human CXCL8. They described the ability to bind to IL-8 and the blocking of IL-8/its receptors interaction by the LY3041658 monoclonal antibody. Therefore, the study has been designed to identify potential small molecules inhibiting interleukin-8 by targeting LY3041658/IL-8 complex structure using an in silico approach. A structure‑based pharmacophore and molecular docking models of the protein active site cavity were generated to identify possible candidates, followed by virtual screening with the ZINC database. ADME analysis of hit compounds was also conducted. Molecular dynamics simulations were then performed to survey the behaviour and stability of the ligand-protein complexes. Furthermore, the MM/PBSA technique has been utilized to evaluate the free binding energy. The final data confirmed that one newly obtained compound, ZINC21882765, may serve as the best potential inhibitor for IL-8.

Deciphering the molecular mechanisms of selective non-covalency demonstrated differentially by 9-Allylnaphtho[1,8-ef]isoindole-7,8,10(9H)-trione (C11) against fibroblast growth factor receptors 1-4

The specific inhibition of aberrant Fibroblast Growth Factor Receptors (FGFRs) has been identified as a feasible strategy to therapeutically ameliorate their respective carcinogenic involvements. High homology among these proteins has however limited efforts towards the discovery of selective small-molecule compounds due to undesirable effects elicited by pan-FGFR inhibitors. A recent study showed the selective activity of a new compound C11 which was >52 times more potent against FGFR1 than FGFR2 and FGFR3, and 4 times than FGFR4. This C11 selective non-covalency was investigated in this study using computational methods since it has remained unresolved. Structural findings revealed that C11 enhanced structural perturbations in FGFR1 with less prominent effects in other FGFRs. High deviations also characterized the C11-bound active pocket of FGFR1 with notable fluctuations across the constituent P-loop, αC helix, hinge region, catalytic, and activation loops. These induced motions were essential for optimal C11 motion an d positioning of its phenalenone ring and prop-2-en-l-yl moiety at the FGFR1 active pocket to interact stably and strongly with A564^FGFR1, L484^FGFR1, Y563^FGFR1, and E562^FGFR1 which as well had high energy contributions. C11 exhibited highly unstable binding in F GFRs2-3 with a more steady interaction with FGFR4. Free binding energy (ΔG_bind) analyses further estimated the highest interaction energy for C11-FGFR1 with favorable desolvation energy that indicated a deep hydrophobic pocket binding for C11 in FGFR1 compared to other FGFRs. We believe rational insights from this study will contribute to the structure-based design of highly specific FGFR1 inhibitors. Communicated by Ramaswamy H. Sarma.

Dipeptidyl peptidase IV inhibition of phytocompounds from Artocarpus champeden (Lour.) Stokes: In silico molecular docking study and ADME-Tox prediction approach

The present study examines the potential activity prediction based on free binding energy (ΔG) and interaction confirmation of phytocompounds from Artocarpus champeden (Lour.) Stokes with macromolecule protein receptor of dipeptidyl peptidase IV (DPP-IV) using in silico molecular docking studies and physicochemical and pharmacokinetic properties (ADME-Tox) prediction approaches. The active subsites of the DPP-IV receptor macromolecule protein Protein Data Bank (ID: 1 × 70) were docked using Autodock v4.2.6 (100 docking runs). A grid box of 52 × 28 × 26 Å points spaced by 0.37 Å was centered on the active site of x = 40.926 Å; y = 50.522 Å; z = 35.031 Å. For ADME-Tox prediction, Swiss ADME online-based application programs were used. The results show that 12 pythocompounds from A. champeden have the potential as DPP-IV inhibitors based on ΔG value and interaction conformation. There are five pythocompounds with lower ΔG values and inhibition constants than the native ligand and seven pythocompounds with ΔG values and inhibition constants close to the native ligand. The 12 compounds form an interaction conformation at the active subsites of the DPP-IV receptor. At the same time, the results of the ADME-Tox prediction analysis showed that the 12 compounds had different physicochemical and pharmacokinetic properties.

Insight of druggable cannabinoids against estrogen receptor β in breast cancer

Breast cancer (BC) is the second most prevalent cancer worldwide. Estrogen receptor beta (ERβ) is an essential protein of breast cells to suppress estrogen-induced uncontrolled proliferation. Thus, small molecules that can modulate and enhance ERβ expression would be an effective agent to suppress BC development. Studies showed that cannabinoid (CB), specifically delta-9-tetrahydrocannabinol (Del9THC), can increase the expression of ERβ and inhibits BC cell proliferation. In this study, less psychoactive and structurally similar analogs of Del9THC were chosen as drug candidates and ERβ was targeted as a therapeutic receptor. Delta-8-tetrahydrocannabinol (Del8THC) and delta-4-isotetrahydrocannabinol (Del4isoTHC) were the drug candidates selected on the basis of literature reports, absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties, medicinal chemistry profile, and physicochemical features. Molecular docking simulations were carried out to determine ligand receptor interactions and binding affinity based on free binding energy. To get a better drug, the structural modification was done on Del8THC and generated a new CB analog called Cannabinoid A. Finally, molecular interaction analysis revealed that two CBs and one of their analog interact with the active site residues of ERβ. Therefore, this study revealed a new way to discover novel drug(s) for BC patients.Communicated by Ramaswamy H. Sarma.