Novel evolutionary models and applications to sequence alignment problems

In this paper, we present a novel graph-theoretical approach for representing a wide variety of sequence analysis problems within a single model. The model allows incorporation of the operations “insertion”, “deletion”, and “substitution”, and various parameters such as relative distances and weights. Conceptually, we refer the problem as the minimum weight common mutated sequence (MWCMS) problem. The MWCMS model has many applications including multiple sequence alignment problem, the phylogenetic analysis, the DNA sequencing problem, and sequence comparison problem, which encompass a core set of very difficult problems in computational biology. Thus the model presented in this paper lays out a mathematical modeling framework that allows one to investigate theoretical and computational issues, and to forge new advances for these distinct, but related problems. Through the introduction of supernodes, and the multi-layer supergraph, we proved that MWCMS is -complete. Furthermore, it was shown that a conflict graph derived from the multi-layer supergraph has the property that a solution to the associated node-packing problem of the conflict graph corresponds to a solution of the MWCMS problem. In this case, we proved that when the number of input sequences is a constant, MWCMS is polynomial-time solvable. We also demonstrated that some well-known combinatorial problems can be viewed as special cases of the MWCMS problem. In particular, we presented theoretical results implied by the MWCMS theory for the minimum weight supersequence problem, the minimum weight superstring problem, and the longest common subsequence problem. Two integer programming formulations were presented and a simple yet elegant decomposition heuristic was introduced. The integer programming instances have proven to be computationally intensive. Consequently, research involving simultaneous column and row generation and parallel computing will be explored. The heuristic algorithm, introduced herein for multiple sequence alignment, overcomes the order-dependent drawbacks of many of the existing algorithms, and is capable of returning good sequence alignments within reasonable computational time. It is able to return the optimal alignment for multiple sequences of length less than 1500 base pairs within 30 minutes. Its algorithmic decomposition nature lends itself naturally for parallel distributed computing, and we continue to explore its flexibility and scalability in a massive parallel environment.     

New results on the melting thermodynamics of a circular DNA chain

We investigate the impact of supercoil period and nonzero supercoil formation energy on the thermal denaturation of a circular DNA. Our analysis is based on a recently proposed generalization of the Poland-Scheraga model that allows the DNA melting to be studied for plasmids with circular topology, where denaturation is accompanied by formation of supercoils. We find that the previously obtained first-order melting transition persists under the generalization discussed. The dependence of the size of the order-parameter jump at the transition point and the associated melting temperature are obtained analytically.     

Dynamics of DNA: Experimental controversies and theoretical insights

Recent experimental advances using fluorescence correlation spectroscopy (FCS) have given unprecedented information about the small-scale kinetics of large biopolymers in solution. However some of the first studies in this direction yielded conflicting results for the mean squared displacement of the tagged end-point of a DNA chain, deviating from traditional theories of polymer dynamics. Spurred by this controversy, we have developed a hydrodynamic mean-field theory for single semiflexible polymers which points to a resolution of the differing experimental observations. The theory precisely captures, without fitting parameters, one set of recent FCS results, reproducing the experimental dynamics over five decades in time and three decades of chain lengths. The success of the theory makes it an excellent candidate for a variety of biophysical contexts where the internal fluctuations of semiflexible polymers play a role.     

用单分子技术研究抗癌药物顺铂对DNA结构的影响

文章作者用磁镊与原子力显微镜研究了抗癌药物顺铂对单个DNA分子结构的影响.当顺铂浓度较低时,DNA链变得柔软,驻留长度从～52 nm显著缩短到～15 nm;当顺铂浓度较高时,DNA表现出凝聚现象.基于单分子拉伸和原子力显微镜(AFM)成像两方面的实验结果,文章作者提出一个顺铂导致的DNA变软(softening)-成环(looping)-缩短(shortening)-凝聚(condensing)模型(简写为SLSC模型)来解释观察到的DNA凝聚,并认为通过远程交联使DNA形成小环结构是铂类抗癌药物作用的重要特征.