Predicting substrate selectivity between UGT1A9 and UGT1A10 using molecular modelling and molecular dynamics approach

Uridine 5′-diphospho-glucuronosyltransferase-1A9 (UGT1A9) expressed in the liver, shows good sequence identity with UGT1A10, expressed in the intestine. Both uridine 5′-diphospho-glucuronosyltransferase (UGT) isoforms show comprehensive overlapping substrate selectivity but there are differences in stereoselectivity, regiospecificity and rate of glucuronidation of the substrates. Multiple sequence alignment analyses of UGT1A9 and UGT1A10 showed that 13% of the residues in N-terminal domain (NTD) are non-identical between them. Herein, authors attempted homology modelling of UGT1A9 and UGT1A10 and validation using software tools and reported mutagenic studies. A molecular docking study of the known substrates is performed on UGT1A9 and UGT1A10 homology models. The non-identical N-terminal residues ranging from 111 to 117 in UGT1A9 and UGT1A10 were identified to play a central role in the substrate selectivity. However, substrate binding is performed by Ser111, Gly115 and Leu117 in UGT1A10 and Gly111, Asp115 and Phe117 in UGT1A9. This study reports new residues in NTD, showing interaction with uridine 5′-diphospho-glucuronic acid which binds with C-terminal domain. Further, molecular dynamics simulations were carried out to study the role of non-identical residues in substrate identification. The study demonstrates the folding of the UGT enzyme, particularly, helix-loop-helix transition and movement of Nα3-2 helix, in response to substrate and co-substrate binding.