Search

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. [ABSTRACT FROM AUTHOR]

Glycoproteins play a central role in the immune response. In this study, we focus on the core of the common O-linked mucin-type glycopeptides. It has been observed that glycosylation stabilizes the protein in a stiffened, extended structure. We provide a unified picture for the conformation and stabilization of the O-glycosidic linkage using O-(2-acetylamino-2-deoxy-α- or -β- d-galacto- or -mannopyranosyl)- N-acetyl- l-serinamide model structures. We have calculated equilibrium geometries of the model structures with the B3LYP/6-31G(2 df, p) method suggested in the Gaussian-4 theory. According to the relative energies, we confirm that the GalNAc-Ser linkage is more stable than its mannose analogues. The natural preference for the α-GalNAc-Ser over the β-GalNAc-Ser anomers can be explained by entropic effects. We explored the hydrogen bonding patterns on the carbohydrate unit calculating highly accurate dRPA@PBE0.25 and dRPA75 energies and found that in some cases, the acetamido group can be fixed by hydrogen bonding from the (O3)H atom, but in most of the cases, it can rotate more freely. The torsion angles in the glycosidic linkage show that the linkage is stiffened more in the α-anomers and the most in the α-GalNAc-Ser structure because of the steric strains in the axial position and by two or three intramolecular hydrogen bonds. We also found that although, in our gas-phase model geometries, the peptide backbone prefers to be in a γ-turn, a structural water molecule can stabilize a polyproline II helix of a proline-rich sequence, a β-sheet, or more likely random coils. [ABSTRACT FROM AUTHOR]