Data Assimilation using a GPU Accelerated Path Integral Monte Carlo Approach

The answers to data assimilation questions can be expressed as path integrals over all possible state and parameter histories. We show how these path integrals can be evaluated numerically using a Markov Chain Monte Carlo method designed to run in parallel on a Graphics Processing Unit (GPU). We demonstrate the application of the method to an example with a transmembrane voltage time series of a simulated neuron as an input, and using a Hodgkin-Huxley neuron model. By taking advantage of GPU computing, we gain a parallel speedup factor of up to about 300, compared to an equivalent serial computation on a CPU, with performance increasing as the length of the observation time used for data assimilation increases.Comment: 5 figures, submitted to Journal of Computational Physic

Structure Function of Polymer Nematic Liquid Crystals: A Monte Carlo Simulation

We present a Monte Carlo simulation of a polymer nematic for varying volume fractions, concentrating on the structure function of the sample. We achieve nematic ordering with stiff polymers made of spherical monomers that would otherwise not form a nematic state. Our results are in good qualitative agreement with theoretical and experimental predictions, most notably the bowtie pattern in the static structure function.Comment: 10 pages, plain TeX, macros included, 3 figures available from archive. Published versio

Hydrophobic and ionic-interactions in bulk and confined water with implications for collapse and folding of proteins

Water and water-mediated interactions determine thermodynamic and kinetics of protein folding, protein aggregation and self-assembly in confined spaces. To obtain insights into the role of water in the context of folding problems, we describe computer simulations of a few related model systems. The dynamics of collapse of eicosane shows that upon expulsion of water the linear hydrocarbon chain adopts an ordered helical hairpin structure with 1.5 turns. The structure of dimer of eicosane molecules has two well ordered helical hairpins that are stacked perpendicular to each other. As a prelude to studying folding in confined spaces we used simulations to understand changes in hydrophobic and ionic interactions in nano droplets. Solvation of hydrophobic and charged species change drastically in nano water droplets. Hydrophobic species are localized at the boundary. The tendency of ions to be at the boundary where water density is low increases as the charge density decreases. Interaction between hydrophobic, polar, and charged residue are also profoundly altered in confined spaces. Using the results of computer simulations and accounting for loss of chain entropy upon confinement we argue and then demonstrate, using simulations in explicit water, that ordered states of generic amphiphilic peptide sequences should be stabilized in cylindrical nanopores

Ligand-Receptor Interactions

The formation and dissociation of specific noncovalent interactions between a variety of macromolecules play a crucial role in the function of biological systems. During the last few years, three main lines of research led to a dramatic improvement of our understanding of these important phenomena. First, combination of genetic engineering and X ray cristallography made available a simultaneous knowledg of the precise structure and affinity of series or related ligand-receptor systems differing by a few well-defined atoms. Second, improvement of computer power and simulation techniques allowed extended exploration of the interaction of realistic macromolecules. Third, simultaneous development of a variety of techniques based on atomic force microscopy, hydrodynamic flow, biomembrane probes, optical tweezers, magnetic fields or flexible transducers yielded direct experimental information of the behavior of single ligand receptor bonds. At the same time, investigation of well defined cellular models raised the interest of biologists to the kinetic and mechanical properties of cell membrane receptors. The aim of this review is to give a description of these advances that benefitted from a largely multidisciplinar approach