Dynamic Investigation of Protein Metal Active Sites: Interplay of XANES and Molecular Dynamics Simulations

Structure-Redox Response Correlation in a Few Select Heme Systems Using X-ray Absorption Spectroelectrochemistry

Heme based biomolecules control some of the most crucial life processes, such as oxygen and electron transport during respiration and energy metabolism, respectively. The active site of the heme, viz., the metal center, plays a key role and attributes functionality to these biomolecules. During the oxygen binding and debinding processes, it is important to note that the oxidation state of iron in hemoglobin (+II in the native form) does not undergo any change. However, the spin states of the metal center change. We present here a comprehensive study of the redox response of such molecules, based on the electronic structure of the active site. The local electronic structure of heme in a few selective molecular systems is studied in operando via synchrotron X-ray absorption spectroscopy (Fe K-edge) and cyclic voltammetry. Our objective is to identify the electronic structural parameters that can effectively be correlated with the redox reversibility. Evolution in these parameters can be followed to trace the overall changes in redox state of the system. Our data indicate that axial coordination and spin state of the iron center are two such parameters that are intimately connected with the redox response.

Ligand pathways in neuroglobin revealed by low-temperature photodissociation and docking experiments

A combined biophysical approach was applied to map gas-docking sites within murine neuroglobin (Ngb), revealing snapshots of events that might govern activity and dynamics in this unique hexacoordinate globin, which is most likely to be involved in gas-sensing in the central nervous system and for which a precise mechanism of action remains to be elucidated. The application of UV-visible microspectroscopy in crystallo, solution X-ray absorption near-edge spectroscopy and X-ray diffraction experiments at 15-40 K provided the structural characterization of an Ngb photolytic intermediate by cryo-trapping and allowed direct observation of the relocation of carbon monoxide within the distal heme pocket after photodissociation. Moreover, X-ray diffraction at 100 K under a high pressure of dioxygen, a physiological ligand of Ngb, unravelled the existence of a storage site for O2 in Ngb which coincides with Xe-III, a previously described docking site for xenon or krypton. Notably, no other secondary sites were observed under our experimental conditions.

Dynamic multiple-scattering treatment of X-ray absorption: Parameterization of a new molecular dynamics force field for myoglobin

We present a detailed analysis of the X-ray absorption near-edge structure (XANES) data on the Fe K-edge of CO Myoglobin based on a combined procedure of Molecular Dynamics (MD) calculations and MXAN (Minuit XANes) data analysis that we call D-MXAN. The ability of performing quantitative XANES data analysis allows us to refine classical force field MD parameters, thus obtaining a reliable tool for the atomic investigation of this important model system for biological macromolecules. The iterative procedure here applied corrects the greatest part of the structural discrepancy between classical MD sampling and experimental determinations. Our procedure, moreover, is able to discriminate between different heme conformational basins visited during the MD simulation, thus demonstrating the necessity of a sampling on the order of tens of nanoseconds, even for an application such X-ray absorption spectroscopy data analysis.

Equilibrium between 5- and 6-Fold Coordination in the First Hydration Shell of Cu(II)

The hydration structure dynamics of Cu(II) ion is characterized by a combination of classical molecular dynamics simulation and X-ray absorption near-edge spectroscopy. Previous experimental data have been analyzed on the basis of 5- or 6-fold first hydration structure, with a quite well-established equatorial structure. This 4-fold equatorial geometry has been our starting point to develop a simple but effective in silico model, which provides ab initio theoretical X-ray absorption spectra in very good agreement with the experimental data. Our results point out two equally populated 6- and 5-fold hydration structures with remarkable different water residence times of 5 and 98 ps, respectively, and a low free energy barrier between first and second hydration shell.

Identification of catalytic sites for oxygen reduction in iron- and nitrogen-doped graphene materials

While platinum has hitherto been the element of choice for catalysing oxygen electroreduction in acidic polymer fuel cells, tremendous progress has been reported for pyrolysed Fe-N-C materials. However, the structure of their active sites has remained elusive, delaying further advance. Here, we synthesized Fe-N-C materials quasi-free of crystallographic iron structures after argon or ammonia pyrolysis. These materials exhibit nearly identical Mössbauer spectra and identical X-ray absorption near-edge spectroscopy (XANES) spectra, revealing the same Fe-centred moieties. However, the much higher activity and basicity of NH3-pyrolysed Fe-N-C materials demonstrates that the turnover frequency of Fe-centred moieties depends on the physico-chemical properties of the support. Following a thorough XANES analysis, the detailed structures of two FeN4 porphyrinic architectures with different O2 adsorption modes were then identified. These porphyrinic moieties are not easily integrated in graphene sheets, in contrast with Fe-centred moieties assumed hitherto for pyrolysed Fe-N-C materials. These new insights open the path to bottom-up synthesis approaches and studies on site-support interactions.

Ultrafast core-loss spectroscopy in four-dimensional electron microscopy

We demonstrate ultrafast core-electron energy-loss spectroscopy in four-dimensional electron microscopy as an element-specific probe of nanoscale dynamics. We apply it to the study of photoexcited graphite with femtosecond and nanosecond resolutions. The transient core-loss spectra, in combination with ab initio molecular dynamics simulations, reveal the elongation of the carbon-carbon bonds, even though the overall behavior is a contraction of the crystal lattice. A prompt energy-gap shrinkage is observed on the picosecond time scale, which is caused by local bond length elongation and the direct renormalization of band energies due to temperature-dependent electron-phonon interactions.

K-edge XANES investigation of octakis(DMSO)lanthanoid(III) complexes in DMSO solution and solid iodides

The potential of high energy XANES (X-ray absorption near edge structure) as a tool for the structural analysis of lanthanoid-containing systems has been explored. The K-edge XANES spectra of La(3+), Gd(3+), and Lu(3+) ions both in DMSO solution and solid octakis(DMSO)lanthanoid(III) iodides have been analysed. Although the K-edges of lanthanoids cover the energy range of 38 (La) to 65 (Lu) keV, the large widths of the core hole states do not appreciably reduce the potential structural information of the XANES data. We show that, for lanthanoid compounds, accurate structural parameters are obtained from the analysis of K-edge XANES signals if a deconvolution procedure is carried out. We found that in solid octakis(DMSO)lanthanoid(III) iodides the Ln(3+) ions are coordinated by eight DMSO ligands arranged in a quite symmetric fashion. In DMSO solution the Ln(3+) ions retain a regular eight-coordination structure and the coordination number does not change along the series. In contrast to when in water the second coordination shell has been found to provide a negligible contribution to the XANES spectra of Ln(3+) ions in DMSO solution.

Effects of the pathological Q212P mutation on human prion protein non-octarepeat copper-binding site

Prion diseases are a class of fatal neurodegenerative disorders characterized by brain spongiosis, synaptic degeneration, microglia and astrocytes activation, neuronal loss and altered redox control. These maladies can be sporadic, iatrogenic and genetic. The etiological agent is the prion, a misfolded form of the cellular prion protein, PrP(C). PrP(C) interacts with metal ions, in particular copper and zinc, through the octarepeat and non-octarepeat binding sites. The physiological implication of this interaction is still unclear, as is the role of metals in the conversion. Since prion diseases present metal dyshomeostasis and increased oxidative stress, we described the copper-binding site located in the human C-terminal domain of PrP-HuPrP(90-231), both in the wild-type protein and in the protein carrying the pathological mutation Q212P. We used the synchrotron-based X-ray absorption fine structure technique to study the Cu(II) and Cu(I) coordination geometries in the mutant, and we compared them with those obtained using the wild-type protein. By analyzing the extended X-ray absorption fine structure and the X-ray absorption near-edge structure, we highlighted changes in copper coordination induced by the point mutation Q212P in both oxidation states. While in the wild-type protein the copper-binding site has the same structure for both Cu(II) and Cu(I), in the mutant the coordination site changes drastically from the oxidized to the reduced form of the copper ion. Copper-binding sites in the mutant resemble those obtained using peptides, confirming the loss of short- and long-range interactions. These changes probably cause alterations in copper homeostasis and, consequently, in redox control.

Identification of 13- and 14-coordinated structures of first hydrated shell of [AuCl4](-) acid aqueous solution by combination of MD and XANES

Molecular dynamics simulations combined with multiple scatting X-ray absorption spectral calculations have been employed to study the local geometry of an [AuCl4](-) cluster in acid aqueous solution. It is found that the previously assumed simple octahedral structure with two H atoms located in the axial position of the [AuCl4](-) plane fails to correctly reproduce the experimental X-ray absorption spectra. Molecular dynamics simulations have revealed a very complicated first hydrated water shell that mainly consists of 13 or 14 ligand water molecules. In general, the first hydrated water molecules are located in two orthogonal quasi-elliptic surfaces around [AuCl4](-) cluster. The water molecules' configuration in the first hydrated shell is further confirmed by the X-ray absorption spectral analysis, which gives a good agreement with the experiment when the 13-hydrated and 14-hydrated complexes are evenly mixed. It shows the power of the combined molecular dynamics simulations and spectra calculations in determining the local structure of molecular ligands around the metal center.

Carbon monoxide binding to the heme group at the dimeric interface modulates structure and copper accessibility in the Cu,Zn superoxide dismutase from Haemophilus ducreyi: in silico and in vitro evidences

X-ray absorption near-edge structure (XANES) spectroscopy and molecular dynamics (MD) simulations have been jointly applied to the study of the Cu,Zn superoxide dismutase from Haemophilus ducreyi (HdSOD) in interaction with the carbon monoxide molecule. The configurational flexibility of the Fe(II)-heme group, intercalated between the two subunits, has been sampled by MD simulations and included in the XANES data analysis without optimization in the structural parameter space. Our results provide an interpretation of the observed discrepancy in the Fe-heme distances as detected by extended X-ray absorption fine structure (EXAFS) spectroscopy and the classical XANES analysis, in which the structural parameters are optimized in a unique structure. Moreover, binding of the CO molecule to the heme induces a long range effect on the Cu,Zn active site, as evidenced by both MD simulations and in vitro experiments. MD simulation of the CO bound system, in fact, highlighted a structural rearrangement of the protein-protein hydrogen bond network in the region of the Cu,Zn active site, correlated with an increase in water accessibility at short distance from the copper atom. In line, in vitro experiments evidenced an increase of copper accessibility to a chelating agent when the CO molecule binds to the heme group, as compared to a heme deprived HdSOD. Altogether, our results support the hypothesis that the HdSOD is a heme-sensor protein, in which binding to small gaseous molecules modulates the enzyme superoxide activity as an adaptive response to the bacterial environment.