Theoretical prediction of optical absorption maxima for photosensory receptor mutants

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.

Fragment molecular orbital calculations on red fluorescent proteins (DsRed and mFruits)

We have performed a series of fragment molecular orbital (FMO) calculations for a family of red fluorescent proteins, DsRed and mFruits. The electronic transition energies were evaluated by the method of configuration interaction singles with perturbative doubles [CIS(D)] including higher-order corrections. The calculated values were in good agreement with the corresponding experimental peak values of spectra. Additionally, the chromophore environment was systematically analyzed in terms of the interaction energies between the pigment moiety and neighboring residues. It was theoretically revealed that the electrostatic interactions play a dominant role in the DsRed chromophore, whereas the color tunings in mFruits are controlled in a more delicate fashion.

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.

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.

Stress tensor analysis of the protein quake of photoactive yellow protein

Immediately after photon absorption, the photoenergy is converted to local stress energy via the ultrafast photoisomerization reaction of the p-coumaric acid (pCA) chromophore in a small water-soluble blue light receptor, photoactive yellow protein (PYP), derived from the halophilic bacterium, Halorhodospira halophila. A series of conformational changes are then induced, which are intimately related with the relaxation process on the energy landscape of PYP. In order to understand the signaling function of PYP in atomic detail, the characterization of the physical mechanism of the protein quake of PYP is important, as is the atomic description of the series of conformational changes associated with the photocycle. Here, we report a theoretical/computational study for the analysis of the intramolecular stress tensor for the dark state and three intermediate states, pR, pB1 and pB2, of PYP. As a result, we found that the magnitude of the stress released during the change from the pR to the pB1 state is significantly large at the hydroxyl oxygen atom of Tyr42, suggesting that this atom is the focus of the protein quake of PYP. This is consistent with previous experimental observations.

Spectral tuning of photoactive yellow protein

We report a theoretical study on the optical properties of a small, water-soluble photosensory receptor, photoactive yellow protein (PYP). A hierarchical ab initio molecular orbital calculation accurately evaluated the optical absorption maximum of the wild-type, as well as the lambda(max) values of 12 mutants. Electronic excitation of the chromophore directly affects the electronic state of nearby atoms in the protein environment. This effect is explicitly considered in the present study. Furthermore, the spectral tuning mechanism of PYP was investigated at the atomic level. The static disorder of a protein molecule is intimately related to the complex nature of its energy landscape. By using molecular dynamics simulation and quantum mechanical structure optimization, we obtained multiple minimum energy conformations of PYP. The statistical distribution of electronic excitation energies of these minima was compared with the hole-burning experiment (Masciangioli, T. [2000] Photochem. Photobiol. 72, 639), a direct observation of the distribution of excitation energies.