Mechanisms of SSBP1 variants in mitochondrial disease: Molecular dynamics simulations reveal stable tetramers with altered DNA binding surfaces.

Several mutations in the gene for the mitochondrial single stranded DNA binding protein (SSBP1) have recently been implicated in human disease, but initial reports are insufficient to explain the molecular mechanism of disease, including the possible role of SSBP1 heterotetramers in heterozygous patients. Here we employed molecular simulations to model the dynamics of wild type and 31 variant SSBP1 tetramer systems, including 7 variant homotetramer and 24 representative heterotetramer systems. Our simulations indicate that all variants are stable and most have stronger intermonomer interactions, reduced solvent accessible surface areas, and a net loss of positive surface charge. We then used structural alignments and phosphate binding simulations to predict DNA binding surfaces on SSBP1. Our models suggest that nearly the entire surface of SSBP1, excluding flexible loops and protruding helices, is available for DNA binding, and we observed several potential DNA binding hotspots. Changes to the protein surface in variant SSBP1 tetramers potentially alter anchor points or wrapping paths, rather than abolishing binding altogether. Overall, our findings disqualify tetramer destabilization or gross disruption of DNA binding as mechanisms of disease. Instead, they are consistent with subtle changes to DNA binding, wrapping, or release that cause rare but consequential failures of mtDNA maintenance, which, in turn, are consistent with the late onset of disease in most of the reported SSBP1 cases.