A BAYESIAN HIERARCHICAL MODEL FOR CONSTRAINED DISTRIBUTED LAG FUNCTIONS: ESTIMATING THE TIME COURSE OF HOSPITALIZATION ASSOCIATED WITH AIR POLLUTION EXPOSURE

Numerous time series studies have provided strong evidence of an association between increased levels of ambient air pollution and increased levels of hospital admissions, typically at 0, 1, or 2 days after an air pollution episode. An important research aim is to extend existing statistical models so that a more detailed understanding of the time course of hospitalization after exposure to air pollution can be obtained. Information about this time course, combined with prior knowledge about biological mechanisms, could provide the basis for hypotheses concerning the mechanism by which air pollution causes disease. Previous studies have identified two important methodological questions: (1) How can we estimate the shape of the distributed lag between increased air pollution exposure and increased mortality or morbidity? and (2) How should we estimate the cumulative population health risk from short-term exposure to air pollution? Distributed lag models are appropriate tools for estimating air pollution health effects that may be spread over several days. However, estimation for distributed lag models in air pollution and health applications is hampered by the substantial noise in the data and the inherently weak signal that is the target of investigation. We introduce an hierarchical Bayesian distributed lag model that incorporates prior information about the time course of pollution effects and combines information across multiple locations. The model has a connection to penalized spline smoothing using a special type of penalty matrix. We apply the model to estimating the distributed lag between exposure to particulate matter air pollution and hospitalization for cardiovascular and respiratory disease using data from a large United States air pollution and hospitalization database of Medicare enrollees in 94 counties covering the years 1999-2002.

Enumerating constrained elementary flux vectors of metabolic networks

The concept of elementary vector is generalised to the case where the steady-state space of the metabolic network is not a flux cone but is a general polyhedron due to further inhomogeneous constraints on the flows through some of the reactions. On one hand, this allows to selectively enumerate elementary modes which satisfy certain optimality criteria and this can yield a large computational gain compared with full enumeration. On the other hand, in contrast to the single optimum found by executing a linear program, this enables a comprehensive description of the set of alternate optima often encountered in flux balance analysis. The concepts are illustrated on a metabolic network model of human cardiac mitochondria.

Discordant responses to potent antiretroviral treatment in previously naive HIV-1-infected adults initiating treatment in resource-constrained countries: the antiretroviral therapy in low-income countries (ART-LINC) collaboration

OBJECTIVES: To assess the frequency of and risk factors for discordant responses at 6 months on highly active antiretroviral therapy (HAART) in previously treatment-naive HIV patients from resource-limited countries. METHODS: The Antiretroviral Therapy in Low-Income Countries Collaboration is a network of clinics providing care and treatment to HIV-infected patients in Africa, Latin America, and Asia. Patients who initiated therapy between 1996 and 2004, were aged 16 years or older, and had a baseline CD4 cell count were included in this analysis. Responses were defined based on plasma viral load (PVL) and CD4 cell count at 6 months as complete virologic and immunologic (VR(+)IR(+)), virologic only (VR(+)IR(-)), immunologic only (VR(-)IR(+)), and nonresponse (VR(-)IR(-)). Multinomial logistic regression was used to assess the association between therapy responses and clinical and demographic variables. RESULTS: Of the 3111 patients eligible for analysis, 1914 had available information at 6 months of therapy: 1074 (56.1%) were VR(+)IR(+), 364 (19.0%) were VR(+)IR(-), 283 (14.8%) were (VR(-)IR(+)), and 193 (10.1%) were VR(-)IR(-). OF THE 3111 patients eligible for analysis, 1914 had available information at 6 months of therapy: 1074 (56.1%) were VRIR, 364 (19.0%) were VRIR, 283 (14.8%) were (VRIR), and 193 (10.1%) were VRIR. Compared with complete responders, virologic-only responders were older, had a higher baseline CD4 cell count, had a lower baseline PVL, and were more likely to have received a nonstandard HAART regimen; immunologic-only responders were younger, had a lower baseline CD4 cell count, had a higher baseline PVL, and were more likely to have received a protease inhibitor-based regimen. CONCLUSIONS: The frequency of and risk factors for discordant responses were comparable to those observed in developed countries. Longer follow-up is needed to assess the long-term impact of discordant responses on mortality in these resource-limited settings.

New sst4/5-selective somatostatin peptidomimetics based on a constrained tryptophan scaffold

The synthesis and biological evaluation of four peptidomimetic analogs of somatostatin based on a constrained Trp residue, 3-amino-indolo[2,3-c]azepin-2-one (Aia), are reported. It is shown that dipeptidomimetics with a D-Aia-Lys sequence, functionalized with N- and C-terminal aromatic substituents, display a good selectivity for both sst4 and sst5. This study allowed us to identify a new highly potent sst5 agonist with good selectivity over the other receptors, except versus sst4.

The effect of matrix bound parathyroid hormone on bone regeneration

INTRODUCTION: Autogenous bone is the most successful bone-grafting material; however, multiple disadvantages continue to drive developments of improved methods for bone regeneration. AIM: The aim of the present study was to test the hypothesis that an arginine-glycine-aspartic acid (RGD) modified polyethylene glycol-based matrix (PEG) containing covalently bound peptides of the parathyroid hormone (PTH(1-34)) enhances bone regeneration to a degree similar to autogenous bone. MATERIAL AND METHODS: Six American foxhounds received a total of 48 cylindrical titanium implants placed in the mandible between the first premolar and the second molar. Five, respectively, 7 months following tooth extraction, implants were placed into the center of surgically created defects. This resulted in a circumferential bone defect simulating an alveolar defect with a circular gap of 1.5 mm. Four treatment modalities were randomly allocated to the four defects per side: (1) PEG-matrix containing 20 microg/ml of PTH(1-34), and 350 microg/ml cys-RGD peptide, (2) PEG alone, (3) autogenous bone and (4) empty defects. Histomorphometric analysis was performed 4 and 12 weeks after implantation. The area fraction of newly formed bone was determined within the former defect and the degree of bone-to-implant contact (BIC) was evaluated both in the defect region and in the apical region of the implant. For statistical analysis ANOVA and subsequent pairwise Student's t-test were applied. RESULTS: Healing was uneventful and all implants were histologically integrated. Histomorphometric analysis after 4 weeks showed an average area fraction of newly formed bone of 41.7+/-1.8% for matrix-PTH, 26.6+/-4.1% for PEG alone, 43.9+/-4.5% for autogenous bone, and 28.9+/-1.5% for empty defects. After 12 weeks, the respective values were 49.4+/-7.0% for matrix-PTH, 39.3+/-5.7% for PEG alone, 50.5+/-3.4% for autogenous bone and 38.7+/-1.9% for empty defects. Statistical analysis after 4 and 12 weeks revealed significantly more newly formed bone in the PTH(1-34) group compared with PEG alone or empty defects, whereas no difference could be detected against autogenous bone. Regarding BIC no significant difference was observed between the four treatment groups neither at 4 nor at 12 weeks. CONCLUSION: It is concluded that an RGD-modified PEG hydrogel containing PTH(1-34) is an effective matrix system to obtain bone regeneration.

Factor XIII activation peptide is released into plasma upon cleavage by thrombin and shows a different structure compared to its bound form

The first step of coagulation factor XIII (FXIII) activation involves cleavage of the FXIII activation peptide (FXIII-AP) by thrombin. However, it is not known whether the FXIII-AP is released into plasma upon cleavage or remains attached to activated FXIII. The aim of the present work was to study the structure of free FXIII-AP, develop an assay for FXIII-AP determination in human plasma, and to answer the question whether FXIII-AP is released into plasma. We used ab-initio modeling and molecular dynamics simulations to study the structure of free FXIII-AP. We raised monoclonal and polyclonal antibodies against FXIII-AP and developed a highly sensitive and specific ELISA method for direct detection of FXIII-AP in human plasma. Structural analysis showed a putative different conformation of the free FXIII-AP compared to FXIII-AP bound to the FXIII protein. We concluded that it might be feasible to develop specific antibodies against the free FXIII-AP. Using our new FXIII-AP ELISA, we found high levels of FXIII-AP in in-vitro activated plasma samples and serum. We showed for the first time that FXIIIAP is detached from activated FXIII and is released into plasma, where it can be directly measured. Our findings may be of major clinical interest in regard to a possible new marker in thrombotic disease.

Phosphatidylethanolamine is the precursor of the ethanolamine phosphoglycerol moiety bound to eukaryotic elongation factor 1A

In addition to its conventional role during protein synthesis, eukaryotic elongation factor 1A is involved in other cellular processes. Several regions of interaction between eukaryotic elongation factor 1A and the translational apparatus or the cytoskeleton have been identified, yet the roles of the different post-translational modifications of eukaryotic elongation factor 1A are completely unknown. One amino acid modification, which so far has only been found in eukaryotic elongation factor 1A, consists of ethanolamine-phosphoglycerol attached to two glutamate residues that are conserved between mammals and plants. We now report that ethanolamine-phosphoglycerol is also present in eukaryotic elongation factor 1A of the protozoan parasite Trypanosoma brucei, indicating that this unique protein modification is of ancient origin. In addition, using RNA-mediated gene silencing against enzymes of the Kennedy pathway, we demonstrate that phosphatidylethanolamine is a direct precursor of the ethanolamine-phosphoglycerol moiety. Down-regulation of the expression of ethanolamine kinase and ethanolamine-phosphate cytidylyltransferase results in inhibition of phosphatidylethanolamine synthesis in T. brucei procyclic forms and, concomitantly, in a block in glycosylphosphatidylinositol attachment to procyclins and ethanolamine-phosphoglycerol modification of eukaryotic elongation factor 1A.

Meprins, membrane-bound and secreted astacin metalloproteinases

The astacins are a subfamily of the metzincin superfamily of metalloproteinases. The first to be characterized was the crayfish enzyme astacin. To date more than 200 members of this family have been identified in species ranging from bacteria to humans. Astacins are involved in developmental morphogenesis, matrix assembly, tissue differentiation and digestion. Family members include the procollagen C-proteinase (BMP1, bone morphogenetic protein 1), tolloid and mammalian tolloid-like, HMP (Hydra vulgaris metalloproteinase), sea urchin BP10 (blastula protein) and SPAN (Strongylocentrotus purpuratus astacin), the 'hatching' subfamily comprising alveolin, ovastacin, LCE, HCE ('low' and 'high' choriolytic enzymes), nephrosin (from carp head kidney), UVS.2 from frog, and the meprins. In the human and mouse genomes, there are six astacin family genes (two meprins, three BMP1/tolloid-like, one ovastacin), but in Caenorhabditis elegans there are 40. Meprins are the only astacin proteinases that function on the membrane and extracellularly by virtue of the fact that they can be membrane-bound or secreted. They are unique in their domain structure and covalent subunit dimerization, oligomerization propensities, and expression patterns. They are normally highly regulated at the transcriptional and post-translational levels, localize to specific membranes or extracellular spaces, and can hydrolyse biologically active peptides, cytokines, extracellular matrix (ECM) proteins and cell-surface proteins. The in vivo substrates of meprins are unknown, but the abundant expression of these proteinases in the epithelial cells of the intestine, kidney and skin provide clues to their functions.

Parallel algorithm for solving integer linear programs

Linear programs, or LPs, are often used in optimization problems, such as improving manufacturing efficiency of maximizing the yield from limited resources. The most common method for solving LPs is the Simplex Method, which will yield a solution, if one exists, but over the real numbers. From a purely numerical standpoint, it will be an optimal solution, but quite often we desire an optimal integer solution. A linear program in which the variables are also constrained to be integers is called an integer linear program or ILP. It is the focus of this report to present a parallel algorithm for solving ILPs. We discuss a serial algorithm using a breadth-first branch-and-bound search to check the feasible solution space, and then extend it into a parallel algorithm using a client-server model. In the parallel mode, the search may not be truly breadth-first, depending on the solution time for each node in the solution tree. Our search takes advantage of pruning, often resulting in super-linear improvements in solution time. Finally, we present results from sample ILPs, describe a few modifications to enhance the algorithm and improve solution time, and offer suggestions for future work.

Vibro-acoustical analysis and design of a multiple-layer constrained viscoelastic damping structure

The goal of this research is to provide a framework for vibro-acoustical analysis and design of a multiple-layer constrained damping structure. The existing research on damping and viscoelastic damping mechanism is limited to the following four mainstream approaches: modeling techniques of damping treatments/materials; control through the electrical-mechanical effect using the piezoelectric layer; optimization by adjusting the parameters of the structure to meet the design requirements; and identification of the damping material’s properties through the response of the structure. This research proposes a systematic design methodology for the multiple-layer constrained damping beam giving consideration to vibro-acoustics. A modeling technique to study the vibro-acoustics of multiple-layered viscoelastic laminated beams using the Biot damping model is presented using a hybrid numerical model. The boundary element method (BEM) is used to model the acoustical cavity whereas the Finite Element Method (FEM) is the basis for vibration analysis of the multiple-layered beam structure. Through the proposed procedure, the analysis can easily be extended to other complex geometry with arbitrary boundary conditions. The nonlinear behavior of viscoelastic damping materials is represented by the Biot damping model taking into account the effects of frequency, temperature and different damping materials for individual layers. A curve-fitting procedure used to obtain the Biot constants for different damping materials for each temperature is explained. The results from structural vibration analysis for selected beams agree with published closed-form results and results for the radiated noise for a sample beam structure obtained using a commercial BEM software is compared with the acoustical results of the same beam with using the Biot damping model. The extension of the Biot damping model is demonstrated to study MDOF (Multiple Degrees of Freedom) dynamics equations of a discrete system in order to introduce different types of viscoelastic damping materials. The mechanical properties of viscoelastic damping materials such as shear modulus and loss factor change with respect to different ambient temperatures and frequencies. The application of multiple-layer treatment increases the damping characteristic of the structure significantly and thus helps to attenuate the vibration and noise for a broad range of frequency and temperature. The main contributions of this dissertation include the following three major tasks: 1) Study of the viscoelastic damping mechanism and the dynamics equation of a multilayer damped system incorporating the Biot damping model. 2) Building the Finite Element Method (FEM) model of the multiple-layer constrained viscoelastic damping beam and conducting the vibration analysis. 3) Extending the vibration problem to the Boundary Element Method (BEM) based acoustical problem and comparing the results with commercial simulation software.

Gender and the use of antiretroviral treatment in resource-constrained settings: findings from a multicenter collaboration

AIMS: To compare the gender distribution of HIV-infected adults receiving highly active antiretroviral treatment (HAART) in resource-constrained settings with estimates of the gender distribution of HIV infection; to describe the clinical characteristics of women and men receiving HAART. METHODS: The Antiretroviral Therapy in Lower-Income Countries, ART-LINC Collaboration is a network of clinics providing HAART in Africa, Latin America, and Asia. We compared UNAIDS data on the gender distribution of HIV infection with the proportions of women and men receiving HAART in the ART-LINC Collaboration. RESULTS: Twenty-nine centers in 13 countries participated. Among 33,164 individuals, 19,989 (60.3%) were women. Proportions of women receiving HAART in ART-LINC centers were similar to, or higher than, UNAIDS estimates of the proportions of HIV-infected women in all but two centers. There were fewer women receiving HAART than expected from UNAIDS data in one center in Uganda and one center in India. Taking into account heterogeneity across cohorts, women were younger than men, less likely to have advanced HIV infection, and more likely to be anemic at HAART initiation. CONCLUSIONS: Women in resource-constrained settings are not necessarily disadvantaged in their access to HAART. More attention needs to be paid to ensuring that HIV-infected men are seeking care and starting HAART.

A REGIONAL ASSESSMENT OF GROUNDWATER AVAILABILITY AS CONSTRAINED BY LOCAL HYDROGEOLOGY AND ENVIRONMENTAL LIMITS TO STREAMFLOW DEPLETION

Groundwater pumping from aquifers in hydraulic connection with nearby streams is known to cause adverse impacts by decreasing flows to levels below those necessary to maintain aquatic ecosystems. The recent passage of the Great Lakes--St. Lawrence River Basin Water Resources Compact has brought attention to this issue in the Great Lakes region. In particular, the legislation requires the Great Lakes states to enact measures for limiting water withdrawals that can cause adverse ecosystem impacts. This study explores how both hydrogeologic and environmental flow limitations constrain groundwater availability in the Great Lakes Basin. A methodology for calculating maximum allowable pumping rates is presented. Groundwater availability across the basin is shown to be constrained by a combination of hydrogeologic yield and environmental flow limitations varying over both local and regional scales. The results are sensitive to factors such as pumping time and streamflow depletion limits as well as streambed conductance. Understanding how these restrictions constrain groundwater usage and which hydrogeologic characteristics and spatial variables have the most influence on potential streamflow depletions has important water resources policy and management implications.