Control of helix formation in vinylogous gamma-peptides by (E)- and (Z)-double bonds: a way to ion channels and monomolecular nanotubes

[structure: see text] A complete overview on the alternative and competitive helices in vinylogous gamma-peptides is given, which was obtained on the basis of a systematic conformational analysis at various levels of ab initio MO theory (HF/6-31G*, DFT/B3LYP/6-31G*, PCM/HF/6-31G*). Contrary to the parent gamma-peptides, there is a strict control of helix formation by the configuration of the double bond between the C(alpha) and C(beta) atoms of the monomer constituents. (E)-Double bonds favor helices with larger pseudocycles beginning with 14- up to 27-membered hydrogen-bonded rings, whereas the (Z)-configuration of the double bonds supports a distinct preference of helices with smaller seven- and nine-membered pseudocycles showing interactions between nearest-neighbor peptide bonds. The rather stable helices of the (E)-vinylogous peptides with 22-, 24-, and 27-membered hydrogen-bonded pseudocycles have inner diameters large enough to let molecules or ions pass. Thus, they could be interesting model compounds for the design of membrane channels and monomolecular nanotubes. Since (E)- and (Z)-vinylogous gamma-amino acids and their oligomers are synthetically accessible, our study may stimulate structure research in this novel field of foldamers.     

Hydrazino peptides as foldamers: an extension of the beta-peptide concept.

Replacing the C(beta) atoms in the beta-amino acid constituents of beta-peptides by nitrogen atoms leads to hydrazino peptides. A systematic conformation analysis of blocked hydrazino peptide oligomers of the general type I at the HF/6-31G, MP2/6-31G, and DFT/B3LYP/6-31G levels of ab initio MO theory and on the basis of molecular mechanics reveals a wide variety of secondary structures, as for instance various helices and sheet- and turnlike conformers. Some of them are closely related to secondary structure types found in beta-peptides; others represent novel types. Thus, a very stable, novel helix with 14-membered hydrogen-bonded pseudocycles, which occupies a conformation space different from that of helices with 14-membered rings found among the most stable conformers in beta-peptides, is indicated. The most important secondary structure elements are characterized by interactions between peptidic NH and CO groups. The additional hydrazino N(alpha)H group takes part in special structuring effects but is of lesser importance for secondary structure formation. The influence of environmental effects on the existence and stability of the various structure types is discussed. Due to the wide variety of structural possibilities, hydrazino peptides might be a useful tool for peptide and protein design.     

Helix Formation and Folding in γ‐Peptides and Their Vinylogues

A complete overview of all possible periodic structures with characteristic H‐bonding patterns is provided for oligomers composed of γ‐amino acids (γ‐peptides) and their vinylogues by a systematic conformational search on hexamer model compounds employing ab initio MO theory at various levels of approximation (HF/6‐31G*, DFT/B3LYP/6‐31G*, SCRF/HF/6‐31G*, PCM//HF/6‐31G*). A wide variety of structures with definite backbone conformations and H‐bonds formed in forward and backward directions along the sequence was found in this class of foldamers. All formally conceivable H‐bonded pseudocycles between 7‐ and 24‐membered rings are predicted in the periodic hexamer structures, which are mostly helices. The backbone elongation in comparison to α‐ and β‐peptides allows several possibilities to realize identical H‐bonding patterns. In good agreement with experimental data, helical structures with 14‐ and 9‐membered pseudocycles are most stable. It is shown that the introduction of an (E)‐double bond into the backbone of the γ‐amino acid constituents, which leads to vinylogous γ‐amino acids, supports the folding into helices with larger H‐bonded pseudocycles in the resulting vinylogous γ‐peptides. Due to the considerable potential for secondary‐structure formation, γ‐peptides and their vinylogues might be useful tools in peptide and protein design and even in material sciences.     

Structural and energetic relations between β turns

A systematic quantum chemical study on the structure and stability of the major types of β-turn structures in peptides and proteins was performed at different levels of ab initio MO theory (MP2/6-31G*, HF/6-31G*, HF/3-21G) considering model turns of the general type Ac(SINGLE BOND)X_aa(SINGLE BOND)Y_aa(SINGLE BOND)NHCH₃ with the amino acids glycine, l - and d -alanine, aminoisobutyric acid, and l -proline. The influence of correlation effects, zero-point vibration energies, thermal energies, and entropies on the turn formation was examined. Solvent effects on the turn stabilities were estimated employing quantum chemical continuum approaches (Onsager's self-consistent reaction field and Tomasi's polarizable continuum models). The results provide insight into the geometry and stability relations between the various β-turn subtypes. They show some characteristic deviations from the widely accepted standard rotation angles of β turns. The stability order of the β-turn subtypes depends strongly on the amino acid type. Thus, the replacement of l -amino acids in the two conformation-determining turn positions by d - or α,α-disubstituted amino acid residues generally increases the turn formation tendency and can be used to favor distinct β-turn subtypes in peptide and protein design. The β-turn subtype preferences, depending on amino acid structure modifications, can be well illustrated by molecular dynamics simulations in the gas phase and in aqueous solution. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18 : 1415–1430, 1997     

Potential energy surface of thionylimide

The potential energy surface (PES) of thionylimide has been searched using ab initio MO and density functional calculations. The electronic structures of the isomers of HNSO have been studied using the HF/6‐31+G*, MP2(full)/6‐31+G*, and B3LYP/6‐31+G* levels. Final energies of these molecules have been calculated at the high‐accuracy G2 and CBS‐Q levels. The probable pathways of isomerization of thionylimide to its isomers (e.g., thiocyanic acid, HONS, nitrosothiols) have been explored by studying the three‐ or four‐membered transition states. This study identified total eight possible isomers ( 1–8 ) of HNSO, of which four ( 1–4 ) have already been realized experimentally. Of the remaining four ( 5–8 ), at least two ( 5, 7 ) can be generated experimentally. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theoretical studies on the S-N interactions in sulfoximine

The potential energy surface of sulfoximines has been searched using ab initio MO and Density Functional Calculations. The electronic structures of the isomers of sulfoximine have been studied using HF/6-31+G*, MP2(full)/6-31+G* and B3LYP/6-31+G* levels. Final energies of these molecules have been calculated at the high accuracy G2 and CBS-Q levels. Though a formal SN double bond is generally considered between sulfur and nitrogen in these systems, theoretical studies do not show any π interaction between them. S-N rotational barriers, bond dissociation energies, atomic charge analysis, and NBO analysis all indicate only a single bond across S-N with a very strong ionic interaction.     

An ab initio MO study on the structures and electronic states of hydrogen-bonded O3- HF and SO2-HF complexes

Ab initio molecular orbital (MO) calculations have been carried out for base-hydrogen fluoride (HF) complexes (base = O3 and SO2) in order to elucidate the structures and energetics of the complexes. The ab initio calculations were performed up to the QCISD(T)/6-311++G(d,p) level of theory. In both complexes, hydrogen-bonded structures where the hydrogen of HF orients toward one of the oxygen atoms of bases were obtained as stable forms. The calculations showed that cis and trans isomers exist in both complexes. All calculations for the SO2-HF complex indicated that the cis form is more stable in energy than the trans form. On the other hand, in O3-HF complexes, the stable structures are changed by the ab initio levels of theory used, and the energies of the cis and trans forms are close to each other. From the most sophisticated calculations (QCISD(T)/6-311++G(d,p)//QCISD/6-311+G(d) level), it was predicted that the complex formation energies for cis SO2-HF, trans SO2-HF, cis O3-HF, and trans O3-HF are 6.1, 5.7, 3.4, and 3.6 kcal/mol, respectively, indicating that the binding energy of HF to SO2 is larger than that of O3. The harmonic vibrational frequencies calculated for cis O3-HF and cis SO2-HF complexes were in good agreement with the experimental values measured by Andrews et al. Also, the calculated rotation constants for cis SO2-HF agreed with the experiment.     

Absolute configurations, predominant conformations, and tautomeric structures of enantiomeric tert-butylphenylphosphinothioic acid.

Vibrational absorption and circular dichroism spectra of dextrorotatory, levorotatory, and racemic mixture of tert-butylphenylphosphinothioic acid have been measured in CCl(4) solutions in the 2000-900 cm(-1) region. The conformations for both tautomeric structures of (S)-tert-butylphenylphosphinothioic acid are investigated using the B3LYP functional with the 6-31G* basis set. For the most stable conformation, the absorption and VCD spectra are predicted ab initio using the B3LYP functional with 6-31G*, 6-311G(2d, 2p), 6-31+G, and 6-311G(3df, 3pd) basis sets. A different functional, B3PW91, was also used with the 6-31G* basis set. The predicted spectra are compared to the experimental spectra. The comparison indicates that (-)-tert-butylphenylphosphinothioic acid is of the (S)-configuration and exists in only one tautomeric structure and one conformation in CCl(4) solution.     

Neutron scattering, infra red, Raman spectroscopy and ab initio study of L-threonine.

Inelastic incoherent neutron scattering spectra (IINS) were obtained for normal and deuterated L-threonine. Raman and infrared spectra were also recorded. Geometries were optimized for the zwitterion form using ab initio Hartree-Fock (HF) levels with 6-31G*, 6-311G*, 6-311G** and 6-311 + + G** basis sets. Force fields and normal modes were calculated and used as basis for an assignment of the spectral features.     

Energetics and structure of glycine and alanine based model peptides: Approximate SCC-DFTB,AM1 and PM3 methods in comparison with DFT,HF and MP2 calculations

We calculate relative energies and geometries of important secondary structural elements for small glycine and alanine based polypeptides containing up to eight residues. We compare the performance of the approximate methods AM1, PM3 and self-consistent charge, density-functional tight-binding (SCC-DFTB) to density-functional theory (DFT), Hartree–Fock (HF) and MP2. The SCC-DFTB is able to reproduce structures and relative energies of various peptide models reliably compared to DFT results. The AM1 and PM3 methods show deficiencies in describing important secondary structure elements like extended, helical or turn structures. The discrepancies between different ab initio (HF, MP2) and DFT (B3LYP) methods for medium sized basis sets (6-31G*) also show the need for higher level calculations, since systematic errors found for small molecules may add up when investigating longer polypeptides.     

Computational study of conformational preferences of thioamide-containing azaglycine peptides

The effect of thioamide substitution on the conformational stability of an azaglycine-containing peptide, For-AzaGly-NH2 (1), was investigated for the sake of finding possible applications by using ab initio and DFT methods. As model compounds, For-[psiCSNH]-AzaGly-NH2 (2), For-AzaGly-[psiCSNH]-NH2 (3), and For-[psiCSNH]-AzaGly-[psiCSNH]-NH2 (4) were used. Two-dimensional phi-psi potential energy surfaces (PESs) for 2-4 were calculated at the B3LYP/6-31G*//HF/6-31G* level in gas (epsilon = 1.0) and in water (epsilon = 78.4) by applying the isodensity polarizable continuum model (IPCM) method. On the basis of these PESs, the minimum energy conformations for 2-4 were characterized at the B3LYP level with 6-31G*, 6-311G**, and 6-31+G** basis sets. The remarkable structural effect of thioamide substitution for 2-4 is that beta-strand structure is observed as a global or local minimum. The minima of 2-4 are also compared with those for glycine and thioamide-containing glycine peptides. Our theoretical results demonstrate that compounds 2-4 would be used to design controllable secondary structures.     

A theoretical study of the effect of basis sets on stationary structures for the addition of carbon dioxide to methylamine: A relation among geometries,energy status,and electronic structures

The carbon dioxide (CO₂) to methylamine (CH₃NH₂) addition reaction is investigated using an ab initio MO method with the aim of obtaining an insight into the relation between transition-state structure and reactivity. The complete potential energy hypersurface was explored and the stationary points representing reactants, transition state, and product were localized. To establish the relevance of different basis sets, calculations were carried out at different levels of theory: STO -3G , 3-21G , 4-21G , 4-31G , 6-31G , 6-311G , 6-31G *, and MP 2/6-31G *//HF /6-31G *. A reaction analysis done by correlation of bond-order indices using the More O'Ferrall-Jencks diagram reveals the reaction mechanism as an asynchronous process. The transition state can be described as a four-membered ring. The lone nitrogen electron pair of the amine appears as the driving force of the reaction, where a positive hydrogen rather than a radical is more likely to be transferred. © 1993 John Wiley & Sons, Inc.     

Helix formation in alpha,gamma- and beta,gamma-hybrid peptides: theoretical insights into mimicry of alpha- and beta-peptides

Alpha,gamma- and beta,gamma-hybrid peptides, which are composed of two different homologous amino acid constituents in alternate order, are suggested as novel classes of peptide foldamers. On the basis of a systematic conformational search employing the methods of ab initio MO theory, the possibilities for the formation of periodic secondary structures in these systems are described. The conformational analysis provides a great number of helix conformers widely differing in energy, which can be arranged into three groups: (i) helices with all hydrogen bonds formed in forward direction along the sequence, (ii) helices with all hydrogen bonds in backward direction, and (iii) helices with alternate hydrogen-bond directions (mixed or beta-helices). Most stable are representatives of beta-helices, but their stability decreases considerably in more polar environments in comparison to helix conformers from the other two classes. There is a great similarity between the overall topology of the most stable hybrid peptide helices and typical helices of peptides which are exclusively composed of a single type of homologous amino acids. Thus, the helices of the beta,gamma-hybrid peptides mimic perfectly those of the native alpha-peptides as, for instance, the well-known alpha-helix, whereas the most stable helix conformers of alpha,gamma-hybrid peptides correspond well to the overall structure of beta-peptide helices. The two suggested novel hybrid peptide classes expand considerably the pool of peptide foldamers and may be promising tools in peptide design and in material sciences.     

六氢吡啶和氨形成的氢键团簇C5H10NH(NH3)n(n=1～3)结构的从头算研究

在RHF/6-31G(d)水平下,对C5H10NH(NH3)n(n=1～3)氢键团簇的平衡构型进行了从头算研究,优化得到各种可能的平衡构型.C5H10NH(NH3)为线型氢键结构,而C5H10NH(NH3)2为三元环结构,C5H10NH(NH3)3为四元环结构.在MP2/6-31G(d)//B3LYP/6-31G(d)水平下,对最稳定构型C5H10NH(NH3)n(Ⅰ)(n=1～3)的分子轨道进行布居分析,并且对相应的占据轨道进行指认.C5H10NH(NH3)n(Ⅰ)(n=1～3)垂直电离势的计算结果表明,形成氢键团簇后,分子的垂直电离势降低.     

Molecular orbital predictions of the vibrational frequencies of some molecular ions

Recent spectroscopic advances have led to the first determinations of infrared vibration-rotation bands of polyatomic molecular ions. These initial detections were guided by ab initio predictions of the vibrational frequencies. The calculations reported here predict the vibrational frequencies of additional ions which are candidates for laboratory analysis. Vibrational frequencies of neutral molecules computed at three levels of theory, HF/3-21G, HF/6-31G*, and MP2/6-31G*, were compared with experiment and the effect of scaling was investigated to determine how accurately vibrational frequencies could be predicted. For 92% of the frequencies examined, uniformly scaled HF/6-31G* vibrational frequencies were within 100 cm-1 of experiment with a mean absolute error of 49 cm-1. This relatively simple theory thus seems suitable for predicting vibrational frequencies to guide laboratory spectroscopic searches for ions in the infrared. Hence, the frequencies of 30 molecular ions, many with astrochemical significance,were computed. They are CH2+, CH3+, CH5+, NH2+, NH4+, H3O+, H2F+, SiH2+, PH4+, H3S+, H2Cl+, C2H+, classical C2H3+, nonclassical C2H3+, nonclassical C2H5+, HCNH+, H2CNH2+, H3CNH3+, HCO+, HOC+, H2CO+, H2COH+, H3COH2+, H3CFH+, HN2+, HO2+, C3H+, HOCO+, HCS+, and HSiO+.     

An Ab Initio Molecular Orbital Study of Isophosphinines

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single-point total energy calculation are reported for phosphinine and 13 isophosphinines 7-19 . Isomers 7-11 with an allenic system are calculated to be 8-18 kcal mol m 1 more stable than structures 12-17 with an acetylenic moiety. The calculated energy difference (66.19 kcal mol m 1 ) between phosphinine and the most stable isophosphinine (1-phospha-1,2,4-cyclohexatriene, 10 ) is smaller than the difference (78.96 kcal mol m 1 ) between benzene and the most stable isobenzene (cyclohexa-1,2,4-triene, 2 ). The isophosphinines 18 and 19 , with a butatriene moiety, are calculated to be the least stable isomers.     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

Synthesis and vibrational study of platinum(II) and palladium(II) complexes of glyoxilic acid oxime

New platinum(II) and palladium(II) complexes of glyoxilic acid oxime (gao) have been prepared and characterised by infrared (4000–150 cm⁻¹) and Raman (4000–200 cm⁻¹) spectra. The gao acts as bidentate ligand bonding through the oxime nitrogen and carboxyl oxygen atoms to form neutral bis-chelate square-planar complexes. The lowest energy conformer of the gao ligand (ectt) was selected among 16 theoretically possible conformers on the basis of ab initio calculations at HF/3-21G*, HF/6-31G* and HF/6-311** levels of the theory from which structural parameters and conformational stabilities have been obtained. A complete vibrational assignment of the gao was performed for the lowest energy ectt conformer on the basis of ab initio optimised parameters and normal coordinate analysis calculations (PED). NCA calculations of the complexes studied were also performed.     

Structures and stability of N9, N9 − and N9 + clusters

 Ab initio molecular orbital calculations for N₉, N⁻ ₉ and N⁺ ₉ isomers were carried out at the HF/ 6-31G*, B3PW91/6-31G*, B3LYP/6-31G* and MP2/ 6-31G* levels of theory. Stable equilibrium geometric structures were determined by harmonic vibrational frequency analyses at the HF/6-31G*, B3PW91/6-31G* and B3LYP/6-31G* levels of theory. The most stable free-radical N₉ cluster is structure 1 with C _{2 v} symmetry and that of anion N⁻ ₉ is structure 3 with C _s symmetry. Only one stable structure of the N⁺ ₉ cation with C _{2 v} symmetry was predicted. Their potential application as high-energy-density materials has been examined. Received: 15 June 1999 / Accepted: 11 October 1999 / Published online: 14 March 2000     

Vibrational spectra and structure of benzil and its 18O- and d10-labelled derivatives: a quantum chemical and experimental study

Geometry and vibrational spectroscopic data of benzil-d0 benzil-d10 and benzil-18O calculated at various levels of theory (RHF/6-31G*, B3LYP/6-31G*, BLYP/6-31G*) are reported. The theoretical results are discussed mainly in terms of the comparisons with infrared (4000-100 cm(-1)) and Raman (4000-50 cm(-1)) spectral data. The calculated isotopic frequency shifts, induced by the 18O- and d10-labeling, are in a good agreement with the measured values. A complete vibrational assignment was made with the help of ab initio force field calculations. The data thus obtained were used for reassigning some vibrational frequencies. The results of the optimized molecular structure obtained on the basis of RHF and the DFT calculations are presented and compared with the experimental X-ray diffraction for the benzil-d0 single crystal. It turns out that the best structural parameters are predicted by the B3LYP/6-31G* method.