It is postulated that the overexpression of Carbonic Anhydrase isozyme IX in some cancers contributes to the acidification of the extracellular matrix. It was proved that this promotes the growth and ...metastasis of the tumor. These observations have made Carbonic Anhydrase IX an attractive drug target. In the light of the findings and importance of the glycoprotein in the cancer treatment, we have employed quantum-chemical approaches to study non-covalent interactions in the binding pocket. As a ligand, the acetazolamide (AZM) molecule was chosen, being known as a potential inhibitor exhibiting anticancer properties. First-Principles Molecular Dynamics was performed to study the chalcogen and other non-covalent interactions in the AZM ligand and its complexes with amino acids forming the binding site. Based on Density Functional Theory (DFT) and post-Hartree-Fock methods, the metric and electronic structure parameters were described. The Non-Covalent Interaction (NCI) index and Atoms in Molecules (AIM) methods were applied for qualitative/quantitative analyses of the non-covalent interactions. Finally, the AZM-binding pocket interaction energy decomposition was carried out. Chalcogen bonding in the AZM molecule is an important factor stabilizing the preferred conformation. Free energy mapping via metadynamics and Path Integral molecular dynamics confirmed the significance of the chalcogen bond in structuring the conformational flexibility of the systems. The developed models are useful in the design of new inhibitors with desired pharmacological properties.
The intra- and intermolecular interactions of selected quinolone carboxylic acid derivatives were studied in monomers, dimers and crystals. The investigated compounds are well-recognized as medicines ...or as bases for further studies in drug design. We employed density functional theory (DFT) in its classical formulation to develop gas-phase and solvent reaction field (PCM) models describing geometric, energetic and electronic structure parameters for monomers and dimers. The electronic structure was investigated based on the atoms in molecules (AIM) and natural bond orbital (NBO) theories. Special attention was devoted to the intramolecular hydrogen bonds (HB) present in the investigated compounds. The characterization of energy components was performed using symmetry-adapted perturbation theory (SAPT). Finally, the time-evolution methods of Car-Parrinello molecular dynamics (CPMD) and path integral molecular dynamics (PIMD) were employed to describe the hydrogen bond dynamics as well as the spectroscopic signatures. The vibrational features of the O-H stretching were studied using Fourier transformation of the autocorrelation function of atomic velocity. The inclusion of quantum nuclear effects provided an accurate depiction of the bridged proton delocalization. The CPMD and PIMD simulations were carried out in the gas and crystalline phases. It was found that the polar environment enhances the strength of the intramolecular hydrogen bonds. The SAPT analysis revealed that the dispersive forces are decisive factors in the intermolecular interactions. In the electronic ground state, the proton-transfer phenomena are not favourable. The CPMD results showed generally that the bridged proton is localized at the donor side, with possible proton-sharing events in the solid-phase simulation of stronger hydrogen bridges. However, the PIMD enabled the quantitative estimation of the quantum effects inclusion-the proton position was moved towards the bridge midpoint, but no qualitative changes were detected. It was found that the interatomic distance between the donor and acceptor atoms was shortened and that the bridged proton was strongly delocalized.
The set of TX3-TrX2 (T = C, Si, Ge; Tr = B, Al, Ga; X = F, Cl, Br) molecules offers a rather unique opportunity to study both σ-hole and π-hole dimerization on the tetrel and triel ends, ...respectively. According to the molecular electrostatic potential (MEP) distribution, the π-hole extrema (acidic sites) were more intense than their σ-hole counterparts. The molecules owning the most (CX3-AlX2) and least (SiX3-BX2) intense π-holes were chosen to evaluate their capacities to attract one and two HCN molecules (Lewis bases). We discovered that the energetic characteristics of π-hole dimers severely conflict with the monomers MEP pattern since the weakest π-hole monomer forms a dimer characterized by interaction energy compared to those created by the monomers with noticeably greater power in the π-hole region. This outcome is due to the deformation of the weakest π-hole donor. Furthermore, the MEP analysis for monomers in the geometry of respective dimers revealed a “residual π-hole” site that was able to drive second ligand attachment, giving rise to the two “unusual trimers” examined further by the NCI and QTAIM analyses. Apart from them, the π-hole/π-hole and σ-hole/π-hole trimers have also been obtained throughout this study and described using energetic and geometric parameters. The SAPT approach revealed details of the bonding in one of the “unusual trimers”. Finally, Born-Oppenheimer Molecular Dynamics (BOMD) simulations were carried out to investigate the time evolution of the interatomic distances of the studied complexes as well as their stability.
Two core-modified porphyrins containing metal atoms, namely platinum(
ii
) or platinum(
iv
) and rhodium(
iii
), in place of two NH units, have been obtained by a post-synthetic modification of the ...21,23-ditelluraporphyrin. The products of the tellurium-to-metal exchange, 21-platina-23-rhodaporphyrins, incorporate rhodacyclopentadiene and platinacyclopentadiene units with the metal atoms facing each other. The two molecules exhibit different degrees of metal-metal interaction depending on the oxidation state of platinum, with the NBO bond order being 0.04 for platinum(
iv
) and 0.15 for platinum(
ii
). Consistently, the Quantum Theory of Atoms in Molecules analysis revealed the presence of the bond determinant, the bond critical point, in the platinum(
ii
) species, in contrast to the platinum(
iv
) congener. The two porphyrinoids are interconvertible in redox reactions. They both exhibit fluxional behaviour in solution, studied by
1
H NMR, involving alteration in the metal ion coordination sphere accompanied by the macrocyclic skeleton conformation change.
The incorporation of two metal atoms, platinum and rhodium, in place of NH groups in the porphyrin macrocycle brings them into close proximity. The weak metal-metal bond is switched on and off depending on the oxidation state of the platinum ion.
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.
Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems ...discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller–Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers’ stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm−1–1700 cm−1 and 2300 cm−1–3400 cm−1 in the gas phase and 600 cm−1–1800 cm−1 and 2200 cm−1–3400 cm−1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm−1–1700 cm−1 and 2300 cm−1–3300 cm−1 for the gas phase and one broad absorption region in the solid state between 700 cm−1 and 3100 cm−1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.
Intra- and intermolecular forces competition was investigated in the 9,10-anthraquinone (
) and its derivatives both in vacuo and in the crystalline phase. The 1,8-dihydroxy-9,10-anthraquinone (
) ...and 1,8-dinitro-4,5-dihydroxy-anthraquinone (
) contain Resonance-Assisted Hydrogen Bonds (RAHBs). The intramolecular hydrogen bonds properties were studied in the electronic ground and excited states employing Møller-Plesset second-order perturbation theory (MP2), Density Functional Theory (DFT) method in its classical formulation as well as its time-dependent extension (TD-DFT). The proton potential functions were obtained via scanning the OH distance and the dihedral angle related to the OH group rotation. The topological analysis was carried out on the basis of theories of Atoms in Molecules (AIM-molecular topology, properties of critical points, AIM charges) and Electron Localization Function (ELF-2D maps showing bonding patterns, calculation of electron populations in the hydrogen bonds). The Symmetry-Adapted Perturbation Theory (SAPT) was applied for the energy decomposition in the dimers. Finally, Car-Parrinello molecular dynamics (CPMD) simulations were performed to shed light onto bridge protons dynamics upon environmental influence. The vibrational features of the OH stretching were revealed using Fourier transformation of the autocorrelation function of atomic velocity. It was found that the presence of OH and NO2 substituents influenced the geometric and electronic structure of the anthraquinone moiety. The AIM and ELF analyses showed that the quantitative differences between hydrogen bonds properties could be neglected. The bridged protons are localized on the donor side in the electronic ground state, but the Excited-State Intramolecular Proton Transfer (ESIPT) was noticed as a result of the TD-DFT calculations. The hierarchy of interactions determined by SAPT method indicated that weak hydrogen bonds play modifying role in the organization of these crystal structures, but primary ordering factor is dispersion. The CPMD crystalline phase results indicated bridged proton-sharing in the compound
.
Substitution is well-known to modulate the physico-chemical properties of molecules. In this study, a combined, multifactor approach was employed to determine a plethora of substitution patterns ...using –Br and –O-H in 1,4-naphthoquinone and its derivatives. On the basis of classical Density Functional Theory (DFT), 25 models divided into three groups were developed. The first group contains 1,4-naphthoquinone and its derivatives substituted only by –Br. The second group consists of compounds substituted by –Br and one –O-H group. As a result of the substitution, an intramolecular hydrogen bond was formed. The third group also contains –Br as a substituent, but two –O-H groups were introduced and two intramolecular hydrogen bonds were established. The simulations were performed at the ωB97XD/6-311++G(2d,2p) level of theory. The presence of substituents influenced the electronic structure of the parent compound and its derivatives by inductive effects, but it also affected the geometry of the 2 and 3 groups, due to the intramolecular hydrogen bonding and the formation of a quasi-ring/rings. The static DFT models were applied to investigate the aromaticity changes in the fused rings based on the Harmonic Oscillator Model of Aromaticity (HOMA). The OH stretching was detected for the compounds from groups 2 and 3 and further used to find correlations with energetic parameters. The evolution of the electronic structure was analyzed using Hirshfeld atomic charges and the Substituent Active Region (cSAR) parameter. The proton reaction path was investigated to provide information on the modulation of hydrogen bridge properties by diverse substitution positions on the donor and acceptor sides. Subsequently, Car–Parrinello Molecular Dynamics (CPMD) was carried out in the double-bridged systems (group 3) to assess the cooperative effects in double –O-H-substituted systems. It was determined that the –O-H influence on the core of the molecule is more significant than that of –Br, but the latter has a major impact on the bridge dynamics. The competitive or synergic effect of two –Br substituents was found to depend on the coupling between the intramolecular hydrogen bridges. Thus, the novel mechanism of a secondary (cooperative) substituent effect was established in the double-bridged systems via DFT and CPMD results comparison, consisting of a mediation of the bromine substitutions’ influence by the cooperative proton transfer events in the hydrogen bridges.
Display omitted
•The renewed and corrected interpretations of experimental vibrational spectra and hydrogen bonds in the L- and D-iso ascorbic acids.•The d(OO) = f(γ(OH)) correlation is a tool for ...verification of the dissimilarity between the resonance-assisted hydrogen bond (RAHB) and the charge assisted hydrogen bond (CAHB).•The difference between the bridged proton dynamics of the intramolecular and intermolecular hydrogen bonds in ascorbic acids shown by CPMD and DFT calculations.
The studies of two isomers of ascorbic acid and their deuteroanalogues, presented in the paper, have been accomplished by vibrational spectroscopy methods and quantum-chemical simulations. The spectroscopic research of L-ascorbic and D-isoascorbic acids have been carried out by the infrared (IR) and Raman (R) techniques. On the basis of the obtained results the spectral interpretation of the hydrogen bonded groups of ascorbic acids has been performed. Car-Parrinello Molecular Dynamics (CPMD) and Density Functional Theory (DFT) have been employed to support spectroscopic experimental findings and shed light onto the bridged proton dynamics in the L- and D- isomers of ascorbic acids. The accurate assignments of the hydrogen bond modes have been accomplished with the application of deuterosubstitution, CPMD-solid state simulations and Potential Energy Distribution (PED) analysis. The spectral and structural results have shown that dependency ν(OH) = f(γ(OH)) is the most common for the OHO hydrogen bond, whereas dependency d(OO) = f(γ(OH)) differs as for the ionic and resonance assisted hydrogen bonds.
This study covers the analysis of isomeric forms of nitrophthalic acids with pyridine. This work dwells on the complementary experimental (X-ray, IR and Raman) and theoretical (Car-Parrinello ...Molecular Dynamics (CPMD) and Density Functional Theory (DFT)) studies of the obtained complexes. The conducted studies showed that steric repulsion between the nitro group in ortho-position and the carboxyl group causes significant isomeric changes. Modeling of the nitrophthalic acid-pyridine complex yielded a short strong intramolecular hydrogen bond (SSHB). The transition energy from the isomeric form with an intermolecular hydrogen bond to the isomeric form with an intramolecular hydrogen bond was estimated.
PDF
