The structure of the gas-phase tyrosine–glycine dipeptide

Biological Functions of Antioxidant Dipeptides

In the history of modern nutritional science, understanding antioxidants is one of the major topics. In many cases, food-derived antioxidants have π conjugate or thiol group in their molecular structures because π conjugate stabilizes radical by its delocalization and two thiol groups form a disulfide bond in its antioxidative process. In recent years, antioxidant peptides have received much attention because for their ability to scavenge free radicals, inhibition of lipid peroxidation, chelation of transition metal ions, as well as their additional nutritional value. Among them, dipeptides are attracting much interest as post-amino acids, which have residues in common with amino acids, but also have different physiological properties and functions from those of amino acids. Especially, dipeptides containing moieties of several amino acid (tryptophan, tyrosine, histidine, cysteine, and methionine) possess potent antioxidant activity. This review summarizes previous details of structural property, radical scavenging activity, and biological activity of antioxidant dipeptide. Hopefully, this review will help provide a new insight into the study of the biological functions of antioxidant dipeptides.

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.

A rotational study of the AlaAla dipeptide

Herein, we present the first rotational study of the AlaAla dipeptide, brought into the gas phase by laser ablation. Two different structures have been unveiled in the isolated environment of a supersonic expansion by Fourier transform microwave spectroscopy. These structures have been identified through their rotational and 14N quadrupole coupling constants. The flexibility of the -NH2 and -COOH ends allows the formation of strong intramolecular interactions giving rise to five- and seven-membered ring configurations.

Sugar-peptidic bond interactions: spectroscopic characterization of a model system

Sugars are small carbohydrates which play numerous roles in living organisms such as storage of energy or as structural components. Modifications of specific sites within the glycan chain can modulate a carbohydrate's overall biological function as it happens with nucleic acids and proteins. Hence, identifying discrete carbohydrate modifications and understanding their biological effects is essential. A study of such processes requires of a deep knowledge of the interaction mechanism at the molecular level. Here, we use a combination of laser spectroscopy in jets and quantum mechanical calculations to characterize the interaction between phenyl-β-d-glucopyranoside and N-methylacetamide as a model to understand the interaction between a sugar and a peptide bond. The most stable structure of the molecular aggregate shows that the main interaction between the peptide fragment and the sugar proceeds via a C[double bond, length as m-dash]OH-O2 hydrogen bond. A second conformer was also found, in which the peptide establishes a C[double bond, length as m-dash]OH-O6 hydrogen bond with the hydroxymethyl substituent of the sugar unit. All the conformers present an additional interaction point with the aromatic ring. This particular preference of the peptide for the hydroxyl close to the aromatic ring could explain why glycogenin uses tyrosine in order to convert glucose into glycogen by exposing the O4H hydroxyl group for the other glucoses for the polymerization to take place.

Conformational Interconversions of Amino Acid Derivatives

Exhaustive conformational interconversions including transition structure analyses of N-acetyl-l-glycine-N-methylamide as well as its alanine, serine, and cysteine analogues have been investigated at the MP2/6-31G** level, yielding a total of 142 transition states. Improved estimates of relative energies were obtained by separately extrapolating the Hartree-Fock and MP2 energies to the basis set limit and adding the difference between CCSD(T) and MP2 results with the cc-pVDZ basis set to the extrapolated MP2 results. The performance of eight empirical force fields (AMBER94, AMBER14SB, MM2, MM3, MMFFs, CHARMM22_CMAP, OPLS_2005, and AMOEBAPRO13) in reproducing ab initio energies of transition states was tested. Our results indicate that commonly used class I force fields employing a fixed partial charge model for the electrostatic interaction provide mean errors in the ∼10 kJ/mol range for energies of conformational transition states for amino acid conformers. Modern reparametrized versions, such as CHARMM22_CMAP, and polarizable force fields, such as AMOEBAPRO13, have slightly lower mean errors, but maximal errors are still in the 35 kJ/mol range. There are differences between the force fields in their ability for reproducing conformational transitions classified according to backbone/side-chain or regions in the Ramachandran angles, but the data set is likely too small to draw any general conclusions. Errors in conformational interconversion barriers by ∼10 kJ/mol suggest that the commonly used force field may bias certain types of transitions by several orders of magnitude in rate and thus lead to incorrect dynamics in simulations. It is therefore suggested that information for conformational transition states should be included in parametrizations of new force fields.

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.

Towards an understanding of the nucleation of alpha-para amino benzoic acid from ethanolic solutions: a multi-scale approach

The molecular assembly and subsequent nucleation of para-amino benzoic acid (PABA) from ethanolic solutions is probed using a multi-scale and multi-technique approach. This is applied by examining and interrelating information regarding the molecular, solution-state, cluster, solid-state and surface structures to understand why the alpha form of PABA is crystallised in preference to its low temperature beta form. Calculations suggest that conformational changes within the solute molecule play little or no role in directing the nucleation of either the alpha or beta crystal forms. Combined ab initio and molecular dynamics calculations of the stability of small clusters in solution suggests that the hydrogen-bonded carboxylic acid dimers, present in the alpha structure, are the most stable in solution and play a major role in the self-assembly and polymorphic expression of the alpha form in ethanol in preference to the beta form. These calculations are in good agreement with X-ray small-angle scattering analysis which reveals the presence of PABA clusters in ethanol which are consistent with the size and shape of a carboxylic acid dimer. SAXS studies also reveal the presence of larger cluster structures in a size range 10-40 nm which appear to grow, perhaps reflecting a change in the balance between monomers and dimers within the solution during the nucleation process. The results of crystallisation-kinetics experiments indicate an instantaneous nucleation mechanism where the number of instantaneously nucleated crystallites is calculated to be 1360-660 nuclei per ml and the subsequent growth is found to be only rate limited by diffusion of the growth unit to the crystallite surface. A linear dependence of growth rate with respect to supersaturation is observed for the (0 1 -1) capping face, which is associated with strong π-π stacking interactions. This is consistent with a solid-on-solid mechanism associated with surface roughened growth and concomitant poor lattice-perfection. Conversely, the side (1 0 -1) surface has a growth mechanism consistent with a 2D nucleation birth and spread mechanism. Hence, these mechanisms result in very fast growth along the b-axis and the needle-like morphology that is observed for alpha-PABA.

Comprehensive analysis of energy minima of the 20 natural amino acids

Energy minima of the 20 natural amino acids (capped by a peptide bond at both the N and the C termini, CH3-C(═O)-N(H)-(H)Cα(R)-C(═O)-N(H)-CH3), were obtained by ab initio geometry optimization. Starting with a large number of minima, quickly generated by MarvinView, geometry optimization at the HF/6-31G(d,p) level of theory reduced the number of minima, followed by further optimization at the B3LYP/apc-1 and MP2/cc-pVDZ levels, which caused some minima to disappear and some stable minima to migrate on the Ramachandran map. There is a relation between the number of minima and the size and the flexibility of the side chain. The energy minima of the 20 amino acids are mainly located in the regions of βL, γL, δL, and αL of the Ramachandran map. Multipole moments of atoms occurring in the fragment [-NH-Cα-C(═O)-] common to all 20 amino acids were calculated at the three levels of theory mentioned above. The near parallelism in behavior of these moments between levels of theory is beneficial toward estimating moments with the more expensive B3LYP and MP2 methods from data calculated with the cheaper HF method. Finally, we explored the transferability of properties between different amino acids: the bond length and angles of the common fragment [-NH-Cα(HαCβ)-C'(═O)-] in all amino acids except Gly and Pro. All bond lengths are highly transferable between different amino acids, and the standard deviations are small.

Mimicking anaesthetic-receptor interaction: a combined spectroscopic and computational study of propofol···phenol

Propofol is a general anaesthetic that exerts its action by interaction with the GABA(A) receptor. Crystallographic studies suggest that there is a direct interaction between propofol and the phenolic residue of a tyrosine in the channel. In this study we create propofol···phenol clusters by their co-expansion in jets. The complex is probed using a set of mass-resolved spectroscopic strategies: 2-color REMPI, UV/UV hole-burning, IR/UV double resonance and the novel technique IR/IR/UV triple resonance. The existence of at least six different isomers in the expansion is demonstrated. All the isomers are stabilized by interactions between their aromatic rings. Additionally, in some conformers the OH moieties form hydrogen bonds in some of the isomers, with propofol and phenol alternating their donor-acceptor roles, while in others the -OH···OH angle points to a dipole-dipole interaction. Interpretation of the data in the light of dispersion-corrected DFT calculations shows that shallow barriers separate all the isomers, both in the ground and excited electronic states. Comparison of the structures of the complex with the X-ray diffraction data is also offered.

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.

Extensive conformational searches of 13 representative dipeptides and an efficient method for dipeptide structure determinations based on amino acid conformers

Conformations of peptides are the basis for their property studies and the predictions of peptide structures are highly important in life science but very complex in practice. Here, thorough searches on the potential energy surfaces of 13 representative dipeptides by considering all possible combinations of the bond rotational degrees of freedom are performed using the density functional theory based methods. Careful analyses of the conformers of the 13 dipeptides and the corresponding amino acids reveal the connections between the structures of dipeptide and amino acids. A method for finding all important dipeptide conformers by optimizing a small number of trial structures generated by suitable superposition of the parent amino acid conformations is thus proposed. Applying the method to another eight dipeptides carefully examined by others shows that the new approach is both highly efficient and reliable by providing the most complete ensembles of dipeptide conformers and much improved agreements between the theoretical and experimental IR spectra. The method opens the door for the determination of the stable structures of all dipeptides with a manageable amount of effort. Preliminary result on the applicability of the method to the tripeptide structure determination is also presented. The results are the first step towards proving Anfinsen's hypothesis by revealing the relationships between the structures of the simplest peptide and its constituting amino acids. It implies that the structures of peptides are not only determined by their amino acid sequences, but also closely linked with the amino acid conformations.

Effect of cation complexation on the structure of a conformationally flexible multiply charged anion: stabilization of excess charge in the Na+ x adenosine 5'-triphosphate dianion ion-pair complex

We report a computational study of the conformationally and tautomerically flexible cation-dianion complex of Na(+) with doubly deprotonated adenosine 5'-triphosphate (ATP) using a hierarchical selection method. The method uses molecular dynamics to generate initial conformeric structures, followed by a classification process that groups conformers into five "families" to ensure that a representative sample of structures is retained for further analysis, while very similar conformational structures are eliminated. Hierarchical ab initio calculations (DFT and MP2) of typical conformers of the families are then performed to identify the lowest-energy conformeric structures. The procedure described should provide a useful methodology for conducting higher-level ab initio calculations of medium-sized gas-phase biological molecules for interpreting contemporary laser spectroscopy measurements. For Na(+) x [ATP-2H](2) (considering tautomers where the phosphate chain of ATP is doubly deprotonated), the calculations reveal that the sodium cation interacts directly with the negatively charged phosphates (maximum distance = 2.54 A) in all of the low-energy conformers, while a number of the structures also display close cation-adenine interactions producing compact ball-like structures. These compact structures generally correspond to the lowest-energy conformers. The structural variation between the bare [ATP-2H](2-) molecular ion (Burke et al. J. Phys. Chem. A 2005 , 109 , 9775-9785) and the Na(+) x [ATP-2H](2-) cluster is discussed in detail, including the effect of sodiation on the intramolecular hydrogen-bonding network within ATP in a gas-phase environment.

Vibrational Spectroscopy of Protonated Imidazole and its Complexes with Water Molecules: Ab Initio Anharmonic Calculations and Experiments

The results of anharmonic frequency calculations on neutral imidazole (C3N2H4, Im), protonated imidazole (ImH+), and its complexes with water (ImH+)(H2O)n, are presented and compared to gas phase infrared photodissociation spectroscopy (IRPD) data. Anharmonic frequencies are obtained via ab initio vibrational self-consistent field (VSCF) calculations taking into account pairwise interactions between the normal modes. The key results are: (1) Prediction of anharmonic vibrational frequencies on an MP2 ab initio potential energy surface show excellent agreement with experiment and outstanding improvement over the harmonic frequencies. For example, the ab initio calculated anharmonic frequency for (ImH+)(H2O)N2 exhibits an overall average percentage error of 0.6% from experiment. (2) Anharmonic vibrational frequencies calculated on a semiempirical potential energy surface fitted to ab initio harmonic data represents spectroscopy well, particularly for water complexes. As an example, anharmonic frequencies for (ImH+)H2O and (ImH+)(H2O)2 show an overall average deviation of 1.02% and 1.05% from experiment, respectively. This agreement between theory and experiment also supports the validity and use of the pairwise approximation used in the calculations. (3) Anharmonic coupling due to hydration effects is found to significantly reduce the vibrational frequencies for the NH stretch modes. The frequency of the NH stretch is observed to increase with the removal of a water molecule or replacement of water with N2. This result also indicates the ability of the VSCF method to predict accurate frequencies in a matrix environment. The calculation provides insights into the nature of anharmonic effects in the potential surface. Analysis of percentage anharmoncity in neutral Im and ImH+ shows a higher percentage anharmonicity in the NH and CH stretch modes of neutral Im. Also, we observe that anharmonicity in the NH stretch modes of ImH+ have some contribution from coupling effects, while that of neutral Im has no contribution whatsoever from mode-mode coupling. It is concluded that the incorporation of anharmonic effects in the calculation brings theory and experiment into much closer agreement for these systems.