Inducing high exo selectivity in Diels-Alder reaction by dimethylborane substituent: a DFT study.

In this work, the role of Lewis acid-base (LAB) interaction on the stereoselectivity of the Diels-Alder (DA) reaction has been studied by DFT in gas and solution (dichloromethane) phases. The calculations were performed at the B3LYP/6-311G++ (d, p) level. Two different series of DA reactions were investigated: (1)-three mono-substituted cyclopentadienes + dimethyl(vinyl)borane; (2)-five α,ß-unsaturated carbonyl compounds + cyclopenta-2,4-dien-1-yldimethylborane. The reacting diene and dienophile pairs were chosen to restrict LAB interaction to the exo reaction pathway. It was found that in some of the examined cases, the favorable LAB interaction is so strong that it can lead to a completely exo-selective DA reaction. Furthermore, a novel multistep synthetic method was hypothesized for preparing exo cycloadduct with near 100% stereoselectivity. Our results can open up new avenues toward the rational design of exo-selective DA reactions for synthesizing novel bioorganic compounds.

Design of ionic liquids containing glucose and choline as drug carriers, finding the link between QM and MD studies.

Designing drug delivery systems for therapeutic compounds whose receptors are located in the cytosol of cells is challenging as a bilayer cell membrane is negatively charged. The newly designed drug delivery systems should assist the mentioned drugs in passing the membrane barriers and achieving their targets. This study concentrated on developing novel ionic liquids (ILs) that interact effectively with cell membranes. These ILs are based on glucose-containing choline and are expected to be non-toxic. The binding energies of the known pharmaceutically active ionic liquids were calculated at the B3LYP/6-311++G(d,p) level in the gas phase and compared with those of our newly designed carbohydrate-based ionic liquids. Subsequently, we employed MD simulations to obtain information about the interactions of these known and designed ILs with the cell membrane. In our approach, we adopted QM and MD studies and illustrated that there could be a link between the QM and MD results.

Design of amino acid- and carbohydrate-based anticancer drugs to inhibit polymerase η.

DNA polymerase η (polη) is of significant value for designing new families of anticancer drugs. This protein takes a role in many stages of the cell cycle, including DNA replication, translesion DNA synthesis, and the repairing process of DNA. According to many studies, a high level of expression of polη in most cases has been associated with low rates of patients' survival, regardless of considering the stage of tumor cells. Thus, the design of new drugs with fewer side effects to inhibit polη in cancerous cells has attracted attention in recent years. This project aims to design and explore the alternative inhibitors for polη, which are based on carbohydrates and amino acids. In terms of physicochemical properties, they are similar to the traditional anticancer drugs such as Cytarabine (cytosine arabinose). These alternative inhibitors are supposed to disrupt the DNA replication process in cancerous cells and prevent the tumor cells from mitosis. These newly designed structures, which are based on natural products, are expected to be non-toxic and to have the same chemotherapeutic impact as the traditional agents. The combinatorial use of quantum mechanics studies and molecular dynamic simulation has enabled us to precisely predict the inhibition mechanism of the newly designed structure, which is based on carbohydrates and amino acids, and compare it with that of the traditional chemotherapeutic drugs such as Cytarabine. Our results suggest that the inhibitors containing the natural building blocks of amino acid and carbohydrate could be considered alternative drugs for Cytarabine to block polη.

New pathways of stability for NHCs derived from azole, di-azole, n-tetrazole, and ab-tetrazole, by DFT.

We have investigated the pathways of stability for NHCs derived from azole, di-azole, n-tetrazole, and ab-tetrazole (1a, 2a, 3a, and 4a, respectively), at the M06/6-311++G** level of theory. Optimization and vibrational frequency calculations of ground states (GS) and transition states (TS) are performed to identify Gibbs free energies and nature of stationary points, respectively. Two possible pathways of stability for 1a-4a are compared and contrasted which entail dimerization through hydrogen bonding (HB) and covalent bonding (CB). The CB pathway comprises head to head (HH) and head to tail (HT) dimerizations. Plausible reaction profiles are illustrated for 1a-4a along with the mechanism of each dimerization. Structures 1a-3a show one possibility for HB while 4a represents two possibilities. Structures 1a and 4a display two HH dimers while 2a and 3a show one. Structures 1a-4a undergo HT dimerizations to yield three possible dimers which include trans, cis, and [2+3] isomers. Interestingly, for all 1a-4a, HB dimerization turns out as the most favorable stability pathway for showing no barrier of reaction. Structures 4b and 4c indicate the highest stability with respect to their initial 4a compared to remaining HB dimers 1b-3b. In addition, the 1,2-H shift appears as a possible rearrangement for 1a-4a to yield their corresponding tautomers (1i, 2h, 3h, and 4k, respectively). The reaction profile of this rearrangement indicates that 1a-4a favor HB dimerization pathway more than 1,2-H shift, in terms of kinetic and thermodynamic. Graphical Abstract.

Comparison of acidity and metal ion affinity of D-Glucosamine and N-acetyl-D-glucosamine, a DFT study.

The derivatives of glucose such as glucosamine (ß-D-GlcN) and N-acetyl-D-ß-glucosamine (GlcNAc) are significant in several biological systems. D-GlcN has been used widely to treat osteoarthritis in humans and animal models as well as GlcNAc has been proposed as a treatment for autoimmune diseases. The DFT/B3LYP/6-311++G (d,p) method as well as QTAIM and NBO analyses were used to the acidity values of D-GlcN and GlcNAc sugars and their complexes with alkali ions in the gas phase. The Li+, Na+ and K+ prefer bi-dentate chelate in these complexes. The computed results indicate that metal ion affinity (MIA) in GlcNAc is higher than that in D-GlcN. There are direct correlations between the MIA values of D-GlcN and GlcNAc sugars and the atomic numbers of Li, Na, and K. The calculated acidity values for GlcNAc at C2-NH and C6-HO6 sites are 331 and 333 kcal mol-1, respectively. Whereas the calculated acidity values for D-GlcN at C2-NH and C6-HO6 sites are 365 and 372 kcal mol-1, respectively. The AIM and NBO analyses indicate the presence of intramolecular H-bonds in GlcNAc sugar in both its neutral form and conjugate base; whereas D-GlcN indicates intramolecular H-bonds in only its conjugate base.

Phenylcyclopropane Energetics and Characterization of Its Conjugate Base: Phenyl Substituent Effects and the C-H Bond Dissociation Energy of Cyclopropane.

The α-C-H bond dissociation energy (BDE) of phenylcyclopropane (1) was experimentally determined using Hess' law. An equilibrium acidity determination of 1 afforded ΔH°acid = 389.1 ± 0.8 kcal mol-1, and isotopic labeling established that the α-position of the three-membered ring is the favored deprotonation site. Interestingly, the structure of the base proved to be a key factor in correctly determining the proper ionization site (i.e., secondary amide ions are needed, and primary ones and OH- lead to incorrect conclusions since they scramble the deuterium label). An experimental measurement of the electron affinity of 1-phenylcyclopropyl radical (EA = 17.5 ± 2.8 kcal mol-1) was combined with the ionization energy of hydrogen (313.6 kcal mol-1) to afford BDE = 93.0 ± 2.9 kcal mol-1. This enabled the effect of the phenyl substituent to be evaluated and compared to other situations where it is attached to an sp3- or sp2-hybridized carbon center. M06-2X, CCSD(T), G4, and W1BD computations were also carried out, and a revised C-H BDE for cyclopropane of 108.9 ± 1.0 kcal mol-1 is recommended.

Stereoelectronic effects: a simple yet powerful tool to manipulate anion affinity.

Different strategies are employed in designing strong and selective anion receptors but stereoelectronic effects have been largely ignored. In this work, the stereo configuration of a non-interacting ether is found to have a large impact of more than two orders of magnitude on the binding of a rigid diol with tetrabutylammonium chloride in acetonitrile-d3. A favorable carbon-oxygen dipole and an intramolecular C-HOH hydrogen bond in an equatorially substituted ether is found to be energetically more important than a stabilizing hydrogen bond in the corresponding axially oriented alcohol. IR spectroscopy is also used to probe the structures of the bound complexes and several binding motifs are identified.

Power of a remote hydrogen bond donor: anion recognition and structural consequences revealed by IR spectroscopy.

Natural and synthetic anion receptors are extensively employed, but the structures of their bound complexes are difficult to determine in the liquid phase. Infrared spectroscopy is used in this work to characterize the solution structures of bound anion receptors for the first time, and surprisingly only two of three hydroxyl groups of the neutral aliphatic triols are found to directly interact with Cl(­). The binding constants of these triols with zero to three CF3 groups were measured in a polar environment, and KCD3CN(Cl(­)) = 1.1 × 10(6) M(­1) for the tris(trifluoromethyl) derivative. This is a remarkably large value, and high selectivity with respect to interfering anions such as, Br(­), NO3(­) and NCS(­) is also displayed. The effects of the third "noninteracting" hydroxyl groups on the structures and binding constants were also explored, and surprisingly they are as large or larger than the OH substituents that hydrogen bond to Cl(­). That is, a remote hydroxyl group can play a larger role in binding than two OH substituents that directly interact with an anionic center.

Electrostatically defying cation-cation clusters: can likes attract in a low-polarity environment?

Like-charge ion pairing is commonly observed in protein structures and plays a significant role in biochemical processes. Density functional calculations combined with the conductor-like polarizable continuum model were employed to study the formation possibilities of doubly charged noncovalently linked complexes of a series of model compounds and amino acids in the gas phase and in solution. Hydrogen bond interactions were found to offset the Coulombic repulsion such that cation-cation clusters are minima on the potential energy surfaces and neither counterions nor solvent molecules are needed to hold them together. In the gas phase the dissociation energies are exothermic, and the separation barriers span from 1.7 to 15.6 kcal mol(-1). Liquid-phase computations indicate that the separation enthalpies of the cation-cation complexes become endothermic in water and nonpolar solvents with dielectric constants of ≥7 (i.e., the value for THF). These results reveal that electrostatically defying noncovalent complexes of like-charged ions can overcome their Coulombic repulsion even in low-polarity environments.

Comparison of gas phase intrinsic properties of cytosine and thymine nucleobases with their O-alkyl adducts: different hydrogen bonding preferences for thymine versus O-alkyl thymine.

In recent years, there has been increasing interest in damaged DNA and RNA nucleobases. These damaged nucleobases can cause DNA mutation, resulting in various diseases such as cancer. Alkylating agents are mutagenic and carcinogenic in a variety of prokaryotic and eukaryotic organisms. The present study employs density functional theory (DFT/B3LYP) with the 6-311++G(d,p) basis set to investigate the effect of chemical damage in O-alkyl pyrimidines such as O(4)-methylthymine, O(2)-methylcytosine and O(2)-methylthymine. We compared the intrinsic properties, such as proton affinities, gas phase acidities, equilibrium tautomerization and nucleobase pair's hydrogen bonding properties, of these molecules with those in the normal nucleobases thymine and cytosine. The results are of interest for chemical reasons and also possibly for biological purposes since biological media can be quite non-polar. Furthermore, we found that N1-H of O(4)-methylthymine is less acidic than N1-H of thymine, suggesting that alkyl DNA glycosylase enzyme cannot discriminate this damaged nucleobase from a normal thymine nucleobase. This result indicates that the conjugated base anion of O(4)-methylthymine would be a worse leaving group and O(4)-methylthymine is repaired in genome by demethylation rather than enzyme-catalyzed excision at N1.

Interactions of coinage metal clusters with histidine and their effects on histidine acidity; theoretical investigation.

Understanding the nature of interaction between metal nanoparticles and biomolecules such as amino acids is important in the development and design of biosensors. In this paper, binding of M(3) clusters (M = Au, Ag and Cu) with neutral and anionic forms of histidine amino acid was studied using density functional theory (DFT-B3LYP). It was found that the interaction of histidine with M(3) clusters is governed by two major bonding factors: (a) the anchoring N-M and O-M bonds and (b) the nonconventional N-H···M and O-H···M hydrogen bonds. The nature of these chemical bonds has been investigated based on quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses. In the next step, the effects of Au, Ag and Cu metal clusters on the gas-phase acidity of weak organic acid (histidine) have been explored. The acidity of isolated histidine was compared with the acidity of its Au(3)-, Ag(3)- and Cu(3)-complexed species. Results indicate that upon complexation with M(3) clusters (at 298 K), the gas-phase acidity (GPA) of histidine varies from 339.5 to 312.3, 315.0, and 313.7 kcal mol(-1) for Au(3)-, Ag(3)- and Cu(3)-His complexes, respectively (i.e., its dissociation becomes much less endothermic). These values indicate that a weak organic acid can be converted to a super acid when it is complexed with metal clusters. Also, in order to investigate the acidity value of the imidazole moiety in histidine, histidine methyl ester (His-OMe) was selected. Similarly, the acidity of this compound was compared with the acidity of their Au(3), Ag(3) and Cu(3)-complexed species. After complexation with M(3) clusters at 298 K, the gas-phase acidity (GPA) of His-OMe varies from 333.0 to 280.0, 304.2 and 291.5 kcal mol(-1), respectively. Moreover, pK(a) values were determined in water for isolated and complexed species of His and His-OMe. The resulting pK(a) values were found to decrease upon complexation with M(3) clusters.

Effect of hydrogen bonds on pKa values: importance of networking.

The pK(a) of an acyclic aliphatic heptaol ((HOCH(2)CH(2)CH(OH)CH(2))(3)COH) was measured in DMSO, and its gas-phase acidity is reported as well. This tertiary alcohol was found to be 10(21) times more acidic than tert-butyl alcohol in DMSO and an order of magnitude more acidic than acetic acid (i.e., pK(a) = 11.4 vs 12.3). This can be attributed to a 21.9 kcal mol(-1) stabilization of the charged oxygen center in the conjugate base by three hydrogen bonds and another 6.3 kcal mol(-1) stabilization resulting from an additional three hydrogen bonds between the uncharged primary and secondary hydroxyl groups. Charge delocalization by both the first and second solvation shells may be used to facilitate enzymatic reactions. Acidity constants of a series of polyols were also computed, and the combination of hydrogen-bonding and electron-withdrawing substituents was found to afford acids that are predicted to be extremely acidic in DMSO (i.e., pK(a) < 0). These hydrogen bond enhanced acids represent an attractive class of Brønsted acid catalysts.

Ionic liquid based on α-amino acid anion and N7,N9-dimethylguaninium cation ([dMG][AA]): theoretical study on the structure and electronic properties.

The interactions between five amino acid based anions ([AA](-) (AA = Gly, Phe, His, Try, and Tyr)) and N7,N9-dimethylguaninium cation ([dMG](+)) have been investigated by the hybrid density functional theory method B3LYP together with the basis set 6-311++G(d,p). The calculated interaction energy was found to decrease in magnitude with increasing side-chain length in the amino acid anion. The interaction between the [dMG](+) cation and [AA](-) anion in the most stable configurations of ion pairs is a hydrogen bonding interaction. These hydrogen bonds (H bonds) were analyzed by the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis. Finally, several correlations between electron densities in bond critical points of hydrogen bonds and interaction energy as well as vibrational frequencies in the most stable configurations of ion pairs have been checked.

Interactions of glutathione tripeptide with gold cluster: influence of intramolecular hydrogen bond on complexation behavior.

Understanding the nature of the interaction between metal nanoparticles and biomolecules has been important in the development and design of sensors. In this paper, structural, electronic, and bonding properties of the neutral and anionic forms of glutathione tripeptide (GSH) complexes with a Au(3) cluster were studied using the DFT-B3LYP with 6-31+G**-LANL2DZ mixed basis set. Binding of glutathione with the gold cluster is governed by two different kinds of interactions: Au-X (X = N, O, and S) anchoring bond and Au···H-X nonconventional hydrogen bonding. The influence of the intramolecular hydrogen bonding of glutathione on the interaction of this peptide with the gold cluster has been investigated. To gain insight on the role of intramolecular hydrogen bonding on Au-GSH interaction, we compared interaction energies of Au-GSH complexes with those of cystein and glycine components. Our results demonstrated that, in spite of the ability of cystein to form highly stable metal-sulfide interaction, complexation behavior of glutathione is governed by its intramolecular backbone hydrogen bonding. The quantum theory of atom in molecule (QTAIM) and natural bond orbital analysis (NBO) have also been applied to interpret the nature of interactions in Au-GSH complexes. Finally, conformational flexibility of glutathione during complexation with the Au(3) cluster was investigated by means of monitoring Ramachandran angles.

Mechanisms and kinetics of thiotepa and tepa hydrolysis: DFT study.

N,N',N″-triethylenethiophosphoramide (Thiotepa) and its oxo analogue (Tepa) as the major metabolite are trifunctional alkylating agents with a broad spectrum of antitumor activity. In vivo and vitro studies show alkylation of DNA by Thiotepa and Tepa can follow two pathways, but it remains unclear which pathway represents the precise mechanism of action. In pathway 1, these agents are capable of forming cross-links with DNA molecules via two different mechanisms. In the first mechanism, the ring opening reaction is initiated by protonating the aziridine, which then becomes the primary target of nucleophilic attack by the N7-Guanine. The second one is a direct nucleophilic ring opening of aziridyl group. Thiotepa and Tepa in pathway 2, act as a cell penetrating carrier for aziridine, which is released via hydrolysis. The released aziridine can form a cross-link with N7-Guanine. In this study, we calculated the activation free energy and kinetic rate constant for hydrolysis of these agents and explored interaction of aziridine with Guanine to predict the most probable mechanism by applying density functional theory (DFT) using B3LYP method. In addition, solvent effect was introduced using the conductor-like polarizable continuum model (CPCM) in water, THF and diethylether. Hyperconjugation stabilization factors that have an effect on stability of generated transition state were investigated by natural bond order (NBO) analysis. Furthermore, quantum theory of atoms in molecules (QTAIM) analysis was performed to extract the bond critical points (BCP) properties, because the electron densities can be considered as a good description of the strength of different types of interactions.

Experimental and computational bridgehead C-H bond dissociation enthalpies.

Bridgehead C-H bond dissociation enthalpies of 105.7 ± 2.0, 102.9 ± 1.7, and 102.4 ± 1.9 kcal mol(-1) for bicyclo[2.2.1]heptane, bicyclo[2.2.2]octane, and adamantane, respectively, were determined in the gas phase by making use of a thermodynamic cycle (i.e., BDE(R-H) = ΔH°(acid)(H-X) - IE(H(·)) + EA(X(·))). These results are in good accord with high-level G3 theory calculations, and the experimental values along with G3 predictions for bicyclo[1.1.1]pentane, bicyclo[2.1.1]hexane, bicyclo[3.1.1]heptane, and bicyclo[4.2.1]nonane were found to correlate with the flexibility of the ring system. Rare examples of alkyl anions in the gas phase are also provided.

H/D exchange kinetics: experimental evidence for formation of different b fragment ion conformers/isomers during the gas-phase peptide sequencing.

Electrospray ionization (ESI) Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) combined with H/D exchange reactions was utilized to explore the existence of different b(5)(+) and b(4)(+) fragment ion conformers/isomers of hexapeptide WHWLQL in the gas phase. Distinct H/D exchange trends for protonated WHWLQL ([M + H](+)) and its b(5)(+) and b(4)(+) fragment ions (with ND(3)) were observed. Isolated (12)C(all) isotopomers of both b(5)(+) and b(4)(+) fragment ions yielded bimodal distributions of H/D exchanged product ions. The H/D exchange reaction kinetics also confirmed that b(5)(+) and b(4)(+) fragment ions exist as combination of slow-exchanging ("s") and fast-exchanging ("f") species. The calculated rate constant for the first labile hydrogen exchange of [M + H](+) (k([M + H](+)) = 3.80 +/- 0.7 x 10(-10) cm(3) mol(-1) s(-1)) was approximately 30 and approximately 5 times greater than those for the "s" and "f" species of b(5)(+), respectively. Data from H/D exchange of isolated "s" species at longer ND(3) reaction times confirmed the existence of different conformers or isomers for b(5)(+) fragment ions. The sustained off-resonance irradiation collision-activated dissociation (SORI-CAD) of WHWLQL combined with the H/D exchange reactions indicate that "s" and "f" species of b(5)(+) and b(4)(+) fragment ions can be produced in the ICR cell as well as the ESI source. The significance of these observations for detailed understanding of protein sequencing and ion fragmentation pathways is discussed.

Single-centered hydrogen-bonded enhanced acidity (SHEA) acids: a new class of Brønsted acids.

Hydrogen bonds are the dominant motif for organizing the three-dimensional structures of biomolecules such as carbohydrates, nucleic acids, and proteins, and serve as templates for proton transfer reactions. Computations, gas-phase acidity measurements, and pK(a) determinations in dimethyl sulfoxide on a series of polyols indicate that multiple hydrogen bonds to a single charged center lead to greatly enhanced acidities. A new class of Brønsted acids, consequently, is proposed.

DFT study of the interaction of cytidine and 2'-deoxycytidine with Li+, Na+, and K+: effects of metal cationization on sugar puckering and stability of the N-glycosidic bond.

Density functional theory (DFT) calculations were performed at the B3LYP level with a 6-311++G(d,p) basis set to systematically explore the geometrical multiplicity and binding strength for complexes formed by Li(+), Na(+), and K(+) with cytidine and 2'-deoxycytidine. All computational studies indicate that the metal ion affinity (MIA) decreases from Li(+) to Na(+) and K(+) for cytosine nucleosides. For example, for cytidine the affinity for the above metal ions are 79.5, 55.2, and 41.8 and for 2'-deoxycytidine, 82.8, 57.4, and 42.2 kcal/mol, respectively. It is also interesting to mention that linear correlations between calculated MIA values and the atomic numbers (Z) of the above metal ions were found. The influence of metal cationization on the coordination modes and the strength of the N-glycosidic bond in cytosine nucleosides have been studied. In all cases, the N1-C1' bond distance changes upon introducing a positive charge in the nucleosides. It has been found that metal binding significantly changes the values of the phase angle of pseudorotation P in the sugar unit of these nucleosides. With respect to the sugar ring, metal binding changes the values of the glycosyl torsion angle and sugar ring conformation. The present calculations in the gas phase provide the first clues on the intrinsic chemistry of these systems and may be of value for studies of the influence of metal cations on the conformational behavior and function of nucleic acids.

Comparison of thermochemistry of aspartame (artificial sweetener) and glucose.

We have compared the gas phase thermochemical properties of aspartame (artificial sweetener) and alpha- and beta-glucose. These parameters include metal ion affinities with Li(+)-, Na(+)-, K(+)-, Mg(+2)-, Ca(+2)-, Fe(+2)-, Zn(+2)-ions, and chloride ion affinity by using DFT calculations. For example, for aspartame, the affinity values for the above described metal ions are, respectively, 86.5, 63.2, 44.2, 255.4, 178.4, 235.4, and 300.4, and for beta-glucose are 65.2, 47.3 32.9, 212.9, 140.2, 190.1, and 250.0 kcal mol(-1), respectively. The study shows differences between the intrinsic chemistry of aspartame and glucose.