The enterobacteria lambda phage is a paradigm temperate bacteriophage. Its lysogenic and lytic life cycles echo competition between the DNA binding λ-repressor (CI) and CRO proteins. Here we scrutinize the structure, stability and folding pathways of the λ-repressor protein, that controls the transition from the lysogenic to the lytic state. We first investigate the super-secondary helix-loophelix composition of its backbone. We use a discrete Frenet framing to resolve the backbone spectrum in terms of bond and torsion angles. Instead of four, there appears to be seven individual loops. We model the putative loops using an explicit soliton Ansatz. It is based on the standard soliton profile of the continuum nonlinear Schrödinger equation. The accuracy of the Ansatz far exceeds the B-factor fluctuation distance accuracy of the experimentally determined protein configuration. We then investigate the folding pathways and dynamics of the λ-repressor protein. We introduce a coarse-grained energy function to model the backbone in terms of the Cα atoms and the side-chains in terms of the relative orientation of the C β atoms. We describe the folding dynamics in terms of relaxation dynamics, and find that the folded configuration can be reached from a very generic initial configuration. We conclude that folding is dominated by the temporal ordering of soliton formation. In particular, the third soliton should appear before the first and second. Otherwise, the DNA binding turn does not acquire its correct structure. We confirm the stability of the folded configuration by repeated heating and cooling simulations.
I: INTRODUCTIONThe transition between the lysogenic and the lytic state in bacteriophage λ infected E. coli cell is the paradigm genetic switch mechanism. It is described in numerous molecular biology textbooks and review articles [1]- [7]. The interplay between the lysogeny maintaining λ-repressor (CI) protein and the CRO regulator protein that controls the transition to the lytic state is a simple model for more complex regulatory networks, including those that can lead to cancer in humans.In the present article we describe the physical properties of the λ-repressor protein, that controls the lysogenic-to-lytic transition. We investigate in detail the stability of its native conformation, the dynamics of the folding process, and the landscape of folding pathways. We find that the folded configuration displays a structure which is unique among all known protein structures. We also conclude that the folding pathways are entirely dominated by the loop regions. In particular, the temporal ordering of loop formation appears to be the decisive factor for the protein's ability to reach its native fold. If solitons form in a wrong order the protein may misfold.Full crystallographic information of the experimental λ-repressor structure that we use in our investigation is available in Protein Data Bank (PDB) [8] under the code 1LMB. This structure is a homo-dimer with 92 residues in each of the two monomers. It maintain...