Research Article / Open Access

DOI: 10.31488/jjm.158

Molecular Similarity Relating to a Peptide Sequence Common toAngiotensin II (ANGII) and SARS-Cov2 S-Protein

Wynford Robert Williams*

  1. Faculty of Life Sciences &Education, University of South Wales, Cardiff, U.K

*Corresponding author: W. R. Williams,School of Care Sciences, Faculty of Life Sciences &Education, University of South Wales,Cardiff, CF371DL,U.K ;Tel: 01443 48311

Abstract

The peptidase entity of angiotensin-converting enzyme 2 (ACE2) accommodates the receptor-binding domain of SARS-COV2spike (S) protein initiating cell infection. Compounds binding to ACE2 or the S-protein have the potential to impair cell entry by the virus. The ACE2 substrate angiotensin II (ANGII)shares common amino acid sequences withthe S-protein.This study focuses on one sequence, the N-terminal valine/arginine/aspartate peptideof ANGII. Computer modeling identifiesrelative molecular similarity within the tripeptidesequence andligands for ACE2, ANGII and sigma receptor proteinswith in vitro toxicity forSARS-CoV2.Molecular similarity is also evident withindrugstructures of several neurotransmitter classes and agents from natural sources. Application of these observationsmay proveuseful for the further development of SARS-CoV2antiviral drugs.

Keywords: angiotensin-converting enzyme 2, angiotensin II, antiviral drugs,molecular structure; spike protein, SARS-CoV2,sigma receptor

Introduction

The affinity of SARS-CoV2 (Cv19) envelope proteins for angiotensin-converting enzyme 2 (ACE2) facilitates infection of human cells[1]. ACE2 cell receptorsarepresentwithin several organ systems and expression in lung tissues(elevated in chronic smokers)increases significantly in airways epithelial cells infected with SARS-CoV[2]. There is80% genome sequence identity between SARS-CoV2 and SARS-CoV1,with the latter’sreceptor-binding motif (RBM) for ACE2 comprising of residues 424-494 [3,4]. The amino acid sequence of the Cv19 Spike S1 receptor-binding domain (RBD)isavailable [5]. ACE enzymes and angiotensin peptides are essential components of the complex renin-angiotensin system (RAS) [1]. Individual susceptibility to Cv19 infection may be partially explained by natural variation within the different components of RAS [6].

ACEcarboxypeptidasesare ectoenzymes that cleave amino acid residues from the C-terminus of proteins. Angiotensin 2 (ANGII)an octapeptide (1-8) formed from angiotensin I(1-10) by ACEissubsequently depleted of amino acids by aminopeptidases to ANGIII and smallerpeptides with various effects on blood pressure, aldosterone and sodium retention [7]. The ACE2 glycoprotein, integral to cell membranes, is anchored at the hydrophobic C-terminus with an extracellular N-terminal region containing the catalytic motif of 5 amino acid residues [8,9]. Catalysis of ANG (1-7) from ANGIIsubstrateby ACE2 is completely inhibited in vitro by high concentrations of the C-terminal dipeptide, Pro.Phe[10]. ACE2 is also susceptible to inhibition or activation by other di- and tripeptides and naturally occurring small molecular weight compounds [9,11,12]. The primary physiologic actions of ANGII are mediated via the AT1 receptor (AT1R).The agonist adopts an extended structure within the binding pocket; N-terminal Val3 interacts with the second extracellular loop, whereas the C-terminus interacts with residues of the 7thtransmembrane domain [13].Asp1, Arg2, Tyr4, Phe8, determine the binding affinity and specificity of ANGII for AT1R, whereas Tyr4, His6, Pro7 and Phe8account for agonist potency [14]. Less is known about the requirements for the recognition of ANGII by ACE2.

Several clinically promising antiviral di- and tripeptide drugs have been developed against HIV and Hepatitis C.Knowledge of the full amino acid sequences of ACE2 and SARS-CoV RBD and the application ofsophisticated computational techniques haveyet to provideblocking agents for clinical use [15]. An alternative approach focuses on drug repurposing, the testingof existing compounds for antiviral properties. AT1R blockers such as losartan provide onestrategy to reduce the initial infective stage of Cv19 andresveratrol, a health supplement and dietary constituent,is another compound of interest [16]. Gordon and co-workers have recently identified a significant number of drugsthat target interaction between viral and human proteins [17]. The drugs, including 4-7 residue peptides and cyclic peptides,have beenassessed for their in vitro effects on the infection, growth and cytotoxicityof SARS-CoV2.

Notwithstandingthe evidence supporting alink between the ACE2 receptorandSARS-CoVinfection, histopathological studies do not all support this mechanism[18].The deficit ofdata on the relationship between ACE2 expression and Cv19 mortality,and the requirement for therapeutic agents to inhibit Cv19 infection informs the need for further research. Search of the S-protein sequence of Cv19 reveals two tripeptidesequences, valine/arginine/aspartate (VRD) and arginine/valine/tyrosine (RVY) each occurring twice within the protein [5]. These four amino acid residuesare also present as the sequence aspartate/arginine/valine/tyrosine (DRVY) within the ANGII structure [19]. This study focuseson the VRD amino acid sequence common to ANGII andthe S- protein of SARS-CoV2. Acomputer-aided modeling approach is applied to the investigation ofmolecular similarity within the structures of VRD,ligands of the ANGII system, sigma (σ) antiviral inhibitors and structurally similarcompounds of potential interest.

Methods

The Nemesis software program (Oxford Molecular version 2.1) is used to build molecular structures from contents of the program fragment file and minimise structures by conformational analysis. The VRD peptide and compound structures are minimum energy conformers in an uncharged form. The computational program fits selected paired molecular structures on a three-point basis. Fitting-points comprise of atoms of similar type and partial charge within compound and peptidestructures, identified in the figures with respect to amino acid residue labels. Compound colour-coded atoms in the figures identify ligand fitting-points: carbon-green, nitrogen-blue, oxygen-red, sulphur-yellow. To improve on presentation, bond order within the molecular structures is not shown and some peptide templates are repositioned with respect to each other to provide a better image of the compounds. The Nemesis program computes goodness-of-fit values, in respect of inter-atomic distance at each fitting point and root mean square (RMS) value.

Results

Table 1 lists the investigated compoundswiththeir primary properties and fitting data. The data,encapsulating the fit of each structure to the val/arg/asp (VRD) residues of ANGII (Figure 1 and 2),demonstrate good fitting values with interatomic distances and RMS values respectively ≤0.16Å and ≤0.0200Å. The fits of the ATR ligands BMIC, losartan and L-162313demonstrate the importance of animidazole/imidazopyridine moiety (Figure 1);(templates 3-6). In contrast to BMIC and losartan, L-162313 is an ATRagonist with a fitting-point that can be shared by itssulphonylcarbamate group and the carbonyl of valine [6]Alternative cyclic ringsystems replace the imidazole moiety in structures other than camostat [7].The fits of resveratrol [11,12] are replicated by steroid hormone structures [9,10,13,14]. Progesterone and dexamethasone provide the beststeroid fitsto VRD peptide and one that is replicated by testosterone (0.02Å, 0.08Å, 0.07Å; RMS 0.0069Å – not shown).

Figure 1.Angiotensin II (1) and N-terminal sequence of angiotensin II (2) with compound fits (3-14) to the N-terminal peptide template (grey): 3) BMIC, 4) losartan, 5) L-162313, 6) L-162313, 7) camostat, 8) lovastatin, 9) dexamethasone, 10) progesterone, 11) resveratrol, 12) resveratrol, 13) 17-β estradiol, 14) 17-β estradiol.

Table 1. Fitting data of compounds and the val/arg/asppeptide.

Compound Target/descriptor Fitting points Interatomic distances (Å) RMS (Å)
5-hydroxy indole acetic acid  metabolite C10C5C9 0.08,0.05,0.05 0.0012
17-βE E C10C5C9 0.03,0.16,0.14 0.0123
17-βE E C10C5C9 0.12,0.16,0.05 0.0175
AEAE ACE2 N1N2O10 0.05,0.03,0.06 0.0137
amiodarone ion channels, σ C10C5C9 0.09,0.15,0.08 0.0143
astemizole H C11C10C9 0.09,0.05,0.10 0.0079
BMIC AT N4C9C4 0.04,0.05,0.01 0.0044
BV-T5182 5-HT C10C5C9 0.09,0.10,0.14 0.0041
CJ-1639 D N5C9C10 0.02,0.02,0.01 0.0021
camostat serine protease C11C10C9 0.06,0.06,0.05 0.0021
cloperastine H, σ C10C5N2 0.03,0.02,0.05 0.0019
curcumin protective agent C10C5C9 0.04,0.11,0.09 0.0065
dexamethasone steroid C10C5C9 0.01,0.04,0.04 0.0034
doxorubicin anti-neoplastic C9C5C10 0.08,0.15,0.08 0.0108
haloperidol D, σ O4C4C9 0.01,0.09,0.09 0.0040
KF-26777 A C10C5C9 0.10,0.11,0.03 0.0114
L-162,313 AT N4C9C4 0.01,0.14,0.13 0.0066
L-162,313 AT C4O12C9 0.12,0.06,0.16 0.0200
L-741626 D, 5-HT N4C9C10 0.03,0.09,0.06 0.0068
lisuride α, D, 5-HT N5C9C5 0.03,0.08,0.10 0.0040
losartan AT N4C9C4 0.08,0.01,0.09 0.0034
lovastatin statin, σ C4C5C9 0.05,0.01,0,04 0.0013
LY-344864 5-HT N4C9N1 0.11,0.06,0.05 0.0060
MPFA antioxidant N4C9C10 0.09,0.06,0.15 0.0017
nicotianamine ACE2 N2C10C5 0.08,0.07,0.06 0.0051
pimozide D, 5-HT, H, σ N4C9C5 0.09,0.03,0.10 0.0048
progesterone E, σ C10C5C9 0.02,0.08,0.07 0.0069
prostaglandin I2 prostanoid N5C9C5 0.14, 0.08,0.20 0.0056
RX-821002 α C9C5C10 0.03,0.09,0.11 0.0052
resveratrol protective agent C10C5C9 0.05,0.13,0.12 0.0069
reveratrol protective agent C10C5C9 0.04,0.12,0.13 0.0057
rimcazole σ C11C10N5 0.04,0.03,0.05 0.0017
silibinin protective agent C10C5C9 0.06,0.06,0.01 0.0030
thalidomide anti-neoplastic C9C5C10 0.03,0.12,0.11 0.0029
verapamil ion channels, σ C10C5C9 0.06,0.06,0.04 0.0018

Compounds: AEAE; N-(2-aminoethyl)-1-aziridineethanamine, BMIC; 5-butyl-methyl-imidazole carboxylate 30, MPFA; 3-methyl-5-pentyl-2-furanpentanoic acid.

Receptors: α - adrenergic, σ - sigma, A - adenosine, AT - angiotensin, D - dopamine, E - estrogen, H - histamine, 5-HT - serotonin.

Compound fits to the peptide template in Figures 2 and 3 generally follow the above pattern. Additional fitting characteristics of the structures (Figure 2) include the C11 fitting point of rimcazole and the carbonyl group (C4O4) of haloperidol. In general, peptide fitting-points of the abovecompounds are found on two distinct cyclic rings. The first ring system has one or two fitting-points relating to the arginine guanidinium groupand the second has fitting-points for peptide C4, C5 or C10. These cyclic rings are variously separated by a third cyclic ring, a short alkyl chain, or form part of a fused ring system insteroid-likestructures.The various ring systems fitting to the arginine guanidinium group include aromatic and non aromatic rings (imidazole, piperazine, piperidine,chromone, furan) bicyclic rings (benzeneimidazo, indole, naphthalene, benzofuran) and larger ring systems (imidazopurine, tetracene).Templates 1-7in Figure 2includethe structures of antiviral sigma (σ) receptor ligands identified by Gordon (daunorubicin is replaced here by doxorubicin) [17]. Thepeptide fitting-points of theσ ligands, including lovastatin and progesterone structures, are quite different indescribing triangular pharmacophores of dimensions (Å) 1.25-2.5, 4.8-7.9, 6.3-8.1, in area 4.4-10.1 Å2.Figure 2 also gives the structures of the prostanoidPgI2 and natural lipid MPFAwith fits based on a furan ring and carboxyl group [10,11].

Figure 2.Compound fits (1-11) to the N-terminal peptide template (grey): 1) rimcazole, 2) cloperazine, 3) astemizole, 4) haloperidol, 5) doxorubicin, 6) pimozide, 7) verapamil, 8) amiodarone, 9) thalidomide, 10) prostaglandin I2, 11) MPFA.

Figure 3 expands on the pharmacological diversity within the compounds investigated.These structures, includingclassical receptor agonists and antagonists and natural compounds with more general therapeutic properties, conform to fitting patterns identified in the previous figures.5-hydroxytryptamine (5-HT) antagonists and the metabolite 5-hydroxyindoleacetic acidprovide several different fits to the VRD peptide template. The ACE2 inhibitor structuresAEAEand nicotianamine[11,12]differin not fitting to the arginine guanidinium group andsuperimposing along the axis of the peptide chain. Of the compounds listed in Table 1, 70% fit solely to the arginine residue and 21% to the arginine-aspartate dipeptide.

Figure 3.Compound fits (1-12) to the N-terminal peptide template (grey): 1) lisuride, 2) RX-821002, 3) KF-26777, 4) CJ-1639, 5) L-741626, 6) LY-344864, 7) silibinin, 8) curcumin, 9) BVT-5182, 10) 5-hydroxyindole acetic acid, 11) nicotianamine, 12) AEAE.

Discussion

The results demonstrate relative molecular similarity within the structures of σ receptor ligands, compounds that target the angiotensin pathway,and a tripeptidesequence that contributes to the structures of ANGII and SARS-CoV2 S-protein. The VRD sequenceis also shared by ANG1-9 and ANG1-7 but not by ANGIII. Molecular similarity within the arginine residue, statin and AT1R compounds is perhaps not surprising, as their therapeutic effects in RAS compromised patients are partly due to nitric oxide promoting properties [20]. AT1R blockers and resveratrol are recognized as potential SARS-CoV2 antivirals. Resveratrol displays inhibitory effects against several viruses, including MERS-CoV, and reduces RAS activity in rodent thoracic aorta and a model of nonalchoholic fatty liver disease[21-23]. AEAEand camostat inhibit the function of SARS-CoVS-proteinandprevent Cv19entry into cells [24,25].

Drugs with in vitro toxicity against SARS-CoV2 have been characterized aseither translation inhibitors orσreceptor ligands [17]. The fits of compounds in the latter class(listed in Table 1)do not identify a commonpharmacophore, aspeptide fitting-points and distances between fitting-points are quite different.Although thesecompounds differ instructure, properties and fitting characteristics, they relate to each other through similarity to the peptide VRD. The σ receptor is described as an enigmatic evolutionary isolate, with poorly defined regulation by a remarkably diverse range of ligands [26].There is a marked contrast betweenclassical pharmacologic receptors andthe protein targets ofσ ligands that include enzyme structures.Pharmacophoremapping of over 8500 compounds in a GPCR-focused chemical library has recently created a model of the σ1 receptor, described byfour key chemical features of the ligand-binding pocket alongside eight exclusion volume spheres [27]. Ofthe other drugs investigated here, verapamil, amiodarone, haloperidol, silibinin, genisteinand thalidomide are active in inhibiting the infective or replicative phases of a variety of viruses [28-32].Nicotianamineis an ACE2 inhibitor isolated from soybean whereas AEAE, a blocker of S-protein-mediated cell fusion, isa synthetic compound[9].

Several 5-HT antagonists and the metabolite 5-hydroxyindoleacetic acid relate to the VRD peptide structure. One surprising public health finding of the Cv19 pandemic is the low susceptibility of children to Cv19, raising the possibility of a natural inhibitor to infection within this age group. Interestingly, 5-HT synthesis is highest in childhood, especially between the ages of 2 and 5, declining to adult levels around 11 years of age [33,34]. Data obtained from rodent and cell culture systems demonstrate significant interaction between the serotonin and angiotensin systems. Serotoninergic dorsal raphe nucleus neurons are targets of ANGII; combined 5-HT/ANGII infusion reduces echocardiogram ejection fractions, causing a more aggressive remodeling of valves; blunting of platelet aggregation responses by ANGII relates to loss of the 5-HT transporter (SERT) [35-37].

Minor changes in ACE2 molecular structure are known to interfere with the binding of SARs-CoV. Whereas natural mutations within the RBD may not eliminate infection by SARS-CoV, a single point mutation at aspartate454 abolishes the association of the S1- protein with ACE2 [3,38]. Arginine and valinehave been identified as variable residues within the RBDs of SARS-CoV and SARS-CoV2 that interface with ACE2 [3]. The S-protein domain contains in excess of 3% arginine residues with 16 valine-arginine sequences.

In conclusion, the data presented demonstrate that many drug structures with antiviral effects on SARS-CoV2 share molecular similarity with a peptide sequence within the S-protein, and the same characteristics are evident within compounds as yet untested for antiviral properties. This may be indicative of a common affinity for a viral receptor-protein butdoesnot necessarily relate to theACE2 binding site for the N-terminal residues of ANGII.

References

  1. Gheblawi M, Wang K, Viveiros A, et al.Angiotensin-converting enzyme 2: SARS-COV-2 receptor and regulator of the renin-angiotensin system. Circ Res.2020; 126: 1457-75.

  2. Li G, He X, Zhang L, et al. Assessing ACE2 expression patterns in lung tissues in the pathogenesis of COVID-19. J Autoimmun. 2020.

  3. Yan R, Zhang Y, Li Y, et al. Structural basis for the recognition of SARS-COV-2 by full-length human ACE2. Science 2020; 367: 1444-8.

  4. Li F, Li W, Farzan M, et al. Structure of SARS coronavirus spike-receptor binding domain complexed with receptor. Science 2005; 309: 1864-8.

  5. Chang TJ, Yang DM, Wang M, et al. Genomic analysis and comparative multiple sequence of SARS-CoV2. J Chin Med Assoc. 2020; 83: 537-43.

  6. Hanff TC, Harhay MO, Brown TS, et al. Is there an association between COVID-19 mortality and the renin-angiotensin system-a call for epidemiologic investigations.Clin Infect Dis.2020.

  7. Forrester SJ, Booz GW, Sigmund CD, et al. Angiotensin II signal transduction: an update on mechanisms of physiology and pathophysiology. Physiol Rev.2018; 98: 1627–1738.

  8. Warner FJ, Smith AI, Hooper NM, et al. Angiotensin-converting enzyme-2: a molecular and cellular perspective. Cell MolLife Sci. 2004; 61: 2704-13.

  9. Takahashi S, Yoshiya T, Yoshizawa-Kumagaye K, et al. Nicotianamine is a novel angiotensin-converting enzyme 2 inhibitor in soy-bean. Biomed Res. 2015; 36: 219-24.

  10. Guy JL, Jackson RM, Acharya KR, et al. Angiotensin-converting-enzyme-2 (ACE2): comparative modeling of the active site, specificity requirements, and chloride dependence. Biochemistry 2003; 42:13185-92.

  11. Liao W, Bhullar KS, Chakrabarti S, et al. Egg white-derived tripeptide IRW (Ileu-Arg-Trp) is an activator of angiotensin converting enzyme 2. J Agric Food Chem.2018; 66: 11330-6.

  12. Mores A, Matziari M, Beau F, et al. Development of potent and selective phosphinic peptide inhibitors of angiotensin-converting enzyme 2. J Med Chem. 2008; 51: 2216-26.

  13. Perodin, J, Deraet M, Auger-Messier M, et al. Residues 293 and 294 are ligand contact points of the human angiotensin type 1 receptor. Biochemistry 2002; 41: 14348-56.

  14. Singh KD, Unal H, Desnoyer R, et al. Mechanism of hormone peptide activation of a GPCR: angiotensin II activated state of AT1R initiated by van der Waals attraction. J ChemInf Model. 2019; 59: 373-85.

  15. Han Y, Kral P. Computational design of ACE2-based peptide inhibitors of SARS-Cov-2. ACS Nano 2020; 14: 5143-7.

  16. Magrone T, Magrone M, Jirillo E. Focus on receptors for coronaviruses with special reference to angiotensin-converting enzyme 2 as a potential drug target- a perspective. EndocrMetab Immune Disord Drug Targets 2020; 20: 807-11.

  17. Gordon DE, Jang GM, Bouhaddou M, et al. A SARS-COV-2 protein interaction map reveals targets for drug repurposing. Nature 2020.

  18. Gu J, Korteweg C. Pathology and pathogenesis of severe acute respiratory syndrome. Am J Pathol. 2007; 170: 1136-47.

  19. Alexander SPH, Mathie A, Peters JA.G Protein-coupled Receptors. Br J Pharm. 2011; 164 (S1): S5-S113.

  20. Förstermann U, Sessa WC. Nitric oxide synthetases. Eur Heart J. 2012; 33: 829-37.

  21. Lin SC, Ho CT, Chuo WH, et al. Effective inhibition of MERS-CoV infection by resveratrol. BMC InfectDis.2017; 17: 144

  22. Kim EN, Kim MY, Lim JH, et al. The protective effect of resveratrol on vascular aging by modulation of the renin-angiotensin system. Atherosclerosis 2018; 270: 123-31.

  23. Tiao MM, Lin YJ, Yu HR, et al. Resveratrol ameliorates maternal and post-weaning high-fat diet-induced nonalcoholic fatty liver disease via renin-angiotensin system. Lipids Health Dis. 2018; 17: 178

  24. Huentelman MJ, Zubcevic J, Prada JAH, et al. Structure-based discovery of a novel angiotensin-converting enzyme 2 inhibitor. Hypertension 2004; 44: 903-6.

  25. Hoffman M, Kleine-Weber H, Schroeder S, et al. SARS-COV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 2020; 181: 271-80.

  26. Schmidt HR, Zheng S, Gurpinar E, et al. Crystal structure of the human σ1 receptor. Nature 2016; 532: 527-30.

  27. Alamri MA, Alamri MA. Pharmacophore and docking based sequential virtual screening for the identification of novel sigma 1receptor ligands. Bioinformation2019; 15: 586-95.

  28. Gehring G, Rohrmann K, Atenchong N et al. The clinically approved drugs amiodarone, dronedarone and verapamil inhibit filovirus cell entry. JAntimicrobChemother.2014; 69: 2123-31.

  29. Wunderlich V, Fey F, Sydow G. Antiviral effect of haloperidol on Raucher Murine Leukemia Virus. Arch Geschwulstforsch 1980; 50: 758-62.

  30. Umetsu T, Inoue J, Kogure T, et al. Inhibitory effect of silibinin on hepatitis B virus entry. BiochemBiophys Rep. 2018; 14: 20-5.

  31. LeCher JC, Diep N, Krug PW, et al. Genistein has antiviral activity against Herpes B virus and acts synergistically with antiviral treatments to reduce effective dose. Viruses 2019; 11:

  32. Makonkawkeyoon S, Limson-Pobre RN, Moreira AL, et al. Thalidomide inhibits the replication of Human Immunodeficiency Virus Type I. ProcNatlAcadSci USA. 1993; 90: 5974-78.

  33. Murrin LC, Sanders JD, Bylund DB. Comparison of the maturation of the adrenergic and serotonergic neurotransmitter systems in the brain: implications for differential drug effects on juveniles and adults. BiochemPharmacol. 2007; 73: 1225-36.

  34. Hyland K (2008). Clinical utility of monoamine neurotransmitter metabolite analysis in cerebrospinal fluid. Clin Chem.2008; 54: 633-41.

  35. Pocari CY, Araujo IG, Urzedo-Rodrigues L. et al. Whole-body sodium depletion modifies AT1 mRNA expression and serotonin content in the dorsal raphe nucleus. J Neuroendocrinol. 2019; 31:

  36. Perez J, Diaz N, Tandon I, et al. Elevated serotonin interacts with angiotensin-II to result in altered valve interstitial cell contractability and remodeling. CardiovascEng Technol. 2018; 9: 168-80.

  37. Singh P, Fletcher TW, Li Y, et al. Serotonin uptake rates in platelets from angiotensin II-induced hypertensive mice. Health 2013; 5: 31-9.

  38. Wong SK, Li W, Moore MJ, et al. A 193-amino acid fragment of the SARS coronavirus S protein efficiently binds angiotensin-converting enzyme 2. J Biol Chem.2004; 279: 3197-201.

Received: February 10, 2020;
Accepted: April 07, 2020;
Published: April 09, 2020.

To cite this article : Williams WR.Molecular Similarity Relating to a Peptide Sequence Common to Angiotensin II (ANGII) and SARS-Cov2 S-Protein.Japan Journal of Medicine. 2020;3:2.

©Williams WR.2020.