Anticooperativity of FH···Cl − hydrogen bonds in [FH) n Cl] − clusters ( n = 1…6)

Dynamics of non-covalent interactions during the P-O bond cleavage reaction by ribonuclease A

In this work, an atomistic-scale investigation of the phosphodiester P-O bond cleavage reaction by the enzyme ribonuclease A was carried out using computer simulation techniques. It is shown that during the reaction the network of non-covalent interactions in the active center of the ribonuclease changes significantly, while the role of these non-covalent interactions is different: coordination of the corresponding groups, electron density transfer, and ligand holding in the active center. It is shown that the process of proton transfer from Asp121 to His119 is the first stage of this reaction; at the same time, the hydrogen bond between the phosphate ligand and the imino group of Arg39 is broken, which, although keeping the ligand in the active center, does not allow the ligand to orient itself more conveniently for subsequent proton transfers. Furthermore, the key step of this reaction occurs: proton transfer with the participation of imidazole rings His12 and His119, in which the guiding role is played by several hydrogen bonds with the participation of Phe120, and the role of an electron density carrier is played by the pnictogen bond between the oxygen of the phosphate ligand and the pyridine-like nitrogen of the imidazole ring His119, which was detected for the first time.

Mutual influence of non-covalent interactions formed by imidazole: A systematic quantum-chemical study

Imidazole is a five-membered heterocycle that is part of a number of biologically important molecules such as the amino acid histidine and the hormone histamine. Imidazole has a unique ability to participate in a variety of non-covalent interactions involving the NH group, the pyridine-like nitrogen atom or the π-system. For many biologically active compounds containing the imidazole moiety, its participation in formation of hydrogen bond NH⋯O/N and following proton transfer is the key step of mechanism of their action. In this work a systematic study of the mutual influence of various paired combinations of non-covalent interactions (e.g., hydrogen bonds and π-interactions) involving the imidazole moiety was performed by means of quantum chemistry (PW6B95-GD3/def2-QZVPD) for a series of model systems constructed based on analysis of available x-ray data. It is shown that for considered complexes formation of additional non-covalent interactions can only enhance the proton-donating ability of imidazole. At the same time, its proton-accepting ability can be both enhanced and weakened, depending on what additional interactions are added to a given system. The mutual influence of non-covalent interactions involving imidazole can be classified as weak geometric and strong energetic cooperativity-a small change in the length of non-covalent interaction formed by imidazole can strongly influence its strength. The latter can be used to develop methods for controlling the rate and selectivity of chemical reactions involving the imidazole fragment in larger systems. It is shown that the strong mutual influence of non-covalent interactions involving imidazole is due to the unique ability of the imidazole ring to effectively redistribute electron density in non-covalently bound systems with its participation.

Evolution of vibrational bands upon gradual protonation/deprotonation of arsinic acid H2As(O)OH in media of different polarity

This computational work is devoted to the investigation (MP2/def2-TZVP) of the geometry and IR parameters of arsinic acid H₂AsOOH and its hydrogen-bonded complexes under vacuum and in media with different polarity. The medium effects were accounted for in two ways: (1) implicitly, using the IEFPCM model, varying the dielectric permittivity (ε) and (2) explicitly, by considering hydrogen-bonded complexes of H₂As(O)OH with various hydrogen bond donors (41 complexes) or acceptors (38 complexes), imitating a gradual transition to the As(OH)₂⁺ or AsO₂^- moiety, respectively. It was shown that the transition from vacuum to a medium with ε > 1 causes the As(O)OH fragment to lose its flatness. The solvent polar medium introduces significant changes in the geometry and IR spectral parameters of hydrogen-bonded complexes too: as the polarity of a medium increases, weak hydrogen bonds become weaker, and strong and medium hydrogen bonds become stronger; in the case of a complex with two hydrogen bonds cooperativity effects are observed. In almost all cases the driving force of these changes appears to be preferential solvation of charge-separated structures. In the limiting case of complete deprotonation (or conversely complete protonation) the vibrational frequencies of ν_AsO and ν_As-O turn into ν_As-O(asym) and ν_As-O(sym), respectively. In the intermediate cases the distance between ν_AsO and ν_As-O is sensitive to both implicit solvation and explicit solvation and the systematic changes of this distance can be used for estimation of the degree of proton transfer within the hydrogen bond.

Estimations of OH·N hydrogen bond length from positions and intensities of IR bands

In this computational work applicability of IR spectral parameters for evaluations of OH···N hydrogen bond length is discussed. For a set of 124 complexes with OH···N hydrogen bond formed by combinations of methanol/acetic acid and pyridine (and their fluorine substituted versions) geometries, energies and IR parameters were calculated at MP2/def2-TZVP level of theory. For a number of IR parameters (the shift of proton donor group stretching vibration Δν_s, increase of its intensity I, the low-frequency hydrogen bond stretching vibration ν_σ, bending in-plane δ and out-of-plane γ vibrations) equations linking them with interatomic distances are proposed, the robustness and accuracy of such equations are discussed. The enthalpy of OH···N hydrogen bond formation ΔH was also linked with electron density parameters in (3; -1) critical point.

Self-association of diphenylpnictoginic acids in solution and solid state: covalent vs. hydrogen bonding

Triphenylpnictogens were oxidized to access diphenylpnictioginic acids Ph₂XOOH (X = P, As, Sb, Bi). It was shown that oxidation with chloramine-T does not lead to the cleavage of a C-pnictogen bond. The preliminary reductive cleavage with sodium in liquid ammonia followed by the oxidation with hydrogen peroxide was successfully utilised for the synthesis of diphenylphosphinic and diphenylarsinic acids. It was shown that in solid state (by means of XRD), all diphenylpnictoginic acids form polymeric chains. Diphenylbismuthinic and diphenylantimonic acids form polymeric covalent adducts, while diphenylphosphinic and diphenylarsinic chains are associated through hydrogen bonding. Unlike diphenylphosphinic acid, diphenilarsinic acid forms two polymorphs of hydrogen-bonded infinite chains. In solution in a polar aprotic solvent diphenylarsinic acid, similarly to dimethylarsinic, forms hydrogen-bonded cyclic dimers together with a small amount of cyclic trimers.

Hydration of selenolate moiety: Ab initio investigation of properties of O-H⋯Se(-) hydrogen bonds in CH3 Se(-)⋯(H2 O)n clusters

This work is devoted to investigations of the influence of O-H···Se(-) hydrogen bonds on the electronic shells of selenolate R-Se(-) fragment (R═CH₃ ). The geometric, energetic and nuclear magnetic resonance (NMR) spectral parameters for various conformers of CH₃ Se(-)⋯(H₂ O)_n clusters with n = 0-6 were calculated at CCSD/aug-cc-pVDZ level of theory. For selenolate anion CH₃ Se(-) solvation free energy was calculated, and for water media it is equal to -71.41 kcal/mol. For O-H···Se(-) hydrogen bonds the proportionality coefficients between QTAIM parameters at (3; -1) bond critical point and the strength of an individual hydrogen bond ∆E were proposed. It was shown, that despite a relative weakness of O-H···Se(-) hydrogen bonds, the outer electronic shell of the selenium atom changes significantly upon formation of each hydrogen bond. This, in turn, cause the dramatic change of NMR parameters of selenium nuclei.

On the influence of water molecules on the outer electronic shells of R-SeH, R-Se(-) and R-SeOH fragments in the selenocysteine amino acid residue

In this computational work (MP2/aug-cc-pVTZ) we investigated the features of the outer electronic shells of R-SeH, R-Se(-) and R-SeOH fragments (R = CH3), which can be considered as simplified models for the forms of the active centres of glutathione peroxidases GPx along their catalytic pathway (reduction of peroxides). It is shown that the preferential direction of a nucleophilic attack on the R-Se(-) fragment by a peroxide molecule is determined by the presence of the electron-depleted region of the selenium atom in front of the C-Se bond and nucleophilic attack can be facilitated by the solvation of R-Se(-) by water molecules. Such solvation does not block the direction of potential nucleophilic attack and also leads to the increase of the maximal value of the molecular electrostatic potential on the selenium atom. It was shown that the 77Se NMR chemical shift is sensitive both to the oxidation state and the hydration state of the selenium-containing fragment.

The explicit role of electron exchange in the hydrogen bonded molecular complexes

We applied a set of advanced bonding descriptors to establish the hidden electron density features and binding energy characteristics of intermolecular DH∙∙∙A hydrogen bonds (OH∙∙∙O, NH∙∙∙O and SH∙∙∙O) in 150 isolated and solvated molecular complexes. The exchange-correlation and Pauli potentials as well as corresponding local one-electron forces allowed us to explicitly ascertain how electron exchange defines the bonding picture in the proximity of the H-bond critical point. The electron density features of DH∙∙∙A interaction are governed by alterations in the electron localization in the H-bond region displaying itself in the exchange hole. At that, they do not depend on the variations in the exchange hole mobility. The electrostatic interaction mainly defines the energy of H-bonds of different types, whereas the strengthening/weakening of H-bonds in complexes with varying substituents depends on the barrier height of the exchange potential near the bond critical point. Energy variations between H-bonds in isolated and solvated systems are also caused the electron exchange peculiarities as follows from the corresponding potential and the interacting quantum atom analyses complemented by electron delocalization index calculations. Our approach is based on the bonding descriptors associated with the characteristics of the observable electron density and can be recommended for in-depth studies of non-covalent bonding.

Actual Symmetry of Symmetric Molecular Adducts in the Gas Phase, Solution and in the Solid State

This review discusses molecular adducts, whose composition allows a symmetric structure. Such adducts are popular model systems, as they are useful for analyzing the effect of structure on the property selected for study since they allow one to reduce the number of parameters. The main objectives of this discussion are to evaluate the influence of the surroundings on the symmetry of these adducts, steric hindrances within the adducts, competition between different noncovalent interactions responsible for stabilizing the adducts, and experimental methods that can be used to study the symmetry at different time scales. This review considers the following central binding units: hydrogen (proton), halogen (anion), metal (cation), water (hydrogen peroxide).

Stretching force constants as descriptors of energy and geometry of F···HF hydrogen bonds

In this work applicability of proton donor group stretching vibration force constants k_s and intermolecular stretching force constants k_σ for evaluations of hydrogen bond strength and geometry are discussed. For a set of 30 complexes with F···HF hydrogen bonds in a wide range 0.5-48 kcal/mol by means of quantum chemical calculations equilibrium geometries, complexation energies, vibrational frequencies and corresponding force constants were calculated (MP2/aug-cc-pVTZ). It is shown, that properties of a hydrogen bond are more strictly correlated with the values of force constants than with vibrational frequencies. Easy-to-use equations for estimations of hydrogen bond energy ∆E and geometry (r_FH, r_FF) based on k_s and k_σ values are proposed.