Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Transverse single spin asymmetry in e+p↑→e+J/ψ+X and transverse momentum dependent evolution of the Sivers function

We extend our analysis of transverse single spin asymmetry in electroproduction of J/ψ to include the effect of the scale evolution of the transverse momentum dependent (TMD) parton distribution functions and gluon Sivers function. We estimate single spin asymmetry for JLab, HERMES, COMPASS, and eRHIC energies using the color evaporation model of charmonium production, using an analytically obtained approximate solution of TMD evolution equations discussed in the literature. We find that there is a reduction in the asymmetry compared with our predictions for the earlier case considered by us, wherein the Q2 dependence came only from DGLAP evolution of the unpolarized gluon densities and a different parametrization of the TMD Sivers function was used.

A FUNCTIONAL MODEL FOR PURE Gamma-CONTRACTIONS

A pair of commuting operators (S,P) defined on a Hilbert space H for which the closed symmetrized bidisc Gamma = {(z(1) + z(2), z(1)z(2)) : vertical bar z(1)vertical bar <= 1, vertical bar z(2)vertical bar <= 1} subset of C-2 is a spectral set is called a Gamma-contraction in the literature. A Gamma-contraction (S, P) is said to be pure if P is a pure contraction, i.e., P*(n) -> 0 strongly as n -> infinity Here we construct a functional model and produce a set of unitary invariants for a pure Gamma-contraction. The key ingredient in these constructions is an operator, which is the unique solution of the operator equation S - S*P = DpXDp, where X is an element of B(D-p), and is called the fundamental operator of the Gamma-contraction (S, P). We also discuss some important properties of the fundamental operator.

Laterally structured ripple and square phases with one and two dimensional thickness modulations in a model bilayer system

Molecular dynamics simulations of bilayers in a surfactant/co-surfactant/water system with explicit solvent molecules show formation of topologically distinct gel phases depending upon the bilayer composition. At low temperatures, the bilayers transform from the tilted gel phase, L beta', to the one dimensional (1D) rippled, P beta' phase as the surfactant concentration is increased. More interestingly, we observe a two dimensional (2D) square phase at higher surfactant concentration which, upon heating, transforms to the gel L beta' phase. The thickness modulations in the 1D rippled and square phases are asymmetric in two surfactant leaflets and the bilayer thickness varies by a factor of similar to 2 between maximum and minimum. The 1D ripple consists of a thinner interdigitated region of smaller extent alternating with a thicker non-interdigitated region. The 2D ripple phase is made up of two superimposed square lattices of maximum and minimum thicknesses with molecules of high tilt forming a square lattice translated from the lattice formed with the thickness minima. Using Voronoi diagrams we analyze the intricate interplay between the area-per-head-group, height modulations and chain tilt for the different ripple symmetries. Our simulations indicate that composition plays an important role in controlling the formation of low temperature gel phase symmetries and rippling accommodates the increased area-per-head-group of the surfactant molecules.

Hybrid and Rogue Kinases Encoded in the Genomes of Model Eukaryotes

The highly modular nature of protein kinases generates diverse functional roles mediated by evolutionary events such as domain recombination, insertion and deletion of domains. Usually domain architecture of a kinase is related to the subfamily to which the kinase catalytic domain belongs. However outlier kinases with unusual domain architectures serve in the expansion of the functional space of the protein kinase family. For example, Src kinases are made-up of SH2 and SH3 domains in addition to the kinase catalytic domain. A kinase which lacks these two domains but retains sequence characteristics within the kinase catalytic domain is an outlier that is likely to have modes of regulation different from classical src kinases. This study defines two types of outlier kinases: hybrids and rogues depending on the nature of domain recombination. Hybrid kinases are those where the catalytic kinase domain belongs to a kinase subfamily but the domain architecture is typical of another kinase subfamily. Rogue kinases are those with kinase catalytic domain characteristic of a kinase subfamily but the domain architecture is typical of neither that subfamily nor any other kinase subfamily. This report provides a consolidated set of such hybrid and rogue kinases gleaned from six eukaryotic genomes-S. cerevisiae, D. melanogaster, C. elegans, M. musculus, T. rubripes and H. sapiens-and discusses their functions. The presence of such kinases necessitates a revisiting of the classification scheme of the protein kinase family using full length sequences apart from classical classification using solely the sequences of kinase catalytic domains. The study of these kinases provides a good insight in engineering signalling pathways for a desired output. Lastly, identification of hybrids and rogues in pathogenic protozoa such as P. falciparum sheds light on possible strategies in host-pathogen interactions.

Combining data sets of satellite-retrieved products for basin-scale water balance study: 2. Evaluation on the Mississippi Basin and closure correction model

In this study, we applied the integration methodology developed in the companion paper by Aires (2014) by using real satellite observations over the Mississippi Basin. The methodology provides basin-scale estimates of the four water budget components (precipitation P, evapotranspiration E, water storage change Delta S, and runoff R) in a two-step process: the Simple Weighting (SW) integration and a Postprocessing Filtering (PF) that imposes the water budget closure. A comparison with in situ observations of P and E demonstrated that PF improved the estimation of both components. A Closure Correction Model (CCM) has been derived from the integrated product (SW+PF) that allows to correct each observation data set independently, unlike the SW+PF method which requires simultaneous estimates of the four components. The CCM allows to standardize the various data sets for each component and highly decrease the budget residual (P - E - Delta S - R). As a direct application, the CCM was combined with the water budget equation to reconstruct missing values in any component. Results of a Monte Carlo experiment with synthetic gaps demonstrated the good performances of the method, except for the runoff data that has a variability of the same order of magnitude as the budget residual. Similarly, we proposed a reconstruction of Delta S between 1990 and 2002 where no Gravity Recovery and Climate Experiment data are available. Unlike most of the studies dealing with the water budget closure at the basin scale, only satellite observations and in situ runoff measurements are used. Consequently, the integrated data sets are model independent and can be used for model calibration or validation.

Transverse single spin asymmetry in e plus p up arrow -> e plus J/psi plus X and Q(2) evolution of Sivers function-II

We present estimates of single spin asymmetry in the electroproduction of J/psi taking into account the transverse momentum-dependent (TMD) evolution of the gluon Sivers function. We estimate single spin asymmetry for JLab, HERMES, COMPASS and eRHIC energies using the color evaporation model of J/psi. We have calculated the asymmetry using recent parameters extracted by Echevarria et al. using the Collins-Soper-Sterman approach to TMD evolution. These recent TMD evolution fits are based on the evolution kernel in which the perturbative part is resummed up to next-to-leading logarithmic accuracy. We have also estimated the asymmetry by using parameters which had been obtained by a fit by Anselmino et al., using both an exact numerical and an approximate analytical solution of the TMD evolution equations. We find that the variation among the different estimates obtained using TMD evolution is much smaller than between these on one hand and the estimates obtained using DGLAP evolution on the other. Even though the use of TMD evolution causes an overall reduction in asymmetries compared to the ones obtained without it, they remain sizable. Overall, upon use of TMD evolution, predictions for asymmetries stabilize.

Colored polydimethylsiloxane micropillar arrays for high throughput measurements of forces applied by genetic model organisms

Measuring forces applied by multi-cellular organisms is valuable in investigating biomechanics of their locomotion. Several technologies have been developed to measure such forces, for example, strain gauges, micro-machined sensors, and calibrated cantilevers. We introduce an innovative combination of techniques as a high throughput screening tool to assess forces applied by multiple genetic model organisms. First, we fabricated colored Polydimethylsiloxane (PDMS) micropillars where the color enhances contrast making it easier to detect and track pillar displacement driven by the organism. Second, we developed a semiautomated graphical user interface to analyze the images for pillar displacement, thus reducing the analysis time for each animal to minutes. The addition of color reduced the Young's modulus of PDMS. Therefore, the dye-PDMS composite was characterized using Yeoh's hyperelastic model and the pillars were calibrated using a silicon based force sensor. We used our device to measure forces exerted by wild type and mutant Caenorhabditis elegans moving on an agarose surface. Wild type C. elegans exert an average force of similar to 1 mu N on an individual pillar and a total average force of similar to 7.68 mu N. We show that the middle of C. elegans exerts more force than its extremities. We find that C. elegans mutants with defective body wall muscles apply significantly lower force on individual pillars, while mutants defective in sensing externally applied mechanical forces still apply the same average force per pillar compared to wild type animals. Average forces applied per pillar are independent of the length, diameter, or cuticle stiffness of the animal. We also used the device to measure, for the first time, forces applied by Drosophila melanogaster larvae. Peristaltic waves occurred at 0.4Hz applying an average force of similar to 1.58 mu N on a single pillar. Our colored microfluidic device along with its displacement tracking software allows us to measure forces applied by multiple model organisms that crawl or slither to travel through their environment. (C) 2015 AIP Publishing LLC.

Computational design of model Re/Ru bearing Ni-base superalloys

It is well established that Re and Ru additions to Ni-base superalloys result in improved creep performance and phase stability. However, the role of Re and Ru and their synergetic effects are not well understood, and the first step in understanding these effects is to design alloys with controlled microstructural parameters. A computational approach was undertaken in the present work for designing model alloys with varying levels of Re and Ru. Thermodynamic and first principles calculations were employed complimentarily to design a set of alloys with varying Re and Ru levels, but which were constrained by constant microstructural parameters, i.e., phase fractions and lattice misfit across the alloys. Three ternary/quaternary alloys of type Ni-Al-xRe-yRu were thus designed. These compositions were subsequently cast, homogenized and aged. Experimental results suggest that while the measured volume fraction matches the predicted value in the Ru containing alloy, volume fraction is significantly higher than the designed value in the Re containing alloys. This is possibly due to errors in the thermodynamic database used to predict phase fraction and composition. These errors are also reflected in the mismatch between predicted and measured values of misfit.

Effects of Patterned Substrate on Thin Films Simulated by Family Model

Patterned substrate growth has been a subject of much interest. In this work, characteristics of some statistical properties of a film grown on triangular and vicinal substrates using the Family model are studied. Substrate size and tilt angle are varied. It is found that the interface width and the correlation function increase as the roughness of the pattern is increased. The new scaling exponents are calculated and anomalous scaling is obtained. The transient persistence probability does not show a power law relation when the initial surface is sufficiently rough. The initial rough surface also causes multifractal behavior in the model.

Dissecting Ubiquitin Folding Using the Self-Organized Polymer Model

Folding of Ubiquitin (Ub), a functionally important protein found in eukaryotic organisms, is investigated at low and neutral pH at different temperatures using simulations of the coarse-grained self-organized-polymer model with side chains (SOP-SC). The melting temperatures (T-m's), identified with the peaks in the heat capacity curves, decrease as pH decreases, in qualitative agreement with experiments. The calculated radius of gyration, showing dramatic variations with pH, is in excellent agreement with scattering experiments. At T-m Ub folds in a two-state manner at low and neutral pH. Clustering analysis of the conformations sampled in equilibrium folding trajectories at T-m with multiple transitions between the folded and unfolded states, shows a network of metastable states connecting the native and unfolded states. At low and neutral pH, Ub folds with high probability through a preferred set of conformations resulting in a pH-dependent dominant folding pathway. Folding kinetics reveal that Ub assembly at low pH occurs by multiple pathways involving a combination of nucleation-collapse and diffusion collision mechanism. The mechanism by which Ub folds is dictated by the stability of the key secondary structural elements responsible for establishing long-range contacts and collapse of Ub. Nucleation collapse mechanism holds if the stability of these elements are marginal, as would be the case at elevated temperatures. If the lifetimes associated with these structured microdomains are on the order of hundreds of microseconds, then Ub folding follows the diffusion collision mechanism with intermediates, many of which coincide with those found in equilibrium. Folding at neutral pH is a sequential process with a populated intermediate resembling that sampled at equilibrium. The transition state structures, obtained using a P-fold analysis, are homogeneous and globular with most of the secondary and tertiary structures being native-like. Many of our findings for both the thermodynamics and kinetics of folding are not only in agreement with experiments but also provide missing details not resolvable in standard experiments. The key prediction that folding mechanism varies dramatically with pH is amenable to experimental tests.

Controlled release kinetics of p-aminosalicylic acid from biodegradable crosslinked polyesters for enhanced anti-mycobacterial activity

Unlike conventional polymeric drug delivery systems, where drugs are entrapped in polymers, this study focuses on the incorporation of the drug into the polymer backbone to achieve higher loading and sustained release. Crosslinked, biodegradable, xylitol based polyesters have been synthesized in this study. The bioactive drug moiety, p-aminosalicylic acid (PAS), was incorporated in xylitol based polyesters to impart its anti-mycobacterial activity. To understand the influence of the monomer chemistry on the incorporation of PAS and its subsequent release from the polymer, different diacids have been used. Controlled release profiles of the drug from these polyesters were studied under normal physiological conditions. The degradation of the polyesters varied from 48% to 76% and the release of PAS ranged from 54% to 65% of its initial loading in 7 days. A new model was developed to explain the release kinetics of PAS from the polymer that accounted for the polymer degradation and drug concentration. The thermal, mechanical, drug release and cytocompatibility properties of the polymers indicate their suitability in biomedical applications. The released products from these polymers were observed to be pharmacologically active against Mycobacteria. The high drug loading and sustained release also ensured enhanced efficacy. These polymers form biocompatible, biodegradable polyesters where the sustained release of PAS may be tailored for potential treatment of mycobacterial infections. Statement of significance In the present work, we report on novel polyesters with p-aminosalicylic acid (PAS) incorporated in the polymer backbone. The current work aims to achieve controlled release of PAS and ensures the delivered PAS is stable and pharmacologically active. The novelty of this work primarily involves the synthetic chemistry of polymerization and detailed analysis and efficacy of active PAS delivery. A new kinetic model has been developed to explain the PAS release profiles. These polymers are biodegradable, cytocompatible and anti-mycobacterial in nature. (C) 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Differential proteomic and genomic profiling of mouse striatal cell model of Huntington's disease and control; probable implications to the disease biology

Huntington's disease (HD) is an autosomal dominant disorder of central nervous system caused by expansion of CAG repeats in exon1 of the huntingtin gene (Htt). Among various dysfunctions originated from the mutation in Htt gene, transcriptional deregulation has been considered to be one of the most important abnormalities. Large numbers of investigations identified altered expressions of genes in brains of HD patients and many models of HD. In this study we employed 2D SDS-PAGE/MALDI-MS coupled with 2D-DIGE and real-time PCR experiments of an array of genes focused to HD pathway to determine altered protein and gene expressions in STHdh(Q111)/Hdh(Q111) cells, a cell model of HD and compared with STHdh(Q7)/Hdh(Q7) cells, its wild type counterpart. We annotated 76 proteins from these cells and observed differential expressions of 31 proteins (by 2D-DIGE) involved in processes like unfolded protein binding, negative regulation of neuron apoptosis, response to superoxides etc. Our PCR array experiments identified altered expressions of 47 genes. Altogether significant alteration of 77 genes/proteins could be identified in this HD cell line with potential relevance to HD biology. Biological significance: In this study we intended to find out differential proteomic and genomic profiles in HD condition. We used the STHdh cells, a cellular model for HD and control. These are mouse striatal neuronal cell lines harboring 7 and 111 knock -in CAG repeats in their two alleles. The 111Q containing cell line (STHdh(Q111)/Hdh(Q111)) mimics diseased condition, whereas the 7Q containing ones (STHdh(Q7)/Hdh(Q7)), serves as the proper control cell line. Proteomic experiments were performed earlier to obtain differential expressions of proteins in R6/2 mice models, Hdh(Q) knock -in mice and in plasma and CSF from HD patients. However, no earlier report on proteomic alterations in these two HD cell lines and control was available in literature. It was, therefore, an important objective to find out differential expressions of proteins in these two cell lines. In this study, we annotated 76 proteins from STHdh(Q7)/Hdh(Q7) and STHdh(Q111)/Hdh(Q111) cells using 2D-gel/mass spectrometry. Next, by performing 2D-DIGE, we observed differential expressions of 31 proteins (16 upregulated and 15 downregulated) between these two cell lines. We also performed customized qRT-PCR array focused to HD pathway and found differential expressions of 47 genes (8 gene exptessions increased and 39 genes were decreased significantly). A total of 77 genes/proteins (Htt downregulated in both the studies) were found to be significantly altered from both the experimental paradigms. We validated the differential expressions of Vim, Hypk, Ran, Dstn, Hspa5 and Sod2 either by qRT-PCR or Western blot analysis or both. Out of these 77, similar trends in alteration of 19 out of 31 and 38 out of 47 proteins/genes were reported in earlier studies. Thus our study confirmed earlier observations on differential gene/protein expressions in HD and are really useful. Additionally, we observed differential expression of some novel genes/proteins. One of this was Hypk, a Htt-interacting chaperone protein with the ability to solubilize mHtt aggregated structures in cell lines. We propose that downregulation of Hypk in STHdh-Qm (Q111)/Hdh(Q111) has a causal effect towards HD pathogenesis. Thus the novel findings from our study need further research and might be helpful to understand the molecular mechanism behind HD pathogenesis. (C) 2015 Elsevier B.V. All rights reserved.