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. © 2014 The Authors Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

[1]  Syma Khalid,et al.  Molecular dynamics simulations of DNA within a nanopore: arginine-phosphate tethering and a binding/sliding mechanism for translocation. , 2011, Biochemistry.

[2]  S. Howorka,et al.  Chemical tags facilitate the sensing of individual DNA strands with nanopores. , 2008, Angewandte Chemie.

[3]  M. Troll,et al.  Determination of RNA orientation during translocation through a biological nanopore. , 2006, Biophysical journal.

[4]  Mark Akeson,et al.  Single-molecule analysis of DNA-protein complexes using nanopores , 2007, Nature Methods.

[5]  C. Chipot,et al.  Overcoming free energy barriers using unconstrained molecular dynamics simulations. , 2004, The Journal of chemical physics.

[6]  D. Branton,et al.  Voltage-driven DNA translocations through a nanopore. , 2001, Physical review letters.

[7]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[8]  J Davies,et al.  Inactivation of antibiotics and the dissemination of resistance genes. , 1994, Science.

[9]  D. Branton,et al.  The potential and challenges of nanopore sequencing , 2008, Nature Biotechnology.

[10]  P. Klotman,et al.  Identification and characterization of a cell membrane nucleic acid channel. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[11]  A. Aksimentiev,et al.  Exploring transmembrane transport through α -hemolysin with grid-steered molecular dynamics , 2007 .

[12]  K. Schulten,et al.  Imaging alpha-hemolysin with molecular dynamics: ionic conductance, osmotic permeability, and the electrostatic potential map. , 2005, Biophysical journal.

[13]  E. Koonin,et al.  Horizontal gene transfer in prokaryotes: quantification and classification. , 2001, Annual review of microbiology.

[14]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[15]  Andre Marziali,et al.  Disentangling steric and electrostatic factors in nanoscale transport through confined space. , 2013, Nano letters (Print).

[16]  Christophe Chipot,et al.  Exploring Multidimensional Free Energy Landscapes Using Time-Dependent Biases on Collective Variables. , 2010, Journal of chemical theory and computation.

[17]  K. Schulten,et al.  Orientation discrimination of single-stranded DNA inside the α-hemolysin membrane channel , 2005 .

[18]  P. Coveney,et al.  Determination of Free Energy Profiles for the Translocation of Polynucleotides through α-Hemolysin Nanopores using Non-Equilibrium Molecular Dynamics Simulations. , 2009, Journal of chemical theory and computation.

[19]  Eric Darve,et al.  Calculating Free Energies Using a Scaled-Force Molecular Dynamics Algorithm , 2002 .

[20]  J. Shendure,et al.  Materials and Methods Som Text Figs. S1 and S2 Tables S1 to S4 References Accurate Multiplex Polony Sequencing of an Evolved Bacterial Genome , 2022 .

[21]  G. Crooks.  Nonequilibrium Measurements of Free Energy Differences for Microscopically Reversible Markovian Systems , 1998 .

[22]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[23]  D. Branton,et al.  Characterization of individual polynucleotide molecules using a membrane channel. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[24]  S. Howorka,et al.  Sequence-specific detection of individual DNA strands using engineered nanopores , 2001, Nature Biotechnology.

[25]  H. Bayley,et al.  Continuous base identification for single-molecule nanopore DNA sequencing. , 2009, Nature nanotechnology.

[26]  Y. Danieli.  Guide , 2005 .

[27]  Simone Marsili,et al.  Self-healing umbrella sampling: a non-equilibrium approach for quantitative free energy calculations. , 2006, The journal of physical chemistry. B.

[28]  D. Branton,et al.  Microsecond time-scale discrimination among polycytidylic acid, polyadenylic acid, and polyuridylic acid as homopolymers or as segments within single RNA molecules. , 1999, Biophysical journal.

[29]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[30]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[31]  D. Branton,et al.  Rapid nanopore discrimination between single polynucleotide molecules. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[32]  H. Bayley,et al.  Enhanced translocation of single DNA molecules through α-hemolysin nanopores by manipulation of internal charge , 2008, Proceedings of the National Academy of Sciences.

[33]  H. Bayley,et al.  Recognizing a single base in an individual DNA strand: a step toward DNA sequencing in nanopores. , 2005, Angewandte Chemie.

[34]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[35]  A. Aksimentiev,et al.  Molecular dynamics study of MspA arginine mutants predicts slow DNA translocations and ion current blockades indicative of DNA sequence. , 2012, ACS nano.

[36]  R. Blaese,et al.  Transfer of the bacterial gene for cytosine deaminase to mammalian cells confers lethal sensitivity to 5-fluorocytosine: a negative selection system. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[37]  K. Lieberman,et al.  Processive replication of single DNA molecules in a nanopore catalyzed by phi29 DNA polymerase. , 2010, Journal of the American Chemical Society.

[38]  Alexander D. MacKerell,et al.  An Improved Empirical Potential Energy Function for Molecular Simulations of Phospholipids , 2000 .

[39]  Christophe Chipot,et al.  Exploring the free-energy landscape of a short peptide using an average force. , 2005, The Journal of chemical physics.