Rapid knot detection and application to protein structure prediction

MOTIVATION: Knots in polypeptide chains have been found in very few proteins, and consequently should be generally avoided in protein structure prediction methods. Most effective structure prediction methods do not model the protein folding process itself, but rather seek only to correctly obtain the final native state. Consequently, the mechanisms that prevent knots from occurring in native proteins are not relevant to the modeling process, and as a result, knots can occur with significantly higher frequency in protein models. Here we describe Knotfind, a simple algorithm for knot detection that is fast enough for structure prediction, where tens or hundreds of thousands of conformations may be sampled during the course of a prediction. We have used this algorithm to characterize knots in large populations of model structures generated for targets in CASP 5 and CASP 6 using the Rosetta homology-based modeling method. RESULTS: Analysis of CASP5 models suggested several possible avenues for introduction of knots into these models, and these insights were applied to structure prediction in CASP 6, resulting in a significant decrease in the proportion of knotted models generated. Additionally, using the knot detection algorithm on structures in the Protein Data Bank, a previously unreported deep trefoil knot was found in acetylornithine transcarbamylase. AVAILABILITY: The Knotfind algorithm is available in the Rosetta structure prediction program at http://www.rosettacommons.org.     

Rapid removal of acetimidoyl groups from proteins and peptides. Applications to primary structure determination.

Methylamine buffers can be used for the rapid quantitative removal of acetimidoyl groups from proteins and peptides modified by treatment with ethyl or methyl acetimidate. The half-life for displacement of acetimidoyl groups from fully amidinated proteins incubated in 3.44 M-methylamine/HCl buffer at pH 11.5 and 25 degrees C was approx. 26 min; this half life is 29 times less than that observed in ammonia/HCl buffer under the same conditions of pH and amine concentration. Incubation of acetimidated proteins with methylamine for 4 h resulted in greater than 95% removal of acetimidoyl groups. No deleterious effects on primary structure were detected by amino acid analysis or by automated Edman degradation. Reversible amidination of lysine residues, in conjunction with tryptic digestion, has been successfully applied to the determination of the amino acid sequence of an acetimidated mouse immunoglobulin heavy chain peptide. The regeneration of amino groups in amidinated proteins and peptides by methylaminolysis makes amidination a valuable alternative to citraconoylation and maleoylation in structural studies.     

Refractive index-based determination of detergent concentration and its application to the study of membrane proteins

The concentration of detergent in membrane protein preparations can have a critical role on protein stability, function, and the potential for crystallization. Unfortunately, dialysis or protein concentration can lead to an unknown amount of detergent in the final membrane protein preparations. Here we present a method for the determination of detergent concentration based on refractive index of the detergent solution. This method was applied to quantitate the amount of detergent remaining in solution after concentration in various concentrators. We found that the ability of the tested detergents to pass through the molecular weight cutoff membrane correlates well with detergent micelle size. Therefore, the micelle size can be used as a rough guide to estimate the retention of a given detergent in various molecular weight cutoff concentrators. The refractive index method is exceptionally informative when coupled with size exclusion chromatography and light scattering, and can be used to determine the oligomeric state of the membrane protein, the size of a protein-associated micelle, as well as the amount and size of the unbound detergent micelle.     

Phospholipid component volumes: determination and application to bilayer structure calculations.

We present a new method for the determination of bilayer structure based on a combination of computational studies and laboratory experiments. From molecular dynamics simulations, the volumes of submolecular fragments of saturated and unsaturated phosphatidylcholines in the liquid crystalline state have been extracted with a precision not available experimentally. Constancy of component volumes, both among different lipids and as a function of membrane position for a given lipid, have been examined. The component volumes were then incorporated into the liquid crystallographic method described by Wiener and White (1992. Biophys. J. 61:434-447, and references therein) for determining the structure of a fluid-phase dioleoylphosphatidylcholine bilayer from x-ray and neutron diffraction experiments.     

Purification of BsuE methyltransferase and its application in genome mapping.

We have used a combination of BsuE methyltransferase (M-BsuE) and NotI restriction enzyme to cut genomic DNA at a subset of NotI sites. The usefulness of this system is shown in a re-examination of the restriction map of the human MHC. Combinations of methylases and restriction enzymes can be used to generate cuts at different frequencies in genomic DNA, such that they generate ends complementary to NotI ends, and can be used in conjunction with NotI linking clones in chromosome jumping experiments. These enzyme combinations have the potential to produce cutting sites in genomic DNA spaced at intervals favorable for extensive mapping, fragment enrichment, and cloning efforts.     

Rapid and accurate structure determination of coiled-coil domains using NMR dipolar couplings: application to cGMP-dependent protein kinase Ialpha

Coiled-coil motifs play essential roles in protein assembly and molecular recognition, and are therefore the targets of many ongoing structural and functional studies. However, owing to the dynamic nature of many of the smaller coiled-coil domains, crystallization for X-ray studies is very challenging. Determination of elongated structures using standard NMR approaches is inefficient and usually yields low-resolution structures due to accumulation of small errors over long distances. Here we describe a solution NMR approach based on residual dipolar couplings (RDCs) for rapid and accurate structure determination of coiled-coil dimers. Using this approach, we were able to determine the high-resolution structure of the coiled-coil domain of cGMP-dependent protein kinase Ialpha, a protein of previously unknown structure that is critical for physiological relaxation of vascular smooth muscle. This approach can be extended to solve coiled-coil structures with higher order assemblies.     

The mitochondrial genome of Saprolegnia ferax: organization, gene content and nucleotide sequence

The mitochondrial genome of the peronosporomycete water mold Saprolegnia ferax has been characterized as a 46 930 bp circle containing an 8618 bp large inverted repeat (LIR). Eighteen reading frames encode identified subunits of respiratory complexes I, III, IV and V; 16 encode polypeptides of small and large mitoribosome subunits; and one encodes a subunit of the sec-independent protein translocation pathway. Of four additional putative reading frames three are homologues of those found in the related Phytophthora infestans genome. Protein encoding loci in the tightly compacted genome typically are arranged in operon-like clusters including three abutting and two overlapping pairs of reading frames. Translational RNAs include the mitochondrial small and large subunit rRNAs and 25 tRNA species. No tRNAs are encoded to enable translation of any threonine or the arginine CGR codons. The LIR separates the molecule into 19 274 bp large and 10 420 bp small single copy regions, and it encodes intact duplicate copies of four reading frames encoding known proteins, both rRNAs, and five tRNAs. Partial 3' sequences of three additional reading frames are duplicated at single copy sequence junctions. Active recombination between LIR elements generates two distinctive gene orders and uses the duplicated 3' sequences to maintain intact copies of the partially duplicated loci.     

A hybrid gene team model and its application to genome analysis

It is well-known that functionally related genes occur in a physically clustered form, especially operons in bacteria. By leveraging on this fact, there has recently been an interesting problem formulation known as gene team model, which searches for a set of genes that co-occur in a pair of closely related genomes. However, many gene teams, even experimentally verified operons, frequently scatter within other genomes. Thus, the gene team model should be refined to reflect this observation. In this paper, we generalized the gene team model, that looks for gene clusters in a physically clustered form, to multiple genome cases with relaxed constraints. We propose a novel hybrid pattern model that combines the set and the sequential pattern models. Our model searches for gene clusters with and/or without physical proximity constraint. This model is implemented and tested with 97 genomes (120 replicons). The result was analyzed to show the usefulness of our model. We also compared the result from our hybrid model to those from the traditional gene team model. We also show that predicted gene teams can be used for various genome analysis: operon prediction, phylogenetic analysis of organisms, contextual sequence analysis and genome annotation. Our program is fast enough to provide a service on the web at http://platcom.informatics.indiana.edu/platcom/. Users can select any combination of 97 genomes to predict gene teams.     

The relationship between protein structure and function: a comprehensive survey with application to the yeast genome.

For most proteins in the genome databases, function is predicted via sequence comparison. In spite of the popularity of this approach, the extent to which it can be reliably applied is unknown. We address this issue by systematically investigating the relationship between protein function and structure. We focus initially on enzymes functionally classified by the Enzyme Commission (EC) and relate these to by structurally classified domains the SCOP database. We find that the major SCOP fold classes have different propensities to carry out certain broad categories of functions. For instance, alpha/beta folds are disproportionately associated with enzymes, especially transferases and hydrolases, and all-alpha and small folds with non-enzymes, while alpha+beta folds have an equal tendency either way. These observations for the database overall are largely true for specific genomes. We focus, in particular, on yeast, analyzing it with many classifications in addition to SCOP and EC (i.e. COGs, CATH, MIPS), and find clear tendencies for fold-function association, across a broad spectrum of functions. Analysis with the COGs scheme also suggests that the functions of the most ancient proteins are more evenly distributed among different structural classes than those of more modern ones. For the database overall, we identify the most versatile functions, i.e. those that are associated with the most folds, and the most versatile folds, associated with the most functions. The two most versatile enzymatic functions (hydro-lyases and O-glycosyl glucosidases) are associated with seven folds each. The five most versatile folds (TIM-barrel, Rossmann, ferredoxin, alpha-beta hydrolase, and P-loop NTP hydrolase) are all mixed alpha-beta structures. They stand out as generic scaffolds, accommodating from six to as many as 16 functions (for the exceptional TIM-barrel). At the conclusion of our analysis we are able to construct a graph giving the chance that a functional annotation can be reliably transferred at different degrees of sequence and structural similarity. Supplemental information is available from http://bioinfo.mbb.yale.edu/genome/foldfunc++ +.     

A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome

A model of DNA sequence evolution applicable to coding regions is presented. This represents the first evolutionary model that accounts for dependencies among nucleotides within a codon. The model uses the codon, as opposed to the nucleotide, as the unit of evolution, and is parameterized in terms of synonymous and nonsynonymous nucleotide substitution rates. One of the model's advantages over those used in methods for estimating synonymous and nonsynonymous substitution rates is that it completely corrects for multiple hits at a codon, rather than taking a parsimony approach and considering only pathways of minimum change between homologous codons. Likelihood-ratio versions of the relative-rate test are constructed and applied to data from the complete chloroplast DNA sequences of Oryza sativa, Nicotiana tabacum, and Marchantia polymorpha. Results of these tests confirm previous findings that substitution rates in the chloroplast genome are subject to both lineage-specific and locus-specific effects. Additionally, the new tests suggest tha the rate heterogeneity is due primarily to differences in nonsynonymous substitution rates. Simulations help confirm previous suggestions that silent sites are saturated, leaving no evidence of heterogeneity in synonymous substitution rates.      

HPLC fluorescence method for the determination of nizatidine in human plasma and its application to pharmacokinetic study

A sensitive high‐performance liquid chromatographic (HPLC) method was developed for the determination of nizatidine in human plasma. Nizatidine was derivatized by 4‐fluoro‐7‐nitrobenzofurazan (NBD‐F). Chromatographic separation was performed on a Inertsil C₁₈ column (150 mm × 4.6 mm, 5 µm) using isocratic elution by a mobile phase consisting of methanol/water (55:45) at a flow rate of 1.2 mL/min. Amlodipine was used as the internal standard (IS). Fluorescence detector was used operated at 461 nm (excitation) and 517 nm (emission), respectively. The calibration curve was linear over the range of 50–2000 ng/mL. This method was successfully applied to a pharmacokinetic study after oral administration of a dose (150 mg) of nizatidine. Copyright © 2013 John Wiley & Sons, Ltd.     

Complete nucleotide sequence and genome organization of bovine parvovirus.

We determined the complete nucleotide sequence of bovine parvovirus (BPV), an autonomous parvovirus. The sequence is 5,491 nucleotides long. The terminal regions contain nonidentical imperfect palindromic sequences of 150 and 121 nucleotides. In the plus strand, there are three large open reading frames (left ORF, mid ORF, and right ORF) with coding capacities of 729, 255, and 685 amino acids, respectively. As with all parvoviruses studied to date, the left ORF of BPV codes for the nonstructural protein NS-1 and the right ORF codes for the major parts of the three capsid proteins. The mid ORF probably encodes the major part of the nonstructural protein NP-1. There are promoterlike sequences at map units 4.5, 12.8, and 38.7 and polyadenylation signals at map units 61.6, 64.6, and 98.5. BPV has little DNA homology with the defective parvovirus AAV, with the human autonomous parvovirus B19, or with the other autonomous parvoviruses sequenced (canine parvovirus, feline panleukopenia virus, H-1, and minute virus of mice). Even though the overall DNA homology of BPV with other parvoviruses is low, several small regions of high homology are observed when the amino acid sequences encoded by the left and right ORFs are compared. From these comparisons, it can be shown that the evolutionary relationship among the parvoviruses is B19 in equilibrium with AAV in equilibrium with BPV in equilibrium with MVM. The highly conserved amino acid sequences observed among all parvoviruses may be useful in the identification and detection of parvoviruses and in the design of a general parvovirus vaccine.     

Combinatorial minimization and secondary structure determination of a nucleotide synthase ribozyme

We previously isolated from random sequences ribozymes able to form a glycosidic linkage between a ribose sugar and 4-thiouracil in a reaction that mimics protein-catalyzed nucleotide synthesis. Here we report on two serial in vitro selection experiments that defined the core motif of one of the nucleotide synthase ribozymes and provided improved versions of this ribozyme. The first selection experiment started from a degenerate sequence pool based on the previously isolated sequence and used a selection-amplification protocol that allowed the sequence requirements at the 3' terminus of the ribozyme to be interrogated. Comparing the active sequences identified in this experiment revealed the complicated secondary structure of the nucleotide synthase ribozyme. A second selection was then performed to remove nonessential sequence from the ribozyme. This selection started with a pool with variation introduced in both the sequence and the length of the nonconserved loops and joining regions. This pool was generated using a partial reblocking/deblocking strategy on a DNA synthesizer, allowing the combinatorial synthesis of both point deletions and point substitutions. The consensus ribozyme motif that emerged was an approximately 71 nt pseudoknot structure with five stems and two important joining segments. Comparative sequence analysis and a cross-linking experiment point to the probable location of nucleotide synthesis. The prototype isolate from the second selection was nearly 35 times more efficient than the initial isolate and at least 10(8) times more efficient than an upper limit of an as-yet undetectable uncatalyzed reaction, supporting the idea that RNA-catalyzed nucleotide synthesis might have been important in an RNA world.     

The nucleotide sequence and restriction enzyme sites of the polyoma genome.

We present the 5295 nucleotide-long sequence of the polyoma genome and the restriction enzyme digestion sites predicted from this sequence.     

Chromatid structure: relationship between DNA content and nucleotide sequence diversity

Models of chromatid structure are based on inferences made from genetic, cytological, and cytochemical observations. An alternative approach can provide limits as to the number of identical subunits present in chromatids. This method is based on the demonstration that nucleotide sequence diversity may be estimated from the kinetics of renaturation of denatured DNA. Measurements of DNA content and renaturation rate constants are given for several eukaryotic DNAs. Control experiments involved measurements of renaturation kinetics of DNAs from bacteria and bacteriophage. These estimates show that most of the nucleotide sequences in mouse, Drosophila, and Ciona DNA are present only once per sperm. Since the reduction of DNA content during meiosis indicates that mouse sperm contain a haploid set of chromatids, it follows that a set of mouse meiotic chromatids contains a single copy of most sequences. Models of chromatid structure which postulate multiple subunits with identical nucleotide sequences are therefore not tenable for mouse meiotic chromatids. This method of analyzing nucleotide sequence diversity may be of general use in designing models of chromatid structure in other organisms.