Graph's Topology and Free Energy of a Spin Model on the Graph.
Ontology highlight
ABSTRACT:
SUBMITTER: Choi JM
PROVIDER: S-EPMC5668130 | biostudies-literature | 2017 Feb
REPOSITORIES: biostudies-literature
Publications
Graph's Topology and Free Energy of a Spin Model on the Graph.
Choi Jeong-Mo JM Gilson Amy I AI Shakhnovich Eugene I EI
Physical review letters 20170224 8
In this Letter we investigate a direct relationship between a graph's topology and the free energy of a spin system on the graph. We develop a method of separating topological and energetic contributions to the free energy, and find that considering the topology is sufficient to qualitatively compare the free energies of different graph systems at high temperature, even when the energetics are not fully known. This method was applied to the metal lattice system with defects, and we found that it ...[more]
Similar Datasets
Project description:Free-energy calculations play an important role in the application of computational chemistry to a range of fields, including protein biochemistry, rational drug design, or materials science. Importantly, the free-energy difference is directly related to experimentally measurable quantities such as partition and adsorption coefficients, water activity, and binding affinities. Among several techniques aimed at predicting free-energy differences, perturbation approaches, involving the alchemical transformation of one molecule into another through intermediate states, stand out as rigorous methods based on statistical mechanics. However, despite the importance of free-energy calculations, the applicability of the perturbation approaches is still largely impeded by a number of challenges, including the definition of the perturbation path, i.e., alchemical changes leading to the transformation of one molecule to the other. To address this, an automatic perturbation topology builder based on a graph-matching algorithm is developed, which can identify the maximum common substructure (MCS) of two or multiple molecules and provide the perturbation topologies suitable for free-energy calculations using the GROMOS and the GROMACS simulation packages. Various MCS search options are presented leading to alternative definitions of the perturbation pathway. Moreover, perturbation topologies generated using the default multistate MCS search are used to calculate the changes in free energy between lysine and its two post-translational modifications, 3-methyllysine and acetyllysine. The pairwise free-energy calculations performed on this test system led to a cycle closure of 0.5 ± 0.3 and 0.2 ± 0.2 kJ mol-1, with GROMOS and GROMACS simulation packages, respectively. The same relative free energies between the three states are obtained by employing the enveloping distribution sampling (EDS) approach when compared to the pairwise perturbations. Importantly, this toolkit is made available online as an open-source Python package (https://github.com/drazen-petrov/SMArt).
| S-EPMC8479811 | biostudies-literature
Project description:Despite its role in spin density functional theory and it being the basic observable for describing and understanding magnetic phenomena, few studies have appeared on the electron spin density subtleties thus far. A systematic full topological analysis of this function is lacking, seemingly in contrast to the blossoming in the last 20 years of many studies on the topological features of other scalar fields of chemical interest. We aim to fill this gap by unveiling the kind of information hidden in the spin density distribution that only its topology can disclose. The significance of the spin density critical points, the 18 different ways in which they can be realized and the peculiar topological constraints on their number and kind, arising from the presence of positive and negative spin density regions, is addressed. The notion of molecular spin graphs, spin maxima (minima) joining paths, spin basins and of their valence is introduced. We show that two kinds of structures are associated with a spin-polarized molecule: the usual one, defined through the electron density gradient, and the magnetic structure, defined through the spin density gradient and composed in general by at least two independent spin graphs, related to spin density maxima and minima. Several descriptors, such as the spin polarization index, are introduced to characterize the properties of spin density critical points and basins. The study on the general features of the spin density topology is followed by the specific example of the water molecule in the 3B1 triplet state, using spin density distributions of increasing accuracy.
| S-EPMC7436107 | biostudies-literature
Project description:MotivationSince the amount of published biological sequencing data is growing exponentially, efficient methods for storing and indexing this data are more needed than ever to truly benefit from this invaluable resource for biomedical research. Labeled de Bruijn graphs are a frequently-used approach for representing large sets of sequencing data. While significant progress has been made to succinctly represent the graph itself, efficient methods for storing labels on such graphs are still rapidly evolving.ResultsIn this article, we present RowDiff, a new technique for compacting graph labels by leveraging expected similarities in annotations of vertices adjacent in the graph. RowDiff can be constructed in linear time relative to the number of vertices and labels in the graph, and in space proportional to the graph size. In addition, construction can be efficiently parallelized and distributed, making the technique applicable to graphs with trillions of nodes. RowDiff can be viewed as an intermediary sparsification step of the original annotation matrix and can thus naturally be combined with existing generic schemes for compressed binary matrices. Experiments on 10 000 RNA-seq datasets show that RowDiff combined with multi-BRWT results in a 30% reduction in annotation footprint over Mantis-MST, the previously known most compact annotation representation. Experiments on the sparser Fungi subset of the RefSeq collection show that applying RowDiff sparsification reduces the size of individual annotation columns stored as compressed bit vectors by an average factor of 42. When combining RowDiff with a multi-BRWT representation, the resulting annotation is 26 times smaller than Mantis-MST.Availability and implementationRowDiff is implemented in C++ within the MetaGraph framework. The source code and the data used in the experiments are publicly available at https://github.com/ratschlab/row_diff.
| S-EPMC8346655 | biostudies-literature
Project description:Identifying the conditions that support cooperation in spatial evolutionary game theory has been the focus of a large body of work. In this paper, the classical Prisoner's Dilemma is adopted as an interaction model; agents are placed on graphs and their interactions are constrained by a graph topology. A simple strategy update mechanism is used where agents copy the best performing strategy of their neighbourhood (including themselves). In this paper, we begin with a fully cooperative population and explore the robustness of the population to the introduction of defectors. We introduce a graph structure that has the property that the initial fully cooperative population is robust to any one perturbation (a change of any cooperator to a defector). We present a proof of this property and specify the necessary constraints on the graph. Furthermore, given the standard game payoffs, we calculate the smallest graph which possesses this property. We present an approach for increasing the size of the graph and we show empirically that this extended graph is robust to an increasing percentage of perturbations. We define a new class of graphs for the purpose of future work.
| S-EPMC8097208 | biostudies-literature
Project description:Knowledge of the transfer free energy of amino acids from aqueous solution to a lipid bilayer is essential for understanding membrane protein folding and for predicting membrane protein structure. Here we report a computational approach that can calculate the folding free energy of the transmembrane region of outer membrane ?-barrel proteins (OMPs) by combining an empirical energy function with a reduced discrete state space model. We quantitatively analyzed the transfer free energies of 20 amino acid residues at the center of the lipid bilayer of OmpLA. Our results are in excellent agreement with the experimentally derived hydrophobicity scales. We further exhaustively calculated the transfer free energies of 20 amino acids at all positions in the TM region of OmpLA. We found that the asymmetry of the Gram-negative bacterial outer membrane as well as the TM residues of an OMP determine its functional fold in vivo. Our results suggest that the folding process of an OMP is driven by the lipid-facing residues in its hydrophobic core, and its NC-IN topology is determined by the differential stabilities of OMPs in the asymmetrical outer membrane. The folding free energy is further reduced by lipid A and assisted by general depth-dependent cooperativities that exist between polar and ionizable residues. Moreover, context-dependency of transfer free energies at specific positions in OmpLA predict regions important for protein function as well as structural anomalies. Our computational approach is fast, efficient and applicable to any OMP.
| S-EPMC4881743 | biostudies-literature
Project description:Neuroimaging studies of functional connectivity using graph theory have furthered our understanding of the network structure in temporal lobe epilepsy (TLE). Brain network effects of anti-epileptic drugs could influence such studies, but have not been systematically studied. Resting-state functional MRI was analyzed in 25 patients with TLE using graph theory analysis. Patients were divided into two groups based on anti-epileptic medication use: those taking carbamazepine/oxcarbazepine (CBZ/OXC) (n=9) and those not taking CBZ/OXC (n=16) as a part of their medication regimen. The following graph topology metrics were analyzed: global efficiency, betweenness centrality (BC), clustering coefficient, and small-world index. Multiple linear regression was used to examine the association of CBZ/OXC with graph topology. The two groups did not differ from each other based on epilepsy characteristics. Use of CBZ/OXC was associated with a lower BC. Longer epilepsy duration was also associated with a lower BC. These findings can inform graph theory-based studies in patients with TLE. The changes observed are discussed in relation to the anti-epileptic mechanism of action and adverse effects of CBZ/OXC.
| S-EPMC4490704 | biostudies-literature
Project description:A gene regulatory network can be described at a high level by a directed graph with signed edges, and at a more detailed level by a system of ordinary differential equations (ODEs). The former qualitatively models the causal regulatory interactions between ordered pairs of genes, while the latter quantitatively models the time-varying concentrations of mRNA and proteins. This paper clarifies the connection between the two types of models. We propose a property, called the constant sign property, for a general class of ODE models. The constant sign property characterizes the set of conditions (system parameters, external signals, or internal states) under which an ODE model is consistent with a signed, directed graph. If the constant sign property for an ODE model holds globally for all conditions, then the ODE model has a single signed, directed graph. If the constant sign property for an ODE model only holds locally, which may be more typical, then the ODE model corresponds to different graphs under different sets of conditions. In addition, two versions of constant sign property are given and a relationship between them is proved. As an example, the ODE models that capture the effect of cis-regulatory elements involving protein complex binding, based on the model in the GeneNetWeaver source code, are described in detail and shown to satisfy the global constant sign property with a unique consistent gene regulatory graph. Even a single gene regulatory graph is shown to have many ODE models of GeneNetWeaver type consistent with it due to combinatorial complexity and continuous parameters. Finally the question of how closely data generated by one ODE model can be fit by another ODE model is explored. It is observed that the fit is better if the two models come from the same graph.
| S-EPMC7326199 | biostudies-literature
Project description:Excitations in magnetic structures of the so-called spin-ice materials generate two different peaks in the specific heat and anomalies in entropy in the temperature interval between 0 and 1?K. These points are due to the existence of two low-energy excited global states which seem to transit from a bosonic condensate towards a magnetic neutral plasma in a narrow temperature interval between 0.05???T???1?K. In this paper, we determine the characteristic features of two states and we analyze the possibilities of existence of a BEC state and its phase transition to the magnetic plasma state from a model of two magnetic charge fluids. From the structural analysis of the many-body excitation states, we obtain theoretical results about entropy and specific heat since these two key physical magnitudes announce the phase transitions. We give criteria for distinguishing if some of these phase transitions is of either first or second order.
| S-EPMC6215019 | biostudies-literature
Project description:This study details the electronic structure of the heme–peroxo–copper adduct {[(F8)Fe(DCHIm)]-O2-[Cu(AN)]}+ (LS(AN)) in which O2(2–) bridges the metals in a ?-1,2 or “end-on” configuration. LS(AN) is generated by addition of coordinating base to the parent complex {[(F8)Fe]-O2-[Cu(AN)]}+ (HS(AN)) in which the O2(2–) bridges the metals in an ?-?2:?2 or “side-on” mode. In addition to the structural change of the O2(2–) bridging geometry, coordination of the base changes the spin state of the heme fragment (from S = 5/2 in HS(AN) to S = 1/2 in LS(AN)) that results in an antiferromagnetically coupled diamagnetic ground state in LS(AN). The strong ligand field of the porphyrin modulates the high-spin to low-spin effect on Fe–peroxo bonding relative to nonheme complexes, which is important in the O–O bond cleavage process. On the basis of DFT calculations, the ground state of LS(AN) is dependent on the Fe–O–O–Cu dihedral angle, wherein acute angles (<~150°) yield an antiferromagnetically coupled electronic structure while more obtuse angles yield a ferromagnetic ground state. LS(AN) is diamagnetic and thus has an antiferromagnetically coupled ground state with a calculated Fe–O–O–Cu dihedral angle of 137°. The nature of the bonding in LS(AN) and the frontier molecular orbitals which lead to this magneto-structural correlation provide insight into possible spin topology contributions to O–O bond cleavage by cytochrome c oxidase.
| S-EPMC3226806 | biostudies-literature
Project description:While ubiquitous, energy redistribution remains a poorly understood facet of the nonequilibrium thermodynamics of biomolecules. At the molecular level, finite-size effects, pronounced nonlinearities, and ballistic processes produce behavior that diverges from the macroscale. Here, we show that transient thermal transport reflects macromolecular energy landscape architecture through the topological characteristics of molecular contacts and the nonlinear processes that mediate dynamics. While the former determines transport pathways via pairwise interactions, the latter reflects frustration within the landscape for local conformational rearrangements. Unlike transport through small-molecule systems, such as alkanes, nonlinearity dominates over coherent processes at even quite short time- and length-scales. Our exhaustive all-atom simulations and novel local-in-time and space analysis, applicable to both theory and experiment, permit dissection of energy migration in biomolecules. The approach demonstrates that vibrational energy transport can probe otherwise inaccessible aspects of macromolecular dynamics and interactions that underly biological function.
| S-EPMC6789131 | biostudies-literature