共查询到20条相似文献，搜索用时 15 毫秒
1.
《Journal of computational chemistry》2017,38(31):2671-2674
Parallel cascade selection molecular dynamics (PaCS‐MD) is an enhanced conformational sampling method for searching structural transition pathways from a given reactant to a product. Recently, a temperature‐aided PaCS‐MD (Vinod et al., Eur. Biophys. J. 2016, 45, 463) has been proposed as its extension, in which the temperatures were introduced as additional parameters in conformational resampling, whereas the temperature is fixed in the original PaCS‐MD. In the present study, temperature‐shuffled PaCS‐MD is proposed as a further extension of temperature‐aided PaCS‐MD in which the temperatures are shuffled among different replicas at the beginning of each cycle of conformational resampling. To evaluate their conformational sampling efficiencies, the original, temperature‐aided, and temperature‐shuffled PaCS‐MD were applied to a protein‐folding process of Trp‐cage, and their minimum computational costs to identify the native state were addressed. Through the evaluation, it was confirmed that temperature‐shuffled PaCS‐MD remarkably accelerated the protein‐folding process of Trp‐cage compared with the other methods. © 2017 Wiley Periodicals, Inc. 相似文献
2.
分子模拟在化学、物理、生物、材料等多学科的发展中起着越来越重要的作用。然而,受到当前计算机处理速度的限制,分子模拟计算所能够达到的时间尺度同实验或实际应用中要求的时间尺度相比还存在着巨大的差距。增强抽样方法的发展和应用可以有效地拓宽分子模拟所能研究体系的时间尺度,极大地提高分子模拟的热力学和动力学计算能力。本文中先简单介绍增强抽样方法的发展以及几类增强抽样方法的优缺点,然后重点介绍了我们研究组所发展的温度积分抽样方法(Integrated Tempering Sampling, ITS)的基本思路及其在蛋白质折叠研究中的应用。文章最后总结了增强抽样方法发展的新需求,同时也对此研究方向的广阔发展前景进行了展望。 相似文献
3.
《Journal of computational chemistry》2017,38(11):790-797
The folding processes of mini‐proteins (FSD‐EY/FBPWW28 domain) were computationally investigated by an enhanced conformational sampling method. Through the analyses of trajectories, these mini‐proteins had multiple folding pathways different from a simple two‐state folding, and the multiple folding processes were initiated by partial formations of secondary structures. These findings can be used to understand the folding of large proteins, that is, which secondary structures are partially folded in the early process, and how the remaining parts are sequentially folded. It is found that FSD‐EY (α/β topology) folds by a simple diffusion‐collision mechanism, while the folding process of the FBPWW28 domain (all‐β topology) requires a modification of the diffusion‐collision theory to adequately treat the coil‐sheet transition for the β sheet formation. © 2017 Wiley Periodicals, Inc. 相似文献
4.
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. 相似文献
5.
The recently developed \"temperature intervals with global exchange of replicas\" (TIGER2) algorithm is an efficient replica-exchange sampling algorithm that provides the freedom to specify the number of replicas and temperature levels independently of the size of the system and temperature range to be spanned, thus making it particularly well suited for sampling molecular systems that are considered to be too large to be sampled using conventional replica exchange methods. Although the TIGER2 method is empirical in nature, when appropriately applied it is able to provide sampling that satisfies the balance condition and closely approximates a Boltzmann-weighted ensemble of states. In this work, we evaluated the influence of factors such as temperature range, temperature spacing, replica number, and sampling cycle design on the accuracy of a TIGER2 simulation based on molecular dynamics simulations of alanine dipeptide in implicit solvent. The influence of these factors is further examined by calculating the properties of a complex system composed of the B1 immunoglobulin-binding domain of streptococcal protein G (protein G) in aqueous solution. The accuracy of a TIGER2 simulation is particularly sensitive to the maximum temperature level selected for the simulation. A method to determine the appropriate maximum temperature level to be used in a TIGER2 simulation is presented. 相似文献
6.
Ryuhei Harada Tomotake Nakamura Yasuteru Shigeta 《Journal of computational chemistry》2016,37(8):724-738
As an extension of the Outlier FLOODing (OFLOOD) method [Harada et al., J. Comput. Chem. 2015, 36, 763], the sparsity of the outliers defined by a hierarchical clustering algorithm, FlexDice, was considered to achieve an efficient conformational search as sparsity‐weighted “OFLOOD.” In OFLOOD, FlexDice detects areas of sparse distribution as outliers. The outliers are regarded as candidates that have high potential to promote conformational transitions and are employed as initial structures for conformational resampling by restarting molecular dynamics simulations. When detecting outliers, FlexDice defines a rank in the hierarchy for each outlier, which relates to sparsity in the distribution. In this study, we define a lower rank (first ranked), a medium rank (second ranked), and the highest rank (third ranked) outliers, respectively. For instance, the first‐ranked outliers are located in a given conformational space away from the clusters (highly sparse distribution), whereas those with the third‐ranked outliers are nearby the clusters (a moderately sparse distribution). To achieve the conformational search efficiently, resampling from the outliers with a given rank is performed. As demonstrations, this method was applied to several model systems: Alanine dipeptide, Met‐enkephalin, Trp‐cage, T4 lysozyme, and glutamine binding protein. In each demonstration, the present method successfully reproduced transitions among metastable states. In particular, the first‐ranked OFLOOD highly accelerated the exploration of conformational space by expanding the edges. In contrast, the third‐ranked OFLOOD reproduced local transitions among neighboring metastable states intensively. For quantitatively evaluations of sampled snapshots, free energy calculations were performed with a combination of umbrella samplings, providing rigorous landscapes of the biomolecules. © 2015 Wiley Periodicals, Inc. 相似文献
7.
Phillips JC Braun R Wang W Gumbart J Tajkhorshid E Villa E Chipot C Skeel RD Kalé L Schulten K 《Journal of computational chemistry》2005,26(16):1781-1802
8.
Daniel Beuve&#x;Mery Mehdi Rahman Philippe Ducarme Robert Brasseur 《Journal of computational chemistry》2001,22(2):172-177
We have built clusters of PC using up to 320 processors in parallel. The first cluster was built in 1996 from 64 computers. The parallel architectures are named RAMSES. In this article, we share our experience and discuss our choices. The architecture is composed of Pentium processors connected by a 100 Mbits/s Ethernet network. Because RAMSES are built of independent computers, the system is easily maintained, and can be modified by adding or removing machines. Two programs, which are dedicated to molecular modeling and were previously developed on single machines (OSIRIS and IMPALA), have been adapted. IMPALA models peptide and protein interaction with membranes, OSIRIS simulates an ab initio folding of globular proteins. After parallelization, each program consists of two parts: one is the human/RAMSES interface, and the second is the simulation per se. The parallel architecture is transparent for the users. Tests of the IMPALA and OSIRIS efficiencies are shown and discussed. © 2000 John Wiley & Sons, Inc. J Comput Chem 22: 172–177, 2001 相似文献
9.
Virtual screening of large libraries of small compounds requires fast and reliable automatic docking methods. In this article we present a parallel implementation of a genetic algorithm (GA) and the implementation of an enhanced genetic algorithm (EGA) with niching that lead to remarkable speedups compared to the original version AutoDock 3.0. The niching concept is introduced naturally by sharing genetic information between evolutions of subpopulations that run independently, each on one CPU. A unique set of additionally introduced search parameters that control this information flow has been obtained for drug‐like molecules based on the detailed study of three test cases of different complexity. The average docking time for one compound is of 8.6 s using eight R10,000 processors running at 200 MHz in an Origin 2000 computer. Different genetic algorithms with and without local search (LS) have been compared on an equal workload basis showing EGA/LS to be superior over all alternatives because it finds lower energy solutions faster and more often, particularly for high dimensionality problems. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1971–1982, 2001 相似文献
10.
Schmidt Am Busch M Lopes A Mignon D Simonson T 《Journal of computational chemistry》2008,29(7):1092-1102
11.
The kinetic and thermodynamic aspects of the helix-coil transition in polyalanine-based peptides have been studied at the ensemble level using a distributed computing network. This study builds on a previous report, which critically assessed the performance of several contemporary force fields in reproducing experimental measurements and elucidated the complex nature of helix-coil systems. Here we consider the effects of modifying backbone torsions and the scaling of noncovalent interactions. Although these elements determine the potential of mean force between atoms separated by three covalent bonds (and thus largely determine the local conformational distributions observed in simulation), we demonstrate that the interplay between these factors is both complex and force field dependent. We quantitatively assess the heliophilicity of several helix-stabilizing potentials as well as the changes in heliophilicity resulting from such modifications, which can \"make or break\" the accuracy of a given force field, and our findings suggests that future force field development may need to better consider effect that vary with peptide length. This report also serves as an example of the utility of distributed computing in analyzing and improving upon contemporary force fields at the level of absolute ensemble equilibrium, the next step in force field development. 相似文献
12.
J. Andrew McCammon 《Journal of computational chemistry》2015,36(20):1536-1549
Folding of four fast‐folding proteins, including chignolin, Trp‐cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred‐of‐microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2–2.1 Å of the native NMR or X‐ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second‐order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein‐folding studies. © 2015 Wiley Periodicals, Inc. 相似文献
13.
We critically examine a recently proposed convective replica exchange (cRE) method for enhanced sampling of protein conformation based on theoretical and numerical analysis. The results demonstrate that cRE and related replica exchange with guided annealing (RE‐GA) schemes lead to unbalanced exchange attempt probabilities and break detailed balance whenever the system undergoes slow conformational transitions (relative to the temperature diffusion timescale). Nonetheless, numerical simulations suggest that approximate canonical ensembles can be generated for systems with small conformational transition barriers. This suggests that RE‐GA maybe suitable for simulating intrinsically disordered proteins, an important class of newly recognized functional proteins. The efficacy of RE‐GA is demonstrated by calculating the conformational ensembles of intrinsically disordered kinase inducible domain protein. The results show that RE‐GA helps the protein to escape nonspecific compact states more efficiently and provides several fold speedups in generating converged and largely correct ensembles compared to the standard temperature RE. © 2014 Wiley Periodicals, Inc. 相似文献
14.
All‐atom sampling is a critical and compute‐intensive end stage to protein structural modeling. Because of the vast size and extreme ruggedness of conformational space, even close to the native structure, the high‐resolution sampling problem is almost as difficult as predicting the rough fold of a protein. Here, we present a combination of new algorithms that considerably speed up the exploration of very rugged conformational landscapes and are capable of finding heretofore hidden low‐energy states. The algorithm is based on a hierarchical workflow and can be parallelized on supercomputers with up to 128,000 compute cores with near perfect efficiency. Such scaling behavior is notable, as with Moore's law continuing only in the number of cores per chip, parallelizability is a critical property of new algorithms. Using the enhanced sampling power, we have uncovered previously invisible deficiencies in the Rosetta force field and created an extensive decoy training set for optimizing and testing force fields. © 2012 Wiley Periodicals, Inc. 相似文献
15.
Shaw DE 《Journal of computational chemistry》2005,26(13):1318-1328
Classical molecular dynamics simulations of biological macromolecules in explicitly modeled solvent typically require the evaluation of interactions between all pairs of atoms separated by no more than some distance R, with more distant interactions handled using some less expensive method. Performing such simulations for periods on the order of a millisecond is likely to require the use of massive parallelism. The extent to which such simulations can be efficiently parallelized, however, has historically been limited by the time required for interprocessor communication. This article introduces a new method for the parallel evaluation of distance-limited pairwise particle interactions that significantly reduces the amount of data transferred between processors by comparison with traditional methods. Specifically, the amount of data transferred into and out of a given processor scales as O(R(3/2)p(-1/2)), where p is the number of processors, and with constant factors that should yield a substantial performance advantage in practice. 相似文献
16.
We have investigated the folding pathway of the 36‐residue villin headpiece subdomain (HP‐36) by action‐derived molecular dynamics simulations. The folding is initiated by hydrophobic collapse, after which the concurrent formation of full tertiary structure and α‐helical secondary structure is observed. The collapse is observed to be associated with a couple of specific native contacts contrary to the conventional nonspecific hydrophobic collapse model. Stable secondary structure formation after the collapse suggests that the folding of HP‐36 follows neither the framework model nor the diffusion‐collision model. The C‐terminal helix forms first, followed by the N‐terminal helix positioned in its native orientation. The short middle helix is shown to form last. Signs for multiple folding pathways are also observed. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010 相似文献
17.
A Hamiltonian Replica‐Exchange Molecular Dynamics (REMD) simulation method has been developed that employs a two‐dimensional backbone and one‐dimensional side chain biasing potential specifically to promote conformational transitions in peptides. To exploit the replica framework optimally, the level of the biasing potential in each replica was appropriately adapted during the simulations. This resulted in both high exchange rates between neighboring replicas and improved occupancy/flow of all conformers in each replica. The performance of the approach was tested on several peptide and protein systems and compared with regular MD simulations and previous REMD studies. Improved sampling of relevant conformational states was observed for unrestrained protein and peptide folding simulations as well as for refinement of a loop structure with restricted mobility of loop flanking protein regions. © 2013 Wiley Periodicals, Inc. 相似文献
18.
By using distributed computing techniques and a supercluster of more than 20,000 processors we simulated folding of a 20-residue Trp Cage miniprotein in atomistic detail with implicit GB/SA solvent at a variety of solvent viscosities (gamma). This allowed us to analyze the dependence of folding rates on viscosity. In particular, we focused on the low-viscosity regime (values below the viscosity of water). In accordance with Kramers' theory, we observe approximately linear dependence of the folding rate on 1/gamma for values from 1-10(-1)x that of water viscosity. However, for the regime between 10(-4)-10(-1)x that of water viscosity we observe power-law dependence of the form k approximately gamma(-1/5). These results suggest that estimating folding rates from molecular simulations run at low viscosity under the assumption of linear dependence of rate on inverse viscosity may lead to erroneous results. 相似文献
19.
Konrad Hinsen Eric Pellegrini Sławomir Stachura Gerald R. Kneller 《Journal of computational chemistry》2012,33(25):2043-2048
We present a new version of the program package nMoldyn, which has been originally developed for a neutron‐scattering oriented analysis of molecular dynamics simulations of macromolecular systems (Kneller et al., Comput. Phys. Commun. 1995, 91, 191) and was later rewritten to include in‐depth time series analyses and a graphical user interface (Rog et al., J. Comput. Chem. 2003, 24, 657). The main improvement in this new version and the focus of this article are the parallelization of all the analysis algorithms for use on multicore desktop computers as well as distributed‐memory computing clusters. The parallelization is based on a task farming approach which maintains a simple program structure permitting easy modification and extension of the code to integrate new analysis methods. © 2012 Wiley Periodicals, Inc. 相似文献
20.