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.protein-DNA interaction ͉ molecular dynamics ͉ elastic rod model ͉ DNA loop ͉ large-scale protein motion T he lac repressor (LacI) is a celebrated DNA-binding protein that regulates the function of the lac operon (1), a set of genes responsible for the lactose metabolism in Escherichia coli. In the absence of lactose, LacI inhibits the expression of the operon by binding to two of three specifically recognized DNA sites, called operators, and causing the DNA between the operators to fold into a loop (2-4). Depending on which operators are bound, the loop may have a length of either 384 or 75 bp (4, 5). The smaller loop contains the promoter of the lac operon, which includes binding sites for RNA polymerase and the activator protein CAP (1).Several x-ray structures of LacI, both alone and in complex with DNA, were reported (4, 6, 7). The repressor is a homotetramer folded into a dimer of dimers, two massive ''arms'' connected at the ends by means of a four-helix bundle (Fig. 1). Each arm consists of a core and a DNA-binding head group domain; the lactose binding sites divide the core domains into two subdomains. This architecture is essential for the function of LacI. The protein binds two operators with its two head groups and holds them close together, enforcing the interoperator loop, while leaving the lactose binding sites inside the core domains open for lactose molecules to enter, disrupt the repressor, and induce the expression of the lac operon. At the same time, the loose connections between the LacI domains imply a great degree of structural flexibility that LacI n...