Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.

New measures for estimating surface complementarity and packing at protein-protein interfaces

A number of methods exist that use different approaches to assess geometric properties like the surface complementarity and atom packing at the protein-protein interface. We have developed two new and conceptually different measures using the Delaunay tessellation and interface slice selection to compute the surface complementarity and atom packing at the protein-protein interface in a straightforward manner. Our measures show a strong correlation among themselves and with other existing measures, and can be calculated in a highly time-efficient manner. The measures are discriminative for evaluating biological, as well as non-biological protein-protein contacts, especially from large protein complexes and large-scale structural studies(http://pallab.serc. iisc.ernet.in/nip_nsc). (C) 201 Federation of European Biochemical Societies. Published by Elsevier B. V. All rights reserved.

Interrelations between Dry Eye Syndrome and Tear Fluid Phospholipid Transfer Protein

The simplified model of human tear fluid (TF) is a three-layered structure composed of a homogenous gel-like layer of hydrated mucins, an aqueous phase, and a lipid-rich outermost layer found in the tear-air interface. It is assumed that amphiphilic phospholipids are found adjacent to the aqueous-mucin layer and externally to this a layer composed of non-polar lipids face the tear-air interface. The lipid layer prevents evaporation of the TF and protects the eye, but excess accumulation of lipids may lead to drying of the corneal epithelium. Thus the lipid layer must be controlled and maintained by some molecular mechanisms. In the circulation, phospholipid transfer protein (PLTP) and cholesteryl ester transfer protein (CETP) mediate lipid transfers. The aim of this thesis was to investigate the presence and molecular mechanisms of lipid transfer proteins in human TF. The purpose was also to study the role of these proteins in the development of dry eye syndrome (DES). The presence of TF PLTP and CETP was studied by western blotting and mass spectrometry. The concentration of these proteins was determined by ELISA. The activities of the enzymes were determined by specific lipid transfer assays. To study the molecular mechanisms involved in PLTP mediated lipid transfer Langmuir monolayers and asymmetrical flow field-flow fractionation (AsFlFFF) was used. Ocular tissue samples were stained with monoclonal antibodies against PLTP to study the secretion route of PLTP. Heparin-Sepharose affinity chromatography was used for PLTP pull-down experiments and co-eluted proteins were identified with MALDI-TOF mass spectrometry or Western blot analysis. To study whether PLTP plays any functional role in TF PLTP-deficient mice were examined. The activity of PLTP was also studied in dry eye patients. PLTP is a component of normal human TF, whereas CETP is not. TF PLTP concentration was about 2-fold higher than that in human plasma. Inactivation of PLTP by heat treatment or immunoinhibition abolished the phospholipid transfer activity in tear fluid. PLTP was found to be secreted from lacrimal glands. PLTP seems to be surface active and is capable of accepting lipid molecules without the presence of lipid-protein complexes. The active movement of radioactively labeled lipids and high activity form of PLTP to acceptor particles suggested a shuttle model of PLTP-mediated lipid transfer. In this model, PLTP physically transports lipids between the donor and acceptor. Protein-protein interaction assays revealed ocular mucins as PLTP interaction partners in TF. In mice with a full deficiency of functional PLTP enhanced corneal epithelial damage, increased corneal permeability to carboxyfluorescein, and decreased corneal epithelial occludin expression was demonstrated. Increased tear fluid PLTP activity was observed among human DES patients. These results together suggest a scavenger property of TF PLTP: if the corneal epithelium is contaminated by hydrophobic material, PLTP could remove them and transport them to the superficial layer of the TF or, alternatively, transport them through the naso-lacrimal duct. Thus, PLTP might play an integral role in tear lipid trafficking and in the protection of the corneal epithelium. The increased PLTP activity in human DES patients suggests an ocular surface protective role for this lipid transfer protein.

Using Correlated Parameters for Improved Ranking of Protein-Protein Docking Decoys

A successful protein-protein docking study culminates in identification of decoys at top ranks with near-native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near-native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein-protein interface. In this work, we describe a novel method combining an interface-size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein-protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having <= 10 angstrom backbone root mean square deviation from true binding geometry. Comparisons with other state-of-art methods confirm the improved ranking power of our method without the use of any experiment-guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less-difficult bound binary protein-protein complexes with <= 35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with <= 5 angstrom backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 32: 787-796, 2011.

DockYard-a repository to assist modeling of protein-protein docking

In the absence of interlogs, building docking models is a time intensive task, involving generation of a large pool of docking decoys followed by refinement and screening to identify near native docking solutions. This limits the researcher interested in building docking methods with the choice of benchmarking only a limited number of protein complexes. We have created a repository called dockYard (http://pallab.serc.iisc.ernet.in/dockYard), that allows modelers interested in protein-protein interaction to access large volume of information on protein dimers and their interlogs, and also download decoys for their work if they are interested in building modeling methods. dockYard currently offers four categories of docking decoys derived from: Bound (native dimer co-crystallized), Unbound (individual subunits are crystallized, as well as the target dimer), Variants (match the previous two categories in at least one subunit with 100% sequence identity), and Interlogs (match the previous categories in at least one subunit with >= 90% or >= 50% sequence identity). The web service offers options for full or selective download based on search parameters. Our portal also serves as a repository to modelers who may want to share their decoy sets with the community.

Pleiotropic Roles of a Ribosomal Protein in Dictyostelium discoideum

The cell cycle phase at starvation influences post-starvation differentiation and morphogenesis in Dictyostelium discoideum. We found that when expressed in Saccharomyces cerevisiae, a D. discoideum cDNA that encodes the ribosomal protein S4 (DdS4) rescues mutations in the cell cycle genes cdc24, cdc42 and bem1. The products of these genes affect morphogenesis in yeast via a coordinated moulding of the cytoskeleton during bud site selection. D. discoideum cells that over-or under-expressed DdS4 did not show detectable changes in protein synthesis but displayed similar developmental aberrations whose intensity was graded with the extent of over-or under-expression. This suggested that DdS4 might influence morphogenesis via a stoichiometric effect - specifically, by taking part in a multimeric complex similar to the one involving Cdc24p, Cdc42p and Bem1p in yeast. In support of the hypothesis, the S. cerevisiae proteins Cdc24p, Cdc42p and Bem1p as well as their D. discoideum cognates could be co-precipitated with antibodies to DdS4. Computational analysis and mutational studies explained these findings: a C-terminal domain of DdS4 is the functional equivalent of an SH3 domain in the yeast scaffold protein Bem1p that is central to constructing the bud site selection complex. Thus in addition to being part of the ribosome, DdS4 has a second function, also as part of a multi-protein complex. We speculate that the existence of the second role can act as a safeguard against perturbations to ribosome function caused by spontaneous variations in DdS4 levels.

An Exported Heat Shock Protein 40 Associates with Pathogenesis-Related Knobs in Plasmodium falciparum Infected Erythrocytes

Cell surface structures termed knobs are one of the most important pathogenesis related protein complexes deployed by the malaria parasite Plasmodium falciparum at the surface of the infected erythrocyte. Despite their relevance to the disease, their structure, mechanisms of traffic and their process of assembly remain poorly understood. In this study, we have explored the possible role of a parasite-encoded Hsp40 class of chaperone, namely PFB0090c/PF3D7_0201800 (KAHsp40) in protein trafficking in the infected erythrocyte. We found the gene coding for PF3D7_0201800 to be located in a chromosomal cluster together with knob components KAHRP and PfEMP3. Like the knob components, KAHsp40 too showed the presence of PEXEL motif required for transport to the erythrocyte compartment. Indeed, sub-cellular fractionation and immunofluorescence analysis (IFA) showed KAHsp40 to be exported in the erythrocyte cytoplasm in a stage dependent manner localizing as punctuate spots in the erythrocyte periphery, distinctly from Maurer's cleft, in structures which could be the reminiscent of knobs. Double IFA analysis revealed co-localization of PF3D7_0201800 with the markers of knobs (KAHRP, PfEMP1 and PfEMP3) and components of the PEXEL translocon (Hsp101, PTEX150). KAHsp40 was also found to be in a complex with KAHRP, PfEMP3 and Hsp101 as confirmed by co-immunoprecipitation assay. Our results suggest potential involvement of a parasite encoded Hsp40 in chaperoning knob assembly in the erythrocyte compartment.

An overview of recent advances in structural bioinformatics of protein-protein interactions and a guide to their principles

Rich data bearing on the structural and evolutionary principles of protein protein interactions are paving the way to a better understanding of the regulation of function in the cell. This is particularly the case when these interactions are considered in the framework of key pathways. Knowledge of the interactions may provide insights into the mechanisms of crucial `driver' mutations in oncogenesis. They also provide the foundation toward the design of protein protein interfaces and inhibitors that can abrogate their formation or enhance them. The main features to learn from known 3-D structures of protein protein complexes and the extensive literature which analyzes them computationally and experimentally include the interaction details which permit undertaking structure-based drug discovery, the evolution of complexes and their interactions, the consequences of alterations such as post-translational modifications, ligand binding, disease causing mutations, host pathogen interactions, oligomerization, aggregation and the roles of disorder, dynamics, allostery and more to the protein and the cell. This review highlights some of the recent advances in these areas, including design, inhibition and prediction of protein protein complexes. The field is broad, and much work has been carried out in these areas, making it challenging to cover it in its entirety. Much of this is due to the fast increase in the number of molecules whose structures have been determined experimentally and the vast increase in computational power. Here we provide a concise overview. (C) 2014 Elsevier Ltd. All rights reserved.

Studies on the non-covalent complexes between oleanolic acid and cyclodextrins using electrospray ionization tandem mass spectrometry

Non-covalent inclusion complexes formed between an anti-inflammatory drug, oleanolic acid (OA), and alpha-, beta- and gamma-cyclodextrins (CDs) were investigated by means of solubility studies and electrospray ionization tandem mass spectrometry (ESI-MSn). The order of calculated association constants (K-1:1) of complexes between OA and different CDs in solution is in good agreement with the order of their relative peak intensities and the relative CID energies of the complexes under the same ESI-MSn conditions. These results indicate a direct correlation between the behaviors of solution- and gas-phase complexes. ESI-MS can thus be used to evaluate solution-phase non-covalent complexes successfully. The experimental results show that the most stable 1:1 inclusion complexes between three CDs and OA can be formed, but 2:1 CD-OA complexes can be formed with beta- and gamma-CDs. Multi-component complexes of alpha-CD-OA-beta-CD (1:1:1), alpha-CD-OA-gamma-CD (1:1:1) and beta-CD-OA-gamma-CD (1:1:1) were found in equimolar CD mixtures with excess OA. The formation of 2:1 and multi-component 1:1:1 non-covalent CD-OA complexes indicates that beta- and gamma-CD are able to form sandwich-type inclusion non-covalent complexes with OA. The above results can be partly supported by the relative sizes of OA and CD cavities by molecular modeling calculations.

Consensus scoring for enriching near-native structures from protein-protein docking decoys

The identification of near native protein-protein complexes among a set of decoys remains highly challenging. A stategy for improving the success rate of near native detection is to enrich near native docking decoys in a small number of top ranked decoys. Recently, we found that a combination of three scoring functions (energy, conservation, and interface propensity) can predict the location of binding interface regions with reasonable accuracy. Here, these three scoring functions are modified and combined into a consensus scoring function called ENDES for enriching near native docking decoys. We found that all individual scores result in enrichment for the majority of 28 targets in ZDOCK2.3 decoy set and the 22 targets in Benchmark 2.0. Among the three scores, the interface propensity score yields the highest enrichment in both sets of protein complexes. When these scores are combined into the ENDES consensus score, a significant increase in enrichment of near-native structures is found. For example, when 2000 dock decoys are reduced to 200 decoys by ENDES, the fraction of near-native structures in docking decoys increases by a factor of about six in average. ENDES was implemented into a computer program that is available for download at http://sparks.informatics.iupui.edu.

A binding site for the cyclic adenosine 3',5'-monophosphate-response element-binding protein as a regulatory element in the grp78 promoter.

The 78-kDa glucose-regulated protein (GRP78) is ubiquitously expressed in many cell types. Its promoter contains multiple protein-binding sites and functional elements. In this study we examined a high affinity protein-binding site spanning bp -198 to -180 of the rat grp78 promoter, using nuclear extracts from both B-lymphoid and HeLa cells. This region contains a sequence TGACGTGA which, with the exception of one base, is identical to the cAMP-response element (CRE). Site-directed mutagenesis reveals that this sequence functions as a major basal level regulatory element in hamster fibroblast cells and is also necessary to maintain high promoter activity under stress-induced conditions. By gel mobility shift analysis, we detect two specific protein complexes. The major specific complex I, while immunologically distinct from the 42-kDa CRE-binding protein (CREB), binds most strongly to the grp site, but also exhibits affinity for the CRE consensus sequence. As such, complex I may consist of other members of the CREB/activating transcription factor protein family. The minor specific complex II consists of CREB or a protein antigenically related to it. A nonspecific complex III consists of the Ku autoantigen, an abundant 70- to 80-kDa protein complex in HeLa nuclear extracts. By cotransfection experiments, we demonstrate that in F9 teratocarcinoma cells, the grp78 promoter can be transactivated by the phosphorylated CREB or when the CREB-transfected cells are treated with the calcium ionophore A23187. The differential regulation of the grp78 gene by cAMP in specific cell types and tissues is discussed.

Analysis of IL6 Signal Transduction Model using Reduced Rank Regression

This work presents the application of reduced rank regression to the field of systems biology. A computational approach is used to investigate the mechanisms of the janus-associated kinases/signal transducers and transcription factors (JAK/STAT) and mitogen activated protein kinases (MAPK) signal transduction pathways in hepatic cells stimulated by interleukin-6. The results obtained identify the contribution of individual reactions to the dynamics of the model. These findings are compared to previously available results from sensitivity analysis of the model which focused on the parameters involved and their effect. This application of reduced rank regression allows for an understanding of the individual reaction terms involved in the modelled signal transduction pathways and has the benefit of being computationally inexpensive. The obtained results complement existing findings and also confirm the importance of several protein complexes in the MAPK pathway which hints at benefits that can be achieved by further refining the model.

BRCA1, a 'complex' protein involved in the maintenance of genomic stability

BRCA1 is a major breast and ovarian cancer susceptibility gene, with mutations in this gene predisposing women to a very high risk of developing breast and ovarian tumours. BRCA1 primarily functions to maintain genomic stability via critical roles in DNA repair, cell cycle checkpoint control, transcriptional regulation, apoptosis and mRNA splicing. As a result, BRCA1 mutations often result in defective DNA repair, genomic instability and sensitivity to DNA damaging agents. BRCA1 carries out these different functions through its ability to interact, and form complexes with, a vast array of proteins involved in multiple cellular processes, all of which are considered to contribute to its function as a tumour suppressor. This review discusses and highlights recent research into the functions of BRCA1-related protein complexes and their roles in maintaining genomic stability and tumour suppression.

Real-time monitoring of protein conformational changes using a nano-mechanical sensor.

Proteins can switch between different conformations in response to stimuli, such as pH or temperature variations, or to the binding of ligands. Such plasticity and its kinetics can have a crucial functional role, and their characterization has taken center stage in protein research. As an example, Topoisomerases are particularly interesting enzymes capable of managing tangled and supercoiled double-stranded DNA, thus facilitating many physiological processes. In this work, we describe the use of a cantilever-based nanomotion sensor to characterize the dynamics of human topoisomerase II (Topo II) enzymes and their response to different kinds of ligands, such as ATP, which enhance the conformational dynamics. The sensitivity and time resolution of this sensor allow determining quantitatively the correlation between the ATP concentration and the rate of Topo II conformational changes. Furthermore, we show how to rationalize the experimental results in a comprehensive model that takes into account both the physics of the cantilever and the dynamics of the ATPase cycle of the enzyme, shedding light on the kinetics of the process. Finally, we study the effect of aclarubicin, an anticancer drug, demonstrating that it affects directly the Topo II molecule inhibiting its conformational changes. These results pave the way to a new way of studying the intrinsic dynamics of proteins and of protein complexes allowing new applications ranging from fundamental proteomics to drug discovery and development and possibly to clinical practice.

Cartographie du réseau d'interactions protéiques de la machinerie de transcription dans les cellules humaines

Les protéines sont les macromolécules les plus polyvalentes de la cellule. Elles jouent un rôle fondamental dans la majorité des processus biologiques à travers la formation de complexes multi-protéiques. Durant la transcription, une multitude de facteurs sont impliquées dans le contrôle de l’activité des complexes ARN polymérases. Notre laboratoire s’est intéressé au réseau d’interaction de la machinerie de transcription des ARN polymérases nucléaires, dans le but de mieux comprendre leurs mécanismes de régulation. Pour ce faire, une procédure protéomique comprenant la purification de complexes protéiques par affinité couplée à la spectrométrie de masse et à l’analyse bioinformatique a été développée. La méthode de purification TAP (Tandem Affinity Purification) a été adaptée pour permettre la purification de complexes protéiques solubles assemblés in vivo à partir de cellules humaines. L’objectif de mon projet de maîtrise était de purifier le complexe de l’ARN Pol I ainsi que de poursuivre l’expansion du réseau d’interactions protéine-protéine de la machinerie de transcription de l’ARN Pol II humaine. À l’aide des protéines POLR1E, TWISTNB, POLR2E, PFDN4, MBD2, XPA, CAND1 et PDCD5 étiquetées (TAP-tag) exprimées dans des lignées cellulaires ECR-293, plusieurs complexes protéiques solubles ont été purifiés et analysés par spectrométrie de masse. Les interactions protéiques ont été triées et validées bioinformatiquement pour donner en final une liste d’interactions ayant un haut degré de confiance à partir de laquelle des réseaux d’interactions protéine-protéine ont été créés. Le réseau créé au cours de ce projet connecte plusieurs composantes de la machinerie transcriptionnelle tels que les ARN Pol I, II et III, les complexes RPAP3/R2TP/prefoldin-like, TRiC/CCT, Mi-2/NuRD et des facteurs de transcription et de réparation de l’ADN. Ce type d’analyse nous a permis d’identifier et de caractériser de nouveaux régulateurs de la machinerie de transcription de l’ARN Pol I et II et de mieux comprendre son fonctionnement.