Time-scale and other invariants of integrative mechanical behavior in living cells.

In dealing with systems as complex as the cytoskeleton, we need organizing principles or, short of that, an empirical framework into which these systems fit. We report here unexpected invariants of cytoskeletal behavior that comprise such an empirical framework. We measured elastic and frictional moduli of a variety of cell types over a wide range of time scales and using a variety of biological interventions. In all instances elastic stresses dominated at frequencies below 300 Hz, increased only weakly with frequency, and followed a power law; no characteristic time scale was evident. Frictional stresses paralleled the elastic behavior at frequencies below 10 Hz but approached a Newtonian viscous behavior at higher frequencies. Surprisingly, all data could be collapsed onto master curves, the existence of which implies that elastic and frictional stresses share a common underlying mechanism. Taken together, these findings define an unanticipated integrative framework for studying protein interactions within the complex microenvironment of the cell body, and appear to set limits on what can be predicted about integrated mechanical behavior of the matrix based solely on cytoskeletal constituents considered in isolation. Moreover, these observations are consistent with the hypothesis that the cytoskeleton of the living cell behaves as a soft glassy material, wherein cytoskeletal proteins modulate cell mechanical properties mainly by changing an effective temperature of the cytoskeletal matrix. If so, then the effective temperature becomes an easily quantified determinant of the ability of the cytoskeleton to deform, flow, and reorganize.

Suboptimal filtering and nonlinear time scale transformation for the analysis of multiexponential decays

Multiexponential decays may contain time-constants differing in several orders of magnitudes. In such cases, uniform sampling results in very long records featuring a high degree of oversampling at the final part of the transient. Here, we analyze a nonlinear time scale transformation to reduce the total number of samples with minimum signal distortion, achieving an important reduction of the computational cost of subsequent analyses. We propose a time-varying filter whose length is optimized for minimum mean square error

The problem of the competitiveness of nuclear energy: a biophysical explanation

In this study I try to explain the systemic problem of the low economic competitiveness of nuclear energy for the production of electricity by carrying out a biophysical analysis of its production process. Given the fact that neither econometric approaches nor onedimensional methods of energy analyses are effective, I introduce the concept of biophysical explanation as a quantitative analysis capable of handling the inherent ambiguity associated with the concept of energy. In particular, the quantities of energy, considered as relevant for the assessment, can only be measured and aggregated after having agreed on a pre-analytical definition of a grammar characterizing a given set of finite transformations. Using this grammar it becomes possible to provide a biophysical explanation for the low economic competitiveness of nuclear energy in the production of electricity. When comparing the various unit operations of the process of production of electricity with nuclear energy to the analogous unit operations of the process of production of fossil energy, we see that the various phases of the process are the same. The only difference is related to characteristics of the process associated with the generation of heat which are completely different in the two systems. Since the cost of production of fossil energy provides the base line of economic competitiveness of electricity, the (lack of) economic competitiveness of the production of electricity from nuclear energy can be studied, by comparing the biophysical costs associated with the different unit operations taking place in nuclear and fossil power plants when generating process heat or net electricity. In particular, the analysis focuses on fossil-fuel requirements and labor requirements for those phases that both nuclear plants and fossil energy plants have in common: (i) mining; (ii) refining/enriching; (iii) generating heat/electricity; (iv) handling the pollution/radioactive wastes. By adopting this approach, it becomes possible to explain the systemic low economic competitiveness of nuclear energy in the production of electricity, because of: (i) its dependence on oil, limiting its possible role as a carbon-free alternative; (ii) the choices made in relation to its fuel cycle, especially whether it includes reprocessing operations or not; (iii) the unavoidable uncertainty in the definition of the characteristics of its process; (iv) its large inertia (lack of flexibility) due to issues of time scale; and (v) its low power level.

Characterization of larval moulting cycles in Maja brachydactyla (Brachyura, Majidae) reared in the laboratory

The moulting cycles of all larval instars (zoea I, zoea II, and megalopa) of the spider crab Maja brachydactyla Balss 1922 were studied in laboratory rearing experiments. Morphological changes in the epidermis and cuticle were photographically documented in daily intervals and assigned to successive stages of the moulting cycle (based on Drach's classification system). Our moult-stage characterizations are based on microscopical examination of integumental modifications mainly in the telson, using epidermal condensation, the degree of epidermal retraction (apolysis), and morphogenesis (mainly setagenesis) as criteria. In the zoea II and megalopa, the formation of new setae was also observed in larval appendages including the antenna, maxillule, maxilla, second maxilliped, pleopods, and uropods. As principal stages within the zoea I moulting cycle, we describe postmoult (Drach's stages A–B combined), intermoult (C), and premoult (D), the latter with three substages (D0, D1, and D2). In the zoea II and megalopa, D0 and D1 had to be combined, because morphogenesis (the main characteristic of D1) was unclear in the telson and did not occur synchronically in different appendices. The knowledge of the course and time scale of successive moult-cycle events can be used as a tool for the evaluation of the developmental state within individual larval instars, providing a morphological reference system for physiological and biochemical studies related to crab aquaculture.

Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

Identification of potential small molecule binding pockets on Rho family GTPases

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (.100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.

What can spike train distances tell us about the neural code?

Time scale parametric spike train distances like the Victor and the van Rossum distancesare often applied to study the neural code based on neural stimuli discrimination.Different neural coding hypotheses, such as rate or coincidence coding,can be assessed by combining a time scale parametric spike train distance with aclassifier in order to obtain the optimal discrimination performance. The time scalefor which the responses to different stimuli are distinguished best is assumed to bethe discriminative precision of the neural code. The relevance of temporal codingis evaluated by comparing the optimal discrimination performance with the oneachieved when assuming a rate code.We here characterize the measures quantifying the discrimination performance,the discriminative precision, and the relevance of temporal coding. Furthermore,we evaluate the information these quantities provide about the neural code. Weshow that the discriminative precision is too unspecific to be interpreted in termsof the time scales relevant for encoding. Accordingly, the time scale parametricnature of the distances is mainly an advantage because it allows maximizing thediscrimination performance across a whole set of measures with different sensitivitiesdetermined by the time scale parameter, but not due to the possibility toexamine the temporal properties of the neural code.

Shot noise anomalies in elastic nondegenerate diffusive conductors

We present a theoretical investigation of shot-noise properties in nondegenerate elastic diffusive conductors. Both Monte Carlo simulations and analytical approaches are used. Two interesting phenomena are found: (i) the display of enhanced shot noise for given energy dependences of the scattering time, and (ii) the recovery of full shot noise for asymptotic high applied bias. The first phenomenon is associated with the onset of negative differential conductivity in energy space that drives the system towards a dynamical electrical instability in excellent agreement with analytical predictions. The enhancement is found to be strongly amplified when the dimensionality in momentum space is lowered from three to two dimensions. The second phenomenon is due to the suppression of the effects of long-range Coulomb correlations that takes place when the transit time becomes the shortest time scale in the system, and is common to both elastic and inelastic nondegenerate diffusive conductors. These phenomena shed different light in the understanding of the anomalous behavior of shot noise in mesoscopic conductors, which is a signature of correlations among different current pulses.

Influence of average size and interface passivation on the spectral emission of Si nanocrystals embedded in SiO2

The correlation between the structural (average size and density) and optoelectronic properties [band gap and photoluminescence (PL)] of Si nanocrystals embedded in SiO2 is among the essential factors in understanding their emission mechanism. This correlation has been difficult to establish in the past due to the lack of reliable methods for measuring the size distribution of nanocrystals from electron microscopy, mainly because of the insufficient contrast between Si and SiO2. With this aim, we have recently developed a successful method for imaging Si nanocrystals in SiO2 matrices. This is done by using high-resolution electron microscopy in conjunction with conventional electron microscopy in dark field conditions. Then, by varying the time of annealing in a large time scale we have been able to track the nucleation, pure growth, and ripening stages of the nanocrystal population. The nucleation and pure growth stages are almost completed after a few minutes of annealing time at 1100°C in N2 and afterward the ensemble undergoes an asymptotic ripening process. In contrast, the PL intensity steadily increases and reaches saturation after 3-4 h of annealing at 1100°C. Forming gas postannealing considerably enhances the PL intensity but only for samples annealed previously in less time than that needed for PL saturation. The effects of forming gas are reversible and do not modify the spectral shape of the PL emission. The PL intensity shows at all times an inverse correlation with the amount of Pb paramagnetic centers at the Si-SiO2 nanocrystal-matrix interfaces, which have been measured by electron spin resonance. Consequently, the Pb centers or other centers associated with them are interfacial nonradiative channels for recombination and the emission yield largely depends on the interface passivation. We have correlated as well the average size of the nanocrystals with their optical band gap and PL emission energy. The band gap and emission energy shift to the blue as the nanocrystal size shrinks, in agreement with models based on quantum confinement. As a main result, we have found that the Stokes shift is independent of the average size of nanocrystals and has a constant value of 0.26±0.03 eV, which is almost twice the energy of the Si¿O vibration. This finding suggests that among the possible channels for radiative recombination, the dominant one for Si nanocrystals embedded in SiO2 is a fundamental transition spatially located at the Si¿SiO2 interface with the assistance of a local Si-O vibration.

Magnetostratigraphy of the Miocene continental deposits of the Montes de Castejón (central Ebro basin, Spain): geochronological and paleoenvironmental implications

A detailed magnetostratigraphic study has been carried out in the early to middle Miocene distal alluvial and lacustrine sediments of the Montes de Castejón (central Ebro Basin). The study was based on the analysis of 196 magnetostratigraphic sites sampled along a stratigraphic interval of about 240 meters. Local magnetostratigraphy yielded a sequence of 12 magnetozones (6 normal and 6 reverse) which could be correlated with the Geomagnetic Polarity Time Scale (GPTS) interval C5Cr to C5AD (between 17 and 14.3 Ma.). The sampled sedimentary sequences include the boundary between two tectosedimentary units (TSU, T5 and T6) already defined in the Ebro Basin. The magnetostratigraphy of the Montes de Castejón allows to date the T5/T6 TSU boundary at 16.14 Ma, within chron C5Cn.1n. This magnetostratigraphy also allows us to analyse in detail as well as to discuss the variations in sedimentation rates through space and time between different lacustrine environments: Outer carbonate lacustrine fringes and distal alluvial plains (Montes de Castejón sections) show higher sedimentation rates than offshore lacustrine areas (San Caprasio section, 50 km east of Montes de Castejón).

Slow dynamis in the 3D gonihedric model

We provide analytical evidence of stochastic resonance in polarization switching vertical-cavity surface-emitting lasers (VCSELs). We describe the VCSEL by a two-mode stochastic rate equation model and apply a multiple time-scale analysis. We were able to reduce the dynamical description to a single stochastic differential equation, which is the starting point of the analytical study of stochastic resonance. We confront our results with numerical simulations on the original rate equations, validating the use of a multiple time-scale analysis on stochastic equations as an analytical tool.

Intensity correlation functions of dye lasers: Comparison of colored-gain-noise and colored-loss-noise models

The intensity correlation functions C(t) for the colored-gain-noise model of dye lasers are analyzed and compared with those for the loss-noise model. For correlation times ¿ larger than the deterministic relaxation time td, we show with the use of the adiabatic approximation that C(t) values coincide for both models. For small correlation times we use a method that provides explicit expressions of non-Markovian correlation functions, approximating simultaneously short- and long-time behaviors. Comparison with numerical simulations shows excellent results simultaneously for short- and long-time regimes. It is found that, when the correlation time of the noise increases, differences between the gain- and loss-noise models tend to disappear. The decay of C(t) for both models can be described by a time scale that approaches the deterministic relaxation time. However, in contrast with the loss-noise model, a secondary time scale remains for large times for the gain-noise model, which could allow one to distinguish between both models.

Dark energy equation of state and anthropic selection

We explore the possibility that the dark energy is due to a potential of a scalar field and that the magnitude and the slope of this potential in our part of the Universe are largely determined by anthropic selection effects. We find that, in some models, the most probable values of the slope are very small, implying that the dark energy density stays constant to very high accuracy throughout cosmological evolution. In other models, however, the most probable values of the slope are such that the slow roll condition is only marginally satisfied, leading to a recollapse of the local universe on a time scale comparable to the lifetime of the Sun. In the latter case, the effective equation of state varies appreciably with the redshift, leading to a number of testable predictions.

Simple model for nonexponential relaxation in complex fluids

The nonexponential relaxation occurring in complex dynamics manifested in a wide variety of systems is analyzed through a simple model of diffusion in phase space. It is found that the inability of the system to find its equilibrium state in any time scale becomes apparent in an effective temperature field, which leads to a hierarchy of relaxation times responsible for the slow relaxation phenomena.

Magnetoestratigrafía de las sucesiones del Mioceno medio y superior del Vallés occidental (Depresión del Vallés-Penedés, NE España): Implicaciones biocronológicas y cronoestratigráficas.

The magnetostratigraphic analysis of the middle to late Miocene continental deposits from the Valles-Penedes basin, combined with its well-documented fossil mammal record, provides a well-resoluted chronology for the upper basin infill. It is based on the biostratigraphic and magnetostratigraphic cross-correlation of 18 sections throughout the alluvial and transitional/shallow marine sequences in the Western Valles area. The biostratigraphic framework consists of 24 mammal localities of upper Aragonian and Vallesian age. Correlation of the studied sections to the geomagnetic polarity time scale (GPTS) is based on the distinctive pattern of local magnetozones, as well as the radiometric age of the late Vallesian fauna from the Bicorp Basin (9.6 + 0.3 Ma) and the known relationship of the late Vallesian assemblages with marine beds belonging to the planktonic forarninifera N16 zone. It has led to an absolute dating of the fauna1 events and a precise chronostratigraphy of the Vallesian marnrnal stage in its type area. The Hipparion First Appearance Datum (FAD) defines the lower Vallesian boundary and is dated at 11.1 Ma, at the base of chron C5r. ln. It is in good agreement with radiometric ages from the early Hipparion bearing sites in the Vienna Basin (1 1.1 * 0.5 Ma) and the classic Howenegg locality in Germany (10.8 * 0.3 Ma). It also agrees with the age of the turkish localities of Yailacilar (1 1.6 + 0.25 Ma) and Yenieskihisar-2 (1 1.1 * 0.2 Ma) with absence of Hipparion. Al1 these support the isochrony of the dispersa1 of Hipparion throughout the Mediterranean region. A possible isochrony at a larger geographical scale (Old World, Mesogea) must await more reliable ages of the Hipparion FAD in Asia and Africa. The Cricetulodon FAD that defines the MN9a/MN9b boundary occurs at the middle part of C5n. Assuming an on average constant sedimentation rate, this datum has an age of approximately 10.4 Ma. The earlyllate Vallesian boundary is marked by one of the most distinct fauna1 events of the late Neogene: the dispersa1 of the muridae Progonomys into Europe and North Africa, which coincides with an important macromarnmal turnover. The first extensive appearance of Progonomys in Europe (MN9íMN10 boundary) is dated at 9.7 Ma (C4Ar3r), showing a remarkable diachrony with the Himalayan region. F9i d lly, the FAD of Rotundomys bressnnus occurs in the upper part of C4Ar.ln (9.2-9.3 Ma). The Vallesian spans 2.4 Myr, from 11.1 Ma (CSr.ln) to 8.7 Ma (C4An), and correlates to the early Tortonian.