Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

On the functional relationship between tariffs and welfare

This paper uses a model of trade in two commodities between two countries to establish the following proposition. If the foreign offer curve has no points of inflection and if for each home rate of duty the equilibrium most favorable to the home country is selected (or else there is only one equilibrium), then as the rate of duty increases from zero, home welfare first rises then declines while foreign welfare steadily falls. © 1975.

Functional neuroimaging of autobiographical memory.

Functional neuroimaging studies of autobiographical memory have grown dramatically in recent years. These studies are important because they can investigate the neural correlates of processes that are difficult to study using laboratory stimuli, including: (i) complex constructive processes, (ii) recollective qualities of emotion and vividness, and (iii) remote memory retrieval. Constructing autobiographical memories involves search, monitoring and self-referential processes that are associated with activity in separable prefrontal regions. The contributions of emotion and vividness have been linked to the amygdala and visual cortex respectively. Finally, there is evidence that recent and remote autobiographical memories might activate the hippocampus equally, which has implications for memory-consolidation theories. The rapid development of innovative methods for eliciting personal memories in the scanner provides the opportunity to delve into the functional neuroanatomy of our personal past.

Functional properties of cell-seeded three-dimensionally woven poly(epsilon-caprolactone) scaffolds for cartilage tissue engineering.

Articular cartilage possesses complex mechanical properties that provide healthy joints the ability to bear repeated loads and maintain smooth articulating surfaces over an entire lifetime. In this study, we utilized a fiber-reinforced composite scaffold designed to mimic the anisotropic, nonlinear, and viscoelastic biomechanical characteristics of native cartilage as the basis for developing functional tissue-engineered constructs. Three-dimensionally woven poly(epsilon-caprolactone) (PCL) scaffolds were encapsulated with a fibrin hydrogel, seeded with human adipose-derived stem cells, and cultured for 28 days in chondrogenic culture conditions. Biomechanical testing showed that PCL-based constructs exhibited baseline compressive and shear properties similar to those of native cartilage and maintained these properties throughout the culture period, while supporting the synthesis of a collagen-rich extracellular matrix. Further, constructs displayed an equilibrium coefficient of friction similar to that of native articular cartilage (mu(eq) approximately 0.1-0.3) over the prescribed culture period. Our findings show that three-dimensionally woven PCL-fibrin composite scaffolds can be produced with cartilage-like mechanical properties, and that these engineered properties can be maintained in culture while seeded stem cells regenerate a new, functional tissue construct.

Structural manifestation of the delocalization error of density functional approximations: C(4N+2) rings and C(20) bowl, cage, and ring isomers.

The ground state structure of C(4N+2) rings is believed to exhibit a geometric transition from angle alternation (N < or = 2) to bond alternation (N > 2). All previous density functional theory (DFT) studies on these molecules have failed to reproduce this behavior by predicting either that the transition occurs at too large a ring size, or that the transition leads to a higher symmetry cumulene. Employing the recently proposed perspective of delocalization error within DFT we rationalize this failure of common density functional approximations (DFAs) and present calculations with the rCAM-B3LYP exchange-correlation functional that show an angle-to-bond-alternation transition between C(10) and C(14). The behavior exemplified here manifests itself more generally as the well known tendency of DFAs to bias toward delocalized electron distributions as favored by Huckel aromaticity, of which the C(4N+2) rings provide a quintessential example. Additional examples are the relative energies of the C(20) bowl, cage, and ring isomers; we show that the results from functionals with minimal delocalization error are in good agreement with CCSD(T) results, in contrast to other commonly used DFAs. An unbiased DFT treatment of electron delocalization is a key for reliable prediction of relative stability and hence the structures of complex molecules where many structure stabilization mechanisms exist.

Reformulating time-dependent density functional theory with non-orthogonal localized molecular orbitals.

Time-dependent density functional theory (TDDFT) has broad application in the study of electronic response, excitation and transport. To extend such application to large and complex systems, we develop a reformulation of TDDFT equations in terms of non-orthogonal localized molecular orbitals (NOLMOs). NOLMO is the most localized representation of electronic degrees of freedom and has been used in ground state calculations. In atomic orbital (AO) representation, the sparsity of NOLMO is transferred to the coefficient matrix of molecular orbitals (MOs). Its novel use in TDDFT here leads to a very simple form of time propagation equations which can be solved with linear-scaling effort. We have tested the method for several long-chain saturated and conjugated molecular systems within the self-consistent charge density-functional tight-binding method (SCC-DFTB) and demonstrated its accuracy. This opens up pathways for TDDFT applications to large bio- and nano-systems.

A novel, non-apoptotic role for Scythe/BAT3: a functional switch between the pro- and anti-proliferative roles of p21 during the cell cycle.

BACKGROUND: Scythe/BAT3 is a member of the BAG protein family whose role in apoptosis has been extensively studied. However, since the developmental defects observed in Bat3-null mouse embryos cannot be explained solely by defects in apoptosis, we investigated whether BAT3 is also involved in cell-cycle progression. METHODS/PRINCIPAL FINDINGS: Using a stable-inducible Bat3-knockdown cellular system, we demonstrated that reduced BAT3 protein level causes a delay in both G1/S transition and G2/M progression. Concurrent with these changes in cell-cycle progression, we observed a reduction in the turnover and phosphorylation of the CDK inhibitor p21, which is best known as an inhibitor of DNA replication; however, phosphorylated p21 has also been shown to promote G2/M progression. Our findings indicate that in Bat3-knockdown cells, p21 continues to be synthesized during cell-cycle phases that do not normally require p21, resulting in p21 protein accumulation and a subsequent delay in cell-cycle progression. Finally, we showed that BAT3 co-localizes with p21 during the cell cycle and is required for the translocation of p21 from the cytoplasm to the nucleus during the G1/S transition and G2/M progression. CONCLUSION: Our study reveals a novel, non-apoptotic role for BAT3 in cell-cycle regulation. By maintaining a low p21 protein level during the G1/S transition, BAT3 counteracts the inhibitory effect of p21 on DNA replication and thus enables the cells to progress from G1 to S phase. Conversely, during G2/M progression, BAT3 facilitates p21 phosphorylation by cyclin A/Cdk2, an event required for G2/M progression. BAT3 modulates these pro- and anti-proliferative roles of p21 at least in part by regulating cyclin A abundance, as well as p21 translocation between the cytoplasm and the nucleus to ensure that it functions in the appropriate intracellular compartment during each phase of the cell cycle.

Functional evolution of mammalian odorant receptors.

The mammalian odorant receptor (OR) repertoire is an attractive model to study evolution, because ORs have been subjected to rapid evolution between species, presumably caused by changes of the olfactory system to adapt to the environment. However, functional assessment of ORs in related species remains largely untested. Here we investigated the functional properties of primate and rodent ORs to determine how well evolutionary distance predicts functional characteristics. Using human and mouse ORs with previously identified ligands, we cloned 18 OR orthologs from chimpanzee and rhesus macaque and 17 mouse-rat orthologous pairs that are broadly representative of the OR repertoire. We functionally characterized the in vitro responses of ORs to a wide panel of odors and found similar ligand selectivity but dramatic differences in response magnitude. 87% of human-primate orthologs and 94% of mouse-rat orthologs showed differences in receptor potency (EC50) and/or efficacy (dynamic range) to an individual ligand. Notably dN/dS ratio, an indication of selective pressure during evolution, does not predict functional similarities between orthologs. Additionally, we found that orthologs responded to a common ligand 82% of the time, while human OR paralogs of the same subfamily responded to the common ligand only 33% of the time. Our results suggest that, while OR orthologs tend to show conserved ligand selectivity, their potency and/or efficacy dynamically change during evolution, even in closely related species. These functional changes in orthologs provide a platform for examining how the evolution of ORs can meet species-specific demands.

Relative growth of the limbs and trunk in sifakas: heterochronic, ecological, and functional considerations.

Limb, trunk, and body weight measurements were obtained for growth series of Milne-Edwards's diademed sifaka, Propithecus diadema edwardsi, and the golden-crowned sifaka, Propithecus tattersalli. Similar measures were obtained also for primarily adults of two subspecies of the western sifaka: Propithecus verreauxi coquereli, Coquerel's sifaka, and Propithecus verreauxi verreauxi, Verreaux's sifaka. Ontogenetic series for the larger-bodied P. d. edwardsi and the smaller-bodied P. tattersalli were compared to evaluate whether species-level differences in body proportions result from the differential extension of common patterns of relative growth. In bivariate plots, both subspecies of P. verreauxi were included to examine whether these taxa also lie along a growth trajectory common to all sifakas. Analyses of the data indicate that postcranial proportions for sifakas are ontogenetically scaled, much as demonstrated previously with cranial dimensions for all three species (Ravosa, 1992). As such, P. d. edwardsi apparently develops larger overall size primarily by growing at a faster rate, but not for a longer duration of time, than P. tattersalli and P. verreauxi; this is similar to results based on cranial data. A consideration of Malagasy lemur ecology suggests that regional differences in forage quality and resource availability have strongly influenced the evolutionary development of body-size variation in sifakas. On one hand, the rainforest environment of P. d. edwardsi imposes greater selective pressures for larger body size than the dry-forest environment of P. tattersalli and P. v. coquereli, or the semi-arid climate of P. v. verreauxi. On the other hand, as progressively smaller-bodied adult sifakas are located in the east, west, and northwest, this apparently supports suggestions that adult body size is set by dry-season constraints on food quality and distribution (i.e., smaller taxa are located in more seasonal habitats such as the west and northeast). Moreover, the fact that body-size differentiation occurs primarily via differences in growth rate is also due apparently to differences in resource seasonality (and juvenile mortality risk in turn) between the eastern rainforest and the more temperate northeast and west. Most scaling coefficients for both arm and leg growth range from slight negative allometry to slight positive allometry. Given the low intermembral index for sifakas, which is also an adaptation for propulsive hindlimb-dominated jumping, this suggests that differences in adult limb proportions are largely set prenatally rather than being achieved via higher rates of postnatal hindlimb growth.(ABSTRACT TRUNCATED AT 400 WORDS)

Functional parcellation of attentional control regions of the brain.

Recently, a number of investigators have examined the neural loci of psychological processes enabling the control of visual spatial attention using cued-attention paradigms in combination with event-related functional magnetic resonance imaging. Findings from these studies have provided strong evidence for the involvement of a fronto-parietal network in attentional control. In the present study, we build upon this previous work to further investigate these attentional control systems. In particular, we employed additional controls for nonattentional sensory and interpretative aspects of cue processing to determine whether distinct regions in the fronto-parietal network are involved in different aspects of cue processing, such as cue-symbol interpretation and attentional orienting. In addition, we used shorter cue-target intervals that were closer to those used in the behavioral and event-related potential cueing literatures. Twenty participants performed a cued spatial attention task while brain activity was recorded with functional magnetic resonance imaging. We found functional specialization for different aspects of cue processing in the lateral and medial subregions of the frontal and parietal cortex. In particular, the medial subregions were more specific to the orienting of visual spatial attention, while the lateral subregions were associated with more general aspects of cue processing, such as cue-symbol interpretation. Additional cue-related effects included differential activations in midline frontal regions and pretarget enhancements in the thalamus and early visual cortical areas.

Resolving response, decision, and strategic control: evidence for a functional topography in dorsomedial prefrontal cortex.

The dorsomedial prefrontal cortex (DMPFC) plays a central role in aspects of cognitive control and decision making. Here, we provide evidence for an anterior-to-posterior topography within the DMPFC using tasks that evoke three distinct forms of control demands--response, decision, and strategic--each of which could be mapped onto independent behavioral data. Specifically, we identify three spatially distinct regions within the DMPFC: a posterior region associated with control demands evoked by multiple incompatible responses, a middle region associated with control demands evoked by the relative desirability of decision options, and an anterior region that predicts control demands related to deviations from an individual's preferred decision-making strategy. These results provide new insight into the functional organization of DMPFC and suggest how recent controversies about its role in complex decision making and response mapping can be reconciled.

Functional imaging of numerical processing in adults and 4-y-old children.

Adult humans, infants, pre-school children, and non-human animals appear to share a system of approximate numerical processing for non-symbolic stimuli such as arrays of dots or sequences of tones. Behavioral studies of adult humans implicate a link between these non-symbolic numerical abilities and symbolic numerical processing (e.g., similar distance effects in accuracy and reaction-time for arrays of dots and Arabic numerals). However, neuroimaging studies have remained inconclusive on the neural basis of this link. The intraparietal sulcus (IPS) is known to respond selectively to symbolic numerical stimuli such as Arabic numerals. Recent studies, however, have arrived at conflicting conclusions regarding the role of the IPS in processing non-symbolic, numerosity arrays in adulthood, and very little is known about the brain basis of numerical processing early in development. Addressing the question of whether there is an early-developing neural basis for abstract numerical processing is essential for understanding the cognitive origins of our uniquely human capacity for math and science. Using functional magnetic resonance imaging (fMRI) at 4-Tesla and an event-related fMRI adaptation paradigm, we found that adults showed a greater IPS response to visual arrays that deviated from standard stimuli in their number of elements, than to stimuli that deviated in local element shape. These results support previous claims that there is a neurophysiological link between non-symbolic and symbolic numerical processing in adulthood. In parallel, we tested 4-y-old children with the same fMRI adaptation paradigm as adults to determine whether the neural locus of non-symbolic numerical activity in adults shows continuity in function over development. We found that the IPS responded to numerical deviants similarly in 4-y-old children and adults. To our knowledge, this is the first evidence that the neural locus of adult numerical cognition takes form early in development, prior to sophisticated symbolic numerical experience. More broadly, this is also, to our knowledge, the first cognitive fMRI study to test healthy children as young as 4 y, providing new insights into the neurophysiology of human cognitive development.

Functional desensitization of the isolated beta-adrenergic receptor by the beta-adrenergic receptor kinase: potential role of an analog of the retinal protein arrestin (48-kDa protein).

The beta-adrenergic receptor kinase is an enzyme, possibly analogous to rhodopsin kinase, that multiply phosphorylates the beta-adrenergic receptor only when it is occupied by stimulatory agonists. Since this kinase may play an important role in mediating the process of homologous, or agonist-specific, desensitization, we investigated the functional consequences of receptor phosphorylation by the kinase and possible analogies with the mechanism of action of rhodopsin kinase. Pure hamster lung beta 2-adrenergic receptor, reconstituted in phospholipid vesicles, was assessed for its ability to mediate agonist-promoted stimulation of the GTPase activity of coreconstituted stimulatory guanine nucleotide-binding regulatory protein. When the receptor was phosphorylated by partially (approximately 350-fold) purified preparations of beta-adrenergic receptor kinase, as much as 80% inactivation of its functional activity was observed. However, the use of more highly purified enzyme preparations led to a dramatic decrease in the ability of phosphorylation to inactivate the receptor such that pure enzyme preparations (approximately 20,000-fold purified) caused only minimal (approximately 1off/- 7%) inactivation. Addition of pure retinal arrestin (48-kDa protein or S antigen), which is involved in enhancing the inactivating effect of rhodopsin phosphorylation by rhodopsin kinase, led to partial restoration of the functional effect of beta-adrenergic receptor kinase-promoted phosphorylation (41 +/- 3% inactivation). These results suggest the possibility that a protein analogous to retinal arrestin may exist in other tissues and function in concert with beta-adrenergic receptor kinase to regulate the activity of adenylate cyclase-coupled receptors.

cDNA for the human beta 2-adrenergic receptor: a protein with multiple membrane-spanning domains and encoded by a gene whose chromosomal location is shared with that of the receptor for platelet-derived growth factor.

We have isolated and sequenced a cDNA encoding the human beta 2-adrenergic receptor. The deduced amino acid sequence (413 residues) is that of a protein containing seven clusters of hydrophobic amino acids suggestive of membrane-spanning domains. While the protein is 87% identical overall with the previously cloned hamster beta 2-adrenergic receptor, the most highly conserved regions are the putative transmembrane helices (95% identical) and cytoplasmic loops (93% identical), suggesting that these regions of the molecule harbor important functional domains. Several of the transmembrane helices also share lesser degrees of identity with comparable regions of select members of the opsin family of visual pigments. We have localized the gene for the beta 2-adrenergic receptor to q31-q32 on chromosome 5. This is the same position recently determined for the gene encoding the receptor for platelet-derived growth factor and is adjacent to that for the FMS protooncogene, which encodes the receptor for the macrophage colony-stimulating factor.

Phosphorylation/dephosphorylation of the beta-adrenergic receptor regulates its functional coupling to adenylate cyclase and subcellular distribution.

Prolonged exposure of cells or tissues to drugs or hormones such as catecholamines leads to a state of refractoriness to further stimulation by that agent, known as homologous desensitization. In the case of the beta-adrenergic receptor coupled to adenylate cyclase, this process has been shown to be intimately associated with the sequestration of the receptors from the cell surface through a cAMP-independent process. Recently, we have shown that homologous desensitization in the frog erythrocyte model system is also associated with increased phosphorylation of the beta-adrenergic receptor. We now provide evidence that the phosphorylation state of the beta-adrenergic receptor regulates its functional coupling to adenylate cyclase, subcellular translocation, and recycling to the cell surface during the process of agonist-induced homologous desensitization. Moreover, we show that the receptor phosphorylation is reversed by a phosphatase specifically associated with the sequestered subcellular compartment. At 23 degrees C, the time courses of beta-adrenergic receptor phosphorylation, sequestration, and adenylate cyclase desensitization are identical, occurring without a lag, exhibiting a t1/2 of 30 min, and reaching a maximum at approximately 3 hr. Upon cell lysis, the sequestered beta-adrenergic receptors can be partially recovered in a light membrane vesicle fraction that is separable from the plasma membranes by differential centrifugation. The increased beta-adrenergic receptor phosphorylation is apparently reversed in the sequestered vesicle fraction as the sequestered receptors exhibit a phosphate/receptor stoichiometry that is similar to that observed under basal conditions. High levels of a beta-adrenergic receptor phosphatase activity appear to be associated with the sequestered vesicle membranes. The functional activity of the phosphorylated beta-adrenergic receptor was examined by reconstituting purified receptor with its biochemical effector the guanine nucleotide regulatory protein (Ns) in phospholipid vesicles and assessing the receptor-stimulated GTPase activity of Ns. Compared to controls, phosphorylated beta-adrenergic receptors, purified from desensitized cells, were less efficacious in activating the Ns GTPase activity. These results suggest that phosphorylation of the beta-adrenergic receptor leads to its functional uncoupling and physical translocation away from the cell surface into a sequestered membrane domain. In the sequestered compartment, the phosphorylation is reversed thus enabling the receptor to recycle back to the cell surface and recouple with adenylate cyclase.