Alan Grossfield

A Lipid Pathway for Ligand Binding Is Necessary for a Cannabinoid G Protein-coupled Receptor

Recent isothiocyanate covalent labeling studies have suggested that a classical cannabinoid, (−)-7′-isothiocyanato-11-hydroxy-1′,1′dimethylheptyl-hexahydrocannabinol (AM841), enters the cannabinoid CB2 receptor via the lipid bilayer (Pei, Y., Mercier, R. W., Anday, J. K., Thakur, G. A., Zvonok, A. M., Hurst, D., Reggio, P. H., Janero, D. R., and Makriyannis, A. (2008) Chem. Biol. 15, 1207–1219). However, the sequence of steps involved in such a lipid pathway entry has not yet been elucidated. Here, we test the hypothesis that the endogenous cannabinoid sn-2-arachidonoylglycerol (2-AG) attains access to the CB2 receptor via the lipid bilayer. To this end, we have employed microsecond time scale all-atom molecular dynamics (MD) simulations of the interaction of 2-AG with CB2 via a palmitoyl-oleoyl-phosphatidylcholine lipid bilayer. Results suggest the following: 1) 2-AG first partitions out of bulk lipid at the transmembrane α-helix (TMH) 6/7 interface; 2) 2-AG then enters the CB2 receptor binding pocket by passing between TMH6 and TMH7; 3) the entrance of the 2-AG headgroup into the CB2 binding pocket is sufficient to trigger breaking of the intracellular TMH3/6 ionic lock and the movement of the TMH6 intracellular end away from TMH3; and 4) subsequent to protonation at D3.49/D6.30, further 2-AG entry into the ligand binding pocket results in both a W6.48 toggle switch change and a large influx of water. To our knowledge, this is the first demonstration via unbiased molecular dynamics that a ligand can access the binding pocket of a class A G protein-coupled receptor via the lipid bilayer and the first demonstration via molecular dynamics of G protein-coupled receptor activation triggered by a ligand binding event.

Convergence of molecular dynamics simulations of membrane proteins.

The central question in evaluating almost any result from a molecular dynamics simulation is whether the calculation has converged. Unfortunately, assessing the ergodicity of a single trajectory is very difficult to do. In this work, we assess the sampling of molecular dynamics simulations of the membrane protein rhodopsin by comparing the results from 26 independent trajectories, each run for 100 ns. By examining principal components and cluster populations, we show that rhodopsins fluctuations are not well described by 100 ns of dynamics, and that the sampling is not fully converged even for individual loops. The results serve as a reminder of the caution required when interpreting molecular dynamics simulations of macromolecules. Proteins 2007.

Internal hydration increases during activation of the G-protein-coupled receptor rhodopsin.

Rhodopsin, the membrane protein responsible for dim-light vision, until recently was the only G-protein-coupled receptor (GPCR) with a known crystal structure. As a result, there is enormous interest in studying its structure, dynamics, and function. Here we report the results of three all-atom molecular dynamics simulations, each at least 1.5 micros, which predict that substantial changes in internal hydration play a functional role in rhodopsin activation. We confirm with (1)H magic angle spinning NMR that the increased hydration is specific to the metarhodopsin-I intermediate. The internal water molecules interact with several conserved residues, suggesting that changes in internal hydration may be important during the activation of other GPCRs. The results serve to illustrate the synergism of long-time-scale molecular dynamics simulations and NMR in enhancing our understanding of GPCR function.

Structural and dynamic effects of cholesterol at preferred sites of interaction with rhodopsin identified from microsecond length molecular dynamics simulations.

An unresolved question about GPCR function is the role of membrane components in receptor stability and activation. In particular, cholesterol is known to affect the function of membrane proteins, but the details of its effect on GPCRs are still elusive. Here, we describe how cholesterol modulates the behavior of the TM1‐TM2‐TM7‐helix 8(H8) functional network that comprises the highly conserved NPxxY(x)5,6F motif, through specific interactions with the receptor. The inferences are based on the analysis of microsecond length molecular dynamics (MD) simulations of rhodopsin in an explicit membrane environment. Three regions on the rhodopsin exhibit the highest cholesterol density throughout the trajectory: the extracellular end of TM7, a location resembling the high‐density sterol area from the electron microscopy data; the intracellular parts of TM1, TM2, and TM4, a region suggested as the cholesterol binding site in the recent X‐ray crystallography data on β2‐adrenergic GPCR; and the intracellular ends of TM2‐TM3, a location that was categorized as the high cholesterol density area in multiple independent 100 ns MD simulations of the same system. We found that cholesterol primarily affects specific local perturbations of the helical TM domains such as the kinks in TM1, TM2, and TM7. These local distortions, in turn, relate to rigid‐body motions of the TMs in the TM1‐TM2‐TM7‐H8 bundle. The specificity of the effects stems from the nonuniform distribution of cholesterol around the protein. Through correlation analysis we connect local effects of cholesterol on structural perturbations with a regulatory role of cholesterol in the structural rearrangements involved in GPCR function. Proteins 2009.

Dependence of ion hydration on the sign of the ion’s charge

The solvation of simple ions in water is studied using molecular dynamics simulations with a polarizable force field. Previous simulations using this potential demonstrated that anions are more favorably solvated in water than cations. The present work is an attempt to explain this result by examining the effects of ions on the surrounding water structure, with particular focus on the first solvation shell and its interactions with the surrounding water. We conclude that while the first solvation shell surrounding cations is frustrated by competition between ion-water and water-water interactions, solvation of anions is compatible with good water-water interactions.

Recent progress in the study of G protein-coupled receptors with molecular dynamics computer simulations

G protein-coupled receptors (GPCRs) are a large, biomedically important family of proteins, and the recent explosion of new high-resolution structural information about them has provided an enormous opportunity for computational modeling to make major contributions. In particular, molecular dynamics simulations have become a driving factor in many areas of GPCR biophysics, improving our understanding of lipid-protein interaction, activation mechanisms, and internal hydration. Given that computers will continue to get faster and more structures will be solved, the importance of computational methods will only continue to grow, particularly as simulation research is more closely coupled to experiment.

Concerted Interconversion between Ionic Lock Substates of the β2 Adrenergic Receptor Revealed by Microsecond Timescale Molecular Dynamics

The recently solved crystallographic structures for the A(2A) adenosine receptor and the beta(1) and beta(2) adrenergic receptors have shown important differences between members of the class-A G-protein-coupled receptors and their archetypal model, rhodopsin, such as the apparent breaking of the ionic lock that stabilizes the inactive structure. Here, we characterize a 1.02 mus all-atom simulation of an apo-beta(2) adrenergic receptor that is missing the third intracellular loop to better understand the inactive structure. Although we find that the structure is remarkably rigid, there is a rapid influx of water into the core of the protein, as well as a slight expansion of the molecule relative to the crystal structure. In contrast to the x-ray crystal structures, the ionic lock rapidly reforms, although we see an activation-precursor-like event wherein the ionic lock opens for approximately 200 ns, accompanied by movements in the transmembrane helices associated with activation. When the lock reforms, we see the structure return to its inactive conformation. We also find that the ionic lock exists in three states: closed (or locked), semi-open with a bridging water molecule, and open. The interconversion of these states involves the concerted motion of the entire protein. We characterize these states and the concerted motion underlying their interconversion. These findings may help elucidate the connection between key local events and the associated global structural changes during activation.

Blue matter: strong scaling of molecular dynamics on blue gene/l

This paper presents strong scaling performance data for the Blue Matter molecular dynamics framework using a novel n-body spatial decomposition and a collective communications technique implemented on both MPI and low level hardware interfaces. Using Blue Matter on Blue Gene/L, we have measured scalability through 16,384 nodes with measured time per time-step of under 2.3 milliseconds for a 43,222 atom protein/lipid system. This is equivalent to a simulation rate of over 76 nanoseconds per day and represents an unprecedented time-to-solution for biomolecular simulation as well as continued speed-up to fewer than three atoms per node. On a smaller, solvated lipid system with 13,758 atoms, we have achieved continued speedups through fewer than one atom per node and less than 2 milliseconds/time-step. On a 92,224 atom system, we have achieved floating point performance of over 1.8 TeraFlops/second on 16,384 nodes. Strong scaling of fixed-size classical molecular dynamics of biological systems to large numbers of nodes is necessary to extend the simulation time to the scale required to make contact with experimental data and derive biologically relevant insights.

Validating and improving elastic network models with molecular dynamics simulations.

Elastic network models (ENMs) are a class of simple models intended to represent the collective motions of proteins. In contrast to all‐atom molecular dynamics simulations, the low computational investment required to use an ENM makes them ideal for speculative hypothesis‐testing situations. Historically, ENMs have been validated via comparison to crystallographic B‐factors, but this comparison is relatively low‐resolution and only tests the predictions of relative flexibility. In this work, we systematically validate and optimize a number of ENM‐type models by quantitatively comparing their predictions to microsecond‐scale all‐atom simulations of three different G protein coupled receptors. We show that, despite their apparent simplicity, well‐optimized ENMs perform remarkably well, reproducing the protein fluctuations with an accuracy comparable to what one would expect from all‐atom simulations run for several hundred nanoseconds. Proteins 2010.

LOOS: An extensible platform for the structural analysis of simulations

We have developed LOOS (Lightweight Object-Oriented Structure-analysis library) as an object-oriented library designed to facilitate the rapid development of tools for the structural analysis of simulations. LOOS supports the native file formats of most common simulation packages including AMBER, CHARMM, CNS, Gromacs, NAMD, Tinker, and X-PLOR. Encapsulation and polymorphism are used to simultaneously provide a stable interface to the programmer and make LOOS easily extensible. A rich atom selection language based on the C expression syntax is included as part of the library. LOOS enables students and casual programmer-scientists to rapidly write their own analytical tools in a compact and expressive manner resembling scripting. LOOS is written in C++ and makes extensive use of the Standard Template Library and Boost, and is freely available under the GNU General Public License (version 3) LOOS has been tested on Linux and MacOS X, but is written to be portable and should work on most Unix-based platforms.