The G protein-coupled receptors (GPCRs) are the largest family in the mammalian genome, and are critical to a number of cell signaling processes. As a result, they are of enormous biomedical importance;by some estimates, as many as 50% of new pharmaceuticals target GPCRs. Unsurprisingly, there has been a huge research investment in understanding their biophysics. However, integral membrane proteins are challenging to work with experimentally, leaving an opportunity for computational methods to make a significant contribution. We will use multiscale modeling techniques, including all-atom molecular dynamics simulations and elastic network models, to explore the behavior of several GPCRs, including rhodopsin (and its retinal-free form, opsin) and the ?2-adrenergic receptor (B2AR). Specifically, we will investigate the role of ligand binding in modulating GPCR function, via two separate all-atom molecular dynamics calculations. Microsecond-scale simulations of opsin will, when contrasted with our previous work on rhodopsin in the dark state and during the early stages of activation, allow us to see which interactions in rhodopsin are determined by the presence of the ligand, while the planned simulations of the full activation process will give the first atomic-level view of the structural changes involved in GPCR activation;this knowledge could be critical to the design of novel inhibitors to other GPCRs. The second goal of this proposal is to clarify the role of internal waters in the activation mechanism of GPCRs;our previous simulations described significant increases in the internal hydration of rhodopsin and B2AR. Here, we propose to pursue those observations more rigorously, using automatic pattern recognition methods to correlate hydration changes with functionally interesting protein motions in simulations of rhodopsin, B2AR, and the cannabinoid-2 receptor (CB2). The third goal of the proposal is to develop elastic network models - a simple, computationally inexpensive approach where the protein's interactions are represented as a network of springs - in order to explore larger scale problems not readily amenable to all-atom molecular dynamics, like the modulation of protein motions by G protein binding and GPCR oligomerization. A number of possible network model implementations will be considered, and the models will be carefully validated by quantitative comparison to extensive molecular dynamics simulations, including those proposed for the first aim. The fourth and final goal of the proposal is to assess the validity of a common assumption, that rhodopsin is a good template for understanding GPCR activation in general. To test this hypothesis we will apply multiple computational methods, including long timescale molecular dynamics and elastic network models, to a series of GPCRs, including rhodopsin, opsin, B2AR, and CB2. We will quantitatively correlate the fluctuations of the different GPCRs, with the hypothesis that motions conserved across multiple GPCRs are likely to be functionally significant.

Public Health Relevance

G protein-coupled receptors are the largest family of proteins in the human genome, and are the most targeted proteins for therapeutics development. We will use multiscale modeling methods to connect their structure and dynamics to function, in order to aid in the design and refinement of novel drugs.

National Institute of Health (NIH)
National Institute of General Medical Sciences (NIGMS)
Research Project (R01)
Project #
Application #
Study Section
Biophysics of Neural Systems Study Section (BPNS)
Program Officer
Lyster, Peter
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
University of Rochester
School of Medicine & Dentistry
United States
Zip Code
Salas-Estrada, Leslie A; Leioatts, Nicholas; Romo, Tod D et al. (2018) Lipids Alter Rhodopsin Function via Ligand-like and Solvent-like Interactions. Biophys J 114:355-367
Sur, Sreyoshi; Romo, Tod D; Grossfield, Alan (2018) Selectivity and Mechanism of Fengycin, an Antimicrobial Lipopeptide, from Molecular Dynamics. J Phys Chem B 122:2219-2226
Aytenfisu, Asaminew H; Spasic, Aleksandar; Grossfield, Alan et al. (2017) Revised RNA Dihedral Parameters for the Amber Force Field Improve RNA Molecular Dynamics. J Chem Theory Comput 13:900-915
Lin, Dejun (2015) Generalized and efficient algorithm for computing multipole energies and gradients based on Cartesian tensors. J Chem Phys 143:114115
Lin, Dejun; Grossfield, Alan (2015) Thermodynamics of Micelle Formation and Membrane Fusion Modulate Antimicrobial Lipopeptide Activity. Biophys J 109:750-9
Leioatts, Nicholas; Romo, Tod D; Danial, Shairy Azmy et al. (2015) Retinal Conformation Changes Rhodopsin's Dynamic Ensemble. Biophys J 109:608-17
Romo, Tod D; Grossfield, Alan (2014) How fast is your camera? Timescales for molecular motion and their role in restraining molecular dynamics. Biophys J 106:2549-51
Leioatts, Nicholas; Suresh, Pooja; Romo, Tod D et al. (2014) Structure-based simulations reveal concerted dynamics of GPCR activation. Proteins 82:2538-51
Kimura, Tomohiro; Vukoti, Krishna; Lynch, Diane L et al. (2014) Global fold of human cannabinoid type 2 receptor probed by solid-state 13C-, 15N-MAS NMR and molecular dynamics simulations. Proteins 82:452-65
Leioatts, Nicholas; Mertz, Blake; Martínez-Mayorga, Karina et al. (2014) Retinal ligand mobility explains internal hydration and reconciles active rhodopsin structures. Biochemistry 53:376-85

Showing the most recent 10 out of 23 publications