Activation of the G-Protein-Coupled Receptor Rhodopsin by Water.

Visual rhodopsin is an important archetype for G-protein-coupled receptors, which are membrane proteins implicated in cellular signal transduction. Herein, we show experimentally that approximately 80 water molecules flood rhodopsin upon light absorption to form a solvent-swollen active state. An influx of mobile water is necessary for activating the photoreceptor, and this finding is supported by molecular dynamics (MD) simulations.

A Usual G-Protein-Coupled Receptor in Unusual Membranes.

G-protein-coupled receptors (GPCRs) are the largest family of membrane-bound receptors and constitute about 50% of all known drug targets. They offer great potential for membrane protein nanotechnologies. We report here a charge-interaction-directed reconstitution mechanism that induces spontaneous insertion of bovine rhodopsin, the eukaryotic GPCR, into both lipid- and polymer-based artificial membranes.

Structural basis of G protein-coupled receptor-Gi protein interaction: formation of the cannabinoid CB2 receptor-Gi protein complex.

In this study, we applied a comprehensive G protein-coupled receptor-Gαi protein chemical cross-linking strategy to map the cannabinoid receptor subtype 2 (CB2)-Gαi interface and then used molecular dynamics simulations to explore the dynamics of complex formation. Three cross-link sites were identified using LC-MS/MS and electrospray ionization-MS/MS as follows: 1) a sulfhydryl cross-link between C3.53(134) in TMH3 and the Gαi C-terminal i-3 residue Cys-351; 2) a lysine cross-link between K6.

Retinal ligand mobility explains internal hydration and reconciles active rhodopsin structures.

Rhodopsin, the mammalian dim-light receptor, is one of the best-characterized G-protein-coupled receptors, a pharmaceutically important class of membrane proteins that has garnered a great deal of attention because of the recent availability of structural information. Yet the mechanism of rhodopsin activation is not fully understood. Here, we use microsecond-scale all-atom molecular dynamics simulations, validated by solid-state (2)H nuclear magnetic resonance spectroscopy, to understand the transition between the dark and metarhodopsin I (Meta I) states.

Global fold of human cannabinoid type 2 receptor probed by solid-state 13C-, 15N-MAS NMR and molecular dynamics simulations.

The global fold of human cannabinoid type 2 (CB2 ) receptor in the agonist-bound active state in lipid bilayers was investigated by solid-state (13)C- and (15)N magic-angle spinning (MAS) NMR, in combination with chemical-shift prediction from a structural model of the receptor obtained by microsecond-long molecular dynamics (MD) simulations. Uniformly (13)C- and (15)N-labeled CB2 receptor was expressed in milligram quantities by bacterial fermentation, purified, and functionally reconstituted into liposomes. (13)C MAS NMR spectra were recorded without sensitivity enhancement for direct comparison of Cα, Cβ, and C=O bands of superimposed resonances with predictions from protein structures generated by MD.

Molecular dynamics simulations reveal specific interactions of post-translational palmitoyl modifications with rhodopsin in membranes.

We present a detailed analysis of the behavior of the highly flexible post-translational lipid modifications of rhodopsin from multiple-microsecond all-atom molecular dynamics simulations. Rhodopsin was studied in a realistic membrane environment that includes cholesterol, as well as saturated and polyunsaturated lipids with phosphocholine and phosphoethanolamine headgroups. The simulation reveals striking differences between the palmitoylations at Cys322 and Cys323 as well as between the palmitoyl chains and the neighboring lipids.

Petascale computation performance of lightweight multiscale cardiac models using hybrid programming models.

Future multiscale and multiphysics models must use the power of high performance computing (HPC) systems to enable research into human disease, translational medical science, and treatment. Previously we showed that computationally efficient multiscale models will require the use of sophisticated hybrid programming models, mixing distributed message passing processes (e.g.

Dynamics of confined reactive water in smectite clay-zeolite composites.

The dynamics of water confined to mesoporous regions in minerals such as swelling clays and zeolites is fundamental to a wide range of resource management issues impacting many processes on a global scale, including radioactive waste containment, desalination, and enhanced oil recovery. Large-scale atomic models of freely diffusing multilayer smectite particles at low hydration confined in a silicalite cage are used to investigate water dynamics in the composite environment with the ReaxFF reactive force field over a temperature range of 300-647 K. The reactive capability of the force field enabled a range of relevant surface chemistry to emerge, including acid/base equilibria in the interlayer calcium hydrates and silanol formation on the edges of the clay and inner surface of the zeolite housing.

Performance of hybrid programming models for multiscale cardiac simulations: preparing for petascale computation.

Future multiscale and multiphysics models that support research into human disease, translational medical science, and treatment can utilize the power of high-performance computing (HPC) systems. We anticipate that computationally efficient multiscale models will require the use of sophisticated hybrid programming models, mixing distributed message-passing processes [e.g.

Orthogonal recursive bisection as data decomposition strategy for massively parallel cardiac simulations.

We present the orthogonal recursive bisection algorithm that hierarchically segments the anatomical model structure into subvolumes that are distributed to cores. The anatomy is derived from the Visible Human Project, with electrophysiology based on the FitzHugh-Nagumo (FHN) and ten Tusscher (TT04) models with monodomain diffusion. Benchmark simulations with up to 16,384 and 32,768 cores on IBM Blue Gene/P and L supercomputers for both FHN and TT04 results show good load balancing with almost perfect speedup factors that are close to linear with the number of cores.

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.

Concerted interconversion between ionic lock substates of the beta(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.

Orthogonal recursive bisection data decomposition for high performance computing in cardiac model simulations: dependence on anatomical geometry.

Orthogonal recursive bisection (ORB) algorithm can be used as data decomposition strategy to distribute a large data set of a cardiac model to a distributed memory supercomputer. It has been shown previously that good scaling results can be achieved using the ORB algorithm for data decomposition. However, the ORB algorithm depends on the distribution of computational load of each element in the data set.

Strong scaling and speedup to 16,384 processors in cardiac electro-mechanical simulations.

High performance computing is required to make feasible simulations of whole organ models of the heart with biophysically detailed cellular models in a clinical setting. Increasing model detail by simulating electrophysiology and mechanical models increases computation demands. We present scaling results of an electro - mechanical cardiac model of two ventricles and compare them to our previously published results using an electrophysiological model only.

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,6)F motif, through specific interactions with the receptor.

Large scale cardiac modeling on the Blue Gene supercomputer.

Multi-scale, multi-physical heart models have not yet been able to include a high degree of accuracy and resolution with respect to model detail and spatial resolution due to computational limitations of current systems. We propose a framework to compute large scale cardiac models. Decomposition of anatomical data in segments to be distributed on a parallel computer is carried out by optimal recursive bisection (ORB).

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.

Dynamic structure of retinylidene ligand of rhodopsin probed by molecular simulations.

Rhodopsin is currently the only available atomic-resolution template for understanding biological functions of the G protein-coupled receptor (GPCR) family. The structural basis for the phenomenal dark state stability of 11-cis-retinal bound to rhodopsin and its ultrafast photoreaction are active topics of research. In particular, the beta-ionone ring of the retinylidene inverse agonist is crucial for the activation mechanism.

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.

Retinal counterion switch mechanism in vision evaluated by molecular simulations.

Photoisomerization of the retinylidene chromophore of rhodopsin is the starting point in the vision cascade. A counterion switch mechanism that stabilizes the retinal protonated Schiff base (PSB) has been proposed to be an essential step in rhodopsin activation. On the basis of vibrational and UV-visible spectroscopy, two counterion switch models have emerged.

Contribution of omega-3 fatty acids to the thermodynamics of membrane protein solvation.

Recent NMR experiments and molecular dynamics simulations have indicated that rhodopsin is preferentially solvated by omega-3 fatty acids compared to saturated chains. However, to date no physical theory has been advanced to explain this phenomenon. The present work presents a novel thermodynamic explanation for this preferential solvation based on statistical analysis of 26,100 ns all-atom molecular dynamics simulations of rhodopsin in membranes rich in polyunsaturated chains.

A role for direct interactions in the modulation of rhodopsin by omega-3 polyunsaturated lipids.

Rhodopsin, the G protein-coupled receptor primarily responsible for sensing light, is found in an environment rich in polyunsaturated lipid chains and cholesterol. Biophysical experiments have shown that lipid unsaturation and cholesterol both have significant effects on rhodopsin's stability and function; omega-3 polyunsaturated chains, such as docosahexaenoic acid (DHA), destabilize rhodopsin and enhance the kinetics of the photocycle, whereas cholesterol has the opposite effect. Here, we use molecular dynamics simulations to investigate the possibility that polyunsaturated chains modulate rhodopsin stability and kinetics via specific direct interactions.

Molecular dynamics investigation of dynamical properties of phosphatidylethanolamine lipid bilayers.

We describe the dynamic behavior of a 1-stearoyl-2-oleoyl-phosphatidylethanolamine (SOPE) bilayer from a 20 ns molecular dynamics simulation. The dynamics of individual molecules are characterized in terms of (2)H spin-lattice relaxation rates, nuclear overhauser enhancement spectroscopy (NOESY) cross-relaxation rates, and lateral diffusion coefficients. Additionally, we describe the dynamics of hydrogen bonding through an analysis of hydrogen bond lifetimes and the time evolution of clusters of hydrogen bonded lipids.

Molecular dynamics investigation of the structural properties of phosphatidylethanolamine lipid bilayers.

We report a 14 ns microcanonical (NVE) molecular dynamics simulation of a fully hydrated bilayer of 1-stearoyl-2-oleoyl-phosphatidyethanolamine. This study describes the structure of the bilayer in terms of NMR order parameters and radial distribution functions, and compares them to experimental results and simulations of other lipids. A focus of this work is the characterization of the lipid-water interface, particularly the hydrogen bonding network of the phosphatidylethanolamine (PE) headgroups.

Role of cholesterol and polyunsaturated chains in lipid-protein interactions: molecular dynamics simulation of rhodopsin in a realistic membrane environment.

We present a 118-ns molecular dynamics simulation of rhodopsin embedded in a bilayer composed of a 2:2:1 mixture of 1-stearoyl-2-docosahexaenoyl-phosphatidylcholine (SDPC), 1-stearoyl-2-docosahexaenoyl-phosphatidylethanolamine (SDPE), and cholesterol, respectively. The simulation demonstrates that the protein breaks the lateral and transverse symmetry of the bilayer. Lipids near the protein preferentially reorient such that their unsaturated chains interact with the protein, while the distribution of cholesterol in the membrane complements the variations in rhodopsin's transverse profile.