A comparison of the behavior of functional/basis set combinations for hydrogen-bonding in the water dimer with emphasis on basis set superposition error

Optical band gap of cross-linked, curved, and radical polyaromatic hydrocarbons

The impact of cross-linking, curvature, and radical character on the optical band gap of polyaromatic hydrocarbons has been investigated.

Tuning receptors for the encapsulation of beryllium2+

A method for evaluating second sphere H-bonding interactions as a means to tune specific characteristics in coordinating ligands has been used to evaluate the binding strength of a series of beryllium complexes.

Energetic, Topological and Electric Field Analyses of Cation-Cation Nucleic Acid Interactions in Watson-Crick Disposition

A theoretical study of the effect of the diprotonation on the nucleic acid bases (A : U, A : T and G : C) in Watson-Crick conformation has been carried out by means of DFT computational methods in vacuum. In addition, the corresponding neutral and monoprotonated binary complexes have been considered. Most of the diprotonated species studied are stable, even though the binding energy is positive due to the overall repulsive electrostatic term. Local electrostatic attractive forces in the regions of hydrogen bonds (HBs) are responsible for equilibrium geometries, as shown by the electric field lines connecting the electrophilic and nucleophilic sites involved in the HB interactions. Secondary electrostatic effects also affect the assembling of the nucleic acid complexes in either neutral or cationic form. In particular, the electric field lines flowing from electrophilic sites in one base to nucleophilic sites in the other reinforce the linking between them. Hence, when the nucleophilic site concerns the free lone pair of the heteroatom involved in the HB interaction as acceptor, the HB distance shortens. However, if the free lone pair of the HB acceptor interacts with an electrophilic site in the same molecule, the HB distance elongates, weakening the HB interaction. The topological analysis of the electron density distribution in HB regions indicates that neutral, monoprotonated and diprotonated complexes show no differences in the nature of their HB's.

Ring-Size Effects on the Stability and Spectral Shifts of Hydrogen Bonded Cyclic Ethers Complexes

In order to investigate ring-size effects on the stability and spectral shifts of hydrogen bonded cyclic ethers complexes, the strength of hydrogen bonds in gas phase complexes formed between 2,2,2-trifluoroethanol (TFE) and selected cyclic ethers were examined using FTIR spectroscopy. TFE was chosen as hydrogen bond donor in these complexes, while trimethylene oxide (TMO), tetrahydrofuran (THF) and tetrahydropyran (THP) were selected as hydrogen bond acceptors. Comparable OH-stretching red shifts were observed in the three kinds of complexes. The difference of red shifts is so small (<7 cm⁻¹) for TFE−TMO/THF/THP complexes that one can conclude that their stabilities and the strength of the hydrogen bonds are nearly similar and do not show any marked dependence with the ring size of the hydrogen bond acceptor. The equilibrium constants for the complexation were determined, and atoms-in-molecules (AIM) and natural bond orbital (NBO) analyses were performed to further investigate the intermolecular interactions. Regardless of the ring size, hydrogen bonds in the complexes showed similar strength, in agreement with the observed OH-stretching red shifts.

Template-free H-dimer and H-aggregate formation by dimeric carbocyanine dyes

Dimeric cyanine dyes self-assemble into H-dimers and H-aggregates, which experience de-aggregation with specific biomolecules.

Why lamivudine assembles into double-stranded helices in crystals: salt heterosynthon versus base-pairing homosynthon

Here we were interested in obtaining a better understanding of the competition between the salt heterosynthon and the base-pairing homosynthon formed by the anti-HIV drug lamivudine in the presence of strong acids.

Fingerprints of Both Watson-Crick and Hoogsteen Isomers of the Isolated (Cytosine-Guanine)H+ Pair

 Gas phase protonated guanine-cytosine (CGH⁺) pair was generated using an electrospray ionization source from solutions at two different pH (5.8 and 3.2). Consistent evidence from MS/MS fragmentation patterns and differential ion mobility spectra (DIMS) point toward the presence of two isomers of the CGH⁺ pair, whose relative populations depend strongly on the pH of the solution. Gas phase infrared multiphoton dissociation (IRMPD) spectroscopy in the 900-1900 cm^-1 spectral range further confirms that the Watson-Crick isomer is preferentially produced (91%) at pH = 5.8, while the Hoogsteen isomer predominates (66%) at pH = 3.2). These fingerprint signatures are expected to be useful for the development of new analytical methodologies and to trigger isomer selective photochemical studies of protonated DNA base pairs.

Testing the CP-correction procedure with different DFT methods on H-bonding complexes of κ-carrabiose with water molecules

Interaction of water molecules with κ-carrabiose disaccharide, within three H-bonding complexes, was investigated. Particular interest was focused on the way with which the BSSE correction has to be performed. Two strategies were used, either performing BSSE correction during or after optimization. For this aim, several DFT-functionals (hybrid GGA and hybrid meta-GGA) and 6-31 + G* basis set were considered. The results demonstrated the uselessness of including of BSSE-CP correction during optimization for all complexes. From a structural point of view, a proper H-bonding description was obtained using the PBE0 functional for all complexes. The basis set effect on the BSSE using B3LYP functional was also investigated. The reliability of B3LYP/6-31 + G** and B3LYP/6-31++G** models for the complexes involving one or two water molecules was reported while the use of B3LYP/6-311 + G** or B3LYP/6-311++G** levels was shown to be more appropriate for larger complexes equivalent to that involving three water molecules. CP-corrected interaction energies were demonstrated to be closer to CBS-4 M interaction energies than the uncorrected ones. Graphical abstract Functional and basis set effects on BSSE.

Synthon trends according to acid strength and geometry in salts of N-heterocyclic bases

The hierarchy and robustness of homosynthons and heterosynthons formed by N-heterocyclic bases were assessed experimentally in salts of aminopyrazine (ampyz) andtrans-1,2-bis(4-pyridyl)ethane (BPE) with common strong, moderate and weak acids, and theoretically at the M06-2X/6-31+G** level of theory.

Comparison of tetrel bonds in neutral and protonated complexes of pyridineTF3and furanTF3(T = C, Si, and Ge) with NH3

Protonation not only changes the primary interaction mode between α/β-furanCF₃/p-PyCF₃and NH₃but also prominently enhances the strength of the Si/Ge⋯N tetrel bond.

Predicting solvent-water partitioning of charged organic species using quantum-chemically estimated Abraham pp-LFER solute parameters

Methods for obtaining accurate predictions of solvent-water partitioning for neutral organic chemicals (e.g., K_ow) are well-established. However, methods that provide comparable accuracy are not available for predicting the solvent-water partitioning of ionic species. Previous methods for addressing charge contributions to solvent-water partitioning rely on charged solute descriptors which are obtained from regressions to neutral species descriptors as well as charged descriptors which are specific to unique charge-functionalities and structural moieties. This paper presents a method for obtaining Abraham poly-parameter linear free energy relationship (pp-LFER) descriptors using quantum chemical calculations and molecular structure, only. The method utilizes a large number of solvent-water systems to overcome large errors in individual quantum chemical computations of ionic solvent-water partition coefficients. The result is a single set of quantum-chemically estimated Abraham solute parameters (QCAP) which are solvent-independent, and can be used to predict the solvent-water partitioning of ionic species. Predictions of solvent-water partition coefficients for ionic species using quantum-chemically estimated Abraham parameters (QCAPs) are shown to provide improved accuracy compared over both existing Absolv-estimated Abraham solute parameters (AAP) as well as direct a priori quantum chemical (QC) calculations for partitioning of anionic solutes in 4 organic solvent-water systems (RMS = 0.740, 2.48 and 0.426 for the Absolv, QC and QCAP methods, respectively). For quaternary amine cations in the octanol-water system the RMS errors of the solvent-water partition coefficients were larger and similar between the two Abraham models (RMSE = 0.997 and 1.16, for the AAP and QCAP methods, respectively). Both methods showed significant improvement over direct QC calculations (RMSE = 2.82).

DFT simulation towards evaluation the molecular structure and properties of the heterogeneous C16Mg8O8 nano-cage as selective nano-sensor for H2 and N2 gases

Adsorption of hydrogen (H₂) and nitrogen (N₂) molecules was analyzed on a new fullerene-like C₁₆Mg₈O₈ nano-cage, composed of magnesium, oxygen, and carbon, using density functional theory. A detailed analysis of the energy, geometry, and electronic structure of various H₂ and N₂ adsorptions on the cluster surface was performed. The adsorption energies of H₂ and N₂ were estimated to ranging from -0.16 to -0.52eV, respectively. The most stable adsorption configurations were those in which the H or N atoms of the adsorbates were located near the Mg atom of the cluster surface at different sides. It was found that the heterogeneous C₁₆Mg₈O₈ nano-cluster selectively act against the H₂ and N₂ gaseous molecules. The electrical conductivity of the cluster, arising from HOMO/LUMO energy gap, was more sensitive to N₂ gaseous molecule rather than H₂ one, indicating that the heterogeneous C₁₆Mg₈O₈ nano-cage may be potential nano-sensor for N₂ molecule. These findings were specified by analyzing the characteristics in the electron density of states (DOS).

Hydrogen bonding in the hydroxysulfinyl radical-formic acid-water system: A theoretical study

Quantum chemical methods have been employed to evaluate the possible configurations of the 1:1 and 1:2 HOSO-formic acid complexes and 1:1:1 HOSO-formic acid-water complexes. The first type of complex involves two H bonds, while the other two types comprise three H bonds in a ring. The complexes are relatively stable, with CBS-QB3 computed binding energies of 14.3 kcal mol(-1) , 23.4 kcal mol(-1) , and 21.1 kcal mol(-1) for the lowest-energy structures of the 1:1, 1:2, and 1:1:1 complexes, respectively. Complex formations induce a large spectral red-shift and an enhancement of the IR intensity for the H-bonded OH stretching modes relative to those in the parent monomers. TDDFT calculations of the low-lying electronic excited states demonstrate that the complexes are photochemically quite stable in the troposphere. Small spectral shifts in comparison to the free HOSO radical suggest that the radical and the complexes would not be easily distinguishable using standard UV/vis absorption spectroscopy. © 2016 Wiley Periodicals, Inc.

Regioselective Acylation of Diols and Triols: The Cyanide Effect

Central topics of carbohydrate chemistry embrace structural modifications of carbohydrates and oligosaccharide synthesis. Both require regioselectively protected building blocks that are mainly available via indirect multistep procedures. Hence, direct protection methods targeting a specific hydroxy group are demanded. Dual hydrogen bonding will eventually differentiate between differently positioned hydroxy groups. As cyanide is capable of various kinds of hydrogen bonding and as it is a quite strong sterically nondemanding base, regioselective O-acylations should be possible at low temperatures even at sterically congested positions, thus permitting formation and also isolation of the kinetic product. Indeed, 1,2-cis-diols, having an equatorial and an axial hydroxy group, benzoyl cyanide or acetyl cyanide as an acylating agent, and DMAP as a catalyst yield at -78 °C the thermodynamically unfavorable axial O-acylation product; acyl migration is not observed under these conditions. This phenomenon was substantiated with 3,4-O-unproteced galacto- and fucopyranosides and 2,3-O-unprotected mannopyranosides. Even for 3,4,6-O-unprotected galactopyranosides as triols, axial 4-O-acylation is appreciably faster than O-acylation of the primary 6-hydroxy group. The importance of hydrogen bonding for this unusual regioselectivity could be confirmed by NMR studies and DFT calculations, which indicate favorable hydrogen bonding of cyanide to the most acidic axial hydroxy group supported by hydrogen bonding of the equatorial hydroxy group to the axial oxygen. Thus, the "cyanide effect" is due to dual hydrogen bonding of the axial hydroxy group which enhances the nucleophilicity of the respective oxygen atom, permitting an even faster reaction for diols than for mono-ols. In contrast, fluoride as a counterion favors dual hydrogen bonding to both hydroxy groups leading to equatorial O-acylation.

Assessing the accuracy of some popular DFT methods for computing harmonic vibrational frequencies of water clusters

A wide range of density functional theory (DFT) methods (37 altogether), including pure, hybrid, range-separated hybrid, double-hybrid, and dispersion-corrected functionals, have been employed to compute the harmonic vibrational frequencies of eight small water clusters ranging in size from the dimer to four different isomers of the hexamer. These computed harmonic frequencies have been carefully compared to recently published benchmark values that are expected to be very close to the CCSD(T) complete basis set limit. Of the DFT methods examined here, ωB97 and ωB97X are the most consistently accurate, deviating from the reference values by less than 20 cm(-1) on average and never more than 60 cm(-1). The performance of double-hybrid methods including B2PLYP and mPW2-PLYP is only slightly better than more economical approaches, such as the M06-L pure functional and the M06-2X hybrid functional. Additionally, dispersion corrections offer very little improvement in computed frequencies.

A study on mediating the crystallization behavior of PBT through intermolecular hydrogen-bonding

The intermolecular hydrogen-bonding formed between PBT and TDP can efficiently mediate the crystallization behavior of PBT.

Isolated glyoxylic acid-water 1:1 complexes in low temperature argon matrices

The 1:1 hydrogen bonded complexes between glyoxylic acid (GA) and water are studied in low temperature argon matrices. Four different complex structures were found in deposited matrices. The lowest energy conformer (T1) of GA was found to form complex, where the water molecule was attached to the opposite side of the intramolecular hydrogen bond in the molecule (T1B). Interestingly, this complex was estimated to be+8.0 kJ mol(-1) higher in energy than the most stable structure (T1A), where the water is inserted into the internal hydrogen bond, and also found in solid argon but in smaller abundance. For the second-lowest energy conformer of GA (T2), the two lowest-energy complex structures were identified, with the most stable complex structure (T2A) also being the most abundant in the matrices. The difference between experiment and computational energetic order of the two complex structures of the same GA conformer is explained by contributions of deformation energy upon complexation and the effect of the environment. The computed BSSE-corrected interaction energies are for the two most stable complexes of the two GA conformers for T1A and T2A -42.11 and -45.03 kJ mol(-1), respectively, at the CCSD(T)/aug-cc-pVTZ//B3LYP/aug-cc-pVTZ level of theory.

Membrane alternatives in worlds without oxygen: Creation of an azotosome

The lipid bilayer membrane, which is the foundation of life on Earth, is not viable outside of biology based on liquid water. This fact has caused astronomers who seek conditions suitable for life to search for exoplanets within the "habitable zone," the narrow band in which liquid water can exist. However, can cell membranes be created and function at temperatures far below those at which water is a liquid? We take a step toward answering this question by proposing a new type of membrane, composed of small organic nitrogen compounds, that is capable of forming and functioning in liquid methane at cryogenic temperatures. Using molecular simulations, we demonstrate that these membranes in cryogenic solvent have an elasticity equal to that of lipid bilayers in water at room temperature. As a proof of concept, we also demonstrate that stable cryogenic membranes could arise from compounds observed in the atmosphere of Saturn's moon, Titan, known for the existence of seas of liquid methane on its surface.

Non-covalent interactions in ionic liquid ion pairs and ion pair dimers: a quantum chemical calculation analysis

Weak non-covalent interactions were studied by means of QTAIM and NCI approaches in ion pairs and ion pair dimers of 1-alkyl-3-methylimidazolium cations coupled with perfluorinated anions.

Alternative chelator for ⁸⁹Zr radiopharmaceuticals: radiolabeling and evaluation of 3,4,3-(LI-1,2-HOPO)

Zirconium-89 is an effective radionuclide for antibody-based positron emission tomography (PET) imaging because its physical half-life (78.41 h) matches the biological half-life of IgG antibodies. Desferrioxamine (DFO) is currently the preferred chelator for ⁸⁹Zr⁴⁺; however, accumulation of ⁸⁹Zr in the bones of mice suggests that ⁸⁹Zr⁴⁺ is released from DFO in vivo. An improved chelator for ⁸⁹Zr⁴⁺ could eliminate the release of osteophilic ⁸⁹Zr⁴⁺ and lead to a safer PET tracer with reduced background radiation dose. Herein, we present an octadentate chelator 3,4,3-(LI-1,2-HOPO) (or HOPO) as a potentially superior alternative to DFO. The HOPO ligand formed a 1:1 Zr-HOPO complex that was evaluated experimentally and theoretically. The stability of ⁸⁹Zr-HOPO matched or surpassed that of ⁸⁹Zr-DFO in every experiment. In healthy mice, ⁸⁹Zr-HOPO cleared the body rapidly with no signs of demetalation. Ultimately, HOPO has the potential to replace DFO as the chelator of choice for ⁸⁹Zr-based PET imaging agents.

Electronic circular dichroism behavior of chiral Phthiobuzone

Phthiobuzone is a bis(thiosemicarbazone) derivative with a single chiral center which has been used as a racemate in the clinical treatment of herpes and trachoma diseases. In this study, its two enantiomers were prepared from chiral amino acids and their absolute configurations were investigated by electronic circular dichroism (ECD) combined with modern quantum-chemical calculations using time-dependent density functional theory. It was found that solvation changed both the conformational distribution and the ECD spectrum of each conformer. The theoretical ECD spectra of the two enantiomers were in good agreement with the experimentally determined spectra of the corresponding isomers in dimethyl sulfoxide. The ECD behavior of the bis(thiosemicarbazone) chromophore in a chiral environment is also discussed. Our results indicate that ECD spectroscopy may be a useful tool for the stereochemical evaluation of chiral drugs.

Capping amyloid β-sheets of the tau-amyloid structure VQIVYK with hexapeptides designed to arrest growth. An ONIOM and density functional theory study

We present ONIOM calculations using density functional theory (DFT) as the high and AM1 as the medium level that explore the abilities of different hexapeptide sequences to terminate the growth of a model for the tau-amyloid implicated in Alzheimer’s disease. We delineate and explore several design principles (H-bonding in the side chains, using antiparallel interactions on the growing edge of a parallel sheet, using all-d residues to form rippled interactions at the edge of the sheet, and replacing the H-bond donor N–H’s that inhibit further growth) that can be used individually and in combination to design such peptides that will have a greater affinity for binding to the parallel β-sheet of acetyl-VQIVYK-NHCH₃ than the natural sequence and will prevent another strand from binding to the sheet, thus providing a cap to the growing sheet that arrests further growth. We found peptides in which the Q is replaced by an acetyllysine (aK) residue to be particularly promising candidates, particularly if the reverse sequence (KYVIaKV) is used to form an antiparallel interaction with the sheet.

The interactions of phenylalanines in β-sheet-like structures from molecular orbital calculations using density functional theory (DFT), MP2, and CCSD(T) methods

We present density functional theory calculations designed to evaluate the importance of π-stacking interactions to the stability of in-register Phe residues within parallel β-sheets, such as amyloids. We have used a model of a parallel H-bonded tetramer of acetylPheNH2 as a model and both functionals that were specifically designed to incorporate dispersion effects (DFs), as well as, several traditional functionals which have not been so designed. None of the functionals finds a global minimum for the π-stacked conformation, although two of the DFs find this to be a local minimum. The stacked phenyls taken from the optimized geometries calculated for each functional have been evaluated using MP2 and CCSD(T) calculations for comparison. The results suggest that π-stacking does not make an important contribution to the stability of this system and (by implication) to amyloid formation.

Capping parallel β-sheets of acetyl(Ala)6NH2 with an acetyl(Ala)5ProNH2 can arrest the growth of the sheet, suggesting a potential for curtailing amyloid growth. An ONIOM and density functional theory study

We present ONIOM calculations using B3LYP/d95(d,p) as the high level and AM1 as the medium level on parallel β-sheets containing four strands of Ac-AAAAAA-NH2 capped with either Ac-AAPAAA-NH2 or Ac-AAAPAA-NH2. Because Pro can form H-bonds from only one side of the peptide linkage (that containing the C═O H-bond acceptor), only one of the two Pro-containing strands can favorably add to the sheet on each side. Surprisingly, when the sheet is capped with AAPAAA-NH2 at one edge, the interaction between the cap and sheet is slightly more stabilizing than that of another all Ala strand. Breaking down the interaction enthalpies into H-bonding and distortion energies shows the favorable interaction to be due to lower distortion energies in both the strand and the four-stranded sheet. Because another strand would be inhibited for attachment to the other side of the capping (Pro-containing) strand, we suggest the possible use of Pro residues in peptides designed to arrest the growth of many amyloids.

DFT STUDY OF THE PROTONATION AND DEPROTONATION ENTHALPIES OF BENZOXAZOLE, 1,2-BENZISOXAZOLE AND 2,1-BENZISOXAZOLE AND IMPLICATIONS FOR THE STRUCTURES AND ENERGIES OF THEIR ADDUCTS WITH EXPLICIT WATER MOLECULES

Benzoxazole, 1,2-benzisoxazole and 2,1-benzisoxazole are biologically active molecules with potential applications in drug design. Their interaction with aqueous medium in biological systems may be simulated by considering their interaction with explicit water molecules. Such studies provide information on the structures, energies and type of interactions stabilizing the resulting geometric systems. The objective of the current study was to utilize theoretical approaches to investigate the structures, stabilization energy and binding energy of benzoxazole–water, 1,2-benzisoxazole–water and 2,1-benzisoxazole–water complexes. The calculations were performed utilizing the density functional theory (DFT)/M06-2X/6-311 ++ G(d,p) method and the DFT/ωB97XD method with both the 6-311 ++ G(d,p) and the aug-cc-pVDZ basis sets. The results suggest that the stability of the different clusters depends on interrelated factors including the rings formed by intermolecular hydrogen bonds and the proton affinity (PA) or acidity of the atoms forming the intermolecular hydrogen bonds with the water molecules. A comparison across methods indicates that the results follow similar trends with different methods.

Contribution of the empirical dispersion correction on the conformation of short alanine peptides obtained by gas-phase QM calculations

In this work we analyze the effect of the inclusion of an empirical dispersion term to standard DFT (DFT-D) in the prediction of the conformational energy of the alanine dipeptide (Ala2) and in assessing the relative stabilities of short polyala-nine peptides in helical conformations, i.e., α and 310 helices, from Ala4 to Ala16. The Ala2 conformational energies obtained with the dispersion-corrected GGA functional B97-D are compared to previously published high level MP2 data. Meanwhile, the B97-D performance on larger polyalanine peptides is compared to MP2, B3LYP and RHF calculations obtained at a lower level of theory. Our results show that electron correlation affects the conformational energies of short peptides with a weight that increases with the peptide length. Indeed, while the contribution of vdW forces is significant for larger peptides, in the case of Ala2 it is negligible when compared to solvent effects. Even for short peptides, the inclusion of an empirical dispersion term greatly improves accuracy of DFT methods, providing results that correlate very well with the MP2 reference at no additional computational cost.

Comparison of some dispersion-corrected and traditional functionals as applied to peptides and conformations of cyclohexane derivatives

We compare the energetic and structural properties of fully optimized α-helical and antiparallel β-sheet polyalanines and the energetic differences between axial and equatorial conformations of three cyclohexane derivatives (methyl, fluoro, and chloro) as calculated using several functionals designed to treat dispersion (B97-D, ωB97x-D, M06, M06L, and M06-2X) with other traditional functionals not specifically parametrized to treat dispersion (B3LYP, X3LYP, and PBE1PBE) and with experimental results. Those functionals developed to treat dispersion significantly overestimate interaction enthalpies of folding for the α-helix and predict unreasonable structures that contain Ramachandran φ and ψ and C = O...N H-bonding angles that are out of the bounds of databases compiled the β-sheets. These structures are consistent with overestimation of the interaction energies. For the cyclohexanes, these functionals overestimate the stabilities of the axial conformation, especially when used with smaller basis sets. Their performance improves when the basis set is improved from D95∗∗ to aug-cc-pVTZ (which would not be possible with systems as large as the peptides).

The folding of acetyl(Ala)28NH2 and acetyl(Ala)40NH2 extended strand peptides into antiparallel β-sheets. A density functional theory study of β-sheets with β-turns

We report ONIOM calculations using B3LYP/D95** and AM1 on β-sheet formation from acetyl(Ala)(N)NH(2) (N = 28 or 40). The sheets contain from one to four β-turns for N = 28 and up to six for N = 40. We have obtained four types of geometrically optimized structures. All contain only β-turns. They differ from each other in the types of β-turns formed. The unsolvated sheets containing two turns are most stable. Aqueous solvation (using the SM5.2 and CPCM methods) reduces the stabilities of the folded structures compared to the extended strands.