Energy stable scheme for random batch molecular dynamics.

The computational bottleneck of molecular dynamics is pairwise additive long-range interactions between particles. The random batch Ewald (RBE) method provides a highly efficient and superscalable solver for long-range interactions, but the stochastic nature of this algorithm leads to unphysical self-heating effect during the simulation. We propose an energy stable scheme (ESS) for particle systems by employing a Berendsen-type energy bath. The scheme removes the notorious energy drift, which exists due to the force error even when a symplectic integrator is employed. Combining the RBE with the ESS, the new method provides a perfect solution to the computational bottleneck of molecular dynamics at the microcanonical ensemble. Numerical results for a primitive electrolyte and all-atom pure water systems demonstrate the attractive performance of the algorithm, including its dramatically high accuracy, linear complexity, and overcoming the energy drift for long-time simulations.

A microfluidic chip-based multivalent DNA walker amplification biosensor for the simultaneous detection of multiple food-borne pathogens.

The rapid, simultaneous, sensitive detection of the targets has important application prospects for disease diagnosis and biomedical studies. However, in practical applications, the content of the targets is usually very low, and signal amplification strategies are often needed to improve the detection sensitivity. DNAzyme-driven DNA walkers are an excellent signal amplification strategy due to their outstanding specificity and sensitivity. Food-borne pathogens have always been a foremost threat to human health, and it is an urgent demand to develop a simple, rapid, sensitive, and portable detection method for food-borne pathogens. In addition, there are various species of pathogens, and it is difficult to simultaneously detect multiple pathogens by a single DNA walker. For this reason, a substrate strand with three rA cleavage sites was cleverly designed, and a multivalent DNA walker sensor combined with the microfluidic chip technology was proposed for the simultaneous, rapid, sensitive analysis of Vibrio parahaemolyticus, Salmonella typhimurium, and Staphylococcus aureus. The developed sensor could be used to detect pathogens simultaneously and efficiently with low detection limits and wide detection ranges. Moreover, the combination of gold stirring rod enrichment and DNA walker achieved double amplification, which greatly improved the detection sensitivity. More importantly, by changing the design of the substrate chain, the sensor was expected to be used to detect other targets, thus broadening the scope of practical applications. Therefore, the sensor can build novel detection tool platforms in the field of biosensing.

Silencing of Circ_0135889 Restrains Proliferation and Tumorigenicity of Human Neuroblastoma Cells.

INTRODUCTION: Neuroblastoma (NB) is the most common extracranial solid tumor in infants and young children. Circular ribonucleic acid (RNA) hsa_circ_0135889 (circ_0135889; hsa_circ argonaute 2 _001) is highly expressed in multiple cancer tissues, including NB. However, its role in tumor progression of NB was unclear. METHODS: Real-time quantitative polymerace chain reaction was used to detect RNA expression, and western blotting, or immunohistochemistry was used to measure protein expression. Functional experiments were performed using 3-(4,5-dimethylthiazol-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium, 5-Ethynyl-2'- deoxyuridine, Annexin V-fluorescein isothiocyanate/propidium iodide, and transwell assays, as well as xenograft tumor model. The intermolecular interaction was predicted by online databases and confirmed by dual-luciferase reporter assay and RNA pull-down assay. RESULTS: Circ_0135889 and neuronal differentiation 1 (NEUROD1) were upregulated whilst microRNA (miR)-127-5p, was downregulated in NB tumors and immortalized NB cells. Silencing of circ_0135889 could suppress cell proliferation, migration and invasion, but enhance apoptosis rate of NB cells in vitro. More importantly, circ_0135889 depletion inhibited xenograft tumor growth of NB cells. Circ_0135889 was a sponge for miR-127-5p, and inhibition of miR-127-5p counteracted the inhibitory impact of circ_0135889 knockdown on the malignant behaviors of NB cells. Moreover, NEUROD1 was a direct target of miR-127-5p, and miR-127-5p exerted the anti-tumor role in NB cells by targeting NEUROD1. Furthermore, circ_0135889 regulated NEUROD1expression by sponging miR-127-5p. CONCLUSIONS: Circ_0135889 promoted the tumorigenicity of NB by regulating miR-127-5p/NEUROD1 axis, which might provide a promising therapeutic target for NB.

Improved Random Batch Ewald Method in Molecular Dynamics Simulations.

The random batch Ewald (RBE) is an efficient and accurate method for molecular dynamics (MD) simulations of physical systems at the nano/microscale. The method shows great potential to solve the computational bottleneck of long-range interactions, motivating a necessity to accelerate short-range components of the nonbonded interactions for a further speedup of MD simulations. In this work, we present an improved RBE method for the nonbonding interactions by introducing the random batch idea to constructing neighbor lists for the treatment of both the short-range part of the Ewald splitting and the Lennard-Jones potential. The efficiency that the novel neighbor list algorithm owes to the stochastic minibatch strategy can significantly reduce the total number of neighbors. We obtain the error estimate and convergence by theoretical analysis and implement the improved RBE method in the LAMMPS package. Benchmark simulations are performed to demonstrate the accuracy and stability of the algorithm. Numerical tests on computer performance by conducting large-scaled MD simulations for systems including up to 0.1 billion water molecules are run on massive clusters with up to 50000 CPU cores, demonstrating the attractive features such as the high parallel scalability and memory-saving of the method in comparison to the existing methods.

A random batch Ewald method for charged particles in the isothermal-isobaric ensemble.

We develop an accurate, highly efficient, and scalable random batch Ewald (RBE) method to conduct molecular dynamics simulations in the isothermal-isobaric ensemble (the NPT ensemble) for charged particles in a periodic box. After discretizing the Langevin equations of motion derived using suitable Lagrangians, the RBE method builds the mini-batch strategy into the Fourier space in the Ewald summation for the pressure and forces such that the computational cost is reduced to O(N) per time step. We implement the method in the Large-scale Atomic/Molecular Massively Parallel Simulator package and report accurate simulation results for both dynamical quantities and statistics for equilibrium for typical systems including all-atom bulk water and a semi-isotropic membrane system. Numerical simulations on massive supercomputing cluster are also performed to show promising central processing unit efficiency of the RBE.

Superscalability of the random batch Ewald method.

Coulomb interaction, following an inverse-square force-law, quantifies the amount of force between two stationary and electrically charged particles. The long-range nature of Coulomb interactions poses a major challenge to molecular dynamics simulations, which are major tools for problems at the nano-/micro-scale. Various algorithms are developed to calculate the pairwise Coulomb interactions to a linear scale, but poor scalability limits the size of simulated systems. Here, we use an efficient molecular dynamics algorithm with the random batch Ewald method on all-atom systems where the complete Fourier components in the Coulomb interaction are replaced by randomly selected mini-batches. By simulating the N-body systems up to 108 particles using 10 000 central processing unit cores, we show that this algorithm furnishes O(N) complexity, almost perfect scalability, and an order of magnitude faster computational speed when compared to the existing state-of-the-art algorithms. Further examinations of our algorithm on distinct systems, including pure water, a micro-phase separated electrolyte, and a protein solution, demonstrate that the spatiotemporal information on all time and length scales investigated and thermodynamic quantities derived from our algorithm are in perfect agreement with those obtained from the existing algorithms. Therefore, our algorithm provides a promising solution on scalability of computing the Coulomb interaction. It is particularly useful and cost-effective to simulate ultra-large systems, which is either impossible or very costly to conduct using existing algorithms, and thus will be beneficial to a broad range of problems at nano-/micro-scales.

The Theoretical Model, Method, and Applications of Scattering Photon Burst Counting Based on an Objective Scanning Technique.

Scattering photon burst counting (SPBC) is a single-particle detection method, which is based on measuring scattering photon bursting of single nanoparticles through a detection volume of <1 fL. Although SPBC has been used for bioassays and analysis of nanoparticles, it is necessary to establish its theoretical model and develop a new detection mode in order to further enhance its sensitivity and enlarge its application fields. In this paper, we proposed a theoretical model for the confocal SPBC method and developed a novel SPBC detection mode using the fast objective scanning technique. The computer simulations and experiments documented that this model well describes the relation between photon counts and experimental parameters (such as nanoparticle concentration and diameter, temperature, and viscosity). Based on this model, we developed a novel SPBC detection mode by using the fast objective scanning technique. Compared to the current confocal SPBC method, the sensitivity of this new method was significantly increased due to the significantly increased photon counts per sampling time, the linear detection range is from 0.9 to 90 pM, and the limit of detection is reduced to 40 fM for 30 nm gold nanoparticles. Furthermore, this new method was successfully applied to determine the enzyme activity of caspase-3 and evaluate the inhibition effectiveness of some inhibitors.

Random-batch list algorithm for short-range molecular dynamics simulations.

We propose a fast method for the calculation of short-range interactions in molecular dynamics simulations. The so-called random-batch list method is a stochastic version of the classical neighbor-list method to avoid the construction of a full Verlet list, which introduces two-level neighbor lists for each particle such that the neighboring particles are located in core and shell regions, respectively. Direct interactions are performed in the core region. For the shell zone, we employ a random batch of interacting particles to reduce the number of interaction pairs. The error estimate of the algorithm is provided. We investigate the Lennard-Jones fluid by molecular dynamics simulations and show that this novel method can significantly accelerate the simulations with a factor of several fold without loss of the accuracy. This method is simple to implement, can be well combined with any linked cell methods to further speed up and scale up the simulation systems, and can be straightforwardly extended to other interactions, such as Ewald short-range part, and thus it is promising for large-scale molecular dynamics simulations.

Effects of Kinetic Dielectric Decrement on Ion Diffusion and Capacitance in Electrochemical Systems.

Diffusion of ionic components in electrolytes not only eliminates the gradients of ionic concentrations but also alters the local dielectric environment, and the coupling effect between kinetic dielectric decrement and ionic concentration gradient on the diffusion dynamics is not well understood. Herein, taking the charging process in electrical double layer systems as a case study, we conduct a multiscale investigation of ion diffusions in aqueous electrolytes by combining the dynamic density functional theory and an ion-concentration-dependent dielectric constant model. By properly considering the time evolutions of local dielectric constant coupled with ion density, we report an interesting phenomenon on the suppression of surface charge density that is not captured by conventional models. In addition, we show that the usage of aqueous electrolyte with small dielectric decrement coefficients promotes the capacitance, in quantitative agreement with experimental measurements.

Harmonic surface mapping algorithm for molecular dynamics simulations of particle systems with planar dielectric interfaces.

We have developed an accurate and efficient method for molecular dynamics simulations of charged particles confined by planar dielectric interfaces. The algorithm combines the image-charge method for near field with the harmonic surface mapping, which converts the contribution of infinite far-field charges into a finite number of charges on an auxiliary spherical surface. We approximate the electrostatic potential of far-field charges via spherical harmonic expansion and determine the coefficients by fitting the Dirichlet-to-Neumann boundary condition, which only requires the potential within the simulation cell. Instead of performing the direct evaluation of spherical harmonic series expansion, we use Green's second identity to transform the series expansion into a spherical integral, which can be accurately represented by discrete charges on the sphere. Therefore, the fast multipole method can be readily employed to sum over all charges within and on the sphere, achieving truly linear O(N) complexity. Our algorithm can be applied to a broad range of charged complex fluids under dielectric confinement.

System-integrated technology-enabled model of care to improve the health of stroke patients in rural China: protocol for SINEMA-a cluster-randomized controlled trial.

BACKGROUND: Despite the significant burden of stroke in rural China, secondary prevention of stroke is suboptimal. This study aims to develop a SINEMA for the secondary prevention of stroke in rural China and to evaluate the effectiveness of the model compared with usual care. METHODS: The SINEMA model is being implemented and evaluated through a 1-year cluster-randomized controlled trial in Nanhe County, Hebei Province in China. Fifty villages from 5 townships are randomized in a 1:1 ratio to either the intervention or the control arm (usual care) with a target to enroll 25 stroke survivors per village. Village doctors in the intervention arm (1) receive systematic cascade training by stroke specialists on clinical guidelines, essential medicines and behavior change; (2) conduct monthly follow-up visits with the support of a mobile phone application designed for this study; (3) participate in virtual group activities with other village doctors; 4) receive performance feedback and payment. Stroke survivors participate in a health education and project briefing session, receive monthly follow-up visits by village doctors and receive a voice message call daily as reminders for medication use and physical activities. Baseline and 1-year follow-up survey will be conducted in all villages by trained staff who are blinded of the randomized allocation of villages. The primary outcome will be systolic blood pressure and the secondary outcomes will include diastolic blood pressure, medication adherence, mobility, physical activity level and quality of life. Process and economic evaluation will also be conducted. DISCUSSION: This study is one of very few that aim to promote secondary prevention of stroke in resource-constrained settings and the first to incorporate mobile technologies for both healthcare providers and patients in China. The SINEMA model is innovative as it builds the capacity of primary healthcare workers in the rural area, uses mobile health technologies at the point of care, and addresses critical health needs for a vulnerable community-dwelling patient group. The findings of the study will provide translational evidence for other resource-constrained settings in developing strategies for the secondary prevention of stroke.

Efficient dynamic simulations of charged dielectric colloids through a novel hybrid method.

Modern particle-based simulations increasingly incorporate polarization charges arising from spatially nonuniform permittivity. For complex dielectric geometries, calculation of these induced many-body effects typically requires numerical solvers based upon boundary-element methods, which very significantly increase the required computational effort. For the special case of dielectric spheres, such as colloids or nanoparticles, we recently proposed a semianalytical spectrally accurate hybrid method that combines the method of moments, the image-charge method, and the fast multipole method. The hybrid method is efficient and accurate even when dielectric spheres are closely packed. Here, we extend the method to the evaluation of direct and induced electrostatic forces and demonstrate how this can be incorporated in molecular dynamics simulations. The choice of the relevant numerical parameters for molecular dynamics simulations is discussed in detail, as well as comparisons to the boundary-element method. As a concrete example, we examine the challenging case of binary crystal structures composed of close-packed dielectric colloidal spheres.

Accurate image-charge method by the use of the residue theorem for core-shell dielectric sphere.

An accurate image-charge method (ICM) is developed for ionic interactions outside a core-shell structured dielectric sphere. Core-shell particles have wide applications for which the theoretical investigation requires efficient methods for the Green's function used to calculate pairwise interactions of ions. The ICM is based on an inverse Mellin transform from the coefficients of spherical harmonic series of the Green's function such that the polarization charge due to dielectric boundaries is represented by a series of image point charges and an image line charge. The residue theorem is used to accurately calculate the density of the line charge. Numerical results show that the ICM is promising in fast evaluation of the Green's function, and thus it is useful for theoretical investigations of core-shell particles. This routine can also be applicable for solving other problems with spherical dielectric interfaces such as multilayered media and Debye-Hückel equations.

Harmonic surface mapping algorithm for fast electrostatic sums.

We propose a harmonic surface mapping algorithm (HSMA) for electrostatic pairwise sums of an infinite number of image charges. The images are induced by point sources within a box due to a specific boundary condition which can be non-periodic. The HSMA first introduces an auxiliary surface such that the contribution of images outside the surface can be approximated by the least-squares method using spherical harmonics as basis functions. The so-called harmonic surface mapping is the procedure to transform the approximate solution into a surface charge and a surface dipole over the auxiliary surface, which becomes point images by using numerical integration. The mapping procedure is independent of the number of the sources and is considered to have a low complexity. The electrostatic interactions are then among those charges within the surface and at the integration points, which are all the forms of Coulomb potential and can be accelerated straightforwardly by the fast multipole method to achieve linear scaling. Numerical calculations of the Madelung constant of a crystalline lattice, electrostatic energy of ions in a metallic cavity, and the time performance for large-scale systems show that the HSMA is accurate and fast, and thus is attractive for many applications.

Ion Structure Near a Core-Shell Dielectric Nanoparticle.

A generalized image charge formulation is proposed for the Green's function of a core-shell dielectric nanoparticle for which theoretical and simulation investigations are rarely reported due to the difficulty of resolving the dielectric heterogeneity. Based on the formulation, an efficient and accurate algorithm is developed for calculating electrostatic polarization charges of mobile ions, allowing us to study related physical systems using the Monte Carlo algorithm. The computer simulations show that a fine-tuning of the shell thickness or the ion-interface correlation strength can greatly alter electric double-layer structures and capacitances, owing to the complicated interplay between dielectric boundary effects and ion-interface correlations.

Self-consistent field model for strong electrostatic correlations and inhomogeneous dielectric media.

Electrostatic correlations and variable permittivity of electrolytes are essential for exploring many chemical and physical properties of interfaces in aqueous solutions. We propose a continuum electrostatic model for the treatment of these effects in the framework of the self-consistent field theory. The model incorporates a space- or field-dependent dielectric permittivity and an excluded ion-size effect for the correlation energy. This results in a self-energy modified Poisson-Nernst-Planck or Poisson-Boltzmann equation together with state equations for the self energy and the dielectric function. We show that the ionic size is of significant importance in predicting a finite self energy for an ion in an inhomogeneous medium. Asymptotic approximation is proposed for the solution of a generalized Debye-Hückel equation, which has been shown to capture the ionic correlation and dielectric self energy. Through simulating ionic distribution surrounding a macroion, the modified self-consistent field model is shown to agree with particle-based Monte Carlo simulations. Numerical results for symmetric and asymmetric electrolytes demonstrate that the model is able to predict the charge inversion at high correlation regime in the presence of multivalent interfacial ions which is beyond the mean-field theory and also show strong effect to double layer structure due to the space- or field-dependent dielectric permittivity.

Simulation of electric double layers around charged colloids in aqueous solution of variable permittivity.

The ion distribution around charged colloids in solution has been investigated intensely during the last decade. However, few theoretical approaches have included the influence of variation in the dielectric permittivity within the system, let alone in the surrounding solvent. In this article, we introduce two relatively new methods that can solve the Poisson equation for systems with varying permittivity. The harmonic interpolation method approximately solves the Green's function in terms of a spherical harmonics series, and thus provides analytical ion-ion potentials for the Hamiltonian of charged systems. The Maxwell equations molecular dynamics algorithm features a local approach to electrostatics, allowing for arbitrary local changes of the dielectric constant. We show that the results of both methods are in very good agreement. We also found that the renormalized charge of the colloid, and with it the effective far field interaction, significantly changes if the dielectric properties within the vicinity of the colloid are changed.

Ionic Size Effects: Generalized Boltzmann Distributions, Counterion Stratification, and Modified Debye Length.

Near a charged surface, counterions of different valences and sizes cluster; and their concentration profiles stratify. At a distance from such a surface larger than the Debye length, the electric field is screened by counterions. Recent studies by a variational mean-field approach that includes ionic size effects and by Monte Carlo simulations both suggest that the counterion stratification is determined by the ionic valence-to-volume ratios. Central in the mean-field approach is a free-energy functional of ionic concentrations in which the ionic size effects are included through the entropic effect of solvent molecules. The corresponding equilibrium conditions define the generalized Boltzmann distributions relating the ionic concentrations to the electrostatic potential. This paper presents a detailed analysis and numerical calculations of such a free-energy functional to understand the dependence of the ionic charge density on the electrostatic potential through the generalized Boltzmann distributions, the role of ionic valence-to-volume ratios in the counterion stratification, and the modification of Debye length due to the effect of ionic sizes.

A Screening Condition Imposed Stochastic Approximation for Long-Range Electrostatic Correlations.

The recent random batch Ewald algorithm, originating from a stochastic approximation, performs 1 order of magnitude faster than the mainstream algorithms such as the particle-particle particle-mesh method for handling of long-range electrostatics in large-scale simulations. However, this algorithm fails to fully capture the long-range electrostatic correlations. Here, we demonstrate that, when incorporating a known screening condition in the stochastic approximation, the algorithm could be simply amended without the loss of any efficiency.

Dual-Mode Biosensor for Simultaneous and Rapid Detection of Live and Whole Salmonella typhimurium Based on Bioluminescence and Fluorescence Detection.

Both live and dead Salmonella typhimurium (S.T) are harmful to human health, but there are differences in pathological mechanism, dosage, and security. It is crucial to develop a rapid and simultaneous assay to distinguish and quantify live and dead S.T in foods. Herein, one dual-mode biosensor for simultaneous detection of live and dead S.T was fabricated based on two phage probes, using portable bioluminescence and fluorescent meter as detectors, respectively. Firstly, a magnetic phage capture probe (M-P1) and a phage signal tag (P2-S) labeled with SYTO 13 fluorescent dye were prepared, respectively. Both M-P1 and P2-S can specifically conjugate with S.T to form a magnetic sandwich complex. After magnetic separation, the isolated complex can emit a fluorescent signal under an excited 365 nm laser, which can reflect the total amount of S.T. Afterwards, the lysozyme was added to decompose the captured live S.T, which can release ATP and produce a bioluminescent signal corresponding to the live S.T amount. The dead S.T concentration can be deduced by the difference between total and live examples. The detection limit of 55 CFU/mL for total S.T and 9 CFU/mL for live ones was within 20 min. The assay was successfully employed in milk samples and prospectively for on-site screening of other dead and live bacteria, while changing the phages for the targets.