Conformational preferences and cis-trans isomerization of azaproline residue

The conformational study of N-acetyl-N'-methylamide of azaproline (Ac-azPro-NHMe, the azPro dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the effects of the replacement of the backbone CHalpha group by the nitrogen atom on the conformational preferences and prolyl cis-trans isomerization in the gas phase and in solution (chloroform and water). The incorporation of the Nalpha atom into the prolyl ring results in the different puckering, backbone population, and barriers to prolyl cis-trans isomerization from those of Ac-Pro-NHMe (the Pro dipeptide). In particular, the azPro dipeptide has a dominant backbone conformation D (beta2) with the cis peptide bond preceding the azPro residue in both the gas phase and solution. This may be ascribed to the favorable electrostatic interaction or intramolecular hydrogen bond between the prolyl nitrogen and the amide hydrogen following the azPro residue and to the absence of the unfavorable interactions between electron lone pairs of the acetyl carbonyl oxygen and the prolyl Nalpha. This calculated higher population of the cis peptide bond is consistent with the results from X-ray and NMR experiments. As the solvent polarity increases, the conformations B and B* with the trans peptide bond become more populated and the cis population decreases more, which is opposite to the results for the Pro dipeptide. The conformation B lies between conformations D and A (alpha) and conformation B* is a mirror image of the conformation B on the phi-psi map. The barriers to prolyl cis-trans isomerization for the azPro dipeptide increase with the increase of solvent polarity, and the cis-trans isomerization proceeds through only the clockwise rotation with omega' approximately +120 degrees about the prolyl peptide bond for the azPro dipeptide in the gas phase and in solution, as seen for the Pro dipeptide. The pertinent distance d(N...H-NNHMe) and the pyramidality of imide nitrogen can describe the role of this hydrogen bond in stabilizing the transition state structure and the lower rotational barriers for the azPro dipeptide than those for the Pro dipeptide in the gas phase and in solution.     

ABEEMσπ/MM力场模型下的多肽构象

在ABEEM/MM蛋白质浮动电荷力场模型的基础上,加入孤对电子和 电荷位点,从而能够体现多肽和蛋白质分子中一些重要的各向异性极化性质,允许非化学键方向的电子转移和极化。利用从头计算数据拟合模型相关参数。计算得到的小分子团簇结合能、偶极矩、氢键键长等性质与从头计算结果符合很好。该经典极化模型力场能够重复量子场下丙氨酸二肽、丙氨酸四肽、甘氨酸三肽的各稳定构象,其稳定性顺序与精密从头计算结果相一致,其结构和能量性质较以往模型有一定提高,并优于其他力场模型。     

Understanding organofluorine chemistry. An introduction to the C-F bond

Fluorine is the most electronegative element in the periodic table. When bound to carbon it forms the strongest bonds in organic chemistry and this makes fluorine substitution attractive for the development of pharmaceuticals and a wide range of speciality materials. Although highly polarised, the C-F bond gains stability from the resultant electrostatic attraction between the polarised C delta+ and F delta- atoms. This polarity suppresses lone pair donation from fluorine and in general fluorine is a weak coordinator. However, the C-F bond has interesting properties which can be understood either in terms of electrostatic/dipole interactions or by considering stereoelectronic interactions with neighbouring bonds or lone pairs. In this tutorial review these fundamental aspects of the C-F bond are explored to rationalise the geometry, conformation and reactivity of individual organofluorine compounds.     

Quantitative study of interactions between oxygen lone pair and aromatic rings: substituent effect and the importance of closeness of contact

Current models describe aromatic rings as polar groups based on the fact that benzene and hexafluorobenzene are known to have large and permanent quadrupole moments. This report describes a quantitative study of the interactions between oxygen lone pair and aromatic rings. We found that even electron-rich aromatic rings and oxygen lone pairs exhibit attractive interactions. Free energies of interactions are determined using the triptycene scaffold and the equilibrium constants were determined by low-temperature 1H NMR spectroscopy. An X-ray structure analysis for one of the model compounds confirms the close proximity between the oxygen and the center of the aromatic ring. Theoretical calculations at the MP2/aug-cc-pVTZ level corroborate the experimental results. The origin of attractive interactions was explored by using aromatic rings with a wide range of substituents. The interactions between an oxygen lone pair and an aromatic ring are attractive at van der Waals' distance even with electron-donating substituents. Electron-withdrawing groups increase the strength of the attractive interactions. The results from this study can be only partly rationalized by using the current models of aromatic system. Electrostatic-based models are consistent with the fact that stronger electron-withdrawing groups lead to stronger attractions, but fail to predict or rationalize the fact that weak attractions even exist between electron-rich arenes and oxygen lone pairs. The conclusion from this study is that aromatic rings cannot be treated as a simple quadrupolar functional group at van der Waals' distance. Dispersion forces and local dipole should also be considered.     

Influence of vibrational energy flow on isomerization of flexible molecules: incorporating non-Rice-Ramsperger-Kassel-Marcus kinetics in the simulation of dipeptide isomerization

The conformational isomerization of a dipeptide, N-acetyl-tryptophan methyl amide (NATMA), is studied computationally by including important dynamical corrections to Rice-Ramsperger-Kassel-Marcus (RRKM) theory for the transition rate between pairs of isomers. The dynamical corrections arise from incomplete or sluggish vibrational energy flow in the dipeptide, a property suggested by the mode-selective chemistry that has been observed by Dian et al. [J. Chem. Phys. 120, 133 (2004)]. We compute the extent and rate of vibrational energy flow in NATMA quantum mechanically using local random matrix theory, which we then use to correct the RRKM theory rates. The latter rates are then introduced into a master equation to study the population dynamics of the dipeptide. Incomplete or slow vibrational energy flow is found to enhance the conformational selectivity of NATMA over RRKM estimates.     

Dipole moments and molecular conformations ofO-vinyl- andO-ethylacetoximes in solutions

Dipole moments of the molecules ofO-vinyl- andO-ethylacetoximes in octane and THF were measured at different temperatures. In octane, the equilibrium mixture of stableO-vinylacetoxime conformers consists of forms with different polarity. Polarity of these conformers, which are due to rotation about the C−O and N−O bonds, is mainly determined by the relative position of the lone electron pairs of nitrogen and oxygen atoms. The fraction of the polar conformer increases as temperature increases, as well as on going from octane to THF. Saturation of the vinyl group in theO-vinylacetoxime molecule nearly halves the molecular dipole moment in the nonpolar medium, which radically differentiates this compound from vinyl ethers. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1685–1690, September, 1999.     

Computed electrostatic potentials and average local ionization energies on the molecular surfaces of some tetracyclines

The electrostatic potentials and average local ionization energies computed on the molecular surfaces of four tetracyclines have been investigated with the objective of identifying common features as well as subtle differences that may be related to their biological activities. The four are the parent molecule tetracycline, chlortetracycline, oxytetracycline, and doxycyline. The calculations were carried out at the HF/STO‐3G*//STO‐3G* level. Our electrostatic potential results show that each molecule has a large negative region that extends along its lower portion, consistent with its ability to complex Ca²⁺ and Mg²⁺ ions. Although the surface electrostatic potentials of the four tetracyclines show many similarities, our statistical measure of local polarity allows us to label doxycycline as the one with the lowest degree of local polarity, consistent with its longer half‐life in vivo. The regions in the tetracyclines with the most reactive electrons are the amide nitrogen lone pairs and certain carbons of the outermost rings. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 160–169, 2001     

Dibenzo[a,c]phenazine: a polarity-insensitive hydrogen-bonding probe

A derivative of phenazine, dibenzo[a,c]phenazine (DBPZ), can be used as a very good hydrogen-bonding probe unlike its parent phenazine molecule. Steady-state absorption and fluorescence studies reveal that DBPZ is completely insensitive to polarity of the medium. However, DBPZ can form a hydrogen bond very efficiently in its first excited singlet state. The extent of this excited-state hydrogen-bond formation depends both on size and on hydrogen-bond donor ability of the solvents. Time-resolved fluorescence studies and theoretical calculations also suggest that this hydrogen-bond formation is much more favorable in the excited state as compared to the ground state. In the excited state, the electron density is pushed toward the nitrogen atoms from the benzene rings, thereby increasing the dipole moment of the DBPZ molecule. Although the dipole moment of DBPZ increases upon photoexcitation, like other polarity probes, the molecule remains fully insensitive to the polarity of the interacting solvent. This unusual behavior of DBPZ as compared to simple phenazine and other polarity probes is due to the structure of the molecule. Hydrogen atoms at the 1 and 8 positions of DBPZ are sterically interacting with a lone pair of electrons on the proximate nitrogen atoms and make both of the nitrogen atoms inaccessible to solvent molecules. For this reason, DBPZ cannot sense the polarity of the medium. However, DBPZ can only sense solvents, those that have hydrogen with some electropositive nature, that is, the hydrogen-bond donating solvents. Hydrogen being the smallest among all elements can only interact with the lone pair of electrons of nitrogen atoms. Thus, DBPZ can act as a sensor for the hydrogen-bond donating solvents irrespective of their dielectrics.     

Advancing beyond charge analysis using the electronic localization function: Chemically intuitive distribution of electrostatic moments

We propose here an evaluation of chemically intuitive distributed electrostatic moments using the topological analysis of the electron localization function (ELF). As this partition of the total charge density provides an accurate representation of the molecular dipole, the distributed electrostatic moments based on the ELF partition (DEMEP) allows computing of local moments located at non atomic centers such as lone pairs, sigma bonds and pi systems. As the local dipole contribution can be decomposed in polarization and charge transfer components, our results indicate that local dipolar polarization of the lone pairs and chemical reactivity are closely related whereas the charge transfer contribution is the key factor driving the local bond dipole. Results on relevant molecules show that local dipole contributions can be used to rationalize inductive polarization effects in alcohols derivatives and typical hydrogen bond interactions. Moreover, bond quadrupole polarization moments being related to a pi character enable to discuss bond multiplicities, and to sort families of molecules according to their bond order. That way, the nature of the C-O bond has been revisited for several typical systems by means of the DEMEP analysis which appears also helpful to discuss aromaticity. Special attention has been given to the carbon monoxide molecule, to the CuCO complex and to a weak intramolecular N|-CO interaction involved in several biological systems. In this latter case, it is confirmed that the bond formation is mainly linked to the CO bond polarization. Transferability tests show that the approach is suitable for the design of advanced force fields.     

Fine Tuning of Heme Reactivity: Hydrogen-Bonding and Dipole Interactions Affecting Ligand Binding to Hemoproteins

Heme reactivity in hemoproteins is governed by the microenvironment near the ligand binding site. In order to quantify polarity effects on heme ligand binding, the kinetics of O₂ and CO binding have been measured for a series of synthetic heme models equipped with a range of groups of varying dipole moments positioned near the heme coordination site. For hemes with polar aprotic groups, both O₂ on (k′) and off rates (k) are found to be dependent on the dipole moment. For model systems containing protic groups, the O₂ off rate is substantially reduced due to hydrogen bonding with the coordinated O₂. The hydrogen-bonding stabilization is estimated to be 0.7 and 1.6 kcal/mol for an alcohol and a secondary amide, respectively. CO binding displays little correlation with a polarity effect; instead it seems to depend upon the size and position of the polar group.     

Solvent polarity induced structural changes in 2,6-diamino-9,10-anthraquinone dye

Photophysical properties of 2,6-diamino-9,10-anthraquinone (2,6-DAAQ) dye have been investigated in different solvents and solvent mixtures. The fluorescence quantum yields, fluorescence lifetimes, radiative rate constants, nonradiative rate constants and absorption and fluorescence spectral characteristics show unusual deviations in the lower polarity aprotic solvents in comparison to those in other aprotic solvents of medium to higher polarities. The results indicate that the dye exists in different structural forms in the lower and in the medium to higher polarity solvents. Drawing an analogy with the results reported for other amino-substituted dyes, it is inferred that 2,6-DAAQ dye adopts a planar intramolecular charge transfer (ICT) structure in medium to higher polarity solvents, where the amino lone pairs are in good resonance with the anthraquinone pi-cloud. In the lower polarity solvents, however, the dye is inferred to exist in a nonplanar structure where the amino lone pairs are not in good resonance with the anthraquinone pi-cloud. Due to these structural differences, the dye displays significantly different photophysical behavior in the lower polarity solvents than in the other solvents of medium to higher polarities. Supportive evidence for the above structural changes has been obtained from ab initio quantum chemical calculations on the structures of the dye under different conditions. Unusual deviations in the photophysical properties of 2,6-DAAQ dye in protic solvents in comparison to those in aprotic solvents of similar polarities are attributed to the intermolecular hydrogen bonding effect involving the OH groups of the protic solvents and the quinonoid oxygens of the dye.     

An interpretation of the enhancement of the water dipole moment due to the presence of other water molecules

The dipole moment of the gas phase water monomer is 1.85 D. When solvated in bulk water, the dipole moment of an individual water molecule is observed to be enhanced to the much larger value of 2.9 +/- 0.6 D. To understand the origin of this dipole moment enhancement, the effective fragment potential (EFP) method is used to solvate an ab initio water molecule to predict the dipole moments for various cluster sizes. The dipole moment as a function of cluster size, nH 2O, is investigated [for n = 6-20 (even n), 26, 32, 41, and 50]. Localized charge distributions are used in conjunction with localized molecular orbitals to interpret the dipole moment enhancement. These calculations suggest that the enhancement of the dipole moment originates from the decrease of the angle between the dipole vectors of the lone pairs on oxygen as the number of hydrogen bonds to that oxygen increases. Thus, the decreased angle, and the consequent increase in water dipole moment, is most likely to occur in environments with a larger number of hydrogen bonds, such as the center of a cluster of water molecules.     

Conformational preferences of N-methoxycarbonyl proline dipeptide

The conformational study on N-methoxycarbonyl-L-proline-N'-methylamide (Moc-Pro-NHMe, prolylcarbamate) is carried out using ab initio HF and density functional B3LYP methods with the self-consistent reaction field method in the gas phase and in solution (chloroform, acetonitrile, and water). The replacement of the N-acetyl group by the N-methoxycarbonyl group results in the changes in conformational preferences, populations for backbone and prolyl puckering, and barriers to cis-trans isomerization of the prolyl residue in the gas phase and in solution, although there are small changes in the geometry of the prolyl peptide bond and the torsion angles of backbone and prolyl ring. The cis population increases with the increase of solvent polarity, as found for Ac-Pro-NHMe (prolylamide), but it is amplified by 9% in the gas phase and about 17% in solution for prolylcarbamate compared with those for prolylamide. It is found that the cis-trans isomerization for prolylcarbamate proceeds through the clockwise rotation with omega' approximately +120 degrees about the prolyl peptide bond in the gas phase and in solution, as found for prolylamide. However, the rotational barriers to the cis-trans isomerization for prolylcarbamate are calculated to be 3.7-4.7 kcal/mol lower than those of prolylamide in the gas phase and in solution, and are found to be less sensitive to the solvent polarity. The calculated rotational barriers for prolylcarbamate in chloroform and water are in good agreement with the observed values. The shorter hydrogen-bond distance between the prolyl nitrogen and the amide H (H(NHMe)) of the NHMe group, the decrease in electron overlap of the prolyl C-N bond, and the favorable electrostatic interaction between the ester oxygen and the amide H(NHMe) for the transition state seem to play a role in lowering the rotational barrier of prolylcarbamate. The smaller molecular dipole moments of the ground- and transition-state structures for prolylcarbamate in the gas phase and in solution seem to be one of factors to make the rotational barrier less sensitive to the solvent polarity. As the solvent polarity increases (i.e., from the gas phase to chloroform to acetonitrile), the value of DeltaH(tc)(double dagger) decreases and the magnitude of DeltaS(tc)(double dagger) increases for prolylcarbamate, which results in a nearly constant value of the rotational barrier.     

Conformational preferences and cis-trans isomerization of L-lactic acid residue

The conformational study on N-acetyl- N'-methylamide of l-lactic acid (Ac-Lac-NHMe, the Lac dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore its backbone conformational preferences and cis-trans isomerization for the depsipeptide with an ester bond in the gas phase and in solution. In the gas phase and in chloroform, the conformation tB with a trans depsipeptide bond is most preferred for the Lac dipeptide, whose backbone torsion angles are phi approximately -150 degrees and psi approximately -5 degrees , juxtaposed to those of the 3 10-helical structure. The larger shift in phi is brought to reduce the repulsion between the two carbonyl carbons of the acetyl and NHMe groups. However, the polyproline II-like tF conformation becomes more populated and the relative stability of conformation tB decreases significantly as the solvent polarity increases. This may be ascribed to weakening a C(5) hydrogen bond between the depsipeptidyl oxygen and the carboxyl amide hydrogen that plays a role in stabilizing the conformation tB in the gas phase and in chloroform. The cis populations about the depsipeptide bond are nearly negligible in the gas phase and in solution. The rotational barriers to the cis-trans isomerization of the depsipeptide bond for the Lac dipeptide are calculated to be about 11 kcal/mol, which is about half of those for the Ala dipeptide, although they increase somewhat with the increase of solvent polarity. The cis-trans isomerization of the depsipeptide bond proceeds through either clockwise or anticlockwise rotations with torsion angles of about +90 degrees or -90 degrees , respectively, in the gas phase and in solution, whereas it has been known that the isomerization proceeds through only the clockwise rotation for alanyl and prolyl peptide bonds. The pertinent distances between the depsipeptidyl oxygen and the carboxyl amide hydrogen can describe the role of this hydrogen bond in stabilizing the transition state structures in the gas phase and in solution.     

Polarization-angle-scanning two-dimensional spectroscopy: application to dipeptide structure determination

Coherent two-dimensional optical spectroscopy based on a heterodyne-detected stimulated photon echo measurement technique requires four ultrashort pulses whose pulse-to-pulse delay times, wavevectors, and frequencies are experimentally controllable variables. In addition, the polarization directions of the four radiations can also be arbitrarily adjusted. We show that the polarization-angle-scanning two-dimensional spectroscopy can be of effective use to selectively suppress either all the diagonal peaks or a cross-peak in a given two-dimensional spectrum. Theoretical relationships between the transition dipole vectors of a given pair of coupled modes or quantum transitions and the polarization angle configuration making the corresponding cross-peak vanish are established. Here, to shed light into the underlying principles of the polarization-angle-scanning two-dimensional spectroscopy, we considered the amide I vibrations of various isotope-labeled dipeptide conformers and show that one can selectively suppress a cross-peak by properly controlling the polarization angle of a chosen beam among them. Once the relative directions of the amide I transition dipole vectors are determined using the polarization-angle-scanning technique theoretically proposed here, they can serve as a set of constraints for determining structures of model peptides. The present work demonstrates that the polarization-controlled two-dimensional vibrational or electronic spectroscopy can provide invaluable information on intricate details of molecular structures.     

Solvent effects on the thioamide rotational barrier: an experimental and theoretical study.

The solvent effect on the C-N rotational barriers of N,N-dimethylthioformamide (DMTF) and N,N-dimethylthioacetamide (DMTA) has been investigated using ab initio theory and NMR spectroscopy. Selective inversion recovery NMR experiments were used to measure rotational barriers in a series of solvents. These data are compared to ab initio results at the G2(MP2) theoretical level. The latter are corrected for large amplitude vibrational motions to give differences in free energy. The calculated gas phase barriers are in very good agreement with the experimental values. Solvation effects were calculated using reaction field theory. This approach has been found to give barriers that are in good agreement with experiment for many aprotic, nonaromatic solvents that do not engage in specific interactions with the solute molecules. The calculated solution-phase barriers for the thioamides using the above solvents are also in good agreement with the observed barriers. The solvent effect on the thioamide rotational barrier is larger than that for the amides because the thioamides have a larger ground-state dipole moment, and there is a larger change in dipole moment with increasing solvent polarity. The transition-state dipole moments for the amides and thioamides are relatively similar. The origin of the C-N rotational barrier and its relation to the concept of amide "resonance" is examined.     

Amide I vibrational circular dichroism of dipeptide: Conformation dependence and fragment analysis

The amide I vibrational circular dichroic response of alanine dipeptide analog (ADA) was theoretically investigated and the density functional theory calculation and fragment analysis results are presented. A variety of vibrational spectroscopic properties, local and normal mode frequencies, coupling constant, dipole, and rotational strengths, are calculated by varying two dihedral angles determining the three-dimensional ADA conformation. Considering two monopeptide fragments separately, we show that the amide I vibrational circular dichroism of the ADA can be quantitatively predicted. For several representative conformations of the model ADA, vibrational circular dichroism spectra are calculated by using both the density functional theory calculation and fragment analysis methods.