How amoeboids self-organize into a fruiting body: multicellular coordination in Dictyostelium discoideum.

When individual amoebae of the cellular slime mold Dictyostelium discoideum are starving, they aggregate to form a multicellular migrating slug, which moves toward a region suitable for culmination. The culmination of the morphogenesis involves complex cell movements that transform a mound of cells into a globule of spores on a slender stalk. The movement has been likened to a "reverse fountain," whereby prestalk cells in the upper part form a stalk that moves downwards and anchors to the substratum, while prespore cells in the lower part move upwards to form the spore head. So far, however, no satisfactory explanation has been produced for this process. Using a computer simulation that we developed, we now demonstrate that the processes that are essential during the earlier stages of the morphogenesis are in fact sufficient to produce the dynamics of the culmination stage. These processes are cAMP signaling, differential adhesion, cell differentiation, and production of extracellular matrix. Our model clarifies the processes that generate the observed cell movements. More specifically, we show that periodic upward movements, caused by chemotactic motion, are essential for successful culmination, because the pressure waves they induce squeeze the stalk downwards through the cell mass. The mechanisms revealed by our model have a number of self-organizing and self-correcting properties and can account for many previously unconnected and unexplained experimental observations.

Release of virus from lymphoid tissue affects human immunodeficiency virus type 1 and hepatitis C virus kinetics in the blood.

Kinetic parameters of human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV) infections have been estimated from plasma virus levels following perturbation of the chronically infected (quasi-) steady state. We extend previous models by also considering the large pool of virus localized in the lymphoid tissue (LT) compartment. The results indicate that the fastest time scale of HIV-1 plasma load decay during therapy probably reflects the clearance rate of LT virus and not, as previously supposed, the clearance rate of virus in plasma. This resolves the discrepancy between the clearance rate estimates during therapy and those based on plasma apheresis experiments. In the extended models plasma apheresis measurements are indeed expected to reflect the plasma decay rate. We can reconcile all current HIV-1 estimates with this model when, on average, the clearance rate of virus in plasma is approximately 20 day(-1), that of LT virus is approximately 3 day(-1), and the death rate of virus-producing cells is approximately 0.5 day(-1). The fast clearance in the LT compartment increases current estimates for total daily virus production. Because HCV is produced in the liver, we let virus be produced into the blood compartment of our model. The results suggest that extending current HCV models with an LT compartment is not likely to affect current estimates for kinetic parameters and virus production. Estimates for treatment efficacy might be affected, however.

Small variations in multiple parameters account for wide variations in HIV-1 set-points: a novel modelling approach.

Steady-state levels of HIV-1 viraemia in the plasma vary more than a 1,000-fold between HIV-positive patients and are thought to be influenced by several different host and viral factors such as host target cell availability, host anti-HIV immune response and the virulence of the virus. Previous mathematical models have taken the form of classical ecological food-chain models and are unable to account for this multifactorial nature of the disease. These models suggest that the steady-state viral load (i.e. the set-point) is determined by immune response parameters only. We have devised a generalized consensus model in which the conventional parameters are replaced by so-called 'process functions'. This very general approach yields results that are insensitive to the precise form of the mathematical model. Here we applied the approach to HIV-1 infections by estimating the steady-state values of several process functions from published patient data. Importantly, these estimates are generic because they are independent of the precise form of the underlying processes. We recorded the variation in the estimated steady-state values of the process functions in a group of HIV-1 patients. We developed a novel model by providing explicit expressions for the process functions having the highest patient-to-patient variation in their estimated values. Small variations from patient to patient for several parameters of the new model collectively accounted for the large variations observed in the steady-state viral burden. The novel model remains in full agreement with previous models and data.

Estimating relative fitness in viral competition experiments.

The relative fitness of viral variants has previously been defined as the slope of the logarithmic ratio of the genotype or phenotype frequencies in time plots of pairwise competition experiments. Developing mathematical models for such experiments by employing the conventional coefficient of selection s, we demonstrate that this logarithmic ratio gives the fitness difference, rather than the relative fitness. This fitness difference remains proportional to the actual replication rate realized in the particular experimental setup and hence cannot be extrapolated to other situations. Conversely, the conventional relative fitness (1 + s) should be more generic. We develop an approach to compute the generic relative fitness in conventional competition experiments. This involves an estimation of the total viral replication during the experiment and requires an estimate of the average lifetime of productively infected cells. The novel approach is illustrated by estimating the relative fitness, i.e., the relative replication rate, of a set of zidovudine-resistant human immunodeficiency virus type 1 variants. A tool for calculating the relative fitness from observed changes in viral load and genotype (or phenotype) frequencies is publically available on the website at http://www-binf.bio.uu.nl/( approximately )rdb/fitness.html.

Responses of complex cells in area 17 of the cat to bi-vectorial transparent motion.

We examined the responses to transparent motion of complex cells in cat area 17 which show directional selectivity to moving random pixel arrays (RPAs). The response to an RPA moving in the cell's preferred direction is inhibited when a second RPA is transparently moving in another direction. The inhibition by the second pattern is quantified as a function of its direction. The response to a pattern moving in the preferred direction is never completely suppressed, not even when a second pattern is moving transparently in the opposite direction. To the extent that supra-spontaneous firing rates signal the presence of the optimal velocity vector, these cells therefore still signal the presence of this line-label stimulus despite additional opposing, or otherwise directed, motion components. The results confirm previous suggestions that, for the computation of motion energy in cat area 17 complex cells, a full opponent stage is not plausible. Furthermore, we show that the response to a combination of two motion vectors can be predicted by the average of the responses to the individual components.