SimShift: identifying structural similarities from NMR chemical shifts

MOTIVATION: An important quantity that arises in NMR spectroscopy experiments is the chemical shift. The interpretation of these data is mostly done by human experts; to our knowledge there are no algorithms that predict protein structure from chemical shift sequences alone. One approach to facilitate this process could be to compare two such sequences, where the structure of one protein has already been resolved. Our claim is that similarity of chemical shifts thereby found implies structural similarity of the respective proteins. RESULTS: We present an algorithm to identify structural similarities of proteins by aligning their associated chemical shift sequences. To evaluate the correctness of our predictions, we propose a benchmark set of protein pairs that have high structural similarity, but low sequence similarity (because with high sequence similarity the structural similarities could easily be detected by a sequence alignment algorithm). We compare our results with those of HHsearch and SSEA and show that our method outperforms both in >50% of all cases.     

Computational identification of proteins for selectivity assays

At the stage of optimization of a chemical series the compounds are normally assayed for binding or inhibition on the target protein as well as on several proteins from a selectivity panel. These proteins are normally identified on the basis of sequence homology to the target protein. Experimental selectivity data are also taken into account if available. Cases when a nonhomologous protein has a significant affinity to the compound series are going to be missed if the selectivity panel is identified by homology. Experimental data is usually either unavailable or limited to a small fraction of proteins that should be considered. We have developed a computational method of identification of selectivity panel proteins. It is based on the evaluation of binding site similarity to the target protein using docking scores of target-selected molecular probes. These probes are obtained by docking a large library of drug-like compounds to the target protein followed by selecting a diverse subset from the best virtual binders. Docking scores of these probes to other proteins measure binding site similarity to the target. Because the method does not require prior knowledge of either affinities or structures of inhibitors for the target, it can be applied to any protein with known 3D structure. Validation of the method includes rediscovery of nonhomologous proteins that bind common ligands: estradiol, tamoxifen, and riboflavin. Given 3D structures, the method can effectively discriminate proteins with similar binding sites from random proteins independent of sequence homology.     

A combination method of chemical with enzyme reactions for identification of membrane proteins

A simple method for effective analysis of various proteins has been developed, including membrane proteins, with LC-MS/MS, using CNBr and acetic acid cleavage in one reaction for the digestion of both the M/ and /D/ positions within the target proteins. This dual chemical reaction has been compared with traditional CNBr or an acid cleavage method using a rat kidney membrane fraction and it showed an advantage of the dual reaction with respect to a high number of peptides detected and a high protein recovery. Furthermore, when this dual chemical reaction was combined with trypsin digestion, the number of proteins surprisingly increased approximately 3.0 times more than in the cases with the trypsin digestion only. It was also 1.9 times more than in cases dealing with Tube-Gel trypsin digestion, which is one of the most efficient digestion methods. In addition, it was shown that this dual chemical reaction could be applied to an in-gel digestion. Using the combination of the chemical and enzyme reaction, 172 proteins including 95 membrane proteins were identified. This indicated that this method is one of the efficient systems in single MS/MS analysis. In particular, many membrane proteins identified in this study were detected by a new combination, but not by a traditional trypsin digestion method.     

Validation of archived chemical shifts through atomic coordinates

The public archives containing protein information in the form of NMR chemical shift data at the BioMagResBank (BMRB) and of 3D structure coordinates at the Protein Data Bank are continuously expanding. The quality of the data contained in these archives, however, varies. The main issue for chemical shift values is that they are determined relative to a reference frequency. When this reference frequency is set incorrectly, all related chemical shift values are systematically offset. Such wrongly referenced chemical shift values, as well as other problems such as chemical shift values that are assigned to the wrong atom, are not easily distinguished from correct values and effectively reduce the usefulness of the archive. We describe a new method to correct and validate protein chemical shift values in relation to their 3D structure coordinates. This method classifies atoms using two parameters: the per‐atom solvent accessible surface area (as calculated from the coordinates) and the secondary structure of the parent amino acid. Through the use of Gaussian statistics based on a large database of 3220 BMRB entries, we obtain per‐entry chemical shift corrections as well as Z scores for the individual chemical shift values. In addition, information on the error of the correction value itself is available, and the method can retain only dependable correction values. We provide an online resource with chemical shift, atom exposure, and secondary structure information for all relevant BMRB entries ( http://www.ebi.ac.uk/pdbe/nmr/vasco ) and hope this data will aid the development of new chemical shift‐based methods in NMR. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Isotopically labeled crosslinking reagents: resolution of mass degeneracy in the identification of crosslinked peptides

Mass spectrometry in three dimensions (MS3D) is a newly developed method for the determination of protein structures involving intramolecular chemical crosslinking of proteins, proteolytic digestion of the resulting adducts, identification of crosslinks by mass spectrometry (MS), peak assignment using theoretical mass lists, and computational reduction of crosslinks to a structure by distance geometry methods. To facilitate the unambiguous identification of crosslinked peptides from proteolytic digestion mixtures of crosslinked proteins by MS, we introduced double 18O isotopic labels into the crosslinking reagent to provide the crosslinked peptides with a characteristic isotope pattern. The presence of doublets separated by 4 Da in the mass spectra of these materials allowed ready discrimination between crosslinked and modified peptides, and uncrosslinked peptides using automated intelligent data acquisition (IDA) of MS/MS data. This should allow ready automation of the method for application to whole expressible proteomes.     

First identification of tannin-binding proteins in saliva of Papio hamadryas using MS/MS mass spectrometry

Hamadryas baboons possess salivary proline-rich proteins (PRP), as indicated by the presence of pink-staining protein bands using 1D SDS gel electrophoresis and Coomassie R250 staining. The ability of these protein bands to interact with tannic acid was further examined. In a tannin-binding assay using 5 μg tannic acid mixed with hamadryas whole saliva, we recently found four distinct protein bands of apparently 72, 55, 20, and 15 kDa that were precipitated during the experiments. In this work, we were able to identify these protein bands in a follow-up analysis using MS/MS mass spectrometry after excising such bands out of air-dried gels. Albumin and α-amylase were present in the tannic acid-protein complexes, with albumin already known to nonspecifically interact with a great diversity of chemical compounds. More interesting, we also identified a basic PRP and a cystatin precursor protein. This was the first successful attempt to identify a PRP from precipitated tannin-protein complexes in hamadryas baboons using MS/MS mass spectrometry. On the other hand, the role of cystatins in tannin binding is not yet well understood. However, there are recent reports on cystatin expression in saliva of rats responding to astringent dietary compounds. In conclusion, the follow-up data on tannin-binding proteins present in salivary secretions from hamadryas baboons adds important knowledge to primate physiology and feeding ecology, in order to shed light on the establishment and development of food adaptations in primates. It also demonstrates that tannin binding is characteristic for PRP, but might not be restricted to this particular group of proteins in primate species.     

Semi-automatic high-throughput determination of plasma protein binding using a 96-well plate filtrate assembly and fast liquid chromatography-tandem mass spectrometry

A semi-automatic, high-throughput method has been developed to rapidly assess plasma protein binding of new chemical entities in drug discovery phase. New chemical entities are mixed with plasma and the unbound fractions are separated from the bound fraction by ultrafiltration in a 96-well filtrate assembly. The unbound fractions are then analyzed by fast liquid chromatography-tandem mass spectrometry (LC-MS/MS). Sample handling is automated by a robotic system. Employing a cocktail approach where multiple new chemical entities are allowed to bind to plasma proteins in the same well has further increased the throughput. We have validated the method with 12 commercially available compounds. The plasma protein binding data obtained by this method are comparable with the literature values. This method enables the determination of protein binding for 32 compounds in one single experiment instead of 1-2 compounds using the conventional methods.     

Rapid and reliable protein structure determination via chemical shift threading

Protein structure determination using nuclear magnetic resonance (NMR) spectroscopy can be both time-consuming and labor intensive. Here we demonstrate how chemical shift threading can permit rapid, robust, and accurate protein structure determination using only chemical shift data. Threading is a relatively old bioinformatics technique that uses a combination of sequence information and predicted (or experimentally acquired) low-resolution structural data to generate high-resolution 3D protein structures. The key motivations behind using NMR chemical shifts for protein threading lie in the fact that they are easy to measure, they are available prior to 3D structure determination, and they contain vital structural information. The method we have developed uses not only sequence and chemical shift similarity but also chemical shift-derived secondary structure, shift-derived super-secondary structure, and shift-derived accessible surface area to generate a high quality protein structure regardless of the sequence similarity (or lack thereof) to a known structure already in the PDB. The method (called E-Thrifty) was found to be very fast (often?<?10 min/structure) and to significantly outperform other shift-based or threading-based structure determination methods (in terms of top template model accuracy)—with an average TM-score performance of 0.68 (vs. 0.50–0.62 for other methods). Coupled with recent developments in chemical shift refinement, these results suggest that protein structure determination, using only NMR chemical shifts, is becoming increasingly practical and reliable. E-Thrifty is available as a web server at http://ethrifty.ca.     

Accurate ab initio prediction of NMR chemical shifts of nucleic acids and nucleic acids/protein complexes

NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients (r²) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.     

A novel peak detection approach with chemical noise removal using short‐time FFT for prOTOF MS data

Peak detection is a pivotal first step in biomarker discovery from MS data and can significantly influence the results of downstream data analysis steps. We developed a novel automatic peak detection method for prOTOF MS data, which does not require a priori knowledge of protein masses. Random noise is removed by an undecimated wavelet transform and chemical noise is attenuated by an adaptive short‐time discrete Fourier transform. Isotopic peaks corresponding to a single protein are combined by extracting an envelope over them. Depending on the S/N, the desired peaks in each individual spectrum are detected and those with the highest intensity among their peak clusters are recorded. The common peaks among all the spectra are identified by choosing an appropriate cut‐off threshold in the complete linkage hierarchical clustering. To remove the 1 Da shifting of the peaks, the peak corresponding to the same protein is determined as the detected peak with the largest number among its neighborhood. We validated this method using a data set of serial peptide and protein calibration standards. Compared with MoverZ program, our new method detects more peaks and significantly enhances S/N of the peak after the chemical noise removal. We then successfully applied this method to a data set from prOTOF MS spectra of albumin and albumin‐bound proteins from serum samples of 59 patients with carotid artery disease compared to vascular disease‐free patients to detect peaks with S/N≥2. Our method is easily implemented and is highly effective to define peaks that will be used for disease classification or to highlight potential biomarkers.     

Generation of 3D templates of active sites of proteins with rigid prosthetic groups

MOTIVATION: With the increasing availability of protein structures, the generation of biologically meaningful 3D patterns from the simultaneous alignment of several protein structures is an exciting prospect: active sites could be better understood, protein functions and protein 3D structures could be predicted more accurately. Although patterns can already be generated at the fold and topological levels, no system produces high-resolution 3D patterns including atom and cavity positions. To address this challenge, our research focuses on generating patterns from proteins with rigid prosthetic groups. Since these groups are key elements of protein active sites, the generated 3D patterns are expected to be biologically meaningful. RESULTS: In this paper, we present a new approach which allows the generation of 3D patterns from proteins with rigid prosthetic groups. Using 237 protein chains representing proteins containing porphyrin rings, our method was validated by comparing 3D templates generated from homologues with the 3D structure of the proteins they model. Atom positions were predicted reliably: 93% of them had an accuracy of 1.00 A or less. Moreover, similar results were obtained regarding chemical group and cavity positions. Results also suggested our system could contribute to the validation of 3D protein models. Finally, a 3D template was generated for the active site of human cytochrome P450 CYP17, the 3D structure of which is unknown. Its analysis showed that it is biologically meaningful: our method detected the main patterns of the cytochrome P450 superfamily and the motifs linked to catalytic reactions. The 3D template also suggested the position of a residue, which could be involved in a hydrogen bond with CYP17 substrates and the shape and location of a cavity. Comparisons with independently generated 3D models comforted these hypotheses. AVAILABILITY: Alignment software (Nestor3D) is available at http://www.kingston.ac.uk/~ku33185/Nestor3D.html     

Four-dimensional heteronuclear correlation experiments for chemical shift assignment of solid proteins

Chemical shift assignment is the first step in all established protocols for structure determination of uniformly labeled proteins by NMR. The explosive growth in recent years of magic-angle spinning (MAS) solid-state NMR (SSNMR) applications is largely attributable to improved methods for backbone and side-chain chemical shift correlation spectroscopy. However, the techniques developed so far have been applied primarily to proteins in the size range of 5–10 kDa, despite the fact that SSNMR has no inherent molecular weight limits. Rather, the degeneracy inherent to many 2D and 3D SSNMR spectra of larger proteins has prevented complete unambiguous chemical shift assignment. Here we demonstrate the implementation of 4D backbone chemical shift correlation experiments for assignment of solid proteins. The experiments greatly reduce spectral degeneracy at a modest cost in sensitivity, which is accurately described by theory. We consider several possible implementations and investigate the CANCOCX pulse sequence in detail. This experiment involves three cross polarization steps, from H to CA[i], CA[i] to N[i], and N[i] to C′[i−1], followed by a final homonuclear mixing period. With short homonuclear mixing times (<20 ms), backbone correlations are observed with high sensitivity; with longer mixing times (>200 ms), long-range correlations are revealed. For example, a single 4D experiment with 225 ms homonuclear mixing time reveals ∼200 uniquely resolved medium and long-range correlations in the 56-residue protein GB1. In addition to experimental demonstrations in the 56-residue protein GB1, we present a theoretical analysis of anticipated improvements in resolution for much larger proteins and compare these results in detail with the experiments, finding good agreement between experiment and theory under conditions of stable instrumental performance.     

A chemical method for investigating disulfide-coupled peptide and protein folding

Investigations of protein folding have largely involved studies using disulfide-containing proteins, as disulfide-coupled folding of proteins permits the folding intermediates to be trapped and their conformations determined. Over the last decade, a combination of new biotechnical and chemical methodology has resulted in a remarkable acceleration in our understanding of the mechanism of disulfide-coupled protein folding. In particular, expressed protein ligation, a combination of native chemical ligation and an intein-based approach, permits specifically labeled proteins to be easily produced for studies of protein folding using biophysical methods, such as NMR spectroscopy and X-ray crystallography. A method for regio-selective formation of disulfide bonds using chemical procedures has also been established. This strategy is particularly relevant for the study of disulfide-coupled protein folding, and provides us not only with the native conformation, but also the kinetically trapped topological isomer with native disulfide bonds. Here we review recent developments and applications of biotechnical and chemical methods to investigations of disulfide-coupled peptide and protein folding. Chemical additives designed to accelerate correct protein folding and to avoid non-specific aggregation are also discussed.     

Liquid chromatography MALDI MS/MS for membrane proteome analysis

Membrane proteins play critical roles in many biological functions and are often the molecular targets for drug discovery. However, their analysis presents a special challenge largely due to their highly hydrophobic nature. We present a surfactant-aided shotgun proteomics approach for membrane proteome analysis. In this approach, membrane proteins were solubilized and digested in the presence of SDS followed by newly developed auto-offline liquid chromatography/matrix-assisted laser desorption ionization (LC/MALDI) tandem MS analysis. Because of high tolerance of MALDI to SDS, one-dimensional (1D) LC separation can be combined with MALDI for direct analysis of protein digests containing SDS, without the need for extensive sample cleanup. In addition, the heated droplet interface used in LC/MALDI can work with high flow LC separations, allowing a relatively large amount of protein digest to be used for 1D LC/MALDI which facilitates the detection of low abundance proteins. The proteome identification results obtained by LC/MALDI are compared to the gel electrophoresis/MS method as well as the shotgun proteomics method using 2D LC/electrospray ionization MS. It is demonstrated that, while LC/MALDI provides more extensive proteome coverage compared to the other two methods, these three methods are complementary to each other and a combination of these methods should provide a more comprehensive membrane proteome analysis.     

Two-dimensional mass spectra generated from the analysis of 15N-labeled and unlabeled peptides for efficient protein identification and de novo peptide sequencing

Protein identification has been greatly facilitated by database searches against protein sequences derived from product ion spectra of peptides. This approach is primarily based on the use of fragment ion mass information contained in a MS/MS spectrum. Unambiguous protein identification from a spectrum with low sequence coverage or poor spectral quality can be a major challenge. We present a two-dimensional (2D) mass spectrometric method in which the numbers of nitrogen atoms in the molecular ion and the fragment ions are used to provide additional discriminating power for much improved protein identification and de novo peptide sequencing. The nitrogen number is determined by analyzing the mass difference of corresponding peak pairs in overlaid spectra of (15)N-labeled and unlabeled peptides. These peptides are produced by enzymatic or chemical cleavage of proteins from cells grown in (15)N-enriched and normal media, respectively. It is demonstrated that, using 2D information, i.e., m/z and its associated nitrogen number, this method can, not only confirm protein identification results generated by MS/MS database searching, but also identify peptides that are not possible to identify by database searching alone. Examples are presented of analyzing Escherichia coli K12 extracts that yielded relatively poor MS/MS spectra, presumably from the digests of low abundance proteins, which can still give positive protein identification using this method. Additionally, this 2D MS method can facilitate spectral interpretation for de novo peptide sequencing and identification of posttranslational or other chemical modifications. We envision that this method should be particularly useful for proteome expression profiling of organelles or cells that can be grown in (15)N-enriched media.     

Review physical and chemical aspects of stabilization of compounds in silk

The challenge of stabilization of small molecules and proteins has received considerable interest. The biological activity of small molecules can be lost as a consequence of chemical modifications, while protein activity may be lost due to chemical or structural degradation, such as a change in macromolecular conformation or aggregation. In these cases, stabilization is required to preserve therapeutic and bioactivity efficacy and safety. In addition to use in therapeutic applications, strategies to stabilize small molecules and proteins also have applications in industrial processes, diagnostics, and consumer products like food and cosmetics. Traditionally, therapeutic drug formulation efforts have focused on maintaining stability during product preparation and storage. However, with growing interest in the fields of encapsulation, tissue engineering, and controlled release drug delivery systems, new stabilization challenges are being addressed; the compounds or protein of interest must be stabilized during: (1) fabrication of the protein or small molecule-loaded carrier, (2) device storage, and (3) for the duration of intended release needs in vitro or in vivo. We review common mechanisms of compound degradation for small molecules and proteins during biomaterial preparation (including tissue engineering scaffolds and drug delivery systems), storage, and in vivo implantation. We also review the physical and chemical aspects of polymer-based stabilization approaches, with a particular focus on the stabilizing properties of silk fibroin biomaterials.     

A Chemical and Enzymatic Approach to Study Site-Specific Sumoylation

A variety of cellular pathways are regulated by protein modifications with ubiquitin-family proteins. SUMO, the Small Ubiquitin-like MOdifier, is covalently attached to lysine on target proteins via a cascade reaction catalyzed by E1, E2, and E3 enzymes. A major barrier to understanding the diverse regulatory roles of SUMO has been a lack of suitable methods to identify protein sumoylation sites. Here we developed a mass-spectrometry (MS) based approach combining chemical and enzymatic modifications to identify sumoylation sites. We applied this method to analyze the auto-sumoylation of the E1 enzyme in vitro and compared it to the GG-remnant method using Smt3-I96R as a substrate. We further examined the effect of smt3-I96R mutation in vivo and performed a proteome-wide analysis of protein sumoylation sites in Saccharomyces cerevisiae. To validate these findings, we confirmed several sumoylation sites of Aos1 and Uba2 in vivo. Together, these results demonstrate that our chemical and enzymatic method for identifying protein sumoylation sites provides a useful tool and that a combination of methods allows a detailed analysis of protein sumoylation sites.