The effect of methylation on the structure of uracil

Adding hydrogen atoms to molecular models via fragment superimposition

BACKGROUND

Most experimentally determined structures of biomolecules lack annotated hydrogen positions due to their low electron density. However, thorough structure analysis and simulations require knowledge about the positions of hydrogen atoms. Existing methods for their prediction are either limited to a certain range of molecules or only work effectively on small compounds.

RESULTS

We present a novel algorithm that compiles fragments of molecules with known hydrogen atom positions into a library. Using this library the method is able to predict hydrogen positions for molecules with similar moieties. We show that the method is able to accurately assign hydrogen atoms to most organic compounds including biomacromolecules, if a sufficiently large library is used.

CONCLUSIONS

We bundled the algorithm into the open-source Python package and command line program Hydride. Since usually no additional parametrization is necessary for the problem at hand, the software works out-of-box for a wide range of molecular systems usually within a few seconds of computation time. Hence, we believe that Hydride could be a valuable tool for structural biologists and biophysicists alike.

Theoretical modeling of DNA electron hole transport through polypyrimidine sequences: a QM/MM study

The phenomenon of DNA hole transport (HT) has attracted of scientists for several decades, mainly due to its potential application in molecular electronics. As electron holes mostly localize on purine bases in DNA, the majority of scientific effort has been invested into chemically modifying the structures of adenine and guanine in order to increase their HT-mediating properties. In this work we examine an alternative, never yet explored, way of affecting the HT efficiency by forcing electron holes to localize on pyrimidine bases and move between them. Using an enhanced and revised version of our previously developed QM/MM model, we perform simulations of HT through polyadenine, polycytosine, polyguanine, and polythymine stacks according to a multistep hopping mechanism. From these simulations, kinetic parameters for HT are obtained. The results indicate a particularly high efficiency of cytosine→cytosine hopping, which is about ten times higher than the G → G hopping. We also discuss possible improvement of cytosine HT by modifying the oxidoreductive properties of complementary guanine residues.