The impact of gastrointestinal dysmotility on the aerodigestive microbiome of pediatric lung transplant recipients

BACKGROUND

Delayed gastric emptying has been associated with increased graft rejection, although the mechanism of this association is not known. This study aims to investigate the interrelationship between delays in gastrointestinal motility and the diversity and composition of gastric, oropharyngeal, and lung microbiomes in pediatric lung transplant recipients.

METHODS

We prospectively recruited 23 pediatric lung transplant recipients and 98 pediatric patients with respiratory symptoms undergoing combined endoscopy and bronchoscopy. Gastric, oropharyngeal, and bronchoalveolar lavage samples were collected for 16S sequencing. Gastric samples were also analyzed for bile composition using liquid chromatography.

RESULTS

Patients who underwent lung transplantation had significantly reduced alpha diversity in gastric and oropharyngeal sites compared with patients with respiratory symptoms. This reduction in alpha diversity was especially evident in gastric samples in patients with delayed gastric emptying defined as abnormal gastric emptying on nuclear scintigraphy or as an elevation in gastric bile concentration (p ≤ 0.05). Whereas monocolonies were seen in the lungs of patients who underwent transplantation, these were not the same microbes seen in the stomach; the microbial overlap between lung and gastric samples within patients was low, and data indicated high individual variation between lung transplant recipients. Other contributors to reduced alpha diversity included antibiotics in combination with proton pump inhibitors, especially in gastric and oropharyngeal samples.

CONCLUSIONS

Lung transplant recipients have reduced microbial diversity in gastric fluid (GF) and oropharynx compared with patients who did not undergo lung transplantation. The decreased alpha diversity in GF may be associated with dysmotility.

Progressive Shifts in the Gut Microbiome Reflect Prediabetes and Diabetes Development in a Treatment-Naive Mexican Cohort

Type 2 diabetes (T2D) is a global epidemic that affects more than 8% of the world's population and is a leading cause of death in Mexico. Diet and lifestyle are known to contribute to the onset of T2D. However, the role of the gut microbiome in T2D progression remains uncertain. Associations between microbiome composition and diabetes are confounded by medication use, diet, and obesity. Here we present data on a treatment-naive cohort of 405 Mexican individuals across varying stages of T2D severity. Associations between gut bacteria and more than 200 clinical variables revealed a defined set of bacterial genera that were consistent biomarkers of T2D prevalence and risk. Specifically, gradual increases in blood glucose levels, beta cell dysfunction, and the accumulation of measured T2D risk factors were correlated with the relative abundances of four bacterial genera. In a cohort of 25 individuals, T2D treatment-predominantly metformin-reliably returned the microbiome to the normoglycemic community state. Deep clinical characterization allowed us to broadly control for confounding variables, indicating that these microbiome patterns were independent of common T2D comorbidities, like obesity or cardiovascular disease. Our work provides the first solid evidence for a direct link between the gut microbiome and T2D in a critically high-risk population. In particular, we show that increased T2D risk is reflected in gradual changes in the gut microbiome. Whether or not these T2D-associated changes in the gut contribute to the etiology of T2D or its comorbidities remains to be seen.

Oil Hydrocarbon Degradation by Caspian Sea Microbial Communities

The Caspian Sea, which is the largest landlocked body of water on the planet, receives substantial annual hydrocarbon input from anthropogenic sources (e.g., industry, agriculture, oil exploration, and extraction) and natural sources (e.g., mud volcanoes and oil seeps). The Caspian Sea also receives substantial amounts of runoff from agricultural and municipal sources, containing nutrients that have caused eutrophication and subsequent hypoxia in the deep, cold waters. The effect of decreasing oxygen saturation and cold temperatures on oil hydrocarbon biodegradation by a microbial community is not well characterized. The purpose of this study was to investigate the effect of oxic and anoxic conditions on oil hydrocarbon biodegradation at cold temperatures by microbial communities derived from the Caspian Sea. Water samples were collected from the Caspian Sea for study in experimental microcosms. Major taxonomic orders observed in the ambient water samples included Flavobacteriales, Actinomycetales, and Oceanospirillales. Microcosms were inoculated with microbial communities from the deepest waters and amended with oil hydrocarbons for 17 days. Hydrocarbon degradation and shifts in microbial community structure were measured. Surprisingly, oil hydrocarbon biodegradation under anoxic conditions exceeded that under oxic conditions; this was particularly evident in the degradation of aromatic hydrocarbons. Important microbial taxa associated with the anoxic microcosms included known oil degraders such as Oceanospirillaceae. This study provides knowledge about the ambient community structure of the Caspian Sea, which serves as an important reference point for future studies. Furthermore, this may be the first report in which anaerobic biodegradation of oil hydrocarbons exceeds aerobic biodegradation.

Pasture flock chicken cecal microbiome responses to prebiotics and plum fiber feed amendments

When prebiotics and other fermentation substrates are delivered to animals as feed supplements, the typical goal is to improve weight gain and feed conversion. In this work, we examined pasture flock chicken cecal contents using next generation sequencing (NGS) to identify and understand the composition of the microbiome when prebiotics and fermentation substrates were supplemented. We generated 16S rRNA sequencing data for 120 separate cecal samples from groups of chickens receiving one of 3 prebiotics or fiber feed additives. The data indicated that respective feed additives enrich for specific bacterial community members and modulate the diversity of the microbiome. We applied synthetic learning in microbial ecology (SLiME) analysis to interpret 16S rRNA microbial community data and identify specific bacterial operational taxonomic units (OTU) that are predictive of the particular feed additives used in these experiments. The results suggest that feed can influence microbiome composition in a predictable way, and thus diet may have indirect effects on weight gain and feed conversion through the microbiome.

Fecal microbiota transplant from a rational stool donor improves hepatic encephalopathy: A randomized clinical trial

Recurrent hepatic encephalopathy (HE) is a leading cause of readmission despite standard of care (SOC) associated with microbial dysbiosis. Fecal microbiota transplantation (FMT) may improve dysbiosis; however, it has not been studied in HE. We aimed to define whether FMT using a rationally derived stool donor is safe in recurrent HE compared to SOC alone. An open-label, randomized clinical trial with a 5-month follow-up in outpatient men with cirrhosis with recurrent HE on SOC was conducted with 1:1 randomization. FMT-randomized patients received 5 days of broad-spectrum antibiotic pretreatment, then a single FMT enema from the same donor with the optimal microbiota deficient in HE. Follow-up occurred on days 5, 6, 12, 35, and 150 postrandomization. The primary outcome was safety of FMT compared to SOC using FMT-related serious adverse events (SAEs). Secondary outcomes were adverse events, cognition, microbiota, and metabolomic changes. Participants in both arms were similar on all baseline criteria and were followed until study end. FMT with antibiotic pretreatment was well tolerated. Eight (80%) SOC participants had a total of 11 SAEs compared to 2 (20%) FMT participants with SAEs (both FMT unrelated; P = 0.02). Five SOC and no FMT participants developed further HE (P = 0.03). Cognition improved in the FMT, but not the SOC, group. Model for End-Stage Liver Disease (MELD) score transiently worsened postantibiotics, but reverted to baseline post-FMT. Postantibiotics, beneficial taxa, and microbial diversity reduction occurred with Proteobacteria expansion. However, FMT increased diversity and beneficial taxa. SOC microbiota and MELD score remained similar throughout.

CONCLUSION

FMT from a rationally selected donor reduced hospitalizations, improved cognition, and dysbiosis in cirrhosis with recurrent HE. (Hepatology 2017;66:1727-1738).

Legionellosis acquired through a dental unit: a case study

In 2012, an elderly immunocompromised man died from legionellosis at a hospital in Uppsala, Sweden. The patient had visited a dental ward at the hospital during the incubation period. Legionella spp. at a concentration of 2000 colony-forming units/L were isolated from the cupfiller outlet providing water for oral rinsing. Isolates from the patient and the dental unit were Legionella pneumophila serogroup 1, subgroup Knoxville and ST9. Pulsed-field gel electrophoresis and whole-genome sequencing strongly suggested that the isolates were of common origin. This report presents one of few documented cases of legionellosis acquired through a dental unit.

The Microbial Olympics 2016

Following the success of the inaugural games, the Microbial Olympics return with a new series of events and microbial competitors. The games may have moved to a new hosting venue, but the dedication to training, fitness, competition (and yes, education and humour) lives on.

Recurrent Clostridium difficile infection associates with distinct bile acid and microbiome profiles

BACKGROUND

The healthy microbiome protects against the development of Clostridium difficile infection (CDI), which typically develops following antibiotics. The microbiome metabolises primary to secondary bile acids, a process if disrupted by antibiotics, may be critical for the initiation of CDI.

AIM

To assess the levels of primary and secondary bile acids associated with CDI and associated microbial changes.

METHODS

Stool and serum were collected from patients with (i) first CDI (fCDI), (ii) recurrent CDI (rCDI) and (iii) healthy controls. 16S rRNA sequencing and bile salt metabolomics were performed. Random forest regression models were constructed to predict disease status. PICRUSt analyses were used to test for associations between predicted bacterial bile salt hydrolase (BSH) gene abundances and bile acid levels.

RESULTS

Sixty patients (20 fCDI, 19 rCDI and 21 controls) were enrolled. Secondary bile acids in stool were significantly elevated in controls compared to rCDI and fCDI (P < 0.0001 and P = 0.0007 respectively). Primary bile acids in stool were significantly elevated in rCDI compared to controls (P < 0.0001) and in rCDI compared to fCDI (P = 0.02). Using random forest regression, we distinguished rCDI and fCDI patients 84.2% of the time using bile acid ratios. Stool deoxycholate to glycoursodeoxycholate ratio was the single best predictor. PICRUSt analyses found significant differences in predicted abundances of bacterial BSH genes in stool samples across the groups.

CONCLUSIONS

Primary and secondary bile acid composition in stool was different in those with rCDI, fCDI and controls. The ratio of stool deoxycholate to glycoursodeoxycholate was the single best predictor of disease state and may be a potential biomarker for recurrence.

"Every Gene Is Everywhere but the Environment Selects": Global Geolocalization of Gene Sharing in Environmental Samples through Network Analysis

The spatial distribution of microbes on our planet is famously formulated in the Baas Becking hypothesis as “everything is everywhere but the environment selects.” While this hypothesis does not strictly rule out patterns caused by geographical effects on ecology and historical founder effects, it does propose that the remarkable dispersal potential of microbes leads to distributions generally shaped by environmental factors rather than geographical distance. By constructing sequence similarity networks from uncultured environmental samples, we show that microbial gene pool distributions are not influenced nearly as much by geography as ecology, thus extending the Bass Becking hypothesis from whole organisms to microbial genes. We find that gene pools are shaped by their broad ecological niche (such as sea water, fresh water, host, and airborne). We find that freshwater habitats act as a gene exchange bridge between otherwise disconnected habitats. Finally, certain antibiotic resistance genes deviate from the general trend of habitat specificity by exhibiting a high degree of cross-habitat mobility. The strong cross-habitat mobility of antibiotic resistance genes is a cause for concern and provides a paradigmatic example of the rate by which genes colonize new habitats when new selective forces emerge.

Detection of low-abundance bacterial strains in metagenomic datasets by eigengenome partitioning

Analyses of metagenomic datasets that are sequenced to a depth of billions or trillions of bases can uncover hundreds of microbial genomes, but naive assembly of these data is computationally intensive, requiring hundreds of gigabytes to terabytes of RAM. We present latent strain analysis (LSA), a scalable, de novo pre-assembly method that separates reads into biologically informed partitions and thereby enables assembly of individual genomes. LSA is implemented with a streaming calculation of unobserved variables that we call eigengenomes. Eigengenomes reflect covariance in the abundance of short, fixed-length sequences, or k-mers. As the abundance of each genome in a sample is reflected in the abundance of each k-mer in that genome, eigengenome analysis can be used to partition reads from different genomes. This partitioning can be done in fixed memory using tens of gigabytes of RAM, which makes assembly and downstream analyses of terabytes of data feasible on commodity hardware. Using LSA, we assemble partial and near-complete genomes of bacterial taxa present at relative abundances as low as 0.00001%. We also show that LSA is sensitive enough to separate reads from several strains of the same species.

Abstract 321: A High-Salt Diet Alters the Composition of Intestinal Microbiota in Mice

Objective: High-salt intake is associated with cardiovascular disease and hypertension. A substantial number of patients are sensitive to salt, but environmental factors influencing this phenomenon are poorly understood. The gut microbiome has not been considered in the context of salt and hypertension research to date. It is known to respond to fluctuations in lifestyle and diet, and is increasingly recognized as an important factor which influences host health and disease. Using next generation sequencing methods we aimed to determine the effect of a high-salt diet on gut microbiome composition.

Methods and Results: Male C57BL6/J mice, 10 weeks of age, were individually housed and either subjected to a purified normal salt (0.5% sodium) or high-salt diet (4% sodium + 1% in drinking water) ad libitum for 14 days. Feces was collected on days 0 and 14 of the diets. DNA was extracted from feces, the V4 region of the 16S rRNA gene amplified with PCR. Samples were multiplexed and sequenced with MiSeq. Reads were quality controlled and operational taxonomic units (OTUs) were called by comparison to the Greengenes database.

The high-salt diet created a distinct gut microbiome composition compared to the normal-salt diet (analysis of Jensen-Shannon divergence). Whereas normal diet’s gut microbiome is composed mostly of OTUs from the Bacteroidetes phylum (70% of reads), high-salt diet causes a shift in relative abundance to other phyla (increased Firmicutes:Bacteroidetes ratio). Moreover, the genus Allobaculum significantly spikes after high-salt feeding. Both phenomena have recently been implicated in obesity studies.

Conclusion: An increase in dietary sodium chloride alters the gut microbiome composition. Our findings indicate new perspectives on how salt impacts the body.

Computational methods for high-throughput comparative analyses of natural microbial communities

One of the most widely employed methods in metagenomics is the amplification and sequencing of the highly conserved ribosomal RNA (rRNA) genes from organisms in complex microbial communities. rRNA surveys, typically using the 16S rRNA gene for prokaryotic identification, provide information about the total diversity and taxonomic affiliation of organisms present in a sample. Greatly enhanced by high-throughput sequencing, these surveys have uncovered the remarkable diversity of uncultured organisms and revealed unappreciated ecological roles ranging from nutrient cycling to human health. This chapter outlines the best practices for comparative analyses of microbial community surveys. We explain how to transform raw data into meaningful units for further analysis and discuss how to calculate sample diversity and community distance metrics. Finally, we outline how to find associations of species with specific metadata and true correlations between species from compositional data. We focus on data generated by next-generation sequencing platforms, using the Illumina platform as a test case, because of its widespread use especially among researchers just entering the field.

The slow:fast substitution ratio reveals changing patterns of natural selection in gamma-proteobacterial genomes

Different microbial species are thought to occupy distinct ecological niches, subjecting each species to unique selective constraints, which may leave a recognizable signal in their genomes. Thus, it may be possible to extract insight into the genetic basis of ecological differences among lineages by identifying unusual patterns of substitutions in orthologous gene or protein sequences. We used the ratio of substitutions in slow- versus fast-evolving sites (nucleotides in DNA, or amino acids in protein sequence) to quantify deviations from the typical pattern of selective constraint observed across bacterial lineages. We propose that elevated S:F in one branch (an excess of slow-site substitutions) can indicate a functionally relevant change, due to either positive selection or relaxed evolutionary constraint. In a genome-wide comparative study of gamma-proteobacterial proteins, we find that cell-surface proteins involved with motility and secretion functions often have high S:F ratios, whereas information-processing genes do not. Change in evolutionary constraints in some species is evidenced by increased S:F ratios within functionally related sets of genes (for example, energy production in Pseudomonas fluorescens), whereas other species apparently evolve mostly by drift (for example, uniformly elevated S:F across most genes in Buchnera spp.). Overall, S:F reveals several species-specific, protein-level changes with potential functional/ecological importance. As microbial genome projects yield more species-rich gene trees, the S:F ratio will become an increasingly powerful tool for uncovering functional genetic differences among species.

The evolution of two-component systems in bacteria reveals different strategies for niche adaptation

Two-component systems including histidine protein kinases represent the primary signal transduction paradigm in prokaryotic organisms. To understand how these systems adapt to allow organisms to detect niche-specific signals, we analyzed the phylogenetic distribution of nearly 5,000 histidine protein kinases from 207 sequenced prokaryotic genomes. We found that many genomes carry a large repertoire of recently evolved signaling genes, which may reflect selective pressure to adapt to new environmental conditions. Both lineage-specific gene family expansion and horizontal gene transfer play major roles in the introduction of new histidine kinases into genomes; however, there are differences in how these two evolutionary forces act. Genes imported via horizontal transfer are more likely to retain their original functionality as inferred from a similar complement of signaling domains, while gene family expansion accompanied by domain shuffling appears to be a major source of novel genetic diversity. Family expansion is the dominant source of new histidine kinase genes in the genomes most enriched in signaling proteins, and detailed analysis reveals that divergence in domain structure and changes in expression patterns are hallmarks of recent expansions. Finally, while these two modes of gene acquisition are widespread across bacterial taxa, there are clear species-specific preferences for which mode is used.

Pathways containing histidine protein kinases (HPKs) represent a key mechanism for signal transduction, especially in bacteria. These systems help cells to sense and respond to their environment by detecting external cues and effecting internal responses such as changes in gene expression. As such, they are believed to play a key role in niche adaptation, yet their evolution is difficult to study due to the large number of paralogous subfamilies. This work extends previous large-scale gene evolution studies by considering complex paralogy relationships, and uncovers an abundance of horizontal transfers, gene duplications, and domain shuffling that have marked the evolutionary history of HPKs. An important finding of this study is qualitative differences between the main strategies for acquiring new HPKs (horizontal gene transfer and gene duplication). Hallmarks of the latter process include domain shuffling and the generation of “orphan” HPKs not co-transcribed with a cognate response regulator.

Transcriptome profiling of Shewanella oneidensis gene expression following exposure to acidic and alkaline pH

The molecular response of Shewanella oneidensis MR-1 to variations in extracellular pH was investigated based on genomewide gene expression profiling. Microarray analysis revealed that cells elicited both general and specific transcriptome responses when challenged with environmental acid (pH 4) or base (pH 10) conditions over a 60-min period. Global responses included the differential expression of genes functionally linked to amino acid metabolism, transcriptional regulation and signal transduction, transport, cell membrane structure, and oxidative stress protection. Response to acid stress included the elevated expression of genes encoding glycogen biosynthetic enzymes, phosphate transporters, and the RNA polymerase sigma-38 factor (rpoS), whereas the molecular response to alkaline pH was characterized by upregulation of nhaA and nhaR, which are predicted to encode an Na+/H+ antiporter and transcriptional activator, respectively, as well as sulfate transport and sulfur metabolism genes. Collectively, these results suggest that S. oneidensis modulates multiple transporters, cell envelope components, and pathways of amino acid consumption and central intermediary metabolism as part of its transcriptome response to changing external pH conditions.

Cellular response of Shewanella oneidensis to strontium stress

The physiology and transcriptome dynamics of the metal ion-reducing bacterium Shewanella oneidensis strain MR-1 in response to nonradioactive strontium (Sr) exposure were investigated. Studies indicated that MR-1 was able to grow aerobically in complex medium in the presence of 180 mM SrCl2 but showed severe growth inhibition at levels above that concentration. Temporal gene expression profiles were generated from aerobically grown, mid-exponential-phase MR-1 cells shocked with 180 mM SrCl2 and analyzed for significant differences in mRNA abundance with reference to data for nonstressed MR-1 cells. Genes with annotated functions in siderophore biosynthesis and iron transport were among the most highly induced (>100-fold [P < 0.05]) open reading frames in response to acute Sr stress, and a mutant (SO3032::pKNOCK) defective in siderophore production was found to be hypersensitive to SrCl2 exposure, compared to parental and wild-type strains. Transcripts encoding multidrug and heavy metal efflux pumps, proteins involved in osmotic adaptation, sulfate ABC transporters, and assimilative sulfur metabolism enzymes also were differentially expressed following Sr exposure but at levels that were several orders of magnitude lower than those for iron transport genes. Precipitate formation was observed during aerobic growth of MR-1 in broth cultures amended with 50, 100, or 150 mM SrCl2 but not in cultures of the SO3032::pKNOCK mutant or in the abiotic control. Chemical analysis of this precipitate using laser-induced breakdown spectroscopy and static secondary ion mass spectrometry indicated extracellular solid-phase sequestration of Sr, with at least a portion of the heavy metal associated with carbonate phases.

Transcriptome analysis of Shewanella oneidensis MR-1 in response to elevated salt conditions

Whole-genomic expression patterns were examined in Shewanella oneidensis cells exposed to elevated sodium chloride. Genes involved in Na(+) extrusion and glutamate biosynthesis were significantly up-regulated, and the majority of chemotaxis/motility-related genes were significantly down-regulated. The data also suggested an important role for metabolic adjustment in salt stress adaptation in S. oneidensis.

Global transcriptome analysis of the heat shock response of Shewanella oneidensis

Shewanella oneidensis is an important model organism for bioremediation studies because of its diverse respiratory capabilities. However, the genetic basis and regulatory mechanisms underlying the ability of S. oneidensis to survive and adapt to various environmentally relevant stresses is poorly understood. To define this organism's molecular response to elevated growth temperatures, temporal gene expression profiles were examined in cells subjected to heat stress by using whole-genome DNA microarrays for S. oneidensis. Approximately 15% (n = 711) of the total predicted S. oneidensis genes (n = 4,648) represented on the microarray were significantly up- or downregulated (P < 0.05) over a 25-min period after shift to the heat shock temperature. As expected, the majority of the genes that showed homology to known chaperones and heat shock proteins in other organisms were highly induced. In addition, a number of predicted genes, including those encoding enzymes in glycolysis and the pentose cycle, serine proteases, transcriptional regulators (MerR, LysR, and TetR families), histidine kinases, and hypothetical proteins were induced. Genes encoding membrane proteins were differentially expressed, suggesting that cells possibly alter their membrane composition or structure in response to variations in growth temperature. A substantial number of the genes encoding ribosomal proteins displayed downregulated coexpression patterns in response to heat stress, as did genes encoding prophage and flagellar proteins. Finally, a putative regulatory site with high conservation to the Escherichia coli sigma32-binding consensus sequence was identified upstream of a number of heat-inducible genes.

Reconstruction of regulatory and metabolic pathways in metal-reducing delta-proteobacteria

A study of the genetic and regulatory factors in several biosynthesis, metal ion homeostasis, stress response, and energy metabolism pathways suggests that phylogenetically diverse δ-proteobacteria have homologous regulatory components.

Background

Relatively little is known about the genetic basis for the unique physiology of metal-reducing genera in the delta subgroup of the proteobacteria. The recent availability of complete finished or draft-quality genome sequences for seven representatives allowed us to investigate the genetic and regulatory factors in a number of key pathways involved in the biosynthesis of building blocks and cofactors, metal-ion homeostasis, stress response, and energy metabolism using a combination of regulatory sequence detection and analysis of genomic context.

Results

In the genomes of δ-proteobacteria, we identified candidate binding sites for four regulators of known specificity (BirA, CooA, HrcA, sigma-32), four types of metabolite-binding riboswitches (RFN-, THI-, B12-elements and S-box), and new binding sites for the FUR, ModE, NikR, PerR, and ZUR transcription factors, as well as for the previously uncharacterized factors HcpR and LysX. After reconstruction of the corresponding metabolic pathways and regulatory interactions, we identified possible functions for a large number of previously uncharacterized genes covering a wide range of cellular functions.

Conclusions

Phylogenetically diverse δ-proteobacteria appear to have homologous regulatory components. This study for the first time demonstrates the adaptability of the comparative genomic approach to de novo reconstruction of a regulatory network in a poorly studied taxonomic group of bacteria. Recent efforts in large-scale functional genomic characterization of Desulfovibrio species will provide a unique opportunity to test and expand our predictions.

Contact order revisited: influence of protein size on the folding rate

Guided by the recent success of empirical model predicting the folding rates of small two-state folding proteins from the relative contact order (CO) of their native structures, by a theoretical model of protein folding that predicts that logarithm of the folding rate decreases with the protein chain length L as L(2/3), and by the finding that the folding rates of multistate folding proteins strongly correlate with their sizes and have very bad correlation with CO, we reexamined the dependence of folding rate on CO and L in attempt to find a structural parameter that determines folding rates for the totality of proteins. We show that the Abs_CO = CO x L, is able to predict rather accurately folding rates for both two-state and multistate folding proteins, as well as short peptides, and that this Abs_CO scales with the protein chain length as L(0.70 +/- 0.07) for the totality of studied single-domain proteins and peptides.

Biological networks

Recent advances in high-throughput methods have provided us with a first glimpse of the overall structure of molecular interaction networks in biological systems. Ultimately, we expect that such information will change how we think about biological systems in a fundamental way. Instead of viewing the genetic parts list of an organism as a loose collection of biochemical activities, in the best case, we anticipate discrete networks of function to bridge the gap between genotype and phenotype, and to do so in a more profound way than the current qualitative classification of linked reactions into familiar pathways, such as glycolysis and the MAPK signal transduction cascades. At the present time, however, we are still far from a complete answer to the most basic question: what can we learn about biology by studying networks? Promising steps in this direction have come from such diverse approaches as mathematical analysis of global network structure, partitioning networks into functionally related modules and motifs, and even de novo design of networks. A complete picture will probably require integrating the data obtained from all of these approaches with modeling efforts at many different levels of detail.

Simple physical models connect theory and experiment in protein folding kinetics

Our understanding of the principles underlying the protein-folding problem can be tested by developing and characterizing simple models that make predictions which can be compared to experimental data. Here we extend our earlier model of folding free energy landscapes, in which each residue is considered to be either folded as in the native state or completely disordered, by investigating the role of additional factors representing hydrogen bonding and backbone torsion strain, and by using a hybrid between the master equation approach and the simple transition state theory to evaluate kinetics near the free energy barrier in greater detail. Model calculations of folding phi-values are compared to experimental data for 19 proteins, and for more than half of these, experimental data are reproduced with correlation coefficients between r=0.41 and 0.88; calculations of transition state free energy barriers correlate with rates measured for 37 single domain proteins (r=0.69). The model provides insight into the contribution of alternative-folding pathways, the validity of quasi-equilibrium treatments of the folding landscape, and the magnitude of the Arrhenius prefactor for protein folding. Finally, we discuss the limitations of simple native-state-based models, and as a more general test of such models, provide predictions of folding rates and mechanisms for a comprehensive set of over 400 small protein domains of known structure.

Critical role of beta-hairpin formation in protein G folding

Comparison of the folding mechanisms of proteins with similar structures but very different sequences can provide fundamental insights into the determinants of protein folding mechanisms. Despite very little sequence similarity, the approximately 60 residue IgG binding domains of protein G and protein L both consist of a single helix packed against a four-stranded sheet formed by two symmetrically disposed beta-hairpins. We demonstrate that, as in the case of protein L, one of the two beta-turns of protein G is formed and the other disrupted in the folding transition state. Unlike protein L, however, in protein G it is the second beta-turn that is formed in the folding transition state ensemble. Substitution of an Asp residue by Ala in protein G that eliminates an i,i+2 side chain-main chain hydrogen bond in the second beta-turn slows the folding rate approximately 20-fold but has virtually no effect on the unfolding rate. Taken together with previous results, these findings suggest that the presence of an intact beta-turn in the folding transition state is a consequence of the overall topology of protein L and protein G, but the particular hairpin that is formed is determined by the detailed interatomic interactions that determine the free energies of formation of the isolated beta-hairpins.

Characterization of the unfolding pathway of the cell-cycle protein p13suc1 by molecular dynamics simulations: implications for domain swapping

BACKGROUND

The p13suc1 gene product is a member of the cks (cyclin-dependent protein kinase subunit) protein family and has been implicated in regulation of the cell cycle. Various crystal structures of suc1 are available, including a globular, monomeric form and a beta-strand exchanged dimer. It has been suggested that conversions between these forms, and perhaps others, may be important in the regulation of the cell cycle.

RESULTS

We have undertaken molecular dynamics simulations of protein unfolding to investigate the conformational properties of suc1. Unfolding transition states were identified for each of four simulations. These states contain some native secondary structure, primarily helix alpha1 and the core of the beta sheet. The hydrophobic core is loosely packed. Further unfolding leads to an intermediate state that is slightly more expanded than the transition state, but with considerably fewer nonlocal, tertiary packing contacts and less secondary structure. The helices are fluctuating but partially formed in the denatured state and beta2 and beta4 remain associated.

CONCLUSIONS

It appears that suc1 folds by a nucleation-condensation mechanism, similar to that observed for two-state folding proteins. However, suc1 forms an intermediate during unfolding and contains considerable residual structure in the denatured state. The stability of the beta2-beta4 residual structure is surprising, because beta4 is the strand involved in domain swapping. This stability suggests that the domain-swapping event, if physiologically relevant, may require the assistance of additional factors in vivo or occur early in the folding process.

Experiment and theory highlight role of native state topology in SH3 folding

We use a combination of experiments, computer simulations and simple model calculations to characterize, first, the folding transition state ensemble of the src SH3 domain, and second, the features of the protein that determine its folding mechanism. Kinetic analysis of mutations at 52 of the 57 residues in the src SH3 domain revealed that the transition state ensemble is even more polarized than suspected earlier: no single alanine substitution in the N-terminal 15 residues or the C-terminal 9 residues has more than a two-fold effect on the folding rate, while such substitutions at 15 sites in the central three-stranded beta-sheet cause significant decreases in the folding rate. Molecular dynamics (MD) unfolding simulations and ab initio folding simulations on the src SH3 domain exhibit a hierarchy of folding similar to that observed in the experiments. The similarity in folding mechanism of different SH3 domains and the similar hierarchy of structure formation observed in the experiments and the simulations can be largely accounted for by a simple native state topology-based model of protein folding energy landscapes.

Prediction of protein-folding mechanisms from free-energy landscapes derived from native structures

Guided by recent experimental results suggesting that protein-folding rates and mechanisms are determined largely by native-state topology, we develop a simple model for protein folding free-energy landscapes based on native-state structures. The configurations considered by the model contain one or two contiguous stretches of residues ordered as in the native structure with all other residues completely disordered; the free energy of each configuration is the difference between the entropic cost of ordering the residues, which depends on the total number of residues ordered and the length of the loop between the two ordered segments, and the favorable attractive interactions, which are taken to be proportional to the total surface area buried by the ordered residues in the native structure. Folding kinetics are modeled by allowing only one residue to become ordered/disordered at a time, and a rigorous and exact method is used to identify free-energy maxima on the lowest free-energy paths connecting the fully disordered and fully ordered configurations. The distribution of structure in these free-energy maxima, which comprise the transition-state ensemble in the model, are reasonably consistent with experimental data on the folding transition state for five of seven proteins studied. Thus, the model appears to capture, at least in part, the basic physics underlying protein folding and the aspects of native-state topology that determine protein-folding mechanisms.

Matching theory and experiment in protein folding

There has been considerable progress made over the past year in linking experimental and theoretical approaches to protein folding. Recent results from several independent lines of investigation suggest that protein folding mechanisms and landscapes are largely determined by the topology of the native state and are relatively insensitive to details of the interatomic interactions. This dependence on low-resolution structural features, rather than high-resolution detail, suggests that it should be possible to describe the fundamental physics of the folding process using relatively low-resolution models. Recent experiments have set benchmarks for testing new models and progress has been made in developing theoretical models for interpreting and predicting experimental results.

Evolutionary relationships among ammonia- and nitrite-oxidizing bacteria

Comparative 16S rRNA sequencing was used to evaluate phylogenetic relationships among selected strains of ammonia- and nitrite-oxidizing bacteria. All characterized strains were shown to be affiliated with the proteobacteria. The study extended recent 16S rRNA-based studies of phylogenetic diversity among nitrifiers by the comparison of eight strains of the genus Nitrobacter and representatives of the genera Nitrospira and Nitrospina. The later genera were shown to be affiliated with the delta subdivision of the proteobacteria but did not share a specific relationship to each other or to other members of the delta subdivision. All characterized Nitrobacter strains constituted a closely related assemblage within the alpha subdivision of the proteobacteria. As previously observed, all ammonia-oxidizing genera except Nitrosococcus oceanus constitute a monophyletic assemblage within the beta subdivision of the proteobacteria. Errors in the 16S rRNA sequences for two strains previously deposited in the databases by other investigators (Nitrosolobus multiformis C-71 and Nitrospira briensis C-128) were corrected. Consideration of physiology and phylogenetic distribution suggested that nitrite-oxidizing bacteria of the alpha and gamma subdivisions are derived from immediate photosynthetic ancestry. Each nitrifier retains the general structural features of the specific ancestor's photosynthetic membrane complex. Thus, the nitrifiers, as a group, apparently are not derived from an ancestral nitrifying phenotype.

Surface histamine receptors in rat peritoneal mast cells. Separation of receptor-bearing cells via binding to a histamine-protein conjugate, Cell-ect

The Seragen Cell-ect Total Histamine KitTM was used to separate histamine receptor-bearing rat peritoneal mast cells from cells lacking these receptors. It was found that approximately 75% of the peritoneal mast cell population carried cell surface histamine receptors. The results further suggest that the mast cell histamine receptors present are mainly of an H1-type, as judged by the capacity of a specific H1-antagonist to reduce the histamine receptor-dependent cell adhesion. Moreover, an H1-agonist is less efficient in this respect and an H2-antagonist does not affect the cell adhesion at all. A possible functional role for these receptors, however, remains to be clarified.

Cyclic nucleotide involvement in histamine release from mast cells--a reevaluation

Secretory events in cells in general are accompanied by increased levels of cyclic AMP. In mast cells, however, the pattern is reversed. Thus histamine release is associated with a fall in cAMP. It has been suggested that the lowered levels of cAMP lead to an increase in membrane permeability towards calcium and that an influx of such ions triggers the release mechanisms. It has further been reported that high levels of cAMP inhibit histamine release by decreasing the permeability. However, evidence has now accumulated indicating that this general concept is far too simplistic. Studies are reviewed which imply that there is little or no correlation between histamine release and intracellular levels of cyclic nucleotides. A new working hypothesis with respect to the role of these nucleotides in mast cell secretion is proposed.