Conditional and Unconditional Simulation of Healthy Patients' Visual Fields

This paper describes a simulation problem, motivated by the study of glaucoma, a very serious and widespread ocular illness. To ascertain whether a patient suffers from glaucoma, a perimetric test is done, but the evolution of the disease is very slow, and large longitudinal sets of tests taken on the same patient are needed to study its evolution, to analyze the efficiency of existing methods to detect the progression of glaucoma and to develop new ones. Simulation can be a very useful procedure to get appropriate data sets to work with. Our aim in this work is to simulate several VFs in a healthy patient to reflect his evolution in time. We use a spatio‐temporal model to simulate from, taking into account the correlation existing between the observed (or simulated) values in space and time. Two different simulation procedures (unconditional and conditional) are studied, and applied to obtain the simulations we are interested in. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Big Moose Basin: simulation of response to acidic deposition

The ILWAS model has been enhanced for application to multiple-lake hydrologic basins. This version of the model has been applied to the Big Moose basin, which includes Big Moose Lake and its tributary streams, lakes, and watersheds. The basin, as defined, includes an area of 96 km², with over 20 lakes and ponds, and 70 km of streams. Hydrologic and chemical calibrations have been made using data from seven sampling stations. When total atmospheric sulfur loading to the basin is halved, the model predicts, after four years of simulation, a decreasing sulfate concentration and to a lesser extent a rising alkalinity at Big Moose Lake outlet. At the end of four years, the results show an increase in pH of 0.1 to 0.5 pH units depending upon season.     

Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base and the associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. With one exception we find that for all globin sequences one of the known globin folds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structural research and future development of our approach.     

Application of a directed conformational search for generating 3-D coordinates for protein structures from alpha-carbon coordinates.

A directed conformational search algorithm using the program CONGEN (ref. 3), which samples backbone conformers, is described. The search technique uses information from the partially built structures to direct the search process and is tested on the problem of generating a full set of backbone Cartesian coordinates given only alpha-carbon coordinates. The method has been tested on six proteins of known structure, varying in size and classification, and was able to generate the original backbone coordinates with RMSs ranging from 0.30-0.87A for the alpha-carbons and 0.5-0.99A RMSs for the backbone atoms. Cis peptide linkages were also correctly identified. The procedure was also applied to two proteins available with only alpha-carbon coordinates in the Brookhaven Protein Data Bank; thioredoxin (SRX) and triacylglycerol acylhydrolase (TGL). All-atom models are proposed for the backbone of both these proteins. In addition, the technique was applied to randomized coordinates of flavodoxin to assess the effects of irregularities in the data on the final RMS. This study represents the first time a deterministic conformational search was used on such a large scale.     

Theory of chaperonin action: inertial model for enhancement of prokaryotic Rubisco assembly.

We have performed a computational simulation of the aggregation and chaperonin-dependent reconstitution of dimeric prokaryotic ribulose bisphosphate carboxylase/oxygenase (Rubisco), based on the data of P. Goloubinoff et al. (1989, Nature 342, 884-889) and P. V. Viitanen et al. (1990, Biochemistry 29, 5665-5671). The aggregation is simulated by a set of 12 differential equations representing the aggregation of the Rubisco folding intermediate, Rubisco-I, with itself and with aggregates of Rubisco-I, leading up to dodecamers. Four rate constants, applying to forward or reverse steps in the aggregation process, were included. Optimal values for these constants were determined using the ellipsoid algorithm as implemented by one of us (Ecker, J.G. & Kupferschmid, M., 1988, Introduction to Operations Research, Wiley, New York, pp. 315-322). Intensive exploration of simpler aggregation models did not identify an alternative that could simulate the data as well as this one. The activity of the chaperonin in this system was simulated by using this aggregation model, combined with a model similar to that proposed by Goloubinoff et al. (1989). The model assumes that the chaperonin can bind the folding intermediate rapidly, and that the chaperonin complex releases the Rubisco molecule slowly, permitting time for its spontaneous folding while interacting with the chaperonin. This is followed by self-association of the folded Rubisco monomer to yield the active dimeric Rubisco. A modification of the model that simulates temperature effects was also constructed. The most important results we obtained indicate that the chaperonin-dependent reconstitution of Rubisco can be simulated adequately without invoking any catalysis of folding by the chaperonin. In addition, the simulations predict values for the association rate constant of Rubisco-I with the chaperonin, and other variables, that are subject to experimental verification.     

Evaluation of two unstructured mathematical models for the penicillin G fedbatch fermentation

The mathematical model for the penicillin G fed-batch fermentation proposed by Heijnen et al. (1979) is compared with the model of Bajpai & Reuß (1980). Although the general structure of these models is similar, the difference in metabolic assumptions and specific growth and production kinetics results in a completely different behaviour towards product optimization. A detailed analysis of both models reveals some physical and biochemical shortcomings. It is shown that it is impossible to make a reliable estimation of the model parameters, only using experimental data of simple constant glucose feed rate fermentations with low initial substrate amount. However, it is demonstrated that some model parameters might be key factors in concluding whether or not altering the substrate feeding strategy has an important influence on the final amount of product.It is illustrated that feeding strategy optimization studies can be a tool in designing experiments for parameter estimation purposes.     

Aluminum geochemistry in peatland waters

The chemical speciation of aluminum was examined in surface water samples from Sphagnum peatlands in north-central Minnesota, from peatlands along the Canadian east coast, and from bogs in the Pennine Mountain area of England. In highly organic ([DOC] 50 mg L^–1 ), low pH waters, 80–90% of total dissolved Al was complexed with organic matter (OM), while in waters with low DOC ([DOC] 5 mg L^–1) 54–86% of total dissolved Al existed as Al⁺³ or other inorganic Al species. Batch titrations of OM with Al revealed a high Al binding capacity, 1.4–2.8 mol (mg DOC)^–1, that generally was unsaturated with Al. Titrations of OM with Al in conjunction with a continuous distribution model were used to determine Al-OM conditional stability constants. Binding capacity (mol Al (mg DOC)^–1) and strength (formation constant) increased from pH 3 to 5 but decreased above pH 5 due to formation of AI-hydroxy species including A1(OH)₃ (s). The high binding capacity of OM in bog waters facilitates metal mobility, especially in low pH (< 5) wetlands where metal solubility is high and OM concentrations are highest. Results showed that the relative degree of organic matter saturation with metal ions was important in modeling AI speciation in bog waters.     

Evolutionary analysis of the multigene pregnancy-specific β1-Glycoprotein family: Separation of historical and nonhistorical signals

The pregnancy-specific ₁-glycoproteins (PSG) form a large family of closely related proteins. Using newly developed methods of sequence analysis, in combination with protein modeling, we provide a framework for investigating the evolution and biological function of genes like the PSG. Evolutionary trees, based on C-terminal sequence, group PSG genes in a manner consistent with their genomic organization. Trees constructed using the N-terminal domain sequences are unreliable as an indicator of phylogeny because of non-neutral processes of sequence change. During duplication of the PSG genes, evolutionary pressures have resulted in a gradient of constrained change across each gene. The N-terminal domains show a nonrandom pattern of amino acid substitutions clustered in the immunoglobulin complementarity-determining region (CDR)-like regions, which appear to be important in the function of the protein.     

Hybridomas in a bioreactor cascade: modeling and determination of growth and death kinetics

Hybridomas were cultured under steady-state conditions in a series of two continuous stirred-tank reactors (CSTRs), using a serum-free medium. The substrate not completely converted in the first CSTR, was transported with the cells to the second one and very low growth rates, high death rates, and lysis of viable cells were observed in this second CSTR. These conditions are hardly accessible in a single vessel, because such experiments would be extremely time-consuming and unstable due to a low viability. In contrast to what is often observed in literature, kinetic parameters could thus be derived without the neccessity for extrapolation to lower growth rates. Good agreement with literature averages for other hybridomas was found. Furthermore, showing that the reactor series is a valuable research tool for kinetic studies under extreme conditions, the possibility to observe cell death under stable and defined steady-state conditions offers interesting opportunities to investigate apoptosis and necrosis. Additionally, a model was developed that describes hybridoma growth and monoclonal antibody production in the bioreactor cascade on the basis of glutamine metabolism. Good agreement between the model and the experiments was found.Abbreviation MAb Monoclonal antibody Nomenclature C _AConcentration of any (mol m^-3) component A D Dilution rate (s^-1) K _dDeath-rate constant (mol m^-3) K _lLysis-rate constant (mol m^-3) K _sMonod constant (mol m^-3) m Maintenance coefficient (mol cell^-1 s^-1) q Specific consumption (mol cell^-1 s^-1) or production rate t Time (s) X Cell concentration (cell m^-3) Y Yield coefficient (cell mol^-1) Greek symbols _d Specific death rate (s^-1) _l Specific lysis rate (s^-1) of viable cells _net Net specific growth (s^-1) rate _true True specific growth (s^-1) rate