Biophysical characterization of tubulin tyrosine ligase-like proteins on microtubule dynamics

Report for the scientific sojourn carried out at the Cell Biology and Biophysics Unit from the National Institutes of Health, from 2010 to 2012.

Distributed all-atom simulations of intrinsically unstructured proteins

The Computational Biophysics Group at the Universitat Pompeu Fabra (GRIB-UPF) hosts two unique computational resources dedicated to the execution of large scale molecular dynamics (MD) simulations: (a) the ACMD molecular-dynamics software, used on standard personal computers with graphical processing units (GPUs); and (b) the GPUGRID. net computing network, supported by users distributed worldwide that volunteer GPUs for biomedical research. We leveraged these resources and developed studies, protocols and open-source software to elucidate energetics and pathways of a number of biomolecular systems, with a special focus on flexible proteins with many degrees of freedom. First, we characterized ion permeation through the bactericidal model protein Gramicidin A conducting one of the largest studies to date with the steered MD biasing methodology. Next, we addressed an open problem in structural biology, the determination of drug-protein association kinetics; we reconstructed the binding free energy, association, and dissaciociation rates of a drug like model system through a spatial decomposition and a Makov-chain analysis. The work was published in the Proceedings of the National Academy of Sciences and become one of the few landmark papers elucidating a ligand-binding pathway. Furthermore, we investigated the unstructured Kinase Inducible Domain (KID), a 28-peptide central to signalling and transcriptional response; the kinetics of this challenging system was modelled with a Markovian approach in collaboration with Frank Noe’s group at the Freie University of Berlin. The impact of the funding includes three peer-reviewed publication on high-impact journals; three more papers under review; four MD analysis components, released as open-source software; MD protocols; didactic material, and code for the hosting group.

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.

Oscillatory magnetic tweezers based on ferromagnetic beads and simple coaxial coils

We report the design and validation of simple magnetic tweezers for oscillating ferromagnetic beads in the piconewton and nanometer scales. The system is based on a single pair of coaxial coils operating in two sequential modes: permanent magnetization of the beads through a large and brief pulse of magnetic field and generation of magnetic gradients to produce uniaxial oscillatory forces. By using this two step method, the magnetic moment of the beads remains constant during measurements. Therefore, the applied force can be computed and varies linearly with the driving signal. No feedback control is required to produce well defined force oscillations over a wide bandwidth. The design of the coils was optimized to obtain high magnetic fields (280 mT) and gradients (2 T/m) with high homogeneity (5% variation) within the sample. The magnetic tweezers were implemented in an inverted optical microscope with a videomicroscopy-based multiparticle tracking system. The apparatus was validated with 4.5 ¿m magnetite beads obtaining forces up to ~2 pN and subnanometer resolution. The applicability of the device includes microrheology of biopolymer and cell cytoplasm, molecular mechanics, and mechanotransduction in living cells.

Nanotecnologies en mecànica molecular i cel·lular

Algunes nanotècniques recents permeten la manipulació de biomolècules i cèl·lules en escala nanomètrica amb la mesura simultània de la força aplicada amb resolució de piconewtons. Aquestes escales de desplaçament i força, i la possibilitat de treballar en medi líquid, fan que siguin eines molt útils per a l'estudi de les propietats mecàniques de molècules i cèl·lules individuals en condicions fisiològiques. Entre les tècniques més utilitzades es troben el microscopi de força atòmica, les trampes de làser i les microesferes magnètiques. En aquest treball es descriuen els principis de funcionament d'aquestes tècniques en aplicacions biològiques i, en particular, en l'estudi de la mecànica molecular i cel·lular.

Model for curvature-driven pearling instability in membranes

A phase-field model for dealing with dynamic instabilities in membranes is presented. We use it to study curvature-driven pearling instability in vesicles induced by the anchorage of amphiphilic polymers on the membrane. Within this model, we obtain the morphological changes reported in recent experiments. The formation of a homogeneous pearled structure is achieved by consequent pearling of an initial cylindrical tube from the tip. For high enough concentration of anchors, we show theoretically that the homogeneous pearled shape is energetically less favorable than an inhomogeneous one, with a large sphere connected to an array of smaller spheres.

Polymer-induced tubulation in lipid vesicles

A mechanism of extraction of tubular membranes from a lipid vesicle is presented. A concentration gradient of anchoring amphiphilic polymers generates tubes from budlike vesicle protrusions. We explain this mechanism in the framework of the Canham-Helfrich model. The energy profile is analytically calculated and a tube with a fixed length, corresponding to an energy minimum, is obtained in a certain regime of parameters. Further, using a phase-field model, we corroborate these results numerically. We obtain the growth of tubes when a polymer source is added, and the budlike shape after removal of the polymer source, in accordance with recent experimental results.

Dynamical properties of discrete breathers in curved chains with first and second neighbor interaction

We present the study of discrete breather dynamics in curved polymerlike chains consisting of masses connected via nonlinear springs. The polymer chains are one dimensional but not rectilinear and their motion takes place on a plane. After constructing breathers following numerically accurate procedures, we launch them in the chains and investigate properties of their propagation dynamics. We find that breather motion is strongly affected by the presence of curved regions of polymers, while the breathers themselves show a very strong resilience and remarkable stability in the presence of geometrical changes. For chains with strong angular rigidity we find that breathers either pass through bent regions or get reflected while retaining their frequency. Their motion is practically lossless and seems to be determined through local energy conservation. For less rigid chains modeled via second neighbor interactions, we find similarly that chain geometry typically does not destroy the localized breather states but, contrary to the angularly rigid chains, it induces some small but constant energy loss. Furthermore, we find that a curved segment acts as an active gate reflecting or refracting the incident breather and transforming its velocity to a value that depends on the discrete breathers frequency. We analyze the physical reasoning behind these seemingly general breather properties.

Interaction of a discrete breather with a lattice junction

We study the scattering of a moving discrete breather (DB) on a junction in a Fermi-Pasta-Ulam chain consisting of two segments with different masses of the particles. We consider four distinct cases: (i) a light-heavy (abrupt) junction in which the DB impinges on the junction from the segment with lighter mass, (ii) a heavy-light junction, (iii) an up mass ramp in which the mass in the heavier segment increases continuously as one moves away from the junction point, and (iv) a down mass ramp. Depending on the mass difference and DB characteristics (frequency and velocity), the DB can either reflect from, or transmit through, or get trapped at the junction or on the ramp. For the heavy-light junction, the DB can even split at the junction into a reflected and a transmitted DB. The latter is found to subsequently split into two or more DBs. For the down mass ramp the DB gets accelerated in several stages, with accompanying radiation (phonons). These results are rationalized by calculating the Peierls-Nabarro barrier for the various cases. We also point out implications of our results in realistic situations such as electron-phonon coupled chains.

Modeling a processive molecular motor: the conventional kinesin

We present a model that allows for the derivation of the experimentally accesible observables: spatial steps, mean velocity, stall force, useful power, efficiency and randomness, etc. as a function of the [adenosine triphosphate] concentration and an external load F. The model presents a minimum of adjustable parameters and the theoretical predictions compare well with the available experimental results.

Noise-enhanced excitability in bistable activator-inhibitor media

We show that external fluctuations induce excitable behavior in a bistable spatially extended system with activator-inhibitor dynamics of the FitzHugh-Nagumo type. This can be understood as a mechanism for sustained signal propagation in bistable media. The phase diagram of the stochastic system is analytically obtained and numerically verified. For small-noise intensities, front propagation becomes unstable, and excitable pulses arise as the only possible spatiotemporal behavior of the system. For large-noise intensities, on the other hand, the system enters an effective regime of oscillatory behavior, where it exhibits spontaneous nucleation of pulses and synchronized firing.

Critical Behavior and Axis Defining Symmetry Breaking in Hydra Embryonic Development

The formation of a hollow cellular sphere is often one of the first steps of multicellular embryonic development. In the case of Hydra, the sphere breaks its initial symmetry to form a foot-head axis. During this process a gene, ks1, is increasingly expressed in localized cell domains whose size distribution becomes scale-free at the axis-locking moment. We show that a physical model based solely on the production and exchange of ks1-promoting factors among neighboring cells robustly reproduces the scaling behavior as well as the experimentally observed spontaneous and temperature-directed symmetry breaking.

Pattern selection in a lattice of pulse-coupled oscillators

We study spatio-temporal pattern formation in a ring of N oscillators with inhibitory unidirectional pulselike interactions. The attractors of the dynamics are limit cycles where each oscillator fires once and only once. Since some of these limit cycles lead to the same pattern, we introduce the concept of pattern degeneracy to take it into account. Moreover, we give a qualitative estimation of the volume of the basin of attraction of each pattern by means of some probabilistic arguments and pattern degeneracy, and show how they are modified as we change the value of the coupling strength. In the limit of small coupling, our estimative formula gives a pefect agreement with numerical simulations.

Mechanisms of synchronization and pattern formation in a lattice of pulse-coupled oscillators

We analyze the physical mechanisms leading either to synchronization or to the formation of spatiotemporal patterns in a lattice model of pulse-coupled oscillators. In order to make the system tractable from a mathematical point of view we study a one-dimensional ring with unidirectional coupling. In such a situation, exact results concerning the stability of the fixed of the dynamic evolution of the lattice can be obtained. Furthermore, we show that this stability is the responsible for the different behaviors.