Potential-Energy and Free-Energy Surfaces of Glycyl-Phenylalanyl-Alanine (GFA) Tripeptide: Experiment and Theory

Structural Properties of Phenylalanine-Based Dimers Revealed Using IR Action Spectroscopy

Peptide segments with phenylalanine residues are commonly found in proteins that are related to neurodegenerative diseases. However, the self-assembly of phenylalanine-based peptides can be also functional. Peptides containing phenylalanine residues with different side caps, composition, and chemical alteration can form different types of nanostructures that find many applications in technology and medicine. Various studies have been performed in order to explain the remarkable stability of the resulting nanostructures. Here, we study the early stages of self-assembly of two phenylalanine derived peptides in the gas phase using IR action spectroscopy. Our focus lies on the identification of the key intra- and intermolecular interactions that govern the formation of the dimers. The far-IR region allowed us to distinguish between structural families and to assign the 2-(2-amino-2-phenylacetamido)-2-phenylacetic acid (PhgPhg) dimer to a very symmetric structure with two intermolecular hydrogen bonds and its aromatic rings folded away from the backbone. By comparison with the phenylalanine-based peptide cyclic L-phenylalanyl-L-phenylalanine (cyclo-FF), we found that the linear FF dimer likely adopts a less ordered structure. However, when one more phenylalanine residue is added (FFF), a more structurally organized dimer is formed with several intermolecular hydrogen bonds.

Determination of vibrational band positions in the E-hook of β-tubulin

Raman spectral characterization of the β-TUBB2A E-hook hexapeptide, EGEDEA, is determined through experimental analysis combined with full geometry optimizations and corresponding harmonic vibrational frequency computations employing DFT methods. The hexapeptide is first broken down into di- and tetrapeptide fragments which are analyzed both quantum chemically and experimentally, and then combined to achieve an energetic minimum of the large EGEDEA hexapeptide. The Raman spectral characterization of EGEDEA band positions are then verified via the literature and comparison to the small fragment's similarly located band positions. The approach employed provides further evidence for the use of fragments as a helpful tool in characterization of the vibrational band positions of large peptides. STATEMENT OF SIGNIFICANCE: To investigate β-TUBB2A E-hook hexapeptide, a unique approach is employed whereby the hexapeptide is broken into fragments, EG, ED, EA, EGED, and EDEA and analyzed via experimental Raman spectroscopy of the crystalline solids. The experimentally observed vibrational band positions are compared to those computed using and scaled from DFT methods and Pople's 6-311+G(2df,2pd) basis set. The reported vibrational band positions are also confirmed by previously reported bands of similar peptides in the literature. This methodology facilitates differentiation between the behaviors of various side chains and their influence on the structure of the hexapeptide, providing insight into not only the nature of the peptide but also defining regions for potential protein and cytoplasmic interactions, without requiring excessive computing resources or overly-sensitive experimental methods.

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.

Site-specific binding of a water molecule to the sulfa drugs sulfamethoxazole and sulfisoxazole: a laser-desorption isomer-specific UV and IR study

Using isomer-specific IR spectroscopy, we show that sulfamethoxazole and sulfisoxazole exhibit distinct site specificities for binding a water molecule.

Computational study on single molecular spectroscopy of tyrosin-glycine, tryptophane-glycine and glycine-tryptophane

Quantum chemistry calculations play a fundamental role in revealing the molecular structures observed in gas-phase spectroscopic measurements. The supersonic jet cooling widely used in single molecular spectroscopy experiment is a non-equilibrium process and often causes confusion on the theoretical and experimental comparison. A computational approach is proposed here to account for the effect of the non-equilibrium cooling on the experimental spectra and applied to the cases of tyrosin-glycine (YG), tryptophane-glycine (WG) and glycine-tryptophane (GW). The low energy conformers of YG, WG and GW are obtained through thorough conformational searches. The structural features and equilibrium distributions of conformations and the energy barriers for conformer conversions are then determined. Three classes of transition energy barriers, high, medium and low, are found for the conversions among conformers with distinctly different, similar and the same structural types, respectively. The final conformation populations are determined by assuming an initial temperature of about 450 K and allowing for only the conformation conversion with a low energy barrier to occur during the rapid cooling process. The results provide a natural explanation for the numbers of YG, WG and GW conformations observed experimentally. The theoretical conformation assignments are also in good agreement with the experimental IR data.

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.

Exploring short intramolecular interactions in alkylaromatic substrates

From proteins and peptides to semiconducting polymers, aliphatic chains on aromatic groups are recurring motifs in macromolecules from very diverse application fields. Fields in which molecular folding and packing determine the macroscopic physical properties that make such advanced materials appealing in the first place. Within each macromolecule, the intrinsic structure of each unit defines how it interacts with its neighbours, ultimately opening up or denying certain backbone conformations. This eventually also determines how macromolecules interact with each other. This account deals specifically with the conformational problem of many common alkylaromatic units, examining the features of an intramolecular interaction involving a side chain with as few as three methylene groups. A set of 23 model compounds featuring an intramolecular interaction between an aliphatic X-H (X = C, N, O, and S) bond and an aromatic ring was considered. Quantitative computational analysis was made possible, thanks to complete basis set extrapolated CCSD(T) calculations and NCI topological analysis, the latter of which revealed an elaborate network of dispersive and steric interactions leading to somewhat unintuitive and unexpected results, such as the higher energetic stability of certain twisted conformational isomers over those with extended side chains. Vicinal covalent effects from polarizing groups and various heteroatoms, along with the occurrence of non-dispersive phenomena, were also investigated. The conclusions drawn from the investigation include a comprehensive set of guidelines intended to aid in the prediction of the most stable conformation for this class of building blocks. Our findings affect a variety of different research fields, including the tailoring of functional materials for organic electronics and photovoltaics, with insights into a rational treatment of conformational disorder, and the study of protein- and peptide-folding preferences, putting an emphasis on peculiar interactions between the backbone and aromatic residues.

Multi-conformer molecules in solutions: an NMR-based DFT/MP2 conformational study of two glucopyranosides of a vitamin E model compound

Overall geometries of both glucosyl derivatives of PMC were found on the basis of their NMR spectra in CDCl₃and relatedδ_H,C/ⁿJ_HHIEF-PCM(UFF,CHCl₃)/DFT calculational results.

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.

A conformation-selective IR-UV study of the dipeptides Ac-Phe-Ser-NH2 and Ac-Phe-Cys-NH2: probing the SH···O and OH···O hydrogen bond interactions

The conformational preferences of peptides are mainly controlled by the stabilizing effect of intramolecular interactions. In peptides with polar side chains, not only the backbone but also the side chain interactions determine the resulting conformations. In this paper, the conformational preferences of the capped dipeptides Ac-Phe-Ser-NH2 (FS) and Ac-Phe-Cys-NH2 (FC) are resolved under laser-desorbed jet cooling conditions using IR-UV ion dip spectroscopy and density functional theory (DFT) quantum chemistry calculations. As serine (Ser) and cysteine (Cys) only differ in an OH (Ser) or SH (Cys) moiety; this subtle alteration allows us to study the effect of the difference in hydrogen bonding for an OH and SH group in detail, and its effect on the secondary structure. IR absorption spectra are recorded in the NH stretching region (3200-3600 cm(-1)). In combination with quantum chemical calculations the spectra provide a direct view of intramolecular interactions. Here, we show that both FS as FC share a singly γ-folded backbone conformation as the most stable conformer. The hydrogen bond strength of OH···O (FS) is stronger than that of SH···O (FC), resulting in a more compact gamma turn structure. A second conformer is found for FC, showing a β turn interaction.

Exploring the conformational preferences of 20-residue peptides in isolation: Ac-Ala19-Lys + H+vs. Ac-Lys-Ala19 + H+ and the current reach of DFT

Using a high-level density functional and an exhaustive search of conformation space, the predicted conformation of a 20-amino acid peptide explains two seemingly contradictory experiments.

Investigation of the molecular mechanisms of electronic decoherence within a quinone cofactor

The notion of decoherence is particularly adapted to discuss the quantum-to-classical transition in the context of chemical reactions. Decoherence can be modeled by computing the time evolution of nuclear wave packets evolving on distinct potential energy surfaces, here using density functional theory (DFT) and Born–Oppenheimer molecular dynamics simulations. We investigate a redox cofactor of biological interest (tryptophan tryptophylquinone, TTQ) found in the enzyme methylamine dehydrogenase. We also report the first systematic comparison of semi-empirical DFT (tight-binding DFT) and classical force field approaches for estimating decoherence in molecular systems. In the TTQ cofactor, we find that decoherence combines structural and dynamical aspects: it is initiated by the divergent motions of few atoms and then propagates dynamically to the remaining atoms. It is the mass effect of all the atoms that leads to decoherence within a few femtosecond.

QM and QM/MM simulations of proteins

Molecular dynamics simulations of biomolecules have matured into powerful tools of structural biology. In addition to the commonly used empirical force field potentials, quantum mechanical descriptions are gaining popularity for structure optimization and dynamic simulations of peptides and proteins. In this chapter, we introduce methodological developments such as the QM/MM framework and linear-scaling QM that make efficient calculations on large biomolecules possible. We identify the most common scenarios in which quantum descriptions of peptides and proteins are employed, such as structural refinement, force field development, treatment of unusual residues, and predicting spectroscopic and exited state properties. The benefits and shortcomings of QM potentials, in comparison to classical force fields, are discussed, with special emphasis on the sampling problems of protein conformational space. Finally, recent examples of QM/MM calculations in light-sensitive membrane proteins illustrate typical applications of the reviewed methods.

Comprehensive conformational studies of five tripeptides and a deduced method for efficient determinations of peptide structures

Thorough searches on the potential energy surfaces of five tripeptides, GGG, GYG, GWG, TGG, and MGG, were performed by considering all possible combinations of the bond rotational degrees of freedom with a semiempirical and ab initio combined computational approach. Structural characteristics of the obtained stable tripeptide conformers were carefully analyzed. Conformers of the five tripeptides were found to be closely connected with conformers of their constituting dipeptides and amino acids. A method for finding all important tripeptide conformers by optimizing a small number of trial structures generated by suitable superposition of the parent amino acid and dipeptide conformers is thus proposed. Applying the method to another five tripeptides, YGG, FGG, WGG, GFA, and GGF, studied before shows that the new approach is both efficient and reliable by providing the most complete ensembles of tripeptide conformers. The method is further generalized for application to larger peptides by introducing the breeding and mutation concepts in a genetic algorithm way. The generalized method is verified to be capable of finding tetrapeptide conformers with secondary structures of strands, helices, and turns, which are highly populated in larger peptides. This show some promise for the proposed method to be applied for the structural determination of larger peptides.

Theoretical calculation of the OH vibrational overtone spectra of 1,5-pentanediol and 1,6-hexanediol

It is well-known that intramolecular hydrogen bonding affects the relative energetics of conformers, as well as the OH stretching peak positions, intensities, and width. In this study we simulated the Δv(OH) = 3, 4 overtone spectra of 1,5-pentanediol (PeD) and 1,6-hexanediol (HD) using the peak positions, intensities, and width calculated from the B3LYP/6-31+G(d,p) method. Furthermore, room temperature free energy calculations were performed using B3LYP/6-31+G(d,p) MP2/6-31+G(d,p), and MP2/6-311++G(3df,3pd) to obtain the relative population of the conformers. From the calculation of 109 and 381 distinct conformers for PeD and HD, respectively, we find that for these long chain diols the intramolecular hydrogen bonded conformers are not the most dominant conformation at room temperature. This is in stark contrast with shorter chain diols such as ethylene glycol for which the hydrogen bonded conformer dominates the population at room temperature. On the other hand, we found that the correlation between the hydrogen bonded OH red shift versus the homogeneous width, Γ = 0.0155(Δω)(1.36), which was derived for shorter chain diols, is valid even for these longer chain diols. We also showed that the intramolecular hydrogen bonded OH initially decays through the CCOH torsion and COH bending mode no matter how long the alkanediol chain length is for 1,n-alkanediols for n up to 6.

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.

Mutual relationship between stacking and hydrogen bonding in DNA. Theoretical study of guanine-cytosine, guanine-5-methylcytosine, and their dimers

The mutual relationship between stacking and hydrogen-bonding and the possible influence of stacking in the different behavior of cytosine (C) and 5-methylcytosine (C') in DNA have been studied through complete DFT optimization of different structures of G-C and G-C' dimers (i.e., G-C/C-G and G-C'/C'-G), using four different functionals. Our results show that stacking leads to an increase of the O(6)...H-N(4) hydrogen bond length and to a simultaneous decrease of the N(2)-H...O(2) one, in such a way that both lengths approach each other and, in some cases, an inversion occurs. These results suggest that stacking can be a factor to explain the disparity between theory and experiment on the relative strength of the two lateral hydrogen bonds. Regarding the difference between cytosine and 5-methylcytosine, we have shown that methylation enhances the stacking interactions, mainly due to the increase of polarizability. Methylation also favors the existence of slid structures which can produce local distortions of DNA.

Effects of the biological backbone on stacking interactions at DNA-protein interfaces: the interplay between the backbone···π and π···π components

The (gas-phase) MP2/6-31G*(0.25) π···π stacking interactions between the five natural bases and the aromatic amino acids calculated using (truncated) monomers composed of conjugated rings and/or (extended) monomers containing the biological backbone (either the protein backbone or deoxyribose sugar) were previously compared. Although preliminary energetic results indicated that the protein backbone strengthens, while the deoxyribose sugar either strengthens or weakens, the interaction calculated using truncated models, the reasons for these effects were unknown. The present work explains these observations by dissecting the interaction energy of the extended complexes into individual backbone···π and π···π components. Our calculations reveal that the total interaction energy of the extended complex can be predicted as a sum of the backbone···π and π···π components, which indicates that the biological backbone does not significantly affect the ring system through π-polarization. Instead, we find that the backbone can indirectly affect the magnitude of the π···π contribution by changing the relative ring orientations in extended dimers compared with truncated dimers. Furthermore, the strengths of the individual backbone···π contributions are determined to be significant (up to 18 kJ mol(-1)). Therefore, the origin of the energetic change upon model extension is found to result from a balance between an additional (attractive) backbone···π component and differences in the strength of the π···π interaction. In addition, to understand the effects of the biological backbone on the stacking interactions at DNA-protein interfaces in nature, we analyzed the stacking interactions found in select DNA-protein crystal structures, and verified that an additive approach can be used to examine the strength of these interactions in biological complexes. Interestingly, although the presence of attractive backbone···π contacts is qualitatively confirmed using the quantum theory of atoms in molecules (QTAIM), QTAIM electron density analysis is unable to quantitatively predict the additive relationship of these interactions. Most importantly, this work reveals that both the backbone···π and π···π components must be carefully considered to accurately determine the overall stability of DNA-protein assemblies.

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.

Quantum Chemical Benchmark Energy and Geometry Database for Molecular Clusters and Complex Molecular Systems (www.begdb.com): A Users Manual and Examples

Our previous benchmark CCSD(T)/ complete basis set limit calculations were collected into a database named begdb - Benchmark Energy and Geometry DataBase. Web-based interface to this database was prepared and is available at www.begdb.com. Users can browse, search and plot the data online or download structures and energy tables.