Screening for functional expression and overexpression of a family of diiron-containing interfacial membrane proteins using the univector recombination system

The large number of uncharacterized genes emerging from genome sequencing projects has resulted in a need for quick and reliable screening methods for protein expression parameters. We have utilized the univector plasmid recombination system (as previously reported) to develop a series of vectors for rapid screening for expression in Escherichia coli. A high level of recombinant protein expression is a requirement for purification of protein for structural determination and other purposes. In other applications, successful complementation of a missing enzyme activity in E. coli, as well as directed evolution studies and metabolic engineering, often require a much lower level of protein expression. In this report we describe the construction of a number of new pHOST vectors that can be screened for both low- and high-level expression. We isolated a mutant vector for MBP fusions that exhibited a more optimal level of expression for complementation of aerobic respiration in hemA(-) E. coli, our functional assay for the alternative oxidase. We then demonstrated the use of our system to rapidly screen for both optimal functional expression and optimal overexpression of the alternative oxidase as well as two other members of a family of membrane-bound diiron carboxylate proteins, the plastid terminal oxidase and 5-demethoxyquinone hydroxylase.     

Estimation and confidence regions for multi-dimensional effective dose

The problem of finding confidence regions for multiple predictor variables corresponding to given expected values of a response variable has not been adequately resolved. Motivated by an example from a study on hyperbaric exposure using a logistic regression model, we develop a conceptual framework for the estimation of the multi-dimensional effective dose for binary outcomes. The k -dimensional effective dose can be determined by conditioning on k - 1 components and solving for the last component as a conditional univariate effective dose. We consider various approaches for calculating confidence regions for the multi-dimensional effective dose and compare them via a simulation study for a range of possible designs. We analyze data related to decompression sickness to illustrate our procedure. Our results provide a practical approach to finding confidence regions for predictor variables for a given response value.     

Mass Ligand Binding Screening for Receptor Antagonists: Prototype New Drugs and Blind Alleys

Abstract

The ligand binding assay is a powerful tool in the search for antagonists for novel receptors, and for identification of novel classes of antagonists for well-known receptors. Ligand binding mass screening can be adapted for very high throughput. In order for mass screening to be useful, it is necessary to strictly define the binding characteristics for a compound to be considered a putative receptor antagonist. In practice, we have found that synthetic pursuit of a compound with a K_i of ± 1 uM is likely to lead down a blind alley unless very good evidence for specificity is available. Even potent competitors for binding should be thoroughly evaluated in assays of biological activity before a synthetic program is initiated in earnest.     

Selecting pedigrees for linkage analysis of a quantitative trait: the expected number of informative meioses.

With evidence of segregation at a major locus for a quantitative trait having been found, a logical next step is to select a subset of the pedigrees to include in a linkage study to map the major locus. Ideally this subset should include much of the linkage information in the sample but include only a fraction of the pedigrees. We previously described a strategy for selecting pedigrees for linkage analysis of a quantitative trait on the basis of a pedigree likelihood-ratio statistic. For quantitative traits controlled by a major locus with a rare dominant allele, the likelihood-ratio strategy extracted nearly all the information for linkage while typically requiring marker data on only about one-third of the pedigrees. Here, we describe a new strategy to select pedigrees for linkage analysis on the basis of the expected number of potentially informative meioses in each pedigree. We demonstrate that this informative-meioses strategy provides an efficient and more general means to select pedigrees for a linkage study of a quantitative trait.     

Hexanoylation of a VPAC2 receptor-preferring ligand markedly increased its selectivity and potency

We synthesized a VIP analog that combines mutations that decrease the affinity for the VPAC1 receptor but maintain a high affinity for the VPAC2 receptor with an amino-terminal hexanoylation that increases the affinity for the VPAC2 receptor with a limited decrease in the affinity of the VPAC1 receptor. The resulting Hexanoyl[A19,K(27,28)]VIP had the expected properties of a high affinity for the VPAC2 receptor and a low affinity for the VPAC1 receptor and also a low affinity for the PAC1 and secretin receptors. With a 1000-fold preference for the VPAC2 receptor and a IC50 value of binding of 1 nM, this compound is the most potent and the most selective agonist presently described.     

The role of phadiatop in the screening of respiratory allergy

Phadiatop is a new in vitro screening test for respiratory allergy. This test, based on the RAST procedure, detects in serum, the presence of specific IgE to a mixture of common inhalent allergens. Among 70 patients (26 children and 44 adults) consulting for respiratory syndrome, Phadiatop was positive in 31 cases. There were a good correlation between this new test and skin tests (59% for adults and 92% for children), total IgE (70% for adults and 65% for children) and RAST (93% for adults and 96% for children). Phadiatop, with a specificity of 100%, a sensitivity of 82% (76% for adults and 92% for children) and an efficiency of 90% (86% for adults and 96% for children), is a more accurate test than total IgE and could be an excellent in vitro screening test for respiratory allergy.     

How large a managed area is needed to protect a threatened animal species? Combining simple dispersal and population models

The size and shape of the managed area for a threatened species to have a stable or growing population is a central issue for conservation management, for example for kiwi, Apteryx sp. Combining geometric probability results for retention of juveniles to breed in the protected area with a standard matrix population model allows the creation of an explicit relationship between the minimum area needed and how far juveniles of the species disperse to establish breeding territories. For a given set of demographic parameters for the population, and a rectangular protected area, there is a quadratic relationship between area and dispersal distance. Extensions for a circular protected area, and for a probability distribution on dispersal distance, are considered. The results are applied to kiwi, using established population parameters, giving results that match closely those from a previously published simulation. This approach can provide simple tools, readily implemented in a spreadsheet, for assessment of the size and shape of protected areas needed in conservation management of animals that disperse before breeding.     

Quantitative analysis of the ion-dependent folding stability of DNA triplexes

A DNA triplex is formed through binding of a third strand to the major groove of a duplex. Due to the high charge density of a DNA triplex, metal ions are critical for its stability. We recently developed the tightly bound ion (TBI) model for ion-nucleic acids interactions. The model accounts for the potential correlation and fluctuations of the ion distribution. We now apply the TBI model to analyze the ion dependence of the thermodynamic stability for DNA triplexes. We focus on two experimentally studied systems: a 24-base DNA triplex and a pair of interacting 14-base triplexes. Our theoretical calculations for the number of bound ions indicate that the TBI model provides improved predictions for the number of bound ions than the classical Poisson-Boltzmann (PB) equation. The improvement is more significant for a triplex, which has a higher charge density than a duplex. This is possibly due to the higher ion concentration around the triplex and hence a stronger ion correlation effect for a triplex. In addition, our analysis for the free energy landscape for a pair of 14-mer triplexes immersed in an ionic solution shows that divalent ions could induce an attractive force between the triplexes. Furthermore, we investigate how the protonated cytosines in the triplexes affect the stability of the triplex helices.     

Integrated allele-specific polymerase chain reaction-capillary electrophoresis microdevice for single nucleotide polymorphism genotyping

An integrated allele-specific (AS) polymerase chain reaction (PCR) and capillary electrophoresis (CE) microdevice has been developed for multiplex single nucleotide polymorphism (SNP) genotyping on a portable instrumentation, which was applied for on-site identification of HANWOO (Korean indigenous beef cattle). Twelve sets of primers were designed for targeting beef cattle's eleven SNP loci for HANWOO verification and one primer set for a positive PCR control, and the success rate for identification of HANWOO was demonstrated statistically. The AS PCR and CE separation for multiplex SNP typing was carried out on a glass-based microchip consisting of four layers: a microchannel plate for microfluidic control, a Pt-electrode plate for a resistance temperature detector (RTD), a poly(dimethylsiloxane) (PDMS) membrane and a manifold glass for microvalve function. The operation of the sample loading, AS PCR, microvalve, and CE on a chip was automated with a portable genetic analyzer, and the laser-induced fluorescence detection was performed on a miniaturized fluorescence detector. The blind samples were correctly identified as a HANWOO by showing one or two amplicon peaks in the electropherogram, while the imported beef cattle revealed more than five peaks. Our genetic analysis platform provides rapid, accurate, and on-site multiplex SNP typing.     

Constructing and counting phylogenetic invariants.

The method of invariants is an approach to the problem of reconstructing the phylogenetic tree of a collection of m taxa using nucleotide sequence data. Models for the respective probabilities of the 4m possible vectors of bases at a given site will have unknown parameters that describe the random mechanism by which substitution occurs along the branches of a putative phylogenetic tree. An invariant is a polynomial in these probabilities that, for a given phylogeny, is zero for all choices of the substitution mechanism parameters. If the invariant is typically non-zero for another phylogenetic tree, then estimates of the invariant can be used as evidence to support one phylogeny over another. Previous work of Evans and Speed showed that, for certain commonly used substitution models, the problem of finding a minimal generating set for the ideal of invariants can be reduced to the linear algebra problem of finding a basis for a certain lattice (that is, a free Z-module). They also conjectured that the cardinality of such a generating set can be computed using a simple "degrees of freedom" formula. We verify this conjecture. Along the way, we explain in detail how the observations of Evans and Speed lead to a simple, computationally feasible algorithm for constructing a minimal generating set.     

Tools for comparative protein structure modeling and analysis

The following resources for comparative protein structure modeling and analysis are described (http://salilab.org): MODELLER, a program for comparative modeling by satisfaction of spatial restraints; MODWEB, a web server for automated comparative modeling that relies on PSI-BLAST, IMPALA and MODELLER; MODLOOP, a web server for automated loop modeling that relies on MODELLER; MOULDER, a CPU intensive protocol of MODWEB for building comparative models based on distant known structures; MODBASE, a comprehensive database of annotated comparative models for all sequences detectably related to a known structure; MODVIEW, a Netscape plugin for Linux that integrates viewing of multiple sequences and structures; and SNPWEB, a web server for structure-based prediction of the functional impact of a single amino acid substitution.     

Timing, memory for intervals, and memory for untimed stimuli: The role of instructional ambiguity

Theories of animal timing have had to account for findings that the memory for the duration of a timed interval appears to be dramatically shorted within a short time of its termination. This finding has led to the subjective shortening hypothesis and it has been proposed to account for the poor memory that animals appear to have for the initial portion of a timed interval when a gap is inserted in the to-be-timed signal. It has also been proposed to account for the poor memory for a relatively long interval that has been discriminated from a shorter interval. I suggest here a simpler account in which ambiguity between the gap or retention interval and the intertrial interval results in resetting the clock, rather than forgetting the interval. The ambiguity hypothesis, together with a signal salience mechanism that determines how quickly the clock is reset at the start of the intertrial interval can account for the results of the reported timing experiments that have used the peak procedure. Furthermore, instructional ambiguity rather than memory loss may account for the results of many animal memory experiments that do not involve memory for time.     

Analysis of the free energy in a stochastic RNA secondary structure model

There are two custom ways for predicting RNA secondary structures: minimizing the free energy of a conformation according to a thermodynamic model and maximizing the probability of a folding according to a stochastic model. In most cases, stochastic grammars are used for the latter alternative applying the maximum likelihood principle for determining a grammar's probabilities. In this paper, building on such a stochastic model, we will analyze the expected minimum free energy of an RNA molecule according to Turner's energy rules. Even if the parameters of our grammar are chosen with respect to structural properties of native molecules only (and therefore, independent of molecules' free energy), we prove formulae for the expected minimum free energy and the corresponding variance as functions of the molecule's size which perfectly fit the native behavior of free energies. This gives proof for a high quality of our stochastic model making it a handy tool for further investigations. In fact, the stochastic model for RNA secondary structures presented in this work has, for example, been used as the basis of a new algorithm for the (nonuniform) generation of random RNA secondary structures.     

POPULATION GENETICS OF INTRAGAMETOPHYTIC SELFING

Intragametophytic selfing is a mode of reproduction occurring in homosporous ferns where two gametes from the same haploid gametophyte form a completely homozygous sporophyte. The inbreeding equilibrium is derived for a population with partial intragametophytic selfing, selfing, and outcrossing. Procedures for directly estimating the extent of intragametophytic selfing and selfing using parent-offspring data are given. The conditions for a stable polymorphism from a heterozygous-advantage fitness model are more restrictive for partial intragametophytic selfing than for selfing. The rate of decay of gametic disequilibrium is slower for partial intragametophytic selfing than for selfing. Based on these findings, one would predict that plants with intragametophytic selfing would have less polymorphism for loci with a heterozygous advantage and more gametic disequilibrium between neutral loci than is expected for populations with an equivalent amount of selfing. Data from several studies are consistent with these predictions.     

Molecular modeling and dynamics simulations of PNP from Streptococcus agalactiae

This work describes for the first time a structural model of purine nucleoside phosphorylase from Streptococcus agalactiae (SaPNP). PNP catalyzes the cleavage of N-ribosidic bonds of the purine ribonucleosides and 2-deoxyribonucleosides in the presence of inorganic orthophosphate as a second substrate. This enzyme is a potential target for the development of antibacterial drugs. We modeled the complexes of SaPNP with 15 different ligands in order to determine the structural basis for the specificity of these ligands against SaPNP. The application of a novel empirical scoring function to estimate the affinity of a ligand for a protein was able to identify the ligands with high affinity for PNPs. The analysis of molecular dynamics trajectory for SaPNP indicates that the functionally important motifs have a very stable structure. This new structural model together with a novel empirical scoring function opens the possibility to explorer larger library of compounds in order to identify the new inhibitors for PNPs in virtual screening projects.     

High-affinity binding of Ca2+ to bovine alpha-lactalbumin in the absence and presence of EGTA.

Literature values for the Kd for Ca2+ in bovine alpha-lactalbumin range over 3 orders of magnitude. There is a difference between two results obtained with EGTA as a metal-ion buffer, partly because different values for the Kd of Ca2+-EGTA were used in the calculations, and a much wider difference between results obtained in the presence and absence of EGTA, which has been attributed to an interaction between EGTA and the protein. Titrations in a flow-dialysis cell showed that Mn2+ competed with Ca2+ for the high-affinity site on the protein, and the results, combined with a Kd for Mn2+ of 2.1 +/- 0.1 microM, which was determined fluorimetrically, gave a Kd for Ca2+ of 1.3 +/- 0.1 nM. When alpha-lactalbumin containing 45Ca2+ was titrated with EGTA in a flow-dialysis cell, and widely accepted metal-chelation data for EGTA were used in the calculations, a Kd for Ca2+ of 1.10 +/- 0.03 nM was obtained. The results from the two methods are so similar as to indicate that the affinity for Ca2+ was unaffected by the presence of EGTA.     

The search for hematologic-immunological parameters with prognostic value in children with various tumors]

In a total of 35 children affected with malignant sufferings (except leucaemias) the haematological-immunological parameters of children surviving for a short time are compared with those surviving for a long time. Compared with children surviving for a short time those surviving for long time had a more "normal" finding initially and during the control examination made after more than 5 years of a survival time which was free of complaints and recidives. The immunologically intact condition of the organism is assumed to be decisive for this course of disease.     

Affinity of corticosteroids for mineralocorticoid and glucocorticoid receptors of the rabbit kidney: effect of steroid substitution

Corticosteroid derivatives coupled in the C3, C7 or C17 position with a long aliphatic chain were synthesized in order to select a suitable ligand for the preparation of a biospecific affinity adsorbent for mineralocorticoid receptor purification. The affinity of these derivatives for mineralocorticoid receptors (MR) and glucocorticoid receptors (GR) was explored in rabbit kidney cytosol. In this model, aldosterone bound to a single class of receptors with high affinity (Kd 1 nM) and mineralocorticoid specificity. RU26988, a highly specific ligand for GR, did not compete for these sites. The C7 and C17 positions were found to be of crucial importance in the steroid's interaction with the mineralocorticoid receptors, since the linkage of a long side chain in these positions induced complete loss of affinity. Hence, deoxycorticosterone no longer bound to MR after 17 beta substitution with a 9-carbon aliphatic chain. This loss of affinity was not observed for glucocorticoids. The 17 beta nonylamide derivative of dexamethasone still competed for GR. Increasing the length of the C7 side of the spirolactone SC26304 suppressed its affinity for MR. Finally, C3 was an appropriate position for steroid substitution. The 3-nonylamide of carboxymethyloxime deoxycorticosterone bound to MR but not to GR, and therefore constitutes a suitable ligand for the preparation of a mineralocorticoid adsorbent.