Mapping of Steroids Binding to 17β-Hydroxysteroid Dehydrogenase Type 1 Using Monte Carlo Energy Minimization Reveals Alternative Binding Modes

Investigation of structure-activity relationships in organophosphates-cholinesterase interaction using docking analysis

It is known than the most potent homologues in various series of O,O-dialkylphosphates are the dibutyl or diamyl derivatives toward mammalian cholinesterases (ChEs) (both Acetyl- and Butyryl-ChEs), and the dimethyl or diethyl ones toward insect AChEs. To investigate the ChE interaction with organophosphorus inhibitors (OPIs) in more detail, we have performed in silico docking of the series of O,O-dialkylfluorophosphates into active center of different ChEs - both from mammals (human and mouse AChEs and horse BChE), and from insects (spring grain aphid AChE belonging to AChE-1 type, and housefly AChE belonging to AChE-2 type). According to the modeling results, one radical is directed to the anionic site W84, another to the acyl pocket. In addition to well-known residues 288 and 290 (Torpedo AChE sequence numbering), we showed an essential influence of residue 400 - a short alkyl residue in mammalian ChEs and phenylalanine in insect ChEs. Phenylalanine in this position creates sterical hindrance for proper orientation of the OPI molecule, which increases the distance between the catalytic serine gamma-oxygen and phosphorus, and decreases the angle of nucleophylic attack. This suggestion was supported by docking of dibutylfluorophosphate into the active center of AChEs with in silico mutations. Thus, we suggest both the angle of nucleophylic attack and the distance between the catalytic serine and phosphorus atom as measures of productivity of OPI binding.

Docking flexible ligands in proteins with a solvent exposure- and distance-dependent dielectric function

Physics-based force fields for ligand-protein docking usually determine electrostatic energy with distance-dependent dielectric (DDD) functions, which do not fully account for the dielectric permittivity variance between approximately 2 in the protein core and approximately 80 in bulk water. Here we propose an atom-atom solvent exposure- and distance-dependent dielectric (SEDDD) function, which accounts for both electrostatic and dehydration energy components. Docking was performed using the ZMM program, the AMBER force field, and precomputed libraries of ligand conformers. At the seeding stage, hundreds of thousands of positions and orientations of conformers from the libraries were sampled within the rigid protein. At the refinement stage, the ten lowest-energy structures from the seeding stage were Monte Carlo-minimized with the flexible ligand and flexible protein. A search was considered a success if the root mean square deviation (RMSD) of the ligand atoms in the apparent global minimum from the x-ray structure was <2 A. Calculations on an examining set of 60 ligand-protein complexes with different DDD functions and solvent-exclusion energy revealed outliers in most of which the ligand-binding site was located at the protein surface. Using a training set of 16 ligand-protein complexes, which did not overlap with the examining set, we parameterized the SEDDD function to minimize the RMSD of the apparent global minima from the x-ray structures. Recalculation of the examining set with the SEDDD function demonstrated a 20% increase in the success rate versus the best-performing DDD function.

Dynamics of the streptavidin-biotin complex in solution and in its crystal lattice: distinct behavior revealed by molecular simulations

We present a 250 ns simulation of the wild-type, biotin-liganded streptavidin tetramer in the solution phase and compare the trajectory to two previously published simulations of the protein in its crystal lattice. By performing both types of simulations, we are able to interpret the protein's behavior in solution in the context of its X-ray structure. We find that the rate of conformational sampling is increased in solution over the lattice environment, although the relevant conformational space in solution is also much larger, as indicated by overall fluctuations in the positions of backbone atoms. We also compare the distributions of chi1 angles sampled by side chains exposed to solvent in the lattice and in the solution phase, obtaining overall good agreement between the distributions obtained in our most rigorous lattice simulation and the crystallographic chi1 angles. We observe changes in the chi1 distributions in the solution phase, and note an apparent progression of the distributions as the environment changes from a tightly packed lattice filled with crystallization media to a bath of pure water. Finally, we examine the interaction of biotin and streptavidin in each simulation, uncovering a possible alternate conformation of the biotin carboxylate tail. We also note that a hydrogen bond observed to break transiently in previous solution-phase simulations is predominantly broken in this much longer solution-phase trajectory; in the lattice simulations, the lattice environment appears to help maintain the hydrogen bond, but more sampling will be needed to confirm whether the simulation model truly gives good agreement with the X-ray data in the lattice simulations. We expect that pairing solution-phase biomolecular simulations with crystal lattice simulations will help to validate simulation models and improve the interpretation of experimentally determined structures.

A 3D QSAR model of 17beta-HSD1 inhibitors based on a thieno[2,3-d]pyrimidin-4(3H)-one core applying molecular dynamics simulations and ligand-protein docking

The 17beta-hydroxysteroid dehydrogenase type 1 (17beta-HSD1) enzyme plays a crucial role in female hormonal regulation by catalysing the NADPH-dependent reduction of the less potent estrone E1 into the biologically active estradiol E2. Because 17beta-HSD1 is a key enzyme in E2 biosynthesis, it has emerged as an attractive drug target for inhibitor development. Herein we report the plausible binding modes and a 3D QSAR model of 17beta-HSD1 inhibitors based on a (di)cycloalkenothieno[2,3-d]pyrimidin-4(3H)-one core. Two generated enzyme complexes with potent inhibitors were subjected to molecular dynamics simulation to mimic the dynamic process of inhibitor binding. A set of 17beta-HSD1 inhibitors based on the thieno[2,3-d]pyrimidin-4(3H)-one core were docked into the resulting active site, and a CoMFA model employing the most extensive training set to date was generated. The model was validated with an external test set. Active site residues involved in inhibitor binding and CoMFA fields for steric and electrostatic interactions were identified. The model will be used to guide structural modifications of 17beta-HSD1 inhibitors based on a thieno[2,3-d]pyrimidin-4(3H)-one core in order to improve the biological activity as well as in the design of novel 17beta-HSD1 inhibitors.

Comparison of crystal structures of human androgen receptor ligand-binding domain complexed with various agonists reveals molecular determinants responsible for binding affinity

Androgens exert their effects by binding to the highly specific androgen receptor (AR). In addition to natural potent androgens, AR binds a variety of synthetic agonist or antagonist molecules with different affinities. To identify molecular determinants responsible for this selectivity, we have determined the crystal structure of the human androgen receptor ligand-binding domain (hARLBD) in complex with two natural androgens, testosterone (Testo) and dihydrotestosterone (DHT), and with an androgenic steroid used in sport doping, tetrahydrogestrinone (THG), at 1.64, 1.90, and 1.75 A resolution, respectively. Comparison of these structures first highlights the flexibility of several residues buried in the ligand-binding pocket that can accommodate a variety of ligand structures. As expected, the ligand structure itself (dimension, presence, and position of unsaturated bonds that influence the geometry of the steroidal nucleus or the electronic properties of the neighboring atoms, etc.) determines the number of interactions it can make with the hARLBD. Indeed, THG--which possesses the highest affinity--establishes more van der Waals contacts with the receptor than the other steroids, whereas the geometry of the atoms forming electrostatic interactions at both extremities of the steroid nucleus seems mainly responsible for the higher affinity measured experimentally for DHT over Testo. Moreover, estimation of the ligand-receptor interaction energy through modeling confirms that even minor modifications in ligand structure have a great impact on the strength of these interactions. Our crystallographic data combined with those obtained by modeling will be helpful in the design of novel molecules with stronger affinity for the AR.

Molecular docking simulations of steroid substrates into human cytosolic hydroxysteroid dehydrogenases (AKR1C1 and AKR1C2): insights into positional and stereochemical preferences

AKR1C1 and AKR1C2 are human cytosolic hydroxysteroid dehydrogenases, which play pivotal roles in the metabolism and action of natural and synthetic steroid hormones. The two enzymes are highly homologous, and have distinct positional and stereochemical preferences with various substrates. We performed molecular docking simulations of three steroid substrates, including an androgen (5alpha-dihydrotestosterone, DHT), a progestin (progesterone, PRO), and a synthetic hormone ([7alpha,17alpha]-17-hydroxy-7-methyl-19-norpregn-5(10)-en-20-yn-3-one or tibolone, TIB), into the active sites of the two enzymes. For each substrate and enzyme pair, the activity inferred by the "productive" docking models (in which the spatial arrangement of the steroid and the cofactor would permit a reaction) matched the experimentally observed positional and stereochemical outcome. These productive conformations were energetically and statistically favored except for TIB and PRO with AKR1C2, where experimentally strong substrate inhibition and low activity were observed, respectively. Results showed that (i) a 3-ketosteroid (DHT) and a 20-ketosteroid (PRO) were reduced by AKR1C1 since the carbonyl groups could occupy the same position by "backwards" binding of steroids; (ii) 3alpha-reduced (DHT) and 3beta-reduced (TIB) products were formed by AKR1C2 since the angular methyl groups of the steroids were inverted by "upside-down" binding of steroids; and (iii) the 3beta- and 3alpha-reduction of DHT by AKR1C1 and AKR1C2, respectively occurred since the steroids employed a "swinging" motion to present opposite faces to the cofactor. Favorable nonproductive modes were observed with all substrates in both enzymes in which the steroid was bound at a "near-entry" position and/or an "in-middle" position, which may influence the reaction coordinate.