One nanosecond molecular dynamics simulation of the N-terminal domain of the lambda repressor protein

We have carried out molecular dynamics simulation of the N-terminal domain of the lambda repressor protein in a surrounding environment including explicit waters and ions. We observe two apparent dynamics substates in the nanosecond protein simulation, the transition occurring around 500 ps. The existence of these two apparent substates results from a high flexibility of the arm in each monomer, a relative flexibility of both arms with respect to each other, and a relative displacement of the recognition helices from 30 to 40 A of interhelical distance. Many amino acid residues, including those involved in DNA recognition, undergo a simultaneous transition in their side-chain conformations, consistent with the relationship between side-chain conformation and secondary structural elements, as observed in protein crystal structures. This result suggests plausible conformational changes experienced by the protein upon DNA binding. On the whole, the non-consensus monomer appears to be more flexible than its consensus counterpart.