Cost-effectiveness and safety of the molecular targeted drugs afatinib, gefitinib and erlotinib as first-line treatments for patients with advanced EGFR mutation-positive non-small-cell lung cancer

Epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKIs), including gefitinib, erlotinib and afatinib are standard first-line treatments for EGFR gene mutation-positive non-small cell lung cancer. The present study aimed to compare the cost-effectiveness of using erlotinib, afatinib or gefitinib. The safety of EGFR-TKIs was also investigated. Expected costs were calculated based on data from patients with advanced EGFR mutation-positive non-small-cell lung cancer who were treated with gefitinib, erlotinib or afatinib. Literature was collected to obtain the necessary clinical information for calculating the probability and the validity of each chemotherapy. Median survival time (MST) was used to evaluate the therapeutic effect of the regimens. The cost-effectiveness ratio was calculated using expected costs and MSTs for the three regimens. The cost-effectiveness ratio per month was JPY 386,859.4/MST for afatinib, JPY 264,788.7/MST for gefitinib and JPY 397,039.9/MST for erlotinib. Significant differences were observed between the three groups (p<0.001). The incremental cost-effectiveness ratio (ICER) of gefitinib compared with afatinib per month was JPY 122,070.7/MST. The ICER of gefitinib compared with erlotinib was JPY -69,605.9/MST. Adverse effects of Grade 3 and higher, including diarrhoea (28.6%) and paronychia (14.3%) were observed in the afatinib treatment group. Paronychia (23.1%) was observed in the erlotinib treatment group, while none were observed in the gefitinib treatment group. These findings demonstrate that gefitinib is more cost effective in comparison with the afatinib and erlotinib regimens, although the afatinib and erlotinib regimens were well-tolerated and produce sufficient effects.

Infection-related complications during treatment for childhood acute lymphoblastic leukemia

Background

Comprehensive studies on neutropenia and infection-related complications in patients with acute lymphoblastic leukemia (ALL) are lacking.

Patients and methods

We evaluated infection-related complications that were grade ≥3 on National Cancer Institute's Common Terminology Criteria for Adverse Events (version 3.0) and their risk factors in 409 children with newly diagnosed ALL throughout the treatment period.

Results

Of the 2420 infection episodes, febrile neutropenia and clinically or microbiologically documented infection were seen in 1107 and 1313 episodes, respectively. Among documented infection episodes, upper respiratory tract was the most common site (n = 389), followed by ear (n = 151), bloodstream (n = 147), and gastrointestinal tract (n = 145) infections. These episodes were more common during intensified therapy phases such as remission induction and reinduction, but respiratory and ear infections, presumably viral in origin, also occurred during continuation phases. The 3-year cumulative incidence of infection-related death was low (1.0±0.9%, n = 4), including 2 from Bacillus cereus bacteremia. There was no fungal infection-related mortality. Age 1-9.9 years at diagnosis was associated with febrile neutropenia (P = 0.002) during induction and febrile neutropenia and documented infection (both P < 0.001) during later continuation. White race was associated with documented infection (P = 0.034) during induction. Compared with low-risk patients, standard- and high-risk patients received more intensive therapy during early continuation and had higher incidences of febrile neutropenia (P < 0.001) and documented infections (P = 0.043). Furthermore, poor neutrophil surge after dexamethasone pulses during continuation, which can reflect the poor bone marrow reserve, was associated with infections (P < 0.001).

Conclusions

The incidence of infection-related death was low. However, young age, white race, intensive chemotherapy, and lack of neutrophil surge after dexamethasone treatment were associated with infection-related complications. Close monitoring for prompt administration of antibiotics and modification of chemotherapy should be considered in these patients.

Safety of an oral anticancer agent (trifluridine/tipiracil combination tablet) in patients with advanced and recurrent colorectal cancer

We retrospectively studied the safety of trifluridine/tipiracil combination tablet (TAS-102) monotherapy in patients with advanced and recurrent colorectal cancer. Adverse events to TAS-102 monotherapy were observed in 22 out of 23 cases (95.7%). The most frequent adverse events were neutropenia (69.6%), nausea (53.2%), and malaise (30.4%). Treatment was postponed in 54 (59.3%) out of 91 courses, and in 34 (66.7%) of these 54 courses, the delay in treatment was due to bone marrow suppression. Seven patients with peritoneal metastases suffered from nausea, whilst none of the patients without peritoneal metastases had nausea (p = 0.0139). Nausea and vomiting during a previous chemotherapy cycle was significantly associated with nausea after TAS-102 treatment (p = 0.0007), and the treatment cycles were significantly longer in patients with grade 3 or 4 neutropenia (p = 0.0061). Our results suggest that the incidence of nausea was higher in patients treated with TAS-102. Therefore, it is important to inform patients of the risk of these toxicities and to provide enhanced supportive care. Moreover, we recommend that, for patients with repeated treatment postponement due to neutropenia, the dosage should be fixed based on therapeutic efficacy and prognosis.

Structural and functional analyses of Barth syndrome-causing mutations and alternative splicing in the tafazzin acyltransferase domain

Tafazzin is a mitochondrial phospholipid transacylase, and its mutations cause Barth syndrome (BTHS). Human tafazzin gene produces four distinct alternatively spliced transcripts. To understand the molecular mechanisms of tafazzin deficiency, we performed an atomic resolution analysis of the influence of the BTHS mutations and of alternative splicing on the structure and function of tafazzin. From the three-dimensional (3D) homology modeling of tafazzin, we identified candidate amino acid residues that contribute to cardiolipin binding and to mitochondrial membrane associations that facilitate acyl-transfer reactions. Primate specific exon 5, which is alternatively spliced, is predicted to correspond to an intrinsically unstructured region in the protein. We proposed that this region should change the substrate-binding affinity and/or contribute to primate-specific molecular interactions. Exon 7, another alternatively spliced exon, encodes a region forming a part of the putative substrate-binding cleft, suggesting that the gene products lacking exon 7 will lose their substrate-binding ability. We demonstrate a clear localization of the BTHS mutations at residues responsible for membrane association, substrate binding, and the conformational stability of tafazzin. These findings provide new insights into the function of defective tafazzin and the pathogenesis of BTHS at the level of protein 3D structure and the evolution of alternatively spliced exons in primates.

AS-EAST: a functional annotation tool for putative proteins encoded by alternatively spliced transcripts

Summary: Alternative Splicing Effects ASsessment Tools (AS-EAST) is an online tool for the functional annotation of putative proteins encoded by transcripts generated by alternative splicing (AS). When provided with a transcript sequence, AS-EAST identifies regions altered by AS events in the putative protein sequence encoded by the transcript. Users can evaluate the predicted function of the putative protein by inspecting whether functional domains are included in the altered regions. Moreover, users can infer the loss of inter-molecular interactions in the protein network according to whether the AS events affect interaction residues observed in the 3D structure of the reference isoform. The information obtained from AS-EAST will help to design experimental analyses for the functional significance of novel splice isoforms.

Availability: The online tool is freely available at http://as-alps.nagahama-i-bio.ac.jp/ASEAST/.

Contact:m_shionyu@nagahama-i-bio.ac.jp

The role of IL-33 and its receptor ST2 in human nasal epithelium with allergic rhinitis

BACKGROUND

Interleukin (IL)-33 is a novel member of the IL-1 cytokine family and a ligand for the orphan IL-1 family receptor ST2. The IL-33 induces T helper 2-type inflammatory responses and is considered to play a crucial rule in allergic inflammations, such as asthma and atopic dermatitis. However, the role of IL-33 and its receptor ST2 in allergic rhinitis remains unknown.

OBJECTIVE

We investigated expression of IL-33 and ST2 in the nasal epithelium of patients with allergic rhinitis and the mechanisms of the production of cytokines/chemokines induced by treatment with IL-33 using normal human nasal epithelial cells (HNECs) in vitro.

METHODS

Expression of IL-33 and ST2 in normal and allergic rhinitis nasal mucosa was evaluated by reverse transcription- and real-time polymerase chain reactions and immunohistochemical methods. The IL-33 in serum, and IL-8 and GM-CSF were measured by ELISA. For in vitro experiments, HNECs in primary culture were used.

RESULTS

The IL-33 levels in the sera of patients with allergic rhinitis were significantly higher than that in normal controls. Expression of IL-33 and ST2 was significantly elevated in the epithelium from patients with allergic rhinitis. The IL-33 mRNA in HNECs in vitro was significantly induced by treatment with IFN-γ and the toll-like receptor 9 ligand ODN2006. The IL-33-induced production of IL-8 and GM-CSF from HNECs in vitro was significantly suppressed by corticosteroid treatment and distinct signal transduction inhibitors of ERK, p38 MAPK, JNK, NF-κB and epidermal growth factor receptor.

CONCLUSIONS AND CLINICAL RELEVANCE

The IL-33 and its receptor ST2 play important roles in allergic rhinitis. The IL-33-mediated inflammatory responses via ST2 are regulated by distinct signalling pathways in HNECs and the IL-33/ST2 pathway may provide new therapeutic targets for allergic rhinitis.

Revisiting gap locations in amino acid sequence alignments and a proposal for a method to improve them by introducing solvent accessibility

In comparative modeling, the quality of amino acid sequence alignment still constitutes a major bottleneck in the generation of high quality models of protein three-dimensional (3D) structures. Substantial efforts have been made to improve alignment quality by revising the substitution matrix, introducing multiple sequences, replacing dynamic programming with hidden Markov models, and incorporating 3D structure information. Improvements in the gap penalty have not been a major focus, however, following the development of the affine gap penalty and of the secondary structure dependent gap penalty. We revisited the correlation between protein 3D structure and gap location in a large protein 3D structure data set, and found that the frequency of gap locations approximated to an exponential function of the solvent accessibility of the inserted residues. The nonlinearity of the gap frequency as a function of accessibility corresponded well to the relationship between residue mutation pattern and residue accessibility. By introducing this relationship into the gap penalty calculation for pairwise alignment between template and target amino acid sequences, we were able to obtain a sequence alignment much closer to the structural alignment. The quality of the alignments was substantially improved on a pair of sequences with identity in the “twilight zone” between 20 and 40%. The relocation of gaps by our new method made a significant improvement in comparative modeling, exemplified here by the Bacillus subtilis yitF protein. The method was implemented in a computer program, ALAdeGAP (ALignment with Accessibility dependent GAp Penalty), which is available at http://cib.cf.ocha.ac.jp/target_protein/. Proteins 2011; © 2011 Wiley-Liss, Inc.

RESOPS: a database for analyzing the correspondence of RNA editing sites to protein three-dimensional structures

Transcripts from mitochondrial and chloroplast DNA of land plants often undergo cytidine to uridine conversion-type RNA editing events. RESOPS is a newly built database that specializes in displaying RNA editing sites of land plant organelles on protein three-dimensional (3D) structures to help elucidate the mechanisms of RNA editing for gene expression regulation. RESOPS contains the following information: unedited and edited cDNA sequences with notes for the target nucleotides of RNA editing, conceptual translation from the edited cDNA sequence in pseudo-UniProt format, a list of proteins under the influence of RNA editing, multiple amino acid sequence alignments of edited proteins, the location of amino acid residues coded by codons under the influence of RNA editing in protein 3D structures and the statistics of biased distributions of the edited residues with respect to protein structures. Most of the data processing procedures are automated; hence, it is easy to keep abreast of updated genome and protein 3D structural data. In the RESOPS database, we clarified that the locations of residues switched by RNA editing are significantly biased to protein structural cores. The integration of different types of data in the database also help advance the understanding of RNA editing mechanisms. RESOPS is accessible at http://cib.cf.ocha.ac.jp/RNAEDITING/.

AS-ALPS: a database for analyzing the effects of alternative splicing on protein structure, interaction and network in human and mouse

We have constructed a database, AS-ALPS (alternative splicing-induced alteration of protein structure), which provides information that would be useful for analyzing the effects of alternative splicing (AS) on protein structure, interactions with other bio-molecules and protein interaction networks in human and mouse. Several AS events have been revealed to contribute to the diversification of protein structure, which results in diversification of interaction partners or affinities, which in turn contributes to regulation of bio-molecular networks. Most AS variants, however, are only known at the sequence level. It is important to determine the effects of AS on protein structure and interaction, and to provide candidates for experimental targets that are relevant to network regulation by AS. For this purpose, the three-dimensional (3D) structures of proteins are valuable sources of information; however, these have not been fully exploited in any other AS-related databases. AS-ALPS is the only AS-related database that describes the spatial relationships between protein regions altered by AS ('AS regions') and both the proteins' hydrophobic cores and sites of inter-molecular interactions. This information makes it possible to infer whether protein structural stability and/or protein interaction are affected by each AS event. AS-ALPS can be freely accessed at http://as-alps.nagahama-i-bio.ac.jp and http://genomenetwork.nig.ac.jp/as-alps/.

Eosinophilic oesophagitis in patients presenting with dysphagia--a prospective analysis

BACKGROUND

Eosinophilic oesophagitis (EoO) may be a common finding in adults presenting with dysphagia.

AIM

To identify the risk factors and prevalence of EoO in an adult population with dysphagia.

METHODS

All patients with dysphagia referred for an upper endoscopy (EGD) were asked to participate. Patients completed a detailed questionnaire followed by EGD with four quadrant biopsies in the distal and mid-oesophagus. Primary endpoint was the prevalence of EoO; secondary endpoints included age, gender, asthma, food allergies, gastro-oesophageal reflux disease/dysphagia score and endoscopic findings.

RESULTS

Two hundred and sixty-one patients enrolled between December 2005 and January 2007. Thirty-one patients (12%) met pathological criteria for EoO. There was no difference in EoO prevalence within each gender. Mean age of EoO patients was 42 +/- 15 vs. 61 +/- 15 for non-EoO patients (P < 0.001). EoO was diagnosed in 35% of patients <50 years of age. EoO was present in 22% of asthmatics vs. 9% non-asthmatics (P < 0.01). EoO was present in 36.8% of patients with self-reported food allergies vs. 9.3% those without allergy (P < 0.001). A 13/31(42%) of EoO patients did not have the classic EGD findings (rings +/- furrows) and would have been missed without oesophageal biopsies.

CONCLUSIONS

Eosinophilic oesophagitis was diagnosed in 12% of the patients presenting with dysphagia with relative risk of 9.5 if age <50 years. Oesophageal biopsies are warranted in patients presenting with dysphagia especially in the younger population. Patients may not present with classic endoscopic findings and EoO can be missed without biopsies.

Correlation between amino acid residues converted by RNA editing and functional residues in protein three-dimensional structures in plant organelles

BACKGROUND

In plant organelles, specific messenger RNAs (mRNAs) are subjected to conversion editing, a process that often converts the first or second nucleotide of a codon and hence the encoded amino acid. No systematic patterns in converted sites were found on mRNAs, and the converted sites rarely encoded residues located at the active sites of proteins. The role and origin of RNA editing in plant organelles remain to be elucidated.

RESULTS

Here we study the relationship between amino acid residues encoded by edited codons and the structural characteristics of these residues within proteins, e.g., in protein-protein interfaces, elements of secondary structure, or protein structural cores. We find that the residues encoded by edited codons are significantly biased toward involvement in helices and protein structural cores. RNA editing can convert codons for hydrophilic to hydrophobic amino acids. Hence, only the edited form of an mRNA can be translated into a polypeptide with helix-preferring and core-forming residues at the appropriate positions, which is often required for a protein to form a functional three-dimensional (3D) structure.

CONCLUSION

We have performed a novel analysis of the location of residues affected by RNA editing in proteins in plant organelles. This study documents that RNA editing sites are often found in positions important for 3D structure formation. Without RNA editing, protein folding will not occur properly, thus affecting gene expression. We suggest that RNA editing may have conferring evolutionary advantage by acting as a mechanism to reduce susceptibility to DNA damage by allowing the increase in GC content in DNA while maintaining RNA codons essential to encode residues required for protein folding and activity.

The H-Invitational Database (H-InvDB), a comprehensive annotation resource for human genes and transcripts

Here we report the new features and improvements in our latest release of the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/), a comprehensive annotation resource for human genes and transcripts. H-InvDB, originally developed as an integrated database of the human transcriptome based on extensive annotation of large sets of full-length cDNA (FLcDNA) clones, now provides annotation for 120 558 human mRNAs extracted from the International Nucleotide Sequence Databases (INSD), in addition to 54 978 human FLcDNAs, in the latest release H-InvDB_4.6. We mapped those human transcripts onto the human genome sequences (NCBI build 36.1) and determined 34 699 human gene clusters, which could define 34 057 (98.1%) protein-coding and 642 (1.9%) non-protein-coding loci; 858 (2.5%) transcribed loci overlapped with predicted pseudogenes. For all these transcripts and genes, we provide comprehensive annotation including gene structures, gene functions, alternative splicing variants, functional non-protein-coding RNAs, functional domains, predicted sub cellular localizations, metabolic pathways, predictions of protein 3D structure, mapping of SNPs and microsatellite repeat motifs, co-localization with orphan diseases, gene expression profiles, orthologous genes, protein-protein interactions (PPI) and annotation for gene families. The current H-InvDB annotation resources consist of two main views: Transcript view and Locus view and eight sub-databases: the DiseaseInfo Viewer, H-ANGEL, the Clustering Viewer, G-integra, the TOPO Viewer, Evola, the PPI view and the Gene family/group.

Retention of local conformational compactness in unfolding of barnase; Contribution of end-to-end interactions within quasi-modules

To understand how protein reduces the conformational space to be searched for the native structure, it is crucial to characterize ensembles of conformations on the way of folding processes, in particular ensembles of relatively long-range structures connecting between an extensively unfolded state and a state with a native-like overall chain topology. To analyze such intermediate conformations, we performed multiple unfolding molecular dynamics simulations of barnase at 498K. Some short-range structures such as part of helix and turn were well sustained while most of the secondary structures and the hydrophobic cores were eventually lost, which is consistent with the results by other experimental and computational studies. The most important novel findings were persistence of long-range relatively compact substructures, which was captured by exploiting the concept of module. Module is originally introduced to describe the hierarchical structure of a globular protein in the native state. Modules are conceptually such relatively compact substructures that are resulted from partitioning the native structure of a globular protein completely into several contiguous segments with the least extended conformations. We applied this concept of module to detect a possible hierarchical structure of each snapshot structure in unfolding processes as well. Along with this conceptual extension, such detected relatively compact substructures are named quasi-modules. We found almost perfect persistence of quasi-module boundaries that are positioned close to the native module boundaries throughout the unfolding trajectories. Relatively compact conformations of the quasi-modules seemed to be retained mainly by hydrophobic interactions formed between residues located at both terminal regions within each module. From these results, we propose a hypothesis that hierarchical folding with the early formation of quasi-modules effectively reduces search space for the native structure.

Alternative splicing in human transcriptome: Functional and structural influence on proteins

Alternative splicing is a molecular mechanism that produces multiple proteins from a single gene, and is thought to produce variety in proteins translated from a limited number of genes. Here we analyzed how alternative splicing produced variety in protein structure and function, by using human full-length cDNAs on the assumption that all of the alternatively spliced mRNAs were translated to proteins. We found that the length of alternatively spliced amino acid sequences, in most cases, fell into a size shorter than that of average protein domain. We evaluated comprehensively the presumptive three-dimensional structures of the alternatively spliced products to assess the impact of alternative splicing on gene function. We found that more than half of the products encoded proteins which were involved in signal transduction, transcription and translation, and more than half of alternatively spliced regions comprised interaction sites between proteins and their binding partners, including substrates, DNA/RNA, and other proteins. Intriguingly, 67% of the alternatively spliced isoforms showed significant alterations to regions of the protein structural core, which likely resulted in large conformational change. Based on those findings, we speculate that there are a large number of cases that alternative splicing modulates protein networks through significant alteration in protein conformation.

Large-scale identification and characterization of alternative splicing variants of human gene transcripts using 56,419 completely sequenced and manually annotated full-length cDNAs

We report the first genome-wide identification and characterization of alternative splicing in human gene transcripts based on analysis of the full-length cDNAs. Applying both manual and computational analyses for 56,419 completely sequenced and precisely annotated full-length cDNAs selected for the H-Invitational human transcriptome annotation meetings, we identified 6877 alternative splicing genes with 18 297 different alternative splicing variants. A total of 37,670 exons were involved in these alternative splicing events. The encoded protein sequences were affected in 6005 of the 6877 genes. Notably, alternative splicing affected protein motifs in 3015 genes, subcellular localizations in 2982 genes and transmembrane domains in 1348 genes. We also identified interesting patterns of alternative splicing, in which two distinct genes seemed to be bridged, nested or having overlapping protein coding sequences (CDSs) of different reading frames (multiple CDS). In these cases, completely unrelated proteins are encoded by a single locus. Genome-wide annotations of alternative splicing, relying on full-length cDNAs, should lay firm groundwork for exploring in detail the diversification of protein function, which is mediated by the fast expanding universe of alternative splicing variants.

Coverage of whole proteome by structural genomics observed through protein homology modeling database

We have been developing FAMSBASE, a protein homology-modeling database of whole ORFs predicted from genome sequences. The latest update of FAMSBASE ( http://daisy.nagahama-i-bio.ac.jp/Famsbase/ ), which is based on the protein three-dimensional (3D) structures released by November 2003, contains modeled 3D structures for 368,724 open reading frames (ORFs) derived from genomes of 276 species, namely 17 archaebacterial, 130 eubacterial, 18 eukaryotic and 111 phage genomes. Those 276 genomes are predicted to have 734,193 ORFs in total and the current FAMSBASE contains protein 3D structure of approximately 50% of the ORF products. However, cases that a modeled 3D structure covers the whole part of an ORF product are rare. When portion of an ORF with 3D structure is compared in three kingdoms of life, in archaebacteria and eubacteria, approximately 60% of the ORFs have modeled 3D structures covering almost the entire amino acid sequences, however, the percentage falls to about 30% in eukaryotes. When annual differences in the number of ORFs with modeled 3D structure are calculated, the fraction of modeled 3D structures of soluble protein for archaebacteria is increased by 5%, and that for eubacteria by 7% in the last 3 years. Assuming that this rate would be maintained and that determination of 3D structures for predicted disordered regions is unattainable, whole soluble protein model structures of prokaryotes without the putative disordered regions will be in hand within 15 years. For eukaryotic proteins, they will be in hand within 25 years. The 3D structures we will have at those times are not the 3D structure of the entire proteins encoded in single ORFs, but the 3D structures of separate structural domains. Measuring or predicting spatial arrangements of structural domains in an ORF will then be a coming issue of structural genomics.

Survey of Conserved Alternative Splicing Events of mRNAs Encoding SR Proteins in Land Plants

The serine/arginine-rich (SR) protein family plays an important role in constitutive and alternative splicing (AS). These proteins regulate AS in a tissue-specific and stress-responsive manner. Pre-mRNAs encoding SR proteins are often alternatively spliced, and these AS events may be important for the regulation of AS events of other pre-mRNAs. In this study, we analyzed AS events of SR proteins in Arabidopsis thaliana and Oryza sativa (rice). We found three sets of AS events conserved between Arabidopsis and rice. These conserved AS events were found in the plant-novel-SR protein, SC35-like (SCL), and two-Zn-knuckles-type 9G8 subfamilies. Each member of these subfamilies has at least one RNA recognition motif (RRM) and at least one intron in the RRM-encoded region. We found that the conserved AS events occurred in these introns and, in each case, the conserved AS events resulted in mature mRNAs encoding proteins with incomplete RRMs. To search for the evolutionary origin of these AS events, we analyzed SR proteins in Physcomitrella patens (moss) in addition to those in Arabidopsis and rice. We found moss homologues of the plant-novel-SR protein, SCL, and the two-Zn-knuckles-type 9G8 subfamilies in silico, and these homologues have long introns at the same location of the conserved AS sites in Arabidopsis and rice. Such long introns are quite specific for alternatively spliced introns concerning the Arabidopsis SR protein genes. The long introns found in the moss SR protein genes strongly suggested that conserved AS events in moss SR protein genes might be similar to those in Arabidopsis and rice. We traced the evolutionary origin of the conserved AS events to 400 MYA, when plants first invaded land. These events are likely important in the regulation of whole AS events and likely contribute to the complicated transcriptome described by AS. The complicated transcriptome created by regulated AS events might have provided plants tolerance against droughts or temperature shifts and given them the ability to live on land.

Disulfide linkages and a three-dimensional structure model of the extracellular ligand-binding domain of guanylyl cyclase C

Guanylyl cyclase C (GC-C) is a single-transmembrane receptor that is specifically activated by endogenous ligands, including guanylin, and the exogenous ligand, heat-stable enterotoxin. Using combined HPLC separation and MS analysis techniques the positions of the disulfide linkages in the extracellular ligand-binding domain (ECD) of GC-C were determined to be between Cys7-Cys94, Cys72-Cys77, Cys101-Cys128 and Cys179-Cys226. Furthermore, a three-dimensional structural model of the ECD was constructed by homology modeling, using the structure of the ECD of GC-A as a template (van den Akker et al., 2000, Nature, 406: 101-104) and the information of the disulfide linkages. Although the GC-C model was similar to the known structure of GC-A, importantly its ligand-binding site appears to be located on the quite different region from that in GC-A.

An empirical approach for detecting nucleotide-binding sites on proteins

Protein structure data in the PDB (Protein Data Bank) were used to construct empirical scores of nucleotide-protein interactions. A simple strategy to evaluate the spatial distribution of protein atoms around the base moieties of nucleotides was applied to categorize adenine, guanine, nicotinamide and flavin nucleotide-binding sites. In addition to the known nucleotide-binding motifs, the empirical scores detected several other features that were shared among proteins with different folds. The empirical scores were also used to predict the binding sites on protein molecules and a comprehensive test of the prediction system was performed. As a result, adenine, guanine, nicotinamide and flavin sites were detected with efficiencies of 31, 29, 32 and 40%, respectively. The predictions were judged to be successful if the predicted base with the best score was located within a 3.0 A r.m.s.d. from the known ligand positions.

Het-PDB Navi.: a database for protein-small molecule interactions

The genomes of more than 100 species have been sequenced, and the biological functions of encoded proteins are now actively being researched. Protein function is based on interactions between proteins and other molecules. One approach to assuming protein function based on genomic sequence is to predict interactions between an encoded protein and other molecules. As a data source for such predictions, knowledge regarding known protein-small molecule interactions needs to be compiled. We have, therefore, surveyed interactions between proteins and other molecules in Protein Data Bank (PDB), the protein three-dimensional (3D) structure database. Among 20,685 entries in PDB (April, 2003), 4,189 types of small molecules were found to interact with proteins. Biologically relevant small molecules most often found in PDB were metal ions, such as calcium, zinc, and magnesium. Sugars and nucleotides were the next most common. These molecules are known to act as cofactors for enzymes and/or stabilizers of proteins. In each case of interactions between a protein and small molecule, we found preferred amino acid residues at the interaction sites. These preferences can be the basis for predicting protein function from genomic sequence and protein 3D structures. The data pertaining to these small molecules were collected in a database named Het-PDB Navi., which is freely available at http://daisy.nagahama-i-bio.ac.jp/golab/hetpdbnavi.html and linked to the official PDB home page.

Role of KaiC phosphorylation in the circadian clock system of Synechococcus elongatus PCC 7942

In the cyanobacterium Synechococcus elongatus PCC 7942, KaiA, KaiB, and KaiC are essential proteins for the generation of a circadian rhythm. KaiC is proposed as a negative regulator of the circadian expression of all genes in the genome, and its phosphorylation is regulated positively by KaiA and negatively by KaiB and shows a circadian rhythm in vivo. To study the functions of KaiC phosphorylation in the circadian clock system, we identified two autophosphorylation sites, Ser-431 and Thr-432, by using mass spectrometry (MS). We generated Synechococcus mutants in which these residues were substituted for alanine by using site-directed mutagenesis. Phosphorylation of KaiC was reduced in the single mutants and was completely abolished in the double mutant, indicating that KaiC is also phosphorylated at these sites in vivo. These mutants lost circadian rhythm, indicating that phosphorylation at each of the two sites is essential for the control of the circadian oscillation. Although the nonphosphorylatable mutant KaiC was able to form a hexamer in vitro, it failed to form a clock protein complex with KaiA, KaiB, and SasA in the Synechococcus cells. When nonphosphorylatable KaiC was overexpressed, the kaiBC promoter activity was only transiently repressed. These results suggest that KaiC phosphorylation regulates its transcriptional repression activity by controlling its binding affinity for other clock proteins.

Integrative annotation of 21,037 human genes validated by full-length cDNA clones

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for non-protein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology.

Structure and function of a family 10 beta-xylanase chimera of Streptomyces olivaceoviridis E-86 FXYN and Cellulomonas fimi Cex

The catalytic domain of xylanases belonging to glycoside hydrolase family 10 (GH10) can be divided into 22 modules (M1 to M22; Sato, Y., Niimura, Y., Yura, K., and Go, M. (1999) Gene (Amst.) 238, 93-101). Inspection of the crystal structure of a GH10 xylanase from Streptomyces olivaceoviridis E-86 (SoXyn10A) revealed that the catalytic domain of GH10 xylanases can be dissected into two parts, an N-terminal larger region and C-terminal smaller region, by the substrate binding cleft, corresponding to the module border between M14 and M15. It has been suggested that the topology of the substrate binding clefts of GH10 xylanases are not conserved (Charnock, S. J., Spurway, T. D., Xie, H., Beylot, M. H., Virden, R., Warren, R. A. J., Hazlewood, G. P., and Gilbert, H. J. (1998) J. Biol. Chem. 273, 32187-32199). To facilitate a greater understanding of the structure-function relationship of the substrate binding cleft of GH10 xylanases, a chimeric xylanase between SoXyn10A and Xyn10A from Cellulomonas fimi (CfXyn10A) was constructed, and the topology of the hybrid substrate binding cleft established. At the three-dimensional level, SoXyn10A and CfXyn10A appear to possess 5 subsites, with the amino acid residues comprising subsites -3 to +1 being well conserved, although the +2 subsites are quite different. Biochemical analyses of the chimeric enzyme along with SoXyn10A and CfXyn10A indicated that differences in the structure of subsite +2 influence bond cleavage frequencies and the catalytic efficiency of xylooligosaccharide hydrolysis. The hybrid enzyme constructed in this study displays fascinating biochemistry, with an interesting combination of properties from the parent enzymes, resulting in a low production of xylose.

Relationship between gastroesophageal reflux disease and myocardial ischemia. Effect of reflux on temporary activity of autonomic nervous system

PURPOSE

Assessment of the gastroesophageal reflux disease (GERD) influence on myocardial ischemia and autonomic nervous system (ANS) activity.

MATERIAL AND METHODS

In 50 patients with angiographically confirmed ischemic heart disease (IHD) in I-III CCS class, simultaneous 24-hour ECG and esophageal pH-metry monitoring was performed. We assessed: (1) GERD occurrence in patients with IHD, (2) influence of pathological reflux (PR) on myocardial ischemia--number and total duration of ST depression episodes in GERD and non-GERD patients, (3) temporary activity of ANS was determined according to the dynamics of spectral HRV (Heart Rate Variability) analysis components (LF, HF, VLF, LF/HF).

RESULTS

23 patients (46%) fullfilled the GERD criteria. Patients with GERD had significantly higher number of ST depression episodes (4.13 vs 2.85, p = 0.013) as well as longer total duration of ischemia (64.73 vs 35.2 min, p = 0.034). Spectral HRV analysis showed the significant decrease of LF/HF ratio (p < 0.035), which indicates the sympathovagal balance shift towards the parasympathetic system caused by PR.

CONCLUSIONS

1. GERD is frequent condition in patients with angiographically confirmed IHD. Coexistence of GERD may predispose to the myocardial ischemia. 2. Gastroesophageal reflux may cause the shift of sympathovagal balance towards its parasympathetic component. This mechanism may induce esophago-cardiac reflex, leading to diminished myocardial perfusion.

Distinct interaction of versican/PG-M with hyaluronan and link protein

The proteoglycan aggregate is the major structural component of the cartilage matrix, comprising hyaluronan (HA), link protein (LP), and a large chondroitin sulfate (CS) proteoglycan, aggrecan. Here, we found that another member of aggrecan family, versican, biochemically binds to both HA and LP. Functional analyses of recombinant looped domains (subdomains) A, B, and B' of the N-terminal G1 domain revealed that the B-B' segment of versican is adequate for binding to HA and LP, whereas A and B-B' of aggrecan bound to LP and HA, respectively. BIAcore trade mark analyses showed that the A subdomain of versican G1 enhances HA binding but has a negligible effect on LP binding. Overlay sensorgrams demonstrated that versican G1 or its B-B' segment forms a complex with both HA and LP. We generated a molecular model of the B-B' segment, in which a deletion and an insertion of B' and B are critical for stable structure and HA binding. These results provide important insights into the mechanisms of formation of the proteoglycan aggregate and HA binding of molecules containing the link module.

Chondroitin sulfate synthase-2. Molecular cloning and characterization of a novel human glycosyltransferase homologous to chondroitin sulfate glucuronyltransferase, which has dual enzymatic activities

Chondroitin sulfate is found in a variety of tissues as proteoglycans and consists of repeating disaccharide units of N-acetylgalactosamine and glucuronic acid residues with sulfate residues at various places. We found a novel human gene (GenBank accession number AB086063) that possesses a sequence homologous with the human chondroitin sulfate glucuronyltransferase gene which we recently cloned and characterized. The full-length open reading frame encodes a typical type II membrane protein comprising 775 amino acids. The protein had a domain containing beta 3-glycosyltransferase motif but lacked a typical beta 4-glycosyltransferase motif, which is the same as chondroitin sulfate glucuronyltransferase, whereas chondroitin synthase had both domains. The putative catalytic domain was expressed in COS-7 cells as a soluble enzyme. Surprisingly, both glucuronyltransferase and N-acetylgalactosaminyltransferase activities were observed when chondroitin, chondroitin sulfate, and their oligosaccharides were used as the acceptor substrates. The reaction products were identified to have the linkage of GlcUA beta 1-3GalNAc and GalNAc beta 1-4GlcUA at the non-reducing terminus of chondroitin for glucuronyltransferase activity and N-acetylgalactosaminyltransferase activity, respectively. Quantitative real time PCR analysis revealed that the transcripts were ubiquitously expressed in various human tissues but highly expressed in the pancreas, ovary, placenta, small intestine, and stomach. These results indicate that this enzyme could synthesize chondroitin sulfate chains as a chondroitin sulfate synthase that has both glucuronyltransferase and N-acetylgalactosaminyltransferase activities. Sequence analysis based on three-dimensional structure revealed the presence of not typical but significant beta 4-glycosyltransferase architecture.

Enlarged FAMSBASE: protein 3D structure models of genome sequences for 41 species

Enlarged FAMSBASE is a relational database of comparative protein structure models for the whole genome of 41 species, presented in the GTOP database. The models are calculated by Full Automatic Modeling System (FAMS). Enlarged FAMSBASE provides a wide range of query keys, such as name of ORF (open reading frame), ORF keywords, Protein Data Bank (PDB) ID, PDB heterogen atoms and sequence similarity. Heterogen atoms in PDB include cofactors, ligands and other factors that interact with proteins, and are a good starting point for analyzing interactions between proteins and other molecules. The data may also work as a template for drug design. The present number of ORFs with protein 3D models in FAMSBASE is 183 805, and the database includes an average of three models for each ORF. FAMSBASE is available at http://famsbase.bio.nagoya-u.ac.jp/famsbase/.

Decreased expression of DNA-dependent protein kinase, a DNA repair protein, during human colon carcinogenesis

DNA-dependent protein kinase (DNA-PK), consisting of a catalytic subunit (DNA-PKcs) and the Ku70 and Ku86 proteins, participates in the repair of DNA double-strand breaks (DSBs). We assessed its expression immunohistochemically in normal human colon tissue, colon adenomas, colon carcinomas, and normal tissue distant from carcinomas. Normal colonocytes expressed all DNA-PK proteins. Compared with the expression in normal tissue [176.62 +/- 18.56 (the intensity of expression x the percentage of cells expressing this protein), mean + SE], the expression of Ku70 was significantly reduced in adenomas (36.62 +/- 11.09; P < 0.001) and carcinomas (85.68 +/- 15.76; P < 0.01), as was the expression of Ku86 [(113.10 +/- 10.22 versus 41.66 +/- 14.71 in adenomas (P < 0.01) or versus 85.68 +/- 15.76 in carcinomas (P < 0.05)]. The expression of DNA-PKcs was not significantly changed. The marked underexpression of Ku70 and Ku86 starting at the adenoma stage may be crucial to the development of colon cancer.

Two KaiA-binding domains of cyanobacterial circadian clock protein KaiC

kaiABC, a gene cluster, encodes KaiA, KaiB and KaiC proteins that are essential to circadian rhythms in the unicellular cyanobacterium Synechococcus sp. strain PCC 7942. Kai proteins can interact with each other in all possible combinations. This study identified two KaiA-binding domains (C(KABD1) and C(KABD2)) in KaiC at corresponding regions of its duplicated structure. Clock mutations on the two domains and kaiA altered the strength of C(KABD)-KaiA interactions assayed by the yeast two-hybrid system. Thus, interaction between KaiA and KaiC through C(KABD1) and C(KABD2) is likely important for circadian timing in the cyanobacterium.

Ecosystem impacts of three sequential hurricanes (Dennis, Floyd, and Irene) on the United States' largest lagoonal estuary, Pamlico Sound, NC

Three sequential hurricanes, Dennis, Floyd, and Irene, affected coastal North Carolina in September and October 1999. These hurricanes inundated the region with up to 1 m of rainfall, causing 50- to 500-year flooding in the watershed of the Pamlico Sound, the largest lagoonal estuary in the United States and a key West Atlantic fisheries nursery. We investigated the ecosystem-level impacts on and responses of the Sound to the floodwater discharge. Floodwaters displaced three-fourths of the volume of the Sound, depressed salinity by a similar amount, and delivered at least half of the typical annual nitrogen load to this nitrogen-sensitive ecosystem. Organic carbon concentrations in floodwaters entering Pamlico Sound via a major tributary (the Neuse River Estuary) were at least 2-fold higher than concentrations under prefloodwater conditions. A cascading set of physical, chemical, and ecological impacts followed, including strong vertical stratification, bottom water hypoxia, a sustained increase in algal biomass, displacement of many marine organisms, and a rise in fish disease. Because of the Sound's long residence time ( approximately 1 year), we hypothesize that the effects of the short-term nutrient enrichment could prove to be multiannual. A predicted increase in the frequency of hurricane activity over the next few decades may cause longer-term biogeochemical and trophic changes in this and other estuarine and coastal habitats.

Capacity of thermomonospora alba XylA to impart thermostability in family F/10 chimeric xylanases

To reveal structure-function relationships of family F/10 glycanases, an in vitro molecular level shuffling experiment was conducted to accumulate useful amino acid residues from two homologous F/10 xylanases, FXYN of Streptomyces olivaceoviridis E-86 and XylA of Thermomonospora alba ULJB1, into a single chimeric xylanase. The parent genes were shuffled by crossovers at selected module borders using self-priming Polymerase Chain Reaction (PCR)s. The shuffled constructs, designated as FXYN-M3/4-XylA, FXYN-M9/10-XylA, and FXYN-M14/15-XylA were cloned and their nucleotide sequences were confirmed. Two chimera, FXYN-M3/4-XylA and FXYN-M14/15-XylA, demonstrated activity against RBB-xylan and were over-expressed as His-tag fusion proteins under control of T5 promoter of pQE60. The homogeneously pure chimeric proteins, FXYN-M3/4-XylA and FXYN-M14/15-XylA showed improved thermal and pH profiles compared to those of one of the parents, FXYN. This was apparently due to the influence of amino acids inherited from thermophilic XylA. Measured K(m) and kcat values were closer to those of the other parent, XylA. Interestingly, a significant level of heat tolerance up to 60 degrees C, was recorded for FXYN-M3/4-XylA in comparison to only 40 degrees C for FXYN-M14/15-XylA though their temperature optima did not correlates with their thermal stability. These results indicated that the amino acid residues of the larger T. alba XylA DNA fragment present in FXYN-M3/4-XylA were responsible for inducing its thermal stability.

Mutational analysis of genes involved in pilus structure, motility and transformation competency in the unicellular motile cyanobacterium Synechocystis sp. PCC 6803

The relevance of pilus-related genes to motility, pilus structure on the cell surface and competency of natural transformation was studied by gene disruption analysis in the unicellular motile cyanobacterium Synechocystis: sp. PCC 6803. The genes disrupted in this study were chosen as related to the pil genes for biogenesis of the type IV pili in a Gram-negative bacterium Pseudomonas aeruginosa. It was found that motility of Synechocystis cells was lost in the mutants of slr0063, slr1274, slr1275, slr1276, slr1277 and sll1694 together with a simultaneous loss of the thick pili on the cell surface. Competency of the natural transformation was lost in the mutants listed above and slr0197-disruptant. The gene slr0197 was previously predicted as a competence gene by a search with sequence-independent DNA-binding structure [Yura et al. (1999) DNA Res. 6: 75]. It was suggested that both DNA uptake for natural transformation and motility are mediated by a specific type IV-like pilus structure, while a putative DNA-binding protein encoded by slr0197 is additionally required for the DNA uptake. Based on the homology with the pil genes in P: aeruginosa, slr0063, slr1274, slr1275, slr1276, slr1277 and sll1694 were designated pilB1, pilM, pilN, pilO, pilQ and pilA1, respectively. The gene slr0197 was designated comA.

Module shuffling of a family F/10 xylanase: replacement of modules M4 and M5 of the FXYN of Streptomyces olivaceoviridis E-86 with those of the Cex of Cellulomonas fimi

To facilitate an understanding of structure-function relationships, chimeric xylanases were constructed by module shuffling between the catalytic domains of the FXYN from Streptomyces olivaceoviridis E-86 and the Cex from Cellulomonas fimi. In the family F/10 xylanases, the modules M4 and M5 relate to substrate binding so that modules M4 and M5 of the FXYN were replaced with those of the Cex and the chimeric enzymes denoted FCF-C4, FCF-C5 and FCF-C4,5 were constructed. The k(cat) value of FCF-C5 for p-nitrophenyl-beta-D-cellobioside was similar to that of the FXYN (2.2 s(-1)); however, the k(cat) value of FCF-C4 for p-nitrophenyl-beta-D-cellobioside was significantly higher (7.0 s(-1)). The loss of the hydrogen bond between E46 and S22 or the presence of the I49W mutation would be expected to change the position of Q88, which plays a pivotal role in discriminating between glucose and xylose, resulting in the increased k(cat) value observed for FCF-C4 acting on p-nitrophenyl-beta-D-cellobioside since module M4 directly interacts with Q88. To investigate the synergistic effects of the different modules, module M10 of the FCF-C4 chimera was replaced with that of the Cex. The effects of replacement of module M4 and M10 were almost additive with regard to the K:(m) and k(cat) values.

Significance of a two-domain structure in subunits of phycobiliproteins revealed by the normal mode analysis

Phycobiliproteins are basic building blocks of phycobilisomes, a supra-molecular assembly for the light-capturing function of photosynthesis in cyanobacteria and red algae. One functional form of phycobiliproteins is a trimeric form consisting of three identical units having C(3) symmetry, with each unit composed of two kinds of subunits, the alpha-subunit and beta-subunit. These subunits have similar chain folds and can be divided into either globin-like or X-Y helices domains. We studied the significance of this two-domain structure for their assembled structures and biological function (light-absorption) using a normal mode analysis to investigate dynamic aspects of their three-dimensional structures. We used C-phycocyanin (C-PC) as an example, and focused on the interactions between the two domains. The normal mode analysis was carried out for the following two cases: 1) the whole subunit, including the two domains; and 2) the globin-like domain alone. By comparing the dynamic properties, such as correlative movements between residues and the fluctuations of individual residues, we found that the X-Y helices domain plays an important role not only in the C(3) symmetry assemblies of the subunits in phycobiliproteins, but also in stabilizing the light absorption property by suppressing the fluctuation of the specific Asp residues near the chromophore. Interestingly, the conformation of the X-Y helices domain corresponds to that of a module in pyruvate phosphate dikinase (PPDK). The module in PPDK is involved in the interactions of two domains, just as the X-Y helices domain is involved in the interactions of two subunits. Finally, we discuss the mechanical construction of the C-PC subunits based on the normal mode analysis.

An investigation of the nature and function of module 10 in a family F/10 xylanase FXYN of Streptomyces olivaceoviridis E-86 by module shuffling with the Cex of Cellulomonas fimi and by site-directed mutagenesis

Although the amino acid homology in the catalytic domain of FXYN xylanase from Streptomyces olivaceoviridis E-86 and Cex xylanase from Cellulomonas fimi is only 50%, an active chimeric enzyme was obtained by replacing module 10 in FXYN with module 10 from Cex. In the family F/10 xylanases, module 10 is an important region as it includes an acid/base catalyst and a substrate binding residue. In FXYN, module 10 consists of 15 amino acid residues, while in Cex it consists of 14 amino acid residues. The Km and kcat values of the chimeric xylanase FCF-C10 for PNP-xylobioside (PNP-X2) were 10-fold less than those for FXYN. CD spectral data indicated that the structure of the chimeric enzyme was similar to that of FXYN. Based on the comparison of the amino acid sequences of FXYN and Cex in module 10, we constructed four mutants of FXYN. When D133 or S135 of FXYN was deleted, the kinetic properties were not changed from those of FXYN. By deletion of both D133 and S135, the Km value for PNP-X2 decreased from the 2.0 mM of FXYN to 0.6 mM and the kcat value decreased from the 20 s(-1) of FXYN to 8.7 s(-1). Insertion of Q140 into the doubly deleted mutant further reduced the Km value to 0.3 mM and the kcat value to 3.8 s(-1). These values are close to those for the chimeric enzyme FCF-C10. These results indicate that module 10 itself is able to accommodate changes in the sequence position of amino acids which are critical for enzyme function. Since changes of the spatial position of these amino acids would be expected to result in enzyme inactivation, module 10 must have some flexibility in its tertiary structure. The structure of module 10 itself also affects the substrate specificity of the enzyme.

Module-intron correlation and intron sliding in family F/10 xylanase genes

Xylanases are classified into two families, numbered F/10 and G/11 according to the similarity of amino acid sequences of their catalytic domain (Henrissat, B., Bairoch, A., 1993. New families in the classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem. J. 293, 781-788). Three-dimensional structure of the catalytic domain of the family F/10 xylanase was reported (White, A., Withers, S.G., Gilkes, N.R., Rose, D.R., 1994. Crystal structure of the catalytic domain of the beta-1,4-glycanase Cex from Cellulomonas fimi. Biochemistry 33, 12546-12552). The domain was decomposed into 22 modules by centripetal profiles (Go, M., Nosaka, M., 1987. Protein architecture and the origin of introns. Cold Spring Harbor Symp. Quant. Biol. 52, 915-924; Noguti, T., Sakakibara, H., Go, M., 1993. Localization of hydrogen-bonds within modules in barnase. Proteins 16, 357-363). A module is a contiguous polypeptide segment of amino acid residues having a compact conformation within a globular domain. Collected 31 intron sites of the family F/10 xylanase genes from fungus were found to be correlated to module boundaries with considerable statistical force (p values <0.001). The relationship between the intron locations and protein structures provides supporting evidence for the ancient origin of introns, because such a relationship cannot be expected by random insertion of introns into eukaryotic genes, but it rather suggests pre-existence of introns in the ancestral genes of prokaryotes and eukaryotes. A phylogenetic tree of the fungal and bacterial xylanase sequences made two clusters; one includes both the bacterial and fungal genes, but the other consists of only fungal genes. The mixed cluster of bacterial genes without introns and the fungal genes with introns further supports the ancient origin of introns. Comparison of the conserved base sequences of introns indicates that sliding of a splice site occurred in Aspergillus kawachii gene by one base from the ancestral position. Substrate-binding sites of xylanase are localized on eight modules, and introns are found at both termini of six out of these functional modules. This result suggests that introns might play a functional role in shuffling the exons encoding the substrate-binding modules.

Putative mechanism of natural transformation as deduced from genome data

Genetic transformation is widely utilized in molecular biology as a tool for gene cloning in Escherichia coli and for gene mapping in Bacillus subtilis. Several strains of eubacteria can naturally take up exogenous DNA and integrate the DNA into their own genomes. Molecular details of natural transformation, however, remained to be elucidated. The complete genome of a cyanobacterium, Synechocystis sp. PCC6803, has been sequenced. This bacterium has been used to examine functions of a particular gene. The genome is considered to carry information on natural transformable characteristics of Synechocystis. The first step in genetic transformation is the uptake of exogenous DNA. Proteins with non-specific DNA binding features are required, because specificity in the exogenous DNA has not been demonstrated. Such proteins have modules interacting with the phosphate backbone of DNA, including helix-turn-helix modules. Using a consensus pattern of the phosphate-binding helix-turn-helix module, we searched through the genome data of Synechocystis for genes or open reading frame (ORF) products with the pattern in primary structures. We found that an ORF, slr0197, has the pattern in duplicate at the C-terminal region. We also found that the ORF product has a hydrophobic segment at the N-terminal region, which is followed by two internal repeats of the endonuclease domain. Based on these observations, we propose a model for the initial stage of genetic transformation. This is apparently the first report on molecular mechanisms of natural transformation.

Repetitive use of a phosphate-binding module in DNA polymerase beta, Oct-1 POU domain and phage repressors

Motifs for sequence specific-protein-DNA interactions, such as helix-turn-helix, zinc finger and leucine zipper, are now better understood as a result of extensive studies of three-dimensional (3D) structures of transcription factors. On the other hand, little attention has been paid to motifs for sequence nonspecific binding, namely DNA-phosphate binding. To address the question whether different transcription factors and DNA manipulation enzymes, that is enzymes that work on DNA, share a similar mode of phosphate binding, we surveyed interactions between DNA and protein module, a structural unit of a globular protein. We analyzed the modular organization of DNA polymerase beta and found that residues making contact with DNA phosphates were localized to five modules. Structural comparison of these phosphate-binding modules against others in transcription factors and DNA manipulation enzymes revealed that DNA polymerase beta, the Oct-1 POU domain, 434 Cro and the Arc repressor have a phosphate-binding module with 3D structures similar to one another. This newly detected module, the phosphate-binding helix-turn-helix (pbHTH) module, named for its function and 3D structure, interacts with DNA by (i) making hydrogen bonds between a DNA phosphodiester oxygen and an amino hydrogen of the main chain located at the N-terminus of a C-terminal alpha-helix, and (ii) making electrostatic interactions between DNA phosphates and side chains of lysine or arginine. Finding structurally and functionally similar phosphate-binding units in different transcription factors and DNA manipulation enzymes suggests that shuffling of modules is not limited to the DNA base-recognition motif. Phosphate-binding modules are apparently also shuffled in DNA-binding proteins.

Molecular analysis of the dnaK locus of Leptospira interrogans serovar Copenhageni

Analysis of the dnaK locus of Leptospira interrogans serovar Copenhageni identified four genes in the order hrcA, grpE, dnaK and dnaJ. This is the first time a homologue of hrcA has been identified in a spirochete. The hrcA gene and a regulatory sequence, designated CIRCE, play a significant role in the regulation of the dnaK locus of several Gram+ organisms. Their presence upstream of dnaK in Leptospira suggested a similiar regulatory mechanism. Transcriptional analysis using reverse transcriptase-PCR demonstrated transcription of all four genes and indicated that hrcA and grpE were co-transcribed, as were grpE and dnaK. Whilst hrcA, grpE and dnaK were closely linked on the chromosome, transcription terminators between dnaK and dnaJ and downstream of dnaJ suggested that this latter gene exists in its own operon. Primer extension analysis located functional promoters upstream of hrcA and grpE; however, no evidence of a functional promoter could be found for dnaJ. Moreover, transcripts encompassing the first three genes or the entire locus could not be demonstrated, suggesting that the four genes are regulated independently at the transcriptional level. These results indicate that the regulation of the dnaK locus of Leptospira differs somewhat from that observed in other organisms.

Structure—activity relationships of some indolo[3,2-c]quinolines with antimalarial activity

The synthesis, physicochemical characterization and in vitro antimalarial activity of a series of indolo[3,2-c]quinolines (9a-f) are described. There is only a poor correlation between the activity and hydrophobicity. In contrast, 33% of the observed variation in antimalarial activity can be attributed to the size of the side chain attached to position 9 of the indoloquinoline ring. An increase in the size of this dibasic side chain generally results in a reduction in activity, suggesting that it is accommodated in a site/cavity of limited size on the receptor. More significantly, the charge on the distal nitrogen (N3) on the side chain, located 10-11 A from the quinoline N, could account for 75% of the observed variation. Since a large charge on N3 is associated with improved antimalarial activity, it is suggested that N3 is protonated and functions as a H bond donor in the drug-receptor interaction.

Mechanical stability of compact modules of barnase

Globular proteins are composed of structural elements such as secondary structures and modules. Modules are compact segments consisting of 10-40 contiguous amino acid residues and are often encoded by exons. Therefore, the view that the modular organization of proteins is a result of exon-shuffling or -fusion is given support. Secondary structures such as alpha-helix and beta-sheet are stabilized by hydrogen bonds and are thus considered to be stable, structural elements of a globular domain. Since module boundaries are often located on alpha-helices or beta-sheets, it is not obvious whether the modules are mechanically stable. We carried out molecular dynamics simulations on modules of barnase, a bacterial RNase from Bacillus amyloliquefaciens, for 1 ns in vacuo and 150 ps in water. Five of six modules (M1, 1-24; M2, 25-52; M3, 53-73; M4, 74-88; M5, 89-98) retained native-like conformations during these simulations. Only the C-terminal module (M6, 99-110) was deformed; it is less compact than the other modules. As the modules are mechanically stable they are suitable as parts combined into proteins. Together with RNase activity of the three isolated modules of barnase, M2, M3 and M6, our study supports the view that modules were indeed original building blocks of proteins.

Ligase chain reaction assay for human mutations: the Sickle Cell by LCR assay

We can detect the beta-globin gene sickle cell mutation by using an assay based on the ligase chain reaction. The simultaneous amplification of the human growth hormone gene in the same reaction serves as a control for the amount of template DNA or amplification efficiency. Ligation products, which are biotinylated at one end and tagged with an arbitrary "tail" sequence at the other, are captured by hybridization to "tail"-complementary oligonucleotides immobilized on polystyrene microwells. The captured ligation products are detected colorimetrically by use of streptavidin-alkaline phosphatase conjugate. In a study of 24 subjects, the assay unequivocally discriminated among normal, carrier, and sickle cell genotypes.

Adaptive amino acid replacements accompanied by domain fusion in reverse transcriptase

Two basic processes are involved in protein evolution: One is amino acid replacement and another is reorganization of structural or functional units of proteins. Multidomain or multifunctional proteins are thought to have evolved by fusion of smaller structural units such as modules or domains. Reverse transcriptase (RT) is one of such fused proteins. The N-terminal part forms of globular domain with polymerase activity and the C-terminal part forms another globular domain with ribonuclease H activity (RNase H domain). There are single-domain enzymes which are homologous with the RNase H domain. The group of enzymes is called type I ribonuclease H (RNase HI). It is most likely that the ancestors of RNase HI and the polymerase domain were fused and became contemporary RT. At fusion, amino acid replacements presumably occurred at the interface of the domains to reinforce the interdomain interactions. Such replaced amino acid residues are conserved during evolution of the fused enzyme. We analyzed the pattern of amino acid replacement at each residue site in the free form, RNase HI group, and the integrated form, RNase H domain group. Then we compared the patterns between the two forms. Drastic fitting replacements of amino acid residues occurred at four of 29 residue sites involved in interdomain contact. Hydrophilic amino acid residues of the free form were substituted with hydrophobic or ambivalent ones in the integrated form. These substitutions aid in stabilizing the fused conformation by hydrophobic interactions at the interface of the domains. These observations imply that domain fusion could have occurred with only a relatively small number of adaptive amino acid substitutions.