Ground-state properties of the retinal molecule: from quantum mechanical to classical mechanical computations of retinal proteins

Thermal isomerization rates in retinal analogues using Ab-Initio molecular dynamics

For a detailed understanding of chemical processes in nature and industry, we need accurate models of chemical reactions in complex environments. While Eyring transition state theory is commonly used for modeling chemical reactions, it is most accurate for small molecules in the gas phase. A wide range of alternative rate theories exist that can better capture reactions involving complex molecules and environmental effects. However, they require that the chemical reaction is sampled by molecular dynamics simulations. This is a formidable challenge since the accessible simulation timescales are many orders of magnitude smaller than typical timescales of chemical reactions. To overcome these limitations, rare event methods involving enhanced molecular dynamics sampling are employed. In this work, thermal isomerization of retinal is studied using tight-binding density functional theory. Results from transition state theory are compared to those obtained from enhanced sampling. Rates obtained from dynamical reweighting using infrequent metadynamics simulations were in close agreement with those from transition state theory. Meanwhile, rates obtained from application of Kramers' rate equation to a sampled free energy profile along a torsional dihedral reaction coordinate were found to be up to three orders of magnitude higher. This discrepancy raises concerns about applying rate methods to one-dimensional reaction coordinates in chemical reactions.

Existence of two substates in the O intermediate of the bacteriorhodopsin photocycle

The proton pumping cycle of bacteriorhodopsin (bR) is initiated when the retinal chromophore with the 13-trans configuration is photo-isomerized into the 13-cis configuration. To understand the recovery processes of the initial retinal configuration that occur in the late stage of the photocycle, we have performed a comprehensive analysis of absorption kinetics data collected at various pH levels and at different salt concentrations. The result of analysis revealed the following features of the late stages of the trans photocycle. i) Two substates occur in the O intermediate. ii) The visible absorption band of the first substate (O1) appears at a much shorter wavelength than that of the late substate (O2). iii) O1 is in rapid equilibrium with the preceding state (N), but O1 becomes less stable than N when an ionizable residue (X₁) with a pKa value of 6.5 (in 2 M KCl) is deprotonated. iv) At a low pH and at a low salt concentration, the decay time constant of O2 is longer than those of the preceding states, but the relationship between these time constants is altered when the medium pH or the salt concentration is increased. On the basis of the present observations and previous studies on the structure of the chromophore in O, we suspect that the retinal chromophore in O1 takes on a distorted 13-cis configuration and the O1-to-O2 transition is accompanied by cis-to-trans isomerization about C13C14 bond.

Two substates in the O intermediate of the light-driven proton pump archaerhodopsin-2

The proton pumping cycle of archaerhodopsin-2 (aR2) was investigated over a wide pH range and at different salt concentrations. We have found that two substates, which are spectroscopically and kinetically distinguishable, occur in the O intermediate. The first O-intermediate (O1) absorbs maximumly at ~580 nm, whereas the late O-intermediate (O2) absorbs maximumly at 605 nm. At neutral pH, O1 is in rapid equilibrium with the N intermediate. When the medium pH is increased, O1 becomes less stable than N and, in proportion to the amount of O1 in the dynamic equilibrium between N and O1, the formation rate of O2 decreases. By contrast, the decay rate of O2 increases ~100 folds when the pH of a low-salt membrane suspension is increased from 5.5 to 7.5 or when the salt concentration is increased to 2 M KCl. Together with our recent study on two substates in the O intermediate of bacteriorhodopsin (bR), the present study suggests that the thermally activated re-isomerization of the retinylidene chromophore into the initial all-trans configuration takes place in the O1-to-O2 transition; that is, O1 contains a distorted 13-cis chromophore. It is also found that the pKa value of the key ionizable residue (Asp101^aR2, Asp96^bR) in the proton uptake channel is elevated in the O1 state of aR2 as compared to the O1 state of bR. This implies that the structural property of O1 in the aR2 photocycle can be investigated over a wide pH range.

Conserved hydrogen-bond motifs of membrane transporters and receptors

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

Catalysis of Ground State cis[Formula: see text] trans Isomerization of Bacteriorhodopsin's Retinal Chromophore by a Hydrogen-Bond Network

For the photocycle of the membrane protein bacteriorhodopsin to proceed efficiently, the thermal 13-cis to all-trans back-isomerization of the retinal chromophore must return the protein to its resting state on a time-scale of milliseconds. Here, we report on quantum mechanical/molecular mechanical energy calculations examining the structural and energetic determinants of the retinal cis-trans isomerization in the protein environment. The results suggest that a hydrogen-bonded network consisting of the retinal Schiff base, active site amino acid residues, and water molecules can stabilize the twisted retinal, thus reducing the intrinsic energy cost of the cis-trans thermal isomerization barrier.

Chlorophyll-Derivative Modulation of Rhodopsin Signaling Properties through Evolutionarily Conserved Interaction Pathways

Retinal is the light-absorbing chromophore that is responsible for the activation of visual pigments and light-driven ion pumps. Evolutionary changes in the intermolecular interactions of the retinal with specific amino acids allow for adaptation of the spectral characteristics, referred to as spectral tuning. However, it has been proposed that a specific species of dragon fish has bypassed the adaptive evolutionary process of spectral tuning and replaced it with a single evolutionary event: photosensitization of rhodopsin by chlorophyll derivatives. Here, by using a combination of experimental measurements and computational modeling to probe retinal-receptor interactions in rhodopsin, we show how the binding of the chlorophyll derivative, chlorin-e6 (Ce6) in the intracellular domain (ICD) of the receptor allosterically excites G-protein coupled receptor class A (GPCR-A) conserved long-range correlated fluctuations that connect distant parts of the receptor. These long-range correlated motions are associated with regulating the dynamics and intermolecular interactions of specific amino acids in the retinal ligand-binding pocket that have been associated with shifts in the absorbance peak maximum (λ_max) and hence, spectral sensitivity of the visual system. Moreover, the binding of Ce6 affects the overall global properties of the receptor. Specifically, we find that Ce6-induced dynamics alter the thermal stability of rhodopsin by adjusting hydrogen-bonding interactions near the receptor active-site that consequently also influences the intrinsic conformational equilibrium of the receptor. Due to the conservation of the ICD residues amongst different receptors in this class and the fact that all GPCR-A receptors share a common mechanism of activation, it is possible that the allosteric associations excited in rhodopsin with Ce6 binding are a common feature in all class A GPCRs.

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.

Insights into the formation of carlactone from in-depth analysis of the CCD8-catalyzed reactions

Strigolactones are a new class of phytohormones synthesized from carotenoids via carlactone. The complex structure of carlactone is not easily deducible from its precursor, a cis-configured β-carotene cleavage product, and is thus formed via a poorly understood series of reactions and molecular rearrangements, all catalyzed by only one enzyme, the carotenoid cleavage dioxygenase 8 (CCD8). Moreover, the reactions leading to carlactone are expected to form a second, yet unidentified product. In this study, we used 13 C and 18 O-labeling to shed light on the reactions catalyzed by CCD8. The characterization of the resulting carlactone by LC-MS and NMR, and the identification of the assumed, less accessible second product allowed us to formulate a minimal reaction mechanism for carlactone generation.