Predicting gene ontology annotations of orphan GWAS genes using protein-protein interactions

Background: The number of genome-wide association studies (GWAS) has increased rapidly in the past couple of years, resulting in the identification of genes associated with different diseases. The next step in translating these findings into biomedically useful information is to find out the mechanism of the action of these genes. However, GWAS studies often implicate genes whose functions are currently unknown; for example, MYEOV, ANKLE1, TMEM45B and ORAOV1 are found to be associated with breast cancer, but their molecular function is unknown. Results: We carried out Bayesian inference of Gene Ontology (GO) term annotations of genes by employing the directed acyclic graph structure of GO and the network of protein-protein interactions (PPIs). The approach is designed based on the fact that two proteins that interact biophysically would be in physical proximity of each other, would possess complementary molecular function, and play role in related biological processes. Predicted GO terms were ranked according to their relative association scores and the approach was evaluated quantitatively by plotting the precision versus recall values and F-scores (the harmonic mean of precision and recall) versus varying thresholds. Precisions of similar to 58% and similar to 40% for localization and functions respectively of proteins were determined at a threshold of similar to 30 (top 30 GO terms in the ranked list). Comparison with function prediction based on semantic similarity among nodes in an ontology and incorporation of those similarities in a k nearest neighbor classifier confirmed that our results compared favorably. Conclusions: This approach was applied to predict the cellular component and molecular function GO terms of all human proteins that have interacting partners possessing at least one known GO annotation. The list of predictions is available at http://severus.dbmi.pitt.edu/engo/GOPRED.html. We present the algorithm, evaluations and the results of the computational predictions, especially for genes identified in GWAS studies to be associated with diseases, which are of translational interest.

Identification of heat shock factor binding protein in Plasmodium falciparum

Background: Heat shock factor binding protein (HSBP) was originally discovered in a yeast two-hybrid screen as an interacting partner of heat shock factor (HSF). It appears to be conserved in all eukaryotes studied so far, with yeast being the only exception. Cell biological analysis of HSBP in mammals suggests its role as a negative regulator of heat shock response as it appears to interact with HSF only during the recovery phase following exposure to heat stress. While the identification of HSF in the malaria parasite is still eluding biologists, this study for the first time, reports the presence of a homologue of HSBP in Plasmodium falciparum. Methods: PfHSBP was cloned and purified as his-tag fusion protein. CD (Circular dichroism) spectroscopy was performed to predict the secondary structure. Immunoblots and immunofluorescence approaches were used to study expression and localization of HSBP in P. falciparum. Cellular fractionation was performed to examine subcellular distribution of PfHSBP. Immunoprecipitation was carried out to identify HSBP interacting partner in P. falciparum. Results: PfHSBP is a conserved protein with a high helical content and has a propensity to form homo-oligomers. PfHSBP was cloned, expressed and purified. The in vivo protein expression profile shows maximal expression in trophozoites. The protein was found to exist in oligomeric form as trimer and hexamer. PfHSBP is predominantly localized in the parasite cytosol, however, upon heat shock, it translocates to the nucleus. This study also reports the interaction of PfHSBP with PfHSP70-1 in the cytoplasm of the parasite. Conclusions: This study emphasizes the structural and biochemical conservation of PfHSBP with its mammalian counterpart and highlights its potential role in regulation of heat shock response in the malaria parasite. Analysis of HSBP may be an important step towards identification of the transcription factor regulating the heat shock response in P. falciparum.

Titanium promoted reduction of imines with Grignards, silanes, and zinc: identification of a new mechanism with silanes

Aldimines react with reducing agents, such as Grignards, phenylsilane or zinc in the presence of titanium(IV) isopropoxide to form amines and reductively coupled imines (diamines). Using deuterium labeled reagents, the mechanism of reduction to form amines is described. Reducing agents, such as the Grignard and zinc result in the formation of low valent titanium (LVT), which in turn reduces the imine. On the other hand, phenylsilane reacts by a distinctly different mechanism and where a hydrogen atom from silicon is directly transferred to the titanium coordinated imine. (c) 2014 Elsevier Ltd. All rights reserved.

Residue specific contributions to stability and activity inferred from saturation mutagenesis and deep sequencing

Multiple methods currently exist for rapid construction and screening of single-site saturation mutagenesis (SSM) libraries in which every codon or nucleotide in a DNA fragment is individually randomized. Nucleotide sequences of each library member before and after screening or selection can be obtained through deep sequencing. The relative enrichment of each mutant at each position provides information on its contribution to protein activity or ligand-binding under the conditions of the screen. Such saturation scans have been applied to diverse proteins to delineate hot-spot residues, stability determinants, and for comprehensive fitness estimates. The data have been used to design proteins with enhanced stability, activity and altered specificity relative to wild-type, to test computational predictions of binding affinity, and for protein model discrimination. Future improvements in deep sequencing read lengths and accuracy should allow comprehensive studies of epistatic effects, of combinational variation at multiple sites, and identification of spatially proximate residues.

Identification of miRNAs from French bean (Phaseolus vulgaris) under low nitrate stress

Objective: In this study, we report the role of miRNAs involved under nitrogen starvation from widely grown vegetable crop, French bean. In recent years, a great deal of attention has been paid to the elucidation of miRNAs involved in low nitrate stress. Methods: To identify miRNAs expressed under stress, cDNA libraries were analyzed. Results: We reported the nine potential miRNAs with 67 targets involved in nutrient transporters and other stress specific genes. Among the miRNA sequences obtained 6 sequences belong to miR172 family, one with miR169. RT-PCR analysis of expression of miR172 family was induced upon low nitrate stress while miR169 family was repressed. In addition, Pvu-SN7b and Pvu-miR16 may be new members of miRNA172 and miR169 families, respectively. Conclusion: The targets of Pvu-SN7b were major protein kinases, one among which is the Protein Kinase CK2. CK2 Kinase is found to involve in transcription-directed signaling, gene control and cell-cycle regulation. Other targets of Pvu-SN7b were involved in DNA-dependent transcription regulation, photo-periodism, calcium-mediated signaling. Pvu-miR16 targets Thymidine kinase, the key enzyme of deoxy-nucleotide synthesis. The cleavage of these targets affects cell proliferation there by affecting nodule formation. Pvu-miR8 inhibits translation of its target protein Pre-protein translocase, a membrane-bound protein transporter involved in trans-membrane protein transportation. Together these results denote the response and role of miRNAs to nitrate-limiting conditions in French bean.

Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

Computational framework to integrate the effect of antigen recognition on disease epidemiology outcome: multi-scale approach

Human Leukocyte Antigen (HLA) plays an important role, in presenting foreign pathogens to our immune system, there by eliciting early immune responses. HLA genes are highly polymorphic, giving rise to diverse antigen presentation capability. An important factor contributing to enormous variations in individual responses to diseases is differences in their HLA profiles. The heterogeneity in allele specific disease responses decides the overall disease epidemiological outcome. Here we propose an agent based computational framework, capable of incorporating allele specific information, to analyze disease epidemiology. This framework assumes a SIR model to estimate average disease transmission and recovery rate. Using epitope prediction tool, it performs sequence based epitope detection for a given the pathogenic genome and derives an allele specific disease susceptibility index depending on the epitope detection efficiency. The allele specific disease transmission rate, that follows, is then fed to the agent based epidemiology model, to analyze the disease outcome. The methodology presented here has a potential use in understanding how a disease spreads and effective measures to control the disease.

Identification of a novel AMPK-PEA15 axis in the anoikis-resistant growth of mammary cells

Introduction: Matrix detachment triggers anoikis, a form of apoptosis, in most normal epithelial cells, while acquisition of anoikis resistance is a prime requisite for solid tumor growth. Of note, recent studies have revealed that a small population of normal human mammary epithelial cells (HMECs) survive in suspension and generate multicellular spheroids termed `mammospheres'. Therefore, understanding how normal HMECs overcome anoikis may provide insights into breast cancer initiation and progression. Methods: Primary breast tissue-derived normal HMECs were grown as adherent monolayers or mammospheres. The status of AMP-activated protein kinase (AMPK) and PEA15 signaling was investigated by immunoblotting. Pharmacological agents and an RNA interference (RNAi) approach were employed to gauge their roles in mammosphere formation. Immunoprecipitation and in vitro kinase assays were undertaken to evaluate interactions between AMPK and PEA15. In vitro sphere formation and tumor xenograft assays were performed to understand their roles in tumorigenicity. Results: In this study, we show that mammosphere formation by normal HMECs is accompanied with an increase in AMPK activity. Inhibition or knockdown of AMPK impaired mammosphere formation. Concomitant with AMPK activation, we detected increased Ser(116) phosphorylation of PEA15, which promotes its anti-apoptotic functions. Inhibition or knockdown of AMPK impaired PEA15 Ser(116) phosphorylation and increased apoptosis. Knockdown of PEA15, or overexpression of the nonphosphorylatable S116A mutant of PEA15, also abrogated mammosphere formation. We further demonstrate that AMPK directly interacts with and phosphorylates PEA15 at Ser(116) residue, thus identifying PEA15 as a novel AMPK substrate. Together, these data revealed that AMPK activation facilitates mammosphere formation by inhibition of apoptosis, at least in part, through Ser(116) phosphorylation of PEA15. Since anoikis resistance plays a critical role in solid tumor growth, we investigated the relevance of these findings in the context of breast cancer. Significantly, we show that the AMPK-PEA15 axis plays an important role in the anchorage-independent growth of breast cancer cells both in vitro and in vivo. Conclusions: Our study identifies a novel AMPK-PEA15 signaling axis in the anchorage-independent growth of both normal and cancerous mammary epithelial cells, suggesting that breast cancer cells may employ mechanisms of anoikis resistance already inherent within a subset of normal HMECs. Thus, targeting the AMPK-PEA15 axis might prevent breast cancer dissemination and metastasis.

Computational Study of the Formation of Short Centrosymmetric N-center dot center dot center dot S Supramolecular Synthon and Related Weak Interactions in Crystalline 1,2,4-Triazoles

A comprehensive analysis of the crystal packing and the energetic features of a series of four biologically active molecules belonging to the family of substituted 4-(benzylideneamino)-3-(4-fluoro-3-phenoxyphenyl)-1H-1,2,4-triazole-5-(4 H)-thione derivatives have been performed based on the molecular conformation and the supramolecular packing. This involves the formation of a short centrosymmetric R-2(2)(8) NH...S supramolecular synthon in the solid state, including the presence of CH...S, CH...O, CH...N, CH...F, CH...Cl, CF...FC, CCl...ClC, and CH...pi intermolecular interactions along with pp stacking to evaluate the role of noncovalent interactions in the crystal. The presence of such synthons has a substantial contribution toward the interaction energy (-18 to -20 kcal/mol) as obtained from the PIXEL calculation, wherein the Coulombic and polarization contribution are more significant than the dispersion contribution. The geometrical characteristics of such synthons favor short distance, and the population of related molecules having these geometries is rare as has been obtained from the Cambridge Structural Database (CSD). Furthermore, their interaction energies have been compared with those present in our molecules in the solid state. The topological characteristics of the NH...S supramolecular synthon, in addition to related weak interactions, CH...N, CH...Cl, CF...FC, and CCl...ClC, have been estimated using the quantum theory of atoms in molecules (QTAIM). In addition, an analysis of the Hirshfeld surface and associated fingerprint plots of these four molecules also have provided a platform for the evaluation of the contribution of different atom...atom contacts, which contribute toward the packing of the molecules in solids.

Effect of alkyl substituents in BODIPYs: a comparative DFT computational investigation

Random changes in the alkyl substitution patterns of fluorescent dyes, e.g. BODIPYs, are often accompanied by significant changes in their photophysical properties. To understand such alterations in properties in closely related molecular systems, a comparative DFT (density functional theory) computational investigation was performed in order to comprehend the effects of alkyl substitution in controlling the structural and electronic nature of BODIPY dyes. In this context, a systematic strategy was utilized, considering all possible outcomes of constitutionally-isomeric molecules to understand the alkyl groups' effects on the BODIPY molecules. Four different computational methods {i.e. B3LYP/631G(d); B3LYP/6-311++ G(d,p); wb97xd/6-311++ G(d,p) and mpw1pw91/6-311++ G(d,p)} were employed to rationalize the agreement of the trends associated with the molecular properties. In line with experimental observations, it was found that alkyl substituents in BODIPY dyes situated at 3/5-positions effectively participate in stabilization as well as planarization of such molecules. Screening of all the possible isomeric molecular systems was used to understand the individual properties and overall effects of the typical alkyl substituents in controlling several basic properties of such BODIPY molecules.

From Workstations to Workbenches: Towards Predicting Physicochemically Viable Protein-Protein Interactions Across a Host and a Pathogen

The understanding of protein-protein interactions is indispensable in comprehending most of the biological processes in a cell. Small-scale experiments as well as large-scale high-throughput techniques over the past few decades have facilitated identification and analysis of protein-protein interactions which form the basis of much of our knowledge on functional and regulatory aspects of proteins. However, such rich catalog of interaction data should be used with caution when establishing protein-protein interactions in silico, as the high-throughput datasets are prone to false positives. Numerous computational means developed to pursue genome-wide studies on protein-protein interactions at times overlook the mechanistic and molecular details, thus questioning the reliability of predicted protein-protein interactions. We review the development, advantages, and shortcomings of varied approaches and demonstrate that by providing a structural viewpoint in terms of shape complementarity and interaction energies at protein-protein interfaces coupled with information on expression and localization of proteins homologous to an interacting pair, it is possible to assess the credibility of predicted interactions in biological context. With a focus on human pathogen Mycobacterium tuberculosis H37Rv, we show that such scrupulous use of details at the molecular level can predict physicochemically viable protein-protein interactions across host and pathogen. Such predicted interactions have the potential to provide molecular basis of probable mechanisms of pathogenesis and hence open up ways to explore their usefulness as targets in the light of drug discovery. (c) 2014 IUBMB Life, 66(11):759-774, 2014

Iterated stochastic filters with additive updates for dynamic system identification: Annealing-type iterations and the filter bank

A nonlinear stochastic filtering scheme based on a Gaussian sum representation of the filtering density and an annealing-type iterative update, which is additive and uses an artificial diffusion parameter, is proposed. The additive nature of the update relieves the problem of weight collapse often encountered with filters employing weighted particle based empirical approximation to the filtering density. The proposed Monte Carlo filter bank conforms in structure to the parent nonlinear filtering (Kushner-Stratonovich) equation and possesses excellent mixing properties enabling adequate exploration of the phase space of the state vector. The performance of the filter bank, presently assessed against a few carefully chosen numerical examples, provide ample evidence of its remarkable performance in terms of filter convergence and estimation accuracy vis-a-vis most other competing filters especially in higher dimensional dynamic system identification problems including cases that may demand estimating relatively minor variations in the parameter values from their reference states. (C) 2014 Elsevier Ltd. All rights reserved.

On Whitespace Identification Using Randomly Deployed Sensors

This work considers the identification of the available whitespace, i.e., the regions that do not contain any existing transmitter within a given geographical area. To this end, n sensors are deployed at random locations within the area. These sensors detect for the presence of a transmitter within their radio range r(s) using a binary sensing model, and their individual decisions are combined to estimate the available whitespace. The limiting behavior of the recovered whitespace as a function of n and r(s) is analyzed. It is shown that both the fraction of the available whitespace that the nodes fail to recover as well as their radio range optimally scale as log(n)/n as n gets large. The problem of minimizing the sum absolute error in transmitter localization is also analyzed, and the corresponding optimal scaling of the radio range and the necessary minimum transmitter separation is determined.

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.

Spatio-temporal correlations in aqueous systems: computational studies of static and dynamic heterogeneity by 2D-IR spectroscopy

Local heterogeneity is ubiquitous in natural aqueous systems. It can be caused locally by external biomolecular subsystems like proteins, DNA, micelles and reverse micelles, nanoscopic materials etc., but can also be intrinsic to the thermodynamic nature of the aqueous solution itself (like binary mixtures or at the gas-liquid interface). The altered dynamics of water in the presence of such diverse surfaces has attracted considerable attention in recent years. As these interfaces are quite narrow, only a few molecular layers thick, they are hard to study by conventional methods. The recent development of two dimensional infra-red (2D-IR) spectroscopy allows us to estimate length and time scales of such dynamics fairly accurately. In this work, we present a series of interesting studies employing two dimensional infra-red spectroscopy (2D-IR) to investigate (i) the heterogeneous dynamics of water inside reverse micelles of varying sizes, (ii) supercritical water near the Widom line that is known to exhibit pronounced density fluctuations and also study (iii) the collective and local polarization fluctuation of water molecules in the presence of several different proteins. The spatio-temporal correlation of confined water molecules inside reverse micelles of varying sizes is well captured through the spectral diffusion of corresponding 2D-IR spectra. In the case of supercritical water also, we observe a strong signature of dynamic heterogeneity from the elongated nature of the 2D-IR spectra. In this case the relaxation is ultrafast. We find remarkable agreement between the different tools employed to study the relaxation of density heterogeneity. For aqueous protein solutions, we find that the calculated dielectric constant of the respective systems unanimously shows a noticeable increment compared to that of neat water. However, the `effective' dielectric constant for successive layers shows significant variation, with the layer adjacent to the protein having a much lower value. Relaxation is also slowest at the surface. We find that the dielectric constant achieves the bulk value at distances more than 3 nm from the surface of the protein.