Dynamical systems for discovering protein complexes and functional modules from biological networks

Recent advances in high throughput experiments and annotations via published literature have provided a wealth of interaction maps of several biomolecular networks, including metabolic, protein-protein, and protein-DNA interaction networks. The architecture of these molecular networks reveals important principles of cellular organization and molecular functions. Analyzing such networks, i.e., discovering dense regions in the network, is an important way to identify protein complexes and functional modules. This task has been formulated as the problem of finding heavy subgraphs, the Heaviest k-Subgraph Problem (k-HSP), which itself is NPhard. However, any method based on the k-HSP requires the parameter k and an exact solution of k-HSP may still end up as a “spurious” heavy subgraph, thus reducing its practicability in analyzing large scale biological networks. We proposed a new formulation, called the rank-HSP, and two dynamical systems to approximate its results. In addition, a novel metric, called the Standard deviation and Mean Ratio (SMR), is proposed for use in “spurious” heavy subgraphs to automate the discovery by setting a fixed threshold. Empirical results on both the simulated graphs and biological networks have demonstrated the efficiency and effectiveness of our proposal.

Interacting amino acid preferences of 3D pattern pairs at the binding sites of transient and obligate protein complexes

To assess the physico-chemical characteristics of protein-protein interactions, protein sequences and overall structural folds have been analyzed previously. To highlight this, discovery and examination of amino acid patterns at the binding sites defined by structural proximity in 3-dimensional (3D) space are essential. In this paper, we investigate the interacting preferences of 3D pattern pairs discovered separately in transient and obligate protein complexes. These 3D pattern pairs are not necessarily sequence-consecutive, but each residue in two groups of amino acids from two proteins in a complex is within certain °A threshold to most residues in the other group. We develop an algorithm called AA-pairs by which every pair of interacting proteins is represented as a bipartite graph, and it discovers all maximal quasi-bicliques from every bipartite graph to form our 3D pattern pairs. From 112 and 2533 highly conserved 3D pattern pairs discovered in the transient and obligate complexes respectively, we observe that Ala and Leu is the highest occuring amino acid in interacting 3D patterns of transient (20.91%) and obligate (33.82%) complexes respectively. From the study on the dipeptide composition on each side of interacting 3D pattern pairs, dipeptides Ala-Ala and Ala-Leu are popular in 3D patterns of both transient and obligate complexes. The interactions between amino acids with large hydrophobicity difference are present more in the transient than in the obligate complexes. On contrary, in obligate complexes, interactions between hydrophobic residues account for the top 5 most occuring amino acid pairings.

Heterotrimeric G protein γ subunits provide functional selectivity in Gβγ dimer signaling in arabidopsis

The Arabidopsis thaliana heterotrimeric G protein complex is encoded by single canonical Galpha and Gbeta subunit genes and two Ggamma subunit genes (AGG1 and AGG2), raising the possibility that the two potential G protein complexes mediate different cellular processes. Mutants with reduced expression of one or both Ggamma genes revealed specialized roles for each Ggamma subunit. AGG1-deficient mutants, but not AGG2-deficient mutants, showed impaired resistance against necrotrophic pathogens, reduced induction of the plant defensin gene PDF1.2, and decreased sensitivity to methyl jasmonate. By contrast, both AGG1- and AGG2-deficient mutants were hypersensitive to auxin-mediated induction of lateral roots, suggesting that Gbetagamma1 and Gbetagamma2 synergistically inhibit auxin-dependent lateral root initiation. However, the involvement of each Ggamma subunit in this root response differs, with Gbetagamma1 acting within the central cylinder, attenuating acropetally transported auxin signaling, while Gbetagamma2 affects the action of basipetal auxin and graviresponsiveness within the epidermis and/or cortex. This selectivity also operates in the hypocotyl. Selectivity in Gbetagamma signaling was also found in other known AGB1-mediated pathways. agg1 mutants were hypersensitive to glucose and the osmotic agent mannitol during seed germination, while agg2 mutants were only affected by glucose. We show that both Ggamma subunits form functional Gbetagamma dimers and that each provides functional selectivity to the plant heterotrimeric G proteins, revealing a mechanism underlying the complexity of G protein-mediated signaling in plants.

Immunostimulatory principles from Chlorella pyrenoidosa – part 1 : isolation and biological assessment in vitro

Our proprietary preparation obtained by extraction of Chlorella pyrenoidosa cells, ONC-107 (Respondin™), was recently found to selectively boost antibody response to the influenza vaccine in a human clinical trial. Respondin™ is a potent stimulator of mouse B cell proliferation and an activator of macrophages. Bioactivity-guided resolution concluded that Respondin™ is composed of a mixture of immunostimulatory principles of different chemical nature. A combination of size exclusion, anion exchange and hydrophobic interaction chromatography revealed that the bulk of the immunostimulatory activity resides in polysaccharide/protein complexes with molecular masses larger than 100 kDa that are composed primarily of galactose, rhamnose and arabinose.

Identification of a novel polymorphism in the 3'UTR of the L-arginine transporter gene SLC7A1 : contribution to hypertension and endothelial dysfunction

Background— Endothelial dysfunction because of reduced nitric oxide bioavailability is a key feature of essential hypertension. We have found that normotensive siblings of subjects with essential hypertension have impaired endothelial function accompanied by altered arginine metabolism.

Methods and Results— We have identified a novel C/T polymorphism in the 3′UTR of the principal arginine transporter, solute carrier family 7 (cationic amino acid transporter, y+ system), member 1 gene (SLC7A1). The minor T allele significantly attenuates reporter gene expression (P<0.01) and is impaired in its capacity to form DNA-protein complexes (P<0.05). In 278 hypertensive subjects the frequency of the T allele was 13.3% compared with 7.6% in 498 normotensive subjects (P<0.001). Moreover, the overall genotype distribution observed in hypertensives differed significantly from that in normotensives (P<0.001). To complement these studies, we generated an endothelial-specific transgenic mouse overexpressing l-arginine transporter SLC7A1. The Slc7A1 transgenic mice exhibited significantly enhanced responses to the endothelium-dependent vasodilator acetylcholine (−log EC50 for wild-type versus Slc7A1 transgenic: 6.87±0.10 versus 7.56±0.13; P<0.001). This was accompanied by elevated production of nitric oxide by isolated aortic endothelial cells.

Conclusions— The present study identifies a key, functionally active polymorphism in the 3′UTR of SLC7A1. As such, this polymorphism may account for the apparent link between altered endothelial function, l-arginine, and nitric oxide metabolism and predisposition to essential hypertension.

The P-Loop domain of yeast Clp1 mediates interactions between CF IA and CPF factors in Pre-mRNA 3′ end formation

 Cleavage factor IA (CF IA), cleavage and polyadenylation factor (CPF), constitute major protein complexes required for pre-mRNA 3' end formation in yeast. The Clp1 protein associates with Pcf11, Rna15 and Rna14 in CF IA but its functional role remained unclear. Clp1 carries an evolutionarily conserved P-loop motif that was previously shown to bind ATP. Interestingly, human and archaean Clp1 homologues, but not the yeast protein, carry 5' RNA kinase activity. We show that depletion of Clp1 in yeast promoted defective 3' end formation and RNA polymerase II termination; however, cells expressing Clp1 with mutant P-loops displayed only minor defects in gene expression. Similarly, purified and reconstituted mutant CF IA factors that interfered with ATP binding complemented CF IA depleted extracts in coupled in vitro transcription/3' end processing reactions. We found that Clp1 was required to assemble recombinant CF IA and that certain P-loop mutants failed to interact with the CF IA subunit Pcf11. In contrast, mutations in Clp1 enhanced binding to the 3' endonuclease Ysh1 that is a component of CPF. Our results support a structural role for the Clp1 P-loop motif. ATP binding by Clp1 likely contributes to CF IA formation and cross-factor interactions during the dynamic process of 3' end formation.

Decreased fatty acit Â-oxidation in riboflavin-responsive, multiple acylocoenzyme a dehydrogenase-deficient patients is associated with an increase in uncoupling protein-3

Riboflavin-responsive, multiple acylcoenzyme A dehydrogenase deficiency (RR-MAD), a lipid storage myopathy, is characterized by, among others, a decrease in fatty acid (FA) ß-oxidation capacity. Muscle uncoupling protein 3 (UCP3) is up-regulated under conditions that either increase the levels of circulating free FA and/or decrease FA ß-oxidation. Using a relatively large cohort of seven RR-MAD patients, we aimed to better characterize the metabolic disturbances of this disease and to explore the possibility that it might increase UCP3 expression. A battery of biochemical and molecular tests were performed, which demonstrated decreases in FA ß-oxidation and in the activities of respiratory chain complexes I and II. These metabolic alterations were associated with increases of 3.1- and 1.7-fold in UCP3 mRNA and protein expression, respectively. All parameters were restored to control values after riboflavin treatment. We postulate that the up-regulation of UCP3 in RR-MAD is due to the accumulation of muscle FA/acylCoA. RR-MAD is an optimal model to support the hypothesis that UCP3 is involved in the outward translocation of an excess of FA from the mitochondria and to show that, in humans, the effects of FA on UCP3 expression are direct and independent of fatty acid ß-oxidation.

Thermodynamics of lipophilic drug binding to intestinal fatty acid binding protein and permeation across membranes

Intestinal fatty acid binding protein (I-FABP) is present at high levels in the absorptive cells of the intestine (enterocytes), where it plays a role in the intracellular solubilization of fatty acids (FA). However, I-FABP has also been shown to bind to a range of non-FA ligands, including some lipophilic drug molecules. Thus, in addition to its central role in FA trafficking, I-FABP potentially serves as an important intracellular carrier of lipophilic drugs. In this study we provide a detailed thermodynamic analysis of the binding and stability properties of I-FABP in complex with a series of fibrate and fenamate drugs to provide an insight into the forces driving drug binding to I-FABP. Drug binding and selectivity for I-FABP are driven by the interplay of protein−ligand interactions and solvent processes. The Gibbs free energies (ΔG°) determined from dissociation constants at 25 °C ranged from −6.2 to −10 kcal/mol. The reaction energetics indicate that drug binding to I-FABP is an enthalpy−entropy driven process. The relationship between I-FABP stability and drug binding affinity was examined by pulse proteolysis. There is a strong coupling between drug binding and I-FABP stability. The effect of an I-FABP protein sink on the kinetics and thermodynamics of tolfenamic acid permeation across an artificial phospholipid membrane were investigated. I-FABP significantly decreased the energy barrier for desorption of tolfenamic acid from the membrane into the acceptor compartment. Taken together, these data suggest that the formation of stable drug−I-FABP complexes is thermodynamically viable under conditions simulating the reactant concentrations likely observed in vivo and maybe a significant biochemical process that serves as a driving force for passive intestinal absorption of lipophilic drugs.

Probing the fibrate binding specificity of rat liver fatty acid binding protein

Liver-fatty acid binding protein (L-FABP) is found in high levels in enterocytes and is involved in cytosolic solubilization of fatty acids. In addition, L-FABP has been shown to bind endogenous and exogenous lipophilic compounds, suggesting that it may also play a role in modulating their absorption and disposition within enterocytes. Previously, we have described binding of L-FABP to a range of drugs, including a series of fibrates. In the present study, we have generated structural models of L-FABP-fibrate complexes and undertaken thermodynamic analysis of the binding of fibrates containing either a carboxylic acid or ester functionality. Analysis of the current data reveals that both the location and the energetics of binding are different for fibrates that contain a carboxylate compared to those that do not. As such, the data presented in this study suggest potential mechanisms that underpin molecular recognition and dictate specificity in the interaction between fibrates and L-FABP.

Irreversible and bivalent ligands for selected G-protein coupled receptors

Derivatives of pharmaceutical compounds that could potentially bind irreversibly to a hypertension controlling receptor were synthesised to provide tools for studying this receptor. Also compounds that simultaneously activate two cardiovascular receptors with opposing cellular mechanisms were synthesised. This resulted in a desired partial reduction in the activity of one receptor.

A critical role for the short intracellular C terminus in receptor activity-modifying protein function

Receptor activity-modifying proteins (RAMPs) interact with and modify the behavior of the calcitonin receptor (CTR) and calcitonin receptor-like receptor (CLR). We have examined the contribution of the short intracellular C terminus, using constructs that delete the last eight amino acids of each RAMP. C-Terminal deletion of individual RAMPs had little effect on the signaling profile induced when complexed with CLR in COS-7 or human embryonic kidney (HEK)293 cells. Likewise, confocal microscopy revealed each of the mutant RAMPs translocated hemagglutinin-tagged CLR to the cell surface. In contrast, a pronounced effect of RAMP C-terminal truncation was seen for RAMP/CTRa complexes, studied in COS-7 cells, with significant attenuation of amylin receptor phenotype induction that was stronger for RAMP1 and -2 than RAMP3. The loss of amylin binding upon C-terminal deletion could be partially recovered with overexpression of Gαs, suggesting an impact of the RAMP C terminus on coupling of G proteins to the receptor complex. In HEK293 cells the c-Myc-RAMP1 C-terminal deletion mutant showed high receptor-independent cell surface expression; however, this construct showed low cell surface expression when expressed alone in COS-7 cells, indicating interaction of RAMPs with other cellular components via the C terminus. This mutant also had reduced cell surface expression when coexpressed with CTR. Thus, this study reveals important functionality of the RAMP C-terminal domain and identifies key differences in the role of the RAMP C terminus for CTR versus CLR-based receptors.

Synthesis and biological evaluation of novel folic acid receptor-targeted, β-cyclodextrin-based drug complexes for cancer treatment

Drug targeting is an active area of research and nano-scaled drug delivery systems hold tremendous potential for the treatment of neoplasms. In this study, a novel cyclodextrin (CD)-based nanoparticle drug delivery system has been assembled and characterized for the therapy of folate receptor-positive [FR(+)] cancer. Water-soluble folic acid (FA)-conjugated CD carriers (FACDs) were successfully synthesized and their structures were confirmed by 1D/2D nuclear magnetic resonance (NMR), matrix-assisted laser desorption ionization time-of-flight mass spectrometer (MALDI-TOF-MS), high performance liquid chromatography (HPLC), Fourier transform infrared spectroscopy (FTIR), and circular dichroism. Drug complexes of adamatane (Ada) and cytotoxic doxorubicin (Dox) with FACD were readily obtained by mixed solvent precipitation. The average size of FACD-Ada-Dox was 1.5-2.5 nm. The host-guest association constant Ka was 1,639 M-1 as determined by induced circular dichroism and the hydrophilicity of the FACDs was greatly enhanced compared to unmodified CD. Cellular uptake and FR binding competitive experiments demonstrated an efficient and preferentially targeted delivery of Dox into FR-positive tumor cells and a sustained drug release profile was seen in vitro. The delivery of Dox into FR(+) cancer cells via endocytosis was observed by confocal microscopy and drug uptake of the targeted nanoparticles was 8-fold greater than that of non-targeted drug complexes. Our docking results suggest that FA, FACD and FACD-Ada-Dox could bind human hedgehog interacting protein that contains a FR domain. Mouse cardiomyocytes as well as fibroblast treated with FACD-Ada-Dox had significantly lower levels of reactive oxygen species, with increased content of glutathione and glutathione peroxidase activity, indicating a reduced potential for Dox-induced cardiotoxicity. These results indicate that the targeted drug complex possesses high drug association and sustained drug release properties with good biocompatibility and physiological stability. The novel FA-conjugated β-CD based drug complex might be promising as an anti-tumor treatment for FR(+) cancer.

PTEX is an essential nexus for protein export in malaria parasites

During the blood stages of malaria, several hundred parasite-encoded proteins are exported beyond the double-membrane barrier that separates the parasite from the host cell cytosol. These proteins have a variety of roles that are essential to virulence or parasite growth. There is keen interest in understanding how proteins are exported and whether common machineries are involved in trafficking the different classes of exported proteins. One potential trafficking machine is a protein complex known as the Plasmodium translocon of exported proteins (PTEX). Although PTEX has been linked to the export of one class of exported proteins, there has been no direct evidence for its role and scope in protein translocation. Here we show, through the generation of two parasite lines defective for essential PTEX components (HSP101 or PTEX150), and analysis of a line lacking the non-essential component TRX2 (ref. 12), greatly reduced trafficking of all classes of exported proteins beyond the double membrane barrier enveloping the parasite. This includes proteins containing the PEXEL motif (RxLxE/Q/D) and PEXEL-negative exported proteins (PNEPs). Moreover, the export of proteins destined for expression on the infected erythrocyte surface, including the major virulence factor PfEMP1 in Plasmodium falciparum, was significantly reduced in PTEX knockdown parasites. PTEX function was also essential for blood-stage growth, because even a modest knockdown of PTEX components had a strong effect on the parasite's capacity to complete the erythrocytic cycle both in vitro and in vivo. Hence, as the only known nexus for protein export in Plasmodium parasites, and an essential enzymic machine, PTEX is a prime drug target.

Complex coacervation between flaxseed protein isolate and flaxseed gum

Flaxseed protein isolate (FPI) and flaxseed gum (FG) were extracted, and the electrostatic complexation between these two biopolymers was studied as a function of pH and FPI-to-FG ratio using turbidimetric and electrophoretic mobility (zeta potential) tests. The zeta potential values of FPI, FG, and their mixtures at the FPI-to-FG ratios of 1:1, 3:1, 5:1, 10:1, 15:1 were measured over a pH range 8.0-1.5. The alteration of the secondary structure of FPI as a function of pH was studied using circular dichroism. The proportion of a-helical structure decreased, whereas both β-sheet structure and random coil structure increased with the lowering of pH from 8.0 to 3.0. The acidic pH affected the secondary structure of FPI and the unfolding of helix conformation facilitated the complexation of FPI with FG. The optimum FPI-to-FG ratio for complex coacervation was found to be 3:1. The critical pH values associated with the formation of soluble (pHc) and insoluble (pHΦ1) complexes at the optimum FPI-to-FG ratio were found to be 6.0 and 4.5, respectively. The optimum pH (pHopt) for the optimum complex coacervation was 3.1. The instability and dissolution of FPI-FG complex coacervates started (pHΦ2) at pH2.1. These findings contribute to the development of FPI-FG complex coacervates as delivery vehicles for unstable albeit valuable nutrients such as omega-3 fatty acids.