Calculation of the relative metastabilities of proteins using the CHNOSZ software package

<p>Abstract</p> <p>Background</p> <p>Proteins of various compositions are required by organisms inhabiting different environments. The energetic demands for protein formation are a function of the compositions of proteins as well as geochemical variables including temperature, pressure, oxygen fugacity and pH. The purpose of this study was to explore the dependence of metastable equilibrium states of protein systems on changes in the geochemical variables.</p> <p>Results</p> <p>A software package called CHNOSZ implementing the revised Helgeson-Kirkham-Flowers (HKF) equations of state and group additivity for ionized unfolded aqueous proteins was developed. The program can be used to calculate standard molal Gibbs energies and other thermodynamic properties of reactions and to make chemical speciation and predominance diagrams that represent the metastable equilibrium distributions of proteins. The approach takes account of the chemical affinities of reactions in open systems characterized by the chemical potentials of basis species. The thermodynamic database included with the package permits application of the software to mineral and other inorganic systems as well as systems of proteins or other biomolecules.</p> <p>Conclusion</p> <p>Metastable equilibrium activity diagrams were generated for model cell-surface proteins from archaea and bacteria adapted to growth in environments that differ in temperature and chemical conditions. The predicted metastable equilibrium distributions of the proteins can be compared with the optimal growth temperatures of the organisms and with geochemical variables. The results suggest that a thermodynamic assessment of protein metastability may be useful for integrating bio- and geochemical observations.</p

Speciation and fate of trace metals in estuarine sediments under reduced and oxidized conditions, Seaplane Lagoon, Alameda Naval Air Station (USA)

We have identified important chemical reactions that control the fate of metal-contaminated estuarine sediments if they are left undisturbed (in situ) or if they are dredged. We combined information on the molecular bonding of metals in solids from X-ray absorption spectroscopy (XAS) with thermodynamic and kinetic driving forces obtained from dissolved metal concentrations to deduce the dominant reactions under reduced and oxidized conditions. We evaluated the in situ geochemistry of metals (cadmium, chromium, iron, lead, manganese and zinc) as a function of sediment depth (to 100 cm) from a 60 year record of contamination at the Alameda Naval Air Station, California. Results from XAS and thermodynamic modeling of porewaters show that cadmium and most of the zinc form stable sulfide phases, and that lead and chromium are associated with stable carbonate, phosphate, phyllosilicate, or oxide minerals. Therefore, there is minimal risk associated with the release of these trace metals from the deeper sediments contaminated prior to the Clean Water Act (1975) as long as reducing conditions are maintained. Increased concentrations of dissolved metals with depth were indicative of the formation of metal HS(- )complexes. The sediments also contain zinc, chromium, and manganese associated with detrital iron-rich phyllosilicates and/or oxides. These phases are recalcitrant at near-neutral pH and do not undergo reductive dissolution within the 60 year depositional history of sediments at this site. The fate of these metals during dredging was evaluated by comparing in situ geochemistry with that of sediments oxidized by seawater in laboratory experiments. Cadmium and zinc pose the greatest hazard from dredging because their sulfides were highly reactive in seawater. However, their dissolved concentrations under oxic conditions were limited eventually by sorption to or co-precipitation with an iron (oxy)hydroxide. About 50% of the reacted CdS and 80% of the reacted ZnS were bonded to an oxide-substrate at the end of the 90-day oxidation experiment. Lead and chromium pose a minimal hazard from dredging because they are bonded to relatively insoluble carbonate, phosphate, phyllosilicate, or oxide minerals that are stable in seawater. These results point out the specific chemical behavior of individual metals in estuarine sediments, and the need for direct confirmation of metal speciation in order to constrain predictive models that realistically assess the fate of metals in urban harbors and coastal sediments

Modeling acid-gas generation from boiling chloride brines

<p>Abstract</p> <p>Background</p> <p>This study investigates the generation of HCl and other acid gases from boiling calcium chloride dominated waters at atmospheric pressure, primarily using numerical modeling. The main focus of this investigation relates to the long-term geologic disposal of nuclear waste at Yucca Mountain, Nevada, where pore waters around waste-emplacement tunnels are expected to undergo boiling and evaporative concentration as a result of the heat released by spent nuclear fuel. Processes that are modeled include boiling of highly concentrated solutions, gas transport, and gas condensation accompanied by the dissociation of acid gases, causing low-pH condensate.</p> <p>Results</p> <p>Simple calculations are first carried out to evaluate condensate pH as a function of HCl gas fugacity and condensed water fraction for a vapor equilibrated with saturated calcium chloride brine at 50-150°C and 1 bar. The distillation of a calcium-chloride-dominated brine is then simulated with a reactive transport model using a brine composition representative of partially evaporated calcium-rich pore waters at Yucca Mountain. Results show a significant increase in boiling temperature from evaporative concentration, as well as low pH in condensates, particularly for dynamic systems where partial condensation takes place, which result in enrichment of HCl in condensates. These results are in qualitative agreement with experimental data from other studies.</p> <p>Conclusion</p> <p>The combination of reactive transport with multicomponent brine chemistry to study evaporation, boiling, and the potential for acid gas generation at the proposed Yucca Mountain repository is seen as an improvement relative to previously applied simpler batch evaporation models. This approach allows the evaluation of thermal, hydrological, and chemical (THC) processes in a coupled manner, and modeling of settings much more relevant to actual field conditions than the distillation experiment considered. The actual and modeled distillation experiments do not represent expected conditions in an emplacement drift, but nevertheless illustrate the potential for acid-gas generation at moderate temperatures (<150°C).</p

A metastable equilibrium model for the relative abundances of microbial phyla in a hot spring

Many studies link the compositions of microbial communities to their environments, but the energetics of organism-specific biomass synthesis as a function of geochemical variables has rarely been assessed. We describe a thermodynamic model that integrates geochemical and metagenomic data for biofilms sampled at five sites along a thermal and chemical gradient in the outflow channel of the hot spring known as ‘‘Bison Pool’’ in Yellowstone National Park. The relative abundances of major phyla in individual communities sampled along the outflow channel are modeled by computing metastable equilibrium among model proteins with amino acid compositions derived from metagenomic sequences. Geochemical conditions are represented by temperature and activities of basis species, including pH and oxidation-reduction potential quantified as the activity of dissolved hydrogen. By adjusting the activity of hydrogen, the model can be tuned to closely approximate the relative abundances of the phyla observed in the community profiles generated from BLAST assignments. The findings reveal an inverse relationship between the energy demand to form the proteins at equal thermodynamic activities and the abundance of phyla in the community.Although the metabolisms used by many members of these communities are driven by chemical disequilibria, the results support the possibility that higher-level patterns of chemotrophic microbial ecosystems are shaped by metastable equilibrium states that depend on both the composition of biomass and the environmental conditions

Evaporative evolution of a Na–Cl–NO(3)–K–Ca–SO(4)–Mg–Si brine at 95°C: Experiments and modeling relevant to Yucca Mountain, Nevada

A synthetic Topopah Spring Tuff water representative of one type of pore water at Yucca Mountain, NV was evaporated at 95°C in a series of experiments to determine the geochemical controls for brines that may form on, and possibly impact upon the long-term integrity of waste containers and drip shields at the designated high-level, nuclear-waste repository. Solution chemistry, condensed vapor chemistry, and precipitate mineralogy were used to identify important chemical divides and to validate geochemical calculations of evaporating water chemistry using a high temperature Pitzer thermodynamic database. The water evolved toward a complex "sulfate type" brine that contained about 45 mol % Na, 40 mol % Cl, 9 mol % NO(3), 5 mol % K, and less than 1 mol % each of SO(4), Ca, Mg, ∑CO(2)(aq), F, and Si. All measured ions in the condensed vapor phase were below detection limits. The mineral precipitates identified were halite, anhydrite, bassanite, niter, and nitratine. Trends in the solution composition and identification of CaSO(4 )solids suggest that fluorite, carbonate, sulfate, and magnesium-silicate precipitation control the aqueous solution composition of sulfate type waters by removing fluoride, calcium, and magnesium during the early stages of evaporation. In most cases, the high temperature Pitzer database, used by EQ3/6 geochemical code, sufficiently predicts water composition and mineral precipitation during evaporation. Predicted solution compositions are generally within a factor of 2 of the experimental values. The model predicts that sepiolite, bassanite, amorphous silica, calcite, halite, and brucite are the solubility controlling mineral phases

Calculation of the relative metastabilities of proteins in subcellular compartments of Saccharomyces cerevisiae

[abridged] Background: The distribution of chemical species in an open system at metastable equilibrium can be expressed as a function of environmental variables which can include temperature, oxidation-reduction potential and others. Calculations of metastable equilibrium for various model systems were used to characterize chemical transformations among proteins and groups of proteins found in different compartments of yeast cells. Results: With increasing oxygen fugacity, the relative metastability fields of model proteins for major subcellular compartments go as mitochondrion, endoplasmic reticulum, cytoplasm, nucleus. In a metastable equilibrium setting at relatively high oxygen fugacity, proteins making up actin are predominant, but those constituting the microtubule occur with a low chemical activity. A reaction sequence involving the microtubule and spindle pole proteins was predicted by combining the known intercompartmental interactions with a hypothetical program of oxygen fugacity changes in the local environment. In further calculations, the most-abundant proteins within compartments generally occur in relative abundances that only weakly correspond to a metastable equilibrium distribution. However, physiological populations of proteins that form complexes often show an overall positive or negative correlation with the relative abundances of proteins in metastable assemblages. Conclusions: This study explored the outlines of a thermodynamic description of chemical transformations among interacting proteins in yeast cells. The results suggest that these methods can be used to measure the degree of departure of a natural biochemical process or population from a local minimum in Gibbs energy.Comment: 32 pages, 7 figures; supporting information is available at http://www.chnosz.net/yeas

Transport of Pb and Zn by carboxylate complexes in basinal ore fluids and related petroleum-field brines at 100°C: the influence of pH and oxygen fugacity

It is well established through field observations, experiments, and chemical models that oxidation (redox) state and pH exert a strong influence on the speciation of dissolved components and the solubility of minerals in hydrothermal fluids. log [Image: see text] –pH diagrams were used to depict the influence of oxygen fugacity and pH on monocarboxylate- and dicarboxylate-transport of Pb and Zn in low-temperature (100°C) hydrothermal ore fluids that are related to diagenetic processes in deep sedimentary basins, and allow a first-order comparison of Pb and Zn transport among proposed model fluids for Mississippi Valley-type (MVT) and red-bed related base metal (RBRBM) deposits in terms of their approximate pH and [Image: see text] conditions. To construct these diagrams, total Pb and Zn concentrations and Pb and Zn speciation were calculated as a function of log [Image: see text] and pH for a composite ore-brine with concentrations of major elements, total sulfur, and total carbonate that approximate the composition of MVT and RBRBM model ore fluids and modern basinal brines. In addition to acetate and malonate complexation, complexes involving the ligands Cl(-), HS(-), H(2)S, and OH(- )were included in the model of calculated total metal concentration and metal speciation. Also, in the model, Zn and Pb are competing with the common-rock forming metals Ca, Mg, Na, Fe, and Al for the same ligands. Calculated total Pb concentration and calculated total Zn concentration are constrained by galena and sphalerite solubility, respectively. Isopleths, in log [Image: see text] –pH space, of the concentration of Pb and concentration of Zn in carboxylate (acetate + malonate) complexes illustrate that the oxidized model fluids of T. H. Giordano (in Organic Acids in Geological Processes, ed. E. D. Pittman and M. D. Lewan, Springer-Verlag, New York, 1994, pp. 319–354) and G. M. Anderson (Econ. Geol., 1975, 70, 937–942) are capable of transporting sufficient amounts of Pb (up to 10 ppm) and Zn (up to 100 ppm) in the form of carboxylate complexes to form economic deposits of these metals. On the other hand, the reduced ore fluid models of D. A. Sverjensky (Econ. Geol., 1984, 79, 23–37) and T. H. Giordano and H. L. Barnes (Econ. Geol., 1981, 76, 2200–2211) can at best transport amounts of Pb and Zn, as carboxylate complexes, that are many orders of magnitude below the 1 to 10 ppm minimum required to form economic deposits. Lead and zinc speciation (mol% of total Pb or Zn) in the model ore fluid was calculated at specific log [Image: see text] –pH conditions along the 100, 0.01, and 0.001 ppm total Pb and total Zn isopleths. Along the 100 ppm isopleth conditions are oxidized (∑SO(4 )>> ∑H(2)S) with Pb and Zn predominantly in the form of chloride complexes under acid to mildly alkaline conditions (pH from 3 to approximately 7.5), while hydroxide complexes dominate Pb and Zn speciation under more alkaline conditions. Sulfide complexes are insignificant under these oxidized conditions. For more reduced conditions along the 0.01 and 0.001 ppm isopleths chloride complexes dominate Pb and Zn speciation in the SO(4)(2- )field and near the SO(4)(2-)-reduced sulfur boundary from pH = 4 to approximately 7.5, while hydroxide complexes dominate Pb and Zn speciation under alkaline conditions above pH = 7.5 in the SO(4)(2- )field. In the most reduced fluids (∑H(2)S >> ∑SO(4)) along the 0.01 and 0.001 isopleths, sulfide complexes account for almost 100% of the Pb and Zn in the model fluid. Acetate (monocarboxylate) complexation is significant only under conditions of chloride and hydroxide complex dominance and its effect is maximized in the pH range 5 to 7, where it complexes 2 to 2.6% of the total Pb and 1 to 1.25% of the total Zn. Malonate (dicarboxylate) complexes are insignificant along all isopleths. The speciation results from this study show that deep formation waters characterized by temperatures near 100°C, high oxidation states and ∑H(2)S < 0.03 mg L(-1 )([Image: see text] < 10(-6)), high chlorinities (~ 100000 mg L(-1)), and high but reasonable concentrations of carboxylate anions can mobilize up to 3% of the total Pb and up to 1.3% of the total Zn as carboxylate complexes. Furthermore, these percentages, under the most favorable conditions, correspond to approximately 1 to 100 ppm of these metals in solution; concentrations that are adequate to form economic deposits of these metals. However, the field evidence suggests that all of these optimum conditions for carboxylate complexation are rarely met at the same time. A comparison of the composite ore fluid compositions from this study and modern brine data shows that the ore brines, corresponding to log [Image: see text] –pH conditions based on the Anderson (1975) and Giordano (1994) model fluids, are similar in many respects to modern, high trace-metal petroleum-field brines. The principal differences between modern high trace-metal brines and the composite ore fluids of Anderson (1975) and Giordano (1994) relate to their carboxylate anion content. The reported concentrations of monocarboxylate anions (∑monocbx) and dicarboxylate anions (Edicbx) in high trace-metal petroleum-field brines (< 1 to 300 mg L(-1 )and < 1 mg L(-1), respectively) are significantly lower than the concentrations assumed in the modelled brines of this study (∑monocbx = 7 700 mg L(-1 )and ∑dicbx = 300 mg L(-1)). There are also major differences in the corresponding total chloride to carboxylate ratio (∑m(Cl)/∑m(cbx)) and monocarboxylate to dicarboxylate ratio (∑m(monocbx)/∑m(dicbx)). Modern high trace-metal brines have much higher ∑m(Cl)/∑m(cbx )values and, therefore, the contribution of carboxylate complexes to the total Pb and Zn content in these modern brines is likely to be significantly less than the 1 to 3 percent for the composite ore fluids of Anderson (1975) and Giordano (1994). The composite ore-brine based on the Giordano and Barnes (1981) MVT ore fluid is comparable to the high salinity (> 170 000 mg L(-1 )TDS) subset of modern brines characterized by low trace-metal content and high total reduced sulfur (∑H(2)S). A comparison of the Sverjensky (1984) composite ore-brine with modern petroleum-field brines in terms of ∑H(2)S and Zn content, reveals that this ore fluid corresponds to a "border-type" brine, between modern high trace-metal brines and those with low trace-metal content and high ∑H(2)S. A brine of this type is characterized by values of ∑H(2)S, ∑Zn, and/or ∑Pb within or near the 1 to 10 mg L(-1 )range. Based on brine-composition data from numerous references cited in this paper, border-type brines do exist but are rare. The model results and field evidence presented in this study are consistent with other chemical simulation studies of carboxylate complexation in modern petroleum-field brines. Thus, it appears that carboxylate complexation plays a minor, if not insignificant, role as a transport mechanism for Pb and Zn in high salinity Na–Cl and Na–Ca–Cl basinal brines and related ore fluids

Activity variations attending tungsten skarn formation, Pine Creek, California

An integrated geochemical analysis of the well-exposed Pine Creek, California tungsten skarn deposit has been undertaken to evaluate changes in chemical gradients across various lithologies. Thermodynamic calculations using available experimental and thermodynamic data allow limits to be assigned to the activities of important chemical components in the metasomatic environment. Quantifiable changes in “non-volatile” component activites (CaO, MgO, Al 2 O 3 , Fe 2 O 3 , WO 3 ) and in fugacities (O 2 , F 2 ) have been traced across the system. The activities of Al 2 O 3 , Fe 2 O 3 and WO 3 generally increase from the marble (<10 2 , <10 −6 , <10 −5 respectively), through the outer skarn zone and into the massive garnet skarn (10 −1.7±0.3 , 10 −3.4±0.4 , 10 −4.8±0.1 ) While CaO and MgO activities decrease for the same traverse from 10 −5 and 10 −2.1±1 respectively, to <10 −5.7 and <10 −3 . Calculated oxygen fugacities are 10 −23.5+1.0 at T =800 K (527° C), about one log unit below QFM, and more reducing than that required by Mt-Py-Po. The high variance of the garnet-pyroxene-quartz assemblages adds sufficient uncertainty to the calculated activities for individual specimens that only the large-scale trends survive the small-scale scatter. None of the chemical variables emerge as major independent or controlling factors for the mineralogy or phase compositions. Changes in the activity of one component may be offset by compensatory changes in another resulting in an environment that, while different from Pine Creek, could still host scheelite mineralization. Mass balance calculations indicate that the exposed endoskarn cannot have supplied the necessary chemical components to convert the country rock to skarn.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/47341/1/410_2004_Article_BF00381557.pd