Comparison between computational alanine scanning and per-residue binding free energy decomposition for protein-protein association using MM-GBSA: application to the TCR-p-MHC complex.

Recognition by the T-cell receptor (TCR) of immunogenic peptides (p) presented by Class I major histocompatibility complexes (MHC) is the key event in the immune response against virus-infected cells or tumor cells. A study of the 2C TCR/SIYR/H-2K(b) system using a computational alanine scanning and a much faster binding free energy decomposition based on the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method is presented. The results show that the TCR-p-MHC binding free energy decomposition using this approach and including entropic terms provides a detailed and reliable description of the interactions between the molecules at an atomistic level. Comparison of the decomposition results with experimentally determined activity differences for alanine mutants yields a correlation of 0.67 when the entropy is neglected and 0.72 when the entropy is taken into account. Similarly, comparison of experimental activities with variations in binding free energies determined by computational alanine scanning yields correlations of 0.72 and 0.74 when the entropy is neglected or taken into account, respectively. Some key interactions for the TCR-p-MHC binding are analyzed and some possible side chains replacements are proposed in the context of TCR protein engineering. In addition, a comparison of the two theoretical approaches for estimating the role of each side chain in the complexation is given, and a new ad hoc approach to decompose the vibrational entropy term into atomic contributions, the linear decomposition of the vibrational entropy (LDVE), is introduced. The latter allows the rapid calculation of the entropic contribution of interesting side chains to the binding. This new method is based on the idea that the most important contributions to the vibrational entropy of a molecule originate from residues that contribute most to the vibrational amplitude of the normal modes. The LDVE approach is shown to provide results very similar to those of the exact but highly computationally demanding method.

Gene Ontology Function prediction in Mollicutes using Protein-Protein Association Networks

Abstract Background: Many complex systems can be represented and analysed as networks. The recent availability of large-scale datasets, has made it possible to elucidate some of the organisational principles and rules that govern their function, robustness and evolution. However, one of the main limitations in using protein-protein interactions for function prediction is the availability of interaction data, especially for Mollicutes. If we could harness predicted interactions, such as those from a Protein-Protein Association Networks (PPAN), combining several protein-protein network function-inference methods with semantic similarity calculations, the use of protein-protein interactions for functional inference in this species would become more potentially useful. Results: In this work we show that using PPAN data combined with other approximations, such as functional module detection, orthology exploitation methods and Gene Ontology (GO)-based information measures helps to predict protein function in Mycoplasma genitalium. Conclusions: To our knowledge, the proposed method is the first that combines functional module detection among species, exploiting an orthology procedure and using information theory-based GO semantic similarity in PPAN of the Mycoplasma species. The results of an evaluation show a higher recall than previously reported methods that focused on only one organism network.

Contribution of transmembrane and cytoplasmic domains to membrane protein association

Membrane proteins are critical to every aspect of cell physiology, with their association mediating important biological functions. The transmembrane and cytoplasmic domains are known to be important for their association. In order to characterize their role in detail, we have applied different biophysical techniques in detergent micelles to two model systems. The first study involves FcRγ, a single transmembrane domain protein existing as a disulfide linked homodimer. We investigated the role of a conserved transmembrane polar residue and the cytoplasmic tail in FcRγ homo-interactions. Our results by various biophysical techniques including SDS-PAGE, circular dichroism and sedimentation equilibrium in detergent micelles indicate importance of both the transmembrane polar residue and cytoplasmic tail in maintaining proper conformation for FcRγ homo-interactions. A contrasting second study was on L-selectin, another single transmembrane domain protein with a large extracellular domain and a short cytoplasmic tail. Previous cross-linking experiments indicate its possible dimerization. However, the purified fragment of L-selectin and corresponding mutants did not dimerize when analyzed by TOXCAT assay, sedimentation equilibrium and fluorescence resonance energy transfer. It was likely that the presence of juxtamembrane positively charged residues led to decreased migrational rates in SDS PAGE. In conclusion, complementary biophysical techniques should be used with care when studying membrane protein association in detergent micelles. As an extension to our study on L-selectin, we also investigated its interaction with Calmodulin (CaM) in detergent micelles. CaM was found to interact with different detergents. We applied fluorescence and NMR spectroscopy to characterize the interaction of both the apo and Ca 2+ bound form of CaM, with commonly used detergents, below and above their respective critical micelle concentrations. The interaction of apo-CaM with detergents was found to vary with the nature of the detergent head group, whereas Ca2+-CaM interacted with individual detergent molecules irrespective of the nature of their head group. NMR titration experiments of CaM with detergents indicated involvement of specific residues from the N-lobe, linker and C-lobe of CaM. ^

Controlling protein association and subcellular localization with a synthetic ligand that induces heterodimerization of proteins.

Extracellular growth and differentiation factors induce changes in gene expression in the nucleus by initiating a series of protein associations that alter the subcellular localization of intracellular signaling proteins. Initial events involve receptor homo- or heterodimerization and subsequent recruitment of cytosolic signaling proteins to the inner leaflet of the plasma membrane. Intermediate events involve the translocation of proteins into the nucleus. Late events involve the recruitment of transcriptional activators to the vicinity of specific genes in the nucleus, resulting in increased gene transcription. The ability to induce signals at each of these three phases of signaling pathways is illustrated by the use of a heterodimeric chemical inducer of dimerization that causes a proximal relationship between two different target proteins.

Molecular computations for reactions and phase transitions: applications to protein stabilization, hydrates and catalysis

In this work we have made significant contributions in three different areas of interest: therapeutic protein stabilization, thermodynamics of natural gas clathrate-hydrates, and zeolite catalysis. In all three fields, using our various computational techniques, we have been able to elucidate phenomena that are difficult or impossible to explain experimentally. More specifically, in mixed solvent systems for proteins we developed a statistical-mechanical method to model the thermodynamic effects of additives in molecular-level detail. It was the first method demonstrated to have truly predictive (no adjustable parameters) capability for real protein systems. We also describe a novel mechanism that slows protein association reactions, called the “gap effect.” We developed a comprehensive picture of methioine oxidation by hydrogen peroxide that allows for accurate prediction of protein oxidation and provides a rationale for developing strategies to control oxidation. The method of solvent accessible area (SAA) was shown not to correlate well with oxidation rates. A new property, averaged two-shell water coordination number (2SWCN) was identified and shown to correlate well with oxidation rates. Reference parameters for the van der Waals Platteeuw model of clathrate-hydrates were found for structure I and structure II. These reference parameters are independent of the potential form (unlike the commonly used parameters) and have been validated by calculating phase behavior and structural transitions for mixed hydrate systems. These calculations are validated with experimental data for both structures and for systems that undergo transitions from one structure to another. This is the first method of calculating hydrate thermodynamics to demonstrate predictive capability for phase equilibria, structural changes, and occupancy in pure and mixed hydrate systems. We have computed a new mechanism for the methanol coupling reaction to form ethanol and water in the zeolite chabazite. The mechanism at 400°C proceeds via stable intermediates of water, methane, and protonated formaldehyde.

Protein Assembly onto Cationic Supported Bilayers

Cationic supported bilayers on latex are useful to isolate and immobilize oppositely charged proteins as a monomolecular layer over a range of low protein concentrations and particle number densities. Cholera toxin (CT) from Vibrio cholerae, an 87 kDa AB(5) hexameric protein and bovine serum albumin (BSA) self-assembled on dioctadecyldimethylammonium bromide (DODAB) supported bilayers with high affinity yielding highly organized and monodisperse particulates at 5 x 10(9) particles/mL, over a range of low protein concentrations (0-0.025 mg/mL BSA or CT). Protein association onto the bilayer-covered polystyrene sulfate (PSS) was determined from adsorption isotherms, dynamic light scattering for size distributions and zeta-potential analysis revealing a monomolecular, thin and highly organized protein layer surrounding each particle with potential for biospecific recognition such as antigen-antibody, receptor-ligand, hybridization of oligonucleotide sequences, all of them important in immunodiagnosis, selective biomolecular chromatographic separations, microarrays design and others.

A protein-binding domain, EH, identified in the receptor tyrosine kinase substrate Eps15 and conserved in evolution.

In this report we structurally and functionally define a binding domain that is involved in protein association and that we have designated EH (for Eps15 homology domain). This domain was identified in the tyrosine kinase substrate Eps15 on the basis of regional conservation with several heterogeneous proteins of yeast and nematode. The EH domain spans about 70 amino acids and shows approximately 60% overall amino acid conservation. We demonstrated the ability of the EH domain to specifically bind cytosolic proteins in normal and malignant cells of mesenchymal, epithelial, and hematopoietic origin. These observations prompted our search for additional EH-containing proteins in mammalian cells. Using an EH domain-specific probe derived from the eps15 cDNA, we cloned and characterized a cDNA encoding an EH-containing protein with overall similarity to Eps15; we designated this protein Eps15r (for Eps15-related). Structural comparison of Eps15 and Eps15r defines a family of signal transducers possessing extensive networking abilities including EH-mediated binding and association with Src homology 3-containing proteins.

The association of Na,K-ATPase subunits studied by circular dichroism, surface tension and dilatational elasticity

Different stoichiometries are observed between alpha and beta subunits of Na,K-ATPase that depend on the method employed to solubilize and purify the enzyme. It is not known whether this variability is due to loss of protein-protein association, or is a result of the replacement of essential phospholipids by detergent molecules. With the aim of understanding the effect of enzyme/surfactant ratio on both the catalytic activity and the enzyme structure, we have investigated the bulk and surface properties of the enzyme. The circular dichroism (CD) spectra, surface tension and dilatational surface elasticity results were compared with the residual ATPase activity of the Na,K-ATPase in different surfactant and protein concentrations. Na,K-ATPase in the (alpha beta)(2) form dissociated to the alpha beta form on dilution, and associated to the (alpha beta)(4) form when concentrated. These different stoichiometries have similar ATPase activities and are in equilibrium at C(12)E(8) concentrations below the CIVIC (0.053 mg mL(-1)). At detergent concentrations above the CIVIC the ATPase activity of all forms was abolished, which is concomitant with the dissociation of the a and subunits. (C) 2008 Elsevier Inc. All rights reserved.

Association between single-nucleotide polymorphisms and milk production traits in buffalo

The aim of this study was to identify single-nucleotide polymorphisms (SNPs) in buffaloes associated with milk yield and content, in addition to somatic cell scores based on the cross-species transferability of SNPs from cattle to buffalo. A total of 15,745 SNPs were analyzed, of which 1562 showed 1% significance and 4742 with 5% significance, which were associated for all traits studied. After application of Bonferroni's correction for multiple tests of the traits analyzed, we found 2 significant SNPs placed on cattle chromosomes BTA15 and BTA20, which are homologous to buffalo chromosomes BBU16 and BBU19, respectively. In this genome association study, we found several significant SNPs affecting buffalo milk production and quality. Furthermore, the use of the high-density bovine BeadChip was suitable for genomic analysis in buffaloes. Although extensive chromosome arm homology was described between cattle and buffalo, the exact chromosomal position of SNP markers associated with these economically important traits in buffalo can be determined only through buffalo genome sequencing.

Bovine sperm cells viability during incubation with or without exogenous DNA

The aim of this study was to assess the effect of exogenous DNA and incubation time on the viability of bovine sperm. Sperm were incubated at a concentration of 5 x 10(6)/ml with or without plasmid pEYFP-NUC. Fluorescent probes, propidium iodide/Hoechst 33342, FITC-PSA and JC-1, were used to assess plasma membrane integrity (PMI), acrosome membrane integrity (AMI) and mitochondrial membrane potential (MMP) respectively at 0, 1, 2, 3 and 4 h of incubation. Exogenous DNA addition did not affect sperm viability; however, incubation time was related to sperm deterioration. Simultaneous assessment of PMI, AMI and MMP showed a reduction in the number of sperm with higher viability (integrity of plasma and acrosome membranes and high mitochondrial membrane potential) from 58.7% at 0 h to 7.5% after 4 h of incubation. Lower viability sperm (damaged plasma and acrosome membranes and low mitochondrial membrane potential) increased from 4.6% at 0 h to 25.99% after 4 h of incubation. When PMI, AMI and MMP were assessed separately we noticed a reduction in plasma and acrosome membrane integrity and mitochondrial membrane potential throughout the incubation period. Therefore, exogenous DNA addition does not affect sperm viability, but the viability is reduced by incubation time.

Distributed all-atom simulations of intrinsically unstructured proteins

The Computational Biophysics Group at the Universitat Pompeu Fabra (GRIB-UPF) hosts two unique computational resources dedicated to the execution of large scale molecular dynamics (MD) simulations: (a) the ACMD molecular-dynamics software, used on standard personal computers with graphical processing units (GPUs); and (b) the GPUGRID. net computing network, supported by users distributed worldwide that volunteer GPUs for biomedical research. We leveraged these resources and developed studies, protocols and open-source software to elucidate energetics and pathways of a number of biomolecular systems, with a special focus on flexible proteins with many degrees of freedom. First, we characterized ion permeation through the bactericidal model protein Gramicidin A conducting one of the largest studies to date with the steered MD biasing methodology. Next, we addressed an open problem in structural biology, the determination of drug-protein association kinetics; we reconstructed the binding free energy, association, and dissaciociation rates of a drug like model system through a spatial decomposition and a Makov-chain analysis. The work was published in the Proceedings of the National Academy of Sciences and become one of the few landmark papers elucidating a ligand-binding pathway. Furthermore, we investigated the unstructured Kinase Inducible Domain (KID), a 28-peptide central to signalling and transcriptional response; the kinetics of this challenging system was modelled with a Markovian approach in collaboration with Frank Noe’s group at the Freie University of Berlin. The impact of the funding includes three peer-reviewed publication on high-impact journals; three more papers under review; four MD analysis components, released as open-source software; MD protocols; didactic material, and code for the hosting group.

The effect of hydrostatic pressure on membrane-bound proteins

Many cellular proteins are bound to the surfaces of membranes and participate in various cell signaling responses. Interactions between this group of proteins are in part controlled by the membrane surface to which the proteins are bound. This review focuses on the effects of pressure on membrane-associated proteins. Initially, the effect of pressure on membrane surfaces and how pressure may perturb the membrane binding of proteins is discussed. Next, the effect of pressure on the activity and lateral association of proteins is considered. We then discuss how pressure can be used to gain insight into these types of proteins.

RP1 is a phosphorylation target of CK2 and is involved in cell adhesion

RP1 (synonym: MAPRE2, EB2) is a member of the microtubule binding EB1 protein family, which interacts with APC, a key regulatory molecule in the Wnt signalling pathway. While the other EB1 proteins are well characterized the cellular function and regulation of RP1 remain speculative to date. However, recently RP1 has been implicated in pancreatic cancerogenesis. CK2 is a pleiotropic kinase involved in adhesion, proliferation and anti-apoptosis. Overexpression of protein kinase CK2 is a hallmark of many cancers and supports the malignant phenotype of tumor cells. In this study we investigate the interaction of protein kinase CK2 with RP1 and demonstrate that CK2 phosphorylates RP1 at Ser(236) in vitro. Stable RP1 expression in cell lines leads to a significant cleavage and down-regulation of N-cadherin and impaired adhesion. Cells expressing a Phospho-mimicking point mutant RP1-ASP(236) show a marked decrease of adhesion to endothelial cells under shear stress. Inversely, we found that the cells under shear stress downregulate endogenous RP1, most likely to improve cellular adhesion. Accordingly, when RP1 expression is suppressed by shRNA, cells lacking RP1 display significantly increased cell adherence to surfaces. In summary, RP1 phosphorylation at Ser(236) by CK2 seems to play a significant role in cell adhesion and might initiate new insights in the CK2 and EB1 family protein association.

Common target genes of palatal and gingival fibroblasts for EMD: the microarray approach.

BACKGROUND AND OBJECTIVE Connective tissue grafts are frequently applied, together with Emdogain(®) , for root coverage. However, it is unknown whether fibroblasts from the gingiva and from the palate respond similarly to Emdogain. The aim of this study was therefore to evaluate the effect of Emdogain(®) on fibroblasts from palatal and gingival connective tissue using a genome-wide microarray approach. MATERIAL AND METHODS Human palatal and gingival fibroblasts were exposed to Emdogain(®) and RNA was subjected to microarray analysis followed by gene ontology screening with Database for Annotation, Visualization and Integrated Discovery functional annotation clustering, Kyoto Encyclopedia of Genes and Genomes pathway analysis and the Search Tool for the Retrieval of Interacting Genes/Proteins functional protein association network. Microarray results were confirmed by quantitative RT-PCR analysis. RESULTS The transcription levels of 106 genes were up-/down-regulated by at least five-fold in both gingival and palatal fibroblasts upon exposure to Emdogain(®) . Gene ontology screening assigned the respective genes into 118 biological processes, six cellular components, eight molecular functions and five pathways. Among the striking patterns observed were the changing expression of ligands targeting the transforming growth factor-beta and gp130 receptor family as well as the transition of mesenchymal epithelial cells. Moreover, Emdogain(®) caused changes in expression of receptors for chemokines, lipids and hormones, and for transcription factors such as SMAD3, peroxisome proliferator-activated receptor gamma and those of the ETS family. CONCLUSION The present data suggest that Emdogain(®) causes substantial alterations in gene expression, with similar patterns observed in palatal and gingival fibroblasts.

Salivary Pellicle vs Whole Saliva: The Response of Oral Fibroblasts Based on a Genome-Wide Microarray.

PURPOSE Whole saliva comprises components of the salivary pellicle that spontaneously forms on surfaces of implants and teeth. However, there are no studies that functionally link the salivary pellicle with a possible change in gene expression. MATERIALS AND METHODS This study examined the genetic response of oral fibroblasts exposed to the salivary pellicle and whole saliva. Oral fibroblasts were seeded onto a salivary pellicle and the respective untreated surface. Oral fibroblasts were also exposed to freshly harvested sterile-filtered whole saliva. A genome-wide microarray of oral fibroblasts was performed, followed by gene ontology screening with DAVID functional annotation clustering, KEGG pathway analysis, and the STRING functional protein association network. RESULTS Exposure of oral fibroblasts to saliva caused 61 genes to be differentially expressed (P < .05). Gene ontology screening assigned the respective genes into 262 biologic processes, 3 cellular components, 13 molecular functions, and 7 pathways. Most remarkable was the enrichment in the inflammatory response. None of the genes regulated by whole saliva was significantly changed when cells were placed onto a salivary pellicle. CONCLUSION The salivary pellicle per se does not provoke a significant inflammatory response of oral fibroblasts in vitro, whereas sterile-filtered whole saliva does produce a strong inflammatory response.