Peptide probes derived from pertuzumab by molecular dynamics modeling for HER2 positive tumor imaging.
PLoS Comput Biol 2017;
13:e1005441. [PMID:
28406988 PMCID:
PMC5390981 DOI:
10.1371/journal.pcbi.1005441]
[Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2016] [Accepted: 03/04/2017] [Indexed: 01/19/2023] Open
Abstract
A high level of HER2 expression in breast cancer correlates with a higher tumor growth rate, high metastatic potential, and a poor long-term patient survival rate. Pertuzumab, a human monoclonal antibody, can reduce the effect of HER2 overexpression by preventing HER2 dimerization. In this study, a combination protocol of molecular dynamics modeling and MM/GBSA binding free energy calculations was applied to design peptides that interact with HER2 based on the HER2/pertuzumab crystal structure. Based on a β hairpin in pertuzumab from Glu46 to Lys65—which plays a key role in interacting with HER2—mutations were carried out in silico to improve the binding free energy of the hairpin that interacts with the Phe256-Lys314 of the HER2 protein. Combined the use of one-bead-one-compound library screening, among all the mutations, a peptide (58F63Y) with the lowest binding free energy was confirmed experimentally to have the highest affinity, and it may be used as a new probe in diagnosing and treating HER2-positive breast cancer.
Many therapeutic approaches, including the human monoclonal antibodies trastuzumab and pertuzumab, target the human epidermal growth factor receptor 2 (HER2) of any breast cancer that features HER2 overexpression. Compared to these antibodies, peptides have many advantages, including lower cost, easier synthesis, high affinity, and lower toxicity. Here, we first designed peptides that interact with HER2 protein based on the HER2/pertuzumab crystal structure (PDB entry: 1S78), using a combination protocol of molecular dynamics modeling, molecular mechanics/generalized Born solvent-accessible surface area (MM/GBSA) binding free energy calculations. Then, combined with the peptide library screening, six peptides were selected for further analysis and experimental validations. The results of ex vivo and in vivo experiments confirmed that one peptide (58F63Y) in particular has a strong affinity and high specificity to HER2-overexpressing tumors. This may due to more paired residues and lower binding free energy in peptide 58F63Y and HER2 complex based on free energy decomposition analysis and distances calculation. While both in silico and in vitro screenings point to the same high-affinity peptide, the findings suggest that in silico screening based on calculated binding free energies is rather reliable. Additionally, based on the calculation of binding free energies among mutants, we can reduce the library capacity of one-bead-one-compound screening. In summary, we present a rather simple and rapid means of deriving a peptide with a clear binding site to its target protein.
Collapse