@article {374, title = {Discrimination of RNA fiber structures using solid-state nanopores}, journal = {Nanoscale}, volume = {14}, year = {2022}, month = {05/2023}, pages = {6866-6875}, abstract = {

RNA fibers are a class of biomaterials that can be assembled using HIV-like kissing loop interactions. Because of the programmability of molecular design and low immunorecognition, these structures present an interesting opportunity to solve problems in nanobiotechnology and synthetic biology. However, the experimental tools to fully characterize and discriminate among different fiber structures in solution are limited. Herein, we utilize solid-state nanopore experiments and Brownian dynamics simulations to characterize and distinguish several RNA fiber structures that differ in their degrees of branching. We found that, regardless of the electrolyte type and concentration, fiber structures that have more branches produce longer and deeper ionic current blockades in comparison to the unbranched fibers. Experiments carried out at temperatures ranging from 20{\textendash}60 {\textdegree}C revealed almost identical distributions of current blockade amplitudes, suggesting that the kissing loop interactions in fibers are resistant to heating within this range.

}, issn = {2040-3364}, doi = {10.1039/D1NR08002D}, url = {http://pubs.rsc.org/en/content/articlepdf/2022/NR/D1NR08002D}, author = {Tripathi, Prabhat and Chandler, Morgan and Maffeo, Christopher Michael and Fallahi, Ali and Makhamreh, Amr and Halman, Justin and Aleksei Aksimentiev and Afonin, Kirill A. and Wanunu, Meni} } @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 {364, title = {A nanoscale reciprocating rotary mechanism with coordinated mobility control.}, journal = {Nature Communications}, volume = {12}, year = {2021}, month = {12/2021}, pages = {7138}, chapter = {7138}, abstract = {

Biological molecular motors transform chemical energy into mechanical work by coupling cyclic catalytic reactions to large-scale structural transitions. Mechanical deformation can be surprisingly efficient in realizing such coupling, as demonstrated by the FoF1 ATP synthase. Here, we describe a synthetic molecular mechanism that transforms a rotary motion of an asymmetric camshaft into reciprocating large-scale transitions in a surrounding stator orchestrated by mechanical deformation. We design the mechanism using DNA origami, characterize its structure via cryo-electron microscopy, and examine its dynamic behavior using single-particle fluorescence microscopy and molecular dynamics simulations. While the camshaft can rotate inside the stator by diffusion, the stator{\textquoteright}s mechanics makes the camshaft pause at preferred orientations. By changing the stator{\textquoteright}s mechanical stiffness, we accelerate or suppress the Brownian rotation, demonstrating an allosteric coupling between the camshaft and the stator. Our mechanism provides a framework for manufacturing artificial nanomachines that function because of coordinated movements of their components.

}, keywords = {Cryoelectron Microscopy, DNA, Molecular Conformation, Molecular Dynamics Simulation, Nanostructures, Nanotechnology, Rotation}, issn = {2041-1723}, doi = {10.1038/s41467-021-27230-7}, author = {Bertosin, Eva and Maffeo, Christopher M and Drexler, Thomas and Honemann, Maximilian N and Aleksei Aksimentiev and Dietz, Hendrik} } @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 {360, title = {Synthetic Macrocycle Nanopore for Potassium-Selective Transmembrane Transport}, journal = {Journal of the American Chemical Society}, volume = {143}, year = {2021}, month = {08/2021}, pages = {15975-15983}, chapter = {15975}, abstract = {

Reproducing the structure and function of biological membrane channels, synthetic nanopores have been developed for applications in membrane filtration technologies and biomolecular sensing. Stable stand-alone synthetic nanopores have been created from a variety of materials, including peptides, nucleic acids, synthetic polymers, and solid-state membranes. In contrast to biological nanopores, however, furnishing such synthetic nanopores with an atomically defined shape, including deliberate placement of each and every chemical group, remains a major challenge. Here, we introduce a chemosynthetic macromolecule,\ extended pillararene macrocycle (EPM) as a chemically defined transmembrane nanopore that exhibits selective transmembrane transport. Our ionic current measurements reveal stable insertion of individual EPM nanopores into a lipid bilayer membrane and remarkable cation type-selective transport, with up to a 21-fold selectivity for potassium over sodium ions. Taken together, direct chemical synthesis offers a path to de novo design of a new class of synthetic nanopores with custom transport functionality imprinted in their atomically defined chemical structure.\ 

}, keywords = {Electrical conductivity, MD Simulations, Nanopores, Potassium Ions Selectivity}, doi = {10.1021/jacs.1c04910}, author = {Qiao, Dan and Joshi, Himanshu and Zhu, Huangtianzhi and Wang, Fushi and Xu, Yang and Gao, Jia and Huang, Feihe and Aleksei Aksimentiev and Feng, Jiandong} } @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 {354, title = {Hydrophobic Interactions between DNA Duplexes and Synthetic and Biological Membranes.}, journal = {Journal of the American Chemical Society}, volume = {143}, year = {2021}, month = {05/2021}, pages = {8305-8313}, chapter = {8305}, abstract = {

Equipping DNA with hydrophobic anchors enables targeted interaction with lipid bilayers for applications in biophysics, cell biology, and synthetic biology. Understanding DNA-membrane interactions is crucial for rationally designing functional DNA. Here we study the interactions of hydrophobically tagged DNA with synthetic and cell membranes using a combination of experiments and atomistic molecular dynamics (MD) simulations. The DNA duplexes are rendered hydrophobic by conjugation to a terminal cholesterol anchor or by chemical synthesis of a charge-neutralized alkyl-phosphorothioate (PPT) belt. Cholesterol-DNA tethers to lipid vesicles of different lipid compositions and charges, while PPT DNA binding strongly depends on alkyl length, belt position, and headgroup charge. Divalent cations in the buffer can also influence binding. Our MD simulations directly reveal the complex structure and energetics of PPT DNA within a lipid membrane, demonstrating that longer alkyl-PPT chains provide the most stable membrane anchoring but may disrupt DNA base paring in solution. When tested on cells, cholesterol-DNA is homogeneously distributed on the cell surface, while alkyl-PPT DNA accumulates in clustered structures on the plasma membrane. DNA tethered to the outside of the cell membrane is distinguished from DNA spanning the membrane by nuclease and sphingomyelinase digestion assays. The gained fundamental insight on DNA-bilayer interactions will guide the rational design of membrane-targeting nanostructures.

}, issn = {1520-5126}, doi = {10.1021/jacs.0c13235}, author = {Jones, Sioned F and Joshi, Himanshu and Terry, Stephen J and Burns, Jonathan R and Aleksei Aksimentiev and Eggert, Ulrike S and Howorka, Stefan} } @article {343, title = {Determining the In-Plane Orientation and Binding Mode of Single Fluorescent Dyes in DNA Origami Structures}, journal = {ACS Nano}, volume = {15}, year = {2021}, month = {03/2021}, pages = {5109-5117}, chapter = {5109}, abstract = {

We present a technique to determine the orientation of single fluorophores attached to DNA origami structures based on two measurements. First, the orientation of the absorption transition dipole of the molecule is determined through a polarization-resolved excitation measurement. Second, the orientation of the DNA origami structure is obtained from a DNA-PAINT nanoscopy measurement. Both measurements are performed consecutively on a fluorescence wide-field microscope. We employed this approach to study the orientation of single ATTO 647N, ATTO 643, and Cy5 fluorophores covalently attached to a 2D rectangular DNA origami structure with different nanoenvironments, achieved by changing both the fluorophores{\textquoteright} binding position and immediate vicinity. Our results show that when fluorophores are incorporated with additional space, for example, by omitting nucleotides in an elsewise doublestranded environment, they tend to stick to the DNA and to adopt a preferred orientation that depends more on the specific molecular environment than on the fluorophore type. With the aid of all-atom molecular dynamics simulations, we rationalized our observations and provide insight into the fluorophores{\textquoteright} probable binding modes. We believe this work constitutes an important step toward manipulating the orientation of single fluorophores in DNA origami structures, which is vital for the development of more efficient and reproducible self-assembled nanophotonic devices.\ 

}, keywords = {DNA Nanotechnology, DNA-PAINT, Single-molecule fluorescence polarization-resolved microscopy, Super-resolution microscopy}, doi = {https://doi.org/10.1021/acsnano.0c10259}, url = {https://pubs.acs.org/doi/10.1021/acsnano.0c10259}, author = {Hübner, Kristina and Joshi, Himanshu and Aleksei Aksimentiev and Stefani, Fernando D and Tinnefeld, Philip and Acuna, Guillermo P} } @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 {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 {293, title = {Dynamic Interactions between Lipid-Tethered DNA and Phospholipid Membranes}, journal = {Langmuir}, volume = {34}, year = {2018}, month = {11/2018}, pages = {15084-15092}, chapter = {15084}, abstract = {

Lipid-anchored DNA can attach functional cargo to bilayer membranes in DNA nanotechnology, synthetic biology, and cell biology research. To optimize DNA anchoring, an understanding of DNA-membrane interactions in terms of binding strength, extent, and structural dynamics is required. Here we use experiments and molecular dynamics (MD) simulations to determine how the membrane binding of cholesterol-modified DNA depends on electrostatic and steric factors involving the lipid headgroup charge, duplexed or single-stranded DNA, and the buffer composition. The experiments distinguish between free and membrane vesicle-bound DNA and thereby reveal the surface density of anchored DNA and its binding affinity, something which had previously not been known. The Kd values range from 8.5 {\textpm} 4.9 to 466 {\textpm} 134 μM whereby negatively charged headgroups led to weak binding due to the electrostatic repulsion with respect to the negatively charged DNA. Atomistic MD simulations explain the findings and elucidate the dynamic nature of anchored DNA such as the mushroom-like conformation of single-stranded DNA hovering over the bilayer surface in contrast to a straight-up conformation of double-stranded DNA. The biophysical insight into the binding strength to membranes as well as the molecular accessibility of DNA for hybridization to molecular cargo is expected to facilitate the\  creation of biomimetic DNA versions of natural membrane nanopores and cytoskeletons for research and nanobiotechnology.

}, keywords = {DNA Nanotechnology, Lipid Bilayer anchored DNA, Lipid-DNA binding affinity}, doi = {10.1021/acs.langmuir.8b02271}, url = {https://doi.org/10.1021/acs.langmuir.8b02271}, author = {Arnott, Patrick M. and Joshi, Himanshu and Aleksei Aksimentiev and Howorka, Stefan} } @article {281, title = {Optical voltage sensing using DNA origami}, journal = {Nano Letters}, volume = {18}, year = {2018}, pages = {1962-1971}, abstract = {

We explore the potential of DNA nanotechnology for developing novel optical voltage sensing nano-devices that convert a local change of electric potential into optical signals. As a proof-of-concept of the sensing mechanism, we assembled voltage responsive DNA origami structures labelled with a single pair of FRET dyes. The DNA structures were reversibly immobilised on a nanocapillary tip and underwent controlled structural changes upon application of an electric field. The applied field was monitored through a change in FRET efficiency. By exchanging the position of a single dye, we could tune the voltage sensitivity of our DNA origami structure, demonstrating the flexibility and versatility of our approach. The experimental studies were complemented by coarse-grained simulations that characterised voltage-dependent elastic deformation of the DNA nanostructures and the associated change in the distance between the FRET pair. Our work opens a novel pathway for determining the mechanical properties of DNA origami structures, and highlights potential applications of dynamic DNA nanostructures as voltage sensors.

}, doi = {10.1021/acs.nanolett.7b05354}, url = {https://doi.org/10.1021/acs.nanolett.7b05354}, author = {Hemmig, Elisa Alina and Fitzgerald, Clare and Maffeo, Christopher and Hecker, Lisa and Ochmann, Sarah Elisabeth and Aleksei Aksimentiev and Tinnefeld, Philip and Keyser, Ulrich F} } @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 {258, title = {Picomolar Fingerprinting of Nucleic Acid Nanoparticles Using Solid-State Nanopores}, journal = {ACS Nano}, volume = {11}, year = {2017}, month = {07/2017}, pages = {9701-9710}, abstract = {

Nucleic acid nanoparticles (NANPs) are an emerging class of programmable structures with tunable shape and function. Their promise as tools for fundamental biophysics studies, molecular sensing, and therapeutic applications necessitates methods for their detection and characterization at the single-particle level. In this work, we study electrophoretic transport of individual ring-shaped and cube-shaped NANPs through solid-state nanopores. In the optimal nanopore size range, the particles must deform to pass through, which considerably increases their residence time within the pore. Such anomalously long residence times permit detection of picomolar amounts of NANPs when nanopore measurements are carried out at a high transmembrane bias. In the case of a NANP mixture, the type of individual particle passing through nanopores can be efficiently determined from analysis of a single electrical pulse. Molecular dynamics simulations provide insight into the mechanical barrier to transport of the NANPs and corroborate the difference in the signal amplitudes observed for the two types of particles. Our study serves as a basis for label-free analysis of soft programmable-shape nanoparticles.

}, doi = {10.1021/acsnano.7b04923}, url = {http://dx.doi.org/10.1021/acsnano.7b04923}, author = {Alibakhshi, Mohammad Amin and Halman, Justin R. and Wilson, James and Aleksei Aksimentiev and Afonin, Kirill A. and Wanunu, Meni} } @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 {235, title = {Graphene Nanopores for Protein Sequencing}, journal = {Advanced Functional Materials}, volume = {26}, year = {2016}, month = {06/2016}, pages = {4830-4838}, chapter = {4830}, abstract = {

An inexpensive, reliable method for protein sequencing is essential to unraveling the biological mechanisms governing cellular behavior and disease. Current protein sequencing methods suffer from limitations associated with the size of proteins that can be sequenced, the time, and the cost of the sequencing procedures. This study reports the results of all-atom molecular dynamics simulations that investigated the feasibility of using graphene nanopores for protein sequencing. The study is focused on the biologically significant phenylalanine-glycine repeat peptides (FG-nups){\textemdash}parts of the nuclear pore transport machinery. Surprisingly, FG-nups are found to behave similarly to single stranded DNA: The peptides adhere to graphene and exhibit stepwise translocation when subject to a transmembrane bias or a hydrostatic pressure gradient. Reducing the peptide{\textquoteright}s charge density or increasing the peptide{\textquoteright}s hydrophobicity is found to decrease the translocation speed. Yet, unidirectional and stepwise translocation driven by a transmembrane bias is observed even when the ratio of charged to hydrophobic amino acids is as low as 1:8. The nanopore transport of the peptides is found to produce stepwise modulations of the nanopore ionic current correlated with the type of amino acids present in the nanopore, suggesting that protein sequencing by measuring ionic current blockades may be possible.

}, keywords = {graphene, ionic current, nanopore sequencing, transmembrane transport}, issn = {1616-3028}, doi = {10.1002/adfm.201601272}, url = {http://dx.doi.org/10.1002/adfm.201601272}, author = {Wilson, James and Sloman, Leila and He, Zhiren and Aleksei Aksimentiev} } @article {211, title = {Direct evidence for sequence-dependent attraction between double-stranded DNA controlled by methylation}, journal = {Nature Communications}, volume = {7}, year = {2016}, month = {03/2016}, pages = {11045}, abstract = {

Although proteins mediate highly ordered DNA organization\ in vivo, theoretical studies suggest that homologous DNA duplexes can preferentially associate with one another even in the absence of proteins. Here we combine molecular dynamics simulations with single-molecule fluorescence resonance energy transfer experiments to examine the interactions between duplex DNA in the presence of spermine, a biological polycation. We find that AT-rich DNA duplexes associate more strongly than GC-rich duplexes, regardless of the sequence homology. Methyl groups of thymine acts as a steric block, relocating spermine from major grooves to interhelical regions, thereby increasing DNA{\textendash}DNA attraction. Indeed, methylation of cytosines makes attraction between GC-rich DNA as strong as that between AT-rich DNA. Recent genome-wide chromosome organization studies showed that remote contact frequencies are higher for AT-rich and methylated DNA, suggesting that direct DNA{\textendash}DNA interactions that we report here may play a role in the chromosome organization and gene regulation.

}, doi = {10.1038/ncomms11045}, author = {Jejoong Yoo and Hajin Kim and Aleksei Aksimentiev and Taekjip Ha} } @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 {210, title = {Effect of Cytosine Modifications on DNA Flexibility and Nucleosome Mechanical Stability}, journal = {Nature Communications}, volume = {7}, year = {2016}, pages = {10813}, chapter = {10813}, abstract = {

Cytosine can undergo modifications, forming 5-methylcytosine (5-mC) and its oxidized products 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine (5-fC) and 5-carboxylcytosine (5-caC). Despite their importance as epigenetic markers and as central players in cellular processes, it is not well understood how these modifications influence physical properties of DNA and chromatin. Here we report a comprehensive survey of the effect of cytosine modifications on DNA flexibility. We find that even a single copy of 5-fC increases DNA flexibility markedly. 5-mC reduces and 5-hmC enhances flexibility, and 5-caC does not have a measurable effect. Molecular dynamics simulations show that these modifications promote or dampen structural fluctuations, likely through competing effects of base polarity and steric hindrance, without changing the average structure. The increase in DNA flexibility increases the mechanical stability of the nucleosome and vice versa, suggesting a gene regulation mechanism where cytosine modifications change the accessibility of nucleosomal DNA through their effects on DNA flexibility.

}, doi = {10.1038/ncomms10813}, author = {Thuy Ngo and Jejoong Yoo and Dai, Q and Zhang, Q and He, C and Aleksei Aksimentiev and Taekjip Ha} } @article {202, title = {Highly permeable artificial water channels that self-assemble into two-dimensional arrays}, journal = {Proceedings of the National Academy of Science}, volume = {112}, year = {2015}, month = {07/2015}, pages = {9810-9815}, chapter = {9810}, abstract = {

Bioinspired artificial water channels aim to combine the high permeability and selectivity of biological aquaporin (AQP) water channels with chemical stability. We characterized a new architecture of artificial water channels, peptide-appended pillar[5]arenes (PAPs). The average single-channel osmotic water permeability for PAPs is 1.0({\textpm} 0.3){\texttimes}10-14\ cm3/s or 3.5({\textpm} 1.0){\texttimes}108\ water molecules/s, which is in the range of AQPs (3.4~40.3{\texttimes}108\ water molecules/s) and their current synthetic analogs, carbon nanotubes (CNTs, 9.0{\texttimes}108\ water molecules/s). This permeability is an order of magnitude higher than first-generation artificial water channels (20 to ~107\ water molecules/s). Furthermore, within lipid bilayers, PAP channels can self-assemble into two-dimensional arrays. Importantly for permeable membrane design, the pore density of PAP channel arrays (~2.6{\texttimes}105\ pores/μm2) is two orders of magnitude higher than CNT membranes (0.1~2.5{\texttimes}103\ pores/μm2). PAP channels thus combine the advantages of biological channels and CNTs and improve upon them through their relatively simple synthesis, chemical stability and propensity to form arrays.

}, keywords = {membrane, nanochannel, permeation}, doi = {10.1073/pnas.1508575112}, author = {Yue-xiao Shen and Wen Si and Mustafa Erbakan and Karl Decker and Rita De Zorzi and Patrick O. Saboe and You Jung Kang and Sheereen Majd and Peter J. Butler and Thomas Walz and Aleksei Aksimentiev and Jun-li Hou and Manish Kumar} } @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 {123, title = {A Coarse-Grained Model of Unstructured Single-Stranded DNA Derived from Atomistic Simulation and Single-Molecule Experiment}, journal = {Journal of Chemical Theory and Computation}, volume = {10}, year = {2014}, month = {08/2014}, pages = {2891{\textendash}2896}, abstract = {

A simple coarse-grained model of single-stranded DNA (ssDNA) was developed, featuring only two sites per nucleotide that represent the centers of mass of the backbone and sugar/base groups. In the model, the interactions between sites are described using tabulated bonded potentials optimized to reproduce the solution structure of DNA observed in atomistic molecular dynamics simulations. Isotropic potentials describe nonbonded interactions, implicitly taking into account the solvent conditions to match the experimentally determined radius of gyration of ssDNA. The model reproduces experimentally measured force{\textendash}extension dependence of an unstructured DNA strand across 2 orders of magnitude of the applied force. The accuracy of the model was confirmed by measuring the end-to-end distance of a dT14 fragment via FRET while stretching the molecules using optical tweezers. The model offers straightforward generalization to systems containing double-stranded DNA and DNA binding proteins.

}, doi = {10.1021/ct500193u}, url = {http://pubs.acs.org/doi/abs/10.1021/ct500193u}, author = {Christopher Maffeo and Ngo, Thuy T. M. and Ha, Taekjip and Aleksei Aksimentiev} } @article {24, title = {Toward detection of DNA-bound proteins using solid-state nanopores: insights from computer simulations}, journal = {Electrophoresis}, volume = {33}, year = {2012}, month = {2012 Dec}, pages = {3466-79}, abstract = {

Through all-atom molecular dynamics simulations, we explore the use of nanopores in thin synthetic membranes for detection and identification of DNA binding proteins. Reproducing the setup of a typical experiment, we simulate electric field driven transport of DNA-bound proteins through nanopores smaller in diameter than the proteins. As model systems, we use restriction enzymes EcoRI and BamHI specifically and nonspecifically bound to a fragment of dsDNA, and streptavidin and NeutrAvidin proteins bound to dsDNA and ssDNA via a biotin linker. Our simulations elucidate the molecular mechanics of nanopore-induced rupture of a protein-DNA complex, the effective force applied to the DNA-protein bond by the electrophoretic force in a nanopore, and the role of DNA-surface interactions in the rupture process. We evaluate the ability of the nanopore ionic current and the local electrostatic potential measured by an embedded electrode to report capture of DNA, capture of a DNA-bound protein, and rupture of the DNA-protein bond. We find that changes in the strain on dsDNA can reveal the rupture of a protein-DNA complex by altering both the nanopore ionic current and the potential of the embedded electrode. Based on the results of our simulations, we suggest a new method for detection of DNA binding proteins that utilizes peeling of a nicked double strand under the electrophoretic force in a nanopore.

}, keywords = {Bacterial Proteins, Deoxyribonuclease BamHI, Deoxyribonuclease EcoRI, DNA, DNA-Binding Proteins, Electrochemical Techniques, Models, Chemical, Molecular Dynamics Simulation, Nanopores, Protein Binding, Spectrum Analysis, Static Electricity}, issn = {1522-2683}, doi = {10.1002/elps.201200164}, author = {Jeffrey Comer and Anthony H. Ho and Aleksei Aksimentiev} } @article {30, title = {Optimization of the molecular dynamics method for simulations of DNA and ion transport through biological nanopores}, journal = {Methods Mol Biol}, volume = {870}, year = {2012}, month = {2012}, pages = {165-86}, abstract = {

Molecular dynamics (MD) simulations have become a standard method for the rational design and interpretation of experimental studies of DNA translocation through nanopores. The MD method, however, offers a multitude of algorithms, parameters, and other protocol choices that can affect the accuracy of the resulting data as well as computational efficiency. In this chapter, we examine the most popular choices offered by the MD method, seeking an optimal set of parameters that enable the most computationally efficient and accurate simulations of DNA and ion transport through biological nanopores. In particular, we examine the influence of short-range cutoff, integration timestep and force field parameters on the temperature and concentration dependence of bulk ion conductivity, ion pairing, ion solvation energy, DNA structure, DNA-ion interactions, and the ionic current through a nanopore.

}, keywords = {Biological Transport, DNA, Electric Conductivity, Electroosmosis, Hemolysin Proteins, Ions, Lipid Bilayers, Molecular Dynamics Simulation, Nanopores, Porins, Solvents, Thermodynamics}, issn = {1940-6029}, doi = {10.1007/978-1-61779-773-6_10}, author = {David B Wells and Bhattacharya, Swati and Carr, Rogan and Christopher Maffeo and Anthony H. Ho and Jeffrey Comer and Aleksei Aksimentiev} } @article {34, title = {Control of Nanoscale Environment to Improve Stability of Immobilized Proteins on Diamond Surfaces}, journal = {Adv Funct Mater}, volume = {21}, year = {2011}, month = {2011 Mar 21}, pages = {1040-1050}, abstract = {

Immunoassays for detection of bacterial pathogens rely on the selectivity and stability of bio-recognition elements such as antibodies tethered to sensor surfaces. The search for novel surfaces that improve the stability of biomolecules and assay performance has been pursued for a long time. However, the anticipated improvements in stability have not been realized in practice under physiological conditions because the surface functionalization layers on commonly used substrates, silica and gold, are themselves unstable on time scales of days. In this paper, we show that covalent linking of antibodies to diamond surfaces leads to substantial improvements in biological activity of proteins as measured by the ability to selectively capture cells of the pathogenic bacterium Escherichia coli O157:H7 even after exposure to buffer solutions at 37 {\textdegree}C for extended periods of time, approaching 2 weeks. Our results from ELISA, XPS, fluorescence microscopy, and MD simulations suggest that by using highly stable surface chemistry and controlling the nanoscale organization of the antibodies on the surface, it is possible to achieve significant improvements in biological activity and stability. Our findings can be easily extended to functionalization of micro and nanodimensional sensors and structures of biomedical diagnostic and therapeutic interest.

}, issn = {1616-301X}, doi = {10.1002/adfm.201002251}, author = {Radadia, Adarsh D and Courtney J Stavis and Carr, Rogan and Zeng, Hongjun and King, William P and Carlisle, John A and Aleksei Aksimentiev and Hamers, Robert J and Rashid Bashir} } @article {40, title = {Surface functionalization of thin-film diamond for highly stable and selective biological interfaces}, journal = {Proc Natl Acad Sci U S A}, volume = {108}, year = {2011}, month = {2011 Jan 18}, pages = {983-8}, abstract = {

Carbon is an extremely versatile family of materials with a wide range of mechanical, optical, and mechanical properties, but many similarities in surface chemistry. As one of the most chemically stable materials known, carbon provides an outstanding platform for the development of highly tunable molecular and biomolecular interfaces. Photochemical grafting of alkenes has emerged as an attractive method for functionalizing surfaces of diamond, but many aspects of the surface chemistry and impact on biological recognition processes remain unexplored. Here we report investigations of the interaction of functionalized diamond surfaces with proteins and biological cells using X-ray photoelectron spectroscopy (XPS), atomic force microscopy, and fluorescence methods. XPS data show that functionalization of diamond with short ethylene glycol oligomers reduces the nonspecific binding of fibrinogen below the detection limit of XPS, estimated as > 97\% reduction over H-terminated diamond. Measurements of different forms of diamond with different roughness are used to explore the influence of roughness on nonspecific binding onto H-terminated and ethylene glycol (EG)-terminated surfaces. Finally, we use XPS to characterize the chemical stability of Escherichia coli K12 antibodies on the surfaces of diamond and amine-functionalized glass. Our results show that antibody-modified diamond surfaces exhibit increased stability in XPS and that this is accompanied by retention of biological activity in cell-capture measurements. Our results demonstrate that surface chemistry on diamond and other carbon-based materials provides an excellent platform for biomolecular interfaces with high stability and high selectivity.

}, keywords = {Antibodies, Bacterial, Avidin, Chemistry, Physical, Diamond, Escherichia coli K12, Ethylene Glycol, Fibrinogen, Fluorescence, Microscopy, Atomic Force, Photoelectron Spectroscopy, Protein Binding, Surface Properties}, issn = {1091-6490}, doi = {10.1073/pnas.1006660107}, author = {Courtney J Stavis and Clare, Tami Lasseter and Butler, James E and Radadia, Adarsh D and Carr, Rogan and Zeng, Hongjun and King, William P and Carlisle, John A and Aleksei Aksimentiev and Rashid Bashir and Hamers, Robert J} } @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 {54, title = {Computer modeling in biotechnology: a partner in development}, journal = {Methods Mol Biol}, volume = {474}, year = {2008}, month = {2008}, pages = {181-234}, abstract = {

Computational modeling can be a useful partner in biotechnology, in particular, in nanodevice engineering. Such modeling guides development through nanoscale views of biomolecules and devices not available through experimental imaging methods. We illustrate the role of computational modeling, mainly of molecular dynamics, through four case studies: development of silicon bionanodevices for single molecule electrical recording, development of carbon nano-tube-biomolecular systems as in vivo sensors, development of lipoprotein nanodiscs for assays of single membrane proteins, and engineering of oxygen tolerance into the enzyme hydrogenase for photosynthetic hydrogen gas production. The four case studies show how molecular dynamics approaches were adapted to the specific technical uses through (i) multi-scale extensions, (ii) fast quantum chemical force field evaluation, (iii) coarse graining, and (iv) novel sampling methods. The adapted molecular dynamics simulations provided key information on device behavior and revealed development opportunities, arguing that the "computational microscope" is an indispensable nanoengineering tool.

}, keywords = {Biosensing Techniques, Biotechnology, Computational Biology, Computer Simulation, Models, Molecular, Nanostructures}, issn = {1064-3745}, doi = {10.1007/978-1-59745-480-3_11}, author = {Aleksei Aksimentiev and Robert Brunner and Cohen, Jordi and Jeffrey Comer and Eduardo Cruz-Ch{\'u} and Hardy, David and Rajan, Aruna and Amy Y Shih and G Sigalov and Yin, Ying and Klaus Schulten} } @article {68, title = {The electromechanics of DNA in a synthetic nanopore}, journal = {Biophys J}, volume = {90}, year = {2006}, month = {2006 Feb 1}, pages = {1098-106}, abstract = {

We have explored the electromechanical properties of DNA on a nanometer-length scale using an electric field to force single molecules through synthetic nanopores in ultrathin silicon nitride membranes. At low electric fields, E < 200 mV/10 nm, we observed that single-stranded DNA can permeate pores with a diameter >/=1.0 nm, whereas double-stranded DNA only permeates pores with a diameter >/=3 nm. For pores <3.0 nm diameter, we find a threshold for permeation of double-stranded DNA that depends on the electric field and pH. For a 2 nm diameter pore, the electric field threshold is approximately 3.1 V/10 nm at pH = 8.5; the threshold decreases as pH becomes more acidic or the diameter increases. Molecular dynamics indicates that the field threshold originates from a stretching transition in DNA that occurs under the force gradient in a nanopore. Lowering pH destabilizes the double helix, facilitating DNA translocation at lower fields.

}, keywords = {Biophysics, Computer Simulation, DNA, Electrochemistry, Electromagnetic Fields, Electrophoresis, Electrophoresis, Agar Gel, Hydrogen-Ion Concentration, Ions, Membrane Potentials, Microscopy, Electron, Transmission, Models, Molecular, Nanostructures, Nanotechnology, Nucleic Acid Conformation, Polymerase Chain Reaction, Porosity, Silicon Compounds, Time Factors}, issn = {0006-3495}, doi = {10.1529/biophysj.105.070672}, author = {Jiunn B Heng and Aleksei Aksimentiev and Ho, C and Patrick Marks and Yelena V Grinkova and S Sligar and Klaus Schulten and Gregory Timp} } @article {69, title = {Stretching DNA using the electric field in a synthetic nanopore}, journal = {Nano Lett}, volume = {5}, year = {2005}, month = {2005 Oct}, pages = {1883-8}, abstract = {

The mechanical properties of DNA over segments comparable to the size of a protein-binding site (3-10 nm) are examined using an electric-field-induced translocation of single molecules through a nanometer diameter pore. DNA, immersed in an electrolyte, is forced through synthetic pores ranging from 0.5 to 1.5 nm in radius in a 10 nm thick Si(3)N(4) membrane using an electric field. To account for the stretching and bending, we use molecular dynamics to simulate the translocation. We have found a threshold for translocation that depends on both the dimensions of the pore and the applied transmembrane bias. The voltage threshold coincides with the stretching transition that occurs in double-stranded DNA near 60 pN.

}, keywords = {DNA, Nanotechnology, Static Electricity, Stress, Mechanical}, issn = {1530-6984}, doi = {10.1021/nl0510816}, author = {Jiunn B Heng and Aleksei Aksimentiev and Chuen Ho and Patrick Marks and Yelena V Grinkova and S Sligar and Klaus Schulten and Gregory Timp} } @article {59, title = {Beyond the Gene Chip}, journal = {Bell Labs Tech J}, volume = {10}, year = {2005}, month = {2005}, pages = {5-22}, abstract = {

We describe a prospective strategy for reading the encyclopedic information encoded in the genome: using a nanopore in a membrane formed from an MOS-capacitor to sense the charge in DNA. In principle, as DNA permeates the capacitor-membrane through the pore, the electrostatic charge distribution characteristic of the molecule should polarize the capacitor and induce a voltage on the electrodes that can be measured. Silicon nanofabrication and molecular dynamic simulations with atomic detail are technological linchpins in the development of this detector. The sub-nanometer precision available through silicon nanotechnology facilitates the fabrication of the detector, and molecular dynamics provides us with a means to design it and analyze the experimental outcomes.

}, issn = {1089-7089}, doi = {10.1002/bltj.20102}, author = {Jiunn B Heng and Aleksei Aksimentiev and Ho, C and Valentin Dimitrov and Sorsch, T and Miner, J and Mansfield, W and Klaus Schulten and Gregory Timp} } @article {73, title = {Microscopic Kinetics of DNA Translocation through synthetic nanopores}, journal = {Biophys J}, volume = {87}, year = {2004}, month = {2004 Sep}, pages = {2086-97}, abstract = {

We have previously demonstrated that a nanometer-diameter pore in a nanometer-thick metal-oxide-semiconductor-compatible membrane can be used as a molecular sensor for detecting DNA. The prospects for using this type of device for sequencing DNA are avidly being pursued. The key attribute of the sensor is the electric field-induced (voltage-driven) translocation of the DNA molecule in an electrolytic solution across the membrane through the nanopore. To complement ongoing experimental studies developing such pores and measuring signals in response to the presence of DNA, we conducted molecular dynamics simulations of DNA translocation through the nanopore. A typical simulated system included a patch of a silicon nitride membrane dividing water solution of potassium chloride into two compartments connected by the nanopore. External electrical fields induced capturing of the DNA molecules by the pore from the solution and subsequent translocation. Molecular dynamics simulations suggest that 20-basepair segments of double-stranded DNA can transit a nanopore of 2.2 x 2.6 nm(2) cross section in a few microseconds at typical electrical fields. Hydrophobic interactions between DNA bases and the pore surface can slow down translocation of single-stranded DNA and might favor unzipping of double-stranded DNA inside the pore. DNA occluding the pore mouth blocks the electrolytic current through the pore; these current blockades were found to have the same magnitude as the blockade observed when DNA transits the pore. The feasibility of using molecular dynamics simulations to relate the level of the blocked ionic current to the sequence of DNA was investigated.

}, keywords = {Biological Transport, Biophysics, DNA, Electrophoresis, Electrophysiology, Hemolysin Proteins, Ions, Kinetics, Metals, Models, Molecular, Models, Statistical, Nanotechnology, Nucleic Acid Conformation, Oxides, Semiconductors, Stress, Mechanical, Time Factors}, issn = {0006-3495}, doi = {10.1529/biophysj.104.042960}, author = {Aleksei Aksimentiev and Jiunn B Heng and Gregory Timp and Klaus Schulten} } @article {74, title = {Sizing DNA using a nanometer-diameter pore}, journal = {Biophys J}, volume = {87}, year = {2004}, month = {2004 Oct}, pages = {2905-11}, abstract = {

Each species from bacteria to human has a distinct genetic fingerprint. Therefore, a mechanism that detects a single molecule of DNA represents the ultimate analytical tool. As a first step in the development of such a tool, we have explored using a nanometer-diameter pore, sputtered in a nanometer-thick inorganic membrane with a tightly focused electron beam, as a transducer that detects single molecules of DNA and produces an electrical signature of the structure. When an electric field is applied across the membrane, a DNA molecule immersed in electrolyte is attracted to the pore, blocks the current through it, and eventually translocates across the membrane as verified unequivocally by gel electrophoresis. The relationship between DNA translocation and blocking current has been established through molecular dynamics simulations. By measuring the duration and magnitude of the blocking current transient, we can discriminate single-stranded from double-stranded DNA and resolve the length of the polymer.

}, keywords = {DNA, Electrochemistry, Nanotechnology, Particle Size, Porosity, Reproducibility of Results, Sensitivity and Specificity, Ultrafiltration}, issn = {0006-3495}, doi = {10.1529/biophysj.104.041814}, author = {Jiunn B Heng and Chuen Ho and Kim, Taekyung and Timp, Rolf and Aleksei Aksimentiev and Yelena V Grinkova and S Sligar and Klaus Schulten and Gregory Timp} } @article {183, title = {Morphology of surfaces in mesoscopic polymers, surfactants, electrons, or reaction-diffusion systems: Methods, simulation and measurements}, journal = {Advances in Chemical Physics}, volume = {121}, year = {2002}, pages = {141{\textendash}239}, abstract = {

In this chapter the authors discuss systems with internal surfaces, ordered surfaces, topological transformations, and dynamical scaling. In Section II the authors show specific examples of mesoscopic systems with special attention devoted to the surfaces in the system{\textemdash}that is, periodic surfaces in surfactant systems, periodic surfaces in diblock copolymers, bicontinuous disordered interfaces in spinodally decomposing blends, ordered charge density wave patterns in electron liquids, and dissipative structures in reaction{\textendash}diffusion systems. Section III covers the theory of morphological measures, including the Euler characteristic and the Gaussian and mean curvatures.

}, doi = {10.1002/0471264318.ch3}, author = {Aleksei Aksimentiev and Robert Ho{\l}yst and Marcin Fialkowski} } @article {79, title = {Scaling of the Euler characteristic, surface area, and curvatures in the phase separating or ordering systems}, journal = {Phys Rev Lett}, volume = {86}, year = {2001}, month = {2001 Jan 8}, pages = {240-3}, abstract = {

We present robust scaling laws for the Euler characteristic and curvatures applicable to any symmetric system undergoing phase separating or ordering kinetics. We apply it to the phase ordering in a system of the nonconserved scalar order parameter and find three scaling regimes. The appearance of the preferred nonzero curvature of an interface separating +/- domains marks the crossover to the late stage regime characterized by the Lifshitz-Cahn-Allen scaling.

}, issn = {0031-9007}, author = {Marcin Fialkowski and Aleksei Aksimentiev and Robert Ho{\l}yst} } @article {182, title = {Application of the Euler characteristic to the study of homopolymer blends and copolymer melts}, journal = {Polimery}, volume = {46}, year = {2001}, pages = {309{\textendash}322}, abstract = {

The Euler characteristic is applied to study the phase separation/ordering phenomena occurring in polymer mixtures. Several methods used to compute Euler characteristics are described and discussed in terms of accuracy and efficiency. A novel, robust and exact algorithm based on simplex decomposition is presented. The simulated morphology of spinodal decomposition and the phase equilibrium diagram of gradient copolymers are analyzed in terms of Euler characteristics. Again, the Euler characteristics are applied to specify the connectivity of bicontinuous patterns, to distinguish the disperse and the bicontinuous morphologies, and to detect the percolarion transition. During the spinodal decomposition, the dynamic transformation between the bicontinuous and the droplet patterns is shown to involve the formation of a transient "cylindrical" morphology.

}, url = {http://en.www.ichp.pl/Application-of-the-Euler-characteristic}, author = {Aleksei Aksimentiev and Robert Ho{\l}yst} } @article {80, title = {Influence of the free-energy functional form on simulated morphology of spinodally decomposing blends}, journal = {Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics}, volume = {62}, year = {2000}, month = {2000 Nov}, pages = {6821-30}, abstract = {

The spinodal decomposition of a binary mixture has been studied within several mesoscopic models. It has been found that the form of the equilibrium free energy has a crucial effect on the morphological development in asymmetric blends. We have shown that the principal quantity that determines the topology of the interface (and type of morphology) is the equilibrium minority phase volume fraction, while the transition from bicontinuous to droplet morphology can be treated as a percolation. The concentration dependence of the square gradient coefficient attributed for the Flory-Huggins-de Gennes free energy has no significant influences on the average domain growth, but can be distinguished experimentally from its constant-coefficient alternative by measuring the maximum wave vector of the scattering intensity as a function of the minority phase volume fraction for spinodally decomposing asymmetric blends. The concentration dependence of the Onsager coefficient has the weak, systematic effect of slowing down the morphological development. The local shape of the interface is not affected considerably by the concentration dependence of the square gradient and Onsager coefficient.

}, issn = {1063-651X}, author = {Aleksei Aksimentiev and Robert Ho{\l}yst} } @article {181, title = {Scaling properties of the morphological measures at the early and intermediate stages of the spinodal decomposition in homopolymer blends}, journal = {The Journal of Chemical Physics}, volume = {112}, year = {2000}, pages = {6049{\textendash}6062}, abstract = {

The spinodal decomposition of the homopolymer blends has been studied by the numerical integration of the Cahn{\textendash}Hilliard{\textendash}Cook equation. We have investigated the time evolution of the morphological measures that characterize quantitatively the interface in the system. For symmetric blends we have found that the Euler characteristic of the interface is negative and increases with time as τ0.75 (connectivity of the domains decreases) regardless of the final quench temperature. The homogeneity index of the interface is constant in this case. This suggests that at the level of the integral geometry quantities (Minkowski functionals), the dynamic scaling hypothesis holds for the evolution of the interface morphology in quenched critical systems. The nonuniversal morphological evolution of the asymmetric blends have been studied. Also, we have shown that the thermal fluctuations can modify significantly the curvature distribution.

}, doi = {10.1063/1.481178}, author = {Aleksei Aksimentiev and Krzysztof Moorthi and Robert Ho{\l}yst} } @article {180, title = {Kinetics of the droplet formation at the early and intermediate stages of the spinodal decomposition in homopolymer blends}, journal = {Macromolecular Theory and Simulations}, volume = {9}, year = {2000}, pages = {661{\textendash}674}, abstract = {

The kinetics of the droplet formation during the spinodal decomposition (SD) of the homopolymer blends has been studied by numerical integration of the Cahn-Hilliard-Cook equation. We have found that the droplet formation and growth occurs when the minority phase volume fraction, fm , approaches the percolation threshold value, fthr = 0.3 {\textpm} 0.01. The time for the formation of the disperse droplet morphology (coarsening time) depends only on the equilibrium minority phase volume fraction, fm . fm approaches its equilibrium value logarithmically at the late SD stages, and, therefore, the coarsening time decreases exponentially as the average volume fraction or the quench depth decrease. Since the temporal evolution of the total interfacial area does not depend on the quench conditions and blend morphology, the average droplet size and the droplet number density is determined by the coarsening time. Within the time scale studied, the droplet number density decreases with time as t{\textendash}0.63{\textpm}0.03; the average mean curvature decreases as t{\textendash}0.35{\textpm}0.05; the average Gaussian curvature decreases as t{\textendash}0.42{\textpm}0.03, and the average droplet compactness \~{}V/S3/2 where S is the surface area and V is the volume) approaches a spherical limit logarithmically with time. The droplets with larger area have lower compactness and in the low compactness limit their area is a parabolic function of compactness. The size and shape distribution functions have been also investigated.

}, doi = {10.1002/1521-3919(20001101)9:8<661::AID-MATS661>3.0.CO;2-6}, author = {Aleksei Aksimentiev and Robert Ho{\l}yst and Krzysztof Moorthi} } @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} } @article {178, title = {Phase behavior of gradient copolymers}, journal = {The Journal of Chemical Physics}, volume = {111}, year = {1999}, pages = {2329{\textendash}2339}, abstract = {

Phase properties of gradient ABcopolymer melts which consist of chains with the specified chemical distribution of A and B monomers have been studied within the Landau{\textendash}Ginzburg model. All the melts with the linear distribution of the monomers exhibit only a direct continuous phase transition from disordered to the lamellar phase. The hexagonal, body-centered-cubic, double-gyroid (G), and lamellar-ordered structures have been found in the melts with the monotonic but nonlinear distribution of the monomers. The G structure has been also found in the gradient copolymer melts with the distribution function of monomers similar to the A{\textendash}B{\textendash}A triblock copolymers.

}, doi = {10.1063/1.479504}, author = {Aleksei Aksimentiev and Robert Ho{\l}yst} } @article {177, title = {Single chain conformations in the system of symmetric and asymmetric diblock copolymers}, journal = {Macromolecular Theory and Simulations}, volume = {8}, year = {1999}, pages = {328{\textendash}342}, abstract = {

A theory which describes a local structure and global properties of a diblock copolymer melt has been developed in the framework of the one-loop self-consistent approximation. We have derived expressions for the sizes of a single diblock macromolecule and its parts. Two different behaviors of single macromolecule conformations in the disordered melt have been obtained depending on the asymmetry of chains and morphologies occurring in ordered states after the order-disorder transition (ODT). In the nearly symmetric melt, 0.35 \< f ⪇ 0.5 (f is a composition), the blocks of both types shrink a little initially as the temperature decreases and then, at some temperature, they begin to swell. In strongly asymmetric melts, f \< 0.35, the block of a macromolecule which consists of the monomers of minority type shrinks monotonically, while the other block monotonically swells. We have found nearly Gaussian behavior of the individual blocks and stretching near the chemical bond joining the blocks. Near the ODT the chains are stretched with a magnitude which is of the order of a few percent of their Gaussian sizes. We have calculated the peak position in the scattering curve as a function of the Flory-Huggins interaction parameter, composition and degree of polymerization. Less then 5\% change in the size of copolymer molecules lead to a 25\% shift of the scattering peak in comparison to the Gaussian limit. We have found a good quantitative agreement of our theoretical results with the experimental neutron scattering data.

}, doi = {10.1002/(SICI)1521-3919(19990701)8:4<328::AID-MATS328>3.0.CO;2-3}, author = {Aleksei Aksimentiev and Robert Ho{\l}yst} } @article {179, title = {Single-chain statistics in polymer systems}, journal = {Progress in Polymer Science}, volume = {24}, year = {1999}, pages = {1045{\textendash}1093}, abstract = {

In this review we study the behavior of a single labelled polymer chain in various polymer systems: polymer blends, diblock copolymers, gradient copolymers, ring copolymers, polyelectrolytes, grafted homopolymers, rigid nematogenic polymers, polymers in bad and good solvents, fractal polymers and polymers in fractal environments. We discuss many phenomena related to the single chain behavior, such as: collapse of polymers in bad solvents, protein folding, stretching of polymer brushes, coil{\textendash}rod transition in nematogenic main-chain polymers, knot formation in homopolymer melts, and shrinking and swelling of polymers at temperatures close to the bulk transition temperatures. Our description is mesoscopic, based on two models of polymer systems: the Edwards model with Fixman delta interactions, and the Landau{\textendash}Ginzburg model of phase transitions applied to polymers. In particular, we show the derivation of the Landau{\textendash}Ginzburg model from the Edwards model in the case of homopolymer blends and diblock copolymer melts. In both models, we calculate the radius of gyration and relate them to the correlation function for a single polymer chain. We discuss theoretical results as well as computer simulations and experiments.

}, doi = {10.1016/S0079-6700(99)00023-4}, author = {Aleksei Aksimentiev and Robert Ho{\l}yst} } @article {176, title = {Swelling and shrinking of polymer chains in homopolymer blends}, journal = {Macromolecular Theory and Simulations}, volume = {7}, year = {1998}, pages = {447{\textendash}456}, abstract = {

The radius of gyration, R, of polymer chains in homopolymer blends is studied in the framework of a self-consistent one-loop approximation. We show that the polymer chains can shrink or swell in comparison to the Gaussian chain. Swelling of the polymer chains in a region far away from the critical point is caused by the steric repulsive forces that were included in the model as the constraint of incompressibility. The chains shrink progressively, as we approach the critical region passing through the Gaussian limit, R0 = sqrt(N/6) l, far away from the critical point (N - degree of polymerization, l - length of monomer). The correction responsible for the swelling and the shrinking is small when the concentrations of components ϕ are comparable (N = 1000, \<ϕ\> = 0.5, (R02-R2)/R02 =\ {\textpm}0.02\%). This effect, although small, leads to a local demixing, a sharp shrinking of chains in both components accompanied by a strong change of the global inter-monomer distance, which should be observable experimentally. When the local demixing occurs there is a sudden increase in the scattering intensity (of the order of 30\% for N = 1000, and ϕA = 0.2). The increase of the degree of polymerization of the same type of chains leads to an increase of the swelling-shrinking effects. In addition, the critical concentration of the shorter chains component is smaller in comparison to the value obtained in the Flory-Huggins theory. The self-consistent determination of the radius of gyration and the upper wave-vector cutoff make our model free from any divergences. In the limit of N\ {\textrightarrow} $\infty$ this theory reduces to the random phase approximation (RPA) of de Gennes.

}, doi = {10.1002/(SICI)1521-3919(19990701)8:4<328::AID-MATS328>3.0.CO;2-3}, author = {Aleksei Aksimentiev and Robert Ho{\l}yst} }