Multi-physics interactions drive VEGFR2 relocation on endothelial cells

Figure 2

Valentina Damioli,1, Alberto Salvadori,2,3, Gian Paolo Beretta,1, Cosetta Ravelli,3,4, Stefania Mitola,3,4

1 Università degli Studi di Brescia, DIMI Department of Mechanical and Industrial Engineering, Brescia, 25123, Italy; 2 Università degli Studi di Brescia, DICATAM, Department of Civil, Environmental, Architectural Engineering and Mathematics, Brescia, 25123, Italy; 3 Laboratory for Preventive and Personalized Medicine (MPP Lab), Università degli Studi di Brescia, Brescia, 25123, Italy; 4 Università degli Studi di Brescia, DMMT, Department of Molecular and Translational Medicine, Brescia, 25123, Italy.

DOI: 10.1038/s41598-017-16786-4

Vascular Endothelial Growth Factor Receptor-2 (VEGFR2) is a pro-angiogenic receptor, expressed on endothelial cells (ECs). Although biochemical pathways that follow the VEGFR2 activation are well established, knowledge about the dynamics of receptors on the plasma membrane remains limited. Ligand stimulation induces the polarization of ECs and the relocation of VEGFR2, either in cell protrusions or in the basal aspect in cells plated on ligand-enriched extracellular matrix (ECM). We develop a mathematical model in order to simulate the relocation of VEGFR2 on the cell membrane during the mechanical adhesion of cells onto a ligand-enriched substrate. Co-designing the in vitro experiments with the simulations allows identifying three phases of the receptor dynamics, which are controlled respectively by the high chemical reaction rate, by the mechanical deformation rate, and by the diffusion of free receptors on the membrane. The identification of the laws that regulate receptor polarization opens new perspectives toward developing innovative anti-angiogenic strategies through the modulation of EC activation.

 

 

 

 

Ultrahigh field MR imaging of a subconjunctival anti-glaucoma drug delivery system in a rabbit model

Figure 2

Franziska Kopp,1, Thomas Eickner,2, Stefan Polei,3, Karen Falke,1, Martin Witt,4, Niels Grabow,2, Oliver Stachs,1, Rudolf F. Guthoff,1, Tobias Lindner,3

1 Department of Ophthalmology, Rostock University Medical Center, Rostock, Germany;  2 Institute for Biomedical Engineering, Rostock University Medical Center, Rostock, Germany; 3 Core Facility Multimodal Small Animal Imaging, Rostock University Medical Center, Rostock, Germany; 4 Institute for Anatomy, Rostock University Medical Center, Rostock, Germany.

DOI: 10.1038/s41598-017-15954-w

Local drug delivery systems (DDS) have become a favourable approach for the treatment of numerous diseases. Biomedical imaging techniques such as ultrahigh field magnetic resonance imaging (UHF-MRI) offer unique insight into DDS biodegradation in vivo. We describe the establishment of a 7 Tesla MRI routine for longitudinal in vivo examinations of a subconjunctival DDS for the treatment of glaucoma in a rabbit model. In initial in vitro examinations the T2-relaxation times of the polymeric DDS components were assessed. Imaging of enzymatically degraded depot samples in vitro did not reveal changes in sample morphology or T2-relaxation time. Ex vivo investigations with an enucleated porcine eye showed good correlation of anatomical MRI and histological data. In longitudinal in vivo studies in rabbits, we repeatedly scanned the depot in the same animal over the course of 5 months with an in-plane resolution of 130 µm at scan times of less than 30 minutes. The degradation was quantified using volumetric analysis showing a volume reduction of 82% between 3 and 21 weeks after depot implantation. We have thereby demonstrated the feasibility of our UHF-MRI protocol as a non-invasive imaging routine for qualitative and quantitative, longitudinal evaluation of biodegradable subconjunctival DDS.

 

A data-driven modeling approach to identify disease-specific multi-organ networks driving physiological dysregulation

Warren D. Anderson,1, Danielle DeCicco,2, James S. Schwaber,2, Rajanikanth Vadigepalli,2

1 Center for Public Health Genomics, Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, VA, USA; 2  Daniel Baugh Institute for Functional Genomics and Computational Biology, Department of Pathology, Anatomy, and Cell Biology, Sidney Kimmel Medical College, Thomas Jefferson University, Philadelphia, PA, USA.

DOI: 10.1371/journal.pcbi.1005627

Multiple physiological systems interact throughout the development of a complex disease. Knowledge of the dynamics and connectivity of interactions across physiological systems could facilitate the prevention or mitigation of organ damage underlying complex diseases, many of which are currently refractory to available therapeutics (e.g., hypertension). We studied the regulatory interactions operating within and across organs throughout disease development by integrating in vivo analysis of gene expression dynamics with a reverse engineering approach to infer data-driven dynamic network models of multi-organ gene regulatory influences. We obtained experimental data on the expression of 22 genes across five organs, over a time span that encompassed the development of autonomic nervous system dysfunction and hypertension. We pursued a unique approach for identification of continuous-time models that jointly described the dynamics and structure of multi-organ networks by estimating a sparse subset of ∼12,000 possible gene regulatory interactions. Our analyses revealed that an autonomic dysfunction-specific multi-organ sequence of gene expression activation patterns was associated with a distinct gene regulatory network. We analyzed the model structures for adaptation motifs, and identified disease-specific network motifs involving genes that exhibited aberrant temporal dynamics. Bioinformatic analyses identified disease-specific single nucleotide variants within or near transcription factor binding sites upstream of key genes implicated in maintaining physiological homeostasis. Our approach illustrates a novel framework for investigating the pathogenesis through model-based analysis of multi-organ system dynamics and network properties. Our results yielded novel candidate molecular targets driving the development of cardiovascular disease, metabolic syndrome, and immune dysfunction.

Glial responses to implanted electrodes in the brain

Fig. 2

Joseph W. Salatino,1,2, Kip A. Ludwig,3, Takashi D. Y. Kozai,4,5,6,7, Erin K. Purcell,1,2,8,9

1 Department of Biomedical Engineering, Michigan State University, East Lansing, MI, USA; 2 Institute for Quantitative Health Science and Engineering, Michigan State University, East Lansing, MI, USA; 3 Department of Neurologic Surgery, Mayo Clinic, Rochester, MN, USA; 4 Department of Bioengineering, University of Pittsburgh, Pittsburgh, PA, USA; 5 Center for the Neural Basis of Cognition, University of Pittsburgh, Pittsburgh, PA, USA; 6 McGowan Institute of Regenerative Medicine, University of Pittsburgh, Pittsburgh, PA, USA; 7 Neurotech Center, University of Pittsburgh Brain Institute, Pittsburgh, PA, USA; 8 Department of Electrical and Computer Engineering, Michigan State University, East Lansing, MI, USA, 9 Neuroscience Program, Michigan State University, East Lansing, MI, USA.

DOI: 10.1038/s41551-017-0154-1

The use of implants that can electrically stimulate or record electrophysiological or neurochemical activity in nervous tissue is rapidly expanding. Despite remarkable results in clinical studies and increasing market approvals, the mechanisms underlying the therapeutic effects of neuroprosthetic and neuromodulation devices, as well as their side effects and reasons for their failure, remain poorly understood. A major assumption has been that the signal-generating neurons are the only important target cells of neural-interface technologies. However, recent evidence indicates that the supporting glial cells remodel the structure and function of neuronal networks and are an effector of stimulation-based therapy. Here, we reframe the traditional view of glia as a passive barrier, and discuss their role as an active determinant of the outcomes of device implantation. We also discuss the implications that this has on the development of bioelectronic medical devices.

There are more connections between neurons in the human brain than there are stars in our galaxy1, and there are at least a dozen specific neuronal subtypes in the brain that are recognized as unique on the basis of their distinctive functional and morphological characteristics2,3. There is also growing recognition that non-neuronal supporting cells are more diverse and dynamic than previously appreciated, with distinct classes and subclasses of glia actively shaping the structure and function of neural circuitry4. Although such complexity is a likely requisite for the ability to internalize, integrate and respond to the continuous streams of information that the brain must process, it also makes the effective treatment of neurological disorders especially challenging. In recent years, the development and design of new implantable-device technologies to read-out and write-in electrical and chemical signals to and from the nervous system have created unprecedented opportunities to understand normal brain function and to ameliorate dysfunction resulting from disease or injury.

Although research and clinical applications of implanted electrode arrays continue to experience rapid growth, their usage has outpaced the clear understanding of the mechanisms underlying their benefits, side effects and modes of failure. Originally a precision academic-research tool to measure and modulate neural circuitry at subsecond and submillimetre resolution, implanted electrode arrays have increasingly been used in the clinic to treat an expanding array of medical conditions. Reports in the late 1980s and early 1990s demonstrated compelling preliminary clinical efficacy of deep brain stimulation (DBS) for tremor as a safer alternative to thalamotomy or pallidotomy in medically intractable Parkinson’s disease5. Although the mechanisms underlying its benefits remain the subject of debate6, DBS has since been approved by the US Food and Drug Administration for Parkinson’s disease, essential tremor, obsessive compulsive disorder, dystonia and refractory epilepsy5. Therapeutic indications currently being pursued in clinical studies are rapidly expanding, and include Alzheimer’s disease, depression, Tourette’s syndrome, deafness, blindness, and strategies to promote plasticity in cases of severe stroke or tinnitus6. Electrophysiological and neurochemical recordings have gained traction as a diagnostic tool, as an enabling technology for brain/machine interfaces in paralysis patients and as biomarkers to inform strategies for closed-loop stimulation devices7.

The successful use of chronically implanted neuroprostheses is predicated on the ability to reliably modulate or record signals from surrounding neurons over time (preferably, for many years). This is true for the broad range of clinical and research applications pursued, and for the variety of methods of read-out or write-in of neural activity employed (such as optical or electrical)8. However, problems arising from small signal amplitudes and from signal instability plague implanted recording arrays, limiting their long-term function9,10,11,12,13. Signal amplitudes typically shift on a daily basis14, compromising the likelihood that spike detection crosses the required threshold. This can, in turn, affect apparent firing rates, contributing to the non-stationarity that burdens the use of these signals for prosthetic control13. Studies across animal models often report progressive losses in signal detection in the weeks following implantation14,15. In recordings taken from human subjects, significant changes in unit amplitudes were observed on an intraday basis13. Many of these shifts seem to be related to device micromotion (based on simultaneous effects observed across electrode sites)13, but the vast majority were attributed to a physiological origin (85%). Likewise, in applications that stimulate the central nervous system, desensitization can occur following chronic microstimulation, and inexplicably large placebo effects can follow implantation of non-functional devices16,17. A variety of factors, both biological and non-biological, have been proposed to contribute to observations of instability in neural recordings and to the variable thresholds of neurostimulation10,18. Among these, suboptimal biocompatibility and suboptimal integration with surrounding tissue remains a significant limitation to reliably transfer information to and from the brain through implanted electrode arrays.

 

 

 

Deformable Cell Model of Tissue Growth

Computation 05 00045 g003b

Nikolai Bessonov 1 and Vitaly Volpert 2
1Institute of Problems of Mechanical Engineering, Russian Academy of Sciences, 199178 Saint Petersburg, Russia; 2Institut Camille Jordan, UMR 5208 CNRS, University Lyon 1, 69622 Villeurbanne, France.
This paper is devoted to modelling tissue growth with a deformable cell model. Each cell represents a polygon with particles located at its vertices. Stretching, bending and pressure forces act on particles and determine their displacement. Pressure-dependent cell proliferation is considered. Various patterns of growing tissue are observed. An application of the model to tissue regeneration is illustrated. Approximate analytical models of tissue growth are developed.

The dynamics of molecular evolution over 60,000 generations

pic1

Benjamin H. Good,1,2,3,4,5, Michael J. McDonald,1,2,6, Jeffrey E. Barrick,7,8, Richard E. Lenski,8,9, Michael M. Desai,1,2,3

1 Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts 02138, USA; 2 FAS Center for Systems Biology, Harvard University, Cambridge, Massachusetts 02138, USA; 3 Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA; 4 Department of Physics, University of California Berkeley, Berkeley, California 94720, USA; 5 Department of Bioengineering, University of California Berkeley, Berkeley, California 94720, USA; 6 Centre for Geometric Biology, School of Biological Sciences, Monash University, Clayton, Victoria 3800, Australia; 7 Department of Molecular Biosciences, The University of Texas, Austin, Texas 78712, USA; 8 BEACON Center for the Study of Evolution in Action, Michigan State University, East Lansing, Michigan 48824, USA; 9 Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, Michigan 48824, USA.

DOI: 10.1038/nature24287

The outcomes of evolution are determined by a stochastic dynamical process that governs how mutations arise and spread through a population. However, it is difficult to observe these dynamics directly over long periods and across entire genomes. Here we analyse the dynamics of molecular evolution in twelve experimental populations of Escherichia coli, using whole-genome metagenomic sequencing at five hundred-generation intervals through sixty thousand generations. Although the rate of fitness gain declines over time, molecular evolution is characterized by signatures of rapid adaptation throughout the duration of the experiment, with multiple beneficial variants simultaneously competing for dominance in each population. Interactions between ecological and evolutionary processes play an important role, as long-term quasi-stable coexistence arises spontaneously in most populations, and evolution continues within each clade. We also present evidence that the targets of natural selection change over time, as epistasis and historical contingency alter the strength of selection on different genes. Together, these results show that long-term adaptation to a constant environment can be a more complex and dynamic process than is often assumed.

 

Hydration number: crucial role in nuclear magnetic relaxivity of Gd(III) chelate-based nanoparticles

Figure 3

Rustem Zairov,1, Gulshat Khakimullina,2, Sergey Podyachev,1, Irek Nizameev,3, Georgy Safiullin,4, Rustem Amirov,2, Alberto Vomiero,5, Asiya Mustafin,1

1 A. E. Arbuzov Institute of Organic and Physical Chemistry, Kazan Scientific Center of Russian Academy of Sciences, Arbuzov str., 8, 420088, Kazan, Russia; 2 Kazan (Volga region) Federal university, Kremlyovskaya str., 18, 420008, Kazan, Russia; 3 Kazan National Research Technical University, K. Marx str., 10, 420211, Kazan, Russia; 4 Zavoisky Physical -Technical Institute of Kazan Scientific Center of Russian Academy of Sciences, Sibirsky tract, 10/7, 420029, Kazan, Russia; 5 Division of Materials Science, Department of Engineering Sciences and Mathematics, Luleå University of Technology, SE-971 87, Luleå, Sweden.

DOI: 10.1038/s41598-017-14409-6

Today, nanostructure-based contrast agents (CA) are emerging in the field of magnetic resonance imaging (MRI). Their sensitivity is reported as greatly improved in comparison to commercially used chelate-based ones. The present work is aimed at revealing the factors governing the efficiency of longitudinal magnetic relaxivity (r1) in aqueous colloids of core-shell Gd(III)-based nanoparticles. We report for the first time on hydration number (q) of gadolinium(III) as a substantial factor in controlling r1 values of polyelectrolyte-stabilized nanoparticles built from water insoluble complexes of Gd(III). The use of specific complex structure enables to reveal the impact of the inner-sphere hydration number on both r1 values for the Gd(III)-based nanoparticles and the photophysical properties of their luminescent Tb(III) and Eu(III) counterparts. The low hydration of TTA-based Gd(III) complexes (q ≈ 1) agrees well with the poor relaxivity values (r1 = 2.82 mM−1s−1 and r2 = 3.95 mM−1s−1), while these values tend to increase substantially (r1 = 12.41 mM−1s−1, r2 = 14.36 mM−1s−1) for aqueous Gd(III)-based colloids, when macrocyclic 1,3-diketonate is applied as the ligand (q ≈ 3). The regularities obtained in this work are fundamental in understanding the efficiency of MRI probes in the fast growing field of nanoparticulate contrast agents.

 

Partial volume correction for PET quantification and its impact on brain network in Alzheimer’s disease

Figure 1

Jiarui Yang,1,2, Chenhui Hu,2, Ning Guo,2, Joyita Dutta,2,3, Lucia M. Vaina,1,2, Keith A. Johnson,2,4, Jorge Sepulcre,2,4, Georges El Fakhri,2,4, Quanzheng Li,2,4

1 Boston University, Department of Biomedical Engineering, Boston, 02215, USA; 2 Massachusetts General Hospital, Department of Radiology, Boston, 02114, USA; 3 University of Massachusetts Lowell, Department of Electrical and Computer Engineering, Lowell, 01854, USA; 4 Harvard Medical School, Department of Radiology, Boston, 02115, USA.

DOI: 10.1038/s41598-017-13339-7

Amyloid positron emission tomography (PET) imaging is a valuable tool for research and diagnosis in Alzheimer’s disease (AD). Partial volume effects caused by the limited spatial resolution of PET scanners degrades the quantitative accuracy of PET image. In this study, we have applied a method to evaluate the impact of a joint-entropy based partial volume correction (PVC) technique on brain networks learned from a clinical dataset of AV-45 PET image and compare network properties of both uncorrected and corrected image-based brain networks. We also analyzed the region-wise SUVRs of both uncorrected and corrected images. We further performed classification tests on different groups using the same set of algorithms with same parameter settings. PVC has sometimes been avoided due to increased noise sensitivity in image registration and segmentation, however, our results indicate that appropriate PVC may enhance the brain network structure analysis for AD progression and improve classification performance.

 

High resolution 3-Dimensional imaging of the human cardiac conduction system from microanatomy to mathematical modeling

Figure 7

Robert S. Stephenson,1,2, Andrew Atkinson,3, Petros Kottas,4, Filip Perde,5, Fatemeh Jafarzadeh,3, Mike Bateman,6, Paul A. Iaizzo,6, Jichao Zhao,7, Henggui Zhang,4, Robert H. Anderson,8, Jonathan C. Jarvis,9, Halina Dobrzynski,3

1 Comparative Medicine Lab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; 2 School of Sport and Exercise Sciences, Liverpool John Moores University, Liverpool, UK; 3 Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK; 4 School of Physics and Astronomy, University of Manchester, Manchester, UK; 5 National Institute of Legal Medicine, Bucharest, Romania; 6 The Visible Heart Laboratory, University of Minnesota, Minneapolis, USA; 7 Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand; 8 Institute of Genetic Medicine, University of Newcastle, Newcastle, UK; 9 School of Sport and Exercise Sciences, Liverpool John Moores University, Liverpool, UK.

DOI: 10.1038/s41598-017-07694-8

Cardiac arrhythmias and conduction disturbances are accompanied by structural remodelling of the specialised cardiomyocytes known collectively as the cardiac conduction system. Here, using contrast enhanced micro-computed tomography, we present, in attitudinally appropriate fashion, the first 3-dimensional representations of the cardiac conduction system within the intact human heart. We show that cardiomyocyte orientation can be extracted from these datasets at spatial resolutions approaching the single cell. These data show that commonly accepted anatomical representations are oversimplified. We have incorporated the high-resolution anatomical data into mathematical simulations of cardiac electrical depolarisation. The data presented should have multidisciplinary impact. Since the rate of depolarisation is dictated by cardiac microstructure, and the precise orientation of the cardiomyocytes, our data should improve the fidelity of mathematical models. By showing the precise 3-dimensional relationships between the cardiac conduction system and surrounding structures, we provide new insights relevant to valvar replacement surgery and ablation therapies. We also offer a practical method for investigation of remodelling in disease, and thus, virtual pathology and archiving. Such data presented as 3D images or 3D printed models, will inform discussions between medical teams and their patients, and aid the education of medical and surgical trainees.

 

The role of backward cell migration in two-hit mutants’ production in the stem cell niche

Audrey Bollas,1, Leili Shahriyari,2

1 Department of Mathematics, The Ohio State University, Columbus, OH, United States of America; 2 Mathematical Biosciences Institute, The Ohio State University, Columbus, OH, United States of America.

DOI: 10.1371/journal.pone.0184651

It has been discovered that there are two stem cell groups in the intestinal crypts: central stem cells (CeSCs), which are at the very bottom of the crypt, and border stem cells (BSCs), which are located between CeSCs and transit amplifying cells (TAs). Moreover, backward cell migration from BSCs to CeSCs has been observed. Recently, a bi-compartmental stochastic model, which includes CeSCs and BSCs, has been developed to investigate the probability of two-hit mutant production in the stem cell niche. In this project, we improve this stochastic model by adding the probability of backward cell migration to the model. The model suggests that the probability of two-hit mutant production increases when the frequency of backward cell migration increases. Furthermore, a small non-zero probability of backward cell migration leads to the largest range of optimal values for the frequency of symmetric divisions and the portion of divisions at each stem cell compartment in terms of delaying 2-hit mutant production. Moreover, the probability of two-hit mutant production is more sensitive to the probability of symmetric divisions than to the rate of backward cell migrations. The highest probability of two-hit mutant production corresponds to the case when all stem cell’s divisions are asymmetric.