Next Article in Journal
A Model Assessment of the Occurrence and Reactivity of the Nitrating/Nitrosating Agent Nitrogen Dioxide (NO2) in Sunlit Natural Waters
Next Article in Special Issue
5-Amino-3-methyl-Isoxazole-4-carboxylic Acid as a Novel Unnatural Amino Acid in the Solid Phase Synthesis of α/β-Mixed Peptides
Previous Article in Journal
Flower-like Highly Open-Structured Binder-Free Zn-Co-Oxide Nanosheet for High-Performance Supercapacitor Electrodes
Previous Article in Special Issue
Tirzepatide, a New Era of Dual-Targeted Treatment for Diabetes and Obesity: A Mini-Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Actions of Novel Angiotensin Receptor Blocking Drugs, Bisartans, Relevant for COVID-19 Therapy: Biased Agonism at Angiotensin Receptors and the Beneficial Effects of Neprilysin in the Renin Angiotensin System

by
Graham J. Moore
1,2,*,
Harry Ridgway
3,4,
Konstantinos Kelaidonis
5,
Christos T. Chasapis
6,7,
Irene Ligielli
8,
Thomas Mavromoustakos
8,
Joanna Bojarska
9 and
John M. Matsoukas
1,5,10,*
1
Department of Physiology and Pharmacology, Cumming School of Medicine, University of Calgary, Calgary, AB T2N 4N1, Canada
2
Pepmetics Inc., 772 Murphy Place, Victoria, BC V8Y 3H4, Canada
3
Institute for Sustainable Industries and Liveable Cities, Victoria University, Melbourne, VIC 8001, Australia
4
AquaMem Consultants, Rodeo, New Mexico, NM 88056, USA
5
NewDrug PC, Patras Science Park, 26504 Patras, Greece
6
NMR Facility, Instrumental Analysis Laboratory, School of Natural Sciences, University of Patras, 26504 Patras, Greece
7
Institute of Chemical Engineering Sciences, Foundation for Research and Technology, Hellas (FORTH/ICE-HT), 26504 Patras, Greece
8
Department of Chemistry, National and Kapodistrian University of Athens, 15784 Athens, Greece
9
Institute of General and Ecological Chemistry, Faculty of Chemistry, Lodz University of Technology, Zeromskiego 116, 90-924 Lodz, Poland
10
Institute for Health and Sport, Victoria University, Melbourne, VIC 3030, Australia
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(15), 4854; https://doi.org/10.3390/molecules27154854
Submission received: 15 June 2022 / Revised: 20 July 2022 / Accepted: 26 July 2022 / Published: 29 July 2022
(This article belongs to the Special Issue Advances in Research of Short Peptides II)

Abstract

:
Angiotensin receptor blockers (ARBs) used in the treatment of hypertension and potentially in SARS-CoV-2 infection exhibit inverse agonist effects at angiotensin AR1 receptors, suggesting the receptor may have evolved to accommodate naturally occurring angiotensin ‘antipeptides’. Screening of the human genome has identified a peptide (EGVYVHPV) encoded by mRNA, complementary to that encoding ANG II itself, which is an inverse agonist. Thus, opposite strands of DNA encode peptides with opposite effects at AR1 receptors. Agonism and inverse agonism at AR1 receptors can be explained by a receptor ‘switching’ between an activated state invoking receptor dimerization/G protein coupling and an inverse agonist state mediated by an alternative/second messenger that is slow to reverse. Both receptor states appear to be driven by the formation of the ANG II charge-relay system involving TyrOH-His/imidazole-Carboxylate (analogous to serine proteases). In this system, tyrosinate species formed are essential for activating AT1 and AT2 receptors. ANGII is also known to bind to the zinc-coordinated metalloprotease angiotensin converting enzyme 2 (ACE2) used by the COVID-19 virus to enter cells. Here we report in silico results demonstrating the binding of a new class of anionic biphenyl-tetrazole sartans (‘Bisartans’) to the active site zinc atom of the endopeptidase Neprilysin (NEP) involved in regulating hypertension, by modulating humoral levels of beneficial vasoactive peptides in the RAS such as vasodilator angiotensin (1–7). In vivo and modeling evidence further suggest Bisartans can inhibit ANG II-induced pulmonary edema and may be useful in combatting SARS-CoV-2 infection by inhibiting ACE2-mediated viral entry to cells.

1. Introduction

Angiotensin II is a peptide hormone that mediates the main actions of the renin–angiotensin system (RAS), playing a crucial role in the control of blood pressure and fluid balance. Next, two forms of angiotensin were recognized, namely inactive decapeptide (D-R-V-Y-I-H-P-F-H-L) angiotensin I (ANG I) and powerful octapeptide (DRVYIHPF) angiotensin II (ANG II), the most potent pressor substance known. The vasoconstrictor ANGII peptide was formed through the rapid cleavage of inactive decapeptide ANG I by angiotensin-converting the enzyme (ACE1), dipeptidyl carboxypeptidase, or chymases [1,2,3], whereas ACE2 produces the vasodilator heptapeptide ANG (1–7) from ANGII. In addition, ANG II metabolites, such as ANG III, ANG IV, and ANG (1–7) should not be overlooked [4]. ANG II can bind to the ANG II receptors type 1 and 2 (AT1 and AT2) belonging to the G-protein-coupled receptors. However, most physiological influence is mediated by AT1. More specifically, ANG II acts at AT1 receptors in smooth muscle and other tissues, causing contraction and blood pressure elevation. Lowering ANG II levels with RAS inhibitors, such as ACE1 inhibitors and ANG II receptor blockers (ARB), by direct inhibition with AT1 receptor blockers, provides effective treatment of hypertension. ARBs demonstrate insurmountable blocking effects at AT1 receptors, which are like the effects seen with certain synthetic analogues of ANG II in which the C-terminal Phe is replaced by an aliphatic amino acid such as Ile.
Notably, ANG II is a multifunctional bioactive short peptide. Its regulatory actions are related to not only vasoconstriction (constriction of the blood vessels), increased renal tubular absorption of sodium, and stimulation of adrenal aldosterone production but also to activation of the sympathetic nervous system or the pathogenesis of cardiovascular diseases [2]. ANG II, by AT1, causes cell growth and cellular phenotypic changes, regulates the gene expression of bio-substances, and activates multiple intracellular signaling cascades in cardiac myocytes and fibroblasts and endothelial and smooth muscle cells [3]. Interestingly, an ANG II vaccine, under clinical studies, can be a novel effective therapy against heart failure and hypertension [5,6]. The broad role of angiotensin II as a potent modulator of the immune system has been demonstrated in inflammation, immunity, rheumatoid arthritis, and multiple sclerosis [7,8]. Pioneer research by Steinman and collaborators has shown that blocking an angiotensin-converting enzyme with angiotensin receptor blockers (ARBs) such as Sartans induces potent regulatory T cells and modulates TH1- and TH17-mediated potency [9]. This study is a guide for investigating ARBs and peptide mimetic drugs as possible regulators in the immunotherapy of autoimmune diseases [10,11]. ANG II is involved in diverse pathological situations in terms of tissue remodeling, cell proliferation, migration, tissue invasion, or angiogenesis [12,13]. Thus, ANG II has also relevant in the metastasis of cancers, since a relationship between hypertension and cancer exists [14].
ANG II, the key pressor short peptide, has been a very well-studied peptide by us, and diverse ANG II analogues have been synthesized towards the design and synthesis of non-peptide mimetics as angiotensin II receptor blockers [15,16,17]. In particular, extensive structure activity, nuclear magnetic resonance, and fluorescence studies revealed the importance of the N-terminal domain and of the proline residue to stabilize the active conformation of angiotensin II [18,19,20,21]. In this study, we investigated the mechanism of action of peptide and nonpeptide ARBs at the molecular and pharmacological level, with the aim of understanding the details of their AT1 receptor interactions. Since ARBs (sartans) are, in essence, ‘ANG II look-alikes’, it is anticipated that they will bind to Neprilysin and the ACE2 enzyme, which is the receptor used by the COVID-19 virion for cell entry. In this regard, we have developed a new class of ARB called bisartans (Bis), which contain two biphenyltetrazole groups, enabling bivalent interactions with critical receptor-based positively charged groups, such as (1) R167 of the angiotensin AT1 receptor and (2) Zinc2+ at the active site of ACE2 and Neprilysin. Accordingly, bisartans could have a dual role against COVID-19 by (1) preventing ANG II-induced toxic effects (pulmonary edema/inflammation/cytokine storm) and (2) blocking ACE2-spike protein interactions, particularly involving R and Zn active sites, thereby inhibiting cell entry of the virion.
The beneficial role of Neprilysin (NEP) in Renin Angiotensin System is depicted in Scheme 1. In particular, Angiotensinogen (1–255) is cleaved by renin in the circulation to generate inactive decapeptide angiotensin1 (Ang1). Angiotensin 1 is cleaved by ACE to yield the octapeptide angiotensin II and by Neprilysin enzyme to yield heptapeptide angiotensin (1–7). In contrast to ACE producing vasoconstrictor angiotensin II, Neprilysin(NEP), also a zinc depended protease, is an activator of the alternative RAS in the murine and kidney, producing beneficial vasodilator angiotensin (1–7) from Ang (1–10) and from Ang (1–9), which could lead to therapeutic strategies. Additionally, heptapeptides, angiotensin (1–7), and angiotensin (Ala1–7) (Alamandine) are produced by ACE2, which degrades inflammatory angiotensin II in the renin angiotensin system, counterbalancing the AT1 receptor/angiotensin II axis. The heptapeptides lack phenylalanine, which is the structural requirement in angiotensin II for vasoconstriction and hypertension, and the roles of ACE2 and Neprilysin in the RAS are beneficial [22,23]. This study investigates in silico the docking of sartans and bisartans with neprilysin, in an effort to understand the molecular mechanisms of these interactions.

2. Methods

2.1. Isolated Tissues

Rat isolated smooth muscle tissue bioassays were conducted, as described previously [24,25,26]. Dose-response data were plotted in the format dose/response vs. dose to identify linearity and then as Hill plots, to obtain levels of cooperativity from the slope (Hill coefficient). These protocols for analyzing and interpreting cooperativity in dose-response curves have been described in detail previously [25,26,27].

2.2. In Silico Studies

Docking of ligands to the neprilysin (NEP) (PDB ID: 5JMY) [28] was also performed using AutoDock Vina [29] with default parameters. Point charges and dihedral barriers were initially assigned according to the YAMBER14 force field [30]. However, YAMBER point charges were damped to mimic less polar Gasteiger charges used to optimize the AutoDock scoring function. Docking was performed using non-periodic (walled) boundaries that effectively confined ligands to an approximately 17 × 17 × 32 Å cuboid volume. The setup was implemented using the YASARA molecular-modeling program [31]. The best hits and ligand conformational poses, expressed as kcal/mol free energy of binding, resulting from a minimum of 100 runs per ligand, were reported.
MD simulations were run with YASARA [31]. The setup included an optimization of the hydrogen bonding network to increase the solute stability, and a pKa prediction to fine-tune the protonation states of protein residues at the chosen pH of 7.4. NaCl ions were added with a physiological concentration of 0.9%, with an excess of either Na or Cl to neutralize the cell. After steepest descent and simulated annealing minimizations to remove clashes, the simulation was run for up to 120 ns using the AMBER 14 force field [32] for the solute, GAFF2 [33] and AM1BCC [34] for ligands, and TIP3P for water. The cutoff was 8 Å for van der Waals forces (the default used by AMBER) [35], and no cutoff was applied to electrostatic forces (using the Mesh Ewald algorithm) [36]. MD simulation of the docked Neprilysin-Bisartan complexes was carried out using an NPT ensemble with periodic boundaries at 311 K, including explicit transferable intermolecular potential with 3 points (TIP3P) solvation at liquid density (0.997 g/cc) using AMBER14 parameters. Sodium and chloride ions were added at physiological concentration (0.9 wt%), with an excess of either ion used to maintain system neutrality. The in-vacuo molecular dynamics relaxation of Bis-A occurred in the presence of ZnCl or CH3-guanidine. The MD simulations were run for 3.5 ns in vacuo at 298 K using AMBER-14 parameters. Binding energies and affinity constants should be viewed qualitatively and in relative terms—they highlight differences between ligands and do not take into account entropy differences.

3. Results and Discussion

3.1. Two-State Receptor and Biased Agonism

The octapeptide ANG II acts at AR1 receptors in smooth muscle and other tissues causing contraction and blood pressure elevation. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection involves a viral cell-entry mechanism via angiotensin-converting enzyme 2 (ACE2), implicating angiotensin as a target for potential therapies. Lowering ANG II levels with angiotensin-converting enzyme (ACE) inhibitors or with AR1 receptor blockers (ARBs) provides effective treatment of hypertension. ARBs demonstrate blocking effects at AR1 receptors, which are similar to the effects seen with certain synthetic analogs of ANG II in which the C-terminal Phe is replaced by an aliphatic amino acid such as Ile.
Initially, [Sar1Ile8]ANG II (sarilesin) was designated as a Type 1 desensitizing antagonist [different from competitive surmountable Type II antagonists such as [Sar1Tyr(Me)4]ANG II (sarmesin)], because its effects were insurmountable and reminiscent of the desensitization/tachyphylaxis brought on by supramaximal doses of agonist [24]. However, detailed investigation of dose-response data, transformed into Hanes–Woolf plots and Hill plots [25], has shown that sarilesin induces negative cooperativity (Hill coefficient nH < 1) signifying negative efficacy, which is synonymous with inverse agonism (Table 1). Activities of ANG II analogues are classified as follows: Sar-Arg-Val-Tyr-Ile-His-Pro-Phe (superagonist), Sar-Arg-Val-Tyr(Me)-Ile-His-Pro-Phe (Sarmesin, surmountable antagonist), and Sar-Arg-Val-Tyr-Ile-His-Pro-Ile (Sarilesin, insurmountable inverse agonist) (Table 2).
Moreover, whereas submaximal doses of ANG II result in positive cooperativity (nH > 1 < 2) invoking receptor dimerization with a consequent increase in affinity for ANG II [26], supramaximal doses of ANG II can invoke negative cooperativity and tachyphylaxis (Figure 1). Hill coefficients also show differences in tissue responsiveness to ANG II, with receptor avidity in the decreasing order: aorta, portal vein, and uterus (Table 2).
The prolonged effect of [Sar1Ile8]ANG II does not appear to be due to a slow dissociation rate of this ligand from the receptor, because binding-displacement curves for 125I-ANG II and 125I-[Sar1Ile8]ANG II are superimposable [26]. It seems likely, therefore, that the receptor is locked down as a result of slow dissociation of the second messenger, which is coupled to the ‘inverted’ state of the receptor (Figure 1).
It is worth noting that the receptor mechanisms shown are not exclusive to angiotensin. Receptor dimerization and negative efficacy (inverse agonism) considerations apply to all G proteins coupled with seven transmembrane domain receptors, including those for nonpeptide ligands [27]. Indeed, the threshold effects and tachyphylaxis effects commonly observed in many pharmacological assays are symptomatic of positive and negative cooperativity, respectively.

3.1.1. Angiotensin Peptides and Antipeptides—Regulators of RAS

The replacement of the C-terminal Phe residue of ANG II with an aliphatic amino acid (Ile, Leu, Val, Ala) converts the molecule to an analog with inverse agonist properties. It was of interest to find out if such a molecule is a naturally occurring product encoded by the human genome and designed to counteract ANG II. The sequence of human mRNA, which is complementary to that for ANG II [37], was found to encode the ANG II analogue EGVTVHPV [and likewise the ANG III analogue GVTVHPV]. The C-terminal Val residue confers inverse agonist properties to these ‘antipeptides’, whereas the absence of the Arg side chain has only modest effects on activity. Accordingly, angiotensin ‘antipeptides’ may be produced as a natural countermeasure to angiotensin, by binding to the inverse agonist site (Figure 1).
ANG II has a fundamentally important role in the regulation of cardiovascular function, and, perhaps not surprisingly, its actions are controlled at multiple levels. Axis control points include genetic expression and enzymatic biosynthesis of components of the RAS, giving rise to ANG II, metabolism by ACE2 to the vasodilator peptide ANG II (1–7) acting at AT2 receptors, and signaling effects at AT1 receptors. The present findings add yet another potential level of control for physiological regulation of the RAS, via the production of endogenous antipeptides acting directly as inverse agonists at AT1 receptors. Angiotensin receptor blockers (ARB sartans) used to treat hypertension (see below) apparently bind to a site on the angiotensin receptor which evolved to accommodate naturally occurring peptide inverse agonists.

3.1.2. Charge-Relay Systems

Evidence has been accumulated from structure-reactivity, NMR, and fluorescence lifetime studies which suggest that activation of AT1 receptors is driven by the formation of a charge-relay system (CRS) in ANG II, encompassing triad TyrOH-His-Phe carboxylate, which generates tyrosinate anion species for activating the receptor (Figure 2A). ANG II CRS mechanism is analogous to that in Serine proteases [38]. Notably, methylation of the Tyr hydroxyl group of the superagonist [Sar1]ANG II, or of the inverse agonist [Sar1Ile8]ANG II, leads in both cases to reversible surmountable antagonists, implying the Tyr hydroxyl has a critical role at both the agonist and the inverse agonist states of the receptor. It is the amino acid at position 8 that ultimately determines which receptor state is preferred for binding [39].
The N-terminal portion of ANG II is thought to have a largely supportive role in maintaining the CRS [40]. Accordingly, the role of Arg2 in ANG II is to help stabilize the CRS and mainly act as a chaperone of the C-terminal business end of the molecule [40,41]. Its role is essentially replaced or supplanted by R167 of the receptor when ANG II binds to its receptor (see below). The role of Asp at position 1 of ANG II is considered non-obligatory, and possibly detrimental, since its replacement in [Sar1]ANG II results in a superagonist. Crystallography studies [42,43] have implicated a CRS-like interaction for the binding of the ARB inverse agonist olmesartan to AR1 receptors, wherein an intermolecular CRS is formed with the Tyr35 residue of the receptor (Figure 2B). For ANG II, CRS formation involves an intramolecular interaction (Figure 2A). In another study by Asada et al. [44] the X-ray structure of ANG II bound to AT2 R revealed an interaction of the peptide with Met128, suggested to be a key for receptor activation. The crystal pose of this structure has shown furthermore a hydrogen bond between Tyr OH and PheCOO within the peptide, indicating the formation of tyrosinate species in line with our ANG II model, where tyrosine hydroxylate is interacting with the receptor [45].

3.1.3. Angiotensin Receptor Blockers (ARB Sartans)

The first nonpeptide ARB to be discovered was the surmountable antagonist losartan, which is metabolized to the inverse agonist losartan carboxylate in the bloodstream. This transformation involves the conversion of an imidazole-based neutral hydroxymethyl group of losartan to the carboxylate found in losartan carboxylate (Figure 2A). Similar considerations apply to the neutral carboxamide and carboxylate versions of olmesartan, which are likewise a surmountable antagonist and an inverse agonist, respectively [42]. As outlined above, the peptide equivalents are the TyrOMe and TyrO species in surmountable sarmesin and insurmountable sarilesin, respectively. In other words, it is the carboxylate group in ARBs that defines the inverse agonist activity (Figure 2A), and it is the tyrosinate species in [Sar1Ile8]ANG II that has the same functional role [40,41]. However, the tyrosinate anion in ANG II has two functions: stimulating contraction (positive cooperativity) at submaximal doses and inducing inverse agonism (negative cooperativity) at supramaximal doses. The charge-relay system in ANG II creates two negative charges, Tyr O and COO, the spacing of which is exactly mimicked by the two acid groups in ARBs (Figure 3). These anion pairs form salt bridges with Arg167 and Lys199 on the receptor [42,43,46] (Figure 3). Salt bridges can provide energy for receptor binding affinity, to both the agonist and inverse agonist binding states of the receptor (Figure 1). For the agonist state, the dimerization process, which results in amplification of the contractile response, requires the presence of the Phe ring of ANG II [which is notably absent in ARBs]. The triggering process possibly involves a quadrupolar interaction of the Phe ring of ANG II with the Trp84 ring on the receptor.

3.1.4. SARS-CoV-2, ARBs, and Bisartans

SARS-CoV-2 (COVID-19) patients often present compromised lung function due to severe pulmonary edema. We have noted when conducting rat pressor assays that bolus injections of EC50 doses of ANG II, carried out at 10–20-min intervals, result after several hours in chronic pulmonary edema, as evidenced by the appearance of fluid backing up the tracheal tube. Pulmonary edema seen in COVID-19 infection should respond to ARBs counteracting the effect of toxic ANG II. In addition, ARBs, which are angiotensin look-alikes, should bind to ACE2 and may disrupt the binding of the COVID-19 spike protein, which is vital for the cell entry of the virus [15,47,48]. A recent study reported the discovery and facile synthesis of a new class of sartan-like arterial antihypertensive drugs (angiotensin receptors blockers [ARBS], referred to as ‘bisartans’ [16]. Bisartans are novel bis-alkylated imidazole sartan derivatives bearing dual symmetric inionic biphenyl tetrazole moieties, which in in silico studies exhibited higher binding affinities for the ACE2/spike protein complex (PDB ID: 6LZG), compared to all other known sartans. Bisartans are classified as BisA, B, C, and D, where BisA and B bear the butyl group of position 2 in the imidazole ring, and BisC and D are at position 4, as in losartan [16]. In another study, we further investigated the structural hydropathy traits and polarity properties of mutants that induce the SARS-CoV-2 dominant mutation, which increase transmissibility and infectivity [16]. Arginine blockers, such as ARBs and bisartans, which bear anionic tetrazoles and carboxylates, were found in in silico studies to be ideal candidate drugs for viral infection by weakening S-protein RBD binding to ACE2 and discouraging hydrolysis of cleavage sites [16]. An example of the strong binding of BisA with methylguanidine is where amino groups of guanidine moiety crab the two tetrazoles of BisA. In a similar manner, the two tetrazoles of BisD are crabbed by Zn+2 (Figure 4 and Figure 5).

3.1.5. Bisartans Bound to AT1R

We have developed a new generation of ARBs, called bisartans [15,49]], the structures of which include two tetrazole groups (Figure 6). Bisartans A and B bear the butyl group at position 4 of the imidazole ring, while bisartans C and D are at position 2 as in losartan. Bisartans are potent insurmountable antagonists at AR1 receptors and contain powerful chelators of the zinc atom at the active site of zinc proteases such as ACE2. Coordination of ACE2 Zn+2 with ANG II and analogues involving four residues D, Y, H, and the C-terminal carboxylate has been reported [50]. Knowing that tyrosinate in sarilesin and carboxylate in EXP3174 were both effective at salt bridging to R167, we wondered if tetrazolate, which is a functional mimetic of carboxylate, might also suffice. By using imidazole as the template for mounting two biphenyl tetrazole (BPT) groups, we conjectured that the resulting positive charge on the imidazole should be at about the right distance from the two tetrazoles to subsume the space normally occupied by the Arg2 guanidino group of ANG II, when the peptide is bound to the receptor (Arg2 is predicted to interact with D263 and D281 of the receptor [51]). It turned out that bisartans are indeed potent insurmountable antagonists at AT1 receptors [52] and are the first ARBs discovered that do not use a carboxylate for binding to Arg167 but instead use a second tetrazolate to form the critical salt bridge with Arg167, which appears to provide for insurmountable activity. Bisartans likely exhibit significant bioactivity at AT1 receptors because of tight binding to R167 due to the particular geometry of the double tetrazole unit surrounding the R167 guanidino group. Essentially, bisartans can enfold the R167 guanidino group in the tight embrace of two tetrazole moieties [52]. Interactions between receptor and ARBs were determined from the crystal structure reported by Zhang et al. [42,43]. The critical interactions of olmesartan are with R167, W84, Y35, and K199 residues of AT1R (Figure 6).

3.1.6. Bisartans Bound to NEP

Here we report new in silico results suggesting bisartans also undergo specific binding to the zinc-dependent endopeptidase NEP, which plays a pivotal role in regulating humoral levels of beneficial vasoactive polypeptides, including atrial and brain natriuretic peptides, bradykinin, adrenomedullin, and endothelin-15 [53,54]. Inhibition of NEP causes increased levels of the natriuretic peptides and, thus, represents an important therapeutic target for controlling blood pressure and heart disease in humans. In the present study, a series of 17 conventional angiotensin receptor blockers (ARBs), including four novel bisartans (BisA–BisD) and two potent NEP antagonists, LBQ657 and Candoxatrilat, were docked to the soluble extracellular domain of NEP (PDB ID: 5JMY). Results indicated that AutoDock Vina accurately reproduced the correct X-ray crystallographic posture of LBQ657 in the main Zn2+ pocket of NEP with a root-mean-square deviation (RMSD) ~0.01Å. Dissociation constants (Kd) for BisA and LBQ657 (32,095.31 and 49,026.0703, respectively) were in agreement with their binding energies, although the per-atom binding efficiency for LBQ657 0.3561 kcal/(mol*Atom) was somewhat greater than that of BisA 0.3561 kcal/(mol*Atom).
Of the 17 ligands evaluated, the bisartan BisA exhibited the strongest binding, with two other bisartans (C, D) ranking in the top five scoring positions (Figure 7B). BisA binding was stabilized by hydrophobic and π–π resonance interactions of the drug with contacting NEP residues, as well as cation–π interactions between a terminal tetrazole group and the zinc ion. Molecular dynamics simulations carried out for approximately 40 ns at 311 K in physiological saline demonstrated the solvated NEP-BisA complex (PDB 5JMY-BisA) was stable (average drug RMSD = 1.20 A) over this period (Figure 8). A nearly identical result (average BisA RMSD = 1.64 Å) was obtained for a longer 90 ns MD simulation run under the same conditions, but using a complex comprised of BisA docked to the zinc pocket of human neprilysin (PDB 6GID) (data not shown).
The initial docked conformation of BisA (at t = 0 ns) underwent a conformational change during the initial 300 ps of an MD simulation, causing one tetrazole group of BisA to approach the zinc ion resulting in a stabilizing ionic (cation–pi) bonding interaction. This tetrazole group also formed additional ionic and hydrogen-bonding interactions with Arg717 (Figure 9).

3.1.7. ARBs: Potential Antiviral Drugs for the Treatment of COVID-19

This study was based on the knowledge that cardiovascular disease is entirely related to COVID-19, in terms of mechanisms that trigger the disease [55]. The storm of cytokines released in COVID-19 patients with pneumonia is related to the overexpression of toxic ANG II in the renin–angiotensin system (RAS) [56,57], and clinical studies have shown that morbidity and mortality rate was lower in hypertensive patients infected by SARS-CoV-2 who are taking angiotensin-converting enzyme inhibitors and ANG II-receptor blockers (ARBs), compared to patients not taking these drugs [58,59,60,61,62,63,64,65,66,67]. ANG II, the major RAS component, and SARS-CoV-2 spike protein cleavage, by furin and 3CL, which initiate infection, operate through similar charge-relay-system mechanisms. Tyrosinate in ANG II, serinate in furin, and cysteinate in 3C-like protease (3CLpro) are anions created through the charge-relay-system mechanism and trigger activity via their nucleophile anions. Similar charge-relay-system mechanisms spark and mediate the action of ANG II [45] and the serine-like proteases furin and 3CL [68,69,70,71]. Bisartans were found to bind strongly to the SARS-CoV-2 RBD/ACE2 complex with higher affinity compared to other common sartans. In a similar mode as in the interaction of ARBs with AT1R, bisartans block the critical aminoacid arginine in the furin cleavage site 681–686, which catalyzes the cleavage of the spike protein and triggers COVID-19. Mutation P681R, as in the Delta variant, increases the basicity of the rich arginine cleavage site and the severity of the infection [47]. Bisartans have also been found to bind in silico to 3CL protease such as Pfizer’s Paxlovid combination drugs (PF-07321332 and Ritonavir) [68,69,70]. However, in the COVID-19 assay, BisD was found to be less effective [16,72]. The catalytic center of furin is the triad Asp–His–Ser, and for the 3CL protease, it is the dyad Cys145–His41 [65]. Bisartans act at the three cell entries (ACE2, Furin, 3CLpro) and are ideal potential antiviral drugs for the treatment of COVID-19.

3.1.8. Effects of ACE, ACE2, and NEP in the Renin–Angiotensin System

The above findings with NEP may be extrapolated to the zinc metalloprotease ACE2, which is the receptor used by the COVID-19 virion to enter host cells. This is justified, as both NEP and ACE2 are zinc proteases. Schematic representation of the renin–angiotensin system shows the same beneficial effect of the two proteases, as both produce vasodilatory Ang(1–7) [73] (Scheme 1). Accordingly, it is anticipated that bisartan could bind to the zinc atom in ACE2 and may allosterically inhibit binding and cell entry of the virion (which binds to a different site on the enzyme). In COVID-19, ARBs are known to alleviate some of the symptoms of infection. Bisartans could have a dual role by (1) blocking ANG II induced toxic effects (pulmonary edema/inflammation/cytokine storm) and (2) preventing ACE2–spike protein interactions, particularly involving R and Zn active sites, thereby inhibiting cell entry of the virion.

4. Conclusions

Two states of the AT1 receptor, agonist and inverse agonist states, interact, respectively, with surmountable and insurmountable ARBs and angiotensin analogues. The inverse agonist state can bind a naturally occurring angiotensin ‘antipeptide’, which is encoded by mRNA complementary to that for human ANG II. Agonist activation of the receptor may involve CRS exchange of the Tyr4 of ANG II with Y35 of the receptor [52], whereas inverse agonism is characterized by salt-bridge formation between the ligand and R167 of the receptor [52]. Accordingly, ARB sartans interact preferentially with R167 cation of the AT1R, resulting in the switching of the receptor into its inverted agonist state, and bisartans have the opportunity for bivalent (2 tetrazoles) bonding and, thereby, a stronger interaction leading to increased antihypertensive activity.
Tetrazole can also bind/chelate divalent cations found at the active site of metalloproteases, and Recife model 1 (RM1) energy calculations for the bivalent interaction (both tetrazoles) of bisartan with Zn2+ (−165 kcal/mol) or with the guanidine group of Arg (−106 kcal/mol) are relative to a C-C covalent bond (−330 kcal/mol). In silico studies reveal that bisartans can chelate the active site zinc atom of the enzyme NEP, although only one of the two tetrazoles was found to bind to zinc (Figure 9B), either due to steric hindrance preventing accessibility of the second tetrazole, to compensating multiple secondary interactions, or to limitations of the modeling software/computing power. Visual inspection of the 3D structure in Figure 9 suggests that there is likely enough space to accommodate two tetrazole groups, which is perhaps not surprising since the zinc pocket is designed to bind elongated peptide sequences of natural substrates. Tetrazole is a weak acid [and a functional mimetic of carboxyl, pKa~4], which will bind to protein basic groups in the order Zn2+ > Arg > Lys > others, so there is considerable energy to be gained from the binding of both tetrazoles to zinc. On the other hand, the positive charge on the central imidazole ring of bisartan also has to be taken into account and is shown here (Figure 9B) bonded to F544 of NEP in a parallel-plate interaction (as would be expected for a positively charged imidazole ring bonding to the Phe ring quadrupole). In the initial docking (Figure 9A), R102 and R110 outcompete the zinc ion, but, after relaxing the system (via MD simulation), this tetrazole reorients to interact with zinc (Figure 9B). Why only one tetrazole interacts with zinc is probably best explained by considering that once one tetrazole is captured by zinc, the other tetrazole has to climb out of an energy well to wrap back around the zinc ion. This may be energetically difficult since the butyl chain and biphenyl of bisartan are stabilized by hydrophobic interactions with I535, Y545, G548, T206, F544, and H587 (Figure 9A). Collectively, these interactions, together with the central imidazole, discourage or preclude the dihedral rotations needed to move the second tetrazole into a position where it might begin to interact with zinc.

Author Contributions

Conceptualization and coordination, G.J.M. and J.M.M.; molecular dynamic calculations, H.R., T.M., I.L. and C.T.C.; experimental work and methodology, K.K.; writing—original draft preparation, G.J.M. and J.M.M.; text related to angiotensin II, manuscript editing, J.B.; writing—review and editing, G.J.M. and J.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

G.J.M. would like to thank PepMetics for its support. J.M.M. would like to thank the General Secretariat for Research and Technology, Patras Science Park, the Region of Western Greece (Research and Technology), and the pharmaceutical companies Eli Lilly Greece and Uni-Pharma for supporting their research in multiple sclerosis, hypertension, and COVID-19. We thank Nikolaos Maniotis for typing and editing this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Sample of the compound Bisartan A (BisA) is available from the authors.

Abbreviations

ACE2Angiotensin Converting Enzyme 2
AlaAlanine
AMBERAssisted model building with energy refinement
ANG IIAngiotensin II
AR1Androgen receptor 1
ARBs Angiotensin II receptor blockers
ArgArginine
AT1(2)ANG II receptors type 1 (2)
Azil Azilsartan
BisBisartan
BPTBiphenyl tetrazole
CandeCandesartan
3CLpro3C-like protease
COVID-19Coronavirus diseases
CRSCharge-relay system
DIZEDeminazene aceturate
DNADeoxyribonucleic acid
EproEprosartan
EXP3174Losartan carboxylic acid
HisHistidine
IDIdentifier
IleIsoleucine
IrbeIrbesartan
LeuLeucine
Lo Losartan
LysLysine
MDMolecular dynamics
mRNAMessenger ribonucleic acid
NEPNeprilysin
Olme Olmesartan
PDBProtein Data Bank
PhePhenylalanine
ProProline
RASRenin–angiotensin system
RM1Recife model 1
RMSDRoot-mean-square-deviation
SARS-CoV-2Severe acute respiratory syndrome coronavirus 2
TelmiTelmisartan
TIP3PTransferable intermolecular potential with 3 points
TrpTryptophan
TyrTyrosine
UHFUnrestricted Hartree–Fock
ValValsartan
YAMBERYet Another Model Building and Energy Refinement
YASARAYet Another Scientific Artificial Reality Application

Nomenclature for Bisartans

BisA (4-Butyl-N,N0-bis{[[20-(2H-tetrazol-5-yl)]biphenyl-4-yl] methylimidazolium bromide), BisB (4-Butyl-2-hydroxymethyl-N,N0-bis{[20-(2H-tetrazol-5-yl)-biphenyl-4-yl]methyimidazolium bromide), BisC (2-Butyl-4-chloro-5-hydroxymethyl-N,N0-bis{[20-(2H-tetrazol-5-yl)biphenyl-4-yl]methyl}imidazolium bromide) (Dialkylated losartan), BisD (2-Butyl-N,N0-bis{[20-(2H-tetrazol-5-yl)biphenyl-4-yl]methyl} imidazolium bromide).

References

  1. Bumpus, F.M.; Smeby, R.R.; Page, I.H. Angiotensin, the Renal Pressor Hormone. Circ. Res. 1961, 9, 762–767. [Google Scholar] [CrossRef] [Green Version]
  2. Reid, I.A.; Morris, B.J.; Ganong, W.F. The Renin-Angiotensin System. Annu. Rev. Physiol. 1978, 40, 377–410. [Google Scholar] [CrossRef] [PubMed]
  3. Izumi, Y.; Iwao, H. Chapter 186—Angiotensin II Peptides. In Handbook of Biologically Active Peptides, 2nd ed.; Academic Press: Cambridge, MA, USA, 2013; pp. 1369–1376. [Google Scholar] [CrossRef]
  4. Ardaillou, R.; Chansel, D. Synthesis and effects of active fragments of angiotensin II. Kidney Int. 1997, 52, 1458–1468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Watanabe, R.; Suzuki, J.-i.; Wakayama, K.; Maejima, Y.; Shimamura, M.; Koriyama, H.; Nakagami, H.; Kumagai, H.; Ikeda, Y.; Akazawa, H.; et al. A peptide vaccine targeting angiotensin II attenuates the cardiac dysfunction induced by myocardial infarction. Sci. Rep. 2017, 7, 43920. [Google Scholar] [CrossRef]
  6. Mogi, M. Clinical study on angiotensin II vaccination—The first big step. Hypertens. Res. 2021, 45, 162–163. [Google Scholar] [CrossRef] [PubMed]
  7. Chang, Y.; Wei, W. Angiotensin II in inflammation, immunity and rheumatoid arthritis. Clin. Exp. Immunol. 2015, 179, 137–145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Stegbauer, J.; Lee, D.-H.; Seubert, S.; Ellrichmann, G.; Manzel, A.; Kvakan, H.; Muller, D.N.; Gaupp, S.; Rump, L.C.; Gold, R.; et al. Role of the renin-angiotensin system in autoimmune inflammation of the central nervous system. Proc. Natl. Acad. Sci. USA 2009, 106, 14942–14947. [Google Scholar] [CrossRef] [Green Version]
  9. Platten, M.; Youssef, S.; Hur, E.M.; Ho, P.P.; Han, M.H.; Lanz, T.V.; Phillips, L.K.; Goldstein, M.J.; Bhat, R.; Raine, C.S.; et al. Blocking angiotensin-converting enzyme induces potent regulatory T cells and modulates TH1- and TH17-mediated autoimmunity. Proc. Natl. Acad. Sci. USA 2009, 106, 14948–14953. [Google Scholar] [CrossRef] [Green Version]
  10. Katsara, M.; Matsoukas, J.; Deraos, G.; Apostolopoulos, V. Towards immunotherapeutic drugs and vaccines against multiple sclerosis. Acta Biochim. Biophys. Sin. 2008, 40, 636–642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Apostolopoulos, V.; Rostami, A.; Matsoukas, J. The Long Road of Immunotherapeutics against Multiple Sclerosis. Brain Sci. 2020, 10, 288. [Google Scholar] [CrossRef]
  12. Hsueh, W.A.; Do, Y.S.; Anderson, P.W.; Law, R.E. Angiotensin II in Cell Growth and Matrix Production. Adv. Exp. Med. Biol. 1995, 377, 217–223. [Google Scholar] [CrossRef] [PubMed]
  13. Suzuki, Y.; Ruiz-Ortega, M.; Lorenzo, O.; Ruperez, M.; Esteban, V.; Egido, J. Inflammation and angiotensin II. Int. J. Biochem. Cell Biol. 2003, 35, 881–900. [Google Scholar] [CrossRef]
  14. Ishikane, S.; Takahashi-Yanaga, F. The role of angiotensin II in cancer metastasis: Potential of renin-angiotensin system blockade as a treatment for cancer metastasis. Biochem. Pharmacol. 2018, 151, 96–103. [Google Scholar] [CrossRef] [PubMed]
  15. Matsoukas, J.; Apostolopoulos, V.; Zulli, A.; Moore, G.; Kelaidonis, K.; Moschovou, K.; Mavromoustakos, T. From Angiotensin II to Cyclic Peptides and Angiotensin Receptor Blockers (ARBs): Perspectives of ARBs in COVID-19 Therapy. Molecules 2021, 26, 618. [Google Scholar] [CrossRef]
  16. Ridgway, H.; Moore, G.J.; Mavromoustakos, T.; Tsiodras, S.; Ligielli, I.; Kelaidonis, K.; Chasapis, C.T.; Gadanec, L.K.; Zulli, A.; Apostolopoulos, V.; et al. Discovery of a new generation of angiotensin receptor blocking drugs: Receptor mechanisms and in silico binding to enzymes relevant to SARS-CoV-2. Comput. Struct. Biotechnol. J. 2022, 20, 2091–2111. [Google Scholar] [CrossRef] [PubMed]
  17. Apostolopoulos, V.; Bojarska, J.; Chai, T.-T.; Elnagdy, S.; Kaczmarek, K.; Matsoukas, J.; New, R.; Parang, K.; Lopez, O.P.; Parhiz, H.; et al. A Global Review on Short Peptides: Frontiers and Perspectives. Molecules 2021, 26, 430. [Google Scholar] [CrossRef]
  18. Matsoukas, J.; Cordopatis, P.; Belte, U.; Goghari, M.H.; Ganter, R.C.; Franklin, K.J.; Moore, G.J. Importance of the N-terminal domain of the type II angiotensin antagonist sarmesin for receptor blockade. J. Med. Chem. 1988, 31, 1418–1421. [Google Scholar] [CrossRef]
  19. Matsoukas, J.M.; Agelis, G.; Hondrelis, J.; Yamdagni, R.; Wu, Q.; Ganter, R.; Moore, D.; Moore, G.J.; Smith, J.R. Synthesis and biological activities of angiotensin II, sarilesin, and sarmesin analogs containing Aze or Pip at position 7. J. Med. Chem. 1993, 36, 904–911. [Google Scholar] [CrossRef]
  20. Hondrelis, J.; Lonergan, G.; Voliotis, S.; Matsoukas, J. One pot synthesis and conformation of N-t-butyloxycarbonyl, O-Phenacyl derivatives of proline and other secondary amino acids. Tetrahedron 1990, 46, 565–576. [Google Scholar] [CrossRef]
  21. Polevaya, L.; Mavromoustakos, T.; Zoumboulakis, P.; Grdadolnik, S.G.; Roumelioti, P.; Giatas, N.; Mutule, I.; Keivish, T.; Vlahakos, D.V.; Iliodromitis, E.K.; et al. Synthesis and study of a cyclic angiotensin II antagonist analogue reveals the role of π*–π* interactions in the C-terminal aromatic residue for agonist activity and its structure resemblance with AT1 non-peptide antagonists. Bioorg. Med. Chem. 2001, 9, 1639–1647. [Google Scholar] [CrossRef]
  22. Kaltenecker, C.C.; Domenig, O.; Kopecky, C.; Antlanger, M.; Poglitsch, M.; Berlakovich, G.; Kain, R.; Stegbauer, J.; Rahman, M.; Hellinger, R.; et al. Critical Role of Neprilysin in Kidney Angiotensin Metabolism. Circ. Res. 2020, 127, 593–606. [Google Scholar] [CrossRef] [PubMed]
  23. Domenig, O.; Manzel, A.; Grobe, N.; Kaltenecker, C.; Kovarik, J.; Stegbauer, J.; Gurley, S.B.; Antlanger, M.; Elased, K.; Saemann, M.; et al. The role of neprilysin in angiotensin 1-7 formation in the kidney. J. Hypertens. 2015, 33, e114–e115. [Google Scholar] [CrossRef]
  24. Moore, G.J.; Franklin, K.J.; Nystrom, D.M.; Goghari, M.H. Structure–desensitization relationships of angiotensin analogues in the rat isolated uterus. Can. J. Physiol. Pharmacol. 1985, 63, 966–971. [Google Scholar] [CrossRef]
  25. Moore, G.J.; Scanlon, M.N. Methods for analyzing and interpreting cooperativity in dose-response curves—I. Antagonist effects on angiotensin receptors in smooth muscle. Gen. Pharmacol. Vasc. Syst. 1989, 20, 193–198. [Google Scholar] [CrossRef]
  26. Scanlon, M.N.; Koziarz, P.; Moore, G.J. The relationship between homotropic and heterotropic cooperativity for angiotensin receptors in smooth muscle. Gen. Pharmacol. Vasc. Syst. 1990, 21, 59–65. [Google Scholar] [CrossRef]
  27. Moore, G.J. Methods for analyzing and interpreting cooperativity in dose-response curves—II. Partial agonists acting on muscarinic receptors in smooth muscle. Gen. Pharmacol. Vasc. Syst. 1989, 20, 199–203. [Google Scholar] [CrossRef]
  28. Burley, S.K.; Bhikadiya, C.; Bi, C.; Bittrich, S.; Chen, L.; Crichlow, G.V.; Christie, C.H.; Dalenberg, K.; Di Costanzo, L.; Duarte, J.M.; et al. RCSB Protein Data Bank: Powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences. Nucleic Acids Res. 2021, 49, D437–D451. [Google Scholar] [CrossRef] [PubMed]
  29. Eberhardt, J.; Santos-Martins, D.; Tillack, A.F.; Forli, S. AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python Bindings. J. Chem. Inf. Modeling 2021, 61, 3891–3898. [Google Scholar] [CrossRef] [PubMed]
  30. Krieger, E.; Nielsen, J.E.; Spronk, C.A.E.M.; Vriend, G. Fast empirical pKa prediction by Ewald summation. J. Mol. Graph. Model. 2006, 25, 481–486. [Google Scholar] [CrossRef] [PubMed]
  31. Krieger, E.; Vriend, G. YASARA View—Molecular graphics for all devices—From smartphones to workstations. Bioinformatics 2014, 30, 2981–2982. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and testing of a general amber force field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef] [PubMed]
  34. Jakalian, A.; Jack, D.B.; Bayly, C.I. Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation. J. Comput. Chem. 2002, 23, 1623–1641. [Google Scholar] [CrossRef] [PubMed]
  35. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins Struct. Funct. Bioinform. 2006, 65, 712–725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577–8593. [Google Scholar] [CrossRef] [Green Version]
  37. Moore, G.J.; Ganter, R.C.; Franklin, K.J. Angiotensin ‘antipeptides’: (−)messenger RNA complementary to human angiotensin II (+)messenger RNA encodes an angiotensin receptor antagonist. Biochem. Biophys. Res. Commun. 1989, 160, 1387–1391. [Google Scholar] [CrossRef]
  38. Blow, D.M.; Birktoft, J.J.; Hartley, B.S. Role of a Buried Acid Group in the Mechanism of Action of Chymotrypsin. Nature 1969, 221, 337–340. [Google Scholar] [CrossRef]
  39. Moore, G.J. Designing peptide mimetics. Trends Pharmacol. Sci. 1994, 15, 124–129. [Google Scholar] [CrossRef]
  40. Matsoukas, J.M.; Hondrelis, J.; Keramida, M.; Mavromoustakos, T.; Makriyannis, A.; Yamdagni, R.; Wu, Q.; Moore, G.J. Role of the NH2-terminal domain of angiotensin II (ANG II) and [Sar1]angiotensin II on conformation and activity. NMR evidence for aromatic ring clustering and peptide backbone folding compared with [des-1,2,3]angiotensin II. J. Biol. Chem. 1994, 269, 5303–5312. [Google Scholar]
  41. Matsoukas, J.M.; Bigam, G.; Zhou, N.; Moore, G.J.I. 1H-NMR studies of [Sar1]angiotensin II conformation by nuclear Overhauser effect spectroscopy in the rotating frame (ROESY): Clustering of the aromatic rings in dimethylsulfoxide. Peptides 1990, 11, 359–366. [Google Scholar] [CrossRef]
  42. Zhang, H.; Unal, H.; Desnoyer, R.; Han, G.W.; Patel, N.; Katritch, V.; Karnik, S.S.; Cherezov, V.; Stevens, R.C. Structural Basis for Ligand Recognition and Functional Selectivity at Angiotensin Receptor. J. Biol. Chem. 2015, 290, 29127–29139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zhang, H.; Unal, H.; Gati, C.; Han, G.W.; Liu, W.; Zatsepin, N.A.; James, D.; Wang, D.; Nelson, G.; Weierstall, U.; et al. Structure of the Angiotensin Receptor Revealed by Serial Femtosecond Crystallography. Cell 2015, 161, 833–844. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Asada, H.; Inoue, A.; Ngako Kadji, F.M.; Hirata, K.; Shiimura, Y.; Im, D.; Shimamura, T.; Nomura, N.; Iwanari, H.; Hamakubo, T.; et al. The Crystal Structure of Angiotensin II Type 2 Receptor with Endogenous Peptide Hormone. Structure 2020, 28, 418–425.e4. [Google Scholar] [CrossRef] [PubMed]
  45. Moore, G.J.; Matsoukas, J.M. Angiotensin as a model for hormone—Receptor interactions. Biosci. Rep. 1985, 5, 407–416. [Google Scholar] [CrossRef] [PubMed]
  46. Takezako, T.; Unal, H.; Karnik, S.S.; Node, K. Current topics in angiotensin II type 1 receptor research: Focus on inverse agonism, receptor dimerization and biased agonism. Pharmacol. Res. 2017, 123, 40–50. [Google Scholar] [CrossRef] [PubMed]
  47. Xia, S.; Liu, M.; Wang, C.; Xu, W.; Lan, Q.; Feng, S.; Qi, F.; Bao, L.; Du, L.; Liu, S.; et al. Inhibition of SARS-CoV-2 (previously 2019-nCoV) infection by a highly potent pan-coronavirus fusion inhibitor targeting its spike protein that harbors a high capacity to mediate membrane fusion. Cell Res. 2020, 30, 343–355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Lan, J.; Ge, J.; Yu, J.; Shan, S.; Zhou, H.; Fan, S.; Zhang, Q.; Shi, X.; Wang, Q.; Zhang, L.; et al. Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor. Nature 2020, 581, 215–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Agelis, G.; Resvani, A.; Koukoulitsa, C.; Tůmová, T.; Slaninová, J.; Kalavrizioti, D.; Spyridaki, K.; Afantitis, A.; Melagraki, G.; Siafaka, A.; et al. Rational design, efficient syntheses and biological evaluation of N, N ′-symmetrically bis-substituted butylimidazole analogs as a new class of potent Angiotensin II receptor blockers. Eur. J. Med. Chem. 2013, 62, 352–370. [Google Scholar] [CrossRef] [Green Version]
  50. Fatouros, P.R.; Roy, U.; Sur, S. Modeling Substrate Coordination to Zn-Bound Angiotensin Converting Enzyme 2. bioRxiv 2021, preprint. [Google Scholar] [CrossRef]
  51. Wingler, L.M.; McMahon, C.; Staus, D.P.; Lefkowitz, R.J.; Kruse, A.C. Distinctive Activation Mechanism for Angiotensin Receptor Revealed by a Synthetic Nanobody. Cell 2019, 176, 479–490.e12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Moore, G.J.; Pires, J.M.; Kelaidonis, K.; Gadanec, L.K.; Zulli, A.; Apostolopoulos, V.; Matsoukas, J.M. Receptor Interactions of Angiotensin II and Angiotensin Receptor Blockers—Relevance to COVID-19. Biomolecules 2021, 11, 979. [Google Scholar] [CrossRef] [PubMed]
  53. Jhund, P.S.; McMurray, J.J.V. The neprilysin pathway in heart failure: A review and guide on the use of sacubitril/valsartan. Heart 2016, 102, 1342–1347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Schiering, N.; D’Arcy, A.; Villard, F.; Ramage, P.; Logel, C.; Cumin, F.; Ksander, G.M.; Wiesmann, C.; Karki, R.G.; Mogi, M. Structure of neprilysin in complex with the active metabolite of sacubitril. Sci. Rep. 2016, 6, 27909. [Google Scholar] [CrossRef] [Green Version]
  55. Guo, J.; Huang, Z.; Lin, L.; Lv, J. Coronavirus Disease 2019 (COVID-19) and Cardiovascular Disease: A Viewpoint on the Potential Influence of Angiotensin-Converting Enzyme Inhibitors/Angiotensin Receptor Blockers on Onset and Severity of Severe Acute Respiratory Syndrome Coronavirus 2 Infection. J. Am. Heart Assoc. 2020, 9. [Google Scholar] [CrossRef]
  56. Moreno, M.; Bataller, R. Cytokines and Renin-Angiotensin System Signaling in Hepatic Fibrosis. Clin. Liver Dis. 2008, 12, 825–852. [Google Scholar] [CrossRef]
  57. Ruiz-Ortega, M.; Ruperez, M.; Lorenzo, O.; Esteban, V.; Blanco, J.; Mezzano, S.; Egido, J. Angiotensin II regulates the synthesis of proinflammatory cytokines and chemokines in the kidney. Kidney Int. 2002, 62, S12–S22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Zhang, P.; Zhu, L.; Cai, J.; Lei, F.; Qin, J.-J.; Xie, J.; Liu, Y.-M.; Zhao, Y.-C.; Huang, X.; Lin, L.; et al. Association of Inpatient Use of Angiotensin-Converting Enzyme Inhibitors and Angiotensin II Receptor Blockers With Mortality Among Patients With Hypertension Hospitalized With COVID-19. Circ. Res. 2020, 126, 1671–1681. [Google Scholar] [CrossRef] [PubMed]
  59. Abassi, Z.A.; Skorecki, K.; Heyman, S.N.; Kinaneh, S.; Armaly, Z. COVID-19 infection and mortality: A physiologist’s perspective enlightening clinical features and plausible interventional strategies. Am. J. Physiol.-Lung Cell. Mol. Physiol. 2020, 318, L1020–L1022. [Google Scholar] [CrossRef] [Green Version]
  60. Abassi, Z.A.; Skorecki, K.; Heyman, S.N.; Kinaneh, S.; Armaly, Z. Reply to Letter to the Editor: “Don’t judge too RAShly: The multifaceted role of the renin-angiotensin system and its therapeutic potential in COVID-19”. Am. J. Physiol.-Lung Cell. Mol. Physiol. 2020, 318, L1029–L1030. [Google Scholar] [CrossRef]
  61. Dambha-Miller, H.; Albasri, A.; Hodgson, S.; Wilcox, C.R.; Khan, S.; Islam, N.; Little, P.; Griffin, S.J. Currently prescribed drugs in the UK that could upregulate or downregulate ACE2 in COVID-19 disease: A systematic review. BMJ Open 2020, 10, e040644. [Google Scholar] [CrossRef]
  62. Meng, J.; Xiao, G.; Zhang, J.; He, X.; Ou, M.; Bi, J.; Yang, R.; Di, W.; Wang, Z.; Li, Z.; et al. Renin-angiotensin system inhibitors improve the clinical outcomes of COVID-19 patients with hypertension. Emerg. Microbes Infect. 2020, 9, 757–760. [Google Scholar] [CrossRef] [PubMed]
  63. Ni, W.; Yang, X.; Yang, D.; Bao, J.; Li, R.; Xiao, Y.; Hou, C.; Wang, H.; Liu, J.; Yang, D.; et al. Role of angiotensin-converting enzyme 2 (ACE2) in COVID-19. Crit. Care 2020, 24, 422. [Google Scholar] [CrossRef]
  64. Rico-Mesa, J.S.; White, A.; Anderson, A.S. Outcomes in Patients with COVID-19 Infection Taking ACEI/ARB. Curr. Cardiol. Rep. 2020, 22, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Sriram, K.; Loomba, R.; Insel, P.A. Targeting the renin−angiotensin signaling pathway in COVID-19: Unanswered questions, opportunities, and challenges. Proc. Natl. Acad. Sci. USA 2020, 117, 29274–29282. [Google Scholar] [CrossRef] [PubMed]
  66. Vaduganathan, M.; Vardeny, O.; Michel, T.; McMurray, J.J.V.; Pfeffer, M.A.; Solomon, S.D. Renin–Angiotensin–Aldosterone System Inhibitors in Patients with COVID-19. N. Engl. J. Med. 2020, 382, 1653–1659. [Google Scholar] [CrossRef]
  67. Warner, F.J.; Rajapaksha, H.; Shackel, N.; Herath, C.B. ACE2: From protection of liver disease to propagation of COVID-19. Clin. Sci. 2020, 134, 3137–3158. [Google Scholar] [CrossRef] [PubMed]
  68. Ahmad, B.; Batool, M.; Ain, Q.u.; Kim, M.S.; Choi, S. Exploring the Binding Mechanism of PF-07321332 SARS-CoV-2 Protease Inhibitor through Molecular Dynamics and Binding Free Energy Simulations. Int. J. Mol. Sci. 2021, 22, 9124. [Google Scholar] [CrossRef] [PubMed]
  69. Berteotti, A.; Vacondio, F.; Lodola, A.; Bassi, M.; Silva, C.; Mor, M.; Cavalli, A. Predicting the Reactivity of Nitrile-Carrying Compounds with Cysteine: A Combined Computational and Experimental Study. ACS Med. Chem. Lett. 2014, 5, 501–505. [Google Scholar] [CrossRef] [PubMed]
  70. Ferreira, J.C.; Fadl, S.; Villanueva, A.J.; Rabeh, W.M. Catalytic Dyad Residues His41 and Cys145 Impact the Catalytic Activity and Overall Conformational Fold of the Main SARS-CoV-2 Protease 3-Chymotrypsin-Like Protease. Front. Chem. 2021, 9. [Google Scholar] [CrossRef]
  71. Owen, D.R.; Allerton, C.M.N.; Anderson, A.S.; Aschenbrenner, L.; Avery, M.; Berritt, S.; Boras, B.; Cardin, R.D.; Carlo, A.; Coffman, K.J.; et al. An oral SARS-CoV-2 M pro inhibitor clinical candidate for the treatment of COVID-19. Science 2021, 374, 1586–1593. [Google Scholar] [CrossRef]
  72. Ridgway, H.; Chasapis, C.T.; Kelaidonis, K.; Ligielli, I.; Moore, G.J.; Gadanec, L.K.; Zulli, A.; Apostolopoulos, V.; Mavromoustakos, T.; Matsoukas, J.M. Understanding the Driving Forces That Trigger Mutations in SARS-CoV-2: Mutational Energetics and the Role of Arginine Blockers in COVID-19 Therapy. Viruses 2022, 14, 1029. [Google Scholar] [CrossRef]
  73. Matsoukas, J.M.; Gadanec, L.K.; Zulli, A.; Apostolopoulos, V.; Kelaidonis, K.; Ligielli, I.; Moschovou, K.; Georgiou, N.; Plotas, P.; Chasapis, C.T.; et al. Diminazene Aceturate Reduces Angiotensin II Constriction and Interacts with the Spike Protein of Severe Acute Respiratory Syndrome Coronavirus 2. Biomedicines 2022, 10, 1731. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Receptor switching from an agonist induced state to a desensitized inverse agonist state. Binding of submaximal doses of ANG II to its receptor induces G protein binding and dimerization of the receptor, a concerted mechanism of positive cooperativity (Hill coefficient nH > 1) leading to an increase in agonist affinity and consequent amplification of the contractile response. Different analogues may invoke different levels of cooperativity (kick up the affinity) more than others, e.g., weak or partial agonists. The maximum response may be limited by the available supply of G protein, without which the mode of receptor binding of ANG II [at supramaximal doses] changes to state#2 and induces negative cooperativity (nH < 1) synonymous with inverse agonism (tachyphylaxis). Thus, at high doses, ANG II itself becomes an inverse agonist [as do many partial agonists (2,4)]. Inverse agonists such as ARB drugs, [Sar1Ile8]ANG II, and angiotensin ‘antipeptides’ are unable to excite the receptor and bind to state#2 forming a salt bridge with Arg167 (shown as -X). The resulting insurmountable ‘inverted’ state of the angiotensin receptor engenders smooth muscle relaxation (vasodilation, via an alternative second messenger) for prolonged periods. This receptor-lockdown effect may be due to a slow rate of dissociation of the second messenger (see text).
Figure 1. Receptor switching from an agonist induced state to a desensitized inverse agonist state. Binding of submaximal doses of ANG II to its receptor induces G protein binding and dimerization of the receptor, a concerted mechanism of positive cooperativity (Hill coefficient nH > 1) leading to an increase in agonist affinity and consequent amplification of the contractile response. Different analogues may invoke different levels of cooperativity (kick up the affinity) more than others, e.g., weak or partial agonists. The maximum response may be limited by the available supply of G protein, without which the mode of receptor binding of ANG II [at supramaximal doses] changes to state#2 and induces negative cooperativity (nH < 1) synonymous with inverse agonism (tachyphylaxis). Thus, at high doses, ANG II itself becomes an inverse agonist [as do many partial agonists (2,4)]. Inverse agonists such as ARB drugs, [Sar1Ile8]ANG II, and angiotensin ‘antipeptides’ are unable to excite the receptor and bind to state#2 forming a salt bridge with Arg167 (shown as -X). The resulting insurmountable ‘inverted’ state of the angiotensin receptor engenders smooth muscle relaxation (vasodilation, via an alternative second messenger) for prolonged periods. This receptor-lockdown effect may be due to a slow rate of dissociation of the second messenger (see text).
Molecules 27 04854 g001
Figure 2. Charge relay networks for ANG II and ARBs. (A) Intramolecular charge-relay system for ANG II and (B) intermolecular charge relay-like binding of ARBs to AT1 receptor [42].
Figure 2. Charge relay networks for ANG II and ARBs. (A) Intramolecular charge-relay system for ANG II and (B) intermolecular charge relay-like binding of ARBs to AT1 receptor [42].
Molecules 27 04854 g002
Figure 3. Spatial equivalence of negative charges in losartan carboxylate (EXP3174) and ANG II charge-relay system.
Figure 3. Spatial equivalence of negative charges in losartan carboxylate (EXP3174) and ANG II charge-relay system.
Molecules 27 04854 g003
Figure 4. In vacuo molecular dynamics relaxation (A,B) of Bis-A [−1e] in the presence of ZnCl [+1e] (C) or CH3-guanidine [+1e] (D). The MD simulations were run for 3.5 ns in vacuo at 298 K using AMBER-14 parameters. (A): binding energies as a function of time [ns] for CH3-guanidine (orange) and ZnCl (blue). (B): binding energies of ligands in kcal/mol averaged over the MD trajectories. (C,D): the most relaxed conformations are shown (over the MD simulations) for the BisA-zinc (−16.42 kcal/mol) and BisA-CH3-guanidine (−9.75 kcal/mol) complexes. In both cases, the ligand was trapped between the anionic tetrazole rings. In the case of the BisA-ZnCl simulation, the chloride ion was oriented toward the net positive imidazole group, while the zinc ion was oriented toward the tetrazole groups. Intermolecular interactions using a 6.0-Å cutoff are indicated as colored lines in (C,D) (refer to embedded table). Atomic charges are color-coded as follows: red = more negative; blue = more positive. Note that higher (less negative) energies correspond to stronger binding. Thus, the order of BisA-ligand binding (from strong to weak) was CH3-guanidine to Chloride ion to Zinc ion.
Figure 4. In vacuo molecular dynamics relaxation (A,B) of Bis-A [−1e] in the presence of ZnCl [+1e] (C) or CH3-guanidine [+1e] (D). The MD simulations were run for 3.5 ns in vacuo at 298 K using AMBER-14 parameters. (A): binding energies as a function of time [ns] for CH3-guanidine (orange) and ZnCl (blue). (B): binding energies of ligands in kcal/mol averaged over the MD trajectories. (C,D): the most relaxed conformations are shown (over the MD simulations) for the BisA-zinc (−16.42 kcal/mol) and BisA-CH3-guanidine (−9.75 kcal/mol) complexes. In both cases, the ligand was trapped between the anionic tetrazole rings. In the case of the BisA-ZnCl simulation, the chloride ion was oriented toward the net positive imidazole group, while the zinc ion was oriented toward the tetrazole groups. Intermolecular interactions using a 6.0-Å cutoff are indicated as colored lines in (C,D) (refer to embedded table). Atomic charges are color-coded as follows: red = more negative; blue = more positive. Note that higher (less negative) energies correspond to stronger binding. Thus, the order of BisA-ligand binding (from strong to weak) was CH3-guanidine to Chloride ion to Zinc ion.
Molecules 27 04854 g004
Figure 5. Semiempirical RM1 calculation of intermolecular binding energies (A) for the BisA-CH3 guanidine and (B) BisA-ZnCl (C) complexes. Electrostatic potential distributions for each complex are shown below the graph. Note that lower (more negative) energies in this calculation indicate stronger binding. The RM1 calculations were carried out in in vacuo systems following energy minimization of the relaxed conformers shown in Figure 1, above. The total system charge for each complex was set to 0e for the Unrestricted Hartree–Fock (UHF) calculations. Geometry optimizations were performed by the Pollack–Ribiere conjugate gradient method with a final RMS gradient of 1.0 kcal/Å-mol and spin multiplicity of 1. The electrostatic potential maps (B,C) suggest the zinc ion (going from +2.0e to +0.776e) functioned as a stronger electron withdrawing group compared to CH3-guanidine (going from +1.0e to +0.827e). This effect was reflected in the total electron charge on the respective BisA molecules of each complex (see black-letter labels).
Figure 5. Semiempirical RM1 calculation of intermolecular binding energies (A) for the BisA-CH3 guanidine and (B) BisA-ZnCl (C) complexes. Electrostatic potential distributions for each complex are shown below the graph. Note that lower (more negative) energies in this calculation indicate stronger binding. The RM1 calculations were carried out in in vacuo systems following energy minimization of the relaxed conformers shown in Figure 1, above. The total system charge for each complex was set to 0e for the Unrestricted Hartree–Fock (UHF) calculations. Geometry optimizations were performed by the Pollack–Ribiere conjugate gradient method with a final RMS gradient of 1.0 kcal/Å-mol and spin multiplicity of 1. The electrostatic potential maps (B,C) suggest the zinc ion (going from +2.0e to +0.776e) functioned as a stronger electron withdrawing group compared to CH3-guanidine (going from +1.0e to +0.827e). This effect was reflected in the total electron charge on the respective BisA molecules of each complex (see black-letter labels).
Molecules 27 04854 g005
Figure 6. Interactions between receptor and olmesartan were determined from the crystal structure. The critical interactions of olmesartan are with R167, W84, Y35, and K199 residues of AT1R (reproduced from Zhang et al. [42,43]).
Figure 6. Interactions between receptor and olmesartan were determined from the crystal structure. The critical interactions of olmesartan are with R167, W84, Y35, and K199 residues of AT1R (reproduced from Zhang et al. [42,43]).
Molecules 27 04854 g006
Figure 7. Docking results for selected bisartans and related ARBs to the zinc-dependent endopeptidase NEP (PDB ID: 5JMY). (A) The Zn2+ pocket is enclosed by the bounding box and the zinc ion is circled (magenta sphere). Docking was performed with AutoDock Vina. AutoDock Vina was able to correctly predict the X-ray pose of the bound inhibitor LBQ657 (see Figure 9A below), with an RMSD value < 0.01Å. (B). Bisartan BisA exhibited the strongest binding, with two other bisartans ranking in the top-five scoring positions.
Figure 7. Docking results for selected bisartans and related ARBs to the zinc-dependent endopeptidase NEP (PDB ID: 5JMY). (A) The Zn2+ pocket is enclosed by the bounding box and the zinc ion is circled (magenta sphere). Docking was performed with AutoDock Vina. AutoDock Vina was able to correctly predict the X-ray pose of the bound inhibitor LBQ657 (see Figure 9A below), with an RMSD value < 0.01Å. (B). Bisartan BisA exhibited the strongest binding, with two other bisartans ranking in the top-five scoring positions.
Molecules 27 04854 g007
Figure 8. Selected frame captures from the MD simulation of the docked Neprilysin-BisA complex, indicating the drug was stably bound. MD was carried out using an NPT ensemble with periodic boundaries at 311 K, including explicit transferable intermolecular potential with 3 points (TIP3P) solvation at liquid density (0.997 g/cc) using AMBER14 parameters (see http://www.yasara.org, accessed on 28 July 2022). Sodium and chloride ions were added at physiological concentration (0.9 wt%) with an excess of either ion used to maintain system neutrality. The average RMSD for the BisA drug was 1.20 Å over the time range of from 1 ns to 39.7 ns, indicating it was stably bound in the zinc pocket. Lower graphic shows the drug RMSD and overall system energy as a function of the MD simulation time.
Figure 8. Selected frame captures from the MD simulation of the docked Neprilysin-BisA complex, indicating the drug was stably bound. MD was carried out using an NPT ensemble with periodic boundaries at 311 K, including explicit transferable intermolecular potential with 3 points (TIP3P) solvation at liquid density (0.997 g/cc) using AMBER14 parameters (see http://www.yasara.org, accessed on 28 July 2022). Sodium and chloride ions were added at physiological concentration (0.9 wt%) with an excess of either ion used to maintain system neutrality. The average RMSD for the BisA drug was 1.20 Å over the time range of from 1 ns to 39.7 ns, indicating it was stably bound in the zinc pocket. Lower graphic shows the drug RMSD and overall system energy as a function of the MD simulation time.
Molecules 27 04854 g008
Figure 9. (A) Superimposition of experimental X-ray pose (blue sticks) and best AutoDock Vina pose (yellow sticks) for the NEP inhibitor LBQ657 (RMSD ~0.01 Å) in the Zn2+ pocket. Zinc ion is shown as magenta sphere. Intermolecular interactions are indicated by colored lines: hydrophobic (green lines), π–π-resonance (red lines), and cation–π (blue lines). (B) Best conformational pose of Bisartan A (BisA) docked within the NEP Zn2+ pocket. This pose was used as the starting (t = 0 ns) conformation for the MD simulation described above (see Figure 9), to monitor ligand stability in the zinc pocket. (C) Last frame (39.7 ns) of the MD simulation of the fully solvated BisA-Neprilysin complex (311 K, NPT), showing rearrangement of the drug to undergo stable cation–π zinc coordination with a terminal tetrazole group.
Figure 9. (A) Superimposition of experimental X-ray pose (blue sticks) and best AutoDock Vina pose (yellow sticks) for the NEP inhibitor LBQ657 (RMSD ~0.01 Å) in the Zn2+ pocket. Zinc ion is shown as magenta sphere. Intermolecular interactions are indicated by colored lines: hydrophobic (green lines), π–π-resonance (red lines), and cation–π (blue lines). (B) Best conformational pose of Bisartan A (BisA) docked within the NEP Zn2+ pocket. This pose was used as the starting (t = 0 ns) conformation for the MD simulation described above (see Figure 9), to monitor ligand stability in the zinc pocket. (C) Last frame (39.7 ns) of the MD simulation of the fully solvated BisA-Neprilysin complex (311 K, NPT), showing rearrangement of the drug to undergo stable cation–π zinc coordination with a terminal tetrazole group.
Molecules 27 04854 g009
Scheme 1. Schematic representation of renin–angiotensin system showing the effects of enzymes ACE, ACE2, and NEP. Resulting peptides, AngII and Ang (1–7), interact, respectively, first with AT1R (vasoconstriction effect) and AT2R (vasodilatation effect) and second with MAS receptor (vasodilatation effect).
Scheme 1. Schematic representation of renin–angiotensin system showing the effects of enzymes ACE, ACE2, and NEP. Resulting peptides, AngII and Ang (1–7), interact, respectively, first with AT1R (vasoconstriction effect) and AT2R (vasodilatation effect) and second with MAS receptor (vasodilatation effect).
Molecules 27 04854 sch001
Table 1. Hill coefficients for angiotensin contraction of rat smooth muscle tissues in the absence and presence* of [Sar1Ile8]ANG II (20nM).
Table 1. Hill coefficients for angiotensin contraction of rat smooth muscle tissues in the absence and presence* of [Sar1Ile8]ANG II (20nM).
Uterus1.40.8 *
Portal vein1.70.8 *
Aorta1.90.9 *
Table 2. Activities of ANG II analogues in smooth muscle.
Table 2. Activities of ANG II analogues in smooth muscle.
SUPERAGONISTSar-Arg-Val-Tyr-Ile-His-Pro-Phe
ANTAGONIST (SARMESIN)Sar-Arg-Val-Tyr(Me)-Ile-His-Pro-Phe
(Surmountable)
INVERSE AGONIST (SARILESIN)Sar-Arg-Val-Tyr-Ile-His-Pro-Ile
(Insurmountable)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Moore, G.J.; Ridgway, H.; Kelaidonis, K.; Chasapis, C.T.; Ligielli, I.; Mavromoustakos, T.; Bojarska, J.; Matsoukas, J.M. Actions of Novel Angiotensin Receptor Blocking Drugs, Bisartans, Relevant for COVID-19 Therapy: Biased Agonism at Angiotensin Receptors and the Beneficial Effects of Neprilysin in the Renin Angiotensin System. Molecules 2022, 27, 4854. https://doi.org/10.3390/molecules27154854

AMA Style

Moore GJ, Ridgway H, Kelaidonis K, Chasapis CT, Ligielli I, Mavromoustakos T, Bojarska J, Matsoukas JM. Actions of Novel Angiotensin Receptor Blocking Drugs, Bisartans, Relevant for COVID-19 Therapy: Biased Agonism at Angiotensin Receptors and the Beneficial Effects of Neprilysin in the Renin Angiotensin System. Molecules. 2022; 27(15):4854. https://doi.org/10.3390/molecules27154854

Chicago/Turabian Style

Moore, Graham J., Harry Ridgway, Konstantinos Kelaidonis, Christos T. Chasapis, Irene Ligielli, Thomas Mavromoustakos, Joanna Bojarska, and John M. Matsoukas. 2022. "Actions of Novel Angiotensin Receptor Blocking Drugs, Bisartans, Relevant for COVID-19 Therapy: Biased Agonism at Angiotensin Receptors and the Beneficial Effects of Neprilysin in the Renin Angiotensin System" Molecules 27, no. 15: 4854. https://doi.org/10.3390/molecules27154854

APA Style

Moore, G. J., Ridgway, H., Kelaidonis, K., Chasapis, C. T., Ligielli, I., Mavromoustakos, T., Bojarska, J., & Matsoukas, J. M. (2022). Actions of Novel Angiotensin Receptor Blocking Drugs, Bisartans, Relevant for COVID-19 Therapy: Biased Agonism at Angiotensin Receptors and the Beneficial Effects of Neprilysin in the Renin Angiotensin System. Molecules, 27(15), 4854. https://doi.org/10.3390/molecules27154854

Article Metrics

Back to TopTop