Bifunctional Homodimeric Triokinase/FMN Cyclase: CONTRIBUTION OF PROTEIN DOMAINS TO THE ACTIVITIES OF THE HUMAN ENZYME AND MOLECULAR DYNAMICS SIMULATION OF DOMAIN MOVEMENTS*

Mammalian triokinase, which phosphorylates exogenous dihydroxyacetone and fructose-derived glyceraldehyde, is neither molecularly identified nor firmly associated to an encoding gene. Human FMN cyclase, which splits FAD and other ribonucleoside diphosphate-X compounds to ribonucleoside monophosphate and cyclic X-phosphodiester, is identical to a DAK-encoded dihydroxyacetone kinase. This bifunctional protein was identified as triokinase. It was modeled as a homodimer of two-domain (K and L) subunits. Active centers lie between K1 and L2 or K2 and L1: dihydroxyacetone binds K and ATP binds L in different subunits too distant (≈14 Å) for phosphoryl transfer. FAD docked to the ATP site with ribityl 4′-OH in a possible near-attack conformation for cyclase activity. Reciprocal inhibition between kinase and cyclase reactants confirmed substrate site locations. The differential roles of protein domains were supported by their individual expression: K was inactive, and L displayed cyclase but not kinase activity. The importance of domain mobility for the kinase activity of dimeric triokinase was highlighted by molecular dynamics simulations: ATP approached dihydroxyacetone at distances below 5 Å in near-attack conformation. Based upon structure, docking, and molecular dynamics simulations, relevant residues were mutated to alanine, and k_cat and K_m were assayed whenever kinase and/or cyclase activity was conserved. The results supported the roles of Thr¹¹² (hydrogen bonding of ATP adenine to K in the closed active center), His²²¹ (covalent anchoring of dihydroxyacetone to K), Asp⁴⁰¹ and Asp⁴⁰³ (metal coordination to L), and Asp⁵⁵⁶ (hydrogen bonding of ATP or FAD ribose to L domain). Interestingly, the His²²¹ point mutant acted specifically as a cyclase without kinase activity.