Peter Terpstra

Prediction of the occurrence of the ADP-binding βαβ-fold in proteins, using an amino acid sequence fingerprint

Abstract An amino acid sequence “fingerprint” has been derived that can be used to test if a particular sequence will fold into aβαβ-unit with ADP-binding properties. It was deduced from a careful analysis of the known three-dimensional structures of ADP-binding βαβ-folds. This fingerprint is in fact a set of 11 rules describing the type of amino acid that should occur at a specific position in a peptide fragment. The total length of this fingerprint varies between 29 and 31 residues. By checking against all possible sequences in a database, it appeared that every peptide, which exactly follows this fingerprint, does indeed fold into an ADP-binding βαβ-unit.

Poor replication of candidate genes for major depressive disorder using genome-wide association data.

Data from the Genetic Association Information Network (GAIN) genome-wide association study (GWAS) in major depressive disorder (MDD) were used to explore previously reported candidate gene and single-nucleotide polymorphism (SNP) associations in MDD. A systematic literature search of candidate genes associated with MDD in case–control studies was performed before the results of the GAIN MDD study became available. Measured and imputed candidate SNPs and genes were tested in the GAIN MDD study encompassing 1738 cases and 1802 controls. Imputation was used to increase the number of SNPs from the GWAS and to improve coverage of SNPs in the candidate genes selected. Tests were carried out for individual SNPs and the entire gene using different statistical approaches, with permutation analysis as the final arbiter. In all, 78 papers reporting on 57 genes were identified, from which 92 SNPs could be mapped. In the GAIN MDD study, two SNPs were associated with MDD: C5orf20 (rs12520799; P=0.038; odds ratio (OR) AT=1.10, 95% CI 0.95–1.29; OR TT=1.21, 95% confidence interval (CI) 1.01–1.47) and NPY (rs16139; P=0.034; OR C allele=0.73, 95% CI 0.55–0.97), constituting a direct replication of previously identified SNPs. At the gene level, TNF (rs76917; OR T=1.35, 95% CI 1.13–1.63; P=0.0034) was identified as the only gene for which the association with MDD remained significant after correction for multiple testing. For SLC6A2 (norepinephrine transporter (NET)) significantly more SNPs (19 out of 100; P=0.039) than expected were associated while accounting for the linkage disequilibrium (LD) structure. Thus, we found support for involvement in MDD for only four genes. However, given the number of candidate SNPs and genes that were tested, even these significant may well be false positives. The poor replication may point to publication bias and false-positive findings in previous candidate gene studies, and may also be related to heterogeneity of the MDD phenotype as well as contextual genetic or environmental factors.

Rubredoxin reductase of Pseudomonas oleovorans: Structural relationship to other flavoprotein oxidoreductases based on one NAD and two FAD fingerprints☆

The oxidation of alkanes to alkanols by Pseudomonas oleovorans involves a three-component enzyme system: alkane hydroxylase, rubredoxin and rubredoxin reductase. Alkane hydroxylase and rubredoxin are encoded by the alkBFGHJKL operon, while previous studies indicated that rubredoxin reductase is most likely encoded on the second alk cluster: the alkST operon. In this study we show that alkT encodes the 41 x 10(3) Mr rubredoxin reductase, on the basis of a comparison of the expected amino acid composition of AlkT and the previously established amino acid composition of the purified rubredoxin reductase. The alkT sequence revealed significant similarities between AlkT and several NAD(P)H and FAD-containing reductases and dehydrogenases. All of these enzymes contain two ADP binding sites, which can be recognized by a common beta alpha beta-fold or fingerprint, derived from known structures of cofactor binding enzymes. By means of this amino acid fingerprint we were able to determine that one ADP binding site in rubredoxin reductase (AlkT) is located at the N terminus and is involved in FAD binding, while the second site is located in the middle of the sequence and is involved in the binding of NAD or NADP. In addition, we derived from the sequences of FAD binding reductases a second amino acid fingerprint for FAD binding, and we used this fingerprint to identify a third amino acid sequence in AlkT near the carboxy terminus for binding of the flavin moiety of FAD. On the basis of the known architecture and relative spatial orientations of the NAD and FAD binding sites in related dehydrogenases, a model for part of the tertiary structure of AlkT was developed.

Sequence analysis and molecular characterization of the temperate lactococcal bacteriophage r1t

The temperate lactococcal bacteriophage r1t was isolated from its lysogenic host and its genome was subjected to nucleotide sequence analysis. The linear r1t genome is composed of 33 350 bp and was shown to possess 3′ staggered cohesive ends. Fifty open reading frames (ORFs) were identified which are, probably, organized in a life‐cycle‐specific manner. Nucleotide sequence comparisons, N‐terminal amino acid sequencing and functional analyses enabled the assignment of possible functions to a number of DNA sequences and ORFs. In this way, ORFs specifying regulatory proteins, proteins involved in DNA replication, structural proteins, a holin, a lysin, an integrase, and a dUTPase were putatively identified. One ORF seems to be contained within a self‐splicing group I intron. In addition, the bacteriophage att site required for site‐specific integration into the host chromosome was determined.

Sequence Polymorphisms Cause Many False cis eQTLs

Many investigations have reported the successful mapping of quantitative trait loci (QTLs) for gene expression phenotypes (eQTLs). Local eQTLs, where expression phenotypes map to the genes themselves, are of especially great interest, because they are direct candidates for previously mapped physiological QTLs. Here we show that many mapped local eQTLs in genetical genomics experiments do not reflect actual expression differences caused by sequence polymorphisms in cis-acting factors changing mRNA levels. Instead they indicate hybridization differences caused by sequence polymorphisms in the mRNA region that is targeted by the microarray probes. Many such polymorphisms can be detected by a sensitive and novel statistical approach that takes the individual probe signals into account. Applying this approach to recent mouse and human eQTL data, we demonstrate that indeed many local eQTLs are falsely reported as “cis-acting” or “cis” and can be successfully detected and eliminated with this approach.

Identification and organization of carbon dioxide fixation genes in Xanthobacter flavus H4-14

SummaryThe genes encoding the large (cfxL) and small (cfxS) subunits of ribulose-1,5-bisphosphate carboxylase (RuBisC/O) from Xanthobacter favus H4-14 were identified and characterized. The RuBisC/O genes are separated by 11 by and cotranscribed in Escherichia coli from the lac promoter in the order cfxLS. Primer extension and R-loop experiments with RNA isolated from autotrophically grown X. favus H4-14 showed that transcription of cfxL and cfxS initiated 22 by upstream from cfxL and resulted in a mRNA of at least 2.3 kb. DNA sequence analysis identified the start of an open reading frame transcribed divergently from cfxL, and displaying significant similarities with genes belonging to the lysR family of transcriptional activators. Downstream from cfxS an additional open reading frame was identified with unknown function. Expression studies showed that the genes encoding fructosebisphosphatase (cfxF) and phosphoribulokinase (cfxP) are located downstream from cfxLS. The cfxF and cfxP genes are cotranscribed in the same direction as cfxLS in the order cfxFP.

CLONING AND SEQUENCING OF PSEUDOMONAS GENES DETERMINING SODIUM DODECYL-SULFATE BIODEGRADATION

The nucleotide sequences of two genes involved in sodium dodecyl sulfate (SDS) degradation, by Pseudomonas, have been determined. One of these, sdsA, codes for an alkyl sulfatase (58,957 Da) and has similarity (31.8% identity over a 201-amino acid stretch) to the N terminus of a predicted protein of unknown function from Mycobacterium tuberculosis. The other gene, sdsB, codes for a positive activator protein (33,600 Da) that has extensive similarity with the lysR family of helix-turn-helix DNA-binding activator proteins.

The restriction fragment map of rat-liver mitochondrial DNA: a reconsideration.

Abstract 1. Rat-liver mitochondrial DNA (mtDNA) contains at least 8 cleavage sites for the restriction endonuclease Eco RI, 6 for the restriction endonuclease Hind III, 2 for the restriction endonuclease Bam HI and 11 for the restriction endonuclease Hap II. 2. The physical map of the restriction fragments of Eco RI, Hind III, Bam HI and Hap II is constructed on the basis of: (a) the analysis of partially restricted fragments; (b) analysis of the double digests of total mtDNA; (c) the digestion of isolated restriction fragments with other restriction endonucleases; (d) the identification of fragments of complete single and double digestions and of partially digested fragments containing the base sequences complementary to the 12-S and 16-S RNAs of rat-liver mitochondrial ribosomes. 3. The genes for the ribosomal RNAs are shown to be closely linked. This result differs from data previously reported (Saccone, C., Pepe, G., Cantatore, P., Terpstra, P. and Kroon, A.M. (1976) in The Genetic Function of Mitochondrial DNA, pp. 27–36, Elsevier/North-Holland Biomedical Press, Amsterdam). 4. The origin of replication (D-loop) is localized in the vicinity of the small ribosomal RNA gene and the direction of replication is distant from this gene. 5. The mitochondrial tRNA genes are scattered over the genome as in other animal mtDNAs. The approximate minimal number of tRNA genes is 16–20. 6. We concluded previously that the Eco RI restriction fragments A and D are not adjacent and failed to show the overlap of the 16 S rRNA gene for the Eco RI fragment D and Hind III fragment A. This misinterpretation was due to the fact that the two smallest Eco RI fragments could not be detected with the methods applied and to the lower specific radioactivity of the ribosomal RNAs used in the first series of hybridization experiments.

Homology of Drosophila yolk proteins and the triacylglycerol lipase family

Drosophila yolk proteins consist of a set of related proteins of 50,000 Mr. They are derived from slightly larger precursors by cleavage of a signal peptide. In this respect, they differ from the yolk proteins of other insects which are proteolytic fragments of precursors of 200,000 Mr or larger, termed vitellogenins and probably homologous to the vitellogenins of other egg-laying species. We report here a comparative amino acid analysis demonstrating that the Drosophila yolk proteins are non-homologous to the vitellogenin group of yolk proteins, but surprisingly are related to the triacylglycerol lipase family.

A statistical multiprobe model for analyzing cis and trans genes in genetical genomics experiments with short-oligonucleotide arrays

Short-oligonucleotide arrays typically contain multiple probes per gene. In genetical genomics applications a statistical model for the individual probe signals can help in separating “true” differential mRNA expression from “ghost” effects caused by polymorphisms, misdesigned probes, and batch effects. It can also help in detecting alternative splicing, start, or termination.