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

Modeling Transport Through Synthetic Nanopores

Nanopores in thin synthetic membranes have emerged as convenient tools for high-throughput single-molecule manipulation and analysis. Because of their small sizes and their ability to selectively transport solutes through otherwise impermeable membranes, nanopores have numerous potential applications in nanobiotechnology. For most applications, properties of the nanopore systems have to be characterize at the atomic level, which is currently beyond the limit of experimental methods. Molecular dynamics (MD) simulations can provide the desired information, however several technical challenges have to be met before this method can be applied to synthetic nanopore systems. Here, we highlight our recent work on modeling synthetic nanopores of the most common types. First, we describe a novel graphical tool for setting up all-atom systems incorporating inorganic materials and biomolecules. Next, we illustrate the application of the MD method for silica, silicon nitride, and polyethylene terephthalate nanopores. Following that, we describe a method for modeling synthetic surfaces using a bias potential. Future directions for tool development and nanopore modeling are briefly discussed at the end of this article.

Computer simulations and theory of protein translocation

The translocation of proteins through pores is central to many biological phenomena, such as mitochondrial protein import, protein degradation, and delivery of protein toxins to their cytosolic targets. Because proteins typically have to pass through constrictions that are too narrow to accommodate folded structures, translocation must be coupled to protein unfolding. The simplest model that accounts for such co-translocational unfolding assumes that both translocation and unfolding are accomplished by pulling on the end of the polypeptide chain mechanically. In this Account, we describe theoretical studies and computer simulations of this model and discuss how the time scales of translocation depend on the pulling force and on the protein structure. Computationally, this is a difficult problem because biologically or experimentally relevant time scales of translocation are typically orders of magnitude slower than those accessible by fully atomistic simulations. For this reason, we explore one-dimensional free energy landscapes along suitably defined translocation coordinates and discuss various approaches to their computation. We argue that the free energy landscape of translocation is often bumpy because confinement partitions the protein's configuration space into distinct basins of attraction separated by large entropic barriers. Favorable protein-pore interactions and nonnative interactions within the protein further contribute to the complexity. Computer simulations and simple scaling estimates show that forces of just 2-6 pN are often sufficient to ensure transport of unstructured polypeptides, whereas much higher forces are typically needed to translocate folded protein domains. The unfolding mechanisms found from simulations of translocation are different from those observed in the much better understood case of atomic force microscopy (AFM) pulling studies, in which proteins are unraveled by stretching them between their N- and C-termini. In contrast to AFM experiments, single-molecule experimental studies of protein translocation have just started to emerge. We describe one example of a collaborative study, in which dwell times of beta-hairpin-forming peptides inside the alpha-hemolysin pore were both measured experimentally and estimated using computer simulations. Analysis of the simulated trajectories has explained the experimental finding that more stable hairpins take, on the average, longer to traverse the pore. Despite the insight we have gained, the general relationship between the structure of proteins and their resistance to mechanically driven co-translocational unfolding remains poorly understood. Future theoretical progress likely will be made in conjunction with single-molecule experiments and will require realistic models to account for specific protein-pore interactions and for solvent effects.

Structural model and excitonic properties of the dimeric RC-LH1-PufX complex from Rhodobacter sphaeroides

The light-harvesting apparatus of the purple bacterial photosynthetic unit consists of a pool of peripheral light-harvesting complexes that transfer excitation energy to a reaction center (RC) via the surrounding pigment-protein complex LH1. Recent electron microscopy and atomic force microscopy studies have revealed that RC-LH1 units of Rhodobacter sphaeroides form membrane-bending dimeric complexes together with the polypeptide PufX. We present a structural model for these RC-LH1-PufX dimeric complexes constructed using the molecular dynamics flexible fitting method based on an EM density map. The arrangement of the LH1 BChls displays a distortion near the proposed location of the PufX polypeptide. The resulting atomic model for BChl arrays is used to compute the excitonic properties of the dimeric RC-LH1 complex. A comparison is presented between the structural and excitonic features of the S-shaped dimeric BChl array of Rhodobacter sphaeroides and the circular BChl arrangement found in other purple bacteria.

Microscopic mechanics of hairpin DNA translocation through synthetic nanopores

Nanoscale pores have proved useful as a means to assay DNA and are actively being developed as the basis of genome sequencing methods. Hairpin DNA (hpDNA), having both double-helical and overhanging coil portions, can be trapped in a nanopore, giving ample time to execute a sequence measurement. In this article, we provide a detailed account of hpDNA interaction with a synthetic nanopore obtained through extensive all-atom molecular dynamics simulations. For synthetic pores with minimum diameters from 1.3 to 2.2 nm, we find that hpDNA can translocate by three modes: unzipping of the double helix and--in two distinct orientations--stretching/distortion of the double helix. Furthermore, each of these modes can be selected by an appropriate choice of the pore size and voltage applied transverse to the membrane. We demonstrate that the presence of hpDNA can dramatically alter the distribution of ions within the pore, substantially affecting the ionic current through it. In experiments and simulations, the ionic current relative to that in the absence of DNA can drop below 10% and rise beyond 200%. Simulations associate the former with the double helix occupying the constriction and the latter with accumulation of DNA that has passed through the constriction.

Nonexponential kinetics of DNA escape from alpha-hemolysin nanopores

Throughput and resolution of DNA sequence detection technologies employing nanometer scale pores hinge on accurate kinetic descriptions of DNA motion in nanopores. We present the first detailed experimental study of DNA escape kinetics from alpha-hemolysin nanopores and show that anomalously long escape times for some events result in nonexponential kinetics. From the distribution of first-passage times, we determine that the energy barrier to escape follows a Poisson-like distribution, most likely due to stochastic weak binding events between the DNA and amino acid residues in the pore.

Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics

A novel method to flexibly fit atomic structures into electron microscopy (EM) maps using molecular dynamics simulations is presented. The simulations incorporate the EM data as an external potential added to the molecular dynamics force field, allowing all internal features present in the EM map to be used in the fitting process, while the model remains fully flexible and stereochemically correct. The molecular dynamics flexible fitting (MDFF) method is validated for available crystal structures of protein and RNA in different conformations; measures to assess and monitor the fitting process are introduced. The MDFF method is then used to obtain high-resolution structures of the E. coli ribosome in different functional states imaged by cryo-EM.

Excursion of a single polypeptide into a protein pore: simple physics, but complicated biology

Despite its fundamental and critical importance in molecular biology and practical medical biotechnology, how a polypeptide interacts with a transmembrane protein pore is not yet comprehensively understood. Here, we employed single-channel electrical recordings to reveal the interactions of short polypeptides and small folded proteins with a robust beta-barrel protein pore. The short polypeptides were approximately 25 residues in length, resembling positively charged targeting presequences involved in protein import. The proteins were consisted of positively charged pre-cytochrome b2 fragments (pb2) fused to the small ribonuclease barnase (approximately 110 residues, Ba). Single-molecule experiments exploring the interaction of a folded pb2-Ba protein with a single beta-barrel pore, which contained negatively charged electrostatic traps, revealed the complexity of a network of intermolecular forces, including driving and electrostatic ones. In addition, the interaction was dependent on other factors, such as the hydrophobic content of the interacting polypeptide, the location of the electrostatic trap, the length of the pb2 presequence and temperature. This single-molecule approach together with protein design of either the interacting polypeptide or the pore lumen opens new opportunities for the exploration of the polypeptide-pore interaction at high temporal resolution. Such future studies are also expected to unravel the advantages and limitations of the nanopore technique for the detection and exploration of individual polypeptides.