Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

High incidence of proviral integrations in the Hoxa locus in a new model of E2a-PBX1-induced B-cell leukemia

Relevant mouse models of E2a-PBX1-induced pre-B cell leukemia are still elusive. We now report the generation of a pre-B leukemia model using E2a-PBX1 transgenic mice, which lack mature and precursor T-cells as a result of engineered loss of CD3epsilon expression (CD3epsilon(-/-)). Using insertional mutagenesis and inverse-PCR, we show that B-cell leukemia development in the E2a-PBX1 x CD3epsilon(-/-) compound transgenic animals is significantly accelerated when compared to control littermates, and document several known and novel integrations in these tumors. Of all common integration sites, a small region of 19 kb in the Hoxa gene locus, mostly between Hoxa6 and Hoxa10, represented 18% of all integrations in the E2a-PBX1 B-cell leukemia and was targeted in 86% of these leukemias compared to 17% in control tumors. Q-PCR assessment of expression levels for most Hoxa cluster genes in these tumors revealed an unprecedented impact of the proviral integrations on Hoxa gene expression, with tumors having one to seven different Hoxa genes overexpressed at levels up to 6600-fold above control values. Together our studies set the stage for modeling E2a-PBX1-induced B-cell leukemia and shed new light on the complexity pertaining to Hox gene regulation. In addition, our results show that the Hoxa gene cluster is preferentially targeted in E2a-PBX1-induced tumors, thus suggesting functional collaboration between these oncogenes in pre-B-cell tumors.

The repertoire of G protein-coupled receptors in the human parasite Schistosoma mansoni and the model organism Schmidtea mediterranea.

Background

G protein-coupled receptors (GPCRs) constitute one of the largest groupings of eukaryotic proteins, and represent a particularly lucrative set of pharmaceutical targets. They play an important role in eukaryotic signal transduction and physiology, mediating cellular responses to a diverse range of extracellular stimuli. The phylum Platyhelminthes is of considerable medical and biological importance, housing major pathogens as well as established model organisms. The recent availability of genomic data for the human blood fluke Schistosoma mansoni and the model planarian Schmidtea mediterranea paves the way for the first comprehensive effort to identify and analyze GPCRs in this important phylum.

Results

Application of a novel transmembrane-oriented approach to receptor mining led to the discovery of 117 S. mansoni GPCRs, representing all of the major families; 105 Rhodopsin, 2 Glutamate, 3 Adhesion, 2 Secretin and 5 Frizzled. Similarly, 418 Rhodopsin, 9 Glutamate, 21 Adhesion, 1 Secretin and 11 Frizzled S. mediterranea receptors were identified. Among these, we report the identification of novel receptor groupings, including a large and highly-diverged Platyhelminth-specific Rhodopsin subfamily, a planarian-specific Adhesion-like family, and atypical Glutamate-like receptors. Phylogenetic analysis was carried out following extensive gene curation. Support vector machines (SVMs) were trained and used for ligand-based classification of full-length Rhodopsin GPCRs, complementing phylogenetic and homology-based classification.

Conclusions

Genome-wide investigation of GPCRs in two platyhelminth genomes reveals an extensive and complex receptor signaling repertoire with many unique features. This work provides important sequence and functional leads for understanding basic flatworm receptor biology, and sheds light on a lucrative set of anthelmintic drug targets.

Discovery of Novel Allosteric Modulators of Metabotropic Glutamate Receptor Subtype 5 Reveals Chemical and Functional Diversity and In Vivo Activity in Rat Behavioral Models of Anxiolytic and Antipsychotic Activity

Modulators of metabotropic glutamate receptor subtype 5 (mGluR5) may provide novel treatments for multiple central nervous system (CNS) disorders, including anxiety and schizophrenia. Although compounds have been developed to better understand the physiological roles of mGluR5 and potential usefulness for the treatment of these disorders, there are limitations in the tools available, including poor selectivity, low potency, and limited solubility. To address these issues, we developed an innovative assay that allows simultaneous screening for mGluR5 agonists, antagonists, and potentiators. We identified multiple scaffolds that possess diverse modes of activity at mGluR5, including both positive and negative allosteric modulators (PAMs and NAMs, respectively). 3-Fluoro-5-(3-(pyridine-2-yl)-1,2,4-oxadiazol-5-yl) benzonitrile (VU0285683) was developed as a novel selective mGluR5 NAM with high affinity for the 2-methyl-6-(phenyl-ethynyl)-pyridine (MPEP) binding site. VU0285683 had anxiolytic-like activity in two rodent models for anxiety but did not potentiate phen-cyclidine-induced hyperlocomotor activity. (4-Hydroxypiperidin-1-yl)(4-phenylethynyl) phenyl) methanone (VU0092273) was identified as a novel mGluR5 PAM that also binds to the MPEP site. VU0092273 was chemically optimized to an orally active analog, N-cyclobutyl-6-((3-fluorophenyl) ethynyl) nicotinamide hydrochloride (VU0360172), which is selective for mGluR5. This novel mGluR5 PAM produced a dose-dependent reversal of amphetamine-induced hyperlocomotion, a rodent model predictive of antipsychotic activity. Discovery of structurally and functionally diverse allosteric modulators of mGluR5 that demonstrate in vivo efficacy in rodent models of anxiety and antipsychotic activity provide further support for the tremendous diversity of chemical scaffolds and modes of efficacy of mGluR5 ligands. In addition, these studies provide strong support for the hypothesis that multiple structurally distinct mGluR5 modulators have robust activity in animal models that predict efficacy in the treatment of CNS disorders.

Functional analysis of predicted coiled-coil regions in the Escherichia coli K-12 O-antigen polysaccharide chain length determinant Wzz

Wzz is a membrane protein that determines the chain length distribution of the O-antigen lipopolysaccharide by an unknown mechanism. Wzz proteins consist of two transmembrane helices separated by a large periplasmic loop. The periplasmic loop of Escherichia coli K-12 Wzz (244 amino acids from K65 to A308) was purified and found to be a monomer with an extended conformation, as determined by gel filtration chromatography and analytical ultracentrifugation. Circular dichroism showed that the loop has a 60% helical content. The Wzz periplasmic loop also contains three regions with predicted coiled coils. To probe the function of the predicted coiled coils, we constructed amino acid replacement mutants of the E. coli K-12 Wzz protein, which were designed so that the coiled coils could be separate without compromising the helicity of the individual molecules. Mutations in one of the regions, spanning amino acids 108 to 130 (region I), were associated with a partial defect in O-antigen chain length distribution, while mutants with mutations in the region spanning amino acids 209 to 223 (region III) did not have an apparent functional defect. In contrast, mutations in the region spanning amino acids 153 to 173 (region II) eliminated the Wzz function. This phenotype was associated with protein instability, most likely due to conformational changes caused by the amino acid replacements, which was confirmed by limited trypsin proteolysis. Additional mutagenesis based on a three-dimensional model of region I demonstrated that the amino acids implicated in function are all located at the same face of a predicted alpha-helix, suggesting that a coiled coil actually does not exist in this region. Together, our results suggest that the regions predicted to be coiled coils are important for Wzz function because they maintain the native conformation of the protein, although the existence of coiled coils could not be demonstrated experimentally.

Functional analysis of the glycero-manno-heptose 7-phosphate kinase domain from the bifunctional HldE protein, which is involved in ADP-L-glycero-D-manno-heptose biosynthesis

The core oligosaccharide component of the lipopolysaccharide can be subdivided into inner and outer core regions. In Escherichia coli, the inner core consists of two 3-deoxy-d-manno-octulosonic acid and three glycero-manno-heptose residues. The HldE protein participates in the biosynthesis of ADP-glycero-manno-heptose precursors used in the assembly of the inner core. HldE comprises two functional domains: an N-terminal region with homology to the ribokinase superfamily (HldE1 domain) and a C-terminal region with homology to the cytidylyltransferase superfamily (HldE2 domain). We have employed the structure of the E. coli ribokinase as a template to model the HldE1 domain and predict critical amino acids required for enzyme activity. Mutation of these residues renders the protein inactive as determined in vivo by functional complementation analysis. However, these mutations did not affect the secondary or tertiary structure of purified HldE1, as judged by fluorescence spectroscopy and circular dichroism. Furthermore, in vivo coexpression of wild-type, chromosomally encoded HldE and mutant HldE1 proteins with amino acid substitutions in the predicted ATP binding site caused a dominant negative phenotype as revealed by increased bacterial sensitivity to novobiocin. Copurification experiments demonstrated that HldE and HldE1 form a complex in vivo. Gel filtration chromatography resulted in the detection of a dimer as the predominant form of the native HldE1 protein. Altogether, our data support the notions that the HldE functional unit is a dimer and that structural components present in each HldE1 monomer are required for enzymatic activity.

Bidirectional, iterative approach to the structural delineation of the functional "Chemoprint" in GPR40 for agonist recognition

GPR40, free fatty acid receptor 1 (FFAR1), is a member of the GPCR superfamily and a possible target for the treatment of type 2 diabetes. In this work, we conducted a bidirectional iterative investigation, including computational modeling and site-directed mutagenesis, aimed at delineating amino acid residues forming the functional "chemoprint" of GPR40 for agonist recognition. The computational and experimental studies revolved around the recognition of the potent synthetic agonist GW9508. Our experimentally supported model suggested that H137(4.56), R183(5.39), N244(6.55), and R258(7.35) are directly involved in interactions with the ligand. We have proposed a polarized NH-pi interaction between H137(4.56) and GW9508 as one of the contributing forces leading to the high potency of GW9508. The modeling approach presented in this work provides a general strategy for the exploration of receptor-ligand interactions in G-protein coupled receptors beginning prior to acquisition of experimental data.

Optimal diet choice for large herbivores: an extended contingency model

1. A more general contingency model of optimal diet choice is developed, allowing for simultaneous searching and handling, which extends the theory to include grazing and browsing by large herbivores.</p><p>2. Foraging resolves into three modes: purely encounter-limited, purely handling-limited and mixed-process, in which either a handling-limited prey type is added to an encounter-limited diet, or the diet becomes handling-limited as it expands.</p><p>3. The purely encounter-limited diet is, in general, broader than that predicted by the conventional contingency model,</p><p>4. As the degree of simultaneity of searching and handling increases, the optimal diet expands to the point where it is handling-limited, at which point all inferior prey types are rejected,</p><p>5. Inclusion of a less profitable prey species is not necessarily independent of its encounter rate and the zero-one rule does not necessarily hold: some of the less profitable prey may be included in the optimal diet. This gives an optimal foraging explanation for herbivores' mixed diets.</p><p>6. Rules are shown for calculating the boundary between encounter-limited and handling-limited diets and for predicting the proportion of inferior prey to be included in a two-species diet,</p><p>7. The digestive rate model is modified to include simultaneous searching and handling, showing that the more they overlap, the more the predicted diet-breadth is likely to be reduced.</p>

Measles immune suppression:lessons from the Macaque Model

Measles remains a significant childhood disease, and is associated with a transient immune suppression. Paradoxically, measles virus (MV) infection also induces robust MV-specific immune responses. Current hypotheses for the mechanism underlying measles immune suppression focus on functional impairment of lymphocytes or antigen-presenting cells, caused by infection with or exposure to MV. We have generated stable recombinant MVs that express enhanced green fluorescent protein, and remain virulent in non-human primates. By performing a comprehensive study of virological, immunological, hematological and histopathological observations made in animals euthanized at different time points after MV infection, we developed a model explaining measles immune suppression which fits with the "measles paradox". Here we show that MV preferentially infects CD45RA - memory T-lymphocytes and follicular B-lymphocytes, resulting in high infection levels in these populations. After the peak of viremia MV-infected lymphocytes were cleared within days, followed by immune activation and lymph node enlargement. During this period tuberculin-specific T-lymphocyte responses disappeared, whilst strong MV-specific T-lymphocyte responses emerged. Histopathological analysis of lymphoid tissues showed lymphocyte depletion in the B- and T-cell areas in the absence of apoptotic cells, paralleled by infiltration of T-lymphocytes into B-cell follicles and reappearance of proliferating cells. Our findings indicate an immune-mediated clearance of MV-infected CD45RA - memory T-lymphocytes and follicular B-lymphocytes, which causes temporary immunological amnesia. The rapid oligoclonal expansion of MV-specific lymphocytes and bystander cells masks this depletion, explaining the short duration of measles lymphopenia yet long duration of immune suppression. © 2012 de Vries et al.

Density Functional Theory Study on the Cleavage Mechanism of the Carbonyl Bond in Amides on Flat and Stepped Ru Surfaces: Hydrogen-Induced or Direct C-O Bond Breaking?

We have performed density functional theory (DFT) calculations to investigate the reaction mechanism of the cleavage of the carbonyl bond in amides on both flat and stepped Ru surfaces. The simplest amide molecule, N,N-dimethylacetamide (DMA), was used as the exemplar model molecule. Through the calculations, the most stable transition states (TSs) in all the pathways on both flat and stepped Ru surfaces are identified. Comparing the energy profiles of different reaction pathways, we find that a direct cleavage mechanism is always energetically favored as compared with an alternative hydrogen-induced mechanism on either the flat or stepped Ru surface. It is easier for the dissociation process to occur on the stepped surface than on the flat surface. However, as compared with the terrace, the superiority of step sites boosting the C-O bond dissociation is not as evident as that on CO dissociation.

Identifying functional stakeholder clusters to maximise communication for the ecosystem approach to fisheries management

Interaction with ecological models can improve stakeholder participation in fisheries management. Problems exist in efficiently communicating outputs to stakeholders and an objective method of structuring stakeholder differences is lacking. This paper aims to inform the design of a multi-user communication interface for fisheries management by identifying functional stakeholder groups. Intuitive categorisation of stakeholders, derived from survey responses, is contrasted with an Evidence-Based method derived from analysis of stakeholder literature. Intuitive categorisation relies on interpretation and professional judgement when categorising stakeholders among conventional stakeholder groups. Evidence-Based categorisation quantitatively characterises each stakeholder with a vector of four management objective interest-strength values (Yield, Employment, Profit and Ecosystem Preservation). Survey respondents agreed little in forming intuitive groups and the groups were poorly defined and heterogeneous in interests. In contrast the Evidence-Based clusters were well defined and largely homogeneous, so more useful for identifying functional relations with model outputs. The categorisations lead to two different clusterings of stakeholders and suggest unhelpful stereotyping of stakeholders may occur with the Intuitive categorisation method. Stakeholder clusters based on literature-evidence show a high degree of common interests among clusters and is encouraging for those seeking to maximise dialogue and consensus forming. © 2013 Elsevier Ltd.

Macroscopic limit of time-dependent density-functional theory for adiabatic local approximations of the exchange-correlation kernel

Time-dependent density-functional theory is a rather accurate and efficient way to compute electronic excitations for finite systems. However, in the macroscopic limit (systems of increasing size), for the usual adiabatic random-phase, local-density, or generalized-gradient approximations, one recovers the Kohn-Sham independent-particle picture, and thus the incorrect band gap. To clarify this trend, we investigate the macroscopic limit of the exchange-correlation kernel in such approximations by means of an algebraical analysis complemented with numerical studies of a one-dimensional tight-binding model. We link the failure to shift the Kohn-Sham spectrum of these approximate kernels to the fact that the corresponding operators in the transition space act only on a finite subspace.

Density Functional study of the photoactive yellow protein's chromophore

We have investigated the structural and electronic properties of p-coumaric acid, the chromophore of the photoactive yellow protein (PYP), by means of first-principles molecular dynamics based on density functional theory (DFT). We have studied the chromophore both in the vacuum and in an extended model which includes the nearest residues in the binding pocket of PYP, as derived from crystallographic data. We have characterized the ground state of the isolated chromophore in its protonated and deprotonated forms and computed the energy barrier involved in the trans to cis isomerization process around the carbon-carbon double bond. A comparison of the optimized structures of the chromophore in the vacuum and in the extended protein model, both in the trans (ground state of PYP in the dark) and cis (first light-activated intermediate) configuration, shows how the protein environment affects the chromophore in the first step of the photocycle. Our model gives an energy storage of 25 kcal/mol associated with the trans-to-cia photoisomerization. Finally, we have elucidated the nature of the electronic excitation relevant for the photochemistry of PYP by means of time-dependent DFT calculations.