Modelos mecânicos para o estudo da dinâmica do DNA

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Markovian versus non-Markovian stochastic quantization of a complex-action model

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Analysis of Lattice Size, Energy Density and Denaturation for a One-Dimensional DNA Model

There are several mechanical models to describe the DNA phenomenology. In this work the DNA denaturation is stu- died under thermodynamical and dynamical point of view using the well known Peyrard-Bishop model. The thermody-namics analysis using the transfer integral operator method is briefly reviewed. In particular, the lattice size is discussed and a conjecture about the minimum energy to denaturation is proposed. In terms of the dynamical aspects of the model, the equations of motion for the system are integrated and the results determine the energy density where the denatura- tion occurs. The behavior of the lattice near the phase transition is analyzed. The relation between the thermodynamical and dynamical results is discussed.

Evaluating Measurement Invariance with Censored Ordinal Data: A Monte Carlo Comparison of Alternative Model Estimators and Scales of Measurement

Evaluations of measurement invariance provide essential construct validity evidence. However, the quality of such evidence is partly dependent upon the validity of the resulting statistical conclusions. The presence of Type I or Type II errors can render measurement invariance conclusions meaningless. The purpose of this study was to determine the effects of categorization and censoring on the behavior of the chi-square/likelihood ratio test statistic and two alternative fit indices (CFI and RMSEA) under the context of evaluating measurement invariance. Monte Carlo simulation was used to examine Type I error and power rates for the (a) overall test statistic/fit indices, and (b) change in test statistic/fit indices. Data were generated according to a multiple-group single-factor CFA model across 40 conditions that varied by sample size, strength of item factor loadings, and categorization thresholds. Seven different combinations of model estimators (ML, Yuan-Bentler scaled ML, and WLSMV) and specified measurement scales (continuous, censored, and categorical) were used to analyze each of the simulation conditions. As hypothesized, non-normality increased Type I error rates for the continuous scale of measurement and did not affect error rates for the categorical scale of measurement. Maximum likelihood estimation combined with a categorical scale of measurement resulted in more correct statistical conclusions than the other analysis combinations. For the continuous and censored scales of measurement, the Yuan-Bentler scaled ML resulted in more correct conclusions than normal-theory ML. The censored measurement scale did not offer any advantages over the continuous measurement scale. Comparing across fit statistics and indices, the chi-square-based test statistics were preferred over the alternative fit indices, and ΔRMSEA was preferred over ΔCFI. Results from this study should be used to inform the modeling decisions of applied researchers. However, no single analysis combination can be recommended for all situations. Therefore, it is essential that researchers consider the context and purpose of their analyses.

Air quality simulations for Londrina region, Brazil

The objective of this work were apply and provide a preliminary evaluation of the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) performance, for Londrina region. We performed comparison with measurements obtained in meteorological stations. The model was configured to run with three domains with 27,9 and 3 km of grid resolution, using the ndown program and also was realized a simulation with the model configured to run with a single domain using a land use file based in a classified image for region of MODIS sensor. The emission files to supply the chemistry run were generated based in the work of Martins et al., 2012. RADM2 chemical mechanism and MADE/SORGAM modal aerosol models were used in the simulations. The results demonstrated that model was able to represent coherently the formation and dispersion of the pollution in Metropolitan Region of Londrina and also the importance of using the appropriate land use file for the region.

Two fluid numerical model for coastal applications

Wave breaking is an important coastal process, influencing hydro-morphodynamic processes such as turbulence generation and wave energy dissipation, run-up on the beach and overtopping of coastal defence structures. During breaking, waves are complex mixtures of air and water (“white water”) whose properties affect velocity and pressure fields in the vicinity of the free surface and, depending on the breaker characteristics, different mechanisms for air entrainment are usually observed. Several laboratory experiments have been performed to investigate the role of air bubbles in the wave breaking process (Chanson & Cummings, 1994, among others) and in wave loading on vertical wall (Oumeraci et al., 2001; Peregrine et al., 2006, among others), showing that the air phase is not negligible since the turbulent energy dissipation involves air-water mixture. The recent advancement of numerical models has given valuable insights in the knowledge of wave transformation and interaction with coastal structures. Among these models, some solve the RANS equations coupled with a free-surface tracking algorithm and describe velocity, pressure, turbulence and vorticity fields (Lara et al. 2006 a-b, Clementi et al., 2007). The single-phase numerical model, in which the constitutive equations are solved only for the liquid phase, neglects effects induced by air movement and trapped air bubbles in water. Numerical approximations at the free surface may induce errors in predicting breaking point and wave height and moreover, entrapped air bubbles and water splash in air are not properly represented. The aim of the present thesis is to develop a new two-phase model called COBRAS2 (stands for Cornell Breaking waves And Structures 2 phases), that is the enhancement of the single-phase code COBRAS0, originally developed at Cornell University (Lin & Liu, 1998). In the first part of the work, both fluids are considered as incompressible, while the second part will treat air compressibility modelling. The mathematical formulation and the numerical resolution of the governing equations of COBRAS2 are derived and some model-experiment comparisons are shown. In particular, validation tests are performed in order to prove model stability and accuracy. The simulation of the rising of a large air bubble in an otherwise quiescent water pool reveals the model capability to reproduce the process physics in a realistic way. Analytical solutions for stationary and internal waves are compared with corresponding numerical results, in order to test processes involving wide range of density difference. Waves induced by dam-break in different scenarios (on dry and wet beds, as well as on a ramp) are studied, focusing on the role of air as the medium in which the water wave propagates and on the numerical representation of bubble dynamics. Simulations of solitary and regular waves, characterized by both spilling and plunging breakers, are analyzed with comparisons with experimental data and other numerical model in order to investigate air influence on wave breaking mechanisms and underline model capability and accuracy. Finally, modelling of air compressibility is included in the new developed model and is validated, revealing an accurate reproduction of processes. Some preliminary tests on wave impact on vertical walls are performed: since air flow modelling allows to have a more realistic reproduction of breaking wave propagation, the dependence of wave breaker shapes and aeration characteristics on impact pressure values is studied and, on the basis of a qualitative comparison with experimental observations, the numerical simulations achieve good results.

Thermal excitations of optical nanofibers measured with a fiber-integrated Fabry-Pérot cavity

Efficient coupling of light to quantum emitters, such as atoms, molecules or quantum dots, is one of the great challenges in current research. The interaction can be strongly enhanced by coupling the emitter to the eva-nescent field of subwavelength dielectric waveguides that offer strong lateral confinement of the guided light. In this context subwavelength diameter optical nanofibers as part of a tapered optical fiber (TOF) have proven to be powerful tool which also provide an efficient transfer of the light from the interaction region to an optical bus, that is to say, from the nanofiber to an optical fiber. rnAnother approach towards enhancing light–matter interaction is to employ an optical resonator in which the light is circulating and thus passes the emitters many times. Here, both approaches are combined by experi-mentally realizing a microresonator with an integrated nanofiber waist. This is achieved by building a fiber-integrated Fabry-Pérot type resonator from two fiber Bragg grating mirrors with a stop-band near the cesium D2-line wavelength. The characteristics of this resonator fulfill the requirements of nonlinear optics, optical sensing, and cavity quantum electrodynamics in the strong-coupling regime. Together with its advantageous features, such as a constant high coupling strength over a large volume, tunability, high transmission outside the mirror stop band, and a monolithic design, this resonator is a promising tool for experiments with nanofiber-coupled atomic ensembles in the strong-coupling regime. rnThe resonator's high sensitivity to the optical properties of the nanofiber provides a probe for changes of phys-ical parameters that affect the guided optical mode, e.g., the temperature via the thermo-optic effect of silica. Utilizing this detection scheme, the thermalization dynamics due to far-field heat radiation of a nanofiber is studied over a large temperature range. This investigation provides, for the first time, a measurement of the total radiated power of an object with a diameter smaller than all absorption lengths in the thermal spectrum at the level of a single object of deterministic shape and material. The results show excellent agreement with an ab initio thermodynamic model that considers heat radiation as a volumetric effect and that takes the emitter shape and size relative to the emission wavelength into account. Modeling and investigating the thermalization of microscopic objects with arbitrary shape from first principles is of fundamental interest and has important applications, such as heat management in nano-devices or radiative forcing of aerosols in Earth's climate system. rnUsing a similar method, the effect of the TOF's mechanical modes on the polarization and phase of the fiber-guided light is studied. The measurement results show that in typical TOFs these quantities exhibit high-frequency thermal fluctuations. They originate from high-Q torsional oscillations that couple to the nanofiber-guided light via the strain-optic effect. An ab-initio opto-mechanical model of the TOF is developed that provides an accurate quantitative prediction for the mode spectrum and the mechanically induced polarization and phase fluctuations. These high-frequency fluctuations may limit the ultimate ideality of fiber-coupling into photonic structures. Furthermore, first estimations show that they may currently limit the storage time of nanofiber-based atom traps. The model, on the other hand, provides a method to design TOFs with tailored mechanical properties in order to meet experimental requirements. rn

Protein buffering in model systems and in whole human saliva

The aim of this study was to quantify the buffer attributes (value, power, range and optimum) of two model systems for whole human resting saliva, the purified proteins from whole human resting saliva and single proteins. Two model systems, the first containing amyloglucosidase and lysozyme, and the second containing amyloglucosidase and alpha-amylase, were shown to provide, in combination with hydrogencarbonate and di-hydrogenphosphate, almost identical buffer attributes as whole human resting saliva. It was further demonstrated that changes in the protein concentration as small as 0.1% may change the buffer value of a buffer solution up to 15 times. Additionally, it was shown that there was a protein concentration change in the same range (0.16%) between saliva samples collected at the time periods of 13:00 and others collected at 9:00 am and 17:00. The mode of the protein expression changed between these samples corresponded to the change in basic buffer power and the change of the buffer value at pH 6.7. Finally, SDS Page and Ruthenium II tris (bathophenantroline disulfonate) staining unveiled a constant protein expression in all samples except for one 50 kDa protein band. As the change in the expression pattern of that 50 kDa protein band corresponded to the change in basic buffer power and the buffer value at pH 6.7, it was reasonable to conclude that this 50 kDa protein band may contain the protein(s) belonging to the protein buffer system of human saliva.

Mechanical modeling of Intraneural Ganglion Cyst

Intraneural Ganglion Cyst is a 200 year old mystery related to nerve injury which is yet to be solved. Current treatments for the above problem are relatively simple procedures related to removal of cystic contents from the nerve. However, these treatments may result into neuropathic pain and recurrence of the cyst. The articular theory proposed by Spinner et al., (Spinner et al. 2003) takes into consideration the neurological deficit in Common Peroneal Nerve (CPN) branch of the sciatic nerve and affirms that in addition to the above treatments, ligation of articular branch results into foolproof eradication of the deficit. Mechanical Modeling of the Affected Nerve Cross Section will reinforce the articular theory (Spinner et al. 2003). As the cyst propagates, it compresses the neighboring fascicles and the nerve cross section appears like a signet ring. Hence, in order to mechanically model the affected nerve cross section; computational methods capable of modeling excessively large deformations are required. Traditional FEM produces distorted elements while modeling such deformations, resulting into inaccuracies and premature termination of the analysis. The methods described in this Master’s Thesis are effective enough to be able to simulate such deformations. The results obtained from the model adequately resemble the MRI image obtained at the same location and shows an appearance of a signet ring. This Master’s Thesis describes the neurological deficit in brief followed by detail explanation of the advanced computational methods used to simulate this problem. Finally, qualitative results show the resemblance of mechanical model to MRI images of the Nerve Cross Section at the same location validating the capability of these methods to study this neurological deficit.

The pollen tube: a soft shell with a hard core

Plant cell expansion is controlled by a fine-tuned balance between intracellular turgor pressure, cell wall loosening and cell wall biosynthesis. To understand these processes, it is important to gain in-depth knowledge of cell wall mechanics. Pollen tubes are tip-growing cells that provide an ideal system to study mechanical properties at the single cell level. With the available approaches it was not easy to measure important mechanical parameters of pollen tubes, such as the elasticity of the cell wall. We used a cellular force microscope (CFM) to measure the apparent stiffness of lily pollen tubes. In combination with a mechanical model based on the finite element method (FEM), this allowed us to calculate turgor pressure and cell wall elasticity, which we found to be around 0.3 MPa and 20–90 MPa, respectively. Furthermore, and in contrast to previous reports, we showed that the difference in stiffness between the pollen tube tip and the shank can be explained solely by the geometry of the pollen tube. CFM, in combination with an FEM-based model, provides a powerful method to evaluate important mechanical parameters of single, growing cells. Our findings indicate that the cell wall of growing pollen tubes has mechanical properties similar to rubber. This suggests that a fully turgid pollen tube is a relatively stiff, yet flexible cell that can react very quickly to obstacles or attractants by adjusting the direction of growth on its way through the female transmitting tissue.

Neotethys closure history of Anatolia: insights from 40Ar–39Ar geochronology and P–T estimation in high pressure metasedimentary rocks

The multiple high-pressure (HP), low-temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction- and continental accretion-related evolution of the eastern limb of the long-lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe-Mg-carpholite in three metasedimentary units of the Gondwana-derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single-continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica 40Ar–39Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi-equilibrium thermobarometry), on carpholite-bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite-bearing assemblages were retrogressed through greenschist-facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post-collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tavşanlı Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian-type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll-back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.

Landmark Detection for Fusion of Fundus and MRI Towards a Patient-Specific Multi-Modal Eye Model

Ophthalmologists typically acquire different image modalities to diagnose eye pathologies. They comprise e.g., Fundus photography, Optical Coherence Tomography (OCT), Computed Tomography (CT) and Magnetic Resonance Imaging (MRI). Yet, these images are often complementary and do express the same pathologies in a different way. Some pathologies are only visible in a particular modality. Thus, it is beneficial for the ophthalmologist to have these modalities fused into a single patient-specific model. The presented article’s goal is a fusion of Fundus photography with segmented MRI volumes. This adds information to MRI which was not visible before like vessels and the macula. This article’s contributions include automatic detection of the optic disc, the fovea, the optic axis and an automatic segmentation of the vitreous humor of the eye.

Measuring the mechanical properties of plant cells by combining micro-indentation with osmotic treatments

Growth in plants results from the interaction between genetic and signalling networks and the mechanical properties of cells and tissues. There has been a recent resurgence in research directed at understanding the mechanical aspects of growth, and their feedback on genetic regulation. This has been driven in part by the development of new micro-indentation techniques to measure the mechanical properties of plant cells in vivo. However, the interpretation of indentation experiments remains a challenge, since the force measures results from a combination of turgor pressure, cell wall stiffness, and cell and indenter geometry. In order to interpret the measurements, an accurate mechanical model of the experiment is required. Here, we used a plant cell system with a simple geometry, Nicotiana tabacum Bright Yellow-2 (BY-2) cells, to examine the sensitivity of micro-indentation to a variety of mechanical and experimental parameters. Using a finite-element mechanical model, we found that, for indentations of a few microns on turgid cells, the measurements were mostly sensitive to turgor pressure and the radius of the cell, and not to the exact indenter shape or elastic properties of the cell wall. By complementing indentation experiments with osmotic experiments to measure the elastic strain in turgid cells, we could fit the model to both turgor pressure and cell wall elasticity. This allowed us to interpret apparent stiffness values in terms of meaningful physical parameters that are relevant for morphogenesis.

Fluid-structure interaction of an elastic pipe immersed inside a coaxial rigid pipe: a model for the Tullio phenomenon

An axisymmetric, elastic pipe is filled with an incompressible fluid and is immersed in a second, coaxial rigid pipe which contains the same fluid. A pressure pulse in the outer fluid annulus deforms the elastic pipe which invokes a fluid motion in the fluid core. It is the aim of this study to investigate streaming phenomena in the core which may originate from such a fluid-structure interaction. This work presents a numerical solver for such a configuration. It was developed in the OpenFOAM software environment and is based on the Arbitrary Lagrangian Eulerian (ALE) approach for moving meshes. The solver features a monolithic integration of the one-dimensional, coupled system between the elastic structure and the outer fluid annulus into a dynamic boundary condition for the moving surface of the fluid core. Results indicate that our configuration may serve as a mechanical model of the Tullio Phenomenon (sound-induced vertigo).

Genetic and functional analyses of cell division proteins FtsZ and FtsA in Escherichia coli

In the current model for bacterial cell division, the FtsZ protein forms a ring that marks the division plane, creating a cytoskeletal framework for the subsequent action of other essential division proteins such as FtsA and ZipA. The putative protein complex ultimately generates the division septum. The essential cell division protein FtsZ is a functional and structural homolog of eukaryotic tubulin, and like tubulin, FtsZ hydrolyzes GTP and self-assembles into protein filaments in a strictly GTP-dependent manner. FtsA shares sequence similarity with members of the ATPase superfamily that include actin, but its actual function remains unknown. To test the division model and elucidate functions of the division proteins, this dissertation primarily focuses on the analysis of FtsZ and FtsA in Escherichia coli. ^ By tagging with green fluorescent protein, we first demonstrated that FtsA also exhibits a ring-like structure at the potential division site. The localization of FtsA was dependent on functional FtsZ, suggesting that FtsA is recruited to the septum by the FtsZ ring. In support of this idea, we showed that FtsA and FtsZ directly interact. Using a novel E. coli in situ assay, we found that the FtsA-FtsZ interaction appears to be species-specific, although an interspecies interaction could occur between FtsA and FtsZ proteins from two closely related organisms. In addition, mutagenesis of FtsA revealed that no single domain is solely responsible for its septal localization or interaction with FtsZ. To explore the function of FtsA, we purified FtsA protein and demonstrated that it has ATPase activity. Furthermore, purified FtsA stimulates disassembly of FtsZ polymers in a sedimentation assay but does not affect GTP hydrolysis of FtsZ. This result suggests that in the cell, FtsA may function similarly in regulating dynamic instability of the FtsZ ring during the cell division process. ^ To elucidate the structure-function relationship of FtsZ, we carried out thorough genetic and functional analyses of the mutagenized FtsZ derivatives. Our results indicate that the conserved N-terminal domain of FtsZ is necessary and sufficient for FtsZ self-assembly and localization. Moreover, we discovered a critical role for an extreme C-terminal domain of FtsZ that consists of only 12 residues. Truncated FtsZ derivatives lacking this domain, though able to polymerize and localize, are defective in ring formation in vivo as well as interaction with FtsA and ZipA. Alanine scanning mutagenesis of this region pinpointed at least five residues necessary for the function of FtsZ. Studies of protein levels and protein-protein interactions suggested that these residues may be involved in regulating protein stability and/or FtsZ-FtsA interactions. Interestingly, two of the point mutants exhibited dominant-negative phenotypes. ^ In summary, results from this thesis work have provided additional support for the division machinery model and will contribute to a better understanding of the coordinate functions of FtsA and FtsZ in the cell division process. ^