Using Brownian dynamics simulation, we studied the effect of histone modifications On conformations of an array of nucleosomes in a segment of chromatin. The simulation demonstrated that the segment of chromatin shows...Using Brownian dynamics simulation, we studied the effect of histone modifications On conformations of an array of nucleosomes in a segment of chromatin. The simulation demonstrated that the segment of chromatin shows the dynamic behaviour that its conformation can switch between a state with nearly all of the histones being wrapped by DNA and a state with nearly all of the histones being unwrapped by DNA, thus involving the "cross-talking" interactions among the histones. Each state can stay for a sufficiently long time. These conformational states are essential for gene expression or gene silence. The simulation also shows that these conformational states can be inherited by the daughter DNAs during DNA replication, giving a theoretical explanation of the epigenetic phenomenon.展开更多
The dynamics of two confined colloidal particles is studied by means of Brownian dynamics simulation. The autocorrelation function and cross-correlation function of the two colloidal spheres are computed by utilizing ...The dynamics of two confined colloidal particles is studied by means of Brownian dynamics simulation. The autocorrelation function and cross-correlation function of the two colloidal spheres are computed by utilizing the formulae of hydrodynamic diffusion matrix expanded to different orders, as well as the accurate tensor through numerical algorithm. Furthermore, the numerical results are compared with the experimental results and the theoretical approximation. It is found that the relatively simple theoretical approximation gives good predictions when two spheres are far away from each other, but fails when the two spheres are very close.展开更多
We investigate the statistics of polymer capture by a nanopore using Brownian dynamics simulations. It is found that when the velocity flux is greater than a critical velocity flux, the capture picture is a random sel...We investigate the statistics of polymer capture by a nanopore using Brownian dynamics simulations. It is found that when the velocity flux is greater than a critical velocity flux, the capture picture is a random selection process, otherwise it tends to a statistical process governed by energetic considerations. In addition, the chain ends capture probability decreases as the chain length increases and satisfies a power-law scaling of P0(N)-N^-0.8.展开更多
Brownian dynamics simulations are employed to explore the effects of chain stiffness and trivalent salt concentration on the conformational behavior of spherical polyelectrolyte brush. The rigid brush adopts bundle-li...Brownian dynamics simulations are employed to explore the effects of chain stiffness and trivalent salt concentration on the conformational behavior of spherical polyelectrolyte brush. The rigid brush adopts bundle-like morphology at a wide range of trivalent salt concentration. The number variation of bundles pinned on the colloid surface shows a non-monotonic profile as a function of the chain stiffness. The radial distributions of monomers and ions and the charge ratio between condensed ions and monomers are calculated. The charge inversion is observed for the high salt concentration regardless of chain rigidity. Furthermore, the pair correlation functions of monomer-monomer and monomer-salt cation are used to elucidate the aggregated mechanism of the bundle-like structure.展开更多
It is well known that potassium ion channels have higher permeability than K ions, and the permeable rate of a single K ion channel is about 108 ions per second. We develop a hierarchical model of potassium ion channe...It is well known that potassium ion channels have higher permeability than K ions, and the permeable rate of a single K ion channel is about 108 ions per second. We develop a hierarchical model of potassium ion channel permeation involving ab initio quantum calculations and Brownian dynamics simulations, which can consistently explain a range of channel dynamics. The results show that the average velocity of K ions, the mean permeable time of K ions and the permeable rate of single channel are about 0.92nm/ns, 4.35ns and 2.30 ×10^8 ions/s, respectively.展开更多
By the so-called wormlike chain (WLC) model in polymer physics envision- ing an isotropic rod that is continuously flexible, the force-extension relations of semi- flexible polymer chains strongly constrained by var...By the so-called wormlike chain (WLC) model in polymer physics envision- ing an isotropic rod that is continuously flexible, the force-extension relations of semi- flexible polymer chains strongly constrained by various confinements are theoretically investigated, including a slab-like confinement where the polymer chains are sandwiched between two parallel impenetrable walls, and a capped nanochannel confinement with a circular or rectangular cross-section where the chains are bounded in three directions. The Brownian dynamics (BD) simulations based on the generalized bead-rod (GBR) model are performed to verify the theoretical predictions.展开更多
We simulate the self-assembly of active colloidal molecules from binary mixtures of spherical particles using a Brownian dynamics algorithm.These particles interact via phoretic interactions,which are determined by tw...We simulate the self-assembly of active colloidal molecules from binary mixtures of spherical particles using a Brownian dynamics algorithm.These particles interact via phoretic interactions,which are determined by two independently tunable parameters,surface activity and surface mobility.In systems composed of equal-size particles,we observe the formation of colloidal molecules with well-defined coordination numbers and spatial arrangement,which also display distinct dynamic functions,such as resting,translating,and rotating.By changing the size ratio to 2:1 between the two species,we further observe the formation of colloidal molecules with new structures arising from breaking the size symmetry.By tuning the mutual interactions between the smaller species via their surface mobility,we are able to control their spacing as well as the coordination number of the colloidal molecules.This study highlights the importance of tuning surface parameters and size asymmetry in controlling the structure and the active dynamics of colloidal molecules.展开更多
We study the conformation of an amphiphilic comb-like copolymer in a selective solvent by using Brownian dynamics simulations.Our results demonstrate that,there is an optimal total degree of polymerization to form sph...We study the conformation of an amphiphilic comb-like copolymer in a selective solvent by using Brownian dynamics simulations.Our results demonstrate that,there is an optimal total degree of polymerization to form spherical micelles.Only when the total degree of polymerization is less than the optimal total degree of polymerization,a unimolecular spherical micelle can be formed,otherwise,the spherical micelle will be transformed into a cylindrical micelle.When the total degree of polymerization is less than the optimal total degree of polymerization,the radius of gyration and the stretching factor meet an exponential relationship,where the power exponent decreases with the increase of the length of hydrophilic side chains.The fundamental principles uncovered here afford insights for the optimization of molecular design of amphiphilic comb-like copolymers in various applications.展开更多
Combining Brownian dynamics simulations and self-consistent field theory,we demonstrate that the architectural parameters of coreshell comb-like chains have a decisive influence on their final self-assembled conformat...Combining Brownian dynamics simulations and self-consistent field theory,we demonstrate that the architectural parameters of coreshell comb-like chains have a decisive influence on their final self-assembled conformations.When the ratio of hydrophilic and hydrophobic segments is N_(g,g)∶N_(g,r)=3∶1,no core-shell structure is observed,and unsegregated chains or clusters are substituted.When N_(g,g)∶N_(g,r)=3∶2or 3∶3,the core-shell comb-like chains can be assembled into single micelles or structures where several small micelles are strung together by the backbone,which is similar to the pearl-necklace structure formed by the polyelectrolytes in poor solvent.With the increase of backbone length or grafting density,the probability of forming single micelles becomes lower,but the structures of two,three or more small micelles strung together are more observed.Our results indicate how to obtain the desired self-assembled structures of core-shell comb-like chains by regulating the architectural parameters,which could provide insights for the optimization of molecular design in various applications.展开更多
We investigate the sedimentation dynamics of a binary mixture,the species of which differ by their Stokes coefficients but are identical otherwise.We analyze the sedimentation dynamics and the morphology of the final ...We investigate the sedimentation dynamics of a binary mixture,the species of which differ by their Stokes coefficients but are identical otherwise.We analyze the sedimentation dynamics and the morphology of the final deposits using Brownian dynamics simulations for mixtures with a range of sedimentation velocities of both species.In addition,we use the lattice Boltzmann method to study hydrodynamic effects.We found a threshold in the difference of the sedimentation velocities above which the species in the final deposit are segregated.The degree of segregation increases with the difference in the Stokes coefficients or the sedimentation velocities above the threshold.We propose a simple analytical model that captures the main features of the simulated deposits.展开更多
We investigated the effect of hydrodynamic interaction(HI) on flow-induced polymer translocation through a nanotube by Brownian dynamics simulations. Whether there is HI in the simulation system is separately contro...We investigated the effect of hydrodynamic interaction(HI) on flow-induced polymer translocation through a nanotube by Brownian dynamics simulations. Whether there is HI in the simulation system is separately controlled by using different diffusion tensors. It is found that HI has no effect on critical velocity flux for long po- lymer chains due to the competition between more drag force and the hindrance of chain stretching from HI, however, HI broadens the transition interval. In addition, for flow-induced polymer translocation with HI, the critical velocity flux firstly slowly decreases with the increase of chain length and then becomes identical to that of it without HI, that is, the critical velocity flux is independent of chain length. At the same time, HI also accelerates the translocation process and makes the relative variation amplitude of single bead translocation time smaller. In fact, HI can enhance the intrachain cooperativity to make the whole chain obtain more drag force from fluid field and hinder chain stret- ching, both of which play an important role in translocation process.展开更多
As an example of a very low friction system,Monte Carlo Brownian dynamics simulations have been used to calculate equilibrium structures of a polyelectrolyte brush grafted onto planes.The polymers were calculated in a...As an example of a very low friction system,Monte Carlo Brownian dynamics simulations have been used to calculate equilibrium structures of a polyelectrolyte brush grafted onto planes.The polymers were calculated in a semi-flexible coarse-grain model that is appropriate to treat the charge density of the polyion.The effect of linear charge density on the polyionx,the surface negative charge,and added salts were studied.In salt-free solution,scaling theories predicted the structure well in the lowxregion.In the highxregion,additional shrinkage was found from the theories due to counterion condensation.The effect of surface charge showed not only the repulsion of the polyion from the surface but also the shrinkage in the highxregion due to the additional counterions required for electrical neutrality.The addition of salts led to the shrinkage of the brush heights,and in the highxregion,additional extension was found.The computational strategy for calculating the friction dynamics of the system is also discussed.展开更多
基金Project supported by the National Natural Science Foundation of China (Grants Nos. 10834014 and 10674173)the National Basic Research Program of China (Grant No. 2009CB930704)
文摘Using Brownian dynamics simulation, we studied the effect of histone modifications On conformations of an array of nucleosomes in a segment of chromatin. The simulation demonstrated that the segment of chromatin shows the dynamic behaviour that its conformation can switch between a state with nearly all of the histones being wrapped by DNA and a state with nearly all of the histones being unwrapped by DNA, thus involving the "cross-talking" interactions among the histones. Each state can stay for a sufficiently long time. These conformational states are essential for gene expression or gene silence. The simulation also shows that these conformational states can be inherited by the daughter DNAs during DNA replication, giving a theoretical explanation of the epigenetic phenomenon.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10174041 and 10334020).
文摘The dynamics of two confined colloidal particles is studied by means of Brownian dynamics simulation. The autocorrelation function and cross-correlation function of the two colloidal spheres are computed by utilizing the formulae of hydrodynamic diffusion matrix expanded to different orders, as well as the accurate tensor through numerical algorithm. Furthermore, the numerical results are compared with the experimental results and the theoretical approximation. It is found that the relatively simple theoretical approximation gives good predictions when two spheres are far away from each other, but fails when the two spheres are very close.
基金financially supported by the National Basic Research Program of China(No.2009CB930100)the National Natural Science Foundation of China(Nos.21234007,21304097 and 51473168)the Joint Research Fund for Overseas Chinese,Hong Kong and Macao Young Scientists of the National Natural Science Foundation of China(No.51028301)
文摘We investigate the statistics of polymer capture by a nanopore using Brownian dynamics simulations. It is found that when the velocity flux is greater than a critical velocity flux, the capture picture is a random selection process, otherwise it tends to a statistical process governed by energetic considerations. In addition, the chain ends capture probability decreases as the chain length increases and satisfies a power-law scaling of P0(N)-N^-0.8.
基金financially supported by the National Natural Science Foundation of China (No. 21474005)the Fundamental Research Funds for the Central Universities (No. 3122016L011)
文摘Brownian dynamics simulations are employed to explore the effects of chain stiffness and trivalent salt concentration on the conformational behavior of spherical polyelectrolyte brush. The rigid brush adopts bundle-like morphology at a wide range of trivalent salt concentration. The number variation of bundles pinned on the colloid surface shows a non-monotonic profile as a function of the chain stiffness. The radial distributions of monomers and ions and the charge ratio between condensed ions and monomers are calculated. The charge inversion is observed for the high salt concentration regardless of chain rigidity. Furthermore, the pair correlation functions of monomer-monomer and monomer-salt cation are used to elucidate the aggregated mechanism of the bundle-like structure.
基金Supported by the National Natural Science Foundation of China under Grant Nos 10474018 and 10375016, the Hebei Provincial Natural Science Foundation under Grant Nos C2005000011 and A2004000005, the Subsidization for Construction Programme of Key Subjects in Universities and Colleges of Hebei Province, and the Hebei Provincial Doctoral Foundation under Grant No 2006148.
文摘It is well known that potassium ion channels have higher permeability than K ions, and the permeable rate of a single K ion channel is about 108 ions per second. We develop a hierarchical model of potassium ion channel permeation involving ab initio quantum calculations and Brownian dynamics simulations, which can consistently explain a range of channel dynamics. The results show that the average velocity of K ions, the mean permeable time of K ions and the permeable rate of single channel are about 0.92nm/ns, 4.35ns and 2.30 ×10^8 ions/s, respectively.
基金supported the National Natural Science Foundation of China(Nos.11032006,11072094,and11121202)the Ph.D.Program Foundation of Ministry of Education of China(No.20100211110022)+1 种基金the National Key Project of Magneto-Constrained Fusion Energy Development Program of China(No.2013GB110002)the Fundamental Research Funds for the Central Universities(No.lzujbky2013-1)
文摘By the so-called wormlike chain (WLC) model in polymer physics envision- ing an isotropic rod that is continuously flexible, the force-extension relations of semi- flexible polymer chains strongly constrained by various confinements are theoretically investigated, including a slab-like confinement where the polymer chains are sandwiched between two parallel impenetrable walls, and a capped nanochannel confinement with a circular or rectangular cross-section where the chains are bounded in three directions. The Brownian dynamics (BD) simulations based on the generalized bead-rod (GBR) model are performed to verify the theoretical predictions.
基金the Innovation Program of Guangdong Provincial Department of Education,China(Grant No.2019KTSCX148)the Science and Technology Innovation Commission of Shenzhen(Grant No.JCYJ20170818141727254).
文摘We simulate the self-assembly of active colloidal molecules from binary mixtures of spherical particles using a Brownian dynamics algorithm.These particles interact via phoretic interactions,which are determined by two independently tunable parameters,surface activity and surface mobility.In systems composed of equal-size particles,we observe the formation of colloidal molecules with well-defined coordination numbers and spatial arrangement,which also display distinct dynamic functions,such as resting,translating,and rotating.By changing the size ratio to 2:1 between the two species,we further observe the formation of colloidal molecules with new structures arising from breaking the size symmetry.By tuning the mutual interactions between the smaller species via their surface mobility,we are able to control their spacing as well as the coordination number of the colloidal molecules.This study highlights the importance of tuning surface parameters and size asymmetry in controlling the structure and the active dynamics of colloidal molecules.
基金financially supported by the Scientific Research Project of Yili Normal University(No.2022YSYB009)。
文摘We study the conformation of an amphiphilic comb-like copolymer in a selective solvent by using Brownian dynamics simulations.Our results demonstrate that,there is an optimal total degree of polymerization to form spherical micelles.Only when the total degree of polymerization is less than the optimal total degree of polymerization,a unimolecular spherical micelle can be formed,otherwise,the spherical micelle will be transformed into a cylindrical micelle.When the total degree of polymerization is less than the optimal total degree of polymerization,the radius of gyration and the stretching factor meet an exponential relationship,where the power exponent decreases with the increase of the length of hydrophilic side chains.The fundamental principles uncovered here afford insights for the optimization of molecular design of amphiphilic comb-like copolymers in various applications.
基金financially supported by the Scientific Research Project of Yili Normal University(No.2022YSYB009)the National Natural Science Foundation of China(No.2019M651340)。
文摘Combining Brownian dynamics simulations and self-consistent field theory,we demonstrate that the architectural parameters of coreshell comb-like chains have a decisive influence on their final self-assembled conformations.When the ratio of hydrophilic and hydrophobic segments is N_(g,g)∶N_(g,r)=3∶1,no core-shell structure is observed,and unsegregated chains or clusters are substituted.When N_(g,g)∶N_(g,r)=3∶2or 3∶3,the core-shell comb-like chains can be assembled into single micelles or structures where several small micelles are strung together by the backbone,which is similar to the pearl-necklace structure formed by the polyelectrolytes in poor solvent.With the increase of backbone length or grafting density,the probability of forming single micelles becomes lower,but the structures of two,three or more small micelles strung together are more observed.Our results indicate how to obtain the desired self-assembled structures of core-shell comb-like chains by regulating the architectural parameters,which could provide insights for the optimization of molecular design in various applications.
基金support from the Portuguese Foundation for Science and Technology(FCT)under the contracts no.EXCL/FIS-NAN/0083/2012,SFRH/BD/119240/2016,UIDB/00618/2020,UIDP/00618/2020,EXPL/FISMAC/0406/2021,PTDC/FIS-MAC/28146/2017(LISBOA-01-0145-FEDER-028146)and PTDC/FIS-MAC/5689/2020.
文摘We investigate the sedimentation dynamics of a binary mixture,the species of which differ by their Stokes coefficients but are identical otherwise.We analyze the sedimentation dynamics and the morphology of the final deposits using Brownian dynamics simulations for mixtures with a range of sedimentation velocities of both species.In addition,we use the lattice Boltzmann method to study hydrodynamic effects.We found a threshold in the difference of the sedimentation velocities above which the species in the final deposit are segregated.The degree of segregation increases with the difference in the Stokes coefficients or the sedimentation velocities above the threshold.We propose a simple analytical model that captures the main features of the simulated deposits.
基金Supported by the National Basic Research Program of China(No.2009CB930100), the National Natural Science Foundation of China(Nos.21234007, 21304097, 51473168) and the Joint Research Fund for Overseas Chinese, Hong Kong and Macao Young Scientists of the National Natural Science Foundation of China(No.51028301).
文摘We investigated the effect of hydrodynamic interaction(HI) on flow-induced polymer translocation through a nanotube by Brownian dynamics simulations. Whether there is HI in the simulation system is separately controlled by using different diffusion tensors. It is found that HI has no effect on critical velocity flux for long po- lymer chains due to the competition between more drag force and the hindrance of chain stretching from HI, however, HI broadens the transition interval. In addition, for flow-induced polymer translocation with HI, the critical velocity flux firstly slowly decreases with the increase of chain length and then becomes identical to that of it without HI, that is, the critical velocity flux is independent of chain length. At the same time, HI also accelerates the translocation process and makes the relative variation amplitude of single bead translocation time smaller. In fact, HI can enhance the intrachain cooperativity to make the whole chain obtain more drag force from fluid field and hinder chain stret- ching, both of which play an important role in translocation process.
基金supported by the MEXT program“Elements Strategy Initiative to Form Core Research Center”(since 2012),MEXTMinistry of Education Culture,Sports,Science and Technology,Japan.
文摘As an example of a very low friction system,Monte Carlo Brownian dynamics simulations have been used to calculate equilibrium structures of a polyelectrolyte brush grafted onto planes.The polymers were calculated in a semi-flexible coarse-grain model that is appropriate to treat the charge density of the polyion.The effect of linear charge density on the polyionx,the surface negative charge,and added salts were studied.In salt-free solution,scaling theories predicted the structure well in the lowxregion.In the highxregion,additional shrinkage was found from the theories due to counterion condensation.The effect of surface charge showed not only the repulsion of the polyion from the surface but also the shrinkage in the highxregion due to the additional counterions required for electrical neutrality.The addition of salts led to the shrinkage of the brush heights,and in the highxregion,additional extension was found.The computational strategy for calculating the friction dynamics of the system is also discussed.