Flexibility of Dicer Studied by Implicit Solvent Molecular Dynamics Simulations

Dicer is an enzyme responsible for processing double-stranded RNAs and plays a key role in an RNAi mechanism. Structural insight into the Dicer is provided by the crystal structure of eukaryotic Dicer from Giardia intestinalis. It has been proposed that the structure has three structurally rigid regions that are connected by the flexible hinges. Flexibility of the Dicer is believed to be a critical feature for its function. Spatial arrangement of the RNA-recognition and the catalytic regions is crucial for producing small RNAs of defined length. It has been suggested that in the Giardia Dicer a Platform domain may help in specific arrangement of these regions. To learn more about the role of the Platform domain in Giardia Dicer, we have performed molecular dynamics (MD) simulations of the whole Dicer (WT Dicer) and the Dicer with a deleted platform domain (ΔPlf Dicer). The MD simulations were carried out in an implicit solvent model with two implementations of analytic Generalized Born (GB) solvation model in CHARMM: GBMV (Generalized Born using Molecular Volume) and GBSW (Generalized Born with simple Switching). To detect the key global motions of the Dicer, a principal component analysis (PCA) of the obtained MD trajectories has been used. To further explore the motion of the Dicer, we performed a domain motion analysis with the DYNDOM program. The simulations show that both WT Dicer and ΔPlf Dicer display flexibility which can be described as a movement of two or three domains. The removal of the Platform substantially changed the flexibility and arrangement of these domains. During the MD simulations of ΔPlf Dicer an large movement of the RNA-recognition domain was observed.

[1]  I. MacRae,et al.  Structure of Dicer and mechanistic implications for RNAi. , 2006, Cold Spring Harbor symposia on quantitative biology.

[2]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[3]  Stephen Neidle,et al.  Molecular dynamics and principal components analysis of human telomeric quadruplex multimers. , 2008, Biophysical journal.

[4]  I. MacRae,et al.  Structural determinants of RNA recognition and cleavage by Dicer , 2007, Nature Structural &Molecular Biology.

[5]  C. Brooks,et al.  Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field. , 2006, Journal of the American Chemical Society.

[6]  Angela N. Brooks,et al.  Structural Basis for Double-Stranded RNA Processing by Dicer , 2006, Science.

[7]  H. Berendsen,et al.  Systematic analysis of domain motions in proteins from conformational change: New results on citrate synthase and T4 lysozyme , 1998, Proteins.

[8]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[9]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[10]  D. Case,et al.  Modeling of flap endonuclease interactions with DNA substrate. , 2003, Journal of molecular biology.

[11]  John Mongan,et al.  Interactive essential dynamics , 2004, J. Comput. Aided Mol. Des..

[12]  Charles L. Brooks,et al.  Generalized born model with a simple smoothing function , 2003, J. Comput. Chem..

[13]  D. Beglov,et al.  Atomic Radii for Continuum Electrostatics Calculations Based on Molecular Dynamics Free Energy Simulations , 1997 .

[14]  Alexander D. MacKerell,et al.  Improved treatment of the protein backbone in empirical force fields. , 2004, Journal of the American Chemical Society.

[15]  Viktor Hornak,et al.  HIV-1 protease flaps spontaneously open and reclose in molecular dynamics simulations. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[16]  J. Chocholousová,et al.  Implicit solvent simulations of DNA and DNA-protein complexes: agreement with explicit solvent vs experiment. , 2006, The journal of physical chemistry. B.

[17]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[18]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[19]  Charles L. Brooks,et al.  New analytic approximation to the standard molecular volume definition and its application to generalized Born calculations , 2003, J. Comput. Chem..

[20]  Mensur Dlakic,et al.  DUF283 domain of Dicer proteins has a double-stranded RNA-binding fold , 2006, Bioinform..

[21]  H. Berendsen,et al.  Model‐free methods of analyzing domain motions in proteins from simulation: A comparison of normal mode analysis and molecular dynamics simulation of lysozyme , 1997, Proteins.

[22]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[23]  C. Brooks,et al.  Novel generalized Born methods , 2002 .