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.

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.

The role of soil in vegetated gravelly river braid plains: more than just a passive response?

This paper reviews the role of alluvial soils in vegetated gravelly river braid plains. When considering decadal time scales of river evolution, we argue that it becomes vital to consider soil development as an emergent property of the developing ecosystem. Soil processes have been relatively overlooked in accounts of the interactions between braided river processes and vegetation, although soils have been observed on vegetated fluvial landforms. We hypothesise that soil development plays a major role in the transition (speed and pathway) from a fresh sediment deposit to a vegetated soil-covered landform. Disturbance (erosion and/or deposition), vertical sediment structure (process history), vegetation succession, biological activity and water table fluctuation are seen as the main controls on early alluvial soil evolution. Erosion and deposition processes may not only act as soil disturbing agents, but also as suppliers of ecosystem resources, because of their role in delivering and changing access (e.g. through avulsion) to fluxes of water, fine sediments and organic matter. In turn, the associated initial ecosystem may influence further fluvial landform development, such as through the trapping of fine-grained sediments (e.g. sand) by the engineering action of vegetation and the deposit stabilisation by the developing above and belowground biomass. This may create a strong feedback between geomorphological processes, vegetation succession and soil evolution which we summarise in a conceptual model. We illustrate this model by an example from the Allondon River (CH) and identify the research questions that follow.

Relación entre anchoveta y ambiente a diferentes escalas temporales

Relación entre anchoveta y ambiente a diferentes escalas temporales. Bol Inst Mar Perú 25(1-2):13-21.- Se estudió la relación del ambiente y la anchoveta peruana (Engraulis ringens) en el Ecosistema de Afloramiento Peruano (EAP). Se hipotetiza que en la zona de afloramiento comprendida entre 5°S - 13°S se presenta una relación negativa entre la temperatura, concentración de oxígeno y desembarques de anchoveta en escalas de tiempo decadales. En series de tiempo mensuales observadas entre 1950 – 2008, se analizó la temperatura superficial, subsuperficial, oxígeno, clorofila-a y desembarques de anchoveta y sardina, aplicando métodos estadísticos y espectrales para obtener modos temporales decadales y caracterizar la modulación a baja frecuencia de sus ciclos estacionales. También se usaron datos de reanálisis para caracterizar los cambios estacionales a partir de los forzamientos (remoto y local) ambientales del EAP. Se encontró una fuerte asociación a escala decadal entre las series ambientales costeras y los desembarques de anchoveta y se sugiere que la covariabilidad ambiental en escalas temporales interanuales, intraestacionales, decadales, seculares frente al EAP se incrementará durante la próxima década.

Immunohistochemical and flow cytometric analysis of long-term label-retaining cells in the adult heart.

Cardiac-resident stem/progenitor cells have been identified based on expression of stem cell-associated antigens. However, no single surface marker allows to identify a definite cardiac stem/progenitor cell entity. Hence, functional stem cell markers have been extensively searched for. In homeostatic systems, stem cells divide infrequently and therefore retain DNA labels such as 5-bromo-2'-deoxyuridine, which are diluted with division. We used this method to analyze long-term label-retaining cells in the mouse heart after 14 days of 5-bromo-2'-deoxyuridine administration. Labeled cells were detected using immunohistochemical and flow-cytometric methods after varying chasing periods up to 12 months. Using mathematical models, the observed label dilution could consistently be described in the context of a 2-population model, whereby a population of rapidly dividing cells accounted for an accelerated early decline, and a population of slowly dividing cells accounted for decelerated dilution on longer time scales. Label-retaining cells were preferentially localized in the atria and apical region and stained negative for markers of the major cell lineages present in the heart. Most cells with long-term label-retention expressed stem cell antigen-1 (Sca-1). Sca-1(+)CD31(-) cells formed cell aggregates in culture, out of which lineage-negative (Lin(-))Sca-1(+)CD31(-) cells emerged, which could be cultured for many passages. These cells formed cardiospheres and showed differentiation potential into mesenchymal cell lineages. When cultured in cardiomyogenic differentiation medium, they expressed cardiac-specific genes. In conclusion, recognition of slow-cycling cells provides functional evidence of stem/progenitor cells in the heart. Lin(-)Sca-1(+)CD31(-) cardiac-derived progenitors have a potential for differentiation into cardiomyogenic and mesenchymal cell lineages.

Scaling the microrheology of living cells

We report a scaling law that governs both the elastic and frictional properties of a wide variety of living cell types, over a wide range of time scales and under a variety of biological interventions. This scaling identifies these cells as soft glassy materials existing close to a glass transition, and implies that cytoskeletal proteins may regulate cell mechanical properties mainly by modulating the effective noise temperature of the matrix. The practical implications are that the effective noise temperature is an easily quantified measure of the ability of the cytoskeleton to deform, flow, and reorganize.

Anomalous crossover between thermal and shot noise in macroscopic diffusive conductors

We predict the existence of an anomalous crossover between thermal and shot noise in macroscopic diffusive conductors. We first show that, besides thermal noise, these systems may also exhibit shot noise due to fluctuations of the total number of carriers in the system. Then we show that at increasing currents the crossover between the two noise behaviors is anomalous, in the sense that the low-frequency current spectral density displays a region with a superlinear dependence on the current up to a cubic law. The anomaly is due to the nontrivial coupling in the presence of the long-range Coulomb interaction among the three time scales relevant to the phenomenon, namely, diffusion, transit, and dielectric relaxation time.

Driving rate effects in avalanche-mediated first-order phase transitions

We study the driving-rate and temperature dependence of the power-law exponents that characterize the avalanche distribution in first-order phase transitions. Measurements of acoustic emission in structural transitions in Cu-Zn-Al and Cu-Al-Ni are presented. We show how the observed behavior emerges within a general framework of competing time scales of avalanche relaxation, driving rate, and thermal fluctuations. We confirm our findings by numerical simulations of a prototype model.

Chaotic, memory, and cooling rate effects in spin glasses: Evaluation of the Edwards-Anderson model

We investigate chaotic, memory, and cooling rate effects in the three-dimensional Edwards-Anderson model by doing thermoremanent (TRM) and ac susceptibility numerical experiments and making a detailed comparison with laboratory experiments on spin glasses. In contrast to the experiments, the Edwards-Anderson model does not show any trace of reinitialization processes in temperature change experiments (TRM or ac). A detailed comparison with ac relaxation experiments in the presence of dc magnetic field or coupling distribution perturbations reveals that the absence of chaotic effects in the Edwards-Anderson model is a consequence of the presence of strong cooling rate effects. We discuss possible solutions to this discrepancy, in particular the smallness of the time scales reached in numerical experiments, but we also question the validity of the Edwards-Anderson model to reproduce the experimental results.

Light scattering from suspensions under external gradients

We analyze the light-scattering spectrum of a suspension in a viscoelastic fluid under density and velocity gradients. When a density gradient is present, the dynamic structure factor exhibits universality in the sense that its expression depends only on the reduced frequency and the reduced density gradient. For a velocity gradient, however, the universality breaks down. In this last case we have found a transition point from one to three characteristic frequencies in the spectrum, which is governed by the value of the external gradient. The presence of the viscoelastic time scales introduces a shift in the ``critical¿¿ point.

Disordered backgammon model

In this paper we consider an exactly solvable model that displays glassy behavior at zero temperature due to entropic barriers. The new ingredient of the model is the existence of different energy scales or modes associated with different relaxational time scales. Low-temperature relaxation takes place by partial equilibration of successive lower-energy modes. An adiabatic scaling solution, defined in terms of a threshold energy scale e*, is proposed. For such a solution, modes with energy ee* are equilibrated at the bath temperature, modes with ee* remain out of equilibrium, and relaxation occurs in the neighborhood of the threshold e~e*. The model is presented as a toy example to investigate the conditions related to the existence of an effective temperature in glassy systems and its possible dependence on the energy sector is probed by the corresponding observable.

Tuning synchronization of integrate-and-fire oscillators through mobility

We analyze the emergence of synchronization in a population of moving integrate-and-fire oscillators. Oscillators, while moving on a plane, interact with their nearest neighbor upon firing time. We discover a nonmonotonic dependence of the synchronization time on the velocity of the agents. Moreover, we find that mechanisms that drive synchronization are different for different dynamical regimes. We report the extreme situation where an interplay between the time scales involved in the dynamical processes completely inhibits the achievement of a coherent state. We also provide estimators for the transitions between the different regimes.

Multiplicative noise effects on relaxations from marginal states

Relaxational processes in bistable potentials close to marginal conditions are studied under the combined effect of additive and multiplicative fluctuations. Characteristic time scales associated with the first-passage-time-distribution are analytically obtained. Multiplicative noise introduces large effects on the characteristic decay times, which is particularly significant when relaxations are mediated by fluctuations, i.e., below marginality and for small noise intensity. The relevance of our approach with respect to realistic chemical bistable systems experimentally operated under external noise influences is mentioned.

Dynamics of transient domain structures in nematic liquid crystals: The splay Fréedericksz transition

We discuss the dynamics of the transient pattern formation process corresponding to the splay Fréedericksz transition. The emergence and subsequent evolution of the spatial periodicity is here described in terms of the temporal dependence of the wave numbers corresponding to the maxima of the structure factor. Situations of perpendicular as well as oblique field-induced stripes relative to the initial orientation of the director are both examined with explicit indications of the time scales needed for their appearance and posterior development.

Modelling hydrodynamics in the Rio Parana, Argentina : an evaluation and inter-comparison of reduced-complexity and physics based models applied to a large sand-bed river

Depth-averaged velocities and unit discharges within a 30 km reach of one of the world's largest rivers, the Rio Parana, Argentina, were simulated using three hydrodynamic models with different process representations: a reduced complexity (RC) model that neglects most of the physics governing fluid flow, a two-dimensional model based on the shallow water equations, and a three-dimensional model based on the Reynolds-averaged Navier-Stokes equations. Row characteristics simulated using all three models were compared with data obtained by acoustic Doppler current profiler surveys at four cross sections within the study reach. This analysis demonstrates that, surprisingly, the performance of the RC model is generally equal to, and in some instances better than, that of the physics based models in terms of the statistical agreement between simulated and measured flow properties. In addition, in contrast to previous applications of RC models, the present study demonstrates that the RC model can successfully predict measured flow velocities. The strong performance of the RC model reflects, in part, the simplicity of the depth-averaged mean flow patterns within the study reach and the dominant role of channel-scale topographic features in controlling the flow dynamics. Moreover, the very low water surface slopes that typify large sand-bed rivers enable flow depths to be estimated reliably in the RC model using a simple fixed-lid planar water surface approximation. This approach overcomes a major problem encountered in the application of RC models in environments characterised by shallow flows and steep bed gradients. The RC model is four orders of magnitude faster than the physics based models when performing steady-state hydrodynamic calculations. However, the iterative nature of the RC model calculations implies a reduction in computational efficiency relative to some other RC models. A further implication of this is that, if used to simulate channel morphodynamics, the present RC model may offer only a marginal advantage in terms of computational efficiency over approaches based on the shallow water equations. These observations illustrate the trade off between model realism and efficiency that is a key consideration in RC modelling. Moreover, this outcome highlights a need to rethink the use of RC morphodynamic models in fluvial geomorphology and to move away from existing grid-based approaches, such as the popular cellular automata (CA) models, that remain essentially reductionist in nature. In the case of the world's largest sand-bed rivers, this might be achieved by implementing the RC model outlined here as one element within a hierarchical modelling framework that would enable computationally efficient simulation of the morphodynamics of large rivers over millennial time scales. (C) 2012 Elsevier B.V. All rights reserved.