Anomalous diffusion evidenced by particle tracking at high frame rates.

Diffusion is a common phenomenon in nature and generally is associated with a system trying to reach a local or a global equilibrium state, as a result of highly irregular individual particle motion. Therefore it is of fundamental importance in physics, chemistry and biology. Particle tracking in complex fluids can reveal important characteristics of its properties. In living cells, we coat the microbead with a peptide (RGD) that binds to integrin receptors at the plasma membrane, which connects to the CSK. This procedure is based on the hypothesis that the microsphere can move only if the structure where it is attached move as well. Then, the observed trajectory of microbeads is a probe of the cytoskeleton (CSK), which is governed by several factors, including thermal diffusion, pressure gradients, and molecular motors. The possibility of separating the trajectories into passive and active diffusion may give information about the viscoelasticity of the cell structure and molecular motors activity. And also we could analyze the motion via generalized Stokes-Einstein relation, avoiding the use of any active techniques. Usually a 12 to 16 Frames Per Second (FPS) system is used to track the microbeads in cell for about 5 minutes. Several factors make this FPS limitation: camera computer communication, light, computer speed for online analysis among others. Here we used a high quality camera and our own software, developed in C++ and Linux, to reach high FPS. Measurements were conducted with samples for 10£ and 20£ objectives. We performed sequentially images with different intervals, all with 2 ¹s exposure. The sequences of intervals are in milliseconds: 4 5 ms (maximum speed) 14, 25, 50 and 100 FPS. Our preliminary results highlight the difference between passive and active diffusion, since the passive diffusion is represented by a Gaussian in the distribution of displacements of the center of mass of individual beads between consecutive frames. However, the active process, or anomalous diffusion, shows as long tails in the distribution of displacements.

Comparison of physical properties of two exogenous surfactants: new parameter.

The pulmonary surfactant has essential physical properties for normal lung function. The most important property is the surface tension. In this work, it was evaluated the surface tension of two commercial exogenous surfactants used in surfactant replacement therapy, poractant alfa (Curosurf, Chiesi Farmaceuticals, Italy) and beractant (Survanta, Abbott Laboratories, USA) using new parameters. A Langmuir film balance (Minitrough, KSV Instruments, Finland) was used to measure surface tension of poractant alfa and beractant samples. For both samples, we prepared a solution of 1 mg/m dissolved in chloroform (100π`), which was applied over a subphase of milli-Q water (175 ml) in the chamber of the balance. The chamber has two moving barriers that can change its surface area between a maximal value of 112.5 cm 2 , and a minimal value of 22.5 cm 2, defining a balance cycle. Each surfactant had its surface tension evaluated during 20 balance cycles for three times. Four quantities were calculated from the experiment: Minimum Surface Tension (MTS), defined as the surface tension at minimal surface area during the first cycle; Mean Work Cycle (MWC), defined as the mean hysteresis area of the measured surface tension curve of the last 16 balance cycles; Critical Active Surface Area in Compression (CASAC) or in Expansion (CASAE), defined as the maximal chamber area where the surfactant is active on the surface in compression or expansion. The t-test was applied to verify for statistical significance of the results. Comproved with the MST is the same reported in literature, the differences between MWC, CASAC, and CASAE were statistically significant (p<0.001). The MWC, CASAC and CASAE were higher for poractant alfa than for beractant. A higher MWC for poractant alfa means higher elastic recoil of the lung in comparison with beractant. Using a different methodology, our results showed that poractant alfa is probably more effective in a surfactant replacement therapy than beractant due the use of poractant alfa in relation to the use of beractant in preterm infants with Respiratory Distress Syndrome (RDS).

Dissociative electron attachment to chloromethane.

A low energy electron may attach to a molecule, forming a metastable resonance, which may dissociate into a stable anion and a neutral radical. Chloromethane has been a good target for dissociative electron attachment studies, since it is a small molecule with a clear dissociative ‘sigma*’ shape resonance. We present potential energy curves for CH3Cl and its anion, as a function of the C-Cl distance. Due to the resonant nature of the anion, a correct description requires a treatment based on scattering calculations. In order to compute elastic cross sections and phase shifts we employed the Schwinger multichannel method, implemented with pseudopotentials of Bachelet, Hamann and Schlüter, at the static-exchange plus polarization approximation. At the equilibrium geometry, the resonance was found arround 3.3 eV, in accordance to experience. The incoming electron is captured by a ‘sigma*’ orbital located at the C-Cl bond, which will relax in the presence of this extra electron. We took this bond as the reaction coordinate, and performed several scattering calculations for a series of nuclear conformations. The phase shift obtained in each calculation was fitted by a two component function, consisting in the usual Breit-Wigner profile, which captures the resonant character, and a second order polynomial in the wave number, which accounts for the background contribution. That way, we obtained position and width of the resonance, which allowed us to build the potential energy curve. For larger distances, the anion becomes stable and usual electronic structure calculations suffice. Furthermore, the existence of a dipole-bound anion state is revealed when we employed a set of very diffuse functions. The knowledge on the behaviour of the neutral and anionic electronic states helps us in elucidating how the dissociation takes place.

Fourier´s law from a chain of coupled anharmonic oscillators under energy conserving shot noise.

We analyze the transport of heat along a chain of particles interacting through anharmonic potentials consisting of quartic terms in addition to harmonic quadratic terms and subject to heat reservoirs at its ends. Each particle is also subject to an impulsive shot noise with exponentially distributed waiting times whose effect is to change the sign of its velocity, thus conserving the energy of the chain. We show that the introduction of this energy conserving stochastic noise leads to Fourier's law. That is for large system size L the heat current J behaves as J ‘approximately’ 1/L, which amounts to say that the conductivity k is constant. The conductivity is related to the current by J = kΔT/L, where ΔT is the difference in the temperatures of the reservoirs. The behavior of heat conductivity k for small intensities¸ of the shot noise and large system sizes L are obtained by assuming a scaling behavior of the type k = ‘L POT a Psi’(L’lambda POT a/b’) where a and b are scaling exponents. For the pure harmonic case a = b = 1, characterizing a ballistic conduction of heat when the shot noise is absent. For the anharmonic case we found values for the exponents a and b smaller then 1 and thus consistent with a superdiffusive conduction of heat without the shot noise. We also show that the heat conductivity is not constant but is an increasing function of temperature.

Heat conduction in a chain of harmonic and anharmonic oscillators under the presence of an energy conserving noise

Reproducing Fourier's law of heat conduction from a microscopic stochastic model is a long standing challenge in statistical physics. As was shown by Rieder, Lebowitz and Lieb many years ago, a chain of harmonically coupled oscillators connected to two heat baths at different temperatures does not reproduce the diffusive behaviour of Fourier's law, but instead a ballistic one with an infinite thermal conductivity. Since then, there has been a substantial effort from the scientific community in identifying the key mechanism necessary to reproduce such diffusivity, which usually revolved around anharmonicity and the effect of impurities. Recently, it was shown by Dhar, Venkateshan and Lebowitz that Fourier's law can be recovered by introducing an energy conserving noise, whose role is to simulate the elastic collisions between the atoms and other microscopic degrees of freedom, which one would expect to be present in a real solid. For a one-dimensional chain this is accomplished numerically by randomly flipping - under the framework of a Poisson process with a variable “rate of collisions" - the sign of the velocity of an oscillator. In this poster we present Langevin simulations of a one-dimensional chain of oscillators coupled to two heat baths at different temperatures. We consider both harmonic and anharmonic (quartic) interactions, which are studied with and without the energy conserving noise. With these results we are able to map in detail how the heat conductivity k is influenced by both anharmonicity and the energy conserving noise. We also present a detailed analysis of the behaviour of k as a function of the size of the system and the rate of collisions, which includes a finite-size scaling method that enables us to extract the relevant critical exponents. Finally, we show that for harmonic chains, k is independent of temperature, both with and without the noise. Conversely, for anharmonic chains we find that k increases roughly linearly with the temperature of a given reservoir, while keeping the temperature difference fixed.

Irreversibility by competition on a Glauber-Ising model by means of the entropy production.

An out of equilibrium Ising model subjected to an irreversible dynamics is analyzed by means of a stochastic dynamics, on a effort that aims to understand the observed critical behavior as consequence of the intrinsic microscopic characteristics. The study focus on the kinetic phase transitions that take place by assuming a lattice model with inversion symmetry and under the influence of two competing Glauber dynamics, intended to describe the stationary states using the entropy production, which characterize the system behavior and clarifies its reversibility conditions. Thus, it is considered a square lattice formed by two sublattices interconnected, each one of which is in contact with a heat bath at different temperature from the other. Analytical and numerical treatments are faced, using mean-field approximations and Monte Carlo simulations. For the one dimensional model exact results for the entropy production were obtained, though in this case the phase transition that takes place in the two dimensional counterpart is not observed, fact which is in accordance with the behavior shared by lattice models presenting inversion symmetry. Results found for the stationary state show a critical behavior of the same class as the equilibrium Ising model with a phase transition of the second order, which is evidenced by a divergence with an exponent µ ¼ 0:003 of the entropy production derivative.

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 understood. 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 methacholine (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 ¯bers 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, incorporates the phenomenology of both MCh and FA and reproduces experimental results observed with in vitro exposure of smooth muscle to FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells in a tissue level model. The model can also be used in different biological scales.

Model for the Dynamics of the Cytoskeleton.

Particle tracking of microbeads attached to the cytoskeleton (CSK) reveals an intermittent dynamic. The mean squared displacement (MSD) is subdiffusive for small Δt and superdiffusive for large Δt, which are associated with periods of traps and periods of jumps respectively. The analysis of the displacements has shown a non-Gaussian behavior, what is indicative of an active motion, classifying the cells as a far from equilibrium material. Using Langevin dynamics, we reconstruct the dynamic of the CSK. The model is based on the bundles of actin filaments that link themself with the bead RGD coating, trapping it in an harmonic potential. We consider a one- dimensional motion of a particle, neglecting inertial effects (over-damped Langevin dynamics). The resultant force is decomposed in friction force, elastic force and random force, which is used as white noise representing the effect due to molecular agitation. These description until now shows a static situation where the bead performed a random walk in an elastic potential. In order to modeling the active remodeling of the CSK, we vary the equilibrium position of the potential. Inserting a motion in the well center, we change the equilibrium position linearly with time with constant velocity. The result found exhibits a MSD versus time ’tau’ with three regimes. The first regime is when ‘tau’ < ‘tau IND 0’, where ‘tau IND 0’ is the relaxation time, representing the thermal motion. At this regime the particle can diffuse freely. The second regime is a plateau, ‘tau IND 0’ < ‘tau’ < ‘tau IND 1’, representing the particle caged in the potential. Here, ‘tau IND 1’ is a characteristic time that limit the confinement period. And the third regime, ‘tau’ > ‘tau IND 1’, is when the particles are in the superdiffusive behavior. This is where most of the experiments are performed, under 20 frames per second (FPS), thus there is no experimental evidence that support the first regime. We are currently performing experiments with high frequency, up to 100 FPS, attempting to visualize this diffusive behavior. Beside the first regime, our simple model can reproduce MSD curves similar to what has been found experimentally, which can be helpful to understanding CSK structure and properties.

Metáteses árabes da metáfora: desvelar/velar

Um dos mais intrigantes fatos semânticos da língua árabe é a metátese, transposição de fonemas dentro de uma palavra, frequentemente com relação de sentido entre as formas metatéticas. Em m-th-l, a raiz triconsonantal para a palavra “metáfora”, as metáteses são imensamente sugestivas: Th-L-M and L- Th-M; a ambígua função da metáfora: mostrar e ocultar. A partir desse fato, o artigo discute esses aspectos da metáfora

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.