Comprehensive computational analysis of the molecular mechanism of self-incompatibility in Brassicaceae using improved structure prediction

Abstract

ABSTRACTPlants employ self-incompatibility (SI) to promote cross-fertilization. In Brassicaceae, this process is regulated by the formation of a complex between the pistil determinantSreceptor kinase (SRK) and the pollen determinantS-locus protein 11 (SP11, also known asS-locus cysteine-rich protein, SCR). In our previous study, we used the crystal structures of two eSRK–SP11 complexes inBrassica rapa S8andS9haplotypes and nine computationally predicted complex models to demonstrate that only the SRK ectodomain (eSRK) and SP11 pairs derived from the sameShaplotype exhibit high binding free energy. However, predicting the eSRK–SP11 complex structures for the other 100+Shaplotypes and genera remains difficult because of SP11 polymorphism in sequence and structure. Although protein structure prediction using AlphaFold2 exhibits considerably high accuracy for most protein monomers and complexes, 46% of the predicted SP11 structures that we tested showed < 75 mean per-residue confidence score (pLDDT). Here, we demonstrate that the use of curated multiple sequence alignment (MSA) for cysteine-rich proteins significantly improved model accuracy for SP11 and eSRK–SP11 complexes. Additionally, we calculated the binding free energies of the predicted eSRK–SP11 complexes using molecular dynamics (MD) simulations and observed that someArabidopsishaplotypes formed a binding mode that was critically different from that ofB. rapa S8andS9. Thus, our computational results provide insights into the haplotype-specific eSRK–SP11 binding modes in Brassicaceae at the residue level. The predicted models are freely available at Zenodo,https://doi.org/10.5281/zenodo.8047768.