- Research article
- Open Access
Molecular docking simulation studies on potent butyrylcholinesterase inhibitors obtained from microbial transformation of dihydrotestosterone
© Zafar et al.; licensee Chemistry Central Ltd. 2013
- Received: 6 May 2013
- Accepted: 29 August 2013
- Published: 8 October 2013
Biotransformation is an effective technique for the synthesis of libraries of bioactive compounds. Current study on microbial transformation of dihydrotestosterone (DHT) (1) was carried out to produce various functionalized metabolites.
Microbial transformation of DHT (1) by using two fungal cultures resulted in potent butyrylcholinesterase (BChE) inhibitors. Biotransformation with Macrophomina phaseolina led to the formation of two known products, 5α-androstan-3β,17β-diol (2), and 5β-androstan-3α,17β-diol (3), while biotransformation with Gibberella fujikuroi yielded six known metabolites, 11α,17β-dihydroxyandrost-4-en-3-one (4), androst-1,4-dien-3,17-dione (5), 11α-hydroxyandrost-4-en-3,17-dione (6), 11α-hydroxyandrost-1,4-dien-3,17-dione (7), 12β-hydroxyandrost-1,4-dien-3,17-dione (8), and 16α-hydroxyandrost-1,4-dien-3,17-dione (9). Metabolites 2 and 3 were found to be inactive, while metabolite 4 only weakly inhibited the enzyme. Metabolites 5–7 were identified as significant inhibitors of BChE. Furthermore, predicted results from docking simulation studies were in complete agreement with experimental data. Theoretical results were found to be helpful in explaining the possible mode of action of these newly discovered potent BChE inhibitors. Compounds 8 and 9 were not evaluated for enzyme inhibition activity both in vitro and in silico, due to lack of sufficient quantities.
Biotransformation of DHT (1) with two fungal cultures produced eight known metabolites. Metabolites 5–7 effectively inhibited the BChE activity. Cholinesterase inhibition is among the key strategies in the management of Alzheimer’s disease (AD). The experimental findings were further validated by in silico inhibition studies and possible modes of action were deduced.
- Microbial transformation
- Dihydrotestosterone (DHT)
- Macrophomina phaseolina
- Gibberella fujikuroi
- Butyrylcholinesterase (BChE) inhibition
- Alzheimer’s disease
- Molecular docking simulation
Microbial transformation of steroids is being studied for decades but the need to develop new structural analogues remain strong due to multiple reasons, including quest for medicinally important novel steroids . In continuation of our recent work on the microbial transformation of important steroids , we investigated the microbial biotransformation of DHT. The rationale was to produce the DHT analogues for studying the structure-activity-relationship (SAR) and to synthesize medicinally important compounds with novel activities.
Butyrylcholine (BCh) is an acetylcholine-like compound (Figure 1), hydrolyzed by acetylcholinesterase (AChE) and BChE (pseudocholinesterase), the latter being more efficient. BCh is a synthetic compound, used as a tool to distinguish between acetyl- and BChEs. BChEs is essential for the catalysis of the rapid breakdown of suxamethonium (succinylcholine), a muscle relaxant, frequently used in surgery and electroshock therapy. Prolonged muscle relaxation may follow the administration of drug if pseudocholinesterase activity is defective or markedly decreased [6, 7]. BChE inactivates the neurotransmitter, acetylcholine (ACh). ACh is an important therapeutic target for the treatment of Alzheimer’s disease, characterized by a cholinergic deficit .
DHT was subjected to microbial transformation with two fungal cultures, Macrophomina phaseolina and Gibberella fujikuroi and eight known metabolites 2–9 were obtained. Substrate 1 and its metabolites 2–7 were subjected to AChE and BChE inhibitory activity evaluation. All compounds were found to be inactive against AChE, while metabolite 5–7 have significantly inhibited the BChE. Compounds 8 and 9 were not subjected to the above mentioned activity due to lack of sufficient quantities.
Protein-ligand docking programs are used for the placement of small molecules within the binding pocket of target proteins (receptors) and to rank them according to their binding affinity [9, 10]. In current study, biotransformed products 2–7 were also evaluated in silico to understand their mode of interaction with the BChE. Resolved crystal structure of BChE was used in molecular docking simulation studies. All biotransformed metabolites were docked within the binding pocket of the crystal structure of human BChE (PDB ID 1P0P: 2.30 Å), revealing structural features, responsible of observed enzyme inhibitory activities . MOE docking software was utilized to perform the molecular docking experiment. The outcome of the docking study helped to understand the binding mechanism of compounds with BChE.
Metabolite 2 (C19H32O2) (M + m/z 292.2434, calcd 292.2402) showed no florescence under UV light. The IR spectrum exhibited an absorption at 3350 cm-1 (OH), but no absorption for the ketone group was observed. This suggested that the 2 atomic mass units (a.m.u.) increase in the molecular weight might be due to the reduction of the ketone group of substrate 1 to a hydroxyl group in 2.
The 1H-NMR spectrum of 2 showed two hydroxyl-bearing methine signals at δ 3.49 (m, H-3), and 3.54 (t, J17a,16a,e = 8.4 Hz, H-17). The second hydroxyl signal was resulted from the reduction of the ketone at C-3. The 13C-NMR spectrum had two hydroxyl-bearing carbon signals at δ 71.8 (C-3) and 82.5 (C-17). H-3 (δ 3.49) showed HMBC correlation with C-5 (δ 46.3). H-5 (δ 1.10, m) also showed HMBC cross peaks with C-3 (δ 71.8) and C-10 (δ 36.7). The NOESY interactions between H-3 (δ 3.49) and H-5 (δ 1.10), as well as between H-5 and H-9 (δ 0.64, ddd, J9a,8a = 14.8 Hz, J9a,11a = 10.8 Hz, J9a,11e = 4.4 Hz), indicated the α-orientation (axial) of H-3. H-17 (δ 3.54) showed NOESY interactions with H-14 (δ 0.92), suggesting its α-orientation (Additional file 1). The metabolite 2 was thus characterized as 5α-androstan-3β,17β-diol. The compound 2 has been earlier derived from DHEA (dehydroepiandrosterone) by the action of 17β-hydroxysteroid dehydrogenase type 7 enzyme . 5α-Androstan-3β,17β-diol (2) has been reported to inhibit lipopolysaccharide induced inflammatory response and tumor necrosis factor in human endothelial cells .
The metabolite 3, C19H32O2 [M + = m/z 292.2434 (calcd 292.2402)], was UV inactive, thus lacked α,β-unsaturated carbonyl system. The IR spectrum indicated the presence of -OH (3349 cm-1), but no ketonic absorption. This could be due to the reduction of the carbonyl group at C-3.
The 1H-NMR spectral analysis indicated two hydroxyl-bearing methine-proton triplets at δ 3.94 (J 3a,2/4 = 2.4 Hz, H-3), and 3.55 (J17a,16a,e = 8.8 Hz, H-17), with their respective carbons resonated at δ 67.2 and 82.6, respectively, in the 13C-NMR spectrum (HMQC). HMBC of H-5 (δ 1.48) with carbon at δ 67.2 indicated the reduction of C-3 carbonyl into an OH. The H-3 (δ 3.94) showed NOESY correlation with H-5 (δ 1.48). Interestingly H-5 did not show any NOESY cross peak with H-9 (δ 0.74), indicating that H-5 is not α-oriented. Comparison of the 1H- and 13C-NMR chemical shift data of C-3 and C–5 of metabolites 2 and 3 showed marked differences which may be due to differences in the stereochemistry at these positions (Additional file 2). The metabolite was finally identified as 5β-androstan–3α,17β-diol, reported earlier as a biotransformation product of testosterone by fungus Exophiala jeanselmei.
The molecular formula C19H28O3 [M + , m/z 304.2058] of metabolite 4 was deduced from the HREI-MS (calcd 304.2038). The presence of hydroxyl (3437 cm-1) and carbonyl (1667 cm-1) groups was inferred from the IR spectrum, while UV spectrum also indicated a conjugated ketone (λmax = 233 nm).
The 1H-NMR analysis of 4 displayed some new signals as compared to the starting material DHT (1). A downfield methine signal at δ 4.02 (br. s, W1/2 = 22.0 Hz) and its respective carbon at δ 68.9, indicated the hydroxylation of one of the methylene carbon atoms. A downfield olefinic proton singlet at δ 5.71 (s) and the corresponding carbon at δ 124.6 indicated the introduction of a C = C bond. The HMBC spectrum displayed long-range couplings of the hydroxyl-bearing methine proton (δ 4.02) with C-9 (δ 59.2), C-10 (δ 39.9), and C-13 (δ 43.7), which suggested the position of the -OH at C-11. The stereochemical assignments were based on NOESY interactions between H-11 (δ 4.02), H-8 (δ 1.56), Me-19 (δ 1.31) and Me-18 (δ 0.81). H-11 was thus deduced as β-oriented. The olefinic proton (δ 5.71) showed HMBC with C-10 (δ 39.9), and C-2 (δ 34.2). The enone system was thus deduced at C-3/C-4-C-5, with the olefinic proton at C-4 (Additional file 3). Metabolite 4 was identified as 11α,17β-dihydroxyandrost-4-en-3-one. Hunter et al. in 2009 reported metabolite 4 from biotransformation of testosterone with fungus Myceliophthora thermophila, in an organotypic culture which represents in vivo situation. The cultures consisted of primary rat, porcine, and human hepatocytes .
Metabolite 5 (C19H24O2M + at m/z 284.1726, calcd 284.1776) showed the presence of ketone (1730 cm-1) and a conjugated enone (1657 cm-1) in IR spectrum. UV spectrum showed a strong absorption for conjugated ketone (λmax = 243 nm).
13 C-NMR Chemical shift data of compounds 1–9, δ in ppm
Metabolite 6 (C19H24O3, M + at m/z 302.1852, calcd 302.1882) showed the UV absorption at 234 nm for a conjugated ketone. The IR spectrum of 6 showed absorptions at 1668 (C = C-C = O), 1730 (C = O), and 3451 cm-1 (OH).
A methine proton at δ 4.04 (m, W1/2 = 21.5 Hz) was observed in the 1H-NMR spectrum of 6, while its corresponding methine carbon resonated at δ 68.7 (C-11) in the 13C-NMR spectrum (Table 1). The COSY spectrum analysis indicated interactions between H-11 (δ 4.04) and H-9 (δ 1.15). The H-11 (δ 4.04), also showed HMBC cross peaks with C-13 (δ 47.9), C-10 (δ 40.0), and C-8 (δ 34.6), while NOESY cross peaks were observed between H-11 (δ 4.04), H-8 (δ 1.68), H-18 (δ 0.92), and H-19 (δ 1.31). Based on the above NOESY interactions, H-11 was assigned a β-orientation (Additional file 5). Complete spectral analysis indicated that the metabolite 6 was 11α-hydroxyandrost-4-en-3,17-dione. Koshimura and coworkers obtain compound 6, by hydroxylation of androstenedione with Gelasinospora retispora. Biotransformation of testosterone with Fusarium lini also afforded the same compound .
Molecular formula C19H24O3 (M + m/z 300.1749) was deduced from the HREI-MS of metabolite 7 (calcd 300.1725). The UV analysis suggested a conjugated enone system (λmax 244 nm), while IR spectrum showed absorptions at 3386 (OH), 1731 (C = O) and 1658 cm-1 (C = C-C = O).
The 1H-NMR spectrum of metabolite 7 showed a downfield hydroxyl-bearing methine proton signal at δ 4.08, which was assigned to H/C-11. C-11 resonated at δ 67.7 in the 13C-NMR spectrum (Table 1) of 7. Three olefinic protons appeared at δ 7.73 (d, J1,2 = 10.2 Hz), 6.14 (d, J2,1 = 10.2 Hz) and 6.08 (s), with respective carbons resonated at δ 158.6, 125.1, and 124.8 (Additional file 6). Based on the above spectral data and comparison with the literature , metabolite 7 was identified as 11α-hydroxy-androst-1,4-dien-3,17-dione, earlier obtained from the biotransformation of testosterone with Fusarium lini. A strain of Trichoderma hamatum was utilized for the production of 7 from testosterone, testosterone propionate, androstenedione, 1-dehydrotestosterone, dianabol and progesterone .
Metabolite 8 (C19H24O3, M + m/z 300.1737, calcd 300.1725) showed UV absorption at 245 nm. The IR spectrum showed absorptions at 3389 (OH), 1730 (C = O) and 1658 cm-1 (C = C-C = O).
Metabolites 7 (Table 1) and 8 (Table 1) have structural similarities as evident from their NMR spectra. A methine proton at δ 3.73 (br. s, W1/2 = 13.0 Hz) was observed in the 1H-NMR of 8 (C at δ 86.0). The enone system in ring A of compound 8 showed the same HMBC interactions as in 7. The C-17 ketone carbon (δ 216.2) showed HMBC correlation with the downfield C-12 methine proton (δ 3.73). The OH-bearing methine C-12 (δ 86.0) showed HMBC correlations with H-9 (δ 1.26), H-14 (δ 1.48), and CH3-18 (δ 0.78). The H-12 (δ 3.73) showed NOESY cross peaks with H-14 (δ 1.48), and H-9 (δ 1.26), and thus supported an α-orientation of H-12 (Additional file 7). The metabolite 8 was finally identified as 12β-hydroxyandrost-1,4-dien-3,17-dione, earlier isolated from Halorrhena wulfsbergii leaves . Deoxycholate and β-sitosterol were previously converted to compound 8 by using immobilized Pseudomonas sp. and a thermophilic bacterium, respectively .
Metabolite 9 (C19H24O3, M + = m/z 300.1749, calcd 300.1725), showed UV of 245 nm, while the IR spectrum displayed the same pattern as 8 [(3388 (OH), 1730, 1657 cm-1 (C = O)) for hydroxyl, ketone and enone functionalities, respectively].
The comparison of 1H-NMR spectra of 9 and 8 showed similarities, except for a methine signal at δ 4.36 (br s, W1/2 = 22.3 Hz) in compound 9 which appeared at δ 3.73 in metabolite 8. The 13C-NMR spectrum (Table 1) of 9 showed this new methine carbon at δ 71.1. The proton, geminal to the -OH group (δ 4.36), was correlated with C-13 (δ 38.7), C-14 (δ 47.3) and C-17 (δ 218.3) in HMBC spectrum. The methine C-15 (δ 71.1) showed HMBC correlations with H-14 (δ 1.50) and H2-15 (δ 1.87, 1.98). Based on the above observations, the hydroxyl-bearing methine carbon was assigned to C-16. The H-16 (δ 4.36) showed NOESY cross peaks with H-18 (δ 1.01), but no interaction with H-14 (δ 1.50). Therefore the C-16 proton was deduced to be β-oriented (Additional file 8). The metabolite 9 was thus identified as 16α-hydroxyandrost-1,4-dien-3,17-dione. Numazwama et. al., prepared 16α-hydroxyandrost-1,4-dien-3,17-dione (9) from androst-1,4-dien-3,17-dione upon reaction with hypervalent iodine .
Enzyme inhibition assay
BChE Inhibitory activities (IC 50 ) of compounds 1-7
Mean ± SEMa[ μM]
109.4 ± 1
Galanthamine (positive control)
4.9 ± 0.3
Molecular docking simulation
Re-docking score of co-crystallized ligand butyrylthiocholine (BCh) inside the binding pocket of BChE (PDB ID 1P0P)
Binding free energy (Kcal/mol)
Molecular docking studies of biotransformed DHT derivatives inside the binding cavity of BChE (PDB ID 1P0P)
IC50 ± SEM [μM]
Binding free energy (Kcal/mol)
109.4 ± 1
4.9 ± 0.3
The Figure 6c clearly reflects that no interactions were found between compound3 and binding residues of BChE, which might be due to the inverted orientation. Compound 3 exhibited no inhibition of the enzyme during the experimental studies. Similarly, compounds 1 and 2 also exhibited weak potency against BChE. By comparison of docked pose orientation of all compounds, it was deduced that compound 2 (Figure 6b) has a better orientation than compounds 1, and 4–7. With this orientation, two hydrogen bonds formation were also possible with the TYR128 (2.89 Å) and TYR332 (2.17 Å) residues of BChE, but the molecule would be disorientated for inhibition. This disturbance could be due to the presence of hydroxyl group (OH), instead of carbonyl moiety, and the absence of a double bond in ring “A” that was required for π-π interactions. Compound 4 was found to be more active as compared to 1–3, while less active than compounds 5–7. It would be due to the absence of hydrogen bonding with TYR332 (Figure 6d) even though it is involved in π-π interaction with TRP82. Compound 7 exhibited the highest inhibitory potency, experimentally and theoretically (-6.82 Kcal/mol), due to the double bond in ring “A”, as well as hydrogen bond interactions with TYR128 (3.19 Å) and TYR332 (3.05 Å) and hydrophobic interactions with key residues.
Docking results of compounds 2–7 with BChE provided valuable information about the nature of the binding interactions that were satisfactorily correlated with the experimental studies. This information could be utilized to design new leads against the BChE.
DHT was purchased from Acros-organics (Belgium). Thin layer chromatography was performed on precoated silica gel plates (PF254), purchased from Merck (Germany). Flash silica (E. Merck, Germany) was used for conducting column chromatography (CC). GS-320 (size exclusion) column was used for purification of compounds on preparative HPLC. 1H- and 13C-NMR experiments were conducted in CDCl3 and CD3OD on BrukerAvance-NMR spectrophotometers. The chemical shifts (δ values) are presented in ppm and the coupling constants (J values) in Hertz. JEOL JMS-600H mass spectrometer (Japan) was used for recording EI-MS in m/z (rel. %).
Microorganisms and culture medium
Macrophomina phaseolina (KUCC 730), obtained from Karachi University Culture Collection (KUCC), was grown on Sabouraud dextrose-agar (SDA) at 25°C and stored at 4°C. The medium for M. phaseolina was prepared by dissolving glycerol (10.0 mL), KH2PO4 (5.0 g), glucose (10.0 g), yeast extract (5.0 g), NaCl (5.0 g), and peptone (5.0 g) per litre of distilled H2O.
Gibberella fujikuroi (ATCC 10704) was purchased from the American Type Culture Collection. One litre of fermentation medium for G. fujikuroi was composed of glucose (20.0 g), NH4NO3 (1.0 g), KH2PO4 (5.0 g), MgSO4.7H2O (1.0 g) and trace element solution (2.0 mL). The trace element solution comprised of Co(NO3)3 (0.01 g), FeSO4.7H2O (0.1 g), CuSO4.5H2O (0.1 g), ZnSO4.7H2O (0.16 g), MnSO4. 7H2O (0.01 g), Mo(NH4)3 (0.01 g) in 100 mL distilled water.
Fermentation of DHT (1)
Fermentation with Macrophomina phaseolina
The fermentation medium for M. phaseolina was prepared by following the procedures presented above. 4.0 Litres of the medium was prepared and distributed (100 mL each) among 40 conical flasks (250 mL). The flasks were plugged with cotton and autoclaved at 121°C. The sterilized medium was innoculated with spores of M. phaseolina, obtained from already prepared slants on SDA. The flasks were placed on rotary shaker, and when enough growth was achieved, the flasks were fed uniformly with 1 g of DHT (1), dissolved in 20 mL acetone under aseptic conditions. The flasks were left over rotary shaker for 6 days. The fungal mass was filtered off. The filtrate was extracted with 12 L of dichloromethane (DCM). The DCM extract was dried over anhydrous Na2SO4, and concentrated on a rotary evaporator to obtain a brown gummy material (1.4 g). This gum was fractionated on a silica gel column with gradient petroleum ether- acetone solvent system. Three main fractions were obtained, which on further column chromatography led to the purification of metabolites 2 and 3.
Microbial transformation of DHT (1) with M. phaseolina
The 13C-NMR chemical shifts of compounds 2 and 3 are shown in Table 1. Other data is presented below:
5α-Androstan-3β,17β-diol (2). Colorless crystalline solid. M.p.: 153–155°C (154–157°C ). [α] 25 D : -72.0, c = 0.034, MeOH (+6.20, c = 0.55, CHCl3). 1 H-NMR (CD3OD, 400 MHz): δ 0.97, m, 1.72a, m (H-1); 1.38, m, 1.73a, m (H-2); 3.49, m (H-3); 1.26, m, 1.50, m (H-4); 1.10, m (H-5); 1.28a, m, 1.31a, m (H-6); 0.89, m, 1.68, m (H-7); 1.40, m (H-8); 0.64, ddd (J9a,8a = 14.8 Hz, J9a,11a = 10.8 Hz, J9a,11e = 4.4 Hz) (H-9); 1.30a, m, 1.56a, m (H-11); 1.02, m, 1.79, m (H-12); 0.92, m (H-14); 1.24, m, 1.57a, m (H-15); 1.45, m, 1.95, m (H-16); 3.54, t (J17a,16a/16e = 8.4 Hz) (H-17); 0.71, s (H-18); 0.83, s (H-19), (a = exchangeable assignments). 13 C-NMR (CD3OD, 150 MHz): See Table 1, HREI-MS m/z (mol. formula, calcd value): 292.2434 (C19H32O2, 292.2402). EI-MS m/z (rel. int., %): 292 [M+] (100), 277 (49), 259 (29), 248 (25), 233 (88), 215 (82), 201 (14), 175 (9), 165 (52), 121 (31), 107 (41), 55 (12).
5β-Androstan-3α,17β-diol (3). Colorless crystalline solid. M.p.: 229–231°C (230–231°C ). [α] 25 D : + 23.0, c = 0.028, MeOH (+ 26.0, ethanol ). 1 H-NMR (CD3OD, 400 MHz): δ 1.66, m, 1.82, dt (J1a,1e = 12 Hz, J1a,2a/2e = 3.6 Hz) (H-1); 1.35, m; 1.69, m (H-2); 3.94, t (J = 2.4 Hz) (H-3); 1.45a, m; 1.59, m (H-4); 1.48a, m (H-5); 1.28a, m; 1.31, m (H-6); 1.27a, m; 1.39, m (H-7); 1.21, m (H-8); 0.74, m (H-9); 1.41, m; 1.47a, m (H-11); 0.99, m; 1.93, m (H-12); 0.89, m (H-14); 1.27a, m; 1.43, m (H-15); 1.46a, m; 1.56, m (H-16); 3.55, t (J17a,16a/16e = 8.8 Hz) (H-17); 0.71, s (H-18); 0.82, s (H-19), (a = exchangeable assignments). 13 C-NMR (CD3OD, 100 MHz): See Table 1. HREI-MS m/z (mol. formula, calcd value): 292.2434 (C19H32O2, 292.2402). EI-MS m/z (rel. int., %): 292 [M+] (66), 277 (51), 259 (33), 241 (20), 217 (54), 215 (100), 201 (14), 173 (11), 165 (45), 147 (33), 107 (47), 93 (48), 55 (37).
Fermentation with Gibberella fujikuroi
Two stage fermentation protocol (discussed above) was used for the fermentation of 1 (1.0 g) with G. fujikuroi. Fermentation medium was prepared as discussed earlier. Fermentation continued for 7 days. All the flasks were filtered, extracted with 12 litres of DCM and evaporated in vacuo. A brown extract (2.5 g) was obtained which was fractionated over a silica gel column with petroleum ether- acetone (P.E.- acetone) gradient solvent system. The fractions were compiled based on thin layer chromatography (TLC) profile to obtain five main fractions, which were again subjected to column chromatography (C.C.) (P.E.-acetone) and size exclusion high performance liquid chromatography (HPLC) by using a GS-320 column, and methanol as solvent system to afford pure metabolites 4–9.
Microbial transformation of DHT (1) with G. fujikuroi
The 13C-NMR chemical shifts of compounds 4–9 are presented in Table 1. Other data is presented below:
11α,17β-Dihydroxyandrost-4-en-3-one (4). Colorless crystalline compound. M.p.: 179–181°C (178–182°C ). [α] 25 D : + 95.0, c = 0.096, MeOH (+94, c = 0.1, CHCl3). 1 H-NMR (CDCl3, 400 MHz): δ 2.62, dt (J1a,1e = 14.0 Hz, J1a,2a/2e = 4.0 Hz), 1.98, td (J1e,2a/2e = 13.6 Hz, J1e,1a = 4.0 Hz); 2.41, m; 2.31, m (H-2); 5.71, s (H-4); 2.37, m; 2.27, m (H-6); 1.81, m; 1.02, m (H-7); 1.56a, m (H-8); 1.05a, m (H-9); 4.02, br. s (W1/2 = 22.0 Hz) (H-11); 1.17, m; 2.13, m (H-12); 1.07a, m (H-14); 1.58a, m; 1.27, m (H-15); 1.45, m; 2.08, m (H-16); 3.77, t (J17a,16a/16e = 8.0 Hz) (H-17); 0.81, s (H-18); 1.31, s (H-19), (a = exchangeable assignments). 13 C-NMR (CDCl3, 150 MHz): See Table 1. HREI-MS m/z (mol. formula, calcd value): 304.2058 (C19H28O3, 304.2038). EI-MS m/z (rel. int., %): 304 [M+] (74), 286 (61), 253 (15), 180 (43), 163 (100), 147 (29), 137 (34), 124 (98), 109 (41), 91 (44), 79 (34), 55 (20).
Androst-1,4-dien-3,17-dione (5). Colorless crystals. M.p.: 140–142°C (141–142°C ). [α] 25 D : + 115.5 (c = 0.028, MeOH) (+117, CHCl3). 1 H-NMR (CDCl3, 400 MHz): δ 7.03, d (J1,2 = 10.4 Hz) (H-1); 6.22, d (J2,1 = 10.4 Hz) (H-2); 6.07, s (H-4); 2.41a, m, 2.48, m (H-6); 1.12, m, 2.06, m (H-7); 1.80, m (H-8); 1.09, m (H-9); 1.67, m, 1.84, m (H-11); 1.27a, m, 1.87, m (H-12); 1.26a, m (H-14); 1.57, m, 1.95, m (H-15); 2.10, m, 2.43a, m (H-16); 0.92, s (H-18); 1.24, s (H-19), (a = exchangeable assignments). 13 C-NMR (CDCl3, 150 MHz): See Table 1. HREI-MS m/z (mol. formula, calcd value): 284.1726 (C19H24O2, 284.1776). EI-MS m/z (rel. int., %): 284 [M+] (45), 266 (8), 227 (6), 171 (8), 159 (66), 150 (26), 135 (31), 122 (100), 107 (36), 91 (46), 79 (19), 67 (10), 55 (10).
11α-Hydroxyandrost-4-en-3,17-dione (6). Colorless crystalline compound. M.p.: 238–239°C (240–241°C ). [α] 25 D : + 162.0, c = 0.034, MeOH (+ 165, CHCl3). 1 H-NMR (CDCl3, 400 MHz): δ 2.05, m; 2.65, dt (J1a,1e = 13.6 Hz, J1a,2a/2e = 4.4 Hz) (H-1); 2.31a, m; 2.39, m (H-2); 5.72, s (H-4); 2.29a, m; 2.37, m (H-6); 1.12, m; 1.97a, m (H-7); 1.68, m (H-8); 1.15, m (H-9); 4.04, m (W1/2 = 21.5 Hz) (H-11); 1.30, m; 2.11a, m (H-12); 1.37, m (H-14); 1.53, m; 1.95a, m (H-15); 2.14a, m; 2.48, m (H-16); 0.92, s (H-18); 1.31, s (H-19), (a = exchangeable assignments). 13 C-NMR (CDCl3, 125 MHz): See Table 1. HREI-MS m/z (mol. formula, calcd value): 302.1852 (C19H26O3, 302.1882). EI-MS m/z (rel. int., %): 302 [M+] (24), 290 (35), 280 (21), 248 (51), 230 (95), 215 (21), 198 (19), 186 (28), 159 (31), 139 (100), 138 (88), 136 (52), 105 (29), 87 (65), 55 (18).
11α-Hydroxyandrost-1,4-dien-3,17-dione (7). Colorless crystalline compound. M.p.: 208–210°C (212–214°C ). [α] 25 D : + 85.1, c = 0.037, MeOH (+ 86.5 ). 1 H-NMR (CDCl3, 300 MHz): δ 7.73, d (J1,2 = 10.2 Hz) (H-1); 6.14, d (J2,1 = 10.2 Hz) (H-2); 6.08, s (H-4); 2.41, m, 2.48, m (H-6); 1.15, m, 2.07, m (H-7); 1.82, ddd (J8a,9a/14a = 22.8 Hz, J8a,7a = 11.1 Hz, J8a,7e = 3.9 Hz) (H-8); 1.12, m (H-9); 4.08, m (W 1/2 = 19.6 Hz) (H-11); 1.25, m, 2.22, m (H-12); 1.35, m (H-14); 1.56, m, 1.95, m (H-15); 2.15, m, 2.50, m (H-16); 0.93, s (H-18); 1.31, s (H-19). 13 C-NMR (CDCl3, 150 MHz): See Table 1. HREI-MS m/z (mol. formula, calcd value): 300.1749 (C19H24O3, 300.1725). EI-MS m/z (rel. int., %): 300 [M+] (36), 282 (18), 231 (80), 161 (20), 124 (37), 109 (32), 84 (58), 55 (23).
12β-Hydroxyandrost-1,4-dien-3,17-dione (8). Colorless crystalline compound. M.p.: 163–165°C (164–166°C ). [α] 25 D : + 84.0, c = 0.085, MeOH (+ 83.0, CHCl3). 1 H-NMR (CDCl3, 400 MHz): δ 7.03, d (J1,2 = 10.0 Hz) (H-1); 6.23, dd (J2,1 = 10.0 Hz, J2,4 = 1.2 Hz) (H-2); 6.03, br s (H-4); 2.39, m, 2.48, td (13.6, 4.4) (H-6); 1.15, m, 1.92, m (H-7); 1.81, m (H-8); 1.26, m (H-9); 1.42, m, 2.03, m (H-11); 3.73, br s (W 1/2 = 13.0 Hz) (H-12); 1.48, m (H-14); 1.75, m, 1.86, m (H-15); 1.90, m, 2.33, m (H-16); 0.78, s (H-18); 1.25, s (H-19). 13 C-NMR (CDCl3, 150 MHz): See Table 1. HREI-MS m/z (mol. formula, calcd value): 300.1737 (C19H24O3, 300.1725). EI-MS m/z (rel. int., %): 300 [M+] (66), 282 (8), 231 (100), 161 (20), 124 (37), 109 (32), 84 (58), 55 (43).
16α-Hydroxyandrost-1,4-dien-3,17-dione (9). Colorless crystals. M.p.: 152–153°C. [α] 25 D : + 89.3 (c = 0.043, MeOH). 1 H-NMR (CDCl3, 400 MHz): δ 7.00, d (J1,2 = 10.0 Hz) (H-1); 6.22, d (J2,1 = 10.0 Hz) (H-2); 6.04, s (H-4); 2.35, m, 2.45, m (H-6); 1.32, m, 2.03, m (H-7); 1.77, m (H-8); 1.10, m (H-9); 1.67, m, 1.85, m (H-11); 1.37, m, 1.65, m (H-12); 1.50, m (H-14); 1.87, m, 1.98, m (H-15); 4.36, br s (W 1/2 = 20.8 Hz) (H-16); 1.01, s (H-18); 1.24, s (H-19). 13 C-NMR (CDCl3, 150 MHz): See Table 1. HREI-MS m/z (mol. formula, calcd value): 300.1749 (C19H24O3, 300.1725). EI-MS m/z (rel. int., %): 300 [M+] (100), 272 (36), 246 (69), 239 (43), 216 (39), 200 (38), 190 (24), 165 (33), 147 (55), 108 (63), 92 (31), 81 (17).
BChE inhibitory assay
All experiments were performed in a 96-well microplate in triplicate on SpectraMax340 (Molecular Devices, CA, U. S. A.).
Equine serum BChE, butyrylcholine chloride, and 5,5′-dithio-bis (2-nitrobenzoic) acid were purchased from Sigma. Analytical grade buffers and other chemicals were used. The BChE inhibition was measured by a modified spectrophotometric method . Butyrylcholine chloride was used as substrate and the BChE inhibitory activity was measured by using 5-5′-dithio bis (2-nitrobenzoic) acid (DTNB).
150 μL of 100 mM sodium phosphate buffer [PBS] (pH 8), 10 μL of test compound solution (0.2 mM), and 20 μL BChE solution were mixed together and incubated at 25°C for 15 minutes. The reaction was initiated by the addition of 10 μL DTNB and 10 μL of BCh. A yellow 5-thio-2-nitrobenzoate anion was formed by hydrolysis of the substrate as a result of the reaction of DTNB with thiocholine which was monitored spectrophotometrically at 412 nm. MeOH was used as a solvent to dissolve the test compounds and the controls.
Determination of IC50
IC50 values were calculated by measuring the effect of varying concentrations of test compounds by using the EZ-Fit Enzyme Kinetic program. Galanthamine was used as a standard inhibitor of BChE.
Molecular modeling studies were conducted on a dual processor Intel® Xeon™ CPU 3.00 GHz LINUX work station running under SUSE 11.4. For the selection of best docking program with respect to the current biological system, GOLD , MOE  and Surflex [43, 44] docking programs were employed. Molecular Operating Environment (MOE) docking software was finally selected to study the assembly pattern of biotransformed DHT derivatives in complex with human BChE system. Among thirty-one X-ray crystal structures of human BChE in the protein data bank (PDB) [11, 45], PDB ID 1P0P: 2.30 Å was selected as the target protein based on suitable resolution and co-crystallized ligand, BCh. The entire system was energy minimized by MMFF94 force field , after adding the missing hydrogen atoms and keeping heavy atoms fixed until a Root-mean-square deviation (RMSD) gradient of 0.05 Kcalmol-1 Å-1 was reached. 3D structures of all metabolites 2–7 were drawn by molecule builder which is incorporated in MOE modeling package and the structures were subjected to MMFF94 for energy minimization. Subsequently for the evaluation of potential energy, partial charges were calculated by MMFF94 force field . Both prepared systems (protein and ligand) were introduced for molecular docking simulation. Docking simulations were performed by using two methods (a) Alpha Triangle placement method, and (b) Triangle matcher placement method. All compounds were ranked with London dG scoring function and re-scored by GBVI/WSA dG with the force field refinement strategy . A total of 30 docking poses were generated for each ligand and the pose with the lowest energy was selected for further studies.
In summary, we report here the microbial transformation of dihydrotestosterone (DHT, 1) for the first time by using suspension cultures of M. phaseolina and G. fujikuroi. Compounds 5–7 were found to be significant and specific inhibitors of the BChE, in comparison to standard drug, galanthamine. Current experimental studies revealed that the DHT derivatives actively participated in the inhibition of BChE. Additionally, experimental studies correlated well with theoretical studies thus indicating the behavior of different metabolites inside the binding pocket of BChE, and supported inhibition activity trend, followed by least active to highly active compounds. From the docking analysis, experimental inhibitory potency relationship of compounds 1–7 exhibited the fundamental role of α,β-unsaturated carbonyl moiety in ring “A”. In the bound state of metabolites with BChE, compounds 5–7 exhibited favorable hydrophilic, hydrophobic and hydrogen bond interactions with active site residues of the receptor protein. These compounds also participated in π-π interaction towards the TRP82.
MIC envisioned the concept of the current study, supervised the practical work, helped in solving spectroscopic data and finalized the manuscript. SZ carried out the microbial transformation experiments, purified all the metabolites, interpreted the spectroscopic data and prepared the preliminary draft of the manuscript. KD conducted the biochemical assays for inhibition of BChE. UM conducted the in silico studies and wrote the simulation part of the manuscript. ZH supervised the computational studies. All authors read and approved the final manuscript.
One of the authors, S. Zafar, acknowledges the enabling role of the Higher Education Commission (HEC), Pakistan, by awarding HEC indigenous Ph. D. scholarship. Co-authors are thankful to the Pakistan Science Foundation (PSF) for their financial support to purchase the MOE software for computational studies. This research study was supported by the Higher Education Commission’s grant for the project entitled, “High resolution x-ray analysis of pharmaceuticals by important enzymes in complex with plant-based inhibitors as a basis for rational drug design” (20-1364/R&D/09) to MIC.
- Mahato SB, Garai S: Advances in microbial steroid biotransformation. Steroids. 1997, 62: 332-345. 10.1016/S0039-128X(96)00251-6.View ArticleGoogle Scholar
- Choudhary MI, Zafar S, Khan NT, Ahmed S, Noreen S, Marisini B, Al-Khedairy A, ur-Rahman A: Biotransformation of dehydroepiandrosterone with Macrophomina phaseolina and β-glucuronidase inhibitory activity of transformed products. J Enzyme Inhib Med Chem. 2012, 27: 348-355. 10.3109/14756366.2011.590804.View ArticleGoogle Scholar
- Anderson KM, Liao S: Selective retention of dihydrotestosterone by prostatic nuclei. Nature. 1968, 219: 277-279. 10.1038/219277a0.View ArticleGoogle Scholar
- Grino PB, Grifein JE, Wilson JD: Testosterone at high concentrations interacts with the human androgen receptor similarly to dihydrotestosterone. Endocrinology. 1990, 126: 1165-1172. 10.1210/endo-126-2-1165.View ArticleGoogle Scholar
- Vierhapper H, Maier H, Nowotny P, Waldhausl W: Production rates of testosterone and of dihydrotestosterone in female pattern hair loss. Metabolism. 2003, 52: 927-929. 10.1016/S0026-0495(03)00060-X.View ArticleGoogle Scholar
- Kalow W: Pharmacogenetics. 1962, Philadelphia: W. B. Saunders, Chapter 6Google Scholar
- Lehmann H, Liddell J: Medical Genetics. 1964, New York: Grune & Stratton, Chapter 3Google Scholar
- Greig NH, Utsuki T, Ingram DK, Wang Y, Pepeu G, Scali C, Yu Q, Mamczarz J, Holloway HW, Giordano T, Chen D, Furukawa K, Sambamurti K, Brossi A, Lahiri DK: Selective butyrylcholinesterase inhibition elevates brain acetylcholine, augments learning and lowers Alzheimer beta-amyloid peptide in rodent. Proc Natl Acad Sci. 2005, 102: 17213-17218. 10.1073/pnas.0508575102.View ArticleGoogle Scholar
- Feher M, Williams CIJ: Effect of input differences on the results of docking calculations. J Chem Inf Model. 2009, 49: 1704-1714. 10.1021/ci9000629.View ArticleGoogle Scholar
- Onodera K, Satou K, Hirota HJ: Evaluation of molecular docking programs for virtual screening. J Chem Inf Model. 2007, 47: 1609-1618. 10.1021/ci7000378.View ArticleGoogle Scholar
- Nicolet Y, Lockridge O, Masson P, Fontecilla-Camps JC, Nachon F: Crystal structure of human butyrylcholinesterase and of its complexes with substrates and products. J Biol Chem. 2003, 278: 41141-41147. 10.1074/jbc.M210241200.View ArticleGoogle Scholar
- Bellavance E, Luu-The V, Poirier DJ: Potent and selective steroidal inhibitors of 17beta-hydroxysteroid dehydrogenase type 7, an enzyme that catalyzes the reduction of the key hormones estrone and dihydrotestosterone. J Med Chem. 2009, 52: 7488-7502. 10.1021/jm900921c.View ArticleGoogle Scholar
- Norata GD, Cattaneo P, Poletti A, Catapano AL: The androgen derivative 5alpha-androstane-3beta,17beta-diol inhibits tumor necrosis factor alpha and lipopolysaccharide induced inflammatory response in human endothelial cells and in mice aorta. Atherosclerosis. 2010, 212: 100-106. 10.1016/j.atherosclerosis.2010.05.015.View ArticleGoogle Scholar
- Porter RBR, Gallimore WA, Reese PB: Steroid transformations with Exophiala jeanselmei var. lecanii-corni and Ceratocystis paradoxa. Steroids. 1999, 64: 770-779. 10.1016/S0039-128X(99)00066-5.View ArticleGoogle Scholar
- Hunter C, Watts KR, Dedi C, Dodd HT: An unusual ring-a opening and other reaction in steroid transformation by the thermophilic fungus Myceliophthora thermophila. J Steroid Biochem Mol Biol. 2009, 116: 171-177. 10.1016/j.jsbmb.2009.05.012.View ArticleGoogle Scholar
- Langsch A, Giri S, Acikgoz A, Jasmund I, Frericks B, Bader A: Interspecies difference in liver-specific functions and biotransformation of testosterone of primary rat, porcine and human hepatocyte in an organotypical sandwich culture. Toxicol Lett. 2009, 188: 173-179. 10.1016/j.toxlet.2009.03.022.View ArticleGoogle Scholar
- Pozo OJ, Marcos J, Ventura R, Fabregat A, Segura J: Testosterone metabolism revisited: discovery of new metabolites. Anal Bioanal Chem. 2010, 398: 1759-1770. 10.1007/s00216-010-4082-0.View ArticleGoogle Scholar
- Fahrbach M, Krauss M, Preiss A, Kohler HE, Hollender J: Anaerobic testosterone degradation in steroidobacter denitrificans-identification of transformation products. Environ Pollut. 2010, 158: 2572-2581. 10.1016/j.envpol.2010.05.017.View ArticleGoogle Scholar
- Wei W, Fan S, Wang F, Wei D: A new steroid-transforming strain of Mycobacterium neoaurum and cloning of 3-ketosteroid 9alpha-hydroxylase. Appl Biochem Biotechnol. 2010, 162: 1446-1456. 10.1007/s12010-010-8919-y.View ArticleGoogle Scholar
- Zhao Y, Gong P: Mutant an optimization of androstenonetranslation strain Mycobacterium sp. SH5. Shipi Yu Shengwu Jishu Xuebao. 2010, 29: 150-Google Scholar
- Chaudhari PN, Chaudhari BL, Chincholkar SB: Cholesterol biotransformation to androsta-1,4-diene-3,17-dione by growing cells of Chryseobacterium gleum. Biotechnol Lett. 2010, 32: 695-699. 10.1007/s10529-010-0206-z.View ArticleGoogle Scholar
- Koshimura M, Utsukihara T, Hara A, Mizobuchi S, Horiuchi CA, Kuniyoshi MJ: Hydroxylation of steroid compounds by Gelasinospora retispora. J Mol Catal B-Eznym. 2010, 67: 72-77. 10.1016/j.molcatb.2010.07.008.View ArticleGoogle Scholar
- Al-Aboudi A, Mohammad MY, Musharraf SG, Choudhary MI, ur-Rahman A: Microbial transformation of testosterone by Rhizopus stolonifer and Fusarium lini. Nat Prod Res. 2008, 22: 1498-1509. 10.1080/14786410802234528.View ArticleGoogle Scholar
- Bartmanska A, Dmochowska-Gładysz J: Transformation of steroids by Trichoderma hamatum. Enzyme Microb Technol. 2007, 40: 1615-1621. 10.1016/j.enzmictec.2006.11.011.View ArticleGoogle Scholar
- Nelle S, Charles G, Cave A, Goutarel R: Sciences Chimiques. 1970, 27: 153-Google Scholar
- Graham DLV, Lowe GM, Phipps DA, Bilton RF: Two novel methods used in biodegradation of bile acids and plant sterols. Biochem Soc Trans. 1991, 19: 52S-View ArticleGoogle Scholar
- Numazwama M, Mutsumi A, Ogata M: Reactions of enolizable steroidal 4-en-3-ones and 17-ones with hypervalent iodine. Chem Pharm Bull. 1988, 36: 3381-3386. 10.1248/cpb.36.3381.View ArticleGoogle Scholar
- Zurad EG: New treatments for Alzheimer’s disease. Drug Benefit Trends. 2001, 13: 27-40.Google Scholar
- Stpankova S, Komers K: Cholinesterases and cholinesterase inhibitors. Curr Enzym Inhib. 2008, 4: 160-171. 10.2174/157340808786733631.View ArticleGoogle Scholar
- De-Ruggieri P, Gandolfi C, Guzzi U: Steroids 28. 10 Beta amino steroids: A new class of steroid derivatives. Tetrahedron Lett. 1967, 23: 2195-2200. 10.1016/0040-4020(67)80055-3.View ArticleGoogle Scholar
- Wang X, Liu H, Yan P, Liu J, Li Y, Sun Q, Wang C: Simultaneously rapid deprotection of 3-acyloxy groups and reduction of D-ring ketones (nitrile) of steroids using DIBAL-H/NiCl2. J Chem Res. 2011, 35: 291-10.3184/174751911X13050949941793.View ArticleGoogle Scholar
- Fieser LB, Huang W: Configuration of steroid bromoketones; a correction. J Am Chem Soc. 1953, 75: 4837-4838. 10.1021/ja01115a509.View ArticleGoogle Scholar
- Kabasakalian P, McGlotten J, Basch A, Yudis MD: The stereoselective electrochemical reduction of nonconjugated steroidal ketones and α-ketols. J Org Chem. 1961, 26: 1738-1744. 10.1021/jo01065a011.View ArticleGoogle Scholar
- Meister PD, Adolph W: US Patent 2,877,162. 11alpha-hydroxylation of steroids by sporotrichum. 1959, Kalamazoo, Mich: Upjohn CoGoogle Scholar
- Choudhary MI, Nasir M, Khan SN, Atif M, Ali RA, Khalil SM, ur-Rahman A: Microbial hydroxylation of hydroxyprogesterones and alpha-glucosidase inhibition activity of their metabolites. Z Naturforsch. 2007, 62: 593-Google Scholar
- Choudhary MI, Sultan S, Khan MTH, Yasin A, Shaheen F, Atta-ur-Rahman : Biotransformation of (+)-androst-4-ene-3,17-dione. Nat Prod Res. 2004, 18: 529-535. 10.1080/14786410310001620628.View ArticleGoogle Scholar
- Faramarzi MI, Yazdi MT, Jahandar H, Amini M, Monsef-Esfahani HR: Studies on the microbial transformation of androst-1,4-dien-3,17-dione with Acremonium strictum. J Ind Microbiol Biotechnol. 2007, 33: 725-733.View ArticleGoogle Scholar
- Bernstein S, Littell R, Williams JH: Steroidal cyclic ketals. iv.1 the conversion of 11-keto- to 11α-hydroxysteroids. the preparation of 11-epi-hydrocortisone, and δ4-androstene-11α-ol-3,17-dione. J Am Chem Soc. 1953, 75: 1481-1484. 10.1021/ja01102a509.View ArticleGoogle Scholar
- Raspe G, Kieslich K: 12β-Hydroxy steroid. Naturwissenschaften. 1961, 48: 479-10.1007/BF00590379.View ArticleGoogle Scholar
- Ellman GL, Courtney KD, Andres V, Featherstone RM: A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem Pharmacol. 1961, 7: 88-95. 10.1016/0006-2952(61)90145-9.View ArticleGoogle Scholar
- Cambridge Crystallographic Database: Gold version 5.1. 2008, Cambridge, U.KGoogle Scholar
- Chemical Computing Group Inc: Molecular Operating Environment (MOE). 2011, 010 Sherbrook St. West, Suite # 910, Montreal, QC, Canada, H3A 2R7: Chemical Computing Group Inc, 10Google Scholar
- Biopharmics LLC: Surflex, version 2.11. 2007, San Mateo, CA, USA: Biopharmics LLCGoogle Scholar
- Jain AN: Surflex: fully automatic flexible molecular docking using a molecular similarity-based search engine. J Med Chem. 2003, 46: 499-511. 10.1021/jm020406h.View ArticleGoogle Scholar
- Berman HM: The protein data bank. Nucleic Acids Research. 2000, 28 (1): 235-242. 10.1093/nar/28.1.235.View ArticleGoogle Scholar
- Halgren TA: MMFF VI. MMFF94s option for energy minimization studies. J Comp Chem. 1999, 20: 720-729. 10.1002/(SICI)1096-987X(199905)20:7<720::AID-JCC7>3.0.CO;2-X.View ArticleGoogle Scholar
- Halgren TA: Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J Comp Chem. 1996, 17: 490-519. 10.1002/(SICI)1096-987X(199604)17:5/6<490::AID-JCC1>3.0.CO;2-P.View ArticleGoogle Scholar
- Corbeil CR, Williams CI, Labute PJ: Variability in docking success rates due to dataset preparation. Comput Aid Mol Des. 2012, 26: 775-786. 10.1007/s10822-012-9570-1.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.