NAMD: biomolecular simulation on thousands of processors

Abstract NAMD is a fully featured, production molecular dynamics program for high performance simulation of large biomolecular systems. We have previously, at SC2000, presented scaling results for simulations with cutoff electrostatics on up to 2048 processors of the ASCI Red machine, achieved with an object-based hybrid force and spatial decomposition scheme and an aggressive measurement-based predictive load balancing framework. We extend this work by demonstrating similar scaling on the much faster processors of the PSC Lemieux Alpha cluster, and for simulations employing efficient (order N log N) particle mesh Ewald full electrostatics. This unprecedented scalability in a biomolecular simulation code has been attained through latency tolerance, adaptation to multiprocessor nodes, and the direct use of the Quadrics Elan library in place of MPI by the Charm++/Converse parallel runtime system

Argobots: A Lightweight Low-Level Threading and Tasking Framework

In the past few decades, a number of user-level threading and tasking models have been proposed in the literature to address the shortcomings of OS-level threads, primarily with respect to cost and flexibility. Current state-of-the-art user-level threading and tasking models, however, either are too specific to applications or architectures or are not as powerful or flexible. In this paper, we present Argobots, a lightweight, low-level threading and tasking framework that is designed as a portable and performant substrate for high-level programming models or runtime systems. Argobots offers a carefully designed execution model that balances generality of functionality with providing a rich set of controls to allow specialization by end users or high-level programming models. We describe the design, implementation, and performance characterization of Argobots and present integrations with three high-level models: OpenMP, MPI, and colocated I/O services. Evaluations show that (1) Argobots, while providing richer capabilities, is competitive with existing simpler generic threading runtimes; (2) our OpenMP runtime offers more efficient interoperability capabilities than production OpenMP runtimes do; (3) when MPI interoperates with Argobots instead of Pthreads, it enjoys reduced synchronization costs and better latency-hiding capabilities; and (4) I/O services with Argobots reduce interference with colocated applications while achieving performance competitive with that of a Pthreads approach

Excitons in a Photosynthetic Light-Harvesting System: A Combined Molecular Dynamics/Quantum Chemistry and Polaron Model Study

The dynamics of pigment-pigment and pigment-protein interactions in light-harvesting complexes is studied with a novel approach which combines molecular dynamics (MD) simulations with quantum chemistry (QC) calculations. The MD simulations of an LH-II complex, solvated and embedded in a lipid bilayer at physiological conditions (with total system size of 87,055 atoms) revealed a pathway of a water molecule into the B800 binding site, as well as increased dimerization within the B850 BChl ring, as compared to the dimerization found for the crystal structure. The fluctuations of pigment (B850 BChl) excitation energies, as a function of time, were determined via ab initio QC calculations based on the geometries that emerged from the MD simulations. From the results of these calculations we constructed a time-dependent Hamiltonian of the B850 exciton system from which we determined the linear absorption spectrum. Finally, a polaron model is introduced to describe quantum mechanically both the excitonic and vibrational (phonon) degrees of freedom. The exciton-phonon coupling that enters into the polaron model, and the corresponding phonon spectral function are derived from the MD/QC simulations. It is demonstrated that, in the framework of the polaron model, the absorption spectrum of the B850 excitons can be calculated from the autocorrelation function of the excitation energies of individual BChls, which is readily available from the combined MD/QC simulations. The obtained result is in good agreement with the experimentally measured absorption spectrum.Comment: REVTeX3.1, 23 pages, 13 (EPS) figures included. A high quality PDF file of the paper is available at http://www.ks.uiuc.edu/Publications/Papers/PDF/DAMJ2001/DAMJ2001.pd

Relating data—parallelism and (and—) parallelism in logic programs

Much work has been done in the áreas of and-parallelism and data parallelism in Logic Programs. Such work has proceeded to a certain extent in an independent fashion. Both types of parallelism offer advantages and disadvantages. Traditional (and-) parallel models offer generality, being able to exploit parallelism in a large class of programs (including that exploited by data parallelism techniques). Data parallelism techniques on the other hand offer increased performance for a restricted class of programs. The thesis of this paper is that these two forms of parallelism are not fundamentally different and that relating them opens the possibility of obtaining the advantages of both within the same system. Some relevant issues are discussed and solutions proposed. The discussion is illustrated through visualizations of actual parallel executions implementing the ideas proposed

Direct observation of topoisomerase IA gate dynamics

Type IA topoisomerases cleave single-stranded DNA and relieve negative supercoils in discrete steps corresponding to the passage of the intact DNA strand through the cleaved strand. Although type IA topoisomerases are assumed to accomplish this strand passage via a protein-mediated DNA gate, opening of this gate has never been observed. We developed a single-molecule assay to directly measure gate opening of the Escherichia coli type IA topoisomerases I and III. We found that after cleavage of single-stranded DNA, the protein gate opens by as much as 6.6 nm and can close against forces in excess of 16 pN. Key differences in the cleavage, ligation, and gate dynamics of these two enzymes provide insights into their different cellular functions. The single-molecule results are broadly consistent with conformational changes obtained from molecular dynamics simulations. These results allowed us to develop a mechanistic model of interactions between type IA topoisomerases and single-stranded DNA

Mechanisms of Hybrid Oligomer Formation in the Pathogenesis of Combined Alzheimer's and Parkinson's Diseases

Background: Misfolding and pathological aggregation of neuronal proteins has been proposed to play a critical role in the pathogenesis of neurodegenerative disorders. Alzheimer’s disease (AD) and Parkinson’s disease (PD) are frequent neurodegenerative diseases of the aging population. While progressive accumulation of amyloid b protein (Ab) oligomers has been identified as one of the central toxic events in AD, accumulation of a-synuclein (a-syn) resulting in the formation of oligomers and protofibrils has been linked to PD and Lewy body Disease (LBD). We have recently shown that Ab promotes a-syn aggregation and toxic conversion in vivo, suggesting that abnormal interactions between misfolded proteins might contribute to disease pathogenesis. However the molecular characteristics and consequences of these interactions are not completely clear. Methodology/Principal Findings: In order to understand the molecular mechanisms involved in potential Ab/a-syn interactions, immunoblot, molecular modeling, and in vitro studies with a-syn and Ab were performed. We showed in vivo in the brains of patients with AD/PD and in transgenic mice, Ab and a-synuclein co-immunoprecipitate and form complexes. Molecular modeling and simulations showed that Ab binds a-syn monomers, homodimers, and trimers, forming hybrid ringlike pentamers. Interactions occurred between the N-terminus of Ab and the N-terminus and C-terminus of a-syn. Interacting a-syn and Ab dimers that dock on the membrane incorporated additional a-syn molecules, leading to th