Published online Jul 27, 2023. doi: 10.4331/wjbc.v14.i4.72
Peer-review started: December 30, 2022
First decision: April 13, 2023
Revised: April 30, 2023
Accepted: July 3, 2023
Article in press: July 3, 2023
Published online: July 27, 2023
Processing time: 206 Days and 8.3 Hours
Coronavirus disease 2019 (COVID-19) has demonstrated several clinical manifestations which include not only respiratory system issues but also liver, kidney, and other organ injuries. One of these abnormalities is coagulopathies, including thrombosis and disseminated intravascular coagulation. Because of this, the administration of low molecular weight heparin is required for patients that need to be hospitalized. In addition, Remdesivir is an antiviral that was used against Middle East Acute Respiratory Syndrome, Ebola, Acute Respiratory Syndrome, and other diseases, showing satisfactory results on recovery. Besides, there is evidence suggesting that this medication can provide a better prognosis for patients with COVID-19.
To investigate in silico the interaction between Remdesivir and clotting factors, pursuing a possibility of using it as medicine.
In this in silico study, the 3D structures of angiotensin-converting enzyme 2 (ACE2), Factor I (fibrinogen), Factor II (prothrombin), Factor III (thromboplastin), Factor V (proaccelerin), Factor VII (proconvertin), Factor VIII (antihemophilic factor A), Factor IX (antihemophilic factor B), Factor X (Stuart-Prower factor), and Factor XI (precursor of thromboplastin (these structures are technically called receptors) were selected from the Protein Data Bank. The structures of the antivirals Remdesivir and Osetalmivir (these structures are called ligands) were selected from the PubChem database, while the structure of Atazanavir was selected from the ZINC database. The software AutoDock Tools (ADT) was used to prepare the receptors for molecular docking. Ions, peptides, water molecules, and other ones were removed from each ligand, and then, hydrogen atoms were added to the structures. The grid box was delimited and calculated using the same software ADT. A physiological environment with pH 7.4 is needed to make the ligands interact with the receptors, and still the software Marvin sketch® (ChemAxon®) was used to forecast the protonation state. To perform molecular docking, ADT and Vina software was connected. Using PyMol® software and Discovery studio® software from BIOVIA, it was possible to analyze the amino acid residues from receptors that were involved in the interactions with the ligands. Ligand tortions, atoms that participated in the interactions, and the type, strength, and duration of the interactions were also analyzed using those software.
Molecular docking analysis showed that Remdesivir and ACE2 had an affinity energy of -8.8 kcal/moL, forming a complex with eight hydrogen bonds involving seven atoms of Remdesivir and five amino acid residues of ACE2. Remdesivir and prothrombin had an interaction with six hydrogen bonds involving atoms of the drug and five amino acid residues of the clotting factor. Similar to that, Remdesivir and thromboplastin presented interactions via seven hydrogen bonds involving five atoms of the drug and four residues of the clotting factor. While Remdesivir and Factor V established a complex with seven hydrogen bonds between six antiviral atoms and six amino acid residues from the factor, and Factor VII connected with the drug by four hydrogen bonds, which involved three atoms of the drug and three residues of amino acids of the factor. The complex between Remdesivir and Factor IX formed an interaction via 11 hydrophilic bonds with seven atoms of the drug and seven residues of the clotting factor, plus one electrostatic bond and three hydrophobic interactions. Factor X and Remdesivir had an affinity energy of -9.6 kcal/moL, and the complex presented 10 hydrogen bonds and 14 different hydrophobic interactions which involved nine atoms of the drug and 16 amino acid residues of the clotting factor. The interaction between Remdesivir and Factor XI formed five hydrogen bonds involving five amino acid residues of the clotting factor and five of the antiviral atoms.
Because of the in silico significant affinity, Remdesivir possibly could act in the severe acute respiratory syndrome coronavirus 2 infection blockade by interacting with ACE2 and concomitantly act in the modulation of the coagulation cascade preventing the hypercoagulable state.
Core Tip: In the initial period of the emergence of coronavirus disease 2019, bioinformatics tools filled the need for rapid knowledge about the severe acute respiratory syndrome coronavirus 2 and research in search of drug treatment, since experimental research takes time that was not available. Bioinformatics continues to be used for the same purpose in this study, remembering that this methodology does not rule out either in vivo or in vitro testing.