Angiotensin-converting enzyme 2 (ACE2)
is the host cellular receptor
that locks onto the surface spike protein of the 2002 SARS coronavirus
(SARS-CoV-1) and of the novel, highly transmissible and deadly 2019
SARS-CoV-2, responsible for the COVID-19 pandemic. One strategy to
avoid the virus infection is to design peptides by extracting the
human ACE2 peptidase domain α1-helix, which would
bind to the coronavirus surface protein, preventing the virus entry
into the host cells. The natural α1-helix peptide
has a stronger affinity to SARS-CoV-2 than to SARS-CoV-1. Another
peptide was designed by joining α1 with the second
portion of ACE2 that is far in the peptidase sequence yet grafted
in the spike protein interface with ACE2. Previous studies have shown
that, among several α1-based peptides, the hybrid
peptidic scaffold is the one with the highest/strongest affinity for
SARS-CoV-1, which is comparable to the full-length ACE2 affinity.
In this work, binding and folding dynamics of the natural and designed
ACE2-based peptides were simulated by the well-known coarse-grained
structure-based model, with the computed thermodynamic quantities
correlating with the experimental binding affinity data. Furthermore,
theoretical kinetic analysis of native contact formation revealed
the distinction between these processes in the presence of the different
binding partners SARS-CoV-1 and SARS-CoV-2 spike domains. Additionally,
our results indicate the existence of a two-state folding mechanism
for the designed peptide en route to bind to the spike proteins, in
contrast to a downhill mechanism for the natural α1-helix peptides. The presented low-cost simulation protocol demonstrated
its efficiency in evaluating binding affinities and identifying the
mechanisms involved in the neutralization of spike-ACE2 interaction
by designed peptides. Finally, the protocol can be used as a computer-based
screening of more potent designed peptides by experimentalists searching
for new therapeutics against COVID-19.