Relationship between structure and entropy contributions in an anthraquinone mercapto derivative

The structural and thermodynamic properties of an anthraquinone derivative were studied by means of quantum-chemical calculations. Conformational analysis using ab initio and density functional theory methods revealed 14 low-energy conformers. In order to discuss similarities and differences in entropy of the conformers, the rotational and vibrational contributions to entropy were correlated with changes in conformer structure. The component of the moment of inertia perpendicular to the molecular plane gives significant input to ΔS _rot , whereas the largest contributions to the ΔS _vib have vibrations associated with the τ _S1C20 coordinate.

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 study on the binding mechanism between N6-methyladenine and natural DNA bases

N6-methyladenine (m⁶A) is a rare base naturally occurring in DNA. It is different from the base adenine due to its N-CH₃. Therefore, the base not only pairs with thymine, but also with other DNA bases (cytosine, adenine and guanine). In this work, Møller-Plesset second-order (MP2) method has been used to investigate the binding mechanism between m⁶A and natural DNA bases in gas phase and in aqueous solution. The results show that N-CH₃ changed the way of N6-methyladenine binding to natural DNA bases. The binding style significantly influences the stability of base pairs. The trans-m⁶A:G and trans-m⁶A:C conformers are the most stable among all the base pairs. The existence of solvent can remarkably reduce the stability of the base pairs, and the DNA bases prefer pairing with trans-m⁶A to cis-m⁶A. Besides, the properties of these hydrogen bonds have been analyzed by atom in molecules (AIM) theory, natural bond orbital (NBO) analysis and Wiberg bond indexes (WBI). In addition, pairing with m⁶A decreases the binding energies compared to the normal Watson-Crick base pairs, it may explain the instability of the N6 site methylated DNA in theory.

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

Simulations on the possibility of formation of complexes between fluorouracil drug and cucurbit[n]urils: ab initio van der Waals DFT study

The binding geometry of fluorouracil/cucurbit[n]urils (CB[n]s) complexes with n?=?5–8 is investigated using the first-principles van der Waals density functional (vdW-DF) method, involving full geometry optimization. Such host-guest complexes are typically calculated using conventional DFT method, which significantly underestimates non-local dispersion forces (or vdW contributions) and therefore affects interactions between respected entities. We address here the role of vdW forces for the fluorouracil and CB[n]s molecules which can form directional hydrogen bonds with each other. It was found that the inclusion of dispersion interactions significantly affects the host-guest binding properties and the hydrogen bonding between the molecules provides the main binding mechanism, while results in the same geometries for the considered complexes. The 0.84 eV binding energy, for the thermodynamically favorable state, reveals that the interaction of fluorouracil with CB[n]s is an exothermic interaction and typical for strong hydrogen bonding. Furthermore, we have investigated the binding nature of these host-guest systems in aqueous solution with ab initio MD simulations adopting vdW-DF method. These findings afford evidence for the applicability of the vdW-DF approach and provide a realistic benchmark for the investigation of the host-guest complexes.

Interplay between halogen bonds and hydrogen bonds in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes

The character of the cooperativity between the HOX···OH/SH halogen bond (XB) and the Y―H···(H)OX hydrogen bond (HB) in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes has been investigated by means of second-order Møller?Plesset perturbation theory (MP2) calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The geometries of the complexes have been determined from the most negative electrostatic potentials (V _S,min) and the most positive electrostatic potentials (V _S,max) on the electron density contours of the individual species. The greater the V _S,max values of HY, the larger the interaction energies of halogen-bonded HOX···OH/SH in the termolecular complexes, indicating that the ability of cooperative effect of hydrogen bond on halogen bond are determined by V _S,max of HY. The interaction energies, binding distances, infrared vibrational frequencies, and electron densities ρ at the BCPs of the hydrogen bonds and halogen bonds prove that there is positive cooperativity between these bonds. The potentiation of hydrogen bonds on halogen bonds is greater than that of halogen bonds on hydrogen bonds. QTAIM studies have shown that the halogen bonds and hydrogen bonds are closed-shell noncovalent interactions, and both have greater electrostatic character in the termolecular species compared with the bimolecular species.

Computational insight into novel molecular recognition mechanism of different bioactive GAs and the Arabidopsis receptor GID1A

Gibberellin (GA) is an essential plant hormone and plays a significant role during the growth and development of the higher plants. The molecular recognition mode between GA and receptor Arabidopsis thaliana GIBBERELLIN INSENSITIVE DWARF1 A (AtGID1A) was investigated by molecular docking and dynamics simulations to clarify the selective perceived mechanism of different bioactive GA molecules to AtGID1A. The 6-COOH group of GA, especially its β configuration, was found to be an indispensable pharmacophore group for GA recognition and binding to AtGID1A. Not only does a strong salt bridge interaction between the 6β-COOH group of GA and Arg244 of AtGID1A play a very important role in the GA recognition of the receptor, but also an indirect water bridge interaction between the pharmacophore group 6β-COOH of GA and the residue Tyr322 of AtGID1A is essential for the GA binding to the receptor. The site-directed residues mutant modeling study on the receptor-binding pocket confirmed that the mutations of Arg244 and Tyr322 decreased the GA binding activity due to the disappearances of the salt bridge and the hydrogen bond interaction. The 3β-OH group of GA was well known to be necessary for the GA bioactivity due to its forming a unique hydrogen bond with Tyr127 of AtGID1A. In addition, the hydrophobic interaction between GA and AtGID1A was considered a necessary factor to lock the GA active conformation and stabilize the GA-GID1A complex structure. The novel molecular recognition mode will be beneficial in elucidating the GA regulation function on the growth and development of the higher plants.

Non-covalent interactions – QTAIM and NBO analysis

MP2(full)/6-311++G(3df,3pd) calculations were carried out on complexes linked through various non-covalent Lewis acid – Lewis base interactions. These are: hydrogen bond, dihydrogen bond, hydride bond and halogen bond. The quantum theory of ´atoms in molecules´ (QTAIM) as well as the natural bond orbitals (NBO) method were applied to analyze properties of these interactions. It was found that for the A-H…B hydrogen bond as well as for the A-X…B halogen bond (X designates halogen) the complex formation leads to the increase of s-character in the A-atom hybrid orbital aimed toward the H or X atom. In opposite, for the A…H-B hydride bond, where the H-atom possesses negative charge, the decrease of s-character in the B-atom orbital is observed. All these changes connected with the redistribution of the electron charge being the effect of the complex formation are in line with Bent´s rule. The numerous correlations between energetic, geometrical, NBO and QTAIM parameters were also found.

Theoretical study of the triangular bonding complex formed by carbon tetrabromide, a halide, and a solvent molecule in the gas phase

MP2(full)/aug-cc-pVDZ(-PP) computations predict that new triangular bonding complexes (where X^? is a halide and H–C refers to a protic solvent molecule) consist of one halogen bond and two hydrogen bonds in the gas phase. Carbon tetrabromide acts as the donor in the halogen bond, while it acts as an acceptor in the hydrogen bond. The halide (which commonly acts as an acceptor) can interact with both carbon tetrabromide and solvent molecule (CH₃CN, CH₂Cl₂, CHCl₃) to form a halogen bond and a hydrogen bond, respectively. The strength of the halogen bond obeys the order CBr₄???Cl^? > CBr₄???Br^? > CBr₄???I^?. For the hydrogen bonds formed between various halides and the same solvent molecule, the strength of the hydrogen bond obeys the order C-H???Cl^? > C-H???Br^? > C-H???I^?. For the hydrogen bonds formed between the same halide and various solvent molecules, the interaction strength is proportional to the acidity of the hydrogen in the solvent molecule. The diminutive effect is present between the hydrogen bonds and the halogen bond in chlorine and bromine triangular bonding complexes. Complexes containing iodide ion show weak cooperative effects.

Synthetic and quantum chemical study on the regioselective addition of amines to methyl maleamate

Synthetic and theoretical studies were performed to gain insight into the regioselectivity in the mechanism of aspartyl-isoaspartyl formation, modeled by additions of ammonia and primary amines to methyl maleamate. Reactions between maleamate and aliphatic, araliphatic amines or O-methyl acetimidate lead to the formation of N-substituted isoasparaginates. The size of the amine and the activating effect of the amide and ester group on the double bond are the determining factors of the site of addition. The formation of both isomers was observed only in the case of ammonia addition. The regioselectivity was predicted on the basis of the charge distribution for low-energy methyl maleamate conformers, calculated at the B3LYP/6-311++G(2df,2pd)//B3LYP/6-31+G(d) level, both in gas phase and in methanol. The methyl isoasparaginate over methyl asparaginate product ratio was computed based on the free energy Boltzmann distribution of their conformers. The calculated 2 : 1 ratio is in agreement with the experimental regioselectivity of the addition of nitrogen nucleophiles.

Modeling the activity of glutathione as a hydroxyl radical scavenger considering its neutral non-zwitterionic form

Glutathione is an immensely important antioxidant, particularly in the central nervous system. The scavenging mechanism of glutathione towards the OH radical was studied theoretically, considering its neutral, non-zwitterionic form relevant to acidic media. Gibbs free barrier and released energies involved in hydrogen abstraction from the different sites of glutathione by an OH radical were studied at the B3LYP/6-31G(d,p), B3LYP/AUG-cc-pVDZ, M06/AUG-cc-pVDZ, M06-2X/AUG-cc-pVDZ levels of density functional theory. Solvation in bulk aqueous media was also studied at all these levels of theory employing the polarizable continuum model. Our study shows that a hydroxyl radical can abstract a hydrogen atom easily from glutathione. Thus, glutathione is shown to be an efficient scavenger of OH radicals, which is in agreement with the results of previous studies.

A B3LYP and MP2(full) theoretical investigation into the strength of the C-NO2 bond upon the formation of the molecule-cation interaction between Na+ and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the molecule-cation interaction between Na⁺ and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced in comparison with that in the isolated nitrotriazole molecule upon the formation of molecule-cation interaction. The increment of the C–NO₂ bond dissociation energy (ΔBDE) correlated well with the molecule-cation interaction energy. Electron density shifts analysis showed that the electron density shifted toward the C-NO₂ bond upon complex formation, leading to the strengthened C-NO₂ bond and the possibly reduced explosive sensitivity.

Why tazobactam and sulbactam have different intermediates population with SHV-1 β-lactamase: a molecular dynamics study

The imine intermediates of tazobactam and sulbactam bound to SHV-1 β-lactamase were investigated by molecular dynamics (MD) simulation respectively. Hydrogen bond networks around active site were found different between tazobactam and sulbactam acyl-enzymes. In tazobactam imine intermediate, it was observed that the triazolyl ring formed stable hydrogen bonds with Asn170 and Thr167. The results suggest that conformation of imine determined the population of intermediates. In imine intermediate of tazobactam, the triazolyl ring is trapped in Thr_Asn pocket, and it restricts the rotation of C5-C6 bond so that tazobactam can only generate trans enamine intermediate. Further, conformational cluster analyses are performed to substantiate the results. These findings provide an explanation for the corresponding experimental results, and will be potentially useful in the development of new inhibitors.

A B3LYP and MP2(full) theoretical investigation into the strength of the C–NO2 bond upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced and the charge of nitro group turned more negative in complex in comparison with those in isolated nitrotriazole molecule. The increment of the C–NO₂ bond dissociation energies correlated well with the intermolecular H-bonding interaction energies. Electron density shifts analyses showed that the electron density shifted toward the C–NO₂ bond upon complex formation, leading to the strengthened C–NO₂ bond and the possibly reduced explosive sensitivity.

Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.

Density functional studies on photophysical properties and chemical reactivities of the triarylboranes: effect of the constraint of planarity

The geometric and electronic structures, absorption spectra, transporting properties, chemical reactivity indices and electrostatic potentials of the planar three-coordinate organoboron compounds 1-2 and twisted reference compound Mes ₃ B, have been investigated by employing density functional theory (DFT) and conceptual DFT methods to shed light on the planarity effects on the photophysical properties and the chemical reactivity. The results show that the planar compounds 1-2 exhibit significantly lower HOMO level than Mes ₃ B, owing to the stronger electronic induction effect of boron centers. This feature conspicuously induces a blue shifted absorption for 1, although 1 seemingly possesses more extended conjugation framework than Mes ₃ B. Importantly, the reactivity strength of the boron atoms in 1-2 is much lower than that in Mes ₃ B, despite the fact that the tri-coordinate boron centers of 1-2 are completely naked. The interesting and abnormal phenomenon is caused by the strong p-π electronic interactions, that is, the empty p-orbital of boron center is partly filled by π-electron of the neighbor carbon atoms in 1-2, which are confirmed by the analysis of Laplacian of the electron density and natural bond orbitals. Furthermore, the negative electrostatic potentials of the boron centers in 1-2 also interpret that they are not the most preferred sites for incoming nucleophiles. Moreover, it is also found that the planar compounds 1-2 can act as promising electron transporting materials since the internal reorganization energies for electron are really small.

Investigation of the antioxidant properties of hyperjovinol A through its Cu(II) coordination ability

Hyperjovinol A (2-methyl-1-(2,4,6-trihydroxy-3-(3-hydroxy-3,7-dimethyloct-6-enyl)phen yl)propan-1-one) is an acylated phloroglucinol isolated from Hypericum Jovis and exhibiting antioxidant properties comparable with those of the most common antioxidant drugs. The study models the compound’s antioxidant ability through its ability to coordinate a Cu²⁺ ion and reduce it to Cu⁺. Complexes with a Cu²⁺ ion were calculated for all the low energy and for representative high energy conformers of hyperjovinol A, placing the ion in turn near each of the electron-rich binding sites. The most stable complexes are those in which Cu²⁺ binds simultaneously to the O of the OH in the geranyl-type chain (R′) and the C═C double bond at the end of R′, or to the O of a phenol OH and the O of the OH in R′. The most stable complexes in which Cu²⁺ binds only to one site are those in which it binds to the C═C double bond at the end of R′ or to the sp² O of the COCH(CH₃)₂ acyl group. Cu²⁺ is reduced to Cu⁺ in all complexes. Comparisons with corresponding complexes of other molecular structures in which one or more of the structural features of hyperjovinol A are modified attempt to elucidate the role, for the antioxidant ability, of relevant features of hyperjovinol A, like the presence and position of the OH or the C═C double bond in R′. Calculations at the DFT/B3LYP/6–31+G(d,p) level were performed for all the structures considered. Calculations utilizing the LANL2DZ pseudopotential for the Cu²⁺ ion were also performed for hyperjovinol A.

Microsolvation of Mg2+, Ca2+: strong influence of formal charges in hydrogen bond networks

A stochastic exploration of the quantum conformational spaces in the microsolvation of divalent cations with explicit consideration of up to six solvent molecules [Mg (H ₂ O) _n )]²⁺, (n?=?3, 4, 5, 6) at the B3LYP, MP2, CCSD(T) levels is presented. We find several cases in which the formal charge in Mg²⁺ causes dissociation of water molecules in the first solvation shell, leaving a hydroxide ion available to interact with the central cation, the released proton being transferred to outer solvation shells in a Grotthus type mechanism; this particular finding sheds light on the capacity of Mg²⁺ to promote formation of hydroxide anions, a process necessary to regulate proton transfer in enzymes with exonuclease activity. Two distinct types of hydrogen bonds, scattered over a wide range of distances (1.35–2.15 Å) were identified. We find that in inner solvation shells, where hydrogen bond networks are severely disturbed, most of the interaction energies come from electrostatic and polarization+charge transfer, while in outer solvation shells the situation approximates that of pure water clusters.