Total sequence decomposition distinguishes functional modules, "molegos" in apurinic/apyrimidinic endonucleases

The metal dependence of single-metal mediated phosphodiester bond cleavage: a QM/MM study of a multifaceted human enzyme

Nucleases catalyze the cleavage of phosphodiester bonds in nucleic acids using a range of metal cofactors. Although it is well accepted that many nucleases rely on two metal ions, the one-metal mediated pathway is debated. Furthermore, one-metal mediated nucleases maintain activity in the presence of many different metals, but the underlying reasons for this broad metal specificity are unknown. The human apurinic/apyrimidinic endonuclease (APE1), which plays a key role in DNA repair, transcription regulation, and gene expression, is a prototypical example of a one-metal dependent nuclease. Although Mg²⁺ is the native metal cofactor, APE1 remains catalytically active in the presence of several metals, with the rate decreasing as Mg²⁺ > Mn²⁺ > Ni²⁺ > Zn²⁺, while Ca²⁺ completely abolished the activity. The present work uses quantum mechanics-molecular mechanics techniques to map APE1-facilitated phosphodiester bond hydrolysis in the presence of these metals. The structural differences in stationary points along the reaction pathway shed light on the interplay between several factors that allow APE1 to remain catalytically active for various metals, with the trend in the barrier heights correlating with the experimentally reported APE1 catalytic activity. In contrast, Ca²⁺ significantly changes the metal coordination and active site geometry, and thus completely inhibits catalysis. Our work thereby provides support for the controversial single-metal mediated phosphodiester bond cleavage and clarifies uncertainties regarding the role of the metal and metal identity in this important reaction. This information is key for future medicinal and biotechnological applications including disease diagnosis and treatment, and protein engineering.

Membrane interaction and functional plasticity of inositol polyphosphate 5-phosphatases

In this issue of Structure, Trésaugues and colleagues determined the interaction of membrane-bound phosphoinositides with three clinically significant human inositol polyphosphate 5-phosphatases (I5Ps). A comparison to the structures determined with soluble substrates revealed differences in the binding mode and suggested how the I5Ps and apurinic endonuclease (APE1) activities evolved from the same metal-binding active center.

AllerTOP v.2--a server for in silico prediction of allergens

Allergy is an overreaction by the immune system to a previously encountered, ordinarily harmless substance--typically proteins--resulting in skin rash, swelling of mucous membranes, sneezing or wheezing, or other abnormal conditions. The use of modified proteins is increasingly widespread: their presence in food, commercial products, such as washing powder, and medical therapeutics and diagnostics, makes predicting and identifying potential allergens a crucial societal issue. The prediction of allergens has been explored widely using bioinformatics, with many tools being developed in the last decade; many of these are freely available online. Here, we report a set of novel models for allergen prediction utilizing amino acid E-descriptors, auto- and cross-covariance transformation, and several machine learning methods for classification, including logistic regression (LR), decision tree (DT), naïve Bayes (NB), random forest (RF), multilayer perceptron (MLP) and k nearest neighbours (kNN). The best performing method was kNN with 85.3% accuracy at 5-fold cross-validation. The resulting model has been implemented in a revised version of the AllerTOP server (http://www.ddg-pharmfac.net/AllerTOP).

Physicochemical property consensus sequences for functional analysis, design of multivalent antigens and targeted antivirals

Background

Analysis of large sets of biological sequence data from related strains or organisms is complicated by superficial redundancy in the set, which may contain many members that are identical except at one or two positions. Thus a new method, based on deriving physicochemical property (PCP)-consensus sequences, was tested for its ability to generate reference sequences and distinguish functionally significant changes from background variability.

Methods

The PCP consensus program was used to automatically derive consensus sequences starting from sequence alignments of proteins from Flaviviruses (from the Flavitrack database) and human enteroviruses, using a five dimensional set of Eigenvectors that summarize over 200 different scalar values for the PCPs of the amino acids. A PCP-consensus protein of a Dengue virus envelope protein was produced recombinantly and tested for its ability to bind antibodies to strains using ELISA.

Results

PCP-consensus sequences of the flavivirus family could be used to classify them into five discrete groups and distinguish areas of the envelope proteins that correlate with host specificity and disease type. A multivalent Dengue virus antigen was designed and shown to bind antibodies against all four DENV types. A consensus enteroviral VPg protein had the same distinctive high pKa as wild type proteins and was recognized by two different polymerases.

Conclusions

The process for deriving PCP-consensus sequences for any group of aligned similar sequences, has been validated for sequences with up to 50% diversity. Ongoing projects have shown that the method identifies residues that significantly alter PCPs at a given position, and might thus cause changes in function or immunogenicity. Other potential applications include deriving target proteins for drug design and diagnostic kits.

Simplifying complex sequence information: a PCP-consensus protein binds antibodies against all four Dengue serotypes

Designing proteins that reflect the natural variability of a pathogen is essential for developing novel vaccines and drugs. Flaviviruses, including Dengue (DENV) and West Nile (WNV), evolve rapidly and can "escape" neutralizing monoclonal antibodies by mutation. Designing antigens that represent many distinct strains is important for DENV, where infection with a strain from one of the four serotypes may lead to severe hemorrhagic disease on subsequent infection with a strain from another serotype. Here, a DENV physicochemical property (PCP)-consensus sequence was derived from 671 unique sequences from the Flavitrack database. PCP-consensus proteins for domain 3 of the envelope protein (EdomIII) were expressed from synthetic genes in Escherichia coli. The ability of the purified consensus proteins to bind polyclonal antibodies generated in response to infection with strains from each of the four DENV serotypes was determined. The initial consensus protein bound antibodies from DENV-1-3 in ELISA and Western blot assays. This sequence was altered in 3 steps to incorporate regions of maximum variability, identified as significant changes in the PCPs, characteristic of DENV-4 strains. The final protein was recognized by antibodies against all four serotypes. Two amino acids essential for efficient binding to all DENV antibodies are part of a discontinuous epitope previously defined for a neutralizing monoclonal antibody. The PCP-consensus method can significantly reduce the number of experiments required to define a multivalent antigen, which is particularly important when dealing with pathogens that must be tested at higher biosafety levels.

MD simulation and experimental evidence for Mg²+ binding at the B site in human AP endonuclease 1

Apurinic/apyrimidinic endonuclease 1 (APE1), a central enzyme in the base excision repair pathway, cleaves damaged DNA in Mg(2+) dependent reaction. Despite characterization of nine X-ray crystallographic structures of human APE1, in some cases, bound to various metal ions and substrate/product, the position of the metal ion and its stoichiometry for the cleavage reaction are still being debated. While a mutation of the active site E96Q was proposed to eliminate Mg(2+) binding at the "A" site, we show experimentally that this mutant still requires Mg(2+) at concentration similar to that for the wild type enzyme to cleave the AP site in DNA. Molecular dynamics simulations of the wild type APE1, E96Q and a double missense mutant E96Q + D210N indicate that Mg(2+) placed at the A-site destabilizes the bound AP site-containing DNA. In these simulations, the H-bond chain D238-H309-AP site oxygen is broken and the substrate DNA is shifted away from its crystal structure position (1DE9). In contrast, simulations with the Mg(2+) at site B or A+B sites leave the substrate DNA at the position shown in the crystal structure (1DE9). Taken together our MD simulations and biochemical analysis suggests that Mg(2+) binding at the B site is involved in the reaction mechanism associated with endonuclease function of APE1.

PCP consensus sequences of flaviviruses: correlating variance with vector competence and disease phenotype

BACKGROUND

Computational methods are needed to design multivalent vaccines against flaviviruses (FVs) such as the West Nile virus or the dengue virus (DENV).

OBJECTIVE

We aimed to use physicochemical property (PCP) consensus sequences of FV strains to delineate conserved motifs, areas of maximum variability, and specific loci that correlate with arthropod vector, serotype, and disease severity.

METHODS

PCP consensus sequences for 27 species were prepared from 928 annotated sequences catalogued in Flavitrack. Alignments of these correlated well with the known structures of the NS3 protease domain and envelope (E) proteins. The PCPMer suite was used to identify motifs common to all FVs. Areas of PCP variability that correlated with phenotype were plotted on the structures.

RESULTS

Despite considerable diversity at the amino acid level, PCPs for both proteins were well conserved throughout the FVs. A series of insertions in E separated tick- from mosquito-borne viruses and all arthropod-borne viruses from isolates with no known vector or directly from insects. Comparison of a PCP consensus sequence of E derived from 600 DENV strains (DENV600) with individual ones for DENV1-DENV4 showed that most major serotype-specific variation occurs near these insertions. The DENV600 differed from one prepared from eight hemorrhagic or fatal strains from four DENV serotypes at only three positions, two of which overlap known escape mutant sites.

CONCLUSIONS

Comparing consensus sequences showed that substantial changes occur in only a few areas of the E protein. PCP consensus sequences can contribute to the design of multivalent vaccines.

Characteristic motifs for families of allergenic proteins

The identification of potential allergenic proteins is usually done by scanning a database of allergenic proteins and locating known allergens with a high sequence similarity. However, there is no universally accepted cut-off value for sequence similarity to indicate potential IgE cross-reactivity. Further, overall sequence similarity may be less important than discrete areas of similarity in proteins with homologous structure. To identify such areas, we first classified all allergens and their subdomains in the Structural Database of Allergenic Proteins (SDAP, http://fermi.utmb.edu/SDAP/) to their closest protein families as defined in Pfam, and identified conserved physicochemical property motifs characteristic of each group of sequences. Allergens populate only a small subset of all known Pfam families, as all allergenic proteins in SDAP could be grouped to only 130 (of 9318 total) Pfams, and 31 families contain more than four allergens. Conserved physicochemical property motifs for the aligned sequences of the most populated Pfam families were identified with the PCPMer program suite and catalogued in the webserver MotifMate (http://born.utmb.edu/motifmate/summary.php). We also determined specific motifs for allergenic members of a family that could distinguish them from non-allergenic ones. These allergen specific motifs should be most useful in database searches for potential allergens. We found that sequence motifs unique to the allergens in three families (seed storage proteins, Bet v 1, and tropomyosin) overlap with known IgE epitopes, thus providing evidence that our motif based approach can be used to assess the potential allergenicity of novel proteins.

Identification and analysis of conserved sequence motifs in cytochrome P450 family 2. Functional and structural role of a motif 187RFDYKD192 in CYP2B enzymes

Using a multiple alignment of 175 cytochrome P450 (CYP) family 2 sequences, 20 conserved sequence motifs (CSMs) were identified with the program PCPMer. Functional importance of the CSM in CYP2B enzymes was assessed from available data on site-directed mutants and genetic variants. These analyses suggested an important role of the CSM 8, which corresponds to(187)RFDYKD(192) in CYP2B4. Further analysis showed that residues 187, 188, 190, and 192 have a very high rank order of conservation compared with 189 and 191. Therefore, eight mutants (R187A, R187K, F188A, D189A, Y190A, K191A, D192A, and a negative control K186A) were made in an N-terminal truncated and modified form of CYP2B4 with an internal mutation, which is termed 2B4dH/H226Y. Function was examined with the substrates 7-methoxy-4-(trifluoromethyl)coumarin (7-MFC), 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC), 7-benzyloxy-4-(trifluoromethyl)coumarin (7-BFC), and testosterone and with the inhibitors 4-(4-chlorophenyl)imidazole (4-CPI) and bifonazole (BIF). Compared with the template and K186A, the mutants R187A, R187K, F188A, Y190A, and D192A showed > or =2-fold altered substrate specificity, k(cat), K(m), and/or k(cat)/K(m) for 7-MFC and 7-EFC and 3- to 6-fold decreases in differential inhibition (IC(50,BIF)/IC(50,4-CPI)). Subsequently, these mutants displayed 5-12 degrees C decreases in thermal stability (T(m)) and 2-8 degrees C decreases in catalytic tolerance to temperature (T(50)) compared with the template and K186A. Furthermore, when R187A and D192A were introduced in CYP2B1dH, the P450 expression and thermal stability were decreased. In addition, R187A showed increased activity with 7-EFC and decreased IC(50,BIF)/IC(50,4-CPI) compared with 2B1dH. Analysis of long range residue-residue interactions in the CYP2B4 crystal structures indicated strong hydrogen bonds involving Glu(149)-Asn(177)-Arg(187)-Tyr(190) and Asp(192)-Val(194), which were significantly-reduced/abolished by the Arg(187)-->Ala and Asp(192)-->Alasubstitutions, respectively.

Characterization of Mg2+ binding to the DNA repair protein apurinic/apyrimidic endonuclease 1 via solid-state 25Mg NMR spectroscopy

Apurinic/apyrimidinic endonuclease 1 (APE1), a member of the divalent cation-dependent phosphoesterase superfamily of proteins that retain the conserved four-layered alpha/beta-sandwich structural core, is an essential protein that functions as part of base excision repair to remove mutagenic and cytotoxic abasic sites from DNA. Using low-temperature solid-state (25)Mg NMR spectroscopy and various mutants of APE1, we demonstrate that Mg(2+) binds to APE1 and a functional APE1-substrate DNA complex with an overall stoichiometry of one Mg(2+) per mole of APE1 as predicted by the X-ray work of Tainer and co-workers (Mol, C. D.; Kuo, C. F.; Thayer, M. M.; Cunningham, R. P.; Tainer, J. A. Nature 1995, 374 , 381-386). However, the NMR spectra show that the single Mg(2+) site is disordered. We discuss the probable reasons for the disorder at the Mg(2+) binding site. The most likely source of this disorder is arrangement of the protein-ligands about the Mg(2+) (cis and trans isomers). The existence of these isomers reinforces the notion of the plasticity of the metal binding site within APE1.

Unusual role of a cysteine residue in substrate binding and activity of human AP-endonuclease 1

The mammalian AP-endonuclease (APE1) repairs apurinic/apyrimidinic (AP) sites and strand breaks with 3' blocks in the genome that are formed both endogenously and as intermediates during base excision repair. APE1 has an unrelated activity as a redox activator (and named Ref-1) for several trans-acting factors. In order to identify whether any of the seven cysteine residues in human APE1 affects its enzymatic function, we substituted these singly or multiply with serine. The repair activity is not affected in any of the mutants except those with C99S mutation. The Ser99-containing mutant lost affinity for DNA and its activity was inhibited by 10 mM Mg(2+). However, the Ser99 mutant has normal activity in 2 mM Mg(2+). Using crystallographic data and molecular dynamics simulation, we have provided a mechanistic basis for the altered properties of the C99S mutant. We earlier predicted that Mg(2+), with potential binding sites A and B, binds at the B site of wild-type APE1-substrate complex and moves to the A site after cleavage occurs, as observed in the crystal structure. The APE1-substrate complex is stabilized by a H bond between His309 and the AP site. We now show that this bond is broken to destabilize the complex in the absence of the Mg(2+). This effect due to the mutation of Cys99, approximately 16 A from the active site, on the DNA binding and activity is surprising. Mg(2+) at the B site promotes stabilization of the C99S mutant complex. At higher Mg(2+) concentration the A site is also filled, causing the B-site Mg(2+) to shift together with the AP site. At the same time, the H bond between His309 and the AP site shifts toward the 5' site of DNA. These shifts could explain the lower activity of the C99S mutant at higher [Mg(2+)]. The unexpected involvement of Cys99 in APE1's substrate binding and catalysis provides an example of involvement of a residue far from the active site.

A "moving metal mechanism" for substrate cleavage by the DNA repair endonuclease APE-1

Apurinic/apyrimidinic endonuclease (APE-1) is essential for base excision repair (BER) of damaged DNA. Here molecular dynamics (MD) simulations of APE1 complexed with cleaved and uncleaved damaged DNA were used to determine the role and position of the metal ion(s) in the active site before and after DNA cleavage. The simulations started from an energy minimized wild-type structure of the metal-free APE1/damaged-DNA complex (1DE8). A grid search with one Mg2+ ion located two low energy clusters of Mg2+ consistent with the experimentally determined metal ion positions. At the start of the longer MD simulations, Mg2+ ions were placed at different positions as seen in the crystal structures and the movement of the ion was followed over the course of the trajectory. Our analysis suggests a "moving metal mechanism" in which one Mg2+ ion moves from the B- (more buried) to the A-site during substrate cleavage. The anticipated inversion of the phosphate oxygens occurs during the in-line cleavage reaction. Experimental results, which show competition between Ca2+ and Mg2+ for catalyzing the reaction, and high concentrations of Mg2+ are inhibitory, indicate that both sites cannot be simultaneously occupied for maximal activity.

Molego-based definition of the architecture and specificity of metal-binding sites

Decomposing proteins into "molegos," building blocks that are conserved in sequence and 3D-structure, can identify functional elements. To demonstrate the specificity of the decomposition method, the PCPMer program suite was used to numerically define physical chemical property motifs corresponding to the molegos that make up the metal-containing active sites of three distinct enzyme families, from the dimetallic phosphatases, DNase 1 related nucleases/phosphatases, and dioxygenases. All three superfamilies bind metal ions in a beta-strand core region but differ in the number and type of ions needed for activity. The motifs were then used to automatically identify proteins in the ASTRAL40 database that contained similar motifs. The proteins with the highest PCPMer score in the database were primarily metal-binding enzymes that were related in function to those in the alignment used to generate the PCPMer motif lists. The proteins that contained motifs similar to the dioxygenases differed from those found with PCP-motifs for phosphatases and nucleases. Relatively few metal-binding enzymes were detected when the search was done with PCP-motifs defined for interleukin-1 related proteins, which have a beta-strand core but do not bind metal ions. While the box architecture was constant in each superfamily, the specificity for the metal ion preferred for enzymatic activity is determined by the pattern of carbonyl, hydroxyl or imadazole groups in key positions in the molegos. These results have implications for the design of metal-binding enzymes, and illustrate the ability of the PCPMer approach to distinguish, at the sequence level, structural and functional elements.

Stereophysicochemical variability plots highlight conserved antigenic areas in Flaviviruses

Background

Flaviviruses, which include Dengue (DV) and West Nile (WN), mutate in response to immune system pressure. Identifying escape mutants, variant progeny that replicate in the presence of neutralizing antibodies, is a common way to identify functionally important residues of viral proteins. However, the mutations typically occur at variable positions on the viral surface that are not essential for viral replication. Methods are needed to determine the true targets of the neutralizing antibodies.

Results

Stereophysicochemical variability plots (SVPs), 3-D images of protein structures colored according to variability, as determined by our PCPMer program, were used to visualize residues conserved in their physical chemical properties (PCPs) near escape mutant positions. The analysis showed 1) that escape mutations in the flavivirus envelope protein are variable residues by our criteria and 2) two escape mutants found at the same position in many flaviviruses sit above clusters of conserved residues from different regions of the linear sequence. Conservation patterns in T-cell epitopes in the NS3- protease suggest a similar mechanism of immune system evasion.

Conclusion

The SVPs add another dimension to structurally defining the binding sites of neutralizing antibodies. They provide a useful aid for determining antigenically important regions and designing vaccines.

The Hotdog fold: wrapping up a superfamily of thioesterases and dehydratases

BACKGROUND

The Hotdog fold was initially identified in the structure of Escherichia coli FabA and subsequently in 4-hydroxybenzoyl-CoA thioesterase from Pseudomonas sp. strain CBS. Since that time structural determinations have shown a number of other apparently unrelated proteins also share the Hotdog fold.

RESULTS

Using sequence analysis we unify a large superfamily of HotDog domains. Membership includes numerous prokaryotic, archaeal and eukaryotic proteins involved in several related, but distinct, catalytic activities, from metabolic roles such as thioester hydrolysis in fatty acid metabolism, to degradation of phenylacetic acid and the environmental pollutant 4-chlorobenzoate. The superfamily also includes FapR, a non-catalytic bacterial homologue that is involved in transcriptional regulation of fatty acid biosynthesis. We have defined 17 subfamilies, with some characterisation. Operon analysis has revealed numerous HotDog domain-containing proteins to be fusion proteins, where two genes, once separate but adjacent open-reading frames, have been fused into one open-reading frame to give a protein with two functional domains. Finally we have generated a Hidden Markov Model library from our analysis, which can be used as a tool for predicting the occurrence of HotDog domains in any protein sequence.

CONCLUSIONS

The HotDog domain is both an ancient and ubiquitous motif, with members found in the three branches of life.

Effects of backbone contacts 3' to the abasic site on the cleavage and the product binding by human apurinic/apyrimidinic endonuclease (APE1)

The mammalian apurinic/apyrimidinic (AP) endonuclease (APE1) is a multifunctional protein that plays essential roles in DNA repair and gene regulation. We decomposed the APEs into 12 blocks of highly conserved sequence and structure (molegos). This analysis suggested that residues in molegos common to all APEs, but not to the less specific nuclease, DNase I, would dictate enhanced binding to damaged DNA. To test this hypothesis, alanine was substituted for N226 and N229, which form hydrogen bonds to the DNA backbone 3' of the AP sites in crystal structures of the APE1/DNA complex. While the cleavage rate at AP sites of both N226A and N229A mutants increased, their ability to bind to damaged DNA decreased. The ability of a double mutant (N226A/N229A) to bind damaged DNA was further decreased, while the V(max) was almost identical to that of the wild-type APE1. A double mutant at N226 and R177, a residue that binds to the same phosphate as N229, had a significantly decreased activity and substrate binding. As the affinity for product DNA was decreased in all the mutants, the enhanced reaction rate of the single mutants could be due to alleviation of product inhibition of the enzyme. We conclude that hydrogen bonds to phosphate groups 3' to the cleavage site is essential for APE1's binding to the product DNA, which may be necessary for efficient functioning of the base excision repair pathway. The results indicate that the molego analysis can aid in the redesign of proteins with altered binding affinity and activity.

Major linear IgE epitopes of mountain cedar pollen allergen Jun a 1 map to the pectate lyase catalytic site

Resolution of the 3D structures and IgE epitopes of allergens may identify common or conserved features of allergens. Jun a 1, the predominant allergen in mountain cedar pollen, was chosen as a model for identifying common structural and functional features among a group of plant allergens. In this study, synthetic, overlapping peptides of Jun a 1 and sera from patients allergic to mountain cedar pollen were used to identify linear epitopes. A 3D model of Jun a 1 was produced using the Bacillus subtiles pectate lyase (PL) as a template and validated with biophysical measurements. This allowed mappings of four IgE binding sites on Jun a 1. Two of the epitopes mapped to turns or loops on the surface of the model structure. The other two epitopes mapped to the beta-sheet region, homologous to the catalytic site of PL. This region of Jun a 1 is highly conserved in the group 1 allergens from other cedar trees as well as microbial PLs. The finding that two out of three major IgE epitopes map to highly conserved catalytic regions of group 1 cedar allergens may help to explain the high degree of cross-reactivity between cedar pollen allergens and might represent a pattern of reactivity common to other allergens with catalytic activity.

Mammalian DNA base excision repair proteins: their interactions and role in repair of oxidative DNA damage

The DNA base excision repair (BER) is a ubiquitous mechanism for removing damage from the genome induced by spontaneous chemical reaction, reactive oxygen species (ROS) and also DNA damage induced by a variety of environmental genotoxicants. DNA repair is essential for maintaining genomic integrity. As we learn more about BER, a more complex mechanism emerges which supersedes the classical, simple pathway requiring only four enzymatic reactions. The key to understand the complete BER process is to elucidate how multiple proteins interact with one another in a coordinated process under specific physiological conditions.