MSMBuilder: Statistical Models for Biomolecular Dynamics

MSMBuilder is a software package for building statistical models of high-dimensional time-series data. It is designed with a particular focus on the analysis of atomistic simulations of biomolecular dynamics such as protein folding and conformational change. MSMBuilder is named for its ability to construct Markov state models (MSMs), a class of models that has gained favor among computational biophysicists. In addition to both well-established and newer MSM methods, the package includes complementary algorithms for understanding time-series data such as hidden Markov models and time-structure based independent component analysis. MSMBuilder boasts an easy to use command-line interface, as well as clear and consistent abstractions through its Python application programming interface. MSMBuilder was developed with careful consideration for compatibility with the broader machine learning community by following the design of scikit-learn. The package is used primarily by practitioners of molecular dynamics, but is just as applicable to other computational or experimental time-series measurements.     

Rates and patterns of mitochondrial DNA sequence evolution in fringilline finches (fringilla spp.) and the greenfinch (Carduelis chloris)

Rates and patterns of evolution in partial sequences of five mitochondrial genes (cytochrome b, ATPase 6, NADH dehydrogenase subunit 5, tRNA(Glu), and the control region) were compared among taxa in the passerine bird genera Fringilla and Carduelis. Rates of divergence do not vary significantly among genes, even in comparisons with the control region. Rate variation among lineages is significant only for the control region and NADH dehydrogenase subunit 5, and patterns of variation are consistent with the expectations of neutral theory. Base composition is biased in all genes but is stationary among lineages, and there is evidence for directional mutation pressure only in the control region. Despite these similarities, patterns of substitution differ among genes, consistent with alternative regimes of selective constraint. Rates of nonsynonymous substitution are higher in NADH dehydrogenase subunit 5 than in other protein-coding genes, and transitions exist in elevated proportions relative to transversions. Transitions appear to accumulate linearly with time in tRNA(Glu), and despite exhibiting the highest overall rate of divergence among species, there are no transversional changes in this gene. Finally, for resolving phylogenetic relationships among Fringilla taxa, the combined protein-coding data are broadly similar to those of the control region in terms of phylogenetic informativeness and statistical support.      

NADH:hydroxypyruvate reductase and NADPH:glyoxylate reductase in algae: partial purification and characterization from Chlamydomonas reinhardtii

Hydroxypyruvate and glyoxylate reductase activities were measured in extracts from the unicellular green algae, Chlamydomonas reinhardtii, Chlorella vulgaris, Chlorella miniata, and Dunaliella tertiolecta. Only trace levels of these activities were detectable in the blue-green algae, Anabaena variabilis and Synechococcus leopoliensis. A NADH-dependent hydroxypyruvate reductase was purified 130-fold from Chlamydomonas to a specific activity of 18 mumol NADH oxidized X min-1 X mg protein-1. The pH optimum was 5.0 to 7.0 in the presence of phosphate and the Km(hydroxypyruvate) was 0.05 mM. Substrate inhibition by hydroxypyruvate could be partially relieved by phosphate. The molecular weight, estimated by gel filtration, was 96,000. NADH-dependent glyoxylate reductase activity copurified with the hydroxypyruvate reductase. The Km(glyoxylate) was 10 mM, and the pH optimum was 4.5 to 8.5. A specific NADPH:glyoxylate reductase was also partially purified which did not reduce hydroxypyruvate or pyruvate. The NADPH:glyoxylate reductase had a Km(glyoxylate) of 0.1 mM and a pH optimum of 5.0 to 9.5. These reductases were compared with the pyruvate reductase of Chlamydomonas which also catalyzes the reduction of both hydroxypyruvate and glyoxylate.     

Magnesium-induced inner membrane aggregation in heart mitochondria: prevention and reversal by carboxyatractyloside and bongkrekic acid

Mg(2+) at an optimal concentration of 2mM (ph 6.5) induces large increases (up to 30 percent) in the optical density of bovine heart mitochondria incubated under conditions of low ionic strength (< approx. 0.01). The increases are associated with aggregation (sticking together) of the inner membranes and are little affected by changes in the energy status of the mitochondria. Virtually all of a number of other polyvalent cations tested and Ag(+) induce increases in mitochondrial optical density similar to those induced by Mg(2+), their approximate order of concentration effectiveness in respect to Mg(2+) being: La(3+) > Pb(2+) = Cu(2+) > Cd(2+) > Zn(2+) > Ag(+) > Mn(2+) > Ca(2+) > Mg(2+). With the exception of Mg(2+), all of these cations appear to induce swelling of the mitochondria concomitant with inner membrane aggregation. The inhibitors of the adenine nucleotide transport reaction carboxyatratyloside and bongkrekic acid are capable of preventing and reversing Mg(2+)-induced aggregation at the same low concentration required for complete inhibition of phosphorylating respiration, suggesting that they inhibit the aggregation by binding to the adenine nucleotide carrier. The findings are interpreted to indicate (a) that the inner mitochondrial membrane is normally prevented from aggregating by virtue of its net negative outer surface change, (b) that the cations induce the membrane to aggregate by binding at its outer surface, decreasing the net negative charge, and (c) that carboxyatractyloside and bongkrekic acid inhibit the aggregation by binding to the outer surface of the membrane, increasing the net negative charge.     

The novel flightless-I gene brings together two gene families, actin- binding proteins related to gelsolin and leucine-rich-repeat proteins involved in Ras signal transduction

The Drosophila melanogaster gene flightless-I, involved in gastrulation and muscle degeneration, has Caenorhabditis elegans and human homologues. In these highly conserved genes, two previously known gene families have been brought together, families encoding the actin- binding proteins related to gelsolin and the leucine-rich-repeat (LRR) group of proteins involved in protein-protein interactions. Both these gene families exhibit characteristics of molecular changes involving replication slippage and exon shuffling. Phylogenetic analyses of 19 amino acid sequences of 6 related protein types indicate that actin- associated proteins related to gelsolin are monophyletic to a common ancestor and include flightless proteins. Conversely, comparison of 24 amino acid sequences of LRR proteins including the flightless proteins indicates that flightless proteins are members of a structurally related subgroup. Included in the flightless cluster are human and mouse rsp-1 proteins involved in suppressing v-Ras transformation of cells and the membrane-associated yeast (Saccharomyces cerevisae) adenylate cyclase whose analogous LRRs are required for interaction with Ras proteins. There is a strong possibility that ligands for this group could be related and that flightless may have a similar role in Ras signal transduction. It is hypothesized that an ancestral monomeric gelsolin precursor protein has undergone at least four independent gene reorganization events to account for the structural diversity of the extant family of gelsolin-related proteins and that gene duplication and exon shuffling events occurred prior to or at the beginning of multicellular life, resulting in the evolution of some members of the family soon after the appearance of actin-type proteins.      

Comparison of linkage disequilibrium and haplotype diversity on macro- and microchromosomes in chicken

Background

Toll like receptors (TLR) play the central role in the recognition of pathogen associated molecular patterns (PAMPs). Mutations in the TLR1, TLR2 and TLR4 genes may change the ability to recognize PAMPs and cause altered responsiveness to the bacterial pathogens.

Results

The study presents association between TLR gene mutations and increased susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection. Novel mutations in TLR genes (TLR1- Ser150Gly and Val220Met; TLR2 – Phe670Leu) were statistically correlated with the hindrance in recognition of MAP legends. This correlation was confirmed subsequently by measuring the expression levels of cytokines (IL-4, IL-8, IL-10, IL-12 and IFN-γ) in the mutant and wild type moDCs (mocyte derived dendritic cells) after challenge with MAP cell lysate or LPS. Further in silico analysis of the TLR1 and TLR4 ectodomains (ECD) revealed the polymorphic nature of the central ECD and irregularities in the central LRR (leucine rich repeat) motifs.

Conclusion

The most critical positions that may alter the pathogen recognition ability of TLR were: the 9^th amino acid position in LRR motif (TLR1–LRR10) and 4^th residue downstream to LRR domain (exta-LRR region of TLR4). The study describes novel mutations in the TLRs and presents their association with the MAP infection.     

The oxidative photosynthetic carbon cycle or C2 cycle

The oxidative photosynthetic carbon cycle (or C₂ cycle) is the metabolic pathway responsible for photosynthetic oxygen uptake and the light‐dependent production of carbon dioxide that is termed photorespiration. The C₂ and reductive C₃ cycles coexist, and combined, represent total photosynthetic carbon metabolism. A brief historical review is presented beginning with the early observations of the oxygen inhibition of photosynthesis up to the discovery of the oxygenase activity associated with ribulose 1,5‐bisphosphate carboxylase/oxygenase. The properties and the role of the compartmentalization of the enzymes involved with the pathway and the transport of C₂ cycle intermediates are reviewed. The relationship of the C₂ cycle to photorespiratory nitrogen metabolism and other associated metabolic pathways and the properties and regulation of the C₂ cycle in diverse photosynthetic organisms are discussed.     

THE EXPRESSION OF APELIN AND ITS RECEPTOR APJ DURING DIFFERENT PHYSIOLOGICAL STAGES IN THE BOVINE OVARY

Recent studies implicate that apelin and its receptor APJ may have important role for the modulation of angiogenesis. The aim of this study was to further characterise the regulation of apelin/APJ system in bovine ovary. Experiment 1: corpora lutea (CL) were assigned to the following stages: days 1-2, 3-4, 5-7, 8-12, 13-16, >18 (after regression) of oestrous cycle and of gravidity (month <3, 3-5, 6-7 and >8). Experiment 2: Follicles during maturation were divided into granulosa cells (GC) and theca interna (TI) and were examined separately. Classification of follicles occurred by follicle size and oestradiol-17β (E2) concentration in the follicular fluid (FF) (<0.5 ng/ml, 0.5-5 ng/ml; 5-40 ng/ml; 40-180 ng/ml; >180 ng/ml). Real-time RT-PCR (qPCR) was applied to investigate mRNA expression of examined factors. In general, the expression level of apelin during the oestrous cycle was significantly higher compared to the one during pregnancy. Apelin mRNA levels were always high during the cycle with a tendency of decrease after CL regression. The APJ mRNA in the CL was significantly up regulated on days 5-7 and 8-12 followed by a decrease on days 13-16, and further on days >18. The expression of APJ does not show any significant regulation in the CL throughout pregnancy. The expression of apelin and APJ was not statistically regulated in GC, but was significantly up regulated in follicles with an E2 concentration of more than 5 ng/ml and showed an increase according to growth and maturation of follicles. In conclusion, our data suggest that apelin/APJ system is involved in the mechanism regulating angiogenesis during follicle maturation as well as during CL formation and function in the bovine ovary.