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.

An intraresidue H-bonding motif in selenocysteine and cysteine, revealed by gas phase laser spectroscopy and quantum chemistry calculations.

Models of protein chains containing a seleno-cysteine (Sec) residue have been investigated by gas phase laser spectroscopy in order to document the effect of the H-bonding properties of the SeH group in the folding of the Sec side chain, by comparison with recent data on Ser- and Cys-containing sequences. Experimental data, complemented by quantum chemistry calculations and natural bonding orbital (NBO) analyses, are interpreted in terms of the formation of a so-called 5γ intra-residue motif, which bridges the acceptor chalcogen atom of the side chain to the NH bond of the same residue. This local structure, in which the O/S/Se atom is close to the plane of the N-terminal side amide, is constrained by local backbone-side chain hyperconjugation effects involving the S and Se atoms. Theoretical investigations of the Cys/Sec side chain show that (i) this 5γ motif is an intrinsic feature of these residues, (ii) the corresponding H-bond is strongly non-linear and intrinsically weak, (iii) but enhanced by γ- and ß-turn secondary structures, which promote a more favorable 5γ H-bonding approach and distance. The resulting H-bonds are slightly stronger in selenocysteine than in cysteine, but nearly inexistent in serine, whose side chain in contrast behaves as a H-bonding donor. The modest spectral shifts of the Cys/Sec NH stretches measured experimentally reflect the moderate strength of the 5γ H-bonding, in agreement with the correlation obtained with a NBO-based H-bond strength indicator. The evolution along the Ser, Cys and Sec series emphasizes the compromise between the several factors that control the H-bonding in a hyperconjugation-constrained geometry, among them the chalcogen van der Waals and covalent radii. It also illustrates the 5γ H-bond enhancements with the Sec and Cys residues favoured by the constraints imposed by the γ- and ß-turn structures of the peptide chain.

CC2 Benchmark for Models of Phenylalanine Protein Chains: 0-0 Transition Energies and IR Signatures of the ππ* Excited State.

Extensive benchmarking calculations are presented to assess the accuracy of the standard approximate coupled cluster singles and doubles method (CC2) in studying ππ* excited states properties of model protein chains containing a phenylalanine residue, namely capped peptides, whose ground state conformers adopt the prototypical secondary structural features of proteins. First, the dependence with the basis set of the CC2 excitation energies, CC2 geometry optimizations, and amide A region frequencies of the lowest ππ* excited state in a reference system, the N-acetylphenylalaninylamide, are investigated, and the results are compared with experimental data. Second, at the best level of theory determined, the CC2/aug(N,O,π)-cc-pVDZ//CC2/cc-pVDZ level, a series of capped peptides of increasing size and containing residues of different nature are investigated. Along the series, compared to the experimental values, a mean absolute error of 0.10 eV is achieved for the 0-0 transition energies with a systematic overestimation. In addition, mode-dependent linear scaling functions for the calculated frequencies of the amide A region have been determined from the set of 95 experimental frequencies available; they lead to a quantitative simulation of the observed shifts of the amide A region frequencies upon ππ* excitation (root-mean-square deviation of 5 cm-1). These results confirm the reliability of the CC2 method to characterize the lowest ππ* excited state of such medium-sized systems, emphasizing this class of theoretical approaches as a relevant spectroscopic tool, including for tasks as difficult as conformational assignment.

Identification of ion pairs in solution by IR spectroscopy: crucial contributions of gas phase data and simulations.

In a context where structure elucidation of ion pairs in solution remains a contemporary challenge, this work explores an original approach where accurate gas phase spectroscopic data are used to refine high level quantum chemistry calculations of ion pairs in solution, resulting in an unprecedented level of accuracy in vibrational frequency prediction. First, gas phase studies focus on a series of isolated contact ion pairs (M+, Ph-CH2-COO-, with M = Li, Na, K, Rb, Cs) for which conformer-selective IR spectra in the CO2- stretch region are recorded. These experiments reveal the interactions at play in isolated contact ion pairs, and provide vibrational frequencies enabling us to assess the accuracy of the theoretical approach used, i.e., mode-dependent scaled harmonic frequency calculations at the RI-B97-D3/dhf-TZVPP level. This level of calculation is then employed on large water clusters embedding either a free acetate ion or its contact or solvent-shared pairs with a sodium cation in order to simulate the individual vibrational spectra of these species in solution. This study shows that the stretching modes of carboxylate are sensitive to both solvent-shared and contact ion pair formation. FTIR spectra of solutions of increasing concentrations indeed reveal several spectral changes consistent with the presence of specific types of solvent-shared and contact ion pairs. By providing relevant guidelines for the interpretation of solution phase IR spectra, this work illustrates the potential of the approach for the elucidation of supramolecular structures in electrolyte solutions.

Intrinsic Folding Proclivities in Cyclic ß-Peptide Building Blocks: Configuration and Heteroatom Effects Analyzed by Conformer-Selective Spectroscopy and Quantum Chemistry.

This work describes the use of conformer-selective laser spectroscopy following supersonic expansion to probe the local folding proclivities of four-membered ring cyclic ß-amino acid building blocks. Emphasis is placed on stereochemical effects as well as on the structural changes induced by the replacement of a carbon atom of the cycle by a nitrogen atom. The amide A IR spectra are obtained and interpreted with the help of quantum chemistry structure calculations. Results provide evidence that the building block with a trans-substituted cyclobutane ring has a predilection to form strong C8 hydrogen bonds. Nitrogen-atom substitution in the ring induces the formation of the hydrazino turn, with a related but distinct hydrogen-bonding network: the structure is best viewed as a bifurcated C8/C5 bond with the N heteroatom lone electron pair playing a significant acceptor role, which supports recent observations on the hydrazino turn structure in solution. Surprisingly, this study shows that the cis-substituted cyclobutane ring derivative also gives rise predominantly to a C8 hydrogen bond, although weaker than in the two former cases, a feature that is not often encountered for this building block.

Secondary Structures in Phe-Containing Isolated Dipeptide Chains: Laser Spectroscopy vs Quantum Chemistry.

The intrinsic conformational landscape of two phenylalanine-containing protein chain models (-Gly-Phe- and -Ala-Phe- sequences) has been investigated theoretically and experimentally in the gas phase. The near UV spectroscopy (first ππ* transition of the Phe ring) is obtained experimentally under jet conditions where the conformational features can be resolved. Single-conformation IR spectroscopy in the NH stretch region is then obtained by IR/UV double resonance in the ground state, leading to resolved vibrational spectra that are assigned in terms of conformation and H-bonding content from comparison with quantum chemistry calculations. For the main conformer, whose UV spectrum exhibits a significant Franck-Condon activity in low frequency modes involving peptide backbone motions relative to the Phe chromophore, excited state IR spectroscopy has also been recorded in a UV/IR/UV experiment. The NH stretch spectral changes observed in such a ππ* labeling experiment enable us to determine those NH bonds that are coupled to the phenyl ring; they are compared to CC2 excited state calculations to quantify the geometry change upon ππ* excitation. The complete and consistent series of data obtained enable us to propose an unambiguous assignment for the gallery of conformers observed and to demonstrate that, in these two sequences, three conceptually important local structural motifs of proteins (ß-strands, 27 ribbons, and ß-turns) are represented. The satisfactory agreement between the experimental conformational distribution and the predicted landscape anticipated from the DFT-D approach demonstrates the capabilities of a theoretical method that accounts for dispersive interactions. It also shows that the flaws, inherent to a resonant two-photon ionization detection scheme, often evoked for aromatic chromophores, do not seem to be significant in the case of Phe.

Direct spectroscopic evidence of hyperconjugation unveils the conformational landscape of hydrazides.

The stereochemistry of hydrazides makes them especially interesting as building blocks for molecular design. An exhaustive conformational analysis of three model hydrazides was conducted in a conformer-selective approach by using a combination of high-level quantum chemistry calculations and vibrational spectroscopy in the gas phase and in solution. The NH stretch frequency was found to be highly sensitive to hyperconjugation, thus making it an efficient probe of the conformation of the neighboring nitrogen atom. This property greatly assisted the identification of the isomers observed experimentally in the conformer pool. A rationalization of the hydrazide conformational landscape is proposed, therefore paving the way for a better characterization of secondary structures in larger systems.

On the near UV photophysics of a phenylalanine residue: conformation-dependent ππ* state deactivation revealed by laser spectroscopy of isolated neutral dipeptides.

The primary step of the near UV photophysics of a phenylalanine residue is investigated in one- and two-color pump-probe R2PI nanosecond experiments carried out on specific conformers of the Ac-Gly-Phe-NH2 molecule and related neutral compounds isolated in a supersonic expansion. Compared to toluene, whose ππ* state photophysics is dominated by intersystem crossing with a lifetime of ∼80 ns at the origin, the first ππ* state of Phe in the peptide environment is systematically found to be shorter-lived. The lifetime at the origin of transition is found to be significantly shortened in the presence of a primary amide (-CONH2) group (20-60 ns, depending on the conformer considered), demonstrating the existence of an additional non-radiative relaxation channel related to this chemical group. The quenching effect induced by the peptide environment is still more remarkable beyond the origin of the ππ* state, since vibronic bands of one of the 4 conformers observed (the 27-ribbon conformation) become barely detectable in the ns R2PI experiment, suggesting a significant conformer-selective lifetime shortening (below 100 ps). These results on dipeptides, which extend previous investigations on shorter Phe-containing molecules (N-Ac-Phe-NH2 and N-Ac-Phe-NH-Me), confirm the existence of conformer-dependent non-radiative deactivation processes, whose characteristic timescales range from tens of ns down to hundreds of ps or below. This dynamics is assigned to two distinct mechanisms: a first one, consistent with an excitation energy transfer from the optically active ππ* state to low-lying amide nπ* excited states accessed through conical intersections, especially in the presence of a C-terminal primary amide group (-CONH2); a second one, responsible for the short lifetimes in 2(7) ribbon structures, would be more specifically triggered by phenyl ring vibrational excitations. Implications in terms of spectroscopic probing of Phe in a peptide environment, especially in the presence of a quenching amide group, are discussed.

Unraveling the mechanisms of nonradiative deactivation in model peptides following photoexcitation of a phenylalanine residue.

The mechanisms of nonradiative deactivation of a phenylalanine residue after near-UV photoexcitation have been investigated in an isolated peptide chain model (N-acetylphenylalaninylamide, NAPA) both experimentally and theoretically. Lifetime measurements at the origin of the first ππ* state of jet-cooled NAPA molecules have shown that (i) among the three most stable conformers of the molecule, the folded conformer NAPA B is ∼50-times shorter lived than the extended major conformer NAPA A and (ii) this lifetime is virtually insensitive to deuteration at the NH(2) and NH sites. Concurrent time-dependent density functional theory (TDDFT) based nonadiabatic dynamics simulations in the full dimensionality, carried out for the NAPA B conformer, provided direct insights on novel classes of ultrafast deactivation mechanisms, proceeding through several conical intersections and leading in fine to the ground state. These mechanisms are found to be triggered either (i) by a stretch of the N(Phe)H bond, which leads to an H-transfer to the ring, or (ii) by specific backbone amide distortions. The potential energy surfaces of the NAPA conformers along these critical pathways have been characterized more accurately using the coupled cluster doubles (CC2) method and shown to exhibit barriers that can be overcome with moderate excess energies. These results analyzed in the light of the experimental findings enabled us to assign the short lifetime of NAPA B conformer to a number of easily accessible exit channels from the initial ππ* surface, most importantly the one involving a transfer of electronic excitation to an nπ* surface, induced by distortions of the backbone peptide bond.

Strength of NH···S Hydrogen Bonds in Methionine Residues Revealed by Gas-Phase IR/UV Spectroscopy.

Despite of being ubiquitous in proteins, NHbackbone···S hydrogen bonds linking the sulfur atom of methionine or cysteine to backbone NH groups remain poorly documented. Here, we report vibrationally resolved IR NH stretch spectra of two methionine-containing dipeptides (Ac-Phe-Met-NH2 and Ac-Met-Phe-NH2). The conformations observed for both molecules, assigned with the help of DFT-D quantum chemistry, provide spectroscopic evidence for the formation of NHbackbone···S H-bonds, surprisingly strong enough to challenge the classical intrabackbone NH···O═C H-bonds. The methionine side chain is found to fold locally, forming a H-bond with the neighboring amide groups (NH(i) or NH(i+1)). Comparison with protein data bank structural information shows that such a local folding is also common in proteins where it concerns 24% of the methionine residues that have a sulfur atom linked to a backbone NH group. This convergence between the strength of these NH···S H-bonds and protein structural data illustrates their contribution to the stability of protein chains.

Isolated monohydrates of a model peptide chain: effect of a first water molecule on the secondary structure of a capped phenylalanine.

The formation of monohydrates of capped phenylalanine model peptides, CH(3)-CO-Phe-NH(2) and CH(3)-CO-Phe-NH-CH(3), in a supersonic expansion has been investigated using laser spectroscopy and quantum chemistry methods. Conformational distributions of the monohydrates have been revealed by IR/UV double-resonance spectroscopy and their structures assigned by comparison with DFT-D calculations. A careful analysis of the final hydrate distribution together with a detailed theoretical investigation of the potential energy surface of the monohydrates demonstrates that solvation occurs from the conformational distribution of the isolated peptide monomers. The distribution of the monohydrates appears to be strongly dependent on both the initial monomer conformation (extended or folded backbone) and the solvation site initially occupied by the water molecule. The solvation processes taking place during the cooling can be categorized as follows: (a) solvation without significant structural changes of the peptide, (b) solvation inducing significant distortions of the backbone but retaining the secondary structure, and (c) solvation triggering backbone isomerizations, leading to a modification of the peptide secondary structure. It is observed that solvation by a single water molecule can fold a ß-strand into a γ-turn structure (type c) or induce a significant opening of a γ-turn characterized by an elongated C(7) hydrogen bond (type b). These structural changes can be considered as a first step toward the polyproline II condensed-phase structure, illustrating the role played by the very first water molecule in the solvation process.

Experimental and theoretical investigation of the aromatic-aromatic interaction in isolated capped dipeptides.

Among the forces responsible for shaping proteins, interactions between side chains of aromatic residues play an important role as they are involved in the secondary and the tertiary structures of proteins contributing to the formation of hydrophobic domains. The purpose of this paper is to document this interaction in two capped dipeptides modeling a segment of a protein chain having two consecutive Phe residues, Ac-Phe-Phe-NH(2) and Ac-Phe-D-Phe-NH(2). These two molecules have been investigated in the gas phase by IR/UV double resonance spectroscopy, and the assignment of the observed conformers has been done by comparison with quantum chemistry calculations. Both peptides are found to adopt a beta-turn type I conformation stabilized by an edge-to-face interaction between the two aromatic rings. Comparison with other dipeptides in the literature demonstrates the impact of this aromatic-aromatic interaction on the shape adopted by the peptide chain, and its role among the other shaping forces (H-bonds, NH-pi interactions) is discussed. As an illustration, the H-bond strength is found to be significantly lower in the beta-turn type I conformer, in which the two rings interact, as compared to the similar conformer where such an interaction does not exist. This structural feature due to the backbone distortion induced by the interaction between the aromatic rings makes this system a good test for evaluating the ability of computational methods to correctly account for the competition between these forces. MP2, SCS-MP2, DFT, and DFT-D methods have been assessed in this respect. Comparison between geometries, energies, and frequency calculations illustrate their respective limitations in describing conformations resulting from a subtle equilibrium between the several interactions at play.

Gas phase folding of an (Ala)4 neutral peptide chain: spectroscopic evidence for the formation of a beta-hairpin H-bonding pattern.

IR and UV laser spectroscopy of an Ala-based 4-residue model peptide recorded under gas phase isolated conditions provides evidence for the intrinsic stability of compact folded structures resembling the extremity of a beta-hairpin, with a C(14) H-bond bridging the two ends of the chain, and enables us to assess the capabilities of new quantum chemistry techniques to account for dispersive interactions in a medium-size molecule.

Spectroscopic evidence for the formation of helical structures in gas-phase short peptide chains.

Aminoisobutyric acid (Aib) is a synthetic amino acid known to favor the formation of 3(10) helical structures in condensed phases, namely, crystals. The intrinsic character of these helicogenic properties has been investigated on the Ac-Aib-Phe-Aib-NH2 molecule under isolated conditions, namely, in the gas phase, both experimentally by double-resonance IR/UV spectroscopy and theoretically by quantum chemistry. A convergent set of evidence, based on energetic, IR, and UV spectroscopic data as well as on analogies with the similar peptide Ac-Ala-Phe-Ala-NH2 previously studied, enables us to conclude the formation of an incipient 310 helix in these isolated systems.

Time-resolved photoelectron and photoion fragmentation spectroscopy study of 9-methyladenine and its hydrates: a contribution to the understanding of the ultrafast radiationless decay of excited DNA bases.

The excited state dynamics of the purine base 9-methyladenine (9Me-Ade) has been investigated by time- and energy-resolved photoelectron imaging spectroscopy and mass-selected ion spectroscopy, in both vacuum and water-cluster environments. The specific probe processes used, namely a careful monitoring of time-resolved photoelectron energy distributions and of photoion fragmentation, together with the excellent temporal resolution achieved, enable us to derive additional information on the nature of the excited states (pipi*, npi*, pisigma*, triplet) involved in the electronic relaxation of adenine. The two-step pathway we propose to account for the double exponential decay observed agrees well with recent theoretical calculations. The near-UV photophysics of 9Me-Ade is dominated by the direct excitation of the pipi* ((1)L(b)) state (lifetime of 100 fs), followed by internal conversion to the npi* state (lifetime in the ps range) via conical intersection. No evidence for the involvement of a pisigma* or a triplet state was found. 9Me-Ade-(H(2)O)(n) clusters have been studied, focusing on the fragmentation of these species after the probe process. A careful analysis of the fragments allowed us to provide evidence for a double exponential decay profile for the hydrates. The very weak second component observed, however, led us to conclude that the photophysics were very different compared with the isolated base, assigned to a competition between (i) a direct one-step decay of the initially excited state (pipi* L(a) and/or L(b), stabilised by hydration) to the ground state and (ii) a modified two-step decay scheme, qualitatively comparable to that occurring in the isolated molecule.

Gas phase formation of a 3(10)-helix in a three-residue peptide chain: role of side chain-backbone interactions as evidenced by IR-UV double resonance experiments.

The first spectroscopic evidence for the gas-phase formation of helical structures in short peptide chains is reported, using the IR-UV double resonance technique and DFT quantum chemistry calculations. The study involves three chemically protected peptides, all based on the same Ac-(Ala)3-NH2, (Ac = acetyl, Ala = alanine) tripeptide, in which one of the Ala residues is substituted by the aromatic phenylalanine residue. For the three molecules, only one main conformer is observed in the supersonic expansion. IR analysis shows that the structure of this conformer is strongly dependent upon the substitution site: the helical 310-type structure is observed only when Phe occupies the central residue of the chain. The present work also emphasizes that the 310-helix formation does compete with other archetypal H-bonding patterns, such as 27-ribbon or mixed structures, whose relative energetics can be greatly influenced by a modest NH-aromatic interaction.

Excited states dynamics of DNA and RNA bases: characterization of a stepwise deactivation pathway in the gas phase.

Radiationless deactivation pathways of excited gas phase nucleobases were investigated using mass-selected femtosecond resolved pump-probe resonant ionization. By comparison between nucleobases and methylated species, in which tautomerism cannot occur, we can access intrinsic mechanisms at a time resolution never reported so far (80 fs). At this time resolution, and using appropriate substitution, real nuclear motion corresponding to active vibrational modes along deactivation coordinates can actually be probed. We provide evidence for the existence of a two-step decay mechanism, following a 267 nm excitation of the nucleobases. The time resolution achieved together with a careful zero time-delay calibration between lasers allow us to show that the first step does correspond to intrinsic dynamics rather than to a laser cross correlation. For adenine and 9-methyladenine a first decay component of about 100 fs has been measured. This first step is radically increased to 200 fs when the amino group hydrogen atoms of adenine are substituted by methyl groups. Our results could be rationalized according to the effect of the highly localized nature of the excitation combined to the presence of efficient deactivation pathway along both pyrimidine ring and amino group out-of-plane vibrational modes. These nuclear motions play a key role in the vibronic coupling between the initially excited pipi* and the dark npi* states. This seems to be the common mechanism that opens up the earlier phase of the internal conversion pathway which then, in consideration of the rather fast relaxation times observed, would probably proceed via conical intersection between the npi* relay state and high vibrational levels of the ground state.

Intrinsic folding of small peptide chains: spectroscopic evidence for the formation of beta-turns in the gas phase.

Laser desorption of model peptides coupled to laser spectroscopic techniques enables the gas-phase observation of genuine secondary structures of biology. Spectroscopic evidence for the formation of beta-turns in gas-phase peptide chains containing glycine and phenylalanine residues establishes the intrinsic stability of these forms and their ability to compete with other stable structures. The precise characterization of local minima on the potential energy surface from IR spectroscopy constitutes an acute assessment for the state-of-the-art quantum mechanical calculations also presented. The observation of different types of beta-turns depending upon the residue order within the sequence is found to be consistent with the residue propensities in beta-turns of proteins, which suggests that the prevalence of glycine in type II and II' turns stems essentially from an energetic origin, already at play under isolated conditions.

The gas-phase dipeptide analogue acetyl-phenylalanyl-amide: a model for the study of side chain/backbone interactions in proteins.

The issue of the influence of the side chain/backbone interaction on the local conformational preferences of a phenylalanine residue in a peptide chain is addressed. A synergetic approach is used, which combines gas-phase UV spectroscopy as well as gas-phase IR/UV double-resonance experiments with DFT and post Hartree-Fock calculations. N-Acetyl-Phe-amide was chosen as a model system for which three different conformers were observed. The most stable conformer has been identified as an extended beta(L) conformation of the peptide backbone. It is stabilized by a weak but significant NH-pi interaction bridging the aromatic ring on the residue (i) with the NH group on residue (i+1), with the aromatic side chain being in an anti conformation. This stable conformation corresponds to the common NH(i+1)-aromatic(i) interaction encountered in proteins for the three aromatic residues (phenylalanine, tyrosine, and tryptophan), which illustrates the relevance of gas-phase investigations to structural biology issues. The two other less abundant conformers have been assigned to two gamma-folded backbone conformations that differ by the orientation of the side chain. In all cases, the IR data provided spectroscopic fingerprints of these interactions. Finally, the strong conformational dependence of the fluorescence yield found for N-acetyl-Phe-amide illustrates the role of the environment on the excited-state dynamics of these species, which is often exploited by biochemists to monitor protein structural changes from tryptophan lifetime measurements.