Source author record

Macoto Kikuchi

Macoto Kikuchi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

9 published item(s)

preprint2022arXiv

Evolution enhances mutational robustness and suppresses the emergence of a new phenotype: A new computational approach for studying evolution

The aim of this paper is two-fold. First, we propose a new computational method to investigate the particularities of evolution. Second, we apply this method to a model of gene regulatory networks (GRNs) and explore the evolution of mutational robustness and bistability. Living systems have developed their functions through evolutionary processes. To understand the particularities of this process theoretically, evolutionary simulation (ES) alone is insufficient because the outcomes of ES depend on evolutionary pathways. We need a reference system for comparison. An appropriate reference system for this purpose is an ensemble of the randomly sampled genotypes. However, generating high-fitness genotypes by simple random sampling is difficult because such genotypes are rare. In this study, we used the multicanonical Monte Carlo method developed in statistical physics to construct a reference ensemble of GRNs and compared it with the outcomes of ES. We obtained the following results. First, mutational robustness was significantly higher in ES than in the reference ensemble at the same fitness level. Second, the emergence of a new phenotype, bistability, was delayed in evolution. Third, the bistable group of GRNs contains many mutationally fragile GRNs compared with those in the non-bistable group. This suggests that the delayed emergence of bistability is a consequence of the mutation-selection mechanism.

preprint2016arXiv

The interplay of intrinsic disorder and macromolecular crowding on α-synuclein fibril formation

α-synuclein (α-syn) is an intrinsically disordered protein which is considered to be one of the causes of Parkinson's disease. This protein forms amyloid fibrils when in a highly concentrated solution. The fibril formation of α-syn is induced not only by increases in α-syn concentration but also by macromolecular crowding. In order to investigate the coupled effect of the intrinsic disorder of α-syn and macromolecular crowding, we construct a lattice gas model of α-syn in contact with a crowding agent reservoir based on statistical mechanics. The main assumption is that α-syn can be expressed as coarse-grained particles with internal states coupled with effective volume; and disordered states are modeled by larger particles with larger internal entropy than other states. Thanks to the simplicity of the model, we can exactly calculate the number of conformations of crowding agents, and this enables us to prove that the original grand canonical ensemble with a crowding agent reservoir is mathematically equivalent to a canonical ensemble without crowding agents. In this expression, the effect of macromolecular crowding is absorbed in the internal entropy of disordered states; it is clearly shown that the crowding effect reduces the internal entropy. Based on Monte Carlo simulation, we provide scenarios of crowding-induced fibril formation. We also discuss the recent controversy over the existence of helically folded tetramers of α-syn, and suggest that macromolecular crowding is the key to resolving the controversy.

preprint2013arXiv

How to estimate the number of self-avoiding walks over 10^100? Use random walks

Counting the number of N-step self-avoiding walks (SAWs) on a lattice is one of the most difficult problems of enumerative combinatorics. Once we give up calculating the exact number of them, however, we have a chance to apply powerful computational methods of statistical mechanics to this problem. In this paper, we develop a statistical enumeration method for SAWs using the multicanonical Monte Carlo method. A key part of this method is to expand the configuration space of SAWs to random walks, the exact number of which is known. Using this method, we estimate a number of N-step SAWs on a square lattice, c_N, up to N=256. The value of c_256 is 5.6(1)*10^108 (the number in the parentheses is the statistical error of the last digit) and this is larger than one googol (10^100).

preprint2013arXiv

Multicanonical simulation of the Domb-Joyce model and the Go model: new enumeration methods for self-avoiding walks

We develop statistical enumeration methods for self-avoiding walks using a powerful sampling technique called the multicanonical Monte Carlo method. Using these methods, we estimate the numbers of the two dimensional N-step self-avoiding walks up to N=256 with statistical errors. The developed methods are based on statistical mechanical models of paths which include self-avoiding walks. The criterion for selecting a suitable model for enumerating self-avoiding walks is whether or not the configuration space of the model includes a set for which the number of the elements can be exactly counted. We call this set a scale fixing set. We selected the following two models which satisfy the criterion: the Go model for lattice proteins and the Domb-Joyce model for generalized random walks. There is a contrast between these two models in the structures of the configuration space. The configuration space of the Go model is defined as the universal set of self-avoiding walks, and the set of the ground state conformation provides a scale fixing set. On the other hand, the configuration space of the Domb-Joyce model is defined as the universal set of random walks which can be used as a scale fixing set, and the set of the ground state conformation is the same as the universal set of self-avoiding walks. From the perspective of enumeration performance, we conclude that the Domb-Joyce model is the better of the two. The reason for the performance difference is partly explained by the existence of the first-order phase transition of the Go model.

preprint2013arXiv

Robustness Leads Close to the Edge of Chaos in Coupled Map Networks: toward the understanding of biological networks

Dynamics in biological networks are in general robust against several perturbations. We investigate a coupled map network as a model motivated by gene regulatory networks and design systems which are robust against phenotypic perturbations (perturbations in dynamics), as well as systems which are robust against mutation (perturbations in network structure). To achieve such a design, we apply a multicanonical Monte Carlo method. Analysis based on the maximum Lyapunov exponent and parameter sensitivity shows that systems with marginal stability, which are regarded as systems at the edge of chaos, emerge when robustness against network perturbations is required. This emergence of the edge of chaos is a self-organization phenomenon and does not need a fine tuning of parameters.

preprint2013arXiv

Structural flexibility of intrinsically disordered proteins induces stepwise target recognition

An intrinsically disordered protein (IDP) lacks a stable three-dimensional structure, while it folds into a specific structure when it binds to a target molecule. In some IDP-target complexes, not all target binding surfaces are exposed on the outside, and intermediate states are observed in their binding processes. We consider that stepwise target recognition via intermediate states is a characteristic of IDP binding to targets with "hidden" binding sites. To investigate IDP binding to hidden target binding sites, we constructed an IDP lattice model based on the HP model. In our model, the IDP is modeled as a chain and the target is modeled as a highly coarse-grained object. We introduced motion and internal interactions to the target to hide its binding sites. In the case of unhidden binding sites, a two-state transition between the free states and a bound state is observed, and we consider that this represents coupled folding and binding. Introducing hidden binding sites, we found an intermediate bound state in which the IDP forms various structures to temporarily stabilize the complex. The intermediate state provides a scaffold for the IDP to access the hidden binding site. We call this process multiform binding. We conclude that structural flexibility of IDPs enables them to access hidden binding sites, and this is a functional advantage of IDPs.

preprint1996arXiv

Statistical Dependence Analysis

We review our recent studies on the dynamical correlations in MC simulations from the view point of the statistical dependence. Attentions are paid to the reduction of the statistical degrees of freedom for correlated data. Possible biases on several cumulants, such as the susceptibility and the Binder number due to finite MC length are discussed. A new method for calculating the equilibrium relaxation time from the analysis of the statistical dependence is presented. We apply it to the critical dynamics of the Ising model to estimate the dynamical critical exponent accurately.

preprint1996arXiv

Universal Finite-Size-Scaling Functions

The idea of universal finite-size-scaling functions of the Ising model is tested by Monte Carlo simulations for various lattices. Not only regular lattices such as the square lattice but quasiperiodic lattices such as the Penrose lattice are treated. We show that the finite-size-scaling functions of the order parameter for various lattices are collapsed on a single curve by choosing two nonuniversal scaling metric factors. We extend the idea of the universal finite-size-scaling functions to the order-parameter distribution function. We pay attention to the effects of boundary conditions. Keywords: Universal Finite-Size-Scaling Function; Ising Model; Order-Parameter Probability Distribution Function.

preprint1993arXiv

Recycle of random sequences

The correlation between a random sequence and its transformed sequences is studied. In the case of a permutation operation or, in other word, the shuffling operation, it is shown that the correlation can be so small that the sequences can be regarded as independent random sequences. The applications to the Monte Carlo simulations are also given. This method is especially useful in the Ising Monte Carlo simulation.