An In-Depth Look at DNA Crystals through the Prism of Molecular Dynamics Simulations

CGeNArate: a sequence-dependent coarse-grained model of DNA for accurate atomistic MD simulations of kb-long duplexes

We present CGeNArate, a new model for molecular dynamics simulations of very long segments of B-DNA in the context of biotechnological or chromatin studies. The developed method uses a coarse-grained Hamiltonian with trajectories that are back-mapped to the atomistic resolution level with extreme accuracy by means of Machine Learning Approaches. The method is sequence-dependent and reproduces very well not only local, but also global physical properties of DNA. The efficiency of the method allows us to recover with a reduced computational effort high-quality atomic-resolution ensembles of segments containing many kilobases of DNA, entering into the gene range or even the entire DNA of certain cellular organelles.

Identification of the probable structure of the sAPPα-GABABR1a complex and theoretical solutions for such cases

Amyloid precursor protein (APP) is the core of the pathogenesis of Alzheimer's disease (AD). Existing studies have shown that the soluble secreted APP (sAPPα) fragment obtained from the hydrolysis of APP by α-secretase has a synaptic function. Thereinto, a nine-residue fragment (APP9mer) of the extension domain region of sAPPα can bind directly and selectively to the N-terminal sushi1 domain (SD1) of the γ-aminobutyric acid type B receptor subunit 1a (GABA_BR1a) protein, which can influence synaptic transmission and plasticity by changing the GABA_BR1a conformation. APP9mer is a highly flexible, disordered region, and as such it is difficult to experimentally determine the optimal APPmer-SD1 binding complex. In this study we constructed two types of APP9mer-SD1 complexes through molecular docking and molecular dynamics simulation, aiming to explore the recognition function and mechanism of the specific binding of APP9mer with SD1, from which the most probable APPmer-SD1 model conformation is predicted. All the data from the analyses of RMSD, RMSF, PCA, DCCM and MM/PBSA binding energy as well as comparison with the experimental dissociation constant K_d suggest that 2NC is the most likely conformation to restore the crystal structure of the experimental APP9mer-SD1 complex. Of note, the key recognition residues of APP9mer are D24, D25, D27, W29 and W30, which mainly act on the 9-45 residue domain of SD1 (consisting of two loops and three short β-chains at the N-terminus of SD1). The mini-model with key residues identified establishes the molecular basis with deep insight into the interaction between APP and GABA_BR1a and provides a target for the development of therapeutic strategies for modulating GABA_BR1a-specific signaling in neurological and psychiatric disorders. More importantly, the study offers a theoretical solution for how to determine a biomolecular structure with a highly flexible, disordered fragment embedded within. The flexible fragment involved in a protein structure has to be deserted usually during the structural determination with experimental methods (e.g. X-ray crystallography, etc.).

Molecular basis of Arginine and Lysine DNA sequence-dependent thermo-stability modulation

We have used a variety of theoretical and experimental techniques to study the role of four basic amino acids-Arginine, Lysine, Ornithine and L-2,4-Diaminobutyric acid-on the structure, flexibility and sequence-dependent stability of DNA. We found that the presence of organic ions stabilizes the duplexes and significantly reduces the difference in stability between AT- and GC-rich duplexes with respect to the control conditions. This suggests that these amino acids, ingredients of the primordial soup during abiogenesis, could have helped to equalize the stability of AT- and GC-rich DNA oligomers, facilitating a general non-catalysed self-replication of DNA. Experiments and simulations demonstrate that organic ions have an effect that goes beyond the general electrostatic screening, involving specific interactions along the grooves of the double helix. We conclude that organic ions, largely ignored in the DNA world, should be reconsidered as crucial structural elements far from mimics of small inorganic cations.

Thermodynamics of DNA Hybridization from Atomistic Simulations

Studying DNA hybridization equilibrium at atomistic length scales, either via molecular dynamics (MD) or through commonly used advanced sampling approaches, is notoriously difficult. In this work, we describe an order-parameter-based advanced sampling technique to calculate the free energy of hybridization, and estimate the melting temperature of DNA oligomers at atomistic resolution. The free energy landscapes are reported as a function of a native-topology-based order parameter for the Drew-Dickerson dodecamer and for a range of DNA decamer sequences of different GC content. Our estimated melting temperatures match the experimental numbers within ±15 °C. As a test of the numerical reliability of the procedures employed, it was verified that the predicted free energy surfaces and melting temperatures of the d- and l-enantiomers of the Drew-Dickerson dodecamer were indistinguishable within numerical accuracy.

Framework for Conducting and Analyzing Crystal Simulations of Nucleic Acids to Aid in Modern Force Field Evaluation

Crystal simulations provide useful tools, along with solution simulations, to test nucleic acid force fields, but should be interpreted with care owing to the difficulty of establishing the environmental conditions needed to reproduce experimental crystal packing. These challenges underscore the need to construct proper protocols for carrying out crystal simulations and analyzing results to identify the origin of deviations from crystallographic data. Toward this end, we introduce a novel framework for B-factor decomposition into additive intramolecular, rotational, and translational atomic fluctuation components and partitioning of each of these components into individual asymmetric unit and lattice contributions. We apply the framework to a benchmark set of A-DNA, Z-DNA, and B-DNA double helix systems of various chain lengths. Overall, the intramolecular deviations from the crystal were quite small (≤1.0 Å), suggesting high accuracy of the force field, whereas crystal packing was not well reproduced. The present work establishes a framework to conduct and analyze crystal simulations that ultimately take on issues of crystal packing and can provide insight into nucleic acid force fields.