next up previous contents
Next: Diatomics-in-Molecules Up: Some Recognised Systems Previous: Binary Lennard-Jones   Contents


BLN Off-Lattice Protein Model

The general three-colour bead protein model is specified by keyword BLN. The potential follows the form described in Proc. Natl. Acad. Sci. USA, 100, 10712, 2003, expect that the coefficients Ai, Bi, Ci and Di include a factor of ε explicitly.

0pt

V   =   $\displaystyle {\frac{{1}}{{2}}}$Kr$\displaystyle \sum_{{i=1}}^{{N-1}}$(Ri, i+1 - Re)2 + $\displaystyle {\frac{{1}}{{2}}}$Kθ$\displaystyle \sum_{i}^{{N-2}}$(θi - θe)2  
    +  ε$\displaystyle \sum_{i}^{{N-3}}$$\displaystyle \Big[$Ai(1 + cos$\displaystyle \varphi_{i}^{}$) + Bi(1 - cos$\displaystyle \varphi_{i}^{}$)  
             + Ci(1 + cos 3$\displaystyle \varphi_{i}^{}$) + Di$\displaystyle \left(\vphantom{1+\cos\left[\varphi_i+\pi/4\right]}\right.$1 + cos$\displaystyle \left[\vphantom{\varphi_i+\pi/4}\right.$$\displaystyle \varphi_{i}^{}$ + π/4$\displaystyle \left.\vphantom{\varphi_i+\pi/4}\right]$$\displaystyle \left.\vphantom{1+\cos\left[\varphi_i+\pi/4\right]}\right)$$\displaystyle \Big]$  
    +  4ε$\displaystyle \sum_{{i=1}}^{{N-2}}$$\displaystyle \sum_{{j=i+2}}^{N}$$\displaystyle \left[\vphantom{S_{12}\left(\frac{\sigma}{R_{ij}}\right)^{\!12}
+S_6\!\left(\frac{\sigma}{R_{ij}}\right)^{\!6}}\right.$S12$\displaystyle \left(\vphantom{\frac{\sigma}{R_{ij}}}\right.$$\displaystyle {\frac{{\sigma}}{{R_{ij}}}}$$\displaystyle \left.\vphantom{\frac{\sigma}{R_{ij}}}\right)^{{\!12}}_{}$ + S6$\displaystyle \left(\vphantom{\frac{\sigma}{R_{ij}}}\right.$$\displaystyle {\frac{{\sigma}}{{R_{ij}}}}$$\displaystyle \left.\vphantom{\frac{\sigma}{R_{ij}}}\right)^{{\!6}}_{}$$\displaystyle \left.\vphantom{S_{12}\left(\frac{\sigma}{R_{ij}}\right)^{\!12}
+S_6\!\left(\frac{\sigma}{R_{ij}}\right)^{\!6}}\right]$, (9)

where Rij is the separation between beads i and j and the units of distance and energy are σ and ε, respectively. The first term represents the bonds linking successive beads in the linear chain, and a value of Kr = 231.2 εσ-2 was used in most of the work on the Honeycutt and Thirumalai frustrated 46-bead model. The second term is a sum over the bond angles, θi, defined by the triplets of atomic positions $\bf R_{i}^{}$ to $\bf R_{{i+2}}^{}$, and values Kθ = 20 ε rad-2 and θe = 105o were used for the 46-bead model. The third term is a sum over the dihedral angles, $\varphi_{i}^{}$, defined by the quartets $\bf R_{i}^{}$ to $\bf R_{{i+3}}^{}$. In the 46-bead model Ai = Ci = 1.2 if the quartet involved no more than one N monomer, generating a preference for the trans conformation ( $\varphi_{i}^{}$ = 180o), whereas if two or three N monomers are involved then Ai = 0 and Ci = 0.2. This choice makes the three neutral segments of the chain flexible and enables them to accommodate turns. A general specification of these parameters is possible in the new BLN framework via the auxiliary file BLNsequence. The last term in ([*]) represents the nonbonded interactions. In the current BLN implementation Re is set equal to σ, i.e. to unity in reduced units.

An appropriate BLNsequence file for the usual 46-bead model contains the following lines:

comment: S12 > 0 and S6 < 0 for B-B, L-L and L-B, N-L and N-B and N-N 1.0D0 -1.0D0 0.33333333333333D0 0.33333333333333D0 1.0D0 0.0D0 comment: coefficients A, B, C, D comment: for Helical, Extended and Turn residues in order, four per line 0.0D0 1.2D0 1.2D0 1.2D0 0.9D0 0.0D0 1.2D0 0.0D0 0.0D0 0.0D0 0.2D0 0.0D0 LBLBLBLBBNNNBBBLBLBBBNNNLLBLLBBLLBNBLBLBLBLNNNLBBLBLBBBL EEEEEETEHTHEEEEEEEEHHEHHHHHHHHHHEHTEEEEEEETTTEEEEEEEE

The penultimate line defines the sequence, and the final line defines which set of Ai, Bi, Ci and Di parameters apply to which parts of the structure.[#!BrownFH03!#]


next up previous contents
Next: Diatomics-in-Molecules Up: Some Recognised Systems Previous: Binary Lennard-Jones   Contents