Molecular dynamics simulations of phenylimidazole inhibitor complexes of cytochrome P450 cam

Molecular dynamics simulations have been performed on three phenylimidazole inhibitor complexes ofP450 _cam, utilizing the X-ray structures and the AMBER suite of programs. Compared to their corresponding optimized X-ray structures, very similar features were observed for the 1-phenylimidazole (1-PI) and 2-phenylimidazole (2-PI) complexes during a 100 ps MD simulation. The 1-PI inhibitor binds as a Type II complex with the imidazole nitrogen as a ligand of the heme iron. Analysis of the inhibitor-enzyme interctions during the MD simulations reveals that electrostatic interactions of the imidazole with the heme and van der Waals interactions of the phenyl ring with nearby hydrophobic residues are dominant. By contrast, 2-PI binds as a Type I inhibitor in the substrate binding pocket, but not as a ligand of the iron. The interactions of this inhibitor are qualitatively different from that of the Type II 1-PI, being mainly electrostatic/H-bonding interactions with a bound water and polar residues. Although the third compound, 4-PI, in common with 1-PI, also binds as a Type II inhibitor, with one nitrogen of the imidazole as a ligand to the iron, the MD average binding orientation deviates significantly from the X-ray structure. The most important changes observed include: (1) the rotation of the imidazole ring of this inhibitor by about 90° to enhance electrostatic interactions of the imidazole NH group with the carbonyl group of LEU244, and (2) the rotation of the carbonyl group of ASP251 to form a H-bond with VAL254. An analysis of the H-bonding network surrounding this substrate in the optimized crystal structure revealed that there is no H-bonding partner either for the free polar NH group in the imidazole ring of 4-phenylimidazole or for the polar carbonyl group of the nearby ASP251 residue. The deviation of the dynamically averaged inhibitor-enzyme structure of the 4-PI complex from the optimized crystal structure can therefore be rationalized as a consequence of the optimization of the electrostatic interactions among the polar groups.     

Computer simulation and interpretation of45Ca efflux profile patterns

Summary Stimulations or inhibitions by various agents of⁴⁵Ca efflux from prelabeled cells or tissues display distinct and reproducible profile patterns when the results are plotted against time as fractional efflux ratios (FER). FER is the fractional efflux of⁴⁵Ca from stimulated cells divided by the fractional efflux from a control unstimulated group. These profile patterns fall into three categories: peak patterns, exponential patterns, and mixed patterns. Each category can be positive (stimulation) or negative (inhibition). The interpretation of these profiles is difficult because⁴⁵Ca efflux depends on three variables: the rate of calcium transport out of the cell, the specific activity of the cell compartment from which the calcium originates, and the concentration of free calcium in this compartment. A computer model based on data obtained by kinetic analyses of⁴⁵Ca desaturation curves and consisting of two distinct intracellular pools was designed to follow the concentration of the traced substance (⁴⁰Ca), the tracer (⁴⁵Ca), and the specific activity of each compartment before, during, and after the stimulation or the inhibition of calcium fluxes at various pool boundaries. The computer model can reproduce all the FER profiles obtained experimentally and bring information which may be helpful to the interpretation of this type of data. Some predictions of the model were tested experimentally, and the results support the views that a peak pattern may reflect a sustained change in calcium transport across the plasma membrane, that an exponential pattern arises from calcium mobilization from an internal subcellular pool, and that a mixed pattern may be caused by a simultaneous change in calcium fluxes at both compartment boundaries.     

Infanticide and fertility among Eskimos: A computer simulation

Until recently, certain Eskimo groups were reported to practice female infanticide in the belief that the time spent suckling a girl would delay the mother's next opportunity to bear a son, males being preferred to females because of their future role as providers in a hunting economy. From sex ratios in census data, rates of female infanticide of up to 66% for some groups have been inferred, leading some ethnographers to conclude that these groups were headed for extinction. Eskimo beliefs regarding the effects of infanticide on fertility, however, are in accord with the results of research on the relation of fertility and lactation: The cessation of lactation following infanticide would significantly shortern the expected interval until the next birth. Given this fact and available field data regarding the parameters of Eskimo population growth, the present computer simulation indicates that Eskimo populations could sustain a rate of 30% female infanticide and still survive. Higher reported rates are explained as the combined result of female infanticide plus the tendency of ethnographers to overestimate to overestimate the ages of juvenile females relative to juvenile males.     

Rate-controlling steps of oxidative phosphorylation in rat liver mitochondria. A synoptic approach of model and experiment

The contribution of different steps to the control of oxidative phosphorylation in isolated rat liver mitochondria was investigated by a combination of experiments and computer simulations. The parameters of the mathematical model of phosphorylating mitochondria were derived from experimental data. The model correctly describes the competition between ATP utilization inside and outside mitochondria for the ATP generated in mitochondria. On the basis of the good agreement between experiments and simulations, the contribution of different steps to the control of respiration was estimated by computing their control strengths, i.e., the influence of their activities on the rate of respiration. The rate-controlling influences vary depending on the load of oxidative phosphorylation. The predominant steps are: in the fully active state (State 3) — the hydrogen supply to the respiratory chain; in the resting state (State 4) — the proton leak of the mitochondrial inner membrane; in states of non-maximum ATP export — the adenine nucleotide translocator. Titrations of respiration with phenylsuccinate, antimycin, oligomycin and carboxyatractyloside completely support these conclusions.     

Simulation of fermentation conditions for vitamin B12 biosynthesis from whey

Vitamin B₁₂ production in fermentation of Propionibacterium shermanii and Propionibacterium arl AKU 1251 in whey permeate medium has been studied. The observed results and simulated expected values obtained by fitting statistical equations to the recorded data showed that 24 h old inoculum, 5 mg iron l^?1 and 4% whey lactose were optimal for vitamin B₁₂ biosynthesis in both strains when fermentation was carried out under anerobic (84 h) and aerobic (84 h) conditions at 30°C. The supplementation of whey medium with 0.5% (NH₄)₂HPO₄ enhanced further the metabolite yield; however, the preference for a mixed carbon source (lactose + d-glucose or lactose + d-fructose) at different levels varied in the strains under study. P. shermanii, under optimal cultural conditions, was found to be a better strain than Propionibacterium arl AKU, 1251 in fermenting whey lactose for product (vitamin B₁₂) formation.     

Periodic operation of immobilized live cell bioreactor for the production of candicidin

A lumped model for cell growth and secondary metabolite production in an immobilized live cell bioreactor has been developed. This model is applied here to simulate the performance of an immobilized bioreactor under steady-state conditions and under conditions of periodically varying concentration of a growth-limiting substrate. The results of the simulation study were experimentally verified in the case of the production of the antibiotic candicidin by Streptomyces griseus in an immobilized bioreactor with forced periodic operation. The results of the studies suggest that periodically operated immobilized live cell bioreactors can provide a potent alternative for the production of non-growth-associated biochemicals, as compared to free cell fermentations, pulsed fermentations with process cycle regeneration, and nonregenerated bioreactors. This work has demonstrated that by frequent pulsing of the growth limiting nutrient, stable extended production can be obtained at high specific cellular productivities.     

Responses to selection on genotypic or phenotypic values in the presence of genes with major effects

Summary Average genotypic responses were compared after selection for genotypic values and for phenotypic values on the basis of single-gene models and multigene models in simulated livestock populations. Single-gene models dealt with single gene control of the genetic differences between animals, while multigene models considered a collection of genes with various magnitudes of effects on a trait. In each case, selection lasted through discrete generations until the fixation of the gene frequencies occurred. Generations to reach fixation were used to compare various models, and the two criteria for selection, for their efficiency in selection. Implications of using these models versus using infinitesimal models for selection in practice are presented.     

Estimates of marker-associated QTL effects in Monte Carlo backcross generations using multiple regression

Summary The decision of whether or not to use QTLassociated markers in breeding programs needs further information about the magnitude of the additive and dominance effects that can be estimated. The objectives of this paper are (1) to apply some of the Moreno-Gonzalez (1993) genetic models to backcross simulation data generated by the Monte Carlo method, and (2) to get simulation information about the number of testing progenies and mapping density in relation to the magnitude of gene effect estimates. Results of the Monte Carlo study show that the stepwise regression analysis was able to detect relatively small additive and dominance effects when the QTL are independently segregating. When testing selfed families derived from backcross individuals, dominance effects had a larger error standard deviation and were estimated at a lower frequency. Linked QTL require a higher marker mapping density on the genome and a larger number of progenies to detect small genetic effects. Reduction of the environmental error variance by evaluating selfed backcross families in replicate experiments increased the power of the test. Expressions of the number of progenies for detecting significant additive effects were developed for some genetic situations. The ratio of the within-backcross genetic variance to the square of a gene effect estimate is a function of the number of progenies, the heritability of the trait, the marker map density and the portion of the genetic variance explained by the model. Different values (from 0 to 1) assigned to (relative position of the QTL in the marker segment) did not cause a large shift in the residual mean square of the model.     

Kin-structured migration and the rate of advance of an advantageous gene

Hemoglobin E, an allele generally considered to confer malarial resistance in heterozygotes, is found in high frequencies across a wide area of Southeast Asia. Apparently it originated as a single-point mutation which was spread by gene flow. The rate of diffusion of this adaptive allele is studied using four computer simulation models. It is shown that in small populations deterministic equations for gene flow may overestimate rates of diffusion. Other aspects of population structure, however, such as kin-structuring of migrant groups, increase the rate of advance. Finally, population growth coupled with the spread of the allele leads to much more rapid diffusion. These results suggest that population structure can be an important factor affecting the diffusion of advantageous genes.