An Atomistic Model of a Precursor State of Light-Induced Channel Opening of Channelrhodopsin

RISMiCal: A software package to perform fast RISM/3D-RISM calculations

Solvent plays an essential role in a variety of chemical, physical, and biological processes that occur in the solution phase. The reference interaction site model (RISM) and its three-dimensional extension (3D-RISM) serve as powerful computational tools for modeling solvation effects in chemical reactions, biological functions, and structure formations. We present the RISM integrated calculator (RISMiCal) program package, which is based on RISM and 3D-RISM theories with fast GPU code. RISMiCal has been developed as an integrated RISM/3D-RISM program that has interfaces with external programs such as Gaussian16, GAMESS, and Tinker. Fast 3D-RISM programs for single- and multi-GPU codes written in CUDA would enhance the availability of these hybrid methods because they require the performance of many computationally expensive 3D-RISM calculations. We expect that our package can be widely applied for chemical and biological processes in solvent. The RISMiCal package is available at https://rismical-dev.github.io.

Unique Vibrational Characteristics and Structures of the Photoexcited Retinal Chromophore in Ion-Pumping Rhodopsins

Photoisomerization of an all-trans-retinal chromophore triggers ion transport in microbial ion-pumping rhodopsins. Understanding chromophore structures in the electronically excited (S1) state provides insights into the structural evolution on the potential energy surface of the photoexcited state. In this study, we examined the structure of the S1-state chromophore in Natronomonas pharaonis halorhodopsin (NpHR), a chloride ion-pumping rhodopsin, using time-resolved resonance Raman spectroscopy. The spectral patterns of the S1-state chromophore were completely different from those of the ground-state chromophore, resulting from unique vibrational characteristics and the structure of the S1 state. Mode assignments were based on a combination of deuteration shifts of the Raman bands and hybrid quantum mechanics-molecular mechanics calculations. The present observations suggest a weakened bond alternation in the π conjugation system. A strong hydrogen-out-of-plane bending band was observed in the Raman spectra of the S1-state chromophore in NpHR, indicating a twisted polyene structure. Similar frequency shifts for the C═N/C═C and C-C stretching modes of the S1-state chromophore in NpHR were observed in the Raman spectra of sodium ion-pumping and proton-pumping rhodopsins, suggesting that these unique features are common to the S1 states of ion-pumping rhodopsins.

Twisting and Protonation of Retinal Chromophore Regulate Channel Gating of Channelrhodopsin C1C2

Channelrhodopsins (ChRs) are light-gated ion channels and central optogenetic tools that can control neuronal activity with high temporal resolution at the single-cell level. Although their application in optogenetics has rapidly progressed, it is unsolved how their channels open and close. ChRs transport ions through a series of interlocking elementary processes that occur over a broad time scale of subpicoseconds to seconds. During these processes, the retinal chromophore functions as a channel regulatory domain and transfers the optical input as local structural changes to the channel operating domain, the helices, leading to channel gating. Thus, the core question on channel gating dynamics is how the retinal chromophore structure changes throughout the photocycle and what rate-limits the kinetics. Here, we investigated the structural changes in the retinal chromophore of canonical ChR, C1C2, in all photointermediates using time-resolved resonance Raman spectroscopy. Moreover, to reveal the rate-limiting factors of the photocycle and channel gating, we measured the kinetic isotope effect of all photoreaction processes using laser flash photolysis and laser patch clamp, respectively. Spectroscopic and electrophysiological results provided the following understanding of the channel gating: the retinal chromophore highly twists upon the retinal Schiff base (RSB) deprotonation, causing the surrounding helices to move and open the channel. The ion-conducting pathway includes the RSB, where inflowing water mediates the proton to the deprotonated RSB. The twisting of the retinal chromophore relaxes upon the RSB reprotonation, which closes the channel. The RSB reprotonation rate-limits the channel closing.

Properties of a Single Amino Acid Residue in the Third Transmembrane Domain Determine the Kinetics of Ambient Light-Sensitive Channelrhodopsin

Channelrhodopsins have been utilized in gene therapy to restore vision in patients with retinitis pigmentosa and their channel kinetics are an important factor to consider in such applications. We investigated the channel kinetics of ComV1 variants with different amino acid residues at the 172nd position. Patch clamp methods were used to record the photocurrents induced by stimuli from diodes in HEK293 cells transfected with plasmid vectors. The channel kinetics (τon and τoff) were considerably altered by the replacement of the 172nd amino acid and was dependent on the amino acid characteristics. The size of amino acids at this position correlated with τon and decay, whereas the solubility correlated with τon and τoff. Molecular dynamic simulation indicated that the ion tunnel constructed by H172, E121, and R306 widened due to H172A variant, whereas the interaction between A172 and the surrounding amino acids weakened compared with H172. The bottleneck radius of the ion gate constructed with the 172nd amino acid affected the photocurrent and channel kinetics. The 172nd amino acid in ComV1 is a key residue for determining channel kinetics as its properties alter the radius of the ion gate. Our findings can be used to improve the channel kinetics of channelrhodopsins.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Insight into the structural dynamics of light sensitive proteins from time-resolved crystallography and quantum chemical calculations

The structural dynamics underlying molecular mechanisms of light-sensitive proteins can be studied by a variety of experimental and computational biophysical techniques. Here we review recent progress in combining time-resolved crystallography at X-ray free electron lasers and quantum chemical calculations to study structural changes in photoenzymes, photosynthetic proteins, photoreceptors, and photoswitchable fluorescent proteins following photoexcitation.

Hybrid QM/MM Free-Energy Evaluation of Drug-Resistant Mutational Effect on the Binding of an Inhibitor Indinavir to HIV-1 Protease

A human immunodeficiency virus-1 (HIV-1) protease is a homodimeric aspartic protease essential for the replication of HIV. The HIV-1 protease is a target protein in drug discovery for antiretroviral therapy, and various inhibitor molecules of transition state analogues have been developed. However, serious drug-resistant mutants have emerged. For understanding the molecular mechanism of the drug resistance, an accurate examination of the impacts of the mutations on ligand binding and enzymatic activity is necessary. Here, we present a molecular simulation study on the ligand binding of indinavir, a potent transition state analogue inhibitor, to the wild-type protein and a V82T/I84V drug-resistant mutant of the HIV-1 protease. We employed a hybrid ab initio quantum mechanical/molecular mechanical (QM/MM) free-energy optimization technique which combines a highly accurate QM description of the ligand molecule and its interaction with statistically ample conformational sampling of the MM protein environment by long-time molecular dynamics simulations. Through the free-energy calculations of protonation states of catalytic groups at the binding pocket and of the ligand-binding affinity changes upon the mutations, we successfully reproduced the experimentally observed significant reduction of the binding affinity upon the drug-resistant mutations and elucidated the underlying molecular mechanism. The present study opens the way for understanding the molecular mechanism of drug resistance through the direct quantitative comparison of ligand binding and enzymatic reaction with the same accuracy.

Time-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin

Channelrhodopsins (ChRs) are microbial light-gated ion channels utilized in optogenetics to control neural activity with light . Light absorption causes retinal chromophore isomerization and subsequent protein conformational changes visualized as optically distinguished intermediates, coupled with channel opening and closing. However, the detailed molecular events underlying channel gating remain unknown. We performed time-resolved serial femtosecond crystallographic analyses of ChR by using an X-ray free electron laser, which revealed conformational changes following photoactivation. The isomerized retinal adopts a twisted conformation and shifts toward the putative internal proton donor residues, consequently inducing an outward shift of TM3, as well as a local deformation in TM7. These early conformational changes in the pore-forming helices should be the triggers that lead to opening of the ion conducting pore.

Ab Initio Evaluation of the Redox Potential of Cytochrome c

Various biochemical activities of metabolism and biosynthesis are fulfilled by redox processes with explicit electron exchange, which furnish redox enzymes with high chemical reactivity. However, theoretical investigation of a redox process, which simultaneously involves a complex electronic change at a redox metal center and conformational reorganization of the surrounding protein environment coupled to the electronic change, requires computationally conflicting approaches, highly accurate quantum chemical calculations, and long-time molecular dynamics (MD) simulations, limiting the physicochemical understanding of biological redox processes. Here, we theoretically examined a redox process of cytochrome c by means of a hybrid molecular simulation technique, which enables one to consistently treat the redox center at the ab initio quantum chemistry level of theory and the protein reorganization with long-time MD simulations on the microsecond timescale. The calculations successfully evaluated a large absolute redox potential, 4.34 eV, with errors of only 0.03 to 0.34 eV to the experimental ones without any problem-specific empirical parameters. Through the long-time MD sampling, large and nonlinear reorganization of the protein environment was unveiled and the molecular determinants for the redox potential were identified. The present ab initio approach significantly expands the applicability of theoretical investigation to biological redox systems with more electronically complicated redox centers such as polynuclear transition metal complexes.

Structure-Function Relationship of Channelrhodopsins

Ion-translocating rhodopsins, especially channelrhodopsins (ChRs), have attracted broad attention as a powerful tool to modulate the membrane potential of cells with light (optogenetics). Because of recent biophysical, spectroscopic, and computational studies, including the structural determination of cation and anion ChRs, our understanding of the molecular mechanism underlying light-gated ion conduction has been greatly advanced. In this chapter, I first describe the background of rhodopsin family proteins including ChR, and how the optogenetics technology has been established from the discovery of first ChR in 2002. I later introduce the recent findings of the structure-function relationship of ChR by comparing the crystal structures of cation and anion ChRs. I further discuss the future goal in the fields of ChR research and optogenetic tool development.

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.

The effect on ion channel of different protonation states of E90 in channelrhodopsin-2: a molecular dynamics simulation

Channelrhodopsin-2 (ChR2) is a cationic channel protein that has been extensively studied in optogenetics. The ion channel is opened via a series of proton transfers and H-bond changes during the photocycle but the detailed mechanism is still unknown. Molecular dynamics (MD) simulations with enhanced sampling were performed on the dark-adapted state (i.e., D470) and two photocycle intermediates (P₁⁵⁰⁰ and P₂³⁹⁰) to study the proton transfer path of the Schiff base and the subsequent conformational changes. The results suggest there are two possible proton transfer pathways from the Schiff base to proton acceptors (i.e., E123 or D253), depending on the protonation of E90. If E90 is protonated in the P₁⁵⁰⁰ state, the proton on the Schiff base will transfer to E123. The polyene chain of 13-cis retinal tilts and opens the channel that detours the blocking central gate (CG) and forms a narrow channel through the transmembrane helices (TM) 2, 3, 6 and 7. In contrast, if E90 deprotonates after retinal isomerization, the primary proton acceptor is D253, and an almost-open channel through TM1, 2, 3 and 7 is generated. The channel diameter is very close to the experimental value. The potential mean force (PMF) suggests that the free energy is extremely low for ions passing through this channel.

With E90 protonated, the proton acceptor of RSBH⁺ is E123 with a narrow channel along TM3; while with E90 deprotonated, proton transfer from RSBH⁺ to D253 generates an approximately open channel along TM2.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.