Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism.

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'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.

The structure and physical properties of a packaged bacteriophage particle.

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles.

Eukaryotic genome and methylome encode DNA fragments' propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between 'omics' 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' local geometry consists of straight segments connected by sharp bends compressing the DNA'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' 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' nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.

Multi-resolution simulation of DNA transport through large synthetic nanostructures.

Modeling and simulation has become an invaluable partner in development of nanopore sensing systems. The key advantage of the nanopore sensing method - the ability to rapidly detect individual biomolecules as a transient reduction of the ionic current flowing through the nanopore - is also its key deficiency, as the current signal itself rarely provides direct information about the chemical structure of the biomolecule. Complementing experimental calibration of the nanopore sensor readout, coarse-grained and all-atom molecular dynamics simulations have been used extensively to characterize the nanopore translocation process and to connect the microscopic events taking place inside the nanopore to the experimentally measured ionic current blockades. Traditional coarse-grained simulations, however, lack the precision needed to predict ionic current blockades with atomic resolution whereas traditional all-atom simulations are limited by the length and time scales amenable to the method. Here, we describe a multi-resolution framework for modeling electric field-driven passage of DNA molecules and nanostructures through to-scale models of synthetic nanopore systems. We illustrate the method by simulating translocation of double-stranded DNA through a solid-state nanopore and a micron-scale slit, capture and translocation of single-stranded DNA in a double nanopore system, and modeling ionic current readout from a DNA origami nanostructure passage through a nanocapillary. We expect our multi-resolution simulation framework to aid development of the nanopore field by providing accurate, to-scale modeling capability to research laboratories that do not have access to leadership supercomputer facilities.

MrDNA: a multi-resolution model for predicting the structure and dynamics of DNA systems.

Although the field of structural DNA nanotechnology has been advancing with an astonishing pace, de novo design of complex 3D nanostructures and functional devices remains a laborious and time-consuming process. One reason for that is the need for multiple cycles of experimental characterization to elucidate the effect of design choices on the actual shape and function of the self-assembled objects. Here, we demonstrate a multi-resolution simulation framework, mrdna, that, in 30 min or less, can produce an atomistic-resolution structure of a self-assembled DNA nanosystem. We demonstrate fidelity of our mrdna framework through direct comparison of the simulation results with the results of cryo-electron microscopy (cryo-EM) reconstruction of multiple 3D DNA origami objects. Furthermore, we show that our approach can characterize an ensemble of conformations adopted by dynamic DNA nanostructures, the equilibrium structure and dynamics of DNA objects constructed using off-lattice self-assembly principles, i.e. wireframe DNA objects, and to study the properties of DNA objects under a variety of environmental conditions, such as applied electric field. Implemented as an open source Python package, our framework can be extended by the community and integrated with DNA design and molecular graphics tools.

Inchworm movement of two rings switching onto a thread by biased Brownian diffusion represent a three-body problem.

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 "three-body problem." 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-thread interactions for kinetics, and ring-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.

Molecular mechanism of DNA association with single-stranded DNA binding protein.

During DNA replication, the single-stranded DNA binding protein (SSB) wraps single-stranded DNA (ssDNA) with high affinity to protect it from degradation and prevent secondary structure formation. Although SSB binds ssDNA tightly, it can be repositioned along ssDNA to follow the advancement of the replication fork. Using all-atom molecular dynamics simulations, we characterized the molecular mechanism of ssDNA association with SSB. Placed in solution, ssDNA-SSB assemblies were observed to change their structure spontaneously; such structural changes were suppressed in the crystallographic environment. Repeat simulations of the SSB-ssDNA complex under mechanical tension revealed a multitude of possible pathways for ssDNA to come off SSB punctuated by prolonged arrests at reproducible sites at the SSB surface. Ensemble simulations of spontaneous association of short ssDNA fragments with SSB detailed a three-dimensional map of local affinity to DNA; the equilibrium amount of ssDNA bound to SSB was found to depend on the electrolyte concentration but not on the presence of the acidic tips of the SSB tails. Spontaneous formation of ssDNA bulges and their diffusive motion along SSB surface was directly observed in multiple 10-µs-long simulations. Such reptation-like motion was confined by DNA binding to high-affinity spots, suggesting a two-step mechanism for SSB diffusion.

Optical Voltage Sensing Using DNA Origami.

We explore the potential of DNA nanotechnology for developing novel optical voltage sensing nanodevices 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 labeled with a single pair of FRET dyes. The DNA structures were reversibly immobilized 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 characterized 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.

De novo reconstruction of DNA origami structures through atomistic molecular dynamics simulation.

The DNA origami method has brought nanometer-precision fabrication to molecular biology labs, offering myriads of potential applications in the fields of synthetic biology, medicine, molecular computation, etc. Advancing the method further requires controlling self-assembly down to the atomic scale. Here we demonstrate a computational method that allows the equilibrium structure of a large, complex DNA origami object to be determined to atomic resolution. Through direct comparison with the results of cryo-electron microscopy, we demonstrate de novo reconstruction of a 4.7 megadalton pointer structure by means of fully atomistic molecular dynamics simulations. Furthermore, we show that elastic network-guided simulations performed without solvent can yield similar accuracy at a fraction of the computational cost, making this method an attractive approach for prototyping and validation of self-assembled DNA nanostructures.

End-to-end attraction of duplex DNA.

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'-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.

Nanopore Translocation Reveals Electrophoretic Force on Noncanonical RNA:DNA Double Helix.

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on the translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of noncanonical RNA:DNA hybrids in electrophoretic transport to the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a noncanonical helix, with distinct transport properties from B-form DD molecules. We find that RD and DD molecules, with the same contour length, move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and the molecular understanding of electrophoretic transport.

Single molecule delivery into living cells.

Controlled manipulation of cultured cells by delivery of exogenous macromolecules is a cornerstone of experimental biology. Here we describe a platform that uses nanopipettes to deliver defined numbers of macromolecules into cultured cell lines and primary cells at single molecule resolution. In the nanoinjection platform, the nanopipette is used as both a scanning ion conductance microscope (SICM) probe and an injection probe. The SICM is used to position the nanopipette above the cell surface before the nanopipette is inserted into the cell into a defined location and to a predefined depth. We demonstrate that the nanoinjection platform enables the quantitative delivery of DNA, globular proteins, and protein fibrils into cells with single molecule resolution and that delivery results in a phenotypic change in the cell that depends on the identity of the molecules introduced. Using experiments and computational modeling, we also show that macromolecular crowding in the cell increases the signal-to-noise ratio for the detection of translocation events, thus the cell itself enhances the detection of the molecules delivered.

A DNA turbine powered by a transmembrane potential across a nanopore.

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

DNA double helix, a tiny electromotor.

Flowing fluid past chiral objects has been used for centuries to power rotary motion in man-made machines. By contrast, rotary motion in nanoscale biological or chemical systems is produced by biasing Brownian motion through cyclic chemical reactions. Here we show that a chiral biological molecule, a DNA or RNA duplex rotates unidirectionally at billions of revolutions per minute when an electric field is applied along the duplex, with the rotation direction being determined by the chirality of the duplex. The rotation is found to be powered by the drag force of the electro-osmotic flow, realizing the operating principle of a macroscopic turbine at the nanoscale. The resulting torques are sufficient to power rotation of nanoscale beads and rods, offering an engineering principle for constructing nanoscale systems powered by electric field.

Super-Resolution Detection of DNA Nanostructures Using a Nanopore.

High-resolution analysis of biomolecules has brought unprecedented insights into fundamental biological processes and dramatically advanced biosensing. Notwithstanding the ongoing resolution revolution in electron microscopy and optical imaging, only a few methods are presently available for high-resolution analysis of unlabeled single molecules in their native states. Here, label-free electrical sensing of structured single molecules with a spatial resolution down to single-digit nanometers is demonstrated. Using a narrow solid-state nanopore, the passage of a series of nanostructures attached to a freely translocating DNA molecule is detected, resolving individual nanostructures placed as close as 6 nm apart and with a surface-to-surface gap distance of only 2 nm. Such super-resolution ability is attributed to the nanostructure-induced enhancement of the electric field at the tip of the nanopore. This work demonstrates a general approach to improving the resolution of single-molecule nanopore sensing and presents a critical advance towards label-free, high-resolution DNA sequence mapping, and digital information storage independent of molecular motors.

Single-molecule biophysics experiments in silico: Toward a physical model of a replisome.

The interpretation of single-molecule experiments is frequently aided by computational modeling of biomolecular dynamics. The growth of computing power and ongoing validation of computational models suggest that it soon may be possible to replace some experiments outright with computational mimics. Here, we offer a blueprint for performing single-molecule studies in silico using a DNA-binding protein as a test bed. We demonstrate how atomistic simulations, typically limited to sub-millisecond durations and zeptoliter volumes, can guide development of a coarse-grained model for use in simulations that mimic single-molecule experiments. We apply the model to recapitulate, in silico, force-extension characterization of protein binding to single-stranded DNA and protein and DNA replacement assays, providing a detailed portrait of the underlying mechanics. Finally, we use the model to simulate the trombone loop of a replication fork, a large complex of proteins and DNA.

Percolation transition prescribes protein size-specific barrier to passive transport through the nuclear pore complex.

Nuclear pore complexes (NPCs) control biomolecular transport in and out of the nucleus. Disordered nucleoporins in the complex's pore form a permeation barrier, preventing unassisted transport of large biomolecules. Here, we combine coarse-grained simulations of experimentally derived NPC structures with a theoretical model to determine the microscopic mechanism of passive transport. Brute-force simulations of protein transport reveal telegraph-like behavior, where prolonged diffusion on one side of the NPC is interrupted by rapid crossings to the other. We rationalize this behavior using a theoretical model that reproduces the energetics and kinetics of permeation solely from statistics of transient voids within the disordered mesh. As the protein size increases, the mesh transforms from a soft to a hard barrier, enabling orders-of-magnitude reduction in permeation rate for proteins beyond the percolation size threshold. Our model enables exploration of alternative NPC architectures and sets the stage for uncovering molecular mechanisms of facilitated nuclear transport.

Leakless end-to-end transport of small molecules through micron-length DNA nanochannels.

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.

Trapping of protein cargo molecules inside DNA origami nanocages.

The development of the DNA origami technique has directly inspired the idea of using three-dimensional DNA cages for the encapsulation and targeted delivery of drug or cargo molecules. The cages would be filled with molecules that would be released at a site of interest upon cage opening triggered by an external stimulus. Though different cage variants have been developed, efficient loading of DNA cages with freely-diffusing cargo molecules that are not attached to the DNA nanostructure and their efficient retention within the cages has not been presented. Here we address these challenges using DNA origami nanotubes formed by a double-layer of DNA helices that can be sealed with tight DNA lids at their ends. In a first step we attach DNA-conjugated cargo proteins to complementary target strands inside the DNA tubes. After tube sealing, the cargo molecules are released inside the cavity using toehold-mediated strand displacement by externally added invader strands. We show that DNA invaders are rapidly entering the cages through their DNA walls. Retention of ∼70 kDa protein cargo molecules inside the cages was, however, poor. Guided by coarse-grained simulations of the DNA cage dynamics, a tighter sealing of the DNA tubes was developed which greatly reduced the undesired escape of cargo proteins. These improved DNA nanocages allow for efficient encapsulation of medium-sized cargo molecules while remaining accessible to small molecules that can be used to trigger reactions, including a controlled release of the cargo via nanocage opening.

Discrimination of RNA fiber structures using solid-state nanopores.

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-60 °C revealed almost identical distributions of current blockade amplitudes, suggesting that the kissing loop interactions in fibers are resistant to heating within this range.