Aneta Jezierska-Mazzarello

Molecular property investigations of an ortho-hydroxy Schiff base type compound with the first-principle molecular dynamics approach.

The structure, proton transfer, and vibrational dynamics under ambient conditions of a selected ortho-hydroxy Schiff base type compound, 2-(N-methyl-alpha-iminoethyl)-4-chlorophenol, containing a very short intramolecular hydrogen bond, were investigated computationally in the gas phase and in the crystal by density functional theory (DFT) based first-principle molecular dynamics (FPMD). It is found that the proton is well localized on the nitrogen side of the O...H...N bridge in the crystal phase, in agreement with X-ray diffraction experiments, while a more labile proton is located most of the time on the oxygen side in a vacuum. Environmental effects on this very strong hydrogen bond thus appear crucial and lead to drastic changes of the infrared (IR) spectrum: The computed gas-phase IR spectrum shows a very broad absorption band that covers frequencies from about 1000 to 3000 cm(-1) assigned to the labile proton. In mere contrast, a much more localized absorption band around 2600-2700 cm(-1) is predicted in the crystal phase. Finally, effects of the quantization of the proton motion on the hydrogen bond structure were estimated in two ways. First, we constructed the one-dimensional (1D) potential energy surface (PES) for the proton along the O...H...N bridge in a vacuum. The 1D Schrodinger equation was then solved. Next, path integral molecular dynamics (PIMD) was performed in the solid state. Inclusion of quantum effects does not affect the observed change of the most probable tautomer, upon going from the gas phase to the crystal.

Hydrogen bonding as a modulator of aromaticity and electronic structure of selected ortho-hydroxybenzaldehyde derivatives.

Properties of hydrogen bonds can induce changes in geometric or electronic structure parameters in the vicinity of the bridge. Here, we focused primarily on the influence of intramolecular H-bonding on the molecular properties in selected ortho-hydroxybenzaldehydes, with additional restricted insight into substituent effects. Static models were obtained in the framework of density functional theory at B3LYP/6-311+G(d,p) level. The electronic structure parameters evolution was analyzed on the basis of Atoms In Molecules (AIM) and Natural Bond Orbitals methods. The aromaticity changes related to the variable proton position and presence of substituents were studied using Harmonic Oscillator Model of Aromaticity (HOMA), Nucleus-Independent Chemical Shift (NICS) and AIM-based parameter of Matta and Hernández-Trujillo. Finally, Car-Parrinello molecular dynamics was applied to study variability of the hydrogen bridge dynamics. The interplay between effects of the substitution and variable position of the bridged proton was discussed. It was found that the hydrogen bond energies are ca. 9-10 kcal/mol, and the bridged proton exhibits some degree of penetration into the acceptor region. The covalent character of the studied hydrogen bond was most observable when the bridged proton reached the middle position between the donor and acceptor regions. The aromaticity indexes showed that the aromaticity of the central phenyl ring is strongly dependent on the bridged proton position. Correlations between these parameters were found and discussed. In the applied time-scale, the analysis of time evolution of geometric parameters showed that the resonance strengthening does not play a crucial role in the studied compounds.

Investigation of 6-fluoroquinolones activity against Mycobacterium tuberculosis using theoretical molecular descriptors: a case study

AbstractA quantitative structure-activity relationship (QSAR) study on a set of 66 structurally-similar 6-fluoroquinolones was performed using a large pool of theoretical molecular descriptors. Ab initio geometry optimizations were carried out to reproduce the geometrical and electronic structure parameters. The resulting molecular structures were confirmed to be minima via harmonic frequency calculations. Obtained atomic charges, HOMO and LUMO energies, orbital electron densities, dipole moment, energy and many other properties served as quantum-chemical descriptors. A multiple linear regression (MLR) technique was applied to generate a linear model for predicting the biological activity, Minimal Inhibitory Concentration (MIC), treated as negative decade logarithm, (pMIC). The heuristic method was used to optimize the model parameters and select the most significant descriptors. The model was tested internally using the CV LOO procedure on the training set and validated against the external validation set. The result (Q2ext = 0.7393), which was obtained on an external, previously excluded validation data set, shows the predictive performances of this model (R2tr = 0.7416, Q2tr = 0.6613) in establishing (Q)SAR of 6-fluoroquinolones. This validated model could be proficiently used to design new 6-fluoroquinolones with possible higher activity.

Interference of H-bonding and substituent effects in nitro- and hydroxy-substituted salicylaldehydes

AbstractTwo intramolecular interactions, i.e., (1) hydrogen bond and (2) substituent effect, were analyzed and compared. For this purpose, the geometry of 4- and 5-X-substituted salicylaldehyde derivatives (X = NO2, H or OH) was optimized by means of B3LYP/6-311 + G(d,p) and MP2/aug-cc-pVDZ methods. The results obtained allowed us to show that substituents (NO2 or OH) in the para or meta position with respect to either OH or CHO in H-bonded systems interact more strongly than in the case of di-substituted species: 4- and 3-nitrophenol or 4- and 3-hydroxybenzaldehyde by ∼31%. The substituent effect due to the intramolecular charge transfer from the para-counter substituent (NO2) to the proton-donating group (OH) is ∼35% greater than for the interaction of para-OH with the proton-accepting group (CHO). The total energy of H-bonding for salicylaldehyde, and its derivatives, is composed of two contributions: ∼80% from the energy of H-bond formation and ∼20% from the energy associated with reorganization of the electron structure of the systems in question. FigureSubstituent effect stabilization energy (SESE) estimation for the salicylaldehyde and its 4- and 5-X-substituted derivatives

Direct observation of the substitution effects on the hydrogen bridge dynamics in selected Schiff bases—A comparative molecular dynamics study

We have studied substituent effects on the properties of the intramolecular hydrogen bond of some ortho-hydroxy Schiff bases using density functional theory (DFT) based first-principle molecular dynamics (FPMD) and path integral molecular dynamics. The studied compounds possess a strong intramolecular hydrogen bond (r((O⋅⋅⋅N)) ≤ 2.6 Å), which can be tuned by substitution to either (i) enhance the basicity of the acceptor moiety by induction effects or (ii) decrease the hydrogen bond length through steric repulsion. DFT calculations and FPMD were employed to investigate structural and dynamical properties of the selected molecules, while quantum effects on the structural properties were assessed using path integral FPMD. The simulations were performed in vacuo and in the solid state to study the influence of the environment on the hydrogen bond and spectroscopic properties. We give computational support to the suggestion that induction effects are less effective to tune the intramolecular hydrogen bond properties of the discussed ortho-hydroxy Schiff bases than the steric or the environmental effects.

Impact of Mercury(II) on proteinase K catalytic center: investigations via classical and Born-Oppenheimer molecular dynamics

Mercury(II) has a strong affinity for the thiol groups in proteins often resulting in the disruption of their biological functions. In this study we present classical and first-principles, DFT-based molecular dynamics (MD) simulations of a complex of Hg(II) and proteinase K, a well-known serine protease with a very broad and diverse enzymatic activity. It contains a catalytic triad formed by Asp39, His69, and Ser224, which is responsible for its biological activity. It was found previously by X-ray diffraction experiments that the presence of Hg(II) inhibits the enzymatic action of proteinase K by affecting the stereochemistry of the triad. Our simulations predict that (i) the overall structure as well as the protein backbone dynamics are only slightly affected by the mercury cation, (ii) depending on the occupied mercury site, the hydrogen bonds of the catalytic triad are either severely disrupted (both bonds for mercury at site 1, and the His69–Ser224 contact for mercury at site 2) or slightly strengthened (the Asp39–His69 bond when mercury is at site 2), (iii) the network of hydrogen bonds of the catalytic triad is not static but undergoes constant fluctuations, which are significantly modified by the presence of the Hg(II) cation, influencing in turn the triad’s ability to carry out the enzymatic function—these facts explain the experimental findings on the inhibition of proteinase K by Hg(II).

Flexibility of a biotinylated ligand in artificial metalloenzymes based on streptavidin—an insight from molecular dynamics simulations with classical and ab initio force fields

In the field of enzymatic catalysis, creating activity from a non catalytic scaffold is a daunting task. Introduction of a catalytically active moiety within a protein scaffold offers an attractive means for the creation of artificial metalloenzymes. With this goal in mind, introduction of a biotinylated d6-piano-stool complex within streptavidin (SAV) affords enantioselective artificial transfer-hydrogenases for the reduction of prochiral ketones. Based on an X-ray crystal structure of a highly selective hybrid catalyst, displaying significant disorder around the biotinylated catalyst [η6-(p-cymene)Ru(Biot-p-L)Cl], we report on molecular dynamics simulations to shed light on the protein–cofactor interactions and contacts. The results of these simulations with classical force field indicate that the SAV-biotin and SAV-catalyst complexes are more stable than ligand-free SAV. The point mutations introduced did not affect significantly the overall behavior of SAV and, unexpectedly, the P64G substitution did not provide additional flexibility to the protein scaffold. The metal-cofactor proved to be conformationally flexible, and the S112K or P64G mutants proved to enhance this effect in the most pronounced way. The network of intermolecular hydrogen bonds is efficient at stabilizing the position of biotin, but much less at fixing the conformation of an extended biotinylated ligand. This leads to a relative conformational freedom of the metal-cofactor, and a poorly localized catalytic metal moiety. MD calculations with ab initio potential function suggest that the hydrogen bonds alone are not sufficient factors for full stabilization of the biotin. The hydrophobic biotin-binding pocket (and generally protein scaffold) maintains the hydrogen bonds between biotin and protein.

Influence of environmental humidity on organization and molecular dynamics of heteromacrocyclic assemblies.

1D and 2D NMR study, Car-Parrinello molecular dynamics, as well as classical molecular dynamics were employed to investigate three derivatives of benzodiazacoronands (achiral compounds which are able to form single crystals with a planar chirality) with intention to explain all subtle effects important during their preorganization, the step anticipating formation of crystals. The experimental study was carried out in two solvents: chloroform and DMSO either containing traces of water (commercial samples) or carefully dried over molecular sieves. Both methods revealed that environmental humidity has a dramatic influence on topology of solute-solvent interactions. Damping of the macrocycle dynamics by its diverse types of interactions with water molecules was shown by computational means. In the most spectacular experiment, we have proved that in chloroform-d during the low temperature measurements traces of water dramatically change the spectral pattern, leading to isochronous NMR signals of the AB spin system of benzodiazacoronand. The temperature of isochronous point (TIP) strongly depends on the benzodiazacoronand/water (BW) ratio. This observation opens a pathway to a new strategy based on variable temperature crystallizations and fitting of BW ratio with hope to optimize conditions for formation of chiral crystals.

Comparison of resonance assisted and charge assisted effects in strengthening of hydrogen bonds in dipyrrins.

This paper deals with the study of two types of hydrogen bonding: a quasi-aromatic hydrogen bonding in dipyrromethene and the ionic one in dipyrromethane. The study focuses on two phenomena-the proton transfer process and tautomeric equilibrium. Metric parameters and spectroscopic assignments have been calculated; this allowed a further comparison of spectral features calculated with four methods (Car-Parrinello molecular dynamics (CPMD), ab initio, density functional theory (DFT), and numerical calculation of anharmonic vibrational levels via a solution of the corresponding 1D Schrödinger equation). A significant dynamics of the bridged proton and bent vibration of pyrrole fragments in dipyrromethane have been exposed by the CPMD calculations. The prevailing of the ionic effect over the π-electronic coupling in the strengthening of the hydrogen bonding has been shown on the basis of the calculated structural, electron-topological, and spectral data as well as potential energy surface (PES). The analysis of the aromaticity and electronic state of pyrrole and chelate moieties depending on the tautomeric equilibrium by the quantum theory of atoms in molecules (QTAIM) method was conducted. The principle divergence in the behavior of aromaticity of the chelate chains in the analyzed compounds was demonstrated.

Impact of proton transfer phenomena on the electronic structure of model Schiff bases: An AIM/NBO/ELF study

Understanding of the electronic structure evolution due to a proton dynamics is a key issue in biochemistry and material science. This paper reports on density functional theory calculations of Schiff bases containing short, strong intramolecular hydrogen bonds where the bridged proton is located: (i) at the donor site, (ii) strongly delocalized, and (iii) at the acceptor site. The mobility of the bridged proton and its influence on the molecular structure and properties of the chosen Schiff base derivatives have been investigated on the basis of Atoms in Molecules, Natural Bond Orbitals, and Electron Localization Function theories. It has been observed that the extent of the bridged proton delocalization is strongly modified by the steric and inductive effects present in the studied compounds introduced by various substituents. It has been shown that: (i) potential energy profiles for the proton motion are extremely dependent on the substitution of the aromatic ring, (ii) the topology of the free electron pairs present at the donor∕acceptor site, as well as their electron populations, are affected qualitatively by the bridged proton position, (iii) the distortion of the molecular structure due to the bridged proton dynamics includes the atomic charge fluctuations, which are in some cases non-monotonic, and (iv) topology of the ELF recognizes events of proton detachment from the donor and attachment to the acceptor. The quantitative and qualitative results shed light onto molecular consequences of the proton transfer phenomena.