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

Coordination of bilayer properties by an inward-rectifier K+ channel is a cooperative process driven by protein-lipid interaction

Physical properties of biological membranes directly or indirectly govern biological processes. Yet, the interplay between membrane and integral membrane proteins is difficult to assess due to reciprocal effects between membrane proteins, individual lipids, and membrane architecture. Using solid-state NMR (SSNMR) we previously showed that KirBac1.1, a bacterial Inward-Rectifier K+ channel, nucleates bilayer ordering and microdomain formation through tethering anionic lipids. Conversely, these lipids cooperatively bind cationic residues to activate the channel and initiate K+ flux. The mechanistic details governing the relationship between cooperative lipid loading and bilayer ordering are, however, unknown. To investigate, we generated KirBac1.1 samples with different concentrations of 13C-lableded phosphatidyl glycerol (PG) lipids and acquired a full suite of SSNMR 1D temperature series experiments using the ordered all-trans (AT) and disordered trans-gauche (TG) acyl conformations as markers of bilayer dynamics. We observed increased AT ordered signal, decreased TG disordered signal, and increased bilayer melting temperature with increased PG concentration. Further, we identified cooperativity between ordering and direct binding of PG lipids, indicating KirBac1.1-driven bilayer ordering and microdomain formation is a classically cooperative Hill-type process driven by and predicated upon direct binding of PG lipids. Our results provide unique mechanistic insight into how proteins and lipids in tandem contribute to supramolecular bilayer heterogeneity in the lipid membrane.

Strategies for acquisition of resonance assignment spectra of highly dynamic membrane proteins: a GPCR case study

In protein nuclear magnetic resonance (NMR), chemical shift assignment provides a wealth of information. However, acquisition of high-quality solid-state NMR spectra depends on protein-specific dynamics. For membrane proteins, bilayer heterogeneity further complicates this observation. Since the efficiency of cross-polarization transfer is strongly entwined with protein dynamics, optimal temperatures for spectral sensitivity and resolution will depend not only on inherent protein dynamics, but temperature-dependent phase properties of the bilayer environment. We acquired 1-, 2-, and 3D homo- and heteronuclear experiments of the chemokine receptor CCR3 in a 7:3 phosphatidylcholine:cholesterol lipid environment. 1D direct polarization, cross polarization (CP), and T2' experiments indicate sample temperatures below - 25 °C facilitate higher CP enhancement and longer-lived transverse relaxation times. T1rho experiments indicate intermediate timescales are minimized below a sample temperature of - 20 °C. 2D DCP NCA experiments indicated optimal CP efficiency and resolution at a sample temperature of - 30 °C, corroborated by linewidth analysis in 3D NCACX at - 30 °C compared to - 5 °C. This optimal temperature is concluded to be directly related the lipid phase transition, measured to be between - 20 and 15 °C based on rINEPT signal of all-trans and trans-gauche lipid acyl conformations. Our results have critical implications in acquisition of SSNMR membrane protein assignment spectra, as we hypothesize that different lipid compositions with different phase transition properties influence protein dynamics and therefore the optimal acquisition temperature.

Cholesterol Biases the Conformational Landscape of the Chemokine Receptor CCR3: A MAS SSNMR-Filtered Molecular Dynamics Study

Cholesterol directs the pathway of ligand-induced G protein-coupled receptor (GPCR) signal transduction. The GPCR C-C motif chemokine receptor 3 (CCR3) is the principal chemotactic receptor for eosinophils, with roles in cancer metastasis and autoinflammatory conditions. Recently, we discovered a direct correlation between bilayer cholesterol and increased agonist-triggered CCR3 signal transduction. However, the allosteric molecular mechanism escalating ligand affinity and G protein coupling is unknown. To study cholesterol-guided CCR3 conformational selection, we implement comparative, objective measurement of protein architectures by scoring shifts (COMPASS) to grade model structures from molecular dynamics simulations. In this workflow, we scored predicted chemical shifts against 2-dimensional solid-state NMR 13C-13C correlation spectra of U-15N,13C-CCR3 samples prepared with and without cholesterol. Our analysis of trajectory model structures uncovers that cholesterol induces site-specific conformational restraint of extracellular loop (ECL) 2 and conserved motion in transmembrane helices and ECL3 not observed in simulations of bilayers with only phosphatidylcholine lipids. PyLipID analysis implicates direct cholesterol agency in CCR3 conformational selection and dynamics. Residue-residue contact scoring shows that cholesterol biases the conformational selection of the orthosteric pocket involving Y411.39, Y1133.32, and E2877.39. Lastly, we observe contact remodeling in activation pathway residues centered on the initial transmission switch, Na+ pocket, and R3.50 in the DRY motif. Our observations have unique implications for understanding of CCR3 ligand recognition and specificity and provide mechanistic insight into how cholesterol functions as an allosteric regulator of CCR3 signal transduction.

In-Cell Solid-State NMR Analysis of Membrane Proteins

Collapse

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.

Specificity of the chromophore-binding site in human cone opsins

The variable composition of the chromophore-binding pocket in visual receptors is essential for vision. The visual phototransduction starts with the cis-trans isomerization of the retinal chromophore upon absorption of photons. Despite sharing the common 11-cis-retinal chromophore, rod and cone photoreceptors possess distinct photochemical properties. Thus, a detailed molecular characterization of the chromophore-binding pocket of these receptors is critical to understanding the differences in the photochemistry of vision between rods and cones. Unlike for rhodopsin (Rh), the crystal structures of cone opsins remain to be determined. To obtain insights into the specific chromophore-protein interactions that govern spectral tuning in human visual pigments, here we harnessed the unique binding properties of 11-cis-6-membered-ring-retinal (11-cis-6mr-retinal) with human blue, green, and red cone opsins. To unravel the specificity of the chromophore-binding pocket of cone opsins, we applied 11-cis-6mr-retinal analog-binding analyses to human blue, green, and red cone opsins. Our results revealed that among the three cone opsins, only blue cone opsin can accommodate the 11-cis-6mr-retinal in its chromophore-binding pocket, resulting in the formation of a synthetic blue pigment (B6mr) that absorbs visible light. A combination of primary sequence alignment, molecular modeling, and mutagenesis experiments revealed the specific amino acid residue 6.48 (Tyr-262 in blue cone opsins and Trp-281 in green and red cone opsins) as a selectivity filter in human cone opsins. Altogether, the results of our study uncover the molecular basis underlying the binding selectivity of 11-cis-6mr-retinal to the cone opsins.

Synthesis of 9-CD3-9-cis-Retinal Cofactor of Isorhodopsin

We report the synthesis of 9-CD₃-9-cis-retinal via a six-step procedure from β-ionone. The steps involve an initial deuteration of the methyl ketone of β-ionone followed by two consecutive Horner-Wadsworth-Emmons (HWE) coupling reactions and their corresponding DIBAL reductions. A final oxidation of the allylic alcohol of the retinol leads to the target compound. This deuterium labeled retinoid is an important cofactor for studying protein-retinoid interactions in isorhodopsin.

Recent advances in biophysical studies of rhodopsins - Oligomerization, folding, and structure

Retinal-binding proteins, mainly known as rhodopsins, function as photosensors and ion transporters in a wide range of organisms. From halobacterial light-driven proton pump, bacteriorhodopsin, to bovine photoreceptor, visual rhodopsin, they have served as prototypical α-helical membrane proteins in a large number of biophysical studies and aided in the development of many cutting-edge techniques of structural biology and biospectroscopy. In the last decade, microbial and animal rhodopsin families have expanded significantly, bringing into play a number of new interesting structures and functions. In this review, we will discuss recent advances in biophysical approaches to retinal-binding proteins, primarily microbial rhodopsins, including those in optical spectroscopy, X-ray crystallography, nuclear magnetic resonance, and electron paramagnetic resonance, as applied to such fundamental biological aspects as protein oligomerization, folding, and structure.

Nuclear magnetic resonance (NMR) applied to membrane–protein complexes

Increasing evidence suggests that most proteins occur and function in complexes rather than as isolated entities when embedded in cellular membranes. Nuclear magnetic resonance (NMR) provides increasing possibilities to study structure, dynamics and assembly of such systems. In our review, we discuss recent methodological progress to study membrane–protein complexes (MPCs) by NMR, starting with expression, isotope-labeling and reconstitution protocols. We review approaches to deal with spectral complexity and limited spectral spectroscopic sensitivity that are usually encountered in NMR-based studies of MPCs. We highlight NMR applications in various classes of MPCs, including G-protein-coupled receptors, ion channels and retinal proteins and extend our discussion to protein–protein complexes that span entire cellular compartments or orchestrate processes such as protein transport across or within membranes. These examples demonstrate the growing potential of NMR-based studies of MPCs to provide critical insight into the energetics of protein–ligand and protein–protein interactions that underlie essential biological functions in cellular membranes.

Proteorhodopsin Activation Is Modulated by Dynamic Changes in Internal Hydration

Proteorhodopsin, a member of the microbial rhodopsin family, is a seven-transmembrane α-helical protein that functions as a light-driven proton pump. Understanding the proton-pumping mechanism of proteorhodopsin requires intimate knowledge of the proton transfer pathway via complex hydrogen-bonding networks formed by amino acid residues and internal water molecules. Here we conducted a series of microsecond time scale molecular dynamics simulations on both the dark state and the initial photoactivated state of blue proteorhodopsin to reveal the structural basis for proton transfer with respect to protein internal hydration. A complex series of dynamic hydrogen-bonding networks involving water molecules exists, facilitated by water channels and hydration sites within proteorhodopsin. High levels of hydration were discovered at each proton transfer site-the retinal binding pocket and proton uptake and release sites-underscoring the critical participation of water molecules in the proton-pumping mechanism. Water-bridged interactions and local water channels were also observed and can potentially mediate long-distance proton transfer between each site. The most significant phenomenon is after isomerization of retinal, an increase in water flux occurs that connects the proton release group, a conserved arginine residue, and the retinal binding pocket. Our results provide a detailed description of the internal hydration of the early photointermediates in the proteorhodopsin photocycle under alkaline pH conditions. These results lay the fundamental groundwork for understanding the intimate role that hydration plays in the structure-function relationship underlying the proteorhodopsin proton-pumping mechanism, as well as providing context for the relationship of hydration in proteorhodopsin to other microbial retinal proteins.

The G protein-coupled receptor rhodopsin: a historical perspective

Rhodopsin is a key light-sensitive protein expressed exclusively in rod photoreceptor cells of the retina. Failure to express this transmembrane protein causes a lack of rod outer segment formation and progressive retinal degeneration, including the loss of cone photoreceptor cells. Molecular studies of rhodopsin have paved the way to understanding a large family of cell-surface membrane proteins called G protein-coupled receptors (GPCRs). Work started on rhodopsin over 100 years ago still continues today with substantial progress made every year. These activities underscore the importance of rhodopsin as a prototypical GPCR and receptor required for visual perception-the fundamental process of translating light energy into a biochemical cascade of events culminating in vision.

Investigation of rhodopsin dynamics in its signaling state by solid-state deuterium NMR spectroscopy

Site-directed deuterium NMR spectroscopy is a valuable tool to study the structural dynamics of biomolecules in cases where solution NMR is inapplicable. Solid-state (2)H NMR spectral studies of aligned membrane samples of rhodopsin with selectively labeled retinal provide information on structural changes of the chromophore in different protein states. Moreover (2)H NMR relaxation time measurements allow one to study the dynamics of the ligand during the transition from the inactive to the active state. Here we describe the methodological aspects of solid-state (2)H NMR spectroscopy for functional studies of rhodopsin, with an emphasis on the dynamics of the retinal cofactor. We provide complete protocols for the preparation of NMR samples of rhodopsin with 11-cis-retinal selectively deuterated at the methyl groups in aligned membranes. In addition we review optimized conditions for trapping the rhodopsin photointermediates; and we address the challenging problem of trapping the signaling state of rhodopsin in aligned membrane films.

Structure-based simulations reveal concerted dynamics of GPCR activation

G protein-coupled receptors (GPCRs) are a vital class of proteins that transduce biological signals across the cell membrane. However, their allosteric activation mechanism is not fully understood; crystal structures of active and inactive receptors have been reported, but the functional pathway between these two states remains elusive. Here, we use structure-based (Gō-like) models to simulate activation of two GPCRs, rhodopsin and the β₂ adrenergic receptor (β₂AR). We used data-derived reaction coordinates that capture the activation mechanism for both proteins, showing that activation proceeds through quantitatively different paths in the two systems. Both reaction coordinates are determined from the dominant concerted motions in the simulations so the technique is broadly applicable. There were two surprising results. First, the main structural changes in the simulations were distributed throughout the transmembrane bundle, and not localized to the obvious areas of interest, such as the intracellular portion of Helix 6. Second, the activation (and deactivation) paths were distinctly nonmonotonic, populating states that were not simply interpolations between the inactive and active structures. These transitions also suggest a functional explanation for β₂AR's basal activity: it can proceed through a more broadly defined path during the observed transitions.

Interactions between fengycin and model bilayers quantified by coarse-grained molecular dynamics

Bacteria, particularly of the genus Bacillus, produce a wide variety of antifungal compounds. They act by affecting the lipid bilayers of fungal membranes, causing curvature-induced strain and eventual permeabilization. One class of these, known as fengycins, has been commercialized for treating agricultural infections and shows some promise as a possible antifungal pharmaceutical. Understanding the mechanism by which fengycins damage lipid bilayers could prove useful to the future development of related antifungal treatments. In this work, we present multi-microsecond-long simulations of fengycin interacting with different lipid bilayer systems. We see fengycin aggregation and uncover a clear aggregation pattern that is partially influenced by bilayer composition. We also quantify some local bilayer perturbations caused by fengycin binding, including curvature of the lipid bilayer and local electrostatic-driven reorganization.

Amino acid conservation and interactions in rhodopsin: probing receptor activation by NMR spectroscopy

Rhodopsin is a classical two-state G protein-coupled receptor (GPCR). In the dark, its 11-cis retinal chromophore serves as an inverse agonist to lock the receptor in an inactive state. Retinal-protein and protein-protein interactions have evolved to reduce the basal activity of the receptor in order to achieve low dark noise in the visual system. In contrast, absorption of light triggers rapid isomerization of the retinal, which drives the conversion of the receptor to a fully active conformation. Several specific protein-protein interactions have evolved that maintain the lifetime of the active state in order to increase the sensitivity of this receptor for dim-light vision in vertebrates. In this article, we review the molecular interactions that stabilize rhodopsin in the dark-state and describe the use of solid-state NMR spectroscopy for probing the structural changes that occur upon light-activation. Amino acid conservation provides a guide for those interactions that are common in the class A GPCRs as well as those that are unique to the visual system. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

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.

Conformational flexibility and structural dynamics in GPCR-mediated G protein activation: a perspective

Structure and dynamics of G proteins and their cognate receptors, both alone and in complex, are becoming increasingly accessible to experimental techniques. Understanding the conformational changes and timelines that govern these changes can lead to new insights into the processes of ligand binding and associated G protein activation. Experimental systems may involve the use of, or otherwise stabilize, non-native environments. This can complicate our understanding of structural and dynamic features of processes such as the ionic lock, tryptophan toggle, and G protein flexibility. While elements in the receptor's transmembrane helices and the C-terminal α5 helix of Gα undergo well-defined structural changes, regions subject to conformational flexibility may be important in fine-tuning the interactions between activated receptors and G proteins. The pairing of computational and experimental approaches will continue to provide powerful tools to probe the conformation and dynamics of receptor-mediated G protein activation.

G-protein-coupled receptor structure, ligand binding and activation as studied by solid-state NMR spectroscopy

GPCRs (G-protein-coupled receptors) are versatile signalling molecules at the cell surface and make up the largest and most diverse family of membrane receptors in the human genome. They convert a large variety of extracellular stimuli into intracellular responses through the activation of heterotrimeric G-proteins, which make them key regulatory elements in a broad range of normal and pathological processes, and are therefore one of the most important targets for pharmaceutical drug discovery. Knowledge of a GPCR structure enables us to gain a mechanistic insight into its function and dynamics, and further aid rational drug design. Despite intensive research carried out over the last three decades, resolving the structural basis of GPCR function is still a major activity. The crystal structures obtained in the last 5 years provide the first opportunity to understand how protein structure dictates the unique functional properties of these complex signalling molecules. However, owing to the intrinsic hydrophobicity, flexibility and instability of membrane proteins, it is still a challenge to crystallize GPCRs, and, when this is possible, it is no longer in its native membrane environment and no longer without modification. Furthermore, the conformational change of the transmembrane α-helices associated with the structure activation increases the difficulty of capturing the activation state of a GPCR to a higher resolution by X-ray crystallography. On the other hand, solid-state NMR may offer a unique opportunity to study membrane protein structure, ligand binding and activation at atomic resolution in the native membrane environment, as well as described functionally significant dynamics. In the present review, we discuss some recent achievements of solid-state NMR for understanding GPCRs, the largest mammalian proteome at ~1% of the total expressed proteins. Structural information, details of determination, details of ligand conformations and the consequences of ligand binding to initiate activation can all be explored with solid-state NMR.

Solid-state NMR spectroscopy of proteins

Solid-state NMR spectroscopy proved to be a versatile tool for characterization of structure and dynamics of complex biochemical systems. In particular, magic angle spinning (MAS) solid-state NMR came to maturity for application towards structural elucidation of biological macromolecules. Current challenges in applying solid-state NMR as well as progress achieved recently will be discussed in the following chapter focusing on conceptual aspects important for structural elucidation of proteins.

Ligand-dependent conformations and dynamics of the serotonin 5-HT(2A) receptor determine its activation and membrane-driven oligomerization properties

From computational simulations of a serotonin 2A receptor (5-HT_2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT_2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization.

The 5-HT_2A receptor for the neurotransmitter serotonin (5-HT) belongs to family A (rhodopsin-like) G-protein coupled receptors (GPCRs), one of the most important classes of membrane proteins that are targeted by an extensive and diverse collection of external stimuli. Recently we learned that different ligands targeting the same GPCR can elicit different biological responses, but the mechanisms remain unknown. We address this fundamental question for the serotonin 5-HT_2A receptor, because it is known to respond to the binding of structurally diverse ligands by producing similar stimuli in the cell, and to the binding of quite similar ligands with dramatically different responses. Molecular dynamics simulations of molecular models of the serotonin 5-HT_2A receptor in complex with pharmacologically distinct ligands show the dynamic rearrangements of the receptor molecule to be different for these ligands, and the nature and extents of the rearrangements reflect the known pharmacological properties of the ligands as full, partial or inverse activators of the receptor. The different rearrangements of the receptor molecule are shown to produce different rearrangements of the surrounding membrane, a remodeling of the environment that can have differential ligand-determined effects on receptor function and association in the cell's membrane.