Unraveling the Mechanisms of Nonradiative Deactivation in Model Peptides Following Photoexcitation of a Phenylalanine Residue

Isomeric and rotational effects in the chemi-ionisation of 1,2-dibromoethene with metastable neon atoms

The specific geometry of a molecule can have a pronounced influence on its chemical reactivity. However, experimental data on reactions of individual molecular isomers are still sparse because they are often difficult to separate and frequently interconvert into one another under ambient conditions. Here, we employ a novel crossed-beam experiment featuring an electrostatically controlled molecular beam combined with a source for radicals and metastables to spatially separate the cis and trans stereoisomers as well as individual rotational states of 1,2-dibromoethene and study their specific reactivities in the chemi-ionisation reaction with excited neon atoms. The experiments reveal pronounced isomeric and rotational specificities in the rates and product branching ratios of the reaction. The present study underlines the importance and combined role of molecular geometry and of rotational motion in the dynamics of chemi-ionisation reactions.

Generalized Energy-Based Fragmentation Approach for the Electronic Emission Spectra of Large Systems

The excited-state (ES) geometry optimization and electronic emission (fluorescence and phosphorescence) spectra and the ES vibrational spectra of large systems are great challenges in quantum chemistry. In this work, we develop a generalized energy-based fragmentation (GEBF) approach to compute the localized ES structures and vibrational frequencies of large systems. In this approach, the ES energy derivatives (gradients or Hessians) for a localized ES of a large system can be obtained by combining the ES energy derivatives of the corresponding active subsystems (including local excitation center) and the ground-state energy derivatives of inactive subsystems. Two strategies are adopted to overcome two difficulties from state-classification and state-tracking for treating specific ESs. First, for state-classification, we develop an improved density-based spatial clustering applied with noise algorithm with a modified transition orbital projection (TOP) algorithm, which allow a certain ES energy and energy derivatives of the whole system to be calculated with different ES energies and energy derivatives of active subsystems. Furthermore, we also employ the TOP algorithm for tracking the ESs in their geometry optimizations at the time-dependent density functional theory (TDDFT) level. Then, the GEBF approach is applied to investigate the optimized ES geometries or ES vibrational frequencies for two typical systems. Our results show that the cost-effective GEBF approach can accurately reproduce the TDDFT fluorescence spectra of the cytosine derivative and the experimental phosphorescence spectra of the β-cyclodextrin derivative. The GEBF approach is expected to be routinely applied to investigate the electronic emission spectra of very large systems with local chromophores.

Excited States Computation of Models of Phenylalanine Protein Chains: TD-DFT and Composite CC2/TD-DFT Protocols

The present benchmark calculations testify to the validity of time-dependent density functional theory (TD-DFT) when exploring the low-lying excited states potential energy surfaces of models of phenylalanine protein chains. Among three functionals suitable for systems exhibiting charge-transfer excited states, LC-ωPBE, CAM-B3LYP, and ωB97X-D, which were tested on a reference peptide system, we selected the ωB97X-D functional, which gave the best results compared to the approximate coupled-cluster singles and doubles (CC2) method. A quantitative agreement for both the geometrical parameters and the vibrational frequencies was obtained for the lowest singlet excited state (a ππ* state) of the series of capped peptides. In contrast, only a qualitative agreement was met for the corresponding adiabatic zero-point vibrational energy (ZPVE)-corrected excitation energies. Two composite protocols combining CC2 and DFT/TD-DFT methods were then developed to improve these calculations. Both protocols substantially reduced the error compared to CC2 and experiment, and the best of both even led to results of CC2 quality at a lower cost, thus providing a reliable alternative to this method for very large systems.

Multiphoton-excited dynamics of the trans or cis structural isomer of 1,2-dibromoethylene

Photofragmentation dynamics of cis and trans isomers of 1,2-dibromoethylene (1,2-DBE) have been investigated by multiphoton excitation using a picosecond (ps) laser pulse. It has been found that the Br₂ ⁺ product ion preferentially originates from the cis isomer rather than from trans. The Boltzmann-type isotropic low kinetic energy components of the Br⁺ and Br₂ ⁺ product state distributions seem to be most likely from the unimolecular reactions of the vibrationally hot cationic ground state generated by the three-photon absorption at the photon energy below ∼38 000 cm^-1. The highly anisotropic kinetic energy components of Br⁺ and Br₂ ⁺ start to appear at the photon energy above ∼38 000 cm^-1, where the D_n (n ≥ 1) - D₀ transition is facilitated within the same ps laser pulse as the parent molecule is efficiently ionized by the two-photon absorption. The transition dipole moment of the D₄ - D₀ transition of the strongest oscillator strength has been theoretically predicted to be parallel to the C-Br bond or C=C bond axis for the trans or cis isomer, respectively. The fast anisotropic with the (β ∼ +2) component in the Br⁺ product distribution is thus likely from the trans isomer, whereas that of Br₂ ⁺ (β ∼ -0.5) should be the consequence of the photodissociation of the cis isomer. The isomer-specific reactivity found here in the picosecond multiphoton excitation of 1,2-DBE provides a nice platform for the better understanding of the structure-reactivity relationship under the harsh condition of the strong or ultrashort optical field.

Model protein excited states: MRCI calculations with large active spaces vs CC2 method

Benchmarking calculations on excited states of models of phenylalanine protein chains are presented to assess the ability of alternative methods to the standard and most commonly used multiconfigurational wave function-based method, the complete active space self-consistent field (CASSCF), in recovering the non-dynamical correlation for systems that become not affordable by the CASSCF. The exploration of larger active spaces beyond the CASSCF limit is benchmarked through three strategies based on the reduction in the number of determinants: the restricted active space self-consistent field, the generalized active space self-consistent field (GASSCF), and the occupation-restricted multiple active space (ORMAS) schemes. The remaining dynamic correlation effects are then added by the complete active space second-order perturbation theory and by the multireference difference dedicated configuration interaction methods. In parallel, the approximate second-order coupled cluster (CC2), already proven to be successful for small building blocks of model proteins in one of our previous works [Ben Amor et al., J. Chem. Phys. 148, 184105 (2018)], is investigated to assess its performances for larger systems. Among the different alternative strategies to CASSCF, our results highlight the greatest efficiency of the GASSCF and ORMAS schemes in the systematic reduction of the configuration interaction expansion without loss of accuracy in both nature and excitation energies of both singlet ππ^* and nπ^* _CO excited states with respect to the equivalent CASSCF calculations. Guidelines for an optimum applicability of this scheme to systems requiring active spaces beyond the complete active space limit are then proposed. Finally, the extension of the CC2 method to such large systems without loss of accuracy is demonstrated, highlighting the great potential of this method to treat accurately excited states, mainly single reference, of very large systems.

Conformer-Specific Tunneling Dynamics Dictated by the Seam Coordinate of the Conical Intersection

The dynamic role of the conical intersection "seam" coordinate has been first revealed in the H fragmentation reaction of ortho(o)-cresol conformers. One of the (3N - 8) dimensional seam coordinates of the S₁(ππ*)/S₂(πσ*) conical intersection has been identified as the CH₃ torsional potential function. The tunneling dynamics of the reactive flux is dictated by its nuclear layout with respect to the CH₃ torsional angle, as the multidimensional tunneling barrier is dynamically shaped along the conical intersection seam. The effective tunneling-barrier weight-averaged over the quantum-mechanical probability along the CH₃ torsional angle perfectly explains the experimental finding: the sharp variation of the tunneling rate ((700-400) ps^-1) with the CH₃ torsional mode excitations within the narrow (0-100 cm^-1) energetic window. The much longer S₁ lifetime of cis compared to trans is ascribed to the higher-lying S₁/S₂ conical intersection of the former. With the use of distinct lifetimes, vibronic bands of each conformer could be completely separated.

Conformational analysis by UV spectroscopy: the decisive contribution of environment-induced electronic Stark effects

UV chromophores are frequently used as probes of the molecular structure. In particular, they are sensitive to the electric field generated by the molecular environment, resulting in the observation of Stark effects on UV spectra. While these environment-induced electronic Stark effects (EI-ESE) are already used for conformational analysis in the condensed phase, this work explores the potential of such an approach when performed at much higher conformational resolution in the gas phase. By investigating model alkali benzylacetate and 4-phenylbutyrate ion pairs, where the electric field applied to the phenyl ring is chemically tuned by changing the nature of the alkali cation, this work demonstrates that precise conformational assignments can be proposed based on the correlation between the conformation-dependent calculated electric fields and the frequency of the electronic transitions observed in the experimental UV spectra. Remarkably, the sole analysis of Stark effects and fragmentation patterns in mass-selected UV spectra provided an accurate and complete conformational analysis, where spectral differences as small as a few cm-1 between electronic transitions were rationalized. This case study illustrates that the identification of EI-ESE together with their interpretation at the modest cost of a ground state electric field calculation qualify UV spectroscopy as a powerful tool for conformational analysis.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

Cryogenic Ion Spectroscopy of a Singly Protonated Peptide DYYVVR: Locating Phosphorylation Sites of a Kinase Domain

Cryogenic ion spectroscopy (CIS) was applied to singly protonated DYYVVR, a tryptic peptide that contains the two active tyrosine residues (Y980 and Y981) of the Janus kinase 3 (JAK3) kinase domain, together with its point mutants (Y980F and Y981F) and phosphorylated peptides (pY980, pY981, and pY980pY981). The two tyrosine chromophores showed distinguishable UV absorption bands at around 35 200 and 35 450 cm^-1, respectively. By comparing with the point mutants, the lower electronic band was assigned to the absorption of Y981, and the higher one was assigned to Y980. When phosphorylated, the UV absorption of the phosphorylated chromophore shifts to higher energy above 36 500 cm^-1 but the unphosphorylated chromophore gives the absorption in the same region. Conformer-specific IR spectroscopy and density functional theory (DFT) calculations were used to tentatively assign the structure of DYYVVR. Two conformations were found, where Y981 is solvated by the protonated side chain of arginine R984, and the orientation of the carboxylic OH of D979 was different between the two. It is demonstrated that CIS can be used to distinguish the two tyrosine chromophores and to locate the phosphorylation site of a kinase domain.

Hydrogen-Bonded Complexes of Fluorophenylacetylenes: To Fluoresce or Not?

The hydrogen-bonded complexes of fluorophenylacetylenesexhibit unusual and interesting fluorescence turn ON/OFF behaviour following excitation to ¹ ππ* (S₁ ) state. The fluorescence switching behaviour can be realized by (i) "change in the intermolecular structure, (ii) change in the position of fluorine substitution and (iii) change in the hydrogen bonding partner or a combination thereof. Experiments indicate that the ≡C-H⋅⋅⋅X (X=O, N) hydrogen bonding with the acetylenic group plays a pivotal role in this switching behaviour. Intriguingly, weaker ≡C-H⋅⋅⋅X hydrogen bonding leads to fluorescence OFF state, which is turned ON by stronger hydrogen bonding. The observed fluorescence this switching behaviour is rationalized on the basis of a phenomenological model which suggests a coupling between the initially excited S₁ state and a dark S_n state in the Franck-Condon region with limited window controlled by the ≡C-H⋅⋅⋅X hydrogen bonding as a crucial parameter. Such fluorescence switching behaviour in hydrogen-bonded complexes is unprecedented and these intriguing results hopefully will stimulate theoreticians to test 'state of the art' theories to explain these observations in a consistent manner.

Infrared Spectroscopy as a Probe of Electronic Energy Transfer

We have combined electronic and vibrational spectroscopy in a cryogenic ion trap to produce highly resolved, conformer-selective spectra for the ground and excited states of a peptide containing two chromophores. These spectra permit us to determine the precise three-dimensional structure of the peptide and give insight into the migration of the electronic excitation from phenylalanine to tyrosine because changes in the excited-state infrared spectra are sensitive to localization of the electronic energy in each chromophore. The well-controlled experimental conditions make this result a stringent test for theoretical methods dealing with electronic energy transfer.

Stereochemistry-dependent structure of hydrogen-bonded protonated dimers: the case of 1-amino-2-indanol

Stereochemistry effects on the structure of molecular aggregates are studied in the prototypical 1-amino-2-indanol. Conformer-selective IR-UV double resonance spectroscopy reveals how stereochemistry shapes its dimers.

Nonradiative Relaxation Mechanisms of UV Excited Phenylalanine Residues: A Comparative Computational Study

The present work is directed toward understanding the mechanisms of excited state deactivation in three neutral model peptides containing the phenylalanine residue. The excited state dynamics of theγL(g+)folded form of N-acetylphenylalaninylamide (NAPA B) and its amide-N-methylated derivative (NAPMA B) is reviewed and compared to the dynamics of the monohydrated structure of NAPA (NAPAH). The goal is to unravel how the environment, and in particular solvation, impacts the photodynamics of peptides. The systems are investigated using reaction path calculations and surface hopping nonadiabatic dynamics based on the coupled cluster doubles (CC2) method and time-dependent density functional theory. The work emphasizes the role that excitation transfer from the phenylππ*to amidenπ*state plays in the deactivation of the three systems and shows how the ease of out-of-plane distortions of the amide group determines the rate of population transfer between the two electronic states. The subsequent dynamics on thenπ*state is barrierless along several pathways and leads to fast deactivation to the ground electronic state.

A computational study on how structure influences the optical properties in model crystal structures of amyloid fibrils

Amyloid fibrils have been shown to have peculiar optical properties since they can exhibit fluorescence in the absence of aromatic residues. In a recent study, we have shown that proton transfer (PT) events along hydrogen bonds (HBs) are coupled to absorption in the near UV range. Here, we gain more insights into the different types of hydrogen bonding interactions that occur in our model systems and the molecular factors that control the susceptibility of the protons to undergo PT and how this couples to the optical properties. In the case of the strong N-C termini interactions, a nearby methionine residue stabilizes the non-zwitterionic NH₂-COOH pair, while zwitterionic NH₃⁺-COO- is stabilized by the proximity of nearby crystallographic water molecules. Proton motion along the hydrogen bonds in the fibril is intimately coupled to the compression of the heavier atoms, similar to what is observed in bulk water. Small changes in the compression of the hydrogen bonds in the protein can lead to significant changes in both the ground and excited state potential energy surfaces associated with PT. Finally, we also reinforce the importance of nuclear quantum fluctuations of protons in the HBs of the amyloid proteins.

Mechanism of ultrafast non-reactive deactivation of the retinal chromophore in non-polar solvents

Counterion sensitive photodynamics of the retinal chromophore in solution.

Nonstatistical UV Fragmentation of Gas-Phase Peptides Reveals Conformers and Their Structural Features

Solving the 3D structure of a biomolecule requires recognition of its conformers and measurements of their individual structural identities, which can be compared with calculations. We employ the phenomenon of nonstatistical photofragmentation, detected by a combination of UV cold ion spectroscopy and high-resolution mass spectrometry, to identify the main conformers of gas-phase peptides and to recover individual UV absorption and mass spectra of all of these conformers in a single laser scan. We first validate this approach with a benchmark dipeptide, Tyr-Ala, and then apply it to a decapeptide, gramicidin S. The revealed characteristic structural difference between the conformers of the latter identifies some of the previously calculated structures of gramicidin S as the most likely geometries of its remaining unsolved conformer.

How does tetraphenylethylene relax from its excited states?

Photocyclization play a key role in the deactivation mechanism of tetraphenylethylene.

Isolated neutral peptides

This chapter examines the structural characterisation of isolated neutral amino-acids and peptides. After a presentation of the experimental and theoretical state-of-the-art in the field, a review of the major structures and shaping interactions is presented. Special focus is made on conformationally-resolved studies which enable one to go beyond simple structural characterisation; probing flexibility and excited-state photophysics are given as examples of promising future directions.

Photo-stability of peptide-bond aggregates: N-methylformamide dimers

The formation of weakly-bound dimers of N-methylformamide (NMF) and the photochemistry of these dimers after irradiation at 248 nm were explored using matrix-isolation spectroscopy. Calculations were used to characterize the diverse isomers and assign their IR spectra; non-adiabatic dynamics was simulated to understand their photo-deactivation mechanism. The most stable dimers, and , were obtained by trans-trans aggregation (N-HO[double bond, length as m-dash]C interactions) and could be identified in the matrix. The main products formed after irradiation are the trans-cis dimers ( and ), also stabilized by N-HO[double bond, length as m-dash]C interactions. In contrast to the photochemistry of the monomers, no dissociative products were observed after 248 nm irradiation of the dimers. The absence of dissociative products can be explained by a proton-transfer mechanism in the excited state that is faster than the photo-dissociative mechanism. The fact that hydrogen bonding has such a significant effect on the photochemical stability of NMF has important implications to understand the stability of peptide-bonded systems to UV irradiation.

On the photostability of peptides after selective photoexcitation of the backbone: prompt versus slow dissociation

Vulnerability of biomolecules to ultraviolet radiation is intimately linked to deexcitation pathways: photostability requires fast internal conversion to the electronic ground state, but also intramolecular vibrational redistribution and cooling on a time scale faster than dissociation. Here we present a protocol to disentangle slow and non-hazardous statistical dissociation from prompt cleavage of peptide bonds by 210 nm light based on experiments on protonated peptides isolated in vacuo and tagged by 18-crown-6 ether (CE). The weakest link in the system is between the charged site and CE, which is remote from the initial site of excitation. Hence loss of CE serves as direct proof that energy has reached the charge-site end, leaving the backbone intact. Our work demonstrates that excitation of tertiary amide moieties (proline linkages) results in both prompt dissociation and statistical dissociation after energy randomisation over all vibrational degrees of freedom.

Solvent effects on the photochemistry of 4-aminoimidazole-5-carbonitrile, a prebiotically plausible precursor of purines

4-Aminoimidazole-5-carbonitrile (AICN) was suggested as a prebiotically plausible precursor of purine nucleobases and nucleotides. Although it can be formed in a sequence of photoreactions, AICN is immune to further irradiation with UV-light. We present state-of-the-art multi-reference quantum-chemical calculations of potential energy surface cuts and conical intersection optimizations to explain the molecular mechanisms underlying the photostability of this compound. We have identified the N-H bond stretching and ring-puckering mechanisms that should be responsible for the photochemistry of AICN in the gas phase. We have further considered the photochemistry of AICN-water clusters, while including up to six explicit water molecules. The calculations reveal charge transfer to solvent followed by formation of an H3O(+) cation, both of which occur on the (1)πσ* hypersurface. Interestingly, a second proton transfer to an adjacent water molecule leads to a (1)πσ*/S0 conical intersection. We suggest that this electron-driven proton relay might be characteristic of low-lying (1)πσ* states in chromophore-water clusters. Owing to its nature, this mechanism might also be responsible for the photostability of analogous organic molecules in bulk water.

Photo-induced reactions from efficient molecular dynamics with electronic transitions using the FIREBALL local-orbital density functional theory formalism

The computational simulation of photo-induced processes in large molecular systems is a very challenging problem. Firstly, to properly simulate photo-induced reactions the potential energy surfaces corresponding to excited states must be appropriately accessed; secondly, understanding the mechanisms of these processes requires the exploration of complex configurational spaces and the localization of conical intersections; finally, photo-induced reactions are probability events, that require the simulation of hundreds of trajectories to obtain the statistical information for the analysis of the reaction profiles. Here, we present a detailed description of our implementation of a molecular dynamics with electronic transitions algorithm within the local-orbital density functional theory code FIREBALL, suitable for the computational study of these problems. As an example of the application of this approach, we also report results on the [2 + 2] cycloaddition of ethylene with maleic anhydride and on the [2 + 2] photo-induced polymerization reaction of two C60 molecules. We identify different deactivation channels of the initial electron excitation, depending on the time of the electronic transition from LUMO to HOMO, and the character of the HOMO after the transition.

Description of Conical Intersections with Density Functional Methods

Conical intersections are perhaps the most significant mechanistic features of chemical reactions occurring through excited states. By providing funnels for efficient non-adiabatic population transfer, conical intersections govern the branching ratio of products of such reactions, similar to what the transition states do for ground-state reactivity. In this regard, intersections between the ground and the lowest excited states play a special role, and the correct description of the potential energy surfaces in their vicinity is crucial for understanding the mechanism and dynamics of excited-state reactions. The methods of density functional theory, such as time-dependent density functional theory, are widely used to describe the excited states of large molecules. However, are these methods suitable for describing the conical intersections or do they lead to artifacts and, consequently, to erroneous description of reaction dynamics? Here we address the first part of this question and analyze the ability of several density functional approaches, including the linear-response time-dependent approach as well as the spin-flip and ensemble formalisms, to provide the correct description of conical intersections and the potential energy surfaces in their vicinity. It is demonstrated that the commonly used linear-response time-dependent theory does not yield a proper description of these features and that one should instead use alternative computational approaches.

Timescales of N–H bond dissociation in pyrrole: a nonadiabatic dynamics study

The excitation wavelength dependent photodynamics of pyrrole are investigated by surface-hopping nonadiabatic dynamics simulations. The results are explained in terms of correct Rydberg–valence interaction in the lowest B₂ state.

A comprehensive study of cold protonated tyramine: UV photodissociation experiments and ab initio calculations

Excited state properties of cold protonated ions are revealed by a combination of laser spectroscopy and ab initio calculations.