Puckering transition of 4-substituted proline residues

The puckering transition of 4-substituted proline residues by electron-withdrawing groups, i.e., 4(R)-hydroxy-L-proline (Hyp) and 4(R)-fluoro-L-proline (Flp) residues, with trans and cis prolyl peptide bonds was studied by adiabatic optimizations along the torsion angle chi1 of the prolyl ring at the HF/6-31+G(d) level. By analyzing the potential energy surface and local minima, it is observed that the puckering transition of the prolyl ring for Hyp and Flp residues proceeds from a down-puckered conformation to an up-puckered one through the transition state with an envelope form having the N atom at the top of envelope and not a planar one for both trans and cis conformers, which is the same as found for the unsubstituted proline residue. At HF/6-31+G(d) and B3LYP/6-311++G(d,p) levels, the structures of the backbone and prolyl ring for local minima of Ac-Hyp-NHMe and Ac-Flp-NHMe are quite similar to those of Ac-Pro-NHMe. However, the relative stability of the up-puckered conformation to the down-puckered one is increased for Ac-Hyp-NHMe with the cis imide bond and for Ac-Flp-NHMe with the trans and cis imide bonds. In particular, the 4(R)-substitution by hydroxy and fluorine groups has brought some structural changes in the prolyl ring of the transition states and the changes in barriers for the puckering transition. The puckering transitions for Ac-Hyp-NHMe and Ac-Flp-NHMe are proven to be predominantly electronically driven by analyzing the electronic and enthalpic contributions to the barriers, as seen for Ac-Pro-NHMe.     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

Theoretical studies on the conformations and properties of phenylazonaphthalenes

The conformations of phenylazo‐2‐naphthalene (I) and phenylazo‐1‐naphthalene (II) have been studied using ab initio methods and density functional theory. The rotational potential energy surfaces about the C N bonds were calculated for both the trans and the cis forms at the PM3 and HF/STO‐3G levels. The PM3 method was found to be incapable for locating the energy minima correctly. The geometries of rotamers obtained from the HF/STO‐3G surface were fully optimized at the HF/6‐31G* and B3LYP/6‐31G* levels. Single‐point MP2/6‐31G* calculations have also been carried out at the HF/6‐31G* geometries. The vibrational spectra, standard thermodynamic functions, heats of formation, and equilibrium molar fractions have been obtained. According to the calculated results, I is comparatively more stable than II. The trans forms have much lower energies than the cis forms, implying that I and II mainly exist in the trans planar forms. The energy gap (ΔE) between the frontier molecular orbitals (MOs) of the cis forms are about 0.7 eV higher than those of the trans ones, indicating that change in geometry from trans to cis form will cause a hypsochromic effect. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 79: 25–33, 2000     

Puckering transition of proline residue in water

The puckering transition of the proline residue with trans and cis prolyl peptide bonds was explored by optimizations along the torsion angle chi1 of the prolyl ring using quantum-chemical methods in water. By analyzing the potential energy surfaces and local minima in water, it is observed that the puckering transition of the proline residue proceeds from a down-puckered conformation to an up-puckered one and vice versa through the transition state with an envelope form having the N atom at the top of the envelope and not a planar one, as seen in the gas phase, although the backbone conformations are different in the gas phase and in water. The barriers to the puckering transition DeltaGup-->down are estimated to be 3.12 and 3.00 kcal/mol for trans and cis conformers at the B3LYP/6-311++G(d,p) level of theory in water, respectively, which are about 1.7 kcal/mol higher than those in the gas phase. Out of 2197 prolines from the 241 high-resolution PDB chains, four transition-state-like structures with the envelope ring puckering are identified. Three of them have the trans prolyl peptide bonds and one has the cis one. The favorable or steric interactions by neighboring residues may be responsible for the stabilization of these transition-state-like ring structures in the proteins.     

Proton affinity of proline and modified prolines using the kinetic method: role of the conformation investigated by ab initio calculations

The proton affinities of proline, cis-3-methylproline and cis-3-ethylproline have been measured by the kinetic method using an ion trap instrument; the values obtained are 936, 940.5, and 943 kJ mol(-1), respectively. The experimental values are consistent with those obtained by high-level ab initio calculations (B3LYP/6-31+G*//B3LYP/6-31G* and B3P86/6-31+G*//B3LYP/6-31G*). Several conformations of neutral and protonated proline were considered, in particular the endo and exo ring structure and the position of the carboxyl group. These results show the importance of the position of the hydrogen atom of the carboxyl group in determining the most stable protonated proline structure.     

Rotational barriers in monomeric CH2=CX-COOH and CH2=CX-CONH2 (X is H or CH3) and vibrational analysis of methacrylic acid and methacrylamide

The internal rotations in acrylic and methacrylic acids CH2=CX-COOH and their amides CH2=CX-CONH2 (X is H or CH3) were investigated by DFT-B3LYP calculations with 6-311+G** basis set. The potential energy curves were consistent with two minima that correspond to planar cis and trans conformation in the case of the acids (or cis and near-trans forms in the case of the amides). Acrylic acid and acrylamide were predicted to have the cis form as the low and predominant conformation of the molecules. In the case of the methacrylic acid and methacrylamide, the conformational relative stability was predicted to reverse as going from the acrylic to the metha compounds. The trans conformer in methacrylic acid or the near-trans in methacrylamide were predicted to be thermodynamically low energy structures of the molecules. The CCCO rotational barrier was calculated to vary from 4 to 6kcal/mol in the four molecules. The OCOH and OCNH torsional barriers were calculated to be about 13 and 22kcal/mol in the acids and the amides, respectively. The vibrational frequencies of methacrylic acid and methacrylamide were computed at the DFT-B3LYP/6-311+G** level and reliable vibrational assignments were made on the basis of normal coordinate analyses and comparison with experimental data of both molecules in their low energy conformations.     

Conformational preference and cis-trans isomerization of 4(R)-substituted proline residues

We report here the conformational preference and prolyl cis-trans isomerization of 4(R)-substituted proline dipeptides, N-acetyl-N'-methylamides of 4(R)-hydroxy-L-proline and 4(R)-fluoro-L-proline (Ac-Hyp-NHMe and Ac-Flp-NHMe, respectively), studied at the HF/6-31+G(d), B3LYP/6-31+G(d), and B3LYP/6-311++G(d,p) levels of theory. The 4(R)-substitution by electron-withdrawing groups did not result in significant changes in backbone torsion angles as well as endocyclic torsion angles of the prolyl ring. However, the small changes in backbone torsion angles phi and psi and the decrease of bond lengths r(Cbeta-Cgamma) or r(Cgamma-Cdelta) appear to induce the increase of the relative stability of the trans up-puckered conformation and to alter the relative stabilities of transition states for prolyl cis-trans isomerization. Solvation free energies of local minima and transition states in chloroform and water were calculated using the conductor-like polarizable continuum model at the HF/6-31+G(d) level of theory. The population of trans up-puckered conformations increases in the order Ac-Pro-NHMe < Ac-Hyp-NHMe < Ac-Flp-NHMe in chloroform and water. The increase in population for trans up-puckered conformations in solution is attributed to the increase in population for the polyproline-II-like conformations with up puckering. The barriers DeltaGct++ to prolyl cis-to-trans isomerization for Ac-Hyp-NHMe and Ac-Flp-NHMe increase as the solvent polarity increases, as seen for Ac-Pro-NHMe. In particular, it was identified that the cis-trans isomerization proceeds through the clockwise rotation about the prolyl peptide bond for Ac-Hyp-NHMe and Ac-Flp-NHMe in chloroform and water, as seen for Ac-Pro-NHMe.     

Density functional computations of proton affinity and gas-phase basicity of proline

The proton affinity and gas-phase basicity of proline were evaluated by using density functional theory coupling the B3-LYP hybrid functional with the extended 6--311++G** basis set. Cis and trans conformations of the carboxyl moiety for both exo and endo ring structures were considered for the neutral proline. The results show that the most stable structure of proline has the endo ring conformation with the carboxyl group in the cis position. The structure at the global minimum is stabilized by an intramolecular hydrogen bond. The nitrogen of the ring in the exo form is the preferred protonation site. The calculated proton affinity (924.3 kJ mol(-1)) and gas-phase basicity (894.4 kJ mol(-1)) are in very good agreement with the experimental counterparts.     

Isomeric structures of protonated ozone: A theoretical study

Ab initio molecular orbital calculations have been used to determine the structure of protonated ozone. Four stable minima were found on the O₃H⁺ singlet potential energy surface. Three forms correspond to ozone protonated at the central oxygen (C_2v) or at the terminal oxygen (two C_s isomers, E and Z). The fourth isomer (C_s) is a derivative of trioxirane. The most stable structure is the planar E form I. The proton affinity of ozone (to give I) is given as 123.6 kcal/mole (MP2/6-31G*//4-31G). The energy difference between I and protonated trioxirane VI is greater than that between ozone and trioxirane.     

A theoretical study of the rotational isomers of nitrosomethanol by semiempirical (AM1) and ab initio methods

The conformational potential energy surface as a function of the two internal torsion angles in C-nitrosomethanol has been obtained using the semiempirical AM1 method. Optimized geometries are reported for the local minima on this surface and also for the corresponding points on the HF/6-31G, 6-31G*, and 6-31G** surfaces. All methods predict cis and trans minima which occur in degenerate pairs, each pair being connected by a transition state of C_s symmetry. The AM1 structures are found to compare well with the corresponding ab initio structures. Ab initio HF/6-31G and HF/6-31G* harmonic vibrational frequencies are reported for the cis and trans forms of nitrosomethanol. When scaled appropriately the calculated frequencies are found to compare well with experimental frequencies. The ab initio calculations predict the energy barrier for cis → trans isomerization to be between 5.8 and 6.5 kcal/mol with the trans → cis isomerization barrier lying between 2.3 and 6.5 kcal/mol. The corresponding AM1 energy barriers are around 1 kcal/mol lower in energy. The ab initio calculations predict the barrier to conversion between the two cis rotamers to be very small with the AM1 value being around 1 kcal/mol. Both AM1 and ab initio calculations predict interconversion between trans rotamers to require between 1.2 and 1.4 kcal/mol.     

UV-laser photoisomerization of fumaryl chloride: the first identification of maleoyl chloride: matrix isolation infrared and ab initio studies

Fourier transform infrared spectra of fumaryl chloride 1 isolated in an argon matrix at 10 K have been analyzed. The comparison between the ab initio HF/6-31G calculated infrared spectra with the experimental ones reveals the existence of three planar conformers, the cis cis 1a, the cis trans 1b and the trans trans 1c. Laser UV irradiation of 1 at lambda = 340 nm yields maleoyl chloride 2 by a carbon carbon double bond photoisomerization process. The first identification of this compound was performed by comparison of the experimental infrared spectra with the calculated ones at the MP2/6-1G** level. AM1 semiempirical and ab initio calculations were used to calculate the structure and the relative stability of the three non planar maleoyl chloride conformers.     

Isomer-specific spectroscopy and conformational isomerization energetics of o-, m-, and p-ethynylstyrenes

The infrared and ultraviolet spectroscopy of o-, m-, and p-ethynylstyrene isomers (oES, mES, and pES) were studied by a combination of methods, including resonance-enhanced two-photon ionization (R2PI), UV-UV hole-burning spectroscopy (UVHB), resonant ion-dip infrared spectroscopy (RIDIRS), and rotationally resolved fluorescence excitation spectroscopy. In addition, the newly developed method of stimulated emission pumping-population transfer spectroscopy (SEP-PTS) was used to determine the energy threshold to conformational isomerization in m-ethynylstyrene. The S(1) <-- S(0) origin transitions of oES and pES occur at 32 369 and 33 407 cm(-1), respectively. In mES, the cis and trans conformations are calculated to be close in energy. In the R2PI spectrum of mES, the two most prominent peaks (32672 and 32926 cm(-1)) were confirmed by UVHB spectroscopy to be S(1) <-- S(0) origins of these two conformers. The red-shifted conformer was identified as the cis structure by least-squares fitting of the rotationally resolved fluorescence excitation spectrum of the origin band. There are also two possible conformations in oES, but transitions due to only one were observed experimentally, as confirmed by UVHB spectroscopy. Density functional theory calculations (B3LYP/6-31+G) predict that the cis-ortho conformer, in which the substituents point toward each other, is about 8 kJ/mol higher in energy than the trans-ortho isomer, and should only be about 5% of the room temperature population of oES. Ground-state infrared spectra in the C-H stretch region (3000-3300 cm(-1)) of each isomer were obtained with RIDIRS. In all three structural isomers, the acetylenic C-H stretch fundamental was split by Fermi resonance. Infrared spectra were also recorded in the excited electronic state, using a UV-IR-UV version of RIDIR spectroscopy. In all three isomers the acetylenic C-H stretch fundamental was unshifted from the ground state, but no Fermi resonance was seen. The first observed and last unobserved transitions in the SEP-PT spectrum were used to place lower and upper bounds on the barrier to cis --> trans isomerization in m-ethynylstyrene of 990-1070 cm(-1). Arguments are given for the lack of a kinetic shift in the measurement. The analogous trans --> cis barrier is in the same range (989-1065 cm(-1)), indicating that the relative energies of the zero-point levels of the two isomers are (E(ZPL)(cis) - E(ZPL)(trans))= -75 to +81 cm(-1). Both the barrier heights and relative energies of the minima are close to those determined by DFT (Becke3LYP/6-31+G) calculations.     

Conformational flexibility in 2,2'-Dioxybiphenyl-chloro-cyclotetraphosphazenes and its relevance to polyphosphazene analogues

The reaction of the cyclotetraphosphazene, [N4P4Cl8], with the difunctional reagent, 2,2'-biphenol, in the presence of potassium carbonate in acetone produced the spiro-substituted derivatives, 2,2'-dioxybiphenylhexachlorocyclotetraphosphazene, bis(2,2'-dioxybiphenyl)tetrachloro-cyclotetraphosphazene, and tris(2,2'-dioxybiphenyl)dichlorocyclotetraphosphazene. Both cis and trans geometrical isomers of the bis compound are observed. Although chromatographic separation of these was unsuccessful, a sample of the trans isomer was obtained by fractional crystallization. The compounds all show non-first-order 31P NMR spectra which were simulated to extract the spectral parameters. Single-crystal X-ray structures of both the trans bis and the tris compounds show that the cyclophosphazene rings exhibit conformational flexibility which gives rise to different crystalline forms being obtained from the same solvent systems. Crystals of trans-bis(2,2'-dioxybiphenyl)tetrachloro-cyclotetraphosphazene were obtained in two different space groups: Pnna (orthorhombic) and P21/n (monoclinic). In the orthorhombic structure, the dominant (72%) conformation of one phosphazene ring is a chair form, and the other (28%) resembles a boat. While for the monoclinic structure, the ring is virtually flat with an oval shape. In both cases the dioxybiphenyl groups are found in R and S configurations in the same molecule and are pi stacked in columns (Pnna) or involved in pi-pi or pi-H interactions (P21/n), thus anchoring the phosphorus atoms of the cyclotetraphosphazenes but still allowing flexibility in the ring conformations. Three crystalline modifications of tris(2,2'-dioxybiphenyl)dichloro-cyclotetraphosphazene were obtained: two in space group P (triclinic), which contained two molecules of dichloromethane in the unit cell, and one solvent-free form in space group P21/n (monoclinic). The cyclophosphazene rings exhibit puckered conformations with the trans-dioxybiphenyl moieties having opposing RS or SR conformations. DFT calculations were carried out on each of the phosphazene ring conformations in trans-bis(2,2'-dioxybiphenyl)tetrachlorocyclotetraphosphazene identified from the X-ray diffraction analysis. It is concluded that intermolecular interactions (i.e., pi-pi or pi-H) between the dioxybiphenyl groups is a factor that modifies the nature of the potential energy surface between the different conformers. The flexibility of the phosphazene ring is supported computationally through the calculated low-energy barriers and experimentally through the highly disordered phosphazene ring conformations observed in the solid state. The results on 2,2'-dioxybiphenyl-substituted cyclotetraphosphazenes provide evidence that microcrystalline domains in their 2,2'-dioxybiphenyl-substituted polyphosphazene analogues will be generated by similar pi-pi and pi-H interactions.     

卤代对苝醌类光敏剂光敏活性影响的量子化学研究

Taking hypomycin B（HMB）as the model compound,HF / 6-31G and TD-B3LYP / 6-31G methods have been employed to explore the effect of chlorine,bromine and iodine substitutions on molecular properties and photosensitization of perylenequinonoid photosensitizer（PQP）. It was found firstly that,the halogen substitutions lowered the EHOMO and ELUMO,and the ΔE. From chlorine,bromine to iodine substitutes,the EHOMO and ELUMO increased,while the corresponding ΔE decreased. Secondly,the halogen substitutions increased the molecular triplet-generating quantum yields and lowered the molecular lowest lying triplet energies,which resulted in the substitutes’similar 1O2 yields with their parent compounds. After halogen substitutions,the molecular adiabatic electron affinities increased,which made the substitutes possess lower O2· - -generating abilities than their parent compounds. Finally,the halogen substitutions lowered the intramolecular hydrogen bond energies,while enhancing the intramolecular proton transfer（IPT）barriers of cis isomers and lowering those of trans isomers on the ground state.     

Substituent effects on the gas-phase fragmentation reactions of sulfonium ion containing peptides

The multistage mass spectrometric (MS/MS and MS3) gas-phase fragmentation reactions of methionine side-chain sulfonium ion containing peptides formed by reaction with a series of para-substituted phenacyl bromide (XBr where X=CH2COC6H4R, and R=--COOH, --COOCH3, --H, --CH3 and --CH2CH3) alkylating reagents have been examined in a linear quadrupole ion trap mass spectrometer. MS/MS of the singly (M+) and multiply ([M++nH](n+1)+) charged precursor ions results in exclusive dissociation at the fixed charge containing side chain, independently of the amino acid composition and precursor ion charge state (i.e., proton mobility). However, loss of the methylphenacyl sulfide side-chain fragment as a neutral versus charged (protonated) species was observed to be highly dependent on the proton mobility of the precursor ion, and the identity of the phenacyl group para-substituent. Molecular orbital calculations were performed at the B3LYP/6-31+G** level of theory to calculate the theoretical proton affinities of the neutral side-chain fragments. The log of the ratio of neutral versus protonated side-chain fragment losses from the derivatized side chain were found to exhibit a linear dependence on the proton affinity of the side-chain fragmentation product, as well as the proton affinities of the peptide product ions. Finally, MS3 dissociation of the nominally identical neutral and protonated loss product ions formed by MS/MS of the [M++H]2+ and [M++2H]3+ precursor ions, respectively, from the peptide GAILM(X)GAILK revealed significant differences in the abundances of the resultant product ions. These results suggest that the protonated peptide product ions formed by gas-phase fragmentation of sulfonium ion containing precursors in an ion trap mass spectrometer do not necessarily undergo intramolecular proton 'scrambling' prior to their further dissociation, in contrast to that previously demonstrated for peptide ions introduced by external ionization sources.     

Conformational behaviour of 2,2′-bipyrrole

The rotational potential around the interannular bond in 2,2-bipyrrole has been calculated making use of standard minimal STO-3G and split valence 4-31G basis sets. Geometrical optimization concerning the most significant interannular internal parameters has been performed with both basis sets. The trans conformer is predicted to be more stable than the cis. The minimal basis set predicts the existence of a cisoid-gauche minimum which after limited optimization becomes very shallow and it seems to be an artifact of the rigid rotor approximation. At 4-31G level, both the trans and cis conformers represent maxima in the potential curve and two gauche minima appear at =46.0° and =147.6°, the latter being the absolute minimum. The absolute maximum of the potential curve corresponds to the cis conformer.     

Ab initio study of the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2) and related C(4)H(10)(2+) isomers.

Ab initio calculations at the MP4(SDTQ)/6-311G//MP2/6-31G level were performed to study the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2), and related C(4)H(10)(2+) isomers. Two doubly hydrogen bridged diborane type trans 1 and cis 2 isomers were located as minima. The trans isomer was found to be more favorable than cis isomer by only 0.6 kcal/mol. Several other minima for C(4)H(10)(2+) were also located. However, the global energy minimum corresponds to C-H (C(4) position) protonated 2-butyl cation 10. Structure 10 was computed to be substantially more stable than 1 by 31.7 kcal/mol. The structure 10 was found to be lower in energy than 2-butyl cation 13 by 34.4 kcal/mol.     

Quantum-chemical calculation of the hydroxylamine molecule and its complex with the proton

The cis and trans conformations of the hydroxylamine molecule were calculated by the SCF LCAO-MO method in the CNDO/2 approximation with optimization of the geometry. It was shown that the cis conformer is more stable in the free state; the barrier to internal rotation of the OH group about the N-O axis was determined. Various schemes for the protonation of hydroxylamine were investigated. The proton affinity of the molecule was calculated (for the N- and O-complexes). It was shown that protonation of NH₂OH at the nitrogen is more favorable, and the trans form of hydroxylamine has the strongest proton affinity.Translated from Teoreticheskaya i éksperimental'naya Khimiya, Vol. 24, No. 1, pp. 122–125, January–February, 1988.     

First-principle computational study on the full conformational space of L-threonine diamide, the energetic stability of cis and trans isomers

First-principle computations were carried out on the conformational space of trans and cis peptide bond isomers of HCO-Thr-NH2. Using the concept of multidimensional conformational analysis (MDCA), geometry optimizations were performed at the B3LYP/6-31G(d) level of theory, and single-point energies as well as thermodynamic functions were calculated at the G3MP2B3 level of theory for the corresponding optimized structures. Two backbone Ramachandran-type potential energy surfaces (PESs) were computed, one each for the cis and trans isomers, keeping the side chain at the fully extended orientation (chi1=chi2=anti). Similarly, two side chain PESs for the cis and trans isomers were generated for the (phi=psi=anti) orientation corresponding to approximately the betaL backbone conformation. Besides correlating the relative Gibbs free energy of the various stable conformations with the number of stabilizing hydrogen bonds, the process of trans-->cis isomerization is discussed in terms of intrinsic stabilities as measured by the computed thermodynamic functions.