Modeling Substrate Coordination to Zn-Bound Angiotensin Converting Enzyme 2

The spike protein in the envelope of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) interacts with the receptor Angiotensin Converting Enzyme 2 (ACE2) on the host cell to facilitate the viral uptake. Angiotensin II (Ang II) peptide, which has a naturally high affinity for ACE2, may be useful in inhibiting this interaction. In this study, we computationally designed several Ang II mutants to find a strong binding sequence to ACE2 receptor and examined the role of ligand substitution in the docking of native as well as mutant Ang II to the ACE2 receptor. The peptide in the ACE2-peptide complex was coordinated to zinc in the ACE2 cleft. Exploratory molecular dynamics (MD) simulations were used to measure the time-based stability of the native and mutant peptides and their receptor complexes. The MD-generated root-mean-square deviation (RMSD) values are mostly similar between the native and seven mutant peptides considered in this work, although the values for free peptides demonstrated higher variation, and often were higher in amplitude than peptides associated with the ACE2 complex. An observed lack of a strong secondary structure in the short peptides is attributed to the latter’s greater flexibility and movement. The strongest binding energies within the ACE2-peptide complexes were observed in the native Ang II and only one of its mutant variants, suggesting ACE2 cleft is designed to provide optimal binding to the native sequence. An examination of the S1 binding site on ACE2 suggests that complex formation alone with these peptides may not be sufficient to allosterically inhibit the binding of SARS-CoV-2 spike proteins. However, it opens up the potential for utilizing AngII-ACE2 binding in the future design of molecular and supramolecular structures to prevent spike protein interaction with the receptor through creation of steric hindrance.