GPCRmd uncovers the dynamics of the 3D-GPCRome.

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.

Publisher Correction: GPCRmd uncovers the dynamics of the 3D-GPCRome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D3 dopamine receptor (D3R) in complex with an antagonist. Close to 60 µs of simulation time was generated and the resulting MD refined models were compared to a D3R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.

Molecular dynamics simulations and NMR spectroscopy studies of trehalose-lipid bilayer systems.

The disaccharide trehalose (TRH) strongly affects the physical properties of lipid bilayers. We investigate interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH using NMR spectroscopy and molecular dynamics (MD) computer simulations. We compare dipolar couplings derived from DMPC/TRH trajectories with those determined (i) experimentally in TRH using conventional high-resolution NMR in a weakly ordered solvent (bicelles), and (ii) by solid-state NMR in multilamellar vesicles (MLV) formed by DMPC. Analysis of the experimental and MD-derived couplings in DMPC indicated that the force field used in the simulations reasonably well describes the experimental results with the exception for the glycerol fragment that exhibits significant deviations. The signs of dipolar couplings, not available from the experiments on highly ordered systems, were determined from the trajectory analysis. The crucial step in the analysis of residual dipolar couplings (RDCs) in TRH determined in a bicelle-environment was access to the conformational distributions derived from the MD trajectory. Furthermore, the conformational behavior of TRH, investigated by J-couplings, in the ordered and isotropic phases is essentially identical, indicating that the general assumptions in the analyses of RDCs are well founded.

Mapping the Interface of a GPCR Dimer: A Structural Model of the A2A Adenosine and D2 Dopamine Receptor Heteromer.

The A2A adenosine (A2AR) and D2 dopamine (D2R) receptors form oligomers in the cell membrane and allosteric interactions across the A2AR-D2R heteromer represent a target for development of drugs against central nervous system disorders. However, understanding of the molecular determinants of A2AR-D2R heteromerization and the allosteric antagonistic interactions between the receptor protomers is still limited. In this work, a structural model of the A2AR-D2R heterodimer was generated using a combined experimental and computational approach. Regions involved in the heteromer interface were modeled based on the effects of peptides derived from the transmembrane (TM) helices on A2AR-D2R receptor-receptor interactions in bioluminescence resonance energy transfer (BRET) and proximity ligation assays. Peptides corresponding to TM-IV and TM-V of the A2AR blocked heterodimer interactions and disrupted the allosteric effect of A2AR activation on D2R agonist binding. Protein-protein docking was used to construct a model of the A2AR-D2R heterodimer with a TM-IV/V interface, which was refined using molecular dynamics simulations. Mutations in the predicted interface reduced A2AR-D2R interactions in BRET experiments and altered the allosteric modulation. The heterodimer model provided insights into the structural basis of allosteric modulation and the technique developed to characterize the A2AR-D2R interface can be extended to study the many other G protein-coupled receptors that engage in heteroreceptor complexes.

Coarse-Grained Molecular Dynamics Simulations of Membrane-Trehalose Interactions.

It is well established that trehalose (TRH) affects the physical properties of lipid bilayers and stabilizes biological membranes. We present molecular dynamics (MD) computer simulations to investigate the interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH. Both atomistic and coarse-grained (CG) interaction models were employed, and the coarse graining of DMPC leads to a reduction in the acyl chain length corresponding to a 1,2-dilauroyl-sn-glycero-3-phosphocholine lipid (DLPC). Several modifications of the Martini interaction model, used for CG simulations, were implemented, resulting in different potentials of mean force (PMFs) for DMPC bilayer-TRH interactions. These PMFs were subsequently used in a simple two-site analytical model for the description of sugar binding at the membrane interface. In contrast to that in atomistic MD simulations, the binding in the CG model was not in agreement with the two-site model. Our interpretation is that the interaction balance, involving water, TRH, and lipids, in the CG systems needs further tuning of the force-field parameters. The area per lipid is only weakly affected by TRH concentration, whereas the compressibility modulus related to the fluctuations of the membrane increases with an increase in TRH content. In agreement with experimental findings, the bending modulus is not affected by the inclusion of TRH. The important aspects of lipid bilayer interactions with biomolecules are membrane curvature generation and sensing. In the present investigation, membrane curvature is generated by artificial buckling of the bilayer in one dimension. It turns out that TRH prefers the regions with the highest curvature, which enables the most favorable situation for lipid-sugar interactions.

Molecular dynamics simulations of membrane-sugar interactions.

It is well documented that disaccharides in general and trehalose (TRH) in particular strongly affect physical properties and functionality of lipid bilayers. We investigate interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH by means of molecular dynamics (MD) computer simulations. Ten different TRH concentrations were studied in the range wTRH = 0-0.20 (w/w). The potential of mean force (PMF) for DMPC bilayer-TRH interactions was determined using two different force fields, and was subsequently used in a simple analytical model for description of sugar binding at the membrane interface. The MD results were in good agreement with the predictions of the model. The net affinities of TRH for the DMPC bilayer derived from the model and MD simulations were compared with experimental results. The area per lipid increases and the membrane becomes thinner with increased TRH concentration, which is interpreted as an intercalation effect of the TRH molecules into the polar part of the lipids, resulting in conformational changes in the chains. These results are consistent with recent experimental observations. The compressibility modulus related to the fluctuations of the membrane increases dramatically with increased TRH concentration, which indicates higher order and rigidity of the bilayer. This is also reflected in a decrease (by a factor of 15) of the lateral diffusion of the lipids. We interpret these observations as a formation of a glassy state at the interface of the membrane, which has been suggested in the literature as a hypothesis for the membrane-sugar interactions.

Molecular dynamics simulations of membranes composed of glycolipids and phospholipids.

Lipid membranes composed of 1,2-di-(9Z,12Z,15Z)-octadecatrienoyl-3-O-ß-D-galactosyl-sn-glycerol or monogalactosyldiacylglycerol (MGDG) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) were studied by means of molecular dynamics (MD) computer simulations. Three lipid compositions were considered: 0%, 20%, and 45% MGDG (by mole) denoted as MG-0, MG-20, and MG-45, respectively. The article is focused on the calculation of NMR dipolar interactions, which were confronted with previously reported experimental couplings. Dynamical processes and orientational distributions relevant for the averaging of dipolar interactions were evaluated. Furthermore, several parameters important for characterization of the bilayer structure, molecular organization, and dynamics were investigated. In general, only a minor change in DMPC properties was observed upon the increased MGDG/DMPC ratio, whereas properties related to MGDG undergo a more pronounced change. This effect was ascribed to the fact that DMPC is a bilayer (L(α)) forming lipid, whereas MGDG prefers a reverse hexagonal (H(II)) arrangement.