Modular Organization in a Cell: Concepts and Applications

Network biology is conceptualized as an interdisciplinary field, lying at the intersection among graph theory, statistical mechanics and biology. Great efforts have been made to promote the concept of network biology and its various applications in life s

CpG islands or CpG clusters: how to identify functional GC-rich regions in a genome?

Background: CpG islands (CGIs), clusters of CpG dinucleotides in GC-rich regions, are often located in the 5' end of genes and considered gene markers. Hackenberg et al. ( 2006) recently developed a new algorithm, CpGcluster, which uses a completely diffe

Phylogenetic analyses of Fringillidae (Aves : Passeriformes) using mitochondrial tRNA gene sequences and secondary structure

Fringillidae is a large and diverse family of Passeriformes. So far, however, Fringillidae relationships deduced from morphological features and by a number of molecular approaches have remained unproven. Recently, much attention has been attracted to mitochondrial tRNA genes, whose sequence and secondary structural characteristics have shown to be useful for Acrodont Lizards and deep-branch phylogenetic studies. In order to identify useful phylogenetic markers and test Fringillidae relationships, we have sequenced three major clusters of mitochondrial tRNA genes from 15 Fringillidae, taxa. A coincident tree, with coturnix as outgroup, was obtained through Maximum-likelihood method using combined dataset of 11 mitochondrial tRNA gene sequences. The result was similar to that through Neighbor-joining but different from Maximum-parsimony methods. Phylogenetic trees constructed with stem-region sequences of 11 genes had many different topologies and lower confidence than with total sequences. On the other hand, some secondary structural characteristics may provide phylogenetic information on relatively short internal branches at under-genus level. In summary, our data indicate that mitochondrial tRNA genes can achieve high confidence on Fringillidae phylogeny at subfamily level, and stem-region sequences may be suitable only at above-family level. Secondary structural characteristics may also be useful to resolve phylogenetic relationship between different genera of Fringillidae with good performance.

A molecular model of braid-like DNA structure

The three-dimensional molecular models of DNA triple helices and triple-stranded brain-like structure were built up by molecular architecture, and their structural features and energy decomposition were examined. The results showed: (i) The base triplet is the element forming braid-like and triple helix DNA; (ii) Under specified conditions, DNA could form the triplet-stranded braid-like structure; (iii) DNA stability of the braid-like structure is less than that of the triple helix structure. (C) 1995 Academic Press Limited.

Comparative studies on parallel and antiparallel duplex and triplex DNA

Parallel strand models for base sequences d(A)(10). d(T)(10), d(AT)(5) . d(TA)(5), d(G(5)C(5)). d(C(5)G(5)), d(GC)(5) . d(CG)(5) and d(CTATAGGGAT). d(GATATCCCTA), where reverse Watson-Crick A-T pairing with two H-bonds and reverse Watson-Crick G-C pairing with one H-bond or with two H-bonds were adopted, and three models of d(T)(14). d(A)(14). d(T)(14) triple helix with different strand orientations were built up by molecular architecture and energy minimization. Comparisons of parallel duplex models with their corresponding B-DNA models and comparisons among the three triple helices showed: (i) conformational energies of parallel AT duplex models were a little lower, while for GC duplex models they were about 8% higher than that of their corresponding B-DNA models; (ii) the energy differences between parallel and B-type duplex models and among the three triple helices arose mainly from base stacking energies, especially for GC base pairing; (iii) the parallel duplexes with one H-bond G-C pairs were less stable than those with two H-bonds G-C pairs. The present paper includes a brief discussion about the effect of base stacking and base sequences on DNA conformations. (C) 1997 Academic Press Limited.

SNPNB: analyzing neighboring-nucleotide biases on single nucleotide polymorphisms (SNPs)

SNPNB is a user-friendly and platform-independent application for analyzing Single Nucleotide Polymorphism NeighBoring sequence context and nucleotide bias patterns, and subsequently evaluating the effective SNP size for the bias patterns observed from the whole data. It was implemented by Java and Perl. SNPNB can efficiently handle genome-wide or chromosome-wide SNP data analysis in a PC or a workstation. It provides visualizations of the bias patterns for SNPs or each type of SNPs.

Quantitative structure-property relationships for octanol-air partition coefficients of polychlorinated naphthalenes, chlorobenzenes and p,p '-DDT

The octanol-air partition coefficient (K-OA) is a key descriptor of chemicals partitioning between the atmosphere and environmental organic phases. Quantitative structure-property relationships (QSPR) are necessary to model and predict KOA from molecular structures. Based on 12 quantum chemical descriptors computed by the PM3 Hamiltonian, using partial least squares (PLS) analysis, a QSPR model for logarithms of K-OA to base 10 (log K-OA) for polychlorinated naphthalenes (PCNs), chlorobenzenes and p,p'-DDT was obtained. The cross-validated Q(cum)(2) value of the model is 0.973, indicating a good predictive ability of the model. The main factors governing log K-OA of the PCNs, chlorobenzenes, and p,p'-DDT are, in order of decreasing importance, molecular size and molecular ability of donating/accepting electrons to participate in intermolecular interactions. The intermolecular dispersive interactions play a leading role in governing log K-OA. The more chlorines in PCN and chlorobenzene molecules, the greater the log K-OA values. Increasing E-LUMO (the energy of the lowest unoccupied molecular orbital) of the molecules leads to decreasing log K-OA values, implying possible intermolecular interactions between the molecules under study and octanol molecules. (C) 2002 Elsevier Science Ltd. All rights reserved.

Substitution Matrices of Residue Triplets Derived from Protein Blocks.

In protein sequence alignment, residue similarity is usually evaluated by substitution matrix, which scores all possible exchanges of one amino acid with another. Several matrices are widely used in sequence alignment, including PAM matrices derived from homologous sequence and BLOSUM matrices derived from aligned segments of BLOCKS. However, most matrices have not addressed the high-order residue-residue interactions that are vital to the bioproperties of protein.With consideration for the inherent correlation in residue triplet, we present a new scoring scheme for sequence alignment. Protein sequence is treated as overlapping and successive 3-residue segments. Two edge residues of a triplet are clustered into hydrophobic or polar categories, respectively. Protein sequence is then rewritten into triplet sequence with 2 · 20 · 2 = 80 alphabets. Using a traditional approach, we construct a new scoring scheme named TLESUMhp (TripLEt SUbstitution Matrices with hydropobic and polar information) for pairwise substitution of triplets, which characterizes the similarity of residue triplets. The applications of this matrix led to marked improvements in multiple sequence alignment and in searching structurally alike residue segments. The reason for the occurrence of the ‘‘twilight zone,’’ i.e., structure explosion of lowidentity sequences, is also discussed.

Potential energy landscape and robustness of a gene regulatory network: Toggle switch

Finding a multidimensional potential landscape is the key for addressing important global issues, such as the robustness of cellular networks. We have uncovered the underlying potential energy landscape of a simple gene regulatory network: a toggle switch. This was realized by explicitly constructing the steady state probability of the gene switch in the protein concentration space in the presence of the intrinsic statistical fluctuations due to the small number of proteins in the cell. We explored the global phase space for the system. We found that the protein synthesis rate and the unbinding rate of proteins to the gene were small relative to the protein degradation rate; the gene switch is monostable with only one stable basin of attraction. When both the protein synthesis rate and the unbinding rate of proteins to the gene are large compared with the protein degradation rate, two global basins of attraction emerge for a toggle switch. These basins correspond to the biologically stable functional states. The potential energy barrier between the two basins determines the time scale of conversion from one to the other. We found as the protein synthesis rate and protein unbinding rate to the gene relative to the protein degradation rate became larger, the potential energy barrier became larger. This also corresponded to systems with less noise or the fluctuations on the protein numbers.

Single-molecule dynamics reveals cooperative binding-folding in protein recognition

The study of associations between two biomolecules is the key to understanding molecular function and recognition. Molecular function is often thought to be determined by underlying structures. Here, combining a single-molecule study of protein binding with an energy-landscape-inspired microscopic model, we found strong evidence that biomolecular recognition is determined by flexibilities in addition to structures. Our model is based on coarse-grained molecular dynamics on the residue level with the energy function biased toward the native binding structure ( the Go model). With our model, the underlying free-energy landscape of the binding can be explored. There are two distinct conformational states at the free-energy minimum, one with partial folding of CBD itself and significant interface binding of CBD to Cdc42, and the other with native folding of CBD itself and native interface binding of CBD to Cdc42. This shows that the binding process proceeds with a significant interface binding of CBD with Cdc42 first, without a complete folding of CBD itself, and that binding and folding are then coupled to reach the native binding state.

Structural analysis of SNARE motifs from sea perch, Lateolabrax japonicus by computerized approaches

Three cDNA sequences encoding four SNARE (N-ethylmaleimide-sensitive fusion protein attachment protein receptors) motifs were cloned from sea perch, and the deduced peptide sequences were analyzed for structural prediction by using 14 different web servers and softwares. The "ionic layer" structure, the three dimensional extension and conformational characters of the SNARE 7S core complex by using bioinformatics approaches were compared respectively with those from mammalian X-ray crystallographic investigations. The result suggested that the formation and stabilization of fish SNARE core complex might be driven by hydrophobic association, hydrogen bond among R group of core amino acids and electrostatic attraction at molecular level. This revealed that the SNARE proteins interaction of the fish may share the same molecular mechanism with that of mammal, indicating the universality and solidity of SNARE core complex theory. This work is also an attempt to get the protein 3D structural information which appears to be similar to that obtained through X-ray crystallography, only by using computerized approaches. (C) 2007 Elsevier Ltd. All rights reserved.

A mechanistic model of algal photoinhibition induced by photodamage to Photosystem-II

Photoinhibition is a central problem for the understanding of plasticity in photosynthesis vs. irradiance response. It effectively reduces the photosynthetic rate. In this contribution, we present a mechanistic model of algal photoinhibition induced by photodamage to photosystem-II. Photosystem-IIs (PSIIs) are assumed to exist in three states: open, closed and inhibited. Photosynthesis is closely associated with the transitions between the three states. The present model is defined by four parameters: effective cross section of PSII, number of PSIIs, turnover time of electron transfer chains and the ratio of rate constant of damage to that of repair of D1 proteins in PSIIs. It gives a photosynthetic response curve of phytoplankton to irradiance (PI-curve). Without photoinhibition, the PI-curve is in hyperbola with the first three parameters. The PI-curve with photoinhibition can be simplified to the same form as the hyperbola by replacing either the number of PSIIs with the number of functional PSIIs or the turnover time of electron transfer chains with the average turnover time.