Revealing the Effect of Stereocontrol on Intermolecular Interactions between Abiotic, Sequence-Defined Polyurethanes and a Ligand.

The development of precision polymer synthesis has facilitated access to a diverse library of abiotic structures wherein chiral monomers are positioned at specific locations within macromolecular chains. These structures are anticipated to exhibit folding characteristics similar to those of biotic macromolecules and possess comparable functionalities. However, the extensive sequence space and numerous variables make selecting a sequence with the desired function challenging. Therefore, revealing sequence-function dependencies and developing practical tools are necessary to analyze their conformations and molecular interactions. In this study, we investigate the effect of stereochemistry, which dictates the spatial location of backbone and pendant groups, on the interaction between sequence-defined oligourethanes and bisphenol A ligands. Various methods are explored to analyze the receptor-like properties of model oligomers and the ligand. The accuracy of molecular dynamics simulations and experimental techniques is assessed to uncover the impact of discrete changes in stereochemical arrangements on the structures of the resulting complexes and their binding strengths. Detailed computational investigations providing atomistic details show that the formed complexes demonstrate significant structural diversity depending on the sequence of stereocenters, thus affecting the oligomer-ligand binding strength. Among the tested techniques, the fluorescence spectroscopy data, fitted to the Stern-Volmer equation, are consistently aligned with the calculations, thus validating the developed simulation methodology. The developed methodology opens a way to engineer the structure of sequence-defined oligomers with receptor-like functionality to explore their practical applications, e.g., as sensory materials.

Experimental Assessment of the Electronic and Geometrical Structure of a Near-Infrared Absorbing and Highly Fluorescent Microbial Rhodopsin.

The recently discovered Neorhodopsin (NeoR) exhibits absorption and emission maxima in the near-infrared spectral region, which together with the high fluorescence quantum yield makes it an attractive retinal protein for optogenetic applications. The unique optical properties can be rationalized by a theoretical model that predicts a high charge transfer character in the electronic ground state (S0) which is otherwise typical of the excited state S1 in canonical retinal proteins. The present study sets out to assess the electronic structure of the NeoR chromophore by resonance Raman (RR) spectroscopy since frequencies and relative intensities of RR bands are controlled by the ground and excited state's properties. The RR spectra of NeoR differ dramatically from those of canonical rhodopsins but can be reliably reproduced by the calculations carried out within two different structural models. The remarkable agreement between the experimental and calculated spectra confirms the consistency and robustness of the theoretical approach.

Benchmarking two-photon absorption strengths of rhodopsin chromophore models with CC3 and CCSD methodologies: An assessment of popular density functional approximations.

This work presents the investigations of the impact of an increasing electron correlation in the hierarchy of coupled-cluster methods, i.e., CC2, CCSD, and CC3, on two-photon absorption (2PA) strengths for the lowest excited state of the minimal rhodopsin's chromophore model-cis-penta-2,4-dieniminium cation (PSB3). For a larger chromophore's model [4-cis-hepta-2,4,6-trieniminium cation (PSB4)], CC2 and CCSD calculations of 2PA strengths were performed. Additionally, 2PA strengths predicted by some popular density functional theory (DFT) functionals differing in HF exchange contribution were assessed against the reference CC3/CCSD data. For PSB3, the accuracy of 2PA strengths increases in the following order: CC2 < CCSD < CC3, with the CC2 deviation from both higher-level methods exceeding 10% at 6-31+G* basis sets and 2% at aug-cc-pVDZ basis set. However, for PSB4, this trend is reversed and CC2-based 2PA strength is larger than the corresponding CCSD value. Among the DFT functionals investigated, CAM-B3LYP and BHandHLYP provide 2PA strengths in best compliance with reference data, however, with the error approaching an order of magnitude.

Quantum chemistry study of the multiphoton absorption in enhanced green fluorescent protein at the single amino acid residue level.

The chromophore (CRO) of fluorescent proteins (FPs) is embedded in a complex environment that is a source of specific interactions with the CRO. Understanding how these interactions influence FPs spectral properties is important for a directed design of novel markers with desired characteristics. In this work, we apply computational chemistry methods to gain insight into one-, two- and three-photon absorption (1PA, 2PA, 3PA) tuning in enhanced green fluorescent protein (EGFP). To achieve this goal, we built EGFP models differing in: i) number and position of hydrogen-bonds (h-bonds) donors to the CRO and ii) the electric field, as approximated by polarizable force field, acting on the CRO. We find that h-bonding to the CRO's phenolate oxygen results in stronger one- and multiphoton absorption. The brighter absorption can be also achieved by creating more positive electric field near the CRO's phenolate moiety. Interestingly, while individual CRO - environment h-bonds usually enhance 1PA and 2PA, it takes a few h-bond donors to enhance 3PA. Clearly, response of the absorption intensity to many-body effects depends on the excitation mechanism. We further employ symmetry-adapted perturbation theory (SAPT) to reveal excellent (2PA) and good (3PA) correlation of multiphoton intensity with electrostatic and induction interaction energies. This points to importance of accounting for mutual CRO - environment polarization in quantitative calculations of absorption spectra in FPs.

Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency.

The activation of rhodopsin, the light-sensitive G-protein-coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. In this study we used hundreds of quantum-classical trajectories to show that, 15 fs after light absorption, a degeneracy between the reactive excited state and a neighbouring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show here that such splitting is modulated by protein electrostatics, thus linking amino acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.

The Role of Hydrogen Bonds and Electrostatic Interactions in Enhancing Two-Photon Absorption in Green and Yellow Fluorescent Proteins.

The spectral properties of fluorescent proteins (FPs) depend on the protein environment of the chromophore (CRO). A deeper understanding of the CRO - environment interactions in terms of FPs' spectral characteristics will allow for a rational design of novel markers with desired properties. Here, we are taking a step towards achieving this important goal. With the time-dependent density functional theory (TDDFT), we calculate one- and two-photon absorption (OPA and TPA) spectra for 5 green FPs (GFPs) and 3 yellow FPs (YFPs) differing in amino acid sequence. The goal is to reveal the roles of: (i) electrostatic interactions, (ii) hydrogen-bonds (h-bonds) and (iii) h-bonds together with distant electrostatic field in absorption spectra tuning. Our results point to design hypothesis towards FPs optimised for TPA-based applications. Both h-bonds and electrostatic interactions co-operate in enhancing TPA cross-section ( σTPA ) for the S0→S1 transition in GFPs. Furthermore, it seems that details of h-bonds network in the CRO's vicinity influences σTPA response to CRO - environment electrostatic interactions in YFPs. We postulate that engineering FPs with more hydrophilic CRO's environment can lead to greater σTPA . We also find that removing h-bonds formed with the CRO's phenolate leads to TPA enhancement for transition to higher excited states than S1 . Particularly Y145 and T203 residues are important in this regard.

Recognition of the True and False Resonance Raman Optical Activity.

Resonance Raman optical activity (RROA) possesses all aspects of a sensitive tool for molecular detection, but its measurement remains challenging. We demonstrate that reliable recording of RROA of chiral colorful compounds is possible, but only after considering the effect of the electronic circular dichroism (ECD) on the ROA spectra induced by the dissolved chiral compound. We show RROA for a number of model vitamin B12 derivatives that are chemically similar but exhibit distinctively different spectroscopic behavior. The ECD/ROA effect is proportional to the concentration and dependent on the optical pathlength of the light propagating through the sample. It can severely alter relative band intensities and signs in the natural RROA spectra. The spectra analyses are supported by computational modeling based on density functional theory. Neglecting the ECD effect during ROA measurement can lead to misinterpretation of the recorded spectra and erroneous conclusions about the molecular structure.

Removing artifacts in polarizable embedding calculations of one- and two-photon absorption spectra of fluorescent proteins.

The multiscale calculations involving excited states may suffer from the electron spill-out (ESO) problem. This seems to be especially the case when the environment of the core region, described with the electronic structure method, is approximated by a polarizable force field. The ESO effect often leads to incorrect physical character of electronic excitations, spreading outside the quantum region, which, in turn, results in erroneous absorption spectra. In this work, we investigate means to remove the artifacts in one-photon absorption (OPA) and two-photon absorption (TPA) spectra of green and yellow fluorescent protein representatives. This includes (i) using different basis sets, (ii) extending the core subsystem beyond the chromophore, (iii) modification of polarization interaction between the core region and its environment, and (iv) including the Pauli repulsion through effective core potentials (ECPs). Our results clearly show that ESO is observed when diffuse functions are used to assemble the multielectron wave function regardless of the exchange-correlation functional used. Furthermore, extending the core region, thus accounting for exchange interactions between the chromophore and its environment, leads to even more spurious excited states. Also, damping the interactions between the core subsystem and the polarizable force field is hardly helpful. In contrast, placing ECPs in the position of sites creating the embedding potential leads to the removal of artificious excited states that presumably should not be observed in the OPA and TPA spectra. We prove that it is a reliable and cost-effective approach for systems where the covalent bond(s) between the core region and its environment must be cut.

What is the Optimal Size of the Quantum Region in Embedding Calculations of Two-Photon Absorption Spectra of Fluorescent Proteins?

We systematically investigate an impact of the size and content of a quantum (QM) region, treated at the density functional theory level, in embedding calculations on one- (OPA) and two-photon absorption (TPA) spectra of the following fluorescent proteins (FPs) models: Aequorea victoria green FP (avGFP) with neutral (avGFP-n) and anionic (avGFP-a) chromophore as well as Citrine FP. We find that amino acid (a.a.) residues as well as water molecules hydrogen-bonded (h-bonded) to the chromophore usually boost both OPA and TPA processes intensity. The presence of hydrophobic a.a. residues in the quantum region also non-negligibly affects both absorption spectra but decreases absorption intensity. We conclude that to reach a quantitative description of OPA and TPA spectra in multiscale modeling of FPs, the quantum region should consist of a chromophore and most of a.a. residues and water molecules in a radius of 0.30-0.35 nm (ca. 200-230 atoms) when the remaining part of the system is approximated by the electrostatic point-charges. The optimal size of the QM region can be reduced to 80-100 atoms by utilizing a more advanced polarizable embedding model. We also find components of the QM region that are specific to a FP under study. We propose that the F165 a.a. residue is important in tuning the TPA spectrum of avGFP-n but not other investigated FPs. In the case of Citrine, Y203 and M69 a.a. residues must definitely be part of the QM subsystem. Furthermore, we find that long-range electrostatic interactions between the QM region and the rest of the protein cannot be neglected even for the most extensive QM regions (ca. 350 atoms).

Flavones' and Flavonols' Antiradical Structure-Activity Relationship-A Quantum Chemical Study.

Flavonoids are known for their antiradical capacity, and this ability is strongly structure-dependent. In this research, the activity of flavones and flavonols in a water solvent was studied with the density functional theory methods. These included examination of flavonoids' molecular and radical structures with natural bonding orbitals analysis, spin density analysis and frontier molecular orbitals theory. Calculations of determinants were performed: specific, for the three possible mechanisms of action-hydrogen atom transfer (HAT), electron transfer-proton transfer (ETPT) and sequential proton loss electron transfer (SPLET); and the unspecific-reorganization enthalpy (RE) and hydrogen abstraction enthalpy (HAE). Intramolecular hydrogen bonding, catechol moiety activity and the probability of electron density swap between rings were all established. Hydrogen bonding seems to be much more important than the conjugation effect, because some structures tends to form more intramolecular hydrogen bonds instead of being completely planar. The very first hydrogen abstraction mechanism in a water solvent is SPLET, and the most privileged abstraction site, indicated by HAE, can be associated with the C3 hydroxyl group of flavonols and C4' hydroxyl group of flavones. For the catechol moiety, an intramolecular reorganization to an o-benzoquinone-like structure occurs, and the ETPT is favored as the second abstraction mechanism.

Assessment of Functionals for TDDFT Calculations of One- and Two-Photon Absorption Properties of Neutral and Anionic Fluorescent Proteins Chromophores.

Performance of DFT functionals with different percentages of exact Hartree-Fock exchange energy (EX) is assessed for recovery of the CC2 reference one- (OPA) and two-photon absorption (TPA) spectra of fluorescent proteins chromophores in vacuo. The investigated DFT functionals, together with their EX contributions are BLYP (0%), B3LYP (20%), B1LYP (25%), BHandHLYP (50%), and CAM-B3LYP (19% at short range and 65% at long range). Our test set consists of anionic and neutral chromophores as naturally occurring in the fluorescent proteins. For the first time, we compare TDDFT and CC2 methods for higher excited states than the S1 state, exhibiting relatively large TPA intensity. Our TDDFT results for neutral chromophores reveal an increase in excitation energies as well as TPA and OPA intensities errors, compared to CC2-derived results, as the DFT functional contains less exact exchange. The long-range-corrected CAM-B3LYP functional performs the best, closely followed by BHandHLYP, while BLYP usually significantly underestimates all investigated spectral properties, hence being the worst in reproducing the reference CC2 results. The hybrid B3LYP and B1LYP functionals can be roughly placed in between. We propose that TDDFT may underestimate the TPA intensities for neutral chromophores of fluorescent proteins due to underestimated oscillator strengths between some excited states. In the case of anionic chromophores, we find that B3LYP and B1LYP functionals overcome others in terms of reproducing CC2 excitation energies. On the other hand, however, TPA intensity is usually significantly underestimated, and in this respect, CAM-B3LYP functional seems to be again superior. In contrast to the case of neutral chromophores, it seems that a large magnitude of excited-state dipole moments or changes in dipole moments upon excitation may be the driving force behind high TPA transition moments.

Initial excited-state relaxation of locked retinal protonated schiff base chromophore. An insight from coupled cluster and multireference perturbation theory calculations.

The initial S1 excited-state relaxation of retinal protonated Schiff base (RPSB) analog with central C11C12 double bond locked by eight-membered ring (locked-11.8) was investigated by means of multireference perturbation theory methods (XMCQDPT2, XMS-CASPT2, MS-CASPT2) as well as single-reference coupled-cluster CC2 method. The analysis of XMCQDPT2-based geometries reveals rather weak coupling between in-plane and out-of-plane structural evolution and minor energetical relaxation of three locked-11.8 conformers. Therefore, a strong coupling between bonds length inversion and backbone out-of-plane deformation resulting in a very steep S1 energy profile predicted by CASSCF/CASPT2 calculations is in clear contradiction with the reference XMCQDPT2 results. Even though CC2 method predicts good quality ground-state structures, the excited-state structures display more advanced torsional deformation leading to ca. 0.2 eV exaggerated energy relaxation and significantly red shifted (0.4-0.7 eV) emission maxima. According to our findings, the initial photoisomerization process in locked-11.8, and possibly in other RPSB analogs, studied fully (both geometries and energies) by multireference perturbation theory may be somewhat slower than predicted by CASSCF/CASPT2 or CC2 methods. © 2018 Wiley Periodicals, Inc.

Is the choice of a standard zeroth-order hamiltonian in CASPT2 ansatz optimal in calculations of excitation energies in protonated and unprotonated schiff bases of retinal?

To account for systematic error of CASPT2 method empirical modification of the zeroth-order Hamiltonian with Ionization Potential-Electron Affinity (IPEA) shift was introduced. The optimized IPEA value (0.25 a.u.), called standard IPEA (S-IPEA), was recommended but due to its unsatisfactory performance in multiple metallic and organic compounds it has been questioned lately as a general parameter working properly for all molecules under CASPT2 study. As we are interested in Schiff bases of retinal, an important question emerging from this conflict of choice, to use or not to use S-IPEA, is whether the introduction of the modified zeroth-order Hamiltonian into CASPT2 ansatz does really improve their energetics. To achieve this goal, we assessed an impact of the IPEA shift value, in a range of 0-0.35 a.u., on vertical excitation energies to low-lying singlet states of two protonated (RPSBs) and two unprotonated (RSBs) Schiff bases of retinal for which experimental data in gas phase are available. In addition, an effect of geometry, basis set, and active space on computed VEEs is also reported. We find, that for these systems, the choice of S-IPEA significantly overestimates both S0 →S1 and S0 →S2 energies and the best theoretical estimate, in reference to the experimental data, is provided with either unmodified zeroth-order Hamiltonian or small value of the IPEA shift in a range of 0.05-0.15 a.u., depending on active space and basis set size, equilibrium geometry, and character of the excited state. © 2018 Wiley Periodicals, Inc.

Excited-state minima and emission energies of retinal chromophore analogues: Performance of CASSCF and CC2 methods as compared with CASPT2.

This study provides gas-phase S1 excited-state geometries along with emission and adiabatic energies for methylated/demethylated and ring-locked analogues of protonated Schiff base retinal models comprising system of five conjugated double bonds (PSB5), using second order multiconfiguration perturbation theory (CASPT2). CASPT2 results serve as reference data to assess the performance of CC2 (second-order approximate coupled cluster singles and doubles) and a commonly used CASSCF/CASPT2 protocol, that is, complete active space self-consistent field (CASSCF) geometry optimization followed by CASPT2 energy calculation. We find that the CASSCF methodology fails to locate planar S1 minimum energy structures for four out of five investigated planar models in contrast to CC2 and CASPT2 methods. However, for those which were found: one planar and two twisted minima, there is an excellent agreement between CASSCF and CASPT2 results in terms of geometrical parameters, one-electron properties, as well as emission and adiabatic energies. CC2 performs well for in-plane S1 minima and their spectroscopic and electronic properties. However, this picture deteriorates for twisted minima. As expected, the CC2 description of the S2 electronic state, with strong multireference and significant double excitation character, is very poor, exhibiting errors in transition energies exceeding 1 eV. They may be substantially diminished by recalculating transition energies with CASPT2 method. Our work shows that CASSCF/CASPT2 and CC2 shortcomings may influence gas-phase retinal analogues' excited state description in a dramatic way. © 2017 Wiley Periodicals, Inc.

Assessing the Accuracy of Various Ab Initio Methods for Geometries and Excitation Energies of Retinal Chromophore Minimal Model by Comparison with CASPT3 Results.

The effect of the quality of the ground-state geometry on excitation energies in the retinal chromophore minimal model (PSB3) was systematically investigated using various single- (within Møller-Plesset and coupled-cluster frameworks) and multiconfigurational [within complete active space self-consistent field (CASSCF) and CASSCF-based perturbative approaches: second-order CASPT2 and third-order CASPT3] methods. Among investigated methods, only CASPT3 provides geometry in nearly perfect agreement with the CCSD(T)-based equilibrium structure. The second goal of the present study was to assess the performance of the CASPT2 methodology, which is popular in computational spectroscopy of retinals, in describing the excitation energies of low-lying excited states of PSB3 relative to CASPT3 results. The resulting CASPT2 excitation energy error is up to 0.16 eV for the S0 → S1 transition but only up to 0.06 eV for the S0 → S2 transition. Furthermore, CASPT3 excitation energies practically do not depend on modification of the zeroth-order Hamiltonian (so-called IPEA shift parameter), which does dramatically and nonsystematically affect CASPT2 excitation energies.

The role of spin-orbit coupling in the photolysis of methylcobalamin.

The photolysis of the methylcobalamin cofactor (MeCbl) in its base-off form was investigated by considering the extent of spin-orbit coupling (SOC). Triplet Co-C photodissociation pathways previously invoked at the density functional theory level using Landau-Zener theory were further validated with ab initio calculations that combine SOC based on multi-state second order perturbation theory. It was determined that SOC is feasible between singlet and triplet states at elongated Co-C distances, leading to photodissociation from the state having dominant σ(dz(2)) character, by either direct coupling with the lowest singlet states or by crossing with SOC mixed triplets.

Assessment of Approximate Coupled-Cluster and Algebraic-Diagrammatic-Construction Methods for Ground- and Excited-State Reaction Paths and the Conical-Intersection Seam of a Retinal-Chromophore Model.

As a minimal model of the chromophore of rhodopsin proteins, the penta-2,4-dieniminium cation (PSB3) poses a challenging test system for the assessment of electronic-structure methods for the exploration of ground- and excited-state potential-energy surfaces, the topography of conical intersections, and the dimensionality (topology) of the branching space. Herein, we report on the performance of the approximate linear-response coupled-cluster method of second order (CC2) and the algebraic-diagrammatic-construction scheme of the polarization propagator of second and third orders (ADC(2) and ADC(3)). For the ADC(2) method, we considered both the strict and extended variants (ADC(2)-s and ADC(2)-x). For both CC2 and ADC methods, we also tested the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) variants. We have explored several ground- and excited-state reaction paths, a circular path centered around the S1/S0 surface crossing, and a 2D scan of the potential-energy surfaces along the branching space. We find that the CC2 and ADC methods yield a different dimensionality of the intersection space. While the ADC methods yield a linear intersection topology, we find a conical intersection topology for the CC2 method. We present computational evidence showing that the linear-response CC2 method yields a surface crossing between the reference state and the first response state featuring characteristics that are expected for a true conical intersection. Finally, we test the performance of these methods for the approximate geometry optimization of the S1/S0 minimum-energy conical intersection and compare the geometries with available data from multireference methods. The present study provides new insight into the performance of linear-response CC2 and polarization-propagator ADC methods for molecular electronic spectroscopy and applications in computational photochemistry.

Impacts of retinal polyene (de)methylation on the photoisomerization mechanism and photon energy storage of rhodopsin.

Ab initio multiconfigurational quantum chemical methodology combined with molecular mechanics (CASPT2//CASSCF/AMBER) was applied to probe impacts of retinal protonated Schiff base (RPSB) polyene methylation and/or demethylation on the mechanism of photochemical isomerization in bovine rhodopsin. We have examined structural and spectroscopic properties of wild-type rhodopsin (with 11-cis-9,13-dimethyl-RPSB) and artificial rhodopsins, hosting four 11-cis-RPSB derivatives, 13-demethyl-, 9-demethyl-, 10-methyl-13-demethyl-, and 10-methyl-RPSB, evolving along the photoisomerization coordinate. It is found that the addition of 10-methyl or/and deletion of 9-/13-methyl groups do not appear to interfere structurally with the photoisomerization pathway in the S1 excited state. Remarkably, the two-mode space-saving mechanism initiated by bond order inversion and followed by asynchronous bicycle-pedal distortion in the RPSB backbone drives the photoreaction in all rhodopsin analogues studied here. However, methylation and/or demethylation is responsible for perturbation of excess energy deposited in the conical intersection structures. The analysis of photon energy stored by bathorhodopsin in synthetic pigments reveals that it is affected by steric crowding of methyl substituents in the RPSB backbone.

Mechanism of Co-C bond photolysis in the base-on form of methylcobalamin.

A mechanism of Co-C bond photodissociation in the base-on form of the methylcobalamin cofactor (MeCbl) has been investigated employing time-dependent density functional theory (TD-DFT), in which the key step involves singlet radical pair generation from the first electronically excited state (S1). The corresponding potential energy surface of the S1 state was constructed as a function of Co-C and Co-Naxial bond distances, and two possible photodissociation pathways were identified on the basis of energetic grounds. These pathways are distinguished by whether the Co-C bond (path A) or Co-Naxial bond (path B) elongates first. Although the final intermediate of both pathways is the same (namely a ligand field (LF) state responsible for Co-C dissociation), the reaction coordinates associated with paths A and B are different. The photolysis of MeCbl is wavelength-dependent, and present TD-DFT analysis indicates that excitation in the visible α/ß band (520 nm) can be associated with path A, whereas excitation in the near-UV region (400 nm) is associated with path B. The possibility of intersystem crossing, and internal conversion to the ground state along path B are also discussed. The mechanism proposed in this study reconciles existing experimental data with previous theoretical calculations addressing the possible involvement of a repulsive triplet state.

Initial excited-state dynamics of an N-alkylated indanylidene-pyrroline (NAIP) rhodopsin analog.

N-Alkylated indanylidene-pyrroline-based molecular switches mimic different aspects of the light-induced retinal chromophore isomerization in rhodopsin: the vertebrate dim-light visual pigment. In particular, they display a similar ultrashort excited-state lifetime, subpicosecond photoproduct appearance time, and photoproduct vibrational coherence. To better understand the early light-induced dynamics of such systems, we measured and modeled the resonance Raman spectra of the Z-isomer of the N-methyl-4-(5'-methoxy-2',2'-dimethyl-indan-1'-ylidene)-5-methyl-2,3-dihydro-2H-pyrrolium (NAIP) switch in methanol solution. It is shown that the data, complemented with a <70 fs excited-state trajectory computation, demonstrate initial excited-state structural dynamics dominated by double-bond expansion and single-bond contraction stretches. This mode subsequently couples with the five-membered ring inversion and double-bond torsion. These results are discussed in the context of the mechanism of the excited-state photoisomerization of NAIP switches in solution and the 11-cis retinal in rhodopsin.