Rhodopsins

We report a multi-scale simulation study of the photocycle of the rhodopsins. The quasi-atomistic representation ("united atoms" UA) of retinal is combined with a minimalist coarse grained (CG, one-bead-per amino acid) representation of the protein, in a hybrid UA/CG approach, which is the homolog of QM/MM, but at lower resolution. An accurate multi-stable parameterization of the model allows simulating each state and transition among them, and the combination of different scale representation allows addressing the entire photocycle. We test the model on bacterial rhodopsin, for which more experimental data are available, and then also report results for mammalian rhodopsins. In particular, the analysis of simulations reveals the spontaneous appearance of meta-stable states in quantitative agreement with experimental data.

Genetically encoded voltage indicators (GEVIs) have found wide applications as molecular tools for visualization of changes in cell membrane potential. Among others, several classes of archaerhodopsin-3-based GEVIs have been developed and have proved themselves promising in various molecular imaging studies. To expand the application range for this type of GEVIs, new variants with absorption band maxima shifted toward the first biological window and enhanced fluorescence signal are required. Here, we integrate computational and experimental strategies to reveal structural factors that distinguish far-red bright archaerhodopsin-3-based GEVIs, Archers, obtained by directed evolution in a previous study (McIsaac et al., PNAS, 2014) and the wild-type archaerhodopsin-3 with an extremely dim fluorescence signal, aiming to use the obtained information in subsequent rational design. We found that the fluorescence can be enhanced by stabilization of a certain conformation of the protein, which, in turn, can be achieved by tuning the pKa value of two titratable residues. These findings were supported further by introducing mutations into wild-type archeorhodopsin-3 and detecting the enhancement of the fluorescence signal. Finally, we came up with a rational design and proposed previously unknown Archers variants with red-shifted absorption bands (λmax up to 640 nm) and potential-dependent bright fluorescence (quantum yield up to 0.97%).

Transitions between different conformational states are ubiquitous in proteins, being involved in signaling, catalysis, and other fundamental activities in cells. However, modeling those processes is extremely difficult, due to the need of efficiently exploring a vast conformational space in order to seek for the actual transition path for systems whose complexity is already high in the stable states. Here we report a strategy that simplifies this task attacking the complexity on several sides. We first apply a minimalist coarse-grained model to Calmodulin, based on an empirical force field with a partial structural bias, to explore the transition paths between the apo-closed state and the Ca-bound open state of the protein. We then select representative structures along the trajectory based on a structural clustering algorithm and build a cleaned-up trajectory with them. We finally compare this trajectory with that produced by the online tool MinActionPath, by minimizing the action integral using a harmonic network model, and with that obtained by the PROMPT morphing method, based on an optimal mass transportation-type approach including physical constraints. The comparison is performed both on the structural and energetic level, using the coarse-grained and the atomistic force fields upon reconstruction. Our analysis indicates that this method returns trajectories capable of exploring intermediate states with physical meaning, retaining a very low computational cost, which can allow systematic and extensive exploration of the multi-stable proteins transition pathways.

Transitions between different conformational states are ubiquitous in proteins, being involved in signaling, catalysis, and other fundamental activities in cells. However, modeling those processes is extremely difficult, due to the need of efficiently exploring a vast conformational space in order to seek for the actual transition path for systems whose complexity is already high in the stable states. Here we report a strategy that simplifies this task attacking the complexity on several sides. We first apply a minimalist coarse-grained model to Calmodulin, based on an empirical force field with a partial structural bias, to explore the transition paths between the apo-closed state and the Ca-bound open state of the protein. We then select representative structures along the trajectory based on a structural clustering algorithm and build a cleaned-up trajectory with them. We finally compare this trajectory with that produced by the online tool MinActionPath, by minimizing the action integral using a harmonic network model, and with that obtained by the PROMPT morphing method, based on an optimal mass transportation-type approach including physical constraints. The comparison is performed both on the structural and energetic level, using the coarse-grained and the atomistic force fields upon reconstruction. Our analysis indicates that this method returns trajectories capable of exploring intermediate states with physical meaning, retaining a very low computational cost, which can allow systematic and extensive exploration of the multi-stable proteins transition pathways.

The swimming behaviour of Halobactetium halobium is influenced by variations in light intensity and wavelength. These stimuli are sensed through two retinal-containing proteins: sensory rhodopsin I and II (SR-I and SR-II), each displaying distinct spectra and photocycles. Light as an energy source is used by H. halobium through two other retinal-containing membrane proteins, acting as photon-driven ionic pumps, bacteriorhodopsin and halorhodopsin. In their ground states, SR-I and SR-II have an absorption maximum at 587 and 487 nm respectively. They have been shown to function as photoreceptors that signal, by a chemical transduction chain, to the motor switch controlling the direction of rotation of a flagellar bundle at the cell pole. Changes from clockwise to counterclockwise rotation or vice versa produce reversals in the swimming direction of H. halobium. Phototaxis occurs because the time interval between reversals is altered by response to light stimuli. SR-I acts as a receptor for attractant (red-orange) light stimuli and displays photochromic properties: when absorption in this wavelength range increases, reversals become less frequent and the pigment is transformed into a blue-absorbing species (SR,&. This intermediate of the SR-I photocycle appears to act both as a signalling state for the reversalsuppressing transduction chain and as a receptor for repellent (blue) light stimuli; these are effective in eliciting more frequent reversals when delivered against a red-orange background. SR-II, also called phoborhodopsin, only senses blue-green light (repellent) stimuli, producing an increased frequency of reversals. This sensor-effecter system has been subjected to various types of analysis and constitutes the second best known example of this kind in prokaryotes, the first being chemotaxis in enterobacteria. Differences and similarities between these two natural models are described. Several mathematical models have been proposed to interpret the mechanisms underlying taxis phenomena. These are typically kinetic schemes for the transduction chain, which are partly hypothetical in H. halobium, due to the limited amount of convincing data on the nature of the (chemical) signal reaching the motor switch. Basically, two types of model have been suggested to account for the time interval distribution of reversals in the swimming behaviour of H. halobium: in the first (of which two variants have been proposed), the stochastic properties of the putative transduction chain signalling to the motor switch are involved; in the second, the hypothesis that an endogeneous oscillator is responsible for the motor switch control is maintained. Structural and biochemical data relevant for a model treatment of the processes occurring at the transduction chain are reported. The available behavioural data in H. halobium are also reported and discussed in detail, comparing the experimental approaches and the different interpretations. The implications of the proposed kinetic models are discussed.

Microbial rhodopsins have recently been discovered in pathogenic fungi and have been postulated to be involved in signaling during the course of an infection. Here, we report on the spectroscopic characterization of a light-driven proton pump rhodopsin (UmRh1) from the smut pathogen Ustilago maydis, the causative agent of tumors in maize plants. Electrophysiology, time-resolved UV/Vis and vibrational spectroscopy indicate a pH-dependent photocycle. We also characterized the impact of the auxin hormone indole-3-acetic acid that was shown to influence the pump activity of UmRh1 on individual photocycle intermediates. A facile pumping activity test was established of UmRh1 expressed in Pichia pastoris cells, for probing proton pumping out of the living yeast cells during illumination. We show similarities and distinct differences to the well-known bacteriorhodopsin from archaea and discuss the putative role of UmRh1 in pathogenesis.

Freshwater ecosystems are critical but fragile environments directly affecting society and its welfare. However, our understanding of genuinely freshwater microbial communities, constrained by our capacity to manipulate its prokaryotic participants in axenic cultures, remains very rudimentary. Even the most abundant components, freshwater Actinobacteria, remain largely unknown. Here, applying deep metagenomic sequencing to the microbial community of a freshwater reservoir, we were able to circumvent this traditional bottleneck and reconstruct de novo seven distinct streamlined actinobacterial genomes. These genomes represent three new groups of photoheterotrophic, planktonic Actinobacteria. We describe for the first time genomes of two novel clades, acMicro (Micrococcineae, related to Luna2,) and acAMD (Actinomycetales, related to acTH1). Besides, an aggregate of contigs belonged to a new branch of the Acidimicrobiales. All are estimated to have small genomes (approximately 1.2 Mb), and their GC content varied from 40 to 61%. One of the Micrococcineae genomes encodes a proteorhodopsin, a rhodopsin type reported for the first time in Actinobacteria. The remarkable potential capacity of some of these genomes to transform recalcitrant plant detrital material, particularly lignin-derived compounds, suggests close linkages between the terrestrial and aquatic realms. Moreover, abundances of Actinobacteria correlate inversely to those of Cyanobacteria that are responsible for prolonged and frequently irretrievable damage to freshwater ecosystems. This suggests that they might serve as sentinels of impending ecological catastrophes.

Freshwater ecosystems are critical but fragile environments directly affecting society and its welfare. However, our understanding of genuinely freshwater microbial communities, constrained by our capacity to manipulate its prokaryotic participants in axenic cultures, remains very rudimentary. Even the most abundant components, freshwater Actinobacteria, remain largely unknown. Here, applying deep metagenomic sequencing to the microbial community of a freshwater reservoir, we were able to circumvent this traditional bottleneck and reconstruct de novo seven distinct streamlined actinobacterial genomes. These genomes represent three new groups of photoheterotrophic, planktonic Actinobacteria. We describe for the first time genomes of two novel clades, acMicro (Micrococcineae, related to Luna2,) and acAMD (Actinomycetales, related to acTH1). Besides, an aggregate of contigs belonged to a new branch of the Acidimicrobiales. All are estimated to have small genomes (approximately 1.2 Mb), and their GC content varied from 40 to 61%. One of the Micrococcineae genomes encodes a proteorhodopsin, a rhodopsin type reported for the first time in Actinobacteria. The remarkable potential capacity of some of these genomes to transform recalcitrant plant detrital material, particularly lignin-derived compounds, suggests close linkages between the terrestrial and aquatic realms. Moreover, abundances of Actinobacteria correlate inversely to those of Cyanobacteria that are responsible for prolonged and frequently irretrievable damage to freshwater ecosystems. This suggests that they might serve as sentinels of impending ecological catastrophes.

The light-dependent ion-transport function of microbial rhodopsin has been widely used in optogenetics for optical control of neural activity. In order to increase the variety of rhodopsin proteins having a wide range of absorption wavelengths, the light absorption properties of various wild-type rhodopsins and their artificially mutated variants were investigated in the literature. Here, we demonstrate that a machine-learning-based (ML-based) data-driven approach is useful for understanding and predicting the light-absorption properties of microbial rhodopsin proteins. We constructed a database of 796 proteins consisting of microbial rhodopsin wildtypes and their variants. We then proposed an ML method that produces a statistical model describing the relationship between amino-acid sequences and absorption wavelengths and demonstrated that the fitted statistical model is useful for understanding colour tuning rules and predicting absorption wavelengths. By applying the ML method to...

Bacteriorhodopsin (bR) is a membrane protein from Halobacterium salinarum that carries out a photoenergetic function-directed proton transport through the cellular membrane. It is a unique object for studying molecular mechanisms of light energy conversion in biological systems [1]. In membranes of halobacteria, bR forms patches called purple membranes in which bR trimers are arranged into a regular hexagonal structure [2]. The stability of this two-dimensional crystal is supported by lipids tightly bound to the protein, such as phosphatidyl glycerophosphate (50%) and sulfo-glycolipids (30%). Acting as a light-dependent proton pump, bR generates electrochemical transmembrane potential used for ATP synthesis [2-5]. bR is a chromoprotein with molecular weight of 26.5 kDa; it consists of apoprotein bacterioopsin and a chromophore group, all-trans retinal, bound covalently to K-216 through a protonated Schiff base (SB). bR consists of 248 amino acids tightly arranged in seven transmembrane α helices, some amino acids of which form a chromophore center for retinal insertion. Absorption of a light quantum causes chromophore isomerization into the 13-cis form and is accompanied by several reversible transformations of the bR molecule

Photoinduced proton pumping in bateriorhodopsin (bR) was extensively studied using multiple experimental methods and utilizing various theoretical modeling approaches. These studies usually refer to well-resolved structural data of bacteriorhodopsin. However, despite obtained results, the origin of the proton pumping force initiated by the electronic excitation of retinal remains questionable. Here we discuss the challenges and approaches to modeling the proton transfer in bR and demonstrate that the process, which starts from the electronic excitation of the retinal molecule, is mainly due to the detailed arrangement of the protein environment. Using quantum chemical calculations, we have revealed that the retinal molecule after its excitation is fixed in the ground state of 13-cis,15-syn configuration, as a result of interaction with specific protein residuals. Thus, reaching this fixed configuration, the proton is first transferred to the Asp-85 residue from the water molecule, w...

We have assigned the mouse rhodopsin gene, Rho, to chromosome 6 using DNA from a set of mousehamster somatic hybrid cell lines and a partial cDNA clone for mouse opsin. This assignment rules out the direct involvement of the rhodopsin gene in the known mouse mutations that produce retinal degeneration, including retinal degeneration slow (rds, chromosome 17), retinal degeneration (rd, chromosome 5), Purkinje cell degeneration (pcd, chromosome 13), and nervous (nr, chromosome 8). Segregation of Rho-specific DNA fragment differences among 50 animals from an interspecific backcross (C57BL/6J X Mus spretus) X C57BL/6J indicates th&$he Rho locus is 4.0 + 2.8 map units distal to the l&s for the protooncogene Raf-2 and 18.0 + 5.4 mapunits proximal to the locus for the proto-oncogene Kr'as-2. Linkage to Raf-1 was confirmed using four sets of recombinant inbred strains. The two loci RAF1 and RHO are also syntenic on human chromosome 3, but on opposite arms.

The sequenced genome of the poly-extremophile Exiguobacterium sp. S17, isolated from modern stromatolites at Laguna Socompa (3,570 m), a High-Altitude Andean Lake (HAAL) in Argentinean Puna revealed a putative proteorhodopsin-encoding gene. The HAAL area is exposed to the highest UV irradiation on Earth, making the microbial community living in the stromatolites test cases for survival strategies under extreme conditions. The heterologous expressed protein E17R from Exiguobacterium (248 amino acids, 85% sequence identity to its ortholog ESR from E. sibiricum) was assembled with retinal displaying an absorbance maximum at 524 nm, which makes it a member of the green-absorbing PR-subfamily. Titration down to low pH values (eventually causing partial protein denaturation) indicated a pK value between two and three. Global fitting of data from laser flash-induced absorption changes gave evidence for an early red-shifted intermediate (its formation being below the experimental resolution...

Rhodopsin is the visual pigment responsible for initiating the phototransduction cascade in vertebrate rod photoreceptors. Although well-characterized in a few model systems, comparative studies of rhodopsin function, particularly for nonmammalian vertebrates are comparatively lacking. Bowerbirds are rare among passerines in possessing a key substitution, D83N, at a site that is otherwise highly conserved among G protein-coupled receptors. While this substitution is present in some dim-light adapted vertebrates, often accompanying another unusual substitution, A292S, its functional relevance in birds is uncertain. To investigate functional effects associated with these two substitutions, we use the rhodopsin gene from the great bowerbird (Ptilonorhynchus nuchalis) as a background for site-directed mutagenesis, in vitro expression and functional characterization. We also mutated these sites in two additional rhodopsins that do not naturally possess N83, chicken and bovine, for comparison. Both sites were found to contribute to spectral blue-shifts, but had opposing effects on kinetic rates. Substitutions at site 83 were found to primarily affect the kinetics of light-activated rhodopsin, while substitutions at site 292 had a larger impact on spectral tuning. The contribution of substitutions at site 83 to spectral tuning in particular depended on genetic background, but overall, the effects of substitutions were otherwise surprisingly additive, and the magnitudes of functional shifts were roughly similar across all three genetic backgrounds. By employing a comparative approach with multiple species, our study provides new insight into the joint impact of sites 83 and 292 on rhodopsin structure-function as well as their evolutionary significance for dim-light vision across vertebrates.

Transducin is the retinal rod outer SCllment (ROS).spcci{'~: G protein couplinlt the photo=xcitcd rhodol~in to cydic GMP-pho~photlicstcra=. Tbc = ,tubunit or transdu¢in is known to ~ ADP.ribosylatod by b,actcrial toxins, Wc inv=stillatcd the possibility that transducin is modified in vitro by an cndOllenous ADP.ribosyltranfferas¢ activity. By usin~ either Reg, cytosoli¢ ¢=tract of Re= or purified transdu¢in in the prera=nca of [=.~P]nicotinamide adenine dinucleotid¢ (NAD'). the ¢t and .~ =ubanitx of tran=ducin wcr¢ found to be rndiolabcled, The labclinil wa= dccrcast.d by snake venom phosphodi=stentsc [ (PDE l). The modification wa,= shown to b¢ mono ADP.ribo=ylahon by analysax ~n thin layer chromatollraphy of the PDE I-hydrolyzed products which rev~led only S'AMP rctidu~. [n addition we robert that =odium nitroprusiide =ctivat¢~ the ADP-ribo=ylation of =ran=due=n,

Transitions between different conformational states are ubiquitous in proteins. A vast class of conformation-changing proteins includes evolutionary switches, which vary their conformation as an effect of few mutations or weak environmental variations. However, modeling those processes is extremely difficult due to the need of efficiently exploring a vast conformational space to look for the actual transition path. In this study, we report a strategy that simplifies this task attacking the complexity on several sides. We first apply a minimalist coarse-grained model to the protein, based on an empirical force field with a partial structural bias toward one or both the reference structures. We then explore the transition paths by means of stochastic molecular dynamics and select representative structures by means of a principal path-based clustering algorithm. We finally compare this trajectory with that produced by independent methods adopting a morphing-oriented approach. Our analysis indicates that the minimalist model returns trajectories capable of exploring intermediate states with physical meaning, retaining a very low computational cost, which can allow systematic and extensive exploration of the multistable proteins transition pathways.

Motivation: Atomistic or coarse grained (CG) potentials derived from statistical distributions of internal variables have recently become popular due to the need of simplified interactions for reaching larger scales in simulations or more efficient conformational space sampling. However, the process of parameterization of accurate and predictive statistics-based force fields requires a huge amount of work and is prone to the introduction of bias and errors. Results: This article introduces SecStAnT, a software for the creation and analysis of protein structural datasets with user-defined primary/secondary structure composition, with a particular focus on the CG representation. In addition, the possibility of managing different resolutions and the primary/secondary structure selectivity allow addressing the mapping-backmapping of atomistic to CG representation and study the secondary to primary structure relations. Sample datasets and distributions are reported, including interpretation of structural features. Availability and implementation: SecStAnT is available free of charge at secstant.sourceforge.net/. Source code is freely available on request, implemented in Java and supported on Linux, MS Windows and OSX.

The last decade has witnessed a renewed interest in the coarse-grained (CG) models for biopolymers, also stimulated by the needs of modern molecular biology, dealing with nano- to micro-sized bio-molecular systems and larger than microsecond timescale. This combination of size and timescale is, in fact, hard to access by atomic-based simulations. Coarse graining the system is a route to be followed to overcome these limits, but the ways of practically implementing it are many and different, making the landscape of CG models very vast and complex.In this paper, the CG models are reviewed and their features, applications and performances compared. This analysis, restricted to proteins, focuses on the minimalist models, namely those reducing at minimum the number of degrees of freedom without losing the possibility of explicitly describing the secondary structures. This class includes models using a single or a few interacting centers (beads) for each amino acid.From this analysis seve...

We have assigned the mouse rhodopsin gene, Rho, to chromosome 6 using DNA from a set of mousehamster somatic hybrid cell lines and a partial cDNA clone for mouse opsin. This assignment rules out the direct involvement of the rhodopsin gene in the known mouse mutations that produce retinal degeneration, including retinal degeneration slow (rds, chromosome 17), retinal degeneration (rd, chromosome 5), Purkinje cell degeneration (pcd, chromosome 13), and nervous (nr, chromosome 8). Segregation of Rho-specific DNA fragment differences among 50 animals from an interspecific backcross (C57BL/6J X Mus spretus) X C57BL/6J indicates th&$he Rho locus is 4.0 + 2.8 map units distal to the l&s for the protooncogene Raf-2 and 18.0 + 5.4 mapunits proximal to the locus for the proto-oncogene Kr'as-2. Linkage to Raf-1 was confirmed using four sets of recombinant inbred strains. The two loci RAF1 and RHO are also syntenic on human chromosome 3, but on opposite arms.

We have assigned the mouse rhodopsin gene, Rho, to chromosome 6 using DNA from a set of mouse-hamster somatic hybrid cell lines and a partial cDNA clone for mouse opsin. This assignment rules out the direct involvement of the rhodopsin gene in the known mouse mutations that produce retinal degeneration, including retinal degeneration slow (rds, chromosome 17), retinal degeneration (rd, chromosome 5), Purkinje cell degeneration (pcd, chromosome 13), and nervous (nr, chromosome 8). Segregation of Rho-specific DNA fragment differences among 50 animals from an interspecific backcross (C57BL/6J X Mus spretus) X C57BL/6J indicates that the Rho locus is 4.0 +/- 2.8 map units distal to the locus for the proto-oncogene Raf-1 and 18.0 +/- 5.4 map units proximal to the locus for the proto-oncogene Kras-2. Linkage to Raf-1 was confirmed using four sets of recombinant inbred strains. The two loci RAF1 and RHO are also syntenic on human chromosome 3, but on opposite arms.

The sequenced genome of the poly-extremophile Exiguobacterium sp. S17, isolated from modern stromatolites at Laguna Socompa (3,570 m), a High-Altitude Andean Lake (HAAL) in Argentinean Puna revealed a putative proteorhodopsin-encoding gene. The HAAL area is exposed to the highest UV irradiation on Earth, making the microbial community living in the stromatolites test cases for survival strategies under extreme conditions. The heterologous expressed protein E17R from Exiguobacterium (248 amino acids, 85% sequence identity to its ortholog ESR from E. sibiricum) was assembled with retinal displaying an absorbance maximum at 524 nm, which makes it a member of the green-absorbing PR-subfamily. Titration down to low pH values (eventually causing partial protein denaturation) indicated a pK value between two and three. Global fitting of data from laser flash-induced absorption changes gave evidence for an early red-shifted intermediate (its formation being below the experimental resolution...

Photochemical reaction dynamics of the primary events in recombinant bacteriorhodopsin (bRrec) was studied by
femtosecond laser absorption spectroscopy with 25-fs time resolution. bRrec was produced in an Escherichia coli expression system. Since bRrec was prepared in a DMPC–CHAPS micelle system in the monomeric form, its comparison with trimeric and monomeric forms of the native bacteriorhodopsin (bRtrim and bRmon, respectively) was carried out. We found that bRrec intermediate I (excited state of bR) was formed in the range of 100 fs, as in the case of bRtrim and bRmon. Further processes, namely the decay of the excited state I and the formation of intermediates J and K of bRrec, occurred more slowly compared to bRtrim, but similarly to bRmon. The lifetime of intermediate I, judging from the signal of ΔAESA(470480 nm), was 0.68 ps (78%) and 4.4 ps (22%) for bRrec, 0.52 ps (73%) and 1.7 ps (27%) for bRmon, and 0.45 ps (90%) and 1.75 ps (10%) for bRtrim. The formation time of intermediate K, judging from the signal of ΔAGSA(625635 nm), was 13.5 ps for bRrec, 9.8 ps for bRmon, and 4.3 ps for bRtrim. In addition, there was a decrease in the photoreaction efficiency of bRrec and bRmon as seen by a decrease in absorbance in the differential spectrum of the intermediate K by ∼14%. Since photochemical properties of bRrec are similar to those of the monomeric form of the native protein, bRrec and its mutants can be considered as a basis for further studies of the mechanism of bacteriorhodopsin functioning.

The sequenced genome of the poly-extremophile Exiguobacterium sp. S17, isolated from modern stromatolites at Laguna Socompa (3,570 m), a High-Altitude Andean Lake (HAAL) in Argentinean Puna revealed a putative proteorhodopsin-encoding gene. The HAAL area is exposed to the highest UV irradiation on Earth, making the microbial community living in the stromatolites test cases for survival strategies under extreme conditions. The heterolo-gous expressed protein E17R from Exiguobacterium (248 amino acids, 85% sequence identity to its ortholog ESR from E. sibiricum) was assembled with retinal displaying an absorbance maximum at 524 nm, which makes it a member of the green-absorbing PR-sub-family. Titration down to low pH values (eventually causing partial protein denaturation) indicated a pK value between two and three. Global fitting of data from laser flash-induced absorption changes gave evidence for an early red-shifted intermediate (its formation being below the experimental resolution) that decayed (τ 1 = 3.5 μs) into another red-shifted intermediate. This species decayed in a two-step process (τ 2 = 84 μs, τ 3 = 11 ms), to which the initial state of E17-PR was reformed with a kinetics of 2 ms. Proton transport capability of the HAAL protein was determined by BLM measurements. Additional blue light irradiation reduced the proton current, clearly identifying a blue light absorbing, M-like intermediate. The apparent absence of this intermediate is explained by closely matching formation and decay kinetics.

Atomistic or coarse grained (CG) potentials derived from statistical distributions of internal variables have recently become popular due to the need of simplified interactions for reaching larger scales in simulations or more efficient conformational space sampling. However, the process of parameterization of accurate and predictive statistics-based force fields requires a huge amount of work and is prone to the introduction of bias and errors. This article introduces SecStAnT, a software for the creation and analysis of protein structural datasets with user-defined primary/secondary structure composition, with a particular focus on the CG representation. In addition, the possibility of managing different resolutions and the primary/secondary structure selectivity allow addressing the mapping-backmapping of atomistic to CG representation and study the secondary to primary structure relations. Sample datasets and distributions are reported, including interpretation of structural features. SecStAnT is available free of charge at secstant.sourceforge.net/. Source code is freely available on request, implemented in Java and supported on Linux, MS Windows and OSX.

The activity within a living cell is based on a complex network of interactions among biomolecules, exchanging information and energy through biochemical processes. These events occur on different scales, from the nano-to the macroscale, spanning about 10 orders of magnitude in the space domain and 15 orders of magnitude in the time domain. Consequently, many different modeling techniquesÑeach proper for a particular time or space scaleÑare commonly used. In addition, a single process often spans more than a single time or space scale. Thus, the necessity arises for combining the modeling techniques in multiscale approaches.

Eukaryotic microbial rhodopsins are widespread bacteriorhodopsin-like proteins found in many lower eukaryotic groups including fungi. Many fungi contain multiple rhodopsins, some significantly diverged from the original bacteriorhodopsin template. Although few fungal rhodopsins have been studied biophysically, both fast-cycling light-driven proton pumps and slow-cycling photosensors have been found. The purpose of this study was to characterize photochemically a new subgroup of fungal rhodopsins, the so-called auxiliary group. The study used the two known rhodopsin genes from the fungal wheat pathogen, Phaeosphaeria nodorum. One of the genes is a member of the auxiliary group while the other is highly similar to previously characterized proton-pumping Leptosphaeria rhodopsin. Auxiliary rhodopsin genes from a range of species form a distinct group with a unique primary structure and are located in carotenoid biosynthesis gene cluster. Amino acid conservation pattern suggests that auxiliary rhodopsins retain the transmembrane core of bacteriorhodopsins, including all residues important for proton transport, but have unique polar intramembrane residues. Spectroscopic characterization of the two yeast-expressed Phaeosphaeria rhodopsins showed many similarities: absorption spectra, conformation of the retinal chromophore, fast photocycling, and carboxylic acid protonation changes. It is likely that both Phaeosphaeria rhodopsins are proton-pumping, at least in vitro. We suggest that auxiliary rhodopsins have separated from their ancestors fairly recently and have acquired the ability to interact with as yet unidentified transducers, performing a photosensory function without changing their spectral properties and basic photochemistry.

Atomistic or coarse grained (CG) potentials derived from statistical distributions of internal variables have recently become popular due to the need of simplified interactions for reaching larger scales in simulations or more efficient conformational space sampling. However, the process of parameterization of accurate and predictive statistics-based force fields requires a huge amount of work and is prone to the introduction of bias and errors. This article introduces SecStAnT, a software for the creation and analysis of protein structural datasets with user-defined primary/secondary structure composition, with a particular focus on the CG representation. In addition, the possibility of managing different resolutions and the primary/secondary structure selectivity allow addressing the mapping-backmapping of atomistic to CG representation and study the secondary to primary structure relations. Sample datasets and distributions are reported, including interpretation of structural features. SecStAnT is available free of charge at secstant.sourceforge.net/. Source code is freely available on request, implemented in Java and supported on Linux, MS Windows and OSX.

Since algal rhodopsins, the eukaryotic seven-transmembrane proteins, are generally difficult to express in Escherichia coli, eukaryotic cells have been used for heterologous expression. Mistic, a membrane-associated protein that was originally discovered in Bacillus subtilis, has been shown to improve the expression levels of many foreign integral membrane proteins in E. coli when used as a fusion partner linked to the N-terminus of cargo proteins.
Here, we expressed two algal rhodopsins with N- and C-terminal Mistic domains in E. coli - Acetabularia rhodopsin I (ARI) and Chlamydomonas sensory rhodopsin B (CSRB, channel rhodopsin 2). UV/VIS spectroscopy, pH titration of proton acceptor residue, laser-induced photolysis and electrophysiological measurement were used for investigating important residues in proton transport and spectroscopic characters of the proteins.
Protein yield of two algal rhodopsins was enhanced, obtaining 0.12 mg of Mistic-ARI and 0.04 mg of Mistic-CSRB per liter of culture. Spheroplasts expression Mistic-ARI had outward proton-pumping activity, indicating protein functionality. Asp89 of ARI changed its protonation state by light absorption, and Asp100 was important for O(600) formation. Electrophysiology revealed that both residues took part in proton transport. The spectroscopic analyses of Mistic-CSRB revealed its characteristics.
Fusion to the membrane-integrating protein Mistic can enhance overexpression of eukaryotic type I rhodopsins in E. coli.
These findings indicate that Mistic fusion and E. coli expression method could be an effective, low cost technique for studying eukaryotic membrane proteins. This may have useful implications, for example, in studying structural characteristics and optogenetics for rhodopsins.

We cloned an intronless, nuclear-encoded opsin gene from an EST library of Acetabularia acetabulum. Acetabularia rhodopsin I (ARI) encodes a protein of 246 amino acids with molecular weight of 27 kDa. ARI was reconstituted in the Xenopus oocyte expression system to characterize its electrophysiological properties utilizing the two-electrode voltage-clamping technique. Oocytes where ARI cRNA was injected displayed outward directed currents in response to light. The maximum action spectrum of ARI was detected at 520 nm green light. Light-stimulated ARI current amplitude was altered by the protons, but not by the other ions in recording solutions, suggesting that the algal rhodopsin is a light-stimulated proton pump. Typical proton-mediated outward current elicited by 520 nm light was characterized with two phases of non-inactivating outward current following initial transient current. Taken together, we here reported cloning of a novel Acetabularia opsin gene which was characterized to be a proton-pump stimulated by light.

Detergent-resistant membrane microdomains in the plasma membrane, known as lipid rafts, have been implicated in various cellular processes. We report here that a low-density Triton X-100-insoluble membrane (detergent-resistant membrane; DRM) fraction is present in bovine rod photoreceptor outer segments (ROS). In dark-adapted ROS, transducin and most of cGMPphosphodiesterase (PDE) were detergent-soluble. When ROS membranes were exposed to light, however, a large portion of transducin localized in the DRM fraction. Furthermore, on addition of guanosine 5-3-O-(thio)triphosphate (GTP␥S) to light-bleached ROS, transducin became detergent-soluble again. PDE was not recruited to the DRM fraction after light stimulus alone, but simultaneous stimulation by light and GTP␥S induced a massive translocation of all PDE subunits to the DRM. A cholesterol-removing reagent, methyl-␤-cyclodextrin, selectively but partially solubilized PDE from the DRM, suggesting that cholesterol contributes, at least in part, to the association of PDE with the DRM. By contrast, transducin was not extracted by the depletion of cholesterol. These data suggest that transducin and PDE are likely to perform their functions in phototransduction by changing their localization between two distinct lipid phases, rafts and surrounding fluid membrane, on disc membranes in an activation-dependent manner.