Oliver Beckstein

MDAnalysis: A Toolkit for the Analysis of Molecular Dynamics Simulations

MDAnalysis is an object‐oriented library for structural and temporal analysis of molecular dynamics (MD) simulation trajectories and individual protein structures. It is written in the Python language with some performance‐critical code in C. It uses the powerful NumPy package to expose trajectory data as fast and efficient NumPy arrays. It has been tested on systems of millions of particles. Many common file formats of simulation packages including CHARMM, Gromacs, Amber, and NAMD and the Protein Data Bank format can be read and written. Atoms can be selected with a syntax similar to CHARMMs powerful selection commands. MDAnalysis enables both novice and experienced programmers to rapidly write their own analytical tools and access data stored in trajectories in an easily accessible manner that facilitates interactive explorative analysis. MDAnalysis has been tested on and works for most Unix‐based platforms such as Linux and Mac OS X. It is freely available under the GNU General Public License from http://mdanalysis.googlecode.com.

Liquid–vapor oscillations of water in hydrophobic nanopores

Water plays a key role in biological membrane transport. In ion channels and water-conducting pores (aquaporins), one-dimensional confinement in conjunction with strong surface effects changes the physical behavior of water. In molecular dynamics simulations of water in short (0.8 nm) hydrophobic pores the water density in the pore fluctuates on a nanosecond time scale. In long simulations (460 ns in total) at pore radii ranging from 0.35 to 1.0 nm we quantify the kinetics of oscillations between a liquid-filled and a vapor-filled pore. This behavior can be explained as capillary evaporation alternating with capillary condensation, driven by pressure fluctuations in the water outside the pore. The free-energy difference between the two states depends linearly on the radius. The free-energy landscape shows how a metastable liquid state gradually develops with increasing radius. For radii > ≈0.55 nm it becomes the globally stable state and the vapor state vanishes. One-dimensional confinement affects the dynamic behavior of the water molecules and increases the self diffusion by a factor of 2–3 compared with bulk water. Permeabilities for the narrow pores are of the same order of magnitude as for biological water pores. Water flow is not continuous but occurs in bursts. Our results suggest that simulations aimed at collective phenomena such as hydrophobic effects may require simulation times >50 ns. For water in confined geometries, it is not possible to extrapolate from bulk or short time behavior to longer time scales.

Molecular Basis of Alternating Access Membrane Transport by the Sodium-Hydantoin Transporter Mhp1

Triangulating to Mechanism Cellular uptake and release of a variety of substrates are mediated by secondary transporters, but no crystal structures are known for all three fundamental states of the transport cycle, which has limited explanations for their proposed mechanisms. Shimamura et al. (p. 470) report a 3.8-angstrom structure of the inward-facing conformation of the bacterial sodium-benzylhydantoin transport protein, Mhp1, complementing the other two available structures. Molecular modeling for the interconversions of these structures shows a simple rigid body rotation of four helices relative to the rest of the structure in which the protein switches reversibly from outward- to inward-facing. Three complementary crystal structures reveal the mechanism of a transport protein in molecular dynamics simulations. The structure of the sodium-benzylhydantoin transport protein Mhp1 from Microbacterium liquefaciens comprises a five-helix inverted repeat, which is widespread among secondary transporters. Here, we report the crystal structure of an inward-facing conformation of Mhp1 at 3.8 angstroms resolution, complementing its previously described structures in outward-facing and occluded states. From analyses of the three structures and molecular dynamics simulations, we propose a mechanism for the transport cycle in Mhp1. Switching from the outward- to the inward-facing state, to effect the inward release of sodium and benzylhydantoin, is primarily achieved by a rigid body movement of transmembrane helices 3, 4, 8, and 9 relative to the rest of the protein. This forms the basis of an alternating access mechanism applicable to many transporters of this emerging superfamily.

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 Å) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (∼1 Å) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge‐bending of the pore‐lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.

Lipidbook: A Public Repository for Force-Field Parameters Used in Membrane Simulations

Lipidbook is a public database for force-field parameters with a special emphasis on lipids, detergents and similar molecules that are of interest when simulating biological membrane systems. It stores parameter files that are supplied by the community. Topologies, parameters and lipid or whole bilayer structures can be deposited in any format for any simulation code, preferably under a license that promotes “open knowledge.” A number of mechanisms are implemented to aid a user in judging the appropriateness of a given parameter set for a project. For instance, parameter sets are versioned, linked to the primary citation via PubMed identifier and digital object identifier (DOI), and users can publicly comment on deposited parameters. Licensing and, hence, the conditions for use and dissemination of academically generated data are often unclear. In those cases it is also possible to provide a link instead of uploading a file. A snapshot of the linked file is then archived using the WebCite® service without further involvement of the user or Lipidbook, thus ensuring a transparent and permanent history of the parameter set. Lipidbook can be accessed freely online at http://lipidbook.bioch.ox.ac.uk. Deposition of data requires online registration.

Simulations of anion transport through OprP reveal the molecular basis for high affinity and selectivity for phosphate

The outer membrane protein OprP from Pseudomonas aeruginosa forms a phosphate selective pore. To understand the mechanism of phosphate permeation and selectivity, we used three simulation techniques [equilibrium molecular dynamics simulations, steered molecular dynamics, and calculation of a potential of mean force (PMF)]. The PMF for phosphate reveals a deep free energy well midway along the OprP channel. Two adjacent phosphate-binding sites (W1 and W2), each with a well depth of ≈8 kT, are identified close to the L3 loop in the most constricted region of the pore. A dissociation constant for phosphate of 6 μM is computed from the PMF, within the range of reported experimental values. The transfer of phosphate between sites W1 and W2 is correlated with changes in conformation of the sidechain of K121, which serves as a “charged brush” to facilitate phosphate passage between the two subsites. OprP also binds chloride, but less strongly than phosphate, as calculated from a Cl− PMF. The difference in affinity and hence selectivity is due to the “Lys-cluster” motif, the positive charges of which interact strongly with a partially dehydrated phosphate ion but are shielded from a Cl− by the hydration shell of the smaller ion. Our simulations suggest that OprP does not conform to the conventional picture of a channel with relatively flat energy landscape for permeant ions, but rather resembles a membrane-inserted binding protein with a high specificity that allows access to a centrally located binding site from both the extracellular and the periplasmic spaces.

Computing ensembles of transitions from stable states: Dynamic importance sampling

There is an increasing dataset of solved biomolecular structures in more than one conformation and increasing evidence that large‐scale conformational change is critical for biomolecular function. In this article, we present our implementation of a dynamic importance sampling (DIMS) algorithm that is directed toward improving our understanding of important intermediate states between experimentally defined starting and ending points. This complements traditional molecular dynamics methods where most of the sampling time is spent in the stable free energy wells defined by these initial and final points. As such, the algorithm creates a candidate set of transitions that provide insights for the much slower and probably most important, functionally relevant degrees of freedom. The method is implemented in the program CHARMM and is tested on six systems of growing size and complexity. These systems, the folding of Protein A and of Protein G, the conformational changes in the calcium sensor S100A6, the glucose–galactose‐binding protein, maltodextrin, and lactoferrin, are also compared against other approaches that have been suggested in the literature. The results suggest good sampling on a diverse set of intermediates for all six systems with an ability to control the bias and thus to sample distributions of trajectories for the analysis of intermediate states.

Flexible Gates Generate Occluded Intermediates in the Transport Cycle of LacY

The major facilitator superfamily (MFS) transporter lactose permease (LacY) alternates between cytoplasmic and periplasmic open conformations to co-transport a sugar molecule together with a proton across the plasma membrane. Indirect experimental evidence suggested the existence of an occluded transition intermediate of LacY, which would prevent leaking of the proton gradient. As no experimental structure is known, the conformational transition is not fully understood in atomic detail. We simulated transition events from a cytoplasmic open conformation to a periplasmic open conformation with the dynamic importance sampling molecular dynamics method and observed occluded intermediates. Analysis of water permeation pathways and the electrostatic free-energy landscape of a solvated proton indicated that the occluded state contains a solvated central cavity inaccessible from either side of the membrane. We propose a pair of geometric order parameters that capture the state of the pathway through the MFS transporters as shown by a survey of available crystal structures and models. We present a model for the occluded state of apo-LacY, which is similar to the occluded crystal structures of the MFS transporters EmrD, PepTSo, NarU, PiPT and XylE. Our simulations are consistent with experimental double electron spin–spin distance measurements that have been interpreted to show occluded conformations. During the simulations, a salt bridge that has been postulated to be involved in driving the conformational transition formed. Our results argue against a simple rigid-body domain motion as implied by a strict “rocker-switch mechanism” and instead hint at an intricate coupling between two flexible gates.

Crystal structure of the sodium–proton antiporter NhaA dimer and new mechanistic insights

A dimeric structure of the sodium–proton antiporter NhaA provides insight into the roles of Asp163 and Lys300 in the transport mechanism.

The alternating access mechanism of transport as observed in the sodium-hydantoin transporter Mhp1

Crystal structures of a membrane protein transporter in three different conformational states provide insights into the transport mechanism.