Structural screening and molecular simulation identify potential ligands against the K700E hot spot variant and functional pockets of SF3B1 to modulate splicing in myelodysplastic syndrome

Abstract

Myelodysplastic syndrome (MDS), a blood disorder with ineffective hematopoiesis and risk of transformation to acute myeloid leukemia, is characterized by recurring cytogenetic and molecular alterations. By chromosome analysis, approximately 60% of patients, carry chromosome 5 and 7 alterations, trisomy of chromosome 8 and may also present with increasingly complex karyotypes, especially in higher grade MDS (MDS with refractory anemia and increased blasts type 1 and 2). Moreover, somatic pathogenic variants in genes associated with aberrant mRNA splicing are frequently mutated with SF3B1 the most frequently mutated. In the setting of SF3B1, the K700E hot-spot mutation is present in approximately 50% of cases. Since recent studies have highlighted modulation of functional dynamics in SF3B1 by mutant splicing factors, the objective of the study was to identify potential small molecule modulators against the frequently mutated RNA splicing factor SF3B1(K700E) and functional allosteric sites by using a molecular structure-based approach and a molecular dynamic simulation. To identify potential SF3B1 modulators, we collected a series of chemical compounds from the Zinc and Enamine database. An initial screen followed by further molecular analysis and simulation using the Schrödinger suite was performed. Parameters used to monitor the stability and binding of the protein-ligand complex included: RMSF, protein-ligand contacts, electrostatic, Van Der Waals forces and binding energies (MMGBSA). A 100-nanosecond simulation showed strong binding between selected compounds and key amino acid residues, including the mutation hot-spot K700E and functional allosteric amino acid residue R630. Ligand binding energies between compounds and key amino acid residues ranged from −50.67 to −58.04 kcal/mol. In brief, small molecule modulators show strong binding to SF3B1 suggesting these compounds may be used against cells harboring the K700E variant or to modulate splicing by targeting functional allosteric sites.