Identifying optimal lipid raft characteristics required to promote nanoscale protein-protein interactions on the plasma membrane

The dynamic lateral segregation of signaling proteins into microdomains is proposed to facilitate signal transduction, but the constraints on microdomain size, mobility, and diffusion that might realize this function are undefined. Here we interrogate a stochastic spatial model of the plasma membrane to determine how microdomains affect protein dynamics. Taking lipid rafts as representative microdomains, we show that reduced protein mobility in rafts segregates dynamically partitioning proteins, but the equilibrium concentration is largely independent of raft size and mobility. Rafts weakly impede small-scale protein diffusion but more strongly impede long-range protein mobility. The long-range mobility of raft-partitioning and raft-excluded proteins, however, is reduced to a similar extent. Dynamic partitioning into rafts increases specific interprotein collision rates, but to maximize this critical, biologically relevant function, rafts must be small (diameter, 6 to 14 nm) and mobile. Intermolecular collisions can also be favored by the selective capture and exclusion of proteins by rafts, although this mechanism is generally less efficient than simple dynamic partitioning. Generalizing these results, we conclude that microdomains can readily operate as protein concentrators or isolators but there appear to be significant constraints on size and mobility if microdomains are also required to function as reaction chambers that facilitate nanoscale protein-protein interactions. These results may have significant implications for the many signaling cascades that are scaffolded or assembled in plasma membrane microdomains.

Exploring the parameter space of a rabbit ventricular action potential model to investigate the effect of variation on action potential and calcium transients

Computational models for cardiomyocyte action potentials (AP) often make use of a large parameter set. This parameter set can contain some elements that are fitted to experimental data independently of any other element, some elements that are derived concurrently with other elements to match experimental data, and some elements that are derived purely from phenomenological fitting to produce the desired AP output. Furthermore, models can make use of several different data sets, not always derived for the same conditions or even the same species. It is consequently uncertain whether the parameter set for a given model is physiologically accurate. Furthermore, it is only recently that the possibility of degeneracy in parameter values in producing a given simulation output has started to be addressed. In this study, we examine the effects of varying two parameters (the L-type calcium current (I(CaL)) and the delayed rectifier potassium current (I(Ks))) in a computational model of a rabbit ventricular cardiomyocyte AP on both the membrane potential (V(m)) and calcium (Ca(2+)) transient. It will subsequently be determined if there is degeneracy in this model to these parameter values, which will have important implications on the stability of these models to cell-to-cell parameter variation, and also whether the current methodology for generating parameter values is flawed. The accuracy of AP duration (APD) as an indicator of AP shape will also be assessed.

Parameter estimation for a phenomenological model of the cardiac action potential

The action potential (ap) of a cardiac cell is made up of a complex balance of ionic currents which flow across the cell membrane in response to electrical excitation of the cell. Biophysically detailed mathematical models of the ap have grown larger in terms of the variables and parameters required to model new findings in subcellular ionic mechanisms. The fitting of parameters to such models has seen a large degree of parameter and module re-use from earlier models. An alternative method for modelling electrically exciteable cardiac tissue is a phenomenological model, which reconstructs tissue level ap wave behaviour without subcellular details. A new parameter estimation technique to fit the morphology of the ap in a four variable phenomenological model is presented. An approximation of a nonlinear ordinary differential equation model is established that corresponds to the given phenomenological model of the cardiac ap. The parameter estimation problem is converted into a minimisation problem for the unknown parameters. A modified hybrid Nelder–Mead simplex search and particle swarm optimization is then used to solve the minimisation problem for the unknown parameters. The successful fitting of data generated from a well known biophysically detailed model is demonstrated. A successful fit to an experimental ap recording that contains both noise and experimental artefacts is also produced. The parameter estimation method’s ability to fit a complex morphology to a model with substantially more parameters than previously used is established.

Increased sedation requirements in patients receiving extracorporeal membrane oxygenation for respiratory and cardiorespiratory failure

Critically ill patients receiving extracorporeal membrane oxygenation (ECMO) are often noted to have increased sedation requirements. However, data related to sedation in this complex group of patients is limited. The aim of our study was to characterise the sedation requirements in adult patients receiving ECMO for cardiorespiratory failure. A retrospective chart review was performed to collect sedation data for 30 consecutive patients who received venovenous or venoarterial ECMO between April 2009 and March 2011. To test for a difference in doses over time we used a regression model. The dose of midazolam received on ECMO support increased by an average of 18 mg per day (95% confidence interval 8, 29 mg, P=0.001), while the dose of morphine increased by 29 mg per day (95% confidence interval 4, 53 mg, P=0.021) The venovenous group received a daily midazolam dose that was 157 mg higher than the venoarterial group (95% confidence interval 53, 261 mg, P=0.005). We did not observe any significant increase in fentanyl doses over time (95% confidence interval 1269, 4337 µg, P=0.94). There is a significant increase in dose requirement for morphine and midazolam during ECMO. Patients on venovenous ECMO received higher sedative doses as compared to patients on venoarterial ECMO. Future research should focus on mechanisms behind these changes and also identify drugs that are most suitable for sedation during ECMO.

An in-vitro 3-D cell culture model for studying pathomechanisms in AMD

Purpose To develop a novel 3-D cell culture model with the view to studying the pathomechanisms underlying the development of age-related macular degeneration (AMD). Our central hypothesis is that the silk structural protein fibroin used in conjunction with cultured human cells can be used to mimic the structural relationships between the RPE and choriocapillaris in health and disease. Methods Co-cultures of human RPE cells (ARPE-19 cells grown in Miller’s medium) and microvascular endothelial cells (HMEC-1 cells grown in endothelial culture medium) were established on opposing sides of a synthetic Bruch’s membrane (3 microns thick) constructed from B mori silk fibroin. Cell attachment was facilitated by pre-coating the fibroin membrane with vitronectin (for ARPE-19 cells) and gelatin (for HMEC-1 cells) respectively. The effects of tropoelastin on attachment of ARPE-19 cells was also examined. Barrier function was examined by measurement of trans-epithelial resistance (TER) using a voltohmmeter (EVOM-2). The phagocytic activity of the synthetic RPE was tested using vitronectin-coated microspheres (2 micron diameter FluoSpheres). In some cultures, membrane defects were created by puncturing within a 24 G needle. The architecture of the synthetic tissue before and after wounding was examined by confocal microscopy after staining for ZO-1 and F-actin. Results The RPE layer of the 3D model developed a cobblestoned morphology (validated by staining for ZO-1 and F-actin), displayed barrier function (validated by measurement of TER) and demonstrated cytoplasmic uptake of vitronectin-coated microspheres. Attachment of ARPE-19 cells to fibroin was unaffected by tropoelastin. Microvascular endothelial cells attached well to the gelatin-coated surface of the fibroin membrane and remained physically separated from the overlaying RPE layer. The fibroin membranes were amenable to puncturing without collapse thus providing the opportunity to study transmembrane migration of the endothelial cells. Conclusions Synthetic Bruch’s membranes constructed from silk fibroin, vitronectin and gelatin, support the co-cultivation of RPE cells and microvascular endothelial cells. The resulting RPE layer displays functions similar to that of native RPE and the entire tri-layered structure displays potential to be used as an in vitro model of choroidal neovascularization.

Molecular and cellular analysis of basement membrane invasion by human breast cancer cells in Matrigel-based in vitro assays

In vitro analyses of basement membrane invasiveness employing Matrigel (a murine tumor extract rich in basement membrane components) have been performed on human breast cancer model systems. Constitutive invasiveness of different human breast cancer (HBC) cell lines has been examined as well as regulation by steroid hormones, growth factors, and oncogenes. Carcinoma cells exhibiting a mesenchymal-like phenotype (vimentin expression, lack of cell border associated uvomorulin) show dramatically increased motility, invasiveness, and metastatic potential in nude mice. These findings support the hypothesis that epithelial to mesenchymal transition (EMT)-like events may be instrumental in the metastatic progression of human breast cancer. The MCF-7 subline MCF-7ADR appears to have undergone such a transition. The importance of such a transition may be reflected in the emergence of vimentin expression as an indicator of poor prognosis in HBC. Matrix degradation and laminin recognition are highlighted as potential targets for antimetastatic therapy, and analyses of laminin attachment and the matrix metalloproteinase (MMP) family in HBC cell lines are summarized. Matrigel-based assays have proved useful in the study of the molecular mechanisms of basement membrane invasiveness, their regulation in HBC cells, and their potential as targets for antimetastatic therapy.

Activation of matrix metalloproteinase-2 (MMP-2) by membrane type 1 matrix metalloproteinase through an artificial receptor for ProMMP-2 generates active MMP-2

The suggested model for pro-matrix metalloproteinase-2 (proMMP-2) activation by membrane type 1 MMP (MT1-MMP) implicates the complex between MT1-MMP and tissue inhibitor of MMP-2 (TIMP-2) as a receptor for proMMP-2. To dissect this model and assess the pathologic significance of MMP-2 activation, an artificial receptor for proMMP-2 was created by replacing the signal sequence of TIMP-2 with cytoplasmic/transmembrane domain of type II transmembrane mosaic serine protease (MSP-T2). Unlike TIMP-2, MSP-T2 served as a receptor for proMMP-2 without inhibiting MT1-MMP, and generated TIMP-2-free active MMP-2 even at a low level of MT1-MMP. Thus, MSP-T2 did not affect direct cleavage of the substrate testican-1 by MT1-MMP, whereas TIMP-2 inhibited it even at the level that stimulates proMMP-2 processing. Expression of MSP-T2 in HT1080 cells enhanced MMP-2 activation by endogenous MT1-MMP and caused intensive hydrolysis of collagen gel. Expression of MSP-T2 in U87 glioma cells, which express a trace level of endogenous MT1-MMP, induced MMP-2 activation and enhanced cell-associated protease activity, activation of extracellular signal-regulated kinase, and metastatic ability into chick embryonic liver and lung. MT1-MMP can exert both maximum MMP-2 activation and direct cleavage of substrates with MSP-T2, which cannot be achieved with TIMP-2. These results suggest that MMP-2 activation by MT1-MMP potentially amplifies protease activity, and combination with direct cleavage of substrate causes effective tissue degradation and enhances tumor invasion and metastasis, which highlights the complex role of TIMP-2. MSP-T2 is a unique tool to analyze physiologic and pathologic roles of MMP-2 and MT1-MMP in comparison with TIMP-2.

Upregulated MT1-MMP/TIMP-2 axis in the TSU-Pr1-B1/B2 model of metastatic progression in transitional cell carcinoma of the bladder

Muscle invasive transitional cell carcinoma (TCC) of the bladder is associated with a high frequency of metastasis, resulting in poor prognosis for patients presenting with this disease. Models that capture and demonstrate step-wise enhancement of elements of the human metastatic cascade on a similar genetic background are useful research tools. We have utilized the transitional cell carcinoma cell line TSU-Pr1 to develop an in vivo experimental model of bladder TCC metastasis. TSU-Pr1 cells were inoculated into the left cardiac ventricle of SCID mice and the development of bone metastases was monitored using high resolution X-ray. Tumor tissue from a single bone lesion was excised and cultured in vitro to generate the TSU-Pr1-B1 subline. This cycle was repeated with the TSU-Pr1-B1 cells to generate the successive subline TSU-Pr1-B2. DNA profiling and karyotype analysis confirmed the genetic relationship of these three cell lines. In vitro, the growth rate of these cell lines was not significantly different. However, following intracardiac inoculation TSU-Pr1, TSU-Pr1-B1 and TSU-Pr1-B2 exhibited increasing metastatic potential with a concomitant decrease in time to the onset of radiologically detectable metastatic bone lesions. Significant elevations in the levels of mRNA expression of the matrix metalloproteases (MMPs) membrane type 1-MMP (MT1-MMP), MT2-MMP and MMP-9, and their inhibitor, tissue inhibitor of metalloprotease-2 (TIMP-2), across the progressively metastatic cell lines, were detected by quantitative PCR. Given the role of MT1-MMP and TIMP-2 in MMP-2 activation, and the upregulation of MMP-9, these data suggest an important role for matrix remodeling, particularly basement membrane, in this progression. The TSU-Pr1-B1/B2 model holds promise for further identification of important molecules.

Myogel, a novel, basement membrane-rich, extracellular matrix derived from skeletal muscle, is highly adipogenic in vivo and in vitro

Background/Aims Biological and synthetic scaffolds play important roles in tissue engineering and are being developed towards human clinical applications. Based on previous work from our laboratory, we propose that extracellular matrices from skeletal muscle could be developed for adipose tissue engineering. Methods Extracellular matrices (Myogels) extracted from skeletal muscle of various species were assessed using biochemical assays including ELISA and Western blotting. Biofunctionality was assessed using an in vitro differentiation assay and a tissue engineering construct model in the rat. Results Myogels were successfully extracted from mice, rats, pigs and humans. Myogels contained significant levels of laminin α4- and α2-subunits and collagen I compared to Matrigel™, which contains laminin 1 (α1β1γ1) and collagen IV. Levels of growth factors such as fibroblast growth factor 2 were significantly higher than Matrigel, vascular endothelial growth factor-A levels were significantly lower and all other growth factors were comparable. Myogels reproducibly stimulated adipogenic differentiation of preadipocytes in vitro and the growth of adipose tissue in the rat. Conclusions We found Myogel induces adipocyte differentiation in vitroand shows strong adipogenic potential in vivo, inducing the growth of well-vascularised adipose tissue. Myogel offers an alternative for current support scaffolds in adipose tissue engineering, allowing the scaling up of animal models towards clinical adipose tissue engineering applications.

Assessment of angiogenic potential--the use of AIDS-KS cell supernatants as an in vitro model

Metastasis, the passage of primary tumour cells throughout the body via the vascular system and their subsequent proliferation into secondary lesions in distant organs, represents a poor prognosis and therefore an understandably feared event for cancer patients. Despite considerable advances in cancer diagnosis and treatment, most deaths are the result of metastases resistant to conventional treatment [1]. Rather than being a random process, metastasis involves a series of organised steps leading to the growth of a secondary tumour. Malignant tumours stimulate the production of new vessels by the host, and this process is a prerequisite for the increase in size of a new tumour [2]. Angiogenesis, not only permits tumour expansion but also allows the entry of tumour cells into the circulation and is probably the most vital event for the metastatic process [3]. Metastasis and angiogenesis [4] have received much attention in recent years. A biological understanding of both phenomena seems to be an urgent priority towards the search for an effective prevention and treatment of tumour progression. Studies in vitro and in vivo have shown that one of the most important barriers to the passage of malignant cells is the basement membrane. The crossing of such barriers is a vital step in the formation of a metastasis [5].

A coupled SPH-DEM model for micro-scale structural deformations of plant cells during drying

A single plant cell was modeled with smoothed particle hydrodynamics (SPH) and a discrete element method (DEM) to study the basic micromechanics that govern the cellular structural deformations during drying. This two-dimensional particle-based model consists of two components: a cell fluid model and a cell wall model. The cell fluid was approximated to a highly viscous Newtonian fluid and modeled with SPH. The cell wall was treated as a stiff semi-permeable solid membrane with visco-elastic properties and modeled as a neo-Hookean solid material using a DEM. Compared to existing meshfree particle-based plant cell models, we have specifically introduced cell wall–fluid attraction forces and cell wall bending stiffness effects to address the critical shrinkage characteristics of the plant cells during drying. Also, a moisture domain-based novel approach was used to simulate drying mechanisms within the particle scheme. The model performance was found to be mainly influenced by the particle resolution, initial gap between the outermost fluid particles and wall particles and number of particles in the SPH influence domain. A higher order smoothing kernel was used with adaptive smoothing length to improve the stability and accuracy of the model. Cell deformations at different states of cell dryness were qualitatively and quantitatively compared with microscopic experimental findings on apple cells and a fairly good agreement was observed with some exceptions. The wall–fluid attraction forces and cell wall bending stiffness were found to be significantly improving the model predictions. A detailed sensitivity analysis was also done to further investigate the influence of wall–fluid attraction forces, cell wall bending stiffness, cell wall stiffness and the particle resolution. This novel meshfree based modeling approach is highly applicable for cellular level deformation studies of plant food materials during drying, which characterize large deformations.

Three-dimensional analysis of budding sites and released virus suggests a revised model for HIV-1 morphogenesis

Current models of HIV-1 morphogenesis hold that newly synthesized viral Gag polyproteins traffic to and assemble at the cell membrane into spherical protein shells. The resulting late-budding structure is thought to be released by the cellular ESCRT machinery severing the membrane tether connecting it to the producer cell. Using electron tomography and scanning transmission electron microscopy, we find that virions have a morphology and composition distinct from late-budding sites. Gag is arranged as a continuous but incomplete sphere in the released virion. In contrast, late-budding sites lacking functional ESCRT exhibited a nearly closed Gag sphere. The results lead us to propose that budding is initiated by Gag assembly, but is completed in an ESCRT-dependent manner before the Gag sphere is complete. This suggests that ESCRT functions early in HIV-1 release-akin to its role in vesicle formation-and is not restricted to severing the thin membrane tether.

A numerical method for the fractional Fitzhugh–Nagumo monodomain model

A fractional FitzHugh–Nagumo monodomain model with zero Dirichlet boundary conditions is presented, generalising the standard monodomain model that describes the propagation of the electrical potential in heterogeneous cardiac tissue. The model consists of a coupled fractional Riesz space nonlinear reaction-diffusion model and a system of ordinary differential equations, describing the ionic fluxes as a function of the membrane potential. We solve this model by decoupling the space-fractional partial differential equation and the system of ordinary differential equations at each time step. Thus, this means treating the fractional Riesz space nonlinear reaction-diffusion model as if the nonlinear source term is only locally Lipschitz. The fractional Riesz space nonlinear reaction-diffusion model is solved using an implicit numerical method with the shifted Grunwald–Letnikov approximation, and the stability and convergence are discussed in detail in the context of the local Lipschitz property. Some numerical examples are given to show the consistency of our computational approach.

Formation of the three-dimensional geometry of the red blood cell membrane

Red blood cells (RBCs) are nonnucleated liquid capsules, enclosed in deformable viscoelastic membranes with complex three dimensional geometrical structures. Generally, RBC membranes are highly incompressible and resistant to areal changes. However, RBC membranes show a planar shear deformation and out of plane bending deformation. The behaviour of RBCs in blood vessels is investigated using numerical models. All the characteristics of RBC membranes should be addressed to develop a more accurate and stable model. This article presents an effective methodology to model the three dimensional geometry of the RBC membrane with the aid of commercial software COMSOL Multiphysics 4.2a and Fortran programming. Initially, a mesh is generated for a sphere using the COMSOL Multiphysics software to represent the RBC membrane. The elastic energy of the membrane is considered to determine a stable membrane shape. Then, the actual biconcave shape of the membrane is obtained based on the principle of virtual work, when the total energy is minimised. The geometry of the RBC membrane could be used with meshfree particle methods to simulate motion and deformation of RBCs in micro-capillaries

A semi-alternating direction method for a 2-D fractional FitzHugh–Nagumo monodomain model on an approximate irregular domain

A FitzHugh-Nagumo monodomain model has been used to describe the propagation of the electrical potential in heterogeneous cardiac tissue. In this paper, we consider a two-dimensional fractional FitzHugh-Nagumo monodomain model on an irregular domain. The model consists of a coupled Riesz space fractional nonlinear reaction-diffusion model and an ordinary differential equation, describing the ionic fluxes as a function of the membrane potential. Secondly, we use a decoupling technique and focus on solving the Riesz space fractional nonlinear reaction-diffusion model. A novel spatially second-order accurate semi-implicit alternating direction method (SIADM) for this model on an approximate irregular domain is proposed. Thirdly, stability and convergence of the SIADM are proved. Finally, some numerical examples are given to support our theoretical analysis and these numerical techniques are employed to simulate a two-dimensional fractional Fitzhugh-Nagumo model on both an approximate circular and an approximate irregular domain.