Homology modeling and molecular dynamics simulation studies of human type 1 3β-hydroxysteroid dehydrogenase: Toward the understanding of cofactor specificity

Engineering the cofactor binding site of 7α-hydroxysteroid dehydrogenase for improvement of catalytic activity, thermostability, and alteration of substrate preference

Hydroxysteroid dehydrogenases (HSDHs) are crucial for bile acid metabolism and influence the size of the bile acid pool and gut microbiota composition. HSDHs with high activity, thermostability, and substrate selectivity are the basis for constructing engineered bacteria for disease treatment. In this study, we designed mutations at the cofactor binding site involving Thr15 and Arg16 residues of HSDH St-2-2. The T15A, R16A, and R16Q mutants exhibited 7.85-, 2.50-, and 4.35-fold higher catalytic activity than the wild type, respectively, while also displaying an altered substrate preference (from taurocholic acid (TCA) to taurochenodeoxycholic acid (TCDCA)). These mutants showed lower Km and higher kcat values, indicating stronger binding to the substrate and resulting in 3190-, 3123-, and 3093-fold higher kcat/Km values for TCDCA oxidation. Furthermore, the Tm values of the T15A, R16A, and R16Q mutants were found to increase by 4.3 °C, 6.0 °C, and 7.0 °C, respectively. Molecular structure analysis indicated that reshaped internal hydrogens and surface mutations could improve catalytic activity and thermostability, and altered interactions among the catalytic triad, cofactor binding sites, and substrates could change substrate preference. This work provides valuable insights into modifying substrate preference as well as enhancing the catalytic activity and thermostability of HSDHs by targeting the cofactor binding site.

Structure-based rational design of hydroxysteroid dehydrogenases for improving and diversifying steroid synthesis

A group of steroidogenic enzymes, hydroxysteroid dehydrogenases are involved in steroid metabolism which is very important in the cell: signaling, growth, reproduction, and energy homeostasis. The enzymes show an inherent function in the interconversion of ketosteroids and hydroxysteroids in a position- and stereospecific manner on the steroid nucleus and side-chains. However, the biocatalysis of steroids reaction is a vital and demanding, yet challenging, task to produce the desired enantiopure products with non-natural substrates or non-natural cofactors, and/or in non-physiological conditions. This has driven the use of protein design strategies to improve their inherent biosynthetic efficiency or activate their silent catalytic ability. In this review, the innate features and catalytic characteristics of enzymes based on sequence-structure-function relationships of steroidogenic enzymes are reviewed. Combining structure information and catalytic mechanisms, progress in protein redesign to stimulate potential function, for example, substrate specificity, cofactor dependence, and catalytic stability are discussed.

Determinants of Cofactor Specificity for the Glucose-6-Phosphate Dehydrogenase from Escherichia coli: Simulation, Kinetics and Evolutionary Studies

Glucose 6-Phosphate Dehydrogenases (G6PDHs) from different sources show varying specificities towards NAD+ and NADP+ as cofactors. However, it is not known to what extent structural determinants of cofactor preference are conserved in the G6PDH family. In this work, molecular simulations, kinetic characterization of site-directed mutants and phylogenetic analyses were used to study the structural basis for the strong preference towards NADP+ shown by the G6PDH from Escherichia coli. Molecular Dynamics trajectories of homology models showed a highly favorable binding energy for residues K18 and R50 when interacting with the 2'-phosphate of NADP+, but the same residues formed no observable interactions in the case of NAD+. Alanine mutants of both residues were kinetically characterized and analyzed with respect to the binding energy of the transition state, according to the kcat/KM value determined for each cofactor. Whereas both residues contribute to the binding energy of NADP+, only R50 makes a contribution (about -1 kcal/mol) to NAD+ binding. In the absence of both positive charges the enzyme was unable to discriminate NADP+ from NAD+. Although kinetic data is sparse, the observed distribution of cofactor preferences within the phylogenetic tree is sufficient to rule out the possibility that the known NADP+-specific G6PDHs form a monophyletic group. While the β1-α1 loop shows no strict conservation of K18, (rather, S and T seem to be more frequent), in the case of the β2-α2 loop, different degrees of conservation are observed for R50. Noteworthy is the fact that a K18T mutant is indistinguishable from K18A in terms of cofactor preference. We conclude that the structural determinants for the strict discrimination against NAD+ in the case of the NADP+-specific enzymes have evolved independently through different means during the evolution of the G6PDH family. We further suggest that other regions in the cofactor binding pocket, besides the β1-α1 and β2-α2 loops, play a role in determining cofactor preference.