Control and controllability of microswimmers by a shearing flow

With the continuing rapid development of artificial microrobots and active particles, questions of microswimmer guidance and control are becoming ever more relevant and prevalent. In both the applications and theoretical study of such microscale swimmers, control is often mediated by an engineered property of the swimmer, such as in the case of magnetically propelled microrobots. In this work, we will consider a modality of control that is applicable in more generality, effecting guidance via modulation of a background fluid flow. Here, considering a model swimmer in a commonplace flow and simple geometry, we analyse and subsequently establish the efficacy of flow-mediated microswimmer positional control, later touching upon a question of optimal control. Moving beyond idealized notions of controllability and towards considerations of practical utility, we then evaluate the robustness of this control modality to sources of variation that may be present in applications, examining in particular the effects of measurement inaccuracy and rotational noise. This exploration gives rise to a number of cautionary observations, which, overall, demonstrate the need for the careful assessment of both policy and behavioural robustness when designing control schemes for use in practice

Turing Instabilities are Not Enough to Ensure Pattern Formation

Symmetry-breaking instabilities play an important role in understanding the mechanisms underlying the diversity of patterns observed in nature, such as in Turing’s reaction–diffusion theory, which connects cellular signalling and transport with the development of growth and form. Extensive literature focuses on the linear stability analysis of homogeneous equilibria in these systems, culminating in a set of conditions for transport-driven instabilities that are commonly presumed to initiate self-organisation. We demonstrate that a selection of simple, canonical transport models with only mild multistable non-linearities can satisfy the Turing instability conditions while also robustly exhibiting only transient patterns. Hence, a Turing-like instability is insufficient for the existence of a patterned state. While it is known that linear theory can fail to predict the formation of patterns, we demonstrate that such failures can appear robustly in systems with multiple stable homogeneous equilibria. Given that biological systems such as gene regulatory networks and spatially distributed ecosystems often exhibit a high degree of multistability and nonlinearity, this raises important questions of how to analyse prospective mechanisms for self-organisation

The effects of rapid yawing on simple swimmer models and planar Jeffery's orbits

Over a sufficiently long period of time, or from an appropriate distance, the motion of many swimmers can appear smooth, with their trajectories appearing almost ballistic in nature and slowly varying in character. These long-time behaviours, however, often mask more complex dynamics, such as the side-to-side snakelike motion exhibited by spermatozoa as they swim, propelled by the frequent and periodic beating of their flagellum. Many models of motion neglect these effects in favour of smoother long-term behaviours, which are often of greater practical interest than the small-scale oscillatory motion. Whilst it may be tempting to ignore any yawing motion, simply assuming that any effects of rapid oscillations cancel out over a period, a precise quantification of the impacts of high-frequency yawing is lacking. In this study, we systematically evaluate the long-term effects of general high-frequency oscillations on translational and angular motion, cast in the context of microswimmers but applicable more generally. Via a multiple-scales asymptotic analysis, we show that rapid oscillations can cause a long-term bias in the average direction of progression. We identify sufficient conditions for an unbiased long-term effect of yawing, and we quantify how yawing modifies the speed of propulsion and the effective hydrodynamic shape when in shear flow. Furthermore, we investigate and justify the long-time validity of the derived leading-order solutions and, by direct computational simulation, we evidence the relevance of the presented results to a canonical microswimmer

Mathematical and computational models of the retina in health, development and disease

The retina confers upon us the gift of vision, enabling us to perceive the world in a manner unparalleled by any other tissue. Experimental and clinical studies have provided great insight into the physiology and biochemistry of the retina; however, there are questions which cannot be answered using these methods alone. Mathematical and computational techniques can provide complementary insight into this inherently complex and nonlinear system. They allow us to characterise and predict the behaviour of the retina, as well as to test hypotheses which are experimentally intractable. In this review, we survey some of the key theoretical models of the retina in the healthy, developmental and diseased states. The main insights derived from each of these modelling studies are highlighted, as are model predictions which have yet to be tested, and data which need to be gathered to inform future modelling work. Possible directions for future research are also discussed. Whilst the present modelling studies have achieved great success in unravelling the workings of the retina, they have yet to achieve their full potential. For this to happen, greater involvement with the modelling community is required, and stronger collaborations forged between experimentalists, clinicians and theoreticians. It is hoped that, in addition to bringing the fruits of current modelling studies to the attention of the ophthalmological community, this review will encourage many such future collaborations

Cardiogenic shock following administration of propofol and fentanyl in a healthy woman: a case report

<p>Abstract</p> <p>Introduction</p> <p>Cardiogenic shock is very uncommon in healthy people. The differential diagnosis for patients with acute heart failure in previously healthy hearts includes acute myocardial infarction and myocarditis. However, many drugs can also depress myocardial function. Propofol and fentanyl are frequently used during different medical procedures. The cardiovascular depressive effect of both drugs has been well established, but the development of cardiogenic shock is very rare when these agents are used.</p> <p>Case presentation</p> <p>After a minor surgical intervention, a 32-year-old Caucasian woman with no significant medical history went into sudden hemodynamic deterioration due to acute heart failure. An urgent echocardiogram showed severe biventricular dysfunction and an estimated left ventricular ejection fraction of 20%. Extracorporeal life support and mechanical ventilation were required. Five days later her ventricular function had fully recovered, which allowed the progressive withdrawal of medical treatment. Prior to her hospital discharge, cardiac MRI showed neither edema nor pathological deposits on the delayed contrast enhancement sequences. At her six-month follow-up examination, the patient was asymptomatic and did not require treatment.</p> <p>Conclusion</p> <p>Although there are many causes of cardiogenic shock, the presence of abrupt hemodynamic deterioration and the absence of a clear cause could be related to the use of propofol and fentanyl.</p

Optimizing Combination Therapies with Existing and Future CML Drugs

Small-molecule inhibitors imatinib, dasatinib and nilotinib have been developed to treat Chromic Myeloid Leukemia (CML). The existence of a triple-cross-resistant mutation, T315I, has been a challenging problem, which can be overcome by finding new inhibitors. Many new compounds active against T315I mutants are now at different stages of development. In this paper we develop an algorithm which can weigh different combination treatment protocols according to their cross-resistance properties, and find the protocols with the highest probability of treatment success. This algorithm also takes into account drug toxicity by minimizing the number of drugs used, and their concentration. Although our methodology is based on a stochastic model of CML microevolution, the algorithm itself does not require measurements of any parameters (such as mutation rates, or division/death rates of cells), and can be used by medical professionals without a mathematical background. For illustration, we apply this algorithm to the mutation data obtained in [1], [2]

Endothelial Cell Capture of Heparin-Binding Growth Factors under Flow

Circulation is an important delivery method for both natural and synthetic molecules, but microenvironment interactions, regulated by endothelial cells and critical to the molecule's fate, are difficult to interpret using traditional approaches. In this work, we analyzed and predicted growth factor capture under flow using computer modeling and a three-dimensional experimental approach that includes pertinent circulation characteristics such as pulsatile flow, competing binding interactions, and limited bioavailability. An understanding of the controlling features of this process was desired. The experimental module consisted of a bioreactor with synthetic endothelial-lined hollow fibers under flow. The physical design of the system was incorporated into the model parameters. The heparin-binding growth factor fibroblast growth factor-2 (FGF-2) was used for both the experiments and simulations. Our computational model was composed of three parts: (1) media flow equations, (2) mass transport equations and (3) cell surface reaction equations. The model is based on the flow and reactions within a single hollow fiber and was scaled linearly by the total number of fibers for comparison with experimental results. Our model predicted, and experiments confirmed, that removal of heparan sulfate (HS) from the system would result in a dramatic loss of binding by heparin-binding proteins, but not by proteins that do not bind heparin. The model further predicted a significant loss of bound protein at flow rates only slightly higher than average capillary flow rates, corroborated experimentally, suggesting that the probability of capture in a single pass at high flow rates is extremely low. Several other key parameters were investigated with the coupling between receptors and proteoglycans shown to have a critical impact on successful capture. The combined system offers opportunities to examine circulation capture in a straightforward quantitative manner that should prove advantageous for biologicals or drug delivery investigations

A Comprehensive Analysis of Shared Loci between Systemic Lupus Erythematosus (SLE) and Sixteen Autoimmune Diseases Reveals Limited Genetic Overlap

In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non–Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10−06) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn's disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non–MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases

Tuning the Mammalian Circadian Clock: Robust Synergy of Two Loops

The circadian clock is accountable for the regulation of internal rhythms in most living organisms. It allows the anticipation of environmental changes during the day and a better adaptation of physiological processes. In mammals the main clock is located in the suprachiasmatic nucleus (SCN) and synchronizes secondary clocks throughout the body. Its molecular constituents form an intracellular network which dictates circadian time and regulates clock-controlled genes. These clock-controlled genes are involved in crucial biological processes including metabolism and cell cycle regulation. Its malfunction can lead to disruption of biological rhythms and cause severe damage to the organism. The detailed mechanisms that govern the circadian system are not yet completely understood. Mathematical models can be of great help to exploit the mechanism of the circadian circuitry. We built a mathematical model for the core clock system using available data on phases and amplitudes of clock components obtained from an extensive literature search. This model was used to answer complex questions for example: how does the degradation rate of Per affect the period of the system and what is the role of the ROR/Bmal/REV-ERB (RBR) loop? Our findings indicate that an increase in the RNA degradation rate of the clock gene Period (Per) can contribute to increase or decrease of the period - a consequence of a non-monotonic effect of Per transcript stability on the circadian period identified by our model. Furthermore, we provide theoretical evidence for a potential role of the RBR loop as an independent oscillator. We carried out overexpression experiments on members of the RBR loop which lead to loss of oscillations consistent with our predictions. These findings challenge the role of the RBR loop as a merely auxiliary loop and might change our view of the clock molecular circuitry and of the function of the nuclear receptors (REV-ERB and ROR) as a putative driving force of molecular oscillations

A Bioinformatics Filtering Strategy for Identifying Radiation Response Biomarker Candidates

The number of biomarker candidates is often much larger than the number of clinical patient data points available, which motivates the use of a rational candidate variable filtering methodology. The goal of this paper is to apply such a bioinformatics filtering process to isolate a modest number (<10) of key interacting genes and their associated single nucleotide polymorphisms involved in radiation response, and to ultimately serve as a basis for using clinical datasets to identify new biomarkers. In step 1, we surveyed the literature on genetic and protein correlates to radiation response, in vivo or in vitro, across cellular, animal, and human studies. In step 2, we analyzed two publicly available microarray datasets and identified genes in which mRNA expression changed in response to radiation. Combining results from Step 1 and Step 2, we identified 20 genes that were common to all three sources. As a final step, a curated database of protein interactions was used to generate the most statistically reliable protein interaction network among any subset of the 20 genes resulting from Steps 1 and 2, resulting in identification of a small, tightly interacting network with 7 out of 20 input genes. We further ranked the genes in terms of likely importance, based on their location within the network using a graph-based scoring function. The resulting core interacting network provides an attractive set of genes likely to be important to radiation response