Parallelizing a DNA simulation code for the Cray MTA-2

The Cray MTA-2 (Multithreaded Architecture) is an unusual parallel supercomputer that promises ease of use and high performance. We describe our experience on the MTA-2 with a molecular dynamics code, SIMU-MD, that we are using to simulate the translocation of DNA through a nanopore in a silicon based ultrafast sequencer. Our sequencer is constructed using standard VLSI technology and consists of a nanopore surrounded by field effect transistors (FETs). We propose to use the FETs to sense variations in charge as a DNA molecule translocates through the pore and thus differentiate between the four building block nucleotides of DNA. We were able to port SIMU-MD, a serial C code, to the MTA with only a modest effort and with good performance. Our porting process needed neither a parallelism support platform nor attention to the intimate details of parallel programming and interprocessor communication, as would have been the case with more conventional supercomputers.

[1]  Thomas Stricker,et al.  Performance characterization of a molecular dynamics code on PC clusters: is there any easy parallelism in CHARMM? , 2002, Proceedings 16th International Parallel and Distributed Processing Symposium.

[2]  A. Parsegian,et al.  Energy of an Ion crossing a Low Dielectric Membrane: Solutions to Four Relevant Electrostatic Problems , 1969, Nature.

[3]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[4]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[5]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[6]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[7]  Yuefan Deng,et al.  The performance of a supercomputer built with commodity components , 2001, Parallel Comput..

[8]  David Fincham,et al.  Optimisation of the Ewald Sum for Large Systems , 1994 .

[9]  Allan Porterfield,et al.  The Tera computer system , 1990 .

[10]  H. G. Petersen,et al.  An algorithm for the simulation of condensed matter which grows as the 3/2 power of the number of particles , 1988 .

[11]  Juhyoun Kwak,et al.  Ion-beam sculpting at nanometre length scales , 2001 .

[12]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[13]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[14]  M G McInnis,et al.  Evidence for genetic linkage of Alzheimer's disease to chromosome 10q. , 2000, Science.

[15]  M. Muthukumar,et al.  Translocation of a confined polymer through a hole. , 2001, Physical review letters.

[16]  M. Glaser.  Atomistic Simulation and Modeling of Smectic Liquid Crystals , 2000 .

[17]  Grant S. Heffelfinger,et al.  Parallel atomistic simulations , 2000 .

[18]  J. Board,et al.  Ewald summation techniques in perspective: a survey , 1996 .

[19]  D. Branton,et al.  Characterization of individual polynucleotide molecules using a membrane channel. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Laxmikant V. Kalé,et al.  Scalable Molecular Dynamics for Large Biomolecular Systems , 2000, ACM/IEEE SC 2000 Conference (SC'00).

[21]  M. Muthukumar,et al.  Polymer translocation through a hole , 1999 .