Hydrogen-bonding studies of amino acid side-chains with DNA base pairs

Unraveling the Nature of Hydrogen Bonds of "Proton Sponges" Based on Car-Parrinello and Metadynamics Approaches

The nature of intra- and intermolecular non-covalent interactions was studied in four naphthalene derivatives commonly referred to as "proton sponges". Special attention was paid to an intramolecular hydrogen bond present in the protonated form of the compounds. The unsubstituted "proton sponge" served as a reference structure to study the substituent influence on the hydrogen bond (HB) properties. We selected three compounds substituted by methoxy, amino, and nitro groups. The presence of the substituents either retained the parent symmetry or rendered the compounds asymmetric. In order to reveal the non-covalent interaction properties, the Hirshfeld surface (HS) was computed for the crystal structures of the studied compounds. Next, quantum-chemical simulations were performed in vacuo and in the crystalline phase. Car-Parrinello molecular dynamics (CPMD), Path Integral Molecular Dynamics (PIMD), and metadynamics were employed to investigate the time-evolution changes of metric parameters and free energy profile in both phases. Additionally, for selected snapshots obtained from the CPMD trajectories, non-covalent interactions and electronic structure were studied. Quantum theory of atoms in molecules (QTAIM) and the Density Overlap Regions Indicator (DORI) were applied for this purpose. It was found based on Hirshfeld surfaces that, besides intramolecular hydrogen bonds, other non-covalent interactions are present and have a strong impact on the crystal structure organization. The CPMD results obtained in both phases showed frequent proton transfer phenomena. The proton was strongly delocalized in the applied time-scale and temperature, especially in the PIMD framework. The use of metadynamics allowed for tracing the free energy profiles and confirming that the hydrogen bonds present in "proton sponges" are Low-Barrier Hydrogen Bonds (LBHBs). The electronic and topological analysis quantitatively described the temperature dependence and time-evolution changes of the electronic structure. The covalency of the hydrogen bonds was estimated based on QTAIM analysis. It was found that strong hydrogen bonds show greater covalency, which is additionally determined by the proton position in the hydrogen bridge.

Understanding the potency of malarial ligand (D44) in plasmodium FKBP35 and modelled halogen atom (Br, Cl, F) functional groups

The present study clearly depicts the understanding of the D44 in Plasmodium FKBP35 around the hinge region. To analyse the binding stability of D44 ligand and to understand the role of halogen bond, hydrogen bond interaction formed between the hinge region amino acids: Isoleucine (Ile74), Phenylalanine (Phe54), Aspartic acid (Asp55) Phenylalanine (Phe64),Tyrosine (Tyr100), Tryptophan (TRP 77) and ligand D44 was portrayed specifically through interaction energy calculations at HF, M062X, MP2 level of theories for different basis set (6-311G**, 6-31+G*, LANL2DZ). The investigation will provide an apparent picture regarding the non-covalent interaction that hold the contact of ligand and amino acids in the hinge region and the implication of modelled functional groups (Br, Cl, F, OSO and NH₂) on ligand, which will help chemist in synthesizing new novel ligands. HOMO, LUMO chart calculated for D44 ligands reveals graphic illustration of orbital's that stimulate for contact. The aim and natural bond orbital analysis identified key contribution of individual hydrogen/halogen bonds that contribute for the binding strength through stabilization energy, ρ and ∇²ρ values. Overall this study finds out that the Stability of D44 in Plasmodium FKBP35 was enhanced by the Halogen atom (Br, Cl, F) functional groups; which provide an innovative pathway for the selection of functional groups that opt for the hinge region side chains on the ligand.

A quantum chemical perspective on the potency of electron donors and acceptors in pnicogen bonds (AS...N, P...N, N...N)

A quantum chemical perspective of 31 structures contains electron acceptors: ASCl₃ (arsenic trichloride), PCl₃ (phosphorous trichloride) and NCl₃ (nitrogen trichloride); forming non-covalent bond with various nitrogen-based electron donors that resulted in pnicogen bonds, AS...N, P...N and N...N were calculated at M062X/def2-QZVP level of theory. Besides the above method, MP2/def2-QZVP and CCSD(T)/def2-QZVP level of theories have also been analysed to have in depth knowledge about the bonds formed. The nature of the bonds was assumed from the electrostatic potential evaluated for all the monomers, where σ hole is positive for all the monomers. The strongest pnicogen bonds are ASCl₃-NF₂H, PCl₃-NCH₃CH₃CH₃ and NCl₃-NCH₃CH₃CH₃ having interaction energies as -4.15, -11.58 and -3.25 kcal/mol, respectively, at MP2/def2-QZVP level of theory. Further at CCSD(T)/def2-QZVP level of theory, ASCl₃-NF₂H and NCl₃-NCH₃CH₃CH₃ are found to be the most stable with interaction energies as -3.53 and -2.45 kcal/mol, respectively. The potential energy surface scan was performed for all the stable complexes in order to confirm the existences of energies are true minima. Moreover to confirm the halogen and pnicogen bonds, AIM analysis was carried out. The results from the above factors of pnicogen bond will help crystal growth, material science and engineering community to explore novel materials, which abide for modernization. Graphical abstract PCl₃-NCH₃CH₃CH₃ complex with 2.61 Å and pnicogen angle of 178.54° is strong, and interaction energy is -11.58 kcal/mol. Electron donors - ASCl₃, PCl₃ and NCl₃ and electron acceptors -NCH₃CH₃CH₃, NH₃C₂ and NHCO have strong electrostatic contribution. High and low values of (ρ) ∇²(ρ) reveal the strong and weak pnicogen bond. Schematic representation of acceptors surrounded by its donors and Electrostatic Potential map.

Hydrogen bond docking preference in furans: OH⋯π vs. OH⋯O

The docking sites of hydrogen bonds in complexes formed between 2,2,2-trifluoroethanol (TFE), furan (Fu), and 2-methyl furan (MF) have been investigated. Using density functional theory (DFT) calculations, gas phase and matrix isolation FTIR spectroscopies, the strengths of OH⋯O and OH⋯π hydrogen bonds in the complexes were compared to find the docking preference. Calculations suggest that the hydrogen bond donor, TFE, is more likely to dock onto the oxygen atom of the aromatic furans ring, and consequently, the OH⋯O type hydrogen bond is relatively stronger than the OH⋯π type. The FTIR spectrum in the OH-stretching fundamental range obtained at room temperatures has been compared with that obtained at extremely low temperatures in the matrix. The fundamental and the red shifts of OH-stretching vibrations were observed in both FTIR spectra, confirming the formation of hydrogen bonded complexes. By assessing the ability of furan and MF to participate in the formation of OH⋯O hydrogen bond, the effect of ring methylation has been highlighted. From the calculated geometric and thermodynamic parameters as well as the frequency shift of the OH-stretching vibrations in complexes, TFE-MF is found to be more stable than TFE-Fu, which suggests that the strength of the OH⋯O hydrogen bond in TFE-MF originates from the high activity of the furan molecule caused by the methylation of the aromatic ring. The present study furthers the knowledge of docking preference in heteroaromatic molecules and is helpful to understand the nature of intermolecular interactions between hydrogen bond donors and acceptors, including both electron-deficient atoms and π cloud.

Energy decomposition scheme based on the generalized Kohn-Sham scheme

In this paper, a new energy decomposition analysis scheme based on the generalized Kohn-Sham (GKS) and the localized molecular orbital energy decomposition analysis (LMO-EDA) scheme, named GKS-EDA, is proposed. The GKS-EDA scheme has a wide range of DFT functional adaptability compared to LMO-EDA. In the GKS-EDA scheme, the exchange, repulsion, and polarization terms are determined by DFT orbitals; the correlation term is defined as the difference of the GKS correlation energy from monomers to supermolecule. Using the new definition, the GKS-EDA scheme avoids the error of LMO-EDA which comes from the separated treatment of EX and EC functionals. The scheme can perform analysis both in the gas and in the condensed phases with most of the popular DFT functionals, including LDA, GGA, meta-GGA, hybrid GGA/meta-GGA, double hybrid, range-separated (long-range correction), and dispersion correction. By the GKS-EDA scheme, the DFT functionals assessment for hydrogen bonding, vdW interaction, symmetric radical cation, charge-transfer, and metal-ligand interaction is performed.