1
|
Choudhury A, Santra S, Ghosh D. Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment. J Chem Theory Comput 2024; 20:4951-4964. [PMID: 38864715 DOI: 10.1021/acs.jctc.4c00027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.
Collapse
Affiliation(s)
- Arpan Choudhury
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Supriyo Santra
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Debashree Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
2
|
Singh NK, Ramamourthy B, Hage N, Kappagantu KM. Optogenetics: Illuminating the Future of Hearing Restoration and Understanding Auditory Perception. Curr Gene Ther 2024; 24:208-216. [PMID: 38676313 DOI: 10.2174/0115665232269742231213110937] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Revised: 10/07/2023] [Accepted: 10/25/2023] [Indexed: 04/28/2024]
Abstract
Hearing loss is a prevalent sensory impairment significantly affecting communication and quality of life. Traditional approaches for hearing restoration, such as cochlear implants, have limitations in frequency resolution and spatial selectivity. Optogenetics, an emerging field utilizing light-sensitive proteins, offers a promising avenue for addressing these limitations and revolutionizing hearing rehabilitation. This review explores the methods of introducing Channelrhodopsin- 2 (ChR2), a key light-sensitive protein, into cochlear cells to enable optogenetic stimulation. Viral- mediated gene delivery is a widely employed technique in optogenetics. Selecting a suitable viral vector, such as adeno-associated viruses (AAV), is crucial in efficient gene delivery to cochlear cells. The ChR2 gene is inserted into the viral vector through molecular cloning techniques, and the resulting viral vector is introduced into cochlear cells via direct injection or round window membrane delivery. This allows for the expression of ChR2 and subsequent light sensitivity in targeted cells. Alternatively, direct cell transfection offers a non-viral approach for ChR2 delivery. The ChR2 gene is cloned into a plasmid vector, which is then combined with transfection agents like liposomes or nanoparticles. This mixture is applied to cochlear cells, facilitating the entry of the plasmid DNA into the target cells and enabling ChR2 expression. Optogenetic stimulation using ChR2 allows for precise and selective activation of specific neurons in response to light, potentially overcoming the limitations of current auditory prostheses. Moreover, optogenetics has broader implications in understanding the neural circuits involved in auditory processing and behavior. The combination of optogenetics and gene delivery techniques provides a promising avenue for improving hearing restoration strategies, offering the potential for enhanced frequency resolution, spatial selectivity, and improved auditory perception.
Collapse
Affiliation(s)
- Namit Kant Singh
- Department of Otorhinolaryngology and Head and Neck Surgery, All India institute of Medical Sciences, Bibinagar, Hyderabad, India
| | - Balaji Ramamourthy
- Department of Otorhinolaryngology and Head and Neck Surgery, All India institute of Medical Sciences, Bibinagar, Hyderabad, India
| | - Neemu Hage
- Department of Otorhinolaryngology and Head and Neck Surgery, All India institute of Medical Sciences, Bibinagar, Hyderabad, India
| | - Krishna Medha Kappagantu
- Department of Otorhinolaryngology and Head and Neck Surgery, All India institute of Medical Sciences, Bibinagar, Hyderabad, India
| |
Collapse
|
3
|
Ricardi N, González-Espinoza CE, Adam S, Church JR, Schapiro I, Wesołowski TA. Embedding Nonrigid Solutes in an Averaged Environment: A Case Study on Rhodopsins. J Chem Theory Comput 2023; 19:5289-5302. [PMID: 37441785 PMCID: PMC10413860 DOI: 10.1021/acs.jctc.3c00285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2023]
Abstract
Many simulation methods concerning solvated molecules are based on the assumption that the solvated species and the solvent can be characterized by some representative structures of the solute and some embedding potential corresponding to this structure. While the averaging of the solvent configurations to obtain an embedding potential has been studied in great detail, this hinges on a single solute structure representation. This assumption is re-examined and generalized for conformationally flexible solutes and tested on 4 nonrigid systems. In this generalized approach, the solute is characterized by a set of representative structures and the corresponding embedding potentials. The representative structures are identified by means of subdividing the statistical ensemble, which in this work is generated by a constant-temperature molecular dynamics simulation. The embedding potential defined in the Frozen-Density Embedding Theory is used to characterize the average effect of the solvent in each subensemble. The numerical examples concern the vertical excitation energies of protonated retinal Schiff bases in protein environments. It is comprehensively shown that subensemble averaging leads to huge computational savings compared with explicit averaging of the excitation energies in the whole ensemble while introducing only minor errors in the case of the systems examined.
Collapse
Affiliation(s)
- Niccolò Ricardi
- Department of Physical Chemistry, University of Geneva, 1205 Geneva, Switzerland
| | | | - Suliman Adam
- Fritz Haber Center for Molecular Dynamics, Hebrew University of Jerusalem Israel, 91904 Jerusalem, Israel
| | - Jonathan R Church
- Fritz Haber Center for Molecular Dynamics, Hebrew University of Jerusalem Israel, 91904 Jerusalem, Israel
| | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics, Hebrew University of Jerusalem Israel, 91904 Jerusalem, Israel
| | | |
Collapse
|
4
|
Xin Q, Zhang W, Yuan S. The Mechanism of the Channel Opening in Channelrhodopsin-2: A Molecular Dynamics Simulation. Int J Mol Sci 2023; 24:ijms24065667. [PMID: 36982741 PMCID: PMC10057421 DOI: 10.3390/ijms24065667] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 03/02/2023] [Accepted: 03/07/2023] [Indexed: 03/18/2023] Open
Abstract
Channelrhodopsin-2 (ChR2) has been one of the most important objects in the study of optogenetics. The retinal chromophore molecule absorbs photons and undergoes an isomerization reaction, which triggers the photocycle, resulting in a series of conformational changes. In this study, a series of intermediate structures (including D470, P500, P390-early, P390-late, and P520 states) of ChR2 in the photocycle were modeled, and molecular dynamics (MD) simulations were performed to elucidate the mechanism of ion channel opening of ChR2. The maximum absorption wavelength of these intermediates calculated by time-dependent density function theory (TD-DFT) is in general agreement with the experimental values, the distribution of water density gradually increases in the process of photocycle, and the radius of the ion channel is larger than 6 Å. All these results indicate that our structural models of the intermediates are reasonable. The evolution of protonation state of E90 during the photocycle is explained. E90 will deprotonate when the P390-early transforms into P390-late, in which the two conformations of P390-early and P390-late obtained from the simulations are consistent with the experimental descriptions. To validate the conductive P520 state, the potential mean force (PMF) of Na+ ions passing through the P520 intermediate was calculated by using steered molecular dynamics (SMD) simulation combined with umbrella sampling. The result shows that the Na+ ions passing through the channel with a very low energy barrier, especially in the central gate, is almost barrierless. This indicates that the channel is open in the P520 state.
Collapse
|
5
|
Gating and ion selectivity of Channelrhodopsins are critical for photo-activated orientation of Chlamydomonas as shown by in vivo point mutation. Nat Commun 2022; 13:7253. [PMID: 36433995 PMCID: PMC9700795 DOI: 10.1038/s41467-022-35018-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Accepted: 11/15/2022] [Indexed: 11/26/2022] Open
Abstract
The green unicellular alga Chlamydomonas reinhardtii with two photoreceptors called channelrhodopsins is a model organism that gave birth to a new scientific field of biomedical studies, optogenetics. Although channelrhodopsins are helping to decipher the activity of the human brain, their functionality has never been extensively studied in the organism of origin, mainly due to the difficulties connected to reverse genetic interventions. In this study, we present a CRISPR-Cas9-based technique that enables a precise in vivo exchange of single amino acids in a selected gene. To shed light on the function of channelrhodopsins ChR1 (C1) and ChR2 (C2) in vivo, we deleted both channelrhodopsins independently in the wild-type strain and introduced point mutations in the remaining channel, causing modified photocycle kinetics and ion selectivity. The mutated strains, ΔC1C2-E123T, ΔC1C2-E90R and ΔC1C2-E90Q, showed about 100-fold decrease in photosensitivity, a reduced photophobic response and faster light adaptation rates due to accelerated photocycle kinetics and reduced Ca2+ conductance. Moreover, the ΔC1C2-E90Q with an additionally reduced H+ permeability produced an electrical response only in the presence of Na+ ions, highlighting a contribution and importance of H+ conductance to photocurrents in the wild-type algae. Finally, in the ΔC1C2-E90R strain with the channelrhodopsin selectivity converted to anions, no photo-responses were detected. We conclude that the precise photocycle kinetics and the particular ion selectivity of channelrhodopsins are the key parameters for efficient phototaxis in low light conditions.
Collapse
|
6
|
Bhatt A, Mukhopadhyaya A, Ali ME. α-Helix in Cystathionine β-Synthase Enzyme Acts as an Electron Reservoir. J Phys Chem B 2022; 126:4754-4760. [PMID: 35687358 DOI: 10.1021/acs.jpcb.2c01657] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The modulation of electron density at the Pyridoxal 5'-phosphate (PLP) catalytic center, because of charge transfer across the α-helix/PLP interface, is the determining factor for the enzymatic activities in the human Cystathionine β-Synthase (hCBS) enzyme. Applying density functional theory calculations, in conjunction with the real space density analysis, we investigated the charge density delocalization across the entire heme-α-helix-PLP electron communication channels. The electron delocalization due to hydrogen bonds at the heme/α-helix and α-helix/PLP interfaces are found to be extended over a very long range, as a result of redistribution of electron densities of the cofactors. Moreover, the internal hydrogen bonds of α-helix that are crucial for its secondary structure also participate in the electron redistribution through the structured hydrogen-bond network. α-Helix is found to accumulate the electron density at the ground state from both of the cofactors and behaves as an electron reservoir for catalytic reaction at the electrophilic center of PLP.
Collapse
Affiliation(s)
- Aashish Bhatt
- Institute of Nano Science and Technology, Sector-81, Mohali, Punjab-140306, India
| | - Aritra Mukhopadhyaya
- Institute of Nano Science and Technology, Sector-81, Mohali, Punjab-140306, India
| | - Md Ehesan Ali
- Institute of Nano Science and Technology, Sector-81, Mohali, Punjab-140306, India
| |
Collapse
|
7
|
Church JR, Amoyal GS, Borin VA, Adam S, Olsen JMH, Schapiro I. Deciphering the Spectral Tuning Mechanism in Proteorhodopsin: The Dominant Role of Electrostatics Instead of Chromophore Geometry. Chemistry 2022; 28:e202200139. [PMID: 35307890 PMCID: PMC9325082 DOI: 10.1002/chem.202200139] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2022] [Indexed: 11/11/2022]
Abstract
Proteorhodopsin (PR) is a photoactive proton pump found in marine bacteria. There are two phenotypes of PR exhibiting an environmental adaptation to the ocean's depth which tunes their maximum absorption: blue‐absorbing proteorhodopsin (BPR) and green‐absorbing proteorhodopsin (GPR). This blue/green color‐shift is controlled by a glutamine to leucine substitution at position 105 which accounts for a 20 nm shift. Typically, spectral tuning in rhodopsins is rationalized by the external point charge model but the Q105L mutation is charge neutral. To study this tuning mechanism, we employed the hybrid QM/MM method with sampling from molecular dynamics. Our results reveal that the positive partial charge of glutamine near the C14−C15 bond of retinal shortens the effective conjugation length of the chromophore compared to the leucine residue. The derived mechanism can be applied to explain the color regulation in other retinal proteins and can serve as a guideline for rational design of spectral shifts.
Collapse
Affiliation(s)
- Jonathan R Church
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, 9190401, Israel
| | - Gil S Amoyal
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, 9190401, Israel
| | - Veniamin A Borin
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, 9190401, Israel
| | - Suliman Adam
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, 9190401, Israel
| | | | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, 9190401, Israel
| |
Collapse
|
8
|
Pedraza-González L, Barneschi L, Padula D, De Vico L, Olivucci M. Evolution of the Automatic Rhodopsin Modeling (ARM) Protocol. Top Curr Chem (Cham) 2022; 380:21. [PMID: 35291019 PMCID: PMC8924150 DOI: 10.1007/s41061-022-00374-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Accepted: 01/29/2022] [Indexed: 10/27/2022]
Abstract
In recent years, photoactive proteins such as rhodopsins have become a common target for cutting-edge research in the field of optogenetics. Alongside wet-lab research, computational methods are also developing rapidly to provide the necessary tools to analyze and rationalize experimental results and, most of all, drive the design of novel systems. The Automatic Rhodopsin Modeling (ARM) protocol is focused on providing exactly the necessary computational tools to study rhodopsins, those being either natural or resulting from mutations. The code has evolved along the years to finally provide results that are reproducible by any user, accurate and reliable so as to replicate experimental trends. Furthermore, the code is efficient in terms of necessary computing resources and time, and scalable in terms of both number of concurrent calculations as well as features. In this review, we will show how the code underlying ARM achieved each of these properties.
Collapse
Affiliation(s)
- Laura Pedraza-González
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy. .,Department of Chemistry and Industrial Chemistry, University of Pisa, Via Moruzzi 13, 56124, Pisa, Italy.
| | - Leonardo Barneschi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy
| | - Daniele Padula
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy.
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy. .,Department of Chemistry, Bowling Green State University, Bowling Green, OH, 43403, USA.
| |
Collapse
|
9
|
Church JR, Olsen JMH, Schapiro I. The Impact of Retinal Configuration on the Protein-Chromophore Interactions in Bistable Jumping Spider Rhodopsin-1. Molecules 2021; 27:71. [PMID: 35011302 PMCID: PMC8746357 DOI: 10.3390/molecules27010071] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2021] [Revised: 12/15/2021] [Accepted: 12/17/2021] [Indexed: 12/02/2022] Open
Abstract
Bistable rhodopsins have two stable forms that can be interconverted by light. Due to their ability to act as photoswitches, these proteins are considered as ideal candidates for applications such as optogenetics. In this work, we analyze a recently crystalized bistable rhodopsin, namely the jumping spider rhodopsin-1 (JSR1). This rhodopsin exhibits identical absorption maxima for the parent and the photoproduct form, which impedes its broad application. We performed hybrid QM/MM simulations to study three isomers of the retinal chromophore: the 9-cis, 11-cis and all-trans configurations. The main aim was to gain insight into the specific interactions of each isomer and their impact on the absorption maximum in JSR1. The absorption spectra were computed using sampled snapshots from QM/MM molecular dynamics trajectories and compared to their experimental counterparts. The chromophore-protein interactions were analyzed by visualizing the electrostatic potential of the protein and projecting it onto the chromophore. It was found that the distance between a nearby tyrosine (Y126) residue plays a larger role in the predicted absorption maximum than the primary counterion (E194). Geometric differences between the isomers were also noted, including a structural change in the polyene chain of the chromophore, as well as changes in the nearby hydrogen bonding network.
Collapse
Affiliation(s)
- Jonathan R. Church
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel;
| | | | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel;
| |
Collapse
|
10
|
Pro219 is an electrostatic color determinant in the light-driven sodium pump KR2. Commun Biol 2021; 4:1185. [PMID: 34645937 PMCID: PMC8514524 DOI: 10.1038/s42003-021-02684-z] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2021] [Accepted: 09/19/2021] [Indexed: 11/13/2022] Open
Abstract
Color tuning in animal and microbial rhodopsins has attracted the interest of many researchers, as the color of their common retinal chromophores is modulated by the amino acid residues forming the chromophore cavity. Critical cavity amino acid residues are often called “color switches”, as the rhodopsin color is effectively tuned through their substitution. Well-known color switches are the L/Q and A/TS switches located in the C and G helices of the microbial rhodopsin structure respectively. Recently, we reported on a third G/P switch located in the F helix of the light-driven sodium pumps of KR2 and JsNaR causing substantial spectral red-shifts in the latter with respect to the former. In order to investigate the molecular-level mechanism driving such switching function, here we present an exhaustive mutation, spectroscopic and computational investigation of the P219X mutant set of KR2. To do so, we study the changes in the absorption band of the 19 possible mutants and construct, semi-automatically, the corresponding hybrid quantum mechanics/molecular mechanics models. We found that the P219X feature a red-shifted light absorption with the only exception of P219R. The analysis of the corresponding models indicate that the G/P switch induces red-shifting variations via electrostatic interactions, while replacement-induced chromophore geometrical (steric) distortions play a minor role. However, the same analysis indicates that the P219R blue-shifted variant has a more complex origin involving both electrostatic and steric changes accompanied by protonation state and hydrogen bond networks modifications. These results make it difficult to extract simple rules or formulate theories for predicting how a switch operates without considering the atomistic details and environmental consequences of the side chain replacement. Nakajima, Pedraza-González et al. provide a comprehensive investigation of amino acid mutations at position 219 of the sodium pump rhodopsin, KR2, and their role in the color tuning of the retinal chromophore. They prepared P219X (X= A, C, D, E, F, G, H, I, K, L, M, N, Q, R, S, T, V, W, Y) mutants of KR2, and find that all mutants are red-shifted, except for P219R, highlighting its role as a color determinant in the light-driven pump KR2.
Collapse
|