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 CHARMM’s 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 Public License from http://mdanalysis.googlecode.com.

MDAnalysis (http://mdanalysis.org) is a library for structural and temporal analysis of molecular dynamics (MD) simulation trajectories and individual protein structures. MD simulations of biological molecules have become an important tool to elucidate the relationship between molecular structure and physiological function. Simulations are performed with highly optimized software packages on HPC resources but most codes generate output trajectories in their own formats so that the development of new trajectory analysis algorithms is confined to specific user communities and widespread adoption and further development is delayed. MDAnalysis addresses this problem by abstracting access to the raw simulation data and presenting a uniform object-oriented Python interface to the user. It thus enables users to rapidly write code that is portable and immediately usable in virtually all biomolecular simulation communities. The user interface and modular design work equally well in complex scripted work flows, as foundations for other packages, and for interactive and rapid prototyping work in IPython / Jupyter notebooks, especially together with molecular visualization provided by nglview and time series analysis with pandas. MDAnalysis is written in Python and Cython and uses NumPy arrays for easy interoperability with the wider scientific Python ecosystem. It is widely used and forms the foundation for more specialized biomolecular simulation tools. MDAnalysis is available under the GNU General Public License v2.

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 freeenergy 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. C hannel and transporter proteins control flow of water, ions, and other solutes across cell membranes. In recent years several channel and pore structures have been solved at near atomic resolution (1-6), which together with three decades of physiological data (7) and theoretical and simulation approaches (8) allow us to describe transport of ions, water, or other small molecules at a molecular level. Water plays a special role here: it either solvates the inner surfaces of the pore and the permeators (for example, ions and small molecules like glycerol), or it is the permeant species itself as in the aquaporin family of water pores (9-11) or in the bacterial peptide channel gramicidin A, whose water transport properties are well studied (12-14). Thus, a better characterization of the behavior of water would improve our understanding of the biological function of a wide range of transporters. The remarkable water transport properties of aquaporins [water is conducted through a long (Ϸ2 nm) and narrow (Ϸ0.3 nm diameter) pore at bulk diffusion rates while at the same time protons are strongly selected against] are the topic of recent simulation studies (15, 16).The shape and dimensions of biological pores and the nature of the pore-lining atoms are recognized as major determinants of function. How the behavior of water depends on these factors is far from understood (17). Water is not a simple liquid because of its strong hydrogen bond network. When confined to narrow geometries like slits or pores ...

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.

Water-filled pores of nanometer dimensions play important roles in chemistry and biology, e.g., as channels through biological membranes. Biological nanopores are frequently gated, i.e., they switch between an open and a closed state. In several ion channel structures the gate is formed by a ring of hydrophobic side chains that do not physically occlude the pore. Here we investigate whether a hydrophobic pore can act as a gate via molecular dynamics simulations of the passage of water through atomistic models of nanopores embedded within a membrane mimetic. Both the geometry of a nanopore and the hydrophilicity vs hydrophobicity of its lining determine whether water enters the channel. For purely hydrophobic pores there is an abrupt transition from a closed state (no water in the pore cavity) to an open state (cavity water at approximately bulk density) once a critical pore radius is exceeded. This critical radius depends on the length of the pore and the radius of the mouth region. Furthermore, a closed hydrophobic nanopore can be opened by adding dipoles to its lining.

A hydrophobic constriction site can act as an efficient barrier to ion and water permeation if its diameter is less than the diameter of an ion's first hydration shell. This hydrophobic gating mechanism is thought to operate in a number of ion channels, e.g. the nicotinic receptor, bacterial mechanosensitive channels (MscL and MscS) and perhaps in some potassium channels (e.g. KcsA, MthK and KvAP). Simplified pore models allow one to investigate the primary characteristics of a conduction pathway, namely its geometry (shape, pore length, and radius), the chemical character of the pore wall surface, and its local flexibility and surface roughness. Our extended (about 0.1 micros) molecular dynamic simulations show that a short hydrophobic pore is closed to water for radii smaller than 0.45 nm. By increasing the polarity of the pore wall (and thus reducing its hydrophobicity) the transition radius can be decreased until for hydrophilic pores liquid water is stable down to a radius comparable to a water molecule's radius. Ions behave similarly but the transition from conducting to non-conducting pores is even steeper and occurs at a radius of 0.65 nm for hydrophobic pores. The presence of water vapour in a constriction zone indicates a barrier for ion permeation. A thermodynamic model can explain the behaviour of water in nanopores in terms of the surface tensions, which leads to a simple measure of 'hydrophobicity' in this context. Furthermore, increased local flexibility decreases the permeability of polar species. An increase in temperature has the same effect, and we hypothesize that both effects can be explained by a decrease in the effective solvent-surface attraction which in turn leads to an increase in the solvent-wall surface free energy.

Sodium/proton (Na+/H+) antiporters, located at the plasma membrane in every cell, are vital for cell homeostasis1. In humans, their dysfunction has been linked to diseases, such as, hypertension, heart failure and epilepsy and they are well-established drug targets2. The best understood model system for Na+/H+ antiport is NhaA from Escherichia coli1,3, where both EM and crystal structures are available4-6. NhaA is made up of two distinct domains, a Core domain and a Dimerisation domain. In the NhaA crystal structure a cavity is located between the two domains providing access to the ion-binding site from the inward-facing surface of the protein1,4. Like many Na+/H+ antiporters, the activity of NhaA is regulated by pH, only becoming active above pH 6.5, where a conformational change is thought to occur7. To date, the only reported NhaA crystal structure is of the low pH inactivated form4. Here, we describe the active-state structure of a Na+/H+ antiporter, NapA from Thermus thermophilus at 3 Å resolution, solved from crystals grown at pH 7.8. In the NapA structure, the Core and Dimerisation domains are in different positions to those seen in NhaA and a negatively charged cavity has now opened to the outside. The extracellular cavity allows access to a strictly conserved aspartate residue thought to directly coordinate ion-binding1,8,9, a role supported here by molecular dynamics simulations. To alternate access to this ion-binding site, however, requires a surprisingly large rotation of the Core domain, some 20° against the Dimerisation interface. We conclude that despite their fast transport rates of up to 1500 ions/sec3, Na+/H+ antiporters operate by a two-domain rocking bundle model, revealing themes relevant to secondary-active transporters in general.

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free energy calculations and our new dynamic importance sampling (DIMS) molecular dynamics (MD) method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular four ionic bonds are identified that open in a sequential, zipper-like fashion and thus dominate the free energy landscape of the transition. Transitions between the closed and open conformations only have to overcome moderate free energy barriers. Unexpectedly, the closed and open state encompass broad free energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental FRET measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS-MD computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.