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

Two 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

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

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)

Structural and electronic structure differences due to the O-H center dot center dot center dot O and O-H center dot center dot center dot S bond formation in selected benzamide derivatives: a first-principles molecular dynamics study

Density functional theory-based methods were employed to obtain static and dynamical descriptions of the molecular properties of 2-hydroxy-N-methylbenzamide and 2-hydroxy-N-methylthiobenzamide; compounds containing O-H center dot center dot center dot O and O-H center dot center dot center dot S strong, intramolecular hydrogen bonds. These compounds are important as analogues of commercial analgesic and antipyretic medicines. In the current study the classical Kohn-Sham method was applied to develop static models describing the geometric parameters and proton potentials. The topological analysis of the electron density was performed via atoms in molecules theory. Subsequently, Car-Parrinello molecular dynamics investigations were performed in vacuo and in the solid state. The geometric and spectroscopic properties were investigated and compared with available experimental data. The influence of quantum effects on the intramolecular hydrogen bond properties were studied via path integral molecular dynamics in the solid state for 2-hydroxy-N-methylbenzamide. We found that the proton behavior depends strongly on the type of acceptor: the sulfur-containing bridge has significantly smaller proton flexibility than the oxygen-bearing analogue, which is reflected in the electronic structure and bridge dynamics