Oncogenic mutations in the DNA-binding domain of FOXO1 disrupt folding: quantitative insights from experiments and molecular simulations

FOXO1, a member of the family of winged-helix motif Forkhead box (FOX) transcription factors, is the most abundantly expressed FOXO member in mature B-cells. Sequencing of diffuse large B-cell lymphoma (DLBCL) tumors and cell lines identified specific mutations in the forkhead domain linked to loss of function. Differential scanning calorimetry and thermal shift assays were used to characterize how eight of these mutations affect the stability of the FOX domain. Mutations L183P and L183R were found to be particularly destabilizing. Electrophoresis mobility shift assays show these same mutations also disrupt FOXO1 binding to their canonical DNA sequences, suggesting the loss of function is due to destabilization of the folded structure. Computational modeling of the effects of mutations on FOXO1 folding was performed using alchemical FEP, and a Markov model of the entire folding reaction was constructed from massively parallel molecular simulations, which predicts folding pathways involving the late folding of helix 3. Although FEP is able to qualitatively predict the destabilization from L183 mutations, we find that a simple hydrophobic transfer model, combined with estimates of unfolded-state solvent accessible surface areas from molecular simulation, is able to more accurately predict changes in folding free energies due to mutations. These results suggest that atomic detail provided from simulation is important for accurate prediction of mutational effects on folding stability. Corresponding disease-associated mutations in other FOX family members support further experimental and computational studies of the folding mechanism of FOX domains.