Theoretical study of the symmetry of the (OH...O)- hydrogen bonds in vinyl alcohol-vinyl alcoholate systems

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Theoretical study on L-H+-L with identical donors: short strong hydrogen bond or not?

Short strong hydrogen bonds (SSHBs) play crucial role in many chemical processes. Recently, as the representative of SSHBs, [F-H-F]- was experimentally observed. [F-H-F]- has a symmetric structure, which can be described as a H+ acid shared by two terminal F- donors (F--H+-F-). To explore whether two identical donors are bound to result in SSHBs, we performed theoretical studies on a series of compounds (L-H+-L) with two identical electron donors (L corresponds to donors containing group 14, 15, 16 and 17 elements). The results show that identical donors do not definitely lead to SSHBs. Instead, typical hydrogen bonds also exist. We found that both electronegativity and basicity contribute to the patterns of hydrogen bonds, where more electronegative and weaker donors benefit to SSHBs. Besides, it was found that zero-point energies also respond to the hydrogen bonding systems. This systemic work is expected to provide more insights into SSHBs.

Synthesis, crystal structures, cytotoxicity and qualitative structure-activity relationship (QSAR) of cis-bis{5-[(E)-2-(aryl)-1-diazenyl]quinolinolato}di-n-butyltin(IV) complexes, (n)Bu2Sn(L)2

A series of cis-bis{5-[(E)-2-(aryl)-1-diazenyl]quinolinolato}di-n-butyltin(IV) complexes has been synthesized and characterized by (1)H-, (13)C-, (119)Sn NMR, ESI-MS (electrospray ionization mass spectrometry), IR and (119m)Sn Mössbauer spectroscopic techniques in combination with elemental analyses. The structures of four di-n-butyltin(IV) complexes, viz., (n)Bu(2)Sn(L(3))(2) (3), (n)Bu(2)Sn(L(4))(2) (4), (n)Bu(2)Sn(L(5))(2) (5) and (n)Bu(2)Sn(L(7))(2).0.5C(6)H(6) (7) (LH=5-[(E)-2-(aryl)-1-diazenyl)quinolin-8-ol) were determined by single crystal X-ray diffraction. In general, the complexes were found to adopt a distorted cis-octahedral arrangement around the tin atom. These complexes retain their solid-state structure in non-coordinating solvent as evidenced by (119)Sn and (13)C NMR spectroscopic results. The in vitro cytotoxicity of di-n-butyltin(IV) complexes (3-8) is reported against seven well characterized human tumour cell lines. The basicity of the two quinolinolato donor N and O atoms of the ligands are discussed in relation to the cytotoxicity data.

Effects of OH Radical Addition on Proton Transfer in the Guanine−Cytosine Base Pair

Double proton transfer (PT) reactions in guanine-cytosine OH radical adducts are studied by the hybrid density functional B3LYP approach. Concerted and stepwise proton-transfer processes are explored between N1(H) on guanine (G) and N3 on cytosine (C), and between N4(H) on C and O6 on G. All systems except GC6OH display a concerted mechanism. 8OHGC has the highest dissociation energy and is 1.2 kcal/mol more stable than the nonradical GC base pair. The origin of the interactions are investigated through the estimation of intrinsic acid-basic properties of the *OH-X monomer (X = G or C). Solvent effects play a significant role in reducing the dissociation energy. The reactions including *OH-C adducts have significantly lower PT barriers than both the nonradical GC pair and the *OH-G adducts. All reactions are endothermic, with the GC6OH --> GC6OHPT reaction has the lowest reaction energy (4.6 kcal/mol). In accordance with earlier results, the estimated NBO charges show that the G moiety carries a slight negative charge (and C a corresponding positive one) in each adduct. The formation of a partial ion pair may be a potential factor leading to the PT reactions being thermodynamically unfavored.

Correlated proton motion in hydrogen bonded systems: tuning proton affinities

The theorem of matching proton affinities (PA) has been widely used in the analysis of hydrogen bonds. However, most experimental and theoretical investigations have to cope with the problem that the variation of the PA of one partner in the hydrogen bond severely affects the properties of the interface between both molecules. The B3LYP/d95+(d,p) analysis of two hydrogen bonds coupled by a 5-methyl-1H-imidazole molecule showed that it is possible to change the PA of one partner of the hydrogen bond while maintaining the properties of the interface. This technique allowed us to correlate various properties of the hydrogen bond directly with the difference in the PAs between both partners: it is possible to tune the potential energy surface of the bonding hydrogen atom from that of an ordinary hydrogen bond (localized hydrogen atom) to that of a low barrier hydrogen bond (LBHB, delocalized hydrogen atom) just by varying the proton affinity of one partner. This correlation shows clearly that matching PAs are of lesser importance for the formation of a LBHB than the relative energy difference between the two tautomers of the hydrogen bond.

Theoretical Study of (CH···C)- Hydrogen Bonds in CH4-nXn(X = F, Cl; n = 0, 1, 2) Systems Complexed with Their Homoconjugate and Heteroconjugate Carbanions

This work deals with a theoretical study of the (CH...C)- hydrogen bonds in CH4, CH3X, and CH2X2 (X = F, Cl) complexed with their homoconjugate and heteroconjugate carbanions. The properties of the complexes are calculated with the B3LYP method using the 6-311++G(d,p) or 6-311++G(2df,2p) basis sets. The deprotonation enthalpies (DPE) of the CH bond or the proton affinities of the carbanions (PA(C-) are calculated as well. All the systems with the exception of the CH4...CHCl2(-) one are characterized by a double minimum potential. In some of the complexes, the (CH(b)...C)- hydrogen bond is linear. In other systems, such as CH3F...CH2F- and CH3F...CHF2(-), there is a large departure from linearity, the systems being stabilized by electrostatic interactions between the nonbonded H of the neutral molecule and the F atom of the carbanion. In the transition state, the (CH(b)...C)- bond is linear, and there is a large contraction of the intermolecular C...C distance. The binding energies vary within a large range, from -1.4 to -11.1 kcal mol(-1) for the stable complexes and -8.6 to -44.1 kcal mol(-1) for the metastable complexes. The energy barriers to proton transfer are between 5 and 20 kcal mol(-1) for the heteroconjugate systems and between 3.8 and 8.3 kcal mol(-1) for the homoconjugate systems. The binding energies of the linear complexes depend exponentially on 1.5DPE - PA(C-), showing that the proton donor is more important than the proton acceptor in determining hydrogen bond strength. The NBO analysis indicates an important electronic reorganization in the two partners. The elongations of the CH bond resulting from the interaction with the carbanion depend on the occupation of the sigma*(CH(b)) antibonding orbitals and on the hybridization of the C bonded to H(b). The frequency shifts of the nu(CH)(A1) stretching vibration range between 15 and 1150 cm(-1). They are linearly correlated to the elongation of the CH(b) bond.

Energetic Preferences for α,β versus β,γ Unsaturation

Density functional theory has been applied at the B3LYP/6-311+G(d,p)//B3LYP/6-31G(d) level to examine the energetics of alpha,beta- versus beta,gamma-unsaturation for some common organic functional groups. Specifically, the relative stabilities of allyl-X (H2C=CHCH2X) and 1-propenyl-X (H3CCH=CHX) isomers have been computed for X = methyl, vinyl, phenyl, formyl, acetyl, methoxy, methylthio, methylsulfinyl, methylsulfonyl, sulfamoyl, and methoxysulfonyl, and the results are compared to available experimental data. The intrinsic preference of 3 kcal/mol for the 1-propenyl isomer when X = CH3 is exceeded by 2-4 kcal/mol for first-row conjugating groups, but it is not met for the sulfur-containing groups. In particular, alpha,beta-unsaturation is favored by less than 1 kcal/mol for the sulfone and sulfonamide analogues, while it is preferred by 8 kcal/mol for the vinyl-substituted case. Detailed structural results and torsional energy profiles are also reported.