Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^

PIM KINASE INHIBITION: SINGLE AGENT AND COMBINATION APPROACHES IN B-CELL LYMPHOMA

Proviral integration site for Moloney murine leukemia virus (Pim) kinases are Ser/Thr/Tyr kinases. They modulate B-cell development but become oncoproteins and promote cancer development once overexpressed. Containing three isoforms, Pim-1, -2 and -3 are known to phosphorylate various substrates that regulate transcription, translation, cell cycle, and survival pathways in both hematological and solid tumors. Mantle cell lymphoma (MCL) is an aggressive B-cell lymphoma. Elevated Pim kinase levels are common in MCL, and it negatively correlates with patient outcome. SGI-1776 is a small molecule inhibitor selective for Pim-1/-3. We hypothesize that SGI-1776 treatment in MCL will inhibit Pim kinase function, and inhibition of downstream substrates phosphorylation will disrupt transcriptional, translational, and cell cycle processes while promoting apoptosis. SGI-1776 treatment induced moderate to high levels of apoptosis in four MCL cell lines (JeKo-1, Mino, SP-53 and Granta-519) and peripheral blood mononuclear cells (PBMCs) from MCL patients. Phosphorylation of transcription and translation regulators, c-Myc and 4E-BP1 declined in both model systems. Additionally, levels of short-lived Mcl-1 mRNA and protein also decreased and correlated with decline of global RNA synthesis. Collectively, our investigations highlight Pim kinases as viable drug targets in MCL and emphasize their roles in transcriptional and translational regulation. We further investigated a combination strategy using SGI-1776 with bendamustine, an FDA-approved DNA-damaging alkylating agent for treating non-Hodgkin’s lymphoma. We hypothesized this combination will enhance SGI-1776-induced transcription and translation inhibition, while promoting bendamustine-triggered DNA damage and inducing additive to synergistic cytotoxicity in B-cell lymphoma. Bendamustine alone resulted in moderate levels of apoptosis induction in MCL cell lines (JeKo-1 and Mino), and in MCL and splenic marginal zone lymphoma (a type of B-cell lymphoma) primary cells. An additive effect in cell killing was observed when combined with SGI-1776. Expectedly, SGI-1776 effectively decreased global RNA and protein synthesis levels, while bendamustine significantly inhibited DNA synthesis and generated DNA damage response. In combination, intensified inhibitory effects in DNA, RNA and protein syntheses were observed. Together, these data suggested feasibility of using Pim kinase inhibitor in combination with chemotherapeutic agents such as bendamustine in B-cell lymphoma, and provided foundation of their mechanism of actions in lymphoma cells.

Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Unveiling Protein Functions through the Dynamics of the Interaction Network

Protein interaction networks have become a tool to study biological processes, either for predicting molecular functions or for designing proper new drugs to regulate the main biological interactions. Furthermore, such networks are known to be organized in sub-networks of proteins contributing to the same cellular function. However, the protein function prediction is not accurate and each protein has traditionally been assigned to only one function by the network formalism. By considering the network of the physical interactions between proteins of the yeast together with a manual and single functional classification scheme, we introduce a method able to reveal important information on protein function, at both micro- and macro-scale. In particular, the inspection of the properties of oscillatory dynamics on top of the protein interaction network leads to the identification of misclassification problems in protein function assignments, as well as to unveil correct identification of protein functions. We also demonstrate that our approach can give a network representation of the meta-organization of biological processes by unraveling the interactions between different functional classes

Purification and Characterization of AsES Protein A Subtilisin Secreted by Acremonium Strictum is a Novel Plant Defense Elicitor.

In this work, the purification and characterization of an extracellular elicitor protein, designated AsES, produced by an avirulent isolate of the strawberry pathogen Acremonium strictum, are reported. The defense eliciting activity present in culture filtrates was recovered and purified by ultrafiltration (cutoff, 30 kDa), anionic exchange (Q-Sepharose, pH 7.5), and hydrophobic interaction (phenyl-Sepharose) chromatographies. Two-dimensional SDS-PAGE of the purified active fraction revealed a single spot of 34 kDa and pI 8.8. HPLC (C2/C18) and MS/MS analysis confirmed purification to homogeneity. Foliar spray with AsES provided a total systemic protection against anthracnose disease in strawberry, accompanied by the expression of defense-related genes (i.e. PR1 and Chi2-1). Accumulation of reactive oxygen species (e.g. H2O2 and O2̇̄) and callose was also observed in Arabidopsis. By using degenerate primers designed from the partial amino acid sequences and rapid amplification reactions of cDNA ends, the complete AsES-coding cDNA of 1167 nucleotides was obtained. The deduced amino acid sequence showed significant identity with fungal serine proteinases of the subtilisin family, indicating that AsES is synthesized as a larger precursor containing a 15-residue secretory signal peptide and a 90-residue peptidase inhibitor I9 domain in addition to the 283-residue mature protein. AsES exhibited proteolytic activity in vitro, and its resistance eliciting activity was eliminated when inhibited with PMSF, suggesting that its proteolytic activity is required to induce the defense response. This is, to our knowledge, the first report of a fungal subtilisin that shows eliciting activity in plants. This finding could contribute to develop disease biocontrol strategies in plants by activating its innate immunity.

Choroid plexus epithelial expression of MDR1 P glycoprotein and multidrug resistance-associated protein contribute to the blood–cerebrospinal-fluid drug-permeability barrier

The blood–brain barrier and a blood–cerebrospinal-fluid (CSF) barrier function together to isolate the brain from circulating drugs, toxins, and xenobiotics. The blood–CSF drug-permeability barrier is localized to the epithelium of the choroid plexus (CP). However, the molecular mechanisms regulating drug permeability across the CP epithelium are defined poorly. Herein, we describe a drug-permeability barrier in human and rodent CP mediated by epithelial-specific expression of the MDR1 (multidrug resistance) P glycoprotein (Pgp) and the multidrug resistance-associated protein (MRP). Noninvasive single-photon-emission computed tomography with 99mTc-sestamibi, a membrane-permeant radiopharmaceutical whose transport is mediated by both Pgp and MRP, shows a large blood-to-CSF concentration gradient across intact CP epithelium in humans in vivo. In rats, pharmacokinetic analysis with 99mTc-sestamibi determined the concentration gradient to be greater than 100-fold. In membrane fractions of isolated native CP from rat, mouse, and human, the 170-kDa Pgp and 190-kDa MRP are identified readily. Furthermore, the murine proteins are absent in CP isolated from their respective mdr1a/1b(−/−) and mrp(−/−) gene knockout littermates. As determined by immunohistochemical and drug-transport analysis of native CP and polarized epithelial cell cultures derived from neonatal rat CP, Pgp localizes subapically, conferring an apical-to-basal transepithelial permeation barrier to radiolabeled drugs. Conversely, MRP localizes basolaterally, conferring an opposing basal-to-apical drug-permeation barrier. Together, these transporters may coordinate secretion and reabsorption of natural product substrates and therapeutic drugs, including chemotherapeutic agents, antipsychotics, and HIV protease inhibitors, into and out of the central nervous system.

Functional identification and reconstitution of an odorant receptor in single olfactory neurons

The olfactory system is remarkable in its capacity to discriminate a wide range of odorants through a series of transduction events initiated in olfactory receptor neurons. Each olfactory neuron is expected to express only a single odorant receptor gene that belongs to the G protein coupled receptor family. The ligand–receptor interaction, however, has not been clearly characterized. This study demonstrates the functional identification of olfactory receptor(s) for specific odorant(s) from single olfactory neurons by a combination of Ca2+-imaging and reverse transcription–coupled PCR analysis. First, a candidate odorant receptor was cloned from a single tissue-printed olfactory neuron that displayed odorant-induced Ca2+ increase. Next, recombinant adenovirus-mediated expression of the isolated receptor gene was established in the olfactory epithelium by using green fluorescent protein as a marker. The infected neurons elicited external Ca2+ entry when exposed to the odorant that originally was used to identify the receptor gene. Experiments performed to determine ligand specificity revealed that the odorant receptor recognized specific structural motifs within odorant molecules. The odorant receptor-mediated signal transduction appears to be reconstituted by this two-step approach: the receptor screening for given odorant(s) from single neurons and the functional expression of the receptor via recombinant adenovirus. The present approach should enable us to examine not only ligand specificity of an odorant receptor but also receptor specificity and diversity for a particular odorant of interest.

Measurement of cytosolic, mitochondrial, and Golgi pH in single living cells with green fluorescent proteins

Many cellular events depend on a tightly compartmentalized distribution of H+ ions across membrane-bound organelles. However, measurements of organelle pH in living cells have been scarce. Several mutants of the Aequorea victoria green fluorescent protein (GFP) displayed a pH-dependent absorbance and fluorescent emission, with apparent pKa values ranging from 6.15 (mutations F64L/S65T/H231L) and 6.4 (K26R/F64L/S65T/Y66W/N146I/M153T/V163A/N164H/H231L) to a remarkable 7.1 (S65G/S72A/T203Y/H231L). We have targeted these GFPs to the cytosol plus nucleus, the medial/trans-Golgi by fusion with galactosyltransferase, and the mitochondrial matrix by using the targeting signal from subunit IV of cytochrome c oxidase. Cells in culture transfected with these cDNAs displayed the expected subcellular localization by light and electron microscopy and reported local pH that was calibrated in situ with ionophores. We monitored cytosolic and nuclear pH of HeLa cells, and mitochondrial matrix pH in HeLa cells and in rat neonatal cardiomyocytes. The pH of the medial/trans-Golgi was measured at steady-state (calibrated to be 6.58 in HeLa cells) and after various manipulations. These demonstrated that the Golgi membrane in intact cells is relatively permeable to H+, and that Cl− serves as a counter-ion for H+ transport and likely helps to maintain electroneutrality. The amenability to engineer GFPs to specific subcellular locations or tissue targets using gene fusion and transfer techniques should allow us to examine pH at sites previously inaccessible.

Agrobacterium VirD2 protein interacts with plant host cyclophilins

Agrobacterium tumefaciens induces crown gall tumors on plants by transferring a nucleoprotein complex, the T-complex, from the bacterium to the plant cell. The T-complex consists of T-DNA, a single-stranded DNA segment of the tumor-inducing plasmid, VirD2, an endonuclease covalently bound to the 5′ end of the T-DNA, and perhaps VirE2, a single-stranded DNA binding protein. The yeast two-hybrid system was used to screen for proteins interacting with VirD2 and VirE2 to identify components in Arabidopsis thaliana that interact with the T-complex. Three VirD2- and two VirE2-interacting proteins were identified. Here we characterize the interactions of VirD2 with two isoforms of Arabidopsis cyclophilins identified by using this analysis. The VirD2 domain interacting with the cyclophilins is distinct from the endonuclease, omega, and the nuclear localization signal domains. The VirD2–cyclophilin interaction is disrupted in vitro by cyclosporin A, which also inhibits Agrobacterium-mediated transformation of Arabidopsis and tobacco. These data strongly suggest that host cyclophilins play a role in T-DNA transfer.

RecA binding to a single double-stranded DNA molecule: A possible role of DNA conformational fluctuations

Most genetic regulatory mechanisms involve protein–DNA interactions. In these processes, the classical Watson–Crick DNA structure sometimes is distorted severely, which in turn enables the precise recognition of the specific sites by the protein. Despite its key importance, very little is known about such deformation processes. To address this general question, we have studied a model system, namely, RecA binding to double-stranded DNA. Results from micromanipulation experiments indicate that RecA binds strongly to stretched DNA; based on this observation, we propose that spontaneous thermal stretching fluctuations may play a role in the binding of RecA to DNA. This has fundamental implications for the protein–DNA binding mechanism, which must therefore rely in part on a combination of flexibility and thermal fluctuations of the DNA structure. We also show that this mechanism is sequence sensitive. Theoretical simulations support this interpretation of our experimental results, and it is argued that this is of broad relevance to DNA–protein interactions.

Protein thermostability above 100°C: A key role for ionic interactions

The discovery of hyperthermophilic microorganisms and the analysis of hyperthermostable enzymes has established the fact that multisubunit enzymes can survive for prolonged periods at temperatures above 100°C. We have carried out homology-based modeling and direct structure comparison on the hexameric glutamate dehydrogenases from the hyperthermophiles Pyrococcus furiosus and Thermococcus litoralis whose optimal growth temperatures are 100°C and 88°C, respectively, to determine key stabilizing features. These enzymes, which are 87% homologous, differ 16-fold in thermal stability at 104°C. We observed that an intersubunit ion-pair network was substantially reduced in the less stable enzyme from T. litoralis, and two residues were then altered to restore these interactions. The single mutations both had adverse effects on the thermostability of the protein. However, with both mutations in place, we observed a fourfold improvement of stability at 104°C over the wild-type enzyme. The catalytic properties of the enzymes were unaffected by the mutations. These results suggest that extensive ion-pair networks may provide a general strategy for manipulating enzyme thermostability of multisubunit enzymes. However, this study emphasizes the importance of the exact local environment of a residue in determining its effects on stability.

Low frequency vibrational modes in proteins: Changes induced by point-mutations in the protein-cofactor matrix of bacterial reaction centers

As a step toward understanding their functional role, the low frequency vibrational motions (<300 cm−1) that are coupled to optical excitation of the primary donor bacteriochlorophyll cofactors in the reaction center from Rhodobacter sphaeroides were investigated. The pattern of hydrogen-bonding interaction between these bacteriochlorophylls and the surrounding protein was altered in several ways by mutation of single amino acids. The spectrum of low frequency vibrational modes identified by femtosecond coherence spectroscopy varied strongly between the different reaction center complexes, including between different mutants where the pattern of hydrogen bonds was the same. It is argued that these variations are primarily due to changes in the nature of the individual modes, rather than to changes in the charge distribution in the electronic states involved in the optical excitation. Pronounced effects of point mutations on the low frequency vibrational modes active in a protein-cofactor system have not been reported previously. The changes in frequency observed indicate a strong involvement of the protein in these nuclear motions and demonstrate that the protein matrix can increase or decrease the fluctuations of the cofactor along specific directions.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Effectors of a developmental mitogen-activated protein kinase cascade revealed by expression signatures of signaling mutants

Despite the importance of mitogen-activated protein kinase (MAPK) signaling in eukaryotic biology, the mechanisms by which signaling yields phenotypic changes are poorly understood. We have combined transcriptional profiling with genetics to determine how the Kss1 MAPK signaling pathway controls dimorphic development in Saccharomyces cerevisiae. This analysis identified dozens of transcripts that are regulated by the pathway, whereas previous work had identified only a single downstream target, FLO11. One of the MAPK-regulated genes is PGU1, which encodes a secreted enzyme that hydrolyzes polygalacturonic acid, a structural barrier to microbial invasion present in the natural plant substrate of S. cerevisiae. A third key transcriptional target is the G1 cyclin gene CLN1, a morphogenetic regulator that we show to be essential for pseudohyphal growth. In contrast, the homologous CLN2 cyclin gene is dispensable for development. Thus, the Kss1 MAPK cascade programs development by coordinately modulating a cell adhesion factor, a secreted host-destroying activity, and a specialized subunit of the Cdc28 cyclin-dependent kinase.

Three functional luciferase domains in a single polypeptide chain

We report a unique case of a gene containing three homologous and contiguous repeat sequences, each of which, after excision, cloning, and expression in Escherichia coli, is shown to code for a peptide catalyzing the same reaction as the native protein, Gonyaulax polyedra luciferase (Mr = 137). This enzyme, which catalyzes the light-emitting oxidation of a linear tetrapyrrole (dinoflagellate luciferin), exhibits no sequence similarities to other luciferases in databases. Sequence analysis also reveals an unusual evolutionary feature of this gene: synonymous substitutions are strongly constrained in the central regions of each of the repeated coding sequences.