Scott E. Feller

Lipid Bilayer Structure Determined by the Simultaneous Analysis of Neutron and X-Ray Scattering Data

Quantitative structures were obtained for the fully hydrated fluid phases of dioleoylphosphatidylcholine (DOPC) and dipalmitoylphosphatidylcholine (DPPC) bilayers by simultaneously analyzing x-ray and neutron scattering data. The neutron data for DOPC included two solvent contrasts, 50% and 100% D(2)O. For DPPC, additional contrast data were obtained with deuterated analogs DPPC_d62, DPPC_d13, and DPPC_d9. For the analysis, we developed a model that is based on volume probability distributions and their spatial conservation. The models design was guided and tested by a DOPC molecular dynamics simulation. The model consistently captures the salient features found in both electron and neutron scattering density profiles. A key result of the analysis is the molecular surface area, A. For DPPC at 50 degrees C A = 63.0 A(2), whereas for DOPC at 30 degrees C A = 67.4 A(2), with estimated uncertainties of 1 A(2). Although A for DPPC agrees with a recently reported value obtained solely from the analysis of x-ray scattering data, A for DOPC is almost 10% smaller. This improved method for determining lipid areas helps to reconcile long-standing differences in the values of lipid areas obtained from stand-alone x-ray and neutron scattering experiments and poses new challenges for molecular dynamics simulations.

Constant surface tension simulations of lipid bilayers: The sensitivity of surface areas and compressibilities

Eight molecular dynamics simulations of a hydrated lipid bilayer have been carried out differing only in the applied surface tension, γ, defining the boundary conditions of the periodic cell. The calculated surface area per molecule and deuterium order parameter profile are found to depend strongly on γ. We present several methods to calculate the area compressibility modulus, KA, from the simulations. Equivalence between the constant area and constant surface tension ensembles is investigated by comparing the present simulations with earlier work from our laboratories and we find simulation results to depend much more strongly on the specified surface area or surface tension than on the ensemble employed.

Molecular dynamics simulations of lipid bilayers

Computer simulation methods are becoming increasingly widespread as tools for studying the structure and dynamics of lipid bilayer membranes. The length scale and time scale accessible to atomic-level molecular dynamics simulations are rapidly increasing, providing insight into the relatively slow motions of molecular reorientation and translation and demonstrating that effects due to the finite size of the simulation cell can influence simulation results. Additionally, significant advances have been made in the complexity of membrane systems studied, including bilayers with cholesterol, small solute molecules, and lipid-protein and lipid-DNA complexes. Especially promising is the progress that continues to be made in the comparison of simulation results with experiment, both to validate the simulation algorithms and to aid in the interpretation of existing experimental data.

Calculation of the dielectric permittivity profile for a nonuniform system: Application to a lipid bilayer simulation

We derive an expression relating the static dielectric permittivity profile for a system nonuniform in one dimension to correlations between the net system dipole moment and the local polarization density. The permittivity profile of a dipalmitoylphosphatidylcholine (DPPC) lipid bilayer in water is calculated from an all-atom 20-ns molecular dynamics simulation. The component of the permittivity parallel to the bilayer shows a nonmonotonic decrease from the value in bulk water to the value in the membrane interior; the interfacial region itself has a very large permittivity, greater than that of bulk water. In high-dielectric regions, obtaining a quantitative estimate of the component normal to the bilayer is not possible because of large numerical uncertainty. However, the calculated correlation function is consistent with a value for the interface at least as large as that of bulk water. In general, the transition to a low-dielectric environment is sharp and is located on the inner border of the region ...

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.

Nuclear Overhauser Enhancement Spectroscopy Cross-Relaxation Rates and Ethanol Distribution across Membranes

Measurement of nuclear Overhauser enhancement spectroscopy cross-relaxation rates between ethanol and palmitoyloleoylphosphatidylcholine bilayers was combined with atomic-level molecular dynamics simulations. The molecular dynamics trajectories yielded autocorrelation functions of proton dipole-dipole interactions, and, consequently, relaxation times and cross-relaxation rates. These analyses allow the measured cross-relaxation rates to be interpreted in terms of relative interaction strengths with the various segments of the lipid molecule. We determined that cross-relaxation between ethanol and specific lipid resonances is primarily determined by the sites of interaction with some modulation due to lipid disorder and to local differences in intramolecular lipid dynamics. The rates scale linearly with the lifetime of temporary ethanol-lipid associations. Ethanol interacts with palmitoyloleoylphosphatidylcholine bilayers primarily via hydrophilic interactions, in particular the formation of hydrogen bonds to the lipid phosphate group. There is a weak contribution to binding from hydrophobic interaction with lipid chain segments near the glycerol. However, the strength of hydrophobic interactions is insufficient to compensate for the energetic loss of locating ethanol in an exclusively hydrophobic environment, resulting in a probability of locating ethanol in the bilayer center that is three orders of magnitude lower than locating ethanol at the lipid/water interface. The low cross-relaxation rates between terminal methyl protons of hydrocarbon chains and ethanol are as much the result of infrequent chain upturns as of brief excursions of ethanol into the region of lipid hydrocarbon chains near the glycerol. The combination of nuclear magnetic resonance measurements and molecular dynamics simulations offers a general pathway to study the interaction of small molecules with the lipid matrix at atomic resolution.

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.

Acyl chain conformations in phospholipid bilayers: a comparative study of docosahexaenoic acid and saturated fatty acids

A variety of experimental methods indicate unique biophysical properties of membranes containing the highly polyunsaturated omega-3 fatty acid, docosahexaenoic acid (DHA). In the following we review the atomically detailed picture of DHA acyl chains structure and dynamics that has emerged from computational studies of this system in our lab. A comprehensive approach, beginning with ab-initio quantum chemical studies of model compounds representing segments of DHA and ending with large scale classical molecular dynamics simulations of DHA-containing bilayers, is described with particular attention paid to contrasting the properties of DHA with those of saturated fatty acids. Connection with experiment is made primarily through comparison with Nuclear Magnetic Resonance (NMR) studies, particularly those that probe details of the chain structure and dynamics. Our computational results suggest that low torsional energy barriers, comparable to kT at physiological conditions, for the rotatable bonds in the DHA chain are the key to the differences observed between polyunsaturated and saturated acyl chains.