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

Pair approximation for a model of vertical and horizontal transmission of parasites.

We apply Stochastic Dynamics method for a differential equations model, proposed by Marc Lipsitch and collaborators (Proc. R. Soc. Lond. B 260, 321, 1995), for which the transmission dynamics of parasites occurs from a parent to its offspring (vertical transmission), and by contact with infected host (horizontal transmission). Herpes, Hepatitis and AIDS are examples of diseases for which both horizontal and vertical transmission occur simultaneously during the virus spreading. Understanding the role of each type of transmission in the infection prevalence on a susceptible host population may provide some information about the factors that contribute for the eradication and/or control of those diseases. We present a pair mean-field approximation obtained from the master equation of the model. The pair approximation is formed by the differential equations of the susceptible and infected population densities and the differential equations of pairs that contribute to the former ones. In terms of the model parameters, we obtain the conditions that lead to the disease eradication, and set up the phase diagram based on the local stability analysis of fixed points. We also perform Monte Carlo simulations of the model on complete graphs and Erdös-Rényi graphs in order to investigate the influence of population size and neighborhood on the previous mean-field results; by this way, we also expect to evaluate the contribution of vertical and horizontal transmission on the elimination of parasite. Pair Approximation for a Model of Vertical and Horizontal Transmission of Parasites.

Quantification of cellular action fibers alignment from confocal microscopic images.

The cellular rheology has recently undergone a rapid development with particular attention to the cytoskeleton mechanical properties and its main components - actin filaments, intermediate filaments, microtubules and crosslinked proteins. However it is not clear what are the cellular structural changes that directly affect the cell mechanical properties. Thus, in this work, we aimed to quantify the structural rearrangement of these fibers that may emerge in changes in the cell mechanics. We created an image analysis platform to study smooth muscle cells from different arteries: aorta, mammary, renal, carotid and coronary and processed respectively 31, 29, 31, 30 and 35 cell image obtained by confocal microscopy. The platform was developed in Matlab (MathWorks) and it uses the Sobel operator to determine the actin fiber image orientation of the cell, labeled with phalloidin. The Sobel operator is used as a filter capable of calculating the pixel brightness gradient, point to point, in the image. The operator uses vertical and horizontal convolution kernels to calculate the magnitude and the angle of the pixel intensity gradient. The image analysis followed the sequence: (1) opens a given cells image set to be processed; (2) sets a fix threshold to eliminate noise, based on Otsu's method; (3) detect the fiber edges in the image using the Sobel operator; and (4) quantify the actin fiber orientation. Our first result is the probability distribution II(Δθ) to find a given fiber angle deviation (Δθ) from the main cell fiber orientation θ0. The II(Δθ) follows an exponential decay II(Δθ) = Aexp(-αΔθ) regarding to its θ0. We defined and determined a misalignment index α of the fibers of each artery kind: coronary αCo = (1.72 ‘+ or =’ 0.36)rad POT -1; renal αRe = (1.43 + or - 0.64)rad POT -1; aorta αAo = (1.42 + or - 0.43)rad POT -1; mammary αMa = (1.12 + or - 0.50)rad POT -1; and carotid αCa = (1.01 + or - 0.39)rad POT -1. The α of coronary and carotid are statistically different (p < 0.05) among all analyzed cells. We discussed our results correlating the misalignment index data with the experimental cell mechanical properties obtained by using Optical Magnetic Twisting Cytometry with the same group of cells.

Rotational diffusion membrane and fluorescence anisotropy.

Structural properties of model membranes, such as lipid vesicles, may be investigated through the addition of fluorescent probes. After incorporation, the fluorescent molecules are excited with linearly polarized light and the fluorescence emission is depolarized due to translational as well as rotational diffusion during the lifetime of the excited state. The monitoring of emitted light is undertaken through the technique of time-resolved fluorescence: the intensity of the emitted light informs on fluorescence decay times, and the decay of the components of the emitted light yield rotational correlation times which inform on the fluidity of the medium. The fluorescent molecule DPH, of uniaxial symmetry, is rather hydrophobic and has collinear transition and emission moments. It has been used frequently as a probe for the monitoring of the fluidity of the lipid bilayer along the phase transition of the chains. The interpretation of experimental data requires models for localization of fluorescent molecules as well as for possible restrictions on their movement. In this study, we develop calculations for two models for uniaxial diffusion of fluorescent molecules, such as DPH, suggested in several articles in the literature. A zeroth order test model consists of a free randomly rotating dipole in a homogeneous solution, and serves as the basis for the study of the diffusion of models in anisotropic media. In the second model, we consider random rotations of emitting dipoles distributed within cones with their axes perpendicular to the vesicle spherical geometry. In the third model, the dipole rotates in the plane of the of bilayer spherical geometry, within a movement that might occur between the monolayers forming the bilayer. For each of the models analysed, two methods are used by us in order to analyse the rotational diffusion: (I) solution of the corresponding rotational diffusion equation for a single molecule, taking into account the boundary conditions imposed by the models, for the probability of the fluorescent molecule to be found with a given configuration at time t. Considering the distribution of molecules in the geometry proposed, we obtain the analytical expression for the fluorescence anisotropy, except for the cone geometry, for which the solution is obtained numerically; (II) numerical simulations of a restricted rotational random walk in the two geometries corresponding to the two models. The latter method may be very useful in the cases of low-symmetry geometries or of composed geometries.

Shape resonance spectrum of cytosine-guanine pairs.

Single and double strand breaks in DNA can be caused by low-energy electrons, the most abundant secondary products of the interaction of ionizing radiation to the biological matter. Attachment of these electrons to biomolecules lead to the formation of transient negative ions (TNIs) [1], often referred to as resonances, a process that may lead to significant vibrational excitation and dissociation. In the present study, we employ the parallel version [2] of the Schwinger Multichannel Method implemented with pseudopotentials [3] to obtain the shape resonance spectrum of cytosine-guanine (CG) pairs, with special attention to π* transient anion states. Recent experimental studies pointed out a quasi-continuum vibrational excitation spectrum for electron collisions against formic acid dimers [4], suggesting that electron attachment into π* valence orbitals could induce proton transfer in these dimers. In addition, our previous studies on the shape resonance spectra of the hydrogen-bonded complexes comprising formic acid and formamide units indicated interesting electron delocalization (localization) effects arising from the presence (absence) of inversion symmetry centers in the complexes [5]. In the present work, we extend the studies on hydrogen-bonded complexes to the CG pair, where localization of ¼¤ anions would be expected, based on the previous results. References [1]. B. Boudaïffa, P. Cloutier, D. Hunting, M. A. Huels, L. Sanche, Science 287, 1658 (2000). [2]. J. S. dos Santos, R. F. da Costa , M. T. do N. Varella, J. Chem. Phys. 136, 084307 (2012). [3]. M. H. F. Bettega, L. G. Ferreira, M. A. P. Lima, Phys. Rev. A 47, 1111 (1993). [4]. M. Allan, Phys. Rev. Lett. 98, 123201 (2007). [5]. T. C. Freitas, S. dA. Sanchez, M. T. do N. Varella, M. H. F. Bettega, Phys. Rev. A 84, 062714 (2011).

Study of surface free energy on a lattice-gas model applied to Liquid bridge

The pulmonary crackling and the formation of liquid bridges are problems that for centuries have been attracting the attention of scientists. In order to study these phenomena, it was developed a canonical cubic lattice-gas­ like model to explain the rupture of liquid bridges in lung airways [A. Alencar et al., 2006, PRE]. Here, we further develop this model and add entropy analysis to study thermodynamic properties, such as free energy and force. The simulations were performed using the Monte Carlo method with Metropolis algorithm. The exchange between gas and liquid particles were performed randomly according to the Kawasaki dynamics and weighted by the Boltzmann factor. Each particle, which can be solid (s), liquid (l) or gas (g), has 26 neighbors: 6 + 12 + 8, with distances 1, √2 and √3, respectively. The energy of a lattice's site m is calculated by the following expression: Em = ∑k=126 Ji(m)j(k) in witch (i, j) = g, l or s. Specifically, it was studied the surface free energy of the liquid bridge, trapped between two planes, when its height is changed. For that, was considered two methods. First, just the internal energy was calculated. Then was considered the entropy. It was fond no difference in the surface free energy between this two methods. We calculate the liquid bridge force between the two planes using the numerical surface free energy. This force is strong for small height, and decreases as the distance between the two planes, height, is increased. The liquid-gas system was also characterized studying the variation of internal energy and heat capacity with the temperature. For that, was performed simulation with the same proportion of liquid and gas particle, but different lattice size. The scale of the liquid-gas system was also studied, for low temperature, using different values to the interaction Jij.

Mechanical alterations in airway smooth muscle after interaction with indoor pollutant - A mathematical model

The viscoelasticity of mammalian lung is determined by the mechanical properties and structural regulation of the airway smooth muscle (ASM). The exposure to polluted air may deteriorate these properties with harmful consequences to individual health. Formaldehyde (FA) is an important indoor pollutant found among volatile organic compounds. This pollutant permeates through the smooth muscle tissue forming covalent bonds between proteins in the extracellular matrix and intracellular protein structure changing mechanical properties of ASM and inducing asthma symptoms, such as airway hyperresponsiveness, even at low concentrations. In the experimental scenario, the mechanical effect of FA is the stiffening of the tissue, but the mechanism behind this effect is not fully w1derstood. Thus, the aim of this study is to reproduce the mechanical behavior of the ASM, such as contraction and stretching, under FA action or not. For this, it was created a two-dimensional viscoelastic network model based on Voronoi tessellation solved using Runge-Kutta method of fourth order. The equilibrium configuration was reached when the forces in different parts of the network were equal. This model simulates the mechanical behavior of ASM through of a network of dashpots and springs. This dashpot-spring mechanical coupling mimics the composition of the actomyosin machinery of ASM through the contraction of springs to a minimum length. We hypothesized that formation of covalent bonds, due to the FA action, can be represented in the model by a simple change in the elastic constant of the springs, while the action of methacholinc (MCh) reduce the equilibrium length of the spring. A sigmoid curve of tension as a function of MCh doses was obtained, showing increased tension when the muscle strip was exposed to FA. Our simulations suggest that FA, at a concentration of 0.1 ppm, can affect the elastic properties of the smooth muscle fibers by a factor of 120%. We also analyze the dynamic mechanical properties, observing the viscous and elastic behavior of the network. Finally, the proposed model, although simple, ir1corporates the phenomenology of both MCh and FA and reproduces experirnental results observed with ir1 vitro exposure of smooth muscle to .FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells ir1 a tissue level model. The model can also be used in different biological scales.

Statistical models for experimental model membranes.

Biological membranes are constituted from lipid bilayers and proteins. Investigation of protein-membrane interaction, essential for biological function of cells, must rest upon solid knowledge of lipid bilayer behavior. Thus, extensive studies of an experimental model for membranes, lipid bilayers in water solution, have been undertaken in the last decades. These systems present structural, thermal and electrical properties which depend on temperature, ionic strength or concentration. In this talk, we shall discuss statistical models for lipid bilayers, as well as the relation between their properties and results for properties of lipid dispersions investigated by the laboratories supervised by Teresa Lamy (IF-USP) and Amando Ito (FFCL-USP).

Stochastic quantization of the spherical model and supersymmetry.

In this work, we reported some results about the stochastic quantization of the spherical model. We started by reviewing some basic aspects of this method with emphasis in the connection between the Langevin equation and the supersymmetric quantum mechanics, aiming at the application of the corresponding connection to the spherical model. An intuitive idea is that when applied to the spherical model this gives rise to a supersymmetric version that is identified with one studied in Phys. Rev. E 85, 061109, (2012). Before investigating in detail this aspect, we studied the stochastic quantization of the mean spherical model that is simpler to implement than the one with the strict constraint. We also highlight some points concerning more traditional methods discussed in the literature like canonical and path integral quantization. To produce a supersymmetric version, grounded in the Nicolai map, we investigated the stochastic quantization of the strict spherical model. We showed in fact that the result of this process is an off-shell supersymmetric extension of the quantum spherical model (with the precise supersymmetric constraint structure). That analysis establishes a connection between the classical model and its supersymmetric quantum counterpart. The supersymmetric version in this way constructed is a more natural one and gives further support and motivations to investigate similar connections in other models of the literature.

The role of asymmetries in rock-paper-scissors biodiversity models.

The maintenance of biodiversity is a long standing puzzle in ecology. It is a classical result that if the interactions of the species in an ecosystem are chosen in a random way, then complex ecosystems can't sustain themselves, meaning that the structure of the interactions between the species must be a central component on the preservation of biodiversity and on the stability of ecosystems. The rock-paper-scissors model is one of the paradigmatic models that study how biodiversity is maintained. In this model 3 species dominate each other in a cyclic way (mimicking a trophic cycle), that is, rock dominates scissors, that dominates paper, that dominates rock. In the original version of this model, this dominance obeys a 'Z IND 3' symmetry, in the sense that the strength of dominance is always the same. In this work, we break this symmetry, studying the effects of the addition of an asymmetry parameter. In the usual model, in a two dimensional lattice, the species distribute themselves according to spiral patterns, that can be explained by the complex Landau-Guinzburg equation. With the addition of asymmetry, new spatial patterns appear during the transient and the system either ends in a state with spirals, similar to the ones of the original model, or in a state where unstable spatial patterns dominate or in a state where only one species survives (and biodiversity is lost).

Low energy elastic and electronically inelastic electron scattering from biomolecules.

Reactions initiated by collisions with low-energy secondary electrons has been found to be the prominent mechanism toward the radiation damage on living tissues through DNA strand breaks. Now it is widely accepted that during the interaction with these secondary species the selective breaking of chemical bonds is triggered by dissociative electron attachment (DEA), that is, the capture of the incident electron and the formation of temporary negative ion states [1,2,3]. One of the approaches largely used toward a deeper understanding of the radiation damage to DNA is through modeling of DEA with its basic constituents (nucleotide bases, sugar and other subunits). We have tried to simplify this approach and attempt to make it comprehensible at a more fundamental level by looking at even simple molecules. Studies involving organic systems such as carboxylic acids, alcohols and simple ¯ve-membered heterocyclic compounds are taken as starting points for these understanding. In the present study we investigate the role played by elastic scattering and electronic excitation of molecules on electron-driven chemical processes. Special attention is focused on the analysis of the in°uence of polarization and multichannel coupling e®ects on the magnitude of elastic and electronically inelastic cross-sections. Our aim is also to investigate the existence of resonances in the elastic and electronically inelastic channels as well as to characterize them with respect to its type (shape, core-excited or Feshbach), symmetry and position. The relevance of these issues is evaluated within the context of possible applications for the modeling of discharge environments and implications in the understanding of mutagenic rupture of DNA chains. The scattering calculations were carried out with the Schwinger multichannel method (SMC) [4] and its implementation with pseudopotentials (SMCPP) [5] at di®erent levels of approximation for impact energies ranging from 0.5 eV to 30 eV. References [1] B. Boudai®a, P. Cloutier, D. Hunting, M. A. Huels and L. Sanche, Science 287, 1658 (2000). [2] X. Pan, P. Cloutier, D. Hunting and L. Sanche, Phys. Rev. Lett. 90, 208102 (2003). [3] F. Martin, P. D. Burrow, Z. Cai, P. Cloutier, D. Hunting and L. Sanche, Phys. Rev. Lett. 93, 068101 (2004). [4] K. Takatsuka and V. McKoy, Phys. Rev. A 24, 2437 (1981); ibid. Phys. Rev. A 30, 1734 (1984). [5] M. H. F. Bettega, L. G. Ferreira and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993).

Atomistic Study of Water Confined in Silica.

In this work, we have used a combined of atomistic simulation methods to explore the effects of confinement of water molecules between silica surfaces. Firstly, the mechanical properties of water severe confined (~3A) between two silica alpha-quartz was determined based on first principles calculations within the density functional theory (DFT). Simulated annealing methods were employed due to the complex potential energry surface, and the difficulties to avoid local minima. Our results suggest that much of the stiffness of the material (46%) remains, even after the insertion of a water monolayer in the silica. Secondly, in order to access typical time scales for confined systems, classical molecular dynamics was used to determine the dynamical properties of water confined in silica cylindrical pores, with diameters varying from 10 to 40A. in this case we have varied the passivation of the silica surface, from 13% to 100% of SiOH, and the other terminations being SiOH2 and SiOH3, the distribution of the different terminations was obtained with a Monte Carlo simulation. The simulations indicates a lowering of the diffusion coefficientes as the diameter decreases, due to the structuration of hydrogen bonds of water molecules; we have also obtained the density profiles of the confined water and the interfacial tension.

Study of the Influence of Localized Vibrational Modes in Charge Transport Properties at Nanoscale Systems.

In molecular and atomic devices the interaction between electrons and ionic vibrations has an important role in electronic transport. The electron-phonon coupling can cause the loss of the electron's phase coherence, the opening of new conductance channels and the suppression of purely elastic ones. From the technological viewpoint phonons might restrict the efficiency of electronic devices by energy dissipation, causing heating, power loss and instability. The state of the art in electron transport calculations consists in combining ab initio calculations via Density Functional Theory (DFT) with Non-Equilibrium Green's Function formalism (NEGF). In order to include electron-phonon interactions, one needs in principle to include a self-energy scattering term in the open system Hamiltonian which takes into account the effect of the phonons over the electrons and vice versa. Nevertheless this term could be obtained approximately by perturbative methods. In the First Born Approximation one considers only the first order terms of the electronic Green's function expansion. In the Self-Consistent Born Approximation, the interaction self-energy is calculated with the perturbed electronic Green's function in a self-consistent way. In this work we describe how to incorporate the electron-phonon interaction to the SMEAGOL program (Spin and Molecular Electronics in Atomically Generated Orbital Landscapes), an ab initio code for electronic transport based on the combination of DFT + NEGF. This provides a tool for calculating the transport properties of materials' specific system, particularly in molecular electronics. Preliminary results will be presented, showing the effects produced by considering the electron-phonon interaction in nanoscale devices.

Spin Orbit Effects in the Electronic Transport Properties of Adsorbed Graphene Nanoribbons

Graphene has received great attention due to its exceptional properties, which include corners with zero effective mass, extremely large mobilities, this could render it the new template for the next generation of electronic devices. Furthermore it has weak spin orbit interaction because of the low atomic number of carbon atom in turn results in long spin coherence lengths. Therefore, graphene is also a promising material for future applications in spintronic devices - the use of electronic spin degrees of freedom instead of the electron charge. Graphene can be engineered to form a number of different structures. In particular, by appropriately cutting it one can obtain 1-D system -with only a few nanometers in width - known as graphene nanoribbon, which strongly owe their properties to the width of the ribbons and to the atomic structure along the edges. Those GNR-based systems have been shown to have great potential applications specially as connectors for integrated circuits. Impurities and defects might play an important role to the coherence of these systems. In particular, the presence of transition metal atoms can lead to significant spin-flip processes of conduction electrons. Understanding this effect is of utmost importance for spintronics applied design. In this work, we focus on electronic transport properties of armchair graphene nanoribbons with adsorbed transition metal atoms as impurities and taking into account the spin-orbit effect. Our calculations were performed using a combination of density functional theory and non-equilibrium Greens functions. Also, employing a recursive method we consider a large number of impurities randomly distributed along the nanoribbon in order to infer, for different concentrations of defects, the spin-coherence length.

First principles calculations for vacancies and antisites in PbSe and PbTe: bulk and nanowire.

The energetic stability and the electronic properties of vacancies (VX) and antisites (XY) in PbSe and PbTe are investigated. PbSe and PbTe are narrow band gap semiconductors and have the potential to be used in infrared detectors, laser, and diodes. They are also of special interest for thermoelectric devices (TE). The calculations are based in the Density Functional Theory (DFT) and the General Gradient Approximation (GGA) for the exchange-correlation term, as implemented in the VASP code. The core and valence electrons are described by the Projected Augmented Wave (PAW) and the Plane Wave (PW) methods, respectively. The defects are studied in the bulk and nanowire (NW) system. Our results show that intrinsec defects (vacancies and antisites) in PbTe have lower formation energies in the NW as compared to the bulk and present a trend in migrate to the surface of the NW. For the PbSe we obtain similar results when compare the formation energy for the bulk and NW. However, the Pb vacancy and the antisites are more stable in the core of the NW. The intrinsec defects are shallow defects for the bulk system. For both PbSe and PbTe VPb is a shallow acceptor defect and VSe and VT e are shallow donor defects for the PbSe and PbTe, respectively. Similar electronic properties are observed for the antisites. For the Pb in the anion site we obtain an n-type semiconductor for both PbSe and PbTe, SeP b is a p-type for the PbSe, and T eP b is a n-type for PbTe. Due the quantum con¯nement effects present in the NW (the band gap open), these defects have different electronic properties for the NW as compared to the bulk. Now these defects give rise to electronic levels in the band gap of the PbTe NW and the VT e present a metallic character. For the PbSe NW a p-type and a n-type semiconductor is obtained for the VP b and P bSe, respectively. On the other hand, deep electronic levels are present in the band gap for the VSe and SePb. These results show that due an enhanced in the electronic density of states (DOS) near the Fermi energy, the defective PbSe and PbTe are candidates for efficient TE devices.