1
|
Lokiarchaeota archaeon schizorhodopsin-2 (LaSzR2) is an inward proton pump displaying a characteristic feature of acid-induced spectral blue-shift. Sci Rep 2020; 10:20857. [PMID: 33257762 PMCID: PMC7704677 DOI: 10.1038/s41598-020-77936-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2020] [Accepted: 10/29/2020] [Indexed: 11/25/2022] Open
Abstract
The photoreactive protein rhodopsin is widespread in microorganisms and has a variety of photobiological functions. Recently, a novel phylogenetically distinctive group named ‘schizorhodopsin (SzR)’ has been identified as an inward proton pump. We performed functional and spectroscopic studies on an uncharacterised schizorhodopsin from the phylum Lokiarchaeota archaeon. The protein, LaSzR2, having an all-trans-retinal chromophore, showed inward proton pump activity with an absorption maximum at 549 nm. The pH titration experiments revealed that the protonated Schiff base of the retinal chromophore (Lys188, pKa = 12.3) is stabilised by the deprotonated counterion (presumably Asp184, pKa = 3.7). The flash-photolysis experiments revealed the presence of two photointermediates, K and M. A proton was released and uptaken from bulk solution upon the formation and decay of the M intermediate. During the M-decay, the Schiff base was reprotonated by the proton from a proton donating residue (presumably Asp172). These properties were compared with other inward (SzRs and xenorhodopsins, XeRs) and outward proton pumps. Notably, LaSzR2 showed acid-induced spectral ‘blue-shift’ due to the protonation of the counterion, whereas outward proton pumps showed opposite shifts (red-shifts). Thus, we can distinguish between inward and outward proton pumps by the direction of the acid-induced spectral shift.
Collapse
|
2
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
3
|
Ernst OP, Lodowski DT, Elstner M, Hegemann P, Brown L, Kandori H. Microbial and animal rhodopsins: structures, functions, and molecular mechanisms. Chem Rev 2014; 114:126-63. [PMID: 24364740 PMCID: PMC3979449 DOI: 10.1021/cr4003769] [Citation(s) in RCA: 759] [Impact Index Per Article: 75.9] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2013] [Indexed: 12/31/2022]
Affiliation(s)
- Oliver P. Ernst
- Departments
of Biochemistry and Molecular Genetics, University of Toronto, 1 King’s College Circle, Medical Sciences Building, Toronto, Ontario M5S 1A8, Canada
| | - David T. Lodowski
- Center
for Proteomics and Bioinformatics, Case
Western Reserve University School of Medicine, 10900 Euclid Avenue, Cleveland, Ohio 44106, United States
| | - Marcus Elstner
- Institute
for Physical Chemistry, Karlsruhe Institute
of Technology, Kaiserstrasse
12, 76131 Karlsruhe, Germany
| | - Peter Hegemann
- Institute
of Biology, Experimental Biophysics, Humboldt-Universität
zu Berlin, Invalidenstrasse
42, 10115 Berlin, Germany
| | - Leonid
S. Brown
- Department
of Physics and Biophysics Interdepartmental Group, University of Guelph, 50 Stone Road East, Guelph, Ontario N1G 2W1, Canada
| | - Hideki Kandori
- Department
of Frontier Materials, Nagoya Institute
of Technology, Showa-ku, Nagoya 466-8555, Japan
| |
Collapse
|
4
|
Pal R, Sekharan S, Batista VS. Spectral Tuning in Halorhodopsin: The Chloride Pump Photoreceptor. J Am Chem Soc 2013; 135:9624-7. [DOI: 10.1021/ja404600z] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Rhitankar Pal
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut
06520-8107, United States
| | - Sivakumar Sekharan
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut
06520-8107, United States
| | - Victor S. Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut
06520-8107, United States
| |
Collapse
|
5
|
Sneskov K, Olsen JMH, Schwabe T, Hättig C, Christiansen O, Kongsted J. Computational screening of one- and two-photon spectrally tuned channelrhodopsin mutants. Phys Chem Chem Phys 2013; 15:7567-76. [DOI: 10.1039/c3cp44350g] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
6
|
Hasegawa JY, Fujimoto KJ, Kawatsu T. A Configuration Interaction Picture for a Molecular Environment Using Localized Molecular Orbitals: The Excited States of Retinal Proteins. J Chem Theory Comput 2012; 8:4452-61. [DOI: 10.1021/ct300510b] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jun-ya Hasegawa
- Fukui Institute for Fundamental
Chemistry, Kyoto University, 34-4 Takano-Nishihiraki-cho,
Sakyo-ku, Kyoto 606-8103, Japan
- Quantum Chemistry Research Institute, Kyodai Katsura Venture Plaza, Goryou
Oohara 1-36, Nishikyo-ku, Kyoto 615-8245, Japan
- Institute for Molecular Science, National Institute of Natural Science, 38 Nishigo-Naka,
Myodaiji, Okazaki 444-8585, Japan
| | - Kazuhiro J. Fujimoto
- Department of Computational Science,
Graduate School of System Informatics, Kobe University, 1-1, Rokkodai, Nada, Kobe 657-8501, Japan
| | - Tsutomu Kawatsu
- Fukui Institute for Fundamental
Chemistry, Kyoto University, 34-4 Takano-Nishihiraki-cho,
Sakyo-ku, Kyoto 606-8103, Japan
- Institute for Molecular Science, National Institute of Natural Science, 38 Nishigo-Naka,
Myodaiji, Okazaki 444-8585, Japan
| |
Collapse
|
7
|
Hayashi T, Matsuura A, Sato H, Sakurai M. Full-Quantum chemical calculation of the absorption maximum of bacteriorhodopsin: a comprehensive analysis of the amino acid residues contributing to the opsin shift. Biophysics (Nagoya-shi) 2012; 8:115-25. [PMID: 27493528 PMCID: PMC4629650 DOI: 10.2142/biophysics.8.115] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2012] [Accepted: 06/14/2012] [Indexed: 12/01/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Tomohiko Hayashi
- Center for Biological Resources and Informatics, Tokyo Institute of Technology, B-62 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan
| | - Azuma Matsuura
- Fujitsu Laboratories, Ltd., 10-1 Morinosato-Wakamiya, Atsugi 243-0197, Japan
| | - Hiroyuki Sato
- Fujitsu Laboratories, Ltd., 10-1 Morinosato-Wakamiya, Atsugi 243-0197, Japan
| | - Minoru Sakurai
- Center for Biological Resources and Informatics, Tokyo Institute of Technology, B-62 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan
| |
Collapse
|
8
|
Orekhov PS, Shaytan AK, Shaitan KV. Calculation of spectral shifts of the mutants of bacteriorhodopsin by QM/MM methods. Biophysics (Nagoya-shi) 2012. [DOI: 10.1134/s0006350912020170] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
9
|
Kaila VRI, Send R, Sundholm D. The effect of protein environment on photoexcitation properties of retinal. J Phys Chem B 2012; 116:2249-58. [PMID: 22166007 DOI: 10.1021/jp205918m] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Retinal is the photon absorbing chromophore of rhodopsin and other visual pigments, enabling the vertebrate vision process. The effects of the protein environment on the primary photoexcitation process of retinal were studied by time-dependent density functional theory (TDDFT) and the algebraic diagrammatic construction through second order (ADC(2)) combined with our recently introduced reduction of virtual space (RVS) approximation method. The calculations were performed on large full quantum chemical cluster models of the bluecone (BC) and rhodopsin (Rh) pigments with 165-171 atoms. Absorption wavelengths of 441 and 491 nm were obtained at the B3LYP level of theory for the respective models, which agree well with the experimental values of 414 and 498 nm. Electrostatic rather than structural strain effects were shown to dominate the spectral tuning properties of the surrounding protein. The Schiff base retinal and a neighboring Glu-113 residue were found to have comparable proton affinities in the ground state of the BC model, whereas in the excited state, the proton affinity of the Schiff base is 5.9 kcal/mol (0.26 eV) higher. For the ground and excited states of the Rh model, the proton affinity of the Schiff base is 3.2 kcal/mol (0.14 eV) and 7.9 kcal/mol (0.34 eV) higher than for Glu-113, respectively. The protein environment was found to enhance the bond length alternation (BLA) of the retinyl chain and blueshift the first absorption maxima of the protonated Schiff base in the BC and Rh models relative to the chromophore in the gas phase. The protein environment was also found to decrease the intensity of the second excited state, thus improving the quantum yield of the photoexcitation process. Relaxation of the BC model on the excited state potential energy surface led to a vanishing BLA around the isomerization center of the conjugated retinyl chain, rendering the retinal accessible for cis-trans isomerization. The energy of the relaxed excited state was found to be 30 kcal/mol (1.3 eV) above the minimum ground state energy, and might be related to the transition state of the thermal activation process.
Collapse
Affiliation(s)
- Ville R I Kaila
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, USA.
| | | | | |
Collapse
|
10
|
Melaccio F, Ferré N, Olivucci M. Quantum chemical modeling of rhodopsin mutants displaying switchable colors. Phys Chem Chem Phys 2012; 14:12485-95. [DOI: 10.1039/c2cp40940b] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
11
|
Hasegawa JY, Fujimoto KJ, Nakatsuji H. Color tuning in photofunctional proteins. Chemphyschem 2011; 12:3106-15. [PMID: 21990164 DOI: 10.1002/cphc.201100452] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2011] [Revised: 08/24/2011] [Indexed: 11/11/2022]
Abstract
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.
Collapse
Affiliation(s)
- Jun-ya Hasegawa
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan.
| | | | | |
Collapse
|
12
|
The application of quantum chemistry and condensed matter theory in studying amino-acids, protein folding and anticancer drug technology. Theor Chem Acc 2009. [DOI: 10.1007/s00214-009-0558-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
13
|
Wanko M, Hoffmann M, Frähmcke J, Frauenheim T, Elstner M. Effect of polarization on the opsin shift in rhodopsins. 2. Empirical polarization models for proteins. J Phys Chem B 2008; 112:11468-78. [PMID: 18729405 DOI: 10.1021/jp802409k] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The explicit treatment of polarization as a many-body interaction in condensed-phase systems represents a current problem in empirical force-field development. Although a variety of efficient models for molecular polarization have been suggested, polarizable force fields are still far from common use nowadays. In this work, we consider interactive polarization models employing Thole's short-range damping scheme and assess them for application on polypeptides. Despite the simplicity of the model, we find mean polarizabilities and anisotropies of amino acid side chains in excellent agreement with MP2/cc-pVQZ benchmark calculations. Combined with restrained electrostatic potential (RESP) derived atomic charges, the models are applied in a quantum-mechanical/molecular-mechanical (QM/MM) approach. An iterative scheme is used to establish a self-consistent mutual polarization between the QM and MM moieties. This ansatz is employed to study the influence of the protein polarizability on calculated optical properties of the protonated Schiff base of retinal in rhodopsin (Rh), bacterio-rhodopsin (bR), and pharaonis sensory rhodopsin II (psRII). The shifts of the excitation energy due to the instantaneous polarization response of the protein to the charge transfer on the retinal chromophore are quantified using the high level ab initio multireference spectroscopy-oriented configuration interaction (SORCI) method. The results are compared with those of previously published QM1/QM2/MM models for bR and psRII.
Collapse
|
14
|
Wanko M, Hoffmann M, Frauenheim T, Elstner M. Effect of polarization on the opsin shift in rhodopsins. 1. A combined QM/QM/MM model for bacteriorhodopsin and pharaonis sensory rhodopsin II. J Phys Chem B 2008; 112:11462-7. [PMID: 18698712 DOI: 10.1021/jp802408g] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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.
Collapse
|
15
|
Braun-Sand S, Sharma PK, Chu ZT, Pisliakov AV, Warshel A. The energetics of the primary proton transfer in bacteriorhodopsin revisited: it is a sequential light-induced charge separation after all. BIOCHIMICA ET BIOPHYSICA ACTA 2008; 1777:441-52. [PMID: 18387356 PMCID: PMC2443747 DOI: 10.1016/j.bbabio.2008.03.001] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/31/2007] [Revised: 02/29/2008] [Accepted: 03/03/2008] [Indexed: 11/26/2022]
Abstract
The light-induced proton transport in bacteriorhodopsin has been considered as a model for other light-induced proton pumps. However, the exact nature of this process is still unclear. For example, it is not entirely clear what the driving force of the initial proton transfer is and, in particular, whether it reflects electrostatic forces or other effects. The present work simulates the primary proton transfer (PT) by a specialized combination of the EVB and the QCFF/PI methods. This combination allows us to obtain sufficient sampling and a quantitative free energy profile for the PT at different protein configurations. The calculated profiles provide new insight about energetics of the primary PT and its coupling to the protein conformational changes. Our finding confirms the tentative analysis of an earlier work (A. Warshel, Conversion of light energy to electrostatic energy in the proton pump of Halobacterium halobium, Photochem. Photobiol. 30 (1979) 285-290) and determines that the overall PT process is driven by the energetics of the charge separation between the Schiff base and its counterion Asp85. Apparently, the light-induced relaxation of the steric energy of the chromophore leads to an increase in the ion-pair distance, and this drives the PT process. Our use of the linear response approximation allows us to estimate the change in the protein conformational energy and provides the first computational description of the coupling between the protein structural changes and the PT process. It is also found that the PT is not driven by twist-modulated changes of the Schiff base's pKa, changes in the hydrogen bond directionality, or other non-electrostatic effects. Overall, based on a consistent use of structural information as the starting point for converging free energy calculations, we conclude that the primary event should be described as a light-induced formation of an unstable ground state, whose relaxation leads to charge separation and to the destabilization of the ion-pair state. This provides the driving force for the subsequent PT steps.
Collapse
Affiliation(s)
- Sonja Braun-Sand
- Department of Chemistry, University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA 90089-1062, USA
- Department of Chemistry, University of Colorado at Colorado Springs (UCCS), Colorado Springs, CO 80918
| | - Pankaz K. Sharma
- Department of Chemistry, University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA 90089-1062, USA
| | - Zhen T. Chu
- Department of Chemistry, University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA 90089-1062, USA
| | - Andrei V. Pisliakov
- Department of Chemistry, University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA 90089-1062, USA
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, 418 SGM Building, 3620 McClintock Avenue, Los Angeles, CA 90089-1062, USA
| |
Collapse
|
16
|
Matsuura A, Sato H, Houjou H, Saito S, Hayashi T, Sakurai M. Accurate evaluation of the absorption maxima of retinal proteins based on a hybrid QM/MM method. J Comput Chem 2007; 27:1623-30. [PMID: 16900496 DOI: 10.1002/jcc.20432] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
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.
Collapse
Affiliation(s)
- Azuma Matsuura
- Fujitsu Laboratories Ltd., 10-1 Morinosato-Wakamiya, Atsugi 243-0197, Japan
| | | | | | | | | | | |
Collapse
|
17
|
Bondar AN, Suhai S, Fischer S, Smith JC, Elstner M. Suppression of the back proton-transfer from Asp85 to the retinal Schiff base in bacteriorhodopsin: A theoretical analysis of structural elements. J Struct Biol 2007; 157:454-69. [PMID: 17189704 DOI: 10.1016/j.jsb.2006.10.007] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2006] [Revised: 09/26/2006] [Accepted: 10/03/2006] [Indexed: 11/15/2022]
Abstract
The transfer of a proton from the retinal Schiff base to the nearby Asp85 protein group is an essential step in the directional proton-pumping by bacteriorhodopsin. To avoid the wasteful back reprotonation of the Schiff base from Asp85, the protein must ensure that, following Schiff base deprotonation, the energy barrier for back proton-transfer from Asp85 to the Schiff base is larger than that for proton-transfer from the Schiff base to Asp85. Here, three structural elements that may contribute to suppressing the back proton-transfer from Asp85 to the Schiff base are investigated: (i) retinal twisting; (ii) hydrogen-bonding distances in the active site; and (iii) the number and location of internal water molecules. The impact of the pattern of bond twisting on the retinal deprotonation energy is dissected by performing an extensive set of quantum-mechanical calculations. Structural rearrangements in the active site, such as changes of the Thr89:Asp85 distance and relocation of water molecules hydrogen-bonding to the Asp85 acceptor group, may participate in the mechanism which ensures that following the transfer of the Schiff base proton to Asp85 the protein proceeds with the subsequent photocycle steps, and not with back proton transfer from Asp85 to the Schiff base.
Collapse
Affiliation(s)
- Ana-Nicoleta Bondar
- Computational Molecular Biophysics, IWR, University of Heidelberg, Im Neuenheimer Feld 368, 69120 Heidelberg, Germany
| | | | | | | | | |
Collapse
|
18
|
Kloppmann E, Becker T, Ullmann GM. Electrostatic potential at the retinal of three archaeal rhodopsins: implications for their different absorption spectra. Proteins 2006; 61:953-65. [PMID: 16247786 DOI: 10.1002/prot.20744] [Citation(s) in RCA: 35] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The color tuning mechanism of the rhodopsin protein family has been in the focus of research for decades. However, the structural basis of the tuning mechanism in general and of the absorption shift between rhodopsins in particular remains under discussion. It is clear that a major determinant for spectral shifts between different rhodopsins are electrostatic interactions between the chromophore retinal and the protein. Based on the Poisson-Boltzmann equation, we computed and compared the electrostatic potential at the retinal of three archaeal rhodopsins: bacteriorhodopsin (BR), halorhodopsin (HR), and sensory rhodopsin II (SRII) for which high-resolution structures are available. These proteins are an excellent test case for understanding the spectral tuning of retinal. The absorption maxima of BR and HR are very similar, whereas the spectrum of SRII is considerably blue shifted--despite the structural similarity between these three proteins. In agreement with their absorption maxima, we find that the electrostatic potential is similar in BR and HR, whereas significant differences are seen for SRII. The decomposition of the electrostatic potential into contributions of individual residues, allowed us to identify seven residues that are responsible for the differences in electrostatic potential between the proteins. Three of these residues are located in the retinal binding pocket and have in fact been shown to account for part of the absorption shift between BR and SRII by mutational studies. One residue is located close to the beta-ionone ring of retinal and the remaining three residues are more than 8 A away from the retinal. These residues have not been discussed before, because they are, partly because of their location, no obvious candidates for the spectral shift among BR, HR, and SRII. However, their contribution to the differences in electrostatic potential is evident. The counterion of the Schiff base, which is frequently discussed to be involved in the spectral tuning, does not contribute to the dissimilarities between the electrostatic potentials.
Collapse
Affiliation(s)
- Edda Kloppmann
- Structural Biology/Bioinformatics, University of Bayreuth, Bayreuth, Germany
| | | | | |
Collapse
|
19
|
López CS, Faza ON, Estévez SL, de Lera AR. Computation of vertical excitation energies of retinal and analogs: Scope and limitations. J Comput Chem 2006; 27:116-23. [PMID: 16273505 DOI: 10.1002/jcc.20305] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
A comprehensive survey of computational methods: semiempirical (ZINDO/S), Time-Dependent Hartree-Fock (TD-HF), Configuration Interaction Singles (CIS), and several approximate functionals within the Time-Dependent Density Functional Theory (TD-DFT) has been carried out for the description of vertical excitation energies and oscillator strengths of retinal and related polyenals. ZINDO and TD-DFT computations showed the best agreement with the experimental data. In particular, hybrid functionals including approximately 25% of exact exchange (B3LYP, B3P86, and PBE0) were found to perform best with these highly conjugated polyenes. A systematic average error of 0.18-0.22 eV has been found after a simple one-parameter correction. Thus, 0.18 eV might be considered the upper limit of accuracy for current one-determinant methods in the computation of vertical excitation energies. The consideration of adiabatic excitations, conformational sampling, solvation, and nondynamic correlation should describe this processes more accurately, but this leads to highly demanding methods beyond feasibility for these large polyenes. The trends observed, particularly the good performance of the ZINDO/S method, should pave the way for the prediction of excited states properties in natural and artificial photoreceptor proteins, thus advancing towards the description of their light-transducing biological role in Nature.
Collapse
Affiliation(s)
- Carlos Silva López
- Departamento de Química Orgánica, Facultade de Química, Universidade de Vigo, Lagoas Marcosende, E-36310, Vigo, Galicia, Spain
| | | | | | | |
Collapse
|
20
|
Affiliation(s)
- Jacopo Tomasi
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa, Italy.
| | | | | |
Collapse
|
21
|
Wanko M, Hoffmann M, Strodel P, Koslowski A, Thiel W, Neese F, Frauenheim T, Elstner M. Calculating Absorption Shifts for Retinal Proteins: Computational Challenges. J Phys Chem B 2005; 109:3606-15. [PMID: 16851399 DOI: 10.1021/jp0463060] [Citation(s) in RCA: 202] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Rhodopsins can modulate the optical properties of their chromophores over a wide range of wavelengths. The mechanism for this spectral tuning is based on the response of the retinal chromophore to external stress and the interaction with the charged, polar, and polarizable amino acids of the protein environment and is connected to its large change in dipole moment upon excitation, its large electronic polarizability, and its structural flexibility. In this work, we investigate the accuracy of computational approaches for modeling changes in absorption energies with respect to changes in geometry and applied external electric fields. We illustrate the high sensitivity of absorption energies on the ground-state structure of retinal, which varies significantly with the computational method used for geometry optimization. The response to external fields, in particular to point charges which model the protein environment in combined quantum mechanical/molecular mechanical (QM/MM) applications, is a crucial feature, which is not properly represented by previously used methods, such as time-dependent density functional theory (TDDFT), complete active space self-consistent field (CASSCF), and Hartree-Fock (HF) or semiempirical configuration interaction singles (CIS). This is discussed in detail for bacteriorhodopsin (bR), a protein which blue-shifts retinal gas-phase excitation energy by about 0.5 eV. As a result of this study, we propose a procedure which combines structure optimization or molecular dynamics simulation using DFT methods with a semiempirical or ab initio multireference configuration interaction treatment of the excitation energies. Using a conventional QM/MM point charge representation of the protein environment, we obtain an absorption energy for bR of 2.34 eV. This result is already close to the experimental value of 2.18 eV, even without considering the effects of protein polarization, differential dispersion, and conformational sampling.
Collapse
Affiliation(s)
- M Wanko
- Department of Theoretical Physics, University of Paderborn, D-33098 Paderborn, Germany
| | | | | | | | | | | | | | | |
Collapse
|
22
|
Bondar AN, Fischer S, Smith JC, Elstner M, Suhai S. Key role of electrostatic interactions in bacteriorhodopsin proton transfer. J Am Chem Soc 2004; 126:14668-77. [PMID: 15521787 DOI: 10.1021/ja047982i] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The first proton transport step following photon absorption in bacteriorhodopsin is from the 13-cis retinal Schiff base to Asp85. Configurational and energetic determinants of this step are investigated here by performing quantum mechanical/molecular mechanical minimum-energy reaction-path calculations. The results suggest that retinal can pump protons when in the 13-cis, 15-anti conformation but not when 13-cis, 15-syn. Decomposition of the proton transfer energy profiles for various possible pathways reveals a conflict between the effect of the intrinsic proton affinities of the Schiff base and Asp85, which favors the neutral, product state (i.e., with Asp85 protonated), with the mainly electrostatic interaction between the protein environment with the reacting partners, which favors the ion pair reactant state (i.e., with retinal protonated). The rate-limiting proton-transfer barrier depends both on the relative orientations of the proton donor and acceptor groups and on the pathway followed by the proton; depending on these factors, the barrier may arise from breaking and forming of hydrogen bonds involving the Schiff base, Asp85, Asp212, and water w402, and from nonbonded interactions involving protein groups that respond to the charge rearrangements in the Schiff base region.
Collapse
Affiliation(s)
- Ana-Nicoleta Bondar
- Computational Molecular Biophysics, IWR, Heidelberg University, Im Neuenheimer Feld 368, Heidelberg, Germany.
| | | | | | | | | |
Collapse
|
23
|
Howard IA, Klein DJ, March NH, Van Alsenoy C, Suhai S, Jánosvalfi Z, Nagy Á. Change in Electronic Structure of Polyenes Due to Interaction with Polyacenes and with Graphitic Strips. J Phys Chem B 2004. [DOI: 10.1021/jp040079k] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- I. A. Howard
- Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium, Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, Oxford University, Oxford, England, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium, Deutsches Krebsforschungszentrum (DKFZ), Molekulare Biophysik, D-60120 Heidelberg, Germany, and Department of Theoretical Physics, University of Debrecen, H-4010 Debrecen, Hungary
| | - D. J. Klein
- Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium, Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, Oxford University, Oxford, England, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium, Deutsches Krebsforschungszentrum (DKFZ), Molekulare Biophysik, D-60120 Heidelberg, Germany, and Department of Theoretical Physics, University of Debrecen, H-4010 Debrecen, Hungary
| | - N. H. March
- Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium, Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, Oxford University, Oxford, England, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium, Deutsches Krebsforschungszentrum (DKFZ), Molekulare Biophysik, D-60120 Heidelberg, Germany, and Department of Theoretical Physics, University of Debrecen, H-4010 Debrecen, Hungary
| | - C. Van Alsenoy
- Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium, Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, Oxford University, Oxford, England, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium, Deutsches Krebsforschungszentrum (DKFZ), Molekulare Biophysik, D-60120 Heidelberg, Germany, and Department of Theoretical Physics, University of Debrecen, H-4010 Debrecen, Hungary
| | - S. Suhai
- Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium, Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, Oxford University, Oxford, England, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium, Deutsches Krebsforschungszentrum (DKFZ), Molekulare Biophysik, D-60120 Heidelberg, Germany, and Department of Theoretical Physics, University of Debrecen, H-4010 Debrecen, Hungary
| | - Z. Jánosvalfi
- Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium, Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, Oxford University, Oxford, England, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium, Deutsches Krebsforschungszentrum (DKFZ), Molekulare Biophysik, D-60120 Heidelberg, Germany, and Department of Theoretical Physics, University of Debrecen, H-4010 Debrecen, Hungary
| | - Á. Nagy
- Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium, Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, Oxford University, Oxford, England, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Antwerp, Belgium, Deutsches Krebsforschungszentrum (DKFZ), Molekulare Biophysik, D-60120 Heidelberg, Germany, and Department of Theoretical Physics, University of Debrecen, H-4010 Debrecen, Hungary
| |
Collapse
|
24
|
Wada M, Sakurai M. A quantum chemical method for rapid optimization of protein structures. J Comput Chem 2004; 26:160-8. [PMID: 15586398 DOI: 10.1002/jcc.20154] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
A quantum chemical method for rapid optimization of protein structures is proposed. In this method, a protein structure is treated as an assembly of amino acid units, and the geometry optimization of each unit is performed with taking the effect of its surrounding environment into account. The optimized geometry of a whole protein is obtained by repeated application of such a local optimization procedure over the entire part of the protein. Here, we implemented this method in the MOPAC program and performed geometry optimization for three different sizes of proteins. Consequently, these results demonstrate that the total energies of the proteins are much efficiently minimized compared with the use of conventional optimization methods, including the MOZYME algorithm (a representative linear-scaling method) with the BFGS routine. The proposed method is superior to the conventional methods in both CPU time and memory requirements.
Collapse
Affiliation(s)
- Mitsuhito Wada
- Makuhari R&D Center, Celestar Lexico-Sciences, Inc., Makuhari Techno Garden D17, 1-3 Nakase, Mihama-ku, Chiba 261-8501, Japan
| | | |
Collapse
|