Effect of Interaction upon Translocation of Confined Polymer Chain Through Nanopore

The effect of the interaction between nanopore and chain monomer on the translocation of a single polymer chain confined in a finite size square through an interacting nanopore to a large space has been studied by two-dimensional bond fluctuation model with Monte Carlo simulation. Results indicate that the free energy barrier before the successful translocation of the chain depends linearly on the chain length as well as the nanopore length for different pore-polymer interaction, and the attractive interaction reduces the free energy barrier, leading to the reduction of the average trapping time.     

Comparative assessment of computational methods for the determination of solvation free energies in alcohol‐based molecules

The determination of differences in solvation free energies between related drug molecules remains an important challenge in computational drug optimization, when fast and accurate calculation of differences in binding free energy are required. In this study, we have evaluated the performance of five commonly used polarized continuum model (PCM) methodologies in the determination of solvation free energies for 53 typical alcohol and alkane small molecules. In addition, the performance of these PCM methods, of a thermodynamic integration (TI) protocol and of the Poisson–Boltzmann (PB) and generalized Born (GB) methods, were tested in the determination of solvation free energies changes for 28 common alkane‐alcohol transformations, by the substitution of an hydrogen atom for a hydroxyl substituent. The results show that the solvation model D (SMD) performs better among the PCM‐based approaches in estimating solvation free energies for alcohol molecules, and solvation free energy changes for alkane‐alcohol transformations, with an average error below 1 kcal/mol for both quantities. However, for the determination of solvation free energy changes on alkane‐alcohol transformation, PB and TI yielded better results. TI was particularly accurate in the treatment of hydroxyl groups additions to aromatic rings (0.53 kcal/mol), a common transformation when optimizing drug‐binding in computer‐aided drug design. © 2013 Wiley Periodicals, Inc.     

Discrimination of single‐stranded DNA homopolymers by sieving out G‐quadruplex using tiny solid‐state nanopores

Nanopore sensor has been developed as a promising technology for DNA sequencing at the single‐base resolution. However, the discrimination of homopolymers composed of guanines from other nucleotides has not been clearly revealed due to the easily formed G‐quadruplex in aqueous buffers. In this work, we report that a tiny silicon nitride nanopore was used to sieve out G tetramers to make sure only homopolymers composed of guanines could translocate through the nanopore, then the 20‐nucleotide long ssDNA homopolymers could be identified and differentiated. It is found that the size of the nucleotide plays a major role in affecting the current blockade as well as the dwell time while DNA is translocating through the nanopore. By the comparison of translocation behavior of ssDNA homopolymers composed of nucleotides with different volumes, it is found that smaller nucleotides can lead to higher translocation speed and lower current blockage, which is also found and validated for the 105‐nucleotide long homopolymers. The studies performed in this work will improve our understanding of nanopore‐based DNA sequencing at single‐base level.     

Accurate pKa Calculation of the Conjugate Acids of Alkanolamines,Alkaloids and Nucleotide Bases by Quantum Chemical Methods

The pK_a of the conjugate acids of alkanolamines, neurotransmitters, alkaloid drugs and nucleotide bases are calculated with density functional methods (B3LYP, M08‐HX and M11‐L) and ab initio methods (SCS‐MP2, G3). Implicit solvent effects are included with a conductor‐like polarizable continuum model (CPCM) and universal solvation models (SMD, SM8). G3, SCS‐MP2 and M11‐L methods coupled with SMD and SM8 solvation models perform well for alkanolamines with mean unsigned errors below 0.20 pK_a units, in all cases. Extending this method to the pK_a calculation of 35 nitrogen‐containing compounds spanning 12 pK_a units showed an excellent correlation between experimental and computational pK_a values of these 35 amines with the computationally low‐cost SM8/M11‐L density functional approach.     

Minimum free energy pathways and free energy profiles for conformational transitions based on atomistic molecular dynamics simulations

An efficient method for the calculation of minimum free energy pathways and free energy profiles for conformational transitions is presented. Short restricted perturbation-targeted molecular dynamics trajectories are used to generate an approximate free energy surface. Approximate reaction pathways for the conformational change are constructed from one-dimensional line segments on this surface using a Monte Carlo optimization. Accurate free energy profiles are then determined along the pathways by means of one-dimensional adaptive umbrella sampling simulations. The method is illustrated by its application to the alanine "dipeptide." Due to the low computational cost and memory demands, the method is expected to be useful for the treatment of large biomolecular systems.     

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. © 2015 Wiley Periodicals, Inc.     

Translocation of single-stranded DNA through the α-hemolysin protein nanopore in acidic solutions

The effect of acidic pH on the translocation of single-stranded DNA through the α-hemolysin pore is investigated. Two significantly different types of events, i.e. deep blockades and shallow blockades, are observed at low pH. The residence times of the shallow blockades are not significantly different from those of the DNA translocation events obtained at or near physiological pH, whereas the deep blockades have much larger residence times and blockage amplitudes. With a decrease in the pH of the electrolyte solution, the percentage of the deep blockades in the total events increases. Furthermore, the mean residence time of these long-lived events is dependent on the length of DNA, and also varies with the nucleotide base, suggesting that they are appropriate for use in DNA analysis. In addition to being used as an effective approach to affect DNA translocation in the nanopore, manipulation of the pH of the electrolyte solution provides a potential means to greatly enhance the sensitivity of nanopore stochastic sensing.     

Tinker‐OpenMM: Absolute and relative alchemical free energies using AMOEBA on GPUs

The capabilities of the polarizable force fields for alchemical free energy calculations have been limited by the high computational cost and complexity of the underlying potential energy functions. In this work, we present a GPU‐based general alchemical free energy simulation platform for polarizable potential AMOEBA. Tinker‐OpenMM, the OpenMM implementation of the AMOEBA simulation engine has been modified to enable both absolute and relative alchemical simulations on GPUs, which leads to a ∼200‐fold improvement in simulation speed over a single CPU core. We show that free energy values calculated using this platform agree with the results of Tinker simulations for the hydration of organic compounds and binding of host–guest systems within the statistical errors. In addition to absolute binding, we designed a relative alchemical approach for computing relative binding affinities of ligands to the same host, where a special path was applied to avoid numerical instability due to polarization between the different ligands that bind to the same site. This scheme is general and does not require ligands to have similar scaffolds. We show that relative hydration and binding free energy calculated using this approach match those computed from the absolute free energy approach. © 2017 Wiley Periodicals, Inc.     

Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations

Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

String method in collective variables: minimum free energy paths and isocommittor surfaces

A computational technique is proposed which combines the string method with a sampling technique to determine minimum free energy paths. The technique only requires to compute the mean force and another conditional expectation locally along the string, and therefore can be applied even if the number of collective variables kept in the free energy calculation is large. This is in contrast with other free energy sampling techniques which aim at mapping the full free energy landscape and whose cost increases exponentially with the number of collective variables kept in the free energy. Provided that the number of collective variables is large enough, the new technique captures the mechanism of transition in that it allows to determine the committor function for the reaction and, in particular, the transition state region. The new technique is illustrated on the example of alanine dipeptide, in which we compute the minimum free energy path for the isomerization transition using either two or four dihedral angles as collective variables. It is shown that the mechanism of transition can be captured using the four dihedral angles, but it cannot be captured using only two of them.     

Spatial blockage of ionic current for electrophoretic translocation of DNA through a graphene nanopore

Graphene nanopore has been promising the ultra‐high resolution for DNA sequencing due to the atomic thickness and excellent electronic properties of the graphene monolayer. The dynamical translocation phenomena and/or behaviors underneath the blocked ionic current, however, have not been well unveiled to date for the translocation of DNA electrophoretically through a graphene nanopore. In this report, the assessment on the sensitivity of ionic current to instantaneous statuses of DNA in a 2.4 nm graphene nanopore was carried out based on the all‐atom molecular dynamics simulations. By filtering out the thermal noise of ionic current, the instantaneous conformational variations of DNA in a graphene nanopore have been unveiled from the fluctuations of ionic current, because of the spatial blockage effect of DNA against ionic current. Interestingly, the neighborhood effect of DNA against ionic current was also observed within a distance of 1.5 nm nearby the graphene nanopore, suggesting the further precise control for DNA translocation through a graphene nanopore in gene sequencing. Moreover, the sensitivity of the blocked ionic current toward the instantaneous conformations of DNA in a graphene nanopore demonstrates the great potential of graphene nanopores in the dynamics analysis of single molecules.     

Sampling free energy surfaces as slices by combining umbrella sampling and metadynamics

Metadynamics (MTD) is a very powerful technique to sample high‐dimensional free energy landscapes, and due to its self‐guiding property, the method has been successful in studying complex reactions and conformational changes. MTD sampling is based on filling the free energy basins by biasing potentials and thus for cases with flat, broad, and unbound free energy wells, the computational time to sample them becomes very large. To alleviate this problem, we combine the standard Umbrella Sampling (US) technique with MTD to sample orthogonal collective variables (CVs) in a simultaneous way. Within this scheme, we construct the equilibrium distribution of CVs from biased distributions obtained from independent MTD simulations with umbrella potentials. Reweighting is carried out by a procedure that combines US reweighting and Tiwary–Parrinello MTD reweighting within the Weighted Histogram Analysis Method (WHAM). The approach is ideal for a controlled sampling of a CV in a MTD simulation, making it computationally efficient in sampling flat, broad, and unbound free energy surfaces. This technique also allows for a distributed sampling of a high‐dimensional free energy surface, further increasing the computational efficiency in sampling. We demonstrate the application of this technique in sampling high‐dimensional surface for various chemical reactions using ab initio and QM/MM hybrid molecular dynamics simulations. Further, to carry out MTD bias reweighting for computing forward reaction barriers in ab initio or QM/MM simulations, we propose a computationally affordable approach that does not require recrossing trajectories. © 2016 Wiley Periodicals, Inc.     

DNA‐Based Nanopore Sensing

Nanopore sensing is an attractive, label‐free approach that can measure single molecules. Although initially proposed for rapid and low‐cost DNA sequencing, nanopore sensors have been successfully employed in the detection of a wide variety of substrates. Early successes were mostly achieved based on two main strategies by 1) creating sensing elements inside the nanopore through protein mutation and chemical modification or 2) using molecular adapters to enhance analyte recognition. Over the past five years, DNA molecules started to be used as probes for sensing rather than substrates for sequencing. In this Minireview, we highlight the recent research efforts of nanopore sensing based on DNA‐mediated characteristic current events. As nanopore sensing is becoming increasingly important in biochemical and biophysical studies, DNA‐based sensing may find wider applications in investigating DNA‐involving biological processes.     

Enhanced conformational sampling method for proteins based on the TaBoo SeArch algorithm: Application to the folding of a mini‐protein,chignolin

The conformational samplings are indispensible for obtaining reliable canonical ensembles, which provide statistical averages of physical quantities such as free energies. However, the samplings of vast conformational space of biomacromolecules by conventional molecular dynamics (MD) simulations might be insufficient, due to their inadequate accessible time‐scales for investigating biological functions. Therefore, the development of methodologies for enhancing the conformational sampling of biomacromolecules still remains as a challenging issue in computational biology. To tackle this problem, we newly propose an efficient conformational search method, which is referred as TaBoo SeArch (TBSA) algorithm. In TBSA, an inverse energy histogram is used to select seeds for the conformational resampling so that states with high frequencies are inhibited, while states with low frequencies are efficiently sampled to explore the unvisited conformational space. As a demonstration, TBSA was applied to the folding of a mini‐protein, chignolin, and automatically sampled the native structure (C_α root mean square deviation < 1.0 Å) with nanosecond order computational costs started from a completely extended structure, although a long‐time 1‐µs normal MD simulation failed to sample the native structure. Furthermore, a multiscale free energy landscape method based on the conformational sampling of TBSA were quantitatively evaluated through free energy calculations with both implicit and explicit solvent models, which enable us to find several metastable states on the folding landscape. © 2015 Wiley Periodicals, Inc.     

Comparison of two simulation methods to compute solvation free energies and partition coefficients

The thermodynamic integration (TI) and expanded ensemble (EE) methods are used here to calculate the hydration free energy in water, the solvation free energy in 1‐octanol, and the octanol‐water partition coefficient for a six compounds of varying functionality using the optimized potentials for liquid simulations (OPLS) all‐atom (AA) force field parameters and atomic charges. Both methods use the molecular dynamics algorithm as a primary component of the simulation protocol, and both have found wide applications in fields such as the calculation of activity coefficients, phase behavior, and partition coefficients. Both methods result in solvation free energies and 1‐octanol/water partition coefficients with average absolute deviations (AAD) from experimental data to within 4 kJ/mol and 0.5 log units, respectively. Here, we find that in simulations the OPLS‐AA force field parameters (with fixed charges) can reproduce solvation free energies of solutes in 1‐octanol with AAD of about half that for the solute hydration free energies using a extended simple point charge (SPC/E) model of water. The computational efficiency of the two simulation methods are compared based on the time (in nanoseconds) required to obtain similar standard deviations in the solvation free energies and 1‐octanol/water partition coefficients. By this analysis, the EE method is found to be a factor of nine more efficient than the TI algorithm. For both methods, solvation free energy calculations in 1‐octanol consume roughly an order of magnitude more CPU hours than the hydration free energy calculations. © 2012 Wiley Periodicals, Inc.     

Steered molecular dynamics simulations of Na+ permeation across the gramicidin A channel

The potential of mean forces (PMF) governing Na+ permeation through gramicidin A (gA) channels with explicit water and membrane was characterized using steered molecular dynamics (SMD) simulations. Constant-force SMD with a steering force parallel to the channel axis revealed at least seven energy wells in each monomer of the channel dimer. Except at the channel dimer interface, each energy well is associated with at least three and at most four backbone carbonyl oxygens and two water oxygens in a pseudo-hexahedral or pseudo-octahedral coordination with the Na+ ion. Repeated constant-velocity SMD by dragging a Na+ ion from each energy well in opposite directions parallel to the channel axis allowed the computation of the PMF across the gA channel, revealing a global minimum corresponding to Na+ binding sites near the entrance of gA at +/-9.3 A from the geometric center of the channel. The effect of volatile anesthetics on the PMF was also analyzed in the presence of halothane molecules. Although the accuracy of the current PMF calculation from SMD simulations is not yet sufficient to quantify the PMF difference with and without anesthetics, the comparison of the overall PMF profiles nevertheless confirms that the anesthetics cause insignificant changes to the structural makeup of the free energy wells along the channel and the overall permeation barrier. On average, the PMF appears less rugged in the outer part of the channel in the presence of anesthetics, consistent with our earlier finding that halothane interaction with anchoring residues makes the gA channel more dynamic. A causal relationship was observed between the reorientation of the coordinating backbone carbonyl oxygen and Na+ transit from one energy well to another, suggesting the possibility that even minute changes in the conformation of pore-lining residues due to dynamic motion could be sufficient to trigger the ion permeation. Because some of the carbonyl oxygens contribute to Na+ coordination in two adjacent energy wells, our SMD results reveal that the atomic picture of ion "hopping" through a gA channel actually involves a Na+ ion being carried in a relay by the coordinating oxygens from one energy well to the next. Steered molecular dynamics complements other computational approaches as an attractive means for the atomistic interpretation of experimental permeation studies.     

Toward robust computational electrochemical predicting the environmental fate of organic pollutants

A number of density functionals was utilized for the calculation of electron attachment free energy for nitrocompounds, quinones and azacyclic compounds. Different solvation models have been tested on the calculation of difference in free energies of solvation of oxidized and reduced forms of nitrocompounds in aqueous solution, quinones in acetonitrile, and azacyclic compounds in dimethylformamide. Gas‐phase free energies evaluated at the mPWB1K/tzvp level and solvation energies obtained using SMD model to compute solvation energies of neutral oxidized forms and PCM(Pauling) to compute solvation energies of anion‐radical reduced forms provide reasonable accuracy of the prediction of electron attachment free energy, difference in free solvation energies of oxidized and reduced forms, and as consequence yield reduction potentials in good agreement with experimental data (mean absolute deviation is 0.15 V). It was also found that SMD/M05‐2X/tzvp method provides reduction potentials with deviation of 0.12 V from the experimental values but in cases of nitrocompounds and quinones this accuracy is achieved due to the cancelation of errors. To predict reduction ability of naturally occurred iron containing species with respect to organic pollutants we exploited experimental data within the framework of Pourbaix (E_h ? pH) diagrams. We conclude that surface‐bound Fe(II) as well as certain forms of aqueous Fe(II)_aq are capable of reducing a variety of nitroaromatic compounds, quinones and novel high energy materials under basic conditions (pH > 8). At the same time, zero‐valent iron is expected to be active under neutral and acidic conditions. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Calculation of free energy landscapes: A histogram reweighted metadynamics approach

We present an efficient method for the calculation of free energy landscapes. Our approach involves a history‐dependent bias potential, which is evaluated on a grid. The corresponding free energy landscape is constructed via a histogram reweighting procedure a posteriori. Because of the presence of the bias potential, it can be also used to accelerate rare events. In addition, the calculated free energy landscape is not restricted to the actual choice of collective variables and can in principle be extended to auxiliary variables of interest without further numerical effort. The applicability is shown for several examples. We present numerical results for the alanine dipeptide and the Met‐Enkephalin in explicit solution to illustrate our approach. Furthermore, we derive an empirical formula that allows the prediction of the computational cost for the ordinary metadynamics variant in comparison with our approach, which is validated by a dimensionless representation. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Cover Picture: Electrophoresis 14'2011

Issue no. 14 is a regular issue consisting of 17 contributions distributed over 3 distinct parts. Part I has 7 articles describing studies on proteins and proteomics including measuring protein mobility using laser Doppler electrophoresis, high sensitivity protein analysis by FESI‐CE‐MALDI‐MS, protocol for SDS‐PAGE separation of myosin heavy chain isoforms, proteomics analysis of a spring wheat cultivar in response to prolonged cold stress, analysis of changes in the 20S proteasome, treatment of acute lymphoblastic leukemia with L‐asparaginase, and improved method for immunostaining of mucin. Part II is on nucleic acids and has 4 contributions on multiplex PCR‐CE analysis, copy number variation, ethnic difference, SNPs, CE‐SSCP, allelic ladder and SNaPshot technique. Part III has 6 articles on various fundamentals and methodologies, including nanofluidics, nanopore sensing, sequential injection setup for CIEF combined with MS detection, determination of imatinib mesylate in chronic myeloid leukemia patients by CE, UV detection of plasma malondialdehyde using CE‐FASI, enantioseparation of basic drugs by CE using clarithromycin lactobionate chiral selector and CE determination of bioactive constituents in Herba Houttuyniae. Featured articles include: Advances in the measurement of protein mobility using laser Doppler electrophoresis – the diffusion barrier technique (( 10.1002/elps.201100108 )) Improved method for immunostaining of mucin separated by supported molecular matrix electrophoresis by optimizing the matrix composition and fixation procedure (( 10.1002/elps.201000608 )) Strategy for high‐fidelity multiplex DNA copy number assay system using capillary electrophoresis devices (( 10.1002/elps.201100093 )) Electrokinetic particle translocation through a nanopore containing a floating electrode (( 10.1002/elps.201100050 )) Sequential injection setup for capillary isoelectric focusing combined with MS detection (( 10.1002/elps.201100012 ))     

BAR-based optimum adaptive steered MD for configurational sampling

The equilibrium and nonequilibrium adaptive alchemical free energy simulation methods optimum Bennett's acceptance ratio and optimum crooks' equation (OCE), based on the statistically optimal bidirectional reweighting estimator named Bennett's Acceptance Ratio or Crooks' equation, perform initial sampling in the staging alchemical transformation and then determine the importance rank of different states via the time-derivative of the variance. The method is proven to give speedups compared with the equal time rule. In the current work, we extend the time derivative of variance guided adaptive sampling method to the configurational space, falling in the term of steered MD (SMD). The SMD approach biasing physically meaningful collective variable (CV) such as one dihedral or one distance to pulling the system from one conformational state to another. By minimizing the variance of the free energy differences along the pathway in an optimized way, a new type of adaptive SMD (ASMD) is introduced. As exhibits in the alchemical case, this adaptive sampling method outperforms the traditional equal-time SMD in nonequilibrium stratification. Also, the method gives much more efficient calculation of potential of mean force than the selection criterion-based ASMD scheme, which is proven to be more efficient than traditional SMD. The OCE workflow is periodicity-of-CV dependent while ASMD is not. The performance is demonstrated in a dihedral flipping case and two distance pulling cases, accounting for periodic and nonperiodic CVs, respectively. © 2019 Wiley Periodicals, Inc.