Free-energy calculations reveal the subtle differences in the interactions of DNA bases with α-hemolysin.

Next generation DNA sequencing methods that utilize protein nanopores have the potential to revolutionize this area of biotechnology. While the technique is underpinned by simple physics, the wild-type protein pores do not have all of the desired properties for efficient and accurate DNA sequencing. Much of the research efforts have focused on protein nanopores, such as α-hemolysin from Staphylococcus aureus. However, the speed of DNA translocation has historically been an issue, hampered in part by incomplete knowledge of the energetics of translocation. Here we have utilized atomistic molecular dynamics simulations of nucleotide fragments in order to calculate the potential of mean force (PMF) through α-hemolysin. Our results reveal specific regions within the pore that play a key role in the interaction with DNA. In particular, charged residues such as D127 and K131 provide stabilizing interactions with the anionic DNA and therefore are likely to reduce the speed of translocation. These regions provide rational targets for pore optimization. Furthermore, we show that the energetic contributions to the protein-DNA interactions are a complex combination of electrostatics and short-range interactions, often mediated by water molecules.

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.

Molecular dynamics simulations of DNA within a nanopore: arginine-phosphate tethering and a binding/sliding mechanism for translocation.

Protein nanopores show great potential as low-cost detectors in DNA sequencing devices. To date, research has largely focused on the staphylococcal pore α-hemolysin (αHL). In the present study, we have developed simplified models of the wild-type αHL pore and various mutants in order to study the translocation dynamics of single-stranded DNA under the influence of an applied electric field. The model nanopores reflect the experimentally measured conductance values in planar lipid bilayers. We show that interactions between rings of cationic amino acids and DNA backbone phosphates result in metastable tethering of nucleic acid molecules within the pore, leading us to propose a "binding and sliding" mechanism for translocation. We also observe folding of DNA into nonlinear conformational intermediates during passage through the confined nanopore environment. Despite adopting nonlinear conformations, the DNA hexamer always exits the pore in the same orientation as it enters (3' to 5') in our simulations. The observations from our simulations help to rationalize experimentally determined trends in residual current and translocation efficiency for αHL and its mutants.