Steric and electronic influences on the torsional energy landscape of retinal

Lipid membrane mimetics and oligomerization tune functional properties of proteorhodopsin

The functional properties of proteorhodopsin (PR) have been found to be strongly modulated by oligomeric distributions and lipid membrane mimetics. This study aims to distinguish and explain their effects by investigating how oligomer formation impacts PR's function of proton transport in lipid-based membrane mimetic environments. We find that PR forms stable hexamers and pentamers in both E. coli membranes and synthetic liposomes. Compared with the monomers, the photocycle kinetics of PR oligomers is ∼2 and ∼4.5 times slower for transitions between the K and M and the M and N photointermediates, respectively, indicating that oligomerization significantly slows PR's rate of proton transport in liposomes. In contrast, the apparent pKa of the key proton acceptor residue D97 (pKaD97) of liposome-embedded PR persists at 6.2-6.6, regardless of cross-protomer modulation of D97, suggesting that the liposome environment helps maintain PR's functional activity at neutral pH. By comparison, when extracted directly from E. coli membranes into styrene-maleic acid lipid particles, the pKaD97 of monomer-enriched E50Q PR drastically increases to 8.9, implying that there is a very low active PR population at neutral pH to engage in PR's photocycle. These findings demonstrate that oligomerization impacts PR's photocycle kinetics, while lipid-based membrane mimetics strongly affect PR's active population via different mechanisms.

Structural insights into light-driven anion pumping in cyanobacteria

Transmembrane ion transport is a key process in living cells. Active transport of ions is carried out by various ion transporters including microbial rhodopsins (MRs). MRs perform diverse functions such as active and passive ion transport, photo-sensing, and others. In particular, MRs can pump various monovalent ions like Na⁺, K⁺, Cl^-, I^-, NO₃^-. The only characterized MR proposed to pump sulfate in addition to halides belongs to the cyanobacterium Synechocystis sp. PCC 7509 and is named Synechocystis halorhodopsin (SyHR). The structural study of SyHR may help to understand what makes an MR pump divalent ions. Here we present the crystal structure of SyHR in the ground state, the structure of its sulfate-bound form as well as two photoreaction intermediates, the K and O states. These data reveal the molecular origin of the unique properties of the protein (exceptionally strong chloride binding and proposed pumping of divalent anions) and sheds light on the mechanism of anion release and uptake in cyanobacterial halorhodopsins. The unique properties of SyHR highlight its potential as an optogenetics tool and may help engineer different types of anion pumps with applications in optogenetics.

Driving torsion scans with wavefront propagation

The parameterization of torsional/dihedral angle potential energy terms is a crucial part of developing molecular mechanics force fields. Quantum mechanical (QM) methods are often used to provide samples of the potential energy surface (PES) for fitting the empirical parameters in these force field terms. To ensure that the sampled molecular configurations are thermodynamically feasible, constrained QM geometry optimizations are typically carried out, which relax the orthogonal degrees of freedom while fixing the target torsion angle(s) on a grid of values. However, the quality of results and computational cost are affected by various factors on a non-trivial PES, such as dependence on the chosen scan direction and the lack of efficient approaches to integrate results started from multiple initial guesses. In this paper, we propose a systematic and versatile workflow called TorsionDrive to generate energy-minimized structures on a grid of torsion constraints by means of a recursive wavefront propagation algorithm, which resolves the deficiencies of conventional scanning approaches and generates higher quality QM data for force field development. The capabilities of our method are presented for multi-dimensional scans and multiple initial guess structures, and an integration with the MolSSI QCArchive distributed computing ecosystem is described. The method is implemented in an open-source software package that is compatible with many QM software packages and energy minimization codes.

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.

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.

Allosteric Effects of the Proton Donor on the Microbial Proton Pump Proteorhodopsin

Proteorhodopsin (PR) is a microbial proton pump that is ubiquitous in marine environments and may play an important role in the oceanic carbon cycle. Photoisomerization of the retinal chromophore in PR leads to a series of proton transfers between specific acidic amino acid residues and the Schiff base of retinal, culminating in a proton motive force to facilitate ATP synthesis. The proton donor in a similar retinal protein, bacteriorhodopsin, acts as a latch to allow the influx of bulk water. However, it is unclear if the proton donor in PR, E108, utilizes the same latch mechanism to become internally hydrated. Here, we used molecular dynamics simulations to model the changes in internal hydration of the blue variant of PR during photoactivation with the proton donor in protonated and deprotonated states. We find that there is a stark contrast in the levels of internal hydration of the cytoplasmic half of PR based on the protonation state of E108. Instead of a latch mechanism, deprotonation of E108 acts as a gate, taking advantage of a nearby polar residue (S61) to promote the formation of a stable water wire from bulk cytoplasm to the retinal-binding pocket over hundreds of nanoseconds. No large-scale conformational changes occur in PR over the microsecond timescale. This subtle yet clear difference in the effect of deprotonation of the proton donor in PR may help explain why the photointermediates that involve the proton donor (i.e., M and N states) have timescales that are orders of magnitude different from the archaeal proton pump, bacteriorhodopsin. In general, our study highlights the importance of understanding how structural fluctuations lead to differences in the way that retinal proteins accomplish the same task.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

The Activation Pathway of Human Rhodopsin in Comparison to Bovine Rhodopsin

Rhodopsin, the photoreceptor of rod cells, absorbs light to mediate the first step of vision by activating the G protein transducin (Gt). Several human diseases, such as retinitis pigmentosa or congenital night blindness, are linked to rhodopsin malfunctions. Most of the corresponding in vivo studies and structure-function analyses (e.g. based on protein x-ray crystallography or spectroscopy) have been carried out on murine or bovine rhodopsin. Because these rhodopsins differ at several amino acid positions from human rhodopsin, we conducted a comprehensive spectroscopic characterization of human rhodopsin in combination with molecular dynamics simulations. We show by FTIR and UV-visible difference spectroscopy that the light-induced transformations of the early photointermediates are very similar. Significant differences between the pigments appear with formation of the still inactive Meta I state and the transition to active Meta II. However, the conformation of Meta II and its activity toward the G protein are essentially the same, presumably reflecting the evolutionary pressure under which the active state has developed. Altogether, our results show that although the basic activation pathways of human and bovine rhodopsin are similar, structural deviations exist in the inactive conformation and during receptor activation, even between closely related rhodopsins. These differences between the well studied bovine or murine rhodopsins and human rhodopsin have to be taken into account when the influence of point mutations on the activation pathway of human rhodopsin are investigated using the bovine or murine rhodopsin template sequences.

Explaining the mobility of retinal in activated rhodopsin and opsin

Computational studies reveal flexibility of the rhodopsin cofactor, retinal, within the protein binding pocket that play a key role in the activated state and regeneration of rhodopsin.

Retinal ligand mobility explains internal hydration and reconciles active rhodopsin structures

Rhodopsin, the mammalian dim-light receptor, is one of the best-characterized G-protein-coupled receptors, a pharmaceutically important class of membrane proteins that has garnered a great deal of attention because of the recent availability of structural information. Yet the mechanism of rhodopsin activation is not fully understood. Here, we use microsecond-scale all-atom molecular dynamics simulations, validated by solid-state (2)H nuclear magnetic resonance spectroscopy, to understand the transition between the dark and metarhodopsin I (Meta I) states. Our analysis of these simulations reveals striking differences in ligand flexibility between the two states. Retinal is much more dynamic in Meta I, adopting an elongated conformation similar to that seen in the recent activelike crystal structures. Surprisingly, this elongation corresponds to both a dramatic influx of bulk water into the hydrophobic core of the protein and a concerted transition in the highly conserved Trp265(6.48) residue. In addition, enhanced ligand flexibility upon light activation provides an explanation for the different retinal orientations observed in X-ray crystal structures of active rhodopsin.

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.

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.

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.