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.

High-resolution UV spectroscopy of 1-indanol

We report on rotationally resolved laser induced fluorescence (LIF) and vibrationally resolved resonance-enhanced multiphoton ionization (REMPI) spectroscopy of the chiral molecule 1-indanol. Spectra of the S₁← S₀ electronic transition are recorded in a jet-cooled, pulsed molecular beam. Using two time-delayed pulsed lasers, the lifetimes of the S₁ state of the two most stable conformers, referred to as eq1 and ax2, have been determined. The S₁← S₀ origin bands of these conformers as well as the transition to a vibrationally excited level in the S₁ state of eq1 are recorded with full rotational resolution (25 MHz observed linewidth) by measuring the LIF intensity following excitation with a tuneable, narrowband cw laser. On selected rotationally resolved electronic transitions, Lamb-dips are measured to confirm the Lorentzian lifetime-contribution to the observed lineshapes. The rotationally resolved S₁← S₀ origin band of a neon-complex of eq1 is measured via LIF as well. The fit of the rotationally resolved LIF spectra of the origin bands to those of an asymmetric rotor yields a standard deviation of about 6 MHz. The resulting spectroscopic parameters are tabulated and compared to the outcome of ab initio calculations. For both conformers as well as for the Ne-eq1 complex, the geometric structures in the S₀ and S₁ states are discussed. For all systems, the transition dipole moment is mainly along the a-axis, the contributions along the b- and c-axes being about one order of magnitude smaller.

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.

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.

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.

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.

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.

Local NH–π interactions involving aromatic residues of proteins: influence of backbone conformation and ππ* excitation on the π H-bond strength, as revealed from studies of isolated model peptides

Gas phase conformer-selective IR spectroscopy combined and relevant quantum chemistry methods document the NH–π interactions in Phe residues.