Exploring the mechanism of a regulatory SNP of KLK3 by molecular dynamics simulation

The SNP -158G>A of KLK3 has been validated as a regulatory SNP (rSNP) by molecular biology assays, but the mechanism of how it affects the binding of an androgen receptor (AR) homodimer with DNA is unclear. In the current study, molecular dynamics simulation was adopted to explain its inner cause. Based on a recent review), three types of intermolecular forces were analyzed, and the differences among them were compared between complexes containing -158 A:T and -158 G:C. Extra hydrophobic contacts caused by the methyl group on the mutated thymine were the most crucial factor to the regulatory effect of this rSNP. Further analysis concerning the relative motion of the two recognition helixes of the AR homodimer indicated that the hydrophobic interactions between the recognition helix B and the major groove containing -158 A:T changed that helix's motion greatly from swaying in a plane at free state to vibrating slightly around an equilibrium position. A relatively full explanation on the occurrence of rSNP -158G>A is presented here.