Structural dynamics of the lac repressor-DNA complex revealed by a multiscale simulation.

A multiscale simulation of a complex between the lac repressor protein (LacI) and a 107-bp-long DNA segment is reported. The complex between the repressor and two operator DNA segments is described by all-atom molecular dynamics; the size of the simulated system comprises either 226,000 or 314,000 atoms. The DNA loop connecting the operators is modeled as a continuous elastic ribbon, described mathematically by the nonlinear Kirchhoff differential equations with boundary conditions obtained from the coordinates of the terminal base pairs of each operator. The forces stemming from the looped DNA are included in the molecular dynamics simulations; the loop structure and the forces are continuously recomputed because the protein motions during the simulations shift the operators and the presumed termini of the loop. The simulations reveal the structural dynamics of the LacI-DNA complex in unprecedented detail. The multiple domains of LacI exhibit remarkable structural stability during the simulation, moving much like rigid bodies. LacI is shown to absorb the strain from the looped DNA mainly through its mobile DNA-binding head groups. Even with large fluctuating forces applied, the head groups tilt strongly and keep their grip on the operator DNA, while the remainder of the protein retains its V-shaped structure. A simulated opening of the cleft of LacI by 500-pN forces revealed the interactions responsible for locking LacI in the V-conformation.

[1]  P. Hagerman Flexibility of DNA. , 1988, Annual review of biophysics and biophysical chemistry.

[2]  B. Müller-Hill,et al.  The three operators of the lac operon cooperate in repression. , 1990, The EMBO journal.

[3]  M. Karplus,et al.  Molecular dynamics simulations in biology , 1990, Nature.

[4]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[5]  T. Steitz,et al.  Crystal structure of lac repressor core tetramer and its implications for DNA looping. , 1995, Science.

[6]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[7]  T Schlick,et al.  Modeling superhelical DNA: recent analytical and dynamic approaches. , 1995, Current opinion in structural biology.

[8]  A. Bensimon,et al.  The Elasticity of a Single Supercoiled DNA Molecule , 1996, Science.

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

[10]  W. Olson,et al.  Simulating DNA at low resolution. , 1996, Current opinion in structural biology.

[11]  G. Chang,et al.  Crystal Structure of the Lactose Operon Repressor and Its Complexes with DNA and Inducer , 1996, Science.

[12]  M. Rief,et al.  Reversible unfolding of individual titin immunoglobulin domains by AFM. , 1997, Science.

[13]  Christian Holm,et al.  How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines , 1998 .

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

[15]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[16]  Klaus Schulten,et al.  Elastic Rod Model of a DNA Loop in the Lac Operon , 1999 .

[17]  C A Spronk,et al.  The solution structure of Lac repressor headpiece 62 complexed to a symmetrical lac operator. , 1999, Structure.

[18]  M. Lewis,et al.  A closer view of the conformation of the Lac repressor bound to operator , 2000, Nature Structural Biology.

[19]  K. Schulten,et al.  Steered molecular dynamics and mechanical functions of proteins. , 2001, Current opinion in structural biology.

[20]  Gregory A. Voth,et al.  Interfacing continuum and molecular dynamics: An application to lipid bilayers , 2001 .

[21]  J. Howard,et al.  Mechanics of Motor Proteins and the Cytoskeleton , 2001 .

[22]  Charles L. Brooks,et al.  Identifying native‐like protein structures using physics‐based potentials , 2002, J. Comput. Chem..

[23]  Klaus Schulten,et al.  Identifying unfolding intermediates of FN-III(10) by steered molecular dynamics. , 2002, Journal of molecular biology.

[24]  Klaus Schulten,et al.  Quasicontinuum Representations of Atomic-Scale Mechanics: From Proteins to Dislocations , 2002 .

[25]  V. Pande,et al.  Absolute comparison of simulated and experimental protein-folding dynamics , 2002, Nature.

[26]  R. Cheong,et al.  Fluorescence Resonance Energy Transfer over ∼130 Basepairs in Hyperstable Lac Repressor-DNA Loops , 2003 .

[27]  K. Schulten,et al.  Structural Basis for Cooperative Dna Binding by Cap and Lac Repressor Each Hand Binds with a High Specificity to a 21 Bp Opera- Tor Dna (lewis Et Test Our Predictions and Further Advance the Study Of , 2022 .

[28]  Klaus Schulten,et al.  Multiscale Method for Simulating Protein-DNA Complexes , 2004, Multiscale Model. Simul..

[29]  K. Schulten,et al.  Modelling DNA loops using continuum and statistical mechanics One contribution of 16 to a Theme ‘The mechanics of DNA’ , 2004, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.