-
A network aggregation model for the dynamics and treatment of neurodegenerative diseases at the brain scale
Authors:
Georgia S. Brennan,
Alain Goriely
Abstract:
Neurodegenerative diseases are associated with the assembly of specific proteins into oligomers and fibrillar aggregates. At the brain scale, these protein assemblies can diffuse through the brain and seed other regions, creating an autocatalytic protein progression. The growth and transport of these assemblies depend on various mechanisms that can be targeted therapeutically. Here, we use spatial…
▽ More
Neurodegenerative diseases are associated with the assembly of specific proteins into oligomers and fibrillar aggregates. At the brain scale, these protein assemblies can diffuse through the brain and seed other regions, creating an autocatalytic protein progression. The growth and transport of these assemblies depend on various mechanisms that can be targeted therapeutically. Here, we use spatially-extended nucleation-aggregation-fragmentation models for the dynamics of prion-like neurodegenerative protein-spreading in the brain to study the effect of different drugs on whole-brain Alzheimer's disease progression.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
Soft cells and the geometry of seashells
Authors:
Gábor Domokos,
Alain Goriely,
Ákos G. Horváth,
Krisztina Regős
Abstract:
A central problem of geometry is the tiling of space with simple structures. The classical solutions, such as triangles, squares, and hexagons in the plane and cubes and other polyhedra in three-dimensional space are built with sharp corners and flat faces. However, many tilings in Nature are characterized by shapes with curved edges, non-flat faces, and few, if any, sharp corners. An important qu…
▽ More
A central problem of geometry is the tiling of space with simple structures. The classical solutions, such as triangles, squares, and hexagons in the plane and cubes and other polyhedra in three-dimensional space are built with sharp corners and flat faces. However, many tilings in Nature are characterized by shapes with curved edges, non-flat faces, and few, if any, sharp corners. An important question is then to relate prototypical sharp tilings to softer natural shapes. Here, we solve this problem by introducing a new class of shapes, the \textit{soft cells}, minimizing the number of sharp corners and filling space as \emph{soft tilings}. We prove that an infinite class of polyhedral tilings can be smoothly deformed into soft tilings and we construct the soft versions of all Dirichlet-Voronoi cells associated with point lattices in two and three dimensions. Remarkably, these ideal soft shapes, born out of geometry, are found abundantly in nature, from cells to shells.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Stability of topological descriptors for neuronal morphology
Authors:
David Beers,
Heather A. Harrington,
Alain Goriely
Abstract:
The topological morphology descriptor of a neuron is a multiset of intervals associated to the shape of the neuron represented as a tree. In practice, topological morphology descriptors are vectorized using persistence images, which can help classify and characterize the morphology of broad groups of neurons. We study the stability of topological morphology descriptors under small changes to neuro…
▽ More
The topological morphology descriptor of a neuron is a multiset of intervals associated to the shape of the neuron represented as a tree. In practice, topological morphology descriptors are vectorized using persistence images, which can help classify and characterize the morphology of broad groups of neurons. We study the stability of topological morphology descriptors under small changes to neuronal morphology. We show that the persistence diagram arising from the topological morphology descriptor of a neuron is stable for the 1-Wasserstein distance against a range of perturbations to the tree. These results guarantee that persistence images of topological morphology descriptors are stable against the same set of perturbations and reliable.
△ Less
Submitted 16 November, 2022;
originally announced November 2022.
-
Brain Chains as Topological Signatures for Alzheimer's Disease
Authors:
Christian Goodbrake,
David Beers,
Travis B. Thompson,
Heather A. Harrington,
Alain Goriely
Abstract:
We propose a topological framework to study the evolution of Alzheimer's disease, the most common neurodegenerative disease. The modeling of this disease starts with the representation of the brain connectivity as a graph and the seeding of a toxic protein in a specific region represented by a vertex. Over time, the accumulation of toxic proteins at vertices and their propagation along edges are m…
▽ More
We propose a topological framework to study the evolution of Alzheimer's disease, the most common neurodegenerative disease. The modeling of this disease starts with the representation of the brain connectivity as a graph and the seeding of a toxic protein in a specific region represented by a vertex. Over time, the accumulation of toxic proteins at vertices and their propagation along edges are modeled by a dynamical system on this graph. These dynamics provide an order on the edges of the graph according to the damage created by high concentrations of proteins. This sequence of edges defines a filtration of the graph. We consider different filtrations given by different disease seeding locations. To study this filtration we propose a new combinatorial and topological method. A filtration defines a maximal chain in the partially ordered set of spanning subgraphs ordered by inclusion. To identify similar graphs, and define a topological signature, we quotient this poset by graph homotopy equivalence, which gives maximal chains in a smaller poset. We provide an algorithm to compute this direct quotient without computing all subgraphs and then propose bounds on the total number of graphs up to homotopy equivalence. To compare the maximal chains generated by this method, we extend Kendall's $d_K$ metric for permutations to more general graded posets and establish bounds for this metric. We then demonstrate the utility of this framework on actual brain graphs by studying the dynamics of tau proteins on the structural connectome. {We show that the proposed topological brain chain equivalence classes distinguish different simulated subtypes of Alzheimer's disease.
△ Less
Submitted 18 September, 2023; v1 submitted 22 August, 2022;
originally announced August 2022.
-
Barcodes distinguish morphology of neuronal tauopathy
Authors:
David Beers,
Despoina Goniotaki,
Diane P. Hanger,
Alain Goriely,
Heather A. Harrington
Abstract:
The geometry of neurons is known to be important for their functions. Hence, neurons are often classified by their morphology. Two recent methods, persistent homology and the topological morphology descriptor, assign a morphology descriptor called a barcode to a neuron equipped with a given function, such as the Euclidean distance from the root of the neuron. These barcodes can be converted into m…
▽ More
The geometry of neurons is known to be important for their functions. Hence, neurons are often classified by their morphology. Two recent methods, persistent homology and the topological morphology descriptor, assign a morphology descriptor called a barcode to a neuron equipped with a given function, such as the Euclidean distance from the root of the neuron. These barcodes can be converted into matrices called persistence images, which can then be averaged across groups. We show that when the defining function is the path length from the root, both the topological morphology descriptor and persistent homology are equivalent. We further show that persistence images arising from the path length procedure provide an interpretable summary of neuronal morphology. We introduce {topological morphology functions}, a class of functions similar to Sholl functions, that can be recovered from the associated topological morphology descriptor. To demonstrate this topological approach, we compare healthy cortical and hippocampal mouse neurons to those affected by progressive tauopathy. We find a significant difference in the morphology of healthy neurons and those with a tauopathy at a postsymptomatic age. We use persistence images to conclude that the diseased group tends to have neurons with shorter branches as well as fewer branches far from the soma.
△ Less
Submitted 7 April, 2022;
originally announced April 2022.
-
Neuronal Oscillations on Evolving Networks: Dynamics, Damage, Degradation, Decline, Dementia, and Death
Authors:
Alain Goriely,
Ellen Kuhl,
Christian Bick
Abstract:
Neurodegenerative diseases, such as Alzheimer's or Parkinson's disease, show characteristic degradation of structural brain networks. This degradation eventually leads to changes in the network dynamics and degradation of cognitive functions. Here, we model the progression in terms of coupled physical processes: The accumulation of toxic proteins, given by a nonlinear reaction-diffusion transport…
▽ More
Neurodegenerative diseases, such as Alzheimer's or Parkinson's disease, show characteristic degradation of structural brain networks. This degradation eventually leads to changes in the network dynamics and degradation of cognitive functions. Here, we model the progression in terms of coupled physical processes: The accumulation of toxic proteins, given by a nonlinear reaction-diffusion transport process, yields an evolving brain connectome characterized by weighted edges on which a neuronal-mass model evolves. The progression of the brain functions can be tested by simulating the resting-state activity on the evolving brain network. We show that while the evolution of edge weights plays a minor role in the overall progression of the disease, dynamic biomarkers predict a transition over a period of 10 years associated with strong cognitive decline.
△ Less
Submitted 18 September, 2020;
originally announced September 2020.
-
Likely cavitation and radial motion of stochastic elastic spheres
Authors:
L. Angela Mihai,
Thomas E. Woolley,
Alain Goriely
Abstract:
The cavitation of solid elastic spheres is a classical problem of continuum mechanics. Here, we study this problem within the context of "stochastic elasticity" where the constitutive parameters are characterised by probability density functions. We consider homogeneous spheres of stochastic neo-Hookean material, composites with two concentric stochastic neo-Hookean phases, and inhomogeneous spher…
▽ More
The cavitation of solid elastic spheres is a classical problem of continuum mechanics. Here, we study this problem within the context of "stochastic elasticity" where the constitutive parameters are characterised by probability density functions. We consider homogeneous spheres of stochastic neo-Hookean material, composites with two concentric stochastic neo-Hookean phases, and inhomogeneous spheres of locally neo-Hookean material with a radially varying parameter. In all cases, we show that the material at the centre determines the critical load at which a spherical cavity forms there. However, while under dead-load traction, a supercritical bifurcation, with stable cavitation, is obtained in a static sphere of stochastic neo-Hookean material, for the composite and radially inhomogeneous spheres, a subcritical bifurcation, with snap cavitation, is also possible. For the dynamic spheres, oscillatory motions are produced under suitable dead-load traction, such that a cavity forms and expands to a maximum radius, then collapses again to zero periodically, but not under impulse traction. Under a surface impulse, a subcritical bifurcation is found in a static sphere of stochastic neo-Hookean material and also in an inhomogeneous sphere, whereas in composite spheres, supercritical bifurcations can occur as well. Given the non-deterministic material parameters, the results can be characterised in terms of probability distributions.
△ Less
Submitted 27 January, 2020; v1 submitted 12 June, 2019;
originally announced June 2019.
-
Likely equilibria of stochastic hyperelastic spherical shells and tubes
Authors:
L. Angela Mihai,
Danielle Fitt,
Thomas E. Woolley,
Alain Goriely
Abstract:
In large deformations, internally pressurised elastic spherical shells and tubes may undergo a limit-point, or inflation, instability manifested by a rapid transition in which their radii suddenly increase. The possible existence of such an instability depends on the material constitutive model. Here, we revisit this problem in the context of stochastic incompressible hyperelastic materials, and a…
▽ More
In large deformations, internally pressurised elastic spherical shells and tubes may undergo a limit-point, or inflation, instability manifested by a rapid transition in which their radii suddenly increase. The possible existence of such an instability depends on the material constitutive model. Here, we revisit this problem in the context of stochastic incompressible hyperelastic materials, and ask the question: what is the probability distribution of stable radially symmetric inflation, such that the internal pressure always increases as the radial stretch increases? For the classic elastic problem, involving isotropic incompressible materials, there is a critical parameter value that strictly separates the cases where inflation instability can occur or not. By contrast, for the stochastic problem, we show that the inherent variability of the probabilistic parameters implies that there is always competition between the two cases. To illustrate this, we draw on published experimental data for rubber, and derive the probability distribution of the corresponding random shear modulus to predict the inflation responses for a spherical shell and a cylindrical tube made of a material characterised by this parameter.
△ Less
Submitted 10 August, 2019; v1 submitted 25 July, 2018;
originally announced August 2018.
-
Likely cavitation in stochastic elasticity
Authors:
L. Angela Mihai,
Danielle Fitt,
Thomas E. Woolley,
Alain Goriely
Abstract:
We revisit the classic problem of elastic cavitation within the framework of stochastic elasticity. For the deterministic elastic problem, involving homogeneous isotropic incompressible hyperelastic spheres under radially symmetric tension, there is a critical dead-load traction at which cavitation can occur for some materials. In addition to the well-known case of stable cavitation post-bifurcati…
▽ More
We revisit the classic problem of elastic cavitation within the framework of stochastic elasticity. For the deterministic elastic problem, involving homogeneous isotropic incompressible hyperelastic spheres under radially symmetric tension, there is a critical dead-load traction at which cavitation can occur for some materials. In addition to the well-known case of stable cavitation post-bifurcation at the critical dead load, we show the existence of unstable snap cavitation for some isotropic materials satisfying Baker-Ericksen inequalities. For the stochastic problem, we derive the probability distribution of the deformations after bifurcation. In this case, we find that, due to the probabilistic nature of the material parameters, there is always a competition between the stable and unstable states. Therefore, at a critical load, stable or unstable cavitation occurs with a given probability, and there is also a probability that the cavity may form under smaller or greater loads than the expected critical value. We refer to these phenomena as `likely cavitation'. Moreover, we provide examples of homogeneous isotropic incompressible materials exhibiting stable or unstable cavitation together with their stochastic equivalent.
△ Less
Submitted 10 August, 2019; v1 submitted 27 July, 2018;
originally announced August 2018.
-
Likely equilibria of the stochastic Rivlin cube
Authors:
L. Angela Mihai,
Thomas E. Woolley,
Alain Goriely
Abstract:
The problem of the Rivlin cube is to determine the stability of all homogeneous equilibria of an isotropic incompressible hyperelastic body under equitriaxial dead loads. Here, we consider the stochastic version of this problem where the elastic parameters are random variables following standard probability laws. Uncertainties in these parameters may arise, for example, from inherent data variatio…
▽ More
The problem of the Rivlin cube is to determine the stability of all homogeneous equilibria of an isotropic incompressible hyperelastic body under equitriaxial dead loads. Here, we consider the stochastic version of this problem where the elastic parameters are random variables following standard probability laws. Uncertainties in these parameters may arise, for example, from inherent data variation between different batches of homogeneous samples, or from different experimental tests. As for the deterministic elastic problem, we consider the following questions: what are the likely equilibria and how does their stability depend on the material constitutive law? In addition, for the stochastic model, the problem is to derive the probability distribution of deformations, given the variability of the parameters.
△ Less
Submitted 6 February, 2019; v1 submitted 24 July, 2018;
originally announced July 2018.
-
Guaranteed upper and lower bounds on the uniform load of contact problems in elasticity
Authors:
L. Angela Mihai,
Alain Goriely
Abstract:
Two mathematical models are developed within the theoretical framework of large strain elasticity for the determination of upper and lower bounds on the total strain energy of a finitely deformed hyperelastic body in unilateral contact with a rigid surface or with an elastic substrate. The model problems take the form of two continuous optimization problems with inequality constraints, the solutio…
▽ More
Two mathematical models are developed within the theoretical framework of large strain elasticity for the determination of upper and lower bounds on the total strain energy of a finitely deformed hyperelastic body in unilateral contact with a rigid surface or with an elastic substrate. The model problems take the form of two continuous optimization problems with inequality constraints, the solutions of which are used to provide an enclosure on the uniform external load acting on the body's surface away from the contact zone.
△ Less
Submitted 15 August, 2016;
originally announced August 2016.
-
Design and stability of a family of deployable structures
Authors:
Thomas Lessinnes,
Alain Goriely
Abstract:
A large family of deployable filamentary structures can be built by connecting two elastic rods along their length. The resulting structure has interesting shapes that can be stabilized by tuning the material properties of each rod. To model this structure and study its stability, we show that the equilibrium equations describing unloaded states can be derived from a variational principle. We then…
▽ More
A large family of deployable filamentary structures can be built by connecting two elastic rods along their length. The resulting structure has interesting shapes that can be stabilized by tuning the material properties of each rod. To model this structure and study its stability, we show that the equilibrium equations describing unloaded states can be derived from a variational principle. We then use a novel geometric method to study the stability of the resulting equilibria. As an example we apply the theory to establish the stability of all possible equilibria of the Bristol ladder.
△ Less
Submitted 19 April, 2016;
originally announced April 2016.
-
Geometric conditions for the positive definiteness of the second variation in one-dimensional problems
Authors:
Thomas Lessinnes,
Alain Goriely
Abstract:
Given a functional for a one-dimensional physical system, a classical problem is to minimize it by finding stationary solutions and then checking the positive definiteness of the second variation. Establishing the positive definiteness is, in general, analytically untractable. However, we show here that a global geometric analysis of the phase-plane trajectories associated with the stationary solu…
▽ More
Given a functional for a one-dimensional physical system, a classical problem is to minimize it by finding stationary solutions and then checking the positive definiteness of the second variation. Establishing the positive definiteness is, in general, analytically untractable. However, we show here that a global geometric analysis of the phase-plane trajectories associated with the stationary solutions leads to generic conditions for minimality. These results provide a straightforward and direct proof of positive definiteness, or lack thereof, in many important cases. In particular, when applied to mechanical systems, the stability or instability of entire classes of solutions can be obtained effortlessly from their geometry in phase-plane, as illustrated on a problem of a mass hanging from an elastic rod with intrinsic curvature.
△ Less
Submitted 9 April, 2016;
originally announced April 2016.
-
Controlled topological transitions in thin film phase separation
Authors:
Matthew G. Hennessy,
Victor M. Burlakov,
Alain Goriely,
Andreas Münch,
Barbara Wagner
Abstract:
In this paper the evolution of a binary mixture in a thin-film geometry with a wall at the top and bottom is considered. By bringing the mixture into its miscibility gap so that no spinodal decomposition occurs in the bulk, a slight energetic bias of the walls towards each one of the constituents ensures the nucleation of thin boundary layers that grow until the constituents have moved into one of…
▽ More
In this paper the evolution of a binary mixture in a thin-film geometry with a wall at the top and bottom is considered. By bringing the mixture into its miscibility gap so that no spinodal decomposition occurs in the bulk, a slight energetic bias of the walls towards each one of the constituents ensures the nucleation of thin boundary layers that grow until the constituents have moved into one of the two layers. These layers are separated by an interfacial region where the composition changes rapidly. Conditions that ensure the separation into two layers with a thin interfacial region are investigated based on a phase-field model. Using matched asymptotic expansions a corresponding sharp-interface problem for the location of the interface is established.
It is then argued that this newly created two-layer system is not at its energetic minimum but destabilizes into a controlled self-replicating pattern of trapezoidal vertical stripes by minimizing the interfacial energy between the phases while conserving their area. A quantitative analysis of this mechanism is carried out via a thin-film model for the free interfaces, which is derived asymptotically from the sharp-interface model.
△ Less
Submitted 27 February, 2014;
originally announced February 2014.
-
The mechanics of a chain or ring of spherical magnets
Authors:
Cameron L. Hall,
Dominic Vella,
Alain Goriely
Abstract:
Strong magnets, such as neodymium-iron-boron magnets, are increasingly being manufactured as spheres. Because of their dipolar characters, these spheres can easily be arranged into long chains that exhibit mechanical properties reminiscent of elastic strings or rods. While simple formulations exist for the energy of a deformed elastic rod, it is not clear whether or not they are also appropriate f…
▽ More
Strong magnets, such as neodymium-iron-boron magnets, are increasingly being manufactured as spheres. Because of their dipolar characters, these spheres can easily be arranged into long chains that exhibit mechanical properties reminiscent of elastic strings or rods. While simple formulations exist for the energy of a deformed elastic rod, it is not clear whether or not they are also appropriate for a chain of spherical magnets. In this paper, we use discrete-to-continuum asymptotic analysis to derive a continuum model for the energy of a deformed chain of magnets based on the magnetostatic interactions between individual spheres. We find that the mechanical properties of a chain of magnets differ significantly from those of an elastic rod: while both magnetic chains and elastic rods support bending by change of local curvature, nonlocal interaction terms also appear in the energy formulation for a magnetic chain. This continuum model for the energy of a chain of magnets is used to analyse small deformations of a circular ring of magnets and hence obtain theoretical predictions for the vibrational modes of a circular ring of magnets. Surprisingly, despite the contribution of nonlocal energy terms, we find that the vibrations of a circular ring of magnets are governed by the same equation that governs the vibrations of a circular elastic ring.
△ Less
Submitted 12 October, 2013;
originally announced October 2013.