Insights into the Binding of Receptor-Binding Domain (RBD) of SARS-CoV-2 Wild Type and B.1.620 Variant with hACE2 Using Molecular Docking and Simulation Approaches

Recently, a new variant, B.1620, with mutations (S477N-E484K) in the spike protein’s receptor-binding domain (RBD) has been reported in Europe. In order to design therapeutic strategies suitable for B.1.620, further studies are required. A detailed investigation of the structural features and variations caused by these substitutions, that is, a molecular level investigation, is essential to uncover the role of these changes. To determine whether and how the binding affinity of ACE2–RBD is affected, we used protein–protein docking and all-atom simulation approaches. Our analysis revealed that B.1.620 binds more strongly than the wild type and alters the hydrogen bonding network. The docking score for the wild type was reported to be −122.6 +/− 0.7 kcal/mol, while for B.1.620, the docking score was −124.9 +/− 3.8 kcal/mol. A comparative binding investigation showed that the wild-type complex has 11 hydrogen bonds and one salt bridge, while the B.1.620 complex has 14 hydrogen bonds and one salt bridge, among which most of the interactions are preserved between the wild type and B.1.620. A dynamic analysis of the two complexes revealed stable dynamics, which corroborated the global stability trend, compactness, and flexibility of the three essential loops, providing a better conformational optimization opportunity and binding. Furthermore, binding free energy revealed that the wild type had a total binding energy of −51.14 kcal/mol, while for B.1.628, the total binding energy was −68.25 kcal/mol. The current findings based on protein complex modeling and bio-simulation methods revealed the atomic features of the B.1.620 variant harboring S477N and E484K mutations in the RBD and the basis for infectivity. In conclusion, the current study presents distinguishing features of B.1.620, which can be used to design structure-based drugs against the B.1.620 variant.

Phytomedicine in Disease Management: In-Silico Analysis of the Binding Affinity of Artesunate and Azadirachtin for Malaria Treatment

In the rural communities of sub-Saharan African (sSA) countries, malaria is being managed using phytocompounds. Artesunate is reported to inhibit Gephyrin E, a central, multi-domain scaffolding protein of inhibitory post-synapses. Neem plant and its metabolites like azadirachtin are being indicated for management of malaria by traditional healers. The present study was aimed to cheminformatically analyse the binding potential of artesunate and azadirachtin with various reactive moieties of Gephyrin E, to reduce malaria scourge. With molecular dynamics (MD), binding free energy estimation and binding affinity of artesunate and azadirachtin to Gephyrin E was done. GRIP docking was done to study the interactions of these test ligands with Gephyrin E (6FGC). MD simulation gave insights to structural changes upon binding of artesunate and azadirachtin in the ligand-binding pocket of Gephyrin E. Root mean square deviation (RMSD) and root mean square fluctuation (RMSF) were calculated. From the estimation, azadirachtin had a total binding energy of −36.97 kcal/mol; artesunate had a binding energy of −35.73 kcal/mol. The GRIP docking results provided a clearer evidence that artesunate has comparatively better binding affinity to Gephyrin E than azadirachtin, and the critical binding sites (in activity order) were cavity 3, 2, 8, and 6 for artesunate while for azadirachtin, it was cavity 6, 3, 8, and 2. The GRIP docking provided detailed interactions at the atomic levels, providing evidence; both compounds have chances to overcome the drug resistance problem, albeit higher for artesunate. Our findings added another piece of evidence that azadirachtin may be effective as an anti-malarial agent. The results herein may provide impetus for more studies into bioactive components of plant origin towards the effective management of malaria disease phenotype.

Computational Evaluation of Abrogation of HBx-Bcl-xL Complex with High-affinity Carbon Nanotubes (Fullerene) to Halt the Hepatitis B Virus Replication

Hepatitis B virus (HBV) is the world’s most prevalent chronic viral infection. More than 350 million individuals are chronic carriers of the virus, with an estimated 2 billion infected persons. For instance, the role of HBx protein in attachment and infection is very obvious and consequently deemed as an important druggable target. Targeting the interface and discovering novel drugs greatly advanced the field of therapeutics development. Therefore, in the current study, HBx to Bcl-xL is abrogated on high-affinity carbon nanotubes using computational structural biology tools. Our analysis revealed that among the total 62 carbon fullerenes, only 13 compounds exhibited inhibitory activity against HBx, which was further confirmed through IFD-based rescoring. Structural dynamics investigation revealed stable binding, compactness, and hydrogen bonds reprogramming. Moreover, the binding free energy calculation results revealed that the top hits1-4 possess the total binding energy of −54.36 kcal/mol (hit1), −50.81 kcal/mol (hit2), −47.09 kcal/mol (hit3), and −45.59 kcal/mol for hit4. In addition, the predicted KD values and bioactivity scores further validated the inhibitory potential of these top hits. The identified compounds need further in vitro and in vivo validation to aid the treatment process of HBV.

Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations

Early Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The estimated total binding energy is ~70 kJ/mol, of which 50–60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem. This leads to a dynamic model whereby the 200 nm long coiled coil attached to the FYVE domain dimer can amplify local hinge-bending motions such that the Rab5-binding domain at the other end of the coiled coil can explore an area of 0.1 μm2 in the search for a second endosome with which to interact.

Structural and Biophysical Investigation of the Key Hotspots on the Surface of Epstein–Barr Nuclear Antigen 1 Essential for DNA Recognition and Pathogenesis

Epstein-Barr Virus (EBV) is considered the most important human pathogen due to its role in infections and cellular malignancies. It has been reported that this Oncolytic virus infects 90% world’s population. EBNA1 is required for DNA binding and survival of the virus and is considered an essential drug target. The biochemical and structural properties of this protein are known, but it is still unclear which residues impart a critical role in the recognition of dsDNA. Intending to disclose only the essential residues in recognition of dsDNA, this study used a computational pipeline to generate an alanine mutant of each interacting residue and determine the impact on the binding. Our analysis revealed that R469A, K514A, Y518A, R521A and R522A are the key hotspots for the recognition of dsDNA by the EBNA1. The dynamics properties, i.e. stability, flexibility, structural compactness, hydrogen bonding frequency, binding affinity, are altered by disrupting the protein-DNA contacts, thereby decreases the binding affinity. In particular, the two arginine substitution, R521A and R522A, significantly affected the total binding energy. Thus, we hypothesize that these residues impart a critical role in the dsDNA recognition and pathogenesis. This study would help to design structure-based drugs against the EBV infections.

A Quantum Chemical Topology Picture of Intermolecular Electrostatic Interactions and Charge Penetration Energy

<div>Basing on the Interacting Quantum Atoms approach, we present herein a conceptual and theoretical framework of short-range electrostatic interactions, whose accurate description is still a challenging problem in molecular modeling. For all the non-covalent complexes in the S66 database, the fragment-based and atomic decomposition of the electrostatic binding energies is performed using both the charge density of the dimers and the unrelaxed densities of the monomers. This energy decomposition together with dispersion corrections gives rise to a pairwise approximation to the total binding energy. It also provides energetic descriptors at varying distance that directly address the atomic and molecular electrostatic interactions as described by point-charge or multipole-based potentials. Additionally, we propose a consistent definition of the charge penetration energy within quantum chemical topology, which is mainly characterized in terms of the intramolecular electrostatic energy. Finally, we discuss some practical implications of our results for the design and validation of electrostatic potentials.</div>

A Quantum Chemical Topology Picture of Intermolecular Electrostatic Interactions and Charge Penetration Energy

<div>Basing on the Interacting Quantum Atoms approach, we present herein a conceptual and theoretical framework of short-range electrostatic interactions, whose accurate description is still a challenging problem in molecular modeling. For all the non-covalent complexes in the S66 database, the fragment-based and atomic decomposition of the electrostatic binding energies is performed using both the charge density of the dimers and the unrelaxed densities of the monomers. This energy decomposition together with dispersion corrections gives rise to a pairwise approximation to the total binding energy. It also provides energetic descriptors at varying distance that directly address the atomic and molecular electrostatic interactions as described by point-charge or multipole-based potentials. Additionally, we propose a consistent definition of the charge penetration energy within quantum chemical topology, which is mainly characterized in terms of the intramolecular electrostatic energy. Finally, we discuss some practical implications of our results for the design and validation of electrostatic potentials.</div>

Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations

Early Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The total binding energy is 70 kJ/mol, of which 50-60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem, allowing binding via a mechanism resembling that of a suction cup connected to a locally rigid stem via a flexible hinge.

Covalent and non-covalent binding free energy calculations for peptidomimetic inhibitors of SARS-CoV-2 main protease

This work employs rigorous absolute binding free energy calculations and QM/MM methods to calculate the total binding energy of two recently crystallized peptidomimetic covalent inhibitors of the SARS-CoV-2 Mpro target.

The Effect of Deuteration on the H2 Receptor Histamine Binding Profile: A Computational Insight into Modified Hydrogen Bonding Interactions

We used a range of computational techniques to reveal an increased histamine affinity for its H2 receptor upon deuteration, which was interpreted through altered hydrogen bonding interactions within the receptor and the aqueous environment preceding the binding. Molecular docking identified the area between third and fifth transmembrane α-helices as the likely binding pocket for several histamine poses, with the most favorable binding energy of −7.4 kcal mol−1 closely matching the experimental value of −5.9 kcal mol−1. The subsequent molecular dynamics simulation and MM-GBSA analysis recognized Asp98 as the most dominant residue, accounting for 40% of the total binding energy, established through a persistent hydrogen bonding with the histamine −NH3+ group, the latter further held in place through the N–H∙∙∙O hydrogen bonding with Tyr250. Unlike earlier literature proposals, the important role of Thr190 is not evident in hydrogen bonds through its −OH group, but rather in the C–H∙∙∙π contacts with the imidazole ring, while its former moiety is constantly engaged in the hydrogen bonding with Asp186. Lastly, quantum-chemical calculations within the receptor cluster model and utilizing the empirical quantization of the ionizable X–H bonds (X = N, O, S), supported the deuteration-induced affinity increase, with the calculated difference in the binding free energy of −0.85 kcal mol−1, being in excellent agreement with an experimental value of −0.75 kcal mol−1, thus confirming the relevance of hydrogen bonding for the H2 receptor activation.