Integrating pharmacogenetics into society: in search of a model

There has been considerable scientific, corporate and policy interest in the more effective use of genetics in both drug development and delivery. Pharmacogenetics - the study of the relationship between an individual's genetic makeup and response to medicinal drugs - has attracted global interest, but will it live up to its promise? Looking beyond the hype that has accompanied much of the commentary in the area, the future of pharmacogenetics will depend on how competing interests and options are resolved.     

Quantitative analysis of mouse urine volatiles: in search of MHC-dependent differences

Genes of the major histocompatibility complex (MHC), which play a critical role in immune recognition, influence mating preference and other social behaviors in mice. Training experiments using urine scent from mice differing only in the MHC complex, from MHC class I mutants or from knock-out mice lacking functional MHC class I molecules (beta2m-deficient), suggest that these behavioral effects are mediated by differences in MHC-dependent volatile components. In search for the physical basis of these behavioral studies, we have conducted a comparison of urinary volatiles in three sub-strains of C57BL/6 mice, a beta2m-deficient mutant lacking functional MHC class I expression and two unrelated inbred strains, using the technique of sorptive extraction with polydimethylsiloxan and subsequent analysis by gas chromatography/mass spectrometry. We show (i) that qualitative differences occur between different inbred strains but not in mice with the C57BL/6 background, (ii) that the individual variability in abundance in the same mouse strain is strongly component-dependent, (iii) that C57BL/6 sub-strains obtained from different provenance show a higher fraction of quantitative differences than a sub-strain and its beta2m-mutant obtained from the same source and (iv) that comparison of the spectra of beta2m mice and the corresponding wild type reveals no qualitative differences in close to 200 major and minor components and only minimal differences in a few substances from an ensemble of 69 selected for quantitative analysis. Our data suggest that odor is shaped by ontogenetic, environmental and genetic factors, and the gestalt of this scent may identify a mouse on the individual and population level; but, within the limits of the ensemble of components analysed, the results do not support the notion that functional MHC class I molecules influence the urinary volatile composition.     

HIV/AIDS: in search of an animal model

AIDS is among the most devastating diseases of our time, claiming the lives of approximately 3 million people per year. The primary cause of AIDS, human immunodeficiency virus type 1 (HIV-1), is a pathogen that is highly specific for humans and generally does not infect or cause disease in other species. This property complicates the generation of animal models that are urgently needed to test new antiretroviral therapies and vaccines. The most practical animal models developed to date consist of infection of rhesus macaques with a simian immunodeficiency virus (SIV) or chimeric HIV/SIV viruses. Although these models are useful for particular applications, the fact that SIV is a distinct virus compared with HIV-1 represents a significant limitation to their use. Here, we discuss the uses and limitations of existing models and recent advances that might lead to better animal models for HIV/AIDS.     

Stochastic search gene suggestion: a Bayesian hierarchical model for gene mapping

Mapping the genes for a complex disease, such as diabetes or rheumatoid arthritis (RA), involves finding multiple genetic loci that may contribute to the onset of the disease. Pairwise testing of the loci leads to the problem of multiple testing. Looking at haplotypes, or linear sets of loci, avoids multiple tests but results in a contingency table with sparse counts, especially when using marker loci with multiple alleles. We propose a hierarchical Bayesian model for case-parent triad data that uses a conditional logistic regression likelihood to model the probability of transmission to a diseased child. We define hierarchical prior distributions on the allele main effects to model the genetic dependencies present in the human leukocyte antigen (HLA) region of chromosome 6. First, we add a hierarchical level for model selection that accounts for both locus and allele selection. This allows us to cast the problem of identifying genetic loci relevant to the disease into a problem of Bayesian variable selection. Second, we attempt to include linkage disequilibrium as a covariance structure in the prior for model coefficients. We evaluate the performance of the procedure with some simulated examples and then apply our procedure to identifying genetic markers in the HLA region that influence risk for RA. Our software is available on the website http://www.epigenetic.org/Linkage/ssgs-public/.     

ESS++: a C++ objected-oriented algorithm for Bayesian stochastic search model exploration

SUMMARY: ESS++ is a C++ implementation of a fully Bayesian variable selection approach for single and multiple response linear regression. ESS++ works well both when the number of observations is larger than the number of predictors and in the 'large p, small n' case. In the current version, ESS++ can handle several hundred observations, thousands of predictors and a few responses simultaneously. The core engine of ESS++ for the selection of relevant predictors is based on Evolutionary Monte Carlo. Our implementation is open source, allowing community-based alterations and improvements. AVAILABILITY: C++ source code and documentation including compilation instructions are available under GNU licence at http://bgx.org.uk/software/ESS.html.     

Concentration compared with total urinary excretion of 11,17-DOA in cynomolgus monkey urine

Strees sensitive molecules exhibit great variation in concentration in the circulation and it may often be advantageous to quantify these in urine or feces rather than in serum or plasma. We advocate that all urine-or feces-should be collected, and that excretion of stress sensitive molecules should be expressed as amounts excreted per time unit per kg body-weight, rather than being expressed as concentrations in samples. Urine and feces excretion varies significantly within and between animals over time, which may render simple concentration measures of molecules of little biological relevance.     

Concentration effects of volatile anesthetics on the properties of model membranes: a coarse-grain approach

To gain insights into the molecular level mechanism of drug action at the membrane site, we have carried out extensive molecular dynamics simulations of a model membrane in the presence of a volatile anesthetic using a coarse-grain model. Six different anesthetic (halothane)/lipid (dimyristoylphosphatidylcholine) ratios have been investigated, going beyond the low doses typical of medical applications. The volatile anesthetics were introduced into a preassembled fully hydrated 512-molecule lipid bilayer and each of the molecular dynamics simulations were carried out at ambient conditions, using the NPT ensemble. The area per lipid increases monotonically with the halothane concentration and the lamellar spacing decreases, whereas the lipid bilayer thickness shows no appreciable differences and only a slight increase upon addition of halothane. The density profiles of the anesthetic molecules display a bimodal distribution along the membrane normal with maxima located close to the lipid-water interface region. We have studied how halothane molecules fluctuate between the two maxima of the bimodal distribution and we observed a different mechanism at low and high anesthetic concentrations. Through the investigation of the reorientational motions of the lipid tails, we found that the anesthetic molecules increase the segmental order of the lipids close to the membrane surface.     

Fast tree search for a triangular lattice model of protein folding

Using a triangular lattice model to study the designability of protein folding, we overcame the parity problem of previous cubic lattice model and enumerated all the sequences and compact structures on a simple two-dimensional triangular lattice model of size 4 5 6 5 4. We used two types of amino acids, hydrophobic and polar, to make up the sequences, and achieved 223W212 different sequences excluding the reverse symmetry sequences. The total string number of distinct compact structures was 219,093, excluding reflection symmetry in the self-avoiding path of length 24 triangular lattice model. Based on this model, we applied a fast search algorithm by constructing a cluster tree. The algorithm decreased the computation by computing the objective energy of non-leaf nodes. The parallel experiments proved that the fast tree search algorithm yielded an exponential speed-up in the model of size 4 5 6 5 4. Designability analysis was performed to understand the search result.     

HaMStR: Profile hidden markov model based search for orthologs in ESTs

Background  

EST sequencing is a versatile approach for rapidly gathering protein coding sequences. They provide direct access to an organism's gene repertoire bypassing the still error-prone procedure of gene prediction from genomic data. Therefore, ESTs are often the only source for biological sequence data from taxa outside mainstream interest. The widespread use of ESTs in evolutionary studies and particularly in molecular systematics studies is still hindered by the lack of efficient and reliable approaches for automated ortholog predictions in ESTs. Existing methods either depend on a known species tree or cannot cope with redundancy in EST data.     

The jigsaw puzzle model: search for conformational specificity in protein interiors

The jigsaw puzzle model postulates that the predominant factor relating primary sequence to three-dimensional fold lies in the stereospecific packing of interdigitating side-chains within densely packed protein interiors. An attempt has been made to check the validity of the model by means of a surface complementarity function. Out of a database of 100 highly resolved protein structures the contacts between buried hydrophobic residues (Leu, Ile, Val, Phe) and their neighbours have been categorized in terms of the extent of side-chain surface area involved in a contact (overlap) and their steric fit (Sm). The results show that the majority of contacts between a buried residue and its immediate neighbours (side-chains) are of high steric fit and in the case of extended overlap at least one of the angular parameters characterizing interresidue geometry to have pronounced deviation from a random distribution, estimated by chi(2). The calculations thus tend to support the "jigsaw puzzle" model in that 75-85% of the contacts involving hydrophobic residues are of high surface complementarity, which, coupled to high overlap, exercise fairly stringent constraints over the possible geometrical orientations between interacting residues. These constraints manifest in simple patterns in the distributions of orientational angles. Approximately 60-80% of the buried side-chain surface packs against neighbouring side-chains, the rest interacting with main-chain atoms. The latter partition of the surface maintains an equally high steric fit (relative to side-chain contacts) emphasizing a non-trivial though secondary role played by main-chain atoms in interior packing. The majority of this class of contacts, though of high complementarity, is of reduced overlap. All residues whether hydrophobic or polar/charged show similar surface complementarity measures upon burial, indicating comparable competence of all amino acids in packing effectively with their atomic environments. The specificity thus appears to be distributed over the entire network of contacts within proteins. The study concludes with a proposal to classify contacts as specific and non-specific (based on overlap and fit), with the former perhaps contributing more to the specificity between sequence and fold than the latter.