Ayori Mitsutake

Replica-exchange multicanonical and multicanonical replica-exchange Monte Carlo simulations of peptides. I. Formulation and benchmark test

The replica-exchange multicanonical algorithm and the multicanonical replica-exchange method for molecular dynamics simulations have recently been developed. In the former method the multicanonical weight factor is determined from a short replica-exchange simulation with the multiple-histogram reweighting techniques. A long multicanonical production run with high statistics is then performed with this weight factor. In this method, the process of determining the multicanonical weight factor is faster and simpler than that in the usual iterative determination. The multicanonical replica-exchange method is a further extension of the first in which a replica-exchange multicanonical simulation is performed with a small number of replicas. In this paper, we give the formulations of these two methods for Monte Carlo simulations and demonstrate the effectiveness of these algorithms for a penta peptide in the gas phase.

Replica-exchange multicanonical and multicanonical replica-exchange Monte Carlo simulations of peptides. II. Application to a more complex system

In Paper I of this series the formulations of the replica-exchange multicanonical algorithm and the multicanonical replica-exchange method for Monte Carlo versions have been presented. The effectiveness of these algorithms were then tested with the system of a penta peptide, Met-enkephalin, in the gas phase. In this article the detailed comparisons of performances of these algorithms together with the regular replica-exchange method are made, taking a more complex system of a 17-residue helical peptide. It is shown that these two new algorithms are more efficient than the regular replica-exchange method.

Multidimensional generalized-ensemble algorithms for complex systems

We give general formulations of the multidimensional multicanonical algorithm, simulated tempering, and replica-exchange method. We generalize the original potential energy function E(0) by adding any physical quantity V of interest as a new energy term. These multidimensional generalized-ensemble algorithms then perform a random walk not only in E(0) space but also in V space. Among the three algorithms, the replica-exchange method is the easiest to perform because the weight factor is just a product of regular Boltzmann-like factors, while the weight factors for the multicanonical algorithm and simulated tempering are not a priori known. We give a simple procedure for obtaining the weight factors for these two latter algorithms, which uses a short replica-exchange simulation and the multiple-histogram reweighting techniques. As an example of applications of these algorithms, we have performed a two-dimensional replica-exchange simulation and a two-dimensional simulated-tempering simulation using an alpha-helical peptide system. From these simulations, we study the helix-coil transitions of the peptide in gas phase and in aqueous solution.

Enhanced sampling algorithms.

In biomolecular systems (especially all-atom models) with many degrees of freedom such as proteins and nucleic acids, there exist an astronomically large number of local-minimum-energy states. Conventional simulations in the canonical ensemble are of little use, because they tend to get trapped in states of these energy local minima. Enhanced conformational sampling techniques are thus in great demand. A simulation in generalized ensemble performs a random walk in potential energy space and can overcome this difficulty. From only one simulation run, one can obtain canonical-ensemble averages of physical quantities as functions of temperature by the single-histogram and/or multiple-histogram reweighting techniques. In this article we review uses of the generalized-ensemble algorithms in biomolecular systems. Three well-known methods, namely, multicanonical algorithm, simulated tempering, and replica-exchange method, are described first. Both Monte Carlo and molecular dynamics versions of the algorithms are given. We then present various extensions of these three generalized-ensemble algorithms. The effectiveness of the methods is tested with short peptide and protein systems.

Replica-exchange extensions of simulated tempering method

In this paper we consider combinations of two well-known generalized-ensemble algorithms, namely, simulated tempering and replica-exchange method. We discuss two examples of such combinations. One is the replica-exchange simulated tempering and the other is the simulated tempering replica-exchange method. In the former method, a short replica-exchange simulation is first performed and the simulated tempering weight factor is obtained by the multiple-histogram reweighting techniques. This process of simulated tempering weight factor determination is faster and simpler than that in the usual iterative process. A long simulated tempering production run is then performed with this weight factor. The latter method is a further extension of the former in which a simulated tempering replica-exchange simulation is performed with a small number of replicas. These algorithms are particularly useful for studying frustrated systems with rough energy landscape. We give the formulations of these two methods in detail and demonstrate their effectiveness taking the example of the system of a 17-residue helical peptide.

Simulated-tempering replica-exchange method for the multidimensional version

In this article, the general formulation of the multidimensional simulated-tempering replica-exchange method is described. In previous works, the one-dimensional replica-exchange simulated-tempering and simulated-tempering replica-exchange methods were developed. For the former method, the weight factor of the one-dimensional simulated tempering is determined by a short replica-exchange simulation and multiple-histogram reweighing techniques. For the latter method, the production run is a replica-exchange simulation with a few replicas not in the canonical ensembles but in the simulated-tempering ensembles which cover wide ranges of temperature. Recently, the general formulation of the multidimensional replica-exchange simulated tempering was presented. In this article, the extension of the simulated-tempering replica-exchange method for the multidimensional version is given. As an example of applications of the algorithm, a two-dimensional replica-exchange simulation and two simulated-tempering replica-exchange simulations have been performed. Here, an alpha-helical peptide system with a model solvent has been used for the applications.

Relaxation mode analysis of a peptide system: Comparison with principal component analysis

This article reports the first attempt to apply the relaxation mode analysis method to a simulation of a biomolecular system. In biomolecular systems, the principal component analysis is a well-known method for analyzing the static properties of fluctuations of structures obtained by a simulation and classifying the structures into some groups. On the other hand, the relaxation mode analysis has been used to analyze the dynamic properties of homopolymer systems. In this article, a long Monte Carlo simulation of Met-enkephalin in gas phase has been performed. The results are analyzed by the principal component analysis and relaxation mode analysis methods. We compare the results of both methods and show the effectiveness of the relaxation mode analysis.

New Computational Approach to Determine Liquid–Solid Phase Equilibria of Water Confined to Slit Nanopores

We devise a new computational approach to compute solid-liquid phase equilibria of confined fluids. Specifically, we extend the multibaric-multithermal ensemble method with an anisotropic pressure control to achieve the solid-liquid phase equilibrium for confined water inside slit nanopores (with slit width h ranging from 5.4 Å to 7.2 Å). A unique feature of this multibaric-multithermal ensemble is that the freezing points of confined water can be determined from the heat-capacity peaks. The new approach has been applied to compute the freezing point of two monolayer ices, namely, a high-density flat rhombic monolayer ice (HD-fRMI) and a high-density puckered rhombic monolayer ice (HD-pRMI) observed in our simulation. We find that the liquid-to-solid transition temperature (or the freezing point) of HD-pRMI is dependent on the slit width h, whereas that of HD-fRMI is nearly independent of the h.

Liquid-solid and solid-solid phase transition of monolayer water: High-density rhombic monolayer ice

Liquid-solid and solid-solid phase transitions of a monolayer water confined between two parallel hydrophobic surfaces are studied by molecular dynamics simulations. The solid phase considered is the high-density rhombic monolayer ice. Based on the computed free energy surface, it is found that at a certain width of the slit nanopore, the monolayer water exhibits not only a high freezing point but also a low energy barrier to crystallization. Moreover, through analyzing the oxygen-hydrogen-oxygen angle distribution and oxygen-hydrogen radial distribution, the high-density monolayer ice is classified as either a flat ice or a puckered ice. The transition between a flat ice and a puckered ice reflects a trade-off between the water-wall interactions and the electrostatic interactions among water molecules.

Molecular dynamics of C-peptide of ribonuclease A studied by replica-exchange Monte Carlo method and diffusion theory

Generalized-ensemble algorithm and diffusion theory have been combined in order to compute the dynamical properties monitored by nuclear magnetic resonance (NMR) experiments from efficient and reliable evaluation of statistical averages. Replica-exchange Monte Carlo simulations have been performed with a C-peptide analogue of ribonuclease A, and the Smoluchowski diffusion equation has been applied. A fairly good agreement between the calculated and measured 1H-NOESY NMR cross peaks has been obtained. The combination of these advanced and continuously improving statistical tools allows the calculation of a wide variety of dynamical properties routinely obtained by experiments.