BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Yuko Okamoto

Yuko Okamoto contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechBiological PhysicsBiomoleculesphysics.comp-phcond-mat.softphysics.chem-phastro-phgr-qc

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2020arXiv

Calculation of the Residual Entropy of Ice Ih by Monte Carlo simulation with the Combination of the Replica-Exchange Wang-Landau algorithm and Multicanonical Replica-Exchange Method

We estimated the residual entropy of ice Ih by the recently developed simulation protocol, namely, the combination of Replica-Exchange Wang-Landau algorithm and Multicanonical Replica-Exchange Method. We employed a model with the nearest neighbor interactions on the three-dimensional hexagonal lattice, which satisfied the ice rules in the ground state. The results showed that our estimate of the residual entropy is found to be within 0.038 % of series expansion estimate by Nagle and within 0.000077 % of PEPS algorithm by Vanderstraeten. In this article, we not only give our latest estimate of the residual entropy of ice Ih but also discuss the importance of the uniformity of a random number generator in MC simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Conformational Search Method for Protein Systems Using Genetic Crossover and Metropolis Criterion

Many proteins carry out their biological functions by forming the characteristic tertiary structures. Therefore, the search of the stable states of proteins by molecular simulations is important to understand their functions and stabilities. However, getting the stable state by conformational search is difficult, because the energy landscape of the system is characterized by many local minima separated by high energy barriers. In order to overcome this difficulty, various sampling and optimization methods for conformations of proteins have been proposed. In this study, we propose a new conformational search method for proteins by using genetic crossover and Metropolis criterion. We applied this method to an $α$-helical protein. The conformations obtained from the simulations are in good agreement with the experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Combination of genetic crossover and replica-exchange method for conformational search of protein systems

We combined the genetic crossover, which is one of the operations of genetic algorithm, and replica-exchange method in parallel molecular dynamics simulations. The genetic crossover and replica-exchange method can search the global conformational space by exchanging the corresponding parts between a pair of conformations of a protein. In this study, we applied this method to an $α$-helical protein, Trp-cage mini protein, which has 20 amino-acid residues. The conformations obtained from the simulations are in good agreement with the experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Predictions of tertiary stuctures of $α$-helical membrane proteins by replica-exchange method with consideration of helix deformations

We propose an improved prediction method of the tertiary structures of $α$-helical membrane proteins based on the replica-exchange method by taking into account helix deformations. Our method allows wide applications because transmembrane helices of native membrane proteins are often distorted. In order to test the effectiveness of the present method, we applied it to the structure predictions of glycophorin A and phospholamban. The results were in accord with experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Optimizations of force-field parameters for protein systems with the secondary-structure stability and instability

We propose a novel method for refining force-field parameters of protein systems. In this method, the agreement of the secondary-structure stability and instability between the protein conformations obtained by experiments and those obtained by molecular dynamics simulations is used as a criterion for the optimization of force-field parameters. As an example of the applications of the present method, we refined the force-field parameter set of the AMBER ff99SB force field by searching the torsion-energy parameter spaces of $ψ$ (N-C$^α$-C-N) and $ζ$ (C$^β$-C$^α$-C-N) of the backbone dihedral angles. We then performed folding simulations of $α$-helical and $β$-hairpin peptides, using the optimized force field. The results showed that the new force-field parameters gave structures more consistent with the experimental implications than the original AMBER ff99SB force field.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Amino-acid-dependent main-chain torsion-energy terms for protein systems

Many commonly used force fields for protein systems such as AMBER, CHARMM, GROMACS, OPLS, and ECEPP have amino-acid-independent force-field parameters of main-chain torsion-energy terms. Here, we propose a new type of amino-acid-dependent torsion-energy terms in the force fields. As an example, we applied this approach to AMBER ff03 force field and determined new amino-acid-dependent parameters for $ψ$ and $ψ'$ angles for each amino acid by using our optimization method, which is one of the knowledge-based approach. In order to test the validity of the new force-field parameters, we then performed folding simulations of $α$-helical and $β$-hairpin peptides, using the optimized force field. The results showed that the new force-field parameters gave structures more consistent with the experimental implications than the original AMBER ff03 force field.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Application of Simulated Tempering and Magnetizing to a Two-Dimensional Potts Model

We applied the simulated tempering and magnetizing (STM) method to the two-dimensional three-state Potts model in an external magnetic field in order to perform further investigations of the STM's applicability. The temperature as well as the external field are treated as dynamical variables updated during the STM simulations. After we obtained adequate information for several lattice sizes $L$ (up to $160\times 160$), we also performed a number of conventional canonical simulations of large lattices, especially in order to illustrate the crossover behavior of the Potts model in external field with increasing $L$. The temperature and external field for larger lattice size simulations were chosen by extrapolation of the detail information obtained by STM. We carefully analyzed the crossover scaling at the phase transitions with respect to the lattice size as well as the temperature and external field. The crossover behavior is clearly observed in the simulations in agreement with theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Free energy calculations of a proton transfer reaction by simulated tempering umbrella sampling first principles molecular dynamics simulations

A new simulated tempering method, which is referred to as simulated tempering umbrella sampling, for calculating the free energy of chemical reactions is proposed. First principles molecular dynamics simulations with this simulated tempering were performed in order to study the intramolecular proton transfer reaction of malonaldehyde in aqueous solution. Conformational sampling in reaction coordinate space can be easily enhanced with this method, and the free energy along a reaction coordinate can be calculated accurately. Moreover, the simulated tempering umbrella sampling provides trajectory data more efficiently than the conventional umbrella sampling method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimizations of protein force fields

In this Chapter we review our works on force fields for molecular simulations of protein systems. We first discuss the functional forms of the force fields and present some extensions of the conventional ones. We then present various methods for force-field parameter optimizations. Finally, some examples of our applications of these parameter optimization methods are given and they are compared with the results from the existing force-fields.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulated Tempering and Magnetizing: An Application of Two-Dimensional Simulated Tempering to Two-Dimensional Ising Model and Its Crossover

We performed two-dimensional simulated tempering (ST) simulations of the two-dimensional Ising model with different lattice sizes in order to investigate the two-dimensional ST's applicability to dealing with phase transitions and to study the crossover of critical scaling behavior. The external field, as well as the temperature, was treated as a dynamical variable updated during the simulations. Thus, this simulation can be referred to as "Simulated Tempering and Magnetizing (STM)." We also performed the "Simulated Magnetizing" (SM) simulations, in which the external field was considered as a dynamical variable and temperature was not. As has been discussed by previous studies, the ST method is not always compatible with first-order phase transitions. This is also true in the magnetizing process. Flipping of the entire magnetization did not occur in the SM simulations under $T_\mathrm{c}$ in large lattice-size simulations. However, the phase changed through the high temperature region in the STM simulations. Thus, the dimensional extension let us eliminate the difficulty of the first-order phase transitions and study wide area of the phase space. We then discuss how frequently parameter-updating attempts should be made for optimal convergence. The results favor frequent attempts. We finally study the crossover behavior of the phase transitions with respect to the temperature and external field. The crossover behavior was clearly observed in the simulations in agreement with the theoretical implications.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Generalized-Ensemble Algorithms for the Isobaric-Isothermal Ensemble

We present generalized-ensemble algorithms for isobaric-isothermal molecular simulations. In addition to the multibaric-multithermal algorithm and replica-exchange method for the isobaric-isothermal ensemble, which have already been proposed, we propose a simulated tempering method for this ensemble. We performed molecular dynamics simulations with these algorithms for an alanine dipeptide system in explicit water molecules to test the effectiveness of the algorithms. We found that these generalized-ensemble algorithms are all useful for conformational sampling of biomolecular systems in the isobaric-isothermal ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Replica-Exchange Molecular Dynamics Simulations for Various Constant Temperature Algorithms

In the replica-exchange molecular dynamics method, where constant-temperature molecular dynamics simulations are performed in each replica, one usually rescales the momentum of each particle after replica exchange. This rescaling method had previously been worked out only for the Gaussian constraint method. In this letter, we present momentum rescaling formulae for four other commonly used constant-temperature algorithms, namely, Langevin dynamics, Andersen algorithm, Nosé-Hoover thermostat, and Nosé-Poincaré thermostat. The effectiveness of these rescaling methods is tested with a small biomolecular system, and it is shown that proper momentum rescaling is necessary to obtain correct results in the canonical ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Replica-Exchange Method in van der Waals Radius Space: Overcoming Steric Restrictions for Biomolecules

We present a new type of the Hamiltonian replica-exchange method, in which not temperatures but the van der Waals radius parameter is exchanged. By decreasing the van der Waals radii that control spatial sizes of atoms, this Hamiltonian replica-exchange method overcomes the steric restrictions and energy barriers. Furthermore, the simulation based on this method escapes from the local-minimum free-energy states and realizes effective sampling in the conformational space. We applied this method to an alanine dipeptide in aqueous solution and showed the effectiveness of the method by comparing the results with those obtained from the conventional canonical method.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Residual Entropy of Ordinary Ice from Multicanonical Simulations

We introduce two simple models with nearest neighbor interactions on 3D hexagonal lattices. Each model allows one to calculate the residual entropy of ice I (ordinary ice) by means of multicanonical simulations. This gives the correction to the residual entropy derived by Linus Pauling in 1935. Our estimate is found to be within less than 0.1% of an analytical approximation by Nagle, which is an improvement of Pauling's result. We pose it as a challenge to experimentalists to improve on the accuracy of a 1936 measurement by Giauque and Stout by about one order of magnitude, which would allow one to identify corrections to Pauling's value unambiguously. It is straightforward to transfer our methods to other crystal systems.

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

Generalized Einstein Theory on Solar and Galactic Scales

We study a generalized Einstein theory with the following two criteria:{\it i}) on the solar scale, it must be consistent with the classical tests of general relativity, {\it ii}) on the galactic scale, the gravitational potential is a sum of Newtonian and Yukawa potentials so that it may explain the flat rotation curves of spiral galaxies. Under these criteria, we find that such a generalized Einstein action must include at least one scalar field and one vector field as well as the quadratic term of the scalar curvature.

0 citations0 reviewsNo code linkedOpen access