Features of the interactions between the methyl-CpG motif and the arginine residues on the surface of MBD proteins

Recognition of the methylated regions of the DNA plays an important role in the epigenetic processes. We analyze the interactions between the methylated DNA and the methyl-CpG-binding proteins using two models. The first model was built from a methylated or non-methylated cytosine, a guanine and an arginine residue in the experimental arrangement. We applied the M06L density functional method with a small, polarized double-ζ basis set for the geometry optimizations, and the MP2 method with polarized triple-ζ basis set for the energy calculations. The second model was built from two methylcytosines, guanines, guanidinium groups plus an additional carboxyl group in the experimental arrangement. We applied the B3LYP method with a small, polarized double-ζ basis set for the geometry optimizations and thermal corrections. The single-point energies were obtained from dual-hybrid dRPA75 and dRPA@PBE0 calculations supplemented by a moderately large polarized triple-ζ basis set. The hydration effects were modeled by adding explicit water molecules. These calculations revealed that the hydrophobic interaction has the largest contribution to the Gibbs interaction energy and turns the arginine side chains into hydrogen bonding position. Our results show that the translation of the protein along the DNA double helix is sterically hindered by the contact of its arginine side chains with the methyl groups of the methyl cytosines. This supports a hopping mechanism for the searching movement of the protein along the DNA.