Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

Molecular mechanisms and evolutionary robustness of a color switch in proteorhodopsins

Proteorhodopsins are widely distributed photoreceptors from marine bacteria. Their discovery revealed a high degree of evolutionary adaptation to ambient light, resulting in blue- and green-absorbing variants that correlate with a conserved glutamine/leucine at position 105. On the basis of an integrated approach combining sensitivity-enhanced solid-state nuclear magnetic resonance (ssNMR) spectroscopy and linear-scaling quantum mechanics/molecular mechanics (QM/MM) methods, this single residue is shown to be responsible for a variety of synergistically coupled structural and electrostatic changes along the retinal polyene chain, ionone ring, and within the binding pocket. They collectively explain the observed color shift. Furthermore, analysis of the differences in chemical shift between nuclei within the same residues in green and blue proteorhodopsins also reveals a correlation with the respective degree of conservation. Our data show that the highly conserved color change mainly affects other highly conserved residues, illustrating a high degree of robustness of the color phenotype to sequence variation.

Quantum Mechanics/Molecular mechanics calculations predict A1, not A2, is present in melanopsin (Opn4m) of red-eared slider turtles (Trachemys scripta elegans)

Melanopsin is a photopigment that plays a role in non-visual, light-driven, cellular processes such as modulation of circadian rhythms, retinal vascular development, and the pupillary light reflex (PLR). In this study, computational methods were used to understand which chromophore is harbored by melanopsin in red-eared slider turtles (Trachemys scripta elegans). In mammals, the vitamin A derivative 11-cis-retinal (A1) is the chromophore, which provides functionality for melanopsin. However, in red-eared slider turtles, a member of the reptilian class, the identity of the chromophore remains unclear. Red-eared slider turtles, similar to other freshwater vertebrates, possess visual pigments that harbor a different vitamin A derivative, 11-cis-3,4-didehydroretinal (A2), making their pigments more sensitive to red-light than blue-light, therefore, suggesting the chromophore to be the A2 derivative instead of the A1. To help resolve the chromophore identity, in this work, computational homology models of melanopsin in red-eared slider turtles were first constructed. Next, quantum mechanics/molecular mechanics (QM/MM) calculations were carried out to compare how A1 and A2 derivatives bind to melanopsin. Time dependent density functional theory (TDDFT) calculations were then used to determine the excitation energy of the pigments. Lastly, calculated excitation energies were compared to experimental spectral sensitivity data from responses by the irises of red-eared sliders. Contrary to what was expected, our results suggest that melanopsin in red-eared slider turtles is more likely to harbor the A1 chromophore than the A2. Furthermore, a glutamine (Q622.56) and tyrosine (Y853.28) residue in the chromophore binding pocket are shown to play a role in the spectral tuning of the chromophore.

Molecular Mechanism of Spectral Tuning by Chloride Binding in Monkey Green Sensitive Visual Pigment

The visual pigments of the cones perceive red, green, and blue colors. The monkey green (MG) pigment possesses a unique Cl^- binding site; however, its relationship to the spectral tuning in green pigments remains elusive. Recently, FTIR spectroscopy revealed the characteristic structural modifications of the retinal binding site by Cl^- binding. Herein, we report the computational structural modeling of MG pigments and quantum-chemical simulation to investigate its spectral redshift and physicochemical relevance when Cl^- is present. Our protein structures reflect the previously suggested structural changes. AlphaFold2 failed to predict these structural changes. Excited-state calculations successfully reproduced the experimental red-shifted absorption energies, corroborating our protein structures. Electrostatic energy decomposition revealed that the redshift results from the His197 protonation state and conformations of Glu129, Ser202, and Ala308; however, Cl^- itself contributes to the blueshift. Site-directed mutagenesis supported our analysis. These modeled structures may provide a valuable foundation for studying cone pigments.

Rhodopsins: An Excitingly Versatile Protein Species for Research, Development and Creative Engineering

The first member and eponym of the rhodopsin family was identified in the 1930s as the visual pigment of the rod photoreceptor cell in the animal retina. It was found to be a membrane protein, owing its photosensitivity to the presence of a covalently bound chromophoric group. This group, derived from vitamin A, was appropriately dubbed retinal. In the 1970s a microbial counterpart of this species was discovered in an archaeon, being a membrane protein also harbouring retinal as a chromophore, and named bacteriorhodopsin. Since their discovery a photogenic panorama unfolded, where up to date new members and subspecies with a variety of light-driven functionality have been added to this family. The animal branch, meanwhile categorized as type-2 rhodopsins, turned out to form a large subclass in the superfamily of G protein-coupled receptors and are essential to multiple elements of light-dependent animal sensory physiology. The microbial branch, the type-1 rhodopsins, largely function as light-driven ion pumps or channels, but also contain sensory-active and enzyme-sustaining subspecies. In this review we will follow the development of this exciting membrane protein panorama in a representative number of highlights and will present a prospect of their extraordinary future potential.

Evolution of the Automatic Rhodopsin Modeling (ARM) Protocol

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.

Efficiently Computing Excitations of Complex Systems: Linear-Scaling Time-Dependent Embedded Mean-Field Theory in Implicit Solvent

Quantum embedding schemes have the potential to significantly reduce the computational cost of first-principles calculations while maintaining accuracy, particularly for calculations of electronic excitations in complex systems. In this work, I combine time-dependent embedded mean field theory (TD-EMFT) with linear-scaling density functional theory and implicit solvation models, extending previous work within the ONETEP code. This provides a way to perform multilevel calculations of electronic excitations on very large systems, where long-range environmental effects, both quantum and classical in nature, are important. I demonstrate the power of this method by performing simulations on a variety of systems, including a molecular dimer, a chromophore in solution, and a doped molecular crystal. This work paves the way for high accuracy calculations to be performed on large-scale systems that were previously beyond the reach of quantum embedding schemes.

Protein Response Effects on Cofactor Excitation Energies from First Principles: Augmenting Subsystem Time-Dependent Density-Functional Theory with Many-Body Expansion Techniques

We investigate the possibility of describing protein response effects on a chromophore excitation by means of subsystem time-dependent density-functional theory (sTDDFT) in combination with a many-body expansion (MBE) approach. While sTDDFT is in principle intrinsically able to include such contributions, addressing cofactor excitations in protein models or entire proteins with full environment-response treatments is currently out of reach. Taking different model structures of the green fluorescent protein (GFP) and bovine rhodopsin as examples, we demonstrate that an embedded-MBE approach based on sTDDFT in its simplest version leads to a good agreement of the predicted protein response effect already at second order. To reproduce reference response effects from nonsubsystem TDDFT calculations quantitatively (error ≤ 5%), however, a third- or even fourth-order MBE may be required. For the latter case, we explore a selective inclusion of fourth-order terms that drastically reduces the computational burden. In addition, we demonstrate how this sTDDFT-MBE treatment can be utilized as an analysis tool to identify residues with dominant response contributions. This, in turn, can be employed to arrive at smaller structural models for light-absorbing proteins, which still feature all of the main characteristics in terms of photoresponse properties.

Electronic Couplings and Electrostatic Interactions Behind the Light Absorption of Retinal Proteins

The photo-functional chromophore retinal exhibits a wide variety of optical absorption properties depending on its intermolecular interactions with surrounding proteins and other chromophores. By utilizing these properties, microbial and animal rhodopsins express biological functions such as ion-transport and signal transduction. In this review, we present the molecular mechanisms underlying light absorption in rhodopsins, as revealed by quantum chemical calculations. Here, symmetry-adapted cluster-configuration interaction (SAC-CI), combined quantum mechanical and molecular mechanical (QM/MM), and transition-density-fragment interaction (TDFI) methods are used to describe the electronic structure of the retinal, the surrounding protein environment, and the electronic coupling between chromophores, respectively. These computational approaches provide successful reproductions of experimentally observed absorption and circular dichroism (CD) spectra, as well as insights into the mechanisms of unique optical properties in terms of chromophore-protein electrostatic interactions and chromophore-chromophore electronic couplings. On the basis of the molecular mechanisms revealed in these studies, we also discuss strategies for artificial design of the optical absorption properties of rhodopsins.

Isatin-based virtual high throughput screening, molecular docking, DFT, QM/MM, MD and MM-PBSA study of novel inhibitors of SARS-CoV-2 main protease

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a rapidly growing health care emergency across the world. One of the viral proteases called main protease or Mpro, plays a crucial role in the replication of SARS-CoV-2. As the structure of Mpro of SARS-CoV-2 is similar to the Mpro of SARS-CoV-1 (responsible for SARS outbreak between 2002 and 2004), we hypothesize that the inhibitors of SARS-CoV-1 Mpro can also inhibit the Mpro of SARS-CoV-2. To test this hypothesis, a total of 79 isatin derivatives, which inhibited Mpro activity under in vitro conditions, were selected from the literature, and then screened through AutoDock Vina. The chemical features of the top 5 isatin derivatives with low binding energies (-8.5 to -8.2 kcal/mol) were used to screen similar types of compounds from several small-molecule libraries containing 15856508 compounds. A total of 1,609 compounds with similarity score ≥ 6 were screened and then subjected to docking as well as ADME analysis. Among the compounds screened, 4 ligands form Zinc drug-like library (ZINC000008848565, ZINC000009513563, ZINC000036759789 and ZINC000046053855) showed good ADMET properties, low binding energy (-8.4 to -8.6 kcal/mol), low interaction energy (-72.62 to -50.01 kcal/mol) and high structural stability with Mpro. Hence, the selected ligands might serve as the lead candidates for further wet laboratory validation, optimization and development.Communicated by Ramaswamy H. Sarma.

Vitamin A1/A2 chromophore exchange: Its role in spectral tuning and visual plasticity

Vertebrate rod and cone photoreceptors detect light via a specialized organelle called the outer segment. This structure is packed with light-sensitive molecules known as visual pigments that consist of a G-protein-coupled, seven-transmembrane protein known as opsin, and a chromophore prosthetic group, either 11-cis retinal ('A₁') or 11-cis 3,4-didehydroretinal ('A₂'). The enzyme cyp27c1 converts A₁ into A₂ in the retinal pigment epithelium. Replacing A₁ with A₂ in a visual pigment red-shifts its spectral sensitivity and broadens its bandwidth of absorption at the expense of decreased photosensitivity and increased thermal noise. The use of vitamin A₂-based visual pigments is strongly associated with the occupation of aquatic habitats in which the ambient light is red-shifted. By modulating the A₁/A₂ ratio in the retina, an organism can dynamically tune the spectral sensitivity of the visual system to better match the predominant wavelengths of light in its environment. As many as a quarter of all vertebrate species utilize A₂, at least during a part of their life cycle or under certain environmental conditions. A₂ utilization therefore represents an important and widespread mechanism of sensory plasticity. This review provides an up-to-date account of the A₁/A₂ chromophore exchange system.

An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima

Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λ_max) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λ_max of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λ_max from the experimental λ_max. Although the major effect of water molecules on λ_max originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λ_max mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.

Multiscale Molecular Modeling in G Protein-Coupled Receptor (GPCR)-Ligand Studies

G protein-coupled receptors (GPCRs) are major drug targets due to their ability to facilitate signal transduction across cell membranes, a process that is vital for many physiological functions to occur. The development of computational technology provides modern tools that permit accurate studies of the structures and properties of large chemical systems, such as enzymes and GPCRs, at the molecular level. The advent of multiscale molecular modeling permits the implementation of multiple levels of theories on a system of interest, for instance, assigning chemically relevant regions to high quantum mechanics (QM) level of theory while treating the rest of the system using classical force field (molecular mechanics (MM) potential). Multiscale QM/MM molecular modeling have far-reaching applications in the rational design of GPCR drugs/ligands by affording precise ligand binding configurations through the consideration of conformational plasticity. This enables the identification of key binding site residues that could be targeted to manipulate GPCR function. This review will focus on recent applications of multiscale QM/MM molecular simulations in GPCR studies that could boost the efficiency of future structure-based drug design (SBDD) strategies.

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.

Adaptive Landscapes in the Age of Synthetic Biology

For nearly a century adaptive landscapes have provided overviews of the evolutionary process and yet they remain metaphors. We redefine adaptive landscapes in terms of biological processes rather than descriptive phenomenology. We focus on the underlying mechanisms that generate emergent properties such as epistasis, dominance, trade-offs and adaptive peaks. We illustrate the utility of landscapes in predicting the course of adaptation and the distribution of fitness effects. We abandon aged arguments concerning landscape ruggedness in favor of empirically determining landscape architecture. In so doing, we transform the landscape metaphor into a scientific framework within which causal hypotheses can be tested.

CpHMD-Then-QM/MM Identification of the Amino Acids Responsible for the Anabaena Sensory Rhodopsin pH-Dependent Electronic Absorption Spectrum

Anabaena Sensory Rhodopsin (ASR), a microbial photoactive protein featuring the retinal chromophore in two different conformations, exhibits a pH-dependent electronic absorption spectrum. Using the recently developed CpHMD-then-QM/MM multiscale protocol applied to ASR embedded in a membrane model, the pH-induced changes in its maximum absorption wavelength have been reproduced and analyzed. While the acidic tiny red-shift is essentially correlated with the deprotonation of an aspartic acid located on the ASR extracellular side, the larger blue-shift experimentally reported at pH values larger than 5 involves a cluster of titrating residues sitting on the cytoplasmic side. The ASR pH-dependent spectrum is the consequence of the competitive stabilization of retinal ground and excited states by the protein electrostatic potential.

a-ARM: Automatic Rhodopsin Modeling with Chromophore Cavity Generation, Ionization State Selection, and External Counterion Placement

The Automatic Rhodopsin Modeling (ARM) protocol has recently been proposed as a tool for the fast and parallel generation of basic hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild type and mutant rhodopsins. However, in its present version, input preparation requires a few hours long user's manipulation of the template protein structure, which also impairs the reproducibility of the generated models. This limitation, which makes model building semiautomatic rather than fully automatic, comprises four tasks: definition of the retinal chromophore cavity, assignment of protonation states of the ionizable residues, neutralization of the protein with external counterions, and finally congruous generation of single or multiple mutations. In this work, we show that the automation of the original ARM protocol can be extended to a level suitable for performing the above tasks without user's manipulation and with an input preparation time of minutes. The new protocol, called a-ARM, delivers fully reproducible (i.e., user independent) rhodopsin QM/MM models as well as an improved model quality. More specifically, we show that the trend in vertical excitation energies observed for a set of 25 wild type and 14 mutant rhodopsins is predicted by the new protocol better than when using the original. Such an agreement is reflected by an estimated (relative to the probed set) trend deviation of 0.7 ± 0.5 kcal mol-1 (0.03 ± 0.02 eV) and mean absolute error of 1.0 kcal mol-1 (0.04 eV).

Assessment of MC-PDFT Excitation Energies for a Set of QM/MM Models of Rhodopsins

A methodology for the automatic production of quantum mechanical/molecular mechanical (QM/MM) models of retinal-binding rhodopsin proteins and subsequent prediction of their spectroscopic properties has been proposed recently by some of the authors. The technology employed for the evaluation of the excitation energies is called Automatic Rhodopsin Modeling (ARM), and it involves the use of the complete active space self-consistent field (CASSCF) method followed by a multiconfiguration second-order perturbation theory (in particular, CASPT2) calculation of external correlation energies. Although it was shown that ARM is capable of successfully reproducing and predicting spectroscopic property trends in chromophore-embedding protein sets, practical applications of such technology are limited by the high computational costs of the multiconfiguration perturbation theory calculations. In the present work we benchmark the more affordable multiconfiguration pair-density functional theory (MC-PDFT) method whose accuracy has been recently validated for retinal chromophores in the gas phase, indicating that MC-PDFT could potentially be used to analyze large (e.g., few hundreds) sets of rhodopsin proteins. Here, we test this theory for a set of rhodopsin QM/MM models whose experimental absorption maxima (λ a max) have been measured. The results indicate that MC-PDFT may be employed to calculate λ a max values for this important class of photoresponsive proteins.

Protonation of the Biliverdin IXα Chromophore in the Red and Far-Red Photoactive States of a Bacteriophytochrome

The tetrapyrrole chromophore biliverdin IXα (BV) in the bacteriophytochrome from Deinococcus radiodurans (DrBphP) is usually assumed to be fully protonated, but this assumption has not been systematically validated by experiments or extensive computations. Here, we use force field molecular dynamics simulations and quantum mechanics/molecular mechanics calculations with density functional theory and XMCQDPT2 methods to investigate the effect of the five most probable protonation forms of BV on structural stability, binding pocket interactions, and absorption spectra in the two photochromic states of DrBphP. While agreement with X-ray structural data and measured UV/vis spectra suggest that in both states the protonated form of the chromophore dominates, we also find that a minor population with a deprotonated D-ring could contribute to the red-shifted tail in the absorption spectra.

Predicting solvatochromic shifts and colours of a solvated organic dye: The example of nile red

The solvatochromic shift, as well as the change in colour of the simple organic dye nile red, is studied in two polar and two non-polar solvents in the context of large-scale time-dependent density-functional theory (TDDFT) calculations treating large parts of the solvent environment from first principles. We show that an explicit solvent representation is vital to resolve absorption peak shifts between nile red in n-hexane and toluene, as well as acetone and ethanol. The origin of the failure of implicit solvent models for these solvents is identified as being due to the strong solute-solvent interactions in form of π-stacking and hydrogen bonding in the case of toluene and ethanol. We furthermore demonstrate that the failures of the computationally inexpensive Perdew-Burke-Ernzerhof (PBE) functional in describing some features of the excited state potential energy surface of the S₁ state of nile red can be corrected for in a straightforward fashion, relying only on a small number of calculations making use of more sophisticated range-separated hybrid functionals. The resulting solvatochromic shifts and predicted colours are in excellent agreement with experiment, showing the computational approach outlined in this work to yield very robust predictions of optical properties of dyes in solution.

Retinal isomerization and water-pore formation in channelrhodopsin-2

Channelrhodopsin-2 (ChR2) is a light-sensitive ion channel widely used in optogenetics. Photoactivation triggers a trans-to-cis isomerization of a covalently bound retinal. Ensuing conformational changes open a cation-selective channel. We explore the structural dynamics in the early photocycle leading to channel opening by classical (MM) and quantum mechanical (QM) molecular simulations. With QM/MM simulations, we generated a protein-adapted force field for the retinal chromophore, which we validated against absorption spectra. In a 4-µs MM simulation of a dark-adapted ChR2 dimer, water entered the vestibules of the closed channel. Retinal all-trans to 13-cis isomerization, simulated with metadynamics, triggered a major restructuring of the charge cluster forming the channel gate. On a microsecond time scale, water penetrated the gate to form a membrane-spanning preopen pore between helices H1, H2, H3, and H7. This influx of water into an ion-impermeable preopen pore is consistent with time-resolved infrared spectroscopy and electrophysiology experiments. In the retinal 13-cis state, D253 emerged as the proton acceptor of the Schiff base. Upon proton transfer from the Schiff base to D253, modeled by QM/MM simulations, we obtained an early-M/P2390-like intermediate. Rapid rotation of the unprotonated Schiff base toward the cytosolic side effectively prevents its reprotonation from the extracellular side. From MM and QM simulations, we gained detailed insight into the mechanism of ChR2 photoactivation and early events in pore formation. By rearranging the network of charges and hydrogen bonds forming the gate, water emerges as a key player in light-driven ChR2 channel opening.

Fine Tuning of Retinal Photoinduced Decay in Solution

Single methylation at position C₁₀ of the all-trans retinal protonated Schiff base switches its excited-state decay in methanol from a slower picosecond into an ultrafast, protein-like subpicosecond process. QM/MM modeling in conjunction with on-the-fly excited-state dynamics provides fundamental understanding of the fine-tuning mechanics that "catalyzes" the photoinduced decay of solvated retinals. Methylation alters the interplay between the ionic S₁ and covalent S₂ states, reducing the excited-state lifetime by favoring the formation of a S₁ transient fluorescent state with fully inverted bond lengths that accounts for the recorded transient spectroscopy and from which a space-saving conical intersection seam is quickly (<1 ps) reached. Minimal and apparently innocent chemical modifications thus affect the characteristic intramolecular charge-transfer of the S₁ state as well as the interaction with the covalent S₂ excited state, eventually providing the high tunability of retinal photophysics and photochemistry and delivering a new concept for the rational design of retinal-based photoactive molecular devices.

pH-dependent absorption spectra of rhodopsin mutant E113Q: On the role of counterions and protein

The absorption spectra of bovine rhodopsin mutant E113Q in solutions were investigated at the molecular level by using a hybrid quantum mechanics/molecular mechanics (QM/MM) method. The calculations suggest the mechanism of the absorption variations of E113Q at different pH values. The results indicate that the polarizations of the counterions in the vicinity of Schiff base under protonation and unprotonation states of the mutant E113Q would be a crucial factor to change the energy gap of the retinal to tune the absorption spectra. Glu-181 residue, which is close to the chromophore, cannot serve as the counterion of the protonated Schiff base of E113Q in dark state. Moreover, the results of the absorption maximum in mutant E113Q with the various anions (Cl^-, Br^-, I^- and NO₃^-) manifested that the mutant E113Q could have the potential for use as a template of anion biosensors at visible wavelength.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling

It was demonstrated in recent studies that some rhodopsins can be used in optogenetics as fluorescent indicators of membrane voltage. One of the promising candidates for these applications is archaerhodopsin-3. While it has already shown encouraging results, there is still a large room for improvement. One of possible directions is increasing the intensity of the protein's fluorescent signal. Rational design of mutants with an improved signal is an important task, which requires both experimental and theoretical studies. Herein, we used a homology-based computational approach to predict the three-dimensional structure of archaerhodopsin-3, and a Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid approach with high-level multireference ab initio methodology (SORCI+Q/AMBER) to model optical properties of this protein. We demonstrated that this methodology allows for reliable prediction of structure and spectral properties of archaerhodopsin-3. The results of this study can be utilized for computational molecular design of efficient fluorescent indicators of membrane voltage for modern optogenetics on the basis of archaerhodopsin-3.

Exploring what prompts ITIC to become a superior acceptor in organic solar cell by combining molecular dynamics simulation with quantum chemistry calculation

A comparison on charge transfer process in OSC between non-fullerene electron acceptor ITIC and PC71BM was taken by microscopic analysis based on the molecular dynamics and quantum chemistry method.

Protein–ligand interaction study of signal transducer smoothened protein with different drugs: molecular docking and QM/MM calculations

A part of signal transducer smoothened (SMO) protein including antitumor agent LY2940680. The site of this antitumor was considered for the docking of 716 ligands.

Assessing the correlation between mutant rhodopsin stability and the severity of retinitis pigmentosa

PURPOSE

Following a previous study that demonstrated a correlation between rhodopsin stability and the severity of retinitis pigmentosa (RP), we investigated whether predictions of severity can be improved with a regional analysis of this correlation. The association between changes to the stability of the protein and the relative amount of rhodopsin reaching the plasma membrane was assessed.

METHODS

Crystallography-based estimations of mutant rhodopsin stability were compared with descriptions in the scientific literature of the visual function of mutation carriers to determine the extent of associations between rhodopsin stability and clinical phenotype. To test the findings of this analysis, three residues of a green fluorescent protein (GFP) tagged rhodopsin plasmid were targeted with site-directed random mutagenesis to generate mutant variants with a range of stability changes. These plasmids were transfected into HEK-293 cells, and then flow cytometry was used to measure rhodopsin on the cells' plasma membrane. The GFP signal was used to measure the ratio between this membrane-bound rhodopsin and total cellular rhodopsin. FoldX stability predictions were then compared with the surface staining data and clinical data from the database to characterize the relationship between rhodopsin stability, the severity of RP, and the expression of rhodopsin at the cell surface.

RESULTS

There was a strong linear correlation between the scale of the destabilization of mutant variants and the severity of retinal disease. A correlation was also seen in vitro between stability and the amount of rhodopsin at the plasma membrane. Rhodopsin is drastically reduced on the surface of cells transfected with variants that differ in their inherent stability from the wild-type by more than 2 kcal/mol. Below this threshold, surface levels are closer to those of the wild-type.

CONCLUSIONS

There is a correlation between the stability of rhodopsin mutations and disease severity and levels of membrane-bound rhodopsin. Measuring membrane-bound rhodopsin with flow cytometry could improve prognoses for poorly characterized mutations and could provide a platform for measuring the effectiveness of treatments.

Comparison of the isomerization mechanisms of human melanopsin and invertebrate and vertebrate rhodopsins

Comparative modeling and ab initio multiconfigurational quantum chemistry are combined to investigate the reactivity of the human nonvisual photoreceptor melanopsin. It is found that both the thermal and photochemical isomerization of the melanopsin 11-cis retinal chromophore occur via a space-saving mechanism involving the unidirectional, counterclockwise twisting of the =C11H-C12H= moiety with respect to its Lys340-linked frame as proposed by Warshel for visual pigments [Warshel A (1976) Nature 260(5553):679-683]. A comparison with the mechanisms documented for vertebrate (bovine) and invertebrate (squid) visual photoreceptors shows that such a mechanism is not affected by the diversity of the three chromophore cavities. Despite such invariance, trajectory computations indicate that although all receptors display less than 100 fs excited state dynamics, human melanopsin decays from the excited state ∼40 fs earlier than bovine rhodopsin. Some diversity is also found in the energy barriers controlling thermal isomerization. Human melanopsin features the highest computed barrier which appears to be ∼2.5 kcal mol(-1) higher than that of bovine rhodopsin. When assuming the validity of both the reaction speed/quantum yield correlation discussed by Warshel, Mathies and coworkers [Weiss RM, Warshel A (1979) J Am Chem Soc 101:6131-6133; Schoenlein RW, Peteanu LA, Mathies RA, Shank CV (1991) Science 254(5030):412-415] and of a relationship between thermal isomerization rate and thermal activation of the photocycle, melanopsin turns out to be a highly sensitive pigment consistent with the low number of melanopsin-containing cells found in the retina and with the extraretina location of melanopsin in nonmammalian vertebrates.

Learning from photobiology how to design molecular devices using a computer

Learning how to model photo-responsive proteins may open the way to the design of lightpowered biomimetic molecular devices.

Modelling retinal chromophores photoisomerization: from minimal models in vacuo to ultimate bidimensional spectroscopy in rhodopsins

Modelling of retinal photoisomerization in different environments is reviewed and ultimate ultrafast electronic spectroscopy is proposed for obtaining new insights.

Generalized QM/MM Force Matching Approach Applied to the 11-cis Protonated Schiff Base Chromophore of Rhodopsin

We extended a previously developed force matching approach to systems with covalent QM/MM boundaries and describe its user-friendly implementation in the publicly available software package CPMD. We applied this approach to the challenging case of the retinal protonated Schiff base in dark state bovine rhodopsin. We were able to develop a highly accurate force field that is able to capture subtle structural changes within the chromophore that have a pronounced influence on the optical properties. The optical absorption spectrum calculated from configurations extracted from a MD trajectory using the new force field is in excellent agreement with QM/MM and experimental references.

Assessment of Density Functional Theory for Describing the Correlation Effects on the Ground and Excited State Potential Energy Surfaces of a Retinal Chromophore Model

In the quest for a cost-effective level of theory able to describe a large portion of the ground and excited potential energy surfaces of large chromophores, promising approaches are rooted in various approximations to the exact density functional theory (DFT). In the present work, we investigate how generalized Kohn-Sham DFT (GKS-DFT), time-dependent DFT (TDDFT), and spin-restricted ensemble-DFT (REKS) methods perform along three important paths characterizing a model retinal chromophore (the penta-2,4-dieniminium cation) in a region of near-degeneracy (close to a conical intersection) with respect to reference high-level multiconfigurational wave function methods. If GKS-DFT correctly describes the closed-shell charge transfer state, only TDDFT and REKS approaches give access to the open-shell diradical, one which sometimes corresponds to the electronic ground state. It is demonstrated that the main drawback of the usual DFT-based methods lies in the absence of interactions between the charge transfer and the diradicaloid configurations. Hence, we test a new computational scheme based on the State-averaged REKS (SA-REKS) approach, which explicitly includes these interactions into account. The State-Interaction SA-REKS (SI-SA-REKS) method significantly improves on the REKS and the SA-REKS results for the target system. The similarities and differences between DFT and wave function-based approaches are analyzed according to (1) the active space dimensions of the wave function-based methods and (2) the relative electronegativities of the allyl and protonated Schiff base moieties.

Selecting DFT methods for use in optimizations of enzyme active sites: applications to ONIOM treatments of DNA glycosylases

When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.

Rhodopsin Absorption from First Principles: Bypassing Common Pitfalls

Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.