Hemeproteins as Targets for Sulfide Species.

Significance: Sulfides are endogenous and ubiquitous signaling species that share the hemeproteins as biochemical targets with O2, nitric oxide, and carbon monoxide. The description of the binding mechanisms is mandatory to anticipate the biochemical relevance of the interaction. Recent Advances: The binding of sulfide to ferric hemeproteins has been described in more than 40 systems, including native proteins, mutants, and model systems. Mechanisms of sulfide binding to ferric hemeproteins have been examined by a combination of kinetic and computational experiments. The distal control of the association process, dissected into the migration of the ligand to the active site and the binding event, reveals that neutral hydrogen sulfide (H2S) reaches the active site and is the predominant binding ligand, while the HS- is excluded by the protein matrix. Experiments with model compounds, devoid of a protein scaffold, reveal that both H2S and HS- can bind the ferric heme if accessing the site. A critical role of the proximal ligand in the prevention of the metal-centered reduction has been experimentally assessed. For metmyoglobin and methemoglobin, the coordination of sulfide leads to noncanonical functions: sulfide storage and its oxidative detoxification have been evidenced under physiological and excess sulfide concentrations, respectively. Critical Issues: The bound species is suggested to predominate in the monoprotonated form, although spectroscopic evidence is pending. Future Directions: A description of the role of hemeproteins as biochemical targets for inorganic sulfide requires understanding the reactivity of bound sulfide, for example: the metal-centered reduction, the reaction with excess sulfide, oxidants, or other gasotransmitters, among other biomolecules.

Conformational Flexibility Drives Cold Adaptation in Pseudoalteromonas haloplanktis TAC125 Globins.

Significance: Temperature is one of the most important drivers in shaping protein adaptations. Many biochemical and physiological processes are influenced by temperature. Proteins and enzymes from organisms living at low temperature are less stable in comparison to high-temperature adapted proteins. The lower stability is generally due to greater conformational flexibility. Recent Advances: Adaptive changes in the structure of cold-adapted proteins may occur at subunit interfaces, distant from the active site, thus producing energy changes associated with conformational transitions transmitted to the active site by allosteric modulation, valid also for monomeric proteins in which tertiary structural changes may play an essential role. Critical Issues: Despite efforts, the current experimental and computational methods still fail to produce general principles on protein evolution, since many changes are protein and species dependent. Environmental constraints or other biological cellular signals may override the ancestral information included in the structure of the protein, thus introducing inaccuracy in estimates and predictions on the evolutionary adaptations of proteins in response to cold adaptation. Future Directions: In this review, we describe the studies and approaches used to investigate stability and flexibility in the cold-adapted globins of the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125. In fact, future research directions will be prescient on more detailed investigation of cold-adapted proteins and the role of fluctuations between different conformational states.

Lessons learned about steered molecular dynamics simulations and free energy calculations.

The calculation of free energy profiles is central in understanding differential enzymatic activity, for instance, involving chemical reactions that require QM-MM tools, ligand migration, and conformational rearrangements that can be modeled using classical potentials. The use of steered molecular dynamics (sMD) together with the Jarzynski equality is a popular approach in calculating free energy profiles. Here, we first briefly review the application of the Jarzynski equality to sMD simulations, then revisit the so-called stiff-spring approximation and the consequent expectation of Gaussian work distributions and, finally, reiterate the practical utility of the second-order cumulant expansion, as it coincides with the parametric maximum-likelihood estimator in this scenario. We illustrate this procedure using simulations of CO, both in aqueous solution and in a carbon nanotube as a model system for biologically relevant nanoheterogeneous environments. We conclude the use of the second-order cumulant expansion permits the use of faster pulling velocities in sMD simulations, without introducing bias due to large dispersion in the non-equilibrium work distribution.

On the accurate estimation of free energies using the jarzynski equality.

The Jarzynski equality is one of the most widely celebrated and scrutinized nonequilibrium work theorems, relating free energy to the external work performed in nonequilibrium transitions. In practice, the required ensemble average of the Boltzmann weights of infinite nonequilibrium transitions is estimated as a finite sample average, resulting in the so-called Jarzynski estimator, ΔF^J . Alternatively, the second-order approximation of the Jarzynski equality, though seldom invoked, is exact for Gaussian distributions and gives rise to the Fluctuation-Dissipation estimator ΔF^FD . Here we derive the parametric maximum-likelihood estimator (MLE) of the free energy ΔF^ML considering unidirectional work distributions belonging to Gaussian or Gamma families, and compare this estimator to ΔF^J . We further consider bidirectional work distributions belonging to the same families, and compare the corresponding bidirectional ΔF^ML∗ to the Bennett acceptance ratio ( ΔF^BAR ) estimator. We show that, for Gaussian unidirectional work distributions, ΔF^FD is in fact the parametric MLE of the free energy, and as such, the most efficient estimator for this statistical family. We observe that ΔF^ML and ΔF^ML∗ perform better than ΔF^J and ΔF^BAR , for unidirectional and bidirectional distributions, respectively. These results illustrate that the characterization of the underlying work distribution permits an optimal use of the Jarzynski equality. © 2018 Wiley Periodicals, Inc.

Sulfide binding properties of truncated hemoglobins.

The truncated hemoglobins from Bacillus subtilis (Bs-trHb) and Thermobifida fusca (Tf-trHb) have been shown to form high-affinity complexes with hydrogen sulfide in their ferric state. The recombinant proteins, as extracted from Escherichia coli cells after overexpression, are indeed partially saturated with sulfide, and even highly purified samples still contain a small but significant amount of iron-bound sulfide. Thus, a complete thermodynamic and kinetic study has been undertaken by means of equilibrium and kinetic displacement experiments to assess the relevant sulfide binding parameters. The body of experimental data indicates that both proteins possess a high affinity for hydrogen sulfide (K = 5.0 x 10(6) and 2.8 x 10(6) M(-1) for Bs-trHb and Tf-trHb, respectively, at pH 7.0), though lower with respect to that reported previously for the sulfide avid Lucina pectinata I hemoglobins (2.9 x 10(8) M(-1)). From the kinetic point of view, the overall high affinity resides in the slow rate of sulfide release, attributed to hydrogen bonding stabilization of the bound ligand by distal residue WG8. A set of point mutants in which these residues have been replaced with Phe indicates that the WG8 residue represents the major kinetic barrier to the escape of the bound sulfide species. Accordingly, classical molecular dynamics simulations of SH(-)-bound ferric Tf-trHb show that WG8 plays a key role in the stabilization of coordinated SH(-) whereas the YCD1 and YB10 contributions are negligible. Interestingly, the triple Tf-trHb mutant bearing only Phe residues in the relevant B10, G8, and CD1 positions is endowed with a higher overall affinity for sulfide characterized by a very fast second-order rate constant and 2 order of magnitude faster kinetics of sulfide release with respect to the wild-type protein. Resonance Raman spectroscopy data indicate that the sulfide adducts are typical of a ferric iron low-spin derivative. In analogy with other low-spin ferric sulfide adducts, the strong band at 375 cm(-1) is tentatively assigned to a Fe-S stretching band. The high affinity for hydrogen sulfide is thought to have a possible physiological significance as H(2)S is produced in bacteria at metabolic steps involved in cysteine biosynthesis and hence in thiol redox homeostasis.