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 transmembrane helix region of the models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode and the second extracellular loop region. 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

Conformational transitions of the serotonin 5-HT3 receptor

International audienceThe serotonin 5-HT3 receptor is a pentameric ligand-gated ion channel (pLGIC). It belongs to a large family of receptors that function as allosteric signal transducers across the plasma membrane1,2; upon binding of neurotransmitter molecules to extracellular sites, the receptors undergo complex conformational transitions that result in transient opening of a pore permeable to ions. 5-HT3 receptors are therapeutic targets for emesis and nausea, irritable bowel syndrome and depression3. In spite of several reported pLGIC structures4-8, no clear unifying view has emerged on the conformational transitions involved in channel gating. Here we report four cryo-electron microscopy structures of the full-length mouse 5-HT3 receptor in complex with the anti-emetic drug tropisetron, with serotonin, and with serotonin and a positive allosteric modulator, at resolutions ranging from 3.2 Å to 4.5 Å. The tropisetron-bound structure resembles those obtained with an inhibitory nanobody5 or without ligand9. The other structures include an 'open' state and two ligand-bound states. We present computational insights into the dynamics of the structures, their pore hydration and free-energy profiles, and characterize movements at the gate level and cation accessibility in the pore. Together, these data deepen our understanding of the gating mechanism of pLGICs and capture ligand binding in unprecedented detail

Combining Structural Data with Computational Methodologies to Investigate Structure–Function Relationships in TRP Channels

Since the emergence of high-resolution three-dimensional structures of membrane proteins, and the increasing availability of state-of-the-art algorithms and high-performance-computing facilities, classical molecular dynamics (MD) simulations have become a routine device to explore the molecular behavior of these proteins. The rise of cryo-electron microscopy (cryo-EM) as a credible experimental tool to resolve structures at an atomic level has revolutionized structural biology in recent years, culminating in the disclosure of the first high-resolution three-dimensional structure of a transient receptor potential (TRP) channel, the vanilloid receptor 1 (TRPV1). As a result, the number of research articles investigating the molecular behavior of TRP channels using macromolecular simulation techniques has proliferated. This review provides an overview of the current state of this field, including our understanding of TRP channel structure, the framework of classical MD simulations, and how to perform such simulations to investigate structure-function relationships in TRP channels