Molecular dynamics (MD) simulations of membrane-embedded G-protein coupled receptors (GPCRs) have rapidly gained popularity among the molecular simulation community in recent years, a trend which has an obvious link to the tremendous pharmaceutical importance of this group of receptors and the increasing availability of crystal structures. In view of the widespread use of this technique, it is of fundamental importance to ensure the reliability and robustness of the methodologies so they yield valid results and enable sufficiently accurate predictions to be made. In this work, 200 ns simulations of the A2a adenosine receptor (A2a AR) have been produced and evaluated in the light of these requirements. The conformational dynamics of the target protein, as obtained from replicate simulations in both the presence and absence of an inverse agonist ligand (ZM241385), have been investigated and compared using principal component analysis (PCA). Results show that, on this time scale, convergence of the replicates is not readily evident and dependent on the types of the protein motions considered. Thus rates of inter- as opposed to intrahelical relaxation and sampling can be different. When studied individually, we find that helices III and IV have noticeably greater stability than helices I, II, V, VI, and VII in the apo form. The addition of the inverse agonist ligand greatly improves the stability of all helices.
A new subfamily of glycosyl hydrolase family GH13 was recently proposed for α-amylases from Anoxybacillus species (ASKA and ADTA), Geobacillus thermoleovorans (GTA, Pizzo, and GtamyII), Bacillus aquimaris (BaqA), and 95 other putative protein homologues. To understand this new GH13 subfamily, we report crystal structures of truncated ASKA (TASKA). ASKA is a thermostable enzyme capable of producing high levels of maltose. Unlike GTA, biochemical analysis showed that Ca(2+) ion supplementation enhances the catalytic activities of ASKA and TASKA. The crystal structures reveal the presence of four Ca(2+) ion binding sites, with three of these binding sites are highly conserved among Anoxybacillus α-amylases. This work provides structural insights into this new GH13 subfamily both in the apo form and in complex with maltose. Furthermore, structural comparison of TASKA and GTA provides an overview of the conformational changes accompanying maltose binding at each subsite.
Crystal structures of holo vitamin D receptor (VDR) revealed a canonical conformation in which the ligand is entrapped in a hydrophobic cavity buried in the ligand-binding domain (LBD). The mousetrap model postulates that helix 12 is positioned away from the domain to expose the interior cavity. However, the extended form of helix 12 is likely due to artifacts during crystallization. In this study, we set out to investigate conformational dynamics of apo VDR using molecular dynamics simulation on microsecond timescale. Here we show the neighboring backbones of helix 2-helix 3n and beta strand 2-helix 6 of LBD, instead of the helix 12, undergo large-scale motion, possibly gating the entrance of ligand to the ligand binding domain. Docking analysis to the simulated open structure of VDR with the estimated free energy of -37.0 kJ/mol, would emphasise the role of H2-H3n and S2-H6 in facilitating the entrance of calcitriol to the LBD of VDR.