Steven J. Plimpton

Bio: Steven J. Plimpton is an academic researcher from Sandia National Laboratories. The author has contributed to research in topics: Parallel algorithm & Direct simulation Monte Carlo. The author has an hindex of 44, co-authored 128 publications receiving 62532 citations.

Noise-induced oscillatory shuttling of NF-{\kappa}B in a two compartment IKK-NF-{\kappa}B-I{\kappa}B-A20 signaling model

Abstract: NF-{\kappa}B is a pleiotropic protein whose nucleo-cytoplasmic trafficking is tightly regulated by multiple negative feedback loops embedded in the NF-{\kappa}B signaling network and contributes to diverse gene expression profiles important in immune cell differentiation, cell apoptosis, and innate immunity. The intracellular signaling processes and their control mechanisms, however, are susceptible to both extrinsic and intrinsic noise. In this article, we present numerical evidence for a universal dynamic behavior of NF-{\kappa}B, namely oscillatory nucleo-cytoplasmic shuttling, due to the fundamentally stochastic nature of the NF-{\kappa}B signaling network. We simulated the effect of extrinsic noise with a deterministic ODE model, using a statistical ensemble approach, generating many copies of the signaling network with different kinetic rates sampled from a biologically feasible parameter space. We modeled the effect of intrinsic noise by simulating the same networks stochastically using the Gillespie algorithm. The results demonstrate that extrinsic noise diversifies the shuttling patterns of NF-{\kappa}B response, whereas intrinsic noise induces oscillatory behavior in many of the otherwise non-oscillatory patterns. We identify two key model parameters which significantly affect the NF-{\kappa}B dynamic response and deduce a two-dimensional phase-diagram of the NF-{\kappa}B response as a function of these parameters. We conclude that if single-cell experiments are performed, a rich variety of NF-{\kappa}B response will be observed, even if population-level experiments, which average response over large numbers of cells, do not evidence oscillatory behavior.

Parallel Tempering Monte Carlo in LAMMPS

Abstract: We present here the details of the implementation of the parallel tempering Monte Carlo technique into a LAMMPS, a heavily used massively parallel molecular dynamics code at Sandia. This technique allows for many replicas of a system to be run at different simulation temperatures. At various points in the simulation, configurations can be swapped between different temperature environments and then continued. This allows for large regions of energy space to be sampled very quickly, and allows for minimum energy configurations to emerge in very complex systems, such as large biomolecular systems. By including this algorithm into an existing code, we immediately gain all of the previous work that had been put into LAMMPS, and allow this technique to quickly be available to the entire Sandia and international LAMMPS community. Finally, we present an example of this code applied to folding a small protein.

Nanoparticle flow, ordering and self-assembly.

Abstract: Nanoparticles are now more than ever being used to tailor materials function and performance in differentiating technologies because of their profound effect on thermo-physical, mechanical and optical properties. The most feasible way to disperse particles in a bulk material or control their packing at a substrate is through fluidization in a carrier, followed by solidification through solvent evaporation/drying/curing/sintering. Unfortunately processing particles as concentrated, fluidized suspensions into useful products remains an art largely because the effect of particle shape and volume fraction on fluidic properties and suspension stability remains unexplored in a regime where particle-particle interaction mechanics is prevalent. To achieve a stronger scientific understanding of the factors that control nanoparticle dispersion and rheology we have developed a multiscale modeling approach to bridge scales between atomistic and molecular-level forces active in dense nanoparticle suspensions. At the largest length scale, two 'coarse-grained' numerical techniques have been developed and implemented to provide for high-fidelity numerical simulations of the rheological response and dispersion characteristics typical in a processing flow. The first is a coupled Navier-Stokes/discrete element method in which the background solvent is treated by finite element methods. The second is a particle based method known as stochastic rotational dynamics. These two methods providemore » a new capability representing a 'bridge' between the molecular scale and the engineering scale, allowing the study of fluid-nanoparticle systems over a wide range of length and timescales as well as particle concentrations. To validate these new methodologies, multi-million atoms simulations explicitly including the solvent have been carried out. These simulations have been vital in establishing the necessary 'subgrid' models for accurate prediction at a larger scale and refining the two coarse-grained methodologies.« less

Fast parallel algorithms for short-range molecular dynamics

Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

Abstract: 抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。