Exploring transmembrane transport through alpha-hemolysin with grid-steered molecular dynamics

Genetically programmed synthetic cells for thermo-responsive protein synthesis and cargo release

Synthetic cells containing genetic programs and protein expression machinery are increasingly recognized as powerful counterparts to engineered living cells in the context of biotechnology, therapeutics and cellular modelling. So far, genetic regulation of synthetic cell activity has been largely confined to chemical stimuli; to unlock their potential in applied settings, engineering stimuli-responsive synthetic cells under genetic regulation is imperative. Here we report the development of temperature-sensitive synthetic cells that control protein production by exploiting heat-responsive mRNA elements. This is achieved by combining RNA thermometer technology, cell-free protein expression and vesicle-based synthetic cell design to create cell-sized capsules able to initiate synthesis of both soluble proteins and membrane proteins at defined temperatures. We show that the latter allows for temperature-controlled cargo release phenomena with potential implications for biomedicine. Platforms like the one presented here can pave the way for customizable, genetically programmed synthetic cells under thermal control to be used in biotechnology.

Protein Deceleration and Sequencing Using Si3N4-CNT Hybrid Nanopores

Protein sequencing is crucial for understanding the complex mechanisms driving biological functions and is of utmost importance in molecular diagnostics and medication development. Nanopores have become an effective tool for single molecule sensing, however, the weak charge and non-uniform charge distribution of protein make capturing and sensing very challenging, which poses a significant obstacle to the development of nanopore-based protein sequencing. In this study, to facilitate capturing of the unfolded protein, highly charged peptide was employed in our simulations, we found that the velocity of unfolded peptide translocating through a hybrid nanopore composed of silicon nitride membrane and carbon nanotube is much slower compared to bare silicon nitride nanopore, it is due to the significant interaction between amino acids and the surface of carbon nanotube. Moreover, by introducing variations in the charge states at the boundaries of carbon nanotube nanopores, the competition and combination of the electrophoretic and electroosmotic flows through the nanopores could be controlled, we then successfully regulated the translocation velocity of unfolded proteins through the hybrid nanopores. The proposed hybrid nanopore effectively retards the translocation velocity of protein through it, facilitates the acquisition of ample information for accurate amino acid identification.

Generating Concentration Gradients across Membranes for Molecular Dynamics Simulations of Periodic Systems

The plasma membrane forms the boundary between a living entity and its environment and acts as a barrier to permeation and flow of substances. Several computational means of calculating permeability have been implemented for molecular dynamics (MD) simulations-based approaches. Except for double bilayer systems, most permeability studies have been performed under equilibrium conditions, in large part due to the challenges associated with creating concentration gradients in simulations utilizing periodic boundary conditions. To enhance the scientific understanding of permeation and complement the existing computational means of characterizing membrane permeability, we developed a non-equilibrium method that enables the generation and maintenance of steady-state gradients in MD simulations. We utilize PBCs advantageously by imposing a directional bias to the motion of permeants so that their crossing of the boundary replenishes the gradient, like a previous study on ions. Under these conditions, a net flow of permeants across membranes may be observed to determine bulk permeability by a direct application of J=PΔc. In the present study, we explore the results of its application to an exemplary O2 and POPC bilayer system, demonstrating accurate and precise permeability measurements. In addition, we illustrate the impact of permeant concentration and the choice of thermostat on the permeability. Moreover, we demonstrate that energetics of permeation can be closely examined by the dissipation of the gradient across the membrane to gain nuanced insights into the thermodynamics of permeability.

Atomistic Characterization of Beta-2-Glycoprotein I Domain V Interaction with Anionic Membranes

Background

Interaction of beta-2-glycoprotein I ( β 2 GPI) with anionic membranes is crucial in antiphospholipid syndrome (APS), implicating the role of it's membrane bind-ing domain, Domain V (DV). The mechanism of DV binding to anionic lipids is not fully understood.

Objectives

This study aims to elucidate the mechanism by which DV of β 2 GPI binds to anionic membranes.

Methods

We utilized molecular dynamics (MD) simulations to investigate the struc-tural basis of anionic lipid recognition by DV. To corroborate the membrane-binding mode identified in the HMMM simulations, we conducted additional simulations using a full mem-brane model.

Results

The study identified critical regions in DV, namely the lysine-rich loop and the hydrophobic loop, essential for membrane association via electrostatic and hydrophobic interactions, respectively. A novel lysine pair contributing to membrane binding was also discovered, providing new insights into β 2 GPI's membrane interaction. Simulations revealed two distinct binding modes of DV to the membrane, with mode 1 characterized by the insertion of the hydrophobic loop into the lipid bilayer, suggesting a dominant mechanism for membrane association. This interaction is pivotal for the pathogenesis of APS, as it facilitates the recognition of β 2 GPI by antiphospholipid antibodies.

Conclusion

The study advances our understanding of the molecular interactions be-tween β 2 GPI's DV and anionic membranes, crucial for APS pathogenesis. It highlights the importance of specific regions in DV for membrane binding and reveals a predominant bind-ing mode. These findings have significant implications for APS diagnostics and therapeutics, offering a deeper insight into the molecular basis of the syndrome.

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.

The carboxy terminus causes interfacial assembly of oleate hydratase on a membrane bilayer

The soluble flavoprotein oleate hydratase (OhyA) hydrates the 9-cis double bond of unsaturated fatty acids. OhyA substrates are embedded in membrane bilayers; OhyA must remove the fatty acid from the bilayer and enclose it in the active site. Here, we show that the positively charged helix-turn-helix motif in the carboxy terminus (CTD) is responsible for interacting with the negatively charged phosphatidylglycerol (PG) bilayer. Super-resolution microscopy of Staphylococcus aureus cells expressing green fluorescent protein fused to OhyA or the CTD sequence shows subcellular localization along the cellular boundary, indicating OhyA is membrane-associated and the CTD sequence is sufficient for membrane recruitment. Using cryo-electron microscopy, we solved the OhyA dimer structure and conducted 3D variability analysis of the reconstructions to assess CTD flexibility. Our surface plasmon resonance experiments corroborated that OhyA binds the PG bilayer with nanomolar affinity and we found the CTD sequence has intrinsic PG binding properties. We determined that the nuclear magnetic resonance structure of a peptide containing the CTD sequence resembles the OhyA crystal structure. We observed intermolecular NOE from PG liposome protons next to the phosphate group to the CTD peptide. The addition of paramagnetic MnCl2 indicated the CTD peptide binds the PG surface but does not insert into the bilayer. Molecular dynamics simulations, supported by site-directed mutagenesis experiments, identify key residues in the helix-turn-helix that drive membrane association. The data show that the OhyA CTD binds the phosphate layer of the PG surface to obtain bilayer-embedded unsaturated fatty acids.

A Marcus-Type Inverted Region in the Translocation Kinetics of a Knotted Protein

Knotted proteins are rare but important species, yet how their complex topologies affect their physical properties is not fully understood. Here we combine single molecule nanopore experiments and all-atom MD simulations to study the electric-field-driven unfolding during the translocation through a model pore of individual protein knots important for methylating tRNA. One of these knots shows an unusual behavior that resembles the behavior of electrons hopping between two potential surfaces: as the electric potential driving the translocation reaction is increased, the rate eventually plateaus or slows back down in the "Marcus inverted regime". Our results shed light on the influence of topology in knotted proteins on their forced translocation through a pore connecting two electrostatic potential wells.

Asymmetric conformations and lipid interactions shape the ATP-coupled cycle of a heterodimeric ABC transporter

Here we used cryo-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy (DEER), and molecular dynamics (MD) simulations, to capture and characterize ATP- and substrate-bound inward-facing (IF) and occluded (OC) conformational states of the heterodimeric ATP binding cassette (ABC) multidrug exporter BmrCD in lipid nanodiscs. Supported by DEER analysis, the structures reveal that ATP-powered isomerization entails changes in the relative symmetry of the BmrC and BmrD subunits that propagates from the transmembrane domain to the nucleotide binding domain. The structures uncover asymmetric substrate and Mg2+ binding which we hypothesize are required for triggering ATP hydrolysis preferentially in one of the nucleotide-binding sites. MD simulations demonstrate that multiple lipid molecules differentially bind the IF versus the OC conformation thus establishing that lipid interactions modulate BmrCD energy landscape. Our findings are framed in a model that highlights the role of asymmetric conformations in the ATP-coupled transport with general implications to the mechanism of ABC transporters.

Ionizable Amino Lipids Distribution and Effects on DSPC/Cholesterol Membranes: Implications for Lipid Nanoparticle Structure

Lipid nanoparticles (LNPs) containing ionizable aminolipids are among the leading platforms for the successful delivery of nucleic-acid-based therapeutics, including messenger RNA (mRNA). The two recently FDA-approved COVID-19 vaccines developed by Moderna and Pfizer/BioNTech belong to this category. Ionizable aminolipids, cholesterol, and DSPC lipids are among the key components of such formulations, crucially modulating physicochemical properties of these formulations and, consequently, the potency of these therapeutics. Despite the importance of these components, the distribution of these molecules in LNPs containing mRNA is not clear. In this study, we used all-atom molecular dynamics (MD) simulations to investigate the distribution and effects of the Lipid-5 (apparent pKa of the lipid nanoparticle = 6.56), a rationally designed and previously reported ionizable aminolipid by Moderna, on lipid bilayers [Mol. Ther. 2018, 26, 1509-1519]. The simulations were conducted with half of the aminolipids charged and half neutral approximately to the expected ionization in the microenvironment of the LNP surface. In all five simulated systems in this work, the cholesterol content was kept constant, whereas the DSPC and Lipid-5 concentrations were changed systematically. We found that at higher concentrations of the ionizable aminolipids, the neutral aminolipids form a disordered aggregate in the membrane interior that preferentially includes cholesterol. The rules underlying the lipid redistribution could be used to rationally choose lipids to optimize the LNP function.

Asymmetric conformations and lipid interactions shape the ATP-coupled cycle of a heterodimeric ABC transporter

To illuminate the structural origin of catalytic asymmetry of heterodimeric ABC transporters and how it shapes the energetics of their conformational cycles, we used cryo-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy (DEER), and molecular dynamics (MD) simulations, to capture and characterize conformational states of the heterodimeric ABC multidrug exporter BmrCD in lipid nanodiscs. In addition to multiple ATP- and substrate-bound inward-facing (IF) conformations, we obtained the structure of an occluded (OC) conformation wherein the unique extracellular domain (ECD) twists to partially open the extracellular gate. In conjunction with DEER analysis of the populations of these conformations, the structures reveal that ATP-powered isomerization entails changes in the relative symmetry of the BmrC and BmrD subunits that propagates from the transmembrane domain (TMD) to the nucleotide binding domain (NBD). The structures uncover asymmetric substrate and Mg 2+ binding which we hypothesize are required for triggering ATP hydrolysis preferentially in one of the nucleotide-binding sites. MD simulations demonstrated that multiple lipid molecules, identified from the cryo-EM density maps, differentially bind the IF versus the OC conformation thus modulating their relative stability. In addition to establishing how lipid interactions with BmrCD modulate the energy landscape, our findings are framed in a distinct transport model that highlights the role of asymmetric conformations in the ATP-coupled cycle with implications to the mechanism of ABC transporters in general.

Targeting lipid-protein interaction to treat Syk-mediated acute myeloid leukemia

Membrane lipids control the cellular activity of kinases containing the Src homology 2 (SH2) domain through direct lipid-SH2 domain interactions. Here we report development of new nonlipidic small molecule inhibitors of the lipid-SH2 domain interaction that block the cellular activity of their host proteins. As a pilot study, we evaluated the efficacy of lipid-SH2 domain interaction inhibitors for spleen tyrosine kinase (Syk), which is implicated in hematopoietic malignancies, including acute myeloid leukemia (AML). An optimized inhibitor (WC36) specifically and potently suppressed oncogenic activities of Syk in AML cell lines and patient-derived AML cells. Unlike ATP-competitive Syk inhibitors, WC36 was refractory to de novo and acquired drug resistance due to its ability to block not only the Syk kinase activity, but also its noncatalytic scaffolding function that is linked to drug resistance. Collectively, our study shows that targeting lipid-protein interaction is a powerful approach to developing new small molecule drugs.

LongBondEliminator: A Molecular Simulation Tool to Remove Ring Penetrations in Biomolecular Simulation Systems

We develop a workflow, implemented as a plugin to the molecular visualization program VMD, that can fix ring penetrations with minimal user input. LongBondEliminator, detects ring piercing artifacts by the long, strained bonds that are the local minimum energy conformation during minimization for some assembled simulation system. The LongBondEliminator tool then automatically treats regions near these long bonds using multiple biases applied through NAMD. By combining biases implemented through the collective variables module, density-based forces, and alchemical techniques in NAMD, LongBondEliminator will iteratively alleviate long bonds found within molecular simulation systems. Through three concrete examples with increasing complexity, a lignin polymer, an viral capsid assembly, and a large, highly glycosylated protein aggrecan, we demonstrate the utility for this method in eliminating ring penetrations from classical MD simulation systems. The tool is available via gitlab as a VMD plugin, and has been developed to be generically useful across a variety of biomolecular simulations.

Harris S, Ramanathan A. Intelligent resolution: Integrating Cryo-EM with AI-driven multi-resolution simulations to observe the severe acute respiratory syndrome coronavirus-2 replication-transcription machinery in action

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) replication transcription complex (RTC) is a multi-domain protein responsible for replicating and transcribing the viral mRNA inside a human cell. Attacking RTC function with pharmaceutical compounds is a pathway to treating COVID-19. Conventional tools, e.g., cryo-electron microscopy and all-atom molecular dynamics (AAMD), do not provide sufficiently high resolution or timescale to capture important dynamics of this molecular machine. Consequently, we develop an innovative workflow that bridges the gap between these resolutions, using mesoscale fluctuating finite element analysis (FFEA) continuum simulations and a hierarchy of AI-methods that continually learn and infer features for maintaining consistency between AAMD and FFEA simulations. We leverage a multi-site distributed workflow manager to orchestrate AI, FFEA, and AAMD jobs, providing optimal resource utilization across HPC centers. Our study provides unprecedented access to study the SARS-CoV-2 RTC machinery, while providing general capability for AI-enabled multi-resolution simulations at scale.

Distinguishing between Similar Miniproteins with Single-Molecule Nanopore Sensing: A Computational Study

A nanopore is a tool in single-molecule sensing biotechnology that offers label-free identification with high throughput. Nanopores have been successfully applied to sequence DNA and show potential in the study of proteins. Nevertheless, the task remains challenging due to the large variability in size, charges, and folds of proteins. Miniproteins have a small number of residues, limited secondary structure, and stable tertiary structure, which can offer a systematic way to reduce complexity. In this computational work, we theoretically evaluated sensing two miniproteins found in the human body using a silicon nitride nanopore. We employed molecular dynamics methods to compute occupied-pore ionic current magnitudes and electronic structure calculations to obtain interaction strengths between pore wall and miniprotein. From the interaction strength, we derived dwell times using a mix of combinatorics and numerical solutions. This latter approach circumvents typical computational demands needed to simulate translocation events using molecular dynamics. We focused on two miniproteins potentially difficult to distinguish owing to their isotropic geometry, similar number of residues, and overall comparable structure. We found that the occupied-pore current magnitudes not to vary significantly, but their dwell times differ by 1 order of magnitude. Together, these results suggest a successful identification protocol for similar miniproteins.

Organic contaminants and atmospheric nitrogen at the graphene-water interface: a simulation study

Ordered nanoscale patterns have been observed by atomic force microscopy at graphene-water and graphite-water interfaces. The two dominant explanations for these patterns are that (i) they consist of self-assembled organic contaminants or (ii) they are dense layers formed from atmospheric gases (especially nitrogen). Here we apply molecular dynamics simulations to study the behavior of dinitrogen and possible organic contaminants at the graphene-water interface. Despite the high concentration of N₂ in ambient air, we find that its expected occupancy at the graphene-water interface is quite low. Although dense (disordered) aggregates of dinitrogen have been observed in previous simulations, our results suggest that they are stable only in the presence of supersaturated aqueous N₂ solutions and dissipate rapidly when they coexist with nitrogen gas near atmospheric pressure. On the other hand, although heavy alkanes are present at only trace concentrations (micrograms per cubic meter) in typical indoor air, we predict that such concentrations can be sufficient to form ordered monolayers that cover the graphene-water interface. For octadecane, grand canonical Monte Carlo suggests nucleation and growth of monolayers above an ambient concentration near 6 μg m^-3, which is less than some literature values for indoor air. The thermodynamics of the formation of these alkane monolayers includes contributions from the hydration free-energy (unfavorable), the free-energy of adsorption to the graphene-water interface (highly favorable), and integration into the alkane monolayer phase (highly favorable). Furthermore, the peak-to-peak distances in AFM force profiles perpendicular to the interface (0.43-0.53 nm), agree with the distances calculated in simulations for overlayers of alkane-like molecules, but not for molecules such as N₂, water, or aromatics. Taken together, these results suggest that ordered domains observed on graphene, graphite, and other hydrophobic materials in water are consistent with alkane-like molecules occupying the interface.

Membrane Mixer: A Toolkit for Efficient Shuffling of Lipids in Heterogeneous Biological Membranes

Molecular dynamics (MD) simulations of biological membranes have achieved such levels of sophistication that are commonly used to predict unresolved structures and various properties of lipids and to substantiate experimental data. While achieving sufficient sampling of lipid dynamics remains a major challenge, a commonly used method to improve lipid sampling, e.g., in terms of specific interactions with membrane-associated proteins, is to randomize the initial arrangement of lipid constituents in multiple replicas of simulations, without changing the overall lipid composition of the membrane of interest. Here, we introduce a method that can rapidly generate multiple replicas of lipid bilayers with different spatial and conformational configurations for any given lipid composition. The underlying algorithm, which allows one to shuffle lipids at any desired level, relies on the application of an external potential, here referred to as the "carving potential", that removes clashes/entanglements before lipid positions are exchanged (shuffled), thereby minimizing the energy penalty due to abrupt lipid repositioning. The method is implemented as "Membrane Mixer Plugin (MMP) 1.0" in VMD, with a convenient graphical user interface that guides the user in setting various options and parameters. The plugin is fully automated and generates new membrane replicas more rapidly and conveniently than other analogous tools. The plugin and its capabilities introduced here can be extended to include additional features in future versions.

Principles of Small-Molecule Transport through Synthetic Nanopores

Synthetic nanopores made from DNA replicate the key biological processes of transporting molecular cargo across lipid bilayers. Understanding transport across the confined lumen of the nanopores is of fundamental interest and of relevance to their rational design for biotechnological applications. Here we reveal the transport principles of organic molecules through DNA nanopores by synergistically combining experiments and computer simulations. Using a highly parallel nanostructured platform, we synchronously measure the kinetic flux across hundreds of individual pores to obtain rate constants. The single-channel transport kinetics are close to the theoretical maximum, while selectivity is determined by the interplay of cargo charge and size, the pores' sterics and electrostatics, and the composition of the surrounding lipid bilayer. The narrow distribution of transport rates implies a high structural homogeneity of DNA nanopores. The molecular passageway through the nanopore is elucidated via coarse-grained constant-velocity steered molecular dynamics simulations. The ensemble simulations pinpoint with high resolution and statistical validity the selectivity filter within the channel lumen and determine the energetic factors governing transport. Our findings on these synthetic pores' structure-function relationship will serve to guide their rational engineering to tailor transport selectivity for cell biological research, sensing, and drug delivery.

Molecular mechanism of capsid disassembly in hepatitis B virus

The disassembly of a viral capsid leading to the release of its genetic material into the host cell is a fundamental step in viral infection. In hepatitis B virus (HBV), the capsid consists of identical protein monomers that dimerize and then arrange themselves into pentamers or hexamers on the capsid surface. By applying atomistic molecular dynamics simulation to an entire solvated HBV capsid subjected to a uniform mechanical stress protocol, we monitor the capsid-disassembly process and analyze the process down to the level of individual amino acids in 20 independent simulation replicas. The strain of an isotropic external force, combined with structural fluctuations, causes structurally heterogeneous cracks to appear in the HBV capsid. Analysis of the monomer-monomer interfaces reveals that, in contrast to the expectation from purely mechanical considerations, the cracks mainly occur within hexameric sites, whereas pentameric sites remain largely intact. Only a small subset of the capsid protein monomers, different in each simulation, are engaged in each instance of disassembly. We identify specific residues whose interactions are most readily lost during disassembly; R127, I139, Y132, N136, A137, and V149 are among the hot spots at the interfaces between dimers that lie within hexamers, leading to disassembly. The majority of these hot-spot residues are conserved by evolution, hinting to their importance for disassembly by avoiding overstabilization of capsids.

Large-Scale Molecular Dynamics Simulations of Cellular Compartments

Molecular dynamics or MD simulation is gradually maturing into a tool for constructing in vivo models of living cells in atomistic details. The feasibility of such models is bolstered by integrating the simulations with data from microscopic, tomographic and spectroscopic experiments on exascale supercomputers, facilitated by the use of deep learning technologies. Over time, MD simulation has evolved from tens of thousands of atoms to over 100 million atoms comprising an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium. In this chapter, we present a step-by-step outline for preparing, executing and analyzing such large-scale MD simulations of biological systems that are essential to life processes. All scripts are provided via GitHub.

Free energy and kinetics of cAMP permeation through connexin26 via applied voltage and milestoning

The connexin family is a diverse group of highly regulated wide-pore channels permeable to biological signaling molecules. Despite the critical roles of connexins in mediating selective molecular signaling in health and disease, the basis of molecular permeation through these pores remains unclear. Here, we report the thermodynamics and kinetics of binding and transport of a second messenger, adenosine-3',5'-cyclophosphate (cAMP), through a connexin26 hemichannel (Cx26). First, inward and outward fluxes of cAMP molecules solvated in KCl solution were obtained from 4 μs of ± 200 mV simulations. These fluxes data yielded a single-channel permeability of cAMP and cAMP/K⁺ permeability ratio consistent with experimentally measured values. The results from voltage simulations were then compared with the potential of mean force (PMF) and the mean first passage times (MFPTs) of a single cAMP without voltage, obtained from a total of 16.5 μs of Voronoi-tessellated Markovian milestoning simulations. Both the voltage simulations and the milestoning simulations revealed two cAMP-binding sites, for which the binding constants K_D and dissociation rates k_off were computed from PMF and MFPTs. The protein dipole inside the pore produces an asymmetric PMF, reflected in unequal cAMP MFPTs in each direction once within the pore. The free energy profiles under opposite voltages were derived from the milestoning PMF and revealed the interplay between voltage and channel polarity on the total free energy. In addition, we show how these factors influence the cAMP dipole vector during permeation, and how cAMP affects the local and nonlocal pore diameter in a position-dependent manner.

Protein Transport through Nanopores Illuminated by Long-Time-Scale Simulations

The transport of molecules through nanoscale confined space is relevant in biology, biosensing, and industrial filtration. Microscopically modeling transport through nanopores is required for a fundamental understanding and guiding engineering, but the short duration and low replica number of existing simulation approaches limit statistically relevant insight. Here we explore protein transport in nanopores with a high-throughput computational method that realistically simulates hundreds of up to seconds-long protein trajectories by combining Brownian dynamics and continuum simulation and integrating both driving forces of electroosmosis and electrophoresis. Ionic current traces are computed to enable experimental comparison. By examining three biological and synthetic nanopores, our study answers questions about the kinetics and mechanism of protein transport and additionally reveals insight that is inaccessible from experiments yet relevant for pore design. The discovery of extremely frequent unhindered passage can guide the improvement of biosensor pores to enhance desired biomolecular recognition by pore-tethered receptors. Similarly, experimentally invisible nontarget adsorption to pore walls highlights how to improve recently developed DNA nanopores. Our work can be expanded to pressure-driven flow to model industrial nanofiltration processes.

Structural-profiling of low molecular weight RNAs by nanopore trapping/translocation using Mycobacterium smegmatis porin A

Folding of RNA can produce elaborate tertiary structures, corresponding to their diverse roles in the regulation of biological activities. Direct observation of RNA structures at high resolution in their native form however remains a challenge. The large vestibule and the narrow constriction of a Mycobacterium smegmatis porin A (MspA) suggests a sensing mode called nanopore trapping/translocation, which clearly distinguishes between microRNA, small interfering RNA (siRNA), transfer RNA (tRNA) and 5 S ribosomal RNA (rRNA). To further profit from the acquired event characteristics, a custom machine learning algorithm is developed. Events from measurements with a mixture of RNA analytes can be automatically classified, reporting a general accuracy of ~93.4%. tRNAs, which possess a unique tertiary structure, report a highly distinguishable sensing feature, different from all other RNA types tested in this study. With this strategy, tRNAs from different sources are measured and a high structural conservation across different species is observed in single molecule.

Electrical unfolding of cytochrome c during translocation through a nanopore constriction

Many small proteins move across cellular compartments through narrow pores. In order to thread a protein through a constriction, free energy must be overcome to either deform or completely unfold the protein. In principle, the diameter of the pore, along with the effective driving force for unfolding the protein, as well as its barrier to translocation, should be critical factors that govern whether the process proceeds via squeezing, unfolding/threading, or both. To probe this for a well-established protein system, we studied the electric-field-driven translocation behavior of cytochrome c (cyt c) through ultrathin silicon nitride (SiN_x) solid-state nanopores of diameters ranging from 1.5 to 5.5 nm. For a 2.5-nm-diameter pore, we find that, in a threshold electric-field regime of ∼30 to 100 MV/m, cyt c is able to squeeze through the pore. As electric fields inside the pore are increased, the unfolded state of cyt c is thermodynamically stabilized, facilitating its translocation. In contrast, for 1.5- and 2.0-nm-diameter pores, translocation occurs only by threading of the fully unfolded protein after it transitions through a higher energy unfolding intermediate state at the mouth of the pore. The relative energies between the metastable, intermediate, and unfolded protein states are extracted using a simple thermodynamic model that is dictated by the relatively slow (∼ms) protein translocation times for passing through the nanopore. These experiments map the various modes of protein translocation through a constriction, which opens avenues for exploring protein folding structures, internal contacts, and electric-field-induced deformability.

Binding mode of SARS-CoV-2 fusion peptide to human cellular membrane

Infection of human cells by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) relies on its binding to a specific receptor and subsequent fusion of the viral and host cell membranes. The fusion peptide (FP), a short peptide segment in the spike protein, plays a central role in the initial penetration of the virus into the host cell membrane, followed by the fusion of the two membranes. Here, we use an array of molecular dynamics simulations that take advantage of the highly mobile membrane mimetic model to investigate the interaction of the SARS-CoV2 FP with a lipid bilayer representing mammalian cellular membranes at an atomic level and to characterize the membrane-bound form of the peptide. Six independent systems were generated by changing the initial positioning and orientation of the FP with respect to the membrane, and each system was simulated in five independent replicas, each for 300 ns. In 73% of the simulations, the FP reaches a stable, membrane-bound configuration, in which the peptide deeply penetrated into the membrane. Clustering of the results reveals three major membrane-binding modes (binding modes 1-3), in which binding mode 1 populates over half of the data points. Taking into account the sequence conservation among the viral FPs and the results of mutagenesis studies establishing the role of specific residues in the helical portion of the FP in membrane association, the significant depth of penetration of the whole peptide, and the dense population of the respective cluster, we propose that the most deeply inserted membrane-bound form (binding mode 1) represents more closely the biologically relevant form. Analysis of FP-lipid interactions shows the involvement of specific residues, previously described as the "fusion-active core residues," in membrane binding. Taken together, the results shed light on a key step involved in SARS-CoV2 infection, with potential implications in designing novel inhibitors.

Mechanical Activation of MscL Revealed by a Locally Distributed Tension Molecular Dynamics Approach

Membrane tension perceived by mechanosensitive (MS) proteins mediates cellular responses to mechanical stimuli and osmotic stresses, and it also guides multiple biological functions including cardiovascular control and development. In bacteria, MS channels function as tension-activated pores limiting excessive turgor pressure, with MS channel of large conductance (MscL) acting as an emergency release valve preventing cell lysis. Previous attempts to simulate gating transitions in MscL by either directly applying steering forces to the protein or by increasing the whole-system tension were not fully successful and often disrupted the integrity of the system. We present a novel, to our knowledge, locally distributed tension molecular dynamics (LDT-MD) simulation method that allows application of forces continuously distributed among lipids surrounding the channel using a specially constructed collective variable. We report reproducible and reversible transitions of MscL to the open state with measured parameters of lateral expansion and conductivity that exactly satisfy experimental values. The LDT-MD method enables exploration of the MscL-gating process with different pulling velocities and variable tension asymmetry between the inner and outer membrane leaflets. We use LDT-MD in combination with well-tempered metadynamics to reconstruct the tension-dependent free-energy landscape for the opening transition in MscL. The flexible definition of the LDT collective variable allows general application of our method to study mechanical activation of any membrane-embedded protein.

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. Here, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMDis 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 biomolecular 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 alchemical 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.

Polynucleotide differentiation using hybrid solid-state nanopore functionalizing with α-hemolysin

We report results from full atomistic molecular dynamics simulations on the properties of biomimetic nanopores. This latter result was obtained through the direct insertion of an α-hemolysin protein inside a hydrophobic solid-state nanopore. Upon translocation of different DNA strands, we demonstrate here that the theoretical system presents the same discrimination properties as the experimental one obtained previously. This opens an interesting way to promote the stability of a specific protein inside a solid nanopore to develop further biomimetic applications for DNA or protein sequencing.

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.

Direct Derivation of Free Energies of Membrane Deformation and Other Solvent Density Variations From Enhanced Sampling Molecular Dynamics

We report a methodology to calculate the free energy of a shape transformation in a lipid membrane directly from a molecular dynamics simulation. The bilayer need not be homogeneous or symmetric and can be atomically detailed or coarse grained. The method is based on a collective variable that quantifies the similarity between the membrane and a set of predefined density distributions. Enhanced sampling of this "Multi-Map" variable re-shapes the bilayer and permits the derivation of the corresponding potential of mean force. Calculated energies thus reflect the dynamic interplay of atoms and molecules, rather than postulated effects. Evaluation of deformations of different shape, amplitude, and range demonstrates that the macroscopic bending modulus assumed by the Helfrich-Canham model is increasingly unsuitable below the 100-Å scale. In this range of major biological significance, direct free-energy calculations reveal a much greater plasticity. We also quantify the stiffening effect of cholesterol on bilayers of different composition and compare with experiments. Lastly, we illustrate how this approach facilitates analysis of other solvent reorganization processes, such as hydrophobic hydration. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Insights into protein sequencing with an α-Hemolysin nanopore by atomistic simulations

Single molecule protein sequencing would represent a disruptive burst in proteomic research with important biomedical impacts. Due to their success in DNA sequencing, nanopore based devices have been recently proposed as possible tools for the sequencing of peptide chains. One of the open questions in nanopore protein sequencing concerns the ability of such devices to provide different signals for all the 20 standard amino acids. Here, using equilibrium all-atom molecular dynamics simulations, we estimated the pore clogging in α-Hemolysin nanopore associated to 20 different homopeptides, one for each standard amino acid. Our results show that pore clogging is affected by amino acid volume, hydrophobicity and net charge. The equilibrium estimations are also supported by non-equilibrium runs for calculating the current blockades for selected homopeptides. Finally, we discuss the possibility to modify the α-Hemolysin nanopore, cutting a portion of the barrel region close to the trans side, to reduce spurious signals and, hence, to enhance the sensitivity of the nanopore.

Sequence-Dependent Unzipping Dynamics of DNA Hairpins in a Nanopore

By applying an electric field to an insulating membrane, movement of charged particles through a nanopore can be induced. The measured ionic current reports on biomolecules passing through the nanopore. In this paper, we explore the sequence-dependent dynamics of DNA unzipping using our recently developed coarse-grained model. We estimated three molecular profiles (the potential of mean force, position-dependent diffusion coefficient, and position-dependent effective charge) for the DNA unzipping of four hairpins with different sequences. We found that the molecular profiles are correlated with the ionic current and molecular events. We also explored the unzipping kinetics using Brownian dynamics. We found that the effect of hairpin structure on the unzipping/translocation times is not only energetic (weaker hairpins unzip more quickly) but also kinetic (different unzipping and translocation pathways play an important role).

Power transformations improve interpolation of grids for molecular mechanics interaction energies

A common strategy for speeding up molecular docking calculations is to precompute nonbonded interaction energies between a receptor molecule and a set of three-dimensional grids. The grids are then interpolated to compute energies for ligand atoms in many different binding poses. Here, I evaluate a smoothing strategy of taking a power transformation of grid point energies and inverse transformation of the result from trilinear interpolation. For molecular docking poses from 85 protein-ligand complexes, this smoothing procedure leads to significant accuracy improvements, including an approximately twofold reduction in the root mean square error at a grid spacing of 0.4 Å and retaining the ability to rank docking poses even at a grid spacing of 0.7 Å. © 2018 Wiley Periodicals, Inc.

Protein sequencing via nanopore based devices: a nanofluidics perspective

Proteins perform a huge number of central functions in living organisms, thus all the new techniques allowing their precise, fast and accurate characterization at single-molecule level certainly represent a burst in proteomics with important biomedical impact. In this review, we describe the recent progresses in the developing of nanopore based devices for protein sequencing. We start with a critical analysis of the main technical requirements for nanopore protein sequencing, summarizing some ideas and methodologies that have recently appeared in the literature. In the last sections, we focus on the physical modelling of the transport phenomena occurring in nanopore based devices. The multiscale nature of the problem is discussed and, in this respect, some of the main possible computational approaches are illustrated.

Single-molecule techniques in biophysics: a review of the progress in methods and applications

Single-molecule biophysics has transformed our understanding of biology, but also of the physics of life. More exotic than simple soft matter, biomatter lives far from thermal equilibrium, covering multiple lengths from the nanoscale of single molecules to up to several orders of magnitude higher in cells, tissues and organisms. Biomolecules are often characterized by underlying instability: multiple metastable free energy states exist, separated by levels of just a few multiples of the thermal energy scale k _B T, where k _B is the Boltzmann constant and T absolute temperature, implying complex inter-conversion kinetics in the relatively hot, wet environment of active biological matter. A key benefit of single-molecule biophysics techniques is their ability to probe heterogeneity of free energy states across a molecular population, too challenging in general for conventional ensemble average approaches. Parallel developments in experimental and computational techniques have catalysed the birth of multiplexed, correlative techniques to tackle previously intractable biological questions. Experimentally, progress has been driven by improvements in sensitivity and speed of detectors, and the stability and efficiency of light sources, probes and microfluidics. We discuss the motivation and requirements for these recent experiments, including the underpinning mathematics. These methods are broadly divided into tools which detect molecules and those which manipulate them. For the former we discuss the progress of super-resolution microscopy, transformative for addressing many longstanding questions in the life sciences, and for the latter we include progress in 'force spectroscopy' techniques that mechanically perturb molecules. We also consider in silico progress of single-molecule computational physics, and how simulation and experimentation may be drawn together to give a more complete understanding. Increasingly, combinatorial techniques are now used, including correlative atomic force microscopy and fluorescence imaging, to probe questions closer to native physiological behaviour. We identify the trade-offs, limitations and applications of these techniques, and discuss exciting new directions.

Iodide Binding in Sodium-Coupled Cotransporters

Several apical iodide translocation pathways have been proposed for iodide efflux out of thyroid follicular cells, including a pathway mediated by the sodium-coupled monocarboxylate transporter 1 (SMCT1), which remains controversial. Herein, we evaluate structural and functional similarities between SMCT1 and the well-studied sodium-iodide symporter (NIS) that mediates the first step of iodide entry into the thyroid. Free-energy calculations using a force field with electronic polarizability verify the presence of a conserved iodide-binding pocket between the TM2, TM3, and TM7 segments in hNIS, where iodide is coordinated by Phe67, Gln72, Cys91, and Gln94. We demonstrate the mutation of residue Gly93 of hNIS to a larger amino acid expels the side chain of a critical tryptophan residue (Trp255) into the interior of the binding pocket, partially occluding the iodide binding site and reducing iodide affinity, which is consistent with previous reports associating mutation of this residue with iodide uptake deficiency and hypothyroidism. Furthermore, we find that the position of Trp255 in this hNIS mutant mirrors that of Trp253 in wild-type hSMCT1, where a threonine (Thr91) occupies the position homologous to that occupied by glycine in wild-type hNIS (Gly93). Correspondingly, mutation of Thr91 to glycine in hSMCT1 makes the pocket structure more like that of wild-type hNIS, increasing its iodide affinity. These results suggest that wild-type hSMCT1 in the inward-facing conformation may bind iodide only very weakly, which may have implications for its ability to transport iodide.

Chasing the open-state structure of pentameric ligand-gated ion channels

Members of the pentameric ligand-gated ion channel family have been crystallized in different conformations, including one in which the transmembrane pore is surprisingly wide. Gonzalez-Gutierrez et al. show that the open-channel conformation of animal members is more similar to the models with narrow pores.

Remarkable advances have been made toward the structural characterization of ion channels in the last two decades. However, the unambiguous assignment of well-defined functional states to the obtained structural models has proved challenging. In the case of the superfamily of nicotinic-receptor channels (also referred to as pentameric ligand-gated ion channels [pLGICs]), for example, two different types of model of the open-channel conformation have been proposed on the basis of structures solved to resolutions better than 4.0 Å. At the level of the transmembrane pore, the open-state models of the proton-gated pLGIC from Gloeobacter violaceus (GLIC) and the invertebrate glutamate-gated Cl^– channel (GluCl) are very similar to each other, but that of the glycine receptor (GlyR) is considerably wider. Indeed, the mean distances between the axis of ion permeation and the Cα atoms at the narrowest constriction of the pore (position −2′) differ by ∼2 Å in these two classes of model, a large difference when it comes to understanding the physicochemical bases of ion conduction and charge selectivity. Here, we take advantage of the extreme open-channel stabilizing effect of mutations at pore-facing position 9′. We find that the I9′A mutation slows down entry into desensitization of GLIC to the extent that macroscopic currents decay only slightly by the end of pH 4.5 solution applications to the extracellular side for several minutes. We crystallize (at pH 4.5) two variants of GLIC carrying this mutation and solve their structures to resolutions of 3.12 Å and 3.36 Å. Furthermore, we perform all-atom molecular dynamics simulations of ion permeation and picrotoxinin block, using the different open-channel structural models. On the basis of these results, we favor the notion that the open-channel structure of pLGICs from animals is much closer to that of the narrow models (of GLIC and GluCl) than it is to that of the GlyR.

SDS-assisted protein transport through solid-state nanopores

Using nanopores for single-molecule sequencing of proteins - similar to nanopore-based sequencing of DNA - faces multiple challenges, including unfolding of the complex tertiary structure of the proteins and enforcing their unidirectional translocation through nanopores. Here, we combine molecular dynamics (MD) simulations with single-molecule experiments to investigate the utility of SDS (Sodium Dodecyl Sulfate) to unfold proteins for solid-state nanopore translocation, while simultaneously endowing them with a stronger electrical charge. Our simulations and experiments prove that SDS-treated proteins show a considerable loss of the protein structure during the nanopore translocation. Moreover, SDS-treated proteins translocate through the nanopore in the direction prescribed by the electrophoretic force due to the negative charge impaired by SDS. In summary, our results suggest that SDS causes protein unfolding while facilitating protein translocation in the direction of the electrophoretic force; both characteristics being advantageous for future protein sequencing applications using solid-state nanopores.

Activation and Desensitization Mechanism of AMPA Receptor-TARP Complex by Cryo-EM

AMPA receptors mediate fast excitatory neurotransmission in the mammalian brain and transduce the binding of presynaptically released glutamate to the opening of a transmembrane cation channel. Within the postsynaptic density, however, AMPA receptors coassemble with transmembrane AMPA receptor regulatory proteins (TARPs), yielding a receptor complex with altered gating kinetics, pharmacology, and pore properties. Here, we elucidate structures of the GluA2-TARP γ2 complex in the presence of the partial agonist kainate or the full agonist quisqualate together with a positive allosteric modulator or with quisqualate alone. We show how TARPs sculpt the ligand-binding domain gating ring, enhancing kainate potency and diminishing the ensemble of desensitized states. TARPs encircle the receptor ion channel, stabilizing M2 helices and pore loops, illustrating how TARPs alter receptor pore properties. Structural and computational analysis suggests the full agonist and modulator complex harbors an ion-permeable channel gate, providing the first view of an activated AMPA receptor.

Modulation of Molecular Flux Using a Graphene Nanopore Capacitor

Modulation of ionic current flowing through nanoscale pores is one of the fundamental biological processes. Inspired by nature, nanopores in synthetic solid-state membranes are being developed to enable rapid analysis of biological macromolecules and to serve as elements of nanofludic circuits. Here, we theoretically investigate ion and water transport through a graphene-insulator-graphene membrane containing a single, electrolyte-filled nanopore. By means of all-atom molecular dynamics simulations, we show that the charge state of such a graphene nanopore capacitor can regulate both the selectivity and the magnitude of the nanopore ionic current. At a fixed transmembrane bias, the ionic current can be switched from being carried by an equal mixture of cations and anions to being carried almost exclusively by either cationic or anionic species, depending on the sign of the charge assigned to both plates of the capacitor. Assigning the plates of the capacitor opposite sign charges can either increase the nanopore current or reduce it substantially, depending on the polarity of the bias driving the transmembrane current. Facilitated by the changes of the nanopore surface charge, such ionic current modulations are found to occur despite the physical dimensions of the nanopore being an order of magnitude larger than the screening length of the electrolyte. The ionic current rectification is accompanied by a pronounced electro-osmotic effect that can transport neutral molecules such as proteins and drugs across the solid-state membrane and thereby serve as an interface between electronic and chemical signals.

Mechanical Trapping of DNA in a Double-Nanopore System

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.

Improved model of hydrated calcium ion for molecular dynamics simulations using classical biomolecular force fields

Calcium ions (Ca(2+) ) play key roles in various fundamental biological processes such as cell signaling and brain function. Molecular dynamics (MD) simulations have been used to study such interactions, however, the accuracy of the Ca(2+) models provided by the standard MD force fields has not been rigorously tested. Here, we assess the performance of the Ca(2+) models from the most popular classical force fields AMBER and CHARMM by computing the osmotic pressure of model compounds and the free energy of DNA-DNA interactions. In the simulations performed using the two standard models, Ca(2+) ions are seen to form artificial clusters with chloride, acetate, and phosphate species; the osmotic pressure of CaAc2 and CaCl2 solutions is a small fraction of the experimental values for both force fields. Using the standard parameterization of Ca(2+) ions in the simulations of Ca(2+) -mediated DNA-DNA interactions leads to qualitatively wrong outcomes: both AMBER and CHARMM simulations suggest strong inter-DNA attraction whereas, in experiment, DNA molecules repel one another. The artificial attraction of Ca(2+) to DNA phosphate is strong enough to affect the direction of the electric field-driven translocation of DNA through a solid-state nanopore. To address these shortcomings of the standard Ca(2+) model, we introduce a custom model of a hydrated Ca(2+) ion and show that using our model brings the results of the above MD simulations in quantitative agreement with experiment. Our improved model of Ca(2+) can be readily applied to MD simulations of various biomolecular systems, including nucleic acids, proteins and lipid bilayer membranes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 752-763, 2016.

Free-energy calculations reveal the subtle differences in the interactions of DNA bases with α-hemolysin

Next generation DNA sequencing methods that utilize protein nanopores have the potential to revolutionize this area of biotechnology. While the technique is underpinned by simple physics, the wild-type protein pores do not have all of the desired properties for efficient and accurate DNA sequencing. Much of the research efforts have focused on protein nanopores, such as α-hemolysin from Staphylococcus aureus. However, the speed of DNA translocation has historically been an issue, hampered in part by incomplete knowledge of the energetics of translocation. Here we have utilized atomistic molecular dynamics simulations of nucleotide fragments in order to calculate the potential of mean force (PMF) through α-hemolysin. Our results reveal specific regions within the pore that play a key role in the interaction with DNA. In particular, charged residues such as D127 and K131 provide stabilizing interactions with the anionic DNA and therefore are likely to reduce the speed of translocation. These regions provide rational targets for pore optimization. Furthermore, we show that the energetic contributions to the protein-DNA interactions are a complex combination of electrostatics and short-range interactions, often mediated by water molecules.

What controls open-pore and residual currents in the first sensing zone of alpha-hemolysin nanopore? Combined experimental and theoretical study

The electrophoretic transport of single-stranded DNA through biological nanopores such as alpha-hemolysin (αHL) is a promising and cost-effective technology with the potential to revolutionize genomics. The rational design of pores with the controlled polymer translocation rates and high contrast between different nucleotides could improve significantly nanopore sequencing applications. Here, we apply a combination of theoretical and experimental methods in an attempt to elucidate several selective modifications in the pore which were proposed to be central for the effective discrimination between purines and pyrimidines. Our nanopore test set includes the wild type αHL and six mutants (E111N/M113X/K147N) in which the cross-section and chemical functionality of the first constriction zone of the pore are modified. Electrophysiological recordings were combined with all-atom Molecular Dynamics simulations (MD) and a recently developed Brownian Dynamics (BROMOC) protocol to investigate residual ion currents and pore-DNA interactions for two homo-polymers e.g. poly(dA)40 or poly(dC)40 blocking the pore. The calculated residual currents and contrast in the poly(dA)40/poly(dC)40 blocked pore are in qualitative agreement with the experimental recordings. We showed that a simple structural metric allows rationalization of key elements in the emergent contrast between purines and pyrimidines in the modified αHL mutants. The shape of the pore and its capacity for hydrogen bonding to a translocated polynucleotide are two essential parameters for contrast optimization. To further probe the impact of these two factors in the ssDNA sensing, we eliminated the effect of the primary constriction using serine substitutions (i.e. E111S/M113S/T145S/K147S) and increased the hydrophobic volume of the central residue in the secondary constriction (L135I). This pore modification sharply increased the contrast between Adenine (A) and Cytosine (C).

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling

Molecular Dynamics Flexible Fitting (MDFF) is an established technique for fitting all-atom structures of molecules into corresponding cryo-electron microscopy (cryo-EM) densities. The practical application of MDFF is simple but requires a user to be aware of and take measures against a variety of possible challenges presented by each individual case. Some of these challenges arise from the complexity of a molecular structure or the limited quality of available structural models and densities to be interpreted, while others stem from the intricacies of MDFF itself. The current article serves as an overview of the strategies that have been developed since MDFF's inception to overcome common challenges and successfully perform MDFF simulations.

Diffusive dynamics of DNA unzipping in a nanopore

When an electric field is applied to an insulating membrane, movement of charged particles through a nanopore is induced. The measured ionic current reports on biomolecules passing through the nanopore. In this work, we explored the kinetics of DNA unzipping in a nanopore using our coarse-grained model (Stachiewicz and Molski, J. Comput. Chem. 2015, 36, 947). Coarse graining allowed a more detailed analysis for a wider range of parameters than all-atom simulations. Dependence of the translocation mode (unzipping or distortion) on the pore diameter was examined, and the threshold voltages were estimated. We determined the potential of mean force, position-dependent diffusion coefficient, and position-dependent effective charge for the DNA unzipping. The three molecular profiles were correlated with the ionic current and molecular events. On the unzipping/translocation force profile, two energy maxima were found, one of them corresponding to the unzipping, and the other to the translocation barriers. The unzipping kinetics were further explored using Brownian dynamics.

Determination of Free Energy Profiles for the Translocation of Polynucleotides through α-Hemolysin Nanopores using Non-Equilibrium Molecular Dynamics Simulations

The translocation of polynucleotides through transmembrane protein pores is a fundamental biological process with important technological and medical relevance. The translocation process is complex, and it is influenced by a range of factors including the diameter and inner surface of the pore, the secondary structure of the polymer, and the interactions between the polymer and protein. In this paper, we perform nonequilibrium constant velocity-steered molecular dynamics simulations of nucleic acid molecule translocation through the protein nanopore α-hemolysin and use Jarzynski's identity to determine the associated free energy profiles. With this approach we are able to explain the observed differences in experimental translocation time through the nanopore between polyadenosine and polydeoxycytidine. The translocation of polynucleotides and single nucleotides through α-hemolysin is investigated. These simulations are computationally intensive as they employ models with atomistic level resolution; in addition to their size, these systems are challenging to study due to the time scales of translocation of large asymmetric molecules. Our simulations provide insight into the role of the interactions between the nucleic acid molecules and the protein pore. Mutated protein pores provide confirmation of residue-specific interactions between nucleotides and the protein pore. By harnessing such molecular dynamics simulations, we gain new physicochemical insight into the translocation process.