@article {389, title = {A DNA turbine powered by a transmembrane potential across a nanopore}, journal = {Nature Nanotechnology}, volume = {19}, year = {2024}, month = {03/2024}, pages = {338--344}, chapter = {338}, abstract = {

Rotary motors play key roles in energy transduction, from macroscale windmills to nanoscale turbines such as ATP synthase in cells. Despite our abilities to construct engines at many scales, developing functional synthetic turbines at the nanoscale has remained challenging. Here, we experimentally demonstrate rationally designed nanoscale DNA origami turbines with three chiral blades. These DNA nanoturbines are 24{\textendash}27 nm in height and diameter and can utilize transmembrane electrochemical potentials across nanopores to drive DNA bundles into sustained unidirectional rotations of up to 10 revolutions s-1. The rotation direction is set by the designed chirality of the turbine. All-atom molecular dynamics simulations show how hydrodynamic flows drive this turbine. At high salt concentrations, the rotation direction of turbines with the same chirality is reversed, which is explained by a change in the anisotropy of the electrophoretic mobility. Our artificial turbines operate autonomously in physiological conditions, converting energy from naturally abundant electrochemical potentials into mechanical work. The results open new possibilities for engineering active robotics at the nanoscale.

}, issn = {1748-3395}, doi = {10.1038/s41565-023-01527-8}, url = {https://doi.org/10.1038/s41565-023-01527-8}, author = {Shi, Xin and Pumm, Anna-Katharina and Maffeo, Christopher and Kohler, Fabian and Feigl, Elija and Zhao, Wenxuan and Verschueren, Daniel and Golestanian, Ramin and Aksimentiev, Aleksei and Dietz, Hendrik and Dekker, Cees} } @article {398, title = {Electro-Osmotic Flow Generation via a Sticky Ion Action}, journal = {ACS nano}, volume = {18}, year = {2024}, pages = {17521-17533}, chapter = {17521-17533}, abstract = {

Selective transport of ions through nanometer-sized pores is fundamental to cell biology and central to many technological processes such as water desalination and electrical energy storage. Conventional methods for generating ion selectivity include placement of fixed electrical charges at the inner surface of a nanopore through either point mutations in a protein pore or chemical treatment of a solid-state nanopore surface, with each nanopore type requiring a custom approach. Here, we describe a general method for transforming a nanoscale pore into a highly selective, anion-conducting channel capable of generating a giant electro-osmotic effect. Our molecular dynamics simulations and reverse potential measurements show that exposure of a biological nanopore to high concentrations of guanidinium chloride renders the nanopore surface positively charged due to transient binding of guanidinium cations to the protein surface. A comparison of four biological nanopores reveals the relationship between ion selectivity, nanopore shape, composition of the nanopore surface, and electro-osmotic flow. Guanidinium ions are also found to produce anion selectivity and a giant electro-osmotic flow in solid-state nanopores via the same mechanism. Our sticky-ion approach to generate electro-osmotic flow can have numerous applications in controlling molecular transport at the nanoscale and for detection, identification, and sequencing of individual proteins.

}, doi = {10.1021/acsnano.4c00829}, author = {Mehrafrooz, Behzad and Yu, Luning and Pandey, Laxmi and Siwy, Zuzanna S and Wanunu, Meni and Aksimentiev, Aleksei} } @article {390, title = {Engineering Biological Nanopore Approaches toward Protein Sequencing}, journal = {ACS Nano}, volume = {17}, year = {2023}, month = {09}, pages = {16369{\textendash}16395}, abstract = {

Biotechnological innovations have vastly improved the capacity to perform large-scale protein studies, while the methods we have for identifying and quantifying individual proteins are still inadequate to perform protein sequencing at the single-molecule level. Nanopore-inspired systems devoted to understanding how single molecules behave have been extensively developed for applications in genome sequencing. These nanopore systems are emerging as prominent tools for protein identification, detection, and analysis, suggesting realistic prospects for novel protein sequencing. This review summarizes recent advances in biological nanopore sensors toward protein sequencing, from the identification of individual amino acids to the controlled translocation of peptides and proteins, with attention focused on device and algorithm development and the delineation of molecular mechanisms with the aid of simulations. Specifically, the review aims to offer recommendations for the advancement of nanopore-based protein sequencing from an engineering perspective, highlighting the need for collaborative efforts across multiple disciplines. These efforts should include chemical conjugation, protein engineering, molecular simulation, machine-learning-assisted identification, and electronic device fabrication to enable practical implementation in real-world scenarios.

}, isbn = {1936-0851}, doi = {10.1021/acsnano.3c05628}, url = {https://doi.org/10.1021/acsnano.3c05628}, author = {Wei, Xiaojun and Penkauskas, Tadas and Reiner, Joseph E. and Kennard, Celeste and Uline, Mark J. and Wang, Qian and Li, Sheng and Aksimentiev, Aleksei and Robertson, Joseph W. F. and Liu, Chang} } @article {366, title = {Expanding the Molecular Alphabet of DNA-Based Data Storage Systems with Neural Network Nanopore Readout Processing}, journal = {Nano Letters}, year = {2022}, month = {02/2022}, abstract = {

DNA is a promising next-generation data storage medium, but challenges remain with synthesis costs and recording latency. Here, we describe a prototype of a DNA data storage system that uses an extended molecular alphabet combining natural and chemically modified nucleotides. Our results show that MspA nanopores can discriminate different combinations and ordered sequences of natural and chemically modified nucleotides in custom-designed oligomers. We further demonstrate singlemolecule sequencing of the extended alphabet using a neural network architecture that classifies raw current signals generated by Oxford Nanopore sequencers with an average accuracy exceeding 60\% (39{\texttimes} larger than random guessing). Molecular dynamics simulations show that the majority of modified nucleotides lead to only minor perturbations of the DNA double helix. Overall, the extended molecular alphabet may potentially offer a nearly 2-fold increase in storage density and potentially the same order of reduction in the recording latency, thereby enabling new implementations of molecular recorders.

}, doi = {10.1021/acs.nanolett.1c04203}, author = {Tabatabaei, S Kasra and Pham, Bach and Pan, Chao and Liu, Jingqian and Chandak, Shubham and Shorkey, Spencer A and Hernandez, Alvaro G and Aleksei Aksimentiev and Chen, Min and Schroeder, Charles M and Milenkovic, Olgica} } @article {369, title = {DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles.}, journal = {Nucleic Acids Res}, volume = {49}, year = {2021}, month = {11/2021}, pages = {11459-11475}, abstract = {

Eukaryotic genome and methylome encode DNA fragments{\textquoteright} propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between {\textquoteright}omics{\textquoteright} and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles{\textquoteright} local geometry consists of straight segments connected by sharp bends compressing the DNA{\textquoteright}s inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles{\textquoteright} inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes{\textquoteright} nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.

}, keywords = {DNA Methylation, DNA, Circular, Molecular Dynamics Simulation, Nucleic Acid Conformation}, issn = {1362-4962}, doi = {10.1093/nar/gkab967}, author = {Yoo, Jejoong and Park, Sangwoo and Maffeo, Christopher and Ha, Taekjip and Aleksei Aksimentiev} } @article {356, title = {The emerging landscape of single-molecule protein sequencing technologies}, journal = {Nature Methods}, volume = {18}, year = {2021}, month = {06/2021}, chapter = {604-617}, abstract = {

Single-cell profiling methods have had a profound impact on the understanding of cellular heterogeneity. While genomes and transcriptomes can be explored at the single-cell level, single-cell profiling of proteomes is not yet established. Here we describe new single-molecule protein sequencing and identification technologies alongside innovations in mass spectrometry that will eventually enable broad sequence coverage in single-cell profiling. These technologies will in turn facilitate biological discovery and open new avenues for ultrasensitive disease diagnostics.

}, doi = {10.1038/s41592-021-01143-1}, url = {https://doi.org/10.1038/s41592-021-01143-1}, author = {Javier Antonio Alfaro and Peggy Bohl{\"a}nder and Mingjie Dai and Mike Filius and Cecil J. Howard and Xander F. van Kooten and Shilo Ohayon and Adam Pomorski and Sonja Schmid and Aleksei Aksimentiev and Eric V. Anslyn and Georges Bedran and Chan Cao and Mauro Chinappi and Etienne Coyaud and Cees Dekker and Gunnar Dittmar and Nicholas Drachman and Rienk Eelkema and David Goodlett and S{\'e}bastien Hentz and Umesh Kalathiya and Neil L. Kelleher and Ryan T. Kelly and Zvi Kelman and Sung Hyun Kim and Bernhard Kuster and David Rodriguez-Larrea and Stuart Lindsay and Giovanni Maglia and Edward M. Marcotte and John P. Marino and Christophe Masselon and Michael Mayer and Patroklos Samaras and Kumar Sarthak and Lusia Sepiashvili and Derek Stein and Meni Wanunu and Mathias Wilhelm and Peng Yin and Amit Meller and Chirlmin Joo} } @article {335, title = {Scalable molecular dynamics on CPU and GPU architectures with NAMD}, journal = {J. Chem. Phys.}, volume = {153}, year = {2020}, month = {07/2020}, chapter = {044130}, abstract = {

NAMD is a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU- based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolec- ular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics\ simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchem- ical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

}, doi = {10.1063/5.0014475}, author = {James C. Phillips and David J. Hardy and Julio D. C. Maia and John E. Stone and Jo{\~a}o V. Ribeiro and Rafael C. Bernardi and Ronak Buch and Giacomo Fiorin and J{\'e}r{\^o}me H{\'e}nin and Wei Jiang and Ryan McGreevy and Marcelo C. R. Melo and Brian K. Radak and Robert D. Skeel and Abhishek Singharoy and Yi Wang and Beno{\^\i}t Roux and Aleksei Aksimentiev and Zaida Luthey-Schulten and Laxmikant V. Kal{\'e} and Klaus Schulten and Christophe Chipot and Emad Tajkhorshid} } @article {317, title = {Electrical recognition of the twenty proteinogenic amino acids using an aerolysin nanopore}, journal = {Nature Biotechnology}, volume = {38}, year = {2020}, month = {02/2020}, pages = {176-181}, chapter = {176}, abstract = {

Efforts to sequence single protein molecules in nanopores have been hampered by the lack of techniques with sufficient sensitivity to discern the subtle molecular differences among all twenty amino acids. Here we report ionic current detection of all twenty proteinogenic amino acids in an aerolysin nanopore with the help of a short polycationic carrier. Application of molecular dynamics simulations revealed that the aerolysin nanopore has a built-in single-molecule trap that fully confines a polycationic carrier-bound amino acid inside the sensing region of the aerolysin. This structural feature means that each amino acid spends sufficient time in the pore for sensitive measurement of the excluded volume of the amino acid. We show that distinct current blockades in wild-type aerolysin can be used to identify 13 of the 20 natural amino acids. Furthermore, we show that chemical modifications, instrumentation advances and nanopore engineering offer a route toward identification of the remaining seven amino acids. These findings may pave the way to nanopore protein sequencing.

Research Highlights

1.\ Tang, Lei. "Sensing proteinogenic amino acids." Nature Methods 17(2):126 (2020) DOI:\ 10.1038/s41592-020-0741-z

2.\ Howorka, Stefan, and Zuzanna S. Siwy. "Reading amino acids in a nanopore." Nature Biotechnology 38(2):159 (2020) DOI:\ 10.1038/s41587-019-0401-y

3.\ Nanopores can identify the amino acids in proteins, the first step to sequencing, Highlight by Illinois news bureau.\ \ 

}, isbn = {1546-1696}, doi = {10.1038/s41587-019-0345-2}, url = {https://doi.org/10.1038/s41587-019-0345-2}, author = {Ouldali, Hadjer and Sarthak, Kumar and Ensslen, Tobias and Piguet, Fabien and Manivet, Philippe and Pelta, Juan and Behrends, Jan C. and Aleksei Aksimentiev and Oukhaled, Abdelghani} } @article {318, title = {Artificial Water Channels Enable Fast and Selective Water Permeation Through Water-wire Networks}, journal = {Nature Nanotechnology.}, volume = {15,}, year = {2020}, month = {01/2020}, pages = {73-79}, chapter = {73}, abstract = {

Artificial water channels are synthetic molecules that aim to mimic the structural and functional features of biological water channels (aquaporins). Here we report on a cluster-forming organic nanoarchitecture, peptide-appended hybrid[4]arene (PAH[4]), as a new class of artificial water channels. Fluorescence experiments and simulations demonstrated that PAH[4]s can form, through lateral diffusion, clusters in lipid membranes that provide synergistic membrane-spanning paths for rapid and selective water permeation through water-wire networks. Quantitative transport studies revealed that PAH[4]s can transport \>109\ water molecules per second per molecule, which is comparable to aquaporin water channels. The performance of these channels exceeds the upper bound limit of current desalination membranes by a factor of ~104, as illustrated by the water/NaCl permeability{\textendash}selectivity trade-off curve. PAH[4]{\textquoteright}s unique properties of a high water/solute permselectivity via cooperative water-wire formation could usher in an alternative design paradigm for permeable membrane materials in separations, energy production, and barrier applications.

}, doi = {https://doi.org/10.1038/s41565-019-0586-8}, author = {Woochul Song and Himanshu Joshi and Ratul Chowdhury and Joseph S. Najem and Yue-xiao Shen and Chao Lang and Codey B. Henderson and Yu-Ming Tu and Megan Farell and Megan E. Pitz and Costas D. Maranas and Paul S. Cremer and Robert J. Hickey and Stephen A. Sarles and Jun-li Hou and Aleksei Aksimentiev and Manish Kumar} } @article {308, title = {Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism}, journal = {Cell}, volume = {179}, year = {2019}, month = {11/2019}, pages = {1098 - 1111.e23}, abstract = {

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore{\textquoteright}s structural integrity and robust energy conversion. Parallels are drawn with the more universal\ mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

}, keywords = {bioenergetics, biological membranes, charge transport, chromatophore, integrative model, MD, mitochondria, Molecular Dynamics Simulation, optical spectroscopy, photosynthesis}, issn = {0092-8674}, doi = {10.1016/j.cell.2019.10.021}, url = {http://www.sciencedirect.com/science/article/pii/S0092867419311717}, author = {Abhishek Singharoy and Christopher Maffeo and Karelia H. Delgado-Magnero and David J.K. Swainsbury and Melih Sener and Ulrich Kleinekath{\"o}fer and John W. Vant and Jonathan Nguyen and Andrew Hitchcock and Barry Isralewitz and Ivan Teo and Danielle E. Chandler and John E. Stone and James C. Phillips and Taras V. Pogorelov and M. Ilaria Mallus and Christophe Chipot and Zaida Luthey-Schulten and D. Peter Tieleman and C. Neil Hunter and Emad Tajkhorshid and Aleksei Aksimentiev and Klaus Schulten} } @article {301, title = {The Oligomerization Landscape of Histones}, journal = {Biophysical Journal}, volume = {116}, year = {2019}, month = {05/2019}, pages = {1845--1855}, doi = {doi.org/10.1016/j.bpj.2019.03.021}, url = {https://www.sciencedirect.com/science/article/pii/S0006349519302784}, author = {Zhao, Haiqing and Winogradoff, David and Dalal, Yamini and Papoian, Garegin A} } @article {290, title = {PoreDesigner for tuning solute selectivity in a robust and highly permeable outer membrane pore}, journal = {Nature Communications}, volume = {9}, year = {2018}, month = {09/2018}, pages = {3661}, abstract = {

Monodispersed angstrom-size pores embedded in a suitable matrix are promising for highly selective membrane-based separations. They can provide substantial energy savings in water treatment and small molecule bioseparations. Such pores present as membrane proteins (chiefly aquaporin-based) are commonplace in biological membranes but difficult to implement in synthetic industrial membranes and have modest selectivity without tunable selectivity. Here we present PoreDesigner, a design workflow to redesign the robust beta-barrel Outer Membrane Protein F as a scaffold to access three specific pore designs that exclude solutes larger than sucrose (\>360 Da), glucose (\>180 Da), and salt (\>58 Da) respectively. PoreDesigner also enables us to design any specified pore size (spanning 3{\textendash}10 {\r A}), engineer its pore profile, and chemistry. These redesigned pores may be ideal for conducting sub-nm aqueous separations with permeabilities exceeding those of classical biological water channels,\ aquaporins, by more than an order of magnitude at over 10 billion water molecules per channel per second.

}, issn = {2041-1723}, doi = {10.1038/s41467-018-06097-1}, url = {https://doi.org/10.1038/s41467-018-06097-1}, author = {Chowdhury, Ratul and Ren, Tingwei and Shankla, Manish and Decker, Karl and Grisewood, Matthew and Prabhakar, Jeevan and Baker, Carol and Golbeck, John H and Aleksei Aksimentiev and Kumar, Manish and Maranas, Costas D} } @article {279, title = {Porphyrin-Assisted Docking of a Thermophage Portal Protein into Lipid Bilayers: Nanopore Engineering and Characterization}, journal = {ACS Nano}, volume = {11}, year = {2017}, month = {12/2017}, pages = {11931-11945}, abstract = {

Nanopore-based sensors for nucleic acid sequencing and single-molecule detection typically employ pore-forming membrane proteins with hydrophobic external surfaces, suitable for insertion into a lipid bilayer. In contrast, hydrophilic pore-containing molecules such as DNA origami, have been shown to require chemical modification to favor insertion into a lipid environment. In this work, we describe a strategy for inserting polar proteins with an inner pore into lipid membranes, focusing here on a circular 12-subunit assembly of the thermophage G20c portal protein. X-ray crystallography, electron microscopy, molecular dynamics and thermal/chaotrope denaturation experiments all find the G20c portal protein to have a highly stable structure, favorable for nanopore sensing applications. Porphyrin conjugation to a cysteine mutant in the protein facilitates the protein{\textquoteright}s insertion into lipid bilayers, allowing us to probe ion transport through the pore. Finally, we probed the portal interior size and shape using a series of cyclodextrins of varying sizes, revealing asymmetric transport that possibly originates from the portal{\textquoteright}s DNA-ratchet function.

}, doi = {10.1021/acsnano.7b06980}, url = {http://pubs.acs.org/doi/abs/10.1021/acsnano.7b06980}, author = {Cressiot*, Benjamin and Greive*, Sandra J. and Si*, Wei and Pascoa, Tomas and Mojtabavi, Mehrnaz and Chechik, Maria and Jenkins, Huw Thomas and Lu, Xueguang and Zhang, Ke and Aleksei Aksimentiev and Antson, Alfred A. and Wanunu, Meni} } @article {274, title = {Quantification of Membrane Protein-Detergent Complex Interactions}, journal = {The Journal of Physical Chemistry B}, volume = {121}, year = {2017}, month = {11/2017}, pages = {10228-10241}, chapter = {10228}, abstract = {

Although fundamentally significant in structural, chemical, and membrane biology, the interfacial protein-detergent complex (PDC) interactions have been modestly examined because of the complicated behavior of both detergents and membrane proteins in aqueous phase. Membrane proteins are prone to unproductive aggregation resulting from poor detergent solvation, but the participating forces in this phenomenon remain ambiguous. Here, we show that using rational membrane protein design, targeted chemical modification, and steady-state fluorescence polarization spectroscopy, the detergent desolvation of membrane proteins can be quantitatively evaluated. We demonstrate that depleting the detergent in the sample well produced a two-state transition of membrane proteins between a fully detergent-solvated state and a detergent-desolvated state, the nature of which depended on the interfacial PDC interactions. Using a panel of six membrane proteins of varying hydrophobic topography, structural fingerprint, and charge distribution on the solvent-accessible surface, we provide direct experimental evidence for the contributions of the electrostatic and hydrophobic interactions to the protein solvation properties. Moreover, all-atom molecular dynamics simulations report the major contribution of the hydrophobic forces exerted at the PDC interface. This semi-quantitative approach might be extended in the future to include studies of the interfacial PDC interactions of other challenging membrane protein systems of unknown structure. This would have practical importance in protein extraction, solubilization, stabilization, and crystallization.

}, doi = {10.1021/acs.jpcb.7b08045}, url = {http://pubs.acs.org/doi/abs/10.1021/acs.jpcb.7b08045}, author = {Wolfe, Aaron John and Si, Wei and Zhang, Zhengqi and Blanden, Adam R and Hsueh, Yi-Ching and Gugel, Jack Fredrick and Pham, Bach and Chen, Min and Loh, Stewart N and Rozovsky, Sharon and Aleksei Aksimentiev and Movileanu, Liviu} } @article {268, title = {Nanoscale Ion Pump Derived from a Biological Water Channel}, journal = {The Journal of Physical Chemistry B }, volume = {121}, year = {2017}, month = {08/2017}, pages = {7899{\textendash}7906 }, abstract = {

Biological molecular machines perform the\ work of supporting life at the smallest of scales, including\ the work of shuttling ions across cell boundaries and against\ chemical gradients. Systems of artifi cial channels at the\ nanoscale can likewise control ionic concentration by way of\ ionic current rectifi cation, species selectivity, and voltage\ gating mechanisms. Here, we theoretically show that a voltagegated,\ ion species-selective, and rectifying ion channel can be\ built using the components of a biological water channel\ aquaporin. Through all-atom molecular dynamics simulations,\ we show that the ionic conductance of a truncated aquaporin\ channel nonlinearly increases with the bias magnitude, depends on the channel{\textquoteright} s orientation, and is highly cation specific but only\ for one polarity of the transmembrane bias. Further, we show that such an unusually complex response of the channel to\ transmembrane bias arises from mechanical motion of a positively charged gate that blocks cation transport. By combining two\ truncated aquaporins, we demonstrate a molecular system that pumps ions against their chemical gradients when subject to an\ alternating transmembrane bias. Our work sets the stage for future biomimicry eff orts directed toward reproducing the function\ of biological ion pumps using synthetic components.

}, doi = {10.1021/acs.jpcb.7b05568}, author = {Karl Decker and Martin Page and Aleksei Aksimentiev} } @article {251, title = {Selective Permeability of Truncated Aquaporin 1 in Silico}, journal = {ACS Biomaterials Science \& Engineering}, volume = {3}, year = {2017}, month = {02/2017}, pages = {342-348}, chapter = {342}, abstract = {

Aquaporin (AQP) proteins function as highly efficient water transport channels that support homeostasis in many types of living cells. Their structure-function relationships have been characterized extensively in fundamental and applied research, primarily via structural analysis, mutational studies, and computational approaches. The present study evaluates the effects of progressive truncations on the permeability and ionic conductivity of AQP-1 (bovine). The use of truncations to determine critical features has not been considered previously, as physical truncation of AQP is likely not technically feasible due to the ornate arrangement of six interwoven alpha helices in a single pore structure. However, structures not obtainable through protein assembly can be realized via synthetic chemistry approaches and studied using molecular dynamics (MD) simulations. Here, we apply the MD method to characterize the permeability of AQP variants truncated along the pore axis from both cytoplasmic and extracellular sides of the channel. The simulation results suggest that AQP-1 can retain its function even after deletion of up to 50\% of the channel{\textquoteright}s length, representing 50\% of proteins{\textquoteright}s molecular mass. Deletions such as these are expected to greatly simplify future biomimicry efforts of reproducing the AQP functionality using synthetic macromolecules. This study demonstrates the potential of in silico approaches to support the creation of streamlined functional analogs of biological molecular machines.

}, doi = {10.1021/acsbiomaterials.6b00583 }, author = {Karl Decker and Martin Page and Ashley Boyd and Irene MacAllister and Mark Ginsberg and Aleksei Aksimentiev} } @article {250, title = {Mechanical Trapping of DNA in a Double-Nanopore System}, journal = {Nano Letters}, volume = {16}, year = {2016}, month = {11/2016}, pages = {8021-8028 }, abstract = {

Nanopores have become ubiquitous components of systems for single-molecule manipulation and detection, in particular DNA sequencing where electric field driven translocation of DNA through a nanopore is used to read out the DNA molecule. Here, we present a double-pore system where two nanopores are drilled in parallel through the same solid-state membrane, which offers new opportunities for DNA manipulation. Our experiments and molecular dynamics simulations show that simultaneous electrophoretic capture of a DNA molecule by the two nanopores mechanically traps the molecule, increasing its residence time within the nanopores by orders of magnitude. Remarkably, by using two unequal-sized nanopores, the pore of DNA entry and exit can be discerned from the ionic current blockades, and the translocation direction can be precisely controlled by small differences in the effective force applied to DNA. The mechanical arrest of DNA translocation using a double-pore system can be straightforwardly integrated into any solid-state nanopore platform, including those using optical or transverse-current readouts.

}, doi = {10.1021/acs.nanolett.6b04642}, author = {Sergii Pud and Shu-Han Chao and Maxim Belkin and Daniel Verschueren and Teun Huijben and Casper van Engelenburg and Cees Dekker and Aleksei Aksimentiev} } @article {316, title = {Promiscuous Histone Mis-Assembly Is Actively Prevented by Chaperones.}, journal = {J Am Chem Soc}, volume = {138}, year = {2016}, month = {10/2016}, pages = {13207-13218}, abstract = {

Histone proteins are essential for the organization, expression, and inheritance of genetic material for eukaryotic cells. A centromere-specific H3 histone variant, centromere protein A (CENP-A), shares about 50\% amino acid sequence identity with H3. CENP-A is required for packaging the centromere and for the proper separation of chromosomes during mitosis. Despite their distinct biological functions, previously reported crystal structures of the CENP-A/H4 and H3/H4 dimers reveal a high degree of similarity. In this work, we characterize the structural dynamics of CENP-A/H4 and H3/H4 dimers based on a dual-resolution approach, using both microsecond-scale explicit-solvent all-atom and coarse-grained (CG) molecular dynamics (MD) simulations. Our data show that the H4 histone is significantly more rigid compared with the H3 histone and its variant CENP-A, hence, serving as a reinforcing structural element within the histone core. We report that the CENP-A/H4 dimer is significantly more dynamic than its canonical counterpart H3/H4, and our results provide a physical explanation for this flexibility. Further, we observe that the centromere-specific chaperone Holliday Junction Recognition Protein (HJURP) stabilizes the CENP-A/H4 dimer by forming a specific electrostatic interaction network. Finally, replacing CENP-A S68 with E68 disrupts the binding interface between CENP-A and HJURP in all-atom MD simulation, and consistently, in vivo experiments demonstrate that replacing CENP-A S68 with E68 disrupts CENP-A{\textquoteright}s localization to the centromere. Based on all our results, we propose that, during the CENP-A/H4 deposition process, the chaperone HJURP protects various substructures of the dimer, serving both as a folding and binding chaperone.

}, keywords = {Centromere Protein A, DNA-Binding Proteins, HeLa Cells, Histones, Humans, Molecular Dynamics Simulation, Protein Folding, Protein Multimerization, Protein Structure, Quaternary, Substrate Specificity}, doi = {10.1021/jacs.6b05355}, author = {Zhao*, Haiqing and Winogradoff*, David and Bui, Minh and Dalal, Yamini and Papoian, Garegin A} } @article {313, title = {Promiscuous histone mis-assembly is actively prevented by chaperones.}, journal = {Journal of the American Chemical Society}, volume = {138}, year = {2016}, pages = {13207--13218}, doi = {10.1021/jacs.6b05355}, author = {Zhao*, Haiqing and Winogradoff*, David and Minh Bui and Yamini Dalal and Garegin Papoian} } @article {208, title = {Mapping fast protein folding with multiple-site fluorescent probes.}, journal = {Proc Natl Acad Sci U S A}, volume = {112}, year = {2015}, month = {06/2015}, pages = {7966-71}, abstract = {

Fast protein folding involves complex dynamics in many degrees of freedom, yet microsecond folding experiments provide only low-resolution structural information. We enhance the structural resolution of the five-helix bundle protein λ6-85 by engineering into it three fluorescent tryptophan-tyrosine contact probes. The probes report on distances between three different helix pairs: 1-2, 1-3, and 3-2. Temperature jump relaxation experiments on these three mutants reveal two different kinetic timescales: a slower timescale for 1-3 and a faster one for the two contacts involving helix 2. We hypothesize that these differences arise from a single folding mechanism that forms contacts on different timescales, and not from changes of mechanism due to adding the probes. To test this hypothesis, we analyzed the corresponding three distances in one published single-trajectory all-atom molecular-dynamics simulation of a similar mutant. Autocorrelation analysis of the trajectory reveals the same "slow" and "fast" distance change as does experiment, but on a faster timescale; smoothing the trajectory in time shows that this ordering is robust and persists into the microsecond folding timescale. Structural investigation of the all-atom computational data suggests that helix 2 misfolds to produce a short-lived off-pathway trap, in agreement with the experimental finding that the 1-2 and 3-2 distances involving helix 2 contacts form a kinetic grouping distinct from 1 to 3. Our work demonstrates that comparison between experiment and simulation can be extended to several order parameters, providing a stronger mechanistic test.

}, keywords = {Fluorescent Dyes, Kinetics, Molecular Dynamics Simulation, Mutation, Protein Folding, Tryptophan, Tyrosine}, issn = {1091-6490}, doi = {10.1073/pnas.1422683112}, url = {http://www.pnas.org/content/112/26/7966.full?sid=da9112c1-35df-4a1b-83b8-cfe4c0c56f7b}, author = {Prigozhin*, Maxim B and Chao*, Shu-Han and Sukenik, Shahar and Pogorelov, Taras V and Gruebele, Martin} } @article {220, title = {Shearing of the CENP-A dimerization interface mediates plasticity in the octameric centromeric nucleosome}, journal = {Scientific reports}, volume = {5}, year = {2015}, doi = {10.1038/srep17038}, author = {Winogradoff, David and Zhao, Haiqing and Dalal, Yamini and Papoian, Garegin A} } @article {219, title = {The acetylation landscape of the H4 histone tail: disentangling the interplay between the specific and cumulative effects}, journal = {J Am Chem Soc}, volume = {137}, year = {2015}, pages = {6245{\textendash}6253}, doi = {10.1021/jacs.5b00235}, author = {Winogradoff, David and Echeverria, Ignacia and Potoyan, Davit A and Papoian, Garegin A} } @article {119, title = {Two Structural Scenarios for Protein Stabilization by PEG.}, journal = {J Phys Chem B}, volume = {118}, year = {2014}, month = {05/2014}, pages = {8388{\textendash}8395}, abstract = {

PEGylation, or addition of poly(ethylene glycol) chains to proteins, is widely used to improve delivery in pharmaceutical applications. Recent studies suggest that stabilization of a protein by PEG, and hence its proteolytic degradability, are sequence-dependent and require only short PEG chains. Here we connect stabilization by short PEG chains directly to the structural dynamics of the protein and PEG chain. We measured the stability of human Pin1 WW domain with PEG-4 at asparagine 19 for a full mutant cycle at two positions thought to influence PEG-protein interaction: Ser16Ala and Tyr23Phe. We then performed explicit solvent molecular dynamics simulations on all PEGylated and PEG-free mutants. The mutant cycle yields a non-additive stabilization effect where the pseudo wild type and double mutant are more stabilized relative to unPEGylated proteins than are the two single mutants. The simulation reveals why: the double mutant suffers loss of beta sheet structure, which PEGylation restores even though the PEG extends as a coil into the solvent. In contrast, in one of the single mutants PEG preferentially interacts with the protein surface while disrupting the interactions of its asparagine host with a nearby methionine side chain. Thus PEG attachment can stabilize a protein differentially depending on the local sequence, and either by interacting with the surface, or by extending into the solvent. A simulation with PEG-45 attached to asparagine 19 shows that PEG even can do both in the same context.

}, issn = {1520-5207}, doi = {10.1021/jp502234s}, author = {Shu-Han Chao and Matthews, Sam S and Paxman, Ryan and Aleksei Aksimentiev and Martin Gruebele and Price, Joshua L} } @article {121, title = {Filling Up the Heme Pocket Stabilizes Apomyoglobin and Speeds Up Its Folding.}, journal = {J Phys Chem B}, year = {2014}, month = {02/2014}, abstract = {

Wild type apomyoglobin folds in at least two steps: the ABGH core rapidly, followed much later by the heme-binding CDEF core. We hypothesize that the evolved heme-binding function of the CDEF core frustrates its folding: it has a smaller contact order and is no more complex topologically than ABGH, and thus, it should be able to fold faster. Therefore, filling up the empty heme cavity of apomyoglobin with larger, hydrophobic side chains should significantly stabilize the protein and increase its folding rate. Molecular dynamics simulations allowed us to design four different mutants with bulkier side chains that increase the native bias of the CDEF region. In vitro thermal denaturation shows that the mutations increase folding stability and bring the protein closer to two-state behavior, as judged by the difference of fluorescence- and circular dichroism-detected protein stability. Millisecond stopped flow measurements of the mutants exhibit refolding kinetics that are over 4 times faster than the wild type{\textquoteright}s. We propose that myoglobin-like proteins not evolved to bind heme are equally stable, and find an example. Our results illustrate how evolution for function can force proteins to adapt frustrated folding mechanisms, despite having simple topologies.

}, issn = {1520-5207}, doi = {10.1021/jp412459z}, author = {Goodman, J S and Shu-Han Chao and Pogorelov, T V and Martin Gruebele} } @article {doi:10.1021/jp501492g, title = {Rectification of Ion Current in Nanopores Depends on the Type of Monovalent Cations: Experiments and Modeling}, journal = {The Journal of Physical Chemistry C}, volume = {118}, number = {18}, year = {2014}, pages = {9809-9819}, abstract = {Rectifying nanopores feature ion currents that are higher for voltages of one polarity compared to the currents recorded for corresponding voltages of the opposite polarity. Rectification of nanopores has been found to depend on the pore opening diameter and distribution of surface charges on the pore walls as well as pore geometry. Very little is known, however, on the dependence of ionic rectification on the type of transported ions of the same charge. We performed experiments with single conically shaped nanopores in a polymer film and recorded current-voltage curves in three electrolytes: LiCl, NaCl, and KCl. Rectification degrees of the pores, quantified as the ratio of currents recorded for voltages of opposite polarities, were the highest for KCl and the lowest for LiCl. The experimental observations could not be explained by a continuum modeling based on the Poisson-Nernst-Planck equations. All-atom molecular dynamics simulations revealed differential binding between Li+, Na+, and K+ ions and carboxyl groups on the pore walls, resulting in changes to both the effective surface charge of the nanopore and cation mobility within the pore.}, doi = {10.1021/jp501492g}, url = {http://dx.doi.org/10.1021/jp501492g}, author = {Gamble, Trevor and Karl Decker and Plett, Timothy S. and Pevarnik, Matthew and Pietschmann, Jan-Frederik and Vlassiouk, Ivan and Aleksei Aksimentiev and Siwy, Zuzanna S.} } @article {26, title = {Molecular dynamics study of MspA arginine mutants predicts slow DNA translocations and ion current blockades indicative of DNA sequence}, journal = {ACS Nano}, volume = {6}, year = {2012}, month = {2012 Aug 28}, pages = {6960-8}, abstract = {

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.

}, keywords = {Arginine, Base Sequence, Biosensing Techniques, Computer Simulation, Conductometry, DNA, Electric Conductivity, Ion Channel Gating, Models, Chemical, Models, Genetic, Models, Molecular, Molecular Sequence Data, Mutation, Porins, Sequence Analysis, DNA, Translocation, Genetic}, issn = {1936-086X}, doi = {10.1021/nn3019943}, author = {Bhattacharya, Swati and Derrington, Ian M and Pavlenok, Mikhail and Niederweis, Michael and Gundlach, Jens H and Aleksei Aksimentiev} } @article {195, title = {Rectification of the current in alpha-hemolysin pore depends on the cation type: the alkali series probed by MD simulations and experiments}, journal = {The Journal of Physical Chemistry C}, volume = {115}, year = {2011}, month = {2011 Feb 21}, pages = {4255-4264}, abstract = {

A striking feature of the alpha-hemolysin channel-a prime candidate for biotechnological applications-is the dependence of its ionic conductance on the magnitude and direction of the applied bias. Through a combination of lipid bilayer single-channel recording and molecular dynamics (MD) simulations, we characterized the current-voltage relationship of alpha-hemolysin for all alkali chloride salts at neutral pH. The rectification of the ionic current was found to depend on the type of cations and increase from Li(+) to Cs(+). Analysis of the MD trajectories yielded a simple quantitative model that related the ionic current to the electrostatic potential, the concentration and effective mobility of ions in the channel. MD simulations reveal that the major contribution to the current asymmetry and rectification properties originates from the cationic contribution to the current that is significantly reduced in a cationic dependent way when the membrane polarity is reversed. The variation of chloride current was found to be less important. We report that the differential affinity of cations for the charged residues positioned at the channel{\textquoteright}s end modulates the number of ions inside the channel stem thus affecting the current properties. Through direct comparison of simulation and experiment, this study evaluates the accuracy of the MD method for prediction of the asymmetric, voltage dependent conductances of a membrane channel.

}, issn = {1932-7447}, doi = {10.1021/jp111441p}, author = {Bhattacharya, Swati and Muzard, L and Payet, Linda and Math{\'e}, J{\'e}r{\^o}me and Bockelmann, Ulrich and Aleksei Aksimentiev and Viasnoff, Virgile} } @article {37, title = {Lipid bilayer coated Al(2)O(3) nanopore sensors: towards a hybrid biological solid-state nanopore}, journal = {Biomed Microdevices}, volume = {13}, year = {2011}, month = {2011 Aug}, pages = {671-82}, abstract = {

Solid-state nanopore sensors are highly versatile platforms for the rapid, label-free electrical detection and analysis of single molecules, applicable to next generation DNA sequencing. The versatility of this technology allows for both large scale device integration and interfacing with biological systems. Here we report on the development of a hybrid biological solid-state nanopore platform that incorporates a highly mobile lipid bilayer on a single solid-state Al(2)O(3) nanopore sensor, for the potential reconstitution of ion channels and biological nanopores. Such a system seeks to combine the superior electrical, thermal, and mechanical stability of Al(2)O(3) solid-state nanopores with the chemical specificity of biological nanopores. Bilayers on Al(2)O(3) exhibit higher diffusivity than those formed on TiO(2) and SiO(2) substrates, attributed to the presence of a thick hydration layer on Al(2)O(3), a key requirement to preserving the biological functionality of reconstituted membrane proteins. Molecular dynamics simulations demonstrate that the electrostatic repulsion between the dipole of the DOPC headgroup and the positively charged Al(2)O(3) surface may be responsible for the enhanced thickness of this hydration layer. Lipid bilayer coated Al(2)O(3) nanopore sensors exhibit excellent electrical properties and enhanced mechanical stability (GΩ seals for over 50 h), making this technology ideal for use in ion channel electrophysiology, the screening of ion channel active drugs and future integration with biological nanopores such as α-hemolysin and MspA for rapid single molecule DNA sequencing. This technology can find broad application in bio-nanotechnology.

}, keywords = {Aluminum Oxide, Biosensing Techniques, DNA, Fluorescence Recovery After Photobleaching, Ion Channels, Lipid Bilayers, Models, Biological, Molecular Dynamics Simulation, Nanopores, Nanotechnology, Phosphatidylcholines, Sequence Analysis, DNA, Silicon Dioxide, Staining and Labeling, Titanium}, issn = {1572-8781}, doi = {10.1007/s10544-011-9537-3}, author = {Venkatesan, Bala Murali and Polans, James and Jeffrey Comer and Sridhar, Supriya and Wendell, David and Aleksei Aksimentiev and Rashid Bashir} } @article {77, title = {Large scale simulation of protein mechanics and function}, journal = {Adv Protein Chem}, volume = {66}, year = {2003}, month = {2003}, pages = {195-247}, keywords = {Amino Acid Motifs, Aquaporins, Catalysis, Cell Membrane, Computer Simulation, Fibronectins, Hydrogen Bonding, Models, Molecular, Protein Conformation, Protein Denaturation, Protein Folding, Protein Structure, Secondary, Proteins, Proton-Translocating ATPases, Substrate Specificity, Time Factors, Water}, issn = {0065-3233}, author = {Tajkhorshid, Emad and Aleksei Aksimentiev and Ilya A Balabin and Gao, Mu and Isralewitz, Barry and Phillips, James C and Zhu, Fangqiang and Klaus Schulten} } @article {78, title = {Diffusion on curved, periodic surfaces}, journal = {Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics}, volume = {60}, year = {1999}, month = {1999 Jul}, pages = {302-7}, abstract = {

We present a simulation algorithm for a diffusion on a curved surface given by the equation phi(r)=0. The algorithm is tested against analytical results known for diffusion on a cylinder and a sphere, and applied to the diffusion on the P, D, and G periodic nodal surfaces. It should find application in an interpretation of two-dimensional exchange NMR spectroscopy data of diffusion on biological membranes.

}, issn = {1063-651X}, author = {Robert Ho{\l}yst and Plewczy{\'n}ski, D and Aleksei Aksimentiev and Burdzy, K} }