Polarizable molecular mechanics studies of Cu(I)/Zn(II) superoxide dismutase: bimetallic binding site and structured
Nohad Gresh1, Krystel El Hage, David Perahia
1Chemistry and Biology, Nucleo(s)tides and Immunology for Therapy (CBNIT), UMR 8601 CNRS, UFR Biomédicale, Paris, France.
Abstract:
The existence of a network of structured waters in the vicinity of the bimetallic site of Cu/Zn-superoxide dismutase (SOD) has been inferred from high-resolution X-ray crystallography. Long-duration molecular dynamics (MD) simulations could enable to quantify the lifetimes and possible interchanges of these waters between themselves as well as with a ligand diffusing toward the bimetallic site. The presence of several charged or polar ligands makes it necessary to resort to second-generation polarizable potentials. As a first step toward such simulations, we benchmark in this article the accuracy of one such potential, sum of interactions between fragments Ab initio computed (SIBFA), by comparisons with quantum mechanics (QM) computations. We first consider the bimetallic binding site of a Cu/Zn-SOD, in which three histidines and a water molecule are bound to Cu(I) and three histidines and one aspartate are bound to Zn(II). The comparisons are made for different His6 complexes with either one or both cations, and either with or without Asp and water. The total net charges vary from zero to three. We subsequently perform preliminary short-duration MD simulations of 296 waters solvating Cu/Zn-SOD. Six representative geometries are selected and energy-minimized. Single-point SIBFA and QM computations are then performed in parallel on model binding sites extracted from these six structures, each of which totals 301 atoms including the closest 28 waters from the Cu metal site. The ranking of their relative stabilities as given by SIBFA is identical to the QM one, and the relative energy differences by both approaches are fully consistent. In addition, the lowest-energy structure, from SIBFA and QM, has a close overlap with the crystallographic one. The SIBFA calculations enable to quantify the impact of polarization and charge transfer in the ranking of the six structures. Five structural waters, which connect Arg141 and Glu131, are endowed with very high dipole moments (2.7-3.0 Debye), equal and larger than the one computed by SIBFA in ice-like arrangements (2.7 D).
More Related Videos
09:22Synthesis and Performance Evaluations of ZnCoS/ZnCdS with Twin Crystal Structure for Multifunctional Redox Photocatalysis in Energy Applications
Published on: July 25, 2025
06:53Author Spotlight: Magnetometric Characterization of Intermediates in the Solid-State Electrochemistry of Redox-Active Metal-Organic Frameworks
Published on: June 9, 2023
Related Concept Videos
02:58Crystal Field Theory - Octahedral Complexes
To explain the observed behavior of transition metal complexes (such as colors), a model involving electrostatic interactions between the electrons from the ligands and the electrons in the unhybridized d orbitals of the central metal atom has been developed. This electrostatic model is crystal field theory (CFT). It helps to understand, interpret, and predict the colors, magnetic behavior, and some structures of coordination compounds of transition metals.
CFT focuses on...
03:02Colors and Magnetism
When atoms or molecules absorb light at the proper frequency, their electrons are excited to higher-energy orbitals. For many main group atoms and molecules, the absorbed photons are in the ultraviolet range of the electromagnetic spectrum, which cannot be detected by the human eye. For coordination compounds, the energy difference between the d orbitals often allows photons in the visible range to be absorbed and emitted, which is seen as colors by the human...
02:46Crystal Field Theory - Tetrahedral and Square Planar Complexes
Crystal field theory (CFT) is applicable to molecules in geometries other than octahedral. In octahedral complexes, the lobes of the dx2−y2 and dz2 orbitals point directly at the ligands. For tetrahedral complexes, the d orbitals remain in place, but with only four ligands located between the axes. None of the orbitals points directly at the tetrahedral ligands. However, the dx2−y2 and dz2 orbitals (along the Cartesian axes) overlap with the ligands less than...
03:45Formation of Complex Ions