# Safe model based optimization balancing exploration and reliability for protein sequence design

**Authors:** Shuuki Takizawa, Keita Mori, Naoto Tanishiki, Dai Yoshimura, Atsushi Ohta, Reiji Teramoto

PMC · DOI: 10.1038/s41598-025-12568-5 · Scientific Reports · 2025-07-29

## TL;DR

This paper introduces a new method for protein design that improves reliability and avoids unreliable predictions.

## Contribution

The novel MD-TPE method penalizes out-of-distribution samples to improve safe optimization in protein sequence design.

## Key findings

- MD-TPE reduces pathological samples compared to TPE in the GFP dataset.
- MD-TPE successfully identifies high-affinity antibody mutants in affinity maturation tasks.

## Abstract

Discovering proteins with desired functionalities using protein engineering is time-consuming. Offline Model-Based Optimization (MBO) accelerates protein sequence design by exploring the vast protein sequence space using a trained proxy model. However, the proxy model often yields excessively good values that are far from the training dataset and causes pathological behavior in the MBO. To address this problem, we propose a mean deviation tree-structured Parzen estimator (MD-TPE) that penalizes unreliable samples located in the out-of-distribution region using the deviation of the predictive distribution of the Gaussian process (GP) model in the objective function to find the solution in the vicinity of the training data, where the proxy model can reliably predict. Upon examining the GFP dataset, compared to TPE, MD-TPE yielded fewer pathological samples. Additionally, it successfully identified mutants with higher binding affinity in the antibody affinity maturation task. Thus, our developed safe optimization approach is useful for protein engineering.

## Linked entities

- **Proteins:** NAL1 (Protein NARROW LEAF 1)

## Full-text entities

- **Genes:** MARVELD3 (MARVEL domain containing 3) [NCBI Gene 91862] {aka MARVD3, MRVLDC3}, CD81 (CD81 molecule) [NCBI Gene 975] {aka CVID6, S5.7, TAPA1, TSPAN28}, RNF130 (ring finger protein 130) [NCBI Gene 55819] {aka G1RP, G1RZFP, GOLIATH, GP}
- **Diseases:** TPE (MESH:D020914), MD (MESH:D010262)
- **Chemicals:** MD (-), PBS (MESH:D007854), Amino acids (MESH:D000596), acids (MESH:D000143)
- **Species:** Homo sapiens (human, species) [taxon 9606]
- **Cell lines:** ESM-1 — Homo sapiens (Human), Transformed cell line (CVCL_XI05), Expi293F — Homo sapiens (Human), Transformed cell line (CVCL_D615)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12307608/full.md

## Figures

4 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12307608/full.md

## References

17 references — full list in the complete paper: https://tomesphere.com/paper/PMC12307608/full.md

---
Source: https://tomesphere.com/paper/PMC12307608