Molecular mechanism of spectral tuning in sensory rhodopsin II

Structural and spectroscopic characterization of the peridinin-chlorophyll a-protein (PCP) complex from Heterocapsa pygmaea (HPPCP)

Light harvesting proteins are optimized to efficiently collect and transfer light energy for photosynthesis. In eukaryotic dinoflagellates these complexes utilize chlorophylls and a special carotenoid, peridinin, and arrange them for efficient excitation energy transfer. At the same time, the carotenoids protect the system by quenching harmful chlorophyll triplet states. Here we use advanced spectroscopic techniques and X-ray structure analysis to investigate excitation energy transfer processes in the major soluble antenna, the peridinin chlorophyll a protein (PCP) from the free living dinoflagellate Heterocapsa pygmaea. We determined the 3D-structure of this complex at high resolution (1.2 Å). For better comparison, we improved the reference structure of this protein from Amphidinium carterae to a resolution of 1.15 Å. We then used fs and ns time-resolved absorption spectroscopy to study the mechanisms of light harvesting, but also of the photoprotective quenching of the chlorophyll triplet state. The photoprotection site was further characterized by Electron Spin Echo Envelope Modulation (ESEEM) spectroscopy to yield information on water molecules involved in triplet-triplet energy transfer. Similar to other PCP complexes, excitation energy transfer from peridinin to chlorophyll is found to be very efficient, with transfer times in the range of 1.6-2.1 ps. One of the four carotenoids, the peridinin 614, is well positioned to quench the chlorophyll triplet state with high efficiency and transfer times in the range of tens of picoseconds. Our structural and dynamic data further support, that the intrinsic water molecule coordinating the chlorophyll Mg ion plays an essential role in photoprotection.

A conserved Trp residue in HwBR contributes to its unique tolerance toward acidic environments

Bacteriorhodopsin (BR) is a light-driven outward proton pump found mainly in halophilic archaea. A BR from an archaeon Haloquadratum walsbyi (HwBR) was found to pump protons under more acidic conditions compared with most known BR proteins. The atomic structural study on HwBR unveiled that a pair of hydrogen bonds between the BC and FG loop in its periplasmic region may be a factor in such improved pumping capability. Here, we further investigated the retinal-binding pocket of HwBR and found that Trp94 contributes to the higher acid tolerance. Through single mutations in a BR from Halobacterium salinarum and HwBR, we examined the conserved tryptophan residues in the retinal-binding pocket. Among these residues of HwBR, mutagenesis at Trp94 facing the periplasmic region caused the most significant disruption to optical stability and proton-pumping capability under acidic conditions. The other tryptophan residues of HwBR exerted little impact on both maximum absorption wavelength and pH-dependent proton pumping. Our findings suggest that the residues from Trp94 to the hydrogen bonds at the BC loop confer both optical stability and functionality on the overall protein in low-pH environments.

Earliest photic zone niches probed by ancestral microbial rhodopsins

For billions of years, life has continuously adapted to dynamic physical conditions near the Earth’s surface. Fossils and other preserved biosignatures in the paleontological record are the most direct evidence for reconstructing the broad historical contours of this adaptive interplay. However, biosignatures dating to Earth’s earliest history are exceedingly rare. Here, we combine phylogenetic inference of primordial rhodopsin proteins with modeled spectral features of the Precambrian Earth environment to reconstruct the paleobiological history of this essential family of photoactive transmembrane proteins. Our results suggest that ancestral microbial rhodopsins likely acted as light-driven proton pumps and were spectrally tuned toward the absorption of green light, which would have enabled their hosts to occupy depths in a water column or biofilm where UV wavelengths were attenuated. Subsequent diversification of rhodopsin functions and peak absorption frequencies was enabled by the expansion of surface ecological niches induced by the accumulation of atmospheric oxygen. Inferred ancestors retain distinct associations between extant functions and peak absorption frequencies. Our findings suggest that novel information encoded by biomolecules can be used as “paleosensors” for conditions of ancient, inhabited niches of host organisms not represented elsewhere in the paleontological record. The coupling of functional diversification and spectral tuning of this taxonomically diverse protein family underscores the utility of rhodopsins as universal testbeds for inferring remotely detectable biosignatures on inhabited planetary bodies.

Electronic Couplings and Electrostatic Interactions Behind the Light Absorption of Retinal Proteins

The photo-functional chromophore retinal exhibits a wide variety of optical absorption properties depending on its intermolecular interactions with surrounding proteins and other chromophores. By utilizing these properties, microbial and animal rhodopsins express biological functions such as ion-transport and signal transduction. In this review, we present the molecular mechanisms underlying light absorption in rhodopsins, as revealed by quantum chemical calculations. Here, symmetry-adapted cluster-configuration interaction (SAC-CI), combined quantum mechanical and molecular mechanical (QM/MM), and transition-density-fragment interaction (TDFI) methods are used to describe the electronic structure of the retinal, the surrounding protein environment, and the electronic coupling between chromophores, respectively. These computational approaches provide successful reproductions of experimentally observed absorption and circular dichroism (CD) spectra, as well as insights into the mechanisms of unique optical properties in terms of chromophore-protein electrostatic interactions and chromophore-chromophore electronic couplings. On the basis of the molecular mechanisms revealed in these studies, we also discuss strategies for artificial design of the optical absorption properties of rhodopsins.

Insights into the Protein Functions and Absorption Wavelengths of Microbial Rhodopsins

Using a quantum mechanical/molecular mechanical approach, the absorption wavelength of the retinal Schiff base was calculated based on 13 microbial rhodopsin crystal structures. The results showed that the protein electrostatic environment decreases the absorption wavelength significantly in the cation-conducting rhodopsin but only slightly in the sensory rhodopsin. Among the microbial rhodopsins with different functions, the differences in the absorption wavelengths are caused by differences in the arrangement of the charged residues at the retinal Schiff base binding moiety, namely, one or two counterions at the three common positions. Among the microbial rhodopsins with similar functions, the differences in the polar residues at the retinal Schiff base binding site are responsible for the differences in the absorption wavelengths. Counterions contribute to an absorption wavelength shift of 50-120 nm, whereas polar groups contribute to a shift of up to ∼10 nm. It seems likely that protein function is directly associated with the absorption wavelength in microbial rhodopsins.

Raman spectroscopy of a near infrared absorbing proteorhodopsin: Similarities to the bacteriorhodopsin O photointermediate

Microbial rhodopsins have become an important tool in the field of optogenetics. However, effective in vivo optogenetics is in many cases severely limited due to the strong absorption and scattering of visible light by biological tissues. Recently, a combination of opsin site-directed mutagenesis and analog retinal substitution has produced variants of proteorhodopsin which absorb maximally in the near-infrared (NIR). In this study, UV-Visible-NIR absorption and resonance Raman spectroscopy were used to study the double mutant, D212N/F234S, of green absorbing proteorhodopsin (GPR) regenerated with MMAR, a retinal analog containing a methylamino modified β-ionone ring. Four distinct subcomponent absorption bands with peak maxima near 560, 620, 710 and 780 nm are detected with the NIR bands dominant at pH <7.3, and the visible bands dominant at pH 9.5. FT-Raman using 1064-nm excitation reveal two strong ethylenic bands at 1482 and 1498 cm-1 corresponding to the NIR subcomponent absorption bands based on an extended linear correlation between λmax and γC = C. This spectrum exhibits two intense bands in the fingerprint and HOOP mode regions that are highly characteristic of the O640 photointermediate from the light-adapted bacteriorhodopsin photocycle. In contrast, 532-nm excitation enhances the 560-nm component, which exhibits bands very similar to light-adapted bacteriorhodopsin and/or the acid-purple form of bacteriorhodopsin. Native GPR and its mutant D97N when regenerated with MMAR also exhibit similar absorption and Raman bands but with weaker contributions from the NIR absorbing components. Based on these results it is proposed that the NIR absorption in GPR-D212N/F234S with MMAR arises from an O-like chromophore, where the Schiff base counterion D97 is protonated and the MMAR adopts an all-trans configuration with a non-planar geometry due to twists in the conjugated polyene segment. This configuration is characterized by extensive charge delocalization, most likely involving nitrogens atoms in the MMAR chromophore.

Purification and Characterization of RhoPDE, a Retinylidene/Phosphodiesterase Fusion Protein and Potential Optogenetic Tool from the Choanoflagellate Salpingoeca rosetta

RhoPDE is a type I rhodopsin/phosphodiesterase gene fusion product from the choanoflagellate Salpingoeca rosetta. The gene was discovered around the time that a similar type I rhodopsin/guanylyl cyclase fusion protein, RhoGC, was shown to control phototaxis of an aquatic fungus through a cGMP signaling pathway. RhoPDE has potential as an optogenetic tool catalyzing the hydrolysis of cyclic nucleotides. Here we provide an expression and purification system for RhoPDE, as well as a crystal structure of the C-terminal phosphodiesterase catalytic domain. We show that RhoPDE contains an even number of transmembrane segments, with N- and C-termini both located on the cytoplasmic surface of the cell membrane. The purified protein exhibits an absorption maximum at 490 nm in the dark state, which shifts to 380 nm upon exposure to light. The protein acts as a cGMP-selective phosphodiesterase. However, the activity does not appear to be modulated by light. The protein is also active with cAMP as a substrate, but with a roughly 5-7-fold lower k_cat. A truncation consisting solely of the phosphodiesterase domain is also active with a k_cat for cGMP roughly 6-9-fold lower than that of the full-length protein. The isolated PDE domain was crystallized, and the X-ray structure showed the protein to be a dimer similar to human PDE9. We anticipate that the purification system introduced here will enable further structural and biochemical experiments to improve our understanding of the function and mechanism of this unique fusion protein.

Theoretical Studies on the Color-Tuning Mechanism in Retinal Proteins

The excited states of the three retinal proteins, bovine rhodopsin (Rh), bacteriorhodopsin (bR), and sensory rhodopsin II (sRII) were studied using the symmetry-adapted cluster-configuration interaction (SAC-CI) and combined quantum mechanical and molecular mechanical (QM/MM) methods. The computed absorption energies are in good agreement with the experimental ones for all three proteins. The spectral tuning mechanism was analyzed in terms of three contributions:  molecular structures of the chromophore in the binding pockets, electrostatic (ES) interaction of the chromophore with the surrounding protein environment, and quantum-mechanical effect between the chromophore and the counterion group. This analysis provided an insight into the mechanism of the large blue-shifts in the absorption peak position of Rh and sRII from that of bR. Protein ES effect is primarily important both in Rh and in sRII, and the structure effect is secondary important in Rh. The quantum-mechanical interaction between the chromophore and the counterion is very important for quantitative reproduction of the excitation energy. These results indicate that the present approach is useful for studying the absorption spectra and the mechanism of the color tuning in the retinal proteins.

Protein Influence on Electronic Spectra Modeled by Multipoles and Polarizabilities

We have developed automatic methods to calculate multipoles and anisotropic polarizabilities for all atoms and bond centers in a protein and to include such a model in the calculation of electronic properties at any level of quantum mechanical theory. This approach is applied for the calculation of the electronic spectra of retinal in rhodopsin at the CASPT2//CASSCF level (second-order multiconfigurational perturbation theory) for the wild-type protein, as well as two mutants and isorhodopsin in QM/MM structures based on two crystal structures. We also perform a detailed investigation of the importance and distance dependence of the multipoles and the polarizabilities for both the absolute and the relative absorption energies. It is shown that the model of the surrounding protein strongly influences the spectrum and that different models give widely different results. For example, the Amber 1994 and 2003 force fields give excitation energies that differ by up to 16 kJ/mol. For accurate excitation energies, multipoles up to quadrupoles and anisotropic polarizabilities are needed. However, interactions with residues more than 10 Å from the chromophore can be treated with a standard polarizable force field without any dipoles or quadrupoles.

Molecular Mechanism of Wide Photoabsorption Spectral Shifts of Color Variants of Human Cellular Retinol Binding Protein II

Color variants of human cellular retinol binding protein II (hCRBPII) created by protein engineering were recently shown to exhibit anomalously wide photoabsorption spectral shifts over ∼200 nm across the visible region. The remarkable phenomenon provides a unique opportunity to gain insight into the molecular basis of the color tuning of retinal binding proteins for understanding of color vision as well as for engineering of novel color variants of retinal binding photoreceptor proteins employed in optogenetics. Here, we report a theoretical investigation of the molecular mechanism underlying the anomalously wide spectral shifts of the color variants of hCRBPII. Computational modeling of the color variants with hybrid molecular simulations of free energy geometry optimization succeeded in reproducing the experimentally observed wide spectral shifts, and revealed that protein flexibility, through which the active site structure of the protein and bound water molecules is altered by remote mutations, plays a significant role in inducing the large spectral shifts.

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.

High efficiency light harvesting by carotenoids in the LH2 complex from photosynthetic bacteria: unique adaptation to growth under low-light conditions

Rhodopin, rhodopinal, and their glucoside derivatives are carotenoids that accumulate in different amounts in the photosynthetic bacterium, Rhodoblastus (Rbl.) acidophilus strain 7050, depending on the intensity of the light under which the organism is grown. The different growth conditions also have a profound effect on the spectra of the bacteriochlorophyll (BChl) pigments that assemble in the major LH2 light-harvesting pigment-protein complex. Under high-light conditions the well-characterized B800-850 LH2 complex is formed and accumulates rhodopin and rhodopin glucoside as the primary carotenoids. Under low-light conditions, a variant LH2, denoted B800-820, is formed, and rhodopinal and rhodopinal glucoside are the most abundant carotenoids. The present investigation compares and contrasts the spectral properties and dynamics of the excited states of rhodopin and rhodopinal in solution. In addition, the systematic differences in pigment composition and structure of the chromophores in the LH2 complexes provide an opportunity to explore the effect of these factors on the rate and efficiency of carotenoid-to-BChl energy transfer. It is found that the enzymatic conversion of rhodopin to rhodopinal by Rbl. acidophilus 7050 grown under low-light conditions results in nearly 100% carotenoid-to-BChl energy transfer efficiency in the LH2 complex. This comparative analysis provides insight into how photosynthetic systems are able to adapt and survive under challenging environmental conditions.

The photochemical determinants of color vision: revealing how opsins tune their chromophore's absorption wavelength

The evolution of a variety of important chromophore-dependent biological processes, including microbial light sensing and mammalian color vision, relies on protein modifications that alter the spectral characteristics of a bound chromophore. Three different color opsins share the same chromophore, but have three distinct absorptions that together cover the entire visible spectrum, giving rise to trichromatic vision. The influence of opsins on the absorbance of the chromophore has been studied through methods such as model compounds, opsin mutagenesis, and computational modeling. The recent development of rhodopsin mimic that uses small soluble proteins to recapitulate the binding and wavelength tuning of the native opsins provides a new platform for studying protein-regulated spectral tuning. The ability to achieve far-red shifted absorption in the rhodopsin mimic system was attributed to a combination of the lack of a counteranion proximal to the iminium, and a uniformly neutral electrostatic environment surrounding the chromophore.

Molecular and evolutionary aspects of microbial sensory rhodopsins

Retinal proteins (~rhodopsins) are photochemically reactive membrane-embedded proteins, with seven transmembrane α-helices which bind the chromophore retinal (vitamin A aldehyde). They are widely distributed through all three biological kingdoms, eukarya, bacteria and archaea, indicating the biological significance of the retinal proteins. Light absorption by the retinal proteins triggers a photoisomerization of the chromophore, leading to the biological function, light-energy conversion or light-signal transduction. This article reviews molecular and evolutionary aspects of the light-signal transduction by microbial sensory receptors and their related proteins. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

The nature of the intramolecular charge transfer state in peridinin

Experimental and theoretical evidence is presented that supports the theory that the intramolecular charge transfer (ICT) state of peridinin is an evolved state formed via excited-state bond-order reversal and solvent reorganization in polar media. The ICT state evolves in <100 fs and is characterized by a large dipole moment (~35 D). The charge transfer character involves a shift of electron density within the polyene chain, and it does not involve participation of molecular orbitals localized in either of the β-rings. Charge is moved from the allenic side of the polyene into the furanic ring region and is accompanied by bond-order reversal in the central portion of the polyene chain. The electronic properties of the ICT state are generated via mixing of the "1(1)Bu(+)" ionic state and the lowest-lying "2(1)Ag(-)" covalent state. The resulting ICT state is primarily (1)Bu(+)-like in character and exhibits not only a large oscillator strength but an unusually large doubly excited character. In most solvents, two populations exist in equilibrium, one with a lowest-lying ICT ionic state and a second with a lowest-lying "2(1)Ag(-)" covalent state. The two populations are separated by a small barrier associated with solvent relaxation and cavity formation.

Full-Quantum chemical calculation of the absorption maximum of bacteriorhodopsin: a comprehensive analysis of the amino acid residues contributing to the opsin shift

Herein, the absorption maximum of bacteriorhodopsin (bR) is calculated using our recently developed method in which the whole protein can be treated quantum mechanically at the level of INDO/S-CIS//ONIOM (B3LYP/6-31G(d,p): AMBER). The full quantum mechanical calculation is shown to reproduce the so-called opsin shift of bR with an error of less than 0.04 eV. We also apply the same calculation for 226 different bR mutants, each of which was constructed by replacing any one of the amino acid residues of the wild-type bR with Gly. This substitution makes it possible to elucidate the extent to which each amino acid contributes to the opsin shift and to estimate the inter-residue synergistic effect. It was found that one of the most important contributions to the opsin shift is the electron transfer from Tyr185 to the chromophore upon excitation. We also indicate that some aromatic (Trp86, Trp182) and polar (Ser141, Thr142) residues, located in the vicinity of the retinal polyene chain and the β-ionone ring, respectively, play an important role in compensating for the large blue-shift induced by both the counterion residues (Asp85, Asp212) and an internal water molecule (W402) located near the Schiff base linkage. In particular, the effect of Trp86 is comparable to that of Tyr185. In addition, Ser141 and Thr142 were found to contribute to an increase in the dipole moment of bR in the excited state. Finally, we provide a complete energy diagram for the opsin shift together with the contribution of the chromophore-protein steric interaction.

meso-arylporpholactones and their reduction products

The rational syntheses of meso-tetraaryl-3-oxo-2-oxaporphyrins 5, known as porpholactones, via MnO(4)(-)-mediated oxidations of the corresponding meso-tetraaryl-2,3-dihydroxychlorins (7) is detailed. Since chlorin 7 is prepared from the parent porphyrin 1, this amounts to a 2-step replacement of a pyrrole moiety in 1 by an oxazolone moiety. The stepwise reduction of the porpholactone 5 results in the formation of chlorin analogues, meso-tetraaryl-3-hydroxy-2-oxachlorin (11) and meso-tetraaryl-2-oxachlorins (12). The reactivity of 11 with respect to nucleophilic substitution by O-, N-, and S-nucleophiles is described. The profound photophysical consequences of the formal replacement of a pyrrole with an oxazolone (porphyrin-like chromophore) or (substituted) oxazole moiety (chlorin-like chromophore with, for the parent oxazolochlorin 12, red-shifted Q(x) band with enhanced oscillator strengths) are detailed and rationalized on the basis of SAC-CI and MNDO-PSDCI molecular orbital theory calculations. The single crystal X-ray structures of the porpholactones point at a minor steric interaction between the carbonyl oxygen and the flanking phenyl group. The essentially planar structures of all chromophores in all oxidation states prove that the observed optical properties originate from the intrinsic electronic properties of the chromophores and are not subject to conformational modulation.

A transporter converted into a sensor, a phototaxis signaling mutant of bacteriorhodopsin at 3.0 Å

Bacteriorhodopsin (BR) and sensory rhodopsin II (SRII), homologous photoactive proteins in haloarchaea, have different molecular functions. BR is a light-driven proton pump, whereas SRII is a phototaxis receptor that transmits a light-induced conformational change to its transducer HtrII. Despite these distinctly different functions, a single residue substitution, Ala215 to Thr215 in the BR retinal-binding pocket, enables its photochemical reactions to transmit signals to HtrII and mediate phototaxis. We pursued a crystal structure of the signaling BR mutant (BR_A215T) to determine the structural changes caused by the A215T mutation and to assess what new photochemistry is likely to be introduced into the BR photoactive site. We crystallized BR_A215T from bicelles and solved its structure to 3.0 Å resolution to enable an atomic-level comparison. The analysis was complemented by molecular dynamics simulation of BR mutated in silico. Three main conclusions regarding the roles of photoactive site residues in signaling emerge from the comparison of BR_A215T, BR, and SRII structures: (i) the Thr215 residue in signaling BR is positioned nearly identically with respect to the retinal chromophore as in SRII, consistent with its role in producing a steric conflict with the retinal C₁₄ group during photoisomerization, proposed earlier to be essential for SRII signaling from vibrational spectroscopy and motility measurements; (ii) Tyr174-Thr204 hydrogen bonding, critical in SRII signaling and mimicked in signaling BR, is likely auxiliary, for example, to maintain Thr204 in the proper position for the steric trigger to occur; and (iii) the primary role of Arg72 in SRII is spectral tuning and not signaling.

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.

The opsin shift and mechanism of spectral tuning in rhodopsin

Molecular dynamics simulations and combined quantum mechanical and molecular mechanical calculations have been performed to investigate the mechanism of the opsin shift and spectral tuning in rhodopsin. A red shift of -980 cm(-1) was estimated in the transfer of the chromophore from methanol solution environment to the protonated Schiff base (PSB)-binding site of the opsin. The conformational change from a 6-s-cis-all-trans configuration in solution to the 6-s-cis-11-cis conformer contributes additional -200 cm(-1), and the remaining effects were attributed to dispersion interactions with the aromatic residues in the binding site. An opsin shift of 2100 cm(-1) was obtained, in reasonable accord with experiment (2730 cm(-1)). Dynamics simulations revealed that the 6-s-cis bond can occupy two main conformations for the β-ionone ring, resulting in a weighted average dihedral angle of about -50°, which may be compared with the experimental estimate of -28° from solid-state NMR and Raman data. We investigated a series of four single mutations, including E113D, A292S, T118A, and A269T, which are located near the PSB, along the polyene chain of retinal and close to the ionone ring. The computational results on absorption energy shift provided insights into the mechanism of spectral tuning, which involves all means of electronic structural effects, including the stabilization or destabilization of either the ground or the electronically excited state of the retinal PSB.

Femtosecond transient absorption spectroscopic study of a carbonyl-containing carotenoid analogue, 2-(all-trans-retinylidene)-indan-1,3-dione

The photophysical properties of a carbonyl-containing carotenoid analogue in an s-cis configuration, relative to the conjugated π system, 2-(all-trans-retinylidene)-indan-1,3-dione (C20Ind), were investigated by femtosecond time-resolved spectroscopy in various solvents. The lifetime of the optically forbidden S(1) state of C20Ind becomes long as solvent polarity increases. This trend is completely opposite to the situation of S(1-ICT) dynamics of carbonyl-containing carotenoids, such as peridinin and fucoxanthin. Excitation energy dependence of the transient absorption measurements shows that the transient absorption spectra in nonpolar solvents were originated from two distinct transient species, while those in polar and protic solvents are due to a single transient species. By referring to the results of MNDO-PSDCI (modified neglect of differential overlap with partial single- and double-configuration interaction) calculations, we conclude: (1) in polar and protic solvents, the S(1) state is generated following excitation up to the S(2) state; (2) in nonpolar solvents, however, both the S(1) and the (1)nπ* states are generated; and (3) C20Ind does not generate the S(1-ICT) state, despite the fact that it has two conjugated carbonyl groups.

Spectral tuning in sensory rhodopsin I from Salinibacter ruber

Organisms utilize light as energy sources and as signals. Rhodopsins, which have seven transmembrane α-helices with retinal covalently linked to a conserved Lys residue, are found in various organisms as distant in evolution as bacteria, archaea, and eukarya. One of the most notable properties of rhodopsin molecules is the large variation in their absorption spectrum. Sensory rhodopsin I (SRI) and sensory rhodopsin II (SRII) function as photosensors and have similar properties (retinal composition, photocycle, structure, and function) except for their λ(max) (SRI, ∼560 nm; SRII, ∼500 nm). An expression system utilizing Escherichia coli and the high protein stability of a newly found SRI-like protein, SrSRI, enables studies of mutant proteins. To determine the residue contributing to the spectral shift from SRI to SRII, we constructed various SRI mutants, in which individual residues were substituted with the corresponding residues of SRII. Three such mutants of SrSRI showed a large spectral blue-shift (>14 nm) without a large alteration of their retinal composition. Two of them, A136Y and A200T, are newly discovered color tuning residues. In the triple mutant, the λ(max) was 525 nm. The inverse mutation of SRII (F134H/Y139A/T204A) generated a spectral-shifted SRII toward longer wavelengths, although the effect is smaller than in the case of SRI, which is probably due to the lack of anion binding in the SRII mutant. Thus, half of the spectral shift from SRI to SRII could be explained by only those three residues taking into account the effect of Cl(-) binding.

The intramolecular charge transfer state in carbonyl-containing polyenes and carotenoids

Numerous femtosecond time-resolved optical spectroscopic experiments have reported that the lifetime of the low-lying S(1) state of carbonyl-containing polyenes and carotenoids decreases with increasing solvent polarity. The effect becomes even more pronounced as the number of double bonds in the conjugated π-electron system decreases. The effect has been attributed to an intramolecular charge transfer (ICT) state coupled to S(1), but it is still not clear what the precise molecular nature of this state is, and how it is able to modulate the spectral and dynamic properties of polyenes and carotenoids. In this work, we examine the nature of the ICT state in three substituted polyenes: crocetindial, which contains two terminal, symmetrically substituted carbonyl groups in conjugation with the π-electron system, 8,8'-diapocarotene-8'-ol-8-al, which has one terminal conjugated carbonyl group and one hydroxyl group, and 8,8'-diapocarotene-8,8'-diol, which has two terminal, symmetrically positioned, hydroxyl groups but no carbonyls. Femtosecond time-resolved optical spectroscopic experiments on these molecules reveal that only the asymmetrically substituted 8,8'-diapocarotene-8'-ol-8-al exhibits any substantial effect of solvent on the excited state spectra and dynamics. The data are interpreted using molecular orbital theory which shows that the ICT state develops via mixing of the low-lying S(1) (2(1)A(g)-like) and S(2) (1(1)B(u)-like) excited singlet states to form a resultant state that preferentially evolves in polar solvent and exhibits a very large (∼25 D) dipole moment. Molecular dynamics calculations demonstrate that the features of the ICT state are present in ∼20 fs.

Stark Absorption Spectroscopy of Peridinin and Allene-Modified Analogues

Stark absorption spectra of peridinin (Per) and five allene-modified analogues and their angular dependence as a function of an externally applied electric field were measured in methyl methacrylate polymer at 77K. In all cases, the energetically lowest absorption band has a significant change of static dipole moment upon photoexcitation (Δμ). In particular, Per has the largest value of |Δμ|. The angles between Δμ and the transition dipole moment of all the analogues were determined. It is suggested that the allene group in Per plays a key role as the electron donor in the charge transfer process following photoexcitation. The results of MNDO-PSDCI calculations support this idea.

Spectral tuning in photoactive yellow protein by modulation of the shape of the excited state energy surface

Protein-chromophore interactions in photoreceptors often shift the chromophore absorbance maximum to a biologically relevant spectral region. A fundamental question regarding such spectral tuning effects is how the electronic ground state S(0) and excited state S(1) are modified by the protein. It is widely assumed that changes in energy gap between S(0) and S(1) are the main factor in biological spectral tuning. We report a generally applicable approach to determine if a specific residue modulates the energy gap, or if it alters the equilibrium nuclear geometry or width of the energy surfaces. This approach uses the effects that changes in these three parameters have on the absorbance and fluorescence emission spectra of mutants. We apply this strategy to a set of mutants of photoactive yellow protein (PYP) containing all 20 side chains at active site residue 46. While the mutants exhibit significant variation in both the position and width of their absorbance spectra, the fluorescence emission spectra are largely unchanged. This provides strong evidence against a major role for changes in energy gap in the spectral tuning of these mutants and reveals a change in the width of the S(1) energy surface. We determined the excited state lifetime of selected mutants and the observed correlation between the fluorescence quantum yield and lifetime shows that the fluorescence spectra are representative of the energy surfaces of the mutants. These results reveal that residue 46 tunes the absorbance spectrum of PYP largely by modulating the width of the S(1) energy surface.

Long-distance proton transfer with a break in the bacteriorhodopsin active site

Bacteriorhodopsin is a proton-pumping membrane protein found in the plasma membrane of the archaeon Halobacterium salinarium. Light-induced isomerization of the retinal chromophore from all-trans to 13-cis leads to a sequence of five conformation-coupled proton transfer steps and the net transport of one proton from the cytoplasmic to the extracellular side of the membrane. The mechanism of the long-distance proton transfer from the primary acceptor Asp85 to the extracellular proton release group during the O --> bR is poorly understood. Experiments suggest that this long-distance transfer could involve a transient state [O] in which the proton resides on the intermediate carrier Asp212. To assess whether the transient protonation of Asp212 participates in the deprotonation of Asp85, we performed hybrid Quantum Mechanics/Molecular Mechanics proton transfer calculations using different protein structures and with different retinal geometries and active site water molecules. The structural models were assessed by computing UV-vis excitation energies and C=O vibrational frequencies. The results indicate that a transient [O] conformer with protonated Asp212 could indeed be sampled during the long-distance proton transfer to the proton release group. Our calculations suggest that, in the starting proton transfer state O, the retinal is strongly twisted and at least three water molecules are present in the active site.

Ultrafast Time-resolved Absorption Spectroscopy of Geometric Isomers of Carotenoids

The structures of a number of stereoisomers of carotenoids have been revealed in three-dimensional X-ray crystallographic investigations of pigment-protein complexes from photosynthetic organisms. Despite these structural elucidations, the reason for the presence of stereoisomers in these systems is not well understood. An important unresolved issue is whether the natural selection of geometric isomers of carotenoids in photosynthetic pigment-protein complexes is determined by the structure of the protein binding site or by the need for the organism to accomplish a specific physiological task. The association of cis isomers of a carotenoid with reaction centers and trans isomers of the same carotenoid with light-harvesting pigment-protein complexes has led to the hypothesis that the stereoisomers play distinctly different physiological roles. A systematic investigation of the photophysics and photochemistry of purified, stable geometric isomers of carotenoids is needed to understand if a relationship between stereochemistry and biological function exists. In this work we present a comparative study of the spectroscopy and excited state dynamics of cis and trans isomers of three different open-chain carotenoids in solution. The molecules are neurosporene (n=9), spheroidene (n=10), and spirilloxanthin (n=13), where n is the number of conjugated pi-electron double bonds. The spectroscopic experiments were carried out on geometric isomers of the carotenoids purified by high performance liquid chromatography (HPLC) and then frozen to 77 K to inhibit isomerization. The spectral data taken at 77 K provide a high resolution view of the spectroscopic differences between geometric isomers. The kinetic data reveal that the lifetime of the lowest excited singlet state of a cis-isomer is consistently shorter than that of its corresponding all-trans counterpart despite the fact that the excited state energy of the cis molecule is typically higher than that of the trans molecule. Quantum theoretical calculations on an n=9 linear polyene were carried out to examine this process. The calculations indicate that the electronic coupling terms are significantly higher for the cis isomer, and when combined with the Franck-Condon factors, predict internal conversion rates roughly double those of the all-trans species. The electronic effects more than offset the decrease in coupling efficiencies associated with the higher system origin energies and explain the observed shorter cis isomer lifetimes.

Water-mediated spectral shifts in rhodopsin and bathorhodopsin

The role of water molecules in spectral tuning of proteins has been left largely unexplored. This topic is important because changing hydrogen bond patterns during the activation process may lead to spectral shifts which can be of diagnostic value for the underlying structures. Arguments put forward in this article are based on spectral shift calculations of the rhodopsin and bathorhodopsin chromophore due to wat2a and 2b in the presence and absence of the counterion and of the amino acids lining the rhodopsin binding pocket. They show, among others, that a single water molecule can shift the absorbance by up to 0.1 eV or 34 nm depending on the environment of the chromophore.

The charge-transfer properties of the S2 state of fucoxanthin in solution and in fucoxanthin chlorophyll-a/c2 protein (FCP) based on stark spectroscopy and molecular-orbital theory

Fucoxanthin chlorophyll-a/c 2 protein (FCP), the membrane-intrinsic light harvesting complex from the diatom Cyclotella meneghiniana, is characterized by Stark spectroscopy to obtain a quantitative measure of the excited-state dipolar properties of the constituent pigments. The electro-optical properties of the carotenoid fucoxanthin (Fx), the primary light harvester in FCP, were determined from the Stark spectrum measured in a MeTHF glass (77 K) and compared to the results from electronic-structure calculations. On photon absorption by Fx, a 17 D change in the static dipole moment (|Delta mu|exp), and a somewhat larger |Delta mu|exp at the red edge, are measured for the S 0 --> S 2 (1 (1)A g (-)-like -->1 (1)B u *+-like) transition. The large change in dipole moment indicates that Fx undergoes photoinduced charge transfer (CT), and underscores the influence of the S 2 state on the polarity-dependent excited-state dynamics of Fx that has so far been attributed to, and discussed in terms of, the S 0 and the S 1/ICT states. MNDO-PSDCI and SACCI-CISD calculations indicate that the 1 (1)B u (*+)-like state intrinsically possesses a dipole moment much smaller than the 2 (1)A g (*-)-like state, suggesting that solvent fields promote the mixing of these two states and could account for the large dipole moments measured here for the S 0 --> S 2 transition. These CT properties of the 1 (1)B u (*+)-like state of Fx are further enhanced in the protein and underpin its photosynthetic capabilities for light harvesting and energy transfer (ET). In FCP, the CT properties of the Fx's vary according to the energetic position: between 450 and 500 nm there appear to be two sets of Fx's that exhibit |Delta mu| exp values on the order of 5 and 15 D, whereas the red-most Fx's, that are very efficient in ET to chlorophyll-a (Chl-a), exhibit strikingly large |Delta mu| exp values on the order of 40 D. Such magnitudes of |Delta mu| exp suggest a mechanism that enhances Coulombic coupling to promote ET from the S 2 state of Fx to Chl-a. These three sets of Fx's, including a fourth red Fx, are identified by fitting the Stark spectrum of FCP with the Stark spectrum of Fx in MeTHF. In contrast to the Fx's in the protein, the electrostatic properties of the Chl's in FCP are comparatively much smaller. Notably, for the Q y band of Chl-a, a |Delta mu| exp of 0.92 D and a change in polarizability ( Delta alpha exp) of 20 A (3), indicate that the Chl-a's are monomeric in nature and decoupled from each other.

Effect of polarization on the opsin shift in rhodopsins. 1. A combined QM/QM/MM model for bacteriorhodopsin and pharaonis sensory rhodopsin II

The optical and IR-spectroscopic properties of the protonated Schiff base of retinal are highly sensitive to the electrostatic environment. This feature makes retinal a useful probe to study structural differences and changes in rhodopsins. It also raises an interest to theoretically predict the spectroscopic response to mutation and structural evolution. Computational models appropriate for this purpose usually combine sophisticated quantum mechanical (QM) methods with molecular mechanics (MM) force fields. In an effort to test and improve the accuracy of these QM/MM models, we consider in this article the effects of polarization and inter-residual charge transfer within the binding pocket of bacteriorhodopsin (bR) and pharaonis sensory rhodopsin II (psRII, also called pharaonis phoborhodopsin, ppR) on the excitation energy using an ab initio QM/QM/MM approach. The results will serve as reference for assessing empirical polarization models in a consecutive article.

Investigating the mechanisms of photosynthetic proteins using continuum electrostatics

Computational methods based on continuum electrostatics are widely used in theoretical biochemistry to analyze the function of proteins. Continuum electrostatic methods in combination with quantum chemical and molecular mechanical methods can help to analyze even very complex biochemical systems. In this article, applications of these methods to proteins involved in photosynthesis are reviewed. After giving a short introduction to the basic concepts of the continuum electrostatic model based on the Poisson-Boltzmann equation, we describe the application of this approach to the docking of electron transfer proteins, to the comparison of isofunctional proteins, to the tuning of absorption spectra, to the analysis of the coupling of electron and proton transfer, to the analysis of the effect of membrane potentials on the energetics of membrane proteins, and to the kinetics of charge transfer reactions. Simulations as those reviewed in this article help to analyze molecular mechanisms on the basis of the structure of the protein, guide new experiments, and provide a better and deeper understanding of protein functions.

Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

We have investigated geometries and excitation energies of bovine rhodopsin and some of its mutants by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in ONIOM scheme, employing B3LYP and BLYP density functionals as well as DFTB method for the QM part and AMBER force field for the MM part. QM/MM geometries of the protonated Schiff-base 11- cis-retinal with B3LYP and DFTB are very similar to each other. TD-B3LYP/MM excitation energy calculations reproduce the experimental absorption maximum of 500 nm in the presence of native rhodopsin environment and predict spectral shifts due to mutations within 10 nm, whereas TD-BLYP/MM excitation energies have red-shift error of at least 50 nm. In the wild-type rhodopsin, Glu113 shifts the first excitation energy to blue and accounts for most of the shift found. Other amino acids individually contribute to the first excitation energy but their net effect is small. The electronic polarization effect is essential for reproducing experimental bond length alternation along the polyene chain in protonated Schiff-base retinal, which correlates with the computed first excitation energy. It also corrects the excitation energies and spectral shifts in mutants, more effectively for deprotonated Schiff-base retinal than for the protonated form. The protonation state and conformation of mutated residues affect electronic spectrum significantly. The present QM/MM calculations estimate not only the experimental excitation energies but also the source of spectral shifts in mutants.

Spectral tuning of deep red cone pigments

Visual pigments are G-protein-coupled receptors that provide a critical interface between organisms and their external environment. Natural selection has generated vertebrate pigments that absorb light from the far-UV (360 nm) to the deep red (630 nm) while using a single chromophore, in either the A1 (11- cis-retinal) or A2 (11- cis-3,4-dehydroretinal) form. The fact that a single chromophore can be manipulated to have an absorption maximum across such an extended spectral region is remarkable. The mechanisms of wavelength regulation remain to be fully revealed, and one of the least well-understood mechanisms is that associated with the deep red pigments. We investigate theoretically the hypothesis that deep red cone pigments select a 6- s- trans conformation of the retinal chromophore ring geometry. This conformation is in contrast to the 6- s- cis ring geometry observed in rhodopsin and, through model chromophore studies, the vast majority of visual pigments. Nomographic spectral analysis of 294 A1 and A2 cone pigment literature absorption maxima indicates that the selection of a 6- s- trans geometry red shifts M/LWS A1 pigments by approximately 1500 cm (-1) ( approximately 50 nm) and A2 pigments by approximately 2700 cm (-1) ( approximately 100 nm). The homology models of seven cone pigments indicate that the deep red cone pigments select 6- s- trans chromophore conformations primarily via electrostatic steering. Our results reveal that the generation of a 6- s- trans conformation not only achieves a significant red shift but also provides enhanced stability of the chromophore within the deep red cone pigment binding sites.

Quantum mechanical/molecular mechanical studies on spectral tuning mechanisms of visual pigments and other photoactive proteins

The protein environments surrounding the retinal tune electronic absorption maximum from 350 to 630 nm. Hybrid quantum mechanical/molecular mechanical (QM/MM) methods can be used in calculating excitation energies of retinal in its native protein environments and in studying the molecular basis of spectral tuning. We hereby review recent QM/MM results on the phototransduction of bovine rhodopsin, bacteriorhodopsin, sensory rhodopsin II, nonretinal photoactive yellow protein and their mutants.

An ion-responsive motif in the second transmembrane segment of rhodopsin-like receptors

A L(M)xxxD(N, E) motif (x=a non-ionic amino acid residue, most frequently A, S, L or F; small capitals indicating a minor representation) is found in the second transmembrane (tm2) segment of most G-protein coupling metazoan receptors of the rhodopsin family (Rh-GPCRs). Changes in signal transduction, agonist binding and receptor cycling are known for numerous receptors bearing evolved or experimentally introduced mutations in this tm2 motif, especially of its aspartate residue. The [Na(+)] sensitivity of the receptor-agonist interaction relates to this aspartate in a number of Rh-GPCRs. Native non-conservative mutations in the tm2 motif only rarely coincide with significant changes in two other ubiquitous features of the rhodopsin family, the seventh transmembrane N(D)PxxY(F) motif and the D(E)RY(W,F) or analogous sequence at the border of the third transmembrane helix and the second intracellular loop. Native tm2 mutations with Rh-GPCRs frequently result in constitutive signaling, and with visual opsins also in shifts to short-wavelength sensitivity. Substitution of a strongly basic residue for the tm2 aspartate in Taste-2 receptors could be connected to a lack of sodium sensing by these receptors. These properties could be consistent with ionic interactions, and even of ion transfer, that involve the tm2 motif. A decrease in cation sensing by this motif is usually connected to an enhanced constitutive interaction of the mutated receptors with cognate G- proteins, and also relates to both the constitutive and the overall activity of the short-wavelength opsins.

Accurate evaluation of the absorption maxima of retinal proteins based on a hybrid QM/MM method

Here we improved our hybrid QM/MM methodology (Houjou et al. J Phys Chem B 2001, 105, 867) for evaluating the absorption maxima of photoreceptor proteins. The renewed method was applied to evaluation of the absorption maxima of several retinal proteins and photoactive yellow protein. The calculated absorption maxima were in good agreement with the corresponding experimental data with a computational error of <10 nm. In addition, our calculations reproduced the experimental gas-phase absorption maxima of model chromophores (protonated all-trans retinal Schiff base and deprotonated thiophenyl-p-coumarate) with the same accuracy. It is expected that our methodology allows for definitive interpretation of the spectral tuning mechanism of retinal proteins.

Ultrafast dynamics and excited state spectra of open-chain carotenoids at room and low temperatures

Many of the spectroscopic features and photophysical properties of carotenoids are explained using a three-state model in which the strong visible absorption of the molecules is associated with an S0 (1(1)Ag-) --> S2 (1(1)Bu+) transition, and the lowest lying singlet state, S1 (2(1)Ag-), is a state into which absorption from the ground state is forbidden by symmetry. However, semiempirical and ab initio quantum calculations have suggested additional excited singlet states may lie either between or in the vicinity of S1 (2(1)Ag-) and S2 (1(1)Bu+), and some ultrafast spectroscopic studies have reported evidence for these states. One such state, denoted S*, has been implicated as an intermediate in the depopulation of S2 (1(1)Bu+) and as a pathway for the formation of carotenoid triplet states in light-harvesting complexes. In this work, we present the results of an ultrafast, time-resolved spectroscopic investigation of a series of open-chain carotenoids derived from photosynthetic bacteria and systematically increasing in their number of pi-electron carbon-carbon double bonds (n). The molecules are neurosporene (n = 9), spheroidene (n = 10), rhodopin glucoside (n = 11), rhodovibrin (n = 12), and spirilloxanthin (n = 13). The molecules were studied in acetone and CS2 solvents at room temperature. These experiments explore the effect of solvent polarity and polarizability on the spectroscopic and kinetic behavior of the molecules. The molecules were also studied in ether/isopentane/ethanol (EPA) glasses at 77 K, in which the spectral resolution is greatly enhanced. Analysis of the data using global fitting techniques has revealed the ultrafast dynamics of the excited states and spectral changes associated with their decay, including spectroscopic features not previously reported. The data are consistent with S* being identified with a twisted conformational structure, the yield of which is increased in molecules having longer pi-electron conjugations. In particular, for the longest molecule in the series, spirilloxanthin, the experiments and a detailed quantum computational analysis reveal the presence of two S* states associated with relaxed S1 (2(1)Ag-) conformations involving nearly planar 6-s-cis and 6-s-trans geometries. We propose that in polar solvents, the ground state of spirilloxanthin takes on a corkscrew conformation that generates a net solute dipole moment while decreasing the cavity formation energy. Upon excitation and relaxation into the S1 (2(1)Ag-) state, the polyene unravels and flattens into a more planar geometry with comparable populations of 6-s-trans and 6-s-cis conformations.

Development of Bacteriorhodopsin Analogues and Studies of Charge Separated Excited States in the Photoprocesses of Linear Polyenes†

Development of bacteriorhodopsin (bR) analogues employing chromophore substitution technique for the purpose of characterizing the binding site of bR and generating bR analogues with novel opto-electronic properties for applications as photoactive element in nanotechnical devices are described. Additionally, the photophysical and photochemical properties of variously substituted diarylpolyenes as models of photobiologically relevant linear polyenes are discussed. The role of charge separated dipolar excited states in the photoprocesses of linear polyenes is highlighted.

Binding, tuning and mechanical function of the 4-hydroxy-cinnamic acid chromophore in photoactive yellow protein

The bacterial photoreceptor protein photoactive yellow protein (PYP) covalently binds the chromophore 4-hydroxy coumaric acid, tuning (spectral) characteristics of this cofactor. Here, we study this binding and tuning using a combination of pointmutations and chromophore analogs. In all photosensor proteins studied to date the covalent linkage of the chromophore to the apoprotein is dispensable for light-induced catalytic activation. We analyzed the functional importance of the covalent linkage using an isosteric chromophore-protein variant in which the cysteine is replaced by a glycine residue and the chromophore by thiomethyl-p-coumaric acid (TMpCA). The model compound TMpCA is shown to weakly complex with the C69G protein. This non-covalent binding results in considerable tuning of both the pKa and the color of the chromophore. The photoactivity of this system, however, was strongly impaired, making PYP the first known photosensor protein in which the covalent linkage of the chromophore is of paramount importance for the functional activity of the protein in vitro. We also studied the influence of chromophore analogs on the color and photocycle of PYP, not only in WT, but especially in the E46Q mutant, to test if effects from both chromophore and protein modifications are additive. When the E46Q protein binds the sinapinic acid chromophore, the color of the protein is effectively changed from yellow to orange. The altered charge distribution in this protein also results in a changed pKa value for chromophore protonation, and a strongly impaired photocycle. Both findings extend our knowledge of the photochemistry of PYP for signal generation.

The color of rhodopsins at the ab initio multiconfigurational perturbation theory resolution

We demonstrate that "brute force" quantum-mechanics/molecular-mechanics computations based on ab initio (i.e., first principles) multiconfigurational perturbation theory can reproduce the absorption maxima of a set of modified bovine rhodopsins with an accuracy allowing for the analysis of the factors determining their colors. In particular, we show that the theory accounts for the changes in excitation energy even when the proteins display the same charge distribution. Three color-tuning mechanisms, leading to changes of close magnitude, are demonstrated to operate in these conditions. The first is based on the change of the conformation of the conjugated backbone of the retinal chromophore. The second operates through the control of the distance between the positive charge residing on the chromophore and the carboxylate counterion. Finally, the third mechanism operates through the changes in orientation of the chromophore relative to the protein. These results offer perspectives for the unbiased computational design of mutants or chemically modified proteins with wanted optical properties.

Computational photochemistry of retinal proteins

High spectral tunability and quantum yield are the striking features of rhodopsin photochemistry. They rely on a strong and complex interaction of their chromophore, the protonated Schiff base of retinal, with its protein environment. In this article, we review the progress in the computational modeling of these systems, focusing on the optical properties and the excited state dynamics. While the earlier success of atomistic theoretical models was based on the breakthrough in X-ray crystallography and combined quantum mechanical molecular mechanical (QM/MM) methodology, recent advances point out the importance of high-level QM methods and the incorporation of effects that are neglected in conventional QM/MM or ONIOM schemes, like polarization and charge transfer.

RM1: A reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I

Twenty years ago, the landmark AM1 was introduced, and has since had an increasingly wide following among chemists due to its consistently good results and time-tested reliability--being presently available in countless computational quantum chemistry programs. However, semiempirical molecular orbital models still are of limited accuracy and need to be improved if the full potential of new linear scaling techniques, such as MOZYME and LocalSCF, is to be realized. Accordingly, in this article we present RM1 (Recife Model 1): a reparameterization of AM1. As before, the properties used in the parameterization procedure were: heats of formation, dipole moments, ionization potentials and geometric variables (bond lengths and angles). Considering that the vast majority of molecules of importance to life can be assembled by using only six elements: C, H, N, O, P, and S, and that by adding the halogens we can now build most molecules of importance to pharmaceutical research, our training set consisted of 1736 molecules, representative of organic and biochemistry, containing C, H, N, O, P, S, F, Cl, Br, and I atoms. Unlike AM1, and similar to PM3, all RM1 parameters have been optimized. For enthalpies of formation, dipole moments, ionization potentials, and interatomic distances, the average errors in RM1, for the 1736 molecules, are less than those for AM1, PM3, and PM5. Indeed, the average errors in kcal x mol(-1) of the enthalpies of formation for AM1, PM3, and PM5 are 11.15, 7.98, and 6.03, whereas for RM1 this value is 5.77. The errors, in Debye, of the dipole moments for AM1, PM3, PM5, and RM1 are, respectively, 0.37, 0.38, 0.50, and 0.34. Likewise, the respective errors for the ionization potentials, in eV, are 0.60, 0.55, 0.48, and 0.45, and the respective errors, in angstroms, for the interatomic distances are 0.036, 0.029, 0.037, and 0.027. The RM1 average error in bond angles of 6.82 degrees is only slightly higher than the AM1 figure of 5.88 degrees, and both are much smaller than the PM3 and PM5 figures of 6.98 degrees and 9.83 degrees, respectively. Moreover, a known error in PM3 nitrogen charges is corrected in RM1. Therefore, RM1 represents an improvement over AM1 and its similar successor PM3, and is probably very competitive with PM5, which is a somewhat different model, and not fully disclosed. RM1 possesses the same analytical construct and the same number of parameters for each atom as AM1, and, therefore, can be easily implemented in any software that already has AM1, not requiring any change in any line of code, with the sole exception of the values of the parameters themselves.

Ground and excited states of retinal schiff base chromophores by multiconfigurational perturbation theory

We have studied the wavelength dependence of retinal Schiff base absorbencies on the protonation state of the chromophore at the multiconfigurational level of theory using second order perturbation theory (CASPT2) within an atomic natural orbital basis set on MP2 optimized geometries. Quantitative agreement between calculated and experimental absorption maxima was obtained for protonated and deprotonated Schiff bases of all-trans- and 11-cis-retinal and intermediate states covering a wavelength range from 610 to 353 nm. These data will be useful as reference points for the calibration of more approximate schemes.