A flexible new computer program for handling DNA sequence data.

A compact new computer program for handling nucleic acid sequence data is presented. It consists of a number of different subsets, which may be used according to a given code system. The program is designed for the determination of restriction enzyme and other recognition sites in correlation with translation patterns, and allows tabulation of codon frequencies and protein molecular weights within specified gene boundaries. The program is especially designed for detection of overlapping genes. The language, is FORTRAN and thus the program may be used on small computers; it may also be used without any prior computer experience. Copies are available on request.     

Construction of restriction maps

A computer program is described, which constructs maps of restrictionendonuclease cleavage sites in linear or circular DNA molecules,given the fragment lengths in single and double digestions withtwo enzymes. The algorithm is based upon a partition methodand a very simple rule to chain fragments. The program is writtenin Prolog II. Received on July 28, 1987; accepted on December 31, 1987     

Definition and identification of homology domains

A method is described for identifying and evaluating regionsof significant similarity between two sequences. The notionof a ‘homology domain’ is employed which definesthe boundaries of a region of sequence homology containing noinsertions or deletions. The relative significance of differentpotential homology domains is evaluated using a non-linear similarityscore related to the probability of finding the observed levelof similarity in the region by chance. The sensitivity of themethod is demonstrated by simulating the evolution of homologydomains and applying the method to their detection. Severalexamples of the use of homology domain identification are given. Received on July 29, 1987; accepted on November 15, 1987     

β-Strand-mediated interactions of protein domains

Protein domains exist by themselves or in combination with other domains to form complex multidomain proteins. Defining domain boundaries in proteins is essential for understanding their evolution and function but is not trivial. More specifically, partitioning domains that interact by forming a single β-sheet is known to be particularly troublesome for automatic structure-based domain decomposition pipelines. Here, we study edge-to-edge β-strand interactions between domains in a protein chain, to help define the boundaries for some more difficult cases where a single β-sheet spanning over two domains gives an appearance of one. We give a number of examples where β-strands belonging to a single β-sheet do not belong to a single domain and highlight the difficulties of automatic domain parsers on these examples. This work can be used as a baseline for defining domain boundaries in homologous proteins or proteins with similar domain interactions in the future.     

Evaluation of Disulfide Bond Position to Enhance the Thermal Stability of a Highly Stable Single Domain Antibody

Single domain antibodies are the small recombinant variable domains derived from camelid heavy-chain-only antibodies. They are renowned for their stability, in large part due to their ability to refold following thermal or chemical denaturation. In addition to refolding after heat denaturation, A3, a high affinity anti-Staphylococcal Enterotoxin B single domain antibody, possesses a melting temperature of ∼84°C, among the highest reported for a single domain antibody. In this work we utilized the recently described crystal structure of A3 to select locations for the insertion of a second disulfide bond and evaluated the impact that the addition of this second bond had on the melting temperature. Four double-disulfide versions of A3 were constructed and each was found to improve the melting temperature relative to the native structure without reducing affinity. Placement of the disulfide bond at a previously published position between framework regions 2 and 3 yielded the largest improvement (>6°C), suggesting this location is optimal, and seemingly provides a universal route to raise the melting temperature of single domain antibodies. This study further demonstrates that even single domain antibodies with extremely high melting points can be further stabilized by addition of disulfide bonds.     

Natural recombination event within the capsid genomic region leading to a chimeric strain of human enterovirus B

Recombination between two strains is a known phenomenon for enteroviruses replicating within a single cell. We describe a recombinant strain recovered from human stools, typed as coxsackievirus B4 (CV-B4) and CV-B3 after partial sequencing of the VP1 and VP2 coding regions, respectively. The strain was neutralized by a polyclonal CV-B3-specific antiserum but not by a CV-B4-specific antiserum. The nucleotide sequence analysis of the whole structural genomic region showed the occurrence of a recombination event at position 1950 within the VP3 capsid gene, in a region coding for the 2b antigenic site previously described for CV-B3. This observation evidences for the first time the occurrence of an interserotypic recombination within the VP2-VP3-VP1 capsid region between two nonpoliovirus enterovirus strains. The neutralization pattern suggests that the major antigenic site is located within the VP2 protein.     

Automatic construction of restriction site maps.

A computer program is described which constructs maps of restriction endonuclease cleavage sites in DNA molecules, given only the fragment lengths. The program utilizes fragment length data from single and double restriction enzyme digests to generate maps for linear or circular molecules. The search for a map can be limited to the unknown (insert) region of a recombinant phage or plasmid. Typical restriction maps with four or five enzymes which cut at three to five unknown sites can be calculated in a few minutes.     

Structure determination of Discoidin II from Dictyostelium discoideum and carbohydrate binding properties of the lectin domain

The social amoeba Dictyostelium discoideum adopts a cohesive stage upon starvation and then produces Discoidin I and II, two proteins able to bind galactose and N-acetyl-galactosamine. The N-terminal domain or discoidin domain (DS) is widely distributed in eukaryotes where it plays a role in extracellular matrix binding while the C-terminal domain displays sequence similarities to invertebrate lectins. We present the first X-ray structures of the wild-type and recombinant Discoidin II in unliganded state and in complex with monosaccharides. The protein forms a homotrimer which presents two binding surfaces situated on the opposite boundaries of the structure. The binding sites of the N-terminal domain contain PEG molecules that could mimics binding of natural ligand. The C-terminal lectin domain interactions with N-acetyl-D-galactosamine and methyl-beta-galactoside are described. The carbohydrate binding sites are located at the interface between monomers. Specificity for galacto configuration can be rationalized since the axial O4 hydroxyl group is involved in several hydrogen bonds with protein side chains. Titration microcalorimetry allowed characterization of affinity and demonstrated the enthalpy-driven character of the interaction. Those results highlight the structural differentiation of the DS domain involved in many cell-adhesion processes from the lectin activity of Dictyostelium discoidins.     

Accessing single nucleotide polymorphisms in genomic DNA by direct multiplex polymerase chain reaction amplification on oligonucleotide microarrays

This study introduces a DNA microarray-based genotyping system for accessing single nucleotide polymorphisms (SNPs) directly from a genomic DNA sample. The described one-step approach combines multiplex amplification and allele-specific solid-phase PCR into an on-chip reaction platform. The multiplex amplification of genomic DNA and the genotyping reaction are both performed directly on the microarray in a single reaction. Oligonucleotides that interrogate single nucleotide positions within multiple genomic regions of interest are covalently tethered to a glass chip, allowing quick analysis of reaction products by fluorescence scanning. Due to a fourfold SNP detection approach employing simultaneous probing of sense and antisense strand information, genotypes can be automatically assigned and validated using a simple computer algorithm. We used the described procedure for parallel genotyping of 10 different polymorphisms in a single reaction and successfully analyzed more than 100 human DNA samples. More than 99% of genotype data were in agreement with data obtained in control experiments with allele-specific oligonucleotide hybridization and capillary sequencing. Our results suggest that this approach might constitute a powerful tool for the analysis of genetic variation.     

A method for identifying splice sites and translation start sites in human genomic sequences

We describe a new method for identifying the sequences that signal the start of translation, and the boundaries between exons and introns (donor and acceptor sites) in human mRNA. According to the mandatory keyword, ORGANISM, and feature key, CDS, a large set of standard data for each signal site was extracted from the ASCII flat file, gbpri.seq, in the GenBank release 108.0. This was used to generate the scoring matrices, which summarize the sequence information for each signal site. The scoring matrices take into account the independent nucleotide frequencies between adjacent bases in each position within the signal site regions, and the relative weight on each nucleotide in proportion to their probabilities in the known signal sites. Using a scoring scheme that is based on the nucleotide scoring matrices, the method has great sensitivity and specificity when used to locate signals in uncharacterized human genomic DNA. These matrices are especially effective at distinguishing true and false sites.     

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.     

Algorithms for identifying local molecular sequence features

Efficient algorithms are described for identifying local molecularsequence features including repeats, dyad symmetry pairingsand aligned matches between sequences, while allowing for errors.Specific applications are given to the genomic sequences ofthe Epstein-Barr virus, Varicella-Zoster virus and the bacteriophages and T7. Received on October 6, 1987; accepted on December 13, 1987     

Functional gene expression domains: defining the functional unit of eukaryotic gene regulation

The term functional domain is often used to describe the region containing the cis acting sequences that regulate a gene locus. "Strong" domain models propose that the domain is a spatially isolated entity consisting of a region of extended accessible chromatin bordered by insulators that have evolved to act as functional boundaries. However, the observation that independently regulated loci can overlap partially or completely raises questions about functional requirements for physically isolated domain structures. An alternative model, the "weak" domain model, proposes that domain structure is determined by the distribution of binding sites for positively acting factors, without a requirement for functional boundaries. The domain would effectively be the region that contains these factor-binding sites. Specificity of promoter-enhancer interactions would play a major role in maintaining the functional autonomy of adjacent genes. Sequences that interfere with these interactions (frequently characterised as insulators) would be selected against if they occurred within the domain but not at the edges, or in the interdomain regions. As a result, insulators would often be found near the borders of domains without necessarily being selected to act as boundaries.     

An integrated software system for microcomputer management of recombinant DNA data.

A database-oriented system (pCP123) is described for the manipulation of recombinant DNA data. This system was developed within the context of an integrated software package with spreadsheet, database, graphing and programming capabilities. The system includes two databases, one of sites and another of regions, coordinately handled by a series of macro-programs operated from four user-define menus. A distinctive feature of the system is the possibility of handling both ends of defined functional or structural regions in situations of simulated deletions or insertions.     

利用酵母人工染色体（YAC）减数分裂同源重组构建人免疫球蛋白κ链基因组大片段

具有同源重叠区的酵母人工染色体(YAC)可以利用酵母细胞减数分裂进行同源重组,从而构建更大的人工染色体基因组,这对生命科学基础研究和生物技术应用研究有着非常重要的意义。本实验以两个含人免疫球蛋白κ链基因簇片段的YAC克隆为材料,通过酵母改型、异型接合、二倍体发孢、单孢子筛选和分子生物学鉴定等技术和方法,利用酵母菌减数分裂同源重组机制,构建了一条包含人的免疫球蛋白κ轻链32个Vκ基因、5个Jκ基因、Cκ基因、Eκ基因和κde基因的YAC重组体,长度约400kb。同时,本实验利用溶壁酶消化法获取单孢子重组体,代替了传统的显微分孢操作。使得利用酵母人工染色体减数分裂同源重组的技术更加简便可行。