COVID-19 liver and gastroenterology findings: An in silico analysis of SARS-CoV-2 interactions with liver molecules

Abstract

BACKGROUND Coronavirus disease 19 (COVID-19) has not only been shown to affect the respiratory system, but has also demonstrated variable clinical presentations including gastrointestinal tract disorders. In addition, abnormalities in liver enzymes have been reported indicating hepatic injury. It is known that severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) might infect cells via the viral receptor angiotensin-converting enzyme 2 (ACE2) which is expressed in several organs including the liver. The viral Spike glycoprotein binds to ACE2 and must be cleaved by Furin and Type 2 Serine Protease to enter the cells. After that, the Akt/mTOR signaling pathway is activated and several COVID-19 changes are triggered. AIM To analyze liver and gastrointestinal symptoms and cell signaling pathways triggered by SARS-CoV-2 infection due to virus-liver interactions in silico . METHODS In this in silico study, the three-dimensional structures of the Akt, mTORC1 and Furin (receptors) were selected from the Protein Data Bank (PDB) and the structures of inhibitors (ligands) MK-2206, CC-223 and Naphthofluorescein were selected from PubChem and ZINC databases. Ligand files were downloaded as 2D structures and converted to optimized 3D structures using ViewerLite 4.2 software. Marvin Sketch ® software was used to calculate prediction of the protonated form of inhibitors in a physiological environment (pH 7.4). AutoDock Tools (ADT) software was used to calculate and delimit the Grid box used in the molecular docking of each structure selected in the PDB. In addition, protonated ligands were prepared for molecular docking using ADT software. Molecular docking was performed using ADT software tools connected to Vina software. Analysis of the amino acid residues involved in ligand interactions, as well as ligand twists, the atoms involved in interactions, bond type and strength of interactions were performed using PyMol ® and Discovery Studio ® (BIOVIA) software. RESULTS Molecular docking analysis showed that the mTORC1/CC-223 complex had affinity energy between the receptor and ligand of -7.7 kcal/moL with interactions ranging from 2.7 to 4.99 Å. There were four significant chemical bonds which involved two of five polypeptide chains that formed the FKBP12–Rapamycin-Binding (FRB) domain. The strongest was a hydrogen bond, the only polar interaction, and Van der Waals interactions shown to be present in 12 residues of mTORC1’s FRB domain. With regard to the Akt/MK-2206 complex there were three Van der Waals interactions and 12 chemical bonds in which seven residues of Akt were involved with all five rings of the MK-2206 structure. In this way, both ASP 388 and GLN 391 bind to the same MK-2206 ring, the smaller one. However, LYS 386 had four chemical bonds with the...

Screening of Small-Molecule Libraries Using SARS-CoV-2-Derived Sequences Identifies Novel Furin Inhibitors

“…In line with our results, permethrin, a Furin inhibitor discovered by in silico screening and validated against synthetic substrates, acts by targeting a likewise Pocket 3 exosite, involving Trp531 and Ala532 [44], much like the anticipated action of P3. The same cavity seems to be able to accommodate other possible Furin inhibitors, such as vitamin B12 and folic acid [45], naphthofluorescein [46], and epicatechin gallate [47]. Therefore, it would be compelling to explore chemical modifications of P3 or analogous compounds.…”