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.

The inability to rapidly generate accurate and robust parameters for novel chemical matter continues to severely limit the application of molecular dynamics (MD) simulations to many biological systems of interest, especially in fields such as drug discovery. Although the release of generalized versions of common classical force fields, e.g., GAFF and CGenFF, have posited guidelines for parameterization of small molecules, many technical challenges remain that have hampered their wide-scale extension. The Force Field Toolkit (ffTK), described herein, minimizes common barriers to ligand parameterization through algorithm and method development, automation of tedious and error-prone tasks, and graphical user interface design. Distributed as a VMD plugin, ffTK facilitates the traversal of a clear and organized workflow resulting in a complete set of CHARMM-compatible parameters. A variety of tools are provided to generate quantum mechanical target data, set up multidimensional optimization routines, and analyze parameter performance. Parameters developed for a small test set of molecules using ffTK were comparable to existing CGenFF parameters in their ability to reproduce experimentally measured values for pure-solvent properties (<15% error from experiment) and free energy of solvation (±0.5 kcal/mol from experiment).

Aquaglyceroporin GlpF selectively conducts water and linear polyalcohols, such as glycerol, across the inner membrane of Escherichia coli. We report steered molecular dynamics simulations of glycerol conduction through GlpF, in which external forces accelerate the transchannel conduction in a manner that preserves the intrinsic conduction mechanism. The simulations reveal channel-glycerol hydrogen bonding interactions and the stereoselectivity of the channel. Employing Jarzynski's identity between free energy and irreversible work, we reconstruct the potential of mean force along the conduction pathway through a time series analysis of molecular dynamics trajectories. This potential locates binding sites and barriers inside the channel; it also reveals a low energy periplasmic vestibule suited for efficient uptake of glycerol from the environment.A quaporins (1), a family of water transporting membrane proteins, are present in all life forms, and defects in their function cause physiological disorders (2). Among more than 150 members identified to date (2), the Escherichia coli glycerol facilitator (GlpF) belongs to the aquaglyceroporin subclass, which is permeable to both water and glycerol. GlpF also stereoselectively conducts longer linear polyalcohols (3, 4). At physiological conditions, all aquaporins exclude charged solutes, including protons, and thereby preserve the electrochemical potential across the cell membrane (2).A 2.2-Å resolution x-ray structure of GlpF revealed a homotetrameric architecture with glycerol and water present inside the channel (5). Each monomeric channel has two halfmembrane spanning repeats related by a quasi-two-fold symmetry. About half of each repeat is ␣-helical; the other half adopts a particular nonhelical structure (5, 6). The N termini of the helical repeats meet at the Asn-Pro-Ala (NPA) motifs located at the channel center. The NPA motifs are conserved among all aquaporins (2), and their spatial arrangement is critical for the biological function of the channel (6-8). The nonhelical repeats expose the backbone carbonyl groups of residues 64-66 and 195-201 toward the channel interior, where they serve as hydrogen acceptors for the substrate (6). The channel diameter measures less than 3.5 Å at its narrowest point, the selectivity filter (SF), lined with residues Trp-48, 9). In the constriction region of the channel, approximately 25 Å long, the substrate is translocated in a single file following a curvilinear pathway (6), a feature critical for excluding proton conduction (10). The hydroxyl groups of glycerol make hydrogen bonds with exposed carbonyl oxygen atoms, polar hydrogen atoms, and water, whereas the aliphatic backbone of glycerol faces the opposite hydrophobic side of the channel (6). The amphipathic channel lining in GlpF ensures the selectivity for linear polyalcohols (5). The amphipathicity of the channel's interior also seems to be important for water transport and is found, although to a less extent, in pure water-conducting aquaporins, such as aquaporin-1 ...

The mechanism of color tuning in the rhodopsin family of proteins has been studied by comparing the optical properties of the light-driven proton pump bacteriorhodopsin (bR) and the light detector sensory rhodopsin II (sRII). Despite a high structural similarity, the maximal absorption is blue-shifted from 568 nm in bR to 497 nm in sRII. The molecular mechanism of this shift is still a matter of debate, and its clarification sheds light onto the general mechanisms of color tuning in retinal proteins. The calculations employ a combined quantum mechanical/molecular mechanical (QM/MM) technique, using a DFT-based method for ground state properties and the semiempirical OM2/MRCI method and ab initio SORCI method for excited state calculations. The high efficiency of the methodology has allowed us to study a wide variety of aspects including dynamical effects. The absorption shift as well as various mutation experiments and vibrational properties have been successfully reproduced. Our results indicate that several sources contribute to the spectral shift between bR and sRII. The main factors are the counterion region at the extracellular side of retinal and the amino acid composition of the binding pocket. Our analysis allows a distinction and identification of the different effects in detail and leads to a clear picture of the mechanism of color tuning, which is in good agreement with available experimental data.