Comparative exploration of hydrogen sulfide and water transmembrane free energy surfaces via orthogonal space tempering free energy sampling

Sulfide oxidation by members of the Sulfolobales

The oxidation of sulfur compounds drives the acidification of geothermal waters. At high temperatures (>80°C) and in acidic conditions (pH <6.0), oxidation of sulfide has historically been considered an abiotic process that generates elemental sulfur (S0) that, in turn, is oxidized by thermoacidophiles of the model archaeal order Sulfolobales to generate sulfuric acid (i.e. sulfate and protons). Here, we describe five new aerobic and autotrophic strains of Sulfolobales comprising two species that were isolated from acidic hot springs in Yellowstone National Park (YNP) and that can use sulfide as an electron donor. These strains significantly accelerated the rate and extent of sulfide oxidation to sulfate relative to abiotic controls, concomitant with production of cells. Yields of sulfide-grown cultures were ∼2-fold greater than those of S0-grown cultures, consistent with thermodynamic calculations indicating more available energy in the former condition than the latter. Homologs of sulfide:quinone oxidoreductase (Sqr) were identified in nearly all Sulfolobales genomes from YNP metagenomes as well as those from other reference Sulfolobales, suggesting a widespread ability to accelerate sulfide oxidation. These observations expand the role of Sulfolobales in the oxidative sulfur cycle, the geobiological feedbacks that drive the formation of acidic hot springs, and landscape evolution.

Physiological roles of hydrogen sulfide in mammalian cells, tissues and organs

H2S belongs to the class of molecules known as gasotransmitters, which also includes nitric oxide (NO) and carbon monoxide (CO). Three enzymes are recognized as endogenous sources of H2S in various cells and tissues: cystathionine g-lyase (CSE), cystathionine β-synthase (CBS) and 3-mercaptopyruvate sulfurtransferase (3-MST). The current article reviews the regulation of these enzymes as well as the pathways of their enzymatic and non-enzymatic degradation and elimination. The multiple interactions of H2S with other labile endogenous molecules (e.g. NO) and reactive oxygen species are also outlined. The various biological targets and signaling pathways are discussed, with special reference to H2S and oxidative posttranscriptional modification of proteins, the effect of H2S on channels and intracellular second messenger pathways, the regulation of gene transcription and translation and the regulation of cellular bioenergetics and metabolism. The pharmacological and molecular tools currently available to study H2S physiology are also reviewed, including their utility and limitations. In subsequent sections, the role of H2S in the regulation of various physiological and cellular functions is reviewed. The physiological role of H2S in various cell types and organ systems are overviewed. Finally, the role of H2S in the regulation of various organ functions is discussed as well as the characteristic bell-shaped biphasic effects of H2S. In addition, key pathophysiological aspects, debated areas, and future research and translational areas are identified A wide array of significant roles of H2S in the physiological regulation of all organ functions emerges from this review.

Intersection of Biotic and Abiotic Sulfur Chemistry Supporting Extreme Microbial Life in Hot Acid

Microbial life on Earth exists within wide ranges of temperature, pressure, pH, salinity, radiation, and water activity. Extreme thermoacidophiles, in particular, are microbes found in hot, acidic biotopes laden with heavy metals and reduced inorganic sulfur species. As chemolithoautotrophs, they thrive in the absence of organic carbon, instead using sulfur and metal oxidation to fuel their bioenergetic needs, while incorporating CO2 as a carbon source. Metal oxidation by these microbes takes place extracellularly, mediated by membrane-associated oxidase complexes. In contrast, sulfur oxidation involves extracellular, membrane-associated, and cytoplasmic biotransformations, which intersect with abiotic sulfur chemistry. This novel lifestyle has been examined in the context of early aerobic life on this planet, but it is also interesting when considering the prospects of life, now or previously, on other solar bodies. Here, extreme thermoacidophily (growth at pH below 4.0, temperature above 55 °C), a characteristic of species in the archaeal order Sulfolobales, is considered from the perspective of sulfur chemistry, both biotic and abiotic, as it relates to microbial bioenergetics. Current understanding of the mechanisms involved are reviewed which are further expanded through recent experimental results focused on imparting sulfur oxidation capacity on a natively nonsulfur oxidizing extremely thermoacidophilic archaeon, Sulfolobus acidocaldarius, through metabolic engineering.

Peptide dimerization-dissociation rates from replica exchange molecular dynamics

We show how accurate rates of formation and dissociation of peptide dimers can be calculated using direct transition counting (DTC) from replica-exchange molecular dynamics (REMD) simulations. First, continuous trajectories corresponding to system replicas evolving at different temperatures are used to assign conformational states. Second, we analyze the entire REMD data to calculate the corresponding rates at each temperature directly from the number of transition counts. Finally, we compare the kinetics extracted directly, using the DTC method, with indirect estimations based on trajectory likelihood maximization using short-time propagators and on decay rates of state autocorrelation functions. For systems with relatively low-dimensional intrinsic conformational dynamics, the DTC method is simple to implement and leads to accurate temperature-dependent rates. We apply the DTC rate-extraction method to all-atom REMD simulations of dimerization of amyloid-forming NNQQ tetrapetides in explicit water. In an assessment of the REMD sampling efficiency with respect to standard MD, we find a gain of more than a factor of two at the lowest temperature.

An In Silico study of TiO2 nanoparticles interaction with twenty standard amino acids in aqueous solution

Titanium dioxide (TiO₂) is probably one of the most widely used nanomaterials, and its extensive exposure may result in potentially adverse biological effects. Yet, the underlying mechanisms of interaction involving TiO₂ NPs and macromolecules, e.g., proteins, are still not well understood. Here, we perform all-atom molecular dynamics simulations to investigate the interactions between TiO₂ NPs and the twenty standard amino acids in aqueous solution exploiting a newly developed TiO₂ force field. We found that charged amino acids play a dominant role during the process of binding to the TiO₂ surface, with both basic and acidic residues overwhelmingly preferred over the non-charged counterparts. By calculating the Potential Mean Force, we showed that Arg is prone to direct binding onto the NP surface, while Lys needs to overcome a ~2 kT free energy barrier. On the other hand, acidic residues tend to form “water bridges” between their sidechains and TiO₂ surface, thus displaying an indirect binding. Moreover, the overall preferred positions and configurations of different residues are highly dependent on properties of the first and second solvation water. These molecular insights learned from this work might help with a better understanding of the interactions between biomolecules and nanomaterials.

Molecular Dynamics Simulation Study of Permeation of Molecules through Skin Lipid Bilayer

Stratum Corneum (SC), the outermost layer of skin, is mainly responsible for skin's barrier function. The complex lipid matrix of SC determines these barrier properties. In this study, the lipid matrix is modeled as an equimolar mixture of ceramide (CER), cholesterol (CHOL), and free fatty acid (FFA). The permeation of water, oxygen, ethanol, acetic acid, urea, butanol, benzene, dimethyl sulfoxide (DMSO), toluene, phenol, styrene, and ethylbenzene across this layer is studied using a constrained MD simulations technique. Several long constrained simulations are performed at a skin temperature of 310 K under NPT conditions. The free energy profiles and diffusion coefficients along the bilayer normal have been calculated for each molecule. Permeability coefficients are also calculated and compared with experimental data. The main resistance for the permeation of hydrophilic and hydrophobic permeants has been found to be in the interior of the lipid bilayer and near the lipid-water interface, respectively. The obtained permeability is found to be a few orders of magnitude higher than experimental values for hydrophilic molecules while for hydrophobic molecules more discrepancy was observed. Overall, the qualitative ranking is consistent with the experiments.

Extended Adaptive Biasing Force Algorithm. An On-the-Fly Implementation for Accurate Free-Energy Calculations

Proper use of the adaptive biasing force (ABF) algorithm in free-energy calculations needs certain prerequisites to be met, namely, that the Jacobian for the metric transformation and its first derivative be available and the coarse variables be independent and fully decoupled from any holonomic constraint or geometric restraint, thereby limiting singularly the field of application of the approach. The extended ABF (eABF) algorithm circumvents these intrinsic limitations by applying the time-dependent bias onto a fictitious particle coupled to the coarse variable of interest by means of a stiff spring. However, with the current implementation of eABF in the popular molecular dynamics engine NAMD, a trajectory-based post-treatment is necessary to derive the underlying free-energy change. Usually, such a posthoc analysis leads to a decrease in the reliability of the free-energy estimates due to the inevitable loss of information, as well as to a drop in efficiency, which stems from substantial read-write accesses to file systems. We have developed a user-friendly, on-the-fly code for performing eABF simulations within NAMD. In the present contribution, this code is probed in eight illustrative examples. The performance of the algorithm is compared with traditional ABF, on the one hand, and the original eABF implementation combined with a posthoc analysis, on the other hand. Our results indicate that the on-the-fly eABF algorithm (i) supplies the correct free-energy landscape in those critical cases where the coarse variables at play are coupled to either each other or to geometric restraints or holonomic constraints, (ii) greatly improves the reliability of the free-energy change, compared to the outcome of a posthoc analysis, and (iii) represents a negligible additional computational effort compared to regular ABF. Moreover, in the proposed implementation, guidelines for choosing two parameters of the eABF algorithm, namely the stiffness of the spring and the mass of the fictitious particles, are proposed. The present on-the-fly eABF implementation can be viewed as the second generation of the ABF algorithm, expected to be widely utilized in the theoretical investigation of recognition and association phenomena relevant to physics, chemistry, and biology.

Generalized Ensemble Sampling of Enzyme Reaction Free Energy Pathways

Free energy path sampling plays an essential role in computational understanding of chemical reactions, particularly those occurring in enzymatic environments. Among a variety of molecular dynamics simulation approaches, the generalized ensemble sampling strategy is uniquely attractive for the fact that it not only can enhance the sampling of rare chemical events but also can naturally ensure consistent exploration of environmental degrees of freedom. In this review, we plan to provide a tutorial-like tour on an emerging topic: generalized ensemble sampling of enzyme reaction free energy path. The discussion is largely focused on our own studies, particularly ones based on the metadynamics free energy sampling method and the on-the-path random walk path sampling method. We hope that this minipresentation will provide interested practitioners some meaningful guidance for future algorithm formulation and application study.

Sampling errors in free energy simulations of small molecules in lipid bilayers

Free energy simulations are a powerful tool for evaluating the interactions of molecular solutes with lipid bilayers as mimetics of cellular membranes. However, these simulations are frequently hindered by systematic sampling errors. This review highlights recent progress in computing free energy profiles for inserting molecular solutes into lipid bilayers. Particular emphasis is placed on a systematic analysis of the free energy profiles, identifying the sources of sampling errors that reduce computational efficiency, and highlighting methodological advances that may alleviate sampling deficiencies. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.