A model for individual and collective cell movement in Dictyostelium discoideum

The cellular slime mold Dictyostelium discoideum is a widely used model system for studying a variety of basic processes in development, including cell–cell signaling, signal transduction, pattern formation, cell motility, and the movement of tissue-like aggregates of cells. Many aspects of cell motion are poorly understood, including how individual cell behavior produces the collective motion of cells observed within the mound and slug. Herein, we describe a biologically realistic model for motile D. discoideum cells that can generate active forces, that interact via surface molecules, and that can detect and respond to chemotactic signals. We model the cells as deformable viscoelastic ellipsoids and incorporate signal transduction and cell–cell signaling by using a previously developed model. The shape constraint restricts the admissible deformations but makes the simulation of a large number of interacting cells feasible. Because the model is based on known processes, the parameters can be estimated or measured experimentally. We show that this model can reproduce the observations on the chemotactic behavior of single cells, streaming during aggregation, and the collective motion of an aggregate of cells driven by a small group of pacemakers. The model predicts that the motion of two-dimensional slugs [Bonner, J. T. (1998) Proc. Natl. Acad. Sci. USA 95, 9355–9359] results from the same behaviors that are exhibited by individual cells; it is not necessary to invoke different mechanisms or behaviors. Our computational experiments also suggest previously uncharacterized phenomena that may be experimentally observable.

The RESID Database of protein structure modifications and the NRL-3D Sequence–Structure Database

The RESID Database is a comprehensive collection of annotations and structures for protein post-translational modifications including N-terminal, C-terminal and peptide chain cross-link modifications. The RESID Database includes systematic and frequently observed alternate names, Chemical Abstracts Service registry numbers, atomic formulas and weights, enzyme activities, taxonomic range, keywords, literature citations with database cross-references, structural diagrams and molecular models. The NRL-3D Sequence–Structure Database is derived from the three-dimensional structure of proteins deposited with the Research Collaboratory for Structural Bioinformatics Protein Data Bank. The NRL-3D Database includes standardized and frequently observed alternate names, sources, keywords, literature citations, experimental conditions and searchable sequences from model coordinates. These databases are freely accessible through the National Cancer Institute–Frederick Advanced Biomedical Computing Center at these web sites: http://www.ncifcrf.gov/RESID, http://www.ncifcrf.gov/ NRL-3D; or at these National Biomedical Research Foundation Protein Information Resource web sites: http://pir.georgetown.edu/pirwww/dbinfo/resid.html, http://pir.georgetown.edu/pirwww/dbinfo/nrl3d.html

Active control of waves in a cochlear model with subpartitions.

Multiscale asymptotic methods developed previously to study macromechanical wave propagation in cochlear models are generalized here to include active control of a cochlear partition having three subpartitions, the basilar membrane, the reticular lamina, and the tectorial membrane. Activation of outer hair cells by stereocilia displacement and/or by lateral wall stretching result in a frequency-dependent force acting between the reticular lamina and basilar membrane. Wavelength-dependent fluid loads are estimated by using the unsteady Stokes' equations, except in the narrow gap between the tectorial membrane and reticular lamina, where lubrication theory is appropriate. The local wavenumber and subpartition amplitude ratios are determined from the zeroth order equations of motion. A solvability relation for the first order equations of motion determines the subpartition amplitudes. The main findings are as follows: The reticular lamina and tectorial membrane move in unison with essentially no squeezing of the gap; an active force level consistent with measurements on isolated outer hair cells can provide a 35-dB amplification and sharpening of subpartition waveforms by delaying dissipation and allowing a greater structural resonance to occur before the wave is cut off; however, previously postulated activity mechanisms for single partition models cannot achieve sharp enough tuning in subpartitioned models.

Comparison of the defensive elicitation activity of cerato-populin and cerato-platanin: a structural model

Plants can defend themselves from potential pathogenic microorganisms relying on a complex interplay of signaling pathways: activation of the MAPK cascade, transcription of defense related genes, production of reactive oxygen species, nitric oxide and synthesis of other defensive compounds such as phytoalexins. These events are triggered by the recognition of pathogen’s effectors (effector-triggered immunity) or PAMPs (PAMP-triggered immunity). The Cerato Platanin Family (CPF) members are Cys-rich proteins secreted and localized on fungal cell walls, involved in several aspects of fungal development and pathogen-host interactions. Although more than hundred genes of the CPF have been identified and analyzed, the structural and functional characterization of the expressed proteins has been restricted only to few members of the family. Interestingly, those proteins have been shown to bind chitin with diverse affinity and after foliar treatment they elicit defensive mechanisms in host and non-host plants. This property turns cerato platanins into interesting candidates, worth to be studied to develop new fungal elicitors with applications in sustainable agriculture. This study focus on cerato-platanin (CP), core member of the family and on the orthologous cerato-populin (Pop1). The latter shows an identity of 62% and an overall homology of 73% with respect to CP. Both proteins are able to induce MAPKs phosphorylation, production of reactive oxygen species and nitric oxide, overexpression of defense’s related genes, programmed cell death and synthesis of phytoalexins. CP, however, when compared to Pop1, induces a faster response and, in some cases, a stronger activity on plane leaves. Aim of the present research is to verify if the dissimilarities observed in the defense elicitation activity of these proteins can be associated to their structural and dynamic features. Taking advantage of the available CP NMR structure, Pop1’s 3D one was obtained by homology modeling. Experimental residual dipolar couplings and 1H, 15N, 13C resonance assignments were used to validate the model. Previous works on CPF members, addressed the highly conserved random coil regions (loops b1-b2 and b2-b3) as sufficient and necessary to induce necrosis in plants’ leaves: that region was investigated in both Pop1 and CP. In the two proteins the loops differ, in their primary sequence, for few mutations and an insertion with a consequent diversification of the proteins’ electrostatic surface. A set of 2D and 3D NMR experiments was performed to characterize both the spatial arrangement and the dynamic features of the loops. NOE data revealed a more extended network of interactions between the loops in Pop1 than in CP. In addition, in Pop1 we identified a salt bridge Lys25/Asp52 and a strong hydrophobic interaction between Phe26/Trp53. These structural features were expected not only to affect the loops’ spatial arrangement, but also to reduce the degree of their conformational freedom. Relaxation data and the order parameter S2 indeed highlighted reduced flexibility, in particular for loop b1-b2 of Pop1. In vitro NMR experiments, where Pop1 and CP were titrated with oligosaccharides, supported the hypothesis that the loops structural and dynamic differences may be responsible for the different chitin-binding properties of the two proteins: CP selectively binds tetramers of chitin in a shallow groove on one side of the barrel defined by loops b1-b2, b2-b3 and b4-b5, Pop1, instead, interacts in a non-specific fashion with oligosaccharides. Because the region involved in chitin-binding is also responsible for the defense elicitation activity, possibly being recognized by plant's receptors, it is reasonable to expect that those structural and dynamic modifications may also justify the different extent of defense elicitation. To test that hypothesis, the initial steps of a protocol aimed to the identify a receptor for CP, in silico, are presented.

Evidence for reciprocal antagonism between motion sensors tuned to coarse and fine features

Early visual processing analyses fine and coarse image features separately. Here we show that motion signals derived from fine and coarse analyses are combined in rather a surprising way: Coarse and fine motion sensors representing the same direction of motion inhibit one another and an imbalance can reverse the motion perceived. Observers judged the direction of motion of patches of filtered two-dimensional noise, centered on 1 and 3 cycles/deg. When both sets of noise were present and only the 3 cycles/deg noise moved, judgments were reversed at short durations. When both sets of noise moved, judgments were correct but sensitivity was impaired. Reversals and impairments occurred both with isotropic noise and with orientation-filtered noise. The reversals and impairments could be simulated in a model of motion sensing by adding a stage in which the outputs of motion sensors tuned to 1 and 3 cycles/deg and the same direction of motion were subtracted from one another. The subtraction model predicted and we confirmed in experiments with orientation-filtered noise that if the 1 cycle/deg noise flickered and the 3 cycles/deg noise moved, the 1 cycle/deg noise appeared to move in the opposite direction to the 3 cycles/deg noise even at long durations.

Well posedness of an integrodifferential kinetic model of Fokker-Planck type for angiogenesis.

Tumor induced angiogenesis processes including the effect of stochastic motion and branching of blood vessels can be described coupling a (nonlocal in time) integrodifferential kinetic equation of Fokker–Planck type with a diffusion equation for the tumor induced ingiogenic factor. The chemotactic force field depends on the flux of blood vessels through the angiogenic factor. We develop an existence and uniqueness theory for this system under natural assumptions on the initial data. The proof combines the construction of fundamental solutions for associated linearized problems with comparison principles, sharp estimates of the velocity integrals and compactness results for this type of kinetic and parabolic operators

Avaliação em 3D da cabeça da mandíbula em crianças portadoras de mordida cruzada posterior unilateral por meio de tomografia computadorizada de feixe cônico

Durante o desenvolvimento da oclusão, a instalação de maloclusões podem resultar em desarmonias dento faciais de natureza e severidade diversas, podendo provocar alterações no desenvolvimento crânio facial, dentre as estruturas envolvidas as Articulações Temporo Mandibulares (ATM), podem sofrer alguma influência, dessa forma a avaliação desta região, no aspecto morfológico e funcional, constituí tema de interesse, sempre que levados em conta os aspectos funcionais da oclusão. A relação entre a forma e a função, tanto das cabeças da mandíbula, bem como o contorno da fossa mandibular com as maloclusões ainda é controversa e não está compreendida por completo, porém a literatura sobre o assunto, demonstra correlação entre a instalação de maloclusões e modificações neste sistema, mesmo que algumas alterações não sejam de ordem estatística e em amostras de indivíduos em tenra idade, as mesmas podem comprometer o desenvolvimento adequado em indivíduos adultos ou mesmo adultos jovens. Tendo como propósito nesse estudo a avaliação das cabeças da mandíbula quanto ao volume e superficíe dos lados direito e esquerdo, cruzado e não cruzado, a amostra selecionada foi de 20 indivíduos com mordida cruzada posterior unilateral, com idades entre 06 e 09 anos de idade, utilizando imagens de tomografia computadorizada por feixe cônico, imagens obtidas por um equipamento modelo i- Cat, sendo utilizado na reformatação e manipulação das imagens o programa computacional - NemoCeph 3D® versão 11.5. Nas medições propostas para esse estudo, utilizou-se o teste t pareado de Student para amostras com distribuição normal. Na observação das tabelas e seus respectivos gráficos, podemos verificar que na comparação entre os lados direito e esquerdo, e cruzado e não cruzado das cabeças da mandíbula, com relação ao volume e superfície, existem diferenças numéricas entre elas, porém não pode ser observado diferenças estatísticas significantes, nessa amostra especifica com a metodologia empregada para esse estudo. Assim foi possível concluir que nas Mordidas Cruzadas Posteriores Unilaterais as cabeças da mandíbula tanto em seu volume como em sua superfície não apresentaram diferenças estatisticamente significantes na amostra estudada.

Finite temperature dynamics of vortices in the two dimensional anisotropic Heisenberg model

We study the effects of finite temperature on the dynamics of non-planar vortices in the classical, two-dimensional anisotropic Heisenberg model with XY- or easy-plane symmetry. To this end, we analyze a generalized Landau-Lifshitz equation including additive white noise and Gilbert damping. Using a collective variable theory with no adjustable parameters we derive an equation of motion for the vortices with stochastic forces which are shown to represent white noise with an effective diffusion constant linearly dependent on temperature. We solve these stochastic equations of motion by means of a Green's function formalism and obtain the mean vortex trajectory and its variance. We find a non-standard time dependence for the variance of the components perpendicular to the driving force. We compare the analytical results with Langevin dynamics simulations and find a good agreement up to temperatures of the order of 25% of the Kosterlitz-Thouless transition temperature. Finally, we discuss the reasons why our approach is not appropriate for higher temperatures as well as the discreteness effects observed in the numerical simulations.

Condensados de Bose-Einstein em redes óticas: a transição superfluido-isolante de Mott em redes hexagonais e a classe de universalidade superfluido-vidro de Bose em 3D

Estudamos transições de fases quânticas em gases bosônicos ultrafrios aprisionados em redes óticas. A física desses sistemas é capturada por um modelo do tipo Bose-Hubbard que, no caso de um sistema sem desordem, em que os átomos têm interação de curto alcance e o tunelamento é apenas entre sítios primeiros vizinhos, prevê a transição de fases quântica superfluido-isolante de Mott (SF-MI) quando a profundidade do potencial da rede ótica é variado. Num primeiro estudo, verificamos como o diagrama de fases dessa transição muda quando passamos de uma rede quadrada para uma hexagonal. Num segundo, investigamos como a desordem modifica essa transição. No estudo com rede hexagonal, apresentamos o diagrama de fases da transição SF-MI e uma estimativa para o ponto crítico do primeiro lobo de Mott. Esses resultados foram obtidos usando o algoritmo de Monte Carlo quântico denominado Worm. Comparamos nossos resultados com os obtidos a partir de uma aproximação de campo médio e com os de um sistema com uma rede ótica quadrada. Ao introduzir desordem no sistema, uma nova fase emerge no diagrama de fases do estado fundamental intermediando a fase superfluida e a isolante de Mott. Essa nova fase é conhecida como vidro de Bose (BG) e a transição de fases quântica SF-BG que ocorre nesse sistema gerou muitas controvérsias desde seus primeiros estudos iniciados no fim dos anos 80. Apesar dos avanços em direção ao entendimento completo desta transição, a caracterização básica das suas propriedades críticas ainda é debatida. O que motivou nosso estudo, foi a publicação de resultados experimentais e numéricos em sistemas tridimensionais [Yu et al. Nature 489, 379 (2012), Yu et al. PRB 86, 134421 (2012)] que violam a lei de escala $\\phi= u z$, em que $\\phi$ é o expoente da temperatura crítica, $z$ é o expoente crítico dinâmico e $ u$ é o expoente do comprimento de correlação. Abordamos essa controvérsia numericamente fazendo uma análise de escalonamento finito usando o algoritmo Worm nas suas versões quântica e clássica. Nossos resultados demonstram que trabalhos anteriores sobre a dependência da temperatura de transição superfluido-líquido normal com o potencial químico (ou campo magnético, em sistemas de spin), $T_c \\propto (\\mu-\\mu_c)^\\phi$, estavam equivocados na interpretação de um comportamento transiente na aproximação da região crítica genuína. Quando os parâmetros do modelo são modificados de maneira a ampliar a região crítica quântica, simulações com ambos os modelos clássico e quântico revelam que a lei de escala $\\phi= u z$ [com $\\phi=2.7(2)$, $z=3$ e $ u = 0.88(5)$] é válida. Também estimamos o expoente crítico do parâmetro de ordem, encontrando $\\beta=1.5(2)$.

Human Motion Analysis Based on Sequential Modeling of Radar Signal and Stereo Image Features

Falls are one of the greatest threats to elderly health in their daily living routines and activities. Therefore, it is very important to detect falls of an elderly in a timely and accurate manner, so that immediate response and proper care can be provided, by sending fall alarms to caregivers. Radar is an effective non-intrusive sensing modality which is well suited for this purpose, which can detect human motions in all types of environments, penetrate walls and fabrics, preserve privacy, and is insensitive to lighting conditions. Micro-Doppler features are utilized in radar signal corresponding to human body motions and gait to detect falls using a narrowband pulse-Doppler radar. Human motions cause time-varying Doppler signatures, which are analyzed using time-frequency representations and matching pursuit decomposition (MPD) for feature extraction and fall detection. The extracted features include MPD features and the principal components of the time-frequency signal representations. To analyze the sequential characteristics of typical falls, the extracted features are used for training and testing hidden Markov models (HMM) in different falling scenarios. Experimental results demonstrate that the proposed algorithm and method achieve fast and accurate fall detections. The risk of falls increases sharply when the elderly or patients try to exit beds. Thus, if a bed exit can be detected at an early stage of this motion, the related injuries can be prevented with a high probability. To detect bed exit for fall prevention, the trajectory of head movements is used for recognize such human motion. A head detector is trained using the histogram of oriented gradient (HOG) features of the head and shoulder areas from recorded bed exit images. A data association algorithm is applied on the head detection results to eliminate head detection false alarms. Then the three dimensional (3D) head trajectories are constructed by matching scale-invariant feature transform (SIFT) keypoints in the detected head areas from both the left and right stereo images. The extracted 3D head trajectories are used for training and testing an HMM based classifier for recognizing bed exit activities. The results of the classifier are presented and discussed in the thesis, which demonstrates the effectiveness of the proposed stereo vision based bed exit detection approach.

First Order Phase Transition in a 3D disordered system

We present a detailed numerical study on the effects of adding quenched impurities to a three dimensional system which in the pure case undergoes a strong first order phase transition (specifically, the ferromagnetic/paramagnetic transition of the site-diluted four states Potts model). We can state that the transition remains first-order in the presence of quenched disorder (a small amount of it) but it turns out to be second order as more impurities are added. A tricritical point, which is studied by means of Finite-Size Scaling, separates the first-order and second-order parts of the critical line. The results were made possible by a new definition of the disorder average that avoids the diverging-variance probability distributions that arise using the standard methodology. We also made use of a recently proposed microcanonical Monte Carlo method in which entropy, instead of free energy, is the basic quantity.

Diseño de un modelo geológico-geotécnico 3D de la Vega Baja del Río Segura (Alicante, SE España)

Los modelos geológico-geotécnicos permiten al ingeniero comprender mejor las condiciones reinantes en un determinado lugar, además de identificar los principales problemas geotécnicos y hacer más realista la estimación de propiedades del suelo. En este trabajo se presenta la metodología empleada para el diseño de un modelo geológico-geotécnico tridimensional de la Vega Baja del Río Segura que consta de cuatro zonas caracterizadas por sus propiedades geotécnicas y su problemática asociada. El modelo resulta fundamentalmente de gran utilidad para la planificación de investigaciones preliminares de obras civiles.

Detección automática con snakes y representación 3D sobre imágenes cerebrales

Comunicación presentada en la VI Conferencia de la Asociación Española para la Inteligencia Artificial (CAEPIA'95), Alicante, 15-17 noviembre 1995.

RASMA: plataforma para la creación de animaciones Stop-motion con asistencia robótica

En este artículo se describe el concepto de plataforma RASMA, Robot-Assisted Stop-Motion Animation, cuya finalidd es facilitar la tarea de generar los fotogramas necesarios para crear una secuencia animada en 2D. Se describen tanto la generación de trayectorias que deben seguir los objetos (en Unity 3D o en Adobe Flash Player), como la exportación/importación de los ficheros de datos en XML, la planificación de las trayectorias del robot, la toma de fotogramas y el ensamblado final de toda la secuencia.

GPGPU implementation of growing neural gas: application to 3D scene reconstruction

Self-organising neural models have the ability to provide a good representation of the input space. In particular the Growing Neural Gas (GNG) is a suitable model because of its flexibility, rapid adaptation and excellent quality of representation. However, this type of learning is time-consuming, especially for high-dimensional input data. Since real applications often work under time constraints, it is necessary to adapt the learning process in order to complete it in a predefined time. This paper proposes a Graphics Processing Unit (GPU) parallel implementation of the GNG with Compute Unified Device Architecture (CUDA). In contrast to existing algorithms, the proposed GPU implementation allows the acceleration of the learning process keeping a good quality of representation. Comparative experiments using iterative, parallel and hybrid implementations are carried out to demonstrate the effectiveness of CUDA implementation. The results show that GNG learning with the proposed implementation achieves a speed-up of 6× compared with the single-threaded CPU implementation. GPU implementation has also been applied to a real application with time constraints: acceleration of 3D scene reconstruction for egomotion, in order to validate the proposal.