Molecular Dynamics Simulations of Rhodopsin Point Mutants at the Cytoplasmic Side of Helices 3 and 6

Abstract

The present work reports on a structural analysis carried out through different computer simulations of a set of rhodopsin mutants with differential functional features in regard to the wild type. Most of these mutants, whose experimental features had previously been reported [Ramon et al. J Biol Chem 282, 14272–14282 (2007)], were designed to perturb a network of electrostatic interactions located at the cytoplasmic sides of transmembrane helices 3 and 6. Geometric and energetic features derived from the detailed analysis of a series of molecular dynamics simulations of the different rhodopsin mutants, involving positions 134(3.49), 247(6.30), and 251(6.34), suggest that the protein structure is sensitive to these mutations through the local changes induced that extend further to the secondary structure of neighboring helices and, ultimately, to the packing of the helical bundle. Overall, the results obtained highlight the complexity of the analyzed network of electrostatic interactions where the effect of each mutation on protein structure can produce rather specific features.     

Critical role of electrostatic interactions of amino acids at the cytoplasmic region of helices 3 and 6 in rhodopsin conformational properties and activation

The cytoplasmic sides of transmembrane helices 3 and 6 of G-protein-coupled receptors are connected by a network of ionic interactions that play an important role in maintaining its inactive conformation. To investigate the role of such a network in rhodopsin structure and function, we have constructed single mutants at position 134 in helix 3 and at positions 247 and 251 in helix 6, as well as combinations of these to obtain double mutants involving the two helices. These mutants have been expressed in COS-1 cells, immunopurified using the rho-1D4 antibody, and studied by UV-visible spectrophotometry. Most of the single mutations did not affect chromophore formation, but double mutants, especially those involving the T251K mutant, resulted in low yield of protein and impaired 11-cis-retinal binding. Single mutants E134Q, E247Q, and E247A showed the ability to activate transducin in the dark, and E134Q and E247A enhanced activation upon illumination, with regard to wild-type rhodopsin. Mutations E247A and T251A (in E134Q/E247A and E134Q/T251A double mutants) resulted in enhanced activation compared with the single E134Q mutant in the dark. A role for Thr(251) in this network is proposed for the first time in rhodopsin. As a result of these mutations, alterations in the hydrogen bond interactions between the amino acid side chains at the cytoplasmic region of transmembrane helices 3 and 6 have been observed using molecular dynamics simulations. Our combined experimental and modeling results provide new insights into the details of the structural determinants of the conformational change ensuing photoactivation of rhodopsin.     

Electrostatic Compensation Restores Trafficking of the Autosomal Recessive Retinitis Pigmentosa E150K Opsin Mutant to the Plasma Membrane

Rhodopsin is the rod photoreceptor G protein-coupled receptor responsible for capturing light. Mutations in the gene encoding this protein can lead to a blinding disease called retinitis pigmentosa, which is inherited frequently in an autosomal dominant manner. The E150K opsin mutant associated with rarely occurring autosomal recessive retinitis pigmentosa localizes to trans-Golgi network membranes rather than to plasma membranes of rod photoreceptor cells. We investigated the molecular mechanisms underlying opsin retention in the Golgi apparatus. Electrostatic calculations reveal that the E150K mutant features an overall accumulation of positive charges between helices H-IV and H-II. Human E150K and several other closely related opsin mutants were then expressed in HEK-293 cells. Spectral characteristics and functional biochemistry of each mutant were analyzed after reconstitution with the cis-retinoid chromophore. UV-visible spectra and rhodopsin/transducin activation assays revealed only minor differences between the purified wild type control and rhodopsin mutants. However, partial restoration of the surface electrostatic charge in the compensatory R69E/E150K double mutant rescues the plasma membrane localization of opsin. These findings emphasize the fundamental importance of electrostatic interactions for appropriate membrane trafficking of opsin and advance our understanding of the pathophysiology of autosomal recessive retinitis pigmentosa due to the E150K mutation.     

Conformational change of the transmembrane helices II and IV of metabotropic glutamate receptor involved in G protein activation

G protein-coupled receptors are classified into several families on the basis of their amino acid sequences and the members of the same family exhibit sequence similarity but those of different families do not. In family 1 GPCRs such as rhodopsin and adrenergic receptor, extensive studies have revealed the stimulus-dependent conformational change of the receptor: the rearrangement of transmembrane helices III and VI is essential for G protein activation. In contrast, in family 3 GPCRs such as metabotropic glutamate receptor (mGluR), the inter-protomer relocation upon ligand binding has been observed but there is much less information about the structural changes of the transmsmbrane helices and the cytoplasmic domains. Here we identified constitutively active mutation sites at the cytoplasmic borders of helices II and IV of mGluR8 and successfully inhibited the G protein activation ability by engineering disulfide cross-linking between these cytoplasmic regions. The analysis of all possible single substitution mutants of these residues revealed that some steric interactions around these sites would be important to keep the receptor protein inactive. These results provided the model that the conformational changes at the cytoplasmic ends of helices II and IV of mGluR are involved in the efficient G protein coupling.     

Uniformity,ideality, and hydrogen bonds in transmembrane alpha-helices

Protein environments substantially influence the balance of molecular interactions that generate structural stability. Transmembrane helices exist in the relatively uniform low dielectric interstices of the lipid bilayer, largely devoid of water and with a very hydrophobic distribution of amino acid residues. Here, through an analysis of bacteriorhodopsin crystal structures and the transmembrane helix structure from M2 protein of influenza A, some helices are shown to be exceptionally uniform in hydrogen bond geometry, peptide plane tilt angle, and backbone torsion angles. Evidence from both the x-ray crystal structures and solid-state NMR structure suggests that the intramolecular backbone hydrogen bonds are shorter than their counterparts in water-soluble proteins. Moreover, the geometry is consistent with a dominance of electrostatic versus covalent contributions to these bonds. A comparison of structure as a function of resolution shows that as the structures become better characterized the helices become much more uniform, suggesting that there is a possibility that many more uniform helices will be observed, even among the moderate resolution membrane protein structures that are currently in the Protein Data Bank that do not show such features.     

Modulation of glycophorin A transmembrane helix interactions by lipid bilayers: molecular dynamics calculations

Starting from the glycophorin A dimer structure determined by NMR, we performed simulations of both dimer and monomer forms in explicit lipid bilayers with constant normal pressure, lateral area, and temperature using the CHARMM potential. Analysis of the trajectories in four different lipids reveals how lipid chain length and saturation modulate the structural and energetic properties of transmembrane helices. Helix tilt, helix-helix crossing angle, and helix accessible volume depend on lipid type in a manner consistent with hydrophobic matching concepts: the most relevant lipid property appears to be the bilayer thickness. Although the net helix-helix interaction enthalpy is strongly attractive, analysis of residue-residue interactions reveals significant unfavorable electrostatic repulsion between interfacial glycine residues previously shown to be critical for dimerization. Peptide volume is nearly conserved upon dimerization in all lipid types, indicating that the monomeric helices pack equally well with lipid as dimer helices do with one another. Enthalpy calculations indicate that the helix-environment interaction energy is lower in the dimer than in the monomer form, when solvated by unsaturated lipids. In all lipid environments there is a marked preference for lipids to interact with peptide predominantly through one rather than both acyl chains. Although our trajectories are not long enough to allow a full thermodynamic treatment, these results demonstrate that molecular dynamics simulations are a powerful method for investigating the protein-protein, protein-lipid, and lipid-lipid interactions that determine the structure, stability and dynamics of transmembrane alpha-helices in membranes.     

Solution structure of the sixth transmembrane helix of the G-protein-coupled receptor, rhodopsin

Low resolution electron density maps have revealed the general orientation of the transmembrane helices of rhodopsin. However, high resolution structural information for the transmembrane domain of the G-protein-coupled receptor, rhodopsin, is as yet unavailable. In this study, a high resolution solution structure is reported for a 15 residue portion of the sixth transmembrane helix of rhodopsin (rhovih) as a free peptide. Helix 6 is one of the transmembrane helices of rhodopsin that contains a proline (amino acid residue 267) and the influence of this proline on the structure of this transmembrane domain was unknown. The structure obtained shows an alpha-helix through most of the sequence. The proline apparently induces only a modest distortion in the helix. Previously, the structure of the intradiskal loop connected to helix 6 was solved. The sequence of this loop contained five residues in common (residues 268-272) with the peptide reported here from the rhovih. The five residues in common between these two structures were superimposed to connect these two structures. The superposition showed a root mean square deviation of 0.2 A. Thus, this five residue sequence formed the same structure in both peptides, indicating that the structure of this region is governed primarily by short range interactions.     

The Effect of Loops on the Structural Organization of α-Helical Membrane Proteins

Loops connecting the transmembrane (TM) α-helices in membrane proteins are expected to affect the structural organization of the thereby connected helices and the helical bundles as a whole. This effect, which has been largely ignored previously, is studied here by analyzing the x-ray structures of 41 α-helical membrane proteins. First we define the loop flexibility ratio, R, and find that 53% of the loops are stretched, where a stretched loop constrains the distance between the two connected helices. The significance of this constraining effect is supported by experiments carried out with bacteriorhodopsin and rhodopsin, in which cutting or eliminating their (predominately stretched) loops has led to a decrease in protein stability, and for rhodopsin, in most cases, also to the destruction of the structure. We show that for nonstretched loops in the extramembranous regions, the fraction of hydrophobic residues is comparable to that for soluble proteins; furthermore (as is also the case for soluble proteins), the hydrophobic residues in these regions are preferentially buried. This is expected to lead to the compact structural organization of the loops, which is transferred to the TM helices, causing them to assemble. We argue that a soluble protein complexed with a membrane protein similarly promotes compactness; other properties of such complexes are also studied. We calculate complementary attractive interactions between helices, including hydrogen bonds and van der Waals interactions of sequential motifs, such as GXXXG. The relative and combined effects of all these factors on the association of the TM helices are discussed and protein structures with only a few of these factors are analyzed. Our study emphasizes the need for classifying membrane proteins into groups according to structural organization. This classification should be considered when procedures for structural analysis or prediction are developed and applied. Detailed analysis of each structure is provided at http://flan.blm.cs.cmu.edu/memloop/     

Evolution of Specificity in Protein-Protein Interactions

Hub proteins are proteins that maintain promiscuous molecular recognition. Because they are reported to play essential roles in cellular control, there has been a special interest in the study of their structural and functional properties, yet the mechanisms by which they evolve to maintain functional interactions are poorly understood. By combining biophysical simulations of coarse-grained proteins and analysis of proteins-complex crystallographic structures, we seek to elucidate those mechanisms. We focus on two types of hub proteins: Multi hubs, which interact with their partners through different interfaces, and Singlish hubs, which do so through a single interface. We show that loss of structural stability is required for the evolution of protein-protein-interaction (PPI) networks, and it is more profound in Singlish hub systems. In addition, different ratios of hydrophobic to electrostatic interfacial amino acids are shown to support distinct network topologies (i.e., Singlish and Multi systems), and therefore underlie a fundamental design principle of PPI in a crowded environment. We argue that the physical nature of hydrophobic and electrostatic interactions, in particular, their favoring of either same-type interactions (hydrophobic-hydrophobic), or opposite-type interactions (negatively-positively charged) plays a key role in maintaining the network topology while allowing the protein amino acid sequence to evolve.     

Evolution of Specificity in Protein-Protein Interactions

Hub proteins are proteins that maintain promiscuous molecular recognition. Because they are reported to play essential roles in cellular control, there has been a special interest in the study of their structural and functional properties, yet the mechanisms by which they evolve to maintain functional interactions are poorly understood. By combining biophysical simulations of coarse-grained proteins and analysis of proteins-complex crystallographic structures, we seek to elucidate those mechanisms. We focus on two types of hub proteins: Multi hubs, which interact with their partners through different interfaces, and Singlish hubs, which do so through a single interface. We show that loss of structural stability is required for the evolution of protein-protein-interaction (PPI) networks, and it is more profound in Singlish hub systems. In addition, different ratios of hydrophobic to electrostatic interfacial amino acids are shown to support distinct network topologies (i.e., Singlish and Multi systems), and therefore underlie a fundamental design principle of PPI in a crowded environment. We argue that the physical nature of hydrophobic and electrostatic interactions, in particular, their favoring of either same-type interactions (hydrophobic-hydrophobic), or opposite-type interactions (negatively-positively charged) plays a key role in maintaining the network topology while allowing the protein amino acid sequence to evolve.     

Molecular dynamics simulations of human prion protein: importance of correct treatment of electrostatic interactions.

Molecular dynamics simulations have been used to investigate the dynamical and structural behavior of a homology model of human prion protein HuPrP(90-230) generated from the NMR structure of the Syrian hamster prion protein ShPrP(90-231) and of ShPrP(<90-231) itself. These PrPs have a large number of charged residues on the protein surface. At the simulation pH 7, HuPrP(90-230) has a net charge of -1 eu from 15 positively and 14 negatively charged residues. Simulations for both PrPs, using the AMBER94 force field in a periodic box model with explicit water molecules, showed high sensitivity to the correct treatment of the electrostatic interactions. Highly unstable behavior of the structured region of the PrPs (127-230) was found using the truncation method, and stable trajectories could be achieved only by including all the long-range electrostatic interactions using the particle mesh Ewald (PME) method. The instability using the truncation method could not be reduced by adding sodium and chloride ions nor by replacing some of the sodium ions with calcium ions. The PME simulations showed, in accordance with NMR experiments with ShPrP and mouse PrP, a flexibly disordered N-terminal part, PrP(90-126), and a structured C-terminal part, PrP(127-230), which includes three alpha-helices and a short antiparallel beta-strand. The simulations showed some tendency for the highly conserved hydrophobic segment PrP(112-131) to adopt an alpha-helical conformation and for helix C to split at residues 212-213, a known disease-associated mutation site (Q212P). Three highly occupied salt bridges could be identified (E146/D144<-->R208, R164<-->D178, and R156<-->E196) which appear to be important for the stability of PrP by linking the stable main structured core (helices B and C) with the more flexible structured part (helix A and strands A and B). Two of these salt bridges involve disease-associated mutations (R208H and D178N). Decreased PrP stability shown by protein unfolding experiments on mutants of these residues and guanidinium chloride or temperature-induced unfolding studies indicating reduced stability at low pH are consistent with stabilization by salt bridges. The fact that electrostatic interactions, in general, and salt bridges, in particular, appear to play an important role in PrP stability has implications for PrP structure and stability at different pHs it may encounter physiologically during normal or abnormal recycling from the pH neutral membrane surface into endosomes or lysomes (acidic pHs) or in NMR experiments (5.2 for ShPrP and 4.5 for mouse PrP).     

Characterization of rhodopsin congenital night blindness mutant T94I

The Thr94 --> Ile mutation in the second transmembrane segment of rhodopsin has been reported to be associated with a congenital night blindness phenotype in a large Irish pedigree. Previously, two other known rhodopsin mutants that cause congenital night blindness, A292E and G90D, have been shown in vitro to constitutively activate the G protein transducin in the absence of a chromophore. The proposed mechanism of constitutive activation of these two mutants is an electrostatic disruption of the active site salt bridge between Glu113 and Lys296 that contributes to stabilization of the protein in the inactive state. Here, the T94I rhodopsin mutant is characterized and compared to the two other known rhodopsin night blindness mutants. The T94I mutant opsin is shown also to constitutively activate transducin. The T94I mutant pigment (with a bound 11-cis-retinal chromophore), like the other known rhodopsin night blindness mutants, is not active in the dark and has wild-type activity upon exposure to light. Similar to the Gly90 --> Asp substitution, position 94 is close enough to the Schiff base nitrogen that an Asp at this position can functionally substitute for the Glu113 counterion. However, in contrast to the other night blindness mutants, the T94I MII intermediate decays with a half-life that is approximately 8-fold slower than in the wild-type MII intermediate. Thus, the one phenotype shared by all congenital night blindness mutants that is different from the wild-type protein is constitutive activation of the apoprotein.     

Mutational analysis demonstrates that specific electrostatic interactions can play a key role in the denatured state ensemble of proteins

The nature of the denatured state ensemble has been controversial for decades owing, in large part, to the difficulty in characterizing the structure and energetics of denatured state interactions. There is increasing evidence for relatively non-specific hydrophobic clustering in the denatured states of some proteins but other types of interactions are much less well characterized. Here, we report the characterization of highly specific electrostatic interactions in the denatured state of a small alpha-beta protein, the N-terminal domain of the ribosomal protein L9 (NTL9). Mutation of Lys12 to Met has been shown to increase the stability of NTL9 significantly through the disruption of denatured state interactions. Here, we describe the analysis of the pH-dependent stability of 13 mutants designed to probe the nature of the Lys12 denatured state interaction. Lys12 is located in a lysine-rich region of the protein but analysis of a set of Lys to Met mutants shows that it plays a unique role in the denatured state. Analysis of mutants of all of the acidic residues in NTL9 shows that Lys12 forms a specific non-native electrostatic interaction with Asp8 in the denatured state ensemble. Thus the distribution of charge-charge interactions in the denatured state ensemble of NTL9 appears to be biased by few key interactions and is very different from that expected in a random coil. We propose that these interactions are not encoded by local sequence effects but rather reflect interactions among residues more distant in sequence. These results demonstrate that electrostatic as well as hydrophobic interactions can play an important role in the denatured state ensemble.     

A model for constitutive lutropin receptor activation based on molecular simulation and engineered mutations in transmembrane helices 6 and 7

Many naturally occurring and engineered mutations lead to constitutive activation of the G protein-coupled lutropin receptor (LHR), some of which also result in reduced ligand responsiveness. To elucidate the nature of interhelical interactions in this heptahelical receptor and changes thereof accompanying activation, we have utilized site-directed mutagenesis on transmembrane helices 6 and 7 of rat LHR to prepare and characterize a number of single, double, and triple mutants. The potent constitutively activating mutants, D556(6.44)H and D556(6.44)Q, were combined with weaker activating mutants, N593(7.45)R and N597(7.49)Q, and the loss-of-responsiveness mutant, N593(7.45)A. The engineered mutants have also been simulated using a new receptor model based on the crystal structure of rhodopsin. The results suggest that constitutive LHR activation by mutations at Asp-556(6.44) is triggered by the breakage or weakening of the interaction found in the wild type receptor between Asp-556(6.44) and Asn-593(7.45). Whereas this perturbation is unique to the activating mutations at Asp-556(6.44), common features to all of the most active LHR mutants are the breakage of the charge-reinforced H-bonding interaction between Arg-442(3.50) and Asp-542(6.30) and the increase in solvent accessibility of the cytosolic extensions of helices 3 and 6, which probably participate in the receptor-G protein interface. Asn-593(7.45) and Asn-597(7.49) also seem to be necessary for the high constitutive activities of D556(6.44)H and D556(6.44)Q and for full ligand responsiveness. The new theoretical model provides a foundation for further experimental work on the molecular mechanism(s) of receptor activation.     

Evolutionary trace of G protein-coupled receptors reveals clusters of residues that determine global and class-specific functions

G protein-coupled receptor (GPCR) activation mediated by ligand-induced structural reorganization of its helices is poorly understood. To determine the universal elements of this conformational switch, we used evolutionary tracing (ET) to identify residue positions commonly important in diverse GPCRs. When mapped onto the rhodopsin structure, these trace residues cluster into a network of contacts from the retinal binding site to the G protein-coupling loops. Their roles in a generic transduction mechanism were verified by 211 of 239 published mutations that caused functional defects. When grouped according to the nature of the defects, these residues sub-divided into three striking sub-clusters: a trigger region, where mutations mostly affect ligand binding, a coupling region near the cytoplasmic interface to the G protein, where mutations affect G protein activation, and a linking core in between where mutations cause constitutive activity and other defects. Differential ET analysis of the opsin family revealed an additional set of opsin-specific residues, several of which form part of the retinal binding pocket, and are known to cause functional defects upon mutation. To test the predictive power of ET, we introduced novel mutations in bovine rhodopsin at a globally important position, Leu-79, and at an opsin-specific position, Trp-175. Both were functionally critical, causing constitutive G protein activation of the mutants and rapid loss of regeneration after photobleaching. These results define in GPCRs a canonical signal transduction mechanism where ligand binding induces conformational changes propagated through adjacent trigger, linking core, and coupling regions.     

Theoretical study of the electrostatically driven step of receptor-G protein recognition

This study proposes a theoretical model describing the electrostatically driven step of the alpha 1 b-adrenergic receptor (AR)-G protein recognition. The comparative analysis of the structural-dynamics features of functionally different receptor forms, i.e., the wild type (ground state) and its constitutively active mutants D142A and A293E, was instrumental to gain insight on the receptor-G protein electrostatic and steric complementarity. Rigid body docking simulations between the different forms of the alpha 1 b-AR and the heterotrimeric G alpha q, G alpha s, G alpha i1, and G alpha t suggest that the cytosolic crevice shared by the active receptor and including the second and the third intracellular loops as well as the cytosolic extension of helices 5 and 6, represents the receptor surface with docking complementarity with the G protein. On the other hand, the G protein solvent-exposed portions that recognize the intracellular loops of the activated receptors are the N-terminal portion of alpha 3, alpha G, the alpha G/alpha 4 loop, alpha 4, the alpha 4/beta 6 loop, alpha 5, and the C-terminus. Docking simulations suggest that the two constitutively active mutants D142A and A293E recognize different G proteins with similar selectivity orders, i.e., G alpha q approximately equal to G alpha s > G alpha i > G alpha t. The theoretical models herein proposed might provide useful suggestions for new experiments aiming at exploring the receptor-G protein interface.     

Brownian dynamics of interactions between glyceraldehyde-3-phosphate dehydrogenase (GAPDH) mutants and F-actin

Brownian dynamics simulations of computer models of GAPDH mutants interacting with F-actin emphasized the electrostatic nature of such interactions, and confirmed the importance of four previously identified lysine residues on the GAPDH structure in these interactions. Mutants were GAPDH models in which one or more of the previously identified lysines had been replaced with alanine. Simulations showed reduced binding of these mutants to F-actin compared to wild-type GAPDH. Binding was significantly reduced by mutating the four lysines; the specific electrostatic interaction energy of the quadruple mutant was -7.3 +/- 1.0 compared to -11.4 +/- 0.5 kcal/mol for the wild enzyme. The BD simulations also reaffirmed the importance of quaternary structure for GAPDH binding F-actin.     

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¹⁶⁷ on H4 and Phe²⁰³ 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.