key: cord-0695275-qgmjwmx7 authors: Payandeh, Zahra; Reza Rahbar, Mohammad; Jahangiri, Abolfazl; Sadat Hashemi, Zahra; Zakeri, Alireza; Jafarisani, Moslem; Javad Rasaee, Mohammad; Khalili, Saeed title: Design of an engineered ACE2 as a novel therapeutics against COVID-19 date: 2020-07-29 journal: J Theor Biol DOI: 10.1016/j.jtbi.2020.110425 sha: aa50ac57f43e0a44ae043d00971eceed123b176a doc_id: 695275 cord_uid: qgmjwmx7 The interaction between the angiotensin-converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of the spike protein from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plays a pivotal role in virus entry into the host cells. Since recombinant ACE2 protein has been suggested as an anti-SARS-CoV-2 therapeutic agent, this study was conducted to design an ACE2 protein with more desirable properties. In this regard, the amino acids with central roles in enzymatic activity of the ACE2 were substituted. Moreover, saturation mutagenesis at the interaction interface between the ACE2 and RBD was performed to increase their interaction affinity. The best mutations to increase the structural and thermal stability of the ACE2 were also selected based on B factors and mutation effects. The obtained resulted revealed that the Arg273Gln and Thr445Gly mutation have drastically reduced the binding affinity of the angiotensin-II into the active site of ACE2. The Thr27Arg mutation was determined to be the most potent mutation to increase the binding affinity. The Asp427Arg mutation was done to decrease the flexibility of the region with high B factor. The Pro451Met mutation along with the Gly448Trp mutation was predicted to increase the thermodynamic stability and thermostability of the ACE2. The designed therapeutic ACE2 would have no enzymatic activity while it could bear stronger interaction with Spike glycoprotein of the SARS-CoV-2. Moreover, decreased in vivo enzymatic degradation would be anticipated due to increased thermostability. This engineered ACE2 could be exploited as a novel therapeutic agent against COVID-19 after necessary evaluations. The SARS-CoV-2 is a non-segmented positive sense RNA virus which belongs to the genus Beta Coronavirus within the Coronaviridae family. Its 30 kilo base pairs genome encodes for multiple structural and non-structural proteins. Like SARS and the Middle East Respiratory Syndrome Coronavirus (MERS-CoV), it mainly causes fever, chronic fatigue, dry cough, shortness of breath, and severe respiratory distress (Cabeça et al., 2013; Corman et al., 2018; Enjuanes et al., 2006; Perlman and Netland, 2009; van Boheemen et al., 2012) . The envelope-anchored trimeric spike (S) protein at the surface of the virus is considered to be the key protein for virus entry into the host cells. The ACE2 (a homologue of ACE) is reported to be the potent human receptor for SARS-CoV-2 (Othman et al., 2020) . This receptor primarily plays physiological roles in the control of vasoconstriction and blood pressure. As a type I membrane protein, ACE2 is found to be expressed in lungs, heart, kidneys and intestine (Donoghue et al., 2000; Zhang et al., 2020; Zhao et al., 2020) . It also has been shown that the cardiovascular diseases are associated with reduced ACE2 expression (Crackower et al., 2002) . Given these properties, ACE2 has been served as a drug target for cardiovascular and kidney diseases. Although the RBD of the SARS-CoV-2 spike protein contains 6 amino acid substitutions compared to the RBD of the SARS spike protein, existence of highly conserved residues at this region corroborates the ACE2 as the receptor for SARS-CoV-2 entry (Boehm and Nabel, 2002; Donoghue et al., 2000; Ge et al., 2013) . To determine whether the SARS-CoV-2 uses ACE2 as the cell entry receptor, Peng Zhou et al. conducted a virus infection study using HeLa cells expressing ACE2 protein and cells without ACE2 expression. They have demonstrated that the SARS-CoV-2 can only infect ACE2-containing cells. They also have demonstrated that SARS-CoV-2 does not use other corona virus receptors including aminopeptidase N, dipeptidyl and peptidase 4 (Wu et al., 2020) . Given the crucial role of ACE2 and RBD interaction in the pathogenesis of SARS-CoV-2, inhibitory moieties to halt this interaction would prevent the virus from infecting the host cells. Targeting the interacting interfaces of ACE2 and RBD via antibodies or other therapeutic small molecules has already been practiced in numerous studies (Chen et al., 2017; Coleman et al., 2014; Der Sarkissian et al., 2008; Pang et al., 2020) . Developing these kind of conventional therapeutics is inherently time consuming and demands recourses. An alternative approach to interfere with the interaction between ACE2 and RBD is the application of molecules which could compete with the cell membrane ACE2 to bind RBD. A recombinant ACE2 molecule which is designed to interact with RBD with higher affinity, without enzymatic activity and higher thermostability could be a suitable alternative. Studies shown that specificity, binding affinity, thermostability and immunogenicity of therapeutic proteins can be altered by amino acid substitutions. Since several amino acids at the binding interface of the interacting molecules are responsible for binding affinity and immunogenicity, substitution of such key amino acids can significantly alter the affinity and immunogenicity (Clark et al., 2006; Xu et al., 2009) . Such effective amino acids could be determined in vitro or using bioinformatics tools and replaced with more efficient amino acids. Bioinformatics approaches provide the opportunity to assess the effect of various sequence alterations on protein affinity, thermostability and enzymatic activity by reduced research costs and time (Payandeh et al., 2018; Payandeh et al., 2019; Sefid et al., 2019) . Given the importance of SARS-CoV-2, its transmission and rapid spread to all countries around the world, it is important to develop new therapeutic approaches to deal with it. Since the ACE2 receptor plays the most important role in transmitting the virus into the pulmonary cells, the purpose of this study is to design a human ACE2 inhibitor that lacks its enzymatic activity and bears higher thermostability and higher affinity towards RBD. The protein sequences for RBD of the COVID-19 spike protein and the ACE2 protein were searched in UniProt (https://www.uniprot.org/) and NCBI (https://www.ncbi.nlm.nih.gov) databases. The 3D structures of these proteins along with their interacting complexes were obtained from PDB database. The latest structures with highest resolution were used for further analyses. Since there were no resolved structures for ACE2 and angiotensin 2 (ANG2) complex, the servers. The complex of ACE2 and ANG2 was used to find the most suitable point mutation candidates. The 2D plot of the interactions between the amino acids of ACE2 and ANG2 was generated by LigPlus software (v.2.1). The complex of ACE2 and ANG2 was fed as the input file and the LigPlot tool was employed. The residue interaction network of ACE2 and ANG2 was extracted from the tertiary structure of the complex of ACE2 and ANG2 by Residue Interaction Network Generator (Piovesan et al., 2016) as provided by http://protein.bio.unipd.it/. Strong and reliable interactions were generated by selecting the strict distance threshold. The network policy for identification of interacting pairs was based on calculating the distance between mass centers of two interacting residues. All types of interactions were considered in this analysis. The generated network further analyzed and visualized by Cytoscape 3.7.2 (Shannon et al., 2003) . To investigate the key residues in correct folding and function of the protein, a centrality analysis was carried out by RINspector (Brysbaert et al., 2018) application of Cytoscape. The centrality analysis was based on the change of average shortest path length under removal of each residue (node in the network) as proposed by (del Sol et al., 2006) and embedded in RINSpector. The shortest path length is the minimum number of edges that are required for connecting the first node to the second one in the network. The mean of these numbers, related to all possible pairs of residues in the network of interaction, is defined as average shortest path. The significance of average shortest path length is determined by the calculation of Z-score. The Z-scores greater than 2 are considered as relevant (Brysbaert et al., 2018) . The complex of ACE2 and ANG2 was fed as input for the mCSM-PPI2 server (http://biosig.unimelb.edu.au/mcsm_ppi2/) to analyze the effects of single-point mutations on protein-protein interactions. All possible mutations were performed for each mutation site to find the most deactivating mutation. The results of performed analyses along with the results of previous reports were considered to select the key amino acids which their mutation would terminate the enzymatic activity of the ACE2. The candidate mutations were exerted using the protein preparation tool from Molegro VIRTUAL DOCKER (v.6.0). The mutated amino acids and the neighboring amino acids were energy minimized to arrive at proper conformations. The Molegro VIRTUAL DOCKER was employed to determine the cavity sizes within the ACE2 structure before and after exertion of the mutations. The structures of the mutated ACE2 and ANG2 were docked again using CABSdock and ClusPro servers to analyze the consequences of the mutations. The structure of ACE2 and RBD complex was extracted from the previously resolved protein complex under the PDB ID of 6M17. This complex was fed as the input for the mCSM-PPI2 server (http://biosig.unimelb.edu.au/mcsm_ppi2/) to analyze the effects of single-point mutations on protein-protein interactions. The Saturation Mutagenesis tool of this server was employed to mutate all residues at the interaction interface to all 19 standard amino acids. The effect of each mutation is reported as predicted affinity change (ΔΔG affinity ) (positive values correspond to affinity increase and vice versa). The LigPlus software was employed to draw the 2D plot of the interactions between the amino acids of ACE2 and RBD at their interaction interface. The mutations with highest possibility to increase the binding affinity between the ACE2 and RBD were exerted in the ACE2 structure by Molegro VIRTUAL DOCKER. The mutated amino acids and the neighboring amino acids were energy minimized to arrive at proper conformations. The effects of the mutation were also analyzed by LigPlus software by drawing the 2D plots of ACE2 and RBD interactions after each mutation. Average B factor for each residue was calculated by Chimera software for both 1R42 and 6M17. The results for 19 possible mutations at the residue with highest B factor were done by iStable2.0 software (http://ncblab.nchu.edu.tw/iStable2/result.php). The PoPMuSiC and HotMuSiC (https://soft.dezyme.com/) were respectively used for the prediction of protein thermodynamic stability changes upon single-site mutations and the prediction of protein thermostability changes upon single-site mutations. A systematic prediction was performed to have the prediction results for 19 mutations per ACE2 residue. The best mutations to enhance the thermal stability were exerted and optimized by Molegro VIRTUAL DOCKER. The CABSflex (http://biocomp.chem.uw.edu.pl/CABSflex2) server was employed for or fast simulation of protein structure flexibility. This tool generates protein dynamics at highly reduced (3 orders of magnitude) computational cost, although with some decrease in resolution. The ACE2 structure with 6 mutations was fed as input for flexibility simulation. The MD analysis was used to investigate the stability of the ACE2 protein following the exerted mutations. This method allows the protein atoms to move freely and find the best location based on the applied force-fields, while placed at a specific temperature and pressure conditions. If the mutated protein continues to maintain its structure and stability after the completion of the molecular dynamics process, the structure could be considered as stable. The GROMACS software version 4.6.5 and the GROMOS96-54 A7 force field were used to perform the MD simulation. The complex was put into a cube box and filled with water using the TIP3 model. If required, the system was neutralized by using Na+ or Cl-ions. Finally, the MD simulation was performed for 30 nanoseconds. The Root Mean Square Deviation (RMSD) and radius of gyration (Rg) was used to evaluate the stability of the mutated ACE2 protein. The predicted mutations were exerted in the structure of the ACE2 from 1R42. The interaction orientation for mutant ACE2 and RBD was reanalyzed by molecular docking. The docking analysis was performed using HADDOCK server (https://wenmr.science.uu.nl/haddock2.4/) (Van Zundert et al., 2016) . The results of the docking analyses were superimposed on the original ACE2 and RBD complex and the RMSD value was calculated using Chimera software (v.1.10.2). The sequences of wild type and mutant ACE2 proteins were analyzed by Bepipred linear epitope prediction tool (http://tools.iedb.org/bcell/) of the IEDB database. The glycosylation sites of the wild type and mutant ACE2 proteins were analyzed by NetNGlyc (http://www.cbs.dtu.dk/services/NetNGlyc/) server. The sequences for RBD of the SARS-CoV-2 spike protein and the ACE2 protein were found under the NCBI accession numbers of QHD43416 and NP_001358344, respectively. The 3D structure for interacting complex of the RBD from spike protein and the ACE2 protein was stored under the PDB ID of 6M17. The ACE2 structure was obtained under the PDB ID of the 1R42. The extra protein chains and the non-protein moieties were removed from both 6M17 and 1R42 structures. The RBD of the SARS-CoV-2 spike protein was also extracted from the 6M17 as a separate structure for further analyses. The molecular docking analyses between the ACE2 and ANG2 proteins indicated that the ANG2 could form a tight interaction with the ACE2 protein within the expected central cleavage of the enzyme. Correct deployment of the ANG2 within the ACE2 active site could be interpreted as the suitability of the employed docking method. The 2D contact, purple lines are either covalent bonds between protein and ligand, or "elastic" bonds within the ligand). The ANG2 is located at the center of the plot, while the ACE2 residues are around it. To determine the functional residues in the protein, a centrality analysis was performed on the residue interaction network of the ACE2-ANG2 complex. As expected, the residue interaction network is somehow a complex network (Figure 2 , left panel). To have a more clear view, the active pocket of the enzyme (ACE2), the interface residues between ACE2 (Chain A) and ANG2 were extracted from the whole network ( Figure 2 , upper panel and Table 1 ). Moreover, the significant residues in the term of centrality were separated and illustrated in Figure 2 lower panel. The Arginine 273 and the Threonine 445 of the ACE2 are among the amino acids attributed with high and significant Z-scores. Moreover, the Threonine 445 is shown to have a hydrogen bond with the Arginine 2 of the ANG2 molecule, while the Arginine 273 is shown to have a hydrogen bond with Aspartate 1 of the ANG2 molecule. The Arginine 2 of ANG2 is determined to be the central residue of ANG2. The increased level of Z-score indicates the increase in the average shortest path length upon node removal. When a significant Z-score assigned to a node (residue), it can be inferred that the related node may play an important role in the communication of the network. Therefore, the subjected residue may be an important structural or functional one. The illustration of active pocket is represented in protein structure to get a better insight about the native structure of the enzyme (Figure 3 ). The Arginine 273 and the Threonine 445 of the ACE2 molecule were mutated to all 19 possible amino acids. The results of the mutagenesis analyses for both selected amino acids have indicated that the mutation to Glycine is the best mutation to decrease the binding affinity between the ACE2 and the ANG2 molecules ( Reanalysis of the molecular docking between the mutated ACE2 and ANG2 molecule indicates that following the exerted mutations the ANG2 molecule loses its ability to make a stable interaction with the active site of the mutated enzyme. (Figure 4) . These results indicate that the performed mutations would most likely terminate the enzymatic activity of the ACE2 by removing highly stabilizing bonds which were provided by original amino acids at positions 273 and 445. contact, purple lines are either covalent bonds between protein and ligand, or "elastic" bonds within the ligand). The ANG2 is located at the center of the plot, while the ACE2 residues are around it (b). The results of the Saturation Mutagenesis analyses were obtained and sorted according to their effects to increase the binding affinity between ACE2 and RBD. The best possible mutations are listed in Table 2 . The mCSM-PPI2 models the effects of variations on the inter-residue noncovalent interaction network using graph kernels, evolutionary information, complex network metrics and energetic terms to better assess the molecular mechanism of the mutation. The employed graph-based structural signatures called mCSM has been shown to be an accurate and high-throughput approach to predict the impact of mutations on protein structure and function, and was one of the first methods capable of assessing the impact of mutations on protein interaction binding affinity (Rodrigues et al., 2019) . The important amino acids which are involved in the interaction between ACE2 and RBD were also plotted by LigPlus software. The amino acids which are involved in hydrogen bonds or hydrophobic interactions are depicted in Figure 5 . Analyses of various single or multiple mutations were performed using their 2D interaction plots. The mutations which were able to conserve the original interacting amino acids at the ACE2 and RBD interface and increasing the number of involved amino acids were considered as suitable mutations. The obtained results demonstrated that the threonine 27 to Arginine mutation is the only mutation which met the selection criteria ( Figure 6 ). Obeying these criteria would grantee the preservation of the ACE2 and RBD interactions in their original orientation. Five ACE2 residues with highest average B factors are listed in Table 3 . The aspartate residue at position 427 has been among top 5 residues with high B factor. The region with highest B factor is depicted in Figure 7 . These high B factor regions are more flexible regions which could be affected by protease activity. The results of mutation analyses have shown that the mutation form aspartate 427 to arginine would increase the thermal stability of the ACE2. The Proline 451 to methionine mutation (PoPMuSiC prediction) and the glycine 448 to tryptophan mutation (HotMuSiC prediction) were predicted to increase the thermodynamic stability and thermostability of the ACE2. These mutations were exerted and optimized using Molegro VIRTUAL DOCKER. Figure 7 . The ACE2 structure colored by B factor change. The structure is colored blue to red as a function of B factor increase. The region with highest B factors is circled. Analyzing the RMSF plot of the perfumed simulation indicated that the flexibility of mutated amino acids during the simulation is not high (Figure 8a ). This means that there are no large-scale conformational transitions of protein system at the mutations sites. The superimposed structures of 10 protein models during the simulation has also indicated the low structural fluctuation ( Figure 8b ) The RMSD diagram of the mutated was plotted as a function of time. The diagram shows the deviation of the protein structure at different simulation times compared to the pre-simulation structure. The results showed that the enzyme reached stability after 10 ns. The fluctuations of the RMSD graph should be due to the amino acid changes. However, at the end of the simulation time the system reached stability (Figure 9a ). The R g is another important factor for protein stability analysis. The lower the R g, the more folded and stable is the protein. As shown, the radius of gyration of the mutant protein decreased during the simulation, which means that its stability is increased during the simulation (Figure 9b ). Figure 9 . The RMSD (a) and Rg (b) plots of the mutant ACE2 as a function of time The results of the molecular docking analyses indicated that the mutant ACE2 is capable of interacting with RBD with proper orientation (Figure 10 ). This interaction orientation could include all original interactions across the ACE2 and RBD interface. The low RMSD (0.8 A) between the original and the mutated ACE2 proteins indicates that the exerted mutation does not make a drastic changes in the protein structure. These results indicate that the mutant ACE2 could compete with the cell membrane ACE2 to bind to the RBD of the Spike protein. The immunological analyses indicated that the mutant ACE2 does not contain any newly formed epitopes as the consequence of the exerted mutations. Moreover, none of the glycosylation sites were omitted or added to the glycosylation pattern of the ACE2 protein following the mutations. The 55, 92, 105, 324, 424 and 548 positions were predicted to be the glycosylation sites at both wild type and mutant ACE2 protein. The intact pattern of the glycosylation would help the mutant ACE2 to be considered as a self-protein upon clinical administration. There is still no definitive drug or effective vaccine to fight against COVID-19. Due to its pivotal role in SARS-CoV-2 infection, the interaction between ACE2 and RBD could be targeted for drug or vaccine development (Prompetchara et al., 2020) . Vaccines are considered to be one of the most cost-effective healthcare tools in the prevention of infectious diseases. As expected, various strategies have been adapted to develop effective vaccines against RBD of SARS-CoV-2 (Chen et al., 2020a) . This is while alternative competitive molecules to interrupt this interaction could be as effective (Molina et al., 2020) . ACE2 could be the alternative target for drug or vaccine design efforts against SARS-CoV-2. This enzyme is expressed in variety of tissues with profound physiological roles (Der Sarkissian et al., 2008; Rentzsch et al., 2008; Tipnis et al., 2000; Wong et al., 2007) . Therefore, ACE2 targeted treatments bear the risk to have negative effects on its inherent physiological roles. Given these circumstances, a recombinant ACE2 molecule which lacks its enzymatic activity with higher thermostability and affinity towards RBD would infections (Monteil et al.) . However, their study does not focus on enhancing the binding affinity and thermostability of the ACE2. It should be noted that recombinant ACE2 has already been tested in phase 1 and phase 2 clinical trials and its administration is reported to be well tolerated by healthy human subjects (Haschke et al., 2013; Khan et al., 2017) . Two mutations were performed at the active site of the ACE2 enzyme according the results of network analysis, 2D interaction plot and previously reported studies. One of the mutations at the active site of the ACE2 protein has already been shown to abolish the enzymatic activity of this protein via site directed mutagenesis (Guy et al., 2005) . This study has corroborated the accuracy of mutation candidates introduced by the employed approach. Although the saturation mutation results proposed the Arginine 273 to Glycine mutation, we have preferred the experimentally approved Arginine to Glutamine mutation. It seems that this mutation could omit the hydrogen bond between Arginine 273 and Aspartate 1 of the ANG2 molecule, while the active site cavity would be occupied by the bulky side chain of the Glutamine. The Threonine 445 to Glycine mutation does the same regarding the hydrogen bond omission and preservation of the occupied volume by an equivalent sidechain. Replacing the selected residues with amino acids of similar size would limit the structural changes to a minimum. These mutations had been shown to have a significant effect on binding orientation and affinity of the ANG2 against the active site of the ACE2. The cavity size of the ACE2 active site has slightly become larger after the introduction of the mutations. This increase in cavity volume seems expected due to the slightly smaller size of the replaced amino acids. These changes would inhibit the ACE2 enzyme to perform its carboxypeptidase function due to the loss of key enzyme and substrate interactions and the lack of proper lock between enzyme and substrate in a correct orientation. Reduced enzymatic activity of ACE2 would bring about reduced physiological consequences upon its prescription. Aside from reduced enzymatic activity, the engineered ACE2 needs to have higher affinity toward the RBD in comparison to the membrane attached ACE2 of the target cells. This property would give the engineered ACE2 to have a winning chance to interact with RBD and compete with membrane attached ACE2. Since the ACE2 and RBD molecules are complexed in a certain orientation (Chen et al., 2020b; Yan et al., 2020) , the mutations to enhance the binding affinity should preserve this orientation. Our results indicate that the threonine 27 to Arginine mutation does not have a drastic effect on the interaction orientation and the interacting amino acids. However, new amino acids with extra hydrogen bonds and hydrophobic interactions are engaged at the interaction interface. These results could be construed as stronger interactions between ACE2 and RBD with similar interaction orientation following the threonine 27 to Arginine mutation. Recombinant ACE2 has already been used for severe acute lung injury, acute ANG2-induced hypertension, ANG2-induced heart hypertrophy, cardiac dysfunction, adverse myocardial remodeling, renal oxidative stress, inflammation, and fibrosis (Chen et al., 2016; Gu et al., 2016; Huentelman et al., 2005; Imai et al., 2005; Zhong et al., 2011; Zou et al., 2014) . However, according to the pharmacokinetic studies the half-life of recombinant ACE2 is only hours in humans and mice due to its fast clearance rates (Haschke et al., 2013; Wysocki et al., 2010) . Resistance against proteolytic degradation is one of the most important factors which could determine the protein half-life. Prior studies involving Lasparaginase, β-galactosidase, λ-repressor, Clostridium difficile toxin A-specific antibodies, and Bacillus subtilis lipase have shown that increasing a protein's thermostability may also increase its resistance to protease digestion (Ahmad et al., 2012; Daniel et al., 1982b; Hussack et al., 2011; Parsell and Sauer, 1989) . Thermal and structural stability of proteins are important determinants of proteolytic stability (Akasako et al., 1995; Daniel et al., 1982a; Parsell and Sauer, 1989) . Proteins are resistant against proteolytic degradation in their native conformation. However, they are highly susceptible to proteases after unfolding. The accessible and flexible loop regions of proteins could be attacked by unspecific proteases even in the native conformation of the highly compact and thermostable proteins (Hubbard, 1998; Sarath et al., 1989) . A number of strategies have been harnessed to extend the low half-life of therapeutic proteins. Amino acid manipulation, post translational modification, bio-conjunction and carrier mediated delivery are the main categories of technologies reported for half-life extension (Zaman et al., 2019) . In order to increase the thermostability of the ACE2 protein we have used the amino acid manipulation strategy. The first approach involved the mutation of amino acids with highest B factors (aspartate 427 to arginine). The diffusion of atomic electron densities with respect to their equilibrium positions due to thermal motion and positional disorder is represented by B factor. Thus, B factor shows the extent of mobility for individual amino acids. Replacing an amino acid with high B factor would increase the protein rigidity (Reetz et al., 2006; Xie et al., 2014) . On the other hands, two other mutations were performed to increase the thermodynamic stability and thermostability of the ACE2. The mutations were selected at the regions which are far from the RDB binding site to avoid any destructive effects on the binding affinity. The obtained MD results have also corroborated the stability of the mutated ACE2 structure. We believe that increased thermodynamic stability and thermostability of the ACE2 would enhance the protein halflife. It should be noted that ACE2 amino acid manipulation could be accompanied with other half-life extension strategies like bio-conjugation. Protein engineering is a powerful tool to design therapeutic proteins with desired properties. Bioinformatics tools are contemporary widely used for protein engineering. These tools could save a lot of time and effort in the case of diseases like COVID-19 which time is essential. Using these tools an engineered ACE2 protein has been designed which is expected to bind RBD with higher affinity while it is more thermostable and lacks enzymatic activity. Since ACE2 is reported to be the main RBD binding receptor, its application as an inhibitory agent provides an evolutionary advantage over other conventionally employed therapeutics. Moreover, using the extracellular domain of ACE2 reduces the possible adverse effects and immunological responses against this therapeutic molecule. In vitro an in vivo examination of the engineered protein would unveil its true potentials. Probing protein stability and proteolytic resistance by loop scanning: a comprehensive mutational analysis High resistance of Escherichia coli ribonuclease HI variant with quintuple thermostabilizing mutations to thermal denaturation, acid denaturation, and proteolytic degradation Angiotensin-converting enzyme 2-a new cardiac regulator RINspector: a Cytoscape app for centrality analyses and DynaMine flexibility prediction Epidemiological and clinical features of human coronavirus infections among different subsets of patients Angiotensin-converting enzyme 2 ameliorates renal fibrosis by blocking the activation of mTOR/ERK signaling in apolipoprotein E-deficient mice The SARS-CoV-2 Vaccine Pipeline: an Overview Optimization of the production process and characterization of the yeast-expressed SARS-CoV recombinant receptor-binding domain (RBD219-N1), a SARS vaccine candidate Structure analysis of the receptor binding of 2019-nCoV Affinity enhancement of an in vivo matured therapeutic antibody using structure-based computational design Purified coronavirus spike protein nanoparticles induce coronavirus neutralizing antibodies in mice Hosts and sources of endemic human coronaviruses Angiotensin-converting enzyme 2 is an essential regulator of heart function A correlation between protein thermostability and resistance to proteolysis A correlation between protein thermostability and resistance to proteolysis Residues crucial for maintaining short paths in network communication mediate signaling in proteins Cardiac overexpression of angiotensin converting enzyme 2 protects the heart from ischemia-induced pathophysiology A novel angiotensin-converting enzyme-related carboxypeptidase (ACE2) converts angiotensin I to angiotensin 1-9 Biochemical aspects of coronavirus replication and virus-host interaction Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor Angiotensinconverting enzyme 2 inhibits lung injury induced by respiratory syncytial virus Identification of critical active-site residues in angiotensin-converting enzyme-2 (ACE2) by site-directed mutagenesis Pharmacokinetics and pharmacodynamics of recombinant human angiotensin-converting enzyme 2 in healthy human subjects The structural aspects of limited proteolysis of native proteins Protection from angiotensin II-induced cardiac hypertrophy and fibrosis by systemic lentiviral delivery of ACE2 in rats Engineered single-domain antibodies with high protease resistance and thermal stability Angiotensin-converting enzyme 2 protects from severe acute lung failure A pilot clinical trial of recombinant human angiotensin-converting enzyme 2 in acute respiratory distress syndrome Potent neutralization of 2019 novel coronavirus by recombinant ACE2-Ig A collection of designed peptides to target SARS-Cov-2-ACE2 interaction: PepI-Covid19 database Inhibition of SARS-CoV-2 infections in engineered human tissues using clinical-grade soluble human ACE2 Interaction of the spike protein RBD from SARS-CoV-2 with ACE2: similarity with SARS-CoV, hot-spot analysis and effect of the receptor polymorphism Potential rapid diagnostics, vaccine and therapeutics for 2019 novel Coronavirus (2019-ncoV): a systematic review The structural stability of a protein is an important determinant of its proteolytic susceptibility in Escherichia coli Ofatumumab monoclonal antibody affinity maturation through in silico modeling Affinity maturation and characterization of the ofatumumab monoclonal antibody Coronaviruses post-SARS: update on replication and pathogenesis The RING 2.0 web server for high quality residue interaction networks Immune responses in COVID-19 and potential vaccines: Lessons learned from SARS and MERS epidemic Iterative saturation mutagenesis on the basis of B factors as a strategy for increasing protein thermostability Transgenic angiotensin-converting enzyme 2 overexpression in vessels of SHRSP rats reduces blood pressure and improves endothelial function mCSM-PPI2: predicting the effects of mutations on protein-protein interactions Proteolytic enzymes: a practical approach In silico engineering towards enhancement of bap-VHH monoclonal antibody binding affinity Cytoscape: a software environment for integrated models of biomolecular interaction networks A human homolog of angiotensin-converting enzyme cloning and functional expression as a captopril-insensitive carboxypeptidase Genomic characterization of a newly discovered coronavirus associated with acute respiratory distress syndrome in humans The HADDOCK2. 2 web server: user-friendly integrative modeling of biomolecular complexes Loss of angiotensin-converting enzyme-2 (Ace2) accelerates diabetic kidney injury Genome composition and divergence of the novel coronavirus (2019-nCoV) originating in China Targeting the Degradation of Angiotensin II With Recombinant Angiotensin-Converting Enzyme 2: Prevention of Angiotensin II-Dependent Hypertension Enhanced enzyme kinetic stability by increasing rigidity within the active site Improving GPX activity of selenium-containing human single-chain Fv antibody by site-directed mutation based on the structural analysis Structural basis for the recognition of the SARS-CoV-2 by full-length human ACE2 Current strategies in extending half-lives of therapeutic proteins The digestive system is a potential route of 2019-nCov infection: a bioinformatics analysis based on single-cell transcriptomes Single-cell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCov Prevention of angiotensin II-mediated renal oxidative stress, inflammation, and fibrosis by angiotensin-converting enzyme 2 Angiotensinconverting enzyme 2 protects from lethal avian influenza A H5N1 infections