Human triokinase/flavin mononucleotide (FMN) cyclase (hTKFC) catalyzes the adenosine triphosphate (ATP)-dependent phosphorylation of D-glyceraldehyde and dihydroxyacetone (DHA), and the cyclizing splitting of flavin adenine dinucleotide (FAD). hTKFC structural models are dimers of identical subunits, each with two domains, K and L, with an L2-K1-K2-L1 arrangement. Two active sites lie between L2-K1 and K2-L1, where triose binds K and ATP binds L, although the resulting ATP-to-triose distance is too large (≈14 Å) for phosphoryl transfer. A 75-ns trajectory of molecular dynamics shows considerable, but transient, ATP-to-DHA approximations in the L2-K1 site (4.83 Å or 4.16 Å). To confirm the trend towards site closure, and its relationship to kinase activity, apo-hTKFC, hTKFC:2DHA:2ATP and hTKFC:2FAD models were submitted to normal mode analysis. The trajectory of hTKFC:2DHA:2ATP was extended up to 160 ns, and 120-ns trajectories of apo-hTKFC and hTKFC:2FAD were simulated. The three systems were comparatively analyzed for equal lengths (120 ns) following the principles of essential dynamics, and by estimating site closure by distance measurements. The full trajectory of hTKFC:2DHA:2ATP was searched for in-line orientations and short distances of DHA hydroxymethyl oxygens to ATP γ-phosphorus. Full site closure was reached only in hTKFC:2DHA:2ATP, where conformations compatible with an associative phosphoryl transfer occurred in L2-K1 for significant trajectory time fractions.
Keywords: FMN cyclase; active-site closure; dihydroxyacetone kinase; essential dynamics; molecular dynamics simulation; normal mode analysis; phosphoryl transfer mechanism; protein domain mobility; triokinase.