Protein Sequence Modules

Abstract

Conserved protein sequence segments are commonly believed to correspond to functional sites in the protein sequence. A novel approach is proposed to profile the changing degree of conservation along the protein sequence, by evaluating the occurrence frequencies of all short oligopeptides of the given sequence in a large proteome database. Thus, a protein sequence conservation profile can be plotted for every protein. The profile indicates where along the sequences the potential functional (conserved) sites are located. The corresponding oligopeptides belonging to the sites are very frequent across many prokaryotic species. Analysis of a representative set of such profiles reveals a common feature of all examined proteins: they consist of sequence modules represented by the peaks of conservation. Typical size of the modules (peak-to-peak distance) is 25–30 amino acid residues.     

Persistence of Dispersal-limited Species in Structured Dynamic Landscapes

Dynamic landscape models have generally assumed random distributions of habitat although real landscapes show spatial organization at many scales. To explore the role of spatial structure in determining the frequency of dispersal-limited forest species, we used a cellular landscape model divided into two zones. Zones were distributed in a random, clustered, or regular spatial pattern. Within each zone habitat cells were randomly destroyed and regenerated, and habitat density and turnover rate were systematically varied. A hypothetical habitat-limited species dispersed between adjacent habitat cells. All trials showed a reduced species frequency relative to a static landscape. Reduction was greater at low habitat density (P = 0.30) than at high density (0.90) suggesting the importance of habitat connectivity in controlling species frequency. The greatest reduction occurred when habitat was concentrated in a small, regularly distributed zone at low habitat density reflecting the enforced isolation of individual habitat cells. Very little reduction was observed when habitat cells were packed into a small clustered zone, a situation promoting connectivity between cells. Moderate–severe frequency reduction occurred when habitat turnover was concentrated in a clustered zone at high habitat density, but little was observed when turnover was widely distributed in a regular or random pattern. These results can be interpreted in terms of a source-sink function in which spatial pattern controlled the degree of contact between landscape zones and determined opportunities for dispersal between habitat cells. We conclude that clustering of forest habitat has the potential to maintain herb species frequency in sparsely forested landscapes. Conversely, clustering of forest disturbance in heavily forested regions, or regular distribution of forest stands (as often occurs in agricultural regions) creates areas which are difficult to colonize, and should be avoided.     

Differential susceptibility epidemic models

We formulate compartmental differential susceptibility (DS) susceptible-infective-removed (SIR) models by dividing the susceptible population into multiple subgroups according to the susceptibility of individuals in each group. We analyze the impact of disease-induced mortality in the situations where the number of contacts per individual is either constant or proportional to the total population. We derive an explicit formula for the reproductive number of infection for each model by investigating the local stability of the infection-free equilibrium. We further prove that the infection-free equilibrium of each model is globally asymptotically stable by qualitative analysis of the dynamics of the model system and by utilizing an appropriately chosen Liapunov function. We show that if the reproductive number is greater than one, then there exists a unique endemic equilibrium for all of the DS models studied in this paper. We prove that the endemic equilibrium is locally asymptotically stable for the models with no disease-induced mortality and the models with contact numbers proportional to the total population. We also provide sufficient conditions for the stability of the endemic equilibrium for other situations. We briefly discuss applications of the DS models to optimal vaccine strategies and the connections between the DS models and predator-prey models with multiple prey populations or host-parasitic interaction models with multiple hosts are also given.This research was partially supported by the Department of Energy under contracts W-7405-ENG-36 and the Applied Mathematical Sciences Program KC-07-01-01.     

Structural and spectroscopic studies of a model for catechol oxidase

A binuclear copper complex, [Cu₂(BPMP)(OAc)₂][ClO₄]·H₂O, has been prepared using the binucleating ligand 2,6-bis[bis(pyridin-2-ylmethylamino)methyl]-4-methylphenol (H-BPMP). The X-ray crystal structure reveals the copper centers to have a five-coordinate square pyramidal geometry, with the acetate ligands bound terminally. The bridging phenolate occupies the apical position of the square-based pyramids and magnetic susceptibility, electron paramagnetic resonance (EPR) and variable-temperature variable-field magnetic circular dichroism (MCD) measurements indicate that the two centers are very weakly antiferromagnetically coupled (J = −0.6 cm⁻¹). Simulation of the dipole–dipole-coupled EPR spectrum showed that in solution the Cu–O–Cu angle was increased from 126° to 160° and that the internuclear distance was larger than that observed crystallographically. The high-resolution spectroscopic information obtained has been correlated with a detailed ligand-field analysis to gain insight into the electronic structure of the complex. Symmetry arguments have been used to demonstrate that the sign of the MCD is characteristic of the tetragonally elongated environment. The complex also displays catecholase activity (k _cat = 15 ± 1.5 min⁻¹, K _M = 6.4 ± 1.8 mM), which is compared with other dicopper catechol oxidase models. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Confidence bands for cumulative incidence curves under the additive risk model

In the context of competing risks, the cumulative incidence function is often used to summarize the cause-specific failure-time data. As an alternative to the proportional hazards model, the additive risk model is used to investigate covariate effects by specifying that the subject-specific hazard function is the sum of a baseline hazard function and a regression function of covariates. Based on such a formulation, we present an approach to constructing simultaneous confidence intervals for the cause-specific cumulative incidence function of patients with given risk factors. A melanoma data set is used for the purpose of illustration.     

Genomic data reveal ancient microendemism in forest scorpions across the California Floristic Province

The California Floristic Province (CFP) in western North America is a globally significant biodiversity hotspot. Elucidating patterns of endemism and the historical drivers of this diversity has been an important challenge of comparative phylogeography for over two decades. We generated phylogenomic data using ddRADseq to examine genetic structure in Uroctonus forest scorpions, an ecologically restricted and dispersal‐limited organism widely distributed across the CFP north to the Columbia River. We coupled our genetic data with species distribution models (SDMs) to determine climatically suitable areas for Uroctonus both now and during the Last Glacial Maximum. Based on our analyses, Uroctonus is composed of two major genetic groups that likely diverged over 2 million years ago. Each of these groups itself contains numerous genetic groups that reveal a pattern of vicariance and microendemism across the CFP. Migration rates among these populations are low. SDMs suggest forest scorpion habitat has remained relatively stable over the last 21 000 years, consistent with the genetic data. Our results suggest tectonic plate rafting, mountain uplift, river drainage formation and climate‐induced habitat fragmentation have all likely played a role in the diversification of Uroctonus. The intricate pattern of genetic fragmentation revealed across a temporal continuum highlights the potential of low‐dispersing species to shed light on small‐scale patterns of biodiversity and the underlying processes that have generated this diversity in biodiversity hotspots.     

The combined effects of energy and disturbance on species richness in protist microcosms

The supply of energy, and the frequency of disturbance can both affect species diversity, often, although not always, producing unimodal diversity patterns along an energy or disturbance gradient. However, it has been suggested that diversity is a combined function of both factors, and the dynamic equilibrium model proposed by Huston (1979) is one formalization of this suggestion. This idea has received little direct testing. Here we carry out such a test using protist microcosms, in a factorial manipulation of six levels of energy (quantity of organic resource) and five levels of disturbance (frequency of temperature shock). Species richness responded to both energy supply and disturbance frequency, although not always unimodally. Patterns of response changed over time, and there was some evidence of an interaction between energy and disturbance at two of six time intervals. However, the specific form of this interaction differs from that predicted by the dynamic equilibrium model, and overall, the results provide little support for the model.     

A High-throughput Method for Measurement of Glomerular Filtration Rate in Conscious Mice

The measurement of glomerular filtration rate (GFR) is the gold standard in kidney function assessment. Currently, investigators determine GFR by measuring the level of the endogenous biomarker creatinine or exogenously applied radioactive labeled inulin (³H or ¹⁴C). Creatinine has the substantial drawback that proximal tubular secretion accounts for ~50% of total renal creatinine excretion and therefore creatinine is not a reliable GFR marker. Depending on the experiment performed, inulin clearance can be determined by an intravenous single bolus injection or continuous infusion (intravenous or osmotic minipump). Both approaches require the collection of plasma or plasma and urine, respectively. Other drawbacks of radioactive labeled inulin include usage of isotopes, time consuming surgical preparation of the animals, and the requirement of a terminal experiment. Here we describe a method which uses a single bolus injection of fluorescein isothiocyanate-(FITC) labeled inulin and the measurement of its fluorescence in 1-2 μl of diluted plasma. By applying a two-compartment model, with 8 blood collections per mouse, it is possible to measure GFR in up to 24 mice per day using a special work-flow protocol. This method only requires brief isoflurane anesthesia with all the blood samples being collected in a non-restrained and awake mouse. Another advantage is that it is possible to follow mice over a period of several months and treatments (i.e. doing paired experiments with dietary changes or drug applications). We hope that this technique of measuring GFR is useful to other investigators studying mouse kidney function and will replace less accurate methods of estimating kidney function, such as plasma creatinine and blood urea nitrogen.     

Genome-scale modeling for Bacillus coagulans to understand the metabolic characteristics

Lactic acid is widely used in many industries, especially in the production of poly-lactic acid. Bacillus coagulans is a promising lactic acid producer in industrial fermentation due to its thermophilic property. In this study, we developed the first genome-scale metabolic model (GEM) of B. coagulans iBag597, together with an enzyme-constrained model ec-iBag597. We measured strain-specific biomass composition and integrated the data into a biomass equation. Then, we validated iBag597 against experimental data generated in this study, including amino acid requirements and carbon source utilization, showing that simulations were generally consistent with the experimental results. Subsequently, we carried out chemostats to investigate the effects of specific growth rate and culture pH on metabolism of B. coagulans. Meanwhile, we used iBag597 to estimate the intracellular metabolic fluxes for those conditions. The results showed that B. coagulans was capable of generating ATP via multiple pathways, and switched among them in response to various conditions. With ec-iBag597, we estimated the protein cost and protein efficiency for each ATP-producing pathway to investigate the switches. Our models pave the way for systems biology of B. coagulans, and our findings suggest that maintaining a proper growth rate and selecting an optimal pH are beneficial for lactate fermentation.     

Lithium Effects on Inositol Phospholipids and Inositol Phosphates: Evaluation of an In Vivo Model for Assessing Polyphosphoinositide Turnover in Brain

Administration of lithium chloride to rats injected intracerebrally with [3H]inositol led to time- and dose-dependent increases in levels of labeled inositol monophosphates in brain. Quantitative analysis of the inositol phosphates by ion chromatography revealed 37- and 20-fold increases in the mass of myo-inositol 1-phosphate and 4-phosphate, respectively, at 4 h intraperitoneal after injections of 6 mEq/kg of lithium chloride. Albeit to a much lesser extent, lithium administration also resulted in an increase in the level of myo-inositol, 1,4-bisphosphate in brain. The lithium-induced increase in content of labeled inositol monophosphates was marked by a concomitant decrease in content of labeled inositol, and after injections of high doses of lithium, e.g., 10 mEq/kg, this was followed by a general decrease in labeling of the inositol phospholipids. In general, animals injected with [3H]inositol but not lithium did not reveal obvious differences in labeling of inositol monophosphates on stimulation by mecamylamine or pilocarpine. However, when animals were injected with [3H]inositol and then lithium, there were large increases in the levels of labeled inositol monophosphates on administration of these compounds. Administration of atropine to the lithium-treated mice led to a partial reduction in the amount of labeled inositol monophosphates accumulated due to the administration of lithium alone. Furthermore, atropine was able to block the pilocarpine-induced increase in level of labeled inositol monophosphates. These results demonstrate the suitable use of the radiotracer technique together with lithium administration for assessing the effects of drugs and receptor agonists on the signaling system involving polyphosphoinositide turnover in brain.