Searching for the simplest structural units to describe the three-dimensional structure of proteins

A comprehensive conformational space analysis of N-formyl-l-tryptophanamide system by using a genetic algorithm for multi-modal search

The conformational space of protected amino acid HCO-Tryptophane-NH₂ was explored by using a new optimization procedure, in order to localize the stable minima on its potential energy surface (PES). The genetic algorithm based on the Multi-Niche Crowding (MNC) technique was used initially to generate a set of optimized structures for title compound. Resulting structures from the genetic algorithm technique will be used hereafter as input conformers at a hierarchy of increasingly more accurate electronic structure calculations (RHF/6-31G+(d) and DFT/B3LYP/6-31G+(d) geometry optimizations). The lowest energy conformer γ_L(g+g+) presents a folded Backbone that is stabilized by strong hydrogen bond noted C₇. This links the carbonyl oxygen of the formyl group and the hydrogen of the amine group. There are further interactions from one hand between the carbonyl oxygen of the formyl group and the neighboring CH group on the pyrrole ring and from other hand between the N-terminus hydrogen and the indole ring in accordance with the experimental results. This work includes also a comparison between the theoretical calculations and the experimental results of X-ray crystallography extracted from protein data bank (PDB).

Accurate structure, thermodynamic and spectroscopic parameters from CC and CC/DFT schemes: the challenge of the conformational equilibrium in glycine

The structures, relative stabilities, and infrared spectra of the six low-energy conformers of glycine have been characterized using a state-of-the-art quantum-mechanical approach allowing the bond distances, conformational enthalpies and vibrational frequencies to be determined well within the chemical accuracy. Transition state structures governing interconversion among the different energy minima have also been characterized. In detail, the gas-phase thermodynamic properties (at 15 K and 410 K) of the glycine conformers considered have been obtained with a 1 kJ mol(-1) accuracy, and it has been shown that the employment of DFT geometries usually reduces such accuracy by at most 0.1 kJ mol(-1). Regarding molecular structures, the use of two different composite schemes allowed us to further confirm the suitability of a rather cost-effective approach and provide geometrical parameters with an overall accuracy better than 0.002 Å for distances and 1 degree for angles. Thanks to a hybrid CC/DFT approach, the infrared spectra of all conformers considered and of several deuterated isotopologues have been reproduced (when experimental data were available) or predicted with an accuracy of 10 cm(-1). Finally, the joint thermodynamic and spectroscopic investigation allowed us to shed some light on the possible observation of elusive conformers. On the whole, the high accuracy of the computational results allows us to draw a fully consistent interpretation of the available experimental data and to obtain a more complete characterization of the potential energy surface of glycine.

Conformational Preferences of N-Acetyl-l-leucine-N‘-methylamide. Gas-Phase and Solution Calculations on the Model Dipeptide

A DFT study of N-acetyl-l-leucine-N'-methylamide conformers in the gas phase and in solution was carried out. The theoretical computational analysis revealed 43 different conformations at the B3LYP/6-31G(d) level of theory in the gas phase. In addition, the effects of three solvents (water, acetonitrile, and chloroform) were included in the calculations using the isodensity polarizable continuum model (IPCM) and the Poisson-Boltzmann self-consistent reaction field (PB-SCRF) method. The stability order of the different conformers in solution has been analyzed. The theoretical results were compared with some experimental data (X-ray, IR, and NMR).

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.

Probing the competition between secondary structures and local preferences in gas phase isolated peptide backbones

Combining laser desorption with a supersonic expansion together with the selectivity of IR/UV double resonance spectroscopy makes it possible to isolate and characterise the gas phase of remarkable backbone conformations of short peptide chains mimicking protein segments. A systematic bottom-up approach involving a conformer-specific IR study of peptide sequences of increasing sizes has enabled us to map the spectral signatures of the intramolecular interactions, which shape the peptide backbone, in particular H-bonds. The precise data collected are directly comparable to the most sophisticated quantum chemistry calculations of these species and therefore constitute a stringent test for the theoretical methods used. One-residue chains reveal the local conformational preference of the backbone and its dependence upon the nature of the residue. The investigation of longer chains provides evidence for a competition between simple successions of local conformational preferences along the chain and more folded structures, in which a new H-bonding network, involving distant H-bonding sites along the backbone, takes place. From three residues, the issue of helical folding can also be addressed. The present review of the gas phase literature data emphasizes the observation of remarkable secondary structures of biology, including short segments of beta-strands, gamma- and beta-turns, combinations of turns, including a 3(10) helix. It also provides evidence for the flexibility of the peptide chains, i.e., a critical influence of rather minor interactions (like side-chain/backbone interactions) on the conformational stability. Finally, the paper will discuss future promising directions of the present approach.

Secondary structures of short peptide chains in the gas phase: Double resonance spectroscopy of protected dipeptides

The conformational structure of short peptide chains in the gas phase is studied by laser spectroscopy of a series of protected dipeptides, Ac-Xxx-Phe-NH(2), Xxx=Gly, Ala, and Val. The combination of laser desorption with supersonic expansion enables us to vaporize the peptide molecules and cool them internally; IR/UV double resonance spectroscopy in comparison to density functional theory calculations on Ac-Gly-Phe-NH(2) permits us to identify and characterize the conformers populated in the supersonic expansion. Two main conformations, corresponding to secondary structures of proteins, are found to compete in the present experiments. One is composed of a doubly gamma-fold corresponding to the 2(7) ribbon structure. Topologically, this motif is very close to a beta-strand backbone conformation. The second conformation observed is the beta-turn, responsible for the chain reversal in proteins. It is characterized by a relatively weak hydrogen bond linking remote NH and CO groups of the molecule and leading to a ten-membered ring. The present gas phase experiment illustrates the intrinsic folding properties of the peptide chain and the robustness of the beta-turn structure, even in the absence of a solvent. The beta-turn population is found to vary significantly with the residues within the sequence; the Ac-Val-Phe-NH(2) peptide, with its two bulky side chains, exhibits the largest beta-turn population. This suggests that the intrinsic stabilities of the 2(7) ribbon and the beta-turn are very similar and that weakly polar interactions occurring between side chains can be a decisive factor capable of controlling the secondary structure.

Intrinsic Folding of Small Peptide Chains: Spectroscopic Evidence for the Formation of β-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.

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF, MP2, DFT, and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.

Peptide model XXVIII: An exploratory ab initio and density functional study on the side-chain-backbone interaction in N-acetyl-L-cysteine- N-methylamide and N-formyl-L-cysteinamide in their γL-backbone conformations

A conformational and electronic study on the energetically preferred conformations (γL) of N- and C-protected L-cysteine (P-CONH-CH(CH2SH)-CONH-Q, where P and Q may be H or Me) was carried out. After restraining the backbone (BB) conformation to its global minimum (γL or C7eq), all nine possible side-chain (SC) conformations were subjected to geometry optimization at the HF/3–21G and the B3LYP/6–31G(d,p) levels of theory. Seven of the nine side-chain conformers were located on the potential-energy surface. All conformers were subjected to an AIM (atoms in molecules) analysis. This study indicates that three of the seven optimized conformers exhibited either or both SC [Formula: see text] BB- or BB [Formula: see text] SC-type intramolecular hydrogen bonding. Five conformers, however, had distances between a proton and a heteroatom that suggested hydrogen bonding.Key words: L-cysteine diamides, side-chain potential-energy surface, ab initio and DFT geometry optimization, AIM analysis, intramolecular hydrogen bonding.

Peptide models XXVII. An exploratory ab initio study on the 21st amino acid side-chain conformations of N-formyl-L-selenocysteinamide (For-L-Sec-NH2) and N-acetyl-L-selenocysteine-N-methylamide (Ac-L-Sec-NHMe) in their γL backbone conformation

Selenocysteine is expected to have 9 × 9 = 81 conformations [3 × 3 = 9 backbone: ψ (g+,a,g-) × ϕ (g+,a,g-) and 3 × 3 = 9 side-chain: χ1 (g+,a,g-) × χ2 (g+,a,g-)]. In the present study, all the torsional modes of the side-chain (χ1: rotation about the Cα-Cβ and χ2: rotation about the Cβ-Se bonds) were investigated in the relaxed γL backbone [(ϕ,ψ); (g-,g+)] conformation. Seven out of the nine expected minima were found at the RHF/3-21G level of theory for N-formyl-L-selenocysteinamide (For-L-Sec-NH2) and N-acetyl-L-selenocysteine-N-methylamide (Ac-L-Sec-NHMe). The stabilization energy exerted by the -CH2-SeH side-chain has been compared with that of -CH2-SH and -CH2-OH. Relative energies of the various conformers were also obtained via single point calculations at the B3LYP/6-31G(d,p) level of theory. Topological analysis of the electron density has been performed by Bader's Atoms in Molecule (AIM) approach using the results. The structures were also optimized at the B3LYP/6-31+G(d,p) level of theory.Key words: selenocysteine side-chain conformations, ab initio MO study, Multidimensional Conformational Analysis (MDCA), Atoms in Molecules (AIM), Bader's electron density analysis.