Probability distributions of the work in the 2D-Ising model

Probability distributions of the magnetic work are computed for the 2D Ising model by means of Monte Carlo simulations. The system is first prepared at equilibrium for three temperatures below, at and above the critical point. A magnetic field is then applied and grown linearly at different rates. Probability distributions of the work are stored and free energy differences computed using the Jarzynski equality. Consistency is checked and the dynamics of the system is analyzed. Free energies and dissipated works are reproduced with simple models. The critical exponent $\delta$ is estimated in an usual manner.Comment: 12 pages, 6 figures. Comments are welcom

Parallel Computing for Probabilistic Response Analysis of High Temperature Composites

The objective of this Phase I research was to establish the required software and hardware strategies to achieve large scale parallelism in solving PCM problems. To meet this objective, several investigations were conducted. First, we identified the multiple levels of parallelism in PCM and the computational strategies to exploit these parallelisms. Next, several software and hardware efficiency investigations were conducted. These involved the use of three different parallel programming paradigms and solution of two example problems on both a shared-memory multiprocessor and a distributed-memory network of workstations

Myostatin signals through miR-34a to regulate Fndc5 expression and browning of white adipocytes

Background/Objectives: Myostatin (Mstn) has a pivotal role in glucose and lipid metabolism. Mstn deficiency leads to the increased browning of white adipose tissue (WAT), which results in the increased energy expenditure and protection against diet-induced obesity and insulin resistance. In this study, we investigated the molecular mechanism(s) through which Mstn regulates browning of white adipocytes. Methods: Quantitative molecular analyses were performed to assess Mstn regulation of miR-34a and Fndc5 expression. miR-34a was overexpressed and repressed to investigate miR-34a regulation of Fndc5. Luciferase reporter analysis verified direct binding between miR-34a and the Fndc5 3′-untranslated region (UTR). The browning phenotype of Mstn−/− adipocytes was assessed through the analysis of brown fat marker gene expression, mitochondrial function and infrared thermography. The role of miR-34a and Fndc5 in this browning phenotype was verified through antibody-mediated neutralization of FNDC5, knockdown of Fndc5 by small interfering RNA and through miR-34a gain-of-function and loss-of-function experiments. Results: Mstn treatment of myoblasts inhibited Fndc5 expression, whereas the loss of Mstn increased Fndc5 levels in muscles and in circulation. Mstn inhibition of Fndc5 is miR-34a dependent. Mstn treatment of C2C12 myoblasts upregulated miR-34a expression, whereas reduced miR-34a expression was noted in Mstn−/− muscle and WAT. Subsequent overexpression of miR-34a inhibited Fndc5 expression, whereas blockade of miR-34a increased Fndc5 expression in myoblasts. Reporter analysis revealed that miR-34a directly suppresses Fndc5 expression through a miR-34a-specific binding site within the Fndc5 3′UTR. Importantly, Mstn-mediated inhibition of Fndc5 was blocked upon miR-34a inhibition. Mstn−/− adipocytes showed reduced miR-34a, enhanced Fndc5 expression and increased thermogenic gene expression, which was reversed upon either neutralization of Fndc5 or Fndc5 knockdown. In agreement, Mstn−/− adipocytes have increased mitochondria, improved mitochondrial function and increased heat production. Conclusions: Mstn regulates Fndc5/Irisin expression and secretion through a novel miR-34a-dependent post-transcriptional mechanism. Loss of Mstn in mice leads to the increased Fndc5/Irisin expression, which contributes to the browning of white adipocytes

Accounting for epistatic interactions improves the functional analysis of protein structures

Motivation: The constraints under which sequence, structure and function coevolve are not fully understood. Bringing this mutual relationship to light can reveal the molecular basis of binding, catalysis and allostery, thereby identifying function and rationally guiding protein redesign. Underlying these relationships are the epistatic interactions that occur when the consequences of a mutation to a protein are determined by the genetic background in which it occurs. Based on prior data, we hypothesize that epistatic forces operate most strongly between residues nearby in the structure, resulting in smooth evolutionary importance across the structure. Methods and Results: We find that when residue scores of evolutionary importance are distributed smoothly between nearby residues, functional site prediction accuracy improves. Accordingly, we designed a novel measure of evolutionary importance that focuses on the interaction between pairs of structurally neighboring residues. This measure that we term pair-interaction Evolutionary Trace yields greater functional site overlap and better structure-based proteome-wide functional predictions. Conclusions: Our data show that the structural smoothness of evolutionary importance is a fundamental feature of the coevolution of sequence, structure and function. Mutations operate on individual residues, but selective pressure depends in part on the extent to which a mutation perturbs interactions with neighboring residues. In practice, this principle led us to redefine the importance of a residue in terms of the importance of its epistatic interactions with neighbors, yielding better annotation of functional residues, motivating experimental validation of a novel functional site in LexA and refining protein function prediction. Contact: [email protected] Supplementary information: Supplementary data are available at Bioinformatics online

A Dual Spring Modeling Approach for Static and Fatigue Failure Assessments of Carbon/Epoxy Composite Sub-Elements

A dual spring model is developed for the static and fatigue damage predictions of laminates interface in composite structures. Stress concentrations can be induced by the defects formed in the fabrication or service process. A conventional S-N based fatigue damage model may not be accurate to predict the fatigue life of a structure with high stress concentration. With the dual spring model, static delamination failure can be simulated using springs of cohesive type material model while fatigue delamination development can be predicted using linear springs, where the crack driving force is computed based on virtual crack closure technique (VCCT). A Paris law type fatigue growth law with its mode mixity is applied for fatigue crack growth prediction. After verified using benchmark examples, including Double Cantilever Beam (DCB), End-Notched Flexure (ENF) and Mix-Mode Bending (MMB), the proposed dual spring model is applied in the static and fatigue damage prediction of NASA/Boeing sub-elements and UTC sub-elements

GWmodelS: a standalone software to train geographically weighted models

With the recent increase in studies on spatial heterogeneity, geographically weighted (GW) models have become an essential set of local techniques, attracting a wide range of users from different domains. In this study, we demonstrate a newly developed standalone GW software, GWmodelS using a community-level house price data set for Wuhan, China. In detail, a number of fundamental GW models are illustrated, including GW descriptive statistics, basic and multiscale GW regression, and GW principle component analysis. Additionally, functionality in spatial data management and batch mapping are presented as essential supplementary activities for GW modeling. The software provides significant advantages in terms of a user-friendly graphical user interface, operational efficiency, and accessibility, which facilitate its usage for users from a wide range of domains

Abundance of unknots in various models of polymer loops

A veritable zoo of different knots is seen in the ensemble of looped polymer chains, whether created computationally or observed in vitro. At short loop lengths, the spectrum of knots is dominated by the trivial knot (unknot). The fractional abundance of this topological state in the ensemble of all conformations of the loop of $N$ segments follows a decaying exponential form, $\sim \exp (-N/N_0)$, where $N_0$ marks the crossover from a mostly unknotted (ie topologically simple) to a mostly knotted (ie topologically complex) ensemble. In the present work we use computational simulation to look closer into the variation of $N_0$ for a variety of polymer models. Among models examined, $N_0$ is smallest (about 240) for the model with all segments of the same length, it is somewhat larger (305) for Gaussian distributed segments, and can be very large (up to many thousands) when the segment length distribution has a fat power law tail.Comment: 13 pages, 6 color figure

Remote exercise testing in pulmonary hypertension (PHRET)

Remote exercise tests for patients with pulmonary hypertension (PH) would improve the telemedicine strategies in this disease. The PHRET study assessed the validity and feasibility of four remote exercise tests performed by PH patients at home. Participants undergoing diagnostic assessment for PH were included. At baseline, patients completed a 6MWT followed by a range of study tests including a Timed Up and Go (TUG) test, a Sit-to-Stand (STS), a Step Test (ST), and a tele-6MWT (T6MWT) performed outside using a GPS-enabled smartphone. Patients performed these tests at home following discharge and at first follow-up. Analysis focused on comparing the results of study tests to the standard 6MWT. The discontinuation rate was 15%. Ninety-seven percent of patients were able to complete a TUG, 92% a STS, 73% a ST, and 49% a T6MWT. At baseline, correlation between the standard 6MWT and study tests, respectively, was T6MWT 0.93, ST 0.78, STS 0.71, and TUG −0.76 (p < 0.001). Direction of change in the study test agreed with the standard 6MWT in 68% of the follow-up ST, 68% of the STS, 71% of the TUG, and 79% of the T6MWT. Patients were able to complete the tests at home, there were no adverse incidents and ≥92% of patients were happy to continue performing home tests. Remote exercise testing is feasible. The T6MWT was a valid remote measure of exercise capacity, but could only be performed by a limited number of patients. The high discontinuation rate may impact the utility of remote tests

The Energy Landscape, Folding Pathways and the Kinetics of a Knotted Protein

The folding pathway and rate coefficients of the folding of a knotted protein are calculated for a potential energy function with minimal energetic frustration. A kinetic transition network is constructed using the discrete path sampling approach, and the resulting potential energy surface is visualized by constructing disconnectivity graphs. Owing to topological constraints, the low-lying portion of the landscape consists of three distinct regions, corresponding to the native knotted state and to configurations where either the N- or C-terminus is not yet folded into the knot. The fastest folding pathways from denatured states exhibit early formation of the N-terminus portion of the knot and a rate-determining step where the C-terminus is incorporated. The low-lying minima with the N-terminus knotted and the C-terminus free therefore constitute an off-pathway intermediate for this model. The insertion of both the N- and C-termini into the knot occur late in the folding process, creating large energy barriers that are the rate limiting steps in the folding process. When compared to other protein folding proteins of a similar length, this system folds over six orders of magnitude more slowly.Comment: 19 page

Residence Time Statistics for Normal and Fractional Diffusion in a Force Field

We investigate statistics of occupation times for an over-damped Brownian particle in an external force field. A backward Fokker-Planck equation introduced by Majumdar and Comtet describing the distribution of occupation times is solved. The solution gives a general relation between occupation time statistics and probability currents which are found from solutions of the corresponding problem of first passage time. This general relationship between occupation times and first passage times, is valid for normal Markovian diffusion and for non-Markovian sub-diffusion, the latter modeled using the fractional Fokker-Planck equation. For binding potential fields we find in the long time limit ergodic behavior for normal diffusion, while for the fractional framework weak ergodicity breaking is found, in agreement with previous results of Bel and Barkai on the continuous time random walk on a lattice. For non-binding potential rich physical behaviors are obtained, and classification of occupation time statistics is made possible according to whether or not the underlying random walk is recurrent and the averaged first return time to the origin is finite. Our work establishes a link between fractional calculus and ergodicity breaking.Comment: 12 page