Free energy-based model of CTCF-mediated chromatin looping in the human genome.

In recent years, high-throughput techniques have revealed considerable structural organization of the human genome with diverse regions of the chromatin interacting with each other in the form of loops. Some of these loops are quite complex and may encompass regions comprised of many interacting chain segments around a central locus. Popular techniques for extracting this information are chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and high-throughput chromosome conformation capture (Hi-C). Here, we introduce a physics-based method to predict the three-dimensional structure of chromatin from population-averaged ChIA-PET data. The approach uses experimentally-validated data from human B-lymphoblastoid cells to generate 2D meta-structures of chromatin using a dynamic programming algorithm that explores the chromatin free energy landscape. By generating both optimal and suboptimal meta-structures we can calculate both the free energy and additionally the relative thermodynamic probability. A 3D structure prediction program with applied restraints then can be used to generate the tertiary structures. The main advantage of this approach for population-averaged experimental data is that it provides a way to distinguish between the principal and the spurious contacts. This study also finds that euchromatin appear to have rather precisely regulated 2D meta-structures compared to heterochromatin. The program source-code is available at https://github.com/plewczynski/looper.

RNA-Puzzles Round III: 3D RNA structure prediction of five riboswitches and one ribozyme.

RNA-Puzzles is a collective experiment in blind 3D RNA structure prediction. We report here a third round of RNA-Puzzles. Five puzzles, 4, 8, 12, 13, 14, all structures of riboswitch aptamers and puzzle 7, a ribozyme structure, are included in this round of the experiment. The riboswitch structures include biological binding sites for small molecules (S-adenosyl methionine, cyclic diadenosine monophosphate, 5-amino 4-imidazole carboxamide riboside 5'-triphosphate, glutamine) and proteins (YbxF), and one set describes large conformational changes between ligand-free and ligand-bound states. The Varkud satellite ribozyme is the most recently solved structure of a known large ribozyme. All puzzles have established biological functions and require structural understanding to appreciate their molecular mechanisms. Through the use of fast-track experimental data, including multidimensional chemical mapping, and accurate prediction of RNA secondary structure, a large portion of the contacts in 3D have been predicted correctly leading to similar topologies for the top ranking predictions. Template-based and homology-derived predictions could predict structures to particularly high accuracies. However, achieving biological insights from de novo prediction of RNA 3D structures still depends on the size and complexity of the RNA. Blind computational predictions of RNA structures already appear to provide useful structural information in many cases. Similar to the previous RNA-Puzzles Round II experiment, the prediction of non-Watson-Crick interactions and the observed high atomic clash scores reveal a notable need for an algorithm of improvement. All prediction models and assessment results are available at http://ahsoka.u-strasbg.fr/rnapuzzles/.

Electron Transport in a Dioxygenase-Ferredoxin Complex: Long Range Charge Coupling between the Rieske and Non-Heme Iron Center.

Dioxygenase (dOx) utilizes stereospecific oxidation on aromatic molecules; consequently, dOx has potential applications in bioremediation and stereospecific oxidation synthesis. The reactive components of dOx comprise a Rieske structure Cys2[2Fe-2S]His2 and a non-heme reactive oxygen center (ROC). Between the Rieske structure and the ROC, a universally conserved Asp residue appears to bridge the two structures forming a Rieske-Asp-ROC triad, where the Asp is known to be essential for electron transfer processes. The Rieske and ROC share hydrogen bonds with Asp through their His ligands; suggesting an ideal network for electron transfer via the carboxyl side chain of Asp. Associated with the dOx is an itinerant charge carrying protein Ferredoxin (Fdx). Depending on the specific cognate, Fdx may also possess either the Rieske structure or a related structure known as 4-Cys-[2Fe-2S] (4-Cys). In this study, we extensively explore, at different levels of theory, the behavior of the individual components (Rieske and ROC) and their interaction together via the Asp using a variety of density function methods, basis sets, and a method known as Generalized Ionic Fragment Approach (GIFA) that permits setting up spin configurations manually. We also report results on the 4-Cys structure for comparison. The individual optimized structures are compared with observed spectroscopic data from the Rieske, 4-Cys and ROC structures (where information is available). The separate pieces are then combined together into a large Rieske-Asp-ROC (donor/bridge/acceptor) complex to estimate the overall coupling between individual components, based on changes to the partial charges. The results suggest that the partial charges are significantly altered when Asp bridges the Rieske and the ROC; hence, long range coupling through hydrogen bonding effects via the intercalated Asp bridge can drastically affect the partial charge distributions compared to the individual isolated structures. The results are consistent with a proton coupled electron transfer mechanism.