On the Oligomeric State of DJ-1 Protein and Its Mutants Associated with Parkinson Disease A COMBINED COMPUTATIONAL AND IN VITRO STUDY

Mutations in the DJ-1 protein are present in patients suffering from familial Parkinson disease. Here we use computational methods and biological assays to investigate the relationship between DJ-1 missense mutations and the protein oligomeric state. Molecular dynamics calculations suggest that: (i) the structure of DJ-1 wild type (WT) in aqueous solution, in both oxidized and reduced forms, is similar to the crystal structure of the reduced form; (ii) the Parkinson disease-causing M26I variant is structurally similar to the WT, consistent with the experimental evidence showing the protein is a dimer as WT; (iii) R98Q is structurally similar to the WT, consistent with the fact that this is a physiological variant; and (iv) the L166P monomer rapidly evolves toward a conformation significantly different from WT, suggesting a change in its ability to oligomerize. Our combined computational and experimental approach is next used to identify a mutant (R28A) that, in contrast to L166P, destabilizes the dimer subunit-subunit interface without significantly changing secondary structure elements

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

Exploring the Dynamical Nature of Intermolecular Hydrogen Bonds in Benzamide, Quinoline and Benzoic Acid Derivatives

The hydrogen bonds properties of 2,6-difluorobenzamide, 5-hydroxyquinoline and 4-hydroxybenzoic acid were investigated by Car–Parrinello and path integral molecular dynamics (CPMD and PIMD), respectively. The computations were carried out in vacuo and in the crystalline phase. The studied complexes possess diverse networks of intermolecular hydrogen bonds (N-H…O, O-H…N and O-H…O). The time evolution of hydrogen bridges gave a deeper insight into bonds dynamics, showing that bridged protons are mostly localized on the donor side; however, the proton transfer phenomenon was registered as well. The vibrational features associated with O-H and N-H stretching were analyzed on the basis of the Fourier transform of the atomic velocity autocorrelation function. The spectroscopic effects of hydrogen bond formation were studied. The PIMD revealed quantum effects influencing the hydrogen bridges providing more accurate free energy sampling. It was found that the N…O or O…O interatomic distances decreased (reducing the length of the hydrogen bridge), while the O-H or N-H covalent bond was elongated, which led to the increase in the proton sharing. Furthermore, Quantum Theory of Atoms in Molecules (QTAIM) was used to give insight into electronic structure parameters. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was employed to estimate the energy contributions to the interaction energy of the selected dimers

Role of Non-Covalent Interactions in Carbonic Anhydrase I—Topiramate Complex Based on QM/MM Approach

Carbonic anhydrase (CA) I with a Topiramate (TPM) complex was investigated on the basis of a Quantum Mechanics/Molecular Mechanics (QM/MM) approach. The QM part was treated using Density Functional Theory (DFT) while the MM was simulated using Amberff14SB and GAFF force fields. In addition, the TIP3P model was applied to reproduce the polar environment influence on the studied complex. Next, three snapshots (after 5 ps, 10 ps, and 15 ps of the simulation time) were taken from the obtained trajectory to provide an insight into the non-covalent interactions present between the ligand and binding pocket of the protein. Our special attention was devoted to the binding site rearrangement, which is known in the literature concerning the complex. This part of the computations was performed using ωB97X functional with Grimme D3 dispersion corrections as well as a Becke–Johnson damping function (D3-BJ). Two basis sets were applied: def2-SVP (for larger models) and def2-TZVPD (for smaller models), respectively. In order to detect and describe non-covalent interactions between amino acids of the binding pocket and the ligand, Independent Gradient Model based on Hirshfeld partitioning (IGMH), Interaction Region Indicator (IRI), Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Orbitals (NBO) methods were employed. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was applied for energy decomposition between the ligand and protein. It was found that during the simulation time, the ligand position in the binding site was preserved. Nonetheless, amino acids interacting with TPM were exchanging during the simulation, thus showing the binding site reorganization. The energy partitioning revealed that dispersion and electrostatics are decisive factors that are responsible for the complex stability

Characterization of intermolecular interactions : from dimers to microsolvation models

Intermodular interactions play an important role in many processes at the molecular level. In the contemporary science, there is a growing interest concerning the characteristics of such interactions. Therefore, the computational chemistry can provide answers to many questions, which could not be answered using experimental methods. The Symmetry-Adapted Perturbation Theory (SAPT) method was applied to characterize the energy partitioning in dimers, trimers and microsolvation models. The investigated complexes belong to various classes of compounds, e.g. • dimers of: NH3 ˑˑˑHX, HF-pyridine, cycloalkanes, hypohalous acids; • trimers of: NH3 ˑˑˑNH3ˑˑˑHF or NH3ˑˑˑHFˑˑˑHF; • microsolvation models (biotin - water molecules). The current study summarizes recent years of our research devoted to the intermolecular interactions