James C. Gumbart

Scalable Molecular Dynamics with NAMD

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.

Structural insight into the biogenesis of β-barrel membrane proteins

β-barrel membrane proteins are essential for nutrient import, signalling, motility and survival. In Gram-negative bacteria, the β-barrel assembly machinery (BAM) complex is responsible for the biogenesis of β-barrel membrane proteins, with homologous complexes found in mitochondria and chloroplasts. Here we describe the structure of BamA, the central and essential component of the BAM complex, from two species of bacteria: Neisseria gonorrhoeae and Haemophilus ducreyi. BamA consists of a large periplasmic domain attached to a 16-strand transmembrane β-barrel domain. Three structural features shed light on the mechanism by which BamA catalyses β-barrel assembly. First, the interior cavity is accessible in one BamA structure and conformationally closed in the other. Second, an exterior rim of the β-barrel has a distinctly narrowed hydrophobic surface, locally destabilizing the outer membrane. And third, the β-barrel can undergo lateral opening, suggesting a route from the interior cavity in BamA into the outer membrane.

Structure of Monomeric Yeast and Mammalian Sec61 Complexes Interacting with the Translating Ribosome

Nascent Chains Revealed Detailed analysis of protein translation and translocation across membranes requires the identification and structural analysis of intermediates involved in these processes (see the Perspective by Kampmann and Blobel). Seidelt et al. (p. 1412, published online 29 October) report the visualization by cryo-electron microscopy of a nascent polypeptide chain in the tunnel of the ribosome at 5.8 angstroms. This resolution allows analysis of the conformation and distinct contacts of the nascent chain within the ribosomal tunnel, which suggests a mechanism by which translational stalling is induced by this peptide. Protein translocation across cellular membranes involves the Sec61 protein, a component of a protein-conducting channel. Whether Sec61 acts as a monomer or as an oligomer during protein translocation has been unclear. Becker et al. (p. 1369, published online 29 October) describe active yeast and mammalian ribosome-Sec61 structures that show the Sec61 complex interacting with the ribosome and a nascent secretory protein signal sequence. The analysis unambiguously reveals that the active protein-conducting channel is a single Sec61 copy with its central pore serving as conduit for the nascent polypeptide. A single copy of a protein-conducting channel molecule provides a conduit for polypeptide translocation across membranes. The trimeric Sec61/SecY complex is a protein-conducting channel (PCC) for secretory and membrane proteins. Although Sec complexes can form oligomers, it has been suggested that a single copy may serve as an active PCC. We determined subnanometer-resolution cryo–electron microscopy structures of eukaryotic ribosome-Sec61 complexes. In combination with biochemical data, we found that in both idle and active states, the Sec complex is not oligomeric and interacts mainly via two cytoplasmic loops with the universal ribosomal adaptor site. In the active state, the ribosomal tunnel and a central pore of the monomeric PCC were occupied by the nascent chain, contacting loop 6 of the Sec complex. This provides a structural basis for the activity of a solitary Sec complex in cotranslational protein translocation.

Cryo-EM structure of the ribosome-SecYE complex in the membrane environment.

The ubiquitous SecY–Sec61 complex translocates nascent secretory proteins across cellular membranes and integrates membrane proteins into lipid bilayers. Several structures of mostly detergent-solubilized Sec complexes have been reported. Here we present a single-particle cryo-EM structure of the SecYEG complex in a membrane environment, bound to a translating ribosome, at subnanometer resolution. Using the SecYEG complex reconstituted in a so-called Nanodisc, we could trace the nascent polypeptide chain from the peptidyltransferase center into the membrane. The reconstruction allowed for the identification of ribosome–lipid interactions. The rRNA helix 59 (H59) directly contacts the lipid surface and appears to modulate the membrane in immediate vicinity to the proposed lateral gate of the protein-conducting channel (PCC). On the basis of our map and molecular dynamics simulations, we present a model of a signal anchor–gated PCC in the membrane.

Rapid parameterization of small molecules using the force field toolkit

The inability to rapidly generate accurate and robust parameters for novel chemical matter continues to severely limit the application of molecular dynamics 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, for example, General Amber Force Field and CHARMM General Force Field, 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, setup 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).

Molecular dynamics simulations of membrane channels and transporters.

Membrane transport constitutes one of the most fundamental processes in all living cells with proteins as major players. Proteins as channels provide highly selective diffusive pathways gated by environmental factors, and as transporters furnish directed, energetically uphill transport consuming energy. X-ray crystallography of channels and transporters furnishes a rapidly growing number of atomic resolution structures, permitting molecular dynamics (MD) simulations to reveal the physical mechanisms underlying channel and transporter function. Ever increasing computational power today permits simulations stretching up to 1 micros, that is, to physiologically relevant time scales. Membrane protein simulations presently focus on ion channels, on aquaporins, on protein-conducting channels, as well as on various transporters. In this review we summarize recent developments in this rapidly evolving field.

Structural basis for iron piracy by pathogenic Neisseria

Neisseria are obligate human pathogens causing bacterial meningitis, septicaemia and gonorrhoea. Neisseria require iron for survival and can extract it directly from human transferrin for transport across the outer membrane. The transport system consists of TbpA, an integral outer membrane protein, and TbpB, a co-receptor attached to the cell surface; both proteins are potentially important vaccine and therapeutic targets. Two key questions driving Neisseria research are how human transferrin is specifically targeted, and how the bacteria liberate iron from transferrin at neutral pH. To address these questions, we solved crystal structures of the TbpA–transferrin complex and of the corresponding co-receptor TbpB. We characterized the TbpB–transferrin complex by small-angle X-ray scattering and the TbpA–TbpB–transferrin complex by electron microscopy. Our studies provide a rational basis for the specificity of TbpA for human transferrin, show how TbpA promotes iron release from transferrin, and elucidate how TbpB facilitates this process.

Lateral opening and exit pore formation are required for BamA function.

The outer membrane of Gram-negative bacteria is replete with a host of β-barrel outer membrane proteins (OMPs). Despite serving a variety of essential functions, including immune response evasion, the exact mechanism of OMP folding and membrane insertion remains largely unclear. The β-barrel assembly machinery complex is required for OMP biogenesis. Crystal structures and molecular dynamics (MD) simulations of the central and essential component, BamA, suggest a mechanism involving lateral opening of its barrel domain. MD simulations reported here reveal an additional feature of BamA: a substrate exit pore positioned above the lateral opening site. Disulfide crosslinks that prevent lateral opening and exit pore formation result in a loss of BamA function, which can be fully rescued by the reductant tris(2-carboxyethyl)phosphine. These data provide strong evidence that lateral opening and exit pore formation are required for BamA function.

Structure of the SecY channel during initiation of protein translocation

Many secretory proteins are targeted by signal sequences to a protein-conducting channel, formed by prokaryotic SecY or eukaryotic Sec61 complexes, and are translocated across the membrane during their synthesis. Crystal structures of the inactive channel show that the SecY subunit of the heterotrimeric complex consists of two halves that form an hourglass-shaped pore with a constriction in the middle of the membrane and a lateral gate that faces the lipid phase. The closed channel has an empty cytoplasmic funnel and an extracellular funnel that is filled with a small helical domain, called the plug. During initiation of translocation, a ribosome–nascent chain complex binds to the SecY (or Sec61) complex, resulting in insertion of the nascent chain. However, the mechanism of channel opening during translocation is unclear. Here we have addressed this question by determining structures of inactive and active ribosome–channel complexes with cryo-electron microscopy. Non-translating ribosome–SecY channel complexes derived from Methanocaldococcus jannaschii or Escherichia coli show the channel in its closed state, and indicate that ribosome binding per se causes only minor changes. The structure of an active E. coli ribosome–channel complex demonstrates that the nascent chain opens the channel, causing mostly rigid body movements of the amino- and carboxy-terminal halves of SecY. In this early translocation intermediate, the polypeptide inserts as a loop into the SecY channel with the hydrophobic signal sequence intercalated into the open lateral gate. The nascent chain also forms a loop on the cytoplasmic surface of SecY rather than entering the channel directly.

The adaptive biasing force method: everything you always wanted to know but were afraid to ask.

In the host of numerical schemes devised to calculate free energy differences by way of geometric transformations, the adaptive biasing force algorithm has emerged as a promising route to map complex free-energy landscapes. It relies upon the simple concept that as a simulation progresses, a continuously updated biasing force is added to the equations of motion, such that in the long-time limit it yields a Hamiltonian devoid of an average force acting along the transition coordinate of interest. This means that sampling proceeds uniformly on a flat free-energy surface, thus providing reliable free-energy estimates. Much of the appeal of the algorithm to the practitioner is in its physically intuitive underlying ideas and the absence of any requirements for prior knowledge about free-energy landscapes. Since its inception in 2001, the adaptive biasing force scheme has been the subject of considerable attention, from in-depth mathematical analysis of convergence properties to novel developments and extensions. The method has also been successfully applied to many challenging problems in chemistry and biology. In this contribution, the method is presented in a comprehensive, self-contained fashion, discussing with a critical eye its properties, applicability, and inherent limitations, as well as introducing novel extensions. Through free-energy calculations of prototypical molecular systems, many methodological aspects are examined, from stratification strategies to overcoming the so-called hidden barriers in orthogonal space, relevant not only to the adaptive biasing force algorithm but also to other importance-sampling schemes. On the basis of the discussions in this paper, a number of good practices for improving the efficiency and reliability of the computed free-energy differences are proposed.