摘要
在这篇文章中,我们构造了一种Level Set模型对生物大分子的分子场进行特征提取和分析。对于3维大分子的分子场上的体数据,我们定义了一个新颖的跟踪几何活动轮廓线的变分方程,这个方程包含两部分:外部能量和内部能量。外部能量用来控制活动轮廓线的移动,而内部能量则作为Level Set函数的惩罚项,以避免该模型在迭代过程中周期性地重新初始化。我们不断改变方程中相关的参数,分析在整个演化过程中这些参数对最后结果的影响。与传统的LevelSet方法相比,我们的方法在分子场上进行特征提取更加精确和有效。
采用上述模型,我们对静态和动态的分子场数据进行分析。对于静态分子场,选取能量跃变区域来揭示分子场特征,同时构造体特征函数来分析LevelSet曲面的拓扑特征区域和整体分子场的能量分布特点。对于动态分子场,我们考察同一蛋白质分子在不同运动状态下的数据,例如观察HIV-1蛋白酶间接受水分子作用改变构象、形成中间通道的整个过程。与此同时,对不同的运动状态,计算分子场多属性Level Set曲面,根据属性的相似程度,寻找分子运动过程中Level Set曲面显著变化的关键时刻。通过上述特征分析方法,结合多种面以及体可视化绘制技术,本文成功地揭示了蛋白质分子场的重要特征。
初步实验表明,本文的方法计算出来的特征与已有生物学结论一致,如HIV-1蛋白酶内部的水分子排除通道,DPS蛋白质与铁离子结合的空腔区域,以及HIV-1蛋白酶在水溶液中进行SMD模拟时的生物活性变化等。生物学家认为这些结果可能对今后研究HIV-1蛋白酶的亲水、疏水性有重要的价值,并可作为“蛋白质-蛋白质”相互作用或“蛋白质-配体”相互作用中首要参考的活性位点,具有非常巨大的研究意义。
In this paper, we exploit a Level Set method to extract and analyze the feature of macromolecular field. For the 3D volume data of the field, we propose a novel variational formulation for geometric active contour, which is composed of an external energy to drive the active contours' motion, and an internal energy which penalizes the Level Set function adaptively to avoid periodically re-initializing. We adjust the parameters of each component and investigate their influence in the whole process. According to the comparison with the traditional frameworks, our approach is shown to be more accurate and efficient on the feature extraction.
Based on the Level Set method constructed in the first part, we give a further analysis for the static and dynamic molecular field data. For the static molecular field, we choose the region where the energy has the quantum leap and analyze the molecular field. Then, we construct the volume feature function to extract the topological region and get the energy distribution of the whole molecular field. For the dynamic molecular field, we consider the same protein in the different states, like observing the process in changing the conformation and forming the middle tunnel. At the mean time, we calculate the multi-attributes of these states and try to find the crucial time during the movement of the molecular. According to the surface and volume rendering techniques, we reveal the significant feature of the macromolecular field.
Preliminary results show that, we successfully find the escape route of water molecules hidden in the HIV-1 protease, the internal cavity of the DPS protein for the iron atom's entry and deposition and biological activity of HIV-1 protease when processing SMD simulation in the aqueous solution. Both of these results are in accordance with the experimental results given by chemists and biologists. Biologists believe that these results may be useful to study the hydrophilic and hydrophobic of HIV-1 protease in the future, and these active cites can be considered as the most important regions in the protein-protein interaction and protein-ligand interaction.
引文
[1]A.R.Leach.Molecular modelling:principles and applications.Pearson Education EMA.Second edition,January 2001.
[2]C.Gibas and P.Jambeck.Developing bioinformatics computer skills.Science Press,2002.
[3]J.C.Whisstock.Lesk am:Prediction of protein function from protein sequence and structure.Quarterly Reviews of Biophysics,36:307-340,2003.
[4]D.L.Wild and M.A.S.Saqi.Structural proteomics:inferring function from protein structure.Current Proteomics,1:59-65,2004.
[5]C.Zhang and S.H.Kim.Overview of structural genomics:from structure to function.Current Opinion in Chemical Biology,7:28-32,2003.
[6]E.E.Hodgkin and W.G.Richards.Molecular similarity based on electrostatic potential and electric field.International Journal of Quantum Chemistry,pages 105-109,1987.
[7]R.Carbo,L.Leyda,and M.Arnau.How similar is a molecule to another? an electron density measure of similarity between two molecular structures.International Journal of Quantum Chemistry,6:1185-1189,1980.
[8]M.V.Kreveld,R.V.Oostrum,C.Bajaj,D.Schikore,and V.Pascucci.Contour trees and small seed set for isosurface traversal.In Proceedings:Thirteen ACM Symposium on Computational Geometry,pages 212-219,1997.
[9]H.Edelsbrunner,D.Letscher,and A.Zomorodian.Topological persistence and simplification.Discrete and Computational Geometry,28(4):511-533,2001.
[10]F.Cazals,F.Chazal,and T.Lewiner.Molecular shape analysis based upon the morse-smale complex and the connolly function.In ACM Symposium on Computational Geometry,pages 351-360,2003.
[11]Y.J.Chiang and X.Lu.Progressive simplification of tetrahedral meshes preserving all isosurface topologies.Computer Graphics Forum,22(3):493-504,2003.
[12]T.Gerstner and R.Pajarola.Topology preserving and controlled topology simplifying multiresolution isosurface extraction.IEEE Visualization:Proceedings of the conference,pages 259-266,2000.
[13]P.Cignoni,D.Constanza,C.Montani,C.Rocchini,and R.Scopigno.Simplification of tetrahedral meshes with accurate error evaluation.IEEE Hsualization:Proceedings of the conference,pages 85-92,2000.
[14]H.Carr,J.Snoeyink,and M.van de Panne.Simplifying flexible isosurfaces using local geometric measures.IEEE Visualization:Proceedings of the conference,pages 497-504,2004.
[15]I.Guskov and Z.Wood.Topological noise removal.In Proceedings of Graphics Interface,pages 19-26,2001.
[16]S.Takahashi,G.M.Nielson,Y.Takeshima,and I.Fujishiro.Topological volume skeletonization using adaptive tetrahedralization.In Proceedings of Geometric Modeling and Processing,pages 227-236,2004.
[17]S.Takahashi,Y.Takeshima,and I.Fujishiro.Topological volume skeletonization and its application to transfer function design.Graphical Models,66(1):24-49,2004.
[18]Z.Wood,H.Hoppe,M.Desbrun,and P.Schroder.Removing excess topology from isosur-faces.A CM Transactions on Graphics,23(2):190-208,2004.
[19]H.Carr,J.Snoeyink,and U.Axen.Computing contour trees in all dimensions.In Symposium on Discrete Algorithms,pages 918-926,2000.
[20]S.Smale.On gradient dynamical systems.Ann.of Math.,74:199-206,1961.
[21]S.Smale.Generalized poincare's conjecture in dimensions greater than four.The Annals of Mathematics,Second Series,74(2):391-406,1961.
[22]H.Edelsbrunner,J.Harer,and A.Zomorodian.Hierarchical morse-smale complexes for piecewise linear 2-manifolds.Discrete and Computational Geometry,30(1):87-107,2003.
[23]P.T.Bremer,H.Edelsbrunner,B.Hamann,and V.Pascucci.A topological hierarchy for functions on triangulated surfaces.IEEE Transactions on Visualization and Computer Graphics,10(4):385-396,2004.
[24]A.Gyulassy,V.Natarajan,V.Pascucci,P.T.Bremer,and B.Hamann.Topology-based simplification for feature extraction from 3d scalar fields.IEEE Visualization:Proceedings of the conference,pages 535-542,2005.
[25]R.Forman.A user's guide to discrete morse theory.2001.
[26]T.Lewiner,H.Lopes,and G.Tavares.Applications of forman's discrete morse theory to topology visualization and mesh compression.IEEE Transactions on Visualization and Computer Graphics,10(5):499-508,2004.
[27]A.Gyulassy,V.Natarajan,V.Pascucci,and B.Hamann.Efficient computation of morsesmale complexes for three-dimensional scalar functions.IEEE Transactions on Visualization and Computer Graphics,13(6):1440-1447,2007.
[28]O.G.Staadt and M.H.Gross.Progressive tetrahedralizations.In Proceedings of IEEE Conference on Visualization,pages 397-402,1998.
[29]V.Natarajan and H.Edelsbrunner.Simplification of three-dimensional density maps.IEEE Transactions on Visualization and Computer Graphics,10(5):587-597,2004.
[30]M.Garland and Y.Zhou.Quadric-based simplification in any dimension.ACM Transactions on Graphics,24(2):209-239,2005.
[31]S.J.Osher and R.Fedkiw.Level set methods and dynamic implicit surfaces.Applied Mathematical Science.Springer,page 153,2002.
[32]J.A.Sethian.Fast Marching Methods and Level Set Methods:Evolving Interfaces in Computational Geometry,Fluid Mechanics.Cambridge University Press:Cambridge,UK,1999.
[33]T.Cootes,A.Hill,C.J.Taylor,and J.Haslam.The use of active shape models for locating structures in medical images.Information Processing in Medical Imaging(IPMI'93) number 687 in Lecture Notes in Computer Science,pages 33-47,1993.
[34]C.Bajaj,Z.Y.Yu,and M.Auer.Volumetric filtering,feature extraction,and visualization of volumetric tomographic molecular imaging.Journal of Structural Biology,2003.
[35]X.C.Tai and T.F.Chan.A survey on multiple level set methods with applications for identifying piecewise constant functions.International Journal of Numerical Analysis and Modelling,1(1):25-48,2004.
[36]R.Whitaker,D.Breen,K.Museth,and N.Soni.Segmentation of biological volume datasets using a level-set.Volume Graphics,pages 249-263,2001.
[37]D.Magee,A.Bulpitt,and E.Berry.Combining 3D deformable models and level set methods for the segmentation of abdominal aortic aneurysms.Proceedings of the British Machine Vision Conference,pages 333-341,2001.
[38]D.Magee,A.Bulpitt,and E.Berry.Level set methods for the 3D segmentation of CT images of abdominal aortic aneurysms.Understanding and Analysis,2001.
[39]L.Zhukov,K.Museth,D.Breen,R.Whitakery,and A.H.Barr.Level set modeling and segmentation of dt-mri brain data.Journal of Electronic Imaging,12(1):125-133,2003.
[40]T.Can,C.Chen,and Y.Wang.Efficient molecular surface generation using Level-set methods.Journal of Molecular Graphics and Modeling,25:442-454,2006.
[41]J.A.Sethian.Curvature and the evolution of fronts.Communications in Mathematical Physics,101:487-499,1985.
[42]J.A.Sethian.Numerical methods for propagating fronts.Variational Methods for Free Surface Interfaces,Concus P,Finn R(eds),pages 155-164,1987.
[43]S.Osher and J.A.Sethian.Fronts propagating with curvature-dependent speed:algorithms based on hamilton-jacobi formulations.Journal of Computational Physics,79:12-49,1988.
[44]H.K.Zhao,S.Osher,B.Merrian,and M.Kang.Irnpficit and nonparametric shape reconstruction from unorganized data using a variational levelset method.Computer Vision and Image Understnading,80:295-319,2000.
[45]H.K.Zhao,S.Osher,and R.Fedkiw.Fast surface reconstruction using the level set method.1st IEEE Workshop on Variational and Levelset Methods,pages 194-202,2001.
[46]T.Tasdizen,R.Whitaker,P.Burchard,and S.Osher.Geometric surface processing via normal maps.Technical Report UUCS-02-03,University of Utah,2002.
[47]M.Kass,A.Witkin,and D.Terzopoulos.Snakes:Active contour models.International Journal of Computer Vision,1(4):321-331,1987.
[48]V.Caselles,R.Kimmel,G.Sapiro,and C.Sbert.Minimal surfaces based object segmentation.Numerische Mathematik,19(4):394-398,1997.
[49]R.Malladi,J.A.Sethian,and B.C.Vemuri.Evolutionary fronts for topology independent shape modeling and recovery.Proceedings of the third European conference on Computer vision,800:3-13,1994.
[50]J.Shah.A common framework for curve evolution,segmentation and anisotropic diffusion.Proceedings of International Conference on Computer Vision and Pattern Recognition,pages 136-142,1996.
[51]N.K.Paragios and R.Deriche.A pde-based levelset approach for detection and tracking of moving objects.Computer Vision,pages 1139-1145,1998.
[52]T.E Chan and L.A.Vese.Active contours without edges.Image Processing,10:266-277,2001.
[53]A.Lefohn,J.Kniss,C.Hansen,and R.Whitaker.A streaming narrow-band algorithm:Interactive deformation and visualization of levelsets.IEEE Transactions on Visualization and Computer Graphics,10(40):422-433,2004.
[54]R.Whitaker,D.Breen,K.Museth,and N.Soni.Segmentation of biological volume datasets using a levelset framework.Volume Graphics,pages 249-263,2001.
[55]D.Breen,R.Whitaker,K.Museth,and L.Zhukov.Levelset Segmentation of Biological Volume Datasets.J.Suri(ed.).Handbook of Medical Image Analysis,Volume I Segmentation Part A,Kluwer,New York,2005.
[56]H.Yamauchi,S.Gumhold,R.Zayer,and H.P.Seidel.Mesh segmentation driven by gaussian curvature.Source Visual Computer,21(8-10 Special Issue SI):659-668,2005.
[57]Y.Shi and W.C.Karl.Real-time tracking using levelsets.Computer Vision and Pattern Recognition,2:20-25,2005.
[58]D.Enright,S.Marschner,and R.Fedkiw.Animation and rendering of complex water surfaces.SIGGRAPH,ACM TOG 21:736-744,2002.
[59]D.Enright,R.Fedkiw,J.Ferziger,and I.Mitchell.A hybrid particle levelset method for improved interface capturing.Journal of Computational Physics,183(1):83-116,2002.
[60]D.Enright,F.Losasso,and R.Fedkiw.A fast and accurate semi-lagrangian particle levelset method.Computers and Structures,83:479-490,2005.
[61]Y.Shi,T.G.Goktekin,A.W.Bargteil,and J.EO'Brien.A method for animating viscoelastic fluids.ACM Transactions on Graphics,23:463-467,2004.
[62]J.Stare.Flows on surfaces of arbitrary topology.Proceedings of SIGGRAPH,pages 724-731,2003.
[63]V.I.Arnold.Geometrical Methods in the Theory of Ordinary Differential Equations.New York:Springer-Verlag,1983.
[64]C.M.Li,C.Y.Xu,C.F.Gui,and M.D.Fox.Level set evolution without re-initialization a new variational formulation.Computer Vision and Pattern Recognition,1(1):430-436,2005.
[65]K.Siddiqi,Y.B.Lauziere,A.Tannenbaum,and S.W.Zucker.Area and length minimizing flows for shape segmentation.Image Processing,IEEE Transactions,7:433-443,1998.
[66]V.Caselles,F.Catte,T.Coll,and F.Dibos.Geometric model for active contours.Numerische Mathematik,66:1-31,1993.
[67]R.Malladi,J.A.Sethian,and B.C.Vemuri.Shape modeling with front propagation:a level set approach.Transaction on Pattern Analysis and Machine Intelligence,17(2):158-175,1995.
[68]L.Evans.Partial Differential Equations.Providence:American Mathematical Society,1998.
[69]T.Michael.Heath.Scientific Computing,An Introductory Survey.McGraw-Hill.Second edition,2002.
[70]D.Eberly,L.Lancaster,and A.Alyassin.On gray scale image measurements.Graphics Models and Image Processing,56(6):550-562,1991.
[71]R.D.Cramer Ⅲ,D.E.Patterson,and J.D.Bunce.Comparative molecular field analysis (CoMFA).1.effect of shape on binding of steroids to carrier proteins.Journal of the American Chemical Society,110:5959-5967,1988.
[72]W.Han,Z.D.Jin,L.Wang,M.Hu,Z.Y.Wang,and Q.S.Peng.The feature extraction of macromolecular field based on level set method.In Proceedings of the CAD/Graphics,pages 108-113,2007.
[73]L.C.Evans.Partial Differential Equations.Providence:American Mathematical Society,1998.
[74]D.Eberly,L.Lancaster,and A.Alyassin.On gray scale image measurements.Graphics Models and Image Processing,56(6):550-562,1991.
[75]B.G.Turner and M.F.Summers.Structural biology of hiv.Journal of Molecular Biology,285(1):1-32,1999.
[76]R.A.Grant,D.J.Filman,S.E.Finkel,R.Kolter,and J.M.Hogle.The crystal structure of dps,a ferritin homolog that binds and protects dna.Nature Structural Biology,5(4):294-303,1998.
[77]P.Ceci,A.Ilari,E.Falvo,and E.Chiancone.The dps protein of agrobacterium tumefaciens does not bind to dna but protects it toward oxidative cleavage.Journal of Biological Chemistry,278(22):20319-20326,2003.
[78]J.C.Phillips,R.Braun,W.Wang,J.Gumbart,E.Tajkhorshid,E.Villa,C.Chipot,R.D.Skeel,L.Kale,and K.Schulten.Scalable molecular dynamics with namd.Journal of Computational Chemistry,26:1781-1802,2005.
[79]A.D.MacKerell,D.Bashford,and M.Bellott et al.All-atom empirical potential for molecular modeling and dynamics studies of proteins.Journal of Computational Chemistry,102:3586-3616,1998.