共查询到20条相似文献，搜索用时 15 毫秒
1.
We present new generalized-ensemble molecular dynamics simulation algorithms, which we refer to as the multibaric-multithermal molecular dynamics. We describe three algorithms based on (1) the Nosé thermostat and the Andersen barostat, (2) the Nosé-Poincaré thermostat and the Andersen barostat, and (3) the Gaussian thermostat and the Andersen barostat. The multibaric-multithermal simulations perform random walks widely both in the potential-energy space and in the volume space. Therefore, one can calculate isobaric-isothermal ensemble averages in wide ranges of temperature and pressure from only one simulation run. We test the effectiveness of the multibaric-multithermal algorithm by applying it to a Lennard-Jones 12-6 potential system. 相似文献
2.
Self‐assembly of molecular systems is an important and general problem that intertwines physics, chemistry, biology, and material sciences. Through understanding of the physical principles of self‐organization, it often becomes feasible to control the process and to obtain complex structures with tailored properties, for example, bacteria colonies of cells or nanodevices with desired properties. Theoretical studies and simulations provide an important tool for unraveling the principles of self‐organization and, therefore, have recently gained an increasing interest. The present article features an extension of a popular code MBN Explorer (MesoBioNano Explorer) aiming to provide a universal approach to study self‐assembly phenomena in biology and nanoscience. In particular, this extension involves a highly parallelized module of MBN Explorer that allows simulating stochastic processes using the kinetic Monte Carlo approach in a three‐dimensional space. We describe the computational side of the developed code, discuss its efficiency, and apply it for studying an exemplary system. © 2014 Wiley Periodicals, Inc. 相似文献
3.
We propose a new type of the Hamiltonian replica‐exchange method (REM) for molecular dynamics (MD) and Monte Carlo simulations, which we refer to as the Coulomb REM (CREM). In this method, electrostatic charge parameters in the Coulomb interactions are exchanged among replicas while temperatures are exchanged in the usual REM. By varying the atom charges, the CREM overcomes free‐energy barriers and realizes more efficient sampling in the conformational space than the REM. Furthermore, this method requires only a smaller number of replicas because only the atom charges of solute molecules are used as exchanged parameters. We performed Coulomb replica‐exchange MD simulations of an alanine dipeptide in explicit water solvent and compared the results with those of the conventional canonical, replica exchange, and van der Waals REMs. Two force fields of AMBER parm99 and AMBER parm99SB were used. As a result, the CREM sampled all local‐minimum free‐energy states more frequently than the other methods for both force fields. Moreover, the Coulomb, van der Waals, and usual REMs were applied to a fragment of an amyloid‐β peptide (Aβ) in explicit water solvent to compare the sampling efficiency of these methods for a larger system. The CREM sampled structures of the Aβ fragment more efficiently than the other methods. We obtained β‐helix, α‐helix, 310‐helix, β‐hairpin, and β‐sheet structures as stable structures and deduced pathways of conformational transitions among these structures from a free‐energy landscape. © 2012 Wiley Periodicals, Inc. 相似文献
4.
An extended system Hamiltonian is proposed to perform molecular dynamics (MD) simulation in the grand canonical ensemble. The Hamiltonian is similar to the one proposed by Lynch and Pettitt (Lynch and Pettitt, J Chem Phys 1997, 107, 8594), which consists of the kinetic and potential energies for real and fractional particles as well as the kinetic and potential energy terms for material and heat reservoirs interacting with the system. We perform a nonlinear scaling of the potential energy parameters of the fractional particle, as well as its mass to vary the number of particles dynamically. On the basis of the equations of motion derived from this Hamiltonian, an algorithm has been proposed for MD simulation at constant chemical potential. The algorithm has been tested for the ideal gas, for the Lennard-Jones fluid over a wide range of temperatures and densities, and for water. The results for the low-density Lennard-Jones fluid are compared with the predictions from a truncated virial equation of state. In the case of the dense Lennard-Jones fluid and water our predicted results are compared with the results reported using other available methods for the calculation of the chemical potential. The method is also applied to the case of vapor-liquid coexistence point predictions. 相似文献
5.
Tetsuro Nagai George A. Pantelopulos Takuya Takahashi John E. Straub 《Journal of computational chemistry》2016,37(21):2017-2028
Simulated tempering (ST) is a generalized‐ensemble algorithm that employs trajectories exploring a range of temperatures to effectively sample rugged energy landscapes. When implemented using the molecular dynamics method, ST can require the use of short time steps for ensuring the stability of trajectories at high temperatures. To address this shortcoming, a mass‐scaling ST (MSST) method is presented in which the particle mass is scaled in proportion to the temperature. Mass scaling in the MSST method leads to velocity distributions that are independent of temperature and eliminates the need for velocity scaling after the accepted temperature updates that are required in conventional ST simulations. The homogeneity in time scales with changing temperature improves the stability of simulations and allows for the use of longer time steps at high temperatures. As a result, the MSST is found to be more efficient than the standard ST method, particularly for cases in which a large temperature range is employed. © 2016 Wiley Periodicals, Inc. 相似文献
6.
《Macromolecular theory and simulations》2018,27(4)
A new molecular modeling algorithm for conducting large‐scale molecular dynamics simulation studies of cross‐linked phenolic resins is developed using a united‐atom model. A phenol–formaldehyde polycondensation system is simulated by a pseudoreaction algorithm taking into consideration (i) the difference in the experimental reaction rate constants at ortho and para positions of phenolic units and (ii) the geometry of the reactants. To avoid formation of locally strained cross‐linked structures that can be generated in a typical cutoff‐distance‐based reaction scheme, a geometrical judgment constraint is applied in the reaction procedure. With this algorithm, cross‐linked network structures of phenolic resins with a maximum conversion (α) of 0.90 are obtained from 10 000 phenols. The density and the tensile modulus of the structure with α of 0.90 at 300 K are 1.2 g cm−3 and 5.4 GPa, respectively. This is in good agreement with experimental values. The strain‐free, highly cross‐linked network structures of phenolic resins exhibit a higher density and tensile modulus compared with structures generated in the absence of the geometrical cutoff. This result demonstrates that the geometrical judgment constraint can effectively avoid the formation of distorted and strained local structures and is necessary for accurate modeling of highly cross‐linked phenolic resins. 相似文献
7.
We propose an improvement of the replica-exchange and replica-permutation methods, which we call the replica sub-permutation method (RSPM). Instead of considering all permutations, this method uses a new algorithm referred to as sub-permutation to perform parameter transition. The RSPM succeeds in reducing the number of combinations between replicas and parameters without the loss of sampling efficiency. For comparison, we applied the replica sub-permutation, replica-permutation, and replica-exchange methods to a β-hairpin mini protein, chignolin, in explicit water. We calculated the transition ratio and number of tunneling events in the parameter space, the number of folding–unfolding events, the autocorrelation function, and the autocorrelation time as measures of sampling efficiency. The results indicate that among the three methods, the proposed RSPM is the most efficient in both parameter and conformational spaces. © 2019 Wiley Periodicals, Inc. 相似文献
8.
The biological function of HIV‐1 integrase (IN) is to integrate viral DNA into the host cell chromosome, and the specific binding of IN with viral DNA is a precondition for IN to function correctly. Beforehand, the binding mode of IN dimer (IN2) with the 27 bp segment of viral DNA before 3′ processing (3′‐P) was obtained via a molecular docking method. Based on the binding mode, the aim of this article was to explore the changes of motive mode and correlative motion for the IN2 and DNA systems after their binding through dynamical cross‐correlation map (DCCM) and principal component analysis (PCA). Finally, solvent effect during the association was analyzed briefly. The results show that there is a significantly increased positive correlation in the interface region between IN2 and viral DNA, and some obvious motive mode changes of the two systems (IN2 and DNA) were also observed after their binding. It was found that water molecules played an important role in the recognition between IN2 and viral DNA through analyzing the water‐mediated hydrogen bonds. 相似文献
9.
Implicit solvent representations, in general, and generalized Born models, in particular, provide an attractive way to reduce the number of interactions and degrees of freedom in a system. The instantaneous relaxation of the dielectric shielding provided by an implicit solvent model can be extremely efficient for high‐throughput and Monte Carlo studies, and a reduced system size can also remove a lot of statistical noise. Despite these advantages, it has been difficult for generalized Born implementations to significantly outperform optimized explicit‐water simulations due to more complex functional forms and the two extra interaction stages necessary to calculate Born radii and the derivative chain rule terms contributing to the force. Here, we present a method that uses a rescaling transformation to make the standard generalized Born expression a function of a single variable, which enables an efficient tabulated implementation on any modern CPU hardware. The total performance is within a factor 2 of simulations in vacuo. The algorithm has been implemented in Gromacs, including single‐instruction multiple‐data acceleration, for three different Born radius models and corresponding chain rule terms. We have also adapted the model to work with the virtual interaction sites commonly used for hydrogens to enable long‐time steps, which makes it possible to achieve a simulation performance of 0.86 μs/day for BBA5 with 1‐nm cutoff on a single quad‐core desktop processor. Finally, we have also implemented a set of streaming kernels without neighborlists to accelerate the non‐cutoff setup occasionally used for implicit solvent simulations of small systems. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010 相似文献
10.
11.
Ikeguchi M 《Journal of computational chemistry》2004,25(4):529-541
A partial rigid-body method of molecular dynamics simulations for proteins and membranes is presented. In this method, the symplectic integrator for rigid bodies is combined with the equations of motion for the NPT ensemble. The standard NPT ensemble is extended to the membrane-specific ensembles, the NPAT (constant normal pressure and lateral surface area of membranes and constant temperature) and NPgammaT (constant normal pressure and lateral surface tension of membranes and constant temperature) ensembles. By more than 30-ns simulations of aqueous proteins and hydrated lipid bilayers, the results of the partial rigid-body method demonstrated excellent conservation of total energy and consistent behavior with the traditional constraint method in terms of structural distribution and fluctuation of proteins and lipids. The efficient implementation of the partial rigid-body method in parallel computation is presented, which is shown to work well in large-scale molecular dynamics simulations. 相似文献
12.
Praprotnik M Hocevar S Hodoscek M Penca M Janezic D 《Journal of computational chemistry》2008,29(1):122-129
A force field of the triclinic framework of AlPO(4)-34, important in methanol-hydrocarbon conversion reactions, was developed using an empirical potential function. Molecular dynamics simulation of an AlPO(4)-34 triclinic framework segment of 1216 atoms, containing the template molecules isopropylamine and water, was performed with explicit consideration of atomic charges. The average RMS difference between instantaneous positions of the framework atoms during 1 ns simulation and their positions in the structure determined from single crystal X-ray diffraction was calculated, and the average structure of the flexible framework was determined. The computed Debye-Waller factors and simulated FTIR spectra are in good agreement with the experimental data. The new force field permits detailed molecular dynamics simulations of flexible, charged aluminophosphate molecular sieves which should lead to a better understanding of the catalytic processes and the crucial role played by templating molecules. 相似文献
13.
F1‐ATPase is an adenosine tri‐phosphate (ATP)‐driven rotary motor enzyme. We investigated the structural fluctuations and concerted motions of subunits in F1‐ATPase using molecular dynamics (MD) simulations. An MD simulation for the α3β3γ complex was carried out for 30 ns. Although large fluctuations of the N‐terminal domain observed in simulations of the isolated βE subunit were suppressed in the complex simulation, the magnitude of fluctuations in the C‐terminal domain was clearly different among the three β subunits (βE, βTP, and βDP). Despite fairly similar conformations of the βTP and βDP subunits, the βDP subunit exhibits smaller fluctuations in the C‐terminal domain than the βTP subunit due to their dissimilar interface configurations. Compared with the βTP subunit, the βDP subunit stably interacts with both the adjacent αDP and αE subunits. This sandwiched configuration in the βDP subunit leads to strongly correlated motions between the βDP and adjacent α subunits. The βDP subunit exhibits an extensive network of highly correlated motions with bound ATP and the γ subunit, as well as with the adjacent α subunits, suggesting that the structural changes occurring in the catalytically active βDP subunit can effectively induce movements of the γ subunit. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010 相似文献
14.
基于分子模拟技术煤焦分子模型构建 《燃料化学学报》2017,45(7):769-779
煤、焦是过程工业的重要原料。因此,有必要深入了解煤、焦分子结构以揭示其反应性。采用Materials Studio 7.0软件,从分子层次研究煤、焦的分子结构。根据已报道的文献,构建煤、焦的初始结构;基于分子力学原理对这些结构进行优化,使得模型物性与煤、焦物性相符;基于退火模拟算法对模型进行优化,从而使得模型密度、元素分析数据与真实值吻合;基于能量最小化原理,对煤、焦模型再次优化,从而获得其稳定、真实的分子构型。由计算结果发现,模型的估算密度、元素组成与已报道一致,说明构建的模型是有效、合理的;在模型优化过程中,相对于其他能量而言,库伦能和范德华能起着重要的作用。因此可以推断在煤、焦热加工过程中,弱键占据主要地位。另外,本文采用分子模拟技术构建煤、焦模型的方法对于构建其他复杂大分子结构有着重要的借鉴作用。 相似文献
15.
We have developed a time-reversible rigid-body (rRB) molecular dynamics algorithm in the isothermal-isobaric (NPT) ensemble. The algorithm is an extension of rigid-body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non-Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state-of-the-art constraint-dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time-step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS-MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL. 相似文献
16.
Andrzej J. Rzepiela Lars V. Schäfer Nicolae Goga H. Jelger Risselada Alex H. De Vries Siewert J. Marrink 《Journal of computational chemistry》2010,31(6):1333-1343
We present an algorithm to reconstruct atomistic structures from their corresponding coarse‐grained (CG) representations and its implementation into the freely available molecular dynamics (MD) program package GROMACS. The central part of the algorithm is a simulated annealing MD simulation in which the CG and atomistic structures are coupled via restraints. A number of examples demonstrate the application of the reconstruction procedure to obtain low‐energy atomistic structural ensembles from their CG counterparts. We reconstructed individual molecules in vacuo (NCQ tripeptide, dipalmitoylphosphatidylcholine, and cholesterol), bulk water, and a WALP transmembrane peptide embedded in a solvated lipid bilayer. The first examples serve to optimize the parameters for the reconstruction procedure, whereas the latter examples illustrate the applicability to condensed‐phase biomolecular systems. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010 相似文献
17.
We have developed an algorithm for sampling the conformational space of large flexible molecules in solution, which combines the molecular dynamics (MD) method and the three-dimensional reference interaction site model (3D-RISM) theory. The solvent-induced force acting on solute atoms was evaluated as the gradient of the solvation free energy with respect to the solute-atom coordinates. To enhance the computation speed, we have applied a multiple timestep algorithm based on the RESPA (Reversible System Propagator Algorithm) to the combined MD/3D-RISM method. By virtue of the algorithm, one can choose a longer timestep for renewing the solvent-induced force compared with that of the conformational update. To illustrate the present MD/3D-RISM simulation, we applied the method to a model of acetylacetone in aqueous solution. The multiple timestep algorithm succeeded in enhancing the computation speed by 3.4 times for this model case. Acetylacetone possesses an intramolecular hydrogen-bonding capability between the hydroxyl group and the carbonyl oxygen atom, and the molecule is significantly stabilized due to this hydrogen bond, especially in gas phase. The intramolecular hydrogen bond was kept intact during almost entire course of the MD simulation in gas phase, while in the aqueous solutions the bond is disrupted in a significant number of conformations. This result qualitatively agrees with the behavior on a free energy barrier lying upon the process for rotating a torsional degree of freedom of the hydroxyl group, where it is significantly reduced in aqueous solution by a cancellation between the electrostatic interaction and the solvation free energy. 相似文献
18.
Stern HA 《Journal of computational chemistry》2004,25(5):749-761
We review principles of non-Hamiltonian statistical mechanics and present a new set of equations and integration algorithm for isothermal-isobaric dynamics. The chief advantage of the present scheme is that it is somewhat simpler than previous methods. We perform numerical simulations to test the accuracy of the algorithm and compare its stability to that of a \"gold standard,\" a symplectic integrator for Hamiltonian dynamics of the same system. The stability of the isothermal-isobaric algorithm is comparable to the stability of the symplectic integrator. 相似文献
19.
Vivek S. Bharadwaj Nathaniel M. Eagan Nicholas M. Wang Prof. Matthew W. Liberatore Prof. C. Mark Maupin 《Chemphyschem》2015,16(13):2810-2817
Despite the importance of fatty‐acid methyl esters (FAMEs) as key components of various green solvents, detergents, plasticizers, and biodiesels, our understanding of these systems at the molecular level is limited. An enhanced molecular‐level perspective of FAMEs will enable a detailed analysis of the polymorph and crystallization phenomena that adversely impact flow properties at low temperatures. Presented here, is the parameterization and validation of a charge‐modified generalized amber force field (GAFF) for eight common FAMEs and two representative biodiesel mixtures. Our simulations accurately reproduce available experimental data (e.g. densities and self‐diffusivity coefficients) and their trends, with respect to temperature and degree of unsaturation. Structural analyses from our simulations provide a more detailed picture of liquid‐phase molecular ordering in FAMEs and confirm recent experimental hypotheses. This study provides a firm foundation to initiate further studies into the mechanisms that drive crystallization phenomena at the molecular level. 相似文献