Global sequencing of proteolytic cleavage sites in apoptosis by specific labeling of protein N termini

The nearly 600 proteases in the human genome regulate a diversity of biological processes, including programmed cell death. Comprehensive characterization of protease signaling in complex biological samples is limited by available proteomic methods. We have developed a general approach for global identification of proteolytic cleavage sites using an engineered enzyme to selectively biotinylate free protein N termini for positive enrichment of corresponding N-terminal peptides. Using this method to study apoptosis, we have sequenced 333 caspase-like cleavage sites distributed among 292 protein substrates. These sites are generally not predicted by in vitro caspase substrate specificity but can be used to predict other physiological caspase cleavage sites. Structural bioinformatic studies show that caspase cleavage sites often appear in surface-accessible loops and even occasionally in helical regions. Strikingly, we also find that a disproportionate number of caspase substrates physically interact, suggesting that these dimeric proteases target protein complexes and networks to elicit apoptosis.     

EVA: continuous automatic evaluation of protein structure prediction servers.

Evaluation of protein structure prediction methods is difficult and time-consuming. Here, we describe EVA, a web server for assessing protein structure prediction methods, in an automated, continuous and large-scale fashion. Currently, EVA evaluates the performance of a variety of prediction methods available through the internet. Every week, the sequences of the latest experimentally determined protein structures are sent to prediction servers, results are collected, performance is evaluated, and a summary is published on the web. EVA has so far collected data for more than 3000 protein chains. These results may provide valuable insight to both developers and users of prediction methods. AVAILABILITY: http://cubic.bioc.columbia.edu/eva. CONTACT: eva@cubic.bioc.columbia.edu     

Statistical potentials for fold assessment

A protein structure model generally needs to be evaluated to assess whether or not it has the correct fold. To improve fold assessment, four types of a residue-level statistical potential were optimized, including distance-dependent, contact, Phi/Psi dihedral angle, and accessible surface statistical potentials. Approximately 10,000 test models with the correct and incorrect folds were built by automated comparative modeling of protein sequences of known structure. The criterion used to discriminate between the correct and incorrect models was the Z-score of the model energy. The performance of a Z-score was determined as a function of many variables in the derivation and use of the corresponding statistical potential. The performance was measured by the fractions of the correctly and incorrectly assessed test models. The most discriminating combination of any one of the four tested potentials is the sum of the normalized distance-dependent and accessible surface potentials. The distance-dependent potential that is optimal for assessing models of all sizes uses both C(alpha) and C(beta) atoms as interaction centers, distinguishes between all 20 standard residue types, has the distance range of 30 A, and is derived and used by taking into account the sequence separation of the interacting atom pairs. The terms for the sequentially local interactions are significantly less informative than those for the sequentially nonlocal interactions. The accessible surface potential that is optimal for assessing models of all sizes uses C(beta) atoms as interaction centers and distinguishes between all 20 standard residue types. The performance of the tested statistical potentials is not likely to improve significantly with an increase in the number of known protein structures used in their derivation. The parameters of fold assessment whose optimal values vary significantly with model size include the size of the known protein structures used to derive the potential and the distance range of the accessible surface potential. Fold assessment by statistical potentials is most difficult for the very small models. This difficulty presents a challenge to fold assessment in large-scale comparative modeling, which produces many small and incomplete models. The results described in this study provide a basis for an optimal use of statistical potentials in fold assessment.     

Combined Effect of Dietary Cadmium and Benzo(a)pyrene on Metallothionein Induction and Apoptosis in the Liver and Kidneys of Bank Voles

Bank voles free living in a contaminated environment have been shown to be more sensitive to cadmium (Cd) toxicity than the rodents exposed to Cd under laboratory conditions. The objective of this study was to find out whether benzo(a)pyrene (BaP), a common environmental co-contaminant, increases Cd toxicity through inhibition of metallothionein (MT) synthesis-a low molecular weight protein that is considered to be primary intracellular component of the protective mechanism. For 6 weeks, the female bank voles were provided with diet containing Cd [less than 0.1 μg/g (control) and 60 μg/g dry wt.] and BaP (0, 5, and 10 μg/g dry wt.) alone or in combination. At the end of exposure period, apoptosis and analyses of MT, Cd, and zinc (Zn) in the liver and kidneys were carried out. Dietary BaP 5 μg/g did not affect but BaP 10 μg/g potentiated rather than inhibited induction of hepatic and renal MT by Cd, and diminished Cd-induced apoptosis in both organs. The hepatic and renal Zn followed a pattern similar to that of MT, attaining the highest level in the Cd + BaP 10-μg/g group. These data indicate that dietary BaP attenuates rather than exacerbates Cd toxicity in bank voles, probably by potentiating MT synthesis and increasing Zn concentration in the liver and kidneys.     

Alignment of protein sequences by their profiles

The accuracy of an alignment between two protein sequences can be improved by including other detectably related sequences in the comparison. We optimize and benchmark such an approach that relies on aligning two multiple sequence alignments, each one including one of the two protein sequences. Thirteen different protocols for creating and comparing profiles corresponding to the multiple sequence alignments are implemented in the SALIGN command of MODELLER. A test set of 200 pairwise, structure-based alignments with sequence identities below 40% is used to benchmark the 13 protocols as well as a number of previously described sequence alignment methods, including heuristic pairwise sequence alignment by BLAST, pairwise sequence alignment by global dynamic programming with an affine gap penalty function by the ALIGN command of MODELLER, sequence-profile alignment by PSI-BLAST, Hidden Markov Model methods implemented in SAM and LOBSTER, pairwise sequence alignment relying on predicted local structure by SEA, and multiple sequence alignment by CLUSTALW and COMPASS. The alignment accuracies of the best new protocols were significantly better than those of the other tested methods. For example, the fraction of the correctly aligned residues relative to the structure-based alignment by the best protocol is 56%, which can be compared with the accuracies of 26%, 42%, 43%, 48%, 50%, 49%, 43%, and 43% for the other methods, respectively. The new method is currently applied to large-scale comparative protein structure modeling of all known sequences.     

PIBASE: a comprehensive database of structurally defined protein interfaces

MOTIVATION: In recent years, the Protein Data Bank (PDB) has experienced rapid growth. To maximize the utility of the high resolution protein-protein interaction data stored in the PDB, we have developed PIBASE, a comprehensive relational database of structurally defined interfaces between pairs of protein domains. It is composed of binary interfaces extracted from structures in the PDB and the Probable Quaternary Structure server using domain assignments from the Structural Classification of Proteins and CATH fold classification systems. RESULTS: PIBASE currently contains 158,915 interacting domain pairs between 105,061 domains from 2125 SCOP families. A diverse set of geometric, physiochemical and topologic properties are calculated for each complex, its domains, interfaces and binding sites. A subset of the interface properties are used to remove interface redundancy within PDB entries, resulting in 20,912 distinct domain-domain interfaces. The complexes are grouped into 989 topological classes based on their patterns of domain-domain contacts. The binary interfaces and their corresponding binding sites are categorized into 18,755 and 30,975 topological classes, respectively, based on the topology of secondary structure elements. The utility of the database is illustrated by outlining several current applications. AVAILABILITY: The database is accessible via the world wide web at http://salilab.org/pibase SUPPLEMENTARY INFORMATION: http://salilab.org/pibase/suppinfo.html.     

Analysis of allosteric signal transduction mechanisms in an engineered fluorescent maltose biosensor

We previously reported the construction of a family of reagentless fluorescent biosensor proteins by the structure-based design of conjugation sites for a single, environmentally sensitive small molecule dye, thus providing a mechanism for the transduction of ligand-induced conformational changes into a macroscopic fluorescence observable. Here we investigate the microscopic mechanisms that may be responsible for the macroscopic fluorescent changes in such Fluorescent Allosteric Signal Transduction (FAST) proteins. As case studies, we selected three individual cysteine mutations (F92C, D95C, and S233C) of Escherichia coli maltose binding protein (MBP) covalently labeled with a single small molecule fluorescent probe, N-((2-iodoacetoxy)ethyl)-N-methyl)amino-7-nitrobenz-2-oxa-1,3-diazole (NBD), each giving rise to a robust FAST protein with a distinct maltose-dependent fluorescence response. The fluorescence emission intensity, anisotropy, lifetime, and iodide-dependent fluorescence quenching were determined for each conjugate in the presence and absence of maltose. Structure-derived solvent accessible surface areas of the three FAST proteins are consistent with experimentally observed quenching data. The D95C protein exhibits the largest fluorescence change upon maltose binding. This mutant was selected for further characterization, and residues surrounding the fluorophore coupling site were mutagenized. Analysis of the resulting mutant FAST proteins suggests that specific hydrogen-bonding interactions between the fluorophore molecule and two tyrosine side-chains, Tyr171 and Tyr176, in the open state but not the closed, are responsible for the dramatic fluorescence response of this construct. Taken together these results provide insights that can be used in future design cycles to construct fluorescent biosensors that optimize signaling by engineering specific hydrogen bonds between a fluorophore and protein.     

Definition of general topological equivalence in protein structures. A procedure involving comparison of properties and relationships through simulated annealing and dynamic programming

A protein is defined as an indexed string of elements at each level in the hierarchy of protein structure: sequence, secondary structure, super-secondary structure, etc. The elements, for example, residues or secondary structure segments such as helices or beta-strands, are associated with a series of properties and can be involved in a number of relationships with other elements. Element-by-element dissimilarity matrices are then computed and used in the alignment procedure based on the sequence alignment algorithm of Needleman & Wunsch, expanded by the simulated annealing technique to take into account relationships as well as properties. The utility of this method for exploring the variability of various aspects of protein structure and for comparing distantly related proteins is demonstrated by multiple alignment of serine proteinases, aspartic proteinase lobes and globins.     

Structural-functional organization of the SV40 virus chromosome. II. Characteristics of viral genome containing deletion in the regulatory region

Structure of SV40 wild type virus genome differing from the well known 776 strain has been characterized. This particular strain (SV40) was used previously for viral chromatin analysis. We show here that SV40 strain differs from 776 strain by deletion in the beginning of the "late" region (enhancer). Restriction nucleases mapping and nucleotide sequencing through this region reveal the absence of one full-length copy of 72 bp repeat.     

Trajectories of ADHD severity over 10 years from childhood into adulthood

The study examined unique trajectories of ADHD severity from childhood (7–16 yo at baseline) through adulthood in a sample of ADHD, bipolar and healthy subjects. Comorbid disorders and temperament were examined as correlates of course of ADHD. N = 81 participants with an ADHD diagnosis, ascertained as a comparison group in a study of bipolar disorder (BP-I), were followed over a 10-year period. Growth mixture modeling (GMM) of ADHD severity was used to investigate trajectories of ADHD severity over 10 years. GMM revealed four trajectories in the N = 251 participants included in these analyses. A persisting high ADHD trajectory had the highest rates of comorbid major depressive disorder and oppositional defiant disorder. This persisting high ADHD group also had higher fantasy and lower persistence and self-directedness compared with those who displayed a pattern of decreasing ADHD symptoms over time. Psychopathologic features that characterize divergent trajectories of ADHD into adulthood are elucidated, and additional, larger studies are warranted.