Agonist binding, agonist affinity and agonist efficacy at G protein-coupled receptors

Measurements of affinity and efficacy are fundamental for work on agonists both in drug discovery and in basic studies on receptors. In this review I wish to consider methods for measuring affinity and efficacy at G protein coupled receptors (GPCRs). Agonist affinity may be estimated in terms of the dissociation constant for agonist binding to a receptor using ligand binding or functional assays. It has, however, been suggested that measurements of affinity are always contaminated by efficacy so that it is impossible to separate the two parameters. Here I show that for many GPCRs, if receptor/G protein coupling is suppressed, experimental measurements of agonist affinity using ligand binding (K-obs) provide quite accurate measures of the agonist microscopic dissociation constant (K-A). Also in pharmacological functional studies, good estimates of agonist dissociation constants are possible. Efficacy can be quantitated in several ways based on functional data ( maximal effect of the agonist (E-max), ratio of agonist dissociation constant to concentration of agonist giving half maximal effect in functional assay ( K-obs/ EC50), a combined parameter EmaxKobs/EC50). Here I show that EmaxKobs/EC50 provides the best assessment of efficacy for a range of agonists across the full range of efficacy for full to partial agonists. Considerable evidence now suggests that ligand efficacy may be dependent on the pathway used to assess it. The efficacy of a ligand may, therefore, be multidimensional. It is still, however, necessary to have accurate measures of efficacy in different pathways.

Towards a model-based integration of co-registered electroencephalography/functional magnetic resonance imaging data with realistic neural population meshes

Brain activity can be measured with several non-invasive neuroimaging modalities, but each modality has inherent limitations with respect to resolution, contrast and interpretability. It is hoped that multimodal integration will address these limitations by using the complementary features of already available data. However, purely statistical integration can prove problematic owing to the disparate signal sources. As an alternative, we propose here an advanced neural population model implemented on an anatomically sound cortical mesh with freely adjustable connectivity, which features proper signal expression through a realistic head model for the electroencephalogram (EEG), as well as a haemodynamic model for functional magnetic resonance imaging based on blood oxygen level dependent contrast (fMRI BOLD). It hence allows simultaneous and realistic predictions of EEG and fMRI BOLD from the same underlying model of neural activity. As proof of principle, we investigate here the influence on simulated brain activity of strengthening visual connectivity. In the future we plan to fit multimodal data with this neural population model. This promises novel, model-based insights into the brain's activity in sleep, rest and task conditions.

Electronic structure of point defects in controlled self-doping of the TiO2 (110) surface: Combined photoemission spectroscopy and density functional theory study

Point defects in metal oxides such as TiO2 are key to their applications in numerous technologies. The investigation of thermally induced nonstoichiometry in TiO2 is complicated by the difficulties in preparing and determining a desired degree of nonstoichiometry. We study controlled self-doping of TiO2 by adsorption of 1/8 and 1/16 monolayer Ti at the (110) surface using a combination of experimental and computational approaches to unravel the details of the adsorption process and the oxidation state of Ti. Upon adsorption of Ti, x-ray and ultraviolet photoemission spectroscopy (XPS and UPS) show formation of reduced Ti. Comparison of pure density functional theory (DFT) with experiment shows that pure DFT provides an inconsistent description of the electronic structure. To surmount this difficulty, we apply DFT corrected for on-site Coulomb interaction (DFT+U) to describe reduced Ti ions. The optimal value of U is 3 eV, determined from comparison of the computed Ti 3d electronic density of states with the UPS data. DFT+U and UPS show the appearance of a Ti 3d adsorbate-induced state at 1.3 eV above the valence band and 1.0 eV below the conduction band. The computations show that the adsorbed Ti atom is oxidized to Ti2+ and a fivefold coordinated surface Ti atom is reduced to Ti3+, while the remaining electron is distributed among other surface Ti atoms. The UPS data are best fitted with reduced Ti2+ and Ti3+ ions. These results demonstrate that the complexity of doped metal oxides is best understood with a combination of experiment and appropriate computations.

Dynamic causal modeling applied to fMRI data shows high reliability

Sensitivity, specificity, and reproducibility are vital to interpret neuroscientific results from functional magnetic resonance imaging (fMRI) experiments. Here we examine the scan–rescan reliability of the percent signal change (PSC) and parameters estimated using Dynamic Causal Modeling (DCM) in scans taken in the same scan session, less than 5 min apart. We find fair to good reliability of PSC in regions that are involved with the task, and fair to excellent reliability with DCM. Also, the DCM analysis uncovers group differences that were not present in the analysis of PSC, which implies that DCM may be more sensitive to the nuances of signal changes in fMRI data.

Modeling chemotaxis reveals the role of reversed phosphotransfer and a bi-functional kinase-phosphatase

Understanding how multiple signals are integrated in living cells to produce a balanced response is a major challenge in biology. Two-component signal transduction pathways, such as bacterial chemotaxis, comprise histidine protein kinases (HPKs) and response regulators (RRs). These are used to sense and respond to changes in the environment. Rhodobacter sphaeroides has a complex chemosensory network with two signaling clusters, each containing a HPK, CheA. Here we demonstrate, using a mathematical model, how the outputs of the two signaling clusters may be integrated. We use our mathematical model supported by experimental data to predict that: (1) the main RR controlling flagellar rotation, CheY6, aided by its specific phosphatase, the bifunctional kinase CheA3, acts as a phosphate sink for the other RRs; and (2) a phosphorelay pathway involving CheB2 connects the cytoplasmic cluster kinase CheA3 with the polar localised kinase CheA2, and allows CheA3-P to phosphorylate non-cognate chemotaxis RRs. These two mechanisms enable the bifunctional kinase/phosphatase activity of CheA3 to integrate and tune the sensory output of each signaling cluster to produce a balanced response. The signal integration mechanisms identified here may be widely used by other bacteria, since like R. sphaeroides, over 50% of chemotactic bacteria have multiple cheA homologues and need to integrate signals from different sources.

A hierarchical Bayesian model for predicting the functional consequences of amino-acid polymorphisms

Genetic polymorphisms in deoxyribonucleic acid coding regions may have a phenotypic effect on the carrier, e.g. by influencing susceptibility to disease. Detection of deleterious mutations via association studies is hampered by the large number of candidate sites; therefore methods are needed to narrow down the search to the most promising sites. For this, a possible approach is to use structural and sequence-based information of the encoded protein to predict whether a mutation at a particular site is likely to disrupt the functionality of the protein itself. We propose a hierarchical Bayesian multivariate adaptive regression spline (BMARS) model for supervised learning in this context and assess its predictive performance by using data from mutagenesis experiments on lac repressor and lysozyme proteins. In these experiments, about 12 amino-acid substitutions were performed at each native amino-acid position and the effect on protein functionality was assessed. The training data thus consist of repeated observations at each position, which the hierarchical framework is needed to account for. The model is trained on the lac repressor data and tested on the lysozyme mutations and vice versa. In particular, we show that the hierarchical BMARS model, by allowing for the clustered nature of the data, yields lower out-of-sample misclassification rates compared with both a BMARS and a frequen-tist MARS model, a support vector machine classifier and an optimally pruned classification tree.

A functional genomics approach reveals novel quantitative trait loci associated with platelet signalling pathways

Platelet response to activation varies widely between individuals but shows interindividual consistency and strong heritability. The genetic basis of this variation has not been properly explored. We therefore systematically measured the effect on function of sequence variation in 97 candidate genes in the collagen and adenosine-diphosphate (ADP) signaling pathways. Resequencing of the genes in 48 European DNA samples nearly doubled the number of known single nucleotide polymorphisms (SNPs) and informed the selection of 1327 SNPs for genotyping in 500 healthy Northern European subjects with known platelet responses to collagen-related peptide (CRP-XL) and ADP. This identified 17 novel associations with platelet function (P < .005) accounting for approximately 46% of the variation in response. Further investigations with platelets of known genotype explored the mechanisms behind some of the associations. SNPs in PEAR1 associated with increased platelet response to CRP-XL and increased PEAR1 protein expression after platelet degranulation. The minor allele of a 3' untranslated region (UTR) SNP (rs2769668) in VAV3 was associated with higher protein expression (P = .03) and increased P-selectin exposure after ADP activation (P = .004). Furthermore the minor allele of the intronic SNP rs17786144 in ITPR1 modified Ca2+ levels after activation with ADP (P < .004). These data provide novel insights into key hubs within platelet signaling networks.

Proteornics analysis unravels the functional repertoire of coronavirus nonstructural protein 3

Severe acute respiratory syndrome (SARS) coronavirus infection and growth are dependent on initiating signaling and enzyme actions upon viral entry into the host cell. Proteins packaged during virus assembly may subsequently form the first line of attack and host manipulation upon infection. A complete characterization of virion components is therefore important to understanding the dynamics of early stages of infection. Mass spectrometry and kinase profiling techniques identified nearly 200 incorporated host and viral proteins. We used published interaction data to identify hubs of connectivity with potential significance for virion formation. Surprisingly, the hub with the most potential connections was not the viral M protein but the nonstructurall protein 3 (nsp3), which is one of the novel virion components identified by mass spectrometry. Based on new experimental data and a bioinformatics analysis across the Coronaviridae, we propose a higher-resolution functional domain architecture for nsp3 that determines the interaction capacity of this protein. Using recombinant protein domains expressed in Escherichia coli, we identified two additional RNA-binding domains of nsp3. One of these domains is located within the previously described SARS-unique domain, and there is a nucleic acid chaperone-like domain located immediately downstream of the papain-like proteinase domain. We also identified a novel cysteine-coordinated metal ion-binding domain. Analyses of interdomain interactions and provisional functional annotation of the remaining, so-far-uncharacterized domains are presented. Overall, the ensemble of data surveyed here paint a more complete picture of nsp3 as a conserved component of the viral protein processing machinery, which is intimately associated with viral RNA in its role as a virion component.

Constrained models of evolution lead to improved prediction of functional linkage from correlated gain and loss of genes

Motivation: We compare phylogenetic approaches for inferring functional gene links. The approaches detect independent instances of the correlated gain and loss of pairs of genes from species' genomes. We investigate the effect on results of basing evidence of correlations on two phylogenetic approaches, Dollo parsminony and maximum likelihood (ML). We further examine the effect of constraining the ML model by fixing the rate of gene gain at a low value, rather than estimating it from the data. Results: We detect correlated evolution among a test set of pairs of yeast (Saccharomyces cerevisiae) genes, with a case study of 21 eukaryotic genomes and test data derived from known yeast protein complexes. If the rate at which genes are gained is constrained to be low, ML achieves by far the best results at detecting known functional links. The model then has fewer parameters but it is more realistic by preventing genes from being gained more than once. Availability: BayesTraits by M. Pagel and A. Meade, and a script to configure and repeatedly launch it by D. Barker and M. Pagel, are available at http://www.evolution.reading.ac.uk .

Functional and molecular characterization of Tachykinins and tachykinin receptors in the mouse uterus

The aim of this study was to analyze the function and expression of tachykinins, tachykinin receptors, and neprilysin (NEP) in the mouse uterus. A previous study showed that the uterotonic effects of substance P (SP), neurokinin A (NKA), and neurokinin B (NKB) in estrogen-treated mice were mainly mediated by the tachykinin NK, receptor. In the present work, further contractility studies were undertaken to determine the nature of the receptors mediating responses to tachykinins in uteri of late pregnant mice. Endpoint and real-time quantitative RTPCR were used to analyze the expression of the genes that encode the tachykinins SP/NKA, NKB, and hemokinin-1 (HK-1) (Tac1, Tac2, and Tac4); and the genes that encode tachykinin NK1 (Tacr1), NK2 (Tacr2), and NK3 (Tacr3) receptors in uteri from pregnant and nonpregnant mice. The data show that the mRNAs of tachykinins (particularly NKB and HK-1), tachykinin receptors, and NEP are locally expressed in the mouse uterus, and their expression changes during the estrous cycle and during pregnancy. The tachykinin INK, receptor is the predominant tachykinin receptor in the nonpregnant and early pregnant mouse and may mediate tachykinin-induced uterine contractions in the nonpregnant mouse. The tachykinin NK, receptor is predominant in the late pregnant mouse and is the main receptor mediating uterotonic responses to tachykinins at late pregnancy. The tachykinin NK, receptor is expressed in considerable amounts only in uteri from nonpregnant diestrous animals, and its physiological significance remains to be clarified.

Diet, immunity and functional foods

Functional foods (specific nutrient and/or food components) should beneficially affect one or more target functions in the body. The use of functional foods as a form of preventive medicine has been the subject of much research over the last two decades. It is well known that nutrition plays a vital role in chronic diseases, but it is only recently that data relating to the effects of specific nutrients or foods on the immune system have become available. This chapter aims to summarize the effects of some functional foods (e.g., prebiotics and micronutrients) on the immune system. It should be noted, however, that studies into the role of functional foods with regard to the human immune system are still in their infancy and a great deal of controversy surrounds the health claims attributed to some functional foods. Consequently, thorough studies are required in human and animal systems if we are to move towards developing a functional diet that provides maximal health benefits.

Evidence that GABA rho subunits contribute to functional ionotropic GABA receptors in mouse cerebellar Purkinje cells

Ionotropic gamma-amino butyric acid (GABA) receptors composed of heterogeneous molecular subunits are major mediators of inhibitory responses in the adult CNS. Here, we describe a novel ionotropic GABA receptor in mouse cerebellar Purkinje cells (PCs) using agents reported to have increased affinity for rho subunit-containing GABA(C) over other GABA receptors. Exogenous application of the GABA(C)-preferring agonist cis-4-aminocrotonic acid (CACA) evoked whole-cell currents in PCs, whilst equimolar concentrations of GABA evoked larger currents. CACA-evoked currents had a greater sensitivity to the selective GABA(C) antagonist (1,2,5,6-tetrahydropyridin-4-yl)methylphosphinic acid (TPMPA) than GABA-evoked currents. Focal application of agonists produced a differential response profile; CACA-evoked currents displayed a much more pronounced attenuation with increasing distance from the PC soma, displayed a slower time-to-peak and exhibited less desensitization than GABA-evoked currents. However, CACA-evoked currents were also completely blocked by bicuculline, a selective agent for GABA(A) receptors. Thus, we describe a population of ionotropic GABA receptors with a mixed GABA(A)/GABA(C) pharmacology. TPMPA reduced inhibitory synaptic transmission at interneurone-Purkinje cell (IN-PC) synapses, causing clear reductions in miniature inhibitory postsynaptic current (mIPSC) amplitude and frequency. Combined application of NO-711 (a selective GABA transporter subtype 1 (GAT-1) antagonist) and SNAP-5114 (a GAT-(2)/3/4 antagonist) induced a tonic GABA conductance in PCs; however, TPMPA had no effect on this current. Immunohistochemical studies suggest that rho subunits are expressed predominantly in PC soma and proximal dendritic compartments with a lower level of expression in more distal dendrites; this selective immunoreactivity contrasted with a more uniform distribution of GABA(A) alpha 1 subunits in PCs. Finally, co-immunoprecipitation studies suggest that rho subunits can form complexes with GABA(A) receptor alpha 1 subunits in the cerebellar cortex. Overall, these data suggest that rho subunits contribute to functional ionotropic receptors that mediate a component of phasic inhibitory GABAergic transmission at IN-PC synapses in the cerebellum.

Efficient modeling and inference for event-related fMRI data

Event-related functional magnetic resonance imaging (efMRI) has emerged as a powerful technique for detecting brains' responses to presented stimuli. A primary goal in efMRI data analysis is to estimate the Hemodynamic Response Function (HRF) and to locate activated regions in human brains when specific tasks are performed. This paper develops new methodologies that are important improvements not only to parametric but also to nonparametric estimation and hypothesis testing of the HRF. First, an effective and computationally fast scheme for estimating the error covariance matrix for efMRI is proposed. Second, methodologies for estimation and hypothesis testing of the HRF are developed. Simulations support the effectiveness of our proposed methods. When applied to an efMRI dataset from an emotional control study, our method reveals more meaningful findings than the popular methods offered by AFNI and FSL. (C) 2008 Elsevier B.V. All rights reserved.

The neural response to emotional prosody, as revealed by functional magnetic resonance imaging

Prosody is an important feature of language, comprising intonation, loudness, and tempo. Emotional prosodic processing forms an integral part of our social interactions. The main aim of this study was to use bold contrast fMRI to clarify the normal functional neuroanatomy of emotional prosody, in passive and active contexts. Subjects performed six separate scanning studies, within which two different conditions were contrasted: (1) "pure" emotional prosody versus rest; (2) congruent emotional prosody versus 'neutral' sentences; (3) congruent emotional prosody versus rest; (4) incongruent emotional prosody versus rest; (5) congruent versus incongruent emotional prosody; and (6) an active experiment in which subjects were instructed to either attend to the emotion conveyed by semantic content or that conveyed by tone of voice. Data resulting from these contrasts were analysed using SPM99. Passive listening to emotional prosody consistently activated the lateral temporal lobe (superior and/or middle temporal gyri). This temporal lobe response was relatively right-lateralised with or without semantic information. Both the separate and direct comparisons of congruent and incongruent emotional prosody revealed that subjects used fewer brain regions to process incongruent emotional prosody than congruent. The neural response to attention to semantics, was left lateralised, and recruited an extensive network not activated by attention to emotional prosody. Attention to emotional prosody modulated the response to speech, and induced right-lateralised activity, including the middle temporal gyrus. In confirming the results of lesion and neuropsychological studies, the current study emphasises the importance of the right hemisphere in the processing of emotional prosody, specifically the lateral temporal lobes. (C) 2003 Elsevier Science Ltd. All rights reserved.

Development and experimental identification of a biomechanical model of the trunk for functional electrical stimulation control in paraplegia

Objectives. Theoretic modeling and experimental studies suggest that functional electrical stimulation (FES) can improve trunk balance in spinal cord injured subjects. This can have a positive impact on daily life, increasing the volume of bimanual workspace, improving sitting posture, and wheelchair propulsion. A closed loop controller for the stimulation is desirable, as it can potentially decrease muscle fatigue and offer better rejection to disturbances. This paper proposes a biomechanical model of the human trunk, and a procedure for its identification, to be used for the future development of FES controllers. The advantage over previous models resides in the simplicity of the solution proposed, which makes it possible to identify the model just before a stimulation session ( taking into account the variability of the muscle response to the FES). Materials and Methods. The structure of the model is based on previous research on FES and muscle physiology. Some details could not be inferred from previous studies, and were determined from experimental data. Experiments with a paraplegic volunteer were conducted in order to measure the moments exerted by the trunk-passive tissues and artificially stimulated muscles. Data for model identification and validation also were collected. Results. Using the proposed structure and identification procedure, the model could adequately reproduce the moments exerted during the experiments. The study reveals that the stimulated trunk extensors can exert maximal moment when the trunk is in the upright position. In contrast, previous studies show that able-bodied subjects can exert maximal trunk extension when flexed forward. Conclusions. The proposed model and identification procedure are a successful first step toward the development of a model-based controller for trunk FES. The model also gives information on the trunk in unique conditions, normally not observable in able-bodied subjects (ie, subject only to extensor muscles contraction).