Exploring transmembrane transport through alpha-hemolysin with grid-steered molecular dynamics

Determination of Free Energy Profiles for the Translocation of Polynucleotides through α-Hemolysin Nanopores using Non-Equilibrium Molecular Dynamics Simulations

The translocation of polynucleotides through transmembrane protein pores is a fundamental biological process with important technological and medical relevance. The translocation process is complex, and it is influenced by a range of factors including the diameter and inner surface of the pore, the secondary structure of the polymer, and the interactions between the polymer and protein. In this paper, we perform nonequilibrium constant velocity-steered molecular dynamics simulations of nucleic acid molecule translocation through the protein nanopore α-hemolysin and use Jarzynski's identity to determine the associated free energy profiles. With this approach we are able to explain the observed differences in experimental translocation time through the nanopore between polyadenosine and polydeoxycytidine. The translocation of polynucleotides and single nucleotides through α-hemolysin is investigated. These simulations are computationally intensive as they employ models with atomistic level resolution; in addition to their size, these systems are challenging to study due to the time scales of translocation of large asymmetric molecules. Our simulations provide insight into the role of the interactions between the nucleic acid molecules and the protein pore. Mutated protein pores provide confirmation of residue-specific interactions between nucleotides and the protein pore. By harnessing such molecular dynamics simulations, we gain new physicochemical insight into the translocation process.

Plasmonic Nanopores for Trapping, Controlling Displacement, and Sequencing of DNA

With the aim of developing a DNA sequencing methodology, we theoretically examine the feasibility of using nanoplasmonics to control the translocation of a DNA molecule through a solid-state nanopore and to read off sequence information using surface-enhanced Raman spectroscopy. Using molecular dynamics simulations, we show that high-intensity optical hot spots produced by a metallic nanostructure can arrest DNA translocation through a solid-state nanopore, thus providing a physical knob for controlling the DNA speed. Switching the plasmonic field on and off can displace the DNA molecule in discrete steps, sequentially exposing neighboring fragments of a DNA molecule to the pore as well as to the plasmonic hot spot. Surface-enhanced Raman scattering from the exposed DNA fragments contains information about their nucleotide composition, possibly allowing the identification of the nucleotide sequence of a DNA molecule transported through the hot spot. The principles of plasmonic nanopore sequencing can be extended to detection of DNA modifications and RNA characterization.

All-Atom Molecular Dynamics Simulation of Protein Translocation through an α-Hemolysin Nanopore

Nanopore sensing is attracting the attention of a large and varied scientific community. One of the main issues in nanopore sensing is how to associate the measured current signals to specific features of the molecule under investigation. This is particularly relevant when the translocating molecule is a protein and the pore is sufficiently narrow to necessarily involve unfolding of the translocating protein. Recent experimental results characterized the cotranslocational unfolding of Thioredoxin (Trx) passing through an α-hemolisin pore, providing evidence for the existence of a multistep process. In this study we report the results of all-atom molecular dynamics simulations of the same system. Our data indicate that Trx translocation involves two main barriers. The first one is an unfolding barrier associated with a translocation intermediate where the N-terminal region of Trx is stuck at the pore entrance in a conformation that strongly resembles the native one. After the abrupt unfolding of the N-terminal region, the Trx enters the α-hemolisin vestibule. During this stage, the constriction is occupied not only by the translocating residue but also by a hairpin-like structure forming a tangle in the constriction. The second barrier is associated with the disentangling of this region.

Dendrimers in Nanoscale Confinement: The Interplay between Conformational Change and Nanopore Entrance

Hyperbranched dendrimers are nanocarriers for drugs, imaging agents, and catalysts. Their nanoscale confinement is of fundamental interest and occurs when dendrimers with bioactive payload block or pass biological nanochannels or when catalysts are entrapped in inorganic nanoporous support scaffolds. The molecular process of confinement and its effect on dendrimer conformations are, however, poorly understood. Here, we use single-molecule nanopore measurements and molecular dynamics simulations to establish an atomically detailed model of pore dendrimer interactions. We discover and explain that electrophoretic migration of polycationic PAMAM dendrimers into confined space is not dictated by the diameter of the branched molecules but by their size and generation-dependent compressibility. Differences in structural flexibility also rationalize the apparent anomaly that the experimental nanopore current read-out depends in nonlinear fashion on dendrimer size. Nanoscale confinement is inferred to reduce the protonation of the polycationic structures. Our model can likely be expanded to other dendrimers and be applied to improve the analysis of biophysical experiments, rationally design functional materials such as nanoporous filtration devices or nanoscale drug carriers that effectively pass biological pores.

Dual mechanism of ion permeation through VDAC revealed with inorganic phosphate ions and phosphate metabolites

In the exchange of metabolites and ions between the mitochondrion and the cytosol, the voltage-dependent anion channel (VDAC) is a key element, as it forms the major transport pathway for these compounds through the mitochondrial outer membrane. Numerous experimental studies have promoted the idea that VDAC acts as a regulator of essential mitochondrial functions. In this study, using a combination of molecular dynamics simulations, free-energy calculations, and electrophysiological measurements, we investigated the transport of ions through VDAC, with a focus on phosphate ions and metabolites. We showed that selectivity of VDAC towards small anions including monovalent phosphates arises from short-lived interactions with positively charged residues scattered throughout the pore. In dramatic contrast, permeation of divalent phosphate ions and phosphate metabolites (AMP and ATP) involves binding sites along a specific translocation pathway. This permeation mechanism offers an explanation for the decrease in VDAC conductance measured in the presence of ATP or AMP at physiological salt concentration. The binding sites occur at similar locations for the divalent phosphate ions, AMP and ATP, and contain identical basic residues. ATP features a marked affinity for a central region of the pore lined by two lysines and one arginine of the N-terminal helix. This cluster of residues together with a few other basic amino acids forms a "charged brush" which facilitates the passage of the anionic metabolites through the pore. All of this reveals that VDAC controls the transport of the inorganic phosphates and phosphate metabolites studied here through two different mechanisms.

Ionic conductivity, structural deformation, and programmable anisotropy of DNA origami in electric field

The DNA origami technique can enable functionalization of inorganic structures for single-molecule electric current recordings. Experiments have shown that several layers of DNA molecules, a DNA origami plate, placed on top of a solid-state nanopore is permeable to ions. Here, we report a comprehensive characterization of the ionic conductivity of DNA origami plates by means of all-atom molecular dynamics (MD) simulations and nanocapillary electric current recordings. Using the MD method, we characterize the ionic conductivity of several origami constructs, revealing the local distribution of ions, the distribution of the electrostatic potential and contribution of different molecular species to the current. The simulations determine the dependence of the ionic conductivity on the applied voltage, the number of DNA layers, the nucleotide content and the lattice type of the plates. We demonstrate that increasing the concentration of Mg(2+) ions makes the origami plates more compact, reducing their conductivity. The conductance of a DNA origami plate on top of a solid-state nanopore is determined by the two competing effects: bending of the DNA origami plate that reduces the current and separation of the DNA origami layers that increases the current. The latter is produced by the electro-osmotic flow and is reversible at the time scale of a hundred nanoseconds. The conductance of a DNA origami object is found to depend on its orientation, reaching maximum when the electric field aligns with the direction of the DNA helices. Our work demonstrates feasibility of programming the electrical properties of a self-assembled nanoscale object using DNA.

Modeling thermophoretic effects in solid-state nanopores

Local modulation of temperature has emerged as a new mechanism for regulation of molecular transport through nanopores. Predicting the effect of such modulations on nanopore transport requires simulation protocols capable of reproducing non-uniform temperature gradients observed in experiment. Conventional molecular dynamics (MD) method typically employs a single thermostat for maintaining a uniform distribution of temperature in the entire simulation domain, and, therefore, can not model local temperature variations. In this article, we describe a set of simulation protocols that enable modeling of nanopore systems featuring non-uniform distributions of temperature. First, we describe a method to impose a temperature gradient in all-atom MD simulations based on a boundary-driven non-equilibrium MD protocol. Then, we use this method to study the effect of temperature gradient on the distribution of ions in bulk solution (the thermophoretic effect). We show that DNA nucleotides exhibit differential response to the same temperature gradient. Next, we describe a method to directly compute the effective force of a thermal gradient on a prototypical biomolecule-a fragment of double-stranded DNA. Following that, we demonstrate an all-atom MD protocol for modeling thermophoretic effects in solid-state nanopores. We show that local heating of a nanopore volume can be used to regulate the nanopore ionic current. Finally, we show how continuum calculations can be coupled to a coarse-grained model of DNA to study the effect of local temperature modulation on electrophoretic motion of DNA through plasmonic nanopores. The computational methods described in this article are expected to find applications in rational design of temperature-responsive nanopore systems.

Designing a hydrophobic barrier within biomimetic nanopores

Nanopores in membranes have a range of potential applications. Biomimetic design of nanopores aims to mimic key functions of biological pores within a stable template structure. Molecular dynamics simulations have been used to test whether a simple β-barrel protein nanopore can be modified to incorporate a hydrophobic barrier to permeation. Simulations have been used to evaluate functional properties of such nanopores, using water flux as a proxy for ionic conductance. The behavior of these model pores has been characterized as a function of pore size and of the hydrophobicity of the amino acid side chains lining the narrow central constriction of the pore. Potential of mean force calculations have been used to calculate free energy landscapes for water and for ion permeation in selected models. These studies demonstrate that a hydrophobic barrier can indeed be designed into a β-barrel protein nanopore, and that the height of the barrier can be adjusted by modifying the number of consecutive rings of hydrophobic side chains. A hydrophobic barrier prevents both water and ion permeation even though the pore is sterically unoccluded. These results both provide insights into the nature of hydrophobic gating in biological pores and channels, and furthermore demonstrate that simple design features may be computationally transplanted into β-barrel membrane proteins to generate functionally complex nanopores.

Structure and response to flow of the glycocalyx layer

The glycocalyx is a sugar-rich layer located at the luminal part of the endothelial cells. It is involved in key metabolic processes and its malfunction is related to several diseases. To understand the function of the glycocalyx, a molecular level characterization is necessary. In this article, we present large-scale molecular-dynamics simulations that provide a comprehensive description of the structure and dynamics of the glycocalyx. We introduce the most detailed, to-date, all-atom glycocalyx model, composed of lipid bilayer, proteoglycan dimers, and heparan sulfate chains with realistic sequences. Our results reveal the folding of proteoglycan ectodomain and the extended conformation of heparan sulfate chains. Furthermore, we study the glycocalyx response under shear flow and its role as a flypaper for binding fibroblast growth factors (FGFs), which are involved in diverse functions related to cellular differentiation, including angiogenesis, morphogenesis, and wound healing. The simulations show that the glycocalyx increases the effective concentration of FGFs, leading to FGF oligomerization, and acts as a lever to transfer mechanical stimulus into the cytoplasmic side of endothelial cells.

xMDFF: molecular dynamics flexible fitting of low-resolution X-ray structures

X-ray crystallography remains the most dominant method for solving atomic structures. However, for relatively large systems, the availability of only medium-to-low-resolution diffraction data often limits the determination of all-atom details. A new molecular dynamics flexible fitting (MDFF)-based approach, xMDFF, for determining structures from such low-resolution crystallographic data is reported. xMDFF employs a real-space refinement scheme that flexibly fits atomic models into an iteratively updating electron-density map. It addresses significant large-scale deformations of the initial model to fit the low-resolution density, as tested with synthetic low-resolution maps of D-ribose-binding protein. xMDFF has been successfully applied to re-refine six low-resolution protein structures of varying sizes that had already been submitted to the Protein Data Bank. Finally, via systematic refinement of a series of data from 3.6 to 7 Å resolution, xMDFF refinements together with electrophysiology experiments were used to validate the first all-atom structure of the voltage-sensing protein Ci-VSP.

GPU-accelerated analysis and visualization of large structures solved by molecular dynamics flexible fitting

Hybrid structure fitting methods combine data from cryo-electron microscopy and X-ray crystallography with molecular dynamics simulations for the determination of all-atom structures of large biomolecular complexes. Evaluating the quality-of-fit obtained from hybrid fitting is computationally demanding, particularly in the context of a multiplicity of structural conformations that must be evaluated. Existing tools for quality-of-fit analysis and visualization have previously targeted small structures and are too slow to be used interactively for large biomolecular complexes of particular interest today such as viruses or for long molecular dynamics trajectories as they arise in protein folding. We present new data-parallel and GPU-accelerated algorithms for rapid interactive computation of quality-of-fit metrics linking all-atom structures and molecular dynamics trajectories to experimentally-determined density maps obtained from cryo-electron microscopy or X-ray crystallography. We evaluate the performance and accuracy of the new quality-of-fit analysis algorithms vis-à-vis existing tools, examine algorithm performance on GPU-accelerated desktop workstations and supercomputers, and describe new visualization techniques for results of hybrid structure fitting methods.

Comparative analysis of nucleotide translocation through protein nanopores using steered molecular dynamics and an adaptive biasing force

The translocation of nucleotide molecules across biological and synthetic nanopores has attracted attention as a next generation technique for sequencing DNA. Computer simulations have the ability to provide atomistic-level insight into important states and processes, delivering a means to develop a fundamental understanding of the translocation event, for example, by extracting the free energy of the process. Even with current supercomputing facilities, the simulation of many-atom systems in fine detail is limited to shorter timescales than the real events they attempt to recreate. This imposes the need for enhanced simulation techniques that expand the scope of investigation in a given timeframe. There are numerous free energy calculation and translocation methodologies available, and it is by no means clear which method is best applied to a particular problem. This article explores the use of two popular free energy calculation methodologies in a nucleotide-nanopore translocation system, using the α-hemolysin nanopore. The first uses constant velocity-steered molecular dynamics (cv-SMD) in conjunction with Jarzynski's equality. The second applies an adaptive biasing force (ABF), which has not previously been applied to the nucleotide-nanpore system. The purpose of this study is to provide a comprehensive comparison of these methodologies, allowing for a detailed comparative assessment of the scientific merits, the computational cost, and the statistical quality of the data obtained from each technique. We find that the ABF method produces results that are closer to experimental measurements than those from cv-SMD, whereas the net errors are smaller for the same computational cost.

A Stabilized Finite Element Method for Modified Poisson-Nernst-Planck Equations to Determine Ion Flow Through a Nanopore

The conventional Poisson-Nernst-Planck equations do not account for the finite size of ions explicitly. This leads to solutions featuring unrealistically high ionic concentrations in the regions subject to external potentials, in particular, near highly charged surfaces. A modified form of the Poisson-Nernst-Planck equations accounts for steric effects and results in solutions with finite ion concentrations. Here, we evaluate numerical methods for solving the modified Poisson-Nernst-Planck equations by modeling electric field-driven transport of ions through a nanopore. We describe a novel, robust finite element solver that combines the applications of the Newton's method to the nonlinear Galerkin form of the equations, augmented with stabilization terms to appropriately handle the drift-diffusion processes. To make direct comparison with particle-based simulations possible, our method is specifically designed to produce solutions under periodic boundary conditions and to conserve the number of ions in the solution domain. We test our finite element solver on a set of challenging numerical experiments that include calculations of the ion distribution in a volume confined between two charged plates, calculations of the ionic current though a nanopore subject to an external electric field, and modeling the effect of a DNA molecule on the ion concentration and nanopore current.

Stretching and controlled motion of single-stranded DNA in locally heated solid-state nanopores

Practical applications of solid-state nanopores for DNA detection and sequencing require the electrophoretic motion of DNA through the nanopores to be precisely controlled. Controlling the motion of single-stranded DNA presents a particular challenge, in part because of the multitude of conformations that a DNA strand can adopt in a nanopore. Through continuum, coarse-grained and atomistic modeling, we demonstrate that local heating of the nanopore volume can be used to alter the electrophoretic mobility and conformation of single-stranded DNA. In the nanopore systems considered, the temperature near the nanopore is modulated via a nanometer-size heater element that can be radiatively switched on and off. The local enhancement of temperature produces considerable stretching of the DNA fragment confined within the nanopore. Such stretching is reversible, so that the conformation of DNA can be toggled between compact (local heating is off) and extended (local heating is on) states. The effective thermophoretic force acting on single-stranded DNA in the vicinity of the nanopore is found to be sufficiently large (4-8 pN) to affect such changes in the DNA conformation. The local heating of the nanopore volume is observed to promote single-file translocation of DNA strands at transmembrane biases as low as 10 mV, which opens new avenues for using solid-state nanopores for detection and sequencing of DNA.

Translocation dynamics of freely jointed Lennard-Jones chains into adsorbing pores

Polymer translocation into adsorbing nanopores is studied by using the Fokker-Planck equation of chain diffusion along the energy landscape calculated with Monte Carlo simulations using the incremental gauge cell method. The free energy profile of a translocating chain was found by combining two independent sub-chains, one free but tethered to a hard wall, and the other tethered inside an adsorbing pore. Translocation dynamics were revealed by application of the Fokker-Planck equation for normal diffusion. Adsorption of polymer chains into nanopores involves a competition of attractive adsorption and repulsive steric hindrance contributions to the free energy. Translocation times fell into two regimes depending on the strength of the adsorbing pore. In addition, we found a non-monotonic dependence of translocation times with increasing adsorption strength, with sharp peak associated with local free energy minima along the translocation coordinate.

Single-stranded DNA within nanopores: conformational dynamics and implications for sequencing; a molecular dynamics simulation study

Engineered protein nanopores, such as those based on α-hemolysin from Staphylococcus aureus have shown great promise as components of next-generation DNA sequencing devices. However, before such protein nanopores can be used to their full potential, the conformational dynamics and translocation pathway of the DNA within them must be characterized at the individual molecule level. Here, we employ atomistic molecular dynamics simulations of single-stranded DNA movement through a model α-hemolysin pore under an applied electric field. The simulations enable characterization of the conformations adopted by single-stranded DNA, and allow exploration of how the conformations may impact on translocation within the wild-type model pore and a number of mutants. Our results show that specific interactions between the protein nanopore and the DNA can have a significant impact on the DNA conformation often leading to localized coiling, which in turn, can alter the order in which the DNA bases exit the nanopore. Thus, our simulations show that strategies to control the conformation of DNA within a protein nanopore would be a distinct advantage for the purposes of DNA sequencing.

Membrane sculpting by F-BAR domains studied by molecular dynamics simulations

Interplay between cellular membranes and their peripheral proteins drives many processes in eukaryotic cells. Proteins of the Bin/Amphiphysin/Rvs (BAR) domain family, in particular, play a role in cellular morphogenesis, for example curving planar membranes into tubular membranes. However, it is still unclear how F-BAR domain proteins act on membranes. Electron microscopy revealed that, in vitro, F-BAR proteins form regular lattices on cylindrically deformed membrane surfaces. Using all-atom and coarse-grained (CG) molecular dynamics simulations, we show that such lattices, indeed, induce tubes of observed radii. A 250 ns all-atom simulation reveals that F-BAR domain curves membranes via the so-called scaffolding mechanism. Plasticity of the F-BAR domain permits conformational change in response to membrane interaction, via partial unwinding of the domains 3-helix bundle structure. A CG simulation covering more than 350 µs provides a dynamic picture of membrane tubulation by lattices of F-BAR domains. A series of CG simulations identified the optimal lattice type for membrane sculpting, which matches closely the lattices seen through cryo-electron microscopy.

To generate organelles, eukaryotic cells sculpt their membranes into compartments, often employing proteins as chaperones, for example, F-BAR domains. The latter induce formation of tubular and vesicular membranes. Functional and structural studies suggest that F-BAR domains sculpt membranes through electrostatic interactions, driving the membrane to match the concave surface of the protein's banana-like shape. Cryo-electron microscopy (cryo-EM) studies provide an average static picture of how F-BAR domains form lattices on the surface of membranes to induce tube formation. Complementing the cryo-EM images, molecular dynamics simulations reported here offer a detailed, dynamic picture of membrane tubulation by a lattice of F-BAR domains and identified lattice types optimally attuned to producing high membrane curvature. The simulations reproduced also a process lasting 350 µs in which lattices of F-BAR domains form a complete tube out of an initially flat membrane. The molecular dynamics study offers, thereby, both a large-scale picture of membrane sculpting by F-BAR domain lattices as well as atomic-level dynamic information about the involvement of the individual F-BAR domain and its interactions with partner F-BAR domains and membrane in the sculpting process.

Polypeptide Translocation Through the Mitochondrial TOM Channel: Temperature-Dependent Rates at the Single-Molecule Level

The TOM protein complex facilitates the transfer of nearly all mitochondrial preproteins across outer mitochondrial membranes. Here we characterized the effect of temperature on facilitated translocation of a mitochondrial presequence peptide pF1β. Ion current fluctuations analysis through single TOM channels revealed thermodynamic and kinetic parameters of substrate binding and allowed determining the energy profile of peptide translocation. The activation energy for the on-rate and off-rate of the presequence peptide into the TOM complex was symmetric with respect to the electric field and estimated to be about 15 and 22 kT per peptide. These values are above that expected for free diffusion of ions in water (6 kT) and reflect the stronger interaction in the channel. Both values are in the range for typical enzyme kinetics and suggest one process without involving large conformational changes within the channel protein.

Molecular dynamics simulations in drug design

This minireview focuses on recent developments in the application of molecular dynamics to drug design. Recent applications of endpoint free-energy computational methods such as molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and generalized Born surface area (MM-GBSA) and linear response methods are described. Recent progress in steered molecular dynamics applied to drug design is reviewed.

Molecular dynamics study of MspA arginine mutants predicts slow DNA translocations and ion current blockades indicative of DNA sequence

The protein nanopore Mycobacteria smegmatis porin A (MspA), can be used to sense individual nucleotides within DNA, potentially enabling a technique known as nanopore sequencing. In this technique, single-stranded DNA electrophoretically moves through the nanopore and results in an ionic current that is nucleotide-specific. However, with a high transport velocity of the DNA within the nanopore, the ionic current cannot be used to distinguish signals within noise. Through extensive (~100 μs in total) all-atom molecular dynamics simulations, we examine the effect of positively charged residues on DNA translocation rate and the ionic current blockades in MspA. Simulation of several arginine mutations show a ~10-30 fold reduction of DNA translocation speed without eliminating the nucleotide induced current blockages. Comparison of our results with similar engineering efforts on a different nanopore (α-hemolysin) reveals a nontrivial effect of nanopore geometry on the ionic current blockades in mutant nanopores.

GPU/CPU Algorithm for Generalized Born/Solvent-Accessible Surface Area Implicit Solvent Calculations

Molecular dynamics methodologies comprise a vital research tool for structural biology. Molecular dynamics has benefited from technological advances in computing, such as multi-core CPUs and graphics processing units (GPUs), but harnessing the full power of hybrid GPU/CPU computers remains difficult. The generalized Born/solvent-accessible surface area implicit solvent model (GB/SA) stands to benefit from hybrid GPU/CPU computers, employing the GPU for the GB calculation and the CPU for the SA calculation. Here, we explore the computational challenges facing GB/SA calculations on hybrid GPU/CPU computers and demonstrate how NAMD, a parallel molecular dynamics program, is able to efficiently utilize GPUs and CPUs simultaneously for fast GB/SA simulations. The hybrid computation principles demonstrated here are generally applicable to parallel applications employing hybrid GPU/CPU calculations.

Slowing and controlling the translocation of DNA in a solid-state nanopore

DNA sequencing methods based on nanopores could potentially represent a low-cost and high-throughput pathway to practical genomics, by replacing current sequencing methods based on synthesis that are limited in speed and cost. The success of nanopore sequencing techniques requires the solution to two fundamental problems: (1) sensing each nucleotide of a DNA strand, in sequence, as it passes through a nanopore; (2) delivering each nucleotide in a DNA strand, in turn, to a sensing site within the nanopore in a controlled manner. It has been demonstrated that a DNA nucleotide can be sensed using electric signals, such as ionic current changes caused by nucleotide blockage at a constriction region in a protein pore or a tunneling current through the nucleotide-bridged gap of two nanoelectrodes built near a solid-state nanopore. However, it is not yet clear how each nucleotide in a DNA strand can be delivered in turn to a sensing site and held there for a sufficient time to ensure high fidelity sensing. This latter problem has been addressed by modifying macroscopic properties, such as a solvent viscosity, ion concentration or temperature. Also, the DNA transistor, a solid state nanopore dressed with a series of metal-dielectric layers has been proposed as a solution. Molecular dynamics simulations provide the means to study and to understand DNA transport in nanopores microscopically. In this article, we review computational studies on how to slow down and control the DNA translocation through a solid-state nanopore.

Predicting the DNA sequence dependence of nanopore ion current using atomic-resolution Brownian dynamics

It has become possible to distinguish DNA molecules of different nucleotide sequences by measuring ion current passing through a narrow pore containing DNA. To assist experimentalists in interpreting the results of such measurements and to improve the DNA sequence detection method, we have developed a computational approach that has both the atomic-scale accuracy and the computational efficiency required to predict DNA sequence-specific differences in the nanopore ion current. In our Brownian dynamics method, the interaction between the ions and DNA is described by three-dimensional potential of mean force maps determined to a 0.03 nm resolution from all-atom molecular dynamics simulations. While this atomic-resolution Brownian dynamics method produces results with orders of magnitude less computational effort than all-atom molecular dynamics requires, we show here that the ion distributions and ion currents predicted by the two methods agree. Finally, using our Brownian dynamics method, we find that a small change in the sequence of DNA within a pore can cause a large change in the ion current, and validate this result with all-atom molecular dynamics.

Parallel Generalized Born Implicit Solvent Calculations with NAMD

Accurate electrostatic descriptions of aqueous solvent are critical for simulation studies of bio-molecules, but the computational cost of explicit treatment of solvent is very high. A computationally more feasible alternative is a generalized Born implicit solvent description which models polar solvent as a dielectric continuum. Unfortunately, the attainable simulation speedup does not transfer to the massive parallel computers often employed for simulation of large structures. Longer cutoff distances, spatially heterogenous distribution of atoms and the necessary three-fold iteration over atom-pairs in each timestep combine to challenge efficient parallel performance of generalized Born implicit solvent algorithms. Here we report how NAMD, a parallel molecular dynamics program, meets the challenge through a unique parallelization strategy. NAMD now permits efficient simulation of large systems whose slow conformational motions benefit most from implicit solvent descriptions due to the inherent low viscosity. NAMD's implicit solvent performance is benchmarked and then illustrated in simulating the ratcheting Escherichia coli ribosome involving ~250,000 atoms.

Characterizing and controlling the motion of ssDNA in a solid-state nanopore

Sequencing DNA in a synthetic solid-state nanopore is potentially a low-cost and high-throughput method. Essential to the nanopore-based DNA sequencing method is the ability to control the motion of a single-stranded DNA (ssDNA) molecule at single-base resolution. Experimental studies showed that the average translocation speed of DNA driven by a biasing electric field can be affected by ionic concentration, solvent viscosity, or temperature. Even though it is possible to slow down the average translocation speed, instantaneous motion of DNA is too diffusive to allow each DNA base to stay in front of a sensor site for its measurement. Using extensive all-atom molecular dynamics simulations, we study the diffusion constant, friction coefficient, electrophoretic mobility, and effective charge of ssDNA in a solid-state nanopore. Simulation results show that the spatial fluctuation of ssDNA in 1 ns is comparable to the spacing between neighboring nucleotides in ssDNA, which makes the sensing of a DNA base very difficult. We demonstrate that the recently proposed DNA transistor could potentially solve this problem by electrically trapping ssDNA inside the DNA transistor and ratcheting ssDNA base-by-base in a biasing electric field. When increasing the biasing electric field, we observed that the translocation of ssDNA changes from ratcheting to steady-sliding. The simulated translocation of ssDNA in the DNA transistor was theoretically characterized using Fokker-Planck analysis.

Control of Nanoscale Environment to Improve Stability of Immobilized Proteins on Diamond Surfaces

Immunoassays for detection of bacterial pathogens rely on the selectivity and stability of bio-recognition elements such as antibodies tethered to sensor surfaces. The search for novel surfaces that improve the stability of biomolecules and assay performance has been pursued for a long time. However, the anticipated improvements in stability have not been realized in practice under physiological conditions because the surface functionalization layers on commonly used substrates, silica and gold, are themselves unstable on time scales of days. In this paper, we show that covalent linking of antibodies to diamond surfaces leads to substantial improvements in biological activity of proteins as measured by the ability to selectively capture cells of the pathogenic bacterium Escherichia coli O157:H7 even after exposure to buffer solutions at 37 °C for extended periods of time, approaching 2 weeks. Our results from ELISA, XPS, fluorescence microscopy, and MD simulations suggest that by using highly stable surface chemistry and controlling the nanoscale organization of the antibodies on the surface, it is possible to achieve significant improvements in biological activity and stability. Our findings can be easily extended to functionalization of micro and nanodimensional sensors and structures of biomedical diagnostic and therapeutic interest.

Formation of salt bridges mediates internal dimerization of myosin VI medial tail domain

The unconventional motor protein, myosin VI, is known to dimerize upon cargo binding to its C-terminal end. It has been shown that one of its tail domains, called the medial tail domain, is a dimerization region. The domain contains an unusual pattern of alternating charged residues and a few hydrophobic residues. To reveal the unknown dimerization mechanism of the medial tail domain, we employed molecular dynamics and single-molecule experimental techniques. Both techniques suggest that the formation of electrostatic-based interhelical salt bridges between oppositely charged residues is a key dimerization factor. For the dimerization to occur, the two identical helices within the dimer do not bind in a symmetric fashion, but rather with an offset of about one helical repeat. Calculations of the dimer-dissociation energy find the contribution of hydrophobic residues to the dimerization process to be minor; they also find that the asymmetric homodimer state is energetically favorable over a state of separate helices.

Modeling Pressure-Driven Transport of Proteins through a Nanochannel

Reducing the size of a nanofluidic channel not only creates new opportunities for high-precision manipulation of biological macromolecules, but also makes the performance of the entire nanofluidic system more susceptible to undesirable interactions between the transported biomolecules and the walls of the channel. In this manuscript, we report molecular dynamics simulations of a pressure-driven flow through a silica nanochannel that characterized, with atomic resolution, adsorption of a model protein to its surface. Although the simulated adsorption of the proteins was found to be nonspecific, it had a dramatic effect on the rate of the protein transport. To determine the relative strength of the protein-silica interactions in different adsorbed states, we simulated flow-induced desorption of the proteins from the silica surface. Our analysis of the protein conformations in the adsorbed states did not reveal any simple dependence of the adsorption strength on the size and composition of the protein-silica contact, suggesting that the heterogeneity of the silica surface may be a important factor.

Modeling nanopores for sequencing DNA

Using nanopores to sequence DNA rapidly and at a low cost has the potential to radically transform the field of genomic research. However, despite all the exciting developments in the field, sequencing DNA using a nanopore has yet to be demonstrated. Among the many problems that hinder development of the nanopore sequencing methods is the inability of current experimental techniques to visualize DNA conformations in a nanopore and directly relate the microscopic state of the system to the measured signal. We have recently shown that such tasks could be accomplished through computation. This chapter provides step-by-step instructions of how to build atomic scale models of biological and solid-state nanopore systems, use the molecular dynamics method to simulate the electric field-driven transport of ions and DNA through the nanopores, and analyze the results of such computational experiments.

Membrane protein dynamics from femtoseconds to seconds

Membrane proteins play a key role in energy conversion, transport, signal recognition, transduction, and other fundamental biological processes. Despite considerable progress in experimental techniques, the determination of structure and dynamics of membrane proteins still represents a great challenge. Computer simulation methods are becoming an increasingly important tool not only in the interpretation of experiments but also in the prediction of membrane protein dynamics. In the present review, we give a brief introduction to molecular modeling techniques currently used to explore protein dynamics on time scales ranging from femtoseconds to microseconds. We then describe a few recent example applications of these techniques to membrane proteins. In conclusion, we also discuss some of the newest developments in simulation methodology that have the potential to further extend the time scale accessible to explore (membrane) protein dynamics.

Mechanical properties of a complete microtubule revealed through molecular dynamics simulation

Microtubules (MTs) are the largest type of cellular filament, essential in processes ranging from mitosis and meiosis to flagellar motility. Many of the processes depend critically on the mechanical properties of the MT, but the elastic moduli, notably the Young's modulus, are not directly revealed in experiment, which instead measures either flexural rigidity or response to radial deformation. Molecular dynamics (MD) is a method that allows the mechanical properties of single biomolecules to be investigated through computation. Typically, MD requires an atomic resolution structure of the molecule, which is unavailable for many systems, including MTs. By combining structural information from cryo-electron microscopy and electron crystallography, we have constructed an all-atom model of a complete MT and used MD to determine its mechanical properties. The simulations revealed nonlinear axial stress-strain behavior featuring a pronounced softening under extension, a possible plastic deformation transition under radial compression, and a distinct asymmetry in response to the two senses of twist. This work demonstrates the possibility of combining different levels of structural information to produce all-atom models suitable for quantitative MD simulations, which extends the range of systems amenable to the MD method and should enable exciting advances in our microscopic knowledge of biology.

Translocation of α-helix chains through a nanopore

The translocation of α-helix chains through a nanopore is studied through Langevin dynamics simulations. The α-helix chains exhibit several different characteristics about their average translocation times and the α-helix structures when they transport through the nanopores under the driving forces. First, the relationship between average translocation times τ and the chain length N satisfies the scaling law, τ∼N(α), and the scaling exponent α depends on the driving force f for the small forces while it is close to the Flory exponent (ν) in the other force regions. For the chains with given chain lengths, it is observed that the dependence of the average translocation times can be expressed as τ∼f(-1/2) for the small forces while can be described as τ∼f in the large force regions. Second, for the large driving force, the average number of α-helix structures N(h) decreases first and then increases in the translocation process. The average waiting time of each bead, especially of the first bead, is also dependent on the driving forces. Furthermore, an elasticity spring model is presented to reasonably explain the change of the α-helix number during the translocation and its elasticity can be locally damaged by the large driving forces. Our results demonstrate the unique behaviors of α-helix chains transporting through the pores, which can enrich our insights into and knowledge on biopolymers transporting through membranes.

Slowing the translocation of double-stranded DNA using a nanopore smaller than the double helix

It is now possible to slow and trap a single molecule of double-stranded DNA (dsDNA), by stretching it using a nanopore, smaller in diameter than the double helix, in a solid-state membrane. By applying an electric force larger than the threshold for stretching, dsDNA can be impelled through the pore. Once a current blockade associated with a translocating molecule is detected, the electric field in the pore is switched in an interval less than the translocation time to a value below the threshold for stretching. According to molecular dynamics (MD) simulations, this leaves the dsDNA stretched in the pore constriction with the base-pairs tilted, while the B-form canonical structure is preserved outside the pore. In this configuration, the translocation velocity is substantially reduced from 1 bp/10 ns to approximately 1 bp/2 ms in the extreme, potentially facilitating high fidelity reads for sequencing, precise sorting, and high resolution (force) spectroscopy.

DNA translocation and unzipping through a nanopore: some geometrical effects

This article explores the role of some geometrical factors on the electrophoretically driven translocations of macromolecules through nanopores. In the case of asymmetric pores, we show how the entry requirements and the direction of translocation can modify the information content of the blocked ionic current as well as the transduction of the electrophoretic drive into a mechanical force. To address these effects we studied the translocation of single-stranded DNA through an asymmetric alpha-hemolysin pore. Depending on the direction of the translocation, we measure the capacity of the pore to discriminate between both DNA orientations. By unzipping DNA hairpins from both sides of the pores we show that the presence of single strand or double strand in the pore can be discriminated based on ionic current levels. We also show that the transduction of the electrophoretic drive into a denaturing mechanical force depends on the local geometry of the pore entrance. Eventually we discuss the application of this work to the measurement of energy barriers for DNA unzipping as well as for protein binding and unfolding.

Computational microscopy of the role of protonable surface residues in nanoprecipitation oscillations

A novel phenomenon has recently been reported in polymeric nanopores. This phenomenon, so-called nanoprecipitation, is characterized by the transient formation of precipitates in the nanopore lumen, producing a sequence of low and high conductance states in the ionic current through the pore. By means of all-atom molecular dynamics simulations, we studied nanoprecipitation for polyethylene terephthalate nanopore immersed in electrolytic solution containing calcium phosphate, covering a total simulation time of 1.24 micros. Our results suggest that protonable surface residues at the nanopore surface, namely carboxyl groups, trigger the formation of precipitates that strongly adhere to the surface, blocking the pore and producing the low conductance state. On the basis of the simulations, we propose a mechanism for the formation of the high conductance state; the mechanism involves detachment of the precipitate from the surface due to reprotonation of carboxyl groups and subsequent translocation of the precipitate out of the pore.

Energy Transfer Dynamics in an RC-LH1-PufX Tubular Photosynthetic Membrane

Light absorption and the subsequent transfer of excitation energy are the first two steps of the photosynthetic process, carried out by protein-bound pigments, mainly bacteriochlorophylls (BChls), in photosynthetic bacteria. BChls are anchored in light-harvesting (LH) complexes, such as light-harvesting complex I (LH1), which directly associates with the reaction center (RC), forming the RC-LH1 core complex. In Rhodobacter sphaeroides, RC-LH1 core complexes contain an additional protein, PufX, and assemble into dimeric RC-LH1-PufX core complexes. In the absence of light-harvesting complexes II, the former complexes can aggregate into a helically ordered tubular photosynthetic membrane. We examined the excitation transfer dynamics in a single RC-LH1-PufX core complex dimer using the hierarchical equations of motion for dissipative quantum dynamics that accurately, yet computationally costly, treat the coupling between BChls and their protein environment. A widely employed description, generalized Förster theory, was also used to calculate the transfer rates of the same excitonic system in order to verify the accuracy of this computationally cheap method. Additionally, in light of the structural uncertainties in the Rhodobacter sphaeroides RC-LH1-PufX core complex, geometrical alterations were introduced in the BChl organization. It is shown that the energy transfer dynamics is not affected by the considered changes in the BChl organization, and that generalized Förster theory provides accurate transfer rates. An all-atom model for a tubular photosynthetic membrane is then constructed on the basis of electron microscopy data, and the overall energy transfer properties of this membrane are computed.

Base-by-base ratcheting of single stranded DNA through a solid-state nanopore

We investigate the base-by-base translocation dynamics of single-stranded DNA (ssDNA) confined in a solid-state nanopore dressed with an electrostatic trap, using all-atom molecular dynamics (MD) simulation. We observe on the simulation time scale of tens of nanoseconds that ssDNA can be driven through the nanopore in a ratchetlike fashion, with a step size equal to the spacing between neighboring phosphate groups in the ssDNA backbone. A 1D-Langevin-like model is derived from atomistic dynamics which can quantitatively describe simulation results and can be used to study dynamics on longer time scales. The controlled ratcheting motion of DNA could potentially enhance the signal-to-noise ratio for nanoelectronic DNA sensing technologies.

Molecular control of ionic conduction in polymer nanopores

Polymeric nanopores show unique transport properties and have attracted a great deal of scientific interest as a test system to study ionic and molecular transport at the nanoscale. By means of all-atom molecular dynamics, we simulated the ion dynamics inside polymeric polyethylene terephthalate nanopores. For this purpose, we established a protocol to assemble atomic models of polymeric material into which we sculpted a nanopore model with the key features of experimental devices, namely a conical geometry and a negative surface charge density. Molecular dynamics simulations of ion currents through the pore show that the protonation state of the carboxyl group of exposed residues have a considerable effect on ion selectivity, by affecting ionic densities and electrostatic potentials inside the nanopores. The role of high concentrations of Ca2+ ions was investigated in detail.

Self-assembly of photosynthetic membranes

Bacterial photosynthetic membranes, also known as chromatophores, are tightly packed with integral membrane proteins that work together to carry out photosynthesis. Chromatophores display a wide range of cellular morphologies; spherical, tubular, and lamellar chromatophores have all been observed in different bacterial species, or with different protein constituents. Through recent computational modeling and simulation, it has been demonstrated that the light-harvesting complexes abundant in chromatophores induce local membrane curvature via multiple mechanisms. These protein complexes assemble to generate a global curvature and sculpt the chromatophores into various cellular-scale architectures.

Continuum electrostatic calculations of the pKa of ionizable residues in an ion channel: dynamic vs. static input structure

We have computed the pK(a)'s of the ionizable residues of a protein ion channel, the Staphylococcus aureus toxin alpha-hemolysin, by using two types of input structures, namely the crystal structure of the heptameric alpha-hemolysin and a set of over four hundred snapshots from a 4.38 ns Molecular Dynamics simulation of the protein inserted in a phospholipid planar bilayer. The comparison of the dynamic picture provided by the Molecular Simulation with the static one based on the X-ray crystal structure of the protein embedded in a lipid membrane allows analyzing the influence of the fluctuations in the protein structure on its ionization properties. We find that the use of the dynamic structure provides interesting information about the sensitivity of the computed pK(a) of a given residue to small changes in the local structure. The calculated pK(a) are consistent with previous indirect estimations obtained from single-channel conductance and selectivity measurements.

Deciphering ionic current signatures of DNA transport through a nanopore

Within just a decade from the pioneering work demonstrating the utility of nanopores for molecular sensing, nanopores have emerged as versatile systems for single-molecule manipulation and analysis. In a typical setup, a gradient of the electrostatic potential captures charged solutes from the solution and forces them to move through a single nanopore, across an otherwise impermeable membrane. The ionic current blockades resulting from the presence of a solute in a nanopore can reveal the type of the solute, for example, the nucleotide makeup of a DNA strand. Despite great success, the microscopic mechanisms underlying the functionality of such stochastic sensors remain largely unknown, as it is not currently possible to characterize the microscopic conformations of single biomolecules directly in a nanopore and thereby unequivocally establish the causal relationship between the observables and the microscopic events. Such a relationship can be determined using molecular dynamics-a computational method that can accurately predict the time evolution of a molecular system starting from a given microscopic state. This article describes recent applications of this method to the process of DNA transport through biological and synthetic nanopores.

Regulation of the protein-conducting channel by a bound ribosome

During protein synthesis, it is often necessary for the ribosome to form a complex with a membrane-bound channel, the SecY/Sec61 complex, in order to translocate nascent proteins across a cellular membrane. Structural data on the ribosome-channel complex are currently limited to low-resolution cryo-electron microscopy maps, including one showing a bacterial ribosome bound to a monomeric SecY complex. Using that map along with available atomic-level models of the ribosome and SecY, we have determined, through molecular dynamics flexible fitting (MDFF), an atomic-resolution model of the ribosome-channel complex. We characterized computationally the sites of ribosome-SecY interaction within the complex and determined the effect of ribosome binding on the SecY channel. We also constructed a model of a ribosome in complex with a SecY dimer by adding a second copy of SecY to the MDFF-derived model. The study involved 2.7-million-atom simulations over altogether nearly 50 ns.

A statistical model for translocation of structured polypeptide chains through nanopores

The translocation process of a globular protein (ubiquitin) across a cylindrical nanopore is studied via molecular dynamics simulations. The ubiquitin is described by a native-centric model on a Calpha carbon backbone to investigate the influence of protein-like structural properties on the translocation mechanism. A thermodynamical and kinetic characterization of the process is obtained by studying the statistics of blockage times, the mobility, and the translocation probability as a function of the pulling force F acting in the pore. The transport dynamics occurs when the force exceeds a threshold Fc depending on a free-energy barrier that ubiquitin has to overcome in order to slide along the channel. Such a barrier results from competition of the unfolding energy and the entropy associated with the confinement effects of the pore. We implement appropriate umbrella sampling simulations to compute the free-energy profile as a function of the position of the ubiquitin center of mass inside of the channel (reaction coordinate). This free energy is then used to construct a phenomenological drift-diffusion model in the reaction coordinate which explains and reproduces the behavior of the observables during the translocation.

Molecular dynamics flexible fitting: a practical guide to combine cryo-electron microscopy and X-ray crystallography

Hybrid computational methods for combining structural data from different sources and resolutions are becoming an essential part of structural biology, especially as the field moves toward the study of large macromolecular assemblies. We have developed the molecular dynamics flexible fitting (MDFF) method for combining high-resolution atomic structures with cryo-electron microscopy (cryo-EM) maps, that results in atomic models representing the conformational state captured by cryo-EM. The method has been applied successfully to the ribosome, a ribonucleoprotein complex responsible for protein synthesis. MDFF involves a molecular dynamics simulation in which a guiding potential, based on the cryo-EM map, is added to the standard force field. Forces proportional to the gradient of the density map guide an atomic structure, available from X-ray crystallography, into high-density regions of a cryo-EM map. In this paper we describe the necessary steps to set up, run, and analyze MDFF simulations and the software packages that implement the corresponding functionalities.

Protein-induced membrane curvature investigated through molecular dynamics flexible fitting

In the photosynthetic purple bacterium Rhodobacter (Rba.) sphaeroides, light is absorbed by membrane-bound light-harvesting (LH) proteins LH1 and LH2. LH1 directly surrounds the reaction center (RC) and, together with PufX, forms a dimeric (RC-LH1-PufX)2 protein complex. In LH2-deficient Rba. sphaeroides mutants, RC-LH1-PufX dimers aggregate into tubular vesicles with a radius of approximately 250-550 A, making RC-LH1-PufX one of the few integral membrane proteins known to actively induce membrane curvature. Recently, a three-dimensional electron microscopy density map showed that the Rba. sphaeroides RC-LH1-PufX dimer exhibits a prominent bend at its dimerizing interface. To investigate the curvature properties of this highly bent protein, we employed molecular dynamics simulations to fit an all-atom structural model of the RC-LH1-PufX dimer within the electron microscopy density map. The simulations reveal how the dimer produces a membrane with high local curvature, even though the location of PufX cannot yet be determined uniquely. The resulting membrane curvature agrees well with the size of RC-LH1-PufX tubular vesicles, and demonstrates how the local curvature properties of the RC-LH1-PufX dimer propagate to form the observed long-range organization of the Rba. sphaeroides tubular vesicles.