Application and comparison of classification algorithms for recognition of Alzheimer's disease in electrical brain activity (EEG)

The early detection of subjects with probable Alzheimer's disease (AD) is crucial for effective appliance of treatment strategies. Here we explored the ability of a multitude of linear and non-linear classification algorithms to discriminate between the electroencephalograms (EEGs) of patients with varying degree of AD and their age-matched control subjects. Absolute and relative spectral power, distribution of spectral power, and measures of spatial synchronization were calculated from recordings of resting eyes-closed continuous EEGs of 45 healthy controls, 116 patients with mild AD and 81 patients with moderate AD, recruited in two different centers (Stockholm, New York). The applied classification algorithms were: principal component linear discriminant analysis (PC LDA), partial least squares LDA (PLS LDA), principal component logistic regression (PC LR), partial least squares logistic regression (PLS LR), bagging, random forest, support vector machines (SVM) and feed-forward neural network. Based on 10-fold cross-validation runs it could be demonstrated that even tough modern computer-intensive classification algorithms such as random forests, SVM and neural networks show a slight superiority, more classical classification algorithms performed nearly equally well. Using random forests classification a considerable sensitivity of up to 85% and a specificity of 78%, respectively for the test of even only mild AD patients has been reached, whereas for the comparison of moderate AD vs. controls, using SVM and neural networks, values of 89% and 88% for sensitivity and specificity were achieved. Such a remarkable performance proves the value of these classification algorithms for clinical diagnostics.

Specific atomic groups and RNA helix geometry in acceptor stem recognition by a tRNA synthetase

Oligonucleotides that recapitulate the acceptor stems of tRNAs are substrates for aminoacylation by many tRNA synthetases in vitro, even though these substrates are missing the anticodon trinucleotides of the genetic code. In the case of tRNAAla a single acceptor stem G⋅U base pair at position 3·70 is essential, based on experiments where the wobble pair has been replaced by alternatives such as I⋅U, G⋅C, and A⋅U, among others. These experiments led to the conclusion that the minor-groove free 2-amino group (of guanosine) of the G⋅U wobble pair is essential for charging. Moreover, alanine-inserting tRNAs (amber suppressors) that replace G⋅U with mismatches such as G⋅A and C⋅A are partially active in vivo and can support growth of an Escherichia coli tRNAAla knockout strain, leading to the hypothesis that a helix irregularity and nucleotide functionalities are important for recognition. Herein we investigate the charging in vitro of oligonucleotide and full-length tRNA substrates that contain mismatches at the position of the G⋅U pair. Although most of these substrates have undetectable activity, G⋅A and C⋅A variants retain some activity, which is, nevertheless, reduced by at least 100-fold. Thus, the in vivo assays are much less sensitive to large changes in aminoacylation kinetic efficiency of 3·70 variants than is the in vitro assay system. Although these functional data do not clarify all of the details, it is now clear that specific atomic groups are substantially more important in determining kinetic efficiency than is a helical distortion. By implication, the activity of mutant tRNAs measured in the in vivo assays appears to be more dependent on factors other than aminoacylation kinetic efficiency.

Novel dimeric interface and electrostatic recognition in bacterial Cu,Zn superoxide dismutase

Eukaryotic Cu,Zn superoxide dismutases (CuZnSODs) are antioxidant enzymes remarkable for their unusually stable β-barrel fold and dimer assembly, diffusion-limited catalysis, and electrostatic guidance of their free radical substrate. Point mutations of CuZnSOD cause the fatal human neurodegenerative disease amyotrophic lateral sclerosis. We determined and analyzed the first crystallographic structure (to our knowledge) for CuZnSOD from a prokaryote, Photobacterium leiognathi, a luminescent symbiont of Leiognathid fish. This structure, exemplifying prokaryotic CuZnSODs, shares the active-site ligand geometry and the topology of the Greek key β-barrel common to the eukaryotic CuZnSODs. However, the β-barrel elements recruited to form the dimer interface, the strategy used to forge the channel for electrostatic recognition of superoxide radical, and the connectivity of the intrasubunit disulfide bond in P. leiognathi CuZnSOD are discrete and strikingly dissimilar from those highly conserved in eukaryotic CuZnSODs. This new CuZnSOD structure broadens our understanding of structural features necessary and sufficient for CuZnSOD activity, highlights a hitherto unrecognized adaptability of the Greek key β-barrel building block in evolution, and reveals that prokaryotic and eukaryotic enzymes diverged from one primordial CuZnSOD and then converged to distinct dimeric enzymes with electrostatic substrate guidance.

Small antibody-like proteins with prescribed ligand specificities derived from the lipocalin fold

We demonstrate that the ligand pocket of a lipocalin from Pieris brassicae, the bilin-binding protein (BBP), can be reshaped by combinatorial protein design such that it recognizes fluorescein, an established immunological hapten. For this purpose 16 residues at the center of the binding site, which is formed by four loops on top of an eight-stranded β-barrel, were subjected to random mutagenesis. Fluorescein-binding BBP variants were then selected from the mutant library by bacterial phage display. Three variants were identified that complex fluorescein with high affinity, exhibiting dissociation constants as low as 35.2 nM. Notably, one of these variants effects almost complete quenching of the ligand fluorescence, similarly as an anti-fluorescein antibody. Detailed ligand-binding studies and site-directed mutagenesis experiments indicated (i) that the molecular recognition of fluorescein is specific and (ii) that charged residues at the center of the pocket are responsible for tight complex formation. Sequence comparison of the BBP variants directed against fluorescein with the wild-type protein and with further variants that were selected against several other ligands revealed that all of the randomized amino acid positions are variable. Hence, a lipocalin can be used for generating molecular pockets with a diversity of shapes. We term this class of engineered proteins “anticalins.” Their one-domain scaffold makes them a promising alternative to antibodies to create a stable receptor protein for a ligand of choice.

Quantitative assessment of enzyme specificity in vivo: P2 recognition by Kex2 protease defined in a genetic system

The specificity of the yeast proprotein-processing Kex2 protease was examined in vivo by using a sensitive, quantitative assay. A truncated prepro-α-factor gene encoding an α-factor precursor with a single α-factor repeat was constructed with restriction sites for cassette mutagenesis flanking the single Kex2 cleavage site (-SLDKR↓EAEA-). All of the 19 substitutions for the Lys (P2) residue in the cleavage site were made. The wild-type and mutant precursors were expressed in a yeast strain lacking the chromosomal genes encoding Kex2 and prepro-α-factor. Cleavage of the 20 sites by Kex2, expressed at the wild-type level, was assessed by using a quantitative-mating assay with an effective range greater than six orders of magnitude. All substitutions for Lys at P2 decreased mating, from 2-fold for Arg to >106-fold for Trp. Eviction of the Kex2-encoding plasmid indicated that cleavage of mutant sites by other cellular proteases was not a complicating factor. Mating efficiencies of strains expressing the mutant precursors correlated well with the specificity (kcat/KM) of purified Kex2 for comparable model peptide substrates, validating the in vivo approach as a quantitative method. The results support the conclusion that KM, which is heavily influenced by the nature of the P2 residue, is a major determinant of cleavage efficiency in vivo. P2 preference followed the rank order: Lys > Arg > Thr > Pro > Glu > Ile > Ser > Ala > Asn > Val > Cys > AsP > Gln > Gly > His > Met > Leu > Tyr > Phe > Trp.

Resected RNA pseudoknots and their recognition by histidyl-tRNA synthetase

Duplexes constituted by closed or open RNA circles paired to single-stranded oligonucleotides terminating with 3′-CCAOH form resected pseudoknots that are substrates of yeast histidyl-tRNA synthetase. Design of this RNA fold is linked to the mimicry of the pseudoknotted amino acid accepting branch of the tRNA-like domain from brome mosaic virus, known to be charged by tyrosyl-tRNA synthetases, with RNA minihelices recapitulating accepting branches of canonical tRNAs. Prediction of the histidylation function of the new family of minimalist tRNA-like structures relates to the geometry of resected pseudoknots that allows proper presentation to histidyl-tRNA synthetase of analogues of the histidine identity determinants N-1 and N73 present in tRNAs. This geometry is such that the analogue of the major N-1 histidine determinant in the RNA circles faces the analogue of the discriminator N73 nucleotide in the accepting oligonucleotides. The combination of identity elements found in tRNAHis species from archaea, eubacteria, and organelles (G-1/C73) is the most efficient for determining histidylation of the duplexes. The inverse combination (C-1/G73) leads to the worst histidine acceptors with charging efficiencies reduced by 2–3 orders of magnitude. Altogether, these findings open new perspectives for understanding evolution of tRNA identity and serendipitous RNA functions.

Protein component of the ribozyme ribonuclease P alters substrate recognition by directly contacting precursor tRNA

The protein component of ribonuclease P (RNase P) binds to the RNA subunit, forming a functional ribonucleoprotein complex in vivo and enhancing the affinity of the precursor tRNA (pre-tRNA) substrate. Photocrosslinking experiments with pre-tRNA bound to RNase P reconstituted with the protein component of Bacillus subtilis ribonuclease P (P protein) site specifically modified with a crosslinking reagent indicate that: (i) the central cleft of P protein directly interacts with the single-stranded 5′ leader sequence of pre-tRNA, and (ii) the orientation and register of the pre-tRNA leader sequence in the central cleft places the protein component in close proximity to the active site. This unique mode of interaction suggests that the catalytic active site in RNase P occurs near the interface of RNA and protein. In contrast to other ribonucleoprotein complexes where the protein mainly stabilizes the active tertiary fold of the RNA, a critical function of the protein component of RNase P is to alter substrate specificity and enhance catalytic efficiency.

Efficient recognition of substrates and substrate analogs by the adenine glycosylase MutY requires the C-terminal domain

The Escherichia coli DNA repair enzyme MutY plays an important role in the prevention of DNA mutations by removing misincorporated adenine residues from 7,8-dihydro-8-oxo-2′-deoxyguanosine:2′-deoxyadenosine (OG:A) mispairs. The N-terminal domain of MutY (Stop 225, Met1–Lys225) has a sequence and structure that is characteristic of a superfamily of base excision repair glycosylases; however, MutY and its homologs contain a unique C-terminal domain. Previous studies have shown that the C-terminal domain confers specificity for OG:A substrates over G:A substrates and exhibits homology to the d(OG)TPase MutT, suggesting a role in OG recognition. In order to provide additional information on the importance of the C-terminal domain in damage recognition, we have investigated the kinetic properties of a form lacking this domain (Stop 225) under multiple- and single-turnover conditions. In addition, the interaction of Stop 225 with a series of non-cleavable substrate and product analogs was evaluated using gel retardation assays and footprinting experiments. Under multiple-turnover conditions Stop 225 exhibits biphasic kinetic behavior with both OG:A and G:A substrates, likely due to rate-limiting DNA product release. However, the rate of turnover of Stop 225 was increased 2-fold with OG:A substrates compared to the wild-type enzyme. In contrast, the intrinsic rate for adenine removal by Stop 225 from both G:A and OG:A substrates is significantly reduced (10- to 25-fold) compared to the wild-type. The affinity of Stop 225 for substrate analogs was dramatically reduced, as was the ability to discriminate between substrate analogs paired with OG over G. Interestingly, similar hydroxyl radical and DMS footprinting patterns are observed for Stop 225 and wild-type MutY bound to DNA duplexes containing OG opposite an abasic site mimic or a non-hydrogen bonding A analog, suggesting that similar regions of the DNA are contacted by both enzyme forms. Importantly, Stop 225 has a reduced ability to prevent DNA mutations in vivo. This implies that the reduced adenine glycosylase activity translates to a reduced capacity of Stop 225 to prevent DNA mutations in vivo.

SfiI endonuclease activity is strongly influenced by the non-specific sequence in the middle of its recognition site

The SfiI endonuclease cleaves DNA at the sequence GGCCNNNN↓NGGCC, where N is any base and ↓ is the point of cleavage. Proteins that recognise discontinuous sequences in DNA can be affected by the unspecified sequence between the specified base pairs of the target site. To examine whether this applies to SfiI, a series of DNA duplexes were made with identical sequences apart from discrete variations in the 5 bp spacer. The rates at which SfiI cleaved each duplex were measured under steady-state conditions: the steady-state rates were determined by the DNA cleavage step in the reaction pathway. SfiI cleaved some of these substrates at faster rates than other substrates. For example, the change in spacer sequence from AACAA to AAACA caused a 70-fold increase in reaction rate. In general, the extrapolated values for kcat and Km were both higher on substrates with inflexible spacers than those with flexible structures. The dinucleotide at the site of cleavage was largely immaterial. SfiI activity is thus highly dependent on conformational variations in the spacer DNA.

Toward rules for 1:1 polyamide:DNA recognition

Polyamides composed of four amino acids, imidazole (Im), pyrrole (Py), hydroxypyrrole (Hp), and β-alanine (β), are synthetic ligands that form highly stable complexes in the minor groove of DNA. Although specific pairing rules within the 2:1 motif can be used to distinguish the four Watson⋅Crick base pairs, a comparable recognition code for 1:1 polyamide:DNA complexes had not been described. To set a quantitative baseline for the field, the sequence specificities of Im, Py, Hp, and β for the four Watson⋅Crick base pairs were determined for two polyamides, Im-β-ImPy-β-Im-β-ImPy-β-Dp (1, for Im, Py, and β) and Im-β-ImHp-β-Im-β-ImPy-β-Dp (2, for Hp), in a 1:1 complex within the DNA sequence context 5′-AAAGAGAAGAG-3′. Im residues do not distinguish G,C from A,T but bind all four base pairs with high affinity. Py and β residues exhibit ≥10-fold preference for A,T over G,C base pairs. The Hp residue displays a unique preference for a single A⋅T base pair with an energetic penalty.

Involvement of Asn-293 in stereospecific agonist recognition and in activation of the beta 2-adrenergic receptor.

To investigate the molecular mechanism for stereospecific binding of agonists to beta 2-adrenergic receptors we used receptor models to identify potential binding sites for the beta-OH-group of the ligand, which defines the chiral center. Ser-165, located in transmembrane helix IV, and Asn-293, situated in the upper half of transmembrane helix VI, were identified as potential binding sites. Mutation of Ser-165 to Ala did not change the binding of either isoproterenol isomer as revealed after transient expression in human embryonic kidney (HEK)-293 cells. In contrast, a receptor mutant in which Asn-293 was replaced by Leu showed substantial loss of stereospecific isoproterenol binding. Adenylyl cyclase stimulation by this mutant after stable expression in CHO cells confirmed the substantial loss of stereospecificity for isoproterenol. In a series of agonists the loss of affinity in the Leu-293 mutant receptor was strongly correlated with the intrinsic activity of the compounds. Full agonists showed a 10-30-fold affinity loss, whereas partial agonists had almost the same affinity for both receptors. Stereospecific recognition of antagonists was unaltered in the Leu-293 mutant receptor. These data indicate a relationship between stereospecificity and intrinsic activity of agonists and suggest that Asn-293 is important for both properties of the agonist-receptor interaction.

Mutational analysis of NM23-H2/NDP kinase identifies the structural domains critical to recognition of a c-myc regulatory element.

NM23-H2, a presumed regulator of tumor metastasis in humans, is a hexameric protein with both enzymatic (NDP kinase) and regulatory (transcriptional activation) activity. While the structure and catalytic mechanisms have been well characterized, the mode of DNA binding is not known. We examined this latter function in a site-directed mutational study and identified residues and domains essential for the recognition of a c-myc regulatory sequence. Three amino acids, Arg-34, Asn-69, and Lys-135, were found among 30 possibilities to be critical for DNA binding. Two of these, Asn-69 and Lys-135, are not conserved between NM23 variants differing in DNA-binding potential, suggesting that DNA recognition resides partly in nonconserved amino acids. All three DNA-binding defective mutant proteins are active enzymatically and appear to be stable hexamers, suggesting that they perform at the level of DNA recognition and that separate functional domains exist for enzyme catalysis and DNA binding. In the context of the known crystal structure of NM23-H2, the DNA-binding residues are located within distinct structural motifs in the monomer, which are exposed to the surface near the 2-fold axis of adjacent subunits in the hexamer. These findings are explained by a model in which NM23-H2 binds DNA with a combinatorial surface consisting of the "outer" face of the dimer. Chemical crosslinking data support a dimeric DNA-binding mode by NM23-H2.

Differential recognition of the type I interferon receptor by interferons tau and alpha is responsible for their disparate cytotoxicities.

Interferon tau (IFN tau), originally identified as a pregnancy recognition hormone, is a type I interferon that is related to the various IFN alpha species (IFN alpha s). Ovine IFN tau has antiviral activity similar to that of human IFN alpha A on the Madin-Darby bovine kidney (MDBK) cell line and is equally effective in inhibiting cell proliferation. In this study, IFN tau was found to differ from IFN alpha A in that is was > 30-fold less toxic to MDBK cells at high concentrations. Excess IFN tau did not block the cytotoxicity of IFN alpha A on MDBK cells, suggesting that these two type I IFNs recognize the type I IFN receptor differently on these cells. In direct binding studies, 125I-IFN tau had a Kd of 3.90 x 10(-10) M for receptor on MDBK cells, whereas that of 125I-IFN alpha A was 4.45 x 10(-11) M. Consistent with the higher binding affinity, IFN alpha A was severalfold more effective than IFN tau in competitive binding against 125I-IFN tau to receptor on MDBK cells. Paradoxically, the two IFNs had similar specific antiviral activities on MDBK cells. However, maximal IFN antiviral activity required only fractional occupancy of receptors, whereas toxicity was associated with maximal receptor occupancy. Hence, IFN alpha A, with the higher binding affinity, was more toxic than IFN tau. The IFNs were similar in inducing the specific phosphorylation of the type I receptor-associated tyrosine kinase Tyk2, and the transcription factors Stat1 alpha and Stat2, suggesting that phosphorylation of these signal transduction proteins is not involved in the cellular toxicity associated with type I IFNs. Experiments using synthetic peptides suggest that differences in the interaction at the N terminal of IFN tau and IFN alpha with the type I receptor complex contribute significantly to differences in high-affinity equilibrium binding of these molecules. It is postulated that such a differential recognition of the receptor is responsible for the similar antiviral but different cytotoxic effects of these IFNs. Moreover, these data imply that receptors are "spare'' with respect to certain biological properties, and we speculate that IFNs may induce a concentration-dependent selective association of receptor subunits.

D1 cap region involved in the receptor recognition and neural cell survival activity of human ciliary neurotrophic factor.

Human ciliary neurotrophic factor (hCNTF), which promotes the cell survival and differentiation of motor and other neurons, is a protein belonging structurally to the alpha-helical cytokine family. hCNTF was subjected to three-dimensional structure modeling and site-directed mutagenesis to analyze its structure-function relationship. The replacement of Lys-155 with any other amino acid residue resulted in abolishment of neural cell survival activity, and some of the Glu-153 mutant proteins had 5- to 10-fold higher biological activity. The D1 cap region (around the boundary between the CD loop and helix D) of hCNTF, including both Glu-153 and Lys-155, was shown to play a key role in the biological activity of hCNTF as one of the putative receptor-recognition sites. In this article, the D1 cap region of the 4-helix-bundle proteins is proposed to be important in receptor recognition and biological activity common to alpha-helical cytokine proteins reactive with gp130, a component protein of the receptors.

Alternate protein frameworks for molecular recognition.

In an effort to determine whether proteins with structures other than the immunoglobulin fold can be used to mimic the ligand binding properties of antibodies, we generated a library from the four-helix bundle protein cytochrome b562 in which the two loops were randomized. Panning of this library against the bovine serum albumin (BSA) conjugate of N-methyl-p-nitrobenzylamine derivative 1 by phage display methods yielded cytochromes in which residues Trp-20, Arg-21, and Ser-22 in loop A and Arg-83 and Trp-84 in loop B were conserved. The individual mutants, which fold into native-like structure, bind selectively to the BSA-1 conjugate with micromolar dissociation constants (Kd), in comparison to a monoclonal antibody that binds selectively to 1 with a Kd of 290 nM. These and other antibody-like receptors may prove useful as therapeutic agents or as reagents for both intra- and extracellular studies.