Zwitterionic conformers of pyrrolysine and their interactions with metal ions—a theoretical study

A total of 16 pyrrolysine conformers in their zwitterionic forms are studied in gas and simulated aqueous phase using a polarizable continuum model (PCM). These conformers are selected on the basis of our study on the intrinsic conformational properties of non-ionic pyrrolysine molecule in gas phase [Das and Mandal (2013) J Mol Model 19:1695?1704]. In aqueous phase, the stable zwitterionic pyrrolysine conformers are characterized by full geometry optimization and vibrational frequency calculations using B3LYP/6-311++G(d,p) level of theory. Single point calculations are also carried out at MP2/6-311++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. The calculated relative energy range of the conformers at B3LYP/6-311++G(d,p) level is 5.19 kcal mol^?1 whereas the same obtained by single point calculations at MP2/6-311++G(d,p) level is 4.58 kcal mol^?1. A thorough analysis reveals that four types of intramolecular H-bonds are present in the conformers; all of which play key roles in determining the energetics and in imparting the observed conformations to the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of the H-bonds. This study also points out that conformers with diverse structural motifs may differ in their thermodynamical stability by a narrow range of relative energy. The effects of metal coordination on the relative stability order and structural features of the conformers are examined by complexing five zwitterionic conformers of pyrrolysine with Cu⁺² through their carboxylate groups. The interaction enthalpies and Gibbs energies, rotational constants, vibrational frequencies and dipole moments of the metal complexes calculated at B3LYP level are also reported. The zwitterionic conformers of pyrrolysine are not stable in gas phase; after geometry optimization they are converted to the non-ionic forms.     

Quantum-mechanical computations on the electronic structure of trans-resveratrol and trans-piceatannol: a theoretical study of the stacking interactions in trans-resveratrol dimers

Accurate quantum-chemical calculations based on the second-order M?ller-Plesset perturbation method (MP2) and density functional theory (DFT) were performed for the first time to investigate the electronic structures of trans-resveratrol and trans-piceatannol, as well as to study the stacking interaction between trans-resveratrol molecules. Ab initio MP2 calculations performed with using standard split-valence Pople basis sets led us to conclude that these compounds have structures that deviate strongly from planarity, whereas the DFT computations for the same basis sets revealed that the equilibrium geometries of these bioactive polyphenols are planar. Furthermore, the results obtained at the MP2(full)/aug-cc-pVTZ and B3LYP/aug-cc-pVTZ levels indicated that the geometries of trans-resveratrol and trans-piceatannol are practically planar at their absolute energy minima. The relative energies of the equilibrium geometries of trans-resveratrol on its potential energy surface were computed at the MP2(full)/aug-cc-pVTZ level. According to the results obtained, a T-shaped (edge-to-phase) conformer of trans-resveratrol dimer is the most stable in vacuum. This T-shaped conformer is mainly stabilized by strong hydrogen bonding and weak C-H...π interactions. Stacked structures with parallel-displaced trans-stilbene skeletons were also found to be energetically stable. The vertical separation and twist angle dependencies of the stacking energy were investigated at the MP2(full)/aug-cc-pVTZ, B3LYP/aug-cc-pVTZ, and HF/aug-cc-pVTZ levels. The standard B3LYP functional and the Hartree-Fock method neglect long-range attractive dispersion interactions. The MP2 computations revealed that the London dispersion energy cannot be neglected at long or short distances. The stacked model considered here may be useful for predicting the quantum nature of the interactions in π-stacked systems of other naturally occurring stilbenoids, and can help to enhance our understanding of the antioxidant and anticancer activities of trans-resveratrol.     

Low-energy conformers of pamidronate and their intramolecular hydrogen bonds: a DFT and QTAIM study

Extensive DFT and ab initio calculations were performed to characterize the conformational space of pamidronate, a typical pharmaceutical for bone diseases. Mono-, di- and tri-protic states of molecule, relevant for physiological pH range, were investigated for both canonical and zwitterionic tautomers. Semiempirical PM6 method were used for prescreening of the single bond rotamers followed by geometry optimizations at the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels. For numerous identified low energy conformers the final electronic energies were determined at the MP2/6-311++G(2df,2p) level and corrected for thermal effects at B3LYP level. Solvation effects were also considered via the COSMO and C-PCM implicit models. Reasonable agreement was found between bond lengths and angle values in comparison with X-ray crystal structures. Relative equilibrium populations of different conformers were determined from molecular partition functions and the role of electronic, vibrational and rotational degrees of freedom on the stability of conformers were analyzed. For no level of theory is a zwitterionic structure stable in the gas-phase while solvation makes them available depending on the protonation state. Geometrically identified intramolecular hydrogen bonds were analyzed by QTAIM approach. All conformers exhibit strong inter-phosphonate hydrogen bonds and in most of them the alkyl-amine side chain is folded on the P-C-P backbone for further hydrogen bond formation.

Theoretical potential functions and vibrational analysis for halocarbonyl azides CXO-NNN (X=F,Cl and Br)

The structural stability of halocarbonyl azides CXO-NNN (X=F, Cl and Br) was investigated by DFT and MP2 calculations using the 6-311++G** basis set. From the calculations, the molecules were found to have an s-cis<--> s-trans conformational equilibrium with cis being the lower -energy form. Full energy optimizations were carried out for the transition states and the minima at the B3LYP/6 -311++G** and MP2/6 -311++G** levels, from which the rotational barriers were calculated to be of the order 8-10 kcal x mol(-1). The vibrational frequencies were computed at the DFT -B3LYP level and the vibrational assignments for the normal modes of the stable conformers were made on the basis of normal coordinate calculations.     

DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

Rational design of carbonitrile-carboxaldehyde cation receptor models: probing the nature of the heteroatom–metal interaction

Amino acids are constituents of proteins and enzymes which take part almost in all metabolic reactions. Glutamic acid, with an ability to form a negatively charged side chain, plays a major role in intra and intermolecular interactions of proteins, peptides, and enzymes. An exhaustive conformational analysis has been performed for all eight possible forms at B3LYP/cc-pVTZ level. All possible neutral, zwitterionic, protonated, and deprotonated forms of glutamic acid structures have been investigated in solution by using polarizable continuum model mimicking water as the solvent. Nine families based on the dihedral angles have been classified for eight glutamic acid forms. The electrostatic effects included in the solvent model usually stabilize the charged forms more. However, the stability of the zwitterionic form has been underestimated due to the lack of hydrogen bonding between the solute and solvent; therefore, it is observed that compact neutral glutamic acid structures are more stable in solution than they are in vacuum. Our calculations have shown that among all eight possible forms, some are not stable in solution and are immediately converted to other more stable forms. Comparison of isoelectronic glutamic acid forms indicated that one of the structures among possible zwitterionic and anionic forms may dominate over the other possible forms. Additional investigations using explicit solvent models are necessary to determine the stability of charged forms of glutamic acid in solution as our results clearly indicate that hydrogen bonding and its type have a major role in the structure and energy of conformers.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

Concerning the solvent effect in the tautomerism of uracil, 5-fluorouracil, and thymine by density-functional theory and ab initio calculations

The tautomerism of uracil, 5-fluorouracil, and thymine has been investigated in the gas phase and in solution. Electron correlation effects were included in ab initio computations at the MP2 level, and DFT calculations were performed using the B3LYP level. Full geometry optimizations were conducted at the HF/6-31G**, HF/6-31+G**, and B3LYP/6-31+G** levels. Single-point MP2/6-31+G** calculations were performed on the HF/6-31+G** optimized geometries. The influence of the solvent was examined from self-consistent reaction field calculations performed with )=2.21 (1,4-dioxane) and )=78.54 (water). The calculated relative free energies ((G) indicate that substitution of uracil at the position group does not change the relative free energy order of the uracil tautomers in the gas phase and in 1,4-dioxane (except at the MP2 level) whereas this ordering changes in water. Attachment of a fluorine atom changes the relative free energy order of uracil tautomers in the gas phase and in solution.     

Accurate micro- and macro- gas phase basicities of hydroxyl-radical-modified pyrimidines estimated by advanced quantum chemistry methods

The gas phase basicities (GPB) of 16 hydroxyl-radical-derived analogues of model pyrimidine nucleosides were analysed at B3LYP/aug-cc-pVDZ and G3MP2B3 levels of theory. All possible tautomeric equilibriums of neutral and protonated forms of pyrimidine analogues were taken into account. The oxidation of pyrimidines usually reduces the GPB values with respect to canonical values. The only exception was observed in the case of 5,6-dihydroxycytidine, which is more basic than model cytidine if both macro- and micro-scopic measures of basicities are used for comparison. In addition, 6-hydroxycytidine is characterised more by the basic character of the O₂ atom compared to the more basic centre of cytidine, despite the fact that the GPB values of both these compounds are almost identical. Although the B3LYP/aug-cc-pVDZ approach seems to be an accurate and robust method for GPB estimation, the microscopic protonation properties are much more sensitive to this method since the difference in energy between some tautomers is often less than 1 kcal/mol with method-dependent succession. In such cases, G3MP2B3 is recommended.     

Calculating the geometry and Raman spectrum of physiological bis(<Emphasis Type="SmallCaps">l</Emphasis>-histidinato)copper(II): an assessment of DFT functionals for aqueous and isolated systems

Reliable density functional theory (DFT) calculations can be performed in conjuction with spectroscopic measurements to elucidate the structural properties of physiologically important bis(amino acidato)copper(II) compounds in solutions. They can provide insight into the influence of intermolecular interactions on the molecular geometry in the crystal lattice or solution when compared with a DFT gas-phase minimum. Our previous paper [Markovi? et al. (2014) Eur J Inorg Chem 198] reported the DFT-determined geometries and Raman spectra for different conformers of physiological bis(l-histidinato)copper(II) with 20 explicit water molecules, as calculated using the B3LYP functional. The present study examined the reliability of those B3LYP results by applying the M06 functional instead, as it should better account for noncovalent interactions. The water molecules were positioned more compactly around the complex by M06 than by B3LYP. The accuracies of the two functionals when compared to relevant experimental data showed that M06 was better at reproducing in-plane Cu?N bond lengths but B3LYP gave more accurate axial Cu?O distances. Both functionals reproduced the experimental Raman spectrum at pH 8 to similar levels of accuracy and provided precise information on the Cu(II) coordination mode and conformation in aqueous solution. Additionally, we assessed several DFT and DFT-D functionals (BP86, B3LYP, B3LYP-D, M06, M06 L, wB97XD, mPW2PLYPD) by using them to model the geometries of experimental bis(l-histidinato)copper(II) crystalline conformations as isolated systems, and then benchmarking the results against those from high-level second-order pertubation Møller–Plesset (MP2) calculations. Although this assessment resulted in an equivocal conclusion because the MP2 results for the isolated complex were inconsistent with the corresponding DFT outcomes, it does provide new information on future benchmark options.     

Aminequinone-hydroxylquinoneimine tautomeric equilibrium revisited: molecular modeling study of the tautomeric equilibrium and substituent effects in 4-(4-R-phenylamino)naphthalene-1,2-diones

Semi-empirical (AM1 and PM3) and DFT (B3LYP/6-31G(d)) calculations were employed to study the tautomeric equilibrium between the aminequinone A and hydroxylquinoneimine B forms of 4-(4-R-phenylamino)naphthalene-1,2-diones. Substituent effects on the tautomeric equilibrium as well as on geometric and electronic parameters were also determined. In the gas phase the hydroxylquinoneimine B form is the most stable, whereas in water the aminequinone A form becomes more stable. The substituents do not modify the relative energies of the two tautomers. These results are in accordance with experimental data reported in the literature.     

The VA, VCD, Raman and ROA spectra of tri-L-serine in aqueous solution

The structures of one conformer of the nonionic neutral and zwitterionic species of L-serinyl L-serinyl L-serine (SSS or tri-L-serine), together with its cationic and anionic species and the capped N-acetyl tri-L-serine N'-methylamide analog were optimized with density functional theory with the Becke 3LYP hybrid exchange correlation (XC) functional and the PW91 GGA XC functional and the 6-31G* and aug-cc-pVDZ basis sets. Subsequently, the vibrational absorption, vibrational circular dichroism, Raman and Raman optical activity spectra were simulated in order to compare them to experimentally measured spectra. In addition, we compare to previously reported studies for both structural determination and spectral simulations and measurements. A comparison of the various ways to treat the effects of the environment and solvation on both the structure and the spectral properties is thoroughly investigated for one conformer, with the goal to determine which level of theory is appropriate to use in the systematic search of the conformational space. In addition, the effects of the counterion, here Cl- anion, are also investigated. Here we present the current state of the art in nanobiology, where the latest methods in experimental and theoretical vibrational spectroscopy are used to gain useful information about the coupling of the nuclear, electronic and magnetic degrees of freedom and structure of tri-L-serine and its capped peptide analog with the environment.     

Variations of the tautomeric preferences and π-electron delocalization for the neutral and redox forms of purine when proceeding from the gas phase (DFT) to water (PCM)

Quantum-chemical calculations were performed for all possible nine neutral tautomers of purine and their oxidized and reduced forms in water {PCM//DFT(B3LYP)/6?311+G(d,p)} and compared to those in the gas phase {DFT(B3LYP)/6?311+G(d,p)}. PCM hydration influences geometries, π-electron delocalization, and relative energies of purine tautomers in different ways. Generally, the harmonic oscillator model of electron delocalization (HOMED) indices increase when proceeding from the gas phase to aequeous solution for the neutral and redox forms of purine. Their changes for the neutral and oxidized tautomers are almost parallel to the relative energies showing that aromaticity plays an important role in the tautomeric preferences. Tautomeric stabilities and tautomeric preferences vary when proceeding from the gas phase to water indicating additionally that intra- and intermolecular interactions affect tautomeric equilibria. The tautomeric mixture of neutral purine in the gas phase consists mainly of the N9H tautomer, whereas two tautomers (N9H and N7H) dominate in water. For oxidized purine, N9H is favored in the gas phase, whereas N1H in water. A gain of one electron dramatically changes the relative stabilities of the CH and NH tautomers that C6H and C8H dominate in the tautomeric mixture in the gas phase, whereas N3H in water. These variations show exceptional sensitivity of the tautomeric purine system on environment in the electron-transfer reactions.     

The influence of the peptide bond on the conformation of amino acids: a theoretical and FT-IR matrix-isolation study of N-acetylproline

A combined experimental matrix-isolation FT-IR and theoretical study has been performed to investigate the conformational behavior of N-acetylproline. The conformational landscape of N-acetylproline was explored using successively higher computational methods, i.e. HF, DFT(B3LYP) and finally MP2. The exploration resulted in 10 conformations with a relative energy difference smaller than 22 kJ.mol(-1) at the HF/3-21G level of theory. These conformations led to six different conformations after DFT(B3LYP) optimizations. Further optimization at the MP2/6-31++G** level of theory resulted in the same six conformations, all of them with an energy difference smaller than 11.5kJ.mol(-1). One conformation with an intramolecular H-bond was found which was energetically the most favorable conformation. The vibrational and thermodynamical features were calculated using the DFT and MP2 methodologies. In the experimental matrix-isolation FT-IR spectrum, the most stable conformation was dominant and at least two non-H-bonded conformations could be identified. An experimental rotamerization constant between the H-bonded and the other non-H-bonded conformations was estimated and appeared to agree reasonably well with the theoretical MP2 predictions. Some new spectral features of N-acetylproline compared to proline were discovered which might be used to discriminate between the acetylated and non-acetylated form.     

Cation recognition of thiacalix[2]thianthrene and p-tert-butylthiacalix[2]thianthrene and their conformers and complexes with Zn(II), Cd(II) and Hg(II): a theoretical investigation

The structures of thiacalix[2]thianthrene, p-tert-butylthiacalix[2]thianthrene and their complexes with Zn²⁺, Cd²⁺ and Hg²⁺ were obtained using B3LYP/LanL2DZ and HF/LanL2DZ calculations. The structures of the most stable conformers of thiacalix[2]thianthrene and p-tert-butylthiacalix[2]thianthrene optimized at either the B3LYP/LanL2DZ or HF/LanL2DZ level are in good agreement with their corresponding X-ray crystallographic structures. The binding energies of cations, Zn²⁺, Cd²⁺ and Hg²⁺ to thiacalix[2]thianthrene and to p-tert-butylthiacalix[2]thianthrene conformers, and the thermodynamic properties of their associations were obtained. The relative selectivities of both thiacalix[2]thianthrene and thiacalix[2]thianthrene conformer are in same order: Zn²⁺ >> Hg²⁺ > Cd²⁺.      

Effect of microsolvation on zwitterionic glycine: an ab initio and density functional theory study

The effect of microsolvation on zwitterionic glycine, considering both (-NH3(+)) as proton donor and (-COO(-)) as proton acceptor at correlated ab initio (MP2) level and density functional methods (B3LYP, PW91, MPW1PW91 and PBE) using 6-311++G** basis set has been reported. DFT methods have been employed so as to compare the performance/quality of different gradient-corrected correlation functionals (PW91, PBE), hybrid functionals (B3LYP, MPW1PW91) and to predict the near quantitative structural and vibrational properties, at reduced computational cost. B3LYP method outperforms among the different DFT methods for the computed hydrogen bond distances and found closer to the value obtained by correlated MP2 level, whereas MPW1PW91 and PBE methods shows very similar values but approximately 0.03 A less, compared to B3LYP method. MP2 calculation and single point CCSD(T)//MP2 calculation have been considered to decompose the interaction energy, including corrections for basis set superposition error (BSSE). Moreover, charge distribution analysis has also been carried out to understand the long raised questions, how and why the two body energies have significant contribution to the total binding energy.     

De(side chain) model of epothilone: bioconformer interconversions DFT study

Using ab initio methods, we have studied conformations of the de(sidechain)de(dioxy)difluoroepothilone model to quantify the effect of stability change between the exo and endo conformers of the epoxy ring. The DFT minimization of the macrolactone ring reveals four low energy conformers, although MP2 predicted five stable structures. The model tested with DFT hybride functional (B3LYP/6–31+G(d,p)) exhibits the global minimum for one of the exo forms (C), experimentally observed in the solid state, but unexpectedly with the MP2 electron correlation method for the virtual endo form (W). Using the QST3 technique, several pathways were found for the conversion of the low energy conformers to the other low energy exo representatives, as well as within the endo analog subset. The potential energy relationships obtained for several exo forms suggest a high conformational mobility between three, experimentally observed, conformers. The high rotational barrier, however, excludes direct equilibrium with experimental EC-derived endo form S. The highest calculated transition state for the conversion of the most stable exo M interligand to the endo S form is approximately a 28 kcal/mol above the energy of the former. The two-step interconversion of the exo H conformer to the endo S requires at least 28 kcal/mol. Surprisingly, we found that the transition state energy of the H form to the virtual endo W has the acceptable value of about 9 kcal/mol and the next energy barrier for free interconversion of endo W to endo S is 13 kcal/mol. Figure DFT Model of Epothilone Interconversions Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ab Initio Calculations of Vibrational Frequencies,Potential Functions of Internal Rotations and Vibrational Infrared and Raman Spectra for 3,3,3-Trifluoropropanal

The conformational behavior and structure of 3,3,3,-trifluoropropal have been investigated by utilizing ab initio calculations with the 6-31G^** basis set (valence double zeta basis with polarization functions on all atoms) at the restricted Hartree Fock (RHF), second-order Møller-Plesset perturbation (MP2), and Density Functional (B3LYP) levels. The molecule is predicted to have a cis Û gauche conformational equilibrium. Full optimization of the transition states was performed and the rotational barriers of both the CHO and CF₃ rotors were calculated. Vibrational frequencies were computed at the three levels and the zero-point corrections were included into the calculated asymmetric CHO rotational barrier. Complete vibrational assignments were made on the basis of normal coordinate calculations for both stable conformers of the molecule.     

Conformational preferences of α,α-trehalose in gas phase and aqueous solution

This work presents an investigation on the conformational preferences of α,α-trehalose in gas phase and aqueous solution. Eighty-one systematically selected structures were studied at the B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) level, giving rise to 40 unique conformers. The 19 lower energy structures and some selected other were further re-optimized at the B3LYP/6-311++G(d,p) level. The main factors accounting for the conformer’s stability were pointed out and discussed. NBO and QTAIM analyses were performed in some selected conformers in order to address the anomeric and exo-anomeric effects as well as intramolecular hydrogen bonding. The effect of solvent water on the relative stability of the conformers was accounted for by applying the conductor-like polarizable continuum model, CPCM.