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. Rhodopsin's 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.

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.

Product formation in rhodopsin by fast hydrogen motions.

The photochemical cis-trans isomerization of retinal in rhodopsin is investigated by structure sampling and excited state QM/MM trajectories with surface hopping. The calculations uncover the motions responsible for photoproduct formation and elucidate the reasons behind the efficient photoisomerization in the primary event of visual transduction.

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.

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 chromophore's 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.

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.

The catalytic activity of proline racemase: a quantum mechanical/molecular mechanical study.

The enzyme proline racemase from the eukaryotic parasite Trypanosoma cruzi (responsible for endemic Chagas disease) catalyzes the reversible stereoinversion of chiral Calpha in proline. We employed a new combined quantum mechanical and molecular mechanical (QM/MM) potential to study the reaction mechanism of the enzyme. Three critical points were found: two almost isoenergetic minima (M1a and M2a), in which the enzyme is bound to L- and D-Pro, respectively, and a transition state (TSCa), unveiling a highly asynchronous concerted process. A systematic analysis was performed on the optimized geometries to point out the key role played by some residues in stabilizing the transition state.

COBRAMM 2.0 - A software interface for tailoring molecular electronic structure calculations and running nanoscale (QM/MM) simulations.

We present a new version of the simulation software COBRAMM, a program package interfacing widely known commercial and academic software for molecular modeling. It allows a problem-driven tailoring of computational chemistry simulations with effortless ground and excited-state electronic structure computations. Calculations can be executed within a pure QM or combined quantum mechanical/molecular mechanical (QM/MM) framework, bridging from the atomistic to the nanoscale. The user can perform all necessary steps to simulate ground state and photoreactions in vacuum, complex biopolymer, or solvent environments. Starting from ground-state optimization, reaction path computations, initial conditions sampling, spectroscopy simulation, and photodynamics with deactivation events, COBRAMM is designed to assist in characterization and analysis of complex molecular materials and their properties. Interpretation of recorded spectra range from steady-state to time-resolved measurements. Various tools help the user to set up the system of interest and analyze the results.

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.

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.

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.

Structure, spectroscopy, and spectral tuning of the gas-phase retinal chromophore: the beta-ionone "handle" and alkyl group effect.

The low-lying singlet states (i.e. S0, S1, and S2) of the chromophore of rhodopsin, the protonated Schiff base of 11-cis-retinal (PSB11), and of its all-trans photoproduct have been studied in isolated conditions by using ab initio multiconfigurational second-order perturbation theory. The computed spectroscopic features include the vertical excitation, the band origin, and the fluorescence maximum of both isomers. On the basis of the S0-->S1 vertical excitation, the gas-phase absorption maximum of PSB11 is predicted to be 545 nm (2.28 eV). Thus, the predicted absorption maximum appears to be closer to that of the rhodopsin pigment (2.48 eV) and considerably red-shifted with respect to that measured in solution (2.82 eV in methanol). In addition, the absorption maxima associated with the blue, green, and red cone visual pigments are tentatively rationalized in terms of the spectral changes computed for PSB11 structures featuring differently twisted beta-ionone rings. More specifically, a blue-shifted absorption maximum is explained in terms of a large twisting of the beta-ionone ring (with respect to the main conjugated chain) in the visual S-cone (blue) pigment chromophore. In contrast, the chromophore of the visual L-cone (red) pigment is expected to have a nearly coplanar beta-ionone ring yielding a six double bond fully conjugated framework. Finally, the M-cone (green) chromophore is expected to feature a twisting angle between 10 and 60 degrees. The spectroscopic effects of the alkyl substituents on the PSB11 spectroscopic properties have also been investigated. It is found that they have a not negligible stabilizing effect on the S1-S0 energy gap (and, thus, cause a red shift of the absorption maximum) only when the double bond of the beta-ionone ring conjugates significantly with the rest of the conjugated chain.

Solvent effects on the vibrational activity and photodynamics of the green fluorescent protein chromophore: a quantum-chemical study.

Vibrational activities in the Raman and resonance Raman spectra of the cationic, neutral, and anionic forms of 4'-hydroxybenzylidene-2,3-dimethyl-imidazolinone, a model compound for the green fluorescent protein chromophore, have been obtained from quantum-chemical calculations in vacuo and with the inclusion of solvent effects through the polarizable continuum model. It is found that inclusion of solvent effects improves slightly the agreement with experimental data for the cationic and neutral forms, whose spectra are qualitatively well-described already by calculations in vacuo. In contrast, inclusion of solvent effects is crucial to reproduce correctly the activities of the anionic form. The structural effects of solvation are remarkable both in the ground and in the lowest excited state of the anionic chromophore and influence not only the vibrational activity but also the photodynamics of the lowest excited state. CASPT2//CASSCF photoreaction paths, computed by including solvent effects at the CASSCF level, indicate a facile torsional deformation around both exocyclic CC bonds. Rotation around the exocyclic CC double bond is shown to lead to a favored radiationless decay channel, more efficient than that in gas phase, and which explains the ultrafast fluorescence decay and ground-state recovery observed in solution. Conversely, rotation around the exocyclic CC single bond accounts for the bottleneck observed in the ground-state recovery cycle. It is also speculated that the ultrafast radiationless decay channel would be hampered in protein for unfavorable electrostatic interactions and steric reasons.