Modelling the coupling of the M-clock and C-clock in lymphatic muscle cells.

Chronic dysfunction of the lymphatic vascular system results in fluid accumulation between cells: lymphoedema. The condition is commonly acquired secondary to diseases such as cancer or the associated therapies. The primary driving force for fluid return through the lymphatic vasculature is provided by contractions of the muscularized lymphatic collecting vessels, driven by electrochemical oscillations. However, there is an incomplete understanding of the molecular and bioelectric mechanisms involved in lymphatic muscle cell excitation, hampering the development and use of pharmacological therapies. Modelling in silico has contributed greatly to understanding the contributions of specific ion channels to the cardiac action potential, but modelling of these processes in lymphatic muscle remains limited. Here, we propose a model of oscillations in the membrane voltage (M-clock) and intracellular calcium concentrations (C-clock) of lymphatic muscle cells. We modify a model by Imtiaz and colleagues to enable the M-clock to drive the C-clock oscillations. This approach differs from typical models of calcium oscillators in lymphatic and related cell types, but is required to fit recent experimental data. We include an additional voltage dependence in the gating variable control for the L-type calcium channel, enabling the M-clock to oscillate independently of the C-clock. We use phase-plane analysis to show that these M-clock oscillations are qualitatively similar to those of a generalised FitzHugh-Nagumo model. We also provide phase plane analysis to understand the interaction of the M-clock and C-clock oscillations. The model and methods have the potential to help determine mechanisms and find targets for pharmacological treatment of lymphoedema.

Correction to: A computational model of a network of initial lymphatics and pre-collectors with permeable interstitium.

In the published paper, we argued that, although there appear to be no data available on the subject, it is inherently unlikely that lymph having traversed a network of initial lymphatics and pre-collectors then encounters an abrupt transition to vessels with all of the typical properties of collecting lymphatics.

A computational model of a network of initial lymphatics and pre-collectors with permeable interstitium.

Initial lymphatic vessels are made up of overlapped endothelial cells that act as unidirectional valves enabling one-way drainage of tissue fluid into the lumen of the initial lymphatics when there is a favourable pressure gradient. Initial lymphatics subsequently drain this fluid into the collecting lymphatics. This paper describes a computational model for a network of passive rat mesenteric lymphatic vessels with sparse secondary valves. The network was simulated with the secondary valves both operational and non-operational. The effects on the cycle-mean outflow-rate from the network of both inflammation and the resistance of the surrounding interstitium were considered. The cycle-mean outflow-rate is sensitive to vessel stiffness. If the influence of primary-valve resistance is reduced relative to that of interstitial resistance and intravascular resistance, there is no absolute advantage of extrinsic pumping, since maximum outflow-rate occurs when vessels are rigid. However, there is relative advantage, in that the outflow-rate at intermediate stiffness is higher with the secondary valves functioning than when they are deactivated. If primary-valve resistance dominates, then extrinsic pumping of non-rigid vessels provides absolute advantage. The nonlinear relation between pressure drop and flow-rate of the endothelial primary valves, combined with downstream compliance and pulsatile external pressure, constitutes a separate mechanism of pumping. By enabling the consideration of interactions between multiple phenomena (primary valves, secondary valves, a real network geometry with multiple branches, deformable vessel walls, interstitial resistance and external pressures), the model offers a perspective for delineating physiological phenomena that have not yet been fully linked to the biomechanics of fluid flow through initial lymphatic networks.

Inhibition of Contraction Strength and Frequency by Wall Shear Stress in a Single-Lymphangion Model.

The phasic contractions of collecting lymphatic vessels are reduced in strength and occur at diminished frequency when a favorable pressure difference and the resulting antegrade flow create large fluid shear stresses at the luminal surface. This paper describes a minimal phenomenological model of this mechanism that is applied to a previously validated numerical model of a phasically contracting lymphangion. The parameters of the inhibition model are quantitatively matched to observations in isolated segments of rat lymphatic vessel, first for mesenteric lymphatics then for thoracic duct, and outcomes from the numerical model are then qualitatively compared with recent observations in isolated segments of rat thoracic duct.

Contraction of collecting lymphatics: organization of pressure-dependent rate for multiple lymphangions.

The paper describes the extension of a previously developed model of pressure-dependent contraction rate to the case of multiple lymphangions. Mechanical factors are key modulators of active lymphatic pumping. As part of the evolution of our lumped-parameter model to match experimental findings, we have designed an algorithm whereby the time until the next contraction depends on lymphangion transmural pressure in the contraction just completed. The functional dependence of frequency on pressure is quantitatively matched to isobaric contraction experiments on isolated lymphatic segments. When each of several lymphangions is given this ability, a scheme for their coordination must be instituted to match the observed synchronization. Accordingly, and in line with an experiment on an isolated lymphatic vessel segment in which we measured contraction sequence and conduction delay, we took the fundamental principle to be that local timing can be overridden by signals to initiate contraction that start in adjacent lymphangions, conducted with a short delay. The scheme leads to retrograde conduction when the lymphangion chain is pumping against an adverse pressure difference, but antegrade conduction when contractions occur with no or a favourable pressure difference. Abolition of these conducted signals leads to chaotic variation of cycle-mean flow-rate from the chain, diastolic duration in each lymphangion, and inter-lymphangion delays. Chaotic rhythm is also seen under other circumstances. Because the model responds to increasing adverse pressure difference by increasing the repetition rate of contractions, it maintains time-average output flow-rate better than one with fixed repetition rate.

Valve-related modes of pump failure in collecting lymphatics: numerical and experimental investigation.

Lymph is transported along collecting lymphatic vessels by intrinsic and extrinsic pumping. The walls have muscle of a type intermediate between blood-vascular smooth muscle and myocardium; a contracting segment between two valves (a lymphangion) constitutes a pump. This intrinsic mechanism is investigated ex vivo in isolated, spontaneously contracting, perfused segments subjected to controlled external pressures. The reaction to varying afterload is probed by slowly ramping up the outlet pressure until pumping fails. Often the failure occurs when the contraction raises intra-lymphangion pressure insufficiently to overcome the outlet pressure, open the outlet valve and cause ejection, but many segments fail by other means, the mechanisms of which are not clear. We here elucidate those mechanisms by resort to a numerical model. Experimental observations are paired with comparable findings from computer simulations, using a lumped-parameter model that incorporates previously measured valve properties, plus new measurements of active contractile and passive elastic properties, and the dependence of contraction frequency on transmural pressure, all taken from isobaric twitch contraction experiments in the same vessel. Surprisingly, the model predicts seven different possible modes of pump failure, each defined by a different sequence of valve events, with their occurrence depending on the parameter values and boundary conditions. Some, but not all, modes were found experimentally. Further model investigation reveals routes by which a vessel exhibiting one mode of failure might under altered circumstances exhibit another.

Pump function curve shape for a model lymphatic vessel.

The transport capacity of a contractile segment of lymphatic vessel is defined by its pump function curve relating mean flow-rate and adverse pressure difference. Numerous system characteristics affect curve shape and the magnitude of the generated flow-rates and pressures. Some cannot be varied experimentally, but their separate and interacting effects can be systematically revealed numerically. This paper explores variations in the rate of change of active tension and the form of the relation between active tension and muscle length, factors not known from experiment to functional precision. Whether the pump function curve bends toward or away from the origin depends partly on the curvature of the passive pressure-diameter relation near zero transmural pressure, but rather more on the form of the relation between active tension and muscle length. A pump function curve bending away from the origin defines a well-performing pump by maximum steady output power. This behaviour is favoured by a length/active-tension relationship which sustains tension at smaller lengths. Such a relationship also favours high peak mechanical efficiency, defined as output power divided by the input power obtained from the lymphangion diameter changes and active-tension time-course. The results highlight the need to pin down experimentally the form of the length/active-tension relationship.

Incorporating measured valve properties into a numerical model of a lymphatic vessel.

An existing lumped-parameter model of multiple lymphangions (lymphatic vascular segments) in series is adapted for the incorporation of recent physiological measurements of lymphatic vascular properties. The new data show very marked nonlinearity of the passive pressure-diameter relation during distension, relative to comparable blood vessels, and complex valve behaviour. Since lymph is transported as a result of either the active contraction or the passive squeezing of vascular segments situated between two one-way valves, the performance of these valves is of primary importance. The valves display hysteresis (the opening and closing pressure drop thresholds differ), a bias to staying open (both state changes occur when the trans-valve pressure drop is adverse) and pressure-drop threshold dependence on transmural pressure. These properties, in combination with the strong nonlinearity that valve operation represents, have in turn caused intriguing numerical problems in the model, and we describe numerical stratagems by which we have overcome the problems. The principal problem is also generalised into a relatively simple mathematical example, for which solution detail is provided using two different solvers.

Development of a model of a multi-lymphangion lymphatic vessel incorporating realistic and measured parameter values.

Our published model of a lymphatic vessel consisting of multiple actively contracting segments between non-return valves has been further developed by the incorporation of properties derived from observations and measurements of rat mesenteric vessels. These included (1) a refractory period between contractions, (2) a highly nonlinear form for the passive part of the pressure-diameter relationship, (3) hysteretic and transmural-pressure-dependent valve opening and closing pressure thresholds and (4) dependence of active tension on muscle length as reflected in local diameter. Experimentally, lymphatic valves are known to be biased to stay open. In consequence, in the improved model, vessel pumping of fluid suffers losses by regurgitation, and valve closure is dependent on backflow first causing an adverse valve pressure drop sufficient to reach the closure threshold. The assumed resistance of an open valve therefore becomes a critical parameter, and experiments to measure this quantity are reported here. However, incorporating this parameter value, along with other parameter values based on existing measurements, led to ineffective pumping. It is argued that the published measurements of valve-closing pressure threshold overestimate this quantity owing to neglect of micro-pipette resistance. An estimate is made of the extent of the possible resulting error. Correcting by this amount, the pumping performance is improved, but still very inefficient unless the open-valve resistance is also increased beyond the measured level. Arguments are given as to why this is justified, and other areas where experimental data are lacking are identified. The model is capable of future adaptation as new experimental data appear.

Simulation of a chain of collapsible contracting lymphangions with progressive valve closure.

The aim of this investigation was to achieve the first step toward a comprehensive model of the lymphatic system. A numerical model has been constructed of a lymphatic vessel, consisting of a short series chain of contractile segments (lymphangions) and of intersegmental valves. The changing diameter of a segment governs the difference between the flows through inlet and outlet valves and is itself governed by a balance between transmural pressure and passive and active wall properties. The compliance of segments is maximal at intermediate diameters and decreases when the segments are subject to greatly positive or negative transmural pressure. Fluid flow is the result of time-varying active contraction causing diameter to reduce and is limited by segmental viscous and valvular resistance. The valves effect a smooth transition from low forward-flow resistance to high backflow resistance. Contraction occurs sequentially in successive lymphangions in the forward-flow direction. The behavior of chains of one to five lymphangions was investigated by means of pump function curves, with variation of valve opening parameters, maximum contractility, lymphangion size gradation, number of lymphangions, and phase delay between adjacent lymphangion contractions. The model was reasonably robust numerically, with mean flow-rate generally reducing as adverse pressure was increased. Sequential contraction was found to be much more efficient than synchronized contraction. At the highest adverse pressures, pumping failed by one of two mechanisms, depending on parameter settings: either mean leakback flow exceeded forward pumping or contraction failed to open the lymphangion outlet valve. Maximum pressure and maximum flow-rate were both sensitive to the contractile state; maximum pressure was also determined by the number of lymphangions in series. Maximum flow-rate was highly sensitive to the transmural pressure experienced by the most upstream lymphangions, suggesting that many feeding lymphatics would be needed to supply one downstream lymphangion chain pumping at optimal transmural pressure.

Unifying scaling theory for vortex dynamics in two-dimensional turbulence.

We present a scaling theory for unforced inviscid two-dimensional turbulence. Our model unifies existing spatial and temporal scaling theories. The theory is based on a self-similar distribution of vortices of different sizes A. Our model uniquely determines the spatial and temporal scaling of the associated vortex number density which allows the determination of the energy spectra and the vortex distributions. We find that the vortex number density scales as n(A,t)-t(-2/3)/A, which implies an energy spectrum E-k(-5), significantly steeper than the classical Batchelor-Kraichnan scaling. High-resolution numerical simulations corroborate the model.

Iterative approach for the numerical simulation of scattering from one- and two-dimensional rough surfaces.

We describe the use of iterative techniques for the solution of the integral equation that arises in an exact treatment of scalar wave scattering from randomly rough surfaces. The surfaces vary in either one or two dimensions, and the special case of a Dirichlet boundary condition is treated. It is found that these techniques, particularly when preconditioning is applied, are much more efficient than direct inversion techniques. Moreover, convergence is obtained for rms roughness of the order of 1, so the techniques have applicability over a wide parameter regime. Convergence is always to the exact solution found by direct inversion, exceptly for cases of extremely large-scaled rms surface heights in which the iterative techniques fail. In addition, by monitoring the residuals in the iteration process, it is immediately clear if the iterative techniques are failing, or performing badly in any given case. Finally, numerical results are compared with existing data in the enhanced backscattering regime.