Abstract
The field of therapeutic peptides is experiencing a surge, fueled by their advantageous features. These include predictable metabolism, enhanced safety profile, high selectivity, and reduced off-target effects compared with small-molecule drugs. Despite progress in addressing limitations associated with peptide drugs, a significant bottleneck remains: the absence of a large-scale in silico screening method for a given protein target structure. Such methods have proven invaluable in accelerating small-molecule drug discovery. The high flexibility of peptide structures and the large diversity of peptide sequences greatly hinder the development of urgently needed computational methods. Here, we report a method called MDockPeP2_VS to address these challenges. It integrates molecular docking with structural conservation between protein folding and protein-peptide binding. Briefly, we discovered that when the interfacial residues are conserved, a sequence fragment derived from a monomeric protein exhibits a high propensity to bind a target protein with a similar conformation. This valuable insight significantly reduces the search space for peptide conformations, resulting in a substantial reduction in computational time and making in silico peptide screening practical. We applied MDockPeP2_VS to develop peptide inhibitors targeting the TEM-1 β-lactamase of Escherichia coli, a key mechanism behind antibiotic resistance in gram-negative bacteria. Among the top 10 peptides selected from in silico screening, TF7 (KTYLAQAAATG) showed significant inhibition of β-lactamase activity with a K (i) value of 1.37 ± 0.37 µM. This fully automated, large-scale structure-based in silico peptide screening software is available for free download at https://zougrouptoolkit.missouri.edu/mdockpep2_vs/download.html.