Increasing synthetic performance of penicillin G acylase from Bacillus megaterium by site-directed mutagenesis.

Site-directed mutagenesis based on predicted modeled structure of pencillin G acylase from Bacillus megaterium (BmPGA) was followed to increase its performance in the kinetically controlled synthesis of cephalexin with high reactant concentrations of 133 mM 7-amino-desaceto-xycephalosporanic acid (7-ADCA) and 267 mM D: -phenylglycine amide (D-PGA). We directed changes in amino acid residues to positions close to the active site that were expected to affect the catalytic performance of penicillin acylase: alpha Y144, alpha F145, and beta V24. Alpha F145 was mutated into tyrosine, alanine, and leucine. Alpha Y144 and beta V24 were mutated into arginine and phenylalanine, respectively. The S/H ratios of three mutants, BmPGAalpha144R, BmPGAbeta24F, and BmPGAbeta24F+alpha144R, were up to 1.3-3.0 times higher values. Compared to the wild-type BmPGA, BmPGAbeta24F+alpha144R showed superior potential of the synthetic performance, allowing the accumulation of up to twofold more cephalexin at significantly higher conversion rates.

Identification and categorization of horizontally transferred genes in prokaryotic genomes.

Horizontal gene transfer (HGT), a process through which genomes acquire genetic materials from distantly related organisms, is believed to be one of the major forces in prokaryotic genome evolution. However, systematic investigation is still scarce to clarify two basic issues about HGT: (1) what types of genes are transferred; and (2) what influence HGT events over the organization and evolution of biological pathways. Genome-scale investigations of these two issues will advance the systematical understanding of HGT in the context of prokaryotic genome evolution. Having investigated 82 genomes, we constructed an HGT database across broad evolutionary timescales. We identified four function categories containing a high proportion of horizontally transferred genes: cell envelope, energy metabolism, regulatory functions, and transport/binding proteins. Such biased function distribution indicates that HGT is not completely random; instead, it is under high selective pressure, required by function restraints in organisms. Furthermore, we mapped the transferred genes onto the connectivity structure map of organism-specific pathways listed in Kyoto Encyclopedia of Genes and Genomes (KEGG). Our results suggest that recruitment of transferred genes into pathways is also selectively constrained because of the tuned interaction between original pathway members. Pathway organization structures still conserve well through evolution even with the recruitment of horizontally transferred genes. Interestingly, in pathways whose organization were significantly affected by HGT events, the operon-like arrangement of transferred genes was found to be prevalent. Such results suggest that operon plays an essential and directional role in the integration of alien genes into pathways.

Searching for potential drug targets in two-component and phosphorelay signal-transduction systems using three-dimensional cluster analysis.

Two-component and phosphorelay signal transduction systems are central components in the virulence and antimicrobial resistance responses of a number of bacterial and fungal pathogens; in some cases, these systems are essential for bacterial growth and viability. Herein, we analyze in detail the conserved surface residue clusters in the phosphotransferase domain of histidine kinases and the regulatory domain of response regulators by using complex structure-based three-dimensional cluster analysis. We also investigate the protein-protein interactions that these residue clusters participate in. The Spo0B-Spo0F complex structure was used as the reference structure, and the multiple aligned sequences of phosphotransferases and response regulators were paired correspondingly. The results show that a contiguous conserved residue cluster is formed around the active site, which crosses the interface of histidine kinases and response regulators. The conserved residue clusters of phosphotransferase and the regulatory domains are directly involved in the functional implementation of two-component signal transduction systems and are good targets for the development of novel antimicrobial agents.

Rational redesign of inhibitors of furin/kexin processing proteases by electrostatic mutations.

AIM: To model the three-dimensional structure and investigate the interaction mechanism of the proprotein convertase furin/kexin and their inhibitors (eglin c mutants). METHODS: The three-dimensional complex structures of furin/kexin with its inhibitors, eglin c mutants, were generated by modeller program using the newly published X-ray crystallographical structures of mouse furin and yeast kexin as templates. The electrostatic interaction energy of each complex was calculated and the results were compared with the experimentally determined inhibition constants to find the correlation between them. RESULTS: High quality models of furin/kexin-eglin c mutants were obtained and used for calculation of the electrostatic interaction energies between the proteases and their inhibitors. The calculated electrostatic energies of interaction showed a linear correlation to the experimental inhibition constants. CONCLUSION: The modeled structures give good explanations of the specificity of eglin c mutants to furin/kexin. The electrostatic interactions play important roles in inhibitory activity of eglin c mutants to furin/kexin. The results presented here provided quantitative structural and functional information concerning the role of the charge-charge interactions in the binding of furin/kexin and their inhibitors.

Learning module networks from genome-wide location and expression data.

We develop a systematic algorithm for discovering network of regulatory modules, which identifies regulatory modules and their regulation program by integrating genome-wide location and expression data. Unlike previous approaches [Eisen, M.B., Spellman, P.T., Brown, P.O. and Botstein, D. (1998) Proc. Natl. Acad. Sci. USA 95, 14863-14868; Tavazoie, S., Hughes, J.D., Campbell, M.J., Cho, R.J. and Church, G.M. (1999) Nat. Genet. 22, 281-285; Ihmels, J., Friedlander, G., Bergmann, S., Sarig, O., Ziv, Y. and Barkai, N. (2002) Nat. Genet. 31, 370-377; Segal, E., Shapira, M., Regev, A., Pe'er, D., Botstein, D., Koller, D. and Friedman, N. (2003) Nat. Genet. 34, 166-176] that relied primarily on gene expression data, our algorithm regards the regulator binding data as prior knowledge that provide direct evidence of physical regulatory interactions. We applied the method to a Saccharomyces cerevisiae genome-wide location data [Lee, T.I., Rinaldi, N.J., Robert, F., Odom, D.T., Bar-Joseph, Z., Gerber, G.K., Hannett, N.M., Harbison, C.T., Thompson, C.M., Simon, I., Zeitlinger, J., Jennings, E.G., Murray, H.L. Gordon, D.B., Ren, B., Wyrick, J.J., Tagne, J.B., Volkert, T.L., Fraenkel, E., Gifford, D.K. and Young, R.A. (2002) Science 298, 799-804] for 106 DNA-binding transcription factors and 250 gene expression experiments under the conditions from the cell cycle to responses to various stress conditions. The results show that our method is able to identify functionally coherent modules and their proper regulators. Supplementary materials are available at http://compbio.sibnet.org/projects/module-network/.

MedBlast: searching articles related to a biological sequence.

UNLABELLED: In the genomic era, researchers often want to know more information about a biological sequence by retrieving its related articles. However, there is no available tool yet to achieve conveniently this goal. Here we developed a new literature-mining tool MedBlast, which uses natural language processing techniques, to retrieve the related articles of a given sequence. An online server of this program is also provided. AVAILABILITY: Both online server and the program are available freely at http://medblast.sibsnet.org

Identification of the key amino acids of glial cell line-derived neurotrophic factor family receptor alpha1 involved in its biological function.

Glial cell line-derived neurotrophic factor (GDNF) plays a critical role in neurodevelopment and survival of midbrain dopaminergic and spinal motor neurons in vitro and in vivo. The biological actions of GDNF are mediated by a two-receptor complex consisting of a glycosylphosphatidylinositol-linked cell surface molecule, the GDNF family receptor alpha1 (GFRalpha1), and receptor protein tyrosine kinase Ret. Although structural analysis of GDNF has been extensively examined, less is known about the structural basis of GFRalpha1 function. In this study, based on evolutionary trace method and relative solvent accessibility prediction of residues, a set of trace residues that are solvent-accessible was selected for site-directed mutagenesis. A series of GFRalpha1 mutations was made, and PC12 cell lines stably expressing different GFRalpha1 mutants were generated. According to the survival and differentiation responses of these stable PC12 cells upon GDNF stimulation and the GDNF-GFRalpha1-Ret interaction assay, residues 152NN153, Arg259, and 316SNS318 in the GFRalpha1 central region were found to be critical for GFRalpha1 binding to GDNF and eliciting downstream signal transduction. The single mutation R259A in the GFRalpha1 molecule simultaneously lost its binding ability to GDNF and Ret. However N152A/N153A or S316A/N317A/S318A mutation in the GFRalpha1 molecule still retained the ability to bind with Ret. These findings suggest that distinct structural elements in GFRalpha1 may be involved in binding to GDNF and Ret.

[Estimating coarse gene networks from yeast gene expression time series].

Gene networks is the collection of gene-gene regulatory relations at the expression level. In this study, a combined approach of the linear transcriptional modeling, identification of promoter elements and gene co-expression clustering is developed to decipher yeast gene networks from expression time series. The cell must reorganize the genomic expression to programs required for growth and survival in each environment. The expression of many genes is regulated by environmental stress. The products of many genes that induced in the environmental stress are involved in metabolism of carbohydrates, structural repairs and even sporulation. Interestingly, it is identified that transcription factors Mcm1 and Dal82 matched their binding sites TT[bond]CC[triple bond]T[double bond]GGAAA and TGAAAAWTTT in cell cycle progression and environmental stress response, respectively. These conclusions agree with the known observations. The results indicate that the approach may be useful for modeling gene networks from microarray data.

[A novel approach for peptide identification by tandem mass spectrometry].

High throughput scoring algorithms that are used to find the match of a tandem mass spectrum to a predicted mass spectrum of a peptide within a database have been applied in shotgun proteomics. However, these algorithms could produce a significant number of incorrect peptide identifications. Here a novel approach was developed to scoring tandem mass spectra against a peptide database, in which fragment ion probabilities, number of enzymatic termini of candidate peptides, matching quality and match pattern between experimental and theoretical spectrum were considered. Benchmarking the novel scorer on a large set of experimental MS/MS spectra, it is demonstrated that PepSearch performs significantly better than the widely used software SEQUEST. The PepSearch software is available at http://compbio.sibsnet.org/projects/pepsearch.

Detecting pathogenicity islands and anomalous gene clusters by iterative discriminant analysis.

We present a simple method to detect pathogenicity islands and anomalous gene clusters in bacterial genomes. The method uses iterative discriminant analysis to define genomic regions that deviate most from the rest of the genome in three compositional criteria: G+C content, dinucleotide frequency and codon usage. Using this method, we identify many virulence-related gene islands, e.g. encoding protein secretion systems, adhesins, toxins, and other anomalous gene clusters, such as prophages. The program and the whole dataset, including the catalogs of genes in the detected anomalous segments, are publicly available at http://compbio.sibsnet.org/projects/pai-ida/. This program can be used in searching for virulence-related factors in newly sequenced bacterial genomes.

[Deducing functional epitopes for GDNF proteins and its specific GFRalpha co-receptors using phylogenetic approach].

Glial cell line-derived neurotrophic factor (GDNF) has received much attention as potential therapeutic agent for the treatment of neurodegenerative diseases. It will be very important to discover the molecular mechanism of this factor and its specific GFRalpha co-receptor. Based on the principle of molecular evolution that site-specific functional importance is relevant to the pressure it undergoes under natural selection, evolutionary trace method was used to identify the functional epitopes in GDNF and GFRalpha families. Some trace residues had been proved to be important in ligand-receptor binding, especially in rat GFRalpha1, where alanine scanning mutagenesis confirmed that sites N(152)N(153), R(259), S(316)N(317)S(318) and Q(247)D(248)S(249) were critical for GFRalpha1 binding to GDNF or Ret and thus affected the formation of GDNF-GFRalpha1-Ret complex.

A structure-function analysis of glial cell-line-derived neurotrophic factor receptor alpha1.

The GFRalpha1 cDNA was amplified by RT-PCR from fetal rat hippocampus. The soluble recombinant GFRalpha1 and its mutants were obtained from an Escherichia coli expression system. The biological activity of soluble GFRalpha1 and its mutants were evaluated in PC12 cells. The results suggest that the central domain of GFRalpha1 is a crucial determinant for ligand binding. This established a solid basis for further study to find the key amino acid mediating the binding of GDNF and GFRalpha1.

[A evolutionary approach to identification of orthologous relationship across proteomes].

How to identify the true orthologous and paralogous relationships among protein families is still a key problem in genome annotation and comparative protemics. Here, a evolutionary approach to ascertainment of the orthologous relationships across the genomes is developed. Forty-four cases of protein families are used in the test for the evolutionary approach. Compared with the method of COG (cluster of orthologous groups of proteins), this approach can generally identify the orthologous relationships and accurately predict the genome function.

Rational Redesign of Inhibitors of Furin/kexin Processing Proteases.

Furin/kexin processing proteases catalyze the proteolysis of large protein precursors involved in many biological processes, such as zymogen activation, peptide hormone synthesis, viral protein processing and receptor maturation, making them potential targets for therapeutic agents. Herein, homology modeling and weighted evolutionary tracing were combined to investigate the interactionmechanism of furin/kex2 with eglin C mutants. The model structures showed that there were many acidic residues in the furin (kex2) binding interface, contributing to specificity for multiple basic residues of their corresponding substrates or inhibitors. Besides, some rational explanations were presented for the different inhibitor/substrate specificity of the furin/kexin members by combining the model structures with results of evolutionary tracing. Based on these analyses,an attempt was made to rationally redesign the eglin C by interface engineering with heterogeneous self-consistent ensemble optimization to improve its inhibitory specificity on furin/kex2. With the model complex structures of furin/kex2 and eglin C variants as structural templates, the P(1), P(2) and P(4) of eglin C were redesigned, respectively. The design results show that both furin and kex2 favored basic residues at P(1), P(2) and P(4) in eglin C, in good agreement with the experimental data. The detection of many specific residues in S' part of furin/kexin sequences made possible designing inhibitors with high specific binding to furin and kex2, respectively. As for furin, the best inhibitor designed was eglin C-P(2)'Glu-P(3)'Asp-P(4)'Arg (only these three positions were shown), while the best eglin C variant for kex2 designed was P(2)'Arg-P(3)'Arg-P(4)'Glu. The structures show that furin and kex2 form distinct interactions with these two eglin C variants. Herein, a strategy was proposed that combine homology modeling, evolutionary tracing and rational interface redesign to investigate enzyme-inhibitor interactions and inhibitor engineering. This computational design gives some rational guidance to further experimental inhibitor engineering.

Delineation of Continuous Domains in Proteins by Differences of Free Energy.

Domain is a protein architecture under proteins' tertiary structure,which can be identified in most of proteins. Different combinations of domains lead to the formation of diverse tertiary structures with diverse function for proteins. The delineation of domains for a protein is important not only conceptually but also practically. Unfortunately, up to now there is not an ideal means to achieve that. This paper proposes a method for domain delineation based on the maximum refolding free energy. The criteria are more or less objective. By using this method, 50 proteins are analyzed. The boundaries for most proteins agree with the data reported in literature. There are a few examples that seem more reasonable, although they are not identical with those in literature.

The Spot Fitting and Expression Quantification for the Nylon Filter cDNA Microarray.

High-density cDNA microarray is important in monitoring large scalegenomic expression profile. With radiolabeled cDNA sample, nylon filter cDNA microarrays have high sensitivity and wide linear range of hybridization signals, so it can be employed to detect many important low-abundance cDNAs. However, the nylon filter tends to be randomly contaminated, and the array spots are prone to dispersion and saturation, resulting in inaccurate expression value. Based on aphysical model, we address a novel treatment with accurate positioning, complement of saturation and correction of interference. Thus the accurate quantificati on of expression can be extracted with remarkably improved reproducibility.

De Novo Interpretation of MS/MS Spectra and Protein Identification via Database Searching.

Peptide sequencing via tandem mass spectrometry(MS/MS)is one of the most powerful tools in proteomics to identify proteins. A new algorithm was developed for de novo interpretation of MS/MS spectra using graph theory and dynamic alignment between real spectra and theoretical spectra. The trustworthy peptides from de novo interpretation were used in protein identification via database searching. A high throughput statistical analysis of SwissProt and TrEMBL protein databases showed that it's enough to identify a protein in database with three sequence tags of four amino acid residues, two sequence tags of five amino acid residues or one sequence tag of eight amino acid residues.

Conservation and Evolution of the "Core Apoptotic Engine" Lesson from the Genome Comparison of Drosophila.

Genome comparison is the main approach to deduce regulatory network from genome sequence. Apoptotic network is one kind of typical regulatory networks. EGL1, CED3, CED4, CED9 and their homologous proteins play essential roles in apoptosis of C.elegans and mammals, and were regarded as the components of the"core apoptotic engine". But in fruit fly, Drosophila melanogaster, this network is incomplete. A series of bioinformatic analyses found the lost chains of "core apoptotic engine" by discovering two homologues of BCL2/CED9 and one of EGL1 in fruit fly genome sequences. These findings proved that the "core apoptotic engine" is indeed widely conserved among multicellular organisms and the evolutionary complexity of this network of Drosophila is between that of C.elegans and mammals.

Reconstruction of ABC Transporter Pathways of Archaea and Comparison of Their Genomes.

Reconstruction and comparison of metabolic pathway and regulatory network is an advanced task in genome function prediction. In this study, many bioinformatic tools were employed to reconstruct all ABC transporter pathways and predict their functional features of an archaeon, Pyrococcus abyssi, on genome scale. The comparison between ABC transporter pathways of P.abyssi and those of another archaeon, M.jannaschii, revealed that there was no peptide uptake ABC transporter system in M.jannaschii. This may result from their different metabolic types.

The Relationship between Yeast Coexpressed Gene Clusters and Their Upstream cis-acting Elements.

DNA microarrays bring biology a new approach to study functions of genes and genomes from their expression pattern on a genomic scale. With its fully sequenced genome and newly published expression patterns available, budding yeast (Saccharomyces cerevisiae) was chosen to carry out an investigation of the relationship between gene 5' upstream cis-acting elements and expression patterns using bioinformatic tools. Results show that genes in the same cluster share common cis-acting element candidates and can be regulated by same transfactors. In the sites found by this study, some sites are corresponding to known cis-acting elements, while others may indicate new ones that can be tested by experiments. The results are helpful to understand more about gene functions, metabolic pathways and genetic networks.