@article {397, title = {Dehydrated Biomimetic Membranes with Skinlike Structure and Function}, journal = {ACS Applied Materials \& Interfaces}, year = {2024}, abstract = {

Novel vapor-permeable materials are sought after for applications in protective wear, energy generation, and water treatment. Current impermeable protective materials effectively block harmful agents but trap heat due to poor water vapor transfer. Here we present a new class of materials, vapor permeable dehydrated nanoporous biomimetic membranes (DBMs), based on channel proteins. This application for biomimetic membranes is unexpected as channel proteins and biomimetic membranes were assumed to be unstable under dry conditions. DBMs mimic human skin{\textquoteright}s structure to offer both high vapor transport and small molecule exclusion under dry conditions. DBMs feature highly organized pores resembling sweat pores in human skin, but at super high densities (\>1012\ pores/cm2). These DBMs achieved exceptional water vapor transport rates, surpassing commercial breathable fabrics by up to 6.2 times, despite containing \>2 orders of magnitude smaller pores (1 nm vs \>700 nm). These DBMs effectively excluded model biological agents and harmful chemicals both in liquid and vapor phases, again in contrast with the commercial breathable fabrics. Remarkably, while hydrated biomimetic membranes were highly permeable to liquid water, they exhibited higher water resistances after dehydration at values \>38 times that of commercial breathable fabrics. Molecular dynamics simulations support our hypothesis that dehydration induced protein hydrophobicity increases which enhanced DBM performance. DBMs hold promise for various applications, including membrane distillation, dehumidification, and protective barriers for atmospheric water harvesting materials.

}, doi = {10.1021/acsami.3c19572}, url = {https://doi.org/10.1021/acsami.3c19572}, author = {Oh, Hyeonji and Tu, Yu-Ming and Samineni, Laximicharan and De Respino, Sophie and Mehrafrooz, Behzad and Joshi, Himanshu and Massenburg, Lynnicia and Lopez-Marques, Horacio and Elessawy, Nada and Song, Woochul and Behera, Harekrushna and Dhiman, Raman and Boorla, Veda Sheersh and Kher, Kartik and Lin, Yi-Chih and Maranas, Costas and Aksimentiev, Aleksei and D. Freeman, Benny and Kumar, Manish} } @article {394, title = {Molecular Determinants of Current Blockade Produced by Peptide Transport Through a Nanopore}, journal = {ACS Nanosci. Au}, year = {2023}, month = {11/2023}, abstract = {

The nanopore sensing method holds the promise of delivering a single molecule technology for identification of biological proteins, direct detection of post-translational modifications, and perhaps de novo determination of a protein{\textquoteright}s amino acid sequence. The key quantity measured in such nanopore sensing experiments is the magnitude of the ionic current passing through a nanopore blocked by a polypeptide chain. Establishing a relationship between the amino acid sequence of a peptide fragment confined within a nanopore and the blockade current flowing through the nanopore remains a major challenge for realizing the nanopore protein sequencing. Using the results of all-atom molecular dynamics simulations, here we compare nanopore sequencing of DNA with nanopore sequencing of proteins. We then delineate the factors affecting the blockade current modulation by the peptide sequence, showing that the current can be determined by (i) the steric footprint of an amino acid, (ii) its interactions with the pore wall, (iii) the local stretching of a polypeptide chain, and (iv) the local enhancement of the ion concentration at the nanopore constriction. We conclude with a brief discussion of the prospects for purely computational prediction of the blockade currents.

}, doi = {10.1021/acsnanoscienceau.3c00046}, author = {Jingqian Liu and Aleksei Aksimentiev} } @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 {393, title = {Long-Range Conductivity in Proteins Mediated by Aromatic Residues}, journal = {ACS Physical Chemistry Au}, volume = {3}, year = {2023}, month = {09}, pages = {444{\textendash}455}, abstract = {

Single-molecule measurements show that many proteins, lacking any redox cofactors, nonetheless exhibit electrical conductance on the order of a nanosiemen over 10 nm distances, implying that electrons can transit an entire protein in less than a nanosecond when subject to a potential difference of less than 1 V. This is puzzling because, for fast transport (i.e., a free energy barrier of zero), the hopping rate is determined by the reorganization energy of approximately 0.8 eV, and this sets the time scale of a single hop to at least 1 μs. Furthermore, the Fermi energies of typical metal electrodes are far removed from the energies required for sequential oxidation and reduction of the aromatic residues of the protein, which should further reduce the hopping current. Here, we combine all-atom molecular dynamics (MD) simulations of non-redox-active proteins (consensus tetratricopeptide repeats) with an electron transfer theory to demonstrate a molecular mechanism that can account for the unexpectedly fast electron transport. According to our MD simulations, the reorganization energy produced by the energy shift on charging (the Stokes shift) is close to the conventional value of 0.8 eV. However, the non-ergodic sampling of molecular configurations by the protein results in reaction-reorganization energies, extracted directly from the distribution of the electrostatic energy fluctuations, that are only \~{}0.2 eV, which is small enough to enable long-range conductivity, without invoking quantum coherent transport. Using the MD values of the reorganization energies, we calculate a current decay with distance that is in agreement with experiment.

}, doi = {10.1021/acsphyschemau.3c00017}, url = {https://doi.org/10.1021/acsphyschemau.3c00017}, author = {Krishnan, Siddharth and Aksimentiev, Aleksei and Lindsay, Stuart and Matyushov, Dmitry} } @article {384, title = {Unidirectional single-file transport of full-length proteins through a nanopore}, journal = {Nature Biotechnology}, year = {2023}, month = {01/2023}, abstract = {

The electrical current blockade of a peptide or protein threading through a nanopore can be used as a fingerprint of the molecule in biosensor applications. However, threading of full-length proteins has only been achieved using enzymatic unfolding and translocation. Here we describe an enzyme-free approach for unidirectional, slow transport of full-length proteins through nanopores. We show that the combination of a chemically resistant biological nanopore, α-hemolysin (narrowest part is \textasciitilde1.4\thinspacenm in diameter), and a high concentration guanidinium chloride buffer enables unidirectional, single-file protein transport propelled by an electroosmotic effect. We show that the mean protein translocation velocity depends linearly on the applied voltage and translocation times depend linearly on length, resembling the translocation dynamics of ssDNA. Using a supervised machine-learning classifier, we demonstrate that single-translocation events contain sufficient information to distinguish their threading orientation and identity with accuracies larger than 90\%. Capture rates of protein are increased substantially when either a genetically encoded charged peptide tail or a DNA tag is added to a protein.

}, issn = {1546-1696}, doi = {10.1038/s41587-022-01598-3}, url = {https://doi.org/10.1038/s41587-022-01598-3}, author = {Yu, Luning and Kang, Xinqi and Li, Fanjun and Mehrafrooz, Behzad and Makhamreh, Amr and Fallahi, Ali and Foster, Joshua C. and Aksimentiev, Aleksei and Chen, Min and Wanunu, Meni} } @article {392, title = {Sulfur-Containing Foldamer-Based Artificial Lithium Channels}, journal = {Angewandte Chemie International Edition}, volume = {62}, year = {2023}, pages = {e202305623}, abstract = {

Abstract Unlike many other biologically relevant ions (Na+, K+, Ca2+, Cl-, etc) and protons, whose cellular concentrations are closely regulated by highly selective channel proteins, Li+ ion is unusual in that its concentration is well tolerated over many orders of magnitude and that no lithium-specific channel proteins have so far been identified. While one naturally evolved primary pathway for Li+ ions to traverse across the cell membrane is through sodium channels by competing with Na+ ions, highly sought-after artificial lithium-transporting channels remain a major challenge to develop. Here we show that sulfur-containing organic nanotubes derived from intramolecularly H-bonded helically folded aromatic foldamers of 3.6\ {\r A}\ in hollow cavity diameter could facilitate highly selective and efficient transmembrane transport of Li+ ions, with high transport selectivity factors of 15.3 and 19.9 over Na+ and K+ ions, respectively.

}, keywords = {Foldamers, H-Bonds, Lithium Channels, Organic Nanotubes, Supramolecular Chemistry}, doi = {https://doi.org/10.1002/anie.202305623}, url = {https://onlinelibrary.wiley.com/doi/abs/10.1002/anie.202305623}, author = {Shen, Jie and R, Deepa and Li, Zhongyan and Oh, Hyeonji and Behera, Harekrushna and Joshi, Himanshu and Kumar, Manish and Aksimentiev, Aleksei and Zeng, Huaqiang} } @article {380, title = {Leakless end-to-end transport of small molecules through micron-length DNA nanochannels}, journal = {Science Advances}, volume = {8}, year = {2022}, month = {09/09/2022}, abstract = {
Designed and engineered protein and DNA nanopores can be used to sense and characterize single molecules and control transmembrane transport of molecular species. However, designed biomolecular pores are less than 100 nm in length and are used primarily for transport across lipid membranes. Nanochannels that span longer distances could be used as conduits for molecules between nonadjacent compartments or cells. Here, we design micrometer-long, 7-nm-diameter DNA nanochannels that small molecules can traverse according to the laws of continuum diffusion. Binding DNA origami caps to channel ends eliminates transport and demonstrates that molecules diffuse from one channel end to the other rather than permeating through channel walls. These micrometer-length nanochannels can also grow, form interconnects, and interface with living cells. This work thus shows how to construct multifunctional, dynamic agents that control molecular transport, opening ways of studying intercellular signaling and modulating molecular transport between synthetic and living cells.
}, doi = {10.1126/sciadv.abq4834}, url = {https://www.science.org/doi/10.1126/sciadv.abq4834}, author = {Li, Yi and Maffeo, Christopher and Joshi, Himanshu and Aleksei Aksimentiev and M{\'e}nard, Brice and Schulman, Rebecca} } @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 {365, title = {Multiple rereads of single proteins at single{\textendash}amino acid resolution using nanopores}, journal = {Science}, volume = {374}, year = {2021}, month = {11/2021}, pages = {1509{\textendash}1513}, abstract = {

A proteomics tool capable of identifying single proteins would be important for cell biology research and applications. Here, we demonstrate a nanopore-based single-molecule peptide reader sensitive to single{\textendash}amino acid substitutions within individual peptides. A DNA-peptide conjugate was pulled through the biological nanopore MspA by the DNA helicase Hel308. Reading the ion current signal through the nanopore enabled discrimination of single{\textendash}amino acid substitutions in single reads. Molecular dynamics simulations showed these signals to result from size exclusion and pore binding. We also demonstrate the capability to {\textquotedblleft}rewind{\textquotedblright} peptide reads, obtaining numerous independent reads of the same molecule, yielding an error rate of \<10-6\ in single amino acid variant identification. These proof-of-concept experiments constitute a promising basis for the development of a single-molecule protein fingerprinting and analysis technology.

}, doi = {10.1021/acs.nanolett.1c04203}, author = {Brinkerhoff, Henry and Kang, Albert SW and Liu, Jingqian and Aleksei Aksimentiev and Dekker, Cees} } @article {351, title = {Foldamer-Based Ultrapermeable and Highly Selective Artificial Water Channels that Exclude Protons}, journal = {Nature Nanotechnology}, volume = {16}, year = {2021}, month = {08/2021}, pages = {911-917}, chapter = {911}, abstract = {

The outstanding capacity of aquaporins (AQPs) for mediating highly selective superfast water transport has inspired the recent development of supramolecular monovalent ion-excluding artificial water channels (AWCs). AWC-based bioinspired membranes are proposed for desalination, water purification, and other separation applications. While some recent progress has been made in synthesizing AWCs that approach the water permeability and ion selectivity of AQPs, a hallmark feature of AQPs{\textemdash}high water transport while excluding protons-has not been reproduced. We report a class of biomimetic, helically folded pore-forming polymeric foldamers that can serve as long-sought-after highly selective ultrafast water-conducting channels with performance exceeding those of AQPs (1.1 {\texttimes} 1010 water molecules per second for AQP1, with high water-over- monovalent-ion transport selectivity (~108 water molecules over Cl-\ ion) conferred by the modularly tunable hydrophobicity of the interior pore surface. The best-performing AWC reported here delivers water transport at an exceptionally high rate, namely, 2.5 times that of AQP1, while concurrently rejecting salts (NaCl and KCl) and even protons.

}, doi = {10.1038/s41565-021-00915-2}, author = {Roy, Arundhati and Shen, Jie and Joshi, Himanshu and Song, Woochul and Tu, Yu-Ming and Chowdhury, Ratul and Ruijuan, Ye and Li, Ning and Ren, Changliang and Kumar, Manish and Aleksei Aksimentiev and Zeng, Huaqiang} } @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 {355, title = {The Manipulation of the Internal Hydrophobicity of FraC Nanopores Augments Peptide Capture and Recognition}, journal = {ACS Nano}, year = {2021}, month = {06/2021}, abstract = {

The detection of analytes and the sequencing of DNA using biological nanopores have seen major advances over recent years. The analysis of proteins and peptides with nanopores, however, is complicated by the complex physicochemical structure of polypeptides and the lack of understanding of the mechanism of capture and recognition of polypeptides by nanopores. In this work, we show that introducing aromatic amino acids at precise positions within the lumen of α-helical fragaceatoxin C (FraC) nanopores increased the capture frequency of peptides and largely improved the discrimination among peptides of similar size. Molecular dynamics simulations determined the sensing region of the nanopore, elucidated the microscopic mechanism enabling accurate characterization of the peptides\ via\ ionic current blockades in FraC, and characterized the effect of the pore modification on peptide discrimination. This work provides insights to improve the recognition and to augment the capture of peptides by nanopores, which is important for developing a real-time and single-molecule size analyzer for peptide recognition and identification.

}, doi = {10.1021/acsnano.0c09958}, url = {https://doi.org/10.1021/acsnano.0c09958}, author = {Florian Leonardus Rudolfus Lucas and Kumar Sarthak and Erica Mariska Lenting and David Coltan and Nieck Jordy van der Heide and Roderick Corstiaan Abraham Versloot and Aleksei Aksimentiev and Giovanni Maglia} } @article {357, title = {A tetrahedral DNA nanorobot with conformational change in response to molecular trigger}, journal = {Nanoscale}, year = {2021}, month = {01/2021}, abstract = {

Dynamic DNA origami nanostructures that respond to external stimuli are promising platforms for cargo delivery and nanoscale sensing. However, the low stability of such nanostructures under physiological conditions presents a major obstacle for their use in biomedical applications. This article describes a stable tetrahedral DNA nanorobot (TDN) programmed to undergo a controlled conformational change in response to epithelial cell adhesion molecule (EpCAM), a molecular biomarker specifically expressed on the circulating tumor cells. Multiresolution molecular dynamics simulations verified the overall stability of the folded TDN design and characterized local distortions in the folded structure. Atomic force microscopy and gel electrophoresis results showed that tetragonal structures are more stable than unfolded DNA origami sheets. Live cell experiments demonstrated the low cytotoxicity and target specificity of TDN. In summary, the proposed TDN can not only effectively resist nuclease catalysis but also has the potential to monitor EpCAM-positive cells precisely.

}, issn = {2040-3364}, doi = {10.1039/D1NR02757C}, url = {http://xlink.rsc.org/?DOI=D1NR02757Chttp://pubs.rsc.org/en/content/articlepdf/2021/NR/D1NR02757Chttp://pubs.rsc.org/en/content/articlepdf/2021/NR/D1NR02757C}, author = {Liu, Fengyu and Liu, Xiaoming and Shi, Qing and Maffeo, Christopher and Kojima, Masaru and Dong, Lixin and Aleksei Aksimentiev and Huang, Qiang and Fukuda, Toshio and Arai, Tatsuo} } @article {340, title = {Chiral Systems Made from DNA}, journal = {Advanced Science}, volume = {8}, year = {2021}, month = {01/2021}, pages = {2003113}, type = {Progress Report}, chapter = {2003113}, abstract = {

The very chemical structure of DNA that enables biological heredity and evolution has non-trivial implications for the self-organization of DNA molecules into larger assemblies and provides limitless opportunities for building functional nanostructures. This progress report discusses the natural organization of DNA into chiral structures and recent advances in creating synthetic chiral systems using DNA as a building material. How nucleic acid chirality naturally comes into play in a diverse array of situations is considered first, at length scales ranging from an individual nucleotide to entire chromosomes. Thereafter, chiral liquid crystal phases formed by dense DNA mixtures are discussed, including the ongoing efforts to understand their origins. The report then summarizes recent efforts directed toward building chiral structures, and other structures of complex topology, using the principle of DNA self-assembly. Discussed last are existing and proposed functional man-made nanostructures designed to either probe or harness DNA{\textquoteright}s chirality, from plasmonics and spintronics to\ biosensing.

}, doi = {10.1002/advs.202003113}, url = {https://doi.org/10.1002/advs.202003113}, author = {D Winogradoff and Pin-Yi Li and H Joshi and L Quednau and C Maffeo and Aleksei Aksimentiev} } @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 {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 {292, title = {Sequence-Dependent DNA Condensation as a Driving Force of DNA Phase Separation}, journal = {Nucleic Acids Research}, volume = {46}, year = {2018}, month = {10/2018}, pages = {9401--9413}, abstract = {

The physical properties of DNA have been suggested to play a central role in spatio-temporal organization of eukaryotic chromosomes. Experimental correlations have been established between the local nucleotide content of DNA and the frequency of inter- and intra-chromosomal contacts but the underlying physical mechanism remains unknown. Here, we combine fluorescence resonance energy transfer (FRET) measurements, precipitation assays, and molecular dynamics simulations to characterize the effect of DNA nucleotide content, sequence, and methylation on inter-DNA association and its correlation with DNA looping. First, we show that the strength of DNA condensation mediated by poly-lysine peptides as a reduced model of histone tails depends on the DNA{\textquoteright}s global nucleotide content but also on the local nucleotide sequence, which turns out to be qualitatively same as the condensation by spermine. Next, we show that the presence and spatial arrangement of C5 methyl groups determines the strength of inter-DNA attraction, partially explaining why RNA resists condensation. Interestingly, multi-color single molecule FRET measurements reveal strong anti-correlation between DNA looping and DNA{\textendash}DNA association, suggesting that a common biophysical mechanism underlies them. We propose that the differential affinity between DNA regions of varying sequence pattern may drive the phase separation of chromatin into chromosomal subdomains.

}, doi = {10.1093/nar/gky639 }, author = {Hyunju Kang and Jejoong Yoo and Hong Soo Lee and Byeong-Kwon Sohn and Seung-Won Lee and Wenjie Ma and Jung-Min Kee and Aleksei Aksimentiev and Hajin Kim} } @article {291, title = {Dynamics of a Molecular Plug Docked onto a Solid-State Nanopore}, journal = {The Journal of Physical Chemistry Letters}, volume = {9}, year = {2018}, month = {08/2018}, pages = {4686-4694}, chapter = {4686}, abstract = {

Docking of a protein- DNA complex onto a nanopore can provide ample observation time, and has enabled collection of analytic applications of biological nanopores, including DNA sequencing. However, the application of the same principle to solid-state nanopores is tempered by poor understanding of the docking process. Here, we elucidate the behavior of individual protein- DNA complexes docked onto a solid-state nanopore by monitoring the nanopore ionic current. Repeat docking of monovalent streptavidin- DNA complexes is found to produce ionic current blockades that fluctuate between discrete levels. We elucidate the roles of the protein plug and the DNA tether in the docking process, finding the docking configurations to determine the multitude of the current blockade levels, whereas the frequency of the current level switching is determined by the interactions between the molecules and the solid-state membrane. Finally, we prove the feasibility of using the nanopore docking principle for single-molecule sensing using solid-state nanopores by detecting conformational changes of a tethered DNA molecule from a random coil to an i-motif state.

}, doi = {10.1021/acs.jpclett.8b01755}, author = {Xin Shi and Qiao Li and Rui Gao and Wei Si and Shao-Chuang Liu and Aleksei Aksimentiev and Yi-Tao Long} } @article {286, title = {A synthetic enzyme built from DNA flips 10$^{7}$ lipids per second in biological membranes}, journal = {Nature Communications}, volume = {9}, year = {2018}, month = {07/2018}, pages = {2426}, abstract = {

Mimicking enzyme function and increasing performance of naturally evolved proteins is one of the most challenging and intriguing aims of nanoscience. Here, we employ DNA nanotechnology to design a synthetic enzyme that substantially outperforms its biological archetypes. Consisting of only eight strands, our DNA nanostructure spontaneously inserts into biological membranes by forming a toroidal pore that connects the membrane{\textquoteright}s inner and outer leaflets. The membrane insertion catalyzes spontaneous transport of lipid molecules between the bilayer leaflets, rapidly equilibrating the lipid composition. Through a combination of microscopic simulations and fluorescence microscopy we find the lipid transport rate catalyzed by the DNA nanostructure exceeds 107 molecules per second, which is three orders of magnitude higher than the rate of lipid transport catalyzed by biological enzymes. Furthermore, we show that our DNA-based enzyme can control the composition of human cell membranes, which opens new avenues for applications of membrane-interacting DNA systems in medicine.

}, keywords = {Computational biophysics, DNA nanostructures, Membrane biophysics, Nanobiotechnology}, doi = {10.1038/s41467-018-04821-5}, url = {http://www.nature.com/articles/s41467-018-04821-5}, author = {Ohmann, Alexander and Li, Chen-Yu and Maffeo, Christopher and Al Nahas, Kareem and Baumann, Kevin N. and G{\"o}pfrich, Kerstin and Yoo, Jejoong and Keyser, Ulrich F. and Aleksei Aksimentiev} } @article {285, title = {Inchworm movement of two rings switching onto a thread by biased Brownian diffusion represent a three-body problem}, journal = {Proceedings of the National Academy of Sciences}, year = {2018}, month = {05/2018}, abstract = {

This article describes a molecular realization of the classical three-body problem, where the motion of three or more bodies is directed by a set of pairwise forces. Surprisingly, motion of the components of the three-body molecular systems is found to be highly choreographed by differences in strength of intercomponent interactions, promoting a rare inchworm-like loading of molecular rings onto a molecular thread. Our work demonstrates the utility of an integrative approach to design and develop functional molecular machines.The coordinated motion of many individual components underpins the operation of all machines. However, despite generations of experience in engineering, understanding the motion of three or more coupled components remains a challenge, known since the time of Newton as the {\textquotedblleft}three-body problem.{\textquotedblright} Here, we describe, quantify, and simulate a molecular three-body problem of threading two molecular rings onto a linear molecular thread. Specifically, we use voltage-triggered reduction of a tetrazine-based thread to capture two cyanostar macrocycles and form a [3]pseudorotaxane product. As a consequence of the noncovalent coupling between the cyanostar rings, we find the threading occurs by an unexpected and rare inchworm-like motion where one ring follows the other. The mechanism was derived from controls, analysis of cyclic voltammetry (CV) traces, and Brownian dynamics simulations. CVs from two noncovalently interacting rings match that of two covalently linked rings designed to thread via the inchworm pathway, and they deviate considerably from the CV of a macrocycle designed to thread via a stepwise pathway. Time-dependent electrochemistry provides estimates of rate constants for threading. Experimentally derived parameters (energy wells, barriers, diffusion coefficients) helped determine likely pathways of motion with rate-kinetics and Brownian dynamics simulations. Simulations verified intercomponent coupling could be separated into ring{\textendash}thread interactions for kinetics, and ring{\textendash}ring interactions for thermodynamics to reduce the three-body problem to a two-body one. Our findings provide a basis for high-throughput design of molecular machinery with multiple components undergoing coupled motion.

}, issn = {0027-8424}, doi = {10.1073/pnas.1719539115}, url = {http://www.pnas.org/content/early/2018/05/03/1719539115}, author = {Benson, Christopher R. and Maffeo, Christopher and Fatila, Elisabeth M. and Liu, Yun and Sheetz, Edward G. and Aleksei Aksimentiev and Singharoy, Abhishek and Flood, Amar H.} } @inbook {277, title = {A Practical Guide to Molecular Dynamics Simulations of DNA Origami Systems}, booktitle = {DNA Nanotechology: Methods and Protocols}, series = {Methods in Molecular Biology}, volume = {1811}, year = {2018}, pages = {209-229}, publisher = {Humana Press}, organization = {Humana Press}, edition = {2}, chapter = {A Practical Guide to Molecular Dynamics Simulations of DNA Origami Systems}, address = {New York}, abstract = {

The DNA origami method exploits the self-assembly property of nucleic acids to build diverse nanoscale systems. The all-atom molecular dynamics (MD) method has emerged as a powerful computational tool for atomic-resolution characterization of the in situ structure and physical properties of DNA origami objects. This chapter provides step-by-step instructions for building atomic-scale models of DNA origami systems, using the MD method to simulate the models, and performing basic analyses of the resulting MD trajectories.\ 

}, issn = {978-1-4939-8581-4}, doi = {10.1007/978-1-4939-8582-1_15}, url = {https://link.springer.com/protocol/10.1007/978-1-4939-8582-1_15}, author = {Jejoong Yoo and Chen-Yu Li and Scott Michael Slone and Christopher Maffeo and Aleksei Aksimentiev}, editor = {Giampaolo Zuccheri} } @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 {237, title = {Large-Conductance Transmembrane Porin Made from DNA Origami}, journal = {ACS Nano}, volume = {10}, year = {2016}, month = {08/2016}, pages = {8207--8214}, abstract = {

DNA nanotechnology allows for the creation of three-dimensional structures at nanometer scale. Here, we use DNA to build the largest synthetic pore in a lipid membrane to date, approaching the dimensions of the nuclear pore complex and increasing the pore-area and the conductance tenfold compared to previous man-made channels. In our design, nineteen cholesterol-tags anchor a megadalton funnel-shaped DNA origami porin in a lipid bilayer membrane (see equilibration trajectory, require login to nanoHUB). Confocal imaging and ionic current recordings reveal spontaneous insertion of the DNA porin into the lipid membrane, creating a transmembrane pore of tens of nanosiemens conductance (see ionic current trajectory). All-atom molecular dynamics simulations characterize the conductance mechanism at the atomic level and independently confirm the DNA porins{\textquoteright} large ionic conductance.

}, doi = {10.1021/acsnano.6b03759}, author = {Kerstin G{\"o}pfrich* and Chen-Yu Li* and Maria Ricci and Satya Prathyusha Bhamidimarri and Jejoong Yoo and Bertalan Gyenes and Alexander Ohmann and Mathias Winterhalter and Aleksei Aksimentiev and Ulrich F Keyser} } @article {236, title = {Ion Channels Made from a Single Membrane-Spanning DNA Duplex}, journal = {Nano Letters}, volume = {16}, year = {2016}, month = {06/2016}, pages = {4665--4669}, abstract = {

Due to their hollow interior, transmembrane channels are capable of opening up pathways for ions across lipid membranes of living cells. Here, we demonstrate ion conduction induced by a single DNA duplex that lacks a hollow central channel. Decorated with six porpyrin-tags, our duplex is designed to span lipid membranes. Combining electrophysiology measurements with all-atom molecular dynamics simulations, we elucidate the microscopic conductance pathway\ (see\ water channel\ trajectory, require login to nanoHUB). Ions flow (see ion transport\ trajectory)\ at the DNA-lipid interface as the lipid head groups tilt towards the amphiphilic duplex forming a toroidal pore filled with water and ions. Ionic current traces produced by the DNA-lipid channel show well-defined insertion steps, closures and gating similar to those observed for traditional protein channels or synthetic pores. Ionic conductances obtained through simulations and experiments are in excellent quantitative agreement. The conductance mechanism realized here, with the smallest possible DNA-based ion channel, offers a route to design a new class of synthetic ion channels with maximum simplicity.

}, doi = {10.1021/acs.nanolett.6b02039}, author = {Kerstin G{\"o}pfrich and Chen-Yu Li and Iwona Mames and Satya Prathyusha Bhamidimarri and Maria Ricci and Jejoong Yoo and Adam Mames and Alexander Ohmann and Mathias Winterhalter and Eugen Stulz and Aleksei Aksimentiev and Ulrich F Keyser} } @article {226, title = {Molecular mechanics of DNA bricks: in situ structure, mechanical properties and ionic conductivity}, journal = {New Journal of Physics}, volume = {18}, year = {2016}, month = {05/2016}, pages = {055012}, abstract = {

The DNA bricks method exploits self-assembly of short DNA fragments to produce custom three-dimensional objects with subnanometer precision. In contrast to DNA origami, the DNA brick method permits a variety of different structures to be realized using the same library of DNA strands. As a consequence of their design, however, assembled DNA brick structures have fewer interhelical connections in comparison to equivalent DNA origami structures. Although the overall shape of the DNA brick objects has been characterized and found to conform to the features of the target designs, the microscopic properties of DNA brick objects remain yet to be determined. Here, we use the all-atom molecular dynamics method to directly compare the structure, mechanical properties and ionic conductivity of DNA brick and DNA origami structures different only by internal connectivity of their consistituent DNA strands. In comparison to equivalent DNA origami structures, the DNA brick structures are found to be less rigid and less dense and have a larger cross-section area normal to the DNA helix direction. At the microscopic level, the junction in the DNA brick structures are found to be right-handed, similar to the structure of individual Holliday junctions (HJ) in solution, which contrasts with the left-handed structure of HJ in DNA origami. Subject to external electric field, a DNA brick plate is more leaky to ions than an equivalent DNA origami plate because of its lower density and larger cross-section area. Overall, our results indicate that the structures produced by the DNA brick method are fairly similar in their overall appearance to those created by the DNA origami method but are more compliant when subject to external forces, which likely is a consequence of their single crossover design.

}, doi = {10.1088/1367-2630/18/5/055012}, url = {http://dx.doi.org/10.1088/1367-2630/18/5/055012}, author = {Scott Michael Slone and Chen-Yu Li and Jejoong Yoo and Aleksei Aksimentiev} } @article {216, title = {The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study}, journal = {PLoS ONE}, volume = {11}, year = {2016}, month = {02/2016}, pages = {e0152002}, abstract = {

LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50\% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.

}, doi = {10.1371/journal.pone.0152002}, url = {http://journals.plos.org/plosone/article?id=10.1371\%2Fjournal.pone.0152002}, author = {Lili Sun and Marcin Tabaka and Sen Hou and Lin Li and Krzysztof Burdzy and Aleksei Aksimentiev and Christopher Maffeo and Xuzhu Zhang and Robert Holyst} } @article {LI2015, title = {Ionic Conductivity, Structural Deformation, and Programmable Anisotropy of DNA Origami in Electric Field}, journal = {ACS Nano}, volume = {9}, year = {2015}, note = {PMID: 25623807}, month = {01/2015}, pages = {1420-1433}, abstract = {

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 Mg2+\ 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.

}, doi = {10.1021/nn505825z}, url = {http://dx.doi.org/10.1021/nn505825z}, author = {Li, Chen-Yu and Hemmig, Elisa A. and Kong, Jinglin and Yoo, Jejoong and Hern{\'a}ndez-Ainsa, Silvia and Keyser, Ulrich F. and Aleksei Aksimentiev} } @article {139, title = {Close encounters with DNA}, journal = {J Phys Condens Matter}, volume = {26}, year = {2014}, month = {10/2014}, pages = {413101}, abstract = {

Over the past ten years, the all-atom molecular dynamics method has grown in the scale of both systems and processes amenable to it and in its ability to make quantitative predictions about the behavior of experimental systems. The field of computational DNA research is no exception, witnessing a dramatic increase in the size of systems simulated with atomic resolution, the duration of individual simulations and the realism of the simulation outcomes. In this topical review, we describe the hallmark physical properties of DNA from the perspective of all-atom simulations. We demonstrate the amazing ability of such simulations to reveal the microscopic physical origins of experimentally observed phenomena. We also discuss the frustrating limitations associated with imperfections of present atomic force fields and inadequate sampling. The review is focused on the following four physical properties of DNA: effective electric charge, response to an external mechanical force, interaction with other DNA molecules and behavior in an external electric field.

}, keywords = {Review}, issn = {1361-648X}, doi = {10.1088/0953-8984/26/41/413101}, author = {Christopher Maffeo and Yoo, J and Jeffrey Comer and Wells, D B and Luan, B and Aleksei Aksimentiev} } @article {15, title = {End-to-end attraction of duplex DNA}, journal = {Nucleic Acids Res}, volume = {40}, year = {2012}, month = {2012 May}, pages = {3812-21}, abstract = {

Recent experiments [Nakata, M. et al., End-to-end stacking and liquid crystal condensation of 6 to 20 basepair DNA duplexes. Science 2007; 318:1276-1279] have demonstrated spontaneous end-to-end association of short duplex DNA fragments into long rod-like structures. By means of extensive all-atom molecular dynamic simulations, we characterized end-to-end interactions of duplex DNA, quantitatively describing the forces, free energy and kinetics of the end-to-end association process. We found short DNA duplexes to spontaneously aggregate end-to-end when axially aligned in a small volume of monovalent electrolyte. It was observed that electrostatic repulsion of 5{\textquoteright}-phosphoryl groups promoted the formation of aggregates in a conformation similar to the B-form DNA double helix. Application of an external force revealed that rupture of the end-to-end assembly occurs by the shearing of the terminal base pairs. The standard binding free energy and the kinetic rates of end-to-end association and dissociation processes were estimated using two complementary methods: umbrella sampling simulations of two DNA fragments and direct observation of the aggregation process in a system containing 458 DNA fragments. We found the end-to-end force to be short range, attractive, hydrophobic and only weakly dependent on the ion concentration. The relation between the stacking free energy and end-to-end attraction is discussed as well as possible roles of the end-to-end interaction in biological and nanotechnological systems.

}, keywords = {DNA, Kinetics, Molecular Dynamics Simulation, Nucleic Acid Conformation}, issn = {1362-4962}, doi = {10.1093/nar/gkr1220}, author = {Christopher Maffeo and Luan, Binquan and Aleksei Aksimentiev} } @inbook {192, title = {Third Generation DNA Sequencing with a Nanopore}, booktitle = {Nanopores: Sensing fundamental biological interactions at the single molecule level}, year = {2011}, pages = {287-311}, publisher = {Springer US}, organization = {Springer US}, chapter = {12}, abstract = {

With the advent of Next-Generation-Sequencing (NGS) technologies, an enormous volume of DNA sequencing data can be generated at low cost, placing genomic science within the grasp of everyday medicine. However, mired in this voluminous data, a new problem has emerged: the assembly of the genome from the short reads. In this chapter we examine the prospects for sequencing DNA using a synthetic nanopore. Nanopore sequencing has the potential for very long reads, reducing the computational burden posed by alignment and genome assembly, while at the same time eliminating logistically challenging and error-prone amplification and library formation due to its exquisite single molecule sensitivity. On the other hand, long high fidelity reads demand stringent control over both the DNA configuration in the pore and the translocation kinetics. We examine the prospects for satisfying these specifications with a synthetic nanopore.

}, keywords = {Comparison between biological and synthetic nanopores, Hydrodynamic focusing, Molecular dynamics simulations, Sequencing DNA, Trapping DNA, λ-DNA}, isbn = {978-1-4419-8251-3}, doi = {10.1007/978-1-4419-8252-0_12}, url = {http://dx.doi.org/10.1007/978-1-4419-8252-0_12}, author = {Gregory Timp and Mirsaidov, Utkur and Timp, Winston and Shim, Jiwook and Wang, Deqiang and Valentin Dimitrov and Scrimgeour, Jan and Lin, Chunchen and Jeffrey Comer and Anthony H. Ho and Zou, Xueqing and Aleksei Aksimentiev and Klaus Schulten}, editor = {Iqbal, Samir M. and Rashid Bashir} } @article {38, title = {Control and reversal of the electrophoretic force on DNA in a charged nanopore}, journal = {J Phys Condens Matter}, volume = {22}, year = {2010}, month = {2010 Nov 17}, pages = {454123}, abstract = {

Electric field driven transport of DNA through solid-state nanopores is the key process in nanopore-based DNA sequencing that promises dramatic reduction of genome sequencing costs. A major hurdle in the development of this sequencing method is that DNA transport through the nanopores occurs too quickly for the DNA sequence to be detected. By means of all-atom molecular dynamics simulations, we demonstrate that the velocity of DNA transport through a nanopore can be controlled by the charge state of the nanopore surface. In particular, we show that the charge density of the nanopore surface controls the magnitude and/or direction of the electro-osmotic flow through the nanopore and thereby can significantly reduce or even reverse the effective electrophoretic force on DNA. Our work suggests a physical mechanism to control DNA transport in a nanopore by chemical, electrical or electrochemical modification of the nanopore surface.

}, keywords = {Base Sequence, Computer Simulation, DNA, Electromagnetic Fields, Electrophoresis, Models, Chemical, Molecular Sequence Data, Nanostructures, Porosity, Sequence Analysis, DNA, Static Electricity, Stress, Mechanical}, issn = {1361-648X}, doi = {10.1088/0953-8984/22/45/454123}, author = {Luan, Binquan and Aleksei Aksimentiev} } @article {44, title = {Electric and Electrophoretic Inversion of the DNA Charge in Multivalent Electrolytes}, journal = {Soft Matter}, volume = {6}, year = {2010}, month = {2010 Jan 1}, pages = {243-246}, abstract = {

Counterion-induced inversion of the DNA charge was characterized through extensive molecular dynamics simulations. We observed reversal of the DNA motion in an external electric field upon increasing the concentration of trivalent or quadrivalent counterions. In the case of a divalent electrolyte, inversion of the DNA{\textquoteright}s electric charge was observed at high concentrations of the electrolyte but not reversal of the DNA{\textquoteright} electrophoretic motion. We demonstrate that inversion of the DNA{\textquoteright}s electrophoretic mobility results from a complex interplay of electrostatics and hydrodynamics.

}, issn = {1744-683X}, doi = {10.1039/b917973a}, author = {Luan, Binquan and Aleksei Aksimentiev} } @article {46, title = {The effect of calcium on the conformation of cobalamin transporter BtuB}, journal = {Proteins}, volume = {78}, year = {2010}, month = {2010 Apr}, pages = {1153-62}, abstract = {

BtuB is a beta-barrel membrane protein that facilitates transport of cobalamin (vitamin B12) from the extracellular medium across the outer membrane of Escherichia coli. It is thought that binding of B12 to BtuB alters the conformation of its periplasm-exposed N-terminal residues (the TonB box), which enables subsequent binding of a TonB protein and leads to eventual uptake of B12 into the cytoplasm. Structural studies determined the location of the B12 binding site at the top of the BtuB{\textquoteright}s beta-barrel, surrounded by extracellular loops. However, the structure of the loops was found to depend on the method used to obtain the protein crystals, which-among other factors-differed in calcium concentration. Experimentally, calcium concentration was found to modulate the binding of the B12 substrate to BtuB. In this study, we investigate the effect of calcium ions on the conformation of the extracellular loops of BtuB and their possible role in B12 binding. Using all-atom molecular dynamics, we simulate conformational fluctuations of several X-ray structures of BtuB in the presence and absence of calcium ions. These simulations demonstrate that calcium ions can stabilize the conformation of loops 3-4, 5-6, and 15-16, and thereby prevent occlusion of the binding site. Furthermore, binding of calcium ions to extracellular loops of BtuB was found to enhance correlated motions in the BtuB structure, which is expected to promote signal transduction. Finally, we characterize conformation dynamics of the TonB box in different X-ray structures and find an interesting correlation between the stability of the TonB box structure and calcium binding.

}, keywords = {Amino Acid Sequence, Bacterial Outer Membrane Proteins, Calcium, Escherichia coli Proteins, Lipid Bilayers, Membrane Transport Proteins, Models, Molecular, Molecular Dynamics Simulation, Molecular Sequence Data, Protein Binding, Protein Conformation, Vitamin B 12, Vitamin B Complex}, issn = {1097-0134}, doi = {10.1002/prot.22635}, author = {Luan, Binquan and Carr, Rogan and Caffrey, Martin and Aleksei Aksimentiev} } @article {55, title = {Strain softening in stretched DNA}, journal = {Phys Rev Lett}, volume = {101}, year = {2008}, month = {2008 Sep 12}, pages = {118101}, abstract = {

The microscopic mechanics of DNA stretching was characterized using extensive molecular dynamics simulations. By employing an anisotropic pressure-control method, realistic force-extension dependences of effectively infinite DNA molecules were obtained. A coexistence of B and S DNA domains was observed during the overstretching transition. The simulations revealed that strain softening may occur in the process of stretching torsionally constrained DNA. The latter observation was qualitatively reconciled with available experimental data using a random-field Ising model.

}, keywords = {Computer Simulation, DNA, Elasticity, Kinetics, Models, Molecular, Nucleic Acid Conformation, Thermodynamics}, issn = {0031-9007}, author = {Luan, Binquan and Aleksei Aksimentiev} } @article {53, title = {DNA attraction in monovalent and divalent electrolytes}, journal = {J Am Chem Soc}, volume = {130}, year = {2008}, month = {2008 Nov 26}, pages = {15754-5}, abstract = {

The dependence of the effective force on the distance between two DNA molecules was directly computed from a set of extensive all-atom molecular dynamics simulations. The simulations revealed that in a monovalent electrolyte the effective force is repulsive at short and long distances but can be attractive in the intermediate range. This attractive force is, however, too weak (approximately 5 pN per turn of a DNA helix) to induce DNA condensation in the presence of thermal fluctuations. In divalent electrolytes, DNA molecules were observed to form a bound state, where Mg(2+) ions bridged minor groves of DNA. The effective force in divalent electrolytes was predominantly attractive, reaching a maximum of 42 pN per one turn of a DNA helix.

}, keywords = {Computer Simulation, DNA, Electrolytes, Models, Molecular, Molecular Conformation}, issn = {1520-5126}, doi = {10.1021/ja804802u}, author = {Luan, Binquan and Aleksei Aksimentiev} } @article {56, title = {Electro-osmotic screening of the DNA charge in a nanopore}, journal = {Phys Rev E Stat Nonlin Soft Matter Phys}, volume = {78}, year = {2008}, month = {2008 Aug}, pages = {021912}, abstract = {

Extensive all-atom molecular dynamics simulations were performed to characterize the microscopic origin of the force experienced by DNA in a bulk electrolyte and a solid-state nanopore when subject to an external electrostatic field E. The effective screening of the DNA charge was found to originate from the hydrodynamic drag of the electro-osmotic flow that is driven by the motion of counterions along the surface of DNA. We show that the effective driving force F in a nanopore obeys the same law as in a bulk electrolyte: F=ximuE , where xi and mu are the friction coefficient and electrophoretic mobility of DNA, respectively. Using this relationship, we suggest a method for determining the effective driving force on DNA in a nanopore that does not require a direct force measurement.

}, keywords = {DNA, Models, Theoretical, Nanostructures, Osmosis, Static Electricity, Water}, issn = {1539-3755}, author = {Luan, Binquan and Aleksei Aksimentiev} } @article {61, title = {In meso crystal structure and docking simulations suggest an alternative proteoglycan binding site in the OpcA outer membrane adhesin}, journal = {Proteins}, volume = {71}, year = {2008}, month = {2008 Apr}, pages = {24-34}, abstract = {

OpcA is an integral outer membrane adhesin protein from Neisseria meningitidis, the causative agent of meningococcal meningitis and septicemia. It binds to sialic acid (SA)-containing polysaccharides on the surface of epithelial cells. The crystal structure of OpcA showed that the protein adopts a 10-stranded beta-barrel structure, with five extensive loop regions on the extracellular side of the membrane. These form a crevice structure, lined with basic residues, which was hypothesized to act as the binding site for polysaccharide ligands. In the current study, a distinctly different OpcA structure has been obtained using crystals grown from a lipidic mesophase. Comparison of the two structures shows that the largest loop (L2), which closes over the end of the beta-barrel in the original crystal form, adopts a much more extended structure by reaching outward and away from the protein. The difference in conformation may be attributable to the absence of zinc ions from the crystallization conditions for the in meso crystal form: in the original structure, two zinc ions were bound to the external loops. Molecular dynamics (MD) simulations performed on the two OpcA models in a lipid bilayer environment demonstrated pronounced loop mobility. These observations support the view that the loop regions of OpcA are capable of a high degree of conformational flexibility. The original binding site for polysaccharide is not present in the in meso crystal form, and is disrupted during MD simulations. Docking analysis suggests a putative alternative location for the SA ligand in the new crystal form of OpcA.

}, keywords = {Adhesins, Bacterial, Bacterial Outer Membrane Proteins, Binding Sites, Computer Simulation, Crystallography, X-Ray, Models, Molecular, Neisseria meningitidis, Protein Binding, Proteoglycans, Spectrum Analysis}, issn = {1097-0134}, doi = {10.1002/prot.21841}, author = {Cherezov, Vadim and Liu, Wei and Derrick, Jeremy P and Luan, Binquan and Aleksei Aksimentiev and Katritch, Vsevolod and Caffrey, Martin} } @article {63, title = {Structure refinement of the OpcA adhesin using molecular dynamics}, journal = {Biophys J}, volume = {93}, year = {2007}, month = {2007 Nov 1}, pages = {3058-69}, abstract = {

OpcA from Neisseria meningitidis, the causative agent of meningococcal meningitis and septicemia, is an integral outer membrane protein that facilitates meningococcal adhesion through binding the proteoglycan receptors of susceptible cells. Two structures of OpcA have been determined by x-ray diffraction to 2 A resolution, revealing dramatically different conformations in the extracellular loops--the protein domain implicated in proteoglycan binding. In the first structure, a positively charged crevice formed by loops 1 and 2 was identified as the site for binding proteoglycans, whereas in the second structure the crevice was not evident as loops 1 and 2 adopted different conformations. To reconcile these results, molecular-dynamics simulations were carried out on both structures embedded in a solvated lipid bilayer membrane. Free of crystal contacts and crystallization agents, the loops were observed to undergo large structural transformations, suggesting that the conformation of the loops in either x-ray structure is affected by crystallization. Subsequent simulations of both structures in their crystal lattices confirmed this conclusion. Based on our molecular-dynamics trajectories, we propose a model for OpcA that combines stable structural features of the available x-ray structures. In this model, all five extracellular loops of OpcA have stable secondary structures. The loops form a funnel that leads to the base of the beta-barrel and that includes Tyr-169 on its exposed surface, which has been implicated in proteoglycan binding.

}, keywords = {Adhesins, Bacterial, Bacterial Outer Membrane Proteins, Computer Simulation, Crystallography, X-Ray, Models, Molecular, Neisseria meningitidis, Protein Structure, Tertiary, Thermodynamics}, issn = {0006-3495}, doi = {10.1529/biophysj.107.106724}, author = {Luan, Binquan and Caffrey, Martin and Aleksei Aksimentiev} } @article {67, title = {The role of molecular modeling in bionanotechnology}, journal = {Phys Biol}, volume = {3}, year = {2006}, month = {2006 Mar}, pages = {S40-53}, abstract = {

Molecular modeling is advocated here as a key methodology for research and development in bionanotechnology. Molecular modeling provides nanoscale images at atomic and even electronic resolution, predicts the nanoscale interaction of unfamiliar combinations of biological and inorganic materials, and evaluates strategies for redesigning biopolymers for nanotechnological uses. The methodology is illustrated in this paper through reviewing three case studies. The first one involves the use of single-walled carbon nanotubes as biomedical sensors where a computationally efficient, yet accurate, description of the influence of biomolecules on nanotube electronic properties through nanotube-biomolecule interactions was developed; this development furnishes the ability to test nanotube electronic properties in realistic biological environments. The second case study involves the use of nanopores manufactured into electronic nanodevices based on silicon compounds for single molecule electrical recording, in particular, for DNA sequencing. Here, modeling combining classical molecular dynamics, material science and device physics, described the interaction of biopolymers, e.g., DNA, with silicon nitrate and silicon oxide pores, furnished accurate dynamic images of pore translocation processes, and predicted signals. The third case study involves the development of nanoscale lipid bilayers for the study of embedded membrane proteins and cholesterol. Molecular modeling tested scaffold proteins, redesigned apolipoproteins found in mammalian plasma that hold the discoidal membranes in the proper shape, and predicted the assembly as well as final structure of the nanodiscs. In entirely new technological areas such as bionanotechnology, qualitative concepts, pictures and suggestions are sorely needed; these three case studies document that molecular modeling can serve a critical role in this respect, even though it may still fall short on quantitative precision.

}, keywords = {Biosensing Techniques, Chemical Engineering, DNA, Lipid Bilayers, Membrane Potentials, Models, Molecular, Nanotechnology, Nanotubes, Carbon, Sequence Analysis, Silicon}, issn = {1478-3975}, doi = {10.1088/1478-3975/3/1/S05}, author = {Lu, Deyu and Aleksei Aksimentiev and Amy Y Shih and Eduardo Cruz-Ch{\'u} and Freddolino, Peter L and Arkhipov, Anton and Klaus Schulten} } @article {191, title = {Simulation of the electric response of DNA translocation through a semiconductor nanopore-capacitor}, journal = {Nanotechnology}, volume = {17}, year = {2006}, pages = {622{\textendash}633}, abstract = {

A multi-scale/multi-material computational model for simulation of the electric signal detected on the electrodes of a metal{\textendash}oxide{\textendash}semiconductor (MOS) capacitor forming a nanoscale artificial membrane, and containing a nanopore with translocating DNA, is presented. The multi-scale approach is based on the incorporation of a molecular dynamics description of a translocating DNA molecule in the nanopore within a three-dimensional Poisson equation self-consistent scheme involving electrolytic and semiconductor charges for the electrostatic potential calculation. The voltage signal obtained from the simulation supports the possibility for single nucleotide resolution with a nanopore device. The electric signal predicted on the capacitor electrodes complements ongoing experiments exploring the use of nanopores in a MOS capacitor membrane for DNA sequencing.

}, doi = {10.1088/0957-4484/17/3/002}, author = {Maria E. Gracheva and A. Xiong and Aleksei Aksimentiev and Klaus Schulten and Gregory Timp and Jean-Pierre Leburton} } @article {190, title = {Electrical signatures of single-stranded DNA with single base mutations in a nanopore capacitor}, volume = {17}, year = {2006}, pages = {3160{\textendash}3165}, abstract = {

In this paper, we evaluate the magnitude of the electrical signals produced by DNA translocation through a 1\ nm diameter nanopore in a capacitor membrane with a numerical multi-scale approach, and assess the possibility of resolving individual nucleotides as well as their types in the absence of conformational disorder. We show that the maximum recorded voltage caused by the DNA translocation is about 35\ mV, while the maximum voltage signal due to the DNA backbone is about 30\ mV, and the maximum voltage of a DNA base is about 8\ mV. Signals from individual nucleotides can be identified in the recorded voltage traces, suggesting a 1\ nm diameter pore in a capacitor can be used to accurately count the number of nucleotides in a DNA strand. Furthermore, we study the effect of a single base substitution on the voltage trace, and calculate the differences among the voltage traces due to a single base mutation for the sequences C3AC7, C3CC7, C3GC7 and C3TC7. The calculated voltage differences are in the 5{\textendash}10\ mV range. The calculated maximum voltage caused by the translocation of individual bases varies from 2 to 9\ mV, which is experimentally detectable.

}, doi = {10.1088/0957-4484/17/13/014}, author = {Maria E. Gracheva and Aleksei Aksimentiev and Jean-Pierre Leburton} }