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

Residue-Residue Mutual Work Analysis of Retinal-Opsin Interaction in Rhodopsin: Implications for Protein-Ligand Binding

Energetic contributions at the single-residue level for retinal-opsin interactions in rhodopsin were studied by combining molecular dynamics simulations, transition path sampling, and a newly developed energy decomposition approach. The virtual work at an infinitesimal time interval was decomposed into the work components on one residue due to its interaction with another residue, which were then averaged over the transition path ensemble along a proposed reaction coordinate. Such residue-residue mutual work analysis on 62 residues within the active center of rhodopsin resulted in a very sparse interaction matrix, which is generally not symmetric but antisymmetric to some extent. Fourteen residues were identified to be major players in retinal relaxation along a plausible pathway from bathorhodopsin to the blue-shifted intermediate, which is in good agreement with an existing NMR study. Based on the matrix of mutual work, a comprehensive network was constructed to provide detailed insights into the chromophore-protein interaction from a viewpoint of energy flow.

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

Recent Insights from Molecular Dynamics Simulations for G Protein-Coupled Receptor Drug Discovery

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. Some ligands that bind to GPCRs activate different downstream signaling pathways. G protein activation, or β-arrestin biased signaling, involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. β-arrestin biased signaling has become a hot target for structure-based drug discovery. However, challenges include that there are few crystal structures available in the Protein Data Bank and that GPCRs are highly dynamic. Hence, molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators’ interactions with receptors and ligands. Here, we highlight recent MD simulation studies and enhanced sampling methods used to study biased G protein-coupled receptor signaling and their conformational dynamics as well as applications to drug discovery.

GPCR structure and function relationship: identification of a biased apelin receptor mutant

Biased ligands of G protein-coupled receptors (GPCRs) may have improved therapeutic benefits and safety profiles. However, the molecular mechanism of GPCR biased signaling remains largely unknown. Using apelin receptor (APJ) as a model, we systematically investigated the potential effects of amino acid residues around the orthosteric binding site on biased signaling. We discovered that a single residue mutation I109A (I109^3.32) in the transmembrane domain 3 (TM3) located in the deep ligand-binding pocket was sufficient to convert a balanced APJ into a G protein signaling biased receptor. APJ I109A mutant receptor retained full capabilities in ligand binding and G protein activation, but was defective in GRK recruitment, β-arrestin recruitment, and downstream receptor-mediated ERK activation. Based on molecular dynamics simulations, we proposed a molecular mechanism for biased signaling of I109A mutant receptor. We postulate that due to the extra space created by I109A mutation, the phenyl group of the last residue (Phe-13) of apelin rotates down and initiates a cascade of conformational changes in TM3. Phe-13 formed a new cluster of hydrophobic interactions with the sidechains of residues in TM3, including F110^3.33 and M113^3.36, which stabilizes the mutant receptor in a conformation favoring biased signaling. Interruption of these stabilizing interactions by double mutation F110A/I109A or M113A/I109A largely restored the β-arrestin-mediated signaling. Taken together, we describe herein the discovery of a biased APJ mutant receptor and provide detailed molecular insights into APJ signaling selectivity, facilitating the discovery of novel therapeutics targeting APJ.

GPCRs: What Can We Learn from Molecular Dynamics Simulations?

Advances in the structural biology of G-protein Coupled Receptors have resulted in a significant step forward in our understanding of how this important class of drug targets function at the molecular level. However, it has also become apparent that they are very dynamic molecules, and moreover, that the underlying dynamics is crucial in shaping the response to different ligands. Molecular dynamics simulations can provide unique insight into the dynamic properties of GPCRs in a way that is complementary to many experimental approaches. In this chapter, we describe progress in three distinct areas that are particularly difficult to study with other techniques: atomic level investigation of the conformational changes that occur when moving between the various states that GPCRs can exist in, the pathways that ligands adopt during binding/unbinding events and finally, the influence of lipids on the conformational dynamics of GPCRs.

Extension and evaluation of the multilevel summation method for fast long-range electrostatics calculations

Several extensions and improvements have been made to the multilevel summation method (MSM) of computing long-range electrostatic interactions. These include pressure calculation, an improved error estimator, faster direct part calculation, extension to non-orthogonal (triclinic) systems, and parallelization using the domain decomposition method. MSM also allows fully non-periodic long-range electrostatics calculations which are not possible using traditional Ewald-based methods. In spite of these significant improvements to the MSM algorithm, the particle-particle particle-mesh (PPPM) method was still found to be faster for the periodic systems we tested on a single processor. However, the fast Fourier transforms (FFTs) that PPPM relies on represent a major scaling bottleneck for the method when running on many cores (because the many-to-many communication pattern of the FFT becomes expensive) and MSM scales better than PPPM when using a large core count for two test problems on Sandia's Redsky machine. This FFT bottleneck can be reduced by running PPPM on only a subset of the total processors. MSM is most competitive for relatively low accuracy calculations. On Sandia's Chama machine, however, PPPM is found to scale better than MSM for all core counts that we tested. These results suggest that PPPM is usually more efficient than MSM for typical problems running on current high performance computers. However, further improvements to MSM algorithm could increase its competitiveness for calculation of long-range electrostatic interactions.

Generalized QM/MM Force Matching Approach Applied to the 11-cis Protonated Schiff Base Chromophore of Rhodopsin

We extended a previously developed force matching approach to systems with covalent QM/MM boundaries and describe its user-friendly implementation in the publicly available software package CPMD. We applied this approach to the challenging case of the retinal protonated Schiff base in dark state bovine rhodopsin. We were able to develop a highly accurate force field that is able to capture subtle structural changes within the chromophore that have a pronounced influence on the optical properties. The optical absorption spectrum calculated from configurations extracted from a MD trajectory using the new force field is in excellent agreement with QM/MM and experimental references.

Retinal conformation governs pKa of protonated Schiff base in rhodopsin activation

We have explored the relationship between conformational energetics and the protonation state of the Schiff base in retinal, the covalently bound ligand responsible for activating the G protein-coupled receptor rhodopsin, using quantum chemical calculations. Guided by experimental structural determinations and large-scale molecular simulations on this system, we examined rotation about each bond in the retinal polyene chain, for both the protonated and deprotonated states that represent the dark and photoactivated states, respectively. Particular attention was paid to the torsional degrees of freedom that determine the shape of the molecule, and hence its interactions with the protein binding pocket. While most torsional degrees of freedom in retinal are characterized by large energetic barriers that minimize structural fluctuations under physiological temperatures, the C6-C7 dihedral defining the relative orientation of the β-ionone ring to the polyene chain has both modest barrier heights and a torsional energy surface that changes dramatically with protonation of the Schiff base. This surprising coupling between conformational degrees of freedom and protonation state is further quantified by calculations of the pKa as a function of the C6-C7 dihedral angle. Notably, pKa shifts of greater than two units arise from torsional fluctuations observed in molecular dynamics simulations of the full ligand-protein-membrane system. It follows that fluctuations in the protonation state of the Schiff base occur prior to forming the activated MII state. These new results shed light on important mechanistic aspects of retinal conformational changes that are involved in the activation of rhodopsin in the visual process.

Molecular dynamics simulations of the adenosine A2a receptor: structural stability, sampling, and convergence

Molecular dynamics (MD) simulations of membrane-embedded G-protein coupled receptors (GPCRs) have rapidly gained popularity among the molecular simulation community in recent years, a trend which has an obvious link to the tremendous pharmaceutical importance of this group of receptors and the increasing availability of crystal structures. In view of the widespread use of this technique, it is of fundamental importance to ensure the reliability and robustness of the methodologies so they yield valid results and enable sufficiently accurate predictions to be made. In this work, 200 ns simulations of the A2a adenosine receptor (A2a AR) have been produced and evaluated in the light of these requirements. The conformational dynamics of the target protein, as obtained from replicate simulations in both the presence and absence of an inverse agonist ligand (ZM241385), have been investigated and compared using principal component analysis (PCA). Results show that, on this time scale, convergence of the replicates is not readily evident and dependent on the types of the protein motions considered. Thus rates of inter- as opposed to intrahelical relaxation and sampling can be different. When studied individually, we find that helices III and IV have noticeably greater stability than helices I, II, V, VI, and VII in the apo form. The addition of the inverse agonist ligand greatly improves the stability of all helices.

Action of molecular switches in GPCRs--theoretical and experimental studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

Model of Abnormal Chromophore-Protein Interaction for Е181К Rhodopsin Mutation: Computer Molecular Dynamics Study

The interaction of the 11-cis-retinal chromophore with the surrounding amino acid residues in the chromophore center of the rhodopsin protein has been investigated for the Е181К mutant form using molecular dynamics simulation. A comparative analysis of the arrangement of the amino acid residues in the chromophore center has been performed for both wild (native) and mutant rhodopsins. It is shown that for the Е181К mutant rhodopsin there is no proper binding of 11-cis-retinal with the surrounding amino acid residues. The distortion of the conformation states in the mutant rhodopsin molecule takes place in both the chromophore center and cytoplasmic domain. Our simulations suggest that a stable covalent linkage of 11-cis-retinal with the protein part (viz. opsin) of the rhodopsin molecule will not form. This, on the other hand, implies that the protein's active site in the cytoplasmic domain, which is responsible for the G-protein binding (so-called transducin), may not be completely blocked.Based on our molecular simulation data, we discuss the possible correlation between retinitis pigmentosa pathogenesis and the structural and functional properties of the rhodopsin protein.

Molecular simulations and solid-state NMR investigate dynamical structure in rhodopsin activation

Rhodopsin has served as the primary model for studying G protein-coupled receptors (GPCRs)-the largest group in the human genome, and consequently a primary target for pharmaceutical development. Understanding the functions and activation mechanisms of GPCRs has proven to be extraordinarily difficult, as they are part of a complex signaling cascade and reside within the cell membrane. Although X-ray crystallography has recently solved several GPCR structures that may resemble the activated conformation, the dynamics and mechanism of rhodopsin activation continue to remain elusive. Notably solid-state ((2))H NMR spectroscopy provides key information pertinent to how local dynamics of the retinal ligand change during rhodopsin activation. When combined with molecular mechanics simulations of proteolipid membranes, a new paradigm for the rhodopsin activation process emerges. Experiment and simulation both suggest that retinal isomerization initiates the rhodopsin photocascade to yield not a single activated structure, but rather an ensemble of activated conformational states. This article is part of a Special Issue entitled: Membrane protein structure and function.

Thermal properties of rhodopsin: insight into the molecular mechanism of dim-light vision

Rhodopsin has developed mechanisms to optimize its sensitivity to light by suppressing dark noise and enhancing quantum yield. We propose that an intramolecular hydrogen-bonding network formed by ∼20 water molecules, the hydrophilic residues, and peptide backbones in the transmembrane region is essential to restrain thermal isomerization, the source of dark noise. We studied the thermal stability of rhodopsin at 55 °C with single point mutations (E181Q and S186A) that perturb the hydrogen-bonding network at the active site. We found that the rate of thermal isomerization increased by 1-2 orders of magnitude in the mutants. Our results illustrate the importance of the intact hydrogen-bonding network for dim-light detection, revealing the functional roles of water molecules in rhodopsin. We also show that thermal isomerization of 11-cis-retinal in solution can be catalyzed by wild-type opsin and that this catalytic property is not affected by the mutations. We characterize the catalytic effect and propose that it is due to steric interactions in the retinal-binding site and increases quantum yield by predetermining the trajectory of photoisomerization. Thus, our studies reveal a balancing act between dark noise and quantum yield, which have opposite effects on the thermal isomerization rate. The acquisition of the hydrogen-bonding network and the tuning of the steric interactions at the retinal-binding site are two important factors in the development of dim-light vision.

Conserved water-mediated hydrogen bond network between TM-I, -II, -VI, and -VII in 7TM receptor activation

Five highly conserved polar residues connected by a number of structural water molecules together with two rotamer micro-switches, TrpVI:13 and TyrVII:20, constitute an extended hydrogen bond network between the intracellular segments of TM-I, -II, -VI, and -VII of 7TM receptors. Molecular dynamics simulations showed that, although the fewer water molecules in rhodopsin were relatively movable, the hydrogen bond network of the beta2-adrenergic receptor was fully loaded with water molecules that were surprisingly immobilized between the two rotamer switches, both apparently being in their closed conformation. Manipulations of the rotamer state of TyrVII:20 and TrpVI:13 demonstrated that these residues served as gates for the water molecules at the intracellular and extracellular ends of the hydrogen bond network, respectively. TrpVI:13 at the bottom of the main ligand-binding pocket was shown to apparently function as a catching trap for water molecules. Mutational analysis of the beta2-adrenergic receptor demonstrated that the highly conserved polar residues of the hydrogen bond network were all important for receptor signaling but served different functions, some dampening constitutive activity (AsnI:18, AspII:10, and AsnVII:13), whereas others (AsnVII:12 and AsnVII:16) located one helical turn apart and sharing a water molecule were shown to be essential for agonist-induced signaling. It is concluded that the conserved water hydrogen bond network of 7TM receptors constitutes an extended allosteric interface between the transmembrane segments being of crucial importance for receptor signaling and that part of the function of the rotamer micro-switches, TyrVII:20 and TrpVI:13, is to gate or trap the water molecules.

Seven transmembrane receptors as shapeshifting proteins: the impact of allosteric modulation and functional selectivity on new drug discovery

It is useful to consider seven transmembrane receptors (7TMRs) as disordered proteins able to allosterically respond to a number of binding partners. Considering 7TMRs as allosteric systems, affinity and efficacy can be thought of in terms of energy flow between a modulator, conduit (the receptor protein), and a number of guests. These guests can be other molecules, receptors, membrane-bound proteins, or signaling proteins in the cytosol. These vectorial flows of energy can yield standard canonical guest allostery (allosteric modification of drug effect), effects along the plane of the cell membrane (receptor oligomerization), or effects directed into the cytosol (differential signaling as functional selectivity). This review discusses these apparently diverse pharmacological effects in terms of molecular dynamics and protein ensemble theory, which tends to unify 7TMR behavior toward cells. Special consideration will be given to functional selectivity (biased agonism and biased antagonism) in terms of mechanism of action and potential therapeutic application. The explosion of technology that has enabled observation of diverse 7TMR behavior has also shown how drugs can have multiple (pluridimensional) efficacies and how this can cause paradoxical drug classification and nomenclatures.

Progress in elucidating the structural and dynamic character of G Protein-Coupled Receptor oligomers for use in drug discovery

G Protein-Coupled Receptors (GPCRs) are the most targeted group of proteins for the development of therapeutic drugs. Until the last decade, structural information about this family of membrane proteins was relatively scarce, and their mechanisms of ligand binding and signal transduction were modeled on the assumption that GPCRs existed and functioned as monomeric entities. New crystal structures of native and engineered GPCRs, together with important biochemical and biophysical data that reveal structural details of the activation mechanism(s) of this receptor family, provide a valuable framework to improve dynamic molecular models of GPCRs with the ultimate goal of elucidating their allostery and functional selectivity. Since the dynamic movements of single GPCR protomers are likely to be affected by the presence of neighboring interacting subunits, oligomeric arrangements should be taken into account to improve the predictive ability of computer-assisted structural models of GPCRs for effective use in drug design.

Light activation of rhodopsin: insights from molecular dynamics simulations guided by solid-state NMR distance restraints

Structural restraints provided by solid-state NMR measurements of the metarhodopsin II intermediate are combined with molecular dynamics simulations to help visualize structural changes in the light activation of rhodopsin. Since the timescale for the formation of the metarhodopsin II intermediate (>1 ms) is beyond that readily accessible by molecular dynamics, we use NMR distance restraints derived from 13C dipolar recoupling measurements to guide the simulations. The simulations yield a working model for how photoisomerization of the 11-cis retinylidene chromophore bound within the interior of rhodopsin is coupled to transmembrane helix motion and receptor activation. The mechanism of activation that emerges is that multiple switches on the extracellular (or intradiscal) side of rhodopsin trigger structural changes that converge to disrupt the ionic lock between helices H3 and H6 on the intracellular side of the receptor.

Retinal conformation and dynamics in activation of rhodopsin illuminated by solid-state H NMR spectroscopy

Solid-state NMR spectroscopy gives a powerful avenue for investigating G protein-coupled receptors and other integral membrane proteins in a native-like environment. This article reviews the use of solid-state (2)H NMR to study the retinal cofactor of rhodopsin in the dark state as well as the meta I and meta II photointermediates. Site-specific (2)H NMR labels have been introduced into three regions (methyl groups) of retinal that are crucially important for the photochemical function of rhodopsin. Despite its phenomenal stability (2)H NMR spectroscopy indicates retinal undergoes rapid fluctuations within the protein binding cavity. The spectral lineshapes reveal the methyl groups spin rapidly about their three-fold (C(3)) axes with an order parameter for the off-axial motion of SC(3) approximately 0.9. For the dark state, the (2)H NMR structure of 11-cis-retinal manifests torsional twisting of both the polyene chain and the beta-ionone ring due to steric interactions of the ligand and the protein. Retinal is accommodated within the rhodopsin binding pocket with a negative pretwist about the C11=C12 double bond. Conformational distortion explains its rapid photochemistry and reveals the trajectory of the 11-cis to trans isomerization. In addition, (2)H NMR has been applied to study the retinylidene dynamics in the dark and light-activated states. Upon isomerization there are drastic changes in the mobility of all three methyl groups. The relaxation data support an activation mechanism whereby the beta-ionone ring of retinal stays in nearly the same environment, without a large displacement of the ligand. Interactions of the beta-ionone ring and the retinylidene Schiff base with the protein transmit the force of the retinal isomerization. Solid-state (2)H NMR thus provides information about the flow of energy that triggers changes in hydrogen-bonding networks and helix movements in the activation mechanism of the photoreceptor.

How to lose a kink and gain a helix: pH independent conformational changes of the fusion domains from influenza hemagglutinin in heterogeneous lipid bilayers

We have simulated two conformations of the fusion domain of influenza hemagglutinin (HA) within explicit water, salt, and heterogeneous lipid bilayers composed of POPC:POPG (4:1). Each conformation has seven different starting points in which the initial peptide structure is the same for each conformation, but the location across the membrane normal and lipid arrangement around the peptide are varied, giving a combined total simulation time of 140 ns. For the HA5 conformation (primary structure from recent NMR spectroscopy at pH = 5), the peptide exhibits a stable and less kinked structure in the lipid bilayer compared to that from the NMR studies. The relative fusogenic behavior of the different conformations has been investigated by calculation of the relative free energy of insertion into the hydrophobic region of lipid bilayer as a function of the depth of immersion. For the HA7 conformations (primary structure from recent NMR spectroscopy at pH = 7.4), while the N-terminal helix preserves its initial structure, the flexible C-terminal chain produces a transient helical motif inside the lipid bilayer. This conformational change is pH-independent, and is closely related to the peptide insertion into the lipid bilayer.

Atomistic insights into rhodopsin activation from a dynamic model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

Molecular simulations of lipid-mediated protein-protein interactions

Recent experimental results revealed that lipid-mediated interactions due to hydrophobic forces may be important in determining the protein topology after insertion in the membrane, in regulating the protein activity, in protein aggregation and in signal transduction. To gain insight into the lipid-mediated interactions between two intrinsic membrane proteins, we developed a mesoscopic model of a lipid bilayer with embedded proteins, which we studied with dissipative particle dynamics. Our calculations of the potential of mean force between transmembrane proteins show that hydrophobic forces drive long-range protein-protein interactions and that the nature of these interactions depends on the length of the protein hydrophobic segment, on the three-dimensional structure of the protein and on the properties of the lipid bilayer. To understand the nature of the computed potentials of mean force, the concept of hydrophilic shielding is introduced. The observed protein interactions are interpreted as resulting from the dynamic reorganization of the system to maintain an optimal hydrophilic shielding of the protein and lipid hydrophobic parts, within the constraint of the flexibility of the components. Our results could lead to a better understanding of several membrane processes in which protein interactions are involved.

Mechanism of signal propagation upon retinal isomerization: insights from molecular dynamics simulations of rhodopsin restrained by normal modes

As one of the best studied members of the pharmaceutically relevant family of G-protein-coupled receptors, rhodopsin serves as a prototype for understanding the mechanism of G-protein-coupled receptor activation. Here, we aim at exploring functionally relevant conformational changes and signal transmission mechanisms involved in its photoactivation brought about through a cis-trans photoisomerization of retinal. For this exploration, we propose a molecular dynamics simulation protocol that utilizes normal modes derived from the anisotropic network model for proteins. Deformations along multiple low-frequency modes of motion are used to efficiently sample collective conformational changes in the presence of explicit membrane and water environment, consistent with interresidue interactions. We identify two highly stable regions in rhodopsin, one clustered near the chromophore, the other near the cytoplasmic ends of transmembrane helices H1, H2, and H7. Due to redistribution of interactions in the neighborhood of retinal upon stabilization of the trans form, local structural rearrangements in the adjoining H3-H6 residues are efficiently propagated to the cytoplasmic end of these particular helices. In the structures obtained by our simulations, all-trans retinal interacts with Cys(167) on H4 and Phe(203) on H5, which were not accessible in the dark state, and exhibits stronger interactions with H5, while some of the contacts made (in the cis form) with H6 are lost.

Activation of G protein-coupled receptors: beyond two-state models and tertiary conformational changes

Transformation of G protein-coupled receptors (GPCRs) from a quiescent to an active state initiates signal transduction. All GPCRs share a common architecture comprising seven transmembrane-spanning alpha-helices, which accommodates signal propagation from a diverse repertoire of external stimuli across biological membranes to a heterotrimeric G protein. Signal propagation through the transmembrane helices likely involves mechanistic features common to all GPCRs. The structure of the light receptor rhodopsin may serve as a prototype for the transmembrane architecture of GPCRs. Early biochemical, biophysical, and pharmacological studies led to the conceptualization of receptor activation based on the context of two-state equilibrium models and conformational changes in protein structure. More recent studies indicate a need to move beyond these classical paradigms and to consider additional aspects of the molecular character of GPCRs, such as the oligomerization and dynamics of the receptor.

Ligand-stabilized conformational states of human beta(2) adrenergic receptor: insight into G-protein-coupled receptor activation

G-protein-coupled receptors (GPCRs) are known to exist in dynamic equilibrium between inactive- and several active-state conformations, even in the absence of a ligand. Recent experimental studies on the beta(2) adrenergic receptor (beta(2)AR) indicate that structurally different ligands with varying efficacies trigger distinct conformational changes and stabilize different receptor conformations. We have developed a computational method to study the ligand-induced rotational orientation changes in the transmembrane helices of GPCRs. This method involves a systematic spanning of the rotational orientation of the transmembrane helices (TMs) that are in the vicinity of the ligand for predicting the helical rotations that occur on ligand binding. The predicted ligand-stabilized receptor conformations are characterized by a simultaneous lowering of the ligand binding energy and a significant gain in interhelical and receptor-ligand hydrogen bonds. Using the beta(2)AR as a model, we show that the receptor conformational state depends on the structure and efficacy of the ligand for a given signaling pathway. We have studied the ligand-stabilized receptor conformations of five different ligands, a full agonist, norepinephrine; a partial agonist, salbutamol; a weak partial agonist, dopamine; a very weak agonist, catechol; and an inverse agonist, ICI-115881. The predicted ligand-stabilized receptor models correlate well with the experimentally observed conformational switches in beta(2)AR, namely, the breaking of the ionic lock between R131(3.50) at the intracellular end of TM3 (part of the DRY motif) and E268(6.30) on TM6, and the rotamer toggle switch on W286(6.48) on TM6. In agreement with trp-bimane quenching experiments, we found that norepinephrine and dopamine break the ionic lock and engage the rotamer toggle switch, whereas salbutamol, a noncatechol partial agonist only breaks the ionic lock, and the weak agonist catechol only engages the rotamer toggle switch. Norepinephrine and dopamine occupy the same binding region, between TM3, TM5, and TM6, whereas the binding site of salbutamol is shifted toward TM4. Catechol binds deeper into the protein cavity compared to the other ligands, making contact with TM5 and TM6. A part of the catechol binding site overlaps with those of dopamine and norepinephrine but not with that of salbutamol. Virtual ligand screening on 10,060 ligands on the norepinephrine-stabilized receptor conformation shows an enrichment of 38% compared to ligand unbound receptor conformation. These results show that ligand-induced conformational changes are important for developing functionally specific drugs that will stabilize a particular receptor conformation. These studies represent the first step toward a more universally applicable computational method for studying ligand efficacy and GPCR activation.

Advances in the Development and Application of Computational Methodologies for Structural Modeling of G-Protein Coupled Receptors

BACKGROUND: Despite the large amount of experimental data accumulated in the past decade on G-protein coupled receptor (GPCR) structure and function, understanding of the molecular mechanisms underlying GPCR signaling is still far from being complete, thus impairing the design of effective and selective pharmaceuticals. OBJECTIVE: Understanding of GPCR function has been challenged even further by more recent experimental evidence that several of these receptors are organized in the cell membrane as homo- or hetero-oligomers, and that they may exhibit unique pharmacological properties. Given the complexity of these new signaling systems, researcher's efforts are turning increasingly to molecular modeling, bioinformatics and computational simulations for mechanistic insights of GPCR functional plasticity. METHODS: We review here current advances in the development and application of computational approaches to improve prediction of GPCR structure and dynamics, thus enhancing current understanding of GPCR signaling. RESULTS/CONCLUSIONS: Models resulting from use of these computational approaches further supported by experiments are expected to help elucidate the complex allosterism that propagates through GPCR complexes, ultimately aiming at successful structure-based rational drug design.

Terahertz spectroscopy of bacteriorhodopsin and rhodopsin: similarities and differences

We studied the low-frequency terahertz spectroscopy of two photoactive protein systems, rhodopsin and bacteriorhodopsin, as a means to characterize collective low-frequency motions in helical transmembrane proteins. From this work, we found that the nature of the vibrational motions activated by terahertz radiation is surprisingly similar between these two structurally similar proteins. Specifically, at the lowest frequencies probed, the cytoplasmic loop regions of the proteins are highly active; and at the higher terahertz frequencies studied, the extracellular loop regions of the protein systems become vibrationally activated. In the case of bacteriorhodopsin, the calculated terahertz spectra are compared with the experimental terahertz signature. This work illustrates the importance of terahertz spectroscopy to identify vibrational degrees of freedom which correlate to known conformational changes in these proteins.

Functional and structural roles of conserved cysteine residues in the carboxyl-terminal domain of the follicle-stimulating hormone receptor in human embryonic kidney 293 cells

The carboxyl-terminal segment of G protein-coupled receptors has one or more conserved cysteine residues that are potential sites for palmitoylation. This posttranslational modification contributes to membrane association, internalization, and membrane targeting of proteins. In contrast to other members of the glycoprotein hormone receptor family (the LH and thyroid-stimulating hormone receptors), it is not known whether the follicle-stimulating hormone receptor (FSHR) is palmitoylated and what are the effects of abolishing its potential palmitoylation sites. In the present study, a functional analysis of the FSHR carboxyl-terminal segment cysteine residues was carried out. We constructed a series of mutant FSHRs by substituting cysteine residues with alanine, serine, or threonine individually and together at positions 629 and 655 (conserved cysteines) and 627 (nonconserved). The results showed that all three cysteine residues are palmitoylated but that only modification at Cys629 is functionally relevant. The lack of palmitoylation does not appear to greatly impair coupling to G(s) but, when absent at position 629, does significantly impair cell surface membrane expression of the partially palmitoylated receptor. All FSHR Cys mutants were capable of binding agonist with the same affinity as the wild-type receptor and internalizing on agonist stimulation. Molecular dynamics simulations at a time scale of approximately 100 nsec revealed that replacement of Cys629 resulted in structures that differed significantly from that of the wild-type receptor. Thus, deviations from wild-type conformation may potentially contribute to the severe impairment in plasma membrane expression and the modest effects on signaling exhibited by the receptors modified in this particular position.

The signaling pathway of rhodopsin

The signal-transduction mechanism of rhodopsin was studied by molecular dynamics (MD) simulations of the high-resolution, inactive structure in an explicit membrane environment. The simulations were employed to calculate equal-time correlations of the fluctuating interaction energy of residue pairs. The resulting interaction-correlation matrix was used to determine a network that couples retinal to the cytoplasmic interface, where transducin binds. Two highly conserved motifs, D(E)RY and NPxxY, were found to have strong interaction correlation with retinal. MD simulations with restraints on each transmembrane helix indicated that the major signal-transduction pathway involves the interdigitating side chains of helices VI and VII. The functional roles of specific residues were elucidated by the calculated effect of retinal isomerization from 11-cis to all-trans on the residue-residue interaction pattern. It is suggested that Glu134 may act as a "signal amplifier" and that Asp83 may introduce a threshold to prevent background noise from activating rhodopsin.

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. We analyzed a total of 23 independent, 100 ns all-atom molecular dynamics simulations of rhodopsin embedded in a lipid bilayer in the microcanonical (N,V,E) ensemble. Analysis of intramolecular fluctuations predicts hydrogen-out-of-plane (HOOP) wagging modes of retinal consistent with those found in Raman vibrational spectroscopy. We show that sampling and ergodicity of the ensemble of simulations are crucial for determining the distribution of conformers of retinal bound to rhodopsin. The polyene chain is rigidly locked into a single, twisted conformation, consistent with the function of retinal as an inverse agonist in the dark state. Most surprisingly, the beta-ionone ring is mobile within its binding pocket; interactions are non-specific and the cavity is sufficiently large to enable structural heterogeneity. We find that retinal occupies two distinct conformations in the dark state, contrary to most previous assumptions. The beta-ionone ring can rotate relative to the polyene chain, thereby populating both positively and negatively twisted 6-s-cis enantiomers. This result, while unexpected, strongly agrees with experimental solid-state (2)H NMR spectra. Correlation analysis identifies the residues most critical to controlling mobility of retinal; we find that Trp265 moves away from the ionone ring prior to any conformational transition. Our findings reinforce how molecular dynamics simulations can challenge conventional assumptions for interpreting experimental data, especially where existing models neglect conformational fluctuations.

Flexibility of ras lipid modifications studied by 2H solid-state NMR and molecular dynamics simulations

Human posttranslationally modified N-ras oncogenes are known to be implicated in numerous human cancers. Here, we applied a combination of experimental and computational techniques to determine structural and dynamical details of the lipid chain modifications of an N-ras heptapeptide in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes. Experimentally, 2H NMR spectroscopy was used to study oriented membranes that incorporated ras heptapeptides with two covalently attached perdeuterated hexadecyl chains. Atomistic molecular dynamics simulations of the same system were carried out over 100 ns including 60 DMPC and 4 ras molecules. Several structural and dynamical experimental parameters could be directly compared to the simulation. Experimental and simulated 2H NMR order parameters for the methylene groups of the ras lipid chains exhibited a systematic difference attributable to the absence of collective motions in the simulation and to geometrical effects. In contrast, experimental 2H NMR spin-lattice relaxation rates for Zeeman order were well reproduced in the simulation. The lack of slower collective motions in the simulation did not appreciably influence the relaxation rates at a Larmor frequency of 115.1 MHz. The experimental angular dependence of the 2H NMR relaxation rates with respect to the external magnetic field was also relatively well simulated. These relaxation rates showed a weak angular dependence, suggesting that the lipid modifications of ras are very flexible and highly mobile in agreement with the low order parameters. To quantify these results, the angular dependence of the 2H relaxation rates was calculated by an analytical model considering both molecular and collective motions. Peptide dynamics in the membrane could be modeled by an anisotropic diffusion tensor with principal values of Dparallel=2.1x10(9) s(-1) and Dperpendicular=4.5x10(5) s(-1). A viscoelastic fitting parameter describing the membrane elasticity, viscosity, and temperature was found to be relatively similar for the ras peptide and the DMPC host matrix. Large motional amplitudes and relatively short correlation times facilitate mixing and dispersal with the lipid bilayer matrix, with implications for the role of the full-length ras protein in signal transduction and oncogenesis.

Molecular dynamics simulations of GlpF in a micelle vs in a bilayer: conformational dynamics of a membrane protein as a function of environment

Octyl glucoside (OG) is a detergent widely employed in structural and functional studies of membrane proteins. To better understand the nature of protein-OG interactions, molecular dynamics simulations (duration 10 ns) have been used to explore an alpha-helical membrane protein, GlpF, in OG micelles and in DMPC bilayers. Greater conformational drift of the extramembraneous protein loops, from the initial X-ray structure, is seen for the GlpF-OG simulations than for the GlpF-DMPC simulation. The mobility of the transmembrane alpha-helices is approximately 1.3x higher in the GlpF-OG than the GlpF-DMPC simulations. The detergent is seen to form an irregular torus around the protein. The presence of the protein leads to a small perturbation in the behavior of the alkyl chains in the OG micelle, namely an approximately 15% increase in the trans-gauche(-)-gauche(+) transition time. Aromatic side chains (Trp, Tyr) and basic side chains (Arg, Lys) play an important role in both protein-detergent (OG) and protein-lipid (DMPC) interactions.

Molecular Dynamics Simulations of Rhodopsin in Different One-Component Lipid Bilayers

Four 20 ns molecular dynamic simulations of rhodopsin embedded in different one-component lipid bilayers have been carried out to ascertain the importance of membrane lipids on the protein structure. Specifically, dimyristoyl phosphatidylcholine (DMPC), dipalmitoyl phosphatidylcholine (DPPC), palmitoyl oleoyl phosphatidylcholine (POPC), and palmitoyl linoleyl phosphatidylcholine (PLPC) lipid bilayers have been considered for the present work. The results reported here provide information on the hydrophobic matching between the protein and the bilayer and about the differential effects of the protein on the thickness of the different membranes. Furthermore, a careful analysis of the individual protein-lipid interactions permits the identification of residues that exhibit permanent interactions with atoms of the lipid environment that may putatively act as hooks of the protein to the membrane. The analysis of the trajectories also provides information about the effect of the bilayer on the protein structure, including secondary structural elements, salt bridges, and rigid-body motions.

Models of beta-amyloid ion channels in the membrane suggest that channel formation in the bilayer is a dynamic process

Here we model the Alzheimer beta-peptide ion channel with the goal of obtaining insight into the mechanism of amyloid toxicity. The models are built based on NMR data of the oligomers, with the universal U-shaped (strand-turn-strand) motif. After 30-ns simulations in the bilayer, the channel dimensions, shapes and subunit organization are in good agreement with atomic force microscopy (AFM). The models use the Abeta(17-42) pentamer NMR-based coordinates. Extension and bending of the straight oligomers lead to two channel topologies, depending on the direction of the curvature: 1), the polar/charged N-terminal beta-strand of Abeta(17-42) faces the water-filled pore, and the hydrophobic C-terminal beta-strand faces the bilayer (CNpNC; p for pore); and 2), the C-terminal beta-strand faces the solvated pore (NCpCN). In the atomistic simulations in a fully solvated DOPC lipid bilayer, the first (CNpNC) channel preserves the pore and conducts solvent; by contrast, hydrophobic collapse blocks the NCpCN channel. AFM demonstrated open pores and collapsed complexes. The final averaged CNpNC pore dimensions (outer diameter 8 nm; inner diameter approximately 2.5 nm) are in the AFM range (8-12 nm; approximately 2 nm, respectively). Further, in agreement with high-resolution AFM images, during the simulations, the channels spontaneously break into ordered subunits in the bilayer; however, we also observe that the subunits are loosely connected by partially disordered inner beta-sheet, suggesting subunit mobility in the bilayer. The cationic channel has strong selective affinity for Ca(2+), supporting experimental calcium-selective beta-amyloid channels. Membrane permeability and consequent disruption of calcium homeostasis were implicated in cellular degeneration. Consequently, the CNpNC channel topology can sign cell death, offering insight into amyloid toxicity via an ion "trap-release" transport mechanism. The observed loosely connected subunit organization suggests that amyloid channel formation in the bilayer is a dynamic, fluid process involving subunit association, dissociation, and channel rearrangements.

Molecular dynamics of rhodopsin and free opsin: computer simulation

Computer simulation was used to perform a comparative study of the molecular dynamics of rhodopsin containing the chromophore group (11-cis-retinal) and free opsin. The molecular dynamics were followed over a time interval of 3000 psec; a total of 3 x 10(6) discrete conformational states of rhodopsin and opsin. The presence of the chromophore group in the chromophore center of opsin was shown to have significant effects on the immediate protein environment of the chromophore and the conformational state of the cytoplasmic domain, but to have virtually no effect on the conformational state of the intradisk domain. The simulation results are used to discuss the possible intramolecular mechanism by which rhodopsin is maintained as a G-protein-coupled receptor in the inactive state, i.e., the function of the chromophore as an effective antagonist ligand.

Unambiguous assignment of intramolecular chemical cross-links in modified mammalian membrane proteins by Fourier transform-tandem mass spectrometry

Fourier transform tandem mass spectrometry (FT-MS/MS) can be used to unambiguously assign intramolecular chemical cross-links to specific amino acid residues even when two or more possible cross-linking sites are adjacent in the cross-linked protein. Bovine rhodopsin (Rho) in its dark-adapted state was intramolecularly cross-linked with lysine-cysteine (K-C) or lysine-lysine (K-K) cross-linkers to obtain interatomic distance information. Large, multiply charged, cross-linked peptide ions containing adjacent lysines, corresponding to Rho(50-86) (K(66) or K(67)) cross-linked to Rho(310-317) (C(316)) or Rho(318-348) (K(325) or K(339)), were fragmented by collision-induced dissociation (CID), infrared multiphoton dissociation (IRMPD), and electron capture dissociation (ECD). Complementary sequence-specific information was obtained by combining cross-link assignments; however, only ECD revealed full palmitoylation of adjacent cysteines (C(322) and C(323)) and cross-linking of K(67) (and not K(66)) to C(316), K(325), and K(339). ECD spectra contained crucial c- and z-ions resulting from cleavage of the bond between K(66) and K(67). To our knowledge, this work also presents the first demonstration that ECD can be used to characterize S-linked fatty acid acylation on cysteines. The comprehensive fragmentation of large peptides by CID, IRMPD, and particularly ECD, in conjunction with the high resolution and mass accuracy of FT-MS/MS, is shown to be a valuable means of characterizing mammalian membrane proteins with both chemical and posttranslational modifications.

Predisposition of the dark state of rhodopsin to functional changes in structure

As the only member of the family of G-protein-coupled receptors for which atomic coordinates are available, rhodopsin is widely studied for insight into the molecular mechanism of G-protein-coupled receptor activation. The currently available structures refer to the inactive, dark state, of rhodopsin, rather than the light-activated metarhodopsin II (Meta II) state. A model for the Meta II state is proposed here by analyzing elastic network normal modes in conjunction with experimental data. Key mechanical features and interactions broken/formed in the proposed model are found to be consistent with the experimental data. The model is further tested by using a set of Meta II fluorescence decay rates measured to empirically characterize the deactivation of rhodopsin mutants. The model is found to correctly predict 93% of the experimentally observed effects in 119 rhodopsin mutants for which the decay rates and misfolding data have been measured, including a systematic analysis of Cys-->Ser replacements reported here. Based on the detailed comparison between model and experiments, a cooperative activation mechanism is deduced that couples retinal isomerization to concerted changes in conformation, facilitated by the intrinsic dynamics of rhodopsin. A global hinge site is identified near the retinal-binding pocket that ensures the efficient propagation of signals from the central transmembrane region to both cytoplasmic and extracellular ends. The predicted activation mechanism opens the transmembrane helices at the critical G-protein binding cytoplasmic domain. This model provides a detailed, mechanistic description of the activation process, extending experimental observations and yielding new insights for further tests.

How a small change in retinal leads to G-protein activation: initial events suggested by molecular dynamics calculations

Rhodopsin is the prototypical G-protein coupled receptor, coupling light activation with high efficiency to signaling molecules. The dark-state X-ray structures of the protein provide a starting point for consideration of the relaxation from initial light activation to conformational changes that may lead to signaling. In this study we create an energetically unstable retinal in the light activated state and then use molecular dynamics simulations to examine the types of compensation, relaxation, and conformational changes that occur following the cis-trans light activation. The results suggest that changes occur throughout the protein, with changes in the orientation of Helices 5 and 6, a closer interaction between Ala 169 on Helix 4 and retinal, and a shift in the Schiff base counterion that also reflects changes in sidechain interactions with the retinal. Taken together, the simulation is suggestive of the types of changes that lead from local conformational change to light-activated signaling in this prototypical system.

Introduction to the membrane protein reviews: the interplay of structure, dynamics, and environment in membrane protein function

In our review, we introduce an organizational scheme for membrane protein function. It is the relationship between structure, dynamics, and environment that endows the membrane and its constituents with remarkable sensitivity and robustness. Our understanding begins with landmark advances like those presented in the following chapters. Membrane proteins are notoriously difficult to study, and so the work presented here on the ADP/ATP carrier [Nury et al. ( 2 )], rhodopsin [Palczewski ( 24 )], and the cytochrome b6f complex [Cramer et al. ( 35 )] represents incredible progress in this now blossoming field.

Effect of different treatments of long-range interactions and sampling conditions in molecular dynamic simulations of rhodopsin embedded in a dipalmitoyl phosphatidylcholine bilayer

The present study analyzes the effect of the simulation conditions on the results of molecular dynamics simulations of G-protein coupled receptors (GPCRs) performed with an explicit lipid bilayer. Accordingly, the present work reports the analysis of different simulations of bovine rhodopsin embedded in a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer using two different sampling conditions and two different approaches for the treatment of long-range electrostatic interactions. Specifically, sampling was carried out either by using the statistical ensembles NVT or NPT (constant number of atoms, a pressure of 1 atm in all directions and fixed temperature), and the electrostatic interactions were treated either by using a twin-cutoff, or the particle mesh Ewald summation method (PME). The results of the present study suggest that the use of the NPT ensemble in combination with the PME method provide more realistic simulations. The use of NPT during the equilibration avoids the need of an a priori estimation of the box dimensions, giving the correct area per lipid. However, once the system is equilibrated, the simulations are irrespective of the sampling conditions used. The use of an electrostatic cutoff induces artifacts on both lipid thickness and the ion distribution, but has no direct effect on the protein and water molecules.

Computational investigations of pore forming peptide assemblies in lipid bilayers

This Letter presents the first application of a three-dimensional numerical molecular theory based modeling approach to study the structure and energetics of assemblies of peptides embedded in lipid bilayers. Coarse-grained models were used for both the peptides and lipids. Both barrel-stave and toroidal pore morphologies for the lipids near the peptide assemblies are found, but at different assembly sizes. The free energy of the assembly is found to have a global free energy minimum for a solution with a membrane-spanning toroidal pore. A pairwise approximation to this free energy is found to underpredict the free energy minimum associated with the membrane-spanning pores.

Interaction of protegrin-1 with lipid bilayers: membrane thinning effect

Protegrins (PG) are important in defending host tissues, preventing infection via an attack on the membrane surface of invading microorganisms. Protegrins have powerful antibiotic abilities, but the molecular-level mechanisms underlying the interactions of their beta-sheet motifs with the membrane are not known. Protegrin-1 (PG-1) is composed of 18 amino acids with a high content of basic residues and two disulfide bonds. Here we focused on the stability of PG-1 at the amphipathic interface in lipid bilayers and on the details of the peptide-membrane interactions. We simulated all-atom models of the PG-1 monomer with explicit water and lipid bilayers composed of both homogeneous POPC (palmitoyl-oleyl-phosphatidylcholine) lipids and a mixture of POPC/POPG (palmitoyl-oleyl-phosphatidylglycerol) (4:1) lipids. We observed that local thinning of the lipid bilayers mediated by the peptide is enhanced in the lipid bilayer containing POPG, consistent with experimental results of selective membrane targeting. The beta-hairpin motif of PG-1 is conserved in both lipid settings, whereas it is highly bent in aqueous solution. The conformational dynamics of PG-1, especially the highly charged beta-hairpin turn region, are found to be mostly responsible for disturbing the membrane. Even though the eventual membrane disruption requires PG-1 oligomers, our simulations clearly show the first step of the monomeric effects. The thinning effects in the bilayer should relate to pore/channel formation in the lipid bilayer and thus be responsible for further defects in the membrane caused by oligomer.

Structure and dynamics of dark-state bovine rhodopsin revealed by chemical cross-linking and high-resolution mass spectrometry

Recent work using chemical cross-linking to define interresidue distance constraints in proteins has shown that these constraints are useful for testing tertiary structural models. We applied this approach to the G-protein-coupled receptor bovine rhodopsin in its native membrane using lysine- and cysteine-targeted bifunctional cross-linking reagents. Cross-linked proteolytic peptides of rhodopsin were identified by combined liquid chromatography and FT-ICR mass spectrometry with automated data-reduction and assignment software. Tandem mass spectrometry was used to verify cross-link assignments and locate the exact sites of cross-link attachment. Cross-links were observed to form between 10 pairs of residues in dark-state rhodopsin. For each pair, cross-linkers with a range of linker lengths were tested to determine an experimental distance-of-closest-approach (DCA) between reactive side-chain atoms. In all, 28 cross-links were identified using seven different cross-linking reagents. Molecular mechanics procedures were applied to published crystal structure data to calculate energetically achievable theoretical DCAs between reactive atoms without altering the position of the protein backbone. Experimentally measured DCAs are generally in good agreement with the theoretical DCAs. However, a cross-link between C316 and K325 in the C-terminal region cannot be rationalized by DCA simulations and suggests that backbone reorientation relative to the crystal coordinates occurs on the timescale of cross-linking reactions. Biochemical and spectroscopic data from other studies have found that the C-terminal region is highly mobile in solution and not fully represented by X-ray crystallography data. Our results show that chemical cross-linking can provide reliable three-dimensional structural information and insight into local conformational dynamics in a membrane protein.

Dynamic models of G-protein coupled receptor dimers: indications of asymmetry in the rhodopsin dimer from molecular dynamics simulations in a POPC bilayer

Based on the growing evidence that G-protein coupled receptors (GPCRs) form homo- and hetero-oligomers, models of GPCR signaling are now considering macromolecular assemblies rather than monomers, with the homo-dimer regarded as the minimal oligomeric arrangement required for functional coupling to the G-protein. The dynamic mechanisms of such signaling assemblies are unknown. To gain some insight into properties of GPCR dimers that may be relevant to functional mechanisms, we study their current structural prototype, rhodopsin. We have carried out nanosecond time-scale molecular dynamics (MD) simulations of a rhodopsin dimer and compared the results to the monomer simulated in the same type of bilayer membrane model composed of an equilibrated unit cell of hydrated palmitoyl-oleoyl-phosphatidyl choline (POPC). The dynamic representation of the homo-dimer reveals the location of structural changes in several regions of the monomeric subunits. These changes appear to be more pronounced at the dimerization interface that had been shown to be involved in the activation process [Proc Natl Acad Sci USA 102:17495, 2005]. The results are consistent with a model of GPCR activation that involves allosteric modulation through a single GPCR subunit per dimer.