Photoactivation Intermediates of a G-Protein Coupled Receptor Rhodopsin Investigated by a Hybrid Molecular Simulation

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.

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.

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.

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

Channelrhodopsins (ChRs) are microbial light-gated ion channels with a retinal chromophore and are widely utilized in optogenetics to precisely control neuronal activity with light. Despite increasing understanding of their structures and photoactivation kinetics, the atomistic mechanism of light gating and ion conduction remains elusive. Here, we present an atomic structural model of a chimeric ChR in a precursor state of the channel opening determined by an accurate hybrid molecular simulation technique and a statistical theory of internal water distribution. The photoactivated structure features extensive tilt of the chromophore accompanied by redistribution of water molecules in its binding pocket, which is absent in previously known photoactivated structures of analogous photoreceptors, and widely agrees with structural and spectroscopic experimental evidence of ChRs. The atomistic model manifests a photoactivated ion-conduction pathway that is markedly different from a previously proposed one and successfully explains experimentally observed mutagenic effects on key channel properties.