STAR: predicting recombination sites from amino acid sequence

Background: Designing novel proteins with site-directed recombination has enormous prospects. By locating effective recombination sites for swapping sequence parts, the probability that hybrid sequences have the desired properties is increased dramatically. The prohibitive requirements for applying current tools led us to investigate machine learning to assist in finding useful recombination sites from amino acid sequence alone. Results: We present STAR, Site Targeted Amino acid Recombination predictor, which produces a score indicating the structural disruption caused by recombination, for each position in an amino acid sequence. Example predictions contrasted with those of alternative tools, illustrate STAR'S utility to assist in determining useful recombination sites. Overall, the correlation coefficient between the output of the experimentally validated protein design algorithm SCHEMA and the prediction of STAR is very high (0.89). Conclusion: STAR allows the user to explore useful recombination sites in amino acid sequences with unknown structure and unknown evolutionary origin. The predictor service is available from http://pprowler.itee.uq.edu.au/star.

The effect of box shape on the dynamic properties of proteins simulated under periodic boundary conditions

The effect of the box shape on the dynamic behavior of proteins simulated under periodic boundary conditions is evaluated. In particular, the influence of simulation boxes defined by the near-densest lattice packing (NDLP) in conjunction with rotational constraints is compared to that of standard box types without these constraints. Three different proteins of varying size, shape, and secondary structure content were examined in the study. The statistical significance of differences in RMSD, radius of gyration, solvent-accessible surface, number of hydrogen bonds, and secondary structure content between proteins, box types, and the application or not of rotational constraints has been assessed. Furthermore, the differences in the collective modes for each protein between different boxes and the application or not of rotational constraints have been examined. In total 105 simulations were performed, and the results compared using a three-way multivariate analysis of variance (MANOVA) for properties derived from the trajectories and a three-way univariate analysis of variance (ANOVA) for collective modes. It is shown that application of roto-translational constraints does not have a statistically significant effect on the results obtained from the different simulations. However, the choice of simulation box was found to have a small (5-10%), but statistically significant effect on the behavior of two of the three proteins included in the study. (c) 2005 Wiley Periodicals, Inc.

Predicting structural disruption of proteins caused by crossover

We present a machine learning model that predicts a structural disruption score from a protein’s primary structure. SCHEMA was introduced by Frances Arnold and colleagues as a method for determining putative recombination sites of a protein on the basis of the full (PDB) description of its structure. The present method provides an alternative to SCHEMA that is able to determine the same score from sequence data only. Circumventing the need for resolving the full structure enables the exploration of yet unresolved and even hypothetical sequences for protein design efforts. Deriving the SCHEMA score from a primary structure is achieved using a two step approach: first predicting a secondary structure from the sequence and then predicting the SCHEMA score from the predicted secondary structure. The correlation coefficient for the prediction is 0.88 and indicates the feasibility of replacing SCHEMA with little loss of precision. ©2005 IEEE

Searching oligo sets of human chromosome 12 using evolutionary strategies

DNA microarray is a powerful tool to measure the level of a mixed population of nucleic acids at one time, which has great impact in many aspects of life sciences research. In order to distinguish nucleic acids with very similar composition by hybridization, it is necessary to design probes with high specificities, i.e. uniqueness, and also sensitivities, i.e., suitable melting temperature and no secondary structure. We make use of available biology tools to gain necessary sequence information of human chromosome 12, and combined with evolutionary strategy (ES) to find unique subsequences representing all predicted exons. The results are presented and discussed.

Direct protein interaction underlies gene-for-gene specificity and coevolution of the flax resistance genes and flax rust avirulence genes

Plant resistance proteins (R proteins) recognize corresponding pathogen avirulence (Avr) proteins either indirectly through detection of changes in their host protein targets or through direct R-Avr protein interaction. Although indirect recognition imposes selection against Avr effector function, pathogen effector molecules recognized through direct interaction may overcome resistance through sequence diversification rather than loss of function. Here we show that the flax rust fungus AvrLS67 genes, whose products are recognized by the L5, L6, and L7 R proteins of flax, are highly diverse, with 12 sequence variants identified from six rust strains. Seven AvrL567 variants derived from Avr alleles induce necrotic responses when expressed in flax plants containing corresponding resistance genes (R genes), whereas five variants from avr alleles do not. Differences in recognition specificity between AvA567 variants and evidence for diversifying selection acting on these genes suggest they have been involved in a gene-specific arms race with the corresponding flax R genes. Yeast two-hybrid assays indicate that recognition is based on direct R-Avr protein interaction and recapitulate the interaction specificity observed in planta. Biochemical analysis of Escherichia coli-produced AvrL567 proteins shows that variants that escape recognition nevertheless maintain a conserved structure and stability, suggesting that the amino acid sequence differences directly affect the R-Avr protein interaction. We suggest that direct recognition associated with high genetic diversity at corresponding R and Avr gene loci represents an alternative outcome of plant-pathogen coevolution to indirect recognition associated with simple balanced polymorphisms for functional and nonfunctional R and Avr genes.

Homomeric ring assemblies of eukaryotic Sm proteins have affinity for both RNA and DNA - Crystal structure of an oligomeric complex of yeast SmF

Sm and Sm-like proteins are key components of small ribonucleoproteins involved in many RNA and DNA processing pathways. In eukaryotes, these complexes contain seven unique Sm or Sm-like (Lsm) proteins assembled as hetero-heptameric rings, whereas in Archaea and bacteria six or seven-membered rings are made from only a single polypeptide chain. Here we show that single Sm and Lsm proteins from yeast also have the capacity to assemble into homo-oligomeric rings. Formation of homo-oligomers by the spliceosomal small nuclear ribonucleoprotein components SmE and SmF preclude hetero-interactions vital to formation of functional small nuclear RNP complexes in vivo. To better understand these unusual complexes, we have determined the crystal structure of the homomeric assembly of the spliceosomal protein SmF. Like its archaeal/bacterial homologs, the SmF complex forms a homomeric ring but in an entirely novel arrangement whereby two heptameric rings form a co-axially stacked dimer via interactions mediated by the variable loops of the individual SmF protein chains. Furthermore, we demonstrate that the homomeric assemblies of yeast Sm and Lsm proteins are capable of binding not only to oligo(U) RNA but, in the case of SmF, also to oligo(dT) single-stranded DNA.

When oceans meet: a teleost shows secondary intergradation at an Indian-Pacific interface

The Indo-West Pacific is characterized by extraordinary marine species diversity. The evolutionary mechanisms responsible for generating this diversity remain puzzling, but are often linked to Pleistocene sea level fluctuations. The impact of these sea level changes on the population genetic architecture of the estuarine fish Lates calcarifer are investigated via a natural experiment in a region of the Indo-West Pacific known to have undergone considerable change during the Pleistocene. L. calcarifer, a coastline-restricted catadromous teleost, provides an excellent model for studying the effects of sea level change as its habitat requirements potentially make it sensitive to the region's physical history. Evidence was found for a large phylogenetic break (4% mtDNA control region; 0.47% ATPase 6 and 8) either side of the Torres Strait, which separates the Western Pacific and Indian Oceans, although some mixing of the clades was evident. This suggests clinal secondary introgression of the clades via contemporary gene flow. Further, populations on Australia's east coast appear to have passed through a bottleneck. This was linked to the historical drying of the Great Barrier Reef coastal lagoon, which resulted in a significant loss of habitat and forced retreat into isolated refugia. These results suggest that historical eustatic changes have left a significant imprint on the molecular diversity within marine species as well as among them in the Indo-West Pacific.

Complexes of 1-thia-4,7,10-triazacyclododecane ([12]aneN(3)S) - Synthesis, structure and stability constants

The 12-membered macrocyclic ligand 1-thia-4,7, 10-triazacyclododecane ([12]aneN(3)S) has been synthesised, although upon crystallization from acetonitrile a product in which carbon dioxide had added to one secondary amine in the macrocyclic ring (H[12]aneN(3)SCO(2). H2O) was isolated and subsequently characterised by X-ray crystallography. The protonation constants for [12]aneN(3)S and stability constants with Zn(II), Pb(II), Cd(II) and Cu(II) have been determined either potentiometrically or spectrophotometrically in aqueous solution, and compared with those measured or reported for the ligands 1-oxa-4,7,10-triazacyclododecane ([12]aneN(3)O) and 1,4,7,10-tetraazacyclododecane ([12]aneN(4)). The magnitudes of the stability constants are consistent with trends observed previously for macrocyclic ligands as secondary amine donors are replaced with oxygen and thioether donors although the stability constant for the [Hg([12]aneN(4))](2+) complex has been estimated from an NMR experiment to be at least three orders of magnitude larger than reported previously. Zinc(II), mercury(II), lead(II), copper(II) and nickel(II) complexes of [12]aneN(3)S have been isolated and characterised by X-ray crystallography. In the case of copper(II), two complexes [Cu([12]aneN(3)S)(H2O)](ClO4)(2) and [Cu-2([12]aneN(3)S)(2)(OH)(2)](ClO4)(2) were isolated, depending on the conditions employed. Molecular mechanics calculations have been employed to investigate the relative metal ion size preferences of the [3333], asym-[2424] and sym-[2424] conformation isomers. The calculations predict that the asym-[2424] conformer is most stable for M-N bond lengths in the range 2.00-2.25 Angstrom whilst for the larger metal ions the [3333] conformer is dominant. The disorder seen in the structure of the [Zn([12]aneN(3)S)(NO3)](+) complex is also explained by the calculations. (C) 1999 Elsevier Science Ltd. All rights reserved.

Reverse transcription of a naturally occurring nonretroviral RNA produces a precise deletion in the majority of its cDNA products

A precise, reproducible deletion made during in vitro reverse transcription of RNA2 from the icosahedral positive-stranded Helicoverpa armigera stunt virus (Tetraviridae) is described. The deletion, located between two hexamer repeats, is a 50-base sequence that includes one copy of the hexamer repeat. Only the Moloney murine leukemia virus reverse transcriptase and its derivative Superscript I, carrying a deletion of the carboxy-terminal RNase H region, showed this response, indicating a template-switching mechanism different from one proposed that involves a RNase H-dependent strand transfer, Superscript II, however, which carries point mutations to reduce RNase H activity, does not cause a deletion. A possible mechanism involves the enzyme pausing at the 3' side of a stem-loop structure and the 3' end of the nascent DNA strand separating from the template and reannealing to the upstream hexamer repeat.

Crystal structure of a heptameric Sm-like protein complex from archaea: Implications for the structure and evolution of snRNPs

The Sm/Lsm proteins associate with small nuclear RNA to form the core of small nuclear ribonucleoproteins, required for processes as diverse as pre-mRNA splicing, mRNA degradation and telomere formation. The Lsm proteins from archaea are likely to represent the ancestral Sm/Lsm domain. Here, we present the crystal structure of the Lsm alpha protein from the thermophilic archaeon Methanobacterium thermoautrophicum at 2.0 Angstrom resolution. The Lsm alpha protein crystallizes as a heptameric ring comprised of seven identical subunits interacting via beta -strand pairing and hydrophobic interactions. The heptamer can be viewed as a propeller-like structure in which each blade consists of a seven-stranded antiparallel beta -sheet formed from neighbouring subunits. There are seven slots on the inner surface of the heptamer ring, each of which is lined by Asp, Asn and Arg residues that are highly conserved in the Sm/Lsm sequences. These conserved slots are likely to form the RNA-binding site. In archaea, the gene encoding Lsm alpha is located next to the L37e ribosomal protein gene in a putative operon, suggesting a role for the Lsm alpha complex in ribosome function or biogenesis. (C) 2001 Academic Press.

Mechanism of ribonuclease inhibition by ribonuclease inhibitor protein based on the crystal structure of its complex with ribonuclease A

We describe the mechanism of ribonuclease inhibition by ribonuclease inhibitor, a protein built of leucine-rich repeats, based on the crystal structure of the complex between the inhibitor and ribonuclease A. The structure was determined by molecular replacement and refined to an R(cryst) of 19.4% at 2.5 Angstrom resolution. Ribonuclease A binds to the concave region of the inhibitor protein comprising its parallel beta-sheet and loops. The inhibitor covers the ribonuclease active site and directly contacts several active-site residues. The inhibitor only partially mimics the RNase-nucleotide interaction and does not utilize the pi phosphate-binding pocket of ribonuclease A, where a sulfate ion remains bound. The 2550 Angstrom(2) of accessible surface area buried upon complex formation may be one of the major contributors to the extremely tight association (K-i = 5.9 x 10(-14) M). The interaction is predominantly electrostatic; there is a high chemical complementarity with 18 putative hydrogen bonds and salt links, but the shape complementarity is lower than in most other protein-protein complexes. Ribonuclease inhibitor changes its conformation upon complex formation; the conformational change is unusual in that it is a plastic reorganization of the entire structure without any obvious hinge and reflects the conformational flexibility of the structure of the inhibitor. There is a good agreement between the crystal structure and other biochemical studies of the interaction. The structure suggests that the conformational flexibility of RI and an unusually large contact area that compensates for a lower degree of complementarity may be the principal reasons for the ability of RI to potently inhibit diverse ribonucleases. However, the inhibition is lost with amphibian ribonucleases that have substituted most residues corresponding to inhibitor-binding residues in RNase A, and with bovine seminal ribonuclease that prevents inhibitor binding by forming a dimer. (C) 1996 Academic Press Limited

Noncoding RNA's in mammals

Transcripts that lack any protein-coding potential represent at least half of the identified elements transcriptome. We review the evidence for the existence of such transcripts in the mammalian transcriptome, and argue that there may be many more noncoding RNAs (ncRNAs) still to be discovered. Relatively few ncRNA “genes” have been ascribed a function based upon mutation analysis. The review discusses possible roles of ncRNAs as cis-acting and trans-acting elements in epigenetic transcriptional control, including monoallelic gene silencing and imprinting. We also consider the evidence that the production of ncRNAs is a common feature of transcriptional enhancers.

Crocodile farms and management of saltwater crocodiles in Northern Territory: Results of a survey of NT crocodile farmers plus analysis of secondary information

After outlining some relevant background information about the NT crocodile farming industry and explaining the purpose of our survey of NT crocodile farmers conducted in the first half of 2005, this paper reports the results of the survey. The information received from the survey is supplemented by secondary data and by information from secondary sources. This report covers the location of respondents; the size of crocodile farms; farmers’ stated knowledge of and attitudes towards the NT Crocodile Management Plan; the involvement of farms in the harvesting of crocodile eggs and the costs involved; views of crocodile farmers about whether the NT Crocodile Management Plan encourages landholders to conserve crocodiles and their perceptions of the benefits to landholders; predicted production trends and trends in the number of farms operating in NT; economic characteristics of crocodile farms producing in NT including the economic advantages and disadvantages of crocodile farming in NT. Concluding comments provide, amongst other things, an overview of the structure of the crocodile farming industry in the NT gleaned from a consideration of data available from the NT Government’s Department of Business, Industry and Resource Development.