
Int. J. Mol. Sci. 2023,24, 13257 2 of 17
sequences with a specific function, like antimicrobial peptide sequences [
20
,
21
], they re-
quire a significant amount of experimental data. Acquiring such extensive data on peptide
sequences interacting with specific PPI-targeting proteins is challenging. Computational
methods exist to generate peptides similar to known PPI inhibitors [
22
], but there are few
small molecules that can inhibit PPIs [7,23].
High-throughput virtual screening (HTVS) technology is a potent and efficient method
for pinpointing drug candidates from a vast compound library [
24
,
25
]. Thus, conducting
HTVS of peptides with the right virtual screening steps and strategies could feasibly yield
bioactive peptide hit candidates [
26
]. The primary technology underpinning HTVS is
molecular docking tools, which encompass both commercial and open-source tools. There
are commercial options available, such as the Molecular Operating Environment suite,
the Schrödinger modeling suite, and the OpenEye Scientific Software suite. Tools like
Schrödinger Glide can perform protein–peptide docking, though they’re essentially crafted
for linear peptides [
27
,
28
]. When these tools are employed for cyclic peptides, the peptide
is treated as a ligand, and protein–ligand docking is undertaken. Glide has introduced a
macrocycle module that accommodates both flexible and rigid docking [
29
]. Nevertheless,
these tools are optimized primarily for conventional small molecule compounds. A thor-
ough understanding of the parameters for conformation generation and scoring functions
is indispensable when docking larger cyclic peptides as ligands.
Conversely, while numerous open-source tools exist for global protein–peptide dock-
ing, such as CABS-dock [
30
,
31
], ClusPro PeptiDock [
32
], and PIPER-FlexPepDock [
33
], and
for local protein–peptide docking, including PEP-FOLD3 [
34
] and DINC 2.0 [
35
], they can-
not handle cyclic peptides. Only two tools, HADDOCK 2.4 [
36
] and AutoDock CrankPep
(ADCP) [
37
,
38
], are specifically tailored for protein–cyclic peptide docking, with ADCP
being recognized as the state-of-the-art. Despite significant advancements in peptide dock-
ing tools, which are crucial for peptide HTVS strategies, the pace of progress in peptide
HTVS remains considerably slower than that for small molecule compounds. Present HTVS
methodologies, even those utilizing cutting-edge techniques, appear insufficient in fully
tapping the potential of peptide libraries. Notably, few reports in the literature highlight
effective use of HTVS, with only a handful of studies showcasing success through both
proprietary and open computational techniques [39].
Given this context, we introduce a novel method in this study to design cyclic pep-
tides targeting PPI. This method harnesses a deep learning (DL) based protein structure
prediction, diverging from the traditional HTVS strategies reliant on molecular docking.
DL techniques, notably AlphaFold2 (AF2) and RoseTTAFold, have made significant strides
in the accurate prediction of 3D protein structures from amino acid sequences [
40
–
42
].
These tools have demonstrated capabilities in predicting individual proteins as well as
more complex assemblies like protein complexes and protein–peptide complexes [
43
–
46
].
AF2 can distinguish between favorable and unfavorable conformational templates. This
suggests that AF2 learns the approximate biophysical energy function of energetically
stable protein backbones, and that DL-based structure prediction methods can link the
sequence space and protein structure space to explore protein structures that satisfy a given
design [
47
]. Deep network hallucinations, a DL-based design method, can address a variety
of protein designs by starting with random and iterative updating of the protein sequence
until the desired properties specified in the loss function are obtained [
48
]. To increase
the possibility of generating more accurate predictions, the loss function converges more
quickly using sequence gradients obtained by inverting the structure prediction network
and using continuous logits, rather than discrete one-hot encodings of the sequence repre-
sentation [
49
–
51
]. In our previous study, we designed linear peptide binders in AfDesign
using a three-step sequence representation of logits, softmax, and one-hot steps [
52
]. AfDe-
sign binder hallucination is based on AF2, so in principle it is not possible to design cyclic
peptides. However, if AF2 can predict the structure of protein–cyclic peptide complexes, it
can design sequences that are likely to form complexes with the target protein, which is
quite efficient methods.