Atmospheric energy for subsurface life on Mars?

The location and density of biologically useful energy sources on Mars will limit the biomass, spatial distribution, and organism size of any biota. Subsurface Martian organisms could be supplied with a large energy flux from the oxidation of photochemically produced atmospheric H2 and CO diffusing into the regolith. However, surface abundance measurements of these gases demonstrate that no more than a few percent of this available flux is actually being consumed, suggesting that biological activity driven by atmospheric H2 and CO is limited in the top few hundred meters of the subsurface. This is significant because the available but unused energy is extremely large: for organisms at 30-m depth, it is 2,000 times previous estimates of hydrothermal and chemical weathering energy and far exceeds the energy derivable from other atmospheric gases. This also implies that the apparent scarcity of life on Mars is not attributable to lack of energy. Instead, the availability of liquid water may be a more important factor limiting biological activity because the photochemical energy flux can only penetrate to 100- to 1,000-m depth, where most H2O is probably frozen. Because both atmospheric and Viking lander soil data provide little evidence for biological activity, the detection of short-lived trace gases will probably be a better indicator of any extant Martian life.

The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.

Energy metabolism, enzymatic flux capacities, and metabolic flux rates in flying honeybees.

Honeybees rely primarily on the oxidation of hexose sugars to provide the energy required for flight. Measurement of VCO2 (equal to VO2, because VCO2/VO2 = 1.0 during carbohydrate oxidation) during flight allowed estimation of steady-state flux rates through pathways of flight muscle energy metabolism. Comparison of Vmax values for flight muscle hexokinase, phosphofructokinase, citrate synthase, and cytochrome c oxidase with rates of carbon and O2 flux during flight reveal that these enzymes operate closer to Vmax in the flight muscles of flying honeybees than in other muscles previously studied. Possible mechanistic and evolutionary implications of these findings are discussed.

Single-pair fluorescence resonance energy transfer on freely diffusing molecules: Observation of Förster distance dependence and subpopulations

Photon bursts from single diffusing donor-acceptor labeled macromolecules were used to measure intramolecular distances and identify subpopulations of freely diffusing macromolecules in a heterogeneous ensemble. By using DNA as a rigid spacer, a series of constructs with varying intramolecular donor-acceptor spacings were used to measure the mean and distribution width of fluorescence resonance energy transfer (FRET) efficiencies as a function of distance. The mean single-pair FRET efficiencies qualitatively follow the distance dependence predicted by Förster theory. Possible contributions to the widths of the FRET efficiency distributions are discussed, and potential applications in the study of biopolymer conformational dynamics are suggested. The ability to measure intramolecular (and intermolecular) distances for single molecules implies the ability to distinguish and monitor subpopulations of molecules in a mixture with different distances or conformational states. This is demonstrated by monitoring substrate and product subpopulations before and after a restriction endonuclease cleavage reaction. Distance measurements at single-molecule resolution also should facilitate the study of complex reactions such as biopolymer folding. To this end, the denaturation of a DNA hairpin was examined by using single-pair FRET.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Equilibrium distributions of topological states in circular DNA: Interplay of supercoiling and knotting

Two variables define the topological state of closed double-stranded DNA: the knot type, K, and ΔLk, the linking number difference from relaxed DNA. The equilibrium distribution of probabilities of these states, P(ΔLk, K), is related to two conditional distributions: P(ΔLk|K), the distribution of ΔLk for a particular K, and P(K|ΔLk) and also to two simple distributions: P(ΔLk), the distribution of ΔLk irrespective of K, and P(K). We explored the relationships between these distributions. P(ΔLk, K), P(ΔLk), and P(K|ΔLk) were calculated from the simulated distributions of P(ΔLk|K) and of P(K). The calculated distributions agreed with previous experimental and theoretical results and greatly advanced on them. Our major focus was on P(K|ΔLk), the distribution of knot types for a particular value of ΔLk, which had not been evaluated previously. We found that unknotted circular DNA is not the most probable state beyond small values of ΔLk. Highly chiral knotted DNA has a lower free energy because it has less torsional deformation. Surprisingly, even at |ΔLk| > 12, only one or two knot types dominate the P(K|ΔLk) distribution despite the huge number of knots of comparable complexity. A large fraction of the knots found belong to the small family of torus knots. The relationship between supercoiling and knotting in vivo is discussed.

Single-molecule protein folding: Diffusion fluorescence resonance energy transfer studies of the denaturation of chymotrypsin inhibitor 2

We report single-molecule folding studies of a small, single-domain protein, chymotrypsin inhibitor 2 (CI2). CI2 is an excellent model system for protein folding studies and has been extensively studied, both experimentally (at the ensemble level) and theoretically. Conformationally assisted ligation methodology was used to synthesize the proteins and site-specifically label them with donor and acceptor dyes. Folded and denatured subpopulations were observed by fluorescence resonance energy transfer (FRET) measurements on freely diffusing single protein molecules. Properties of these subpopulations were directly monitored as a function of guanidinium chloride concentration. It is shown that new information about different aspects of the protein folding reaction can be extracted from such subpopulation properties. Shifts in the mean transfer efficiencies are discussed, FRET efficiency distributions are translated into potentials, and denaturation curves are directly plotted from the areas of the FRET peaks. Changes in stability caused by mutation also are measured by comparing pseudo wild-type CI2 with a destabilized mutant (K17G). Current limitations and future possibilities and prospects for single-pair FRET protein folding investigations are discussed.

Dynamics and folding of single two-stranded coiled-coil peptides studied by fluorescent energy transfer confocal microscopy

We report single-molecule measurements on the folding and unfolding conformational equilibrium distributions and dynamics of a disulfide crosslinked version of the two-stranded coiled coil from GCN4. The peptide has a fluorescent donor and acceptor at the N termini of its two chains and a Cys disulfide near its C terminus. Thus, folding brings the two N termini of the two chains close together, resulting in an enhancement of fluorescent resonant energy transfer. End-to-end distance distributions have thus been characterized under conditions where the peptide is nearly fully folded (0 M urea), unfolded (7.4 M urea), and in dynamic exchange between folded and unfolded states (3.0 M urea). The distributions have been compared for the peptide freely diffusing in solution and deposited onto aminopropyl silanized glass. As the urea concentration is increased, the mean end-to-end distance shifts to longer distances both in free solution and on the modified surface. The widths of these distributions indicate that the molecules are undergoing millisecond conformational fluctuations. Under all three conditions, these fluctuations gave nonexponential correlations on 1- to 100-ms time scale. A component of the correlation decay that was sensitive to the concentration of urea corresponded to that measured by bulk relaxation kinetics. The trajectories provided effective intramolecular diffusion coefficients as a function of the end-to-end distances for the folded and unfolded states. Single-molecule folding studies provide information concerning the distributions of conformational states in the folded, unfolded, and dynamically interconverting states.

Global air-sea flux of CO2: An estimate based on measurements of sea–air pCO2 difference

Approximately 250,000 measurements made for the pCO2 difference between surface water and the marine atmosphere, ΔpCO2, have been assembled for the global oceans. Observations made in the equatorial Pacific during El Nino events have been excluded from the data set. These observations are mapped on the global 4° × 5° grid for a single virtual calendar year (chosen arbitrarily to be 1990) representing a non-El Nino year. Monthly global distributions of ΔpCO2 have been constructed using an interpolation method based on a lateral advection–diffusion transport equation. The net flux of CO2 across the sea surface has been computed using ΔpCO2 distributions and CO2 gas transfer coefficients across sea surface. The annual net uptake flux of CO2 by the global oceans thus estimated ranges from 0.60 to 1.34 Gt-C⋅yr−1 depending on different formulations used for wind speed dependence on the gas transfer coefficient. These estimates are subject to an error of up to 75% resulting from the numerical interpolation method used to estimate the distribution of ΔpCO2 over the global oceans. Temperate and polar oceans of the both hemispheres are the major sinks for atmospheric CO2, whereas the equatorial oceans are the major sources for CO2. The Atlantic Ocean is the most important CO2 sink, providing about 60% of the global ocean uptake, while the Pacific Ocean is neutral because of its equatorial source flux being balanced by the sink flux of the temperate oceans. The Indian and Southern Oceans take up about 20% each.

Increased tricarboxylic acid cycle flux in rat brain during forepaw stimulation detected with 1H[13C]NMR.

NMR spectroscopy was used to test recent proposals that the additional energy required for brain activation is provided through nonoxidative glycolysis. Using localized NMR spectroscopic methods, the rate of C4-glutamate isotopic turnover from infused [1-(13)C]glucose was measured in the somatosensory cortex of rat brain both at rest and during forepaw stimulation. Analysis of the glutamate turnover data using a mathematical model of cerebral glucose metabolism showed that the tricarboxylic acid cycle flux [(V(TCA)] increased from 0.49 +/- 0.03 at rest to 1.48 +/- 0.82 micromol/g/min during stimulation (P < 0.01). The minimum fraction of C4-glutamate derived from C1-glucose was approximately 75%, and this fraction was found in both the resting and stimulated rats. Hence, the percentage increase in oxidative cerebral metabolic rate of glucose use (CMRglc) equals the percentage increases in V(TCA) and cerebral metabolic rate of oxygen consumption (CMRO2). Comparison with previous work for the same rat model, which measured total CMRglc [Ueki, M., Linn, F. & Hossman, K. A. (1988) J. Cereb. Blood Flow Metab. 8, 486-4941, indicates that oxidative CMRglc supplies the majority of energy during sustained brain activation.