Exploration of the Substrate Preference of Lysine Methyltransferase SMYD3 by Molecular Dynamics Simulations

Abstract

SMYD3, a SET and MYND domain containing lysine methyltransferase, catalyzes the transfer of the methyl group from a methyl donor onto the Nϵ group of a lysine residue in the substrate protein. Methylation of MAP3 kinase kinase (MAP3K2) by SMYD3 has been implicated in Ras-driven tumorigenesis. The crystal structure of SMYD3 in complex with MAP3K2 peptide reveals a shallow hydrophobic pocket (P-2), which accommodates the binding of a phenylalanine residue at the -2 position of the substrate (F258) is a crucial determinant of substrate specificity of SMYD3. To better understand the substrate preference of SMYD3 at the -2 position, molecular dynamics (MD) simulations and the MM/GBSA method were performed on the crystal structure of SMYD3-MAP3K2 complex (PDB: 5EX0) after substitution of F258 residue of MAP3K2 to each of the other 19 natural residues, respectively. Binding free energy calculations reveal that the P-2 pocket prefers an aromatic hydrophobic group and none of the substitutions behave better than the wild-type phenylalanine residue does. Furthermore, we investigated the structure-activity relationships (SAR) of a series of non-natural phenylalanine derivative substitutions at the -2 position and found that quite a few modifications on the sidechain of F258 residue could strengthen its binding to the P-2 pocket of SMYD3. These explorations provide insights into developing novel SMYD3 inhibitors with high potency and high selectivity against MAP3K2 and cancer.