Eugene Melamud

SNPs3D: Candidate gene and SNP selection for association studies

BackgroundThe relationship between disease susceptibility and genetic variation is complex, and many different types of data are relevant. We describe a web resource and database that provides and integrates as much information as possible on disease/gene relationships at the molecular level.DescriptionThe resource http://www.SNPs3D.org has three primary modules. One module identifies which genes are candidates for involvement in a specified disease. A second module provides information about the relationships between sets of candidate genes. The third module analyzes the likely impact of non-synonymous SNPs on protein function. Disease/candidate gene relationships and gene-gene relationships are derived from the literature using simple but effective text profiling. SNP/protein function relationships are derived by two methods, one using principles of protein structure and stability, the other based on sequence conservation. Entries for each gene include a number of links to other data, such as expression profiles, pathway context, mouse knockout information and papers. Gene-gene interactions are presented in an interactive graphical interface, providing rapid access to the underlying information, as well as convenient navigation through the network. Use of the resource is illustrated with aspects of the inflammatory response and hypertension.ConclusionThe combination of SNP impact analysis, a knowledge based network of gene relationships and candidate genes, and access to a wide range of data and literature allow a user to quickly assimilate available information, and so develop models of gene-pathway-disease interaction.

Completeness in structural genomics

Structural genomics has the goal of obtaining useful, three-dimensional models of all proteins by a combination of experimental structure determination and comparative model building. We evaluate different strategies for optimizing information return on effort. The strategy that maximizes structural coverage requires about seven times fewer structure determinations compared with the strategy in which targets are selected at random. With a choice of reasonable model quality and the goal of 90% coverage, we extrapolate the estimate of the total effort of structural genomics. It would take ∼16,000 carefully selected structure determinations to construct useful atomic models for the vast majority of all proteins. In practice, unless there is global coordination of target selection, the total effort will likely increase by a factor of three. The task can be accomplished within a decade provided that selection of targets is highly coordinated and significant funding is available.

Stochastic noise in splicing machinery

The number of known alternative human isoforms has been increasing steadily with the amount of available transcription data. To date, over 100 000 isoforms have been detected in EST libraries, and at least 75% of human genes have at least one alternative isoform. In this paper, we propose that most alternative splicing events are the result of noise in the splicing process. We show that the number of isoforms and their abundance can be predicted by a simple stochastic noise model that takes into account two factors: the number of introns in a gene and the expression level of a gene. The results strongly support the hypothesis that most alternative splicing is a consequence of stochastic noise in the splicing machinery, and has no functional significance. The results are also consistent with error rates tuned to ensure that an adequate level of functional product is produced and to reduce the toxic effect of accumulation of misfolding proteins. Based on simulation of sampling of virtual cDNA libraries, we estimate that error rates range from 1 to 10% depending on the number of introns and the expression level of a gene.

Biochemical Characterization of the Methanothermobacter thermautotrophicus Minichromosome Maintenance (MCM) Helicase N-terminal Domains

Minichromosome maintenance helicases are ring-shaped complexes that play an essential role in archaeal and eukaryal DNA replication by separating the two strands of chromosomal DNA to provide the single-stranded substrate for the replicative polymerases. For the archaeal protein it was shown that the N-terminal portion of the protein, which is composed of domains A, B, and C, is involved in multimer formation and single-stranded DNA binding and may also play a role in regulating the helicase activity. Here, a detailed biochemical characterization of the N-terminal region of the Methanothermobacter thermautotrophicus minichromosome maintenance helicase is described. Using biochemical and biophysical analyses it is shown that domain C of the N-terminal portion, located adjacent to the helicase catalytic domains, is required for protein multimerization and that domain B is the main contact region with single-stranded DNA. It is also shown that although oligomerization is not essential for single-stranded DNA binding and ATPase activity, the presence of domain C is essential for helicase activity.

From fold to function.

A number of recent advances have been made in deriving function information from protein structure. A fold relationship to an already characterized protein will often allow general information about function to be deduced. More detailed information can be obtained using sequence relationships to already studied proteins. Methods of deducing function directly from structure, without the use of evolutionary relationships, are developing rapidly. All such methods may be used with models of protein structure, rather than with experimentally determined ones, but model accuracy imposes limitations. The rapid expansion of the structural genomics field has created a new urgency for improved methods of structure-based annotation of function.

Characterization, Kinetics, and Crystal Structures of Fructose-1,6-bisphosphate Aldolase from the Human Parasite, Giardia lamblia

Class I and class II fructose-1,6-bisphosphate aldolases (FBPA), glycolytic pathway enzymes, exhibit no amino acid sequence homology and utilize two different catalytic mechanisms. The mammalian class I FBPA employs a Schiff base mechanism, whereas the human parasitic protozoan Giardia lamblia class II FBPA is a zinc-dependent enzyme. In this study, we have explored the potential exploitation of the Giardia FBPA as a drug target. First, synthesis of FBPA was demonstrated in Giardia trophozoites by using an antibody-based fluorescence assay. Second, inhibition of FBPA gene transcription in Giardia trophozoites suggested that the enzyme is necessary for the survival of the organism under optimal laboratory growth conditions. Third, two crystal structures of FBPA in complex with the transition state analog phosphoglycolohydroxamate (PGH) show that the enzyme is homodimeric and that its active site contains a zinc ion. In one crystal form, each subunit contains PGH, which is coordinated to the zinc ion through the hydroxamic acid hydroxyl and carbonyl oxygen atoms. The second crystal form contains PGH only in one subunit and the active site of the second subunit is unoccupied. Inspection of the two states of the enzyme revealed that it undergoes a conformational transition upon ligand binding. The enzyme cleaves d-fructose-1,6-bisphosphate but not d-tagatose-1,6-bisphosphate, which is a tight binding competitive inhibitor. The essential role of the active site residue Asp-83 in catalysis was demonstrated by amino acid replacement. Determinants of catalysis and substrate recognition, derived from comparison of the G. lamblia FBPA structure with Escherichia coli FBPA and with a closely related enzyme, E. coli tagatose-1,6-bisphosphate aldolase (TBPA), are described.

Structural implication of splicing stochastics

Even though nearly every human gene has at least one alternative splice form, very little is so far known about the structure and function of resulting protein products. It is becoming increasingly clear that a significant fraction of all isoforms are products of noisy selection of splice sites and thus contribute little to actual functional diversity, and may potentially be deleterious. In this study, we examine the impact of alternative splicing on protein sequence and structure in three datasets: alternative splicing events conserved across multiple species, alternative splicing events in genes that are strongly linked to disease and all observed alternative splicing events. We find that the vast majority of all alternative isoforms result in unstable protein conformations. In contrast to that, the small subset of isoforms conserved across species tends to maintain protein structural integrity to a greater extent. Alternative splicing in disease-associated genes produces unstable structures just as frequently as all other genes, indicating that selection to reduce the effects of alternative splicing on this set is not especially pronounced. Overall, the properties of alternative spliced proteins are consistent with the outcome of noisy selection of splice sites by splicing machinery.

The HI0073/HI0074 protein pair from Haemophilus influenzae is a member of a new nucleotidyltransferase family: Structure, sequence analyses, and solution studies

The crystal structure of HI0074 from Haemophilus influenzae, a protein of unknown function, has been determined at a resolution of 2.4 Å. The molecules form an up–down, four‐helix bundle, and associate into homodimers. The fold is most closely related to the substrate‐binding domain of KNTase, yet the amino acid sequences of the two proteins exhibit no significant homology. Sequence analyses of completely and incompletely sequenced genomes reveal that the two adjacent genes, HI0074 and HI0073, and their close relatives comprise a new family of nucleotidyltransferases, with 15 members at the time of writing. The analyses also indicate that this is one of eight families of a large nucleotidyltransferase superfamily, whose members were identified based on the proximity of the nucleotide‐ and substrate‐binding domains on the respective genomes. Both HI0073 and HI0074 were annotated “hypothetical” in the original genome sequencing publication. HI0073 was cloned, expressed, and purified, and was shown to form a complex with HI0074 by polyacrylamide gel electrophoresis under nondenaturing conditions, analytic size exclusion chromatography, and dynamic light scattering. Double‐ and single‐stranded DNA binding assays showed no evidence of DNA binding to HI0074 or to HI0073/HI0074 complex despite the suggestive shape of the putative binding cleft formed by the HI0074 dimer. Proteins 2003;50:249–260.

Research Article: Cloning and partial characterization of the Methanococcoides burtonii minichromosome maintenance (MCM) helicase

Abstract Minichromosome maintenance (MCM) helicases play an essential role in chromosomal DNA replication in archaea and eukarya. Here the isolation and partial characterization of the MCM homologue from the psychophilic archaeaon Methanococcoides burtonii is described. The protein was expressed and purified to near homogeneity from Escherichia coli and, like MCM homologues from other archaeons, the protein was shown to form hexamers.