Molecular mechanism of long-range synergetic color tuning between multiple amino acid residues in conger rhodopsin

Structural Factors Determining the Absorption Spectrum of Channelrhodopsins: A Case Study of the Chimera C1C2

Channelrhodopsins are photosensitive proteins that trigger flagella motion in single-cell algae and have been successfully utilized in optogenetic applications. In optogenetics, light is used to activate neural cells in living organisms, which can be achieved by exploiting the ion channel signaling of channelrhodopsins. Tailoring channelrhodopsins for such applications includes the tuning of the absorption maximum. In order to establish rational design and to obtain a desired spectral shift, a basic understanding of the absorption spectrum is required. We have studied the chimera C1C2 as a representative of this protein family and the first member with an available crystal structure. For this purpose, we sampled the conformations of C1C2 using quantum mechanical/molecular mechanical molecular dynamics and subjected the resulting snapshots of the trajectory to excitation energy calculations using ADC(2) and simplified time-dependent density functional theory. In contrast to previous reports, we found that different hydrogen-bonding networks-involving the retinal protonated Schiff base, the putative counterions E162 and D292, and water molecules-had only a small impact on the absorption spectrum. However, in the case of deprotonated E162, increasing the distance to the Schiff base hydrogen-bonding partner led to a systematic blue shift. The β-ionone ring rotation was identified as another important contributor. Yet the most important factors were found to be the bond length alternation and bond order alternation that were linearly correlated to the absorption maximum by up to 62 and 82%, respectively. We ascribe this novel insight into the structural basis of the absorption spectrum to our enhanced protein setup that includes membrane embedding as well as long and extensive sampling.

Color-tuning of natural variants of heliorhodopsin

Microbial rhodopsins are distributed through many microorganisms. Heliorhodopsins are newly discovered but have an unclear function. They have seven transmembrane helices similar to type-I and type-II rhodopsins, but they are different in that the N-terminal region of heliorhodopsin is cytoplasmic. We chose 13 representative heliorhodopsins from various microorganisms, expressed and purified with an N-terminal His tag, and measured the absorption spectra. The 13 natural variants had an absorption maximum (λmax) in the range 530–556 nm similar to proteorhodopsin (λmax = 490–525 nm). We selected several candidate residues that influence rhodopsin color-tuning based on sequence alignment and constructed mutants via site-directed mutagenesis to confirm the spectral changes. We found two important residues located near retinal chromophore that influence λmax. We also predict the 3D structure via homology-modeling of Thermoplasmatales heliorhodopsin. The results indicate that the color-tuning mechanism of type-I rhodopsin can be applied to understand the color-tuning of heliorhodopsin.

Normal mode analysis and beyond

Normal mode analysis provides a powerful tool in biophysical computations. Particularly, we shed light on its application to protein properties because they directly lead to biological functions. As a result of normal mode analysis, the protein motion is represented as a linear combination of mutually independent normal mode vectors. It has been widely accepted that the large amplitude motions throughout the entire protein molecule can be well described with a few low-frequency normal modes. Furthermore, it is possible to represent the effect of external perturbations, e.g., ligand binding, hydrostatic pressure, as the shifts of normal mode variables. Making use of this advantage, we are able to explore mechanical properties of proteins such as Young's modulus and compressibility. Within thermally fluctuating protein molecules under physiological conditions, tightly packed amino acid residues interact with each other through heat and energy exchanges. Since the structure and dynamics of protein molecules are highly anisotropic, the flow of energy and heat should also be anisotropic. Based on the harmonic approximation of the heat current operator, it is possible to analyze the communication map of a protein molecule. By using this method, the energy transfer pathways of photoactive yellow protein were calculated. It turned out that these pathways are similar to those obtained via the Green-Kubo formalism with equilibrium molecular dynamics simulations, indicating that normal mode analysis captures the intrinsic nature of the transport properties of proteins.

Rhodopsin gene copies in Japanese eel originated in a teleost-specific genome duplication

BACKGROUND

Gene duplication is considered important to increasing the genetic diversity in animals. In fish, visual pigment genes are often independently duplicated, and the evolutionary significance of such duplications has long been of interest. Eels have two rhodopsin genes (rho), one of which (freshwater type, fw-rho) functions in freshwater and the other (deep-sea type, ds-rho) in marine environments. Hence, switching of rho expression in retinal cells is tightly linked with eels' unique life cycle, in which they migrate from rivers or lakes to the sea. These rho genes are apparently paralogous, but the timing of their duplication is unclear due to the deep-branching phylogeny. The aim of the present study is to elucidate the evolutionary origin of the two rho copies in eels using comparative genomics methods.

RESULTS

In the present study, we sequenced the genome of Japanese eel Anguilla japonica and reconstructed two regions containing rho by de novo assembly. We found a single corresponding region in a non-teleostean primitive ray-finned fish (spotted gar) and two regions in a primitive teleost (Asian arowana). The order of ds-rho and the neighboring genes was highly conserved among the three species. With respect to fw-rho, which was lost in Asian arowana, the neighboring genes were also syntenic between Japanese eel and Asian arowana. In particular, the pattern of gene losses in ds-rho and fw-rho regions was the same as that in Asian arowana, and no discrepancy was found in any of the teleost genomes examined. Phylogenetic analysis supports mutual monophyly of these two teleostean synteny groups, which correspond to the ds-rho and fw-rho regions.

CONCLUSIONS

Syntenic and phylogenetic analyses suggest that the duplication of rhodopsin gene in Japanese eel predated the divergence of eel (Elopomorpha) and arowana (Osteoglossomorpha). Thus, based on the principle of parsimony, it is most likely that the rhodopsin paralogs were generated through a whole genome duplication in the ancestor of teleosts, and have remained till the present in eels with distinct functional roles. Our result indicates, for the first time, that teleost-specific genome duplication may have contributed to a gene innovation involved in eel-specific migratory life cycle.

Atomistic design of microbial opsin-based blue-shifted optogenetics tools

Microbial opsins with a bound chromophore function as photosensitive ion transporters and have been employed in optogenetics for the optical control of neuronal activity. Molecular engineering has been utilized to create colour variants for the functional augmentation of optogenetics tools, but was limited by the complexity of the protein-chromophore interactions. Here we report the development of blue-shifted colour variants by rational design at atomic resolution, achieved through accurate hybrid molecular simulations, electrophysiology and X-ray crystallography. The molecular simulation models and the crystal structure reveal the precisely designed conformational changes of the chromophore induced by combinatory mutations that shrink its π-conjugated system which, together with electrostatic tuning, produce large blue shifts of the absorption spectra by maximally 100 nm, while maintaining photosensitive ion transport activities. The design principle we elaborate is applicable to other microbial opsins, and clarifies the underlying molecular mechanism of the blue-shifted action spectra of microbial opsins recently isolated from natural sources.

Chimeric proton-pumping rhodopsins containing the cytoplasmic loop of bovine rhodopsin

G-protein-coupled receptors (GPCRs) transmit stimuli to intracellular signaling systems. Rhodopsin (Rh), which is a prototypical GPCR, possesses an 11-cis retinal. Photoisomerization of 11-cis to all-trans leads to structural changes in the protein of cytoplasmic loops, activating G-protein. Microbial rhodopsins are similar heptahelical membrane proteins that function as bacterial sensors, light-driven ion-pumps, or light-gated channels. They possess an all-trans retinal, and photoisomerization to 13-cis triggers structural changes in protein. Despite these similarities, there is no sequence homology between visual and microbial rhodopsins, and microbial rhodopsins do not activate G-proteins. In this study, new chimeric proton-pumping rhodopsins, proteorhodopsin (PR) and Gloeobacter rhodopsin (GR) were designed by replacing cytoplasmic loops with bovine Rh loops. Although G-protein was not activated by the PR chimeras, all 12 GR chimeras activated G-protein. The GR chimera containing the second cytoplasmic loop of bovine Rh did not activate G-protein. However, the chimera with a second and third double-loop further enhanced G-protein activation. Introduction of an E132Q mutation slowed the photocycle 30-fold and enhanced activation. The highest catalytic activity of the GR chimera was still 3,200 times lower than bovine Rh but only 64 times lower than amphioxus Go-rhodopsin. This GR chimera showed a strong absorption change of the amide-I band on a light-minus-dark difference FTIR spectrum which could represent a larger helical opening, important for G-protein activation. The light-dependent catalytic activity of this GR chimera makes it a potential optogenetic tool for enzymatic activation by light.

Rhodopsin Absorption from First Principles: Bypassing Common Pitfalls

Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.

Gas-Phase Retinal Spectroscopy: Temperature Effects Are But a Mirage

We employ state-of-the-art first-principle approaches to investigate whether temperature effects are responsible for the unusually broad and flat spectrum of protonated Schiff base retinal observed in photodissociation spectroscopy, as has recently been proposed. We first carefully calibrate how to construct a realistic geometrical model of retinal and show that the exchange-correlation M06-2X functional yields an accurate description while the commonly used complete active space self-consistent field method (CASSCF) is not adequate. Using modern multiconfigurational perturbative methods (NEVPT2) to compute the excitations, we then demonstrate that conformations with different orientations of the β-ionone ring are characterized by similar excitations. Moreover, other degrees of freedom identified as active in room-temperature molecular dynamics simulations do not yield the shift required to explain the anomalous spectral shape. Our findings indicate that photodissociation experiments are not representative of the optical spectrum of retinal in the gas phase and call for further experimental characterization of the dissociation spectra.

Structural model of channelrhodopsin

Channelrhodopsins (ChRs) are light-gated cation channels that mediate ion transport across membranes in microalgae (vectorial catalysis). ChRs are now widely used for the analysis of neural networks in tissues and living animals with light (optogenetics). For elucidation of functional mechanisms at the atomic level, as well as for further engineering and application, a detailed structure is urgently needed. In the absence of an experimental structure, here we develop a structural ChR model based on several molecular computational approaches, capitalizing on characteristic patterns in amino acid sequences of ChR1, ChR2, Volvox ChRs, Mesostigma ChR, and the recently identified ChR of the halophilic alga Dunaliella salina. In the present model, we identify remarkable structural motifs that may explain fundamental electrophysiological properties of ChR2, ChR1, and their mutants, and in a crucial validation of the model, we successfully reproduce the excitation energy predicted by absorption spectra.

Color tuning in photofunctional proteins

Depending on protein environment, a single photofunctional chromophore shows a wide variation of photoabsorption/emission energies. This photobiological phenomenon, known as color tuning, is observed in human visual cone pigments, firefly luciferase, and red fluorescent protein. We investigate the origin of color tuning by quantum chemical calculations on the excited states: symmetry-adapted cluster-configuration interaction (SAC-CI) method for excited states and a combined quantum mechanical (QM)/molecular mechanical (MM) method for protein environments. This Minireview summarizes our theoretical studies on the above three systems and explains a common feature of their color-tuning mechanisms. It also discuss the possibility of artificial color tuning toward a rational design of photoabsorption/emission properties.