Abstract
Inherited retinal dystrophies (IRDs) are clinically and genetically heterogeneous disorders. Most IRDs follow a monogenic inheritance pattern. However, an increasing number of unresolved cases suggest the possible contribution of oligogenic or digenic mechanisms. Here, we report two ultra-rare missense variants-AIPL1 R302L and BBS2 P134R-that co-segregate with early-onset nonsyndromic retinal degeneration in affected individuals from a non-consanguineous family. We performed a multi-level computational investigation to assess whether these variants may act through a convergent pathogenic mechanism. Using AlphaFold2-predicted structures, we modeled both wild-type and mutant proteins, introduced point mutations, and performed energy minimization and validation. FoldX, DynaMut2, and DUET all predicted destabilizing effects at the variant sites, corroborated by local disruption of secondary structure and altered surface electrostatics. Comparative docking (via HDOCK and ClusPro) identified a putative interaction interface between the TPR domain of AIPL1 and the β-sheet face of BBS2. This interface was destabilized in the double-mutant model. At the systems level, transcriptomic profiling confirmed co-expression of AIPL1 and BBS2 in human retina and fetal eye, while functional enrichment analysis highlighted overlapping involvement in ciliary and proteostasis pathways. Network propagation suggested that the two proteins may converge on shared interactors relevant to photoreceptor maintenance. Collectively, these in silico results provide structural and systems-level support for a candidate digenic mechanism involving AIPL1 and BBS2. While experimental validation remains necessary, our study proposes a testable mechanistic hypothesis and underscores the value of computational approaches in uncovering complex genetic contributions to IRDs.