Acetylcholine receptor pore permeability studied by molecular dynamics simulation

A dynamic model of the closed-state pore of an acetylcholine receptor (five M2 α-helices stabilized with a (CH₂)₁₀₅ ring) is used to examine the migration of uncharged and charged probe particles equivalent to a hexahydrated sodium ion (van der Waals diameter 7.27 Å) propelled by varied external force along the channel axis. Ion movement through the pore is hindered by steric constraints and electrostatic interactions. The van der Waals gate is formed by helix residues 13′ (A-Val255, B-Val261, C-Val269, D-Val255, and E-Ile264), whereas the negatively charged residues in the upper part of the channel are important for ion selectivity.     

Structure and dynamics of the full-length M1 muscarinic acetylcholine receptor studied by molecular dynamics simulations

The three-dimensional structure of full-length structure of the M₁ muscarinic receptor was obtained through the fragmental homology modeling procedure. A 10-ns molecular dynamics (MD) simulation of the protein imbedded in a lipid slab and surrounded by water molecules was further used to relax the model. It was found that the homology model corresponded to the conformation in the ground state, since no significant motions of the backbone of transmembrane domains were observed. Furthermore, the reliability of the model was validated by analyzing key inter-helical contacts, sidechain-sidechain interactions, the formation of stable aromatic microdomains (clusters) and the docking of acetylcholine to its binding site. Moreover, a few conserved interactions observed in the X-ray structure of rhodopsin, such as inter-helical sidechain-sidechain hydrogen bonds were accurately reproduced in the MD simulation. The coupling of ACh to its binding site was found to be dominated by π-cation and salt bridge interactions, while its conformational space was restrained through van der Waals and hydrogen bond interactions. In general, such features were in very good agreement with the available experimental as well as with theoretical data. Considering the above, the structural information obtained in this study can be used a starting point to investigate the activation mechanism of the receptor and the ability to develop selective agonists and allosteric modulators which could be used for the treatment of Alzheimer’s disease.     

Constitutive activity and ligand-dependent activation of the nuclear receptor CAR-insights from molecular dynamics simulations

The constitutive androstane receptor (CAR) possesses, unlike most other nuclear receptors, a pronounced basal activity in vitro whose structural basis is still not fully understood. Using comparative molecular dynamics simulations of CAR X-ray crystal structures, we evaluated the molecular basis for constitutive activity and ligand-dependent receptor activation. Our results suggest that a combination of van der Waals interactions and hydrogen bonds is required to maintain the activation helix in the active conformation also in absence of a ligand. Furthermore, we identified conformational rearrangements within the ligand-binding pocket upon agonist binding and an influence of CAR inducers pregnanedione and CITCO on the helical conformation of the activation helix. Based on the results a model for ligand-dependent CAR activation is suggested.     

CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations

Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450(cam)(camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450(cam)(cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics.     

DNA condensation by TmHU studied by optical tweezers,AFM and molecular dynamics simulations

The compaction of DNA by the HU protein from Thermotoga maritima (TmHU) is analysed on a single-molecule level by the usage of an optical tweezers-assisted force clamp. The condensation reaction is investigated at forces between 2 and 40 pN applied to the ends of the DNA as well as in dependence on the TmHU concentration. At 2 and 5 pN, the DNA compaction down to 30% of the initial end-to-end distance takes place in two regimes. Increasing the force changes the progression of the reaction until almost nothing is observed at 40 pN. Based on the results of steered molecular dynamics simulations, the first regime of the length reduction is assigned to a primary level of DNA compaction by TmHU. The second one is supposed to correspond to the formation of higher levels of structural organisation. These findings are supported by results obtained by atomic force microscopy.     

Hsc70-Auxilin复合物的拉伸分子动力学模拟

Hsc70与auxilin蛋白组成的系统是Hsp70/Hsp40分子伴侣系统家族的一员,在热休克反应中发挥重要作用。本文为得出auxilin蛋白J结构域的关键氨基酸,首先采用由二硫键交联的Hsc70 R171C与auxilin D876C的复合物结晶结构作为初始模型,进行分子动力学模拟,通过比较平衡后的结合部位发现,将形成二硫键的氨基酸突变为原来的氨基酸结构在结合位点上与生化结果较为相近,之后利用此结构通过拉伸动力学模拟分析了auxilin蛋白J结构域与Hsc70的ATPase功能域的解离过程,并探讨了Hsc70与auxilin蛋白之间的相互作用力。结果表明位于HPD loop上的His874,Asp876,Thr879,螺旋Ⅲ上的Glu884,Asn895,Asp896,Ser899,Glu902,Asn903为关键氨基酸,这些数据符合之前核磁共振实验证实的T抗原J结构域的HPD基序和螺旋Ⅲ与Hsc70的ATPase功能域之间的相互作用。     

Structure-activity relationship studies on vitamin D lactam derivatives as vitamin D receptor antagonist

Structure–activity relationship studies on 1α,25-dihydroxyvitamin D₃-26,23-lactams (DLAMs), antagonists of vitamin D, were conducted, focusing on the substituents of the phenyl group. One of the derivatives (23S,25S)-DLAM-1P-3,5(OEt)₂, showed potent antagonistic activity with an IC₅₀ of 90 nM.     

Development of novel lithocholic acid derivatives as vitamin D receptor agonists

Lithocholic acid (2) was identified as the second endogenous ligand of vitamin D receptor (VDR), though its binding affinity to VDR and its vitamin D activity are very weak compared to those of the active metabolite of vitamin D₃, 1α,25-dihydroxyvitamin D₃ (1). 3-Acylated lithocholic acids were reported to be slightly more potent than lithocholic acid (2) as VDR agonists. Here, aiming to develop more potent lithocholic acid derivatives, we synthesized several derivatives bearing a 3-sulfonate/carbonate or 3-amino/amide substituent, and examined their differentiation-inducing activity toward human promyelocytic leukemia HL-60 cells. Introduction of a nitrogen atom at the 3-position of lithocholic acid (2) decreased the activity, but compound 6 bearing a 3-methylsulfonate group showed more potent activity than lithocholic acid (2) or its acylated derivatives. The binding of 6 to VDR was confirmed by competitive binding assay and X-ray crystallographic analysis of the complex of VDR ligand-binding domain (LBD) with 6.     

Atomistic molecular dynamics simulations of typical and atypical antipsychotic drugs at the dopamine D2 receptor (D2R) elucidates their inhibition mechanism

Dopamine D2 receptor (D2R) plays a pivotal role in nervous systems. Its dysfunction leads to the schizophrenia, Parkinson’s diseases and drug addiction. Since the crystal structure of the D2R was not solved yet, discovering of potent and highly selective anti-psychotic drugs carry challenges for different neurodegenerative diseases. In the current study, we modeled the three-dimensional (3D) structure of the D2R based on a recently crystallized structure of the dopamine D3 receptor. These two receptors share a high amino acid sequence homology (>70%). The interaction of the modeled receptor with well-known atypical and typical anti-psychotic drugs and the inhibition mechanisms of drugs at the catalytic domain were studied via atomistic molecular dynamics simulations. Our results revealed that, class-I and class-II forms of atypical and typical D2R antagonists follow different pathways in the inhibition of the D2Rs.     

Ligand-based pharmacophore modeling and docking studies on vitamin D receptor inhibitors

In recent years, pharmacophore modeling and molecular docking approaches have been extensively used to characterize the structural requirements and explore the conformational space of a ligand in the binding pocket of the selected target protein. Herein, we report a pharmacophore modeling and molecular docking of 45 compounds comprising of the indole scaffold as vitamin D receptor (VDR) inhibitors. Based on the selected best hypothesis (DRRRR.61), an atom-based three-dimensional quantitative structure-activity relationships model was developed to rationalize the structural requirement of biological activity modulating components. The developed model predicted the binding affinity for the training set and test set with R²_(training) = 0.8869 and R²_(test) = 0.8139, respectively. Furthermore, molecular docking and dynamics simulation were performed to understand the underpinning of binding interaction and stability of selected VDR inhibitors in the binding pocket. In conclusion, the results presented here, in the form of functional and structural data, agreed well with the proposed pharmacophores and provide further insights into the development of novel VDR inhibitors with better activity.     

Structural development of p-carborane-based potent non-secosteroidal vitamin D analogs

Non-secosteroidal vitamin D receptor (VDR) ligands are promising candidates for many clinical applications. We recently developed novel non-secosteroidal VDR agonists based on p-carborane (an icosahedral carbon-containing boron cluster) as a hydrophobic core structure. Here, we report the design, synthesis and biological evaluation of carborane-based vitamin D analogs bearing various substituents at the diol moiety. Among the synthesized compounds, methylene derivative 31 exhibited the most potent vitamin D activity, which was comparable to that of the natural hormone, 1α,25(OH)₂D₃. This compound is one of the most potent non-secosteroidal VDR agonists reported to date, and is a promising lead for development of novel drug candidates.     

Premature aging in vitamin D receptor mutant mice

Hypervitaminosis vitamin D₃ has been recently implicated in premature aging through the regulation of 1alpha hydroxylase expression by klotho and fibroblast growth factor-23 (Fgf-23). Here we examined whether the lack of hormonal function of vitamin D₃ in mice is linked to aging phenomena. For this, we used vitamin D₃ receptor (VDR) “Tokyo” knockout (KO) mice (fed with a special rescue diet) and analyzed their growth, skin and cerebellar morphology, as well as overall motor performance. We also studied the expression of aging-related genes, such as Fgf-23, nuclear factor kappaB (NF-kappaB), p53, insulin like growth factor 1 (IGF1) and IGF1 receptor (IGF1R), in liver, as well as klotho in liver, kidney and prostate tissues. Overall, VDR KO mice showed several aging related phenotypes, including poorer survival, early alopecia, thickened skin, enlarged sebaceous glands and development of epidermal cysts. There was no difference either in the structure of cerebellum or in the number of Purkinje cells. Unlike the wildtype controls, VDR KO mice lose their ability to swim after 6 months of age. Expression of all the genes was lower in old VDR KO mice, but only NF-kappaB, Fgf-23, p53 and IGF1R were significantly lower. Since the phenotype of aged VDR knockout mice is similar to mouse models with hypervitaminosis D₃, our study suggests that VDR genetic ablation promotes premature aging in mice, and that vitamin D₃ homeostasis regulates physiological aging.     

Phosphorylation of the arginine/serine repeats of lamin B receptor by SRPK1—Insights from molecular dynamics simulations

Background

Arginine/serine (RS) repeats are found in several proteins in metazoans with a wide variety of functions, many of which are regulated by SR protein kinase 1 (SRPK1)-mediated phosphorylation. Lamin B receptor (LBR) is such a protein implicated in chromatin anchorage to the nuclear envelope.

Methods

Molecular dynamics simulations were used to investigate the conformation of two LBR peptides containing four (human-) and five (turkey-orthologue) consecutive RS dipeptides, in their unphosphorylated and phosphorylated forms and of a conserved peptide, in isolation and in complex with SRPK1. GST pull-down assays were employed to study LBR interactions.

Results

Unphosphorylated RS repeats adopt short, transient helical conformations, whereas serine phosphorylation induces Arginine-claw-like structures. The SRSRSRSPGR peptide, overlapping with the LBR RS repeats, docks into the known, acidic docking groove of SRPK1, in an extended conformation. Phosphorylation by SRPK1 is necessary for the association of LBR with histone H3.

Conclusions

The C-terminal region of the LBR RS domain constitutes a recognition platform for SRPK1, which uses the same recognition mechanism for LBR as for substrates with long RS domains. This docking may promote unfolding of the RS repeats destined to be phosphorylated. Phosphorylation induces Arginine-claw-like conformations, irrespective of the RS-repeat length, that may facilitate interactions with basic partners.

General significance

Our results shed light on the conformational preferences of an important class of repeats before and after their phosphorylation and support the idea that even short RS domains may be constituents of recognition platforms for SRPK1, thus adding to knowledge towards a full understanding of their phosphorylation mechanism.     

Ribonucleotide reductase inhibition by p-alkoxyphenols studied by molecular docking and molecular dynamics simulations

Ribonucleotide reductase (RNR) is necessary for production of the precursor deoxyribonucleotides for DNA synthesis. Class Ia RNR functions via a stable free radical in one of the two components protein R2. The enzyme mechanism involves long range (proton coupled) electron transfer between protein R1 and the tyrosyl radical in protein R2. Earlier experimental studies showed that p-alkoxyphenols inhibit RNR. Here, molecular docking and molecular dynamics simulations involving protein R2 suggest an inhibition mechanism for p-alkoxyphenols . A low energy binding pocket is identified in protein R2. The preferred configuration provides a structural basis explaining their specific binding to the Escherichia coli and mouse R2 proteins. Trp48 (E. coli numbering), on the electron transfer pathway, is involved in the interactions with the inhibitors. The relative order of the binding energies calculated for the phenol derivatives to protein R2 is correlated with earlier experimental data on inhibition efficiency, in turn related to increasing size of the hydrophobic alkyl substituents. Using the configuration identified by molecular docking as a starting point for molecular dynamics simulations, we find that the p-allyloxyphenol interrupts the catalytic electron transfer pathway of the R2 protein by forming hydrogen bonds with Trp48 and Asp237, thus explaining the inhibitory activity of p-alkoxyphenols.