Molecular Mechanism of Wide Photoabsorption Spectral Shifts of Color Variants of Human Cellular Retinol Binding Protein II

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.

Control of Protonated Schiff Base Excited State Decay within Visual Protein Mimics: A Unified Model for Retinal Chromophores

Artificial biomimetic chromophore-protein complexes inspired by natural visual pigments can feature color tunability across the full visible spectrum. However, control of excited state dynamics of the retinal chromophore, which is of paramount importance for technological applications, is lacking due to its complex and subtle photophysics/photochemistry. Here, ultrafast transient absorption spectroscopy and quantum mechanics/molecular mechanics simulations are combined for the study of highly tunable rhodopsin mimics, as compared to retinal chromophores in solution. Conical intersections and transient fluorescent intermediates are identified with atomistic resolution, providing unambiguous assignment of their ultrafast excited state absorption features. The results point out that the electrostatic environment of the chromophore, modified by protein point mutations, affects its excited state properties allowing control of its photophysics with same power of chemical modifications of the chromophore. The complex nature of such fine control is a fundamental knowledge for the design of bio-mimetic opto-electronic and photonic devices.

Accurate Prediction of Absorption Spectral Shifts of Proteorhodopsin Using a Fragment-Based Quantum Mechanical Method

Many experiments have been carried out to display different colors of Proteorhodopsin (PR) and its mutants, but the mechanism of color tuning of PR was not fully elucidated. In this study, we applied the Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps (EE-GMFCC) method to the prediction of excitation energies of PRs. Excitation energies of 10 variants of Blue Proteorhodopsin (BPR-PR105Q) in residue 105GLN were calculated with the EE-GMFCC method at the TD-B3LYP/6-31G* level. The calculated results show good correlation with the experimental values of absorption wavelengths, although the experimental wavelength range among these systems is less than 50 nm. The ensemble-averaged electric fields along the polyene chain of retinal correlated well with EE-GMFCC calculated excitation energies for these 10 PRs, suggesting that electrostatic interactions from nearby residues are responsible for the color tuning. We also utilized the GMFCC method to decompose the excitation energy contribution per residue surrounding the chromophore. Our results show that residues ASP97 and ASP227 have the largest contribution to the absorption spectral shift of PR among the nearby residues of retinal. This work demonstrates that the EE-GMFCC method can be applied to accurately predict the absorption spectral shifts for biomacromolecules.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

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.

NeoR, a near-infrared absorbing rhodopsin

The Rhizoclosmatium globosum genome encodes three rhodopsin-guanylyl cyclases (RGCs), which are predicted to facilitate visual orientation of the fungal zoospores. Here, we show that RGC1 and RGC2 function as light-activated cyclases only upon heterodimerization with RGC3 (NeoR). RGC1/2 utilize conventional green or blue-light-sensitive rhodopsins (λmax = 550 and 480 nm, respectively), with short-lived signaling states, responsible for light-activation of the enzyme. The bistable NeoR is photoswitchable between a near-infrared-sensitive (NIR, λmax = 690 nm) highly fluorescent state (QF = 0.2) and a UV-sensitive non-fluorescent state, thereby modulating the activity by NIR pre-illumination. No other rhodopsin has been reported so far to be functional as a heterooligomer, or as having such a long wavelength absorption or high fluorescence yield. Site-specific mutagenesis and hybrid quantum mechanics/molecular mechanics simulations support the idea that the unusual photochemical properties result from the rigidity of the retinal chromophore and a unique counterion triad composed of two glutamic and one aspartic acids. These findings substantially expand our understanding of the natural potential and limitations of spectral tuning in rhodopsin photoreceptors.

The Effective Conjugation Length Is Responsible for the Red/Green Spectral Tuning in the Cyanobacteriochrome Slr1393g3

The origin of the spectral shift from a red- to a green-absorbing form in a cyanobacteriochrome, Slr1393g3, was identified by combined quantum mechanics/molecular mechanics simulations. This protein, related to classical phytochromes, carries the open-chain tetrapyrrole chromophore phycocyanobilin. Our calculations reveal that the effective conjugation length in the chromophore becomes shorter upon conversion from the red to the green form. This is related to the planarity of the entire chromophore. A large distortion was found for the terminal pyrrole rings A and D; however, the D ring contributes more strongly to the photoproduct tuning, despite a larger change in the twist of the A ring. Our findings implicate that the D ring twist can be exploited to regulate the absorption of the photoproduct. Hence, mutations that affect the D ring twist can lead to rational tuning of the photoproduct absorption, allowing the tailoring of cyanobacteriochromes for biotechnological applications such as optogenetics and bioimaging.

The Dynamic Origin of Color Tuning in Proteins Revealed by a Carotenoid Pigment

Understanding the microscopic origin of the color tuning in pigment-protein complexes is a challenging yet fundamental issue in photoactive biological systems. Here, we propose a possible interpretation by using a state-of-the-art multiscale strategy based on the integration of quantum chemistry and polarizable atomistic embeddings into a dynamic description. By means of such a strategy we are able to resolve the long-standing dispute over the coloration mechanism in the crustacyanin protein. It is shown that the combination of the dynamical flexibility of the carotenoid pigments (astaxanthin) with the responsive protein environment is essential to obtain quantitative predictions of the spectral tuning. The strong linear correlation between the excitation energies and the bond length alternation in the long-chain carotenoids modulated by the dynamical protein environment is a novel finding explaining the high color tunability in crustacyanin.

Theoretical approaches for dynamical ordering of biomolecular systems

BACKGROUND

Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies.

SCOPE OF REVIEW

We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified.

MAJOR CONCLUSIONS

The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function.

GENERAL SIGNIFICANCE

These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

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

Rhodopsin is a G-protein coupled receptor functioning as a photoreceptor for vision through photoactivation of a covalently bound ligand of a retinal protonated Schiff base chromophore. Despite the availability of structural information on the inactivated and activated forms of the receptor, the transition processes initiated by the photoabsorption have not been well understood. Here we theoretically examined the photoactivation processes by means of molecular dynamics (MD) simulations and ab initio quantum mechanical/molecular mechanical (QM/MM) free energy geometry optimizations which enabled accurate geometry determination of the ligand molecule in ample statistical conformational samples of the protein. Structures of the intermediate states of the activation process, blue-shifted intermediate and Lumi, as well as the dark state first generated by MD simulations and then refined by the QM/MM free energy geometry optimizations were characterized by large displacement of the β-ionone ring of retinal along with change in the hydrogen bond of the protonated Schiff base. The ab initio calculations of vibrational and electronic spectroscopic properties of those states well reproduced the experimental observations and successfully identified the molecular origins underlying the spectroscopic features. The structural evolution in the formation of the intermediates provides a molecular insight into the efficient activation processes of the receptor.

Electrostatic versus Resonance Interactions in Photoreceptor Proteins: The Case of Rhodopsin

Light sensing in photoreceptor proteins is subtly modulated by the multiple interactions between the chromophoric unit and its binding pocket. Many theoretical and experimental studies have tried to uncover the fundamental origin of these interactions but reached contradictory conclusions as to whether electrostatics, polarization, or intrinsically quantum effects prevail. Here, we select rhodopsin as a prototypical photoreceptor system to reveal the molecular mechanism underlying these interactions and regulating the spectral tuning. Combining a multireference perturbation method and density functional theory with a classical but atomistic and polarizable embedding scheme, we show that accounting for electrostatics only leads to a qualitatively wrong picture, while a responsive environment can successfully capture both the classical and quantum dominant effects. Several residues are found to tune the excitation by both differentially stabilizing ground and excited states and through nonclassical "inductive resonance" interactions. The results obtained with such a quantum-in-classical model are validated against both experimental data and fully quantum calculations.

Space and Time Evolution of the Electrostatic Potential During the Activation of a Visual Pigment

Animal and microbial retinal proteins employ the Schiff base of retinal as their chromophore. Here, the possible consequences of the charge translocation associated with the light-induced dynamics of the chromophore of a visual opsin are investigated along a representative semiclassical trajectory. We show that the evolution of the electrostatic potential projected by the chromophore onto the surrounding protein displays intense but topographically localized sudden variations in proximity of the decay region. pKa calculations carried out on selected snapshots used as probes, indicate that the only residue which may be sensitive to the electrostatic potential shift is Glu181. Accordingly, our results suggest that the frail Tyr191/268-Glu181-Wat2-Ser186 hydrogen bond network may be perturbed by the transient variations of the electrostatic potential.