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 article, 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. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also 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.

/pdf/scalable-molecular-dynamics-with-namd-51o1z7rwp8.pdf

Scalable molecular dynamics with NAMD

/pdf/scalable-molecular-dynamics-with-namd-3j7xyc70g4.pdf

Scalable Molecular Dynamics with NAMD

One of the most pervasive problems afflicting people throughout the world is inadequate access to clean water and sanitation. Problems with water are expected to grow worse in the coming decades, with water scarcity occurring globally, even in regions currently considered water-rich. Addressing these problems calls out for a tremendous amount of research to be conducted to identify robust new methods of purifying water at lower cost and with less energy, while at the same time minimizing the use of chemicals and impact on the environment. Here we highlight some of the science and technology being developed to improve the disinfection and decontamination of water, as well as efforts to increase water supplies through the safe re-use of wastewater and efficient desalination of sea and brackish water.

Science and technology for water purification in the coming decades

Unlike other methods for docking ligands to the rigid 3D structure of a known protein receptor, Glide approximates a complete systematic search of the conformational, orientational, and positional space of the docked ligand In this search, an initial rough positioning and scoring phase that dramatically narrows the search space is followed by torsionally flexible energy optimization on an OPLS-AA nonbonded potential grid for a few hundred surviving candidate poses The very best candidates are further refined via a Monte Carlo sampling of pose conformation; in some cases, this is crucial to obtaining an accurate docked pose Selection of the best docked pose uses a model energy function that combines empirical and force-field-based terms Docking accuracy is assessed by redocking ligands from 282 cocrystallized PDB complexes starting from conformationally optimized ligand geometries that bear no memory of the correctly docked pose Errors in geometry for the top-ranked pose are less than 1 A in nearly ha

/pdf/glide-a-new-approach-for-rapid-accurate-docking-and-scoring-lnxrmd5zbm.pdf

Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy.

A novel scoring function to estimate protein-ligand binding affinities has been developed and implemented as the Glide 4.0 XP scoring function and docking protocol. In addition to unique water desolvation energy terms, protein-ligand structural motifs leading to enhanced binding affinity are included: (1) hydrophobic enclosure where groups of lipophilic ligand atoms are enclosed on opposite faces by lipophilic protein atoms, (2) neutral-neutral single or correlated hydrogen bonds in a hydrophobically enclosed environment, and (3) five categories of charged-charged hydrogen bonds. The XP scoring function and docking protocol have been developed to reproduce experimental binding affinities for a set of 198 complexes (RMSDs of 2.26 and 1.73 kcal/mol over all and well-docked ligands, respectively) and to yield quality enrichments for a set of fifteen screens of pharmaceutical importance. Enrichment results demonstrate the importance of the novel XP molecular recognition and water scoring in separating active and inactive ligands and avoiding false positives.

/pdf/extra-precision-glide-docking-and-scoring-incorporating-a-2o3e7pl9pg.pdf

Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein-Ligand Complexes

Confinement of matter on the nanometre scale can induce phase transitions not seen in bulk systems1. In the case of water, so-called drying transitions occur on this scale2,3,4,5 as a result of strong hydrogen-bonding between water molecules, which can cause the liquid to recede from nonpolar surfaces to form a vapour layer separating the bulk phase from the surface6. Here we report molecular dynamics simulations showing spontaneous and continuous filling of a nonpolar carbon nanotube with a one-dimensionally ordered chain of water molecules. Although the molecules forming the chain are in chemical and thermal equilibrium with the surrounding bath, we observe pulse-like transmission of water through the nanotube. These transmission bursts result from the tight hydrogen-bonding network inside the tube, which ensures that density fluctuations in the surrounding bath lead to concerted and rapid motion along the tube axis7,8,9. We also find that a minute reduction in the attraction between the tube wall and water dramatically affects pore hydration, leading to sharp, two-state transitions between empty and filled states on a nanosecond timescale. These observations suggest that carbon nanotubes, with their rigid nonpolar structures10,11, might be exploited as unique molecular channels for water and protons, with the channel occupancy and conductivity tunable by changes in the local channel polarity and solvent conditions.

Water conduction through the hydrophobic channel of a carbon nanotube

Simulations of simple solutes (charged and uncharged spheres) in model water have been performed in order to elucidate aspects of solvation in water at ambient and supercriti- cal states. The variation of solvation entropy as a function of solute charge has been used to investigate hydrophobic and hydrophilic ordering and the structure-making and structure- breaking effects of ions. Simulations with model solvents, which differ from water in certain features, have been used to try to identify the particular properties of water that are associat- ed with these phenomena.

/pdf/using-simulation-to-study-solvation-in-water-159q54mzo5.pdf

J. P. Noworyta

Papers

Water conduction through the hydrophobic channel of a carbon nanotube

Using simulation to study solvation in water