Aqueous solvent effects on the conformational space of tryptamine. Structural and electronic analysis

The TRA (3-[2-aminoethyl]indole) is an important neurotransmitter with a close structural and chemical similarity to the neurotransmitter serotonin (5-hydroxytryptamine), and to melatonin (5-methoxy-N-acetyltryptamine), which plays a key role in daily human behavior. Moreover, TRA, and other indolic compounds are very efficient antioxidants. In this work the conformational space of TRA was scanned in aqueous solution, simulating the solvent by the polarizable continuum model. Geometry optimizations were performed at B3LYP/6-31+G** level. Electronic distributions were analyzed at a better calculation level, thus improving the basis set (6-311++G**). A topological study based on Bader’s theory (atoms in molecules) and natural bond orbital (NBO) framework was performed. Structural changes found in solution were related with charge delocalization mechanisms, which explained the changes in the conformational relative population in aqueous phase. Solvent effects on molecular electrostatic potential (MEPs) were also quantified and rationalized through charge delocalization mechanisms, thus connecting changes in MEPs with changes in structure, bond polarization, orbital bonding populations, natural charges, and bond topological properties. Moreover, polarizabilities and dipolar moments were calculated. All conformers were taken into account. Our results are the first prediction of TRA polarizabilities. The results reported contribute to the understanding of the structure, stability and reactivity of TRA and other indole derivatives.     

Electronic structure and conformational properties of <Emphasis Type="Italic">1H</Emphasis>-indole-3-acetic acid

The conformational space of 1H-Indole-3-Acetic Acid (IAA) was scanned using molecular dynamics at semiempirical level, and complemented with functional density calculations at B3LYP/6-31G** level, 14 conformers of lowest energy were obtained. Electronic distributions were analyzed at a higher calculation level, thus improving the basis set (B3LYP/6-311++G**). A topological study based on Bader’s theory (AIM: atoms in molecules) and natural bond orbital (NBO) framework performed with the aim to analyze the stability and reactivity of the conformers allowed the understanding of electronic aspects relevant in the study of the antioxidant properties of IAA. Intramolecular hydrogen bonds were found and were characterized as blue-shifting hydrogen bonding interactions. Furthermore, molecular electrostatic potential maps (MEPs) were obtained and analyzed in the light of AIM and NBO results, thus showing subtle but essential features related not only to reactivity but also with intramolecular weak interactions, charge delocalization and structure stabilization.     

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.     

Theoretical studies on sulfanilamide and derivatives with antibacterial activity: conformational and electronic analysis

Quantum chemical methods have been used to study the conformational and electronic properties of sulfanilamide and derivatives with antibacterial activity. Calculations at B3LYP/6-311++G(3df,2p) level of theory predict the existence of four conformers for sulfanilamide depending on the orientation of p-amino and amide groups. Focusing on the sulfonamide moiety, amide NH₂ and SO₂ groups could exist either in an eclipsed or staggered arrangement. Gas-phase results predict the eclipsed conformer to be most stable but opposite to what has been rationalized previously, no stabilizing hydrogen bonds between those groups has been found through NBO analysis. When solvent effect is taken into account through the IEF-PCM method, staggered conformer is preferred; in fact, eclipsed conformation changed when explicit solvent molecules were included. Conformational analysis of all derivatives has shown two global minima which are specular images. Five out of the seven derivatives studied adopted a particular minimum energy conformation with very similar geometries.     

Experimental and theoretical DFT (B3LYP,X3LYP,CAM-B3LYP and M06-2X) study on electronic structure,spectral features,hydrogen bonding and solvent effects of 4-methylthiadiazole-5-carboxylic acid

Detailed structural, electronic and spectroscopic study of 4-methylthiadiazole-5-carboxylic acid, one of the simplest 1,2,3-thiadiazole derivatives has been performed using density functional theory at four different functionals (B3LYP, X3LYP, CAM-B3LYP and M06-2X). The two possible conformers and their dimeric forms have been investigated for the stability and hence for the calculation of molecular properties of the title compound. Vibrational analysis has been performed with the help of experimental FT-IR and FT-Raman spectra. NBO analysis has been performed to estimate the N–H—O=C hydrogen bond strength and to evaluate the intra and inter molecular charge transfer in the system. Intermolecular hydrogen-bond strength has also been computed using Atoms in Molecules (AIM) theory. To visualise spatial domain, key sites of electron transitions and electron density difference between ground as well as excited states, and their 2D and 3D plots have been computed. Solvent effect on the intermolecular hydrogen bonding have also been investigated using solvents of different polarities. Non-linear optical properties, molecular electrostatic potential surface map (MESP), thermodynamic potentials at different temperatures have also been computed and plotted.     

Structural and electronic properties of Z isomers of (4α→6´´,2α→O→1´´)-phenylflavans substituted with R = H, OH and OCH₃ calculated in aqueous solution with PCM solvation model

In the search for new antioxidants, flavan structures called our attention, as substructures of many important natural compounds, including catechins (flavan-3-ols), simple and dimeric proanthocyanidins, and condensed tannins. In this work the conformational space of the Z-isomers of (4α→6′′, 2α→O→1′′)-phenylflavans substituted with R = H, OH and OCH₃ was scanned in aqueous solution, simulating the solvent by the polarizable continuum model (PCM). Geometry optimizations were performed at B3LYP/6-31 G** level. Electronic distributions were analyzed at a better calculation level, thus improving the basis set (6-311++G**). A topological study based on Bader′s theory (atoms in molecules) and natural bond orbital (NBO) framework was performed. Furthermore, molecular electrostatic potential maps (MEPs) were obtained and thoroughly analyzed. The stereochemistry was discussed, and the effect of the solvent was addressed. Moreover, intrinsic properties were identified, focusing on factors that may be related to their antioxidant properties. Hyperconjugative and inductive effects were described. The coordinated NBO/AIM analysis allowed us to rationalize the changes of MEPs in a polar solvent. To investigate the molecular and structural properties of these compounds in biological media, the polarizabilities and dipolar moments were predicted which were further used to enlighten stability and reactivity properties. All conformers were taken into account. Relevant stereoelectronic aspects were described for understanding the stabilization and antioxidant function of these structures.     

Quantum mechanical studies of lincosamides

Lincosamides are a class of antibiotics used both in clinical and veterinary practice for a wide range of pathogens. This group of drugs inhibits the activity of the bacterial ribosome by binding to the 23S RNA of the large ribosomal subunit and blocking protein synthesis. Currently, three X-ray structures of the ribosome in complex with clindamycin are available in the Protein Data Bank, which reveal that there are two distinct conformations of the pyrrolidinyl propyl group of the bound clindamycin. In this work, we used quantum mechanical methods to investigate the probable conformations of clindamycin in order to explain the two binding modes in the ribosomal 23S RNA. We studied three lincosamide antibiotics: clindamycin, lincomycin, and pirlimycin at the B3LYP level with the 6-31G** basis set. The focus of our work was to connect the conformational landscape and electron densities of the two clindamycin conformers found experimentally with their physicochemical properties. For both functional conformers, we applied natural bond orbital (NBO) analysis and the atoms in molecules (AIM) theory, and calculated the NMR parameters. Based on the results obtained, we were able to show that the structure with the intramolecular hydrogen bond C=O…H-O is the most stable conformer of clindamycin. The charge transfer between the pyrrolidine-derivative ring and the six-atom sugar (methylthiolincosamide), which are linked via an amide bond, was found to be the dominant factor influencing the high stability of this conformer.     

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.

Quantum chemical investigation of intramolecular thione-thiol tautomerism of 1,2,4-triazole-3-thione and its disubstituted derivatives

The one step intramolecular thione-thiol tautomerism of 1,2,4-triazole-3-thione and its disubstituted derivatives has been studied through the use of electronic structure methods. Due to the absence of experimental data for the parent molecule of 1,2,4-triazole-3-thione the structure and energetics of aforementioned tautomers were derived using various basis sets and levels including HF, B3LYP, and MP2 methods. The gas phase results show that in all different levels of theory the most stable tautomer is the thione form. It has also been revealed that B3LYP/6-31G(d,p) level is quite well suited and reliable to investigate these kinds of tautomerism. To account the influence of substituents on the mentioned tautomerization, the tautomerism and conformational properties as well as vibrational analysis of 20 halophenyl and isopyridyl derivatives were investigated using B3LYP/6-31G(d,p) calculations. In all cases the calculations indicate that substituents have no considerable effects on relative stabilities and energy barriers for the thione-thiol proton transfer and the thione forms are the predominant species in the gas phase. In order to figure out the relative stabilities of the species involved in the tautomerism, geometrical and natural bond orbital (NBO) analyses have been employed. It has also been shown that the computed vibrational frequencies of tautomers with different scaling factors could be used to interpret the vibrational frequencies in IR spectrum of similar species.     

Fluoro-substitution effects in deoxyfluoro-D-glucose derivatives: random conformational search and quantum chemical calculation

The effect of substitution by the fluorine atom at different positions of D-glucose was investigated by quantum chemical calculation of the low-energy conformers. These were obtained through the Random conformational search method. The geometries of conformers were optimized at the RHF/6-31(d) level, then reoptimization and vibrational analysis were performed at the B3LYP/6-31+G(d) level. Single-point energies were calculated at the B3LYP/6-311++G(2d,2p) level. The free energies of solvation in water were calculated utilizing the AM1-SM5.4 solvation model. For all substitution positions, the ring conformation does not change much, and the pyranoid 4C1 conformers are dominant, while variations in the substitution site result in different effects in the network of hydrogen bonds, anomeric effect, the solvation free energy, and the ratio of alpha- and beta-anomers.     

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.     

Phenylisoserine in the gas-phase and water: Ab initio studies on neutral and zwitterion conformers

The conformational landscape of phenylisoserine (PhIS) was studied. Trial structures were generated by allowing for all combinations of single-bond rotamers. Based on the B3LYP/aug-cc-pVDZ calculations 54 conformers were found to be stable in the gas phase. The six most stable conformers were further optimized at the B3LYP/aug-cc-pVTZ and MP2/aug-cc-pVDZ levels for which characteristic intramolecular hydrogen bond types were classified. To estimate the influence of water on PhIS conformation, the IEF-PCM/B3LYP/aug-cc-pVDZ calculations were carried out and showed 51 neutral and six zwitterionic conformers to be stable in water solution. According to DFT calculations, the conformer equilibrium in the gas phase is dominated by one conformer, whereas the MP2 calculations suggest three PhIS structures to be significantly populated. Comparison of DFT and MP2 energies of all 57 structures stable in water indicates that, in practice, one zwitterionic and one neutral conformer determine the equilibrium in water. Based on the AIM calculations, we found that for the neutral conformers in vacuum and in water, d(H...B) is linearly correlated with Laplacian at the H-bond critical point.     

Conformational analysis of perezone and dihydroperezone using vibrational circular dichroism

The vibrational circular dichroism (VCD) spectra of perezone and dihydroperezone measured from CDCl₃ solutions were quite similar, suggesting analogous conformations for both molecules. Their absolute configurations were confirmed by comparison of the experimental VCD spectrum of each compound with curves generated from theoretical calculations using density functional theory (DFT) at the B3LYP/DGDZVP level of theory taking into account their conformational mobility. Conformational analysis of the 8-(R) enantiomer showed 19 low energy conformers in a 2.4 kcal/mol energy range, while for 8-(R), with the saturated side alkyl chain, 34 conformers were considered in the first 2 kcal/mol. Initial analyses were carried out using a Monte Carlo searching with the MMFF94 molecular mechanics force field, all MMFF94 conformers were geometrically optimized using DFT at the B3LYP/6-31G(d) level of theory, followed by reoptimization and calculations of their vibrational frequencies at the B3LYP/DGDZVP level. Good agreement between the theoretical 8-(R) enantiomers and experimental VCD curves were observed for both.     

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.     

Conformational study of erythritol and threitol in the gas state by density functional theory calculations

Density functional theory calculations using the B3LYP functional and the 6-311++G(d,p) basis set were carried out on the isolated molecules of erythritol and L-threitol. For the meso isomer, a relatively large number of conformers have to be considered to describe the gas state structure. The lowest energy conformer is characterized by the establishment of a strong intramolecular H-bond between the two terminal hydroxyl groups, giving rise to a seven-membered ring and two additional weaker H-bonds between vicinal OH groups. In the case of L-threitol, two conformers are predominant in the gas state, and both are stabilized by the formation of a cyclic system of four intramolecular hydrogen bonds involving all OH groups. The conformational stability in both diastereomers is discussed in terms of the electronic energy and of the Gibbs energy. The weighted mean enthalpy of both diastereomers in the gas state at 298.15 K was obtained from the thermodynamic data and Boltzmann populations of the low-energy conformers.     

Potential transition-state analogs for glycosyltransferases. Design and DFT calculations of conformational behavior

The structure of a previously calculated transition state (TS) was used to design the [tetrahydro-2-(methylthio)furan-2-yl]methyl phosphate dianion (1) as a new scaffold for transition-state analogs of reactions catalyzed by the inverting glycosyltransferases. This scaffold contains relevant features of the donor and acceptor and represents a new type of potential inhibitors for these enzymes. Available conformational space of 1 was explored using DFT quantum chemical methods by means of two-dimensional potential-energy maps calculated as a function of Phi, Psi, and omega dihedral angles at the B3LYP/6-31+G* level. The calculated potential energy surfaces revealed the existence of several low-energy domains. Structures from these regions were refined at the 6-311++G** level and led to 14 conformers. The stability of conformers is influenced by their environment, and in aqueous solution two conformers dominate the equilibrium. A superposition of calculated conformers with the predicted TS structure revealed that the preferred conformers in solution nicely mimic structural features of the TS. These results imply that 1 has structural properties required to mimic the TS and therefore can be used as a scaffold for further development of TS-analog inhibitors for retaining glycosyltransferases.     

Spectroscopic investigation on glutamic acid by Coulomb-attenuating and double hybrid density functional theory methods

This study deals with the identification of glutamic acid by means of quantum chemical approach. FT-IR, FT-Raman and UV–vis spectra were recorded in the region 4000–400, 4000–50 cm^? 1 and 200–600 nm, respectively. CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p) calculations were performed to obtain the optimised molecular structures, vibrational frequencies and corresponding vibrational assignment, thermodynamic properties and natural bonding orbital (NBO) analysis. The results show that the obtained optimised geometric parameters (bond lengths, bond angles and bond dihedrals) and vibrational frequencies were found to be in good agreement with the experimental results. The calculations of the electronic spectra were compared with the experimental ones. Furthermore, highest occupied molecular orbital and lowest unoccupied molecular orbital analyses and UV–vis spectral analysis were also performed to determine the energy band gaps and transition states. NBO analysis, calculated using density functional theory methods (CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p)), was induced to find inter-molecular atoms. ¹³C and ¹H NMR isotropic chemical shifts were calculated and the assignments made were compared with the ChemDraw Ultra values.     

Complete conformational family of 2',3'-didehydro-2',3'-dideoxyguanosine: quantum chemical and electron density topological study

Comprehensive conformational analysis of the biologically active nucleoside 2',3'-didehydro-2',3'-dideoxyaguanosine (d4G) has been performed at the MP2/6-311++G(d,p)//DFT B3LYP/6-31G(d,p) level of theory. The energetic, geometrical and polar characteristics of twenty d4G conformers as well as their conformational equilibrium were investigated. The electron density topological analysis allowed us to establish that the d4G molecule is stabilized by nine types of intramolecular interactions: O5'H...N3, O5'H...C8, C8H...O5', C2'H...N3, C5'H1...N3, C5'H2...N3, C8H...H1C5', C8H...H2'C5' and N2H1...O5'. The obtained results of conformational analysis permit us to think that d4G may be a terminator of the DNA chain synthesis in the 5'-3' direction. Thus it can be inferred that d4G competes with canonical 2'-deoxyaguanosine in binding an active site of the corresponding enzyme.     

Electronic structure simulations of 2,6-dimethyl-2,5-heptadien-4-one by FTIR,FT-Raman,NMR, UV–vis,NBO and density functional theory

Density functional theory (DFT) (B3LYP and B3PW91) calculations have been carried out for 2,6-dimethyl-2,5-heptadien-4-one (DMHD4O) using 6–311++ G** basis set. Complete vibrational assignment and analysis of the fundamental modes of the compound were carried out from the FTIR and FT-Raman spectral data. The theoretical electronic absorption has been calculated by using time-dependent DFT (TD-DFT) methods and compared with the experimental spectra. The theoretically computed Frontier energy gaps and TD-DFT calculations are in good agreement with the experimental UV–vis spectral absorption. The chemical hardness measured from the Frontier molecular orbital energies of DMHD4O is 0.0693 eV. Electronic stability of the compound arising from hyperconjugative interactions and charge delocalisation were also investigated based on the natural bond orbital (NBO) analysis. Effective stabilisation energy E ⁽²⁾ associated with the interactions of the π and the lone pair of electrons was determined by the NBO analysis. ¹³C and ¹H NMR chemical shifts of the compound have been calculated by means of Gauge-Invariant Atomic Orbital using B3LYP/6–311++ G** method. The partial ionic character of the carbonyl group due to resonance render a partially positive charge to the carbonyl carbon, and thus C4 chemical shift lie in the very downfield 191.6 ppm. Comparison between the experimental and the theoretical results indicates that B3LYP method is able to provide satisfactory results for predicting vibrational, electronic and NMR properties.