Biomedical subjects

P K Maini

Publications and source records attributed to P K Maini.

At least 37 records · Page 2Linked to original sources

PubMed2000

A model for colour pattern formation in the butterfly wing of Papilio dardanus.

The butterfly Papilio dardanus is well known for the spectacular phenotypic polymorphism in the female of the species. We show that numerical simulations of a reaction diffusion model on a geometrically accurate wing domain produce spatial patterns that are consistent with many of those observed on the butterfly. Our results suggest that the wing coloration is due to a simple underlying stripe-like pattern of some pigment-inducing morphogen. We focus on the effect of key factors such as parameter values for mode selection, threshold values which determine colour, wing shape and boundary conditions. The generality of our approach should allow us to investigate other butterfly species. The relationship between these key factors and gene activities is discussed in the context of recent biological advances.

Animals↗

PubMed2000

Clock and induction model for somitogenesis.

After many years of research, somitogenesis is still one of the major unresolved problems in developmental biology. Recent experimental findings show a novel type of pattern formation in which a signal sweeps along the presomitic mesoderm and narrows simultaneously as a new somite is formed. The signal then residues in the posterior half of the new somite, and another wave begins to sweep up from the caudal end. This behaviour is not easily explained by the existing theoretical models. We present a new model for somitogenesis that can account for this behaviour and is consistent with previous experimental observations.

Animals↗

PubMed2000

Enzyme kinetics at high enzyme concentration.

We re-visit previous analyses of the classical Michaelis-Menten substrate-enzyme reaction and, with the aid of the reverse quasi-steady-state assumption, we challenge the approximation d[C]/dt approximately 0 for the basic enzyme reaction at high enzyme concentration. For the first time, an approximate solution for the concentrations of the reactants uniformly valid in time is reported. Numerical simulations are presented to verify this solution. We show that an analytical approximation can be found for the reactants for each initial condition using the appropriate quasi-steady-state assumption. An advantage of the present formalism is that it provides a new procedure for fitting experimental data to determine reaction constants. Finally, a new necessary criterion is found that ensures the validity of the reverse quasi-steady-state assumption. This is verified numerically.

Enzymes↗

PubMed2000

A chemotactic model for the advance and retreat of the primitive streak in avian development.

The formation of the primitive streak in early avian development marks the onset of gastrulation, during which large scale cell movement leads to a trilaminar blastoderm comprising prospective endodermal, mesodermal and ectodermal tissue. During streak formation a specialized group of cells first moves anteriorly as a coherent column, beginning from the posterior end of the prospective anterior-posterior axis (a process called progression), and then reverses course and returns to the most posterior point on the axis (a process called regression). To date little is known concerning the mechanisms controlling either progression or regression. Here we develop a model in which chemotaxis directs the cell movement and which is capable of reproducing the principal features connected with progression and regression of the primitive streak. We show that this model exhibits a number of experimentally-observed features of normal and abnormal streak development, and we propose a number of experimental tests which may serve to illuminate the mechanisms. This paper represents the first attempt to model the global features of primitive streak formation, and provides an initial stage in the development of a more biologically-realistic discrete cell model that will allow for variation of properties between cells and control over movement of individual cells.

Animals↗

PubMed2000

Development and applications of a model for cellular response to multiple chemotactic cues.

The chemotactic response of a cell population to a single chemical species has been characterized experimentally for many cell types and has been extensively studied from a theoretical standpoint. However, cells frequently have multiple receptor types and can detect and respond chemotactically to more than one chemical. How these signals are integrated within the cell is not known. and we therefore adopt a macroscopic phenomenological approach to this problem. In this paper we derive and analyze chemotactic models based on partial differential (chemotaxis) equations for cell movement in response to multiple chemotactic cues. Our derivation generalizes the approach of Othmer and Stevens [29], who have recently developed a modeling framework for studying different chemotactic responses to a single chemical species. The importance of such a generalization is illustrated by the effect of multiple chemical cues on the chemotactic sensitivity and the spatial pattern of cell densities in several examples. We demonstrate that the model can generate the complex patterns observed on the skin of certain animal species and we indicate how the chemotactic response can be viewed as a form of positional indicator.

Animals↗

PubMed2000

Turing instabilities in general systems.

We present necessary and sufficient conditions on the stability matrix of a general n(> or = 2)-dimensional reaction-diffusion system which guarantee that its uniform steady state can undergo a Turing bifurcation. The necessary (kinetic) condition, requiring that the system be composed of an unstable (or activator) and a stable (or inhibitor) subsystem, and the sufficient condition of sufficiently rapid inhibitor diffusion relative to the activator subsystem are established in three theorems which form the core of our results. Given the possibility that the unstable (activator) subsystem involves several species (dimensions), we present a classification of the analytically deduced Turing bifurcations into p (1 < or = p < or = (n - 1)) different classes. For n = 3 dimensions we illustrate numerically that two types of steady Turing pattern arise in one spatial dimension in a generic reaction-diffusion system. The results confirm the validity of an earlier conjecture [12] and they also characterise the class of so-called strongly stable, matrices for which only necessary conditions have been known before [23, 24]. One of the main consequences of the present work is that biological morphogens, which have so far been expected to be single chemical species [1-9], may instead be composed of two or more interacting species forming an unstable subsystem.

Diffusion↗

PubMed1999

Mathematical modelling of extracellular matrix dynamics using discrete cells: fiber orientation and tissue regeneration.

Matrix orientation plays a crucial role in determining the severity of scar tissue after dermal wounding. We present a model framework which allows us to examine the interaction of many of the factors involved in orientation and alignment. Within this framework, cells are considered as discrete objects, while the matrix is modelled as a continuum. Using numerical simulations, we investigate the effect on alignment of changing cell properties and of varying cell interactions with collagen and fibrin.

Animals↗

PubMed1999

Stripe formation in juvenile Pomacanthus explained by a generalized turing mechanism with chemotaxis.

Current interest in pattern formation can be traced to a seminal paper by Turing, who demonstrated that a system of reacting and diffusing chemicals, called morphogens, can interact so as to produce stable nonuniform concentration patterns in space. Recently, a Turing model has been suggested to explain the development of pigmentation patterns on species of growing angelfish such as Pomacanthus semicirculatus, which exhibit readily observed changes in the number, size, and orientation of colored stripes during development of juvenile and adult stages, but the model fails to predict key features of the observations on stripe formation. Here we develop a generalized Turing model incorporating cell growth and movement, we analyze the effects of these processes on patterning, and we demonstrate that the model can explain important features of pattern formation in a growing system such as Pomacanthus. The applicability of classical Turing models to biological pattern formation is limited by virtue of the sensitivity of patterns to model parameters, but here we show that the incorporation of growth results in robustly generated patterns without strict parameter control. In the model, chemotaxis in response to gradients in a morphogen distribution leads to aggregation of one type of pigment cell into a striped spatial pattern.

Animals↗

PubMed1999

The mathematical modelling of cell kinetics in corneal epithelial wound healing.

This paper considers the comparison of experimental spatial and temporal data of mitotic rates measured during corneal epithelial wound healing (CEWH) of a rat model with the predictions of a computer modelling framework. We begin by briefly showing that previous models, used in the study of corneal epithelial wound healing speeds, are inadequate for the study of cell kinetics. We proceed to formulate a new modelling framework more suited to such a study. This framework is simulated in its simplest form, and the results from this motivate a new realisation of the modelling framework, including a caricature of age structuring. Finally, a model with a simple representation of juxtacrine signalling is considered. The final model captures many, though not all, of the trends of the experimental data. This paper thus lays a foundation for the modelling of the cell kinetics of corneal epithelial wound healing, and yields valuable insight regarding the important mechanisms a model should consider in order to reproduce the observed experimental trends.

Animals↗

PubMed1999

Mathematical modelling of anisotropy in fibrous connective tissue.

We present two modelling frameworks for studying dynamic anistropy in connective tissue, motivated by the problem of fibre alignment in wound healing. The first model is a system of partial differential equations operating on a macroscopic scale. We show that a model consisting of a single extracellular matrix material aligned by fibroblasts via flux and stress exhibits behaviour that is incompatible with experimental observations. We extend the model to two matrix types and show that the results of this extended model are robust and consistent with experiment. The second model represents cells as discrete objects in a continuum of ECM. We show that this model predicts patterns of alignment on macroscopic length scales that are lost in a continuum model of the cell population.

Anisotropy↗

PubMed1999

Modelling corneal epithelial wound closure in the presence of physiological electric fields via a moving boundary formalism.

A new framework for the modelling of corneal epithelial wound healing is presented, which can include the presence of a physiological electric field. The difficulty inherent in the inclusion of this biological phenomenon motivates our use of a moving boundary formalism. A key conclusion is that the model predicts a linear relation between the wound healing speed and the physiological electric field strengths over a physiologically large range of electric field strength. Another key point is that this linear relationship between electric field strength and wound healing speed is robust to variations in critical parameters that are difficult to estimate. The linearity is also robust to different realizations of the modelling framework presented.

Animals↗

PubMed1998

Spatially varying equilibria of mechanical models: application to dermal wound contraction.

Mechanochemical models based on the Oster-Murray continuum framework have been applied to a variety of biological settings to obtain an understanding of the morphogenesis of living tissues. Wound-healing in mammalian skin is an important example, because a complex sequence of biochemical and biomechanical responses are orchestrated to close a wound by a combination of new tissue formation and wound contraction. Mechanical interactions between dermal fibroblastic cells and the collagen-rich extracellular matrix are crucial in the development of a contracted wound state. We and others have previously proposed mechanochemical models for wound repair to gain a greater understanding of both normal and abnormal healing. In the present work, the existence of spatially varying equilibria of these models is investigated by using a small-stain approximation and phase-plane techniques, with numerical simulations to confirm the analytical predictions. These results are sources of novel insight into the roles of key biological parameters in determining the mechanical properties of a contracted wound. These methods may also be relevant to other morphogenetic scenarios for which similar mechanochemical models have been proposed.

Animals↗

PubMed1998

Pattern formation in a generalized chemotactic model.

Many models have been proposed for spatial pattern formation in embryology and analyzed for the standard case of zero-flux boundary conditions. However, relatively little attention has been paid to the role of boundary conditions on the form of the final pattern. Here we investigate numerically, the effect of nonstandard boundary conditions on a model pattern generator, which we choose to be of a cell-chemotactic type. We specifically focus on the role of boundary conditions and the effects of scale and aspect ratio, and study the spatiotemporal dynamics of pattern formation. We illustrate the properties of the model by application to the spatiotemporal sequence of skeletal development.

Animals↗

PubMed1997

Hierarchically coupled ultradian oscillators generating robust circadian rhythms.

Ensembles of mutually coupled ultradian cellular oscillators have been proposed by a number of authors to explain the generation of circadian rhythms in mammals. Most mathematical models using many coupled oscillators predict that the output period should vary as the square root of the number of participating units, thus being inconsistent with the well-established experimental result that ablation of substantial parts of the suprachiasmatic nuclei (SCN), the main circadian pacemaker in mammals, does not eliminate the overt circadian functions, which show no changes in the phases or periods of the rhythms. From these observations, we have developed a theoretical model that exhibits the robustness of the circadian clock to changes in the number of cells in the SCN, and that is readily adaptable to include the successful features of other known models of circadian regulation, such as the phase response curves and light resetting of the phase.

Activity Cycles↗

PubMed1997

Role of fibroblast migration in collagen fiber formation during fetal and adult dermal wound healing.

Adult dermal wounds, in contrast to fetal wounds, heal with the formation of scar tissue. A crucial factor in determining the degree of scarring is the ratio of types I and III collagen, which regulates the diameter of the combined fibers. We developed a reaction-diffusion model which focuses on the control of collagen synthesis by different isoforms of the polypeptide transforming growth factor-beta (TGF beta). We used the model to investigate the current controversy as to whether the fibroblasts migrate into the wound from the surrounding unwounded dermis or from the underlying subcutaneous tissue. Numerical simulations of a spatially independent, temporal model led to a value of the collagen ratio consistent with that of healthy tissue for the fetus, but corresponding to scarring in the adult. We investigated the effect of topical application of TGF beta and show that addition of isoform 3 reduces scar tissue formation, in agreement with the experiment. However, numerical solutions of the reaction-diffusion system do not exhibit this sensitivity to growth factor application. Mathematically, this corresponds to the observation that behind healing wave-front solutions, a particular healed state is always selected independent of transients, even though there is a continuum of possible positive steady states. We explain this phenomenon using a caricature system of equations, which reflects the key qualitative features of the full model but has a much simpler mathematical form. Biologically, our results suggest that the migration into a wound of fibroblasts and TGF beta from the surrounding dermis alone cannot account for the essential features of the healing process, and that fibroblasts entering from the underlying subcutaneous tissue are crucial to the healing process.

Animals↗

PubMed1997

Bones, feathers, teeth and coat markings: a unified model.

A first necessary criterion to be met by any model for pattern formation is that it must be able to reproduce the patterns it purports to model. The above examples show how simple ideas from self-organisation can produce spatial patterns of varying complexity that are consistent with those observed experimentally. Second, the model must be consistent with the results of experimental manipulation. Third, the model must make experimentally testable predictions. In this way, a mathematical model can help to elucidate the underlying biochemical and biophysical mechanisms of pattern formation. I have tried to illustrate these ideas with the above examples. These examples also show the breadth of patterning phenomena that can be captured by these models. In animal coat markings, the patterns are laid down simultaneously. In limb development, the skeletal elements are laid down sequentially. The application to feather germ formation illustrates complex sequential regular pattern formation, while the application to tooth primordium formation illustrates sequential irregular pattern formation. Many other patterning phenomena have been studied (see, for example, ref. 7). Pattern formation is one of the central issues in developmental biology and intense interdisciplinary research involving experimentalists and theoreticians is beginning to help us understand how this phenomenon occurs. A detailed understanding of normal development is a necessary first step to the understanding of abnormal development and, hopefully, will help medical science combat developmental defects.

Alligators and Crocodiles↗

PubMed1997

A mathematical model for the capillary endothelial cell-extracellular matrix interactions in wound-healing angiogenesis.

Angiogenesis, the process by which new blood capillaries grow into a tissue from surrounding parent vessels, is a key event in dermal wound healing, malignant-tumour growth, and other pathologic conditions. In wound healing, new capillaries deliver vital metabolites such as amino acids and oxygen to the cells in the wound which are involved in a complex sequence of repair processes. The key cellular constituents of these new capillaries are endothelial cells: their interactions with soluble biochemical and insoluble extracellular matrix (ECM) proteins have been well documented recently, although the biological mechanisms underlying wound-healing angiogenesis are incompletely understood. Considerable recent research, including some continuum mathematical models, have focused on the interactions between endothelial cells and soluble regulators (such as growth factors). In this work, a similar modelling framework is used to investigate the roles of the insoluble ECM substrate, of which collagen is the predominant macromolecular protein. Our model consists of a partial differential equation for the endothelial-cell density (as a function of position and time) coupled to an ordinary differential equation for the ECM density. The ECM is assumed to regulate cell movement (both random and directed) and proliferation, whereas the cells synthesize and degrade the ECM. Analysis and numerical solutions of these equations highlights the roles of these processes in wound-healing angiogenesis. A nonstandard approximation analysis yields insight into the travelling-wave structure of the system. The model is extended to two spatial dimensions (parallel and perpendicular to the plane of the skin), for which numerical simulations are presented. The model predicts that ECM-mediated random motility and cell proliferation are key processes which drive angiogenesis and that the details of the functional dependence of these processes on the ECM density, together with the rate of ECM remodelling, determine the qualitative nature of the angiogenic response. These predictions are experimentally testable, and they may lead towards a greater understanding of the biological mechanisms involved in wound-healing angiogenesis.

Animals↗

PubMed1996

Pattern formation by lateral inhibition with feedback: a mathematical model of delta-notch intercellular signalling.

In many developing tissues, adjacent cells diverge in character so as to create a fine-grained pattern of cells in contrasting states of differentiation. It has been proposed that such patterns can be generated through lateral inhibition--a type of cell-cell interaction whereby a cell that adopts a particular fate inhibits its immediate neighbors from doing likewise. Lateral inhibition is well documented in flies, worms and vertebrates. In all of these organisms, the transmembrane proteins Notch and Delta (or their homologues) have been identified as mediators of the interaction--Notch as receptor, Delta as its ligand on adjacent cells. However, it is not clear under precisely what conditions the Delta-Notch mechanism of lateral inhibition can generate the observed types of pattern, or indeed whether this mechanism is capable of generating such patterns by itself. Here we construct and analyse a simple and general mathematical model of such contact-mediated lateral inhibition. In accordance with experimental data, the model postulates that receipt of inhibition (i.e. activation of Notch) diminished the ability to deliver inhibition (i.e. to produce active Delta). This gives rise to a feedback loop that can amplify differences between adjacent cells. We investigate the pattern-forming potential and temporal behaviour of this model both analytically and through numerical simulation. Inhomogeneities are self-amplifying and develop without need of any other machinery, provided the feedback is sufficiently strong. For a wide range of initial and boundary conditions, the model generates fine-grained patterns similar to those observed in living systems.

Animals↗