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.

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.

Jarzynski's equality is applied to free energy calculations from steered molecular dynamics simulations of biomolecules. The helix-coil transition of deca-alanine in vacuum is used as an example. With about ten trajectories sampled, the second order cumulant expansion, among the various averaging schemes examined, yields the most accurate estimates. We compare umbrella sampling and the present method, and find that their efficiencies are comparable.

Aquaporins are transmembrane channels found in cell membranes of all life forms. We examine their apparently paradoxical property, facilitation of efficient permeation of water while excluding protons, which is of critical importance to preserving the electrochemical potential across the cell membrane. We have determined the structure of the Escherichia coli aquaglyceroporin GlpF with bound water, in native (2.7 angstroms) and in W48F/F200T mutant (2.1 angstroms) forms, and carried out 12-nanosecond molecular dynamics simulations that define the spatial and temporal probability distribution and orientation of a single file of seven to nine water molecules inside the channel. Two conserved asparagines force a central water molecule to serve strictly as a hydrogen bond donor to its neighboring water molecules. Assisted by the electrostatic potential generated by two half-membrane spanning loops, this dictates opposite orientations of water molecules in the two halves of the channel, and thus prevents the formation of a "proton wire," while permitting rapid water diffusion. Both simulations and observations revealed a more regular distribution of channel water and an increased water permeability for the W48F/F200T mutant.

Plants counteract fluctuations in water supply by regulating all aquaporins in the cell plasma membrane. Channel closure results either from the dephosphorylation of two conserved serine residues under conditions of drought stress, or from the protonation of a conserved histidine residue following a drop in cytoplasmic pH due to anoxia during flooding. Here we report the X-ray structure of the spinach plasma membrane aquaporin SoPIP2;1 in its closed conformation at 2.1 A resolution and in its open conformation at 3.9 A resolution, and molecular dynamics simulations of the initial events governing gating. In the closed conformation loop D caps the channel from the cytoplasm and thereby occludes the pore. In the open conformation loop D is displaced up to 16 A and this movement opens a hydrophobic gate blocking the channel entrance from the cytoplasm. These results reveal a molecular gating mechanism which appears conserved throughout all plant plasma membrane aquaporins.