A Monte Carlo method for the simulation of kinetie models

The purpose of this note is to illustrate the feasibility of simulating kinetic systems, such as commonly encountered in photosynthesis research, using the Monte Carlo (MC) method. In this approach, chemical events are considered at the molecular level where they occur randomly and the macroscopic kinetic evolution results from averaging a large number of such events. Their repeated simulation is easily accomplished using digital computing. It is shown that the MC approach is well suited to the capabilities and resources of modern microcomputers. A software package is briefly described and discussed, allowing a simple programming of any kinetic model system and its resolution. The execution is reasonably fast and accurate; it is not subject to such instabilities as found with the conventional analytical approach.Abbreviations MC Monte Carlo - RN random number - PSU photosynthetic unit Dedicated to Prof. L.N.M. Duysens on the occasion of his retirement.     

Manganese-histidine cluster as the functional center of the water oxidation complex in photosynthesis

The recent model of Kambara and Govindjee for water oxidation [Kambara T. and Govindjee (1985) Proc. Natl. Acad. Sci. U.S.A., 82:6119–6123] has been extended in this paper by examining all the data in order to identify the most likely candidate for the redox-active ligand (RAL), suggested to operate between the water oxidizing complex (WOC) and Z, the electron donor to the reaction center P680. We have concluded that a very suitable candidate for RAL is the imidazole moiety of a histidine residue. The electrochemical data available on imidazole derivatives play heavily in this identification of RAL. Thus, we suggest that histidine might play the role of an electron mediator between the WOC and Z. A model of S-states in terms of their plausible chemical identity is presented here.Abbreviations J electronic spin of ion - P680 reaction center chlorophyll - RAL Redox active ligand - S_n state of the oxygen-evolving system - WOC water oxidation complex - Z electron donor to P680 Dedicated to Prof. L.N.M. Duysens on the occasion of his retirement     

Cycling of soil carbon in a Japanese red pine forest II. Changes occurring in the first year after a clear-felling

Cycling of soil carbon in the first year after a clear-felling was compared with that before the felling in a Japanese red pine forest in Hiroshima Prefecture, west Japan. The daily mean temperature at the soil surface in summer was increased after the felling in comparison to that before felling, and the water content of both the A₀ layer and the surface mineral soil was decreased due to the loss of the forest canopy. The rate of weight loss of the A₀ layer was reduced after felling. However, accumulation of the A₀ layer rapidly decreased because of the lack of litter supply to the forest floor. Low soil respiration after felling was mainly caused by the cessation of root respiration. Analysis of annual soil carbon cycling was then conducted using a compartment model. The relative decomposition rate of the A₀ layer decreased whereas that of humus and dead roots in mineral soil increased to some extent after felling. The accumulation of carbon in mineral soil, however, increased slightly due to the supply of humus from roots killed by the felling.     

A soybean growth model based on compensatory growth following defoliation

On the basis of an artificial defoliation experiment, a new growth model of soybean was formulated through a modification of Rudd's (1980) model with regard to his equations for dry matter allocation. Compensatory growth for leaf damage was modelled by a single process in which the dry matter allocation changes dynamically according to the severity of leaf damage. The sums of squared differences between simulations and experimental soybean yields were much smaller in our modified model than in Rudd's original model. The modified model gave a better simulation of yield loss due to defoliation that varied in time and intensity. The relationship between various times and intensities of defoliation and yield loss was shown, which is essential for establishing the dynamic economic injury level in IPM.     

Isolation of a brain peptide identical to the intestinal PHI (peptide HI)

The isolation of a brain peptide identical to the intestinal peptide PHI (peptide HI) is described. The peptide was isolated from porcine brain extract using a chemical assay method based on its C-terminal isoleucine amide structure. The complete amino acid sequence of the peptide was found to be: His-Ala-Asp-Gly-Val-Phe-Thr-Ser-Asp-Phe-Ser-Arg-Leu-Leu-Gly-Gln-Leu-Ser-Ala- Lys-Lys-Tyr-Leu-Glu-Ser-Leu-Ile-NH2. This sequence is identical to the intestinal peptide thus demonstrating PHI to be a brain-gut peptide. The role of PHI in the central nervous system as a neurotransmitter or neuromodulator is discussed.     

Theoretical conformational analysis of phospholipids II. Role of hydration in the gel to liquid crystal transition of phospholipids

To obtain a satisfactory agreement between computed transition temperatures and those determined experimentally, we introduce explicitly water molecules which hydrate the polar headgroup of dipalmitoylphosphatidylethanolamine molecules. The calculated free energy curves as a function of the intermolecular interchain distance and the degree of hydration of the polar groups permit the determination of the transition of the phospholipid system from the gel to the liquid crystalline phase. The detailed structure of the hydration shell is defined using the supermolecular approach.     

Dopamine Uptake by Rat Striatal Synaptosomes: A Compartmental Analysis

Abstract: Dopamine (DA) uptake into synaptosomes from rat corpus striatum was studied in the presence of a monoamine oxidase (MAO) inhibitor and dithiothreitol, by means of a filtration technique. Under these conditions a steady state develops rapidly in which the synaptosomal DA content remains constant while the continuing DA uptake is counterbalanced by DA efflux from the synaptosome. Exchange of synaptosomal [³H]DA and [¹⁴C]DA was measured under these conditions. In timecourse experiments it was found that exchange could be described significantly better by a three-compartment model than by a two-compartment model. However, if synaptosomes from reserpine-pretreated animals were used, analysis according to a three-compartment model did not result in a significantly better fit compared with a two-compartment model. Subsequently, kinetic transfer parameters describing DA fluxes between compartments at different DA concentrations were calculated from the fitted exchange curves. A Michaelis-Menten kinetic analysis indicated that only the in-series three-compartment configuration, in which DA is taken up from the medium into one synaptosomal compartment, from which it can subsequently be transferred to a second compartment without direct access to the medium, gave kinetically acceptable results. Transfer parameters in synaptosomes from reserpine-treated rats were comparable to those parameters describing DA transport between the medium and the first intrasynaptosomal compartment as measured under control conditions. Morover, it was found that potassium depolarization of synaptosomes resulted in a release of DA in a quantity similar to that found in the second intrasynaptosomal compartment. It is suggested that the two intrasynaptosomal compartments found correspond to a cytoplasmatic and vesicular DA pool, respectively. The functional significance of these findings is discussed in terms of the regulation of DA levels within the nerve terminal.     

Modelling Mediterranean pasture dynamies

Pasture vegetation in an open woodland of Quercus rotundifolia subjected to periodic ploughing was sampled in spring during 8 consecutive years. The frequency of herbaceous species was recorded in a total of 69 permanent plots located on 5 adjacent sites with similar lithology, slope and orientation but differing in age since previous ploughing.Vegetation dynamics expressed as trajectories of permanent plots in a non-metric multidimensional scaling space has been modelled in terms of evironmental variables. By fitting a generalized linear model, the dynamics are shown to be related to years since last ploughing, geographical location of plots, total annual rainfall and November rainfall. Meteorological patterns of the sampling period are also described.Abbreviations GLM = Generalized Linear Model - NMDS = Non-metric Multi-Dimensional Scaling - UPGMA = Unweighted Pair-Group Method using Arithmetic Averages     

Growth and production of O+ bream (Abramis brama), O+ roach (Rutilus rutilus) and O+ carp (Cyprinus carpio) in 10 drainable 0.1 Ha ponds

In ten 0.1 ha drainable ponds O⁺ cyprinids were stocked to assess their impact on the aquatic community. Bream and roach were stocked in early June, 1986, and carp in early July, 1986. The ponds were drained mid November, 1986. The relationship between the fish community and its food resource was estimated in terms of production and consumption with use of the data on growth and mortality of the fish. The average gross production of the fish species in the ten ponds (48.4 kg during the experiment) was dominated by the carp (58.5%). A minimum estimate of the food consumption of the O⁺ cyprinids is 760 kJ.m^–2. The average production of bream and roach was much lower than that of carp due to a higher mortality and a more specialized feeding pattern. The gross productions of carp, bream and roach in the ten ponds were not significantly related, although stomach analysis showed diet overlap between the carp and the bream.