Numeric reconstruction of cytoskeleton with finite element method and topology optimization method

The importance of mechanical aspects related to cell activity and its environment is becoming more evident due to their influence in stem cell differentiation and in the development of diseases such as atherosclerosis. The mechanical tension homeostasis is related to normal tissue behavior and its lack may be related to the formation of cancer, which shows a higher mechanical tension. Due to the complexity of cellular activity, the application of simplified models may elucidate which factors are really essential and which have a marginal effect. The development of a systematic method to reconstruct the elements involved in the perception of mechanical aspects by the cell may accelerate substantially the validation of these models. This work proposes the development of a routine capable of reconstructing the topology of focal adhesions and the actomyosin portion of the cytoskeleton from the displacement field generated by the cell on a flexible substrate. Another way to think of this problem is to develop an algorithm to reconstruct the forces applied by the cell from the measurements of the substrate displacement, which would be characterized as an inverse problem. For these kind of problems, the Topology Optimization Method (TOM) is suitable to find a solution. TOM is consisted of an iterative application of an optimization method and an analysis method to obtain an optimal distribution of material in a fixed domain. One way to experimentally obtain the substrate displacement is through Traction Force Microscopy (TFM), which also provides the forces applied by the cell. Along with systematically generating the distributions of focal adhesion and actin-myosin for the validation of simplified models, the algorithm also represents a complementary and more phenomenological approach to TFM. As a first approximation, actin fibers and flexible substrate are represented through two-dimensional linear Finite Element Method. Actin contraction is modeled as an initial stress of the FEM elements. Focal adhesions connecting actin and substrate are represented by springs. The algorithm was applied to data obtained from experiments regarding cytoskeletal prestress and micropatterning, comparing the numerical results to the experimental ones

Simultaneous aligning and smoothing of Surface triangulations

[EN]In this work we develop a procedure to deform a given surface triangulation to obtain its alignment with interior curves. These curves are defined by splines in a parametric space and, subsequently, mapped to the surface triangulation. We have restricted our study to orthogonal mapping, so we require the curves to be included in a patch of the surface that can be orthogonally projected onto a plane (our parametric space). For example, the curves can represent interfaces between different materials or boundary conditions, internal boundaries or feature lines. Another setting in which this procedure can be used is the adaption of a reference mesh to changing curves in the course of an evolutionary process...

Performance evaluation of a parallel algorithm for simultaneous untangling and smoothing of tetrahedral meshes

[EN]A new parallel algorithm for simultaneous untangling and smoothing of tetrahedral meshes is proposed in this paper. We provide a detailed analysis of its performance on shared-memory many-core computer architectures. This performance analysis includes the evaluation of execution time, parallel scalability, load balancing, and parallelism bottlenecks. Additionally, we compare the impact of three previously published graph coloring procedures on the performance of our parallel algorithm. We use six benchmark meshes with a wide range of sizes. Using these experimental data sets, we describe the behavior of the parallel algorithm for different data sizes. We demonstrate that this algorithm is highly scalable when it runs on two different high-performance many-core computers with up to 128 processors...

Combinatorial and Robust Optimisation Models and Algorithms for Railway Applications

This thesis deals with an investigation of combinatorial and robust optimisation models to solve railway problems. Railway applications represent a challenging area for operations research. In fact, most problems in this context can be modelled as combinatorial optimisation problems, in which the number of feasible solutions is finite. Yet, despite the astonishing success in the field of combinatorial optimisation, the current state of algorithmic research faces severe difficulties with highly-complex and data-intensive applications such as those dealing with optimisation issues in large-scale transportation networks. One of the main issues concerns imperfect information. The idea of Robust Optimisation, as a way to represent and handle mathematically systems with not precisely known data, dates back to 1970s. Unfortunately, none of those techniques proved to be successfully applicable in one of the most complex and largest in scale (transportation) settings: that of railway systems. Railway optimisation deals with planning and scheduling problems over several time horizons. Disturbances are inevitable and severely affect the planning process. Here we focus on two compelling aspects of planning: robust planning and online (real-time) planning.

The C-4-Dicarboxylate carriers DcuB and DctA of Escherichia coli: function as cosensors and topology

Das fakultativ anaerobe Enterobakterium Escherichia coli nutzt C4-Dicarboxylate sowohl unter aeroben als auch anaeroben Bedingungen als Kohlenstoff- und Energiequelle. Die Aufnahme der C4-Dicarboxylaten und die Energiekonservierung mittels Fumaratatmung wird durch das Zweikomponentensystem DcuSR reguliert. Die Sensorhistidinkinase DcuS und der nachgeschaltete Responseregulator DcuR aktivieren bei Verfügbarkeit von C4-Dicarboxylaten die Expression der Gene für den Succinat Transporter DctA, den anaeroben Fumarat/Succinat Antiporter DcuB, die Fumarase B sowie die Fumaratreduktase FrdABCD. Die Transportproteine DctA und DcuB wiederum regulieren die Expression der DcuSR-abhängigen Gene negativ. Fehlen von DctA oder DcuB resultiert bereits ohne Effektor in einer maximalen Expression von dctA bzw. dcuB. Durch gerichtete und ungerichtete Mutagenese wurde gezeigt, dass die Transportfunktion des Carriers DcuB unabhängig von seiner regulatorischen Funktion ist. DcuB kann daher als Cosensor des DcuSR Systems angesehen werden.rnUnter Verwendung von Reportergenfusionen von C-terminal verkürzten Konstrukten von DcuB mit der Alkalischen Phosphatase und der β-Galactosidase wurde die Topologie des Multitransmembranproteins DcuB bestimmt. Zusätzlich wurde die Zugänglichkeit bestimmter Aminosäurereste durch chemische Modifikation mit membran-durchlässigen und membran-undurchlässigen Thiolreagenzien untersucht. Die erhaltenen Ergebnisse deuten auf die Existenz eines tief in die Membran reichenden, hydrophilen Kanal hin, welcher zum Periplasma hin geöffnet ist. Mit Hilfe der Topologie-Studien, des Hydropathie-Blots und der Sekundärstruktur-Vorhersage wurde ein Modell des Carriers erstellt. DcuB besitzt kurze, periplasmatisch liegende Proteinenden, die durch 12 Transmembranhelices und zwei große hydrophile Schleifen jeweils zwischen TM VII/VIII und TM XI/XII verbunden sind. Die regulatorisch relevanten Reste K353, T396 und D398 befinden sich innerhalb von TM XI sowie auf der angrenzenden cytoplasmatischen Schleife XI-XII. Unter Berücksichtigung der strukturellen und funktionellen Aspekte wurde ein Regulationsmodell erstellt, welches die gemeinsam durch DcuB und DcuS kontrollierte C4-Dicarboxylat-abhängige Genexpression darstellt. rnDer Effekt von DctA und DcuSR auf die Expression einer dctA´-´lacZ Reportergenfusion und auf die aerobe C4-Dicarboxylat-Aufnahme wurde untersucht. In-vivo FRET-Messungen weisen auf eine direkte Wechselwirkung zwischen dem Carrier DctA und dem Sensor DcuS hin. Dieses Ergebnis stützt die Theorie der Regulation von DcuS durch C4-Dicarboxylate und durch die Cosensoren DctA bzw. DcuB mittels direkter Protein-Protein Interaktion.rn

Hydrodynamical simulations of early-type galaxies: effects of galaxy shape and stellar dynamics on hot coronae

Early-Type galaxies (ETGs) are embedded in hot (10^6-10^7 K), X-ray emitting gaseous haloes, produced mainly by stellar winds and heated by Type Ia supernovae explosions, by the thermalization of stellar motions and occasionally by the central super-massive black hole (SMBH). In particular, the thermalization of the stellar motions is due to the interaction between the stellar and the SNIa ejecta and the hot interstellar medium (ISM) already residing in the ETG. A number of different astrophysical phenomena determine the X-ray properties of the hot ISM, such as stellar population formation and evolution, galaxy structure and internal kinematics, Active Galactic Nuclei (AGN) presence, and environmental effects. With the aid of high-resolution hydrodynamical simulations performed on state-of-the-art galaxy models, in this Thesis we focus on the effects of galaxy shape, stellar kinematics and star formation on the evolution of the X-ray coronae of ETGs. Numerical simulations show that the relative importance of flattening and rotation are functions of the galaxy mass: at low galaxy masses, adding flattening and rotation induces a galactic wind, thus lowering the X-ray luminosity; at high galaxy masses the angular momentum conservation keeps the central regions of rotating galaxies at low density, whereas in non-rotating models a denser and brighter atmosphere is formed. The same dependence from the galaxy mass is present in the effects of star formation (SF): in light galaxies SF contributes to increase the spread in Lx, while at high galaxy masses the halo X-ray properties are marginally sensitive to SF effects. In every case, the star formation rate at the present epoch quite agrees with observations, and the massive, cold gaseous discs are partially or completely consumed by SF on a time-scale of few Gyr, excluding the presence of young stellar discs at the present epoch.

Active regulation of cerebral venous tone: simultaneous arterial and venous transcranial Doppler sonography during a Valsalva manoeuvre

The aim of this study was to analyse the cerebral venous outflow in relation to the arterial inflow during a Valsalva manoeuvre (VM). In 19 healthy volunteers (mean age 24.1 +/- 2.6 years), the middle cerebral artery (MCA) and the straight sinus (SRS) were insonated by transcranial Doppler sonography. Simultaneously the arterial blood pressure was recorded using a photoplethysmographic method. Two VM of 10 s length were performed per participant. Tracings of the variables were then transformed to equidistantly re-sampled data. Phases of the VM were analysed regarding the increase of the flow velocities and the latency to the peak. The typical four phases of the VM were also found in the SRS signal. The relative flow velocity (FV) increase was significantly higher in the SRS than in the MCA for all phases, particularly that of phase IV (p < 0.01). Comparison of the time latency of the VM phases of the MCA and SRS only showed a significant difference for phase I (p < 0.01). In particular, there was no significant difference for phase IV (15.8 +/- 0.29 vs. 16.0 +/- 0.28 s). Alterations in venous outflow in phase I are best explained by a cross-sectional change of the lumen of the SRS, while phases II and III are compatible with a Starling resistor. However, the significantly lager venous than the arterial overshoot in phase IV may be explained by the active regulation of the venous tone.

Dimorphism in Shell Shape and Strength in Two Species of Emydid Turtle

This study investigates the mechanical implications of shell shape differences between males and females of two North American turtle species: Chrysemys picta and Glyptemys insculpta. These species show patterns of sexual dimorphism that are common to many species of turtle. Females have wider and more highly domed shells, whereas males tend to have flatter, more streamlined shells. In addition, the males of many terrestrial species have concave plastra, most likely to accommodate the domed shells of the females while mating. The purpose of this study was to determine whether the known morphological differences in male and female turtle shells are also associated with differences in shell strength. Landmark coordinate data were collected from the shells of males and females of both species. These data were used to create digital models of each shell for finite-element (FE) analysis. FE models were generated by transforming a single base model of a turtle shell to match the shapes of each specimen examined in this study. All models were assigned the same material properties and restraints. Twelve load cases, each representing a predator’s bite at a different location on the carapace, were applied separately to the models. Subsequently, Von Mises stresses were extracted for each element of each model. Overall, the shells of females of both species exhibited significantly lower maximum and average stresses for a given load than those of their male counterparts. Male G. insculpta exhibited significant increases in stresses because of the concave shape of their plastra. We suggest that the mechanical implications of shell shape differences between males and females may have a large impact on many aspects of the biology of these turtle species.

Is there a role for simultaneous hepatic and colorectal resections? A contemporary view from NSQIP

The optimal timing of primary and metastatic tumor management in patients with synchronous hepatic colorectal metastases remains controversial. We aimed to compare perioperative outcomes of simultaneous colorectal/liver resection (SCLR) with isolated resections utilizing a national clinical database.

Statistical finite element model for bone shape and biomechanical properties

We present a framework for statistical finite element analysis combining shape and material properties, and allowing performing statistical statements of biomechanical performance across a given population. In this paper, we focus on the design of orthopaedic implants that fit a maximum percentage of the target population, both in terms of geometry and biomechanical stability. CT scans of the bone under consideration are registered non-rigidly to obtain correspondences in position and intensity between them. A statistical model of shape and intensity (bone density) is computed by means of principal component analysis. Afterwards, finite element analysis (FEA) is performed to analyse the biomechanical performance of the bones. Realistic forces are applied on the bones and the resulting displacement and bone stress distribution are calculated. The mechanical behaviour of different PCA bone instances is compared.