Molecular dynamics simulations have been used to study the differences between two DNA and RNA 14-mer quadruplexes of analogous sequences. Their structures present a completely different fold: DNA forms a bimolecular quadruplex containing antiparallel strands and diagonal loops; RNA forms an intrastrand parallel quadruplex containing a G-tetrad and an hexad, which dimerizes by hexad stacking. We used a multiscale computational approach combining classical Molecular dynamics simulations and density functional theory calculations to elucidate the difference in stability of the 2-folds and their ability in coordinating cations. The presence of 2'-OH groups in the RNA promotes the formation of a large number of intramolecular hydrogen bonds that account for the difference in fold and stability of the two 14-mers. We observe that the adenines in the RNA quadruplex play a key role in conserving the geometry ...
A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on deceptively small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents extensive Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular "puzzles" should serve as ...