OptZyme: Computational Enzyme Redesign Using Transition State Analogues

OptZyme is a new computational procedure for designing improved enzymatic activity (i.e., k_cat or k_cat/K_M) with a novel substrate. The key concept is to use transition state analogue compounds, which are known for many reactions, as proxies for the typically unknown transition state structures. Mutations that minimize the interaction energy of the enzyme with its transition state analogue, rather than with its substrate, are identified that lower the transition state formation energy barrier. Using Escherichia coli β-glucuronidase as a benchmark system, we confirm that K_M correlates (R² = 0.960) with the computed interaction energy between the enzyme and the para-nitrophenyl- β, D-glucuronide substrate, k_cat/K_M correlates (R² = 0.864) with the interaction energy of the transition state analogue, 1,5-glucarolactone, and k_cat correlates (R² = 0.854) with a weighted combination of interaction energies with the substrate and transition state analogue. OptZyme is subsequently used to identify mutants with improved K_M, k_cat, and k_cat/K_M for a new substrate, para-nitrophenyl- β, D-galactoside. Differences between the three libraries reveal structural differences that underpin improving K_M, k_cat, or k_cat/K_M. Mutants predicted to enhance the activity for para-nitrophenyl- β, D-galactoside directly or indirectly create hydrogen bonds with the altered sugar ring conformation or its substituents, namely H162S, L361G, W549R, and N550S.     

Conserved Prosegment Residues Stabilize a Late-Stage Folding Transition State of Pepsin Independently of Ground States

The native folding of certain zymogen-derived enzymes is completely dependent upon a prosegment domain to stabilize the folding transition state, thereby catalyzing the folding reaction. Generally little is known about how the prosegment accomplishes this task. It was previously shown that the prosegment catalyzes a late-stage folding transition between a stable misfolded state and the native state of pepsin. In this study, the contributions of specific prosegment residues to catalyzing pepsin folding were investigated by introducing individual Ala substitutions and measuring the effects on the bimolecular folding reaction between the prosegment peptide and pepsin. The effects of mutations on the free energies of the individual misfolded and native ground states and the transition state were compared using measurements of prosegment-pepsin binding and folding kinetics. Five out of the seven prosegment residues examined yielded relatively large kinetic effects and minimal ground state perturbations upon mutation, findings which indicate that these residues form strengthened and/or non-native contacts in the transition state. These five residues are semi- to strictly conserved, while only a non-conserved residue had no kinetic effect. One conserved residue was shown to form native structure in the transition state. These results indicated that the prosegment, which is only 44 residues long, has evolved a high density of contacts that preferentially stabilize the folding transition state over the ground states. It is postulated that the prosegment forms extensive non-native contacts during the process of catalyzing correct inter- and intra-domain contacts during the final stages of folding. These results have implications for understanding the folding of multi-domain proteins and for the evolution of prosegment-catalyzed folding.     

Comparative Analysis of Light-Harvesting Antennae and State Transition in chlorina and cpSRP Mutants

State transitions in photosynthesis provide for the dynamic allocation of a mobile fraction of light-harvesting complex II (LHCII) to photosystem II (PSII) in state I and to photosystem I (PSI) in state II. In the state I-to-state II transition, LHCII is phosphorylated by STN7 and associates with PSI to favor absorption cross-section of PSI. Here, we used Arabidopsis (Arabidopsis thaliana) mutants with defects in chlorophyll (Chl) b biosynthesis or in the chloroplast signal recognition particle (cpSRP) machinery to study the flexible formation of PS-LHC supercomplexes. Intriguingly, we found that impaired Chl b biosynthesis in chlorina1-2 (ch1-2) led to preferentially stabilized LHCI rather than LHCII, while the contents of both LHCI and LHCII were equally depressed in the cpSRP43-deficient mutant (chaos). In view of recent findings on the modified state transitions in LHCI-deficient mutants (Benson et al., 2015), the ch1-2 and chaos mutants were used to assess the influence of varying LHCI/LHCII antenna size on state transitions. Under state II conditions, LHCII-PSI supercomplexes were not formed in both ch1-2 and chaos plants. LHCII phosphorylation was drastically reduced in ch1-2, and the inactivation of STN7 correlates with the lack of state transitions. In contrast, phosphorylated LHCII in chaos was observed to be exclusively associated with PSII complexes, indicating a lack of mobile LHCII in chaos. Thus, the comparative analysis of ch1-2 and chaos mutants provides new evidence for the flexible organization of LHCs and enhances our understanding of the reversible allocation of LHCII to the two photosystems.In oxygenic photosynthesis, PSII and PSI function in series to convert light energy into the chemical energy that fuels multiple metabolic processes. Most of this light energy is captured by the chlorophyll (Chl) and carotenoid pigments in the light-harvesting antenna complexes (LHCs) that are peripherally associated with the core complexes of both photosystems (Wobbe et al., 2016). However, since the two photosystems exhibit different absorption spectra (Nelson and Yocum, 2006; Nield and Barber, 2006; Qin et al., 2015), PSI or PSII is preferentially excited under naturally fluctuating light intensities and qualities. To optimize photosynthetic electron transfer, the excitation state of the two photosystems must be rebalanced in response to changes in lighting conditions. To achieve this, higher plants and green algae require rapid and precise acclimatory mechanisms to adjust the relative absorption cross-sections of the two photosystems.To date, the phenomenon of state transitions is one of the well-documented short-term acclimatory mechanisms. It allows a mobile portion of the light-harvesting antenna complex II (LHCII) to be allocated to either photosystem, depending on the spectral composition and intensity of the ambient light (Allen and Forsberg, 2001; Rochaix, 2011; Goldschmidt-Clermont and Bassi, 2015; Gollan et al., 2015). State transitions are driven by the redox state of the plastoquinone (PQ) pool (Vener et al., 1997; Zito et al., 1999). When PSI is preferentially excited (by far-red light), the PQ pool is oxidized and all the LHCII is associated with PSII. This allocation of antenna complexes is defined as state I. When light conditions (blue/red light or low light) favor exciton trapping of PSII, the transition from state I to state II occurs. The over-reduced PQ pool triggers the activation of the membrane-localized Ser-Thr kinase STN7, which phosphorylates an N-terminal Thr on each of two major LHCII proteins, LHCB1 and LHCB2 (Allen, 1992; Bellafiore et al., 2005; Shapiguzov et al., 2016). Phosphorylation of LHCII results in the dissociation of LHCII from PSII and triggers its reversible relocation to PSI (Allen, 1992; Rochaix, 2011). Conversely, when the PQ pool is reoxidized, STN7 is inactivated and the constitutively active, thylakoid-associated phosphatase TAP38/PPH1 dephosphorylates LHCII, which then reassociates with PSII (Pribil et al., 2010; Shapiguzov et al., 2010). The physiological significance of state transitions has been demonstrated by the reduction in growth rate seen in the stn7 knock-out mutant under fluctuating light conditions (Bellafiore et al., 2005; Tikkanen et al., 2010).The canonical state transitions model implies spatial and temporal regulation of the allocation of LHC between the two spatially segregated photosystems (Dekker and Boekema, 2005). PSII-LHCII supercomplexes are organized in a tightly packed form in the stacked grana regions of thylakoid membranes, while PSI-LHCI supercomplexes are mainly localized in the nonstacked stromal lamellae and grana margin regions (Dekker and Boekema, 2005; Haferkamp et al., 2010). It has been proposed that, in the grana margin regions, which harbor LHCII and both photosystems, LHCII can migrate rapidly between them (Albertsson et al., 1990; Albertsson, 2001). This idea is supported by the recent discovery of mega complexes containing both photosystems in the grana margin regions (Yokono et al., 2015). Furthermore, phosphorylation of LHCII was found to increase not only the amount of PSI found in the grana margin region of thylakoid membranes (Tikkanen et al., 2008a), but also to modulate the pattern of PSI-PSII megacomplexes under changing light conditions (Suorsa et al., 2015). Nonetheless, open questions remain in relation to the physiological significance of the detection of phosphorylated LHCII in all thylakoid regions, even under the constant light conditions (Grieco et al., 2012; Leoni et al., 2013; Wientjes et al., 2013), although LHCII phosphorylation has been shown to modify the stacking of thylakoid membranes (Chuartzman et al., 2008; Pietrzykowska et al., 2014).State I-to-state II transition is featured by the formation of LHCII-PSI-LHCI supercomplexes, in which LHCII favors the light-harvesting capacity of PSI. Recently, LHCII-PSI-LHCI supercomplexes have been successfully isolated and purified using various detergents (Galka et al., 2012; Drop et al., 2014; Crepin and Caffarri, 2015) or a styrene-maleic acid copolymer (Bell et al., 2015). These findings yielded further insights into the reorganization of supercomplexes associated with state transitions, and it was suggested that phosphorylation of LHCB2 rather than LHCB1 is the essential trigger for the formation of state transition supercomplexes (Leoni et al., 2013; Pietrzykowska et al., 2014; Crepin and Caffarri, 2015; Longoni et al., 2015). Furthermore, characterization of mutants deficient in individual PSI core subunits indicates that PsaH, L, and I are required for docking of LHCII at PSI (Lunde et al., 2000; Zhang and Scheller, 2004; Kouril et al., 2005; Plöchinger et al., 2016).Recently, the state transition capacity has been characterized in the Arabidopsis (Arabidopsis thaliana) mutants with missing LHCI components. Although the Arabidopsis knock-out mutants lacking one of the four LHCI proteins (LHCA1-4) showed enhanced accumulation of LHCII-PSI complexes, the absorption cross-section of PSI under state II conditions was still compromised in the lhca1-4 mutants, and it is suggested that LHCI mediates the detergent-sensitive interaction between ‘extra LHCII’ and PSI (Benson et al., 2015; Grieco et al., 2015). Furthermore, the Arabidopsis mutant ΔLhca lacking all LHCA1-4 proteins was shown to be compensated for the deficiency of LHCI by binding LHCII under state II conditions (Bressan et al., 2016). In spite of this finding, the significant reduction in the absorption cross-section of PSI was still observed in the ΔLhca mutant, suggesting a substantial role of LHCI in light absorption under canopy conditions (Bressan et al., 2016). However, these findings emphasize the acclimatory function of state transitions in balancing light absorption capacity between the two photosystems by modifying their relative antenna size and imply the dynamic and variable organization of PS-LHC supercomplexes.LHC proteins are encoded by the nuclear Lhc superfamily (Jansson, 1994). The biogenesis of LHCs includes the cytoplasmic synthesis of the LHC precursor proteins, their translocation into chloroplasts via the TOC/TIC complex, and their posttranslational targeting and integration into the thylakoid membranes by means of the chloroplast signal recognition particle (cpSRP) machinery (Jarvis and Lopez-Juez, 2013). The posttranslational cpSRP-dependent pathway for the final translocation of LHC proteins into the thylakoid membrane includes interaction of cpSRP43 with LHC apo-proteins and recruitment of cpSRP54 to form a transit complex. Then binding of this tripartite cpSRP transit complex to the SRP receptor cpFtsY follows, which supports docking of the transit complex to thylakoid membranes and its association with the LHC translocase ALB3. Ultimately, ALB3 inserts LHC apo-proteins into the thylakoid membrane (Richter et al., 2010). Importantly, stoichiometric amounts of newly synthesized Chl a and Chl b as well as carotenoid are inserted into the LHC apo-proteins by unknown mechanisms to form the functional LHCs that associate with the core complexes of both photosystems in the thylakoid membranes (Dall’Osto et al., 2015; Wang and Grimm, 2015).The first committed steps in Chl synthesis occur in the Mg branch of the tetrapyrrole biosynthesis pathway. 5-Aminolevulinic acid synthesis provides the precursor for the formation of protoporphyrin IX, which is directed into the Mg branch (Tanaka and Tanaka, 2007; Brzezowski et al., 2015). Chl synthesis ends with the conversion of Chl a to Chl b catalyzed by Chl a oxygenase (CAO; Tanaka et al., 1998; Tomitani et al., 1999). It has been hypothesized that coordination between Chl synthesis and the posttranslational cpSRP pathway is a prerequisite for the efficient integration of Chls into LHC apo-proteins.In this study, we intend to characterize the assembly of LHCs when the availability of Chl molecules or the integration of LHC apo-proteins into thylakoid membranes is limiting. To this end, we compared the assembly of LHCs and the organization of PS-LHC complexes in two different sets of Arabidopsis mutants. Firstly, we used the chlorina1-2 (ch1-2) mutant, which is defective in the CAO gene. The members of the second set of mutants carry knock-out mutations in genes involved in the chloroplast SRP pathway (Richter et al., 2010).Our studies revealed distinct accumulation of PS-LHC supercomplexes between the two sets of mutant relative to wild-type plants. In spite of the defect in synthesis of Chl b, ch1-2 retains predominantly intact PSI-LHCI supercomplexes but has strongly reduced amounts of LHCII. In contrast, the chaos (cpSRP43) mutant exhibits synchronously reduced contents of both LHCI and LHCII, which results in the accumulation of PS core complexes without accompanying LHCs. Thus, the distribution of LHCs in the thylakoid membranes of the two mutants, ch1-2 and chaos, were explored under varying light conditions with the aim of elucidating the influence of modified LHCI/LHCII antenna size on state transitions. Our results contribute to an expanding view on the variety of photosynthetic complexes, which can be observed in Arabidopsis plants with specified mutations in LHC biogenesis.     

Capturing Transition Paths and Transition States for Conformational Rearrangements in the Ribosome

To reveal the molecular determinants of biological function, one seeks to characterize the interactions that are formed in conformational and chemical transition states. In other words, what interactions govern the molecule’s energy landscape? To accomplish this, it is necessary to determine which degrees of freedom can unambiguously identify each transition state. Here, we perform simulations of large-scale aminoacyl-transfer RNA (aa-tRNA) rearrangements during accommodation on the ribosome and project the dynamics along experimentally accessible atomic distances. From this analysis, we obtain evidence for which coordinates capture the correct number of barrier-crossing events and accurately indicate when the aa-tRNA is on a transition path. Although a commonly used coordinate in single-molecule experiments performs poorly, this study implicates alternative coordinates along which rearrangements are accurately described as diffusive movements across a one-dimensional free-energy profile. From this, we provide the theoretical foundation required for single-molecule techniques to uncover the energy landscape governing aa-tRNA selection by the ribosome.     

A Semiempirical Transition State Structure for the First Step in the Alkaline Hydrolysis of Cocaine. Comparison between the Transition State Structure,the Phosphonate Monoester Transition State Analog,and a Newly Designed Thiophosphonate Transition State Analog

Semiempirical molecular orbital calculations have been performed for the first step in the alkaline hydrolysis of the neutral benzoylester of cocaine. Successes, failures, and limitations of these calculations are reviewed. A PM3 calculated transition state structure is compared with the PM3 calculated structure for the hapten used to induce catalytic antibodies for the hydrolysis of cocaine. Implications of these calculations for the computer–aided design of transition state analogs for the induction of catalytic antibodies are discussed.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s0089460020062     

Analysis of Discrete Bioregulatory Networks Using Symbolic Steady States

A discrete model of a biological regulatory network can be represented by a discrete function that contains all available information on interactions between network components and the rules governing the evolution of the network in a finite state space. Since the state space size grows exponentially with the number of network components, analysis of large networks is a complex problem. In this paper, we introduce the notion of symbolic steady state that allows us to identify subnetworks that govern the dynamics of the original network in some region of state space. We state rules to explicitly construct attractors of the system from subnetwork attractors. Using the results, we formulate sufficient conditions for the existence of multiple attractors resp. a cyclic attractor based on the existence of positive resp. negative feedback circuits in the graph representing the structure of the system. In addition, we discuss approaches to finding symbolic steady states. We focus both on dynamics derived via synchronous as well as asynchronous update rules. Lastly, we illustrate the results by analyzing a model of T helper cell differentiation.     

Structural Differentiation of Graphs Using Hosoya-Based Indices

In this paper, we introduce the Hosoya-Spectral indices and the Hosoya information content of a graph. The first measure combines structural information captured by partial Hosoya polynomials and graph spectra. The latter is a graph entropy measure which is based on blocks consisting of vertices with the same partial Hosoya polynomial. We evaluate the discrimination power of these quantities by interpreting numerical results.     

Analyzing and Synthesizing Phylogenies Using Tree Alignment Graphs

Phylogenetic trees are used to analyze and visualize evolution. However, trees can be imperfect datatypes when summarizing multiple trees. This is especially problematic when accommodating for biological phenomena such as horizontal gene transfer, incomplete lineage sorting, and hybridization, as well as topological conflict between datasets. Additionally, researchers may want to combine information from sets of trees that have partially overlapping taxon sets. To address the problem of analyzing sets of trees with conflicting relationships and partially overlapping taxon sets, we introduce methods for aligning, synthesizing and analyzing rooted phylogenetic trees within a graph, called a tree alignment graph (TAG). The TAG can be queried and analyzed to explore uncertainty and conflict. It can also be synthesized to construct trees, presenting an alternative to supertrees approaches. We demonstrate these methods with two empirical datasets. In order to explore uncertainty, we constructed a TAG of the bootstrap trees from the Angiosperm Tree of Life project. Analysis of the resulting graph demonstrates that areas of the dataset that are unresolved in majority-rule consensus tree analyses can be understood in more detail within the context of a graph structure, using measures incorporating node degree and adjacency support. As an exercise in synthesis (i.e., summarization of a TAG constructed from the alignment trees), we also construct a TAG consisting of the taxonomy and source trees from a recent comprehensive bird study. We synthesized this graph into a tree that can be reconstructed in a repeatable fashion and where the underlying source information can be updated. The methods presented here are tractable for large scale analyses and serve as a basis for an alternative to consensus tree and supertree methods. Furthermore, the exploration of these graphs can expose structures and patterns within the dataset that are otherwise difficult to observe.     

Woodland Restoration in the Western Australian Wheatbelt: A Conceptual Framework Using a State and Transition Model

Woodlands dominated by Eucalyptus spp. in temperate southeastern and southwestern Australia have been extensively cleared for agriculture and are often badly degraded by livestock grazing. This has resulted in the loss of biodiversity and widespread land degradation. The continuing decline of these woodlands has become a concern for the conservation of biodiversity, and there is a growing interest among farmers, land managers, and researchers in developing techniques for restoring them. Currently few scientific guidelines exist for undertaking woodland restoration programs. We use a state and transition model to develop hypotheses on restoration strategies for salmon gum (Eucalyptus salmonophloia) woodlands. We consider that this approach provides a suitable framework for organizing knowledge and identifying areas where further information is needed, and hence provides a useful starting point for a restoration program. The model has the potential to provide a tool for land managers with which they can assess the action and effort needed to undertake woodland restoration in agricultural landscapes.     

Identifying Movement States From Location Data Using Cluster Analysis

ABSTRACT Animal movement studies regularly use movement states (e.g., slow and fast) derived from remotely sensed locations to make inferences about strategies of resource use. However, the number of movement state categories used is often arbitrary and rarely inferred from the data. Identifying groups with similar movement characteristics is a statistical problem. We present a framework based on k-means clustering and gap statistic for evaluating the number of movement states without making a priori assumptions about the number of clusters. This allowed us to distinguish 4 movement states using turning angle and step length derived from Global Positioning System locations and head movements derived from tip switches in a neck collar of free-ranging elk (Cervus elaphus) in west central Alberta, Canada. Based on movement characteristics and on the linkage between each state and landscape features, we were able to identify inter-patch movements, intra-patch foraging, rest, and inter-patch foraging movements. Linking behavior to environment (e.g., state-dependent habitat use) can inform decisions on landscape management for wildlife.     

Transition Confidence and Modified Mean Values: Confidence Measures for Hypotheses of Character State Transition between Nodes and Ancestral State Optimizations

A measure is introduced to quantify the reliability of character state transitions and character state stasis (lack of character state change between nodes). The mean value of a character state at an internal node is altered to obtain a measure termed the modified mean. Modified mean values at connecting nodes are employed to determine the support for hypotheses of character state transition/stasis between the nodes. The lowest of the two neighboring modified mean values implied by an hypothesis of character state transition/stasis is used as a direct measure of support for that hypotheses of character state transformation or stasis between the respective nodes.     

Free-Energy Perturbation Simulation on Transition States and Redesign of Butyrylcholinesterase

It is recognized that an ideal anti-cocaine treatment is to accelerate cocaine metabolism by producing biologically inactive metabolites via a route similar to the primary cocaine-metabolizing pathway, i.e., butyrylcholinesterase (BChE)-catalyzed hydrolysis of cocaine. BChE mutants with a higher catalytic activity against (-)-cocaine are highly desired for use as an exogenous enzyme in humans. To develop a rational design for high-activity mutants, we carried out free-energy perturbation (FEP) simulations on various mutations of the transition-state structures in addition to the corresponding free-enzyme structures by using an extended FEP procedure. The FEP simulations on the mutations of both the free-enzyme and transition-state structures allowed us to calculate the mutation-caused shift of the free-energy change from the free enzyme (BChE) to the transition state, and thus to theoretically predict the mutation-caused shift of the catalytic efficiency (k_cat/K_M). The computational predictions are supported by the kinetic data obtained from the wet experiments, demonstrating that the FEP-based computational design approach is promising for rational design of high-activity mutants of an enzyme. One of the BChE mutants designed and discovered in this study has an ∼1800-fold improved catalytic efficiency against (-)-cocaine compared to wild-type BChE. The high-activity mutant may be therapeutically valuable.     

Effect of Solvation on Ozonolysis Reaction Intermediates and Transition States

Electrostatic solvent effects on the ozonolysis of ethylene have been investigated using correlated ab initio and density functional approaches. We use a simple polarizable continuum model for the solvent. It allows us to evaluate the medium effect on both the electronic and nuclear structure of the chemical species involved in the reaction. The computations confirm that basically the reaction proceeds through the Criegee mechanism. However, formation of the van der Waals complexes ethyl-ene/ozone and carbonyl oxide/formaldehyde also appears to play a role. All the calculated species are stabilized with respect to the reactants except the transition state corresponding to the primary ozonide formation. In general, electrostatic solvent effects are relatively small for activation barriers of single reaction steps and more substantial for the corresponding reaction energies. Moreover, the medium significantly modifies the structure of some species for which polarization effects are crucial.     

Vaginal Microbiome Characterization of Nellore Cattle Using Metagenomic Analysis

Understanding of microbial communities inhabiting cattle vaginal tract may lead to a better comprehension of bovine physiology and reproductive health being of great economic interest. Up to date, studies involving cattle microbiota are focused on the gastrointestinal tract, and little is known about the vaginal microbiota. This study aimed to investigate the vaginal microbiome in Nellore cattle, heifers and cows, pregnant and non-pregnant, using a culture independent approach. The main bacterial phyla found were Firmicutes (~40–50%), Bacteroidetes (~15–25%) and Proteobacteria (~5–25%), in addition to ~10–20% of non-classified bacteria. 45–55% of the samples were represented by only ten OTUs: Aeribacillus, Bacteroides, Clostridium, Ruminococcus, Rikenella, Alistipes, Bacillus, Eubacterium, Prevotella and non-classified bacteria. Interestingly, microbiota from all 20 animals could be grouped according to the respiratory metabolism of the main OTUs found, creating three groups of vaginal microbiota in cattle. Archaeal samples were dominated by the Methanobrevibacter genus (Euryarchaeota, ~55–70%). Ascomycota was the main fungal phylum (~80–95%) and Mycosphaerella the most abundant genus (~70–85%). Hormonal influence was not clear, but a tendency for the reduction of bacterial and increase of archaeal populations in pregnant animals was observed. Eukaryotes did not vary significantly between pregnant and non-pregnant animals, but tended to be more abundant on cows than on heifers. The present work describes a great microbial variability in the vaginal community among the evaluated animals and groups (heifers and cows, pregnant and non-pregnant), which is significantly different from the findings previously reported using culture dependent methods, pointing out the need for further studies on this issue. The microbiome found also indicates that the vaginal colonization appears to be influenced by the gastrointestinal community.     

Extracting Transition Rates in Particle Tracking Using Analytical Diffusion Distribution Analysis

Single-particle tracking is an important technique in the life sciences to understand the kinetics of biomolecules. The analysis of apparent diffusion coefficients in vivo, for example, enables researchers to determine whether biomolecules are moving alone, as part of a larger complex, or are bound to large cellular components such as the membrane or chromosomal DNA. A remaining challenge has been to retrieve quantitative kinetic models, especially for molecules that rapidly switch between different diffusional states. Here, we present analytical diffusion distribution analysis (anaDDA), a framework that allows for extracting transition rates from distributions of apparent diffusion coefficients calculated from short trajectories that feature less than 10 localizations per track. Under the assumption that the system is Markovian and diffusion is purely Brownian, we show that theoretically predicted distributions accurately match simulated distributions and that anaDDA outperforms existing methods to retrieve kinetics, especially in the fast regime of 0.1–10 transitions per imaging frame. AnaDDA does account for the effects of confinement and tracking window boundaries. Furthermore, we added the option to perform global fitting of data acquired at different frame times to allow complex models with multiple states to be fitted confidently. Previously, we have started to develop anaDDA to investigate the target search of CRISPR-Cas complexes. In this work, we have optimized the algorithms and reanalyzed experimental data of DNA polymerase I diffusing in live Escherichia coli. We found that long-lived DNA interaction by DNA polymerase are more abundant upon DNA damage, suggesting roles in DNA repair. We further revealed and quantified fast DNA probing interactions that last shorter than 10 ms. AnaDDA pushes the boundaries of the timescale of interactions that can be probed with single-particle tracking and is a mathematically rigorous framework that can be further expanded to extract detailed information about the behavior of biomolecules in living cells.     

Densification and State Transition Across the Missouri Ozarks Landscape

World-wide, some biomes are densifying, or increasing in dense woody vegetation, and shifting to alternative stable states. We quantified densification and state transition between forests ecosystems in historical (ca. 1815–1850) and current (2004–2008) surveys of the Missouri Ozark Highlands, a 5-million ha landscape in southern Missouri, USA. To estimate density of historical forests, we used the Morisita plotless density estimator and applied corrections for surveyor bias. For contemporary forests, we used known densities at plots to predict continuous densities with random forests, an ensemble regression tree method. We also calculated basal area and percent stocking to determine changes in wood volume. Historical forests had densities ranging from about 75 to 320 trees/ha. Current forest densities were about 2.3 times greater and more uniform, at about 300–400 trees/ha (DBH ≥ 12.7 cm). Not all forests have increased in basal area and percent stocking because trees in contemporary forests are smaller in diameter than historical forests. Although oak species still are dominant (as defined by ≥10% composition), oak dominance is being replaced by many fire-sensitive species, of which only eastern redcedar and maples have become dominant. Densification and community changes in functional traits have produced a state transition from open oak forest ecosystems to predominantly closed eastern broadleaf forests in the Missouri Ozarks. This shift is not at equilibrium, as fire-sensitive species are continuing to increase and turnover in long-lived oaks is slow.     

大鼠慢波睡眠和睡眠过渡期脑电信号的小波变换分析

EEG信号经常包含许多快速的时变信息 ,将较长时间段的EEG信号近似看作平稳信号 ,进行FFT谱估计 ,存在其局限性。应用多分辨率小波变换方法 ,在频域和时域上可以同时定位分析大鼠慢波睡眠和睡眠过渡期脑电的动态变化特性。采用慢性埋植电极记录自由活动大鼠的皮层脑电 ,将信号用小波变换分解成δ、θ、α和 β四个分量 ,求各分量的功率和功率百分比的时间变化曲线 ,并与FFT功率谱分析结果进行比较。结果表明 :慢波睡眠期EEG中有 2 6 .2 %± 7.7%的时间段上δ分量功率小于总功率的 5 0 % ,且δ分量较大时 ,其他分量较小 ;δ分量较小时 ,其他分量较大 ,差别显著。此结果揭示了δ节律与θ和α节律之间的一种互补关系。而传统的FFT功率谱分析方法只能显示δ分量为主 (占总功率 70 .6 %± 6 .4 % )的功率谱 ,不能提供时变信息。对于睡眠过渡期的非稳态EEG信号 ,利用小波变换分解得到的θ和α分量可以鉴别出睡眠纺锤波 ,计算睡眠纺锤波的平均持续时间 ,并比较纺锤波和非纺锤波时期各个频谱分量的变化情况。由此可见 ,小波变换可用于计算新的EEG时频定量分析指标用于分析生理、病理和药理作用引起的睡眠EEG的变化过程 ,以弥补传统FFT功率谱分析的不足之处