Bayesian latent variable models for median regression on multiple outcomes

Often a response of interest cannot be measured directly and it is necessary to rely on multiple surrogates, which can be assumed to be conditionally independent given the latent response and observed covariates. Latent response models typically assume that residual densities are Gaussian. This article proposes a Bayesian median regression modeling approach, which avoids parametric assumptions about residual densities by relying on an approximation based on quantiles. To accommodate within-subject dependency, the quantile response categories of the surrogate outcomes are related to underlying normal variables, which depend on a latent normal response. This underlying Gaussian covariance structure simplifies interpretation and model fitting, without restricting the marginal densities of the surrogate outcomes. A Markov chain Monte Carlo algorithm is proposed for posterior computation, and the methods are applied to single-cell electrophoresis (comet assay) data from a genetic toxicology study.     

Hydrogen-bonding propensities of sphingomyelin in solution and in a bilayer assembly: a molecular dynamics study

Sphingomyelin is enriched within lipid microdomains of the cell membrane termed lipid rafts. These microdomains play a part in regulating a variety of cellular events. Computer simulations of the hydrogen-bonding properties of sphingolipids, believed to be central to the organization of these domains, can delineate the possible molecular interactions that underlie this lipid structure. We have therefore used molecular dynamics simulations to unravel the hydrogen-bonding behavior of palmitoylsphingomyelin (PSM). A series of eight simulations of 3 ns each of a single PSM molecule in water showed that the sphingosine OH and NH groups can form hydrogen bonds with the phosphate oxygens of their own polar head, in agreement with NMR data. Simulations of PSM in a bilayer assembly were carried out for 8 ns with three different force field parameterizations. The major physico-chemical parameters of the simulated bilayer agree with those established experimentally. The sphingosine OH group was mainly involved in intramolecular hydrogen bonds, in contrast to the almost exclusive intermolecular hydrogen bonds formed by the amide NH moiety. During the bilayer simulations the intermolecular hydrogen bonds among lipids formed a dynamic network characterized by the presence of hydrogen-bonded lipid clusters of up to nine PSM molecules.     

Asymmetric structural motions of the homomeric alpha7 nicotinic receptor ligand binding domain revealed by molecular dynamics simulation

A homology model of the ligand binding domain of the alpha7 nicotinic receptor is constructed based on the acetylcholine-binding protein crystal structure. This structure is refined in a 10 ns molecular dynamics simulation. The modeled structure proves fairly resilient, with no significant changes at the secondary or tertiary structural levels. The hypothesis that the acetylcholine-binding protein template is in the activated or desensitized state, and the absence of a bound agonist in the simulation suggests that the structure may also be relaxing from this state to the activatable state. Candidate motions that take place involve not only the side chains of residues lining the binding sites, but also the subunit positions that determine the overall shape of the receptor. In particular, two nonadjacent subunits move outward, whereas their partners counterclockwise to them move inward, leading to a marginally wider interface between themselves and an overall asymmetric structure. This in turn affects the binding sites, producing two that are more open and characterized by distinct side-chain conformations of W54 and L118, although motions of the side chains of all residues in every binding site still contribute to a reduction in binding site size, especially the outward motion of W148, which hinders acetylcholine binding. The Cys loop at the membrane interface also displays some flexibility. Although the short simulation timescale is unlikely to sample adequately all the conformational states, the pattern of observed motions suggests how ligand binding may correlate with larger-scale subunit motions that would connect with the transmembrane region that controls the passage of ions. Furthermore, the shape of the asymmetry with binding sites of differing affinity for acetylcholine, characteristic of other nicotinic receptors, may be a natural property of the relaxed, activatable state of alpha7.     

Structure-based design of novel nonpeptide inhibitors of the Src SH2 domain:phosphotyrosine mimetics exploiting multifunctional group replacement chemistry

A series of novel nonpeptide inhibitors of the pp60(c-Src) (Src) SH2 domain is described that exploit multifunctional group replacement of the phenylphosphate moiety of phosphotyrosine (pTyr). Relative to an x-ray structure of citrate complexed to the pTyr binding site of the Src SH2 domain, these nonpeptide ligands illustrate the systematic replacement of the phosphate group by multiple nonhydrolyzable, mono- or dianionic functionalities. Specifically, several phenylalanine (Phe) analogs incorporating key 4' and 3' substituents were synthesized and incorporated into a bicyclic benzamide template previously reported (W. C. Shakespeare et al., Proceedings of the National Academy of Science USA, 2000, Vol. 97, pp. 9373-9378). These pTyr mimetics included 4',3'-diphosphono-Phe (Dpp), 4',3'-dicarboxymethyloxy-Phe (Dcp), and 4'-phosphono-3'-carboxymethyloxy-Phe (Cpp). Noteworthy were nonpeptide inhibitors 8-11 that were 5- to 10-fold more potent than the cognate tetrapeptide ligand Ac-pTyr-Glu-Glu-Ile-NH(2) in binding to the Src SH2 domain.     

M142.2 (cut-6), a novel Caenorhabditis elegans matrix gene important for dauer body shape

The cuticle of the nematode Caenorhabditis elegans is a collagenous extracellular matrix which forms the exoskeleton and defines the shape of the worm. We have characterized the C. elegans gene M142.2, and we show that this is a developmentally regulated gene important for cuticle structure. Transgenic worms expressing M142.2 promoter fused to green fluorescent protein showed that M142.2 is expressed in late embryos and L2d predauers, in the hypodermal cells which synthesize the cuticle. The same temporal pattern was seen by RT-PCR using RNA purified from specific developmental stages. A recombinant fragment of M142.2 was expressed in Escherichia coli and used to raise an antiserum. Immunohistochemistry using the antiserum localized M142.2 to the periphery of the alae of L1 and dauers, forming two longitudinal ribbons over the hypodermal cells. Loss-of-function of M142.2 by RNAi resulted in a novel phenotype: dumpy dauers which lacked alae. M142.2 therefore plays a major role in the assembly of the alae and the morphology of the dauer cuticle; because of its similarity to the other cut genes of the cuticle, we have named the gene cut-6.     

Neurovascular congruence results from a shared patterning mechanism that utilizes Semaphorin3A and Neuropilin-1

Peripheral nerves and blood vessels have similar patterns in quail forelimb development. Usually, nerves extend adjacent to existing blood vessels, but in a few cases, vessels follow nerves. Nerves have been proposed to follow vascular smooth muscle, endothelium, or their basal laminae. Focusing on the major axial blood vessels and nerves, we found that when nerves grow into forelimbs at E3.5-E5, vascular smooth muscle was not detectable by smooth muscle actin immunoreactivity. Additionally, transmission electron microscopy at E5.5 confirmed that early blood vessels lacked smooth muscle and showed that the endothelial cell layer lacks a basal lamina, and we did not observe physical contact between peripheral nerves and these endothelial cells. To test more generally whether lack of nerves affected blood vessel patterns, forelimb-level neural tube ablations were performed at E2 to produce aneural limbs; these had completely normal vascular patterns up to at least E10. To test more generally whether vascular perturbation affected nerve patterns, VEGF(165), VEGF(121), Ang-1, and soluble Flt-1/Fc proteins singly and in combination were focally introduced via beads implanted into E4.5 forelimbs. These produced significant alterations to the vascular patterns, which included the formation of neo-vessels and the creation of ectopic avascular spaces at E6, but in both under- and overvascularized forelimbs, the peripheral nerve pattern was normal. The spatial distribution of semaphorin3A protein immunoreactivity was consistent with a negative regulation of neural and/or vascular patterning. Semaphorin3A bead implantations into E4.5 forelimbs caused failure of nerves and blood vessels to form and to deviate away from the bead. Conversely, semaphorin3A antibody bead implantation was associated with a local increase in capillary formation. Furthermore, neural tube electroporation at E2 with a construct for the soluble form of neuropilin-1 caused vascular malformations and hemorrhage as well as altered nerve trajectories and peripheral nerve defasciculation at E5-E6. These results suggest that neurovascular congruency does not arise from interdependence between peripheral nerves and blood vessels, but supports the hypothesis that it arises by a shared patterning mechanism that utilizes semaphorin3A.     

A multilocus genealogical approach to phylogenetic species recognition in the model eukaryote Neurospora

To critically examine the relationship between species recognized by phylogenetic and reproductive compatibility criteria, we applied phylogenetic species recognition (PSR) to the fungus in which biological species recognition (BSR) has been most comprehensively applied, the well-studied genus Neurospora. Four independent anonymous nuclear loci were characterized and sequenced from 147 individuals that were representative of all described outbreeding species of Neurospora. We developed a consensus-tree approach that identified monophyletic genealogical groups that were concordantly supported by the majority of the loci, or were well supported by at least one locus but not contradicted by any other locus. We recognized a total of eight phylogenetic species, five of which corresponded with the five traditional biological species, and three of which were newly discovered. Not only were phylogenetic criteria superior to traditional reproductive compatibility criteria in revealing the full species diversity of Neurospora, but also significant phylogenetic subdivisions were detected within some species. Despite previous suggestions of hybridization between N. crassa and N. intermedia in nature, and the fact that several putative hybrid individuals were included in this study, no molecular evidence in support of recent interspecific gene flow or the existence of true hybrids was observed. The sequence data from the four loci were combined and used to clarify how the species discovered by PSR were related. Although species-level clades were strongly supported, the phylogenetic relationships among species remained difficult to resolve, perhaps due to conflicting signals resulting from differential lineage sorting.     

The influence of history and contemporary stream hydrology on the evolution of genetic diversity within species: an examination of microsatellite DNA variation in bull trout,Salvelinus confluentus (Pisces: Salmonidae)

An understanding of the relative roles of historical and contemporary factors in structuring genetic variation is a fundamental, but understudied aspect of geographic variation. We examined geographic variation in microsatellite DNA allele frequencies in bull trout (Salvelinus confluentus, Salmonidae) to test hypotheses concerning the relative roles of postglacial dispersal (historical) and current landscape features (contemporary) in structuring genetic variability and population differentiation. Bull trout exhibit relatively low intrapopulation microsatellite variation (average of 1.9 alleles per locus, average He = 0.24), but high levels of interpopulation divergence (F(ST) = 0.39). We found evidence of historical influences on microsatellite variation in the form of a decrease in the number of alleles and heterozygosities in populations on the periphery of the range relative to populations closer to putative glacial refugia. In addition, one region of British Columbia that was colonized later during deglaciation and by more indirect watershed connections showed less developed and more variable patterns of isolation by distance than a similar region colonized earlier and more directly from refugia. Current spatial and drainage interconnectedness among sites and the presence of migration barriers (falls and cascades) within individual streams were found to be important contemporary factors influencing historical patterns of genetic variability and interpopulation divergence. Our work illustrates the limited utility of equilibrium models to delineate population structure and patterns of genetic diversity in recently founded populations or those inhabiting highly heterogeneous environments, and it highlights the need for approaches incorporating a landscape context for population divergence. Substantial microsatellite DNA divergence among bull trout populations may also signal divergence in traits important to population persistence in specific environments.     

Control of Ras cycling by Ca2+

Ras GTPases are binary switches, cycling between an inactive GDP-bound form and an active GTP-bound form at the membrane. They transduce signals into the cytoplasm via effector pathways that regulate cell growth, differentiation and apoptosis. Ras activation is enhanced by guanine nucleotide exchange factors (GEFs); deactivation is accelerated by GTPase-activating proteins (GAPs). Recently, new roles for Ca(2+) and diacylglycerol (DAG) in the control of Ras cycling have emerged with the discovery of a series of novel GEFs and GAPs. These regulators of Ras cycling are likely to play a key role in the information processing of Ca(2+) and DAG signals.     

Checkpoint arrest signaling in response to UV damage is independent of nucleotide excision repair in Saccharomyces cerevisiae

The recognition of DNA double-stranded breaks or single-stranded DNA gaps as a precondition for cell cycle checkpoint arrest has been well established. However, how bulky base damage such as UV-induced pyrimidine dimers elicits a checkpoint response has remained elusive. Nucleotide excision repair represents the main pathway for UV dimer removal that results in strand interruptions. However, we demonstrate here that Rad53p hyperphosphorylation, an early event of checkpoint signaling in Saccharomyces cerevisiae, is independent of nucleotide excision repair (NER), even if replication as a source of secondary DNA damage is excluded. Thus, our data hint at primary base damage or at UV damage (primary or secondary) that does not need to be processed by NER as the relevant substrate of damage-sensing checkpoint proteins.