No abstract

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This paper, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Next, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, e.g., the Tcl scripting language. Finally, the paper provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

No abstract

NAMD is a molecular dynamics program designed for high-performance simulations of very large biological objects on central processing unit (CPU)-and graphics processing unit (GPU)-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics, controlling the temperature, pressure and pH, applying external potentials on tailored grids, leveraging massively parallel resources in multiple-copy simulations, as well as hybrid QM/MM descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations, and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts towards achieving optimal performance on GPUbased architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

A large variety of animals has the ability to sense the geomagnetic field and utilize it as a source of directional (compass) information. It is not known by which biophysical mechanism this magnetoreception is achieved. We investigate the possibility that magnetoreception involves radical-pair processes that are governed by anisotropic hyperfine coupling between (unpaired) electron and nuclear spins. We will show theoretically that fields of geomagnetic field strength and weaker can produce significantly different reaction yields for different alignments of the radical pairs with the magnetic field. As a model for a magnetic sensory organ we propose a system of radical pairs being 1) orientationally ordered in a molecular substrate and 2) exhibiting changes in the reaction yields that affect the visual transduction pathway. We evaluate three-dimensional visual modulation patterns that can arise from the influence of the geomagnetic field on radical-pair systems. The variations of these patterns with orientation and field strength can furnish the magnetic compass ability of birds with the same characteristics as observed in behavioral experiments. We propose that the recently discovered photoreceptor cryptochrome is part of the magnetoreception system and suggest further studies to prove or disprove this hypothesis.