Sequence variation in ligand binding sites in proteins

BACKGROUND: The recent explosion in the availability of complete genome sequences has led to the cataloging of tens of thousands of new proteins and putative proteins. Many of these proteins can be structurally or functionally categorized from sequence conservation alone. In contrast, little attention has been given to the meaning of poorly-conserved sites in families of proteins, which are typically assumed to be of little structural or functional importance. RESULTS: Recently, using statistical free energy analysis of tetratricopeptide repeat (TPR) domains, we observed that positions in contact with peptide ligands are more variable than surface positions in general. Here we show that statistical analysis of TPRs, ankyrin repeats, Cys(2)His(2 )zinc fingers and PDZ domains accurately identifies specificity-determining positions by their sequence variation. Sequence variation is measured as deviation from a neutral reference state, and we present probabilistic and information theory formalisms that improve upon recently suggested methods such as statistical free energies and sequence entropies. CONCLUSION: Sequence variation has been used to identify functionally-important residues in four selected protein families. With TPRs and ankyrin repeats, protein families that bind highly diverse ligands, the effect is so pronounced that sequence "hypervariation" alone can be used to predict ligand binding sites

Fe-EDTA-Bisamide and Fe-ADR-925, The Iron-Bound Hydrolysis Product of the Cardioprotective Agent Dexrazoxane, Cleave DNA Via the Hydroxyl Radical

Use of the antitumor drug doxorubicin is limited by cardiomyopathic side-effects which are believed to be due to iron-mediated hydroxyl radical generation. Dexrazoxane reduces this cardiotoxicity, possibly by removal of iron from doxorubicin by the EDTA-like hydrolysis product of dexrazoxane, ADR-925. However, EDTA-diimides like dexrazoxane, previously used as antitumor agents, are themselves carcinogenic, and recent studies have found that Fe-ADR-925 can also promote hydroxyl radical production. This study demonstrates that, like Fe-EDTA, Fe-ADR-925 and a related desmethyl complex can cleave plasmid DNA under Fenton conditions, and suggests by radical scavenger study that this cleavage is probably via the hydroxyl radical. Differences in DNA cleavage dependence upon concentrations of Fe-EDTA, Fe-ADR-925 and Fe-EDTA-bisamide can be explained by differences in the solution chemistry of the complexes

Targeting the absence: Homozygous DNA deletions as immutable signposts for cancer therapy

Many cancers harbor homozygous DNA deletions (HDs). In contrast to other attributes of cancer cells, their HDs are immutable features that cannot change during tumor progression or therapy. I describe an approach, termed deletion-specific targeting (DST), that employs HDs (not their effects on RNA/protein circuits, but deletions themselves) as the targets of cancer therapy. The DST strategy brings together both existing and new methodologies, including the ubiquitin fusion technique, the split-ubiquitin assay, zinc-finger DNA-recognizing proteins and split restriction nucleases. The DST strategy also employs a feedback mechanism that receives input from a circuit operating as a Boolean OR gate and involves the activation of split nucleases, which destroy DST vector in normal (nontarget) cells. The logic of DST makes possible an incremental and essentially unlimited increase in the selectivity of therapy. If DST strategy can be implemented in a clinical setting, it may prove to be curative and substantially free of side effects

Addressing the unmet need for visualizing Conditional Random Fields in Biological Data

Background: The biological world is replete with phenomena that appear to be ideally modeled and analyzed by one archetypal statistical framework - the Graphical Probabilistic Model (GPM). The structure of GPMs is a uniquely good match for biological problems that range from aligning sequences to modeling the genome-to-phenome relationship. The fundamental questions that GPMs address involve making decisions based on a complex web of interacting factors. Unfortunately, while GPMs ideally fit many questions in biology, they are not an easy solution to apply. Building a GPM is not a simple task for an end user. Moreover, applying GPMs is also impeded by the insidious fact that the complex web of interacting factors inherent to a problem might be easy to define and also intractable to compute upon. Discussion: We propose that the visualization sciences can contribute to many domains of the bio-sciences, by developing tools to address archetypal representation and user interaction issues in GPMs, and in particular a variety of GPM called a Conditional Random Field(CRF). CRFs bring additional power, and additional complexity, because the CRF dependency network can be conditioned on the query data. Conclusions: In this manuscript we examine the shared features of several biological problems that are amenable to modeling with CRFs, highlight the challenges that existing visualization and visual analytics paradigms induce for these data, and document an experimental solution called StickWRLD which, while leaving room for improvement, has been successfully applied in several biological research projects.Comment: BioVis 2014 conferenc

Consensus mutagenesis reveals that non-helical regions influence thermal stability of horseradish peroxidase

The enzyme horseradish peroxidase has many uses in biotechnology but a stabilized derivative would have even wider applicability. To enhance thermal stability, we applied consensus mutagenesis (used successfully with other proteins) to recombinant horseradish peroxidase and generated five single-site mutants. Unexpectedly, these mutations had greater effects on steady-state kinetics than on thermal stability. Only two mutants (T102A, T110V) marginally exceeded the wild type's thermal stability (4% and 10% gain in half-life at 50 °C respectively); the others (Q106R, Q107D, I180F) were less stable than wild type. Stability of a five-fold combination mutant matched that of Q106R, the least-stable single mutant. These results were perplexing: the Class III plant peroxidases display wide differences in thermal stability, yet the consensus mutations failed to reflect these natural variations. We examined the sequence content of Class III peroxidases to determine if there are identifiable molecular reasons for the stability differences observed. Bioinformatic analysis validated our choice of sites and mutations and generated an archetypal peroxidase sequence for comparison with extant sequences. It seems that both genetic variation and differences in protein stability are confined to non-helical regions due to the presence of a highly conserved alpha-helical structural scaffold in these enzymes

Solubilization and Humanization of Paraoxonase-1

Paraoxonase-1 (PON1) is a serum protein, the activity of which is related to susceptibility to cardiovascular disease and intoxication by organophosphorus (OP) compounds. It may also be involved in innate immunity, and it is a possible lead molecule in the development of a catalytic bioscavenger of OP pesticides and nerve agents. Human PON1 expressed in E. coli is mostly found in the insoluble fraction, which motivated the engineering of soluble variants, such as G2E6, with more than 50 mutations from huPON1. We examined the effect on the solubility, activity, and stability of three sets of mutations designed to solubilize huPON1 with fewer overall changes: deletion of the N-terminal leader, polar mutations in the putative HDL binding site, and selection of the subset of residues that became more polar in going from huPON1 to G2E6. All three sets of mutations increase the solubility of huPON1; the HDL-binding mutant has the largest effect on solubility, but it also decreases the activity and stability the most. Based on the G2E6 polar mutations, we “humanized” an engineered variant of PON1 with high activity against cyclosarin (GF) and found that it was still very active against GF with much greater similarity to the human sequence

ConPlex: a server for the evolutionary conservation analysis of protein complex structures

Evolutionary conservation analyses are important for the identification of protein–protein interactions. For protein complex structures, sequence conservation has been applied to determine protein oligomerization states, to characterize native interfaces from non-specific crystal contacts, and to discriminate near-native structures from docking artifacts. However, a user-friendly web-based service for evolutionary conservation analysis of protein complexes has not been available. Therefore, we developed ConPlex (http://sbi.postech.ac.kr/ConPlex/) a web application that enables evolutionary conservation analyses of protein interactions within protein quaternary structures. Users provide protein complex structures; ConPlex automatically identifies protein interfaces and carries out evolutionary conservation analyses for the interface regions. Moreover, ConPlex allows the results of the residue-specific conservation analysis to be displayed on the protein complex structure and provides several options to customize the display output to fit each user’s needs. We believe that ConPlex offers a convenient platform to analyze protein complex structures based on evolutionary conservation of protein–protein interface residues

Characterization of an Italian founder mutation in the RING-finger domain of BRCA1

The identification of founder mutations in cancer predisposing genes is important to improve risk assessment in geographically defined populations, since it may provide specific targets resulting in cost-effective genetic testing. Here, we report the characterization of the BRCA1 c.190T>C (p.Cys64Arg) mutation, mapped to the RING-finger domain coding region, that we detected in 43 hereditary breast/ovarian cancer (HBOC) families, for the large part originating from the province of Bergamo (Northern Italy). Haplotype analysis was performed in 21 families, and led to the identification of a shared haplotype extending over three BRCA1-associated marker loci (0.4 cM). Using the DMLE+2.2 software program and regional population demographic data, we were able to estimate the age of the mutation to vary between 3,100 and 3,350 years old. Functional characterization of the mutation was carried out at both transcript and protein level. Reverse transcriptase-PCR analysis on lymphoblastoid cells revealed expression of full length mRNA from the mutant allele. A green fluorescent protein (GFP)-fragment reassembly assay showed that the p.Cys64Arg substitution prevents the binding of the BRCA1 protein to the interacting protein BARD1, in a similar way as proven deleterious mutations in the RING-domain. Overall, 55 of 83 (66%) female mutation carriers had a diagnosis of breast and/or ovarian cancer. Our observations indicate that the BRCA1 c.190T>C is a pathogenic founder mutation present in the Italian population. Further analyses will evaluate whether screening for this mutation can be suggested as an effective strategy for the rapid identification of at-risk individuals in the Bergamo area

Reduced amino acid alphabets exhibit an improved sensitivity and selectivity in fold assignment

Motivation: Many proteins with vastly dissimilar sequences are found to share a common fold, as evidenced in the wealth of structures now available in the Protein Data Bank. One idea that has found success in various applications is the concept of a reduced amino acid alphabet, wherein similar amino acids are clustered together. Given the structural similarity exhibited by many apparently dissimilar sequences, we undertook this study looking for improvements in fold recognition by comparing protein sequences written in a reduced alphabet. Results: We tested over 150 of the amino acid clustering schemes proposed in the literature with all-versus-all pairwise sequence alignments of sequences in the Distance matrix alignment (DALI) database. We combined several metrics from information retrieval popular in the literature: mean precision, area under the Receiver Operating Characteristic curve and recall at a fixed error rate and found that, in contrast to previous work, reduced alphabets in many cases outperform full alphabets. We find that reduced alphabets can perform at a level comparable to full alphabets in correct pairwise alignment of sequences and can show increased sensitivity to pairs of sequences with structural similarity but low-sequence identity. Based on these results, we hypothesize that reduced alphabets may also show performance gains with more sophisticated methods such as profile and pattern searches