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.

Probing protein stability: towards a computational atomistic, reliable, affordable, and improvable model.

We present an improved application of a recently proposed computational method designed to evaluate the change of free energy as a function of the average value of a suitably chosen collective variable in proteins. The method is based on a full atomistic description of the protein and its environment. The goal is to understand how the protein melting temperature changes upon single-point mutations, because the sign of the temperature variation will allow us to discriminate stabilizing vs. destabilizing mutations in protein sequences. In this refined application the method is based on altruistic well-tempered metadynamics, a variant of multiple-walkers metadynamics. The resulting metastatistics is then modulated by the maximal constrained entropy principle. The latter turns out to be especially helpful in free-energy calculations as it is able to alleviate the severe limitations of metadynamics in properly sampling folded and unfolded configurations. In this work we apply the computational strategy outlined above in the case of the bovine pancreatic trypsin inhibitor, a well-studied small protein, which is a reference for computer simulations since decades. We compute the variation of the melting temperature characterizing the folding-unfolding process between the wild-type protein and two of its single-point mutations that are seen to have opposite effect on the free energy changes. The same approach is used for free energy difference calculations between a truncated form of frataxin and a set of five of its variants. Simulation data are compared to in vitro experiments. In all cases the sign of the change of melting temperature is reproduced, under the further approximation of using an empirical effective mean-field to average out protein-solvent interactions.

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.

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.

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.

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.

SARS-CoV-2 Virion Stabilization by Zn Binding.

Zinc plays a crucial role in the process of virion maturation inside the host cell. The accessory Cys-rich proteins expressed in SARS-CoV-2 by genes ORF7a and ORF8 are likely involved in zinc binding and in interactions with cellular antigens activated by extensive disulfide bonds. In this report we provide a proof of concept for the feasibility of a structural study of orf7a and orf8 proteins. A conceivable hypothesis is that lack of cellular zinc, or substitution thereof, might lead to a significant slowing down of viral maturation.

Dealing with Cu reduction in X-ray absorption spectroscopy experiments.

In this paper we prove in the exemplary case of the amyloid-ß peptide in complex with Cu(ii) that at the current low temperatures employed in XAS experiments, the time needed for collecting a good quality XAS spectrum is significantly shorter than the time after which structural damage becomes appreciable. Our method takes advantage of the well-known circumstance that the transition of Cu from the oxidized to the reduced form under ionizing radiation can be quantified by monitoring a characteristic peak in the pre-edge region. We show that there exists a sufficiently large time window in which good XAS spectra can be acquired before the structure around the oxidized Cu(ii) ion reorganizes itself into the reduced Cu(i) "resting" structure. We suggest that similar considerations apply to other cases of biological interest, especially when dealing with macromolecules in complex with transition metal ions.

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.

Cu(II)-Zn(II) Cross-Modulation in Amyloid-Beta Peptide Binding: An X-ray Absorption Spectroscopy Study.

In this work we analyze at a structural level the mechanism by which Cu(II) and Zn(II) ions compete for binding to the Aß peptides that is involved in the etiology of Alzheimer's disease. We collected X-ray absorption spectroscopy data on samples containing Aß with Cu and Zn at different concentration ratios. We show that the order in which metals are added to the peptide solution matters and that, when Zn is added first, it prevents Cu from binding. On the contrary, when Cu is added first, it does not (completely) prevent Zn binding to Aß peptides. Our analysis suggests that Cu and Zn ions are coordinated to different numbers of histidine residues depending on the [ion]:[peptide] concentration ratio.

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.

Computational and experimental studies on ß-sheet breakers targeting Aß1-40 fibrils.

In this work we present and compare the results of extensive molecular dynamics simulations of model systems comprising an Aß1-40 peptide in water in interaction with short peptides (ß-sheet breakers) mimicking the 17-21 region of the Aß1-40 sequence. Various systems differing in the customized ß-sheet breaker structure have been studied. Specifically we have considered three kinds of ß-sheet breakers, namely Ac-LPFFD-NH2 and two variants thereof, one obtained by substituting the acetyl group with the sulfonic amino acid taurine (Tau-LPFFD-NH2) and a second novel one in which the aspartic acid is substituted by an asparagine (Ac-LPFFN-NH2). Thioflavin T fluorescence, circular dichroism, and mass spectrometry experiments have been performed indicating that ß-sheet breakers are able to inhibit in vitro fibril formation and prevent the ß sheet folding of portions of the Aß1-40 peptide. We show that molecular dynamics simulations and far UV circular dichroism provide consistent evidence that the new Ac-LPFFN-NH2 ß-sheet breaker is more effective than the other two in stabilizing the native α-helix structure of Aß1-40. In agreement with these results thioflavin T fluorescence experiments confirm the higher efficiency in inhibiting Aß1-40 aggregation. Furthermore, mass spectrometry data and molecular dynamics simulations consistently identified the 17-21 Aß1-40 portion as the location of the interaction region between peptide and the Ac-LPFFN-NH2 ß-sheet breaker.

Mesoscopic behavior from microscopic Markov dynamics and its application to calcium release channels.

A major challenge in biology is to understand how molecular processes determine phenotypic features. We address this fundamental problem in a class of model systems by developing a general mathematical framework that allows the calculation of mesoscopic properties from the knowledge of microscopic Markovian transition probabilities. We show how exact analytic formulae for the first and second moments of resident time distributions in mesostates can be derived from microscopic resident times and transition probabilities even for systems with a large number of microstates. We apply our formalism to models of the inositol trisphosphate receptor, which plays a key role in generating calcium signals triggering a wide variety of cellular responses. We demonstrate how experimentally accessible quantities, such as opening and closing times and the coefficient of variation of inter-spike intervals, and other, more elaborated, quantities can be analytically calculated from the underlying microscopic Markovian dynamics. A virtue of our approach is that we do not need to follow the detailed time evolution of the whole system, as we derive the relevant properties of its steady state without having to take into account the often extremely complicated transient features. We emphasize that our formulae fully agree with results obtained by stochastic simulations and approaches based on a full determination of the microscopic system's time evolution. We also illustrate how experiments can be devised to discriminate between alternative molecular models of the inositol trisphosphate receptor. The developed approach is applicable to any system described by a Markov process and, owing to the analytic nature of the resulting formulae, provides an easy way to characterize also rare events that are of particular importance to understand the intermittency properties of complex dynamic systems.

Comparison between SAPS II and SAPS 3 in predicting hospital mortality in a cohort of 103 Italian ICUs. Is new always better?

PURPOSE: More recent severity scores should be more reliable than older ones because they account for the improvement in medical care over time. To provide more insight into this issue, we compared the predictive ability of the Simplified Acute Physiology Score (SAPS) II and SAPS 3 (originally developed from data collected in 1991-1992 and 2002, respectively) on a sample of critically ill patients. METHODS: This was a prospective observational study on 3,661 patients from 103 Italian intensive care units. Standardized mortality ratios (SMRs) were calculated. Assessment of calibration across risk classes was performed using the GiViTI calibration belt. Discrimination was evaluated by means of the area under the receiver operating characteristic analysis. RESULTS: Both scores were shown to discriminate fairly. SAPS 3 largely overpredicted mortality, more than SAPS II (SMR 0.63, 95 % CI 0.60-0.66 vs. 0.87, 95 % CI 0.83-0.91). This result was consistent and statistically significant across all risk classes for SAPS 3. SAPS II did not show relevant deviations from ideal calibration in the first two deciles of risk, whereas in higher-risk classes it overpredicted mortality. CONCLUSIONS: Both scores provided unreliable predictions, but unexpectedly the newer SAPS 3 turned out to overpredict mortality more than the older SAPS II.

Zn induced structural aggregation patterns of ß-amyloid peptides by first-principle simulations and XAS measurements.

We show in this paper that in the presence of Zn ions a peculiar structural aggregation pattern of ß-amyloid peptides in which metal ions are sequentially coordinated to either three or four histidines of nearby peptides is favored. To stabilize this configuration a deprotonated imidazole ring from one of the histidines forms a bridge connecting two adjacent Zn ions. Though present in zeolite imidazolate frameworks, remarkably in biological compounds this peculiar Zn-imidazolate-Zn topology is only found in enzymes belonging to the Cu,Zn-superoxide dismutase family in the form of an imidazolate bridging Cu and Zn. The results we present are obtained by combining X-ray absorption spectroscopy experimental data with detailed first-principle molecular dynamics simulations.

Modeling the interplay of glycine protonation and multiple histidine binding of copper in the prion protein octarepeat subdomains.

The octarepeat region of the prion protein can bind Cu(2+) ions up to full occupancy (one ion per octarepeat) at neutral pH. While crystallographic data show that the HGGG octarepeat subdomain is the basic binding unit, multiple histidine coordination at lower Cu occupancy has been reported by X-ray absorption spectroscopy, EPR, and potentiometric experiments. In this paper we investigate, with first principles Car-Parrinello simulations, the first step for the formation of the Cu low-level binding mode, where four histidine side chains are coordinated to the same Cu(2+) ion. This step involves the further binding of a second histidine to an already HGGG domain bonded Cu(2+) ion. The influence of the pH on the ability of Cu to bind two histidine side chains was taken into account by simulating different protonation states of the amide N atoms of the two glycines lying nearest to the first histidine. Multiple histidine coordination is also seen to occur when glycine deprotonation occurs and the presence of the extra histidine stabilizes the Cu-peptide complex. Though the stabilization effect slightly decreases with the number of deprotonated glycines (reaching a minimum when both N atoms of the two nearest glycines are available as Cu ligands), the system is still capable of binding the second histidine in a 4N tetrahedral (though slightly distorted) coordination, whose energy is very near to that of the crystallographic square-planar 3N1O coordination. This result suggests that at low metal concentration the reorganization energy associated with Cu(II)/Cu(I) reduction is small also at pH approximately 7, when glycines are deprotonated.

Variable costs of ICU patients: a multicenter prospective study.

OBJECTIVE: To analyze the costs of treating critically ill patients. DESIGN AND SETTING: Multicenter, observational, prospective, cohort, bottom-up study on variable costs in 51 ICUs. PATIENTS AND PARTICIPANTS: A total of 1,034 patients aged over 14 years who either spent less than 48 h in the ICU or had multiple trauma, major abdominal surgery, ischemic stroke, chronic obstructive pulmonary disease, cardiac failure, isolated head injury, acute lung injury/adult respiratory distress syndrome (ALI/ARDS), nontraumatic intracranial hemorrhage or coronary surgery. INTERVENTIONS: Data recorded for each patient: length of ICU stay, and cost in euros of all diagnostic and therapeutic procedures, drugs and equipment used, and consultations by physicians from other units. To express cost-efficiency we calculated for each diagnostic group the cost per surviving patient (expenditure for all patients/number of surviving patients) and money loss per patient (expenditure for patients who died/total number of patients). MEASUREMENTS AND RESULTS: Median costs for a multiple trauma patient were euro 4076 and for coronary surgery patient euro 380. The variability is largely due to different lengths of ICU stay. Cost per surviving patient was higher for ALI/ARDS, nontraumatic intracranial hemorrhage, multiple trauma, and emergency abdominal surgery. Money loss per patient was higher for ALI/ARDS and lower for multiple trauma. Planned coronary and major abdominal surgery and short-stay patients were treated most cost-efficiently. CONCLUSIONS: Cost of treatment in an ICU varies widely for different types of patients. Strategies are needed to contain the major determinants of high costs and low cost-efficiency.

Costs of the COPD. Differences between intensive care unit and respiratory intermediate care unit.

INTRODUCTION: To assess whether respiratory intermediate care units (RICUs) are cost effective alternatives to intensive care units (ICUs) for patients with exacerbation of chronic obstructive pulmonary disease (COPD). PATIENTS AND METHODS: Multi-centre, prospective, bottom-up cost study performed in 15 ICUs and 6 RICUs. COPD patients staying longer than 48 h were recruited; those coming from other ICUs/RICUs, with immune-deficiency or stroke, were excluded. After the ICU sample was standardised to the RICU distribution of the reason-for-admission and infusion of a vasoactive drug on admission, 60 ICU patients and 65 RICU patients remained, of the original 164 recruited. For each patient, besides clinical data on admission and discharge, daily information about the resources consumed were recorded and analysed in terms of their costs. RESULTS: Total cost per patient was lower in RICUs than in ICUs (754 vs. 1507 Euro; P < 0.0001). In all items, except drugs and nutrition, we found a significant lower cost in RICUs. Dead patients were noticeably different in terms of disease severity between ICUs and RICUs, while surviving ones were not. CONCLUSIONS: Our study suggests that some COPD patients, less severe and with pure respiratory failure, could be successfully and less costly treated in RICUs.

The relationship between labour cost per patient and the size of intensive care units: a multicentre prospective study.

OBJECTIVE: We examined the relationship between major ICU characteristics and labour cost per patient. DESIGN: Four-week prospective data collection, in which the hours spent by each physician and nurse on both in-ICU and extra-ICU activities were collected. SETTING: Eighty Italian adult ICUs. MEASUREMENTS AND RESULTS: The cost of the time actually spent by ICU staff on ICU patients (labour cost) was computed for each participating unit, by applying to the average annual salaries the proportions of in-ICU activity working time for physicians and nurses. Multiple regression analysis was used to identify ICU characteristics that predict labour costs per patient. Labour cost per patient was positively correlated with ICU mortality and patients average length of stay (slopes =0.67, p =0.048 and 0.09, p <0.0001, respectively). Labour cost per patient decreases almost linearly as the number of beds increases up to about eight, and it remains nearly constant above about twelve beds. The number of patients admitted per physician (not per nurse) increases with the number of beds (Spearman correlation coefficient =0.567, p <0.0001). CONCLUSIONS: Our findings suggest that ICUs with less than about 12 beds are not cost-effective.