Piero Altoè

Conical intersection dynamics of the primary photoisomerization event in vision

Ever since the conversion of the 11-cis retinal chromophore to its all-trans form in rhodopsin was identified as the primary photochemical event in vision, experimentalists and theoreticians have tried to unravel the molecular details of this process. The high quantum yield of 0.65 (ref. 2), the production of the primary ground-state rhodopsin photoproduct within a mere 200 fs (refs 3–7), and the storage of considerable energy in the first stable bathorhodopsin intermediate all suggest an unusually fast and efficient photoactivated one-way reaction. Rhodopsins unique reactivity is generally attributed to a conical intersection between the potential energy surfaces of the ground and excited electronic states enabling the efficient and ultrafast conversion of photon energy into chemical energy. But obtaining direct experimental evidence for the involvement of a conical intersection is challenging: the energy gap between the electronic states of the reacting molecule changes significantly over an ultrashort timescale, which calls for observational methods that combine high temporal resolution with a broad spectral observation window. Here we show that ultrafast optical spectroscopy with sub-20-fs time resolution and spectral coverage from the visible to the near-infrared allows us to follow the dynamics leading to the conical intersection in rhodopsin isomerization. We track coherent wave-packet motion from the photoexcited Franck–Condon region to the photoproduct by monitoring the loss of reactant emission and the subsequent appearance of photoproduct absorption, and find excellent agreement between the experimental observations and molecular dynamics calculations that involve a true electronic state crossing. Taken together, these findings constitute the most compelling evidence to date for the existence and importance of conical intersections in visual photochemistry.

Electrostatic Control of the Photoisomerization Efficiency and Optical Properties in Visual Pigments: On the Role of Counterion Quenching

Hybrid QM(CASPT2//CASSCF/6-31G*)/MM(Amber) computations have been used to map the photoisomerization path of the retinal chromophore in Rhodopsin and explore the reasons behind the photoactivity efficiency and spectral control in the visual pigments. It is shown that while the electrostatic environment plays a central role in properly tuning the optical properties of the chromophore, it is also critical in biasing the ultrafast photochemical event: it controls the slope of the photoisomerization channel as well as the accessibility of the S(1)/S(0) crossing space triggering the ultrafast decay. The roles of the E113 counterion, the E181 residue, and the other amino acids of the protein pocket are explicitly analyzed: it appears that counterion quenching by the protein environment plays a key role in setting up the chromophores optical properties and its photochemical efficiency. A unified scenario is presented that discloses the relationship between spectroscopic and mechanistic properties in rhodopsins and allows us to draw a solid mechanism for spectral tuning in color vision pigments: a tunable counterion shielding appears as the elective mechanism for L<-->M spectral modulation, while a retinal conformational control must dictate S absorption. Finally, it is suggested that this model may contribute to shed new light into mutations-related vision deficiencies that opens innovative perspectives for experimental biomolecular investigations in this field.

Aborted double bicycle-pedal isomerization with hydrogen bond breaking is the primary event of bacteriorhodopsin proton pumping

Quantum mechanics/molecular mechanics calculations based on ab initio multiconfigurational second order perturbation theory are employed to construct a computer model of Bacteriorhodopsin that reproduces the observed static and transient electronic spectra, the dipole moment changes, and the energy stored in the photocycle intermediate K. The computed reaction coordinate indicates that the isomerization of the retinal chromophore occurs via a complex motion accounting for three distinct regimes: (i) production of the excited state intermediate I, (ii) evolution of I toward a conical intersection between the excited state and the ground state, and (iii) formation of K. We show that, during stage ii, a space-saving mechanism dominated by an asynchronous double bicycle-pedal deformation of the C10═C11─C12═C13─C14═N moiety of the chromophore dominates the isomerization. On this same stage a N─H/water hydrogen bond is weakened and initiates a breaking process that is completed during stage iii.

Adenine deactivation in DNA resolved at the CASPT2//CASSCF/AMBER level

We have employed hybrid CASPT2//CASSCF/AMBER calculations to map the (1)L(a)(1pipi*) deactivation path of a single quantum mechanical adenine in a d(A)(10).d(T)(10) double strand in water that is treated at the molecular mechanics level. We find that (a) the L(a) relaxation route is flatter in DNA than in vacuo and (b) the L(a) relaxation energy in DNA is much larger than the stabilization energy of the corresponding L(a) excimer. An intra-monomer relaxation process is found to be compatible with the multiexponential decay recorded in DNA, possibly including the longer (4100 ps) lifetime component.

Multistate Photo-Induced Relaxation and Photoisomerization Ability of Fumaramide Threads: A Computational and Experimental Study

Fumaric and maleic amides are the photoactive units of an important and widely investigated class of photocontrollable rotaxanes as they trigger ring shuttling via a cis-trans photoisomerization. Here, ultrafast decay and photoinduced isomerization in isolated fumaramide and solvated nitrogen-substituted fumaramides (that are employed as threads in those rotaxanes) have been investigated by means of CASPT2//CASSCF computational and time-resolved spectroscopic techniques, respectively. A complex multistate network of competitive deactivation channels, involving both internal conversion and intersystem crossing (ISC) processes, has been detected and characterized that accounts for the picosecond decay and photochemical/photophysical properties observed in the singlet as well as triplet (photosensitized) photochemistry of fumaramides threads. Interestingly, singlet photochemistry appears to follow a non-Kasha rule model, where nonequilibrium dynamical factors control the outcome of the photochemical process: accessible high energy portions of extended crossing seams turn out to drive the deactivation process and ground-state recovery. Concurrently, extended singlet/triplet degenerate regions of twisted molecular structures with significant spin-orbit-coupling values account for ultrafast (picosecond time scale) ISC processes that lead to higher photoisomerization efficiencies. This model discloses the principles behind the intrinsic photochemical reactivity of fumaramide and its control.

Deciphering Intrinsic Deactivation/Isomerization Routes in a Phytochrome Chromophore Model

High level ab initio correlated (CASPT2) computations have been used to elucidate the details of the photoinduced molecular motion and decay mechanisms of a realistic phytochrome chromophore model in vacuo and to explore the reasons underneath its photophysical/photochemical properties. Competitive deactivation routes emerge that unveil the primary photochemical event and the intrinsic photoisomerization ability of this system. The emerged in vacuo based static (i.e., nondynamical) reactivity model accounts for the formation of different excited state intermediates and suggests a qualitative rationale for the short (picosecond) excited state lifetime and ultrafast decay of the emission, its small quantum yield, and the multiexponential decay observed in both solvent and phytochromes. It is thus tentatively suggested that this is a more general deactivation scheme for photoexcited phytochrome chromophores that is independent of the surrounding environment. Spectroscopic properties have also been simulated in both isolated conditions and the protein that satisfactorily match experimental data. For this purpose, preliminary hybrid QM/MM computations at the correlated (CASPT2) level have been used in the protein and are reported here for the first time.

Catalytic Mechanism of Diaminopimelate Epimerase: A QM/MM Investigation.

A QM/MM investigation, based on a DFT(B3LYP)//Amber-ff99 potential, has been carried out to elucidate the mechanism of diaminopimelate epimerase. This enzyme catalyzes the reversible stereoconversion of one of the two stereocenters of diaminopimelate and represents a promising target for rational drug design aimed to develop new selective antibacterial therapeutic agents. The QM/MM computations show that the reaction proceeds through a highly asynchronous mechanism where the side-chain of a negatively charged Cys-73 (thiolate) deprotonates the α-carbon substrate. Simultaneously, the Cys-217 thiolic proton moves toward the same carbon atom on the opposite face, thus determining the configuration inversion. A fingerprint analysis provides a detailed description of the influence of the various residues surrounding the active site and clearly shows the electrostatic nature of the most important contributions to the catalysis.

Spectroscopic Properties of Formaldehyde in Aqueous Solution: Insights from Car-Parrinello and TDDFT/CASPT2 Calculations.

We present Car-Parrinello and Car-Parrinello/molecular mechanics simulations of the structural, vibrational, and electronic properties of formaldehyde in water. The calculated properties of the molecule reproduce experimental values and previous calculations. The n → π* excitation energy, calculated with TDDFT and CASPT2, agrees with experimental data. In particular, it shows a blue shift on going from the gas phase to aqueous solution. Temperature and wave function polarization contributions have been disentangled.

COBRAMM: A Tunable QM/MM Approach to Complex Molecular Architectures. Modelling the Excited and Ground State Properties of Sized Molecular Systems

This contribution describes a new implementation of a general hybrid approach with a modular structure (called COBRAMM: Computations in Bologna Relating Ab‐initio and Molecular Mechanics Methods) that is able to integrate some specialized softwares and acts as a flexible computational environment, thus increasing the flexibility/efficiency of both QM, and MM, and QM/MM calculations. Specifically, QM/MM ground and excited states geometry optimizations, frequency calculations, conical intersection searches and adiabatic/non‐adiabatic molecular dynamics can be performed on a large molecular system, that can be split up to three different layers corresponding to different levels of accuracy. Here we report, together with a description of the method and its implementation, some test examples on very different chemical problems, which span the wide and diversified area of chemistry (from ground to excited states topics) and show the flexibility, general applicability and accuracy of the presented hybrid approac...

Rhodopsin and GFP Chromophores: QM/MM Absorption Spectra in Solvent and Protein

This work presents a QM/MM investigation of the spectral properties of the 11‐cis‐retinal and the GFP chromophore in polar solvents and protein. The results of the computations are in surprising agreement with the experimental values indicating the accuracy of our computational approach. In addition it has been demonstrated the key role of the solvent to create a “virtual conter‐ion”. This effect is due to the reorientation (i.e. polarization) of the polar solvent close to the chromophore: the polarized permanent dipoles of the solvent act similarly to the counter ion, stabilizing the ground state respect to the charge transfer excited state.