Molecular Dynamics Simulation of Dark-adapted Rhodopsin in an Explicit Membrane Bilayer: Coupling between Local Retinal and Larger Scale Conformational Change

3D modeling of the activated states of constitutively active mutants of rhodopsin

The activated (R*) states in constitutively active mutants (CAMs) of G-protein-coupled receptors (GPCRs) are presumably characterized by lower energies than the resting (R) states. If specific configurations of TM helices differing by rotations along the long transmembrane axes possess energies lower than that in the R state for pronounced CAMs, but not for non-CAMs, these particular configurations of TM helices are candidate 3D models for the R* state. The hypothesis was studied in the case of rhodopsin, the only GPCR for which experimentally determined 3D models of the R and R* states are currently available. Indeed, relative energies of the R* state were significantly lower than that of the R state for the rhodopsin mutants G90D/M257Y and E113Q/M257Y (strong CAMs), but not for G90D, E113Q, and M257Y (not CAMs). Next, the developed build-up procedure successfully identified few similar configurations of the TM helical bundle of G90D/M257Y and E113Q/M257Y as possible candidates for the 3D model of the R* state of rhodopsin, all of them being in good agreement with the model suggested by experiment. Since constitutively active mutants are known for many of GPCRs belonging to the large rhodopsin-like family, this approach provides a way for predicting possible 3D structures corresponding to the activated states of the TM regions of many GPCRs for which CAMs have been identified.

Conformational sampling and dynamics of membrane proteins from 10-nanosecond computer simulations

In the current report, we provide a quantitative analysis of the convergence of the sampling of conformational space accomplished in molecular dynamics simulations of membrane proteins of duration in the order of 10 nanoseconds. A set of proteins of diverse size and topology is considered, ranging from helical pores such as gramicidin and small beta-barrels such as OmpT, to larger and more complex structures such as rhodopsin and FepA. Principal component analysis of the C(alpha)-atom trajectories was employed to assess the convergence of the conformational sampling in both the transmembrane domains and the whole proteins, while the time-dependence of the average structure was analyzed to obtain single-domain information. The membrane-embedded regions, particularly those of small or structurally simple proteins, were found to achieve reasonable convergence. By contrast, extra-membranous domains lacking secondary structure are often markedly under-sampled, exhibiting a continuous structural drift. This drift results in a significant imprecision in the calculated B-factors, which detracts from any quantitative comparison to experimental data. In view of such limitations, we suggest that similar analyses may be valuable in simulation studies of membrane protein dynamics, in order to attach a level of confidence to any biologically relevant observations.

Origin and Consequences of Steric Strain in the Rhodopsin Binding Pocket

To study the origin and the effects of steric strain on the chromophore conformation in rhodopsin, we have performed quantum-mechanical calculations on the wild-type retinal chromophore and four retinal derivatives, 13-demethyl-, 10-methyl-13-demethyl-, 10-methyl-, and 9-demethylretinal. For the dynamics of the whole protein, a combined quantum mechanics/molecular mechanics method (DFTB/CHARMM) was used and for the calculation of excited-state properties the nonempirical CASSCF/CASPT2 method. After relaxation inside the protein, all chromophores show significant nonplanar distortions from C10 to C13, most strongly for 10-methylretinal and least pronounced for 9-demethylretinal. In all five cases, the dihedral angle of the C10-C11=C12-C13 bond is negative which attests to the strong chiral discrimination exerted by the protein pocket. The calculations show that the nonplanar distortion of the chromophore, including the sense of rotation, is caused by a combination of two effects: the fitting of both ends to the protein matrix which imposes a distance constraint and the bonding arrangement at the Schiff base terminus. With both the counterion Glu113 and Lys296 displaced off the plane of the chromophore, their binding to N16 exerts a torque on the chromophore. As a result, the polyene chain, from N16 to C13, is twisted in a clockwise manner against the remaining part of the chromophore, leading to a C11=C12 bond with the observed negative dihedral angle. Shifts of the absorption maxima are reproduced correctly, in particular, the red shift of the 10-methyl and the strong blue shift of the 9-demethyl analogue relative to the wild type. Calculated positive rotatory strengths of the alpha-CD bands are in agreement with the calculated absolute conformation of the mutant chromophores.

Molecular dynamics simulations of proteins in lipid bilayers

With recent advances in X-ray crystallography of membrane proteins promising many new high-resolution structures, molecular dynamics simulations will become increasingly valuable for understanding membrane protein function, as they can reveal the dynamic behavior concealed in the static structures. Dramatic increases in computational power, in synergy with more efficient computational methodologies, now allow us to carry out molecular dynamics simulations of any structurally known membrane protein in its native environment, covering timescales of up to 0.1 micros. At the frontiers of membrane protein simulations are ion channels, aquaporins, passive and active transporters, and bioenergetic proteins.

Molecular dynamics simulations of retinal in rhodopsin: from the dark-adapted state towards lumirhodopsin

The formation of photointermediates and conformational changes observed in the retinal chromophore of bilayer-embedded rhodopsin during the early steps of the protein activation have been studied by molecular dynamics (MD) simulation. In particular, the lysine-bound retinal has been examined, focusing on its conformation in the dark-adapted state (10 ns) and on the early steps after the isomerization of the 11-cis bond to trans (up to 10 ns). The parametrization for the chromophore is based on a recent quantum study [Sugihara, M., Buss, V., Entel, P., Elstner, M., and Frauenheim, T. (2002) Biochemistry 41, 15259-15266] and shows good conformational agreement with recent experimental results. The isomerization, induced by switching the function governing the dihedral angle for the C11=C12 bond, was repeated with several different starting conformations. From the repeated simulations, it is shown that the retinal model exhibits a conserved activation pattern. The conformational changes are sequential and propagate outward from the C11=C12 bond, starting with isomerization of the C11=C12 bond, then a rotation of methyl group C20, and followed by increased fluctuations at the beta-ionone ring. The dynamics of these changes suggest that they are linked with photointermediates observed by spectroscopy. The exact moment when these events occur after the isomerization is modulated by the starting conformation, suggesting that retinal isomerizes through multiple pathways that are slightly different. The amplitudes of the structural fluctuations observed for the protein in the dark-adapted state and after isomerization of the retinal are similar, suggesting a subtle mechanism for the transmission of information from the chromophore to the protein.

The second intracellular loop of the calcitonin gene-related peptide receptor provides molecular determinants for signal transduction and cell surface expression

The calcitonin gene-related peptide (CGRP) receptor is a heterodimer of a family B G-protein-coupled receptor, calcitonin receptor-like receptor (CLR), and the accessory protein receptor activity modifying protein 1. It couples to G(s), but it is not known which intracellular loops mediate this. We have identified the boundaries of this loop based on the relative position and length of the juxtamembrane transmembrane regions 3 and 4. The loop has been analyzed by systematic mutagenesis of all residues to alanine, measuring cAMP accumulation, CGRP affinity, and receptor expression. Unlike rhodopsin, ICL2 of the CGRP receptor plays a part in the conformational switch after agonist interaction. His-216 and Lys-227 were essential for a functional CGRP-induced cAMP response. The effect of (H216A)CLR is due to a disruption to the cell surface transport or surface stability of the mutant receptor. In contrast, (K227A)CLR had wild-type expression and agonist affinity, suggesting a direct disruption to the downstream signal transduction mechanism of the CGRP receptor. Modeling suggests that the loop undergoes a significant shift in position during receptor activation, exposing a potential G-protein binding pocket. Lys-227 changes position to point into the pocket, potentially allowing it to interact with bound G-proteins. His-216 occupies a position similar to that of Tyr-136 in bovine rhodopsin, part of the DRY motif of the latter receptor. This is the first comprehensive analysis of an entire intracellular loop within the calcitonin family of G-protein-coupled receptor. These data help to define the structural and functional characteristics of the CGRP-receptor and of family B G-protein-coupled receptors in general.

Molecular dynamics simulations of bovine rhodopsin: influence of protonation states and different membrane-mimicking environments

G-protein coupled receptors (GPCRs) are a protein family of outstanding pharmaceutical interest. GPCR homology models, based on the crystal structure of bovine rhodopsin, have been shown to be valuable tools in the drug-design process. The initial model is often refined by molecular dynamics (MD) simulations, a procedure that has been recently discussed controversially. We therefore analyzed MD simulations of bovine rhodopsin in order to identify contacts that could serve as constraints in the simulation of homology models. Additionally, the effect of an N-terminal truncation, the nature of the membrane mimic, the influence of varying protonation states of buried residues and the importance of internal water molecules was analyzed. All simulations were carried out using the program-package GROMACS. While N-terminal truncation negatively influenced the overall protein stability, a stable simulation was possible in both solvent environments. As regards the protonation state of titratable sites, the experimental data could be reproduced by the program UHBD (University of Houston Brownian Dynamics), suggesting its application for studying homology models of GPCRs. A high flexibility was observed for internal water molecules at some sites. Finally, interhelical hydrogen-bonding interactions could be derived, which can now serve as constraints in the simulations of GPCR homology models.

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. The latter phenomenon suggests a molecular-level mechanism for the experimental finding that cholesterol stabilizes the native dark-adapted state of rhodopsin without binding directly to the protein.

How environment supports a state: molecular dynamics simulations of two states in bacteriorhodopsin suggest lipid and water compensation

The light-driven proton pump bacteriorhodopsin (bR) is a transmembrane protein that uses large conformational changes for proton transfer from the cytoplasmic to the extracellular regions. Crystal structures, due to their solvent conditions, do not resolve the effect of lipid molecules on these protein conformational changes. To begin to understand the molecular details behind such large conformational changes, we simulated two conformations of wild-type bacteriorhodopsin, one of the dark-adapted state and the second of an intermediate (M(O)) state, each within an explicit dimyristoyl-phosphatidylcholine (DMPC) lipid bilayer. The simulations included all-hydrogen and all-atom representations of protein, lipid, and water and were performed for 20 ns. We investigate the equilibrium properties and the dynamic motions of the two conformations in the lipid setting. We note that the conformational state of the M(O) intermediate bR remains markedly different from the dark-adapted bR state in that the M(O) intermediate shows rearrangement of the cytoplasmic portions of helices C, F, and G, and nearby loops. This difference in the states remained throughout the simulations, and the results are stable on the molecular dynamics timescale and provide an illustration of the changes in both lipid and water that help to stabilize a particular state. Our analysis focuses on how the environment adjusts to these two states and on how the dynamics of the helices, loops, and water molecules can be related to the pump mechanism of bacteriorhodopsin. For example, water generally behaves in the same manner on the extracellular sides of both simulations but is decreased in the cytoplasmic region of the M(O) intermediate. We suspect that the different water behavior is closely related to the fluctuations of microcavities volume in the protein interior, which is strongly coupled to the collective motion of the protein. Our simulation result suggests that experimental observation can be useful to verify a decreased number of waters in the cytoplasmic regions of the late-intermediate stages by measuring the rate of water exchange with the interior of the protein.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

How lipids and proteins interact in a membrane: a molecular approach

Membrane proteins in a biological membrane are surrounded by a shell or annulus of 'solvent' lipid molecules. These lipid molecules in general interact rather non-specifically with the protein molecules, although a few 'hot-spots' may be present on the protein where anionic lipids bind with high affinity. Because of the low structural specificity of most of the annular sites, the composition of the lipid annulus will be rather similar to the bulk lipid composition of the membrane. The structures of the solvent lipid molecules are important in determining the conformational state of a membrane protein, and hence its activity, through charge and hydrogen bonding interactions between the lipid headgroups and residues in the protein, and through hydrophobic matching between the protein and the surrounding lipid bilayer. Evidence is also accumulating for the presence of 'co-factor' lipid molecules binding with high specificity to membrane proteins, often between transmembrane alpha-helices, and often being essential for activity.

Tools for channels: moving towards molecular calculations of gating and permeation in ion channel biophysics

Recent X-ray structures of voltage gated potassium channels provide an exciting opportunity to connect molecular structures with measured biological function. Two of the most important connections for these channels are: first, to the molecular basis behind selectivity and the associated free energy profile underlying ionic current flow and, second, to a true molecular understanding of the large-scale conformational transitions that underlie voltage dependent gating. But, existing computational tools need to be further developed to reach these goals. In this contribution to the symposia on sampling methods we outline our dynamic importance sampling method for sampling large-scale conformational transitions as well as our studies with non-equilibrium work events and equilibrium overlap sampling (OS) methods for sampling events related to the calculation of relative free energies.

The retinal conformation and its environment in rhodopsin in light of a new 2.2 A crystal structure

A new high-resolution structure is reported for bovine rhodopsin, the visual pigment in rod photoreceptor cells. Substantial improvement of the resolution limit to 2.2 A has been achieved by new crystallization conditions, which also reduce significantly the probability of merohedral twinning in the crystals. The new structure completely resolves the polypeptide chain and provides further details of the chromophore binding site including the configuration about the C6-C7 single bond of the 11-cis-retinal Schiff base. Based on both an earlier structure and the new improved model of the protein, a theoretical study of the chromophore geometry has been carried out using combined quantum mechanics/force field molecular dynamics. The consistency between the experimental and calculated chromophore structures is found to be significantly improved for the 2.2 A model, including the angle of the negatively twisted 6-s-cis-bond. Importantly, the new crystal structure refinement reveals significant negative pre-twist of the C11-C12 double bond and this is also supported by the theoretical calculation although the latter converges to a smaller value. Bond alternation along the unsaturated chain is significant, but weaker in the calculated structure than the one obtained from the X-ray data. Other differences between the experimental and theoretical structures in the chromophore binding site are discussed with respect to the unique spectral properties and excited state reactivity of the chromophore.

Lipid-protein interactions of integral membrane proteins: a comparative simulation study

The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.