Deep mutational scanning reveals a de novo disulfide bond and combinatorial mutations for engineering thermostable myoglobin

Engineering protein stability is a critical challenge in biotechnology. Here, we used massively parallel deep mutational scanning (DMS) to comprehensively explore the mutational stability landscape of human myoglobin (hMb) and identify key mutations that enhance stability. Our DMS approach involved screening over 10,000 hMb variants by yeast surface display, single-cell sorting, and high-throughput DNA sequencing. We show how surface display levels serve as a proxy for thermostability of soluble hMb variants and report strong correlations between DMS-derived display levels and top-performing machine learning stability prediction algorithms. This approach led to the discovery of a variant with a de novo disulfide bond between residues R32C and C111, which increased thermostability by >12°C compared with wild-type hMb. By combining single stabilizing mutations with R32C, we engineered combinatorial variants that exhibited predominantly additive effects on stability with minimal epistasis. The most stable combinatorial variant exhibited a denaturation temperature exceeding 89°C, representing a >17°C improvement over wild-type hMb. Our findings demonstrate the capabilities in DMS-assisted combinatorial protein engineering to guide the discovery of thermostable variants and highlight the potential of massively parallel mutational analysis for the development of proteins for industrial and biomedical applications.