Machine Learning-Based Development of Gadolinium Binding Peptides.
Authors
Affiliations (7)
Affiliations (7)
- Department of Chemical Engineering & Material Science, Michigan State University, East Lansing, Michigan, USA.
- Michigan State University College of Human Medicine, Grand Rapids, Michigan, USA.
- Department of Computer Science and Engineering, Michigan state university, East Lansing, Michigan, USA.
- Department of Chemistry, Michigan State University, East Lansing, Michigan, USA.
- Department of Radiology, Michigan State University, East Lansing, Michigan, USA.
- Faculty of Life Sciences, Tel Aviv University, Tel-Aviv, Israel.
- School of Biomedical Engineering, Tel Aviv University, Tel Aviv, Israel.
Abstract
Gadolinium-based contrast agents (GBCAs) are indispensable tools in magnetic resonance imaging (MRI), yet their clinical use is limited by non-specific tissue accumulation, low molecular specificity, and safety concerns. Protein and peptide scaffolds provide a promising alternative because they can bind metal ions with high selectivity and enable precise molecular targeting. However, identifying short peptide motifs with optimal gadolinium (Gd<sup>3+</sup>) coordination and high relaxivity remains a major challenge. Here, we used a machine-learning-driven peptide evolution platform, the Protein Optimization Engineering Tool (POET), to design and optimize short Gd-binding motifs that enhance longitudinal relaxivity (r<sub>1</sub>). Two algorithmic strategies were tested: motif-based and regular-expression-based representations. Both algorithms were trained on an initial set of 74 twelve-amino-acid peptides derived from natural EF-hand scaffolds. Through two rounds of directed evolution and experimental screening, POET predicted peptides with up to a 24% increase in r<sub>1</sub> ratio compared with the best natural EF-hand. Further analysis revealed that peptides exhibiting higher relaxivity generally possessed a more negative net charge and lower isoelectric point than the buffer pH, indicating stronger electrostatic stabilization of Gd<sup>3+</sup>. Sequence enrichment analysis showed that acidic and small polar residues, particularly aspartic acid, glycine, and threonine, were selectively favored during evolution, while bulky hydrophobic and basic residues were depleted. These compositional trends align with improved solubility and enhanced metal coordination. Together, these results demonstrate a generalizable framework that integrates computational evolution with biophysical screening to discover new biologically derived Gd-binding motifs. This approach provides a scalable route to engineer responsive, tunable, and biocompatible MRI contrast tags for precision imaging and molecular diagnostics.