Ab Initio Quantum Mechanical Study on the Origin of the pK(a) Differences of the Proton Sponges 1,8-Bis(dimethylamino)naphthalene, 1,8-Bis(dimethylamino)-2,7-dimethoxynaphthalene, 1,6-Dimethyl-1,6-diazacyclodecane, and 1,6-Diazabicyclo[4.4.4]tetradecane

Ab initio quantum mechanical calculations were used in studying the origin of the exceptionally high basicities of four diamines (13-16) with pK(a1) values ranging from 12.1 to 25. The computational approach involved the calculation of the gas-phase proton affinities of the molecules studied at the MP2/6-31G//HF/6-31G level and the solvation energies with the polarizable continuum model at the HF/6-31G level. The calculated gas-phase and aqueous-phase proton affinities of a structurally diverse series of amines were compared with the corresponding experimental gas-phase proton affinities and pK(a1) values. The calculated values were found to be in reasonable agreement with the experimental ones. The high basicities of the studied diamines were found to originate from the nitrogen lone-pair repulsion, solvation effects, and strong intramolecular hydrogen bonds. Each of these factors were found to be able to increase the pK(a1) values of the high-basicity diamines by 2-6 pK(a) units. The relative contributions of the factors varied between the compounds. The nitrogen lone-pair repulsion was estimated to be the most important factor in increasing the pK(a1) values. In addition, barriers for proton transfers between the nitrogens of selected diamines were calculated, and comparison was made between the barrier heights and the geometries of the diamines.

Critical role of helix 12 of the vitamin D(3) receptor for the partial agonism of carboxylic ester antagonists

The carboxy-terminal alpha-helix of a nuclear receptor ligand-binding domain (LBD), helix 12, contains a critical, ligand-modulated interface for the interaction with coactivator proteins. In this study, using the example of the vitamin D receptor (VDR) and the partial antagonist ZK159222, the role of helix 12 (residues 417-427) for both antagonistic and agonistic receptor actions was investigated. Amino acid residue G423 was demonstrated to be critical for partial agonism of ZK159222, but not for the activity of the natural VDR agonist, 1alpha,25-dihydroxyvitamin D(3) (1alpha,25(OH)(2)D(3)). The amount of partial agonism of ZK159222 increased when helix 12 was truncated by the last four amino acid residues (Delta424-27) and augmented even more, when in addition helix 12 of VDR's dimerization partner, retinoid X receptor (RXR), was truncated. In contrast, the low agonism of a structural derivative of ZK159222, ZK168281, was not affected comparably, whereas other close structural relatives of ZK159222 even demonstrated the same agonistic activity as that of 1alpha,25(OH)(2)D(3). The amount of agonism of ZK159222 and ZK168281 at different variations of helix 12 correlated well with VDR's ability to complex with coactivator proteins and inversely correlated with the strength of the compound's antagonistic action on 1alpha,25(OH)(2)D(3) signalling. Molecular dynamics simulations of the LBD complexed with the two antagonists could explain their different action by demonstrating a more drastic displacement of helix 12 through ZK168281 than through ZK159222. Moreover, the modelling could indicate a kink of helix 12 at amino acid residue G423, which provides the last four amino acid residues of helix 12 with a modulatory role for the partial agonism of some VDR antagonists, such as ZK159222. In conclusion, partial agonism of a VDR antagonist is lower the more it disturbs helix 12 in taking the optimal position for coactivator interaction.

Amino acids important for ligand specificity of the human constitutive androstane receptor

The human constitutive androstane receptor (CAR, NR1I3) is an important ligand-activated regulator of oxidative and conjugative enzymes and transport proteins. Because of the lack of a crystal structure of the ligand-binding domain (LBD), wide species differences in ligand specificity and the scarcity of well characterized ligands, the factors that determine CAR ligand specificity are not clear. To address this issue, we developed highly defined homology models of human CAR LBD to identify residues lining the ligand-binding pocket and to perform molecular dynamics simulations with known human CAR modulators. The roles of 22 LBD residues for basal activity, ligand selectivity, and interactions with co-regulators were studied using site-directed mutagenesis, mammalian co-transfection, and yeast two-hybrid assays. These studies identified several amino acids within helices 3 (Asn(165)), 5 (Val(199)), 11 (Tyr(326), Ile(330), and Gln(331)), and 12 (Leu(343) and Ile(346)) that contribute to the high basal activity of human CAR. Unique residues within helices 3 (Ile(164) and Asn(165)), 5 (Cys(202) and His(203)), and 7 (Phe(234) and Phe(238)) were found control the selectivity for CAR activators and inhibitors. A single residue in helix 7 (Phe(243)) appears to explain the human/mouse species difference in response of CAR to 17alpha-ethynyl-3,17beta-estradiol.

Birch PR-10c interacts with several biologically important ligands

PR-10c is a unique member of PR-10 proteins in birch, since it is the only one known to be post-translationally modified by glutathione and is not constitutively expressed in pollen. Both reduced and S-glutathiolated forms of PR-10c show low ribonuclease activity. However, the major function of the protein is apparently not yet resolved. Our protein-ligand interaction studies with saturation transfer difference (STD) NMR revealed that PR-10c interacts with several biologically important molecules, including cytokinin, flavonoid glycosides, sterols and emodin. Competition study with deoxycholate and kinetin revealed no statistically significant binding interference, indicating that these ligands have different binding sites in PR-10c. Ligand docking studies with a molecular model of PR-10c support the STD NMR results of ligand binding and binding epitopes, suggesting that there are three potential binding sites in PR-10c: two in the hydrophobic cavity and one in the glycine-rich loop. Our docking calculations suggested that only kinetin interacts with the glycine-rich loop, the binding occurring through its adenine moiety. Clear ligand specificity could be observed in the binding of nucleotide derivatives. S-glutathiolation of PR-10c did not affect kinetin binding. The present results suggest that birch PR-10c is a multifunctional protein, which has diverse roles in plant stress responses.

Vitamin D receptor agonists specifically modulate the volume of the ligand-binding pocket

Existing crystal structure data has indicated that 1alpha,25-dihydroxyvitamin D(3) (1alpha,25(OH)(2) D(3)) and its analogues bind the ligand-binding pocket (LBP) of the human vitamin D receptor in a very similar fashion. Because docking of a ligand into the LBP is a more flexible process than crystallography can monitor, we analyzed 1alpha,25(OH)(2)D(3), its 20-epi derivative MC1288, the two side-chain analogues Gemini and Ro43-83582 (a hexafluoro-derivative) by molecular dynamics simulations in a complex with the vitamin D receptor ligand-binding domain and a co-activator peptide. Superimposition of the structures showed that the side chain of MC1288, the first side chain of the conformation II of Gemini, the second side chain of Ro43-83582 in conformation I and the first side chain of Ro43-83582 in conformation II take the same agonistic position as the side chain of 1alpha,25(OH)(2)D(3). Compared with the LBP of the natural hormone MC1288 reduced the volume by 17%, and Gemini expanded it by 19%. The shrinking of the LBP of MC1288 and its expansion to accommodate the second side chain of Gemini or Ro43-83582 is the combined result of minor movements of more than 30 residues and major movements of a few critical amino acids. The agonist-selective recognition of anchoring OH groups by the conformational flexible residues Ala-303, Leu-309, and His-397 was confirmed by in vitro assays. In summary, variations in the volume of agonists lead to adaptations in the volume of the LBP and alternative contacts of anchoring OH-groups.

Gene Regulatory Potential of Nonsteroidal Vitamin D Receptor Ligands

The seco-steroid 1α,25-dihydroxyvitamin D3 [1α,25(OH)2D3] is a promising drug candidate due to its pleiotropic function including the regulation of calcium homeostasis, bone mineralization and cellular proliferation, differentiation, and apoptosis. We report here a novel class of nonsteroidal compounds, represented by the bis-aromatic molecules CD4409, CD4420, and CD4528, as ligands of the 1α,25(OH)2D3 receptor (VDR). Taking the known diphenylmethane derivative LG190178 as a reference, this study provides molecular evaluation of the interaction of nonsteroidal ligands with the VDR. All four nonsteroidal compounds were shown to induce VDR-retinoid X receptor heterodimer complex formation on a 1α,25(OH)2D3 response element, stabilize the agonistic conformation of the VDR ligand-binding domain, enable the interaction of VDR with coactivator proteins and contact with their three hydroxyl groups the same residues within the ligand-binding pocket of the VDR as 1α,25(OH)2D3. Molecular dynamics simulations demonstrated that all four nonsteroidal ligands take a shape within the ligand-binding pocket of the VDR that is very similar to that of the natural ligand. CD4528 is mimicking the natural hormone best and was found to be in vitro at least five times more potent than LG190178. In living cells, CD4528 was only two times less potent than 1α,25(OH)2D3 and induced mRNA expression of the VDR target gene CYP24 in a comparable fashion. At a noncalcemic dose of 150 μg/kg, CD4528 showed in vivo a clear induction of CYP24 expression and therefore may be used as a lead compound for the development of therapeutics against psoriasis, osteoporosis, and cancer.

Antagonist- and inverse agonist-driven interactions of the vitamin D receptor and the constitutive androstane receptor with corepressor protein

Ligand-dependent signal transduction by nuclear receptors (NRs) includes dynamic exchanges of coactivator (CoA) and corepressor (CoR) proteins. Here we focused on the structural determinants of the antagonist- and inverse agonist-enhanced interaction of the endocrine NR vitamin D receptor (VDR) and the adopted orphan NR constitutive androstane receptor (CAR) from two species with the CoR NR corepressor. We found that the pure VDR antagonist ZK168281 and the human CAR inverse agonist clotrimazole are both effective inhibitors of the CoA interaction of their respective receptors, whereas ZK168281 resembled more the mouse CAR inverse agonist androstanol in its ability to recruit CoR proteins. Molecular dynamics simulations resulted in comparable models for the CoR receptor interaction domain peptide bound to VDR/antagonist or CAR/inverse agonist complexes. A salt bridge between the CoR and a conserved lysine in helix 4 of the NR is central to this interaction, but also helix 12 was stabilized by direct contacts with residues of the CoR. Fixation of helix 12 in the antagonistic/inverse agonistic conformation prevents an energetically unfavorable free floatation of the C terminus. The comparable molecular mechanisms that explain the similar functional profile of antagonist and inverse agonists are likely to be extended from VDR and CAR to other members of the NR superfamily and may lead to the design of even more effective ligands.

Evaluation of a novel electronic eigenvalue (EEVA) molecular descriptor for QSAR/QSPR studies: validation using a benchmark steroid data set

A novel electronic eigenvalue (EEVA) descriptor of molecular structure for use in the derivation of predictive QSAR/QSPR models is described. Like other spectroscopic QSAR/QSPR descriptors, EEVA is also invariant as to the alignment of the structures concerned. Its performance was tested with respect to the CBG (corticosteroid binding globulin) affinity of 31 benchmark steroids. It appeared that the electronic structure of the steroids, i.e., the "spectra" derived from molecular orbital energies, is directly related to the CBG binding affinities. The predictive ability of EEVA is compared to other QSAR approaches, and its performance is discussed in the context of the Hammett equation. The good performance of EEVA is an indication of the essential quantum mechanical nature of QSAR. The EEVA method is a supplement to conventional 3D QSAR methods, which employ fields or surface properties derived from Coulombic and van der Waals interactions.

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.

A structural basis for the species-specific antagonism of 26,23-lactones on vitamin D signaling

The 26,23-lactone derivative of 1alpha,25-dihydroxyvitamin D3, TEI-9647, is a partial antagonist of the of human vitamin D receptor (VDR). However, we found that TEI-9647 in rat cells behaves as a weak VDR agonist. This behavior could be mimicked in human cells by the double mutagenesis of human VDR (specifically C403S and C410N). The increased agonistic action of TEI-9647 correlates to a gain in the interaction of the VDR with coactivator protein and a decreased stabilization of the antagonistic conformation of the receptor. Molecular dynamics simulations indicated that TEI-9647 acts as antagonist of human VDR by reducing the stability of helix 12 of the ligand binding domain. In contrast, N410 of the rat VDR stabilized, via backbone contacts, the interaction between helices 11 and 12. This results in TEI-9647 becoming a weak agonist in this organism.

Corepressor excess shifts the two-side chain vitamin D analog Gemini from an agonist to an inverse agonist of the vitamin D receptor

The vitamin D receptor (VDR) is an endocrine nuclear receptor that binds with high affinity its natural ligand 1alpha,25-dihydroxyvitamin D3. Gemini is a 1alpha,25-dihydroxyvitamin D3 analog with two identical side chains that, despite its significantly increased volume, binds to the VDR and can function as a potent agonist. This study demonstrates that, at excess corepressor (CoR) levels, Gemini shifts from an agonist to an inverse agonist that actively recruits CoR proteins to the VDR and mediates superrepression. Under these conditions Gemini stabilizes the VDR into a silent conformation, in which helix 12 of the ligand-binding domain is repositioned and thus unable to contribute to coactivator interaction. Amino acid F422 has been described as the lock of helix 12 and seems to be the most critical VDR residue in the inverse agonistic action of Gemini. Molecular dynamics simulations of the Gemini-VDR complex support these observation by indicating that the second side chain of Gemini induces tension to the receptor structure that can be released by a shift of helix 12. Taken together, Gemini is the first described (conditional) inverse agonist to an endocrine nuclear receptor and may function as a sensor for the cell-specific coactivator/CoR ratio.

The Effects of Sequence Variation on Genome-wide NRF2 Binding--New Target Genes and Regulatory SNPs

Transcription factor binding specificity is crucial for proper target gene regulation. Motif discovery algorithms identify the main features of the binding patterns, but the accuracy on the lower affinity sites is often poor. Nuclear factor E2-related factor 2 (NRF2) is a ubiquitous redox-activated transcription factor having a key protective role against endogenous and exogenous oxidant and electrophile stress. Herein, we decipher the effects of sequence variation on the DNA binding sequence of NRF2, in order to identify both genome-wide binding sites for NRF2 and disease-associated regulatory SNPs (rSNPs) with drastic effects on NRF2 binding. Interactions between NRF2 and DNA were studied using molecular modelling, and NRF2 chromatin immunoprecipitation-sequence datasets together with protein binding microarray measurements were utilized to study binding sequence variation in detail. The binding model thus generated was used to identify genome-wide binding sites for NRF2, and genomic binding sites with rSNPs that have strong effects on NRF2 binding and reside on active regulatory elements in human cells. As a proof of concept, miR-126–3p and -5p were identified as NRF2 target microRNAs, and a rSNP (rs113067944) residing on NRF2 target gene (Ferritin, light polypeptide, FTL) promoter was experimentally verified to decrease NRF2 binding and result in decreased transcriptional activity.