General method for the derivation and numerical solution of epithelial transport models

Summary A general method is presented for the formulation and numerical evaluation of mathematical models describing epithelial transport. The method is based on the principles of conservation of mass, and maintenance of electroneutrality within the cells and bathing solutions. It is therefore independent of the specific membrane transport mechanisms, and can be used to evaluate different models describing arbitrary transport processes (including passive, active and cotransport processes). Detailed numerical methods are presented that allow computation of steady-state and transient responses under open-circuit, current-clamp and voltage-clamp conditions, using a general-purpose laboratory minicomputer. To evaluate the utility of this approach, a specific model is presented that is consistent with the Koefoed-Johnson and Ussing hypothesis of sodium transport in tight epithelia (Acta Physiol. Scand. 42:298–308, 1958). This model considers passive transport of an arbitrary number of permeant solutes, active transport of sodium and potassium, and osmotically induced water transport across the apical and basolateral membranes. Results of the model are compared to published experimental measurements in rabbit urinary bladder epithelium.     

Growth and development of immature life stages ofPropylaea 14-punctata L. andCoccinella 7-punctata L. [Col.: Coccinellidae] simulated by the metabolic pool model

The conversion of aphid prey tissue (Acyrthosiphon pisum Harris) into predator biomass (immature life stages ofPropylaea 14-punctata L. andCoccinella 7-punctata L.) is calculated by plotting weight gain against assimilation (i.e. consumption minus egestion). This concept is added to the metabolic pool model byGutierrez et al. (1981) that enables the simulation of growth and development of a predator on a physiological basis. Physiological time is expressed in daydegrees above lower development thresholds for both species. Visual examination of observed and calculated values showed that the model satisfactorily describes the growth patterns of the above predators.      

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     

Comparison of the ribosomal RNA genes in four closely relatedCucurbitaceae

Cucurbitaceae are characterized by a high copy number for nuclear ribosomal RNA genes. We have investigated the genomic ribosomal DNA (rDNA) of four closely related species of this family with respect to structure, length heterogeneity, and evolution. InCucumis melo (melon) there are two main length variants of rDNA repeats with 10.7 and 10.55kb.Cucumis sativus (cucumber) shows at least three repeat types with 11.5, 10.5, and 10.2kb.Cucurbita pepo (zucchini) has two different repeat types with 10.0 and 9.3kb. There are also two different repeat types inCucurbita maxima (pumpkin) of about 11.2 and 10.5kb. Restriction enzyme mapping of the genomic rDNA of these four plants and of cloned repeats ofC. sativus shows further heterogeneities which are due to methylation or point mutations. By comparison of the restriction enzyme maps it was possible to trace some evolutionary events in the family ofCucurbitaceae. Some aspects of regulation and function of the middle repetitive rRNA genes (here between 2000 and 10000 copies) are discussed.     

The versatility of proteolytic enzymes

The growing realization of their physiological importance has generated renewed interest in the study of proteolytic enzymes. Modern methods of protein chemistry and molecular biology have revealed new insights into the protein and gene structure of a variety of protein precursors and their processing by limited proteolysis. Examples are given in this review for transmembrane processes and the role of signal peptidases of both eukaryotic and prokaryotic origin, the processing of prohormones and precursors of growth factors, protein components of blood coagulation, fibrinolysis, and of the complement system, and a group of granulocyte proteases, including the mast cell serine proteases. The relationship of homologous domains found in many of these proteases and their zymogens to protein evolution is a recurrent theme of this discussion.     

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.     

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.