Characterization of ACE2 naturally occurring missense variants: impact on subcellular localization and trafficking
Human Genomics volume 16, Article number: 35 (2022)
Human angiotensin-converting enzyme 2 (ACE2), a type I transmembrane receptor physiologically acting as a carboxypeptidase enzyme within the renin-angiotensin system (RAS), is a critical mediator of infection by several severe acute respiratory syndrome (SARS) corona viruses. For instance, it has been demonstrated that ACE2 is the primary receptor for the SARS-CoV-2 entry to many human cells through binding to the viral spike S protein. Consequently, genetic variability in ACE2 gene has been suggested to contribute to the variable clinical manifestations in COVID-19. Many of those genetic variations result in missense variants within the amino acid sequence of ACE2. The potential effects of those variations on binding to the spike protein have been speculated and, in some cases, demonstrated experimentally. However, their effects on ACE2 protein folding, trafficking and subcellular targeting have not been established.
In this study we aimed to examine the potential effects of 28 missense variants (V801G, D785N, R768W, I753T, L731F, L731I, I727V, N720D, R710H, R708W, S692P, E668K, V658I, N638S, A627V, F592L, G575V, A501T, I468V, M383I, G173S, N159S, N149S, D38E, N33D, K26R, I21T, and S19P) distributed across the ACE2 receptor domains on its subcellular trafficking and targeting through combinatorial approach involving in silico analysis and experimental subcellular localization analysis. Our data show that none of the studied missense variants (including 3 variants predicted to be deleterious R768W, G575V, and G173S) has a significant effect on ACE2 intracellular trafficking and subcellular targeting to the plasma membrane.
Although the selected missense variants display no significant change in ACE2 trafficking and subcellular localization, this does not rule out their effect on viral susceptibility and severity. Further studies are required to investigate the effect of ACE2 variants on its expression, binding, and internalization which might explain the variable clinical manifestations associated with the infection.
The global infection and mortality rates of the Coronavirus Disease 2019 (COVID-19), caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), continue to increase and thus inflecting unprecedent economic and health burden worldwide . The highly contagious SARS-CoV-2 displays a wide spectrum of clinical presentations ranging from asymptomatic infection to severe cardiorespiratory failure that acquire hospitalization, mechanical respiratory, and cardiac support for some cases . Understanding the basis of the extreme interindividual clinical variability may aid clinicians and researchers to assign the most suitable and personalized supportive treatments to patients, improve the efficacy of the available vaccines, and potentially aid in the development of novel effective therapies. Unhealthy habits like smoking or individuals affected by chronic conditions such as obesity, hypertension, cardiovascular diseases, and diabetes are more prone to severe illness and worse prognosis [3, 4]. In addition to the viral genome variability, other risk factors may influence disease clinical severity and presentation including patient’s gender, race, age, and his/her genetic make-up .
Multiple studies have linked genetic variations in angiotensin-converting enzyme 2 (ACE2) to the clinical heterogeneity in infected patients as the receptor is primarily utilized via SARS-CoV-2 for cellular entry to initiate the infection process . Various studies have investigated the influence of naturally occurring single nucleotide polymorphisms (SNPs) on infection variable susceptibility and severity through its effect on binding affinity to the viral S protein targeting ACE2 for COVID-19 therapy [6, 7]. For example, it has been speculated that the East Asian populations would be more susceptible to the severe form of the infection due to certain high allele frequency variants that may lead to higher ACE2 expression . However, multiple ACE2 variants were found to exert protective effects against COVID-19 through impairing ACE2 expression and/or function [9, 10]. ACE2 is a metalloproteinase type 1 transmembrane protein made of 805 amino acids and mainly plays a role in balancing the renin-angiotensin system (RAS) through the conversion of angiotensin II (Ang II) to angiotensin 1–7 . ACE2 receptor is trafficked to the cell surface through the secretory pathway where it is initially synthesized in the endoplasmic reticulum (ER). Once properly folded and post-translationally modified, ACE2 is transported to the Golgi apparatus for further complex posttranslational modifications and folding and then transported to the plasma membrane by vesicular transport . In the ER, tagged misfolded proteins are directed for proteasomal degradation via the endoplasmic reticulum associated degradation (ERAD) .
By the beginning of the spread of COVID-19 infection in 2019, Cao et al. and other groups have demonstrated that ACE2 expression levels and genetic variation may influence its interaction with the SARS-CoV-2 spike (S) protein [8, 12, 13]. Therefore, the highlighted variable outcomes may explain some of the interindividual variability of the infection onset, susceptibility, and severity . Different genetic variants were accounted to mainly affect ACE2 binding affinity to the viral S protein, its expression level, or internalization but little is known regarding its effect on the receptor processing and trafficking through the secretory pathway to the plasma membrane . Some missense variants in secretory proteins like plasma membrane proteins and receptors are known for their deleterious effects and disease causation at various levels [14, 15]. Amino acid substitutions lying away from critically and functionally important protein domains may indirectly result in a loss of function effect due to total or partial retention of the protein in the ER and thus mis-trafficking [14, 16,17,18]. Despite their possible intact biological function, mis-localized membranous proteins lose their function due to their quantitative or partial loss from their distinct functional cellular location. In fact, the Q1069R missense mutation in the ACE2 homologue, angiotensin-converting enzyme (ACE), was found to be sequestered by the ER quality control machinery and prevented from trafficking to the cell surface . Therefore, we hypothesize that some missense variants in ACE2 receptor might exert trafficking defects on this receptor and its levels at the plasma membrane. Moreover, considering ACE2 critical biological functions, partial ER retention or delay might potentially explain some of the interindividual COVID-19 clinical variability and may provide a drug target for SARS-CoV-2 and other coronaviruses infections [1, 19].
In their comparative genetic analysis at the time of initiating this study, Cao et al. have pointed out 30 missense variants out of total 62 genetic variations in the ACE2 coding region for their potential effect on the protein amino acid sequence (Fig. 1A) . We decided then to generate all those variants and evaluate their effects on the subcellular localization and N-glycosylation profile of the ACE2 receptor. Our findings indicate very limited or no detectable effects of these variations on the subcellular localization of ACE2 which may augment the notion that the biological functions of this receptor are essential, and its partial loss cannot be tolerated.
In silico prediction of the structural effects of ACE2 variants
ACE2 reference SNP cluster ID (rsID) with their global allele frequency was retrieved from Ensembl database (https://asia.ensembl.org/Homo_sapiens/Transcript/Variation_Transcript/Table?db=core;g=ENSG00000130234;r=X:15494566-15607236;t=ENST00000427411). To assess the effect of 28 ACE2 nonsynonymous variants on protein function, different in silico prediction tools have been used. SIFT (Sorting Intolerant from Tolerant) algorithm (https://sift.bii.a-star.edu.sg) was utilized to predict whether the studied ACE2 variants affect ACE2 protein function. If SIFT’s score < 0.05 the variant is considered tolerated, and if score > 0.05 the variant is considered to affect protein function . PolyPhen-2 (Polymorphism Phenotyping v2) algorithm (http://genetics.bwh.harvard.edu/pph2/index.shtml) was also used to predict the impact of ACE2 missense variants on its structure and function using different sequence and structure-based predictive features. Two pairs of datasets, HumanDiv and HumanVar, were used to evaluate the consequent damage. PolyPhen-2 scores range between 0 and 1. Benign variants have scores in the range of 0 and 0.15, possibly damaging variants have a score in the range of 0.15 and 0.9 and confidently damaging variants have a score between 0.9 and 1 . PROVEAN (Protein Variation Effect Analyzer) software (http://provean.jcvi.org/seq_submit.php) was also used to predict whether ACE2 amino acid substitutions influence its biological function , where a score less than − 2.5 correspond to a deleterious variant. Finally, to evaluate whether these variants might be disease causing, Mutation Taster has been applied (https://www.mutationtaster.org/MutationTaster69/index.html) . If a variant was found to affect protein function in at least one of the prediction tools, it is classified as possibly deleterious, otherwise it is benign.
Analysis of the effect of missense variants on protein stability and their impact on protein structure using protein modeling
I-Mutant tool was employed to predict the effect of selected missense variants on the protein stability . For this, the protein sequence of ACE2 and missense variants data was submitted in FASTA format. The primary amino acid sequence of ACE2 obtained from the UniProt (Accession No. Q9BYF1) and 28 missense variants generated manually were used as input. Moreover, the coordinates of ACE2 crystal structure were obtained from the protein data bank (PDB), PDB ID 6M17, excluding the coordinates of receptor binding domain (RBD) of SARS-CoV-2 from the co-complexed BoAT1 dimer. I-TASSER was used to model the missing C terminal residues (769–805) and N terminal residues (1–19) of the ACE2 transmembrane helices . SWISS-MODEL was used to generate the three-dimensional homology models of mutant proteins using modeled full-length structure of ACE2 as a template .
Mutagenesis primers design and generation of ACE2 missense variants by site-directed mutagenesis
FLAG-tagged Human ACE2 wild type plasmid (NM_021804) was purchased from OriGene Inc. (RC208442). Among Cao et al.’s reported variants, we have generated 28 missense variants using the wild type construct as a template, via Quick-Change site-directed mutagenesis kit with the Pfu Ultra High-Fidelity DNA polymerase (Stratagene). The primers used for the mutagenesis were designed (Additional file 1: Table S1) using PrimerX software (https://www.bioinformatics.org/primerx/) and purchased from Metabion International AG (https://www.metabion.com/). The generation of the desired variants was confirmed by the dideoxy Sanger DNA sequencing using the ABI 3130xl automated fluorescent Genetic Analyzer (Applied Biosystems). Clustal Omega software (https://www.ebi.ac.uk/Tools/msa/clustalo/) was used for sequence alignments.
HeLa and HEK293 cells were cultured in Dulbecco’s Modified Eagle Medium (Gibco) supplemented with 10% fetal bovine serum (Gibco), antibiotic–antimycotic (Gibco) at 37 °C and 5% CO2 as previously described .
Immunofluorescence and confocal microscopy
HeLa cells were seeded on sterilized cover slips for imaging. Cells were co-transfected with WT or mutated ACE2 plasmid and GFP-tagged HRas, a plasma membrane marker. The methodology has been described previously . Twenty-four hours later, cells were washed three times with phosphate-buffered saline (PBS) and then fixed with methanol at − 20 °C for 5 min. Fixed cells were then blocked with 3% bovine serum albumin (Sigma-Aldrich) for 30 min at room temperature. Fixed cells were co-stained with Anti-Flag primary antibody (1:100 Cell Signaling) and anti-Calnexin (1:50 Santa Cruz Biotechnology) for 1 Hrs in the dark at room temperature. Cells were then washed three times with PBS and then incubated with the respective secondary antibody (Thermo Fischer Scientific) for 45 min in the dark at room temperature. Afterward, cells were then washed and mounted with immmunofluor medium (ICN Biomedicals) and images were acquired using the 100 × objective Nikon confocal Eclipse 80I microscope (Nikon Instruments Inc.). Images were further analyzed and merged using ImageJ software .
Forty-eight hours post transfection, HEK293 cells seeded in 6 well-plates were lysed according to manufacturers’ instructions in RIPA lysis buffer (Pierce Inc.) along with protease inhibitor cocktail (Pierce). Total proteins were quantified by the colorimetric bicinchoninic acid protein assay (BCA kit, Pierce). 20 µg total protein lysate was resolved on 4–12% SDS-PAGE precast gradient gel (GeneScript) followed by transfer into PVDF membrane. Membrane was then incubated with primary antibodies: anti-ACE2 (1:1000 Santa Cruz, cat# sc-390851) anti-Flag antibody (1:1000 Cell Signaling, cat # 8146S), anti-GFP (1:1000 Cell Signaling, cat# 2955S) and anti-actin (1:1000 Santa Cruz Biotechnology, cat# sc-47778) and their corresponding secondary antibodies (Sigma-Aldrich). Membranes were then incubated with Enhanced Chemiluminescence Plus reagent (Pierce) and developed using the Typhoon FLA 9500 imager (GE Healthcare Biosciences, Piscataway, NJ, USA). Blot analysis quantification was then performed using ImageJ software .
Glycosylation sensitivity and resistance assays
HEK293 cells seeded in 6-well plates were co-transfected with wild type or mutant ACE2 and GFP plasmids. Forty-eight hours later, cells were harvested, and lysates were then denatured in denaturation buffer at 100 °C for 10 min based on manufacturer’s protocol. Equal amounts of the proteins were incubated at 37 °C for 3 Hrs in presence or absence of 10U of endoglycosidase H (Sigma-Aldrich). Samples were then resolved on 4–12% SDS-PAGE gel and processed for western blotting as previously described.
N-linked oligosaccharides were removed by PNGase F treatment (New England Biolabs). Cell lysates were denatured at 100 °C for 10 min and equal amounts of glycoproteins were then incubated at 37 °C for 1 Hrs in presence and absence of PNGase F enzyme. Samples were then resolved on 4–12% SDS-PAGE gradient gels and proceeded for western blotting.
Protein stability analysis and half-life determination
Twenty-four hours post transfection, HEK293 cells seeded in 6 well-plates were treated with 100 μg/ml cycloheximide (CHX) (Sigma-Aldrich) to stop new protein translation for different time intervals (4, 8, 12, 18, and 24 Hrs). DMSO-treated cells at the same time intervals were taken as control. Cells were then harvested and proceeded by western blotting.
In silico analysis of ACE2 naturally occurring missense variants
All 28 ACE2 missense variants were tested using different bioinformatic predictive tools to identify the possible functional effects of these variants on ACE2. Altogether, only 3 variants (R768W, G575V, and G173S) were predicted to be deleterious by all the evaluated algorithms (Table 1). Twelve variants were detected to affect protein function by SIFT algorithm. Only 4 variants were found to be deleterious by PROVEAN analysis and 13 were predicted to be damaging by PolyPhen-2 HumanVar, and PolyPhen-2 HumanDiv. Predictions by Mutation Taster show that 10 of the studied substitutions might be disease causing where all others have no predicted phenotypic effect and are considered polymorphic. In total, 10 variants (V801G, N720D, E668K, V658I, F592L, N159S, N149S, D38E, K26R, and I21T) were found to have no effect on the function of ACE2 by all the tested algorithms, classified as benign (Additional file 1: Table S2).
Analysis of the structural stability and effect of missense variants on ACE2 protein
I-mutant Suite was performed to evaluate the effect of missense variants on the overall protein stability. Results demonstrated that all studied variants decreased the ACE2 stability, compared to wild type ACE2 except for S19P, D38Q and S692P variants (Table 2). Among the 28 studied variants, I21T, N33D, F592L, D785N, and V801G were found to introduce the highest instability in ACE2 with a Gibbs free energy change value (ΔΔG) of − 2.16, − 2.10, − 2.97, − 2.16, and − 3.20, respectively.
Moreover, protein modeling was performed to determine the effect of missense variants on ACE2 structure and function. The full-length ACE2 protein was modeled, and missense variants were generated on the protein structure (Fig. 1B). Among 28 studied variants, 13 are located in peptidase domain (PD) of ACE2 which is a potential binding site for RBD of the S protein of SARS-CoV-2. Although some variants are present distant from the ACE2 interface and do not form direct interactions with the RBD of SARS-CoV-2, they might exert structural deformities on ACE2 protein. The evaluation of S19P indicated that the substitution of conserved serine into a hydrophobic proline is likely to disturb the interactions with other ACE2 residues important for protein structure (Fig. 2A1, A2). The substitution of non-polar isoleucine to a polar threonine at 21 amino acid position is likely to produce additional hydrogen bonds with A25 and E87 residues of ACE2 which could disrupt protein conformation (Fig. 2B1, B2). Moreover, the physiochemical changes of K26R and N33D variants are significant and have been shown to enhance the binding affinity of SARS-CoV-2 (Fig. 2C1–D2). However, the substitution of charged residues at amino acid position 38 is less likely to cause significant changes in protein conformation and SARS-CoV-2 binding affinity (Fig. 2E1, E2). Asparagine for serine substitutions at positions 149 and 159 are physiochemically insignificant and may have no effect on protein structure or function (Fig. 2F1–G2). G173 is buried in the core of ACE2 and is also important in forming the catalytic site. G173S variant is likely to disrupt the catalytic site formation as well as catalysis and substrate specificity (Fig. 2H1, H2). The physiochemical properties of M383I, I468V and G575l variants are insignificant and are unlikely to disrupt protein structure and function (Fig. 2I1–L2). A501T variant is expected to disturb the polar interaction of alanine with Glu181 (Fig. 2K1, K2).
Out of 11 variants present in the collectrin domain, five variants (A627V, N638S, V658I, I727V, and L731I) are physiochemically insignificant and are not observed to affect protein structure and function (Fig. 2N1–W2). E668K variant is likely to disrupt the polar interaction of wild type glutamic acid with E667 which could be important for protein conformation (Fig. 2Q1, Q2). The substitution of polar serine with a hydrophobic proline residue at 692 position is likely to disrupt the hydrogen bonds of serine with N159 and N690 (Fig. 2R1, R2). The substitution of charged arginine with neutral tryptophan at 708 position is expected to disrupt the hydrogen bonds and salt bridges formed by R708 with D719, L722, E723, and I727 residues of ACE2 (Additional file 2: Fig. S1, S2). Such changes of S692P and R708W could lead to improper protein folding. Residues 636–658 and 708–717 present at the neck region are observed to form stable polar interactions. These interactions are important for stable dimer formation. The possible loss of charge in case of R710H variant is likely to disturb these interactions, which could result in improper dimer formation (Fig. 2T1, T2). I753T, R763W, D785N, and V801G variants are located at the cytoplasmic end of ACE2 (Fig. 2Y1–Z6). The substitutions are physiochemically significant, particularly R763W, but given their location and the lack of obvious intramolecular interactions, these variants are unlikely to affect protein structure or function.
Exogenously expressed FLAG-tagged wild type ACE2 localizes to the plasma membrane and has a half-life of about 12 Hrs
ACE2 endogenous protein expression in the human HEK293 and HeLa cell lines was assessed by immunoblotting assay using anti-ACE2 monoclonal antibody (Santa Cruz; 1:1000 dilution, cat# sc-390851). Our blot shows that ACE2 protein was not detectable in either HEK293 or HeLa cell lines even when high amounts of cell lysates (50 μg) were analyzed (Additional file 2: Fig. S1A). However, overexpressed Flag-tagged ACE2 in HEK293 cells was detected with the same anti-ACE2 specific antibody at ~ 120 KDa molecular weight.
To establish the intracellular localization of wild type overexpressed ACE2, HeLa cells were co-transfected with Flag-tagged ACE2 and GFP-tagged HRas plasmids. Confocal microscopy images of stained cells displayed a plasma membranal profile for WT ACE2 which overlapped with GFP-tagged HRas as shown in Fig. 3A. Red staining of WT ACE2 by anti-Flag antibody, colocalizes with the green GFP-tagged HRas, staining the plasma membrane, where no colocalization with the blue stained ER marker, Calnexin, was displayed. Similarly, overexpressed ACE2 displays a similar pattern in HEK293 (Additional file 2: Fig. S1B). To confirm the trafficking of ACE2 across the secretory pathway, the N-glycosylation profile of WT ACE2 was evaluated. Glycoprotein N-glycosylation is an enzyme directed process that occurs at specific asparagine residue in N-X-T/S sequence motif, in the ER that is then further modified in the Golgi apparatus. ACE2 is reported to have 7N-glycosylation sites as shown in Fig. 1A [29, 30]. Digestion of protein lysates from HEK293 cells overexpressing WT ACE2 with Endoglycosidase H (Endo H) enzyme, which cleaves off the immature N-glycans only, showed that around 2.8% ± 0.36 of the protein is digestible to a lower molecular weight band (~ 100 KDa) suggestion high level of maturation of the WT protein. The remaining 97.2% did not change their molecular weight protein band of ~ 120KDa suggesting that it has acquired the complex N-glycans that usually take place in the Golgi complex before trafficking to the plasma membrane (Fig. 3B). On the other hand, treatment of N-glycans of ACE2 WT protein by the PNGase F enzyme, a peptide N-glycosidase F that cleaves off N-glycans regardless of their glycans maturation stage, resulted in shift from a high molecular weight band (120 KDa) to a lower molecular weight protein band (~ 100KDa). These results demonstrate that WT ACE2 acquires fully mature N-glycans quantitatively and presumably has high maturation rate.
To gain insight on the half-life, stability, and turnover of the overexpressed WT ACE2, we overexpressed the protein for 24 Hrs then added the protein synthesis inhibitor cycloheximide at 100 µg/ml and quantified the remaining ACE2 by Western blotting at several subsequent time points up to 18 Hrs. As shown in Fig. 3C, D, overexpressed WT ACE2 has an approximate half-life of 12 Hrs suggesting a relatively slow turn over.
All studied ACE2 variants traffic normally and localize to the plasma membrane resembling the WT protein
To assess the effects of the studied ACE2 missense variants, we expressed them individually in HeLa and HEK293 cell lines. HEK293 cells overexpressing ACE2 mutants display normal protein expression pattern by western blotting compared to WT ACE2 (Additional file 2: Fig. S2). HeLa cells were co-transfected with an ACE2 construct (WT or mutants) and GFP-tagged HRas plasmids then the subcellular localization of the proteins was evaluated by confocal microscopy. Interestingly, as shown in Fig. 4 and in Additional file 2: Fig. S3, none of these variants (including the possibly deleterious variants) affected the apparent subcellular localization of the ACE2 protein. For all the studied variants, ACE2 appeared at the plasma membrane and co-localized with GFP-tagged HRas.
To confirm these results, the N- glycosylation profiles of all the studied variants were assessed by Endo H digestion sensitivity and resistance assay and compared to WT ACE2. Similarly, all mutants have shown no significant change in the electrophoretic behaviors compared to WT ACE2, where the vast majority of the expressed ACE2 appears to be resistant to Endo H treatment suggesting that they are fully glycosylated and have normal maturation levels. The Endo H resistant band at ~ 120 KDa accounts for over 96% of the expressed protein for all studied variants. Whereas the lower molecular weight band account to less than ~ 2.8% of total expressed ACE2 as shown in Fig. 5.
The current ongoing COVID-19 pandemic created an urgent need for deciphering the interlink between SARS-CoV-2 and its primary cellular receptor, ACE2. Several studies have extensively investigated ACE2 polymorphic footprint and its associated effects on its structure, binding, and stability [8, 9, 31, 32]. Genetic variations in ACE2 gene are regarded as a potential risk factor in COVID-19 patients . In this context, different predictive studies based on bioinformatics and simulation tools have generated a bulk of data that helped identify major residues in ACE2 and their consequent effect on SARS-CoV-2 binding [34,35,36].
In the current study, we have evaluated the impact of several ACE2 coding missense variants using different predictive algorithms and investigated their effect on ACE2 protein subcellular localization, trafficking, and membrane availability. ACE2 gene displays a unique polymorphic profile in the human population in which 332 missense variants were reported in Ensembl database. Interestingly, in comparison to its homologue the angiotensin-converting enzyme (ACE), ACE2 displays a lower probability of losing its function by genetic mutations, where probability of being loss-of-function intolerant of ACE2 is pLI = 0.998 by gnomAD database, noting that if pLI > 0.9 gene is considered extremely intolerant [31, 37]. Unlikely, a missense variant in ACE protein (Q1069R) was reported to be associated with renal tubular disease resulting in premature death caused by improper localization of ACE and loss of its function . Although different global or conditioned knockout mouse models have been generated for ACE2, all have reported serious associated phenotypes including developmental, cardiovascular, renal and respiratory clinical manifestations .
Overexpressed ACE2 in HEK293 cells in our study was detected at a similar molecular weight ~ 120 KDa compared to endogenous ACE2 detected in other cell lines like HuH7, Caco-2, Calu-3, and HepG2 cells . Moreover, our generated missense variants displayed similar expression profile compared to WT as shown in Additional file 2: Fig. S2. It is worth noting that other ACE2 human variants studied by Shukla et al. and Bhattacharjee et al. (including S19P assessed in our analysis) have displayed no significant change in ACE2 expression level in a similar cellular model [40, 41]. While the surface expression level hasn’t changed, this doesn’t mean that these variants don’t affect the susceptibility of the virus. The binding affinity to the RBD of SARS-CoV-2 was varying significantly compared to WT ACE2 in Bhattacharjee et al. study , considering that the affinity-infectivity relationship is still debatable in various studies .
All studied genetic variants in this paper had a low population distribution with minor allele frequency less than 1%. Although eight of the variants (K26R, I468V, A627V, N638S, S692P, N720D, and L731F/L731I) were considered major hotspot by Cao et al., and were distributed in different populations, two of them, K26R and N720D, were shown in our analysis to be benign by the different evaluated tools. Stability analysis of K26R and N720D variants was decreased by − 0.34 kcal/mol and − 1.38 kcal/mol (Table 2), respectively. Additionally, glycosylation analysis and immunofluorescence assays display no effect on the cellular trafficking and cell surface localization of these two variants (Additional file 2: Fig S2, Fig. S4 and S5). Substitutions of ACE2 residues involved in SARS-CoV-2 binding have been shown to alter the binding affinity of SARS-CoV-2 S protein [6, 43]. K26 residue is present in the binding domain with SARS-CoV-2 receptor RBD, several studies suggest that residual substitution of lysine with arginine at this site enhances ACE2-SARS-CoV-2 binding and could contribute to higher susceptibility to SARS-CoV-2 [6, 34, 40, 44]. Unlikely, N720 amino acid residue is located in the peptidase domain in a proximity to type II transmembrane serine protease (TMPRSS2) cleavage site. Aside from having lower stable structure compared to WT ACE2, substituting asparagine with aspartic acid at this residue is found to weaken TMPRSS2-ACE2 complex and consequently augmenting SARS-CoV-2 viral entry [31, 45]. Although our analysis displays benign effect of these two variants, and no change in the trafficking and expression profiles, it still can affect their susceptibility to bind to the spike protein. Our analysis has also revealed three deleterious variants (R768W, G575V, and G173S) that were found to affect ACE2 functionality by all the prediction tools in our study. Similarly, intracellular localization and cell surface availability of these variants were not affected in our studied model. Knowing that G575 and G173 are located in the peptidase domain of ACE2 protein, G173S variant is predicted to stabilize ACE2-SARS-CoV-2 complex , noting that this substitution might affect ACE2 catalytic activity, while the effect of G575V on the binding affinity to SARS-CoV-2 RBD is not investigated yet. Our results come in line with another in silico analysis showing that R768W is a high risk ACE2 variant and might exert a deleterious effect on ACE2 structure . Additionally, previous studies have shown that the 43 long amino acid topological domain of ACE2 doesn’t affect significantly the cell surface expression and SARS-CoVs mediated viral entry [47, 48]. These data fit with our reported results showing that all mutants present in the cytoplasmic tail of ACE2 (V801G, D785N and R768W) display no change in ACE2 trafficking and transmembrane localization.
In this context, little has been reported in literature about ACE2 biosynthesis and intracellular trafficking. Evident role of N-linked glycosylation on membrane proteins stability, compartmental trafficking, and cell surface expression has been largely reported [49,50,51]. Among the seven glycosylation sites of ACE2, N90, N322 and N546 form glycan-mediated interactions between ACE2 and SARS-CoV-2 RBD complex , noting that N90 and N322 have opposing effects on SARS-CoV-2 binding . ACE2 deprived of all N-glycans was found to accumulate in the endoplasmic reticulum with no significant effect on its enzymatic carboxypeptidase activity . N-linked glycans are usually attached to an asparagine residue of N-X-S/T motif, where X is any amino acid except proline and S/T are serine/threonine amino acids. Among the variants we were interested in studying is S692P, two amino acids proximal to N690 glycosylation site. Substituting serine for proline in the S/T site would consequently lead to the loss of N-glycosylation at N690. We show in our analysis that this substitution increases ACE2 stability (ΔΔG = 1.41 kcal/mol), however, immunofluorescence data show no change in the trafficking of this variant compared to the WT ACE2.
In summary, throughout this study, we show that none of the coding variants included in our work display an effect on ACE2 protein subcellular trafficking and its cell surface availability. This might be due to the high importance of this gene and therefore its intolerance to the loss of its function. Noting that the deleterious effect reported by the computational tools might be affecting the carboxypeptidase activity and binding affinity of ACE2 that were not evaluated in our paper, rather than its trafficking and SARS-CoVs related effect. Apart from ACE2 trafficking modulation, the receptor membranous expression and availability might be influenced at different levels . Considering that ACE2 gene lies on the X chromosome, it is strongly evident that females have an advantage due to the greater chance to form structurally variable ACE2 dimers compared to males who are hemizygous for the ACE2 gene . ACE2 expression is highly affected by sex hormones where estrogen strongly elevates its expression explaining the variable effect of COVID-19 in both sexes . Moreover, ACE2 expression might also be affected by different epigenetic changes providing some evidence to explain the difference in the interindividual susceptibility and clinical variability in infected patients [56, 57]. It is worth noting that a recent Mendelian Randomization (MR) study suggests that ACE2 expression is a causative player in the susceptibility and severity of SARS-CoV-2 infection . Although ACE2 genetic variants didn’t display any significant effect on ACE2 cell surface availability through its trafficking pathway, it still might affect the pharmacokinetics of ACE2 including its decay and consequently SARS-CoV-2 entry. Collectively, these data have important implication on COVID-19 progression and to face these outcomes more detailed knowledge is needed to understand the normal mechanisms controlling ACE2 expression and trafficking at the cellular level as well as its role in the pathogenesis of SARS-CoVs infection.
Availability of data and materials
All data generated or analyzed during the study are included in this article and its additional files. Any further requirements are available from the corresponding author upon request.
Angiotensin-converting enzyme 2
- Ang II:
Corona Virus Disease-2019
- Endo H:
Endoplasmic reticulum-associated degradation
Probability of being loss-of-function intolerant
- PNGase F:
Peptide N-glycosidase F
Receptor binding domain
Severe acute respiratory syndrome corona viruses
Single nucleotide polymorphism
Type II transmembrane serine protease
Badawi S, Ali BR. ACE2 nascence, trafficking, and SARS-CoV-2 pathogenesis: the saga continues. Hum Genomics. 2021;15(1):8.
Bai Y, Yao L, Wei T, Tian F, Jin DY, Chen L, et al. Presumed asymptomatic carrier transmission of COVID-19. JAMA. 2020;323(14):1406–7.
Zhang J, Wang X, Jia X, Li J, Hu K, Chen G, et al. Risk factors for disease severity, unimprovement, and mortality in COVID-19 patients in Wuhan, China. Clin Microbiol Infect. 2020;26(6):767–72.
Richardson S, Hirsch JS, Narasimhan M, Crawford JM, McGinn T, Davidson KW, et al. Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the New York City area. JAMA. 2020;323(20):2052–9.
Pereira NL, Ahmad F, Byku M, Cummins NW, Morris AA, Owens A, et al. COVID-19: understanding inter-individual variability and implications for precision medicine. Mayo Clin Proc. 2021;96(2):446–63.
Suryamohan K, Diwanji D, Stawiski EW, Gupta R, Miersch S, Liu J, et al. Human ACE2 receptor polymorphisms and altered susceptibility to SARS-CoV-2. Commun Biol. 2021;4(1):1–11.
Jia H, Neptune E, Cui H. Targeting ACE2 for COVID-19 therapy: opportunities and challenges. Am J Respir Cell Mol Biol. 2021;64(4):416–25.
Cao Y, Li L, Feng Z, Wan S, Huang P, Sun X, et al. Comparative genetic analysis of the novel coronavirus (2019-nCoV/SARS-CoV-2) receptor ACE2 in different populations. Cell Discov. 2020;6(1):1–4.
Benetti E, Tita R, Spiga O, Ciolfi A, Birolo G, Bruselles A, et al. ACE2 gene variants may underlie interindividual variability and susceptibility to COVID-19 in the Italian population. Eur J Hum Genet. 2020;28(11):1602–14.
Novelli A, Biancolella M, Borgiani P, Cocciadiferro D, Colona VL, D’Apice MR, et al. Analysis of ACE2 genetic variants in 131 Italian SARS-CoV-2-positive patients. Hum Genomics. 2020;14(1):29.
Liu X, Yang N, Tang J, Liu S, Luo D, Duan Q, et al. Downregulation of angiotensin-converting enzyme 2 by the neuraminidase protein of influenza A (H1N1) virus. Virus Res. 2014;24(185):64–71.
Saih A, Baba H, Bouqdayr M, Ghazal H, Hamdi S, Kettani A, et al. In Silico analysis of high-risk missense variants in human ACE2 gene and susceptibility to SARS-CoV-2 infection. BioMed Res Int. 2021;9(2021):e6685840.
Chen F, Zhang Y, Li X, Li W, Liu X, Xue X. The impact of ACE2 polymorphisms on COVID-19 disease: susceptibility, severity, and therapy. Front Cell Infect Microbiol. 2021. https://doi.org/10.3389/fcimb.2021.753721.
Kizhakkedath P, John A, Al-Sawafi BK, Al-Gazali L, Ali BR. Endoplasmic reticulum quality control of LDLR variants associated with familial hypercholesterolemia. FEBS Open Bio. 2019;9(11):1994–2005.
Fukuda R, Okiyoneda T. Cystic fibrosis transmembrane conductance regulator (CFTR) ubiquitylation as a novel pharmaceutical target for cystic fibrosis. Pharmaceuticals. 2020;13(4):75.
Mohamed FE, Al Sorkhy M, Ghattas MA, Al-Gazali L, Al-Dirbashi O, Al-Jasmi F, et al. The pharmacological chaperone N-n-butyl-deoxygalactonojirimycin enhances β-galactosidase processing and activity in fibroblasts of a patient with infantile GM1-gangliosidosis. Hum Genet. 2020;139(5):657–73.
Needham PG, Guerriero CJ, Brodsky JL. Chaperoning endoplasmic reticulum-associated degradation (ERAD) and protein conformational diseases. Cold Spring Harb Perspect Biol. 2019;11(8):a033928.
Gariballa N, Ali BR. Endoplasmic reticulum associated protein degradation (ERAD) in the pathology of diseases related to TGFβ signaling pathway: future therapeutic perspectives. Front Mol Biosci. 2020;7:575608.
de Oliveira RM, Marijanovic Z, Carvalho F, Miltényi GM, Matos JE, Tenreiro S, et al. Impaired proteostasis contributes to renal tubular dysgenesis. PLoS ONE. 2011;6(6):e20854.
Vaser R, Adusumalli S, Leng SN, Sikic M, Ng PC. SIFT missense predictions for genomes. Nat Protoc. 2016;11(1):1–9.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9.
Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and Indels. PLoS ONE. 2012;7(10):e46688.
Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat Methods. 2014;11(4):361–2.
Capriotti E, Fariselli P, Casadio R. I-Mutant20: predicting stability changes upon mutation from the protein sequence or structure. Nucl Acids Res. 2005;1(33):W306-310.
Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008;23(9):40.
Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46(W1):W296-303.
Hume AN, Buttgereit J, Al-Awadhi AM, Al-Suwaidi SS, John A, Bader M, et al. Defective cellular trafficking of missense NPR-B mutants is the major mechanism underlying acromesomelic dysplasia-type Maroteaux. Hum Mol Genet. 2009;18(2):267–77.
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of Image Analysis. Nat Methods. 2012;9(7):671–5.
Rowland R, Brandariz-Nuñez A. Analysis of the role of N-linked glycosylation in cell surface expression, function, and binding properties of SARS-CoV-2 receptor ACE2. Microbiol Spectr. 2021;9(2):e01199-21.
Mehdipour AR, Hummer G. Dual nature of human ACE2 glycosylation in binding to SARS-CoV-2 spike. Proc Natl Acad Sci USA. 2021;118(19):e2100425118.
Al-Mulla F, Mohammad A, Al Madhoun A, Haddad D, Ali H, Eaaswarkhanth M, et al. ACE2 and FURIN variants are potential predictors of SARS-CoV-2 outcome: a time to implement precision medicine against COVID-19. Heliyon. 2021;7(2):e06133.
Ortiz-Fernández L, Sawalha AH. Genetic variability in the expression of the SARS-CoV-2 host cell entry factors across populations. Genes Immun. 2020;21(4):269–72.
Hou Y, Zhao J, Martin W, Kallianpur A, Chung MK, Jehi L, et al. New insights into genetic susceptibility of COVID-19: an ACE2 and TMPRSS2 polymorphism analysis. BMC Med. 2020;15(18):216.
Ren W, Zhu Y, Lan J, Chen H, Wang Y, Shi H, et al. Susceptibilities of human ACE2 genetic variants in coronavirus infection. J Virol. 2022;96(1):e0149221.
Sorokina M, Teixeira J, Barrera-Vilarmau S, Paschke R, Papasotiriou I, Rodrigues JPGLM, et al. Structural models of human ACE2 variants with SARS-CoV-2 Spike protein for structure-based drug design. Sci Data. 2020;7(1):309.
Heinzelman P, Romero PA. Discovery of human ACE2 variants with altered recognition by the SARS-CoV-2 spike protein. PLoS ONE. 2021;16(5):e0251585.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581(7809):434–43.
Jia H, Yue X, Lazartigues E. ACE2 mouse models: a toolbox for cardiovascular and pulmonary research. Nat Commun. 2020;11(1):5165.
Sherman EJ, Emmer BT. ACE2 protein expression within isogenic cell lines is heterogeneous and associated with distinct transcriptomes. Sci Rep. 2021;11(1):15900.
Shukla N, Roelle SM, Suzart VG, Bruchez AM, Matreyek KA. Mutants of human ACE2 differentially promote SARS-CoV and SARS-CoV-2 spike mediated infection. PLOS Pathog. 2021;17(7):e1009715.
Bhattacharjee MJ, Lin JJ, Chang CY, Chiou YT, Li TN, Tai CW, et al. Identifying primate ACE2 variants that confer resistance to SARS-CoV-2. Mol Biol Evol. 2021;38(7):2715–31.
MacGowan SA, Barton MI, Kutuzov M, Dushek O, van der Merwe PA, Barton GJ. Missense variants in human ACE2 strongly affect binding to SARS-CoV-2 Spike providing a mechanism for ACE2 mediated genetic risk in Covid-19: case study in affinity predictions of interface variants. PLoS Comput Biol. 2022;18(3):e1009922.
Chan KK, Dorosky D, Sharma P, Abbasi SA, Dye JM, Kranz DM, et al. Engineering human ACE2 to optimize binding to the spike protein of SARS coronavirus 2. Science. 2020;369(6508):1261–5.
Barton MI, MacGowan SA, Kutuzov MA, Dushek O, Barton GJ, Van der Merwe PA. Effects of common mutations in the SARS-CoV-2 Spike RBD and its ligand, the human ACE2 receptor on binding affinity and kinetics. eLife. 2021;10:e70658.
Mohammad A, Marafie SK, Alshawaf E, Abu-Farha M, Abubaker J, Al-Mulla F. Structural analysis of ACE2 variant N720D demonstrates a higher binding affinity to TMPRSS2. Life Sci. 2020;15(259):118219.
Hadi-Alijanvand H, Rouhani M. Studying the effects of ACE2 mutations on the stability, dynamics, and dissociation process of SARS-CoV-2 S1/hACE2 complexes. J Proteome Res. 2020;19(11):4609–23.
Inoue Y, Tanaka N, Tanaka Y, Inoue S, Morita K, Zhuang M, et al. Clathrin-dependent entry of severe acute respiratory syndrome coronavirus into target cells expressing ACE2 with the cytoplasmic tail deleted. J Virol. 2007;81(16):8722–9.
Karthika T, Joseph J, Das VRA, Nair N, Charulekha P, Roji MD, et al. SARS-CoV-2 cellular entry is independent of the ACE2 cytoplasmic domain signaling. Cells. 2021;10(7):1814.
Singh R, Almutairi MM, Pacheco-Andrade R, Almiahuob MYM, Di Fulvio M. Impact of hybrid and complex N-glycans on cell surface targeting of the endogenous chloride cotransporter Slc12a2. Int J Cell Biol. 2015;2015:505294.
Weng TY, Chiu WT, Liu HS, Cheng HC, Shen MR, Mount DB, et al. Glycosylation regulates the function and membrane localization of KCC4. Biochim Biophys Acta. 2013;1833(5):1133–46.
Wang T, Nakagawa S, Miyake T, Setsu G, Kunisue S, Goto K, et al. Identification and functional characterisation of N-linked glycosylation of the orphan G protein-coupled receptor Gpr176. Sci Rep. 2020;10(1):4429.
Zhao P, Praissman JL, Grant OC, Cai Y, Xiao T, Rosenbalm KE, et al. Virus-receptor interactions of glycosylated SARS-CoV-2 spike and human ACE2 receptor. Cell Host Microbe. 2020;28(4):586-601.e6.
Beyerstedt S, Casaro EB, Rangel ÉB. COVID-19: angiotensin-converting enzyme 2 (ACE2) expression and tissue susceptibility to SARS-CoV-2 infection. Eur J Clin Microbiol Infect Dis. 2021;3:1–15.
Gemmati D, Bramanti B, Serino ML, Secchiero P, Zauli G, Tisato V. COVID-19 and individual genetic susceptibility/receptivity: role of ACE1/ACE2 genes, immunity, inflammation and coagulation: might the double X-chromosome in females be protective against SARS-CoV-2 compared to the single X-chromosome in males? Int J Mol Sci. 2020;21(10):3474.
Chen J, Jiang Q, Xia X, Liu K, Yu Z, Tao W, et al. Individual variation of the SARS-CoV-2 receptor ACE2 gene expression and regulation. Aging Cell. 2020;19(7):e13168.
Lima RS, Rocha LPC, Moreira PR. Genetic and epigenetic control of ACE2 expression and its possible role in COVID-19. Cell Biochem Funct. 2021. https://doi.org/10.1002/cbf.3648.
Novelli G, Biancolella M, Mehrian-Shai R, Colona VL, Brito AF, Grubaugh ND, et al. COVID-19 one year into the pandemic: from genetics and genomics to therapy, vaccination, and policy. Hum Genomics. 2021;15(1):27.
Zheng M. ACE2 and COVID-19 susceptibility and severity. Aging Dis. 2022;13(2):360.
We would like to thank Dr Saeed Tariq for his assistance in confocal microscopy imaging.
This work was funded by the Abu Dhabi ASPIRE Award for Research Excellence (AARE) (grant AARE19-086, 21M136).
Ethics approval and consent to participate
Consent for publication
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be constructed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Badawi, S., Mohamed, F.E., Alkhofash, N.R. et al. Characterization of ACE2 naturally occurring missense variants: impact on subcellular localization and trafficking. Hum Genomics 16, 35 (2022). https://doi.org/10.1186/s40246-022-00411-1
- Angiotensin-converting enzyme 2 (ACE2)
- Subcellular trafficking