Patch dynamics and metapopulation theory: the case of successional species

We present a mathematical framework that combines extinction-colonization dynamics with the dynamics of patch succession. We draw an analogy between the epidemiological categorization of individuals (infected, susceptible, latent and resistant) and the patch structure of a spatially heterogeneous landscape (occupied-suitable, empty-suitable, occupied-unsuitable and empty-unsuitable). This approach allows one to consider life-history attributes that influence persistence in patchy environments (e.g., longevity, colonization ability) in concert with extrinsic processes (e.g., disturbances, succession) that lead to spatial heterogeneity in patch suitability. It also allows the incorporation of seed banks and other dormant life forms, thus broadening patch occupancy dynamics to include sink habitats. We use the model to investigate how equilibrium patch occupancy is influenced by four critical parameters: colonization rate? extinction rate, disturbance frequency and the rate of habitat succession. This analysis leads to general predictions about how the temporal scaling of patch succession and extinction-colonization dynamics influences long-term persistence. We apply the model to herbaceous, early-successional species that inhabit open patches created by periodic disturbances. We predict the minimum disturbance frequency required far viable management of such species in the Florida scrub ecosystem. (C) 2001 Academic Press.

Diffusion modeling of percutaneous absorption kinetics: 2. Finite vehicle volume and solvent deposited solids

The diffusion model for percutaneous absorption is developed for the specific case of delivery to the skin being limited by the application of a finite amount of solute. Two cases are considered; in the first, there is an application of a finite donor (vehicle) volume, and in the second, there are solvent-deposited solids and a thin vehicle with a high partition coefficient. In both cases, the potential effect of an interfacial resistance at the stratum corneum surface is also considered. As in the previous paper, which was concerned with the application of a constant donor concentration, clearance limitations due to the viable eqidermis, the in vitro sampling rate, or perfusion rate in vivo are included. Numerical inversion of the Laplace domain solutions was used for simulations of solute flux and cumulative amount absorbed and to model specific examples of percutaneous absorption of solvent-deposited solids. It was concluded that numerical inversions of the Laplace domain solutions for a diffusion model of the percutaneous absorption, using standard scientific software (such as SCIENTIST, MicroMath Scientific software) on modern personal computers, is a practical alternative to computation of infinite series solutions. Limits of the Laplace domain solutions were used to define the moments of the flux-time profiles for finite donor volumes and the slope of the terminal log flux-time profile. The mean transit time could be related to the diffusion time through stratum corneum, viable epidermal, and donor diffusion layer resistances and clearance from the receptor phase. Approximate expressions for the time to reach maximum flux (peak time) and maximum flux were also derived. The model was then validated using reported amount-time and flux-time profiles for finite doses applied to the skin. It was concluded that for very small donor phase volume or for very large stratum corneum-vehicle partitioning coefficients (e.g., for solvent deposited solids), the flux and amount of solute absorbed are affected by receptor conditions to a lesser extent than is obvious for a constant donor constant donor concentrations. (C) 2001 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 90:504-520, 2001.

NMR structure and backbone dynamics of a concatemer of epidermal growth factor homology modules of the human low-density lipoprotein receptor

The ligand-binding region of the low-density lipoprotein (LDL) receptor is formed by seven N-terminal, imperfect, cysteine-rich (LB) modules. This segment is followed by an epidermal growth factor precursor homology domain with two N-terminal, tandem, EGF-like modules that are thought to participate in LDL binding and recycling of the endocytosed receptor to the cell surface. EGF-A and the concatemer, EGF-AB, of these modules were expressed in Escherichia coli. Correct protein folding of EGF-A and the concatemer EGF-AB was achieved in the presence or absence of calcium ions, in contrast to the LB modules, which require them for correct folding. Homonuclear and heteronuclear H-1-N-15 NMR spectroscopy at 17.6 T was used to determine the three-dimensional structure of the concatemer. Both modules are formed by two pairs of short, anti-parallel beta -strands. In the concatemer, these modules have a fixed relative orientation, stabilized by calcium ion-binding and hydrophobic interactions at the interface. N-15 longitudinal and transverse relaxation rates, and {H-1}-N-15 heteronuclear NOEs were used to derive a model-free description of the backbone dynamics of the molecule. The concatemer appears relatively rigid, particularly near the calcium ion-binding site at the module interface, with an average generalized order parameter of 0.85 +/- 0.11. Some mutations causing familial hypercholesterolemia may now be rationalized. Mutations of D41, D43 and E44 in the EGF-B calcium ion-binding region may affect the stability of the linker and thus the orientation of the tandem modules. The diminutive core also provides little structural stabilization, necessitating the presence of disulfide bonds. The structure and dynamics of EGF-AB contrast with the N-terminal LB modules, which require calcium ions both for folding to form the correct disulfide connectivities and for maintenance of the folded structure, and are connected by highly mobile linking peptides. (C) 2001 Academic Press.

Identifying trajectories of adolescent smoking: An application of latent growth mixture modeling

The goal of the current study was to identify discrete longitudinal patterns of change in adolescent smoking using latent growth mixture modeling. Five distinct longitudinal patterns were identified. A group of early rapid escalators was characterized by early escalation (at age 13) that rapidly increased to heavy smoking. A pattern characterized by occasional puffing up until age 15, at which time smoking escalated to moderate levels was also identified (late moderate escalators). Another group included adolescents who, after age 15, began to escalate slowly in their smoking to light (0.5 cigarettes per month) levels (late slow escalators). Finally, a group of stable light smokers (those who smoked 1-2 cigarettes per month) and a group of stable puffers (those. who smoked only a few puffs per month) were also identified. The stable puffer group was the largest group and represented 25% of smokers.

Solid-state NMR structure determination of melittin in a lipid environment

Solid-state C-13 NMR spectroscopy was used to investigate the three-dimensional structure of melittin as lyophilized powder and in ditetradecylphosphatidylcholine (DTPC) membranes. The distance between specifically labeled carbons in analogs [1-C-13]Gly3-[2-C-13]Ala4, [1-C-13]Gly3-[2-C-13]Leu6, [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 was measured by rotational resonance. As expected, the internuclear distances measured in [1-C-13]Gly3-[2-C-13]Ala4 and [1-C-13]Gly3-[2-C-13]Leu6 were consistent with alpha -helical structure in the N-terminus irrespective of environment. The Internuclear distances measured in [1-C-13]Leu13-[2-C-13]Ala15, [2-C-13]Leu13-[1-C-13]Ala15, and [1-C-13]Leu13-[2-C-13]Leu16 revealed, via molecular modeling, some dependence upon environment for conformation in the region of the bend in helical structure induced by Pro14. A slightly larger interhelical angle between the N- and C-terminal helices was indicated for peptide in dry or hydrated gel state DTPC (139 degrees -145 degrees) than in lyophilized powder (121 degrees -139 degrees) or crystals (129 degrees). The angle, however, is not as great as deduced for melittin in aligned bilayers of DTPC in the liquid-crystalline state (similar to 160 degrees) (R. Smith, F. Separovic, T. J. Milne, A. Whittaker, F. M. Bennett, B. A. Cornell, and A. Makriyannis, 1994, J. Mol, Biol 241:456-466). The study illustrates the utility of rotational resonance in determining local structure within peptide-lipid complexes.

Enzymology, structure, and dynamics of acetohydroxy acid isomeroreductase

Acetohydroxy acid isomeroreductase is a key enzyme involved in the biosynthetic pathway of the amino acids isoleucine, valine, and leucine. This enzyme is of great interest in agrochemical research because it is present only in plants and microorganisms, making it a potential target for specific herbicides and fungicides. Moreover, it catalyzes an unusual two-step reaction that is of great fundamental interest. With a view to characterizing both the mechanism of inhibition by potential herbicides and the complex reaction mechanism, various techniques of enzymology, molecular biology, mass spectrometry, X-ray crystallography, and theoretical simulation have been used. The results and conclusions of these studies are described briefly in this paper.

A formal representation of assumptions in process modelling

In this work, we present a systematic approach to the representation of modelling assumptions. Modelling assumptions form the fundamental basis for the mathematical description of a process system. These assumptions can be translated into either additional mathematical relationships or constraints between model variables, equations, balance volumes or parameters. In order to analyse the effect of modelling assumptions in a formal, rigorous way, a syntax of modelling assumptions has been defined. The smallest indivisible syntactical element, the so called assumption atom has been identified as a triplet. With this syntax a modelling assumption can be described as an elementary assumption, i.e. an assumption consisting of only an assumption atom or a composite assumption consisting of a conjunction of elementary assumptions. The above syntax of modelling assumptions enables us to represent modelling assumptions as transformations acting on the set of model equations. The notion of syntactical correctness and semantical consistency of sets of modelling assumptions is defined and necessary conditions for checking them are given. These transformations can be used in several ways and their implications can be analysed by formal methods. The modelling assumptions define model hierarchies. That is, a series of model families each belonging to a particular equivalence class. These model equivalence classes can be related to primal assumptions regarding the definition of mass, energy and momentum balance volumes and to secondary and tiertinary assumptions regarding the presence or absence and the form of mechanisms within the system. Within equivalence classes, there are many model members, these being related to algebraic model transformations for the particular model. We show how these model hierarchies are driven by the underlying assumption structure and indicate some implications on system dynamics and complexity issues. (C) 2001 Elsevier Science Ltd. All rights reserved.

Using a cartographic modeling language to manipulate spectral satellite imagery

Land related information about the Earth's surface is commonIJ found in two forms: (1) map infornlation and (2) satellite image da ta. Satellite imagery provides a good visual picture of what is on the ground but complex image processing is required to interpret features in an image scene. Increasingly, methods are being sought to integrate the knowledge embodied in mop information into the interpretation task, or, alternatively, to bypass interpretation and perform biophysical modeling directly on derived data sources. A cartographic modeling language, as a generic map analysis package, is suggested as a means to integrate geographical knowledge and imagery in a process-oriented view of the Earth. Specialized cartographic models may be developed by users, which incorporate mapping information in performing land classification. In addition, a cartographic modeling language may be enhanced with operators suited to processing remotely sensed imagery. We demonstrate the usefulness of a cartographic modeling language for pre-processing satellite imagery, and define two nerv cartographic operators that evaluate image neighborhoods as post-processing operations to interpret thematic map values. The language and operators are demonstrated with an example image classification task.

Structural Dynamics in (ND4)2[Cu(D2O)6](SO4)2

The extended X-ray absorption fine structure spectroscopy (EXAFS) of (ND4)(2)[CU(D2O)(6)](SO4)(2) at 5, 14,100, 200, and 298 K is reported. This indicates that the Cu-O bond lengths of the Cu(D2O)(6)(2+) ion do not change significantly within this temperature range, which contrasts with EPR results and X-ray and neutron diffraction experiments, which imply that two of the Cu-(D2O) bonds converge in length as the temperature is raised. The EXAFS measurements thus confirm that the bond distances yielded by the diffraction experiments refer to the average positions of ligands involved in a dynamic equilibrium in which the directions of the long and intermediate bonds of the Jahn-Teller distorted Cu(D2O)(6)(2+) ion are interchanged in the crystal lattice. Analysis of the displacement parameters is consistent with this interpretation, as are the wave functions calculated using a model involving Jahn-Teller vibronic coupling and the influence of lattice strain interactions.

Stochastic gauges in quantum dynamics for many-body simulations

Quantum dynamics simulations can be improved using novel quasiprobability distributions based on non-orthogonal Hermitian kernel operators. This introduces arbitrary functions (gauges) into the stochastic equations. which can be used to tailor them for improved calculations. A possible application to full quantum dynamic simulations of BEC's is presented. (C) 2001 Elsevier Science B.V. All rights reserved.

Haptic information stabilizes and destabilizes coordination dynamics

Goal-directed, coordinated movements in humans emerge from a variety of constraints that range from 'high-level' cognitive strategies based oil perception of the task to 'low-level' neuromuscular-skeletal factors such as differential contributions to coordination from flexor and extensor muscles. There has been a tendency in the literature to dichotomize these sources of constraint, favouring one or the other rather than recognizing and understanding their mutual interplay. In this experiment, subjects were required to coordinate rhythmic flexion and extension movements with an auditory metronome, the rate of which was systematically increased. When subjects started in extension on the beat of the metronome, there was a small tendency to switch to flexion at higher rates, but not vice versa. When subjects: were asked to contact a physical stop, the location of which was either coincident with or counterphase to the auditor) stimulus, two effects occurred. When haptic contact was coincident with sound, coordination was stabilized for both flexion and extension. When haptic contact was counterphase to the metronome, coordination was actually destabilized, with transitions occurring from both extension to flexion on the beat and from flexion to extension on the beat. These results reveal the complementary nature of strategic and neuromuscular factors in sensorimotor coordination. They also suggest the presence of a multimodal neural integration process-which is parametrizable by rate and context - in which intentional movement, touch and sound are bound into a single, coherent unit.