Within just a decade from the pioneering work demonstrating the utility of nanopores for molecular sensing, nanopores have emerged as versatile systems for single-molecule manipulation and analysis. ...In a typical setup, a gradient of the electrostatic potential captures charged solutes from the solution and forces them to move through a single nanopore, across an otherwise impermeable membrane. The ionic current blockades resulting from the presence of a solute in a nanopore can reveal the type of the solute, for example, the nucleotide makeup of a DNA strand. Despite great success, the microscopic mechanisms underlying the functionality of such stochastic sensors remain largely unknown, as it is not currently possible to characterize the microscopic conformations of single biomolecules directly in a nanopore and thereby unequivocally establish the causal relationship between the observables and the microscopic events. Such a relationship can be determined using molecular dynamics-a computational method that can accurately predict the time evolution of a molecular system starting from a given microscopic state. This article describes recent applications of this method to the process of DNA transport through biological and synthetic nanopores.
In contrast to ordinary polymers, the vast majority of biological macromolecules adopt highly ordered three-dimensional structures that define their functions. The key to folding of a biopolymer into ...a unique 3D structure or to an assembly of several biopolymers into a functional unit is a delicate balance between the attractive and repulsive forces that also makes such self-assembly reversible under physiological conditions. The all-atom molecular dynamics (MD) method has emerged as a powerful tool for studies of individual biomolecules and their functional assemblies, encompassing systems of ever increasing complexity. However, advances in parallel computing technology have outpaced the development of the underlying theoretical models-the molecular force fields, pushing the MD method into an untested territory. Recent tests of the MD method have found the most commonly used molecular force fields to be out of balance, overestimating attractive interactions between charged and hydrophobic groups, which can promote artificial aggregation in MD simulations of multi-component protein, nucleic acid, and lipid systems. One route towards improving the force fields is through the NBFIX corrections method, in which the intermolecular forces are calibrated against experimentally measured quantities such as osmotic pressure by making atom pair-specific adjustments to the non-bonded interactions. In this article, we review development of the NBFIX (Non-Bonded FIX) corrections to the AMBER and CHARMM force fields and discuss their implications for MD simulations of electrolyte solutions, dense DNA systems, Holliday junctions, protein folding, and lipid bilayer membranes.
Efforts to sequence single protein molecules in nanopores
have been hampered by the lack of techniques with sufficient sensitivity to discern the subtle molecular differences among all twenty amino ...acids. Here we report ionic current detection of all twenty proteinogenic amino acids in an aerolysin nanopore with the help of a short polycationic carrier. Application of molecular dynamics simulations revealed that the aerolysin nanopore has a built-in single-molecule trap that fully confines a polycationic carrier-bound amino acid inside the sensing region of the aerolysin. This structural feature means that each amino acid spends sufficient time in the pore for sensitive measurement of the excluded volume of the amino acid. We show that distinct current blockades in wild-type aerolysin can be used to identify 13 of the 20 natural amino acids. Furthermore, we show that chemical modifications, instrumentation advances and nanopore engineering offer a route toward identification of the remaining seven amino acids. These findings may pave the way to nanopore protein sequencing.
Over the past decades, molecular dynamics (MD) simulations of biomolecules have become a mainstream biophysics technique. As the length and time scales amenable to the MD method increase, ...shortcomings of the empirical force fields, which have been developed and validated using relatively short simulations of small molecules, become apparent. One common artifact is aggregation of water-soluble biomolecules driven by artificially strong charge-charge interactions. Here, we report a systematic atom pair-specific refinement of Lennard-Jones parameters (NBFIX) describing amine-carboxylate and amine-phosphate interactions, which bring MD simulations of basic peptide-mediated nucleic acid assemblies and lipid bilayer membranes into better agreement with experimental data. As our refinement does not affect the existing parametrization of bonded interactions or alter the solvation free energies, it improves the realism of an MD simulation without introducing additional artifacts.
The DNA origami method permits folding of long single-stranded DNA into complex 3D structures with subnanometer precision. Transmission electron microscopy, atomic force microscopy, and recently ...cryo-EM tomography have been used to characterize the properties of such DNA origami objects, however their microscopic structures and dynamics have remained unknown. Here, we report the results of all-atom molecular dynamics simulations that characterized the structural and mechanical properties of DNA origami objects in unprecedented microscopic detail. When simulated in an aqueous environment, the structures of DNA origami objects depart from their idealized targets as a result of steric, electrostatic, and solvent-mediated forces. Whereas the global structural features of such relaxed conformations conform to the target designs, local deformations are abundant and vary in magnitude along the structures. In contrast to their free-solution conformation, the Holliday junctions in the DNA origami structures adopt a left-handed antiparallel conformation. We find the DNA origami structures undergo considerable temporal fluctuations on both local and global scales. Analysis of such structural fluctuations reveals the local mechanical properties of the DNA origami objects. The lattice type of the structures considerably affects global mechanical properties such as bending rigidity. Our study demonstrates the potential of all-atom molecular dynamics simulations to play a considerable role in future development of the DNA origami field by providing accurate, quantitative assessment of local and global structural and mechanical properties of DNA origami objects.
Meters of DNA wrap around histone proteins to form nucleosomes and fit inside the micron-diameter nucleus. For the genetic information encoded in the DNA to become available for transcription, ...replication, and repair, the DNA–histone assembly must be disrupted. Experiment has indicated that the outer stretches of nucleosomal DNA “breathe” by spontaneously detaching from and reattaching to the histone core. Here, we report direct observation of spontaneous DNA breathing in atomistic molecular dynamics simulations, detailing a microscopic mechanism of the DNA breathing process. According to our simulations, the outer stretches of nucleosomal DNA detach in discrete steps involving 5 or 10 base pairs, with the detachment process being orchestrated by the motion of several conserved histone residues. The inner stretches of nucleosomal DNA are found to be more stably associated with the histone core by more abundant nonspecific DNA–protein contacts, providing a microscopic interpretation of nucleosome unraveling experiments. The CG content of nucleosomal DNA is found to anticorrelate with the extent of unwrapping, supporting the possibility that AT-rich segments may signal the start of transcription by forming less stable nucleosomes.
Display omitted
•Outer stretches of DNA spontaneously and reversibly detach in all-atom MD simulations.•Each turn of DNA detaches after the sequential breaking of two salt bridges.•Several highly conserved histone residues orchestrate the process of DNA detachment.•Less CG content was found to correlate with a greater nucleosome unraveling.
Control over interactions with biomolecules holds the key to applications of graphene in biotechnology. One such application is nanopore sequencing, where a DNA molecule is electrophoretically driven ...through a graphene nanopore. Here we investigate how interactions of single-stranded DNA and a graphene membrane can be controlled by electrically biasing the membrane. The results of our molecular dynamics simulations suggest that electric charge on graphene can force a DNA homopolymer to adopt a range of strikingly different conformations. The conformational response is sensitive to even very subtle nucleotide modifications, such as DNA methylation. The speed of DNA motion through a graphene nanopore is strongly affected by the graphene charge: a positive charge accelerates the motion, whereas a negative charge arrests it. As a possible application of the effect, we demonstrate stop-and-go transport of DNA controlled by the charge of graphene. Such on-demand transport of DNA is essential for realizing nanopore sequencing.
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–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–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 “rewind” peptide reads, obtaining numerous independent reads of the same molecule, yielding an error rate of <10
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.
Spontaneous assembly of DNA molecules into compact structures is ubiquitous in biological systems. Experiment has shown that polycations can turn electrostatic self-repulsion of DNA into attraction, ...yet the physical mechanism of DNA condensation has remained elusive. Here, we report the results of atomistic molecular dynamics simulations that elucidated the microscopic structure of dense DNA assemblies and the physics of interactions that makes such assemblies possible. Reproducing the setup of the DNA condensation experiments, we measured the internal pressure of DNA arrays as a function of the DNA-DNA distance, showing a quantitative agreement between the results of our simulations and the experimental data. Analysis of the MD trajectories determined the DNA-DNA force in a DNA condensate to be pairwise, the DNA condensation to be driven by electrostatics of polycations and not hydration, and the concentration of bridging cations, not adsorbed cations, to determine the magnitude and the sign of the DNA-DNA force. Finally, our simulations quantitatively characterized the orientational correlations of DNA in DNA arrays as well as diffusive motion of DNA and cations.
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.
PDF
