Predicting the Structure of Enzymes with Metal Cofactors: The Example of [FeFe] Hydrogenases.

The advent of deep learning algorithms for protein folding opened a new era in the ability of predicting and optimizing the function of proteins once the sequence is known. The task is more intricate when cofactors like metal ions or small ligands are essential to functioning. In this case, the combined use of traditional simulation methods based on interatomic force fields and deep learning predictions is mandatory. We use the example of [FeFe] hydrogenases, enzymes of unicellular algae promising for biotechnology applications to illustrate this situation. [FeFe] hydrogenase is an iron-sulfur protein that catalyzes the chemical reduction of protons dissolved in liquid water into molecular hydrogen as a gas. Hydrogen production efficiency and cell sensitivity to dioxygen are important parameters to optimize the industrial applications of biological hydrogen production. Both parameters are related to the organization of iron-sulfur clusters within protein domains. In this work, we propose possible three-dimensional structures of Chlorella vulgaris 211/11P [FeFe] hydrogenase, the sequence of which was extracted from the recently published genome of the given strain. Initial structural models are built using: (i) the deep learning algorithm AlphaFold; (ii) the homology modeling server SwissModel; (iii) a manual construction based on the best known bacterial crystal structure. Missing iron-sulfur clusters are included and microsecond-long molecular dynamics of initial structures embedded into the water solution environment were performed. Multiple-walkers metadynamics was also used to enhance the sampling of structures encompassing both functional and non-functional organizations of iron-sulfur clusters. The resulting structural model provided by deep learning is consistent with functional [FeFe] hydrogenase characterized by peculiar interactions between cofactors and the protein matrix.

Chasing the Elusive "In-Between" State of the Copper-Amyloid ß Complex by X-ray Absorption through Partial Thermal Relaxation after Photoreduction.

The redox activity of Cu ions bound to the amyloid-ß (Aß) peptide is implicated as a source of oxidative stress in the context of Alzheimer's disease. In order to explain the efficient redox cycling between CuII -Aß (distorted square-pyramidal) and CuI -Aß (digonal) resting states, the existence of a low-populated "in-between" state, prone to bind Cu in both oxidation states, has been postulated. Here, we exploited the partial X-ray induced photoreduction at 10 K, followed by a thermal relaxation at 200 K, to trap and characterize by X-ray Absorption Spectroscopy (XAS) a partially reduced Cu-Aß1-16 species different from the resting states. Remarkably, the XAS spectrum is well-fitted by a previously proposed model of the "in-between" state, hence providing the first direct spectroscopic characterization of an intermediate state. The present approach could be used to explore and identify the catalytic intermediates of other relevant metal complexes.

Modelling Protein Plasticity: The Example of Frataxin and Its Variants.

Frataxin (FXN) is a protein involved in storage and delivery of iron in the mitochondria. Single-point mutations in the FXN gene lead to reduced production of functional frataxin, with the consequent dyshomeostasis of iron. FXN variants are at the basis of neurological impairment (the Friedreich's ataxia) and several types of cancer. By using altruistic metadynamics in conjunction with the maximal constrained entropy principle, we estimate the change of free energy in the protein unfolding of frataxin and of some of its pathological mutants. The sampled configurations highlight differences between the wild-type and mutated sequences in the stability of the folded state. In partial agreement with thermodynamic experiments, where most of the analyzed variants are characterized by lower thermal stability compared to wild type, the D104G variant is found with a stability comparable to the wild-type sequence and a lower water-accessible surface area. These observations, obtained with the new approach we propose in our work, point to a functional switch, affected by single-point mutations, of frataxin from iron storage to iron release. The method is suitable to investigate wide structural changes in proteins in general, after a proper tuning of the chosen collective variable used to perform the transition.

Serial time-resolved crystallography of photosystem II using a femtosecond X-ray laser.

Photosynthesis, a process catalysed by plants, algae and cyanobacteria converts sunlight to energy thus sustaining all higher life on Earth. Two large membrane protein complexes, photosystem I and II (PSI and PSII), act in series to catalyse the light-driven reactions in photosynthesis. PSII catalyses the light-driven water splitting process, which maintains the Earth's oxygenic atmosphere. In this process, the oxygen-evolving complex (OEC) of PSII cycles through five states, S0 to S4, in which four electrons are sequentially extracted from the OEC in four light-driven charge-separation events. Here we describe time resolved experiments on PSII nano/microcrystals from Thermosynechococcus elongatus performed with the recently developed technique of serial femtosecond crystallography. Structures have been determined from PSII in the dark S1 state and after double laser excitation (putative S3 state) at 5 and 5.5 Å resolution, respectively. The results provide evidence that PSII undergoes significant conformational changes at the electron acceptor side and at the Mn4CaO5 core of the OEC. These include an elongation of the metal cluster, accompanied by changes in the protein environment, which could allow for binding of the second substrate water molecule between the more distant protruding Mn (referred to as the 'dangler' Mn) and the Mn3CaOx cubane in the S2 to S3 transition, as predicted by spectroscopic and computational studies. This work shows the great potential for time-resolved serial femtosecond crystallography for investigation of catalytic processes in biomolecules.

Atomic structure of granulin determined from native nanocrystalline granulovirus using an X-ray free-electron laser.

To understand how molecules function in biological systems, new methods are required to obtain atomic resolution structures from biological material under physiological conditions. Intense femtosecond-duration pulses from X-ray free-electron lasers (XFELs) can outrun most damage processes, vastly increasing the tolerable dose before the specimen is destroyed. This in turn allows structure determination from crystals much smaller and more radiation sensitive than previously considered possible, allowing data collection from room temperature structures and avoiding structural changes due to cooling. Regardless, high-resolution structures obtained from XFEL data mostly use crystals far larger than 1 µm3 in volume, whereas the X-ray beam is often attenuated to protect the detector from damage caused by intense Bragg spots. Here, we describe the 2 Å resolution structure of native nanocrystalline granulovirus occlusion bodies (OBs) that are less than 0.016 µm3 in volume using the full power of the Linac Coherent Light Source (LCLS) and a dose up to 1.3 GGy per crystal. The crystalline shell of granulovirus OBs consists, on average, of about 9,000 unit cells, representing the smallest protein crystals to yield a high-resolution structure by X-ray crystallography to date. The XFEL structure shows little to no evidence of radiation damage and is more complete than a model determined using synchrotron data from recombinantly produced, much larger, cryocooled granulovirus granulin microcrystals. Our measurements suggest that it should be possible, under ideal experimental conditions, to obtain data from protein crystals with only 100 unit cells in volume using currently available XFELs and suggest that single-molecule imaging of individual biomolecules could almost be within reach.

Multi-scale theoretical approach to X-ray absorption spectra in disordered systems: an application to the study of Zn(ii) in water.

We develop a multi-scale theoretical approach aimed at calculating from first principles X-ray absorption spectra of liquid solutions and disordered systems. We test the method by considering the paradigmatic case of Zn(ii) in water which, besides being relevant in itself, is also of interest for biology. With the help of classical molecular dynamics simulations we start by producing bunches of configurations differing for the Zn(ii)-water coordination mode. Different coordination modes are obtained by making use of the so-called dummy atoms method. From the collected molecular dynamics trajectories, snapshots of a more manageable subsystem encompassing the metal site and two solvation layers are cut out. Density functional theory is used to optimize and relax these reduced system configurations employing a uniform dielectric to mimic the surrounding bulk liquid water. On the resulting structures, fully quantum mechanical X-ray absorption spectra calculations are performed by including core-hole effects and core-level shifts. The proposed approach does not rely on any guessing or fitting of the force field or of the atomic positions of the system. The comparison of the theoretically computed spectrum with the experimental Zn K-edge XANES data unambiguously demonstrates that among the different a priori possible geometries, Zn(ii) in water lives in an octahedral coordination mode.

Visualizing a protein quake with time-resolved X-ray scattering at a free-electron laser.

We describe a method to measure ultrafast protein structural changes using time-resolved wide-angle X-ray scattering at an X-ray free-electron laser. We demonstrated this approach using multiphoton excitation of the Blastochloris viridis photosynthetic reaction center, observing an ultrafast global conformational change that arises within picoseconds and precedes the propagation of heat through the protein. This provides direct structural evidence for a 'protein quake': the hypothesis that proteins rapidly dissipate energy through quake-like structural motions.

In vivo protein crystallization opens new routes in structural biology.

Protein crystallization in cells has been observed several times in nature. However, owing to their small size these crystals have not yet been used for X-ray crystallographic analysis. We prepared nano-sized in vivo-grown crystals of Trypanosoma brucei enzymes and applied the emerging method of free-electron laser-based serial femtosecond crystallography to record interpretable diffraction data. This combined approach will open new opportunities in structural systems biology.

Conformation sequence recovery of a non-periodic object from a diffraction-before-destruction experiment.

Knowledge of the sequence of different conformational states of a protein molecule is key to better understanding its biological function. A diffraction pattern from a single conformational state can be captured with an ultrafast X-ray Free-Electron Laser (XFEL) before the target is completely annihilated by the radiation. In this paper, we report the first experimental demonstration of conformation sequence recovery using diffraction patterns from randomly ordered conformations of a non-periodic object using the dimensional reduction technique Isomap and coherent diffraction imaging.

Explosion dynamics of sucrose nanospheres monitored by time of flight spectrometry and coherent diffractive imaging at the split-and-delay beam line of the FLASH soft X-ray laser.

We use a Mach-Zehnder type autocorrelator to split and delay XUV pulses from the FLASH soft X-ray laser for triggering and subsequently probing the explosion of aerosolised sugar balls. FLASH was running at 182 eV photon energy with pulses of 70 fs duration. The delay between the pump-probe pulses was varied between zero and 5 ps, and the pulses were focused to reach peak intensities above 10¹6W/cm² with an off-axis parabola. The direct pulse triggered the explosion of single aerosolised sucrose nano-particles, while the delayed pulse probed the exploding structure. The ejected ions were measured by ion time of flight spectrometry, and the particle sizes were measured by coherent diffractive imaging. The results show that sucrose particles of 560-1000 nm diameter retain their size for about 500 fs following the first exposure. Significant sample expansion happens between 500 fs and 1 ps. We present simulations to support these observations.

Copper-zinc cross-modulation in prion protein binding.

In this paper we report a systematic XAS study of a set of samples in which Cu(II) was progressively added to complexes in which Zn(II) was bound to the tetra-octarepeat portion of the prion protein. This work extends previous EPR and XAS analysis in which, in contrast, the effect of adding Zn(II) to Cu(II)-tetra-octarepeat complexes was investigated. Detailed structural analysis of the XAS spectra taken at both the Cu and Zn K-edge when the two metals are present at different relative concentrations revealed that Zn(II) and Cu(II) ions compete for binding to the tetra-octarepeat peptide by cross-regulating their relative binding modes. We show that the specific metal-peptide coordination mode depends not only, as expected, on the relative metal concentrations, but also on whether Zn(II) or Cu(II) was first bound to the peptide. In particular, it seems that the Zn(II) binding mode in the absence of Cu(II) is able to promote the formation of small peptide clusters in which triplets of tetra-octarepeats are bridged by pairs of Zn ions. When Cu(II) is added, it starts competing with Zn(II) for binding, disrupting the existing peptide cluster arrangement, despite the fact that Cu(II) is unable to completely displace Zn(II). These results may have a bearing on our understanding of peptide-aggregation processes and, with the delicate cross-regulation balancing we have revealed, seem to suggest the existence of an interesting, finely tuned interplay among metal ions affecting protein binding, capable of providing a mechanism for regulation of metal concentration in cells.

Probing the Dynamic Landscape: From Static to Time-Resolved X-Ray Absorption Spectroscopy to Investigate Copper Redox Chemistry in Neurodegenerative Disorders.

Copper (Cu), with its ability to exist in various oxidation states, notably Cu(I) and Cu(II), plays a crucial role in diverse biological redox reactions. This includes its involvement in pathways associated with oxidative stress in neurodegenerative disorders such as Alzheimer's disease, Parkinson's disease, and Transmissible Spongiform Encephalopathies. This paper offers an overview of X-ray Absorption Spectroscopy (XAS) studies designed to elucidate the interactions between Cu ions and proteins or peptides associated with these neurodegenerative diseases. The emphasis lies on XAS specificity, revealing the local coordination environment, and on its sensitivity to Cu oxidation states. Furthermore, the paper focuses on XAS applications targeting the characterization of intermediate reaction states and explores the opportunities arising from recent advancements in time-resolved XAS at ultrabright synchrotron and Free Electron Laser radiation sources.

Crystallographic data processing for free-electron laser sources.

A processing pipeline for diffraction data acquired using the `serial crystallography' methodology with a free-electron laser source is described with reference to the crystallographic analysis suite CrystFEL and the pre-processing program Cheetah. A detailed analysis of the nature and impact of indexing ambiguities is presented. Simulations of the Monte Carlo integration scheme, which accounts for the partially recorded nature of the diffraction intensities, are presented and show that the integration of partial reflections could be made to converge more quickly if the bandwidth of the X-rays were to be increased by a small amount or if a slight convergence angle were introduced into the incident beam.

The INFN-LNF present and future accelerator-based light facilities.

The INFN-Frascati National Laboratory (LNF) is nowadays running a 0.51 GeV electron-positron collider, DA Φ NE, that also represents the synchrotron radiation source of the beamlines of the DA Φ NE-Light facility. Not being DA Φ NE dedicated to synchrotron radiations activities, the DA Φ NE-Light facility can use it mainly in parasitic mode. Particle accelerators and high energy physics (HEP) have been and are the main core of the LNF research activity, but like other HEP international laboratories also LNF is now moving in the direction of developing a dedicated free electron laser (FEL) user facility, EuPRAXIA@SPARC_Lab, based on plasma acceleration. This new facility in the framework of the EuPRAXIA (European Plasma Research Accelerator with eXcellence in Applications) EU project should produce FEL radiation beams for a wide range of applications using a smaller accelerator compared to actual radio frequency-based accelerator sources dimensions.

Revisiting the pro-oxidant activity of copper: interplay of ascorbate, cysteine, and glutathione.

Copper (Cu) is essential for most organisms, but it can be poisonous in excess, through mechanisms such as protein aggregation, trans-metallation, and oxidative stress. The latter could implicate the formation of potentially harmful reactive oxygen species (O2•-, H2O2, and HO•) via the redox cycling between Cu(II)/Cu(I) states in the presence of dioxygen and physiological reducing agents such as ascorbate (AscH), cysteine (Cys), and the tripeptide glutathione (GSH). Although the reactivity of Cu with these reductants has been previously investigated, the reactions taking place in a more physiologically relevant mixture of these biomolecules are not known. Hence, we report here on the reactivity of Cu with binary and ternary mixtures of AscH, Cys, and GSH. By measuring AscH and thiol oxidation, as well as HO• formation, we show that Cu reacts preferentially with GSH and Cys, halting AscH oxidation and also HO• release. This could be explained by the formation of Cu-thiolate clusters with both GSH and, as we first demonstrate here, Cys. Moreover, we observed a remarkable acceleration of Cu-catalyzed GSH oxidation in the presence of Cys. We provide evidence that both thiol-disulfide exchange and the generated H2O2 contribute to this effect. Based on these findings, we speculate that Cu-induced oxidative stress may be mainly driven by GSH depletion and/or protein disulfide formation rather than by HO• and envision a synergistic effect of Cys on Cu toxicity.

Femtosecond free-electron laser x-ray diffraction data sets for algorithm development.

We describe femtosecond X-ray diffraction data sets of viruses and nanoparticles collected at the Linac Coherent Light Source. The data establish the first large benchmark data sets for coherent diffraction methods freely available to the public, to bolster the development of algorithms that are essential for developing this novel approach as a useful imaging technique. Applications are 2D reconstructions, orientation classification and finally 3D imaging by assembling 2D patterns into a 3D diffraction volume.

Is styrene competitive for dopamine receptor binding?

The potential role of styrene oxide in altering the dopaminergic pathway in the ear is investigated by means of molecular docking and molecular dynamics simulations. We estimate the binding affinity of both styrene oxide and dopamine to the dopaminergic receptor DrD2 by computing the free-energy difference, ∆G, between the configuration where the ligand is bound to the receptor and the situation in which it is "infinitely" far away from it. The results show that the styrene oxide has a somewhat lower affinity for binding with respect to dopamine, which, however, may not be enough to prevent exogenous high concentration styrene oxide to compete with endogenous dopamine for DrD2 binding.

Zinc modulates copper coordination mode in prion protein octa-repeat subdomains.

In this work we present and analyse XAS measurements carried out on various portions of Prion-protein tetra-octa-repeat peptides in complexes with Cu(II) ions, both in the presence and in the absence of Zn(II). Because of the ability of the XAS technique to provide detailed local structural information, we are able to demonstrate that Zn acts by directly interacting with the peptide, in this way competing with Cu for binding with histidine. This finding suggests that metal binding competition can be important in the more general context of metal homeostasis.

Zn-Induced Interactions Between SARS-CoV-2 orf7a and BST2/Tetherin.

We present in this work a first X-ray Absorption Spectroscopy study of the interactions of Zn with human BST2/tetherin and SARS-CoV-2 orf7a proteins as well as with some of their complexes. The analysis of the XANES region of the measured spectra shows that Zn binds to BST2, as well as to orf7a, thus resulting in the formation of BST2-orf7a complexes. This structural information confirms the the conjecture, recently put forward by some of the present Authors, according to which the accessory orf7a (and possibly also orf8) viral protein are capable of interfering with the BST2 antiviral activity. Our explanation for this behavior is that, when BST2 gets in contact with Zn bound to the orf7a Cys15 ligand, it has the ability of displacing the metal owing to the creation of a new disulfide bridge across the two proteins. The formation of this BST2-orf7a complex destabilizes BST2 dimerization, thus impairing the antiviral activity of the latter.

An XAS study of the sulfur environment in human neuromelanin and its synthetic analogs.

Neuromelanin is a complex molecule accumulating in the catecholaminergic neurons that undergo a degenerative process in Parkinson's disease. It has been shown to play either a protective or a toxic role depending on whether it is present in the intraneuronal or extraneuronal milieu. Understanding its structure and synthesis mechanisms is mandatory to clarify the reason for this remarkable dual behavior. In the present study, X-ray absorption spectroscopy is employed to investigate the sulfur binding mode in natural human neuromelanin, synthetic neuromelanins, and in certain structurally known model compounds, namely cysteine and decarboxytrichochrome C. Based on comparative fits of human and synthetic neuromelanin spectra in terms of those of model compounds, the occurrence of both cysteine- and trichochrome-like sulfur coordination modes is recognized, and the relative abundance of these two types of structural arrangement is determined. Data on the amount of cysteine- and trichochrome-like sulfur measured in this way indicate that among the synthetic neuromelanins those produced by enzymatic oxidation are the most similar ones to natural neuromelanin. The interest of the method described here lies in the fact that it allows the identification of different sulfur coordination environments in a physically nondestructive way.