Exploring the impact of CYP11A1's missense SNPs on the interaction between CYP11A1 and cholesterol: A comprehensive structural analysis and MD simulation study

Abstract

The process of steroidogenesis plays a vital role in human physiology as it governs the biosynthesis of mineralocorticoids, glucocorticoids, and androgens. These three classes of steroid hormones are primarily produced in the adrenal and gonadal glands through steroidogenesis pathways. Initiated by the side chain cleavage of cholesterol (CLR), this process leads to the conversion of cholesterol into pregnenolone and isocaproic aldehyde. The enzyme CYP11A1, encoded by the CYP11A1 gene, plays a key role in catalyzing the side chain cleavage of CLR. Several single nucleotide polymorphisms (SNPs) have been identified in the CYP11A1 gene, which may predispose carriers to disorders associated with abnormal steroidogenesis. Specifically, missense SNPs in the CYP11A1 gene have the potential to negatively impact the interaction between CYP11A1 and CLR, thus affecting the overall metabolome of steroid hormones. In this computational study, we focused on a specific set of missense SNPs reported in the CYP11A1 gene, aiming to identify variants that directly impact the interaction between CYP11A1 and CLR. The three-dimensional structure of the CYP11A1-CLR complex was obtained from the RCSB Protein Data Bank, while missense SNPs in the CYP11A1 gene were retrieved from Ensembl. To predict the most deleterious variants, we utilized the ConSurf server, SIFT, and PolyPhen. Furthermore, we assessed the impact of induced amino acid (AA) substitutions on the CYP11A1-CLR interaction using the PRODIGY server, PyMol, and Ligplot programs. Additionally, molecular dynamics (MD) simulations were conducted to analyze the effects of deleterious variants on the structural dynamics of the CYP11A1-CLR complex. Among the 8096 retrieved variants, we identified ten missense SNPs (E91K, W147G, R151W, R151Q, S391C, V392M, Q395K, Q416E, R460W, and R460Q) as deleterious for the interaction between CYP11A1 and CLR. MD simulations of the CYP11A1-CLR complexes carrying these deleterious AA substitutions revealed that Q416E, W147G, R460Q, and R460W had the most pronounced impacts on the structural dynamics of the complex. Consequently, these missense SNPs were considered the most deleterious ones. Further functional tests are recommended to assess the impact of these four missense SNPs on the enzymatic activity of CYP11A1. Moreover, Genome-Wide Association Studies (GWAS) should be conducted to determine the significance of their association with abnormal steroidogenesis diseases in various patient groups.

Keywords: CYP11A1; In silico; Missense SNPs; Molecular dynamics; P450scc; Steroidogenesis.