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.

Fast parallel algorithms for short-range molecular dynamics

http://nano-bio.ehu.es/files/articles/Gonze_CPC_2009_590.pdf

ABINIT: First-principles approach to material and nanosystem properties

Chemokines are chemotactic cytokines that control the migratory patterns and positioning of all immune cells. Although chemokines were initially appreciated as important mediators of acute inflammation, we now know that this complex system of approximately 50 endogenous chemokine ligands and 20 G protein–coupled seven-transmembrane signaling receptors is also critical for the generation of primary and secondary adaptive cellular and humoral immune responses. Recent studies demonstrate important roles for the chemokine system in the priming of naive T cells, in cell fate decisions such as effector and memory cell differentiation, and in regulatory T cell function. In this review, we focus on recent advances in understanding how the chemokine system orchestrates immune cell migration and positioning at the organismic level in homeostasis, in acute inflammation, and during the generation and regulation of adoptive primary and secondary immune responses in the lymphoid system and peripheral nonlymphoid tissue.

Chemokines and Chemokine Receptors: Positioning Cells for Host Defense and Immunity

Significance When the human genome folds up inside the cell nucleus, it is spatially partitioned into numerous loops and contact domains. How these structures form is unknown. Here, we show that data from high-resolution spatial proximity maps are consistent with a model in which a complex, including the proteins CCCTC-binding factor (CTCF) and cohesin, mediates the formation of loops by a process of extrusion. Contact domains form as a byproduct of this process. The model accurately predicts how the genome will fold, using only information about the locations at which CTCF is bound. We demonstrate the ability to reengineer loops and domains in a predictable manner by creating highly targeted mutations, some as small as a single base pair, at CTCF sites. We recently used in situ Hi-C to create kilobase-resolution 3D maps of mammalian genomes. Here, we combine these maps with new Hi-C, microscopy, and genome-editing experiments to study the physical structure of chromatin fibers, domains, and loops. We find that the observed contact domains are inconsistent with the equilibrium state for an ordinary condensed polymer. Combining Hi-C data and novel mathematical theorems, we show that contact domains are also not consistent with a fractal globule. Instead, we use physical simulations to study two models of genome folding. In one, intermonomer attraction during polymer condensation leads to formation of an anisotropic “tension globule.” In the other, CCCTC-binding factor (CTCF) and cohesin act together to extrude unknotted loops during interphase. Both models are consistent with the observed contact domains and with the observation that contact domains tend to form inside loops. However, the extrusion model explains a far wider array of observations, such as why loops tend not to overlap and why the CTCF-binding motifs at pairs of loop anchors lie in the convergent orientation. Finally, we perform 13 genome-editing experiments examining the effect of altering CTCF-binding sites on chromatin folding. The convergent rule correctly predicts the affected loops in every case. Moreover, the extrusion model accurately predicts in silico the 3D maps resulting from each experiment using only the location of CTCF-binding sites in the WT. Thus, we show that it is possible to disrupt, restore, and move loops and domains using targeted mutations as small as a single base pair.

Chromatin extrusion explains key features of loop and domain formation in wild-type and engineered genomes

Tissue-Resident Memory T Cells

Chemokines have a central role in regulating processes essential to the immune function of T cells, such as their migration within lymphoid tissues and targeting of pathogens in sites of inflammation. Here we track T cells using multi-photon microscopy to demonstrate that the chemokine CXCL10 enhances the ability of CD8+ T cells to control the pathogen Toxoplasma gondii in the brains of chronically infected mice. This chemokine boosts T-cell function in two different ways: it maintains the effector T-cell population in the brain and speeds up the average migration speed without changing the nature of the walk statistics. Notably, these statistics are not Brownian; rather, CD8+ T-cell motility in the brain is well described by a generalized Levy walk. According to our model, this unexpected feature enables T cells to find rare targets with more than an order of magnitude more efficiency than Brownian random walkers. Thus, CD8+ T-cell behaviour is similar to Levy strategies reported in organisms ranging from mussels to marine predators and monkeys, and CXCL10 aids T cells in shortening the average time taken to find rare targets.

/pdf/generalized-levy-walks-and-the-role-of-chemokines-in-kx5jlbglcb.pdf

Generalized Levy walks and the role of chemokines in migration of effector CD8+ T cells

The cell nucleus must continually resist and respond to intercellular and intracellular mechanical forces to transduce mechanical signals and maintain proper genome organization and expression. Altered nuclear mechanics is associated with many human diseases, including heart disease, progeria, and cancer. Chromatin and nuclear envelope A-type lamin proteins are known to be key nuclear mechanical components perturbed in these diseases, but their distinct mechanical contributions are not known. Here we directly establish the separate roles of chromatin and lamin A/C and show that they determine two distinct mechanical regimes via micromanipulation of single isolated nuclei. Chromatin governs response to small extensions (<3 μm), and euchromatin/heterochromatin levels modulate the stiffness. In contrast, lamin A/C levels control nuclear strain stiffening at large extensions. These results can be understood through simulations of a polymeric shell and cross-linked polymer interior. Our results provide a framework for understanding the differential effects of chromatin and lamin A/C in cell nuclear mechanics and their alterations in disease.

Chromatin and lamin A determine two different mechanical response regimes of the cell nucleus

Chromatin decompaction via increasing euchromatin or decreasing heterochromatin results in a softer nucleus and abnormal nuclear blebbing, independent of lamin perturbations. Conversely, increasing...

Edward J. Banigan

Papers

Generalized Levy walks and the role of chemokines in migration of effector CD8+ T cells

Chromatin and lamin A determine two different mechanical response regimes of the cell nucleus

Chromatin histone modifications and rigidity affect nuclear morphology independent of lamins

Chromatin’s physical properties shape the nucleus and its functions

Loop extrusion: theory meets single-molecule experiments.