Showing papers in "arXiv: Biological Physics in 2009"

Folding@Home and Genome@Home: Using distributed computing to tackle previously intractable problems in computational biology

Abstract: For decades, researchers have been applying computer simulation to address problems in biology. However, many of these "grand challenges" in computational biology, such as simulating how proteins fold, remained unsolved due to their great complexity. Indeed, even to simulate the fastest folding protein would require decades on the fastest modern CPUs. Here, we review novel methods to fundamentally speed such previously intractable problems using a new computational paradigm: distributed computing. By efficiently harnessing tens of thousands of computers throughout the world, we have been able to break previous computational barriers. However, distributed computing brings new challenges, such as how to efficiently divide a complex calculation of many PCs that are connected by relatively slow networking. Moreover, even if the challenge of accurately reproducing reality can be conquered, a new challenge emerges: how can we take the results of these simulations (typically tens to hundreds of gigabytes of raw data) and gain some insight into the questions at hand. This challenge of the analysis of the sea of data resulting from large-scale simulation will likely remain for decades to come.

The Geant4-DNA project

Abstract: The Geant4-DNA project proposes to develop an open-source simulation software based and fully included in the general-purpose Geant4 Monte Carlo simulation toolkit. The main objective of this software is to simulate biological damages induced by ionising radiation at the cellular and sub-cellular scale. This project was originally initiated by the European Space Agency for the prediction of deleterious effects of radiation that may affect astronauts during future long duration space exploration missions. In this paper, the Geant4-DNA collaboration presents an overview of the whole ongoing project, including its most recent developments already available in the last Geant4 public release (9.3 BETA), as well as an illustration example simulating the direct irradiation of a chromatin fibre. Expected extensions involving several research domains, such as particle physics, chemistry and cellular and molecular biology, within a fully interdiciplinary activity of the Geant4 collaboration are also discussed.

Metastable Chimera States in Community-Structured Oscillator Networks

Abstract: A system of symmetrically coupled identical oscillators with phase lag is presented, which is capable of generating a large repertoire of transient (metastable) "chimera" states in which synchronisation and desynchronisation co-exist. The oscillators are organised into communities, such that each oscillator is connected to all its peers in the same community and to a subset of the oscillators in other communities. Measures are introduced for quantifying metastability, the prevalence of chimera states, and the variety of such states a system generates. By simulation, it is shown that each of these measures is maximised when the phase lag of the model is close, but not equal, to pi/2. The relevance of the model to a number of fields is briefly discussed, with particular emphasis on brain dynamics.

Theoretical Study of Physisorption of Nucleobases on Boron Nitride Nanotubes: A New Class of Hybrid Nano-Bio Materials

Abstract: We investigate the adsorption of the nucleic acid bases, adenine (A), guanine (G), cytosine (C), thymine (T) and uracil (U) on the outer wall of a high curvature semiconducting single-walled boron nitride nanotube (BNNT) by first principles density functional theory calculations. The calculated binding energy shows the order: G>A\approxC\approxT\approxU implying that the interaction strength of the (high-curvature) BNNT with the nucleobases, G being an exception, is nearly the same. A higher binding energy for the G-BNNT conjugate appears to result from a stronger hybridization of the molecular orbitals of G and BNNT, since the charge transfer involved in the physisorption process is insignificant. A smaller energy gap predicted for the G-BNNT conjugate relative to that of the pristine BNNT may be useful in application of this class of biofunctional materials to the design of the next generation sensing devices.

Shear-Induced Unfolding Activates von Willebrand Factor A2 Domain for Proteolysis

Abstract: To avoid pathological platelet aggregation by von Willebrand factor (VWF), VWF multimers are regulated in size and reactivity for adhesion by ADAMTS13-mediated proteolysis in a shear flow dependent manner. We examined if tensile stress in VWF under shear flow activates the VWF A2 domain for cleavage by ADAMTS13 using molecular dynamics simulations. We indeed observed stepwise unfolding of A2 and exposure of its deeply buried ADAMTS13 cleavage site. Interestingly, disulfide bonds in the adjacent and highly homologous VWF A1 and A3 domains obstruct their mechanical unfolding. We generated a full length mutant VWF featuring a homologous disulfide bond in A2 (N1493C and C1670S), in an attempt to lock A2 against unfolding. We find this mutant to feature ADAMTS13-resistant behavior in vitro. Our results yield molecular-detail evidence for the force-sensoring function of VWF A2, by revealing how tension in VWF due to shear flow selectively exposes the A2 proteolysis site to ADAMTS13 for cleavage while keeping the folded remainder of A2 intact and functional. We find the unconventional knotted Rossman fold of A2 to be the key to this mechanical response, tailored for regulating VWF size and activity. Based on our model we can explain the pathomechanism of some natural mutations in the VWF A2 domain that significantly increase the cleavage by ADAMTS13 without shearing or chemical denaturation, and provide with the cleavage-activated A2 conformation a structural basis for the design of inhibitors for VWF type 2 diseases.

Generalized Voronoi Tessellation as a Model of Two-dimensional Cell Tissue Dynamics

Abstract: Voronoi tessellations have been used to model the geometric arrangement of cells in morphogenetic or cancerous tissues, however so far only with flat hypersurfaces as cell-cell contact borders. In order to reproduce the experimentally observed piecewise spherical boundary shapes, we develop a consistent theoretical framework of multiplicatively weighted distance functions, defining generalized finite Voronoi neighborhoods around cell bodies of varying radius, which serve as heterogeneous generators of the resulting model tissue. The interactions between cells are represented by adhesive and repelling force densities on the cell contact borders. In addition, protrusive locomotion forces are implemented along the cell boundaries at the tissue margin, and stochastic perturbations allow for non-deterministic motility effects. Simulations of the emerging system of stochastic differential equations for position and velocity of cell centers show the feasibility of this Voronoi method generating realistic cell shapes. In the limiting case of a single cell pair in brief contact, the dynamical nonlinear Ornstein-Uhlenbeck process is analytically investigated. In general, topologically distinct tissue conformations are observed, exhibiting stability on different time scales, and tissue coherence is quantified by suitable characteristics. Finally, an argument is derived pointing to a tradeoff in natural tissues between cell size heterogeneity and the extension of cellular lamellae.

Infrared Surface-Plasmon-Resonance -- a novel biophysical tool for studying living cell

Abstract: We discuss the Surface-Plasmon-Resonance (SPR) technique based on Fourier -Transform - InfraRed (FTIR) spectrometry. We explore the potential of the infrared surface plasmon resonance technique for biological studies in aqueous solutions and compare it to the conventional surface plasmon technique operating in the visible range. We demonstrate that the sensitivity of the SPR technique in the infrared range is not lower and in fact is even higher. We show several examples of applying FTIR-SPR for biological studies: (i) monitoring D-glucose concentration in solution, and (ii) measuring D-glucose uptake by erythrocytes in suspension. We emphasize the advantages of infrared SPR for studying living cell cultures and show how this technique can be used for characterization of (i) cholesterol penetration into plasma membrane, and (ii) transferrin-induced clathrin-mediated endocytosis.

Mapping the conformations of biological assemblies

Abstract: Mapping conformational heterogeneity of macromolecules presents a formidable challenge to X-ray crystallography and cryo-electron microscopy, which often presume its absence. This has severely limited our knowledge of the conformations assumed by biological systems and their role in biological function, even though they are known to be important. We propose a new approach to determining to high resolution the three-dimensional conformations of biological entities such as molecules, macromolecular assemblies, and ultimately cells, with existing and emerging experimental techniques. This approach may also enable one to circumvent current limits due to radiation damage and solution purification.

The effects of fluid viscosity on the kinematics and material properties of C. elegans swimming at low Reynolds number

Abstract: The effects of fluid viscosity on the kinematics of a small swimmer at low Reynolds number are investigated in both experiments and in a simple model. The swimmer is the nematode Caenorhabditis elegans, which is an undulating roundworm approximately 1 mm long. Experiments show that the nematode maintains a highly periodic swimming behavior as the fluid viscosity is varied from 1.0 mPa-s to 12 mPa-s. Surprisingly, the nematode's swimming speed (~0.35 mm/s) is nearly insensitive to the range of fluid viscosities investigated here. However, the nematode's beating frequency decreases to an asymptotic value (~1.7 Hz) with increasing fluid viscosity. A simple model is used to estimate the nematode's Young's modulus and tissue viscosity. Both material properties increase with increasing fluid viscosity. It is proposed that the increase in Young's modulus may be associated with muscle contraction in response to larger mechanical loading while the increase in effective tissue viscosity may be associated with the energy necessary to overcome increased fluid drag forces.

On axoplasmic pressure waves and their possible role in nerve impulse propagation

Abstract: It is suggested that the propagation of the action potential is accompanied by an axoplasmic pressure pulse propagating in the axoplasm along the axon length. The pressure pulse stretch-modulates voltage-gated Na (Nav) channels embedded in the axon membrane, causing their accelerated activation and inactivation and increasing peak channel conductance. As a result, the action potential propagates due to mechano-electrical activation of Nav channels by straggling ionic currents and the axoplasmic pressure pulse. The velocity of such propagation is higher than in the classical purely electrical Nav activation mechanism, and it may be close to the velocity of propagation of pressure pulses in the axoplasm. Extracellular Ca ions influxing during the voltage spike, or Ca ions released from intracellular stores, may trigger a mechanism that generates and augments the pressure pulse, thus opposing its viscous decay. The model can potentially explain a number of phenomena that are not contained within the purely electrical Hodgkin-Huxley-type framework: the Meyer-Overton rule for the effectiveness of anesthetics, as well as various mechanical, optical and thermodynamic phenomena accompanying the action potential. It is shown that the velocity of propagation of axoplasmic pressure pulses is close to the measured velocity of the nerve impulse, both in absolute magnitude and in dependence on axon diameter, degree of myelination and temperature.

Heterogeneous path ensembles for conformational transitions in semi-atomistic models of adenylate kinase

Abstract: We performed "weighted ensemble" path-sampling simulations of adenylate kinase, using several semi-atomistic protein models. Our study investigated both the biophysics of conformational transitions as well as the possibility of increasing model accuracy without sacrificing good sampling. Biophysically, the path ensembles show significant heterogeneity and the explicit possibility of two principle pathways in the Open-Closed transition. We recently showed, under certain conditions, a "symmetry of hetereogeneity" is expected between the forward and the reverse transitions: the fraction of transitions taking a specific pathway/channel will be the same in both the directions. Our path ensembles are analyzed in the light of the symmetry relation and its conditions. In the realm of modeling, we employed an all-atom backbone with various levels of residue interactions. Because reasonable path sampling required only a few weeks of single-processor computing time with these models, the addition of further chemical detail should be feasible.

Measurement of Cytoplasmic Streaming in Chara Corallina by Magnetic Resonance Velocimetry

Abstract: In aquatic plants such as the Characean algae, the force generation that drives cyclosis is localized within the cytoplasm, yet produces fluid flows throughout the vacuole. For this to occur the tonoplast must transmit hydrodynamic shear efficiently. Here, using magnetic resonance velocimetry, we present the first whole-cell measurements of the cross-sectional longitudinal velocity field in Chara corallina and show that it is in quantitative agreement with a recent theoretical analysis of rotational cytoplasmic streaming driven by bidirectional helical forcing in the cytoplasm, with direct shear transmission by the tonoplast.

Influence of Local and Residual Structures on the Scaling Behavior and Dimensions of Unfolded Proteins

Abstract: Although recent spectroscopic studies of chemically denatured proteins hint at significant nonrandom residual structure, the results of extensive small angle X-ray scattering studies suggest random coil behavior, calling for a coherent understanding of these seemingly contradicting observations. Here, we report the results of a Monte Carlo study of the effects of two types of local structures, a helix and Polyproline II (PPII) helix, on the dimensions of random coil polyalanine chains viewed as a model of highly denatured proteins. With an alpha helix content of 20%, corresponding to the Ramachandran probability of being in the helical basin, experimentally observed radii of gyration are recovered. Experimental radii are similarly recovered at an a helix content of 87%, providing an explanation for the previously puzzling experimental finding that the dimensions of the highly helical methanol-induced unfolded state are experimentally indistinguishable from those of the helix-poor urea-unfolded state. In contrast, the radius of gyration increases monotonically with increasing PPII content, and is always more expanded than the dimensions observed experimentally. These results suggest that PPII is unlikely the sole, dominant preferred conformation for unfolded proteins.

Molecular Dynamics Simulations of Cytochrome c un-folding in AOT Reverse Micelles: the first steps

Abstract: This paper explores the reduced form of horse cytochrome c confined in reverse micelles (RM) of so-dium bis-(2-ethylhexyl) sulfosuccinate (AOT) in isooctane by molecular dynamics simulation. RMs of two sizes were constructed at a water content of Wo = [H2O]/[AOT] = 5.5 and 9.1. Our results show that the protein secondary structure and the heme conformation both depend on micellar hydration. At low hydration, the protein structure and the heme moiety remain stable, whereas at high water content the protein becomes unstable and starts to unfold. At Wo = 9.1, according to the X-ray structure, conforma-tional changes are mainly localized on protein loops and around the heme moiety, where we observe a partial opening of the heme crevice. These findings suggest that within our time window (10 ns), the structural changes observed at the heme level are the first steps of the protein denaturation process, pre-viously described experimentally in micellar solutions. In addition, a specific binding of AOT molecules to a few lysine residues of the protein was found only in the small-sized RM.

Scaling of human body mass with height: the Body Mass Index revisited

Abstract: We adapt a biomechanical argument of Rashevsky, which places limits on the stress experienced by a torso supported by the legs, to deduce that body mass $m$ of growing children should scale as the $p$th power of height $h$ with $7/3

Modelling light-driven proton pumps in artificial photosynthetic reaction centers

Abstract: We study a model of a light-induced proton pump in artificial reaction centers. The model contains a molecular triad with four electron states (i.e., one donor state, two photosensitive group states, and one acceptor state) as well as a molecular shuttle having one electron and one proton-binding sites. The shuttle diffuses between the sides of the membrane and translocates protons energetically uphill: from the negative side to the positive side of the membrane, harnessing for this purpose the energy of the electron-charge-separation produced by light. Using methods of quantum transport theory we calculate the range of light intensity and transmembrane potentials that maximize both the light-induced proton current and the energy transduction efficiency. We also study the effect of temperature on proton pumping. The light-induced proton pump in our model gives a quantum yield of proton translocation of about 55 %. Thus, our results explain previous experiments on these artificial photosynthetic reaction centers.

How proteins squeeze through polymer networks: a Cartesian lattice study

Abstract: In this paper a lattice model for the diffusional transport of particles in the interphase cell nucleus is proposed. The dynamic behaviour of single chains on the lattice is investigated and Rouse scaling is verified. Dynamical dense networks are created by a combined version of the bond fluctuation method and a Metropolis Monte Carlo algorithm. Semidilute behaviour of the dense chain networks is shown. By comparing diffusion of particles in a static and a dynamical chain network, we demonstrate that chain diffusion does not alter the diffusion process of small particles. However, we prove that a dynamical network facilitates the transport of large particles. By weighting the mean square displacement trajectories of particles in the static chain network data from the dynamical network can be reconstructed. Additionally, it is shown that subdiffusive behaviour of particles on short time scales results from trapping processes in the crowded environment of the chain network. In the presented model a protein with 30 nm diameter has an effective diffusion coefficient of 1.24E-11 m^2/s in a chromatin fiber network.

Dynamical model of DNA-protein interaction: effect of protein charge distribution and mechanical properties

Abstract: The mechanical model based on beads and springs, which we recently proposed to study non-specific DNA-protein interactions [J. Chem. Phys. 130, 015103 (2009)], was improved by describing proteins as sets of interconnected beads instead of single beads. In this paper, we first compare the results obtained with the updated model with those of the original one and then use it to investigate several aspects of the dynamics of DNA sampling, which could not be accounted for by the original model. These aspects include the effect on the speed of DNA sampling of the regularity and/or randomness of the protein charge distribution, the charge and location of the search site, and the shape and deformability of the protein. We also discuss the efficiency of facilitated diffusion, that is, the extent to which the combination of 1D sliding along the DNA and 3D diffusion in the cell can lead to faster sampling than pure 3D diffusion of the protein.

The intrinsic load-resisting capacity of kinesin

Abstract: Conventional kinesin is a homodimeric motor protein that is capable of walking unidirectionally along a cytoskeletal filament. While previous experiments indicated unyielding unidirectionality against an opposing load up to the so-called stall force, recent experiments also observed limited processive backwalking under superstall loads. This theoretical study seeks to elucidate the molecular mechanical basis for kinesin's steps over the full range of external loads that can possibly be applied to the dimer. We found that kinesin's load-resisting capacity is largely determined by a synergic ratchet-and-pawl mechanism inherent in the dimer. Load susceptibility of this inner molecular mechanical mechanism underlies kinesin's response to various levels of external loads. Computational implementation of the mechanism enabled us to rationalize major trends observed experimentally in kinesin's stalemate and consecutive back steps. The study also predicts several distinct features of kinesin's load-affected motility, which are seemingly counterintuitive but readily verifiable by future experiment.

Steady-state simulations using weighted ensemble path sampling

Abstract: We extend the weighted ensemble (WE) path sampling method to perform rigorous statistical sampling for systems at steady state. The straightforward steady-state implementation of WE is directly practical for simple landscapes, but not when significant metastable intermediates states are present. We therefore develop an enhanced WE scheme, building on existing ideas, which accelerates attainment of steady state in complex systems. We apply both WE approaches to several model systems confirming their correctness and efficiency by comparison with brute-force results. The enhanced version is significantly faster than the brute force and straightforward WE for systems with WE bins that accurately reflect the reaction coordinate(s). The new WE methods can also be applied to equilibrium sampling, since equilibrium is a steady state.

Collective fluctuations in networks of noisy components

Abstract: Collective dynamics result from interactions among noisy dynamical components. Examples include heartbeats, circadian rhythms, and various pattern formations. Because of noise in each component, collective dynamics inevitably involve fluctuations, which may crucially affect functioning of the system. However, the relation between the fluctuations in isolated individual components and those in collective dynamics is unclear. Here we study a linear dynamical system of networked components subjected to independent Gaussian noise and analytically show that the connectivity of networks determines the intensity of fluctuations in the collective dynamics. Remarkably, in general directed networks including scale-free networks, the fluctuations decrease more slowly with the system size than the standard law stated by the central limit theorem. They even remain finite for a large system size when global directionality of the network exists. Moreover, such nontrivial behavior appears even in undirected networks when nonlinear dynamical systems are considered. We demonstrate it with a coupled oscillator system.

Enzyme-Based Logic Systems for Information Processing

Abstract: We review enzymatic systems which involve biocatalytic reactions utilized for information processing (biocomputing). Extensive ongoing research in biocomputing, mimicking Boolean logic gates has been motivated by potential applications in biotechnology and medicine. Furthermore, novel sensor concepts have been contemplated with multiple inputs processed biochemically before the final output is coupled to transducing "smart-material" electrodes and other systems. These applications have warranted recent emphasis on networking of biocomputing gates. First few-gate networks have been experimentally realized, including coupling, for instance, to signal-responsive electrodes for signal readout. In order to achieve scalable, stable network design and functioning, considerations of noise propagation and control have been initiated as a new research direction. Optimization of single enzyme-based gates for avoiding analog noise amplification has been explored, as were certain network-optimization concepts. We review and exemplify these developments, as well as offer an outlook for possible future research foci. The latter include design and uses of non-Boolean network elements, e.g., filters, as well as other developments motivated by potential novel sensor and biotechnology applications.

Detrended Fluctuation Analysis of Systolic Blood Pressure Control Loop

Abstract: We use detrended fluctuation analysis (DFA) to study the dynamics of blood pressure oscillations and its feedback control in rats by analyzing systolic pressure time series before and after a surgical procedure that interrupts its control loop. We found, for each situation, a crossover between two scaling regions characterized by exponents that reflect the nature of the feedback control and its range of operation. In addition, we found evidences of adaptation in the dynamics of blood pressure regulation a few days after surgical disruption of its main feedback circuit. Based on the paradigm of antagonistic, bipartite (vagal and sympathetic) action of the central nerve system, we propose a simple model for pressure homeostasis as the balance between two nonlinear opposing forces, successfully reproducing the crossover observed in the DFA of actual pressure signals.

Binary Fingerprints at Fluctuation-Enhanced Sensing

Abstract: We developed a simple way to generate binary patterns based on spectral slopes in different frequency ranges at fluctuation-enhanced sensing. Such patterns can be considered as binary "fingerprints" of odors. The method has experimentally been demonstrated with a commercial semiconducting metal oxide (Taguchi) sensor exposed to bacterial odors (Escherichia coli and Anthrax-surrogate Bacillus subtilis) and processing their stochastic signals. With a single Taguchi sensor, the situations of empty chamber, tryptic soy agar (TSA) medium, or TSA with bacteria could be distinguished with 100% reproducibility. The bacterium numbers were in the range of 25 thousands to 1 million. To illustrate the relevance for ultra-low power consumption, we show that this new type of signal processing and pattern recognition task can be implemented by a simple analog circuitry and a few logic gates with total power consumption in the microWatts range.

Transport by molecular motors in the presence of static defects

Abstract: The transport by molecular motors along cytoskeletal filaments is studied theoretically in the presence of static defects. The movements of single motors are described as biased random walks along the filament as well as binding to and unbinding from the filament. Three basic types of defects are distinguished, which differ from normal filament sites only in one of the motors' transition probabilities. Both stepping defects with a reduced probability for forward steps and unbinding defects with an increased probability for motor unbinding strongly reduce the velocities and the run lengths of the motors with increasing defect density. For transport by single motors, binding defects with a reduced probability for motor binding have a relatively small effect on the transport properties. For cargo transport by motors teams, binding defects also change the effective unbinding rate of the cargo particles and are expected to have a stronger effect.

A simple self-organized swimmer driven by molecular motors

Abstract: We investigate a self-organized swimmer at low Reynolds numbers. The microscopic swimmer is composed of three spheres that are connected by two identical active linker arms. Each linker arm contains molecular motors and elastic elements and can oscillate spontaneously. We find that such a system immersed in a viscous fluid can self-organize into a state of directed swimming. The swimmer provides a simple system to study important aspects of the swimming of micro-organisms.

Segregation of Polymers in Confined Spaces

Abstract: We investigate the motion of two overlapping polymers with self-avoidance confined in a narrow 2d box. A statistical model is constructed using blob free-energy arguments. We find spontaneous segregation under the condition: $L > R_{//}$, and mixing under $L < R_{//}$, where L is the length of the box, and $R_{//}$ the polymer extension in an infinite slit. Segregation time scales are determined by solving a mean first-passage time problem, and by performing Monte Carlo simulations. Predictions of the two methods show good agreement. Our results may elucidate a driving force for chromosomes segregation in bacteria.

Self-Propelled Rods near Surfaces

Abstract: We study the behavior of self-propelled nano- and micro-rods in three dimensions, confined between two parallel walls, by simulations and scaling arguments. Our simulations include thermal fluctuations and hydrodynamic interactions, which are both relevant for the dynamical behavior at nano- to micrometer length scales. In order to investigate the importance hydrodynamic interactions, we also perform Brownian-dynamics-like simulations. In both cases, we find that self-propelled rods display a strong surface excess in confined geometries. An analogy with semi-flexible polymers is employed to derive scaling laws for the dependence on the wall distance, the rod length, and the propulsive force. The simulation data confirm the scaling predictions.

Active transport and cluster formation on 2D networks

Abstract: We introduce a model for active transport on inhomogeneous networks embedded in a diffusive environment and investigate the formation of particle clusters. In the presence of a hard-core interaction, cluster sizes exhibit an algebraically decaying distribution in a large parameter regime, indicating the existence of clusters on all scales. The results are compared with a diffusion limited aggregation model and active transport on a regular network. For both models we observe aggregation of particles to clusters which are characterized by a finite size-scale if the relevant time-scales and particle densities are considered.

Absolute free energies estimated by combining pre-calculated molecular fragment libraries

Abstract: The absolute free energy -- or partition function, equivalently -- of a molecule can be estimated computationally using a suitable reference system. Here, we demonstrate a practical method for staging such calculations by growing a molecule based on a series of fragments. Significant computer time is saved by pre-calculating fragment configurations and interactions for re-use in a variety of molecules. We employ such fragment libraries and interaction tables for amino acids and capping groups to estimate free energies for small peptides. Equilibrium ensembles for the molecules are generated at no additional computational cost, and are used to check our results by comparison to standard dynamics simulation.