Molecular modeling results indicate that the VDR contains two overlapping ligand binding pockets (LBP). Differential ligand stability and fractional occupancy of the two LBP has been physiochemically linked to the regulation of VDR-dependent genomic and non-genomic cellular responses. The purpose of this report is to develop an unbiased molecular modeling protocol that serves as a good starting point in simulating the dynamic interaction between 1α,25(OH) 2 -vitamin D 3 (1,25D3) and the VDR LBP. To accomplish this goal, the flexible docking protocol developed allowed for flexibility in the VDR ligand and the VDR atoms that form the surfaces of the VDR LBP. This approach blindly replicated the 1,25D3 conformation and side-chain dynamics observed in the VDR x-ray structure. The results are also consistent with the previously published tenants of the vitamin D sterol (VDS)-VDR conformational ensemble model. Furthermore, we used flexible docking in combination with whole cell patch clamp electrophysiology and steroid competition assays to demonstrate that a) new non-vitamin D VDR ligands show a different pocket selectivity when compared to 1,25D3 that is qualitatively consistent with their ability to stimulate chloride channels and b) a new route of ligand binding provides a novel hypothesis describing the structural nuances that underlie hypercalceamia.

A series of monoprotonated aliphatic diamines has been examined, which crystallize in three general motifs: salt-bridged, cyclic, or clustered. The monoprotonated triflic acid salt of Me2N(CH2)4NMe2 forms a proton-bridged cyclic cation. The internal N-N distance is 2.66 A, with the bridging proton in the middle, having an NHN angle >/=172 degrees. The triflate oxygens lie more than 4 A away from the midpoint between the nitrogen atoms, indicating that a salt bridge does not form. The average NH distance in a solid sample was determined by measuring the 15N-H dipolar coupling in the triflic acid salt of the completely deuterated diamine (CD3)2N(CD2)4N(CD3)2. The value of the dipolar coupling constant, 5250 +/- 90 Hz, corresponds to an average NH distance of 1.32 A, nearly half-the NN distance. That result agrees with DFT calculations, which give a double-well potential minimum for proton transit between the two amino groups, having a zero-point vibrational level close to the barrier top. Theory predicts that the maximum value of the zero point vibrational wave function is almost coincident with a local potential energy maximum, consistent with the experimental findings.

Vibrational spectra of the conjugate acid of Me(2)NCH(2)CH(2)CH(2)CH(2)NMe(2) (N,N,N',N'-tetramethylputrescine) have been examined in the gaseous and crystalline phases using Infrared Multiple Photon Dissociation (IRMPD) spectroscopy, Inelastic Neutron Scattering (INS), and high pressure Raman spectroscopy. A band observed near 530 cm(-1) is assigned to the asymmetric stretch of the bridging proton between the two nitrogens, based on deuterium substitution and pressure dependence. The NN distance measured by X-ray crystallography gives a good match to DFT calculations, and the experimental band position agrees with the value predicted from theory using a 2-dimensional potential energy surface. The reduced dimensionality potential energy surface, which treats the ion as though it possesses a linear NHN geometry, shows low barriers to proton transit from one nitrogen to the other, with zero point levels close to the barrier tops. In contrast, two other related systems containing strong hydrogen bonds do not exhibit the same spectroscopic signature of a low barrier hydrogen bond (LBHB). On the one hand, the IRMPD spectra of the conjugate acid ions of the amino acid N,N,N',N'-tetramethylornithine (in which the two nitrogens have different basicities) show fewer bands and no comparable isotopic shifts in the low frequency domain. On the other hand, the IRMPD spectrum of the shorter homologue Me(2)NCH(2)CH(2)CH(2)NMe(2) (N,N,N',N'-tetramethyl-1,3-propanediamine), for which the NHN bond angle deviates substantially from linearity, displays more than one band in the 1100-1400 cm(-1) domain, which vanish as a consequence of deuteration.

PurposeThe thermodynamically favored complex between the nuclear vitamin D receptor (VDR) and 1α,25(OH)2-vitamin D3 (1,25D3) triggers a shift in equilibrium to favor VDR binding to DNA, heterodimerization with the nuclear retinoid x receptor (RXR) and subsequent regulation of gene transcription. The key amino acids and structural requirements governing VDR binding to nuclear coactivators (NCoA) are well defined. Yet very little is understood about the internal changes in amino acid flexibility underpinning the control of ligand affinity, helix 12 conformation and function. Herein, we use molecular dynamics (MD) to study how the backbone and side-chain flexibility of the VDR differs when a) complexed to 1α,25(OH)2-vitamin D3 (1,25D3, agonist) and (23S),25-dehydro-1α(OH)-vitamin D3-26,23-lactone (MK, antagonist); b) residues that form hydrogen bonds with the C25-OH (H305 and H397) of 1,25D3 are mutated to phenylalanine; c) helix 12 conformation is changed and ligand is removed; and d) x-ray water near the C1- and C3-OH groups of 1,25D3 are present or replaced with explicit solvent.MethodsWe performed molecular dynamic simulations on the apo- and holo-VDRs and used T-Analyst to monitor the changes in the backbone and side-chain flexibility of residues that form regions of the VDR ligand binding pocket (LBP), NCoA surface and control helix 12 conformation.ResultsThe VDR-1,25D3 and VDR-MK MD simulations demonstrate that 1,25D3 and MK induce highly similar changes in backbone and side-chain flexibility in residues that form the LBP. MK however did increase the backbone and side-chain flexibility of L404 and R274 respectively. MK also induced expansion of the VDR charge clamp (i.e. NCoA surface) and weakened the intramolecular interaction between H305---V418 (helix 12) and TYR401 (helix 11). In VDR_FF, MK induced a generally more rigid LBP and stronger interaction between F397 and F422 than 1,25D3, and reduced the flexibility of the R274 side-chain. Lastly the VDR MD simulations indicate that R274 can sample multiple conformations in the presence of ligand. When the R274 is extended, the β-OH group of 1,25D3 lies proximal to the backbone carbonyl oxygen of R274 and the side-chain forms H-bonds with hinge domain residues. This differs from the x-ray, kinked geometry, where the side-chain forms an H-bond with the 1α-OH group. Furthermore, 1,25D3, but not MK was observed to stabilize the x-ray geometry of R274 during the > 30 ns MD runs.ConclusionsThe MD methodology applied herein provides an in silico foundation to be expanded upon to better understand the intrinsic flexibility of the VDR and better understand key side-chain and backbone movements involved in the bimolecular interaction between the VDR and its’ ligands.Electronic supplementary materialThe online version of this article (doi:10.1186/2193-9616-1-2) contains supplementary material, which is available to authorized users.