Generation of dynamic three-dimensional genome structure through phase separation of chromatin.

Three-dimensional genome structure and dynamics play critical roles in regulating DNA functions. Flexible chromatin structure and movements suggested that the genome is dynamically phase separated to form A (active) and B (inactive) compartments in interphase nuclei. Here, we examine this hypothesis by developing a polymer model of the whole genome of human cells and assessing the impact of phase separation on genome structure. Upon entry to the G1 phase, the simulated genome expanded according to heterogeneous repulsion among chromatin chains, which moved chromatin heterogeneously, inducing phase separation of chromatin. This repulsion-driven phase separation quantitatively reproduces the experimentally observed chromatin domains, A/B compartments, lamina-associated domains, and nucleolus-associated domains, consistently explaining nuclei of different human cells and predicting their dynamic fluctuations. We propose that phase separation induced by heterogeneous repulsive interactions among chromatin chains largely determines dynamic genome organization.

Role of the reaction-structure coupling in temperature compensation of the KaiABC circadian rhythm.

When the mixture solution of cyanobacterial proteins, KaiA, KaiB, and KaiC, is incubated with ATP in vitro, the phosphorylation level of KaiC shows stable oscillations with the temperature-compensated circadian period. Elucidating this temperature compensation is essential for understanding the KaiABC circadian clock, but its mechanism has remained a mystery. We analyzed the KaiABC temperature compensation by developing a theoretical model describing the feedback relations among reactions and structural transitions in the KaiC molecule. The model showed that the reduced structural cooperativity should weaken the negative feedback coupling among reactions and structural transitions, which enlarges the oscillation amplitude and period, explaining the observed significant period extension upon single amino-acid residue substitution. We propose that an increase in thermal fluctuations similarly attenuates the reaction-structure feedback, explaining the temperature compensation in the KaiABC clock. The model explained the experimentally observed responses of the oscillation phase to the temperature shift or the ADP-concentration change and suggested that the ATPase reactions in the CI domain of KaiC affect the period depending on how the reaction rates are modulated. The KaiABC clock provides a unique opportunity to analyze how the reaction-structure coupling regulates the system-level synchronized oscillations of molecules.

Organization of fast and slow chromatin revealed by single-nucleosome dynamics.

Understanding chromatin organization and dynamics is important, since they crucially affect DNA functions. In this study, we investigate chromatin dynamics by statistically analyzing single-nucleosome movement in living human cells. Bimodal nature of the mean square displacement distribution of nucleosomes allows for a natural categorization of the nucleosomes as fast and slow. Analyses of the nucleosome-nucleosome correlation functions within these categories along with the density of vibrational modes show that the nucleosomes form dynamically correlated fluid regions (i.e., dynamic domains of fast and slow nucleosomes). Perturbed nucleosome dynamics by global histone acetylation or cohesin inactivation indicate that nucleosome-nucleosome interactions along with tethering of chromatin chains organize nucleosomes into fast and slow dynamic domains. A simple polymer model is introduced, which shows the consistency of this dynamic domain picture. Statistical analyses of single-nucleosome movement provide rich information on how chromatin is dynamically organized in a fluid manner in living cells.

The Structural Rule Distinguishing a Superfold: A Case Study of Ferredoxin Fold and the Reverse Ferredoxin Fold.

Superfolds are folds commonly observed among evolutionarily unrelated multiple superfamilies of proteins. Since discovering superfolds almost two decades ago, structural rules distinguishing superfolds from the other ordinary folds have been explored but remained elusive. Here, we analyzed a typical superfold, the ferredoxin fold, and the fold which reverses the N to C terminus direction from the ferredoxin fold as a case study to find the rule to distinguish superfolds from the other folds. Though all the known structural characteristics for superfolds apply to both the ferredoxin fold and the reverse ferredoxin fold, the reverse fold has been found only in a single superfamily. The database analyses in the present study revealed the structural preferences of αß- and ßα-units; the preferences separate two α-helices in the ferredoxin fold, preventing their collision and stabilizing the fold. In contrast, in the reverse ferredoxin fold, the preferences bring two helices near each other, inducing structural conflict. The Rosetta folding simulations suggested that the ferredoxin fold is physically much more realizable than the reverse ferredoxin fold. Therefore, we propose that minimal structural conflict or minimal frustration among secondary structures is the rule to distinguish a superfold from ordinary folds. Intriguingly, the database analyses revealed that a most stringent structural rule in proteins, the right-handedness of the ßαß-unit, is broken in a set of structures to prevent the frustration, suggesting the proposed rule of minimum frustration among secondary structural units is comparably strong as the right-handedness rule of the ßαß-unit.

Effects of a single dose of tablets containing lactononadecapeptide on cognitive function in healthy adults: a randomized, double-blind, cross-over, placebo-controlled trial.

Lactononadecapeptide (LNDP; NIPPLTQTPVVVPPFLQPE) is a memory-improving peptide. The current study aimed to determine the effects of a single dose of tablets containing LNDP on cognitive function in healthy Japanese men aged 30-59 years. A randomized, double-blind, cross-over, placebo-controlled trial was conducted in participants randomly assigned to receive LNDP or placebo tablets. The Uchida-Kraepelin test was used to induce cognitive load in participants as a model of work load. Cognitive function was evaluated using the Japanese version of the CNS Vital Signs. Composite memory and verbal memory were significantly higher following consumption of LNDP than placebo tablets. Carryover effects were observed in attention and concentration domains so that period 1 data was analyzed. LNDP consumption led to higher processing speed, executive function, and cognitive flexibility than placebo. Thus, supplementation with a single dose of LNDP tablets may improve cognitive functions including memory, attention, concentration, and information processing in daily life.

Orally active anti-hypertensive peptides found based on enteroendocrine cell responses to a dipeptide library.

We previously reported that an orally administered dipeptide, Arg-Phe (RF), which causes enteroendocrine cell responses, lowered blood pressure in spontaneously hypertensive rats (SHRs). In this study, we found that Phe-Trp (FW), induced the most potent enteroendocrine cell responses out of total 338 dipeptides. An FW analogue, Phe-Trp-Gly-Lys (FWGK), which was effectively produced by tryptic digestion of bovine serum albumin, decreased blood pressure after oral administration. The minimum effective dose of FWGK (50 µg/kg) was 1/300 of that of RF (15 mg/kg). FWGK stimulated cholecystokinin (CCK) secretion in the enteroendocrine cells and exhibited vasorelaxing and antihypertensive effects via the CCK1 system.

Heterogeneous Spatial Distribution of Transcriptional Activity in Budding Yeast Nuclei.

Recent microscopic and simulation studies have shown that the genome structure fluctuates dynamically in the nuclei of budding yeast Saccharomyces cerevisiae. This genome-wide movement should lead to the fluctuations of individual genes in their territorial regions. This raises an intriguing question of whether the resulting distribution of genes is correlated to their transcriptional activity. An effective method for examining this correlation is to analyze how the spatial distribution of genes and their transcriptional activity are modified by mutation. In this study, we analyzed the modification observed in a budding yeast mutant in which genes necessary for anchoring telomeres to the nuclear envelope, yku70 and esc1, are silenced. Taddei et al. reported that 60 genes are clearly misregulated by this mutation, with 28 and 32 genes downregulated and upregulated, respectively. We calculated the probability density maps of the misregulated genes using a model of dynamical movement of the yeast genome in both wild-type (WT) and yku70 esc1 mutant and showed that the density of downregulated genes is larger near the nucleolus, whereas the density of upregulated genes is larger at the opposite side of the nucleus. By comparing these genes with those highly (top 200 of transcriptome) and lowly (bottom 200) expressed, we showed that the simulated distribution of 28 downregulated (12 out of 32 upregulated) genes has a distinctly larger overlap with the distribution of lowly (highly) expressed genes in the mutant than in the WT. The remaining 20 upregulated genes are localized near the nuclear envelope both in the WT and in the mutant. These results showed that the transcriptional level of genes is affected by their spatial distribution, thus highlighting the importance of the structural regulation in the yeast genome.

Folding pathway of a multidomain protein depends on its topology of domain connectivity.

How do the folding mechanisms of multidomain proteins depend on protein topology? We addressed this question by developing an Ising-like structure-based model and applying it for the analysis of free-energy landscapes and folding kinetics of an example protein, Escherichia coli dihydrofolate reductase (DHFR). DHFR has two domains, one comprising discontinuous N- and C-terminal parts and the other comprising a continuous middle part of the chain. The simulated folding pathway of DHFR is a sequential process during which the continuous domain folds first, followed by the discontinuous domain, thereby avoiding the rapid decrease in conformation entropy caused by the association of the N- and C-terminal parts during the early phase of folding. Our simulated results consistently explain the observed experimental data on folding kinetics and predict an off-pathway structural fluctuation at equilibrium. For a circular permutant for which the topological complexity of wild-type DHFR is resolved, the balance between energy and entropy is modulated, resulting in the coexistence of the two folding pathways. This coexistence of pathways should account for the experimentally observed complex folding behavior of the circular permutant.

Eddy current and coupled landscapes for nonadiabatic and nonequilibrium complex system dynamics.

Physical and biological systems are often involved with coupled processes of different time scales. In the system with electronic and atomic motions, for example, the interplay between the atomic motion along the same energy landscape and the electronic hopping between different landscapes is critical: the system behavior largely depends on whether the intralandscape motion is slower (adiabatic) or faster (nonadiabatic) than the interlandscape hopping. For general nonequilibrium dynamics where Hamiltonian or energy function is unknown a priori, the challenge is how to extend the concepts of the intra- and interlandscape dynamics. In this paper we establish a theoretical framework for describing global nonequilibrium and nonadiabatic complex system dynamics by transforming the coupled landscapes into a single landscape but with additional dimensions. On this single landscape, dynamics is driven by gradient of the potential landscape, which is closely related to the steady-state probability distribution of the enlarged dimensions, and the probability flux, which has a curl nature. Through an example of a self-regulating gene circuit, we show that the curl flux has dramatic effects on gene regulatory dynamics. The curl flux and landscape framework developed here are easy to visualize and can be used to guide further investigation of physical and biological nonequilibrium systems.

Coupling of lever arm swing and biased Brownian motion in actomyosin.

An important unresolved problem associated with actomyosin motors is the role of Brownian motion in the process of force generation. On the basis of structural observations of myosins and actins, the widely held lever-arm hypothesis has been proposed, in which proteins are assumed to show sequential structural changes among observed and hypothesized structures to exert mechanical force. An alternative hypothesis, the Brownian motion hypothesis, has been supported by single-molecule experiments and emphasizes more on the roles of fluctuating protein movement. In this study, we address the long-standing controversy between the lever-arm hypothesis and the Brownian motion hypothesis through in silico observations of an actomyosin system. We study a system composed of myosin II and actin filament by calculating free-energy landscapes of actin-myosin interactions using the molecular dynamics method and by simulating transitions among dynamically changing free-energy landscapes using the Monte Carlo method. The results obtained by this combined multi-scale calculation show that myosin with inorganic phosphate (Pi) and ADP weakly binds to actin and that after releasing Pi and ADP, myosin moves along the actin filament toward the strong-binding site by exhibiting the biased Brownian motion, a behavior consistent with the observed single-molecular behavior of myosin. Conformational flexibility of loops at the actin-interface of myosin and the N-terminus of actin subunit is necessary for the distinct bias in the Brownian motion. Both the 5.5-11 nm displacement due to the biased Brownian motion and the 3-5 nm displacement due to lever-arm swing contribute to the net displacement of myosin. The calculated results further suggest that the recovery stroke of the lever arm plays an important role in enhancing the displacement of myosin through multiple cycles of ATP hydrolysis, suggesting a unified movement mechanism for various members of the myosin family.

VS-APPLE: A Virtual Screening Algorithm Using Promiscuous Protein-Ligand Complexes.

As the number of structurally resolved protein-ligand complexes increases, the ligand-binding pockets of many proteins have been found to accommodate multiple different compounds. Effective use of these structural data is important for developing virtual screening (VS) methods that identify bioactive compounds. Here, we introduce a VS method, VS-APPLE (Virtual Screening Algorithm using Promiscuous Protein-Ligand complExes), based on promiscuous protein-ligand binding structures. In VS-APPLE, multiple ligands bound to a pocket are combined into a query template for screening. Both the structural match between a test compound and the multiple-ligand template and the possible collisions between the test compound and the target protein are evaluated by an efficient geometric hashing method. The performance of VS-APPLE was examined on a filtered, clustered version of the Directory of Useful Decoys data set. In Area Under the Curve analyses of this data set, VS-APPLE outperformed several popular screening programs. Judging from the performance of VS-APPLE, the structural data of promiscuous protein-ligand bindings could be further analyzed and exploited for developing VS methods.

Multiple coupled landscapes and non-adiabatic dynamics with applications to self-activating genes.

Many physical, chemical and biochemical systems (e.g. electronic dynamics and gene regulatory networks) are governed by continuous stochastic processes (e.g. electron dynamics on a particular electronic energy surface and protein (gene product) synthesis) coupled with discrete processes (e.g. hopping among different electronic energy surfaces and on and off switching of genes). One can also think of the underlying dynamics as the continuous motion on a particular landscape and discrete hoppings among different landscapes. The main difference of such systems from the intra-landscape dynamics alone is the emergence of the timescale involved in transitions among different landscapes in addition to the timescale involved in a particular landscape. The adiabatic limit when inter-landscape hoppings are fast compared to continuous intra-landscape dynamics has been studied both analytically and numerically, but the analytical treatment of the non-adiabatic regime where the inter-landscape hoppings are slow or comparable to continuous intra-landscape dynamics remains challenging. In this study, we show that there exists mathematical mapping of the dynamics on 2(N) discretely coupled N continuous dimensional landscapes onto one single landscape in 2N dimensional extended continuous space. On this 2N dimensional landscape, eddy current emerges as a sign of non-equilibrium non-adiabatic dynamics and plays an important role in system evolution. Many interesting physical effects such as the enhancement of fluctuations, irreversibility, dissipation and optimal kinetics emerge due to non-adiabaticity manifested by the eddy current illustrated for an N = 1 self-activator. We further generalize our theory to the N-gene network with multiple binding sites and multiple synthesis rates for discretely coupled non-equilibrium stochastic physical and biological systems.

Time scales in epigenetic dynamics and phenotypic heterogeneity of embryonic stem cells.

A remarkable feature of the self-renewing population of embryonic stem cells (ESCs) is their phenotypic heterogeneity: Nanog and other marker proteins of ESCs show large cell-to-cell variation in their expression level, which should significantly influence the differentiation process of individual cells. The molecular mechanism and biological implication of this heterogeneity, however, still remain elusive. We address this problem by constructing a model of the core gene-network of mouse ESCs. The model takes account of processes of binding/unbinding of transcription factors, formation/dissolution of transcription apparatus, and modification of histone code at each locus of genes in the network. These processes are hierarchically interrelated to each other forming the dynamical feedback loops. By simulating stochastic dynamics of this model, we show that the phenotypic heterogeneity of ESCs can be explained when the chromatin at the Nanog locus undergoes the large scale reorganization in formation/dissolution of transcription apparatus, which should have the timescale similar to the cell cycle period. With this slow transcriptional switching of Nanog, the simulated ESCs fluctuate among multiple transient states, which can trigger the differentiation into the lineage-specific cell states. From the simulated transitions among cell states, the epigenetic landscape underlying transitions is calculated. The slow Nanog switching gives rise to the wide basin of ESC states in the landscape. The bimodal Nanog distribution arising from the kinetic flow running through this ESC basin prevents transdifferentiation and promotes the definite decision of the cell fate. These results show that the distribution of timescales of the regulatory processes is decisively important to characterize the fluctuation of cells and their differentiation process. The analyses through the epigenetic landscape and the kinetic flow on the landscape should provide a guideline to engineer cell differentiation.

Conformational flexibility of loops of myosin enhances the global bias in the actin-myosin interaction landscape.

A long-standing controversy on the mechanism of an actomyosin motor is the role of the Brownian motion of the myosin head in force generation. In order to shed light on this problem, we calculate free-energy landscapes of interaction between an actin filament and the head (S1) of myosin II by using a coarse-grained model of actomyosin. The results show that the free-energy landscape has a global gradient toward the strong-binding site on actin filament, which explains the biased Brownian motion of myosin S1 observed in a single-molecule experiment [Kitamura et al., Nature, 1999, 397, 129 and Biophysics, 2005, 1, 1]. The distinct global gradient in the landscape is brought about only when the conformation of loop 2 at the actin interface of myosin S1 is flexible. The conformational flexibility of loop 3 also contributes to the gradient in the landscape by compensating the role of loop 2. Though the structure of loop 2 is expanded in the weak-binding state, loop 2 shows the larger fluctuation of compaction and expansion due to the actin-myosin interactions as myosin S1 moves toward the strong-binding site on actin filament. Hence, the increase in the compaction-expansion fluctuation of loop 2, the stronger binding of myosin to actin, and the biased Brownian motion of myosin S1 are coupled with each other and should take place in a concurrent way. This predicted coupling should provide opportunities to further test the hypothesis of the biased Brownian motion in actomyosin.

Investigating the Transepithelial Transport and Enzymatic Stability of Lactononadecapeptide (NIPPLTQTPVVVPPFLQPE), a 19-Amino Acid Casein-Derived Peptide in Caco-2 Cells.

Lactononadecapeptide (LNDP; NIPPLTQTPVVVPPFLQPE), a casein-derived peptide comprising 19 residues, is known for its capacity to enhance cognitive function. This study aimed to explore the transepithelial transport and stability of LNDP. Results showed that LNDP retained over 90% stability after 2 h of treatment with gastrointestinal enzymes. The stability of LNDP on Caco-2 cell monolayers ranged from 93.4% ± 0.9% to 101.1% ± 1.2% over a period of 15-60 min, with no significant differences at each time point. The permeability of LNDP across an artificial lipid membrane was very low with the effective permeability of 3.6 × 10-11 cm/s. The Caco-2 assay demonstrated that LNDP could traverse the intestinal epithelium, with an apparent permeability of 1.22 × 10-6 cm/s. Its transport was significantly inhibited to 67.9% ± 5.0% of the control by Gly-Pro, a competitor of peptide transporter 1 (PEPT1). Furthermore, PEPT1 knockdown using siRNA significantly inhibited LNDP transport by 77.6% ± 1.9% in Caco-2 cell monolayers. The LNDP uptake in PEPT1-expressing HEK293 cells was significantly higher (54.5% ± 14.6%) than that in mock cells. These findings suggest that PEPT1 plays a crucial role in LNDP transport, and LNDP exhibits good resistance to gastrointestinal enzymes.

Entropic mechanism of large fluctuation in allosteric transition.

A statistical mechanical model of allosteric transitions in proteins is developed by extending the structure-based model of protein folding to cases of multiple native conformations. The partition function is calculated exactly within the model and the free-energy surface reflecting allostery is derived. This approach is applied to an example protein, the receiver domain of the bacterial enhancer-binding protein NtrC. The model predicts the large entropy associated with a combinatorial number of preexisting transition routes. This large entropy lowers the free-energy barrier of the allosteric transition, which explains the large structural fluctuation observed in the NMR data of NtrC. The global allosteric transformation of NtrC is explained by the shift of preexisting distribution of conformations upon phosphorylation, but the local structural adjustment around the phosphorylation site is explained by the complementary induced-fit mechanism. Structural disordering accompanied by fluctuating interactions specific to two allosteric conformations underlies a large number of routes of allosteric transition.

Unidirectional Brownian motion observed in an in silico single molecule experiment of an actomyosin motor.

The actomyosin molecular motor, the motor composed of myosin II and actin filament, is responsible for muscle contraction, converting chemical energy into mechanical work. Although recent single molecule and structural studies have shed new light on the energy-converting mechanism, the physical basis of the molecular-level mechanism remains unclear because of the experimental limitations. To provide a clue to resolve the controversy between the lever-arm mechanism and the Brownian ratchet-like mechanism, we here report an in silico single molecule experiment of an actomyosin motor. When we placed myosin on an actin filament and allowed myosin to move along the filament, we found that myosin exhibits a unidirectional Brownian motion along the filament. This unidirectionality was found to arise from the combination of a nonequilibrium condition realized by coupling to the ATP hydrolysis and a ratchet-like energy landscape inherent in the actin-myosin interaction along the filament, indicating that a Brownian ratchet-like mechanism contributes substantially to the energy conversion of this molecular motor.

Effects of Prebiotic Yeast Mannan on Gut Health and Sleep Quality in Healthy Adults: A Randomized, Double-Blind, Placebo-Controlled Study.

Human gut health is closely related to sleep. We aimed to evaluate the efficacy of yeast mannan (YM) in improving bowel habits and sleep quality, along with metabolomics in fecal samples. A total of 40 healthy adults (age range, 22-64 years) with discomfort in defecation were enrolled and randomly allocated to receive either YM (n = 20; 1.1 g/day) or placebo (n = 20) for four weeks. Participants recorded their defecation habits throughout the test periods. Sleep electroencephalogram (EEG) recording using an EEG device and fecal sampling were performed pre- and post-treatment. The YM group significantly increased defecation frequency and stool volumes compared to the placebo group. After 4 weeks of treatment, the non-REM sleep stage 3 (N3) duration in the YM group was significantly higher than that in the placebo group. YM ingestion significantly lengthened total time in bed (TIB) and significantly shortened N3 latency compared to placebo intake during the trial. The metabolomics analysis found a total of 20 metabolite differences between the YM and placebo groups. As a result of stepwise linear regression, changes in fecal propionate and gamma-aminobutyric acid (GABA) levels were identified as the primary factors explaining changes in TIB and N3 latency, respectively. Our findings suggest that the prebiotic YM could be beneficial to gut health and sleep quality.

Condensed but liquid-like domain organization of active chromatin regions in living human cells.

In eukaryotes, higher-order chromatin organization is spatiotemporally regulated as domains, for various cellular functions. However, their physical nature in living cells remains unclear (e.g., condensed domains or extended fiber loops; liquid-like or solid-like). Using novel approaches combining genomics, single-nucleosome imaging, and computational modeling, we investigated the physical organization and behavior of early DNA replicated regions in human cells, which correspond to Hi-C contact domains with active chromatin marks. Motion correlation analysis of two neighbor nucleosomes shows that nucleosomes form physically condensed domains with ~150-nm diameters, even in active chromatin regions. The mean-square displacement analysis between two neighbor nucleosomes demonstrates that nucleosomes behave like a liquid in the condensed domain on the ~150 nm/~0.5 s spatiotemporal scale, which facilitates chromatin accessibility. Beyond the micrometers/minutes scale, chromatin seems solid-like, which may contribute to maintaining genome integrity. Our study reveals the viscoelastic principle of the chromatin polymer; chromatin is locally dynamic and reactive but globally stable.

Dynamical modeling of three-dimensional genome organization in interphase budding yeast.

Eukaryotic genome is organized in a set of chromosomes each of which consists of a chain of DNA and associated proteins. Processes involving DNA such as transcription, duplication, and repair, therefore, should be intrinsically related to the three-dimensional organization of the genome. In this article, we develop a computational model of the three-dimensional organization of the haploid genome of interphase budding yeast by regarding chromosomes as chains moving under the constraints of nuclear structure and chromatin-chromatin interactions. The simulated genome structure largely fluctuates with the diffusive movement of chromosomes. This fluctuation, however, is not completely random, as parts of chromosomes distribute in characteristic ways to form "territories" in the nucleus. By suitably taking account of constraints arising from the data of the chromosome-conformation-capture measurement, the model explains the observed fluorescence data of chromosome distributions and motions.