Stoichiometry of chromatin-associated protein complexes revealed by label-free quantitative mass spectrometry-based proteomics

Many cellular proteins assemble into macromolecular protein complexes. The identification of protein–protein interactions and quantification of their stoichiometry is therefore crucial to understand the molecular function of protein complexes. Determining the stoichiometry of protein complexes is usually achieved by mass spectrometry-based methods that rely on introducing stable isotope-labeled reference peptides into the sample of interest. However, these approaches are laborious and not suitable for high-throughput screenings. Here, we describe a robust and easy to implement label-free relative quantification approach that combines the detection of high-confidence protein–protein interactions with an accurate determination of the stoichiometry of the identified protein–protein interactions in a single experiment. We applied this method to two chromatin-associated protein complexes for which the stoichiometry thus far remained elusive: the MBD3/NuRD and PRC2 complex. For each of these complexes, we accurately determined the stoichiometry of the core subunits while at the same time identifying novel interactors and their stoichiometry.     

Coupling dynamics and evolutionary information with structure to identify protein regulatory and functional binding sites

Binding sites in proteins can be either specifically functional binding sites (active sites) that bind specific substrates with high affinity or regulatory binding sites (allosteric sites), that modulate the activity of functional binding sites through effector molecules. Owing to their significance in determining protein function, the identification of protein functional and regulatory binding sites is widely acknowledged as an important biological problem. In this work, we present a novel binding site prediction method, Active and Regulatory site Prediction (AR-Pred), which supplements protein geometry, evolutionary, and physicochemical features with information about protein dynamics to predict putative active and allosteric site residues. As the intrinsic dynamics of globular proteins plays an essential role in controlling binding events, we find it to be an important feature for the identification of protein binding sites. We train and validate our predictive models on multiple balanced training and validation sets with random forest machine learning and obtain an ensemble of discrete models for each prediction type. Our models for active site prediction yield a median area under the curve (AUC) of 91% and Matthews correlation coefficient (MCC) of 0.68, whereas the less well-defined allosteric sites are predicted at a lower level with a median AUC of 80% and MCC of 0.48. When tested on an independent set of proteins, our models for active site prediction show comparable performance to two existing methods and gains compared to two others, while the allosteric site models show gains when tested against three existing prediction methods. AR-Pred is available as a free downloadable package at https://github.com/sambitmishra0628/AR-PRED_source .     

Thermodynamic and electrostatic properties of ternary Oxytricha nova TEBP-DNA complex

Telomeres constitute the nucleoprotein ends of eukaryotic chromosomes which are essential for their proper function. Telomere end binding protein (TEBP) from Oxytricha nova was among the first telomeric proteins, which were well characterized biologically. TEBP consists of two protein subunits (alpha, beta) and forms a ternary complex with single stranded telomeric DNA containing tandem repeats TTTTGGGG. This work presents the characterization of the thermodynamic and electrostatic properties of this complex by computational chemistry methods (continuum Poisson-Boltzmann and solvent accessible surface calculations). Our calculations give a new insight into molecular properties of studied system. Based on the thermodynamic analysis we provide a rationale for the experimental observation that alpha and ssDNA forms a binary complex and the beta subunit joins alpha:ssDNA complex only after the latter is formed. Calculations of distribution of the molecular electrostatic potential for protein subunits alone and for all possible binary complexes revealed the important role of the "guiding funnel" potential generated by alpha:ssDNA complex. This potential may help the beta subunit to dock to the already formed alpha:DNA intermediate in highly steric and electrostatic favorable manner. Our pK(a) calculations of TEBP are able to explain the experimental mobility shifts of the complex in electrophoretic non-denaturating gels.     

Localization of the high affinity calcium-binding site on tubulin molecule

Tubulin is a calcium-binding protein. Two different modes of interaction of calcium with tubulin have been described: a high affinity interaction to one or two binding sites and lower affinity interactions to several other binding sites. In the present study, we have used limited proteolysis of tubulin with trypsin, chymotrypsin, and subtilisin to localize the high affinity calcium-binding sites. Our results indicate that two sites are located in the carboxyl-terminal region of both tubulin subunits, and that tubulin deprived of its carboxyl-terminal region is able to polymerize in the presence of 0.5 mM calcium.     

Structural characterization of the human eukaryotic initiation factor 3 protein complex by mass spectrometry

Protein synthesis in mammalian cells requires initiation factor eIF3, an approximately 800-kDa protein complex that plays a central role in binding of initiator methionyl-tRNA and mRNA to the 40 S ribosomal subunit to form the 48 S initiation complex. The eIF3 complex also prevents premature association of the 40 and 60 S ribosomal subunits and interacts with other initiation factors involved in start codon selection. The molecular mechanisms by which eIF3 exerts these functions are poorly understood. Since its initial characterization in the 1970s, the exact size, composition, and post-translational modifications of mammalian eIF3 have not been rigorously determined. Two powerful mass spectrometric approaches were used in the present study to determine post-translational modifications that may regulate the activity of eIF3 during the translation initiation process and to characterize the molecular structure of the human eIF3 protein complex purified from HeLa cells. In the first approach, the bottom-up analysis of eIF3 allowed for the identification of a total of 13 protein components (eIF3a-m) with a sequence coverage of approximately 79%. Furthermore 29 phosphorylation sites and several other post-translational modifications were unambiguously identified within the eIF3 complex. The second mass spectrometric approach, involving analysis of intact eIF3, allowed the detection of a complex with each of the 13 subunits present in stoichiometric amounts. Using tandem mass spectrometry four eIF3 subunits (h, i, k, and m) were found to be most easily dissociated and therefore likely to be on the periphery of the complex. It is noteworthy that none of these four subunits were found to be phosphorylated. These data raise interesting questions about the function of phosphorylation as it relates to the core subunits of the complex.     

Identification and characterization of fragment binding sites for allosteric ligand design using the site identification by ligand competitive saturation hotspots approach (SILCS-Hotspots)

BackgroundFragment-based ligand design is used for the development of novel ligands that target macromolecules, most notably proteins. Central to its success is the identification of fragment binding sites that are spatially adjacent such that fragments occupying those sites may be linked to create drug-like ligands. Current experimental and computational approaches that address this problem typically identify only a limited number of sites as well as use a limited number of fragment types.MethodsThe site-identification by ligand competitive saturation (SILCS) approach is extended to the identification of fragment bindings sites, with the method termed SILCS-Hotspots. The approach involves precomputation of the SILCS FragMaps following which the identification of Hotspots, performed by identifying of all possible fragment binding sites on the full 3D structure of the protein followed by spatial clustering.ResultsThe SILCS-Hotspots approach identifies a large number of sites on the target protein, including many sites not accessible in experimental structures due to low binding affinities and binding sites on the protein interior. The identified sites are shown to recapitulate the location of known drug-like molecules in both allosteric and orthosteric binding sites on seven proteins including the androgen receptor, the CDK2 and Erk5 kinases, PTP1B phosphatase and three GPCRs; the β2-adrenergic, GPR40 fatty-acid binding and M2-muscarinic receptors. Analysis indicates the importance of considering all possible fragment binding sites, and not just those accessible to experimental methods, when identifying novel binding sites and performing ligand design versus just considering the most favorable sites. The approach is shown to identify a larger number of known binding sites of drug-like molecules versus the commonly used FTMap and Fpocket methods.General significanceThe present results indicate the potential utility of the SILCS-Hotspots approach for fragment-based rational design of ligands, including allosteric modulators.     

The 90S preribosome is a multimodular structure that is assembled through a hierarchical mechanism

The 90S preribosomal particle is required for the production of the 18S rRNA from a pre-rRNA precursor. Despite the identification of the protein components of this particle, its mechanism of assembly and structural design remain unknown. In this work, we have combined biochemical studies, proteomic techniques, and bioinformatic analyses to shed light into the rules of assembly of the yeast 90S preribosome. Our results indicate that several protein subcomplexes work as discrete assembly subunits that bind in defined steps to the 35S pre-rRNA. The assembly of the t-UTP subunit is an essential step for the engagement of at least five additional subunits in two separate, and mutually independent, assembling routes. One of these routes leads to the formation of an assembly intermediate composed of the U3 snoRNP, the Pwp2p/UTP-B, subunit and the Mpp10p complex. The other assembly route involves the stepwise binding of Rrp5p and the UTP-C subunit. We also report the use of a bioinformatic approach that provides a model for the topological arrangement of protein components within the fully assembled particle. Together, our data identify the mechanism of assembly of the 90S preribosome and offer novel information about its internal architecture.     

Genome-wide identification of protein binding sites on RNAs in mammalian cells

RNA-binding proteins (RBPs) are proteins that bind to the RNA and participate in forming ribonucleoprotein complexes. They have crucial roles in various biological processes such as RNA splicing, editing, transport, maintenance, degradation, intracellular localization and translation. The RBPs bind RNA with different RNA-sequence specificities and affinities, thus, identification of protein binding sites on RNAs (R-PBSs) will deeper our understanding of RNA-protein interactions. Currently, high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP, also known as CLIP-Seq) is one of the most powerful methods to map RNA-protein binding sites or RNA modification sites. However, this method is only used for identification of single known RBPs and antibodies for RBPs are required. Here we developed a novel method, called capture of protein binding sites on RNAs (RPBS-Cap) to identify genome-wide protein binding sites on RNAs without using antibodies. Double click strategy is used for the RPBS-Cap assay. Proteins and RNAs are UV-crosslinked in vivo first, then the proteins are crosslinked to the magnetic beads. The RNA elements associated with proteins are captured, reverse transcribed and sequenced. Our approach has potential applications for studying genome-wide RNA-protein interactions.     

Crystal structure of human vacuolar protein sorting protein 29 reveals a phosphodiesterase/nuclease-like fold and two protein-protein interaction sites

Vacuolar protein sorting protein 29 (Vps29p), which is involved in retrograde trafficking from prevacuolar endosomes to the trans-Golgi network, performs its biological functions by participating in the formation of a "retromer complex." In human cells, this complex comprises four conserved proteins: hVps35p, hVps29p, hVps26p, and sorting nexin 1 protein (SNX1). Here, we report the crystal structure of hVps29p at 2.1 Angstroms resolution, the first three-dimensional structure of the retromer subunits. This novel structure adopts a four-layered alpha-beta-beta-alpha sandwich fold. hVps29p contains a metal-binding site that is very similar to the active sites of some proteins of the phosphodiesterase/nuclease protein family, indicating that hVps29p may carry out chemically similar functions. Structure and sequence conservation analysis suggests that hVps29p contains two protein-protein interaction sites. One site, which potentially serves as the interface between hVps29p and hVps35p, comprises 5 conserved hydrophobic and 8 hydrophilic residues. The other site is relatively more hydrophilic and may serve as a binding interface with hVps26p, SNX1, or other target proteins.     

A switch in the mechanism of communication between the two DNA-binding sites in the SfiI restriction endonuclease

While many Type II restriction enzymes are dimers with a single DNA-binding cleft between the subunits, SfiI is a tetramer of identical subunits. Two of its subunits (a dimeric unit) create one DNA-binding cleft, and the other two create a second cleft on the opposite side of the protein. The two clefts bind specific DNA cooperatively to give a complex of SfiI with two recognition sites. This complex is responsible for essentially all of the DNA-cleavage reactions by SfiI: virtually none is due to the complex with one site. The communication between the DNA-binding clefts was examined by disrupting one of the very few polar interactions in the otherwise hydrophobic interface between the dimeric units: a tyrosine hydroxyl was removed by mutation to phenylalanine. The mutant protein remained tetrameric in solution and could bind two DNA sites. But instead of being activated by binding two sites, like wild-type SfiI, it showed maximal activity when bound to a single site and had a lower activity when bound to two sites. This interaction across the dimer interface thus enforces in wild-type SfiI a cooperative transition between inactive and active states in both dimers, but without this interaction as in the mutant protein, a single dimer can undergo the transition to give a stable intermediate with one inactive dimer and one active dimer.     

"Half of the sites" binding of D-glyceraldehyde-3-phosphate dehydrogenase folding intermediate with GroEL

Two D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) folding intermediate subunits bind with chaperonin 60 (GroEL) to form a stable complex, which can no longer bind with additional GAPDH intermediate subunits, but does bind with one more lysozyme folding intermediate or one chaperonin 10 (GroES) molecule, suggesting that the two GAPDH subunits bind at one end of the GroEL molecule displaying a "half of the sites" binding profile. For lysozyme, GroEL binds with either one or two folding intermediates to form a stable 1:1 or 1:2 complex with one substrate on each end of the GroEL double ring for the latter. The 1:1 complex of GroEL.GroES binds with one lysozyme or one dimeric GAPDH folding intermediate to form a stable ternary complex. Both complexes of GroEL.lysozyme1 and GroEL.GAPDH2 bind with one GroES molecule only at the other end of the GroEL molecule forming a trans ternary complex. According to the stoichiometry of GroEL binding with the GAPDH folding intermediate and the formation of ternary complexes containing GroEL.GAPDH2, it is suggested that the folding intermediate of GAPDH binds, very likely in the dimeric form, with GroEL at one end only.     

Mass spectrometric characterization of the affinity-purified human 26S proteasome complex

The 26S proteasome is a multisubunit complex responsible for degradation of ubiquitinated substrates, which plays a critical role in regulating various biological processes. To fully understand the function and regulation of the proteasome complex, an important step is to elucidate its subunit composition and posttranslational modifications. Toward this goal, a new affinity purification strategy has been developed using a derivative of the HB tag for rapid isolation of the human 26S proteasome complex for subsequent proteomic analysis. The purification of the complex is achieved from stable 293 cell lines expressing a HB-tagged proteasome subunit and by high-affinity streptavidin binding with TEV cleavage elution. The complete composition of the 26S proteasome complex, including recently assigned new subunits, is identified by LC-MS/MS. In addition, all known proteasome activator proteins and components involved in the ubiquitin-proteasome degradation pathway are identified. Aside from the subunit composition, the N-terminal modification and phosphorylation of the proteasome subunits have been characterized. Twelve novel phosphorylation sites from eight subunits have been identified, and N-terminal modifications are determined for 25 subunits, 12 of which have not been previously reported in mammals. We also observe different N-terminal processing of subunit Rpn2, which results in identification of two different N-termini of the protein. This work presents the first comprehensive characterization of the human 26S proteasome complex by affinity purification and tandem mass spectrometry. The detailed proteomic profiling obtained here is significant to future studies aiming at a complete understanding of the structure-function relationship of the human 26S proteasome complex.     

Subunit composition of Rhodothermus marinus respiratory complex I

The basic structural characterization of complex I is still not trivial due to its complexity, not only in the number of its protein constituents but especially because of the different properties of the several subunits. Bacterial complex I generally contains 14 subunits: 7 hydrophilic proteins located in the peripheral arm and 7 hydrophobic proteins present in the membrane arm. It is the identification of the hydrophobic proteins that makes the characterization of complex I, and of membrane proteins in general, very difficult. In this article, we report the identification of the subunits of complex I from Rhodothermus marinus. The original approach, presented here, combined several protein and peptides separation strategies (different reversed phase materials, high-performance liquid chromatography, and gel electrophoresis) with different identification methods (electrospray ionization-tandem mass spectrometry, matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry, and Edman degradation analysis) and represents a step forward in the characterization of membrane proteins that studies are still technically highly challenging. The combination of the different methodologies allowed the identification of complex I canonical subunits and also a possible novel subunit, namely a pterin-4α-carbinolamine dehydratase (PCD). This was the first time that a PCD was suggested to be part of complex I, and its possible regulatory role is discussed.