Maurizio Tomaiuolo

A systems approach to hemostasis: 3. Thrombus consolidation regulates intrathrombus solute transport and local thrombin activity

Hemostatic thrombi formed after a penetrating injury have a distinctive structure in which a core of highly activated, closely packed platelets is covered by a shell of less-activated, loosely packed platelets. We have shown that differences in intrathrombus molecular transport emerge in parallel with regional differences in platelet packing density and predicted that these differences affect thrombus growth and stability. Here we test that prediction in a mouse vascular injury model. The studies use a novel method for measuring thrombus contraction in vivo and a previously characterized mouse line with a defect in integrin αIIbβ3 outside-in signaling that affects clot retraction ex vivo. The results show that the mutant mice have a defect in thrombus consolidation following vascular injury, resulting in an increase in intrathrombus transport rates and, as predicted by computational modeling, a decrease in thrombin activity and platelet activation in the thrombus core. Collectively, these data (1) demonstrate that in addition to the activation state of individual platelets, the physical properties of the accumulated mass of adherent platelets is critical in determining intrathrombus agonist distribution and platelet activation and (2) define a novel role for integrin signaling in the regulation of intrathrombus transport rates and localization of thrombin activity.

A systems approach to hemostasis: 1. The interdependence of thrombus architecture and agonist movements in the gaps between platelets.

Hemostatic thrombi develop a characteristic architecture in which a core of highly activated platelets is covered by a shell of less-activated platelets. Here we have used a systems biology approach to examine the interrelationship of this architecture with transport rates and agonist distribution in the gaps between platelets. Studies were performed in mice using probes for platelet accumulation, packing density, and activation plus recently developed transport and thrombin activity probes. The results show that intrathrombus transport within the core is much slower than within the shell. The region of slowest transport coincides with the region of greatest packing density and thrombin activity, and appears prior to full platelet activation. Deleting the contact-dependent signaling molecule, Sema4D, delays platelet activation, but not the emergence of the low transport region. Collectively, these results suggest a timeline in which initial platelet accumulation and the narrowing gaps between platelets create a region of reduced transport that facilitates local thrombin accumulation and greater platelet activation, whereas faster transport rates within the shell help to limit thrombin accumulation and growth of the core. Thus, from a systems perspective, platelet accumulation produces an altered microenvironment that shapes thrombus architecture, which in turn affects agonist distribution and subsequent thrombus growth.

A systems approach to hemostasis: 2. Computational analysis of molecular transport in the thrombus microenvironment

Hemostatic thrombi formed after a penetrating injury have a heterogeneous architecture in which a core of highly activated, densely packed platelets is covered by a shell of less-activated, loosely packed platelets. In the first manuscript in this series, we show that regional differences in intrathrombus protein transport rates emerge early in the hemostatic response and are preserved as the thrombus develops. Here, we use a theoretical approach to investigate this process and its impact on agonist distribution. The results suggest that hindered diffusion, rather than convection, is the dominant mechanism responsible for molecular movement within the thrombus. The analysis also suggests that the thrombus core, as compared with the shell, provides an environment for retaining soluble agonists such as thrombin, affecting the extent of platelet activation by establishing agonist-specific concentration gradients radiating from the site of injury. This analysis accounts for the observed weaker activation and relative instability of platelets in the shell and predicts that a failure to form a tightly packed thrombus core will limit thrombin accumulation, a prediction tested by analysis of data from mice with a defect in clot retraction.

Fast-activating voltage- and calcium-dependent potassium (BK) conductance promotes bursting in pituitary cells: a dynamic clamp study.

The electrical activity pattern of endocrine pituitary cells regulates their basal secretion level. Rat somatotrophs and lactotrophs exhibit spontaneous bursting and have high basal levels of hormone secretion, while gonadotrophs exhibit spontaneous spiking and have low basal hormone secretion. It has been proposed that the difference in electrical activity between bursting somatotrophs and spiking gonadotrophs is due to the presence of large conductance potassium (BK) channels on somatotrophs but not on gonadotrophs. This is one example where the role of an ion channel type may be clearly established. We demonstrate here that BK channels indeed promote bursting activity in pituitary cells. Blocking BK channels in bursting lacto-somatotroph GH4C1 cells changes their firing activity to spiking, while further adding an artificial BK conductance via dynamic clamp restores bursting. Importantly, this burst-promoting effect requires a relatively fast BK activation/deactivation, as predicted by computational models. We also show that adding a fast-activating BK conductance to spiking gonadotrophs converts the activity of these cells to bursting. Together, our results suggest that differences in BK channel expression may underlie the differences in electrical activity and basal hormone secretion levels among pituitary cell types and that the rapid rate of BK channel activation is key to its role in burst promotion.

A THEORETICAL INVESTIGATION OF SYMPATRIC EVOLUTION OF TEMPORAL REPRODUCTIVE ISOLATION AS ILLUSTRATED BY MARINE BROADCAST SPAWNERS

Abstract Recent theory suggests that frequency-dependent disruptive selection in combination with assortative mating can lead to the establishment of reproductive isolation in sympatry. Here we explore how temporal variation in reproduction might simultaneously generate both disruptive selection and assortative mating, and result in sympatric speciation. The conceptual framework of the model may be applicable to biological systems with negative frequency-dependent selection, such as marine broadcast spawners or systems with pollinator limitation. We present a model that is motivated by recent findings in marine broadcast spawners and is parameterized with data from the Montastraea annularis species complex. Broadcast spawners reproduce via external fertilization and synchronous spawning is required to increase the probability of successful fertilization, but empirical evidence shows that as density increases, so does the risk of polyspermy. Polyspermy is the fusion of multiple sperm with an egg at fertilization, a process that makes the embryo unviable. Synchrony can therefore also act as a source of negative density-dependent disruptive selection. Model analysis shows that the interaction between polyspermy and spawning synchrony can lead to temporal reproductive isolation in sympatry and that, more generally, increased density promotes maintenance of genetic variation.

Modeling how reproductive ecology can drive protein diversification and result in linkage disequilibrium between sperm and egg proteins.

Gamete‐recognition proteins determine whether sperm and eggs are compatible at fertilization, and they often evolve rapidly. The source of selection driving the evolution of these proteins is still debated. It has been suggested that sexual conflict can result in proliferation of genetic variation and possibly linkage disequilibrium between sperm and egg proteins. Empirical evidence suggests that both male and female reproductive success can be predicted by their sperm ligand genotype, but why female success can be predicted by a protein expressed only in males is unknown. Here we use mathematical modeling to investigate the interaction between reproductive behavior and sperm availability on the evolution of sperm ligands and egg receptors. We consider haploid and diploid expression in gametes in two possible ecological scenarios, monogamous spawning and competitive spawning. Reproductive behavior plays an important role in determining possible outcomes resulting from sexual conflict. Sperm limitation selects for common genotypes regardless of mating behavior. Under conditions of sperm abundance, competitive spawning provides conditions for the persistence of allelic variation and gametic disequilibrium. With monogamous spawning, such conditions are more restrictive.

Harnessing the Platelet Signaling Network to Produce an Optimal Hemostatic Response

Once released into the circulation by megakaryocytes, circulating platelets can undergo rapid activation at sites of vascular injury and resist unwarranted activation, which can lead to heart attacks and strokes. Historically, the signaling mechanisms underlying the regulation of platelet activation have been approached as a collection of individual pathways unique to agonist. This review takes a different approach, casting platelet activation as the product of a signaling network, in which activating and restraining mechanisms interact in a flexible network that regulates platelet adhesiveness, cohesion between platelets, granule secretion, and the formation of a stable hemostatic thrombus.

Two Types of Burst Firing in Gonadotrophin-Releasing Hormone Neurones

Gonadotrophin‐releasing hormone (GnRH) neurones fire spontaneous bursts of action potentials, although little is understood about the underlying mechanisms. In the present study, we report evidence for two types of bursting/oscillation driven by different mechanisms. Properties of these different types are clarified using mathematical modelling and a recently developed active‐phase/silent‐phase correlation technique. The first type of GnRH neurone (1–2%) exhibits slow (∼0.05 Hz) spontaneous oscillations in membrane potential. Action potential bursts are often observed during oscillation depolarisation, although some oscillations were entirely subthreshold. Oscillations persist after blockade of fast sodium channels with tetrodotoxin (TTX) and blocking receptors for ionotropic fast synaptic transmission, indicating that they are intrinsically generated. In the second type of GnRH neurone, bursts were irregular and TTX caused a stable membrane potential. The two types of bursting cells exhibited distinct active‐phase/silent‐phase correlation patterns, which is suggestive of distinct mechanisms underlying the rhythms. Further studies of type 1 oscillating cells revealed that the oscillation period was not affected by current or voltage steps, although amplitude was sometimes damped. Oestradiol, an important feedback regulator of GnRH neuronal activity, acutely and markedly altered oscillations, specifically depolarising the oscillation nadir and initiating or increasing firing. Blocking calcium‐activated potassium channels, which are rapidly reduced by oestradiol, had a similar effect on oscillations. Kisspeptin, a potent activator of GnRH neurones, translated the oscillation to more depolarised potentials, without altering period or amplitude. These data show that there are at least two distinct types of GnRH neurone bursting patterns with different underlying mechanisms.

Models of Electrical Activity: Calibration and Prediction Testing on the Same Cell

Mathematical models are increasingly important in biology, and testability is becoming a critical issue. One limitation is that one model simulation tests a parameter set representing one instance of the biological counterpart, whereas biological systems are heterogeneous in their properties and behavior, and a model often is fitted to represent an ideal average. This is also true for models of a cells electrical activity; even within a narrowly defined population there can be considerable variation in electrophysiological phenotype. Here, we describe a computational experimental approach for parameterizing a model of the electrical activity of a cell in real time. We combine the inexpensive parallel computational power of a programmable graphics processing unit with the flexibility of the dynamic clamp method. The approach involves 1), recording a cells electrical activity, 2), parameterizing a model to the recording, 3), generating predictions, and 4), testing the predictions on the same cell used for the calibration. We demonstrate the experimental feasibility of our approach using a cell line (GH4C1). These cells are electrically active, and they display tonic spiking or bursting. We use our approach to predict parameter changes that can convert one pattern to the other.

Investigating Heterogeneity of Intracellular Calcium Dynamics in Anterior Pituitary Lactotrophs Using a Combined Modelling ⁄Experimental Approach

Cell responses are commonly heterogeneous, even within a subpopulation. In the present study, we investigate the source of heterogeneity in the Ca2+ response of anterior pituitary lactotrophs to a Ca2+ mobilisation agonist, thyrotrophin‐releasing hormone. This response is characterised by a sharp increase of cytosolic Ca2+ concentration as a result of mobilisation of Ca2+ from intracellular stores, followed by a decrease to an elevated plateau level that results from Ca2+ influx. We focus on heterogeneity of the evoked Ca2+ spike under extracellular Ca2+ free conditions. We introduce a method that uses the information provided by a mathematical model to characterise the source of heterogeneity. This method compares scatter plots of features of the Ca2+ response obtained experimentally with those made from the mathematical model. The model scatter plots reflect random variation of parameters over different ranges, and matching the experimental and model scatter plots allows us to predict which parameters are most variable. We find that a large degree of variation in Ca2+ efflux is a likely key contributor to the heterogeneity of Ca2+ responses to thyrotrophin‐releasing hormone in lactotrophs. This technique is applicable to any situation in which the heterogeneous biological response is described by a mathematical model.