A Novel Signal Processing Measure to Identify Exact and Inexact Tandem Repeat Patterns in DNA Sequences

The identification and analysis of repetitive patterns are active areas of biological and computational research. Tandem repeats in telomeres play a role in cancer and hypervariable trinucleotide tandem repeats are linked to over a dozen major neurodegenerative genetic disorders. In this paper, we present an algorithm to identify the exact and inexact repeat patterns in DNA sequences based on orthogonal exactly periodic subspace decomposition technique. Using the new measure our algorithm resolves the problems like whether the repeat pattern is of period Open image in new window or its multiple (i.e., 2 Open image in new window , 3 Open image in new window , etc.), and several other problems that were present in previous signal-processing-based algorithms. We present an efficient algorithm of Open image in new window , where Open image in new window is the length of DNA sequence and Open image in new window is the window length, for identifying repeats. The algorithm operates in two stages. In the first stage, each nucleotide is analyzed separately for periodicity, and in the second stage, the periodic information of each nucleotide is combined together to identify the tandem repeats. Datasets having exact and inexact repeats were taken up for the experimental purpose. The experimental result shows the effectiveness of the approach.     

The evolution of domain-content in bacterial genomes

Background  

Across all sequenced bacterial genomes, the number of domains n _c in different functional categories c scales as a power-law in the total number of domains n, i.e. , with exponents α _c that vary across functional categories. Here we investigate the implications of these scaling laws for the evolution of domain-content in bacterial genomes and derive the simplest evolutionary model consistent with these scaling laws.     

An ant colony optimization algorithm for phylogenetic estimation under the minimum evolution principle

Background  

Distance matrix methods constitute a major family of phylogenetic estimation methods, and the minimum evolution (ME) principle (aiming at recovering the phylogeny with shortest length) is one of the most commonly used optimality criteria for estimating phylogenetic trees. The major difficulty for its application is that the number of possible phylogenies grows exponentially with the number of taxa analyzed and the minimum evolution principle is known to belong to the -hard class of problems.     

A replica exchange Monte Carlo algorithm for protein folding in the HP model

Background  

The ab initio protein folding problem consists of predicting protein tertiary structure from a given amino acid sequence by minimizing an energy function; it is one of the most important and challenging problems in biochemistry, molecular biology and biophysics. The ab initio protein folding problem is computationally challenging and has been shown to be -hard even when conformations are restricted to a lattice. In this work, we implement and evaluate the replica exchange Monte Carlo (REMC) method, which has already been applied very successfully to more complex protein models and other optimization problems with complex energy landscapes, in combination with the highly effective pull move neighbourhood in two widely studied Hydrophobic Polar (HP) lattice models.     

Variation in the Correlation of G + C Composition with Synonymous Codon Usage Bias among Bacteria

G + C composition at the third codon position (GC3) is widely reported to be correlated with synonymous codon usage bias. However, no quantitative attempt has been made to compare the extent of this correlation among different genomes. Here, we applied Shannon entropy from information theory to measure the degree of GC3 bias and that of synonymous codon usage bias of each gene. The strength of the correlation of GC3 with synonymous codon usage bias, quantified by a correlation coefficient, varied widely among bacterial genomes, ranging from Open image in new window 0.07 to 0.95. Previous analyses suggesting that the relationship between GC3 and synonymous codon usage bias is independent of species are thus inconsistent with the more detailed analyses obtained here for individual species.     

16S rRNA gene based analysis of <Emphasis Type="Italic">Enterobacter sakazakii</Emphasis> strains from different sources and development of a PCR assay for identification

Background  

E. sakazakii is considered to be an opportunistic pathogen, implicated in food borne diseases causing meningitis or enteritis especially in neonates and infants. Cultural standard identification procedures for E. sakazakii include the observation of yellow pigmentation of colonies and a positive glucosidase activity. Up to now, only one PCR system based on a single available 16S rRNA gene sequence has been published for E. sakazakii identification. However, in our hands a preliminary evaluation of this system to a number of target and non-target strains showed significant specificity problems of this system. In this study full-length 16S rRNA genes of thirteen E. sakazakii strains from food, environment and human origin as well as the type strain ATCC 51329 were sequenced. Based on this sequence data a new specific PCR system for E. sakazakii was developed and evaluated.     

Oxygen cost of dynamic or isometric exercise relative to recruited muscle mass

Background

Oxygen cost of different muscle actions may be influenced by different recruitment and rate coding strategies. The purpose of this study was to account for these strategies by comparing the oxygen cost of dynamic and isometric muscle actions relative to the muscle mass recruited via surface electrical stimulation of the knee extensors.

Methods

Comparisons of whole body pulmonary Δ $\dot{V}$O₂ were made in seven young healthy adults (1 female) during 3 minutes of dynamic or isometric knee extensions, both induced by surface electrical stimulation. Recruited mass was quantified in T₂ weighted spin echo magnetic resonance images.

Results

The Δ $\dot{V}$O₂ for dynamic muscle actions, 242 ± 128 ml ? min^-1 (mean ± SD) was greater (p = 0.003) than that for isometric actions, 143 ± 99 ml ? min^-1. Recruited muscle mass was also greater (p = 0.004) for dynamic exercise, 0.716 ± 282 versus 0.483 ± 0.139 kg. The rate of oxygen consumption per unit of recruited muscle (${\dot{V}\text{O}}_{2^{\text{RM}}}$) was similar in dynamic and isometric exercise (346 ± 162 versus 307 ± 198 ml ? kg^-1 ? min^-1; p = 0.352), but the ${\dot{V}\text{O}}_{2^{\text{RM}}}$ calculated relative to initial knee extensor torque was significantly greater during dynamic exercise 5.1 ± 1.5 versus 3.6 ± 1.6 ml ? kg^-1 ? Nm^-1 ? min^-1 (p = 0.019).

Conclusion

These results are consistent with the view that oxygen cost of dynamic and isometric actions is determined by different circumstances of mechanical interaction between actin and myosin in the sarcomere, and that muscle recruitment has only a minor role.

     

Gene Systems Network Inferred from Expression Profiles in Hepatocellular Carcinogenesis by Graphical Gaussian Model

Hepatocellular carcinoma (HCC) in a liver with advanced-stage chronic hepatitis C (CHC) is induced by hepatitis C virus, which chronically infects about 170 million people worldwide. To elucidate the associations between gene groups in hepatocellular carcinogenesis, we analyzed the profiles of the genes characteristically expressed in the CHC and HCC cell stages by a statistical method for inferring the network between gene systems based on the graphical Gaussian model. A systematic evaluation of the inferred network in terms of the biological knowledge revealed that the inferred network was strongly involved in the known gene-gene interactions with high significance Open image in new window , and that the clusters characterized by different cancer-related responses were associated with those of the gene groups related to metabolic pathways and morphological events. Although some relationships in the network remain to be interpreted, the analyses revealed a snapshot of the orchestrated expression of cancer-related groups and some pathways related with metabolisms and morphological events in hepatocellular carcinogenesis, and thus provide possible clues on the disease mechanism and insights that address the gap between molecular and clinical assessments.     

Results of the oxygen Fick method in a closed blood circulation model including “total arteriovenous diffusive shunt of oxygen”

It is considered that arteriovenous diffusive shunts of oxygen may cause inaccuracy of the oxygen Fick method as where is the pulmonary oxygen uptake, is the cardiac output, and CaO₂ and CvO₂ are the arterial and venous oxygen contents, respectively. A simple circulation model, including the whole circulation with nine well-mixed compartments (C1, ... C9), is constructed: the is assigned as constant as 6000 ml min⁻¹; the blood portions of 60 ml move at an interval of 600 ms. C1 and C2 compartments, each having 60 ml volume, represent the blood of pulmonary microcirculation, C3 represents the arterial blood with a volume of 1500 ml, C4, ..., C8, each also having a volume of 60 ml, represent the blood of peripheral microcirculation, whereas C9 represents the venous blood with a volume of 3000 ml. The pulmonary oxygen uptake , related to C1 and C2, the oxygen release , related to C4,...,C8, as well as a “total arteriovenous diffusive shunt of oxygen” , from the arterial blood (C3) to the venous blood (C9), are calculated simultaneously. The alveolar gas has a constant oxygen partial pressure, and the pulmonary diffusion capacity is also constant; similar to modeling the pulmonry, oxygen diffusion, constant partial oxygen pressures for all peripheral tissues as well as constant diffusion capacities for all peripheral oxygen diffusion are also assigned. The diffusion capacities for the (between C3 and C9) are arbitrarily assigned. The Fick method gives incorrect results depending on the total arteriovenous diffusive shunt of oxygen . But the mechanism determining the magnitude of remains unclear.     

An evolutionary analytical model of a complementary circular code simulating the protein coding genes, the 5′ and 3′ regions

The self-complementary subset ∪{AAA,TTT} with = {AAC, AAT, ACC, ATC, ATT, CAG, CTC, CTG, GAA, GAC, GAG, GAT, GCC, GGC, GGT, GTA, GTC, GTT, TAC, TTC} of 22 trinucleotides has a preferential occurrence in the frame 0 (reading frame established by the ATG start trinucleotide) of protein (coding) genes of both prokaryotes and eukaryotes. The subsets ∪{CCC} and ∪{GGG} of 21 trinucleotides have a preferential occurrence in the shifted frames 1 and 2 respectively (frame 0 shifted by one and two nucleotides respectively in the 5′-3′ direction). and are complementary to each other. The subset contains the subset which has the rarity property (6 × 10⁻⁸) to be a complementary maximal circular code with two permutated maximal circular codes and in the frames 1 and 2 respectively. is called a C³ code. A quantitative study of these three subsets in the three frames 0, 1, 2 of protein genes, and the 5′ and 3′ regions of eukaryotes, shows that their occurrence frequencies are constant functions of the trinucleotide positions in the sequences. The frequencies of in the frame 0 of protein genes are 49, 28.5 and 22.5% respectively. In contrast, the frequencies of in the 5′ and 3′ regions of eukaryotes, are independent of the frame. Indeed, the frequency of in the three frames of 5′ (respectively 3′) regions is equal to 35.5% (respectively 38%) and is greater than the frequencies and , both equal to 32.25% (respectively 31%) in the three frames. Several frequency asymmetries unexpectedly observed (e.g. the frequency difference between and in the frame 0), are related to a new property of the subset involving substitutions. An evolutionary analytical model at three parameters (p, q, t) based on an independent mixing of the 22 codons (trinucleotides in frame 0) of with equiprobability (1/22) followed by t ≈ 4 substitutions per codon according to the proportions p ≈ 0.1; q ≈ 0.1 and r = 1 − p − q ≈ 0.8 in the three codon sites respectively, retrieves the frequencies of observed in the three frames of protein genes and explains these asymmetries. Furthermore, the same model (0.1, 0.1, t) after t ≈ 22 substitutions per codon, retrieves the statistical properties observed in the three frames of the 5′ and 3′ regions. The complex behaviour of these analytical curves is totally unexpected and a priori difficult to imagine.     

'Genome order index' should not be used for defining compositional constraints in nucleotide sequences - a case study of the Z-curve

Background  

The Z-curve is a three dimensional representation of DNA sequences proposed over a decade ago and has been extensively applied to sequence segmentation, horizontal gene transfer detection, and sequence analysis. Based on the Z-curve, a "genome order index," was proposed, which is defined as S = a ²+ c ²+t ²+g ², where a, c, t, and g are the nucleotide frequencies of A, C, T, and G, respectively. This index was found to be smaller than 1/3 for almost all tested genomes, which was taken as support for the existence of a constraint on genome composition. A geometric explanation for this constraint has been suggested. Each genome was represented by a point P whose distance from the four faces of a regular tetrahedron was given by the frequencies a, c, t, and g. They claimed that an inscribed sphere of radius r = 1/ contains almost all points corresponding to various genomes, implying that S <r ². The distribution of the points P obtained by S was studied using the Z-curve.     

Conformational analysis of cyclic dipeptides: I. , , , and

Conformational analysis of two pairs of synthetic cyclodipeptides formed by interaction of both side chain functional groups ( , and ) and of the main and side chains ( , and ) was achieved by the method of molecular mechanics. The energetically optimal conformational states of the molecules under study were determined. It was shown that the conformational motility of cyclic system of the compounds under study depends on the relative arrangement of the amide groups and the number of atoms in the cycle.     

High groundwater nitrate concentrations inhibit eutrophication of sulphate-rich freshwater wetlands

During the last 60 years, pollution of the groundwater with has greatly increased in many parts of Europe, as a consequence of excessive use of manure and synthetic fertilisers. Monitoring of groundwater-fed wetlands indicated that sediments with high concentrations had the lowest Fe and concentrations in the pore water. A comparison of two restored open water fens, differing in supply via the groundwater, indicated that the redox potential and the sulphate ( ) reduction rate were lower when the groundwater contained not only but also high concentrations. The lower reduction rates in the -rich open water fen were associated with lower concentrations and the presence of plant species characteristic of clear water. In contrast, the higher reduction rates in the -poor open water fen were associated with very high concentrations and massive development of plant species characteristic of eutrophic environments. Investigations at -rich seepage sites in black alder carrs, showed that high concentrations in the pore water caused chlorosis in the alder carr vegetation, due to lower availability of Fe in the pore water and less Fe uptake by the plants. Experimental desiccation of sediments proved that the -rich seepage sites contained no oxidisable FeS _x , contrary to -poor locations, which became acidified and mobilised extremely high amounts of due to FeS _x oxidation. A laboratory experiment showed that addition to sediments led to reduced releases of Fe, and S^2–, very likely due to the oxidation of reduced Fe and S compounds. Overall, the results confirmed that is an energetically more favourable electron acceptor in anaerobic sediments than Fe and , and that high loads function as a redox buffer, preventing reduction of Fe and . Limited reduction prevents S^2– -mediated mobilisation of from Fe- complexes. At a higher redox potential, reduced Fe, including FeS _x , was oxidised, increasing the content of Fe(III) capable of binding . This prevented increased availability and the concomitant massive development of plant species characteristic of eutrophic environments.     

Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of various benzodiazepine analogues of γ-secretase inhibitors

A 3D QSAR analysis has been performed on a series of 67 benzodiazepine analogues reported as γ-secretase inhibitors using molecular field analysis (MFA), with G/PLS to predict steric and electrostatic molecular field interaction for the activity. The MFA study was carried out using a training set of 54 compounds. The predictive ability of model developed was assessed using a test set of 13 compounds ( as high as 0.729). The analyzed MFA model has demonstrated a good fit, having r² value of 0.858 and cross validated coefficient, value as 0.790. The analysis of the best MFA model provided insight into possible modification of the molecules for better activity.   Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The humus layer determines SO<Subscript>4</Subscript><Superscript>2−</Superscript> isotope values in the mineral soil

Biocycling of sulfur (S) has been proposed to play an important role in the recovery of ecosystems following anthropogenic S deposition. Here, we investigated the importance of the humus layer in the biocycling of S in three forested catchments in the Gårdsjön area of southwestern Sweden with differing S inputs and S isotope signature values. These experimental sites consisted of two reference catchments and the Gårdsjön roof experiment catchment (G1), where anthropogenic deposition was intercepted from 1991 until May 2002 by a roof placed over the entire catchment area. Under the roof, controlled levels of deposition were applied, using a sprinkler system, and the only form of S added was marine SO₄²⁻ with a δ of +19.5‰.We installed ion exchange resin bags at the interface between the humus layer and mineral soil at each of the catchments to collect SO₄²⁻ passing through the humus. The resin bags were installed on four occasions, in 1999 and 2000, covering two summer and two winter periods. The ions collected by each bag during these sampling periods were then eluted and their δ values and SO₄²⁻ concentrations determined. The most striking result is that the average δ value in the resin bags was more than 12‰ lower compared to that of the sprinkler water in the G1 roof catchment. There was no increasing trend in the isotope value in the resin bag SO₄²⁻ despite that the roof treatment has been on-going for almost 10 years; the average value for all resin bags was +7.1‰. The highest δ values found in the G1 roof catchment were between +11‰ and +12‰. However, these values were all obtained from resin bags installed at a single sampling location. Throughfall and resin bag δ values were more similar in the two reference catchments: about +7.5‰ in both cases. There was, however, an increase in resin bag δ values during the first winter period, from about +7‰ to +9‰. The resin bag δ value was linearly and positively related (r² = 0.26, p < 0.001) to the amount of SO₄²⁻ extracted from the resin bags, if relatively high amounts (>50 mmol m⁻²) were excluded. High amounts of resin bag SO₄²⁻ seemed to be related to groundwater inputs, as indicated by the δ value. Our results suggest that rapid immobilization of SO₄²⁻ into a large organic S pool may alter the S isotope value and affect the δ values measured in the mineral soil and runoff.     

Characterization of microbial NO production, N2O production and CH4 oxidation initiated by aeration of anoxic rice field soil

Intermittent drainage of rice fields isdiscussed as an option to mitigate emission ofCH₄, an important greenhouse gas. HoweverN₂O, a potentially more effective greenhouse gas,may be emitted during the aeration phase. Therefore,the metabolism of NO, N₂O, NH ,NO and NO and the kinetics ofCH₄ oxidation were measured after aeration ofmethanogenic rice field soil. Before aeration, thesoil contained NH in relatively highconcentrations (about 4 mM), while NO andNO were almost undetectable. Immediatelyafter aeration both NO and N₂O were produced withrates of about 15 pmol h^-1 gdw^-1 and 5 pmolh^-1 gdw^-1, respectively. Simultaneously,NH decreased while NO accumulated. Later on, NO was depletedwhile NO concentrations increased.Characteristic phases of nitrogen turnover wereassociated with the activities of ammonium oxidizers,nitrite oxidizers and denitrifiers. Oxidation ofNH and production of NO and N₂O wereinhibited by 10 Pa acetylene demonstrating thatnitrification was obligatory for the initiation ofnitrogen turnover and production of NO and N₂O.Ammonium oxidation was not limited by the availableNH and thus, concomittant production of NOand N₂O was not stimulated by addition ofNH . However, addition of NO stimulated production of NO and N₂O in bothanoxic and aerated rice soil slurries. In this case,10 Pa acetylene did not inhibit the production of NOand N₂O demonstrating that it was due todenitrification which was obviously limited by theavailability of NO . In the aerated soilslurries CH₄ was only oxidized if present atelevated concentrations >50 ppmv CH₄). Atatmospheric CH₄ concentrations (1.7 ppmv)CH₄ was not consumed, but was even slightly produced.CH₄ oxidation activity increased afterpreincubation at 20% CH₄, and then CH₄was also oxidized at atmospheric concentrations. CH₄oxidation kinetics exhibited sigmoid characteristicsat low CH₄ concentrations presumably because ofinhibition of CH₄ oxidation by NH .     

Probing the interactions of the solvated electron with DNA by molecular dynamics simulations: bromodeoxyuridine substituted DNA

Solvated electrons () are produced during water radiolysis and can interact with biological substrates, including DNA. To augment DNA damage, radiosensitizers such as bromo-deoxyuridine (BUdR), often referred to as an “electron affinic radiosensitizer”, are incorporated in place of isosteric thymidine. However, little is known about the primary interactions of with DNA. In the present study we addressed this problem by applying molecular modeling and molecular dynamics (MD) simulations to a system of normal (BUdR·A)-DNA and a hydrated electron, where the excess electron was modeled as a localized (H₂O)₆ anionic cluster. Our goals were to evaluate the suitability of the MD simulations for this application; to characterize the motion of around DNA (e.g., diffusion coefficients); to identify and describe configurational states of close localization to DNA; and to evaluate the structural dynamics of DNA in the presence of . The results indicate that has distinct space-preferences for forming close contacts with DNA and is more likely to interact directly with nucleotides other than BUdR. Several classes of DNA - contact sites, all within the major groove, were distinguished depending on the structure of the intermediate water layer H-bonding pattern (or its absence, i.e., a direct H-bonding of with DNA bases). Large-scale structural perturbations were identified during and after the approached the DNA from the major groove side, coupled with deeper penetration of sodium counterions in the minor groove. Figure A rare configuration showing direct interaction between the solvated electron and DNA, where (yellow) and N7(A16) are H-bonded. The close approach from the major groove side invokes deep Na⁺ (magenta) penetration into the minor DNA groove (Fig. 7a).     

A Novel Approach for Bioremediation of a Coastal Marine Wastewater Effluent Based on Artificial Microbial Mats

A novel alternative for wastewater effluent bioremediation was developed using constructed microbial mats on low-density polyester. This biotechnology showed high removal efficiencies for nitrogen and phosphorous in a short retention time (48 h): 94% for orthophosphate (7.78 g m³ d⁻¹), 79% for ammonium (11.30 g m⁻³ d⁻¹), 78% for nitrite (7.46 g m⁻³ d⁻¹), and 83% for nitrate (8.55 g m⁻³ d⁻¹). The microbial mats were dominated by Cyanobacteria genera such as Chroococcus sp., Lyngbya sp., and bacteria of the subclass Proteobacteria representative of the Eubacteria Domain. Nitzschia sp. was the dominant Eukaryote Domain. Various N and P substrates in the wastewater permit the growth of self-forming and self-sustaining bacterial, microalgal, and cyanobacterial communities on a polyester support. The result is the continuous, self-sufficient growth of microbial mats. This is an innovative, economical, and environmentally safe alternative for the treatment of wastewater effluents in coastal marine environments.