NAMD is a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and 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; and hybrid quantum-mechanical/molecular-mechanical 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 toward achieving optimal performance on GPU-based 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.

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.

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.