Absorption wavelength along chromophore low-barrier hydrogen bonds

Molecular origins of absorption wavelength variation among phycocyanobilin-binding proteins

Phycocyanobilin (PCB)-binding proteins, including cyanobacteriochromes and phytochromes, function as photoreceptors and exhibit a wide range of absorption maximum wavelengths. To elucidate the color-tuning mechanisms among these proteins, we investigated seven crystal structures of six PCB-binding proteins: Anacy_2551g3, AnPixJg2, phosphorylation-responsive photosensitive histidine kinase, RcaE, Sb.phyB(PG)-PCB, and Slr1393g3. Employing a quantum chemical/molecular mechanical approach combined with a polarizable continuum model, our analysis revealed that differences in absorption wavelengths among PCB-binding proteins primarily arise from variations in the shape of the PCB molecule itself, accounting for a ∼150 nm difference. Remarkably, calculated excitation energies sufficiently reproduced the absorption wavelengths of these proteins spanning ∼200 nm, including 728 nm for Anacy_2551g3. However, assuming the hypothesized lactim conformation resulted in a significant deviation from the experimentally measured absorption wavelength for Anacy_2551g3. The significantly red-shifted absorption wavelength of Anacy_2551g3 can unambiguously be explained by the significant overlap of molecular orbitals between the two pyrrole rings at both edges of the PCB chromophore without the need to hypothesize lactim formation.

Stretching vibrational frequencies and pKa differences in H-bond networks of protein environments

The experimentally measured stretching vibrational frequencies of O-D [νO-D(donor)] and C=O [νC=O(donor)] H-bond donor groups can provide valuable information about the H-bonds in proteins. Here, using a quantum mechanical/molecular mechanical approach, the relationship between these vibrational frequencies and the difference in pKa values between H-bond donor and acceptor groups [ΔpKa(donor … acceptor)] in bacteriorhodopsin and photoactive yellow protein environments was investigated. The results show that νO-D(donor) is correlated with ΔpKa(donor … acceptor), regardless of the specific protein environment. νC=O(donor) is also correlated with ΔpKa(donor … acceptor), although the correlation is weak because the C=O bond does not have a proton. Importantly, the shifts in νO-D(donor) and νC=O(donor) are not caused by changes in pKa(donor) alone, but rather by changes in ΔpKa(donor … acceptor). Specifically, a decrease in ΔpKa(donor … acceptor) can lead to proton release from the H-bond donor group toward the acceptor group, resulting in shifts in the vibrational frequencies of the protein environment. These findings suggest that changes in the stretching vibrational frequencies, in particular νO-D(donor), can be used to monitor proton transfer in protein environments.

Mechanism of Absorption Wavelength Shift Depending on the Protonation State of the Acrylate Group in Chlorophyll c

Diatoms can use light in the blue-green region because they have chlorophyll c (Chlc) in light-harvesting antenna proteins, fucoxanthin and chlorophyll a/c-binding protein (FCP). Chlc has a protonatable acrylate group (-CH═CH-COOH/COO-) conjugated to the porphyrin ring. As the absorption wavelength of Chlc changes upon the protonation of the acrylate group, Chlc is a candidate component that is responsible for photoprotection in diatoms, which switches the FCP function between light-harvesting and energy-dissipation modes depending on the light intensity. Here, we investigate the mechanism by which the absorption wavelength of Chlc changes owing to the change in the protonation state of the acrylate group, using a quantum mechanical/molecular mechanical approach. The calculated absorption wavelength of the Soret band of protonated Chlc is ∼25 nm longer than that of deprotonated Chlc, which is due to the delocalization of the lowest (LUMO) and second lowest (LUMO+1) unoccupied molecular orbitals toward the acrylate group. These results suggest that in FCP, the decrease in pH on the lumenal side under high-light conditions leads to protonation of Chlc and thereby a red shift in the absorption wavelength.

Mechanism of Absorption Wavelength Shift of Bacteriorhodopsin During Photocycle

Bacteriorhodopsin, a light-driven proton pump, alters the absorption wavelengths in the range of 410-617 nm during the photocycle. Here, we report the absorption wavelengths, calculated using 12 bacteriorhodopsin crystal structures (including the BR, BR_13-cis, J, K₀, K_E, K_L, L, M, N, and O state structures) and a combined quantum mechanical/molecular mechanical/polarizable continuum model (QM/MM/PCM) approach. The QM/MM/PCM calculations reproduced the experimentally measured absorption wavelengths with a standard deviation of 4 nm. The shifts in the absorption wavelengths can be explained mainly by the following four factors: (i) retinal Schiff base deformation/twist induced by the protein environment, leading to a decrease in the electrostatic interaction between the protein environment and the retinal Schiff base; (ii) changes in the protonation state of the protein environment, directly altering the electrostatic interaction between the protein environment and the retinal Schiff base; (iii) changes in the protonation state; or (iv) isomerization of the retinal Schiff base, where the absorption wavelengths of the isomers originally differ.

Proton transfer and conformational changes along the hydrogen bond network in heliorhodopsin

Heliorhodopsin releases a proton from the Schiff base during the L-state to M-state transition but not toward the protein bulk surface. Here we investigate proton transfer and induced structural changes along the H-bond network in heliorhodopsin using a quantum mechanical/molecular mechanical approach and molecular dynamics simulations. Light-induced proton transfer could occur from the Schiff base toward Glu107, reorienting Ser76, followed by subsequent proton transfer toward His80. His80 protonation induces the reorientation of Trp246 on the extracellular surface, originating from the electrostatic interaction that propagates along the transmembrane H-bond network [His80…His23…H₂O_[H23/Q26]…Gln26…Trp246] over a distance of 15 Å. Furthermore, it induces structural fluctuation on the intracellular side in the H-bond network [His80…Asn16…Tyr92…Glu230…Arg104…Glu149], opening the inner cavity at the Tyr92 moiety. These may be a basis of how light-induced proton transfer causes conformational changes during the M-state to O-state transition.

Proton transfer reactions in donor site mutants of ESR, a retinal protein from Exiguobacterium sibiricum

Light-driven proton transport by microbial retinal proteins such as archaeal bacteriorhodopsin involves carboxylic residues as internal proton donors to the catalytic center which is a retinal Schiff base (SB). The proton donor, Asp96 in bacteriorhodopsin, supplies a proton to the transiently deprotonated Schiff base during the photochemical cycle. Subsequent proton uptake resets the protonated state of the donor. This two step process became a distinctive signature of retinal based proton pumps. Similar steps are observed also in many natural variants of bacterial proteorhodopsins and xanthorhodopsins where glutamic acid residues serve as a proton donor. Recently, however, an exception to this rule was found. A retinal protein from Exiguobacterium sibiricum, ESR, contains a Lys residue in place of Asp or Glu, which facilitates proton transfer from the bulk to the SB. Lys96 can be functionally replaced with the more common donor residues, Asp or Glu. Proton transfer to the SB in the mutants containing these replacements (K96E and K96D/A47T) is much faster than in the proteins lacking the proton donor (K96A and similar mutants), and in the case of K96D/A47T, comparable with that in the wild type, indicating that carboxylic residues can replace Lys96 as proton donors in ESR. We show here that there are important differences in the functioning of these residues in ESR from the way Asp96 functions in bacteriorhodopsin. Reprotonation of the SB and proton uptake from the bulk occur almost simultaneously during the M to N transition (as in the wild type ESR at neutral pH), whereas in bacteriorhodopsin these two steps are well separated in time and occur during the M to N and N to O transitions, respectively.

Release of a Proton and Formation of a Low-Barrier Hydrogen Bond between Tyrosine D and D2-His189 in Photosystem II

In photosystem II (PSII), the second-lowest oxidation state (S₁) of the oxygen-evolving Mn₄CaO₅ cluster is the most stable, as the radical form of the redox-active D2-Tyr160 is considered to be a candidate that accepts an electron from the lowest oxidation state (S₀) in the dark. Using quantum mechanical/molecular mechanical calculations, we investigated the redox potential (E _m) of TyrD and its H-bond partner, D2-His189. The potential energy profile indicates that the release of a proton from the TyrD...D2-His189 pair leads to the formation of a low-barrier H-bond. The E _m depends on the H⁺ position along the low-barrier H-bond, e.g., 680 mV when the H⁺ is at the D2-His189 moiety and 800 mV when the H⁺ is at the TyrD moiety, which can explain why TyrD mediates both the S₀ to S₁ oxidation and the S₂ to S₁ reduction.