13C tracer experiments and metabolite balancing for metabolic flux analysis: comparing two approaches

Conventional metabolic flux analysis uses the information gained from determination of measurable fluxes and a steady-state assumption for intracellular metabolites to calculate the metabolic fluxes in a given metabolic network. The determination of intracellular fluxes depends heavily on the correctness of the assumed stoichiometry including the presence of all reactions with a noticeable impact on the model metabolite balances. Determination of fluxes in complex metabolic networks often requires the inclusion of NADH and NADPH balances, which are subject to controversial debate. Transhydrogenation reactions that transfer reduction equivalents from NADH to NADPH or vice versa can usually not be included in the stoichiometric model, because they result in singularities in the stoichiometric matrix. However, it is the NADPH balance that, to a large extent, determines the calculated flux through the pentose phosphate pathway. Hence, wrong assumptions on the presence or activity of transhydrogenation reactions will result in wrong estimations of the intracellular flux distribution. Using 13C tracer experiments and NMR analysis, flux analysis can be performed on the basis of only well established stoichiometric equations and measurements of the labeling state of intracellular metabolites. Neither NADH/NADPH balancing nor assumptions on energy yields need to be included to determine the intracellular fluxes. Because metabolite balancing methods and the use of 13C labeling measurements are two different approaches to the determination of intracellular fluxes, both methods can be used to verify each other or to discuss the origin and significance of deviations in the results. Flux analysis based entirely on metabolite balancing and flux analysis, including labeling information, have been performed independently for a wild-type strain of Aspergillus oryzae producing alpha-amylase. Two different nitrogen sources, NH4+ and NO3-, have been used to investigate the influence of the NADPH requirements on the intracellular flux distribution. The two different approaches to the calculation of fluxes are compared and deviations in the results are discussed. Copyright 1998 John Wiley & Sons, Inc.     

The crystal structure of a quercetin 2,3-dioxygenase from Bacillus subtilis suggests modulation of enzyme activity by a change in the metal ion at the active site(s)

Common structural motifs, such as the cupin domains, are found in enzymes performing different biochemical functions while retaining a similar active site configuration and structural scaffold. The soil bacterium Bacillus subtilis has 20 cupin genes (0.5% of the total genome) with up to 14% of its genes in the form of doublets, thus making it an attractive system for studying the effects of gene duplication. There are four bicupins in B. subtilis encoded by the genes yvrK, yoaN, yxaG, and ywfC. The gene products of yvrK and yoaN function as oxalate decarboxylases with a manganese ion at the active site(s), whereas YwfC is a bacitracin synthetase. Here we present the crystal structure of YxaG, a novel iron-containing quercetin 2,3-dioxygenase with one active site in each cupin domain. Yxag is a dimer, both in solution and in the crystal. The crystal structure shows that the coordination geometry of the Fe ion is different in the two active sites of YxaG. Replacement of the iron at the active site with other metal ions suggests modulation of enzymatic activity in accordance with the Irving-Williams observation on the stability of metal ion complexes. This observation, along with a comparison with the crystal structure of YvrK determined recently, has allowed for a detailed structure-function analysis of the active site, providing clues to the diversification of function in the bicupin family of proteins.     

Intramolecular cooperativity in a protein binding site assessed by combinatorial shotgun scanning mutagenesis

Combinatorial shotgun alanine-scanning was used to assess intramolecular cooperativity in the high affinity site (site 1) of human growth hormone (hGH) for binding to its receptor. A total of 19 side-chains were analyzed and statistically significant data were obtained for 145 of the 171 side-chain pairs. The analysis revealed that 90% of the side-chain pairs exhibited no statistically significant pair interactions, and the remaining 10% of side-chain pairs exhibited only small interactions corresponding to cooperative interaction energies with magnitudes less than 0.4 kcal/mol. The statistical predictions were tested by measuring affinities for purified mutant proteins and were found to be accurate for five of six side-chain pairs tested. The results reveal that hGH site 1 behaves in a highly additive manner and suggest that shotgun scanning should be useful for assessing cooperative effects in other protein-protein interactions.     

Crystal structures of subtilisin BPN' variants containing disulfide bonds and cavities: concerted structural rearrangements induced by mutagenesis

The X-ray structures of four genetically engineered disulfide variants of subtilisin have been analyzed to determine the energetic and structural constraints involved in inserting disulfide bonds into proteins. Each of the engineered disulfides exhibited atypical sets of dihedral angles compared with known structures of natural disulfide bridges in proteins and affected its local structural environment to a different extent. The disulfides located in buried regions, Cys26-Cys232 and Cys29-Cys119, induced larger changes than did Cys24-Cys87 and Cys22-Cys87, which are located on the surface of the molecule. An analysis of the concerted changes in secondary structure units such as alpha-helices and beta-sheets indicated systematic long-range effects. The observed changes in the mutants were largely distributed asymmetrically around the inserted disulfides, reflecting different degrees of inherent flexibility of neighboring secondary structure types. The disulfide substitution in each variant molecule created some invaginations or cavities, causing a reorganization of the surrounding water structure. These changes are described, as well as the changes in side chain positions of groups that border the cavities.     

Analysis of solvent structure in proteins using neutron D2O-H2O solvent maps: pattern of primary and secondary hydration of trypsin.

A method of determining the water structure in protein crystals is described using neutron solvent difference maps. These maps are obtained by comparing the changes in diffracted intensities between two data sets, one in which H2O is the major solvent constituent, and a second in which D2O is the solvent medium. To a good first approximation, the protein atom contributions to the scattering intensities in both data sets are equal and cancel, but since H2O and D2O have very different neutron-scattering properties, their differences are accentuated to reveal an accurate representation of the solvent structure. The method also employs a series of density modification steps that impose known physical constraints on the density distribution function in the unit cell by making real space modifications directly to the density maps. Important attributes of the method are that (1) it is less subjective in the assignment of water positions than X-ray analysis; (2) there is threefold improvement in the signal-to-noise ratio for the solvent density; and (3) the iterative density modification produces a low-biased representation of the solvent density. Tests showed that water molecules with as low as 10% occupancy could be confidently assigned. About 300 water sites were assigned for trypsin from the refined solvent density; 140 of these sites were defined in the maps as discrete peaks, while the remaining were found within less-ordered channels of density. There is a very good correspondence between the sites in the primary hydration layer and waters found in the X-ray structure. Most water sites are clustered into H-bonding networks, many of which are found along intermolecular contact zones. The bound water is equally distributed between contacting apolar and polar atoms at the protein interface. A common occurrence at hydrophobic surfaces is that apolar atoms are circumvented by one or more waters that are part of a larger water network. When the effects on surface accessibility by neighboring molecules in the crystal lattice are taken into consideration, only about 29% of the surface does not interface ordered water. About 25% of the ordered water is found in the second hydration sphere. In many instances these waters bridge larger clusters of primary layer waters. It is apparent that, in certain regions of the crystal, the organization of ordered water reflects the characteristics of the crystal environment more than those of trypsin's surface alone.     

Variability of conformations at crystal contacts in BPTI represent true low-energy structures: correspondence among lattice packing and molecular dynamics structures.

The structures of five basic pancreatic trypsin inhibitor (BPTI) molecules are compared to establish the extent and nature of the conformational variability resulting from crystal packing effects. BPTI is an ideal system to evaluate such factors because of the availability of high resolution X-ray models of five different BPTI structures, each in a different crystal packing environment. Differences observed among the structures are found to be distributed throughout the molecule, although the regions that display most variability are associated with the loop structures (residues 14-17 and 24-29). The regions of structure that show the largest rms deviations from the mean of the five packing motifs correlate well with the presence of intermolecular contacts in the crystal lattice. For most of the molecules there is also a correspondence between a larger number of intermolecular contacts and systematically higher B-factors, although it is not apparent whether this is induced by the crystal contact or results from the fact that the contacts are made predominantly through surface loops. The conformational differences seen among the X-ray models constitute more than local shifts at the lattice contact surfaces, and in fact involve in some cases the making and breaking of intramolecular H-bonds. The magnitudes of the differences among packing models are significantly larger than those usually associated with changes induced by mutagenesis; for instance; the structural differences at the site of mutation observed on removing an internal disulfide from the molecule are significantly less than those associated with lattice contact effects. The crystal packing conformations are compared to representative structures of BPTI generated during a 96-psec molecular dynamics (MD) simulation. This comparison shows a high level of correspondence between the protein flexibility indicated by the X-ray and MD analyses, and specifically between those regions that are most variable. This suggests that the regions that show most variability among the crystal packing models are not artifacts of crystallization, but rather represent true low-energy conformers that have been preferentially selected by crystallization factors.     

Population parameters and exploitation rate of Sarotherodon galilaeus galilaeus (Cichlidae) in Lakes Doukon and Togbadji,Benin

Growth, mortality, recruitment and relative yield per recruit of Sarotherodon galilaeus galilaeus from Lakes Doukon and Togbadji were studied. Data on total length, total weight and sex were recorded on a monthly basis between January and December 2013 for S. g. galilaeus captured by local fishers. The estimated asymptotic lengths L_∞ were 26.2 and 23.6?cm for Lakes Doukon and Togbadji, respectively, while the growth rate K was 0.73 in Lake Doukon and 0.87 in Lake Togbadji. Estimates of fishing mortality, 0.27 and 0.47 y^?1 for Doukon and Togbadji, respectively, were low relative to natural mortality, 1.51 and 1.74 y^?1, respectively. Sizes at first sexual maturity were 12.8 and 13.2?cm for females and males, respectively, in Lake Doukon, and 11.5 and 12.4?cm for females and males, respectively, in Lake Togbadji. The size at first capture was estimated at 13.3 and 12.7?cm for Lakes Doukon and Togbadji, respectively, which, in the light of the size at maturity estimates, indicates that fish spawn at least once before capture. The current exploitation rates of 0.15 for Lake Doukon and 0.21 for Lake Togbadji suggest that their stocks of S. g. galilaeus are not overexploited in either lake.     

Histopathological biomarkers and genotoxicity in gill and liver tissues of Nile tilapia Oreochromis niloticus from a polluted part of the Nile River,Egypt

Fish health is affected by water pollution. Oreochromis niloticus collected during summer 2014 from El-Serw, a polluted site on the Nile River, were compared with fish from a reference site, El-Zamalek. Histopathological changes were detected in gill and liver tissue samples using light and electron microscopy. In addition, the degree of DNA damage was measured using the comet assay. To indicate the severity of water pollution at the two sites, physico-chemical properties and heavy metal concentrations were investigated. Gill damage, including lamellar cell hyperplasia and aneurysm, was observed in the fish samples from the polluted site. The livers of fish from the polluted area showed necrosis and an increase in melanomacrophage centres. Histochemical results confirmed a marked rise of gill mucopolysaccharides and a reduction of carbohydrate stored in hepatocytes. Electron microscopy revealed clear alterations in gill and liver tissue of fish from the polluted site. The comet assay showed highly significant DNA damage in tilapia collected from the polluted site, compared to those from the reference site. Histopathological biomarkers and the comet assay may therefore be sensitive indicators of exposure to mixtures of aquatic pollutants in Nile tilapia.