Using Pilot Systems to Execute Many Task Workloads on Supercomputers

High performance computing systems have historically been designed to support applications comprised of mostly monolithic, single-job workloads. Pilot systems decouple workload specification, resource selection, and task execution via job placeholders and late-binding. Pilot systems help to satisfy the resource requirements of workloads comprised of multiple tasks. RADICAL-Pilot (RP) is a modular and extensible Python-based pilot system. In this paper we describe RP's design, architecture and implementation, and characterize its performance. RP is capable of spawning more than 100 tasks/second and supports the steady-state execution of up to 16K concurrent tasks. RP can be used stand-alone, as well as integrated with other application-level tools as a runtime system

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein–DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50–100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA

The effect of S-substitution at the O6-guanine site on the structure and dynamics of a DNA oligomer containing a G:T mismatch

The effect of S-substitution on the O6 guanine site of a 13-mer DNA duplex containing a G:T mismatch is studied using molecular dynamics. The structure, dynamic evolution and hydration of the S-substituted duplex are compared with those of a normal duplex, a duplex with Ssubstitution on guanine, but no mismatch and a duplex with just a G:T mismatch. The S-substituted mismatch leads to cell death rather than repair. One suggestion is that the G:T mismatch recognition protein recognises the S-substituted mismatch (GS:T) as G:T. This leads to a cycle of futile repair ending in DNA breakage and cell death. We find that some structural features of the helix are similar for the duplex with the G:T mismatch and that with the S-substituted mismatch, but differ from the normal duplex, notably the helical twist. These differences arise from the change in the hydrogen-bonding pattern of the base pair. However a marked feature of the S-substituted G:T mismatch duplex is a very large opening. This showed considerable variability. It is suggested that this enlarged opening would lend support to an alternative model of cell death in which the mismatch protein attaches to thioguanine and activates downstream damage-response pathways. Attack on the sulphur by reactive oxygen species, also leading to cell death, would also be aided by the large, variable opening

The Role of Methylation in the Intrinsic Dynamics of B- and Z-DNA

Methylation of cytosine at the 5-carbon position (5mC) is observed in both prokaryotes and eukaryotes. In humans, DNA methylation at CpG sites plays an important role in gene regulation and has been implicated in development, gene silencing, and cancer. In addition, the CpG dinucleotide is a known hot spot for pathologic mutations genome-wide. CpG tracts may adopt left-handed Z-DNA conformations, which have also been implicated in gene regulation and genomic instability. Methylation facilitates this B-Z transition but the underlying mechanism remains unclear. Herein, four structural models of the dinucleotide d(GC)5 repeat sequence in B-, methylated B-, Z-, and methylated Z-DNA forms were constructed and an aggregate 100 nanoseconds of molecular dynamics simulations in explicit solvent under physiological conditions was performed for each model. Both unmethylated and methylated B-DNA were found to be more flexible than Z-DNA. However, methylation significantly destabilized the BII, relative to the BI, state through the Gp5mC steps. In addition, methylation decreased the free energy difference between B- and Z-DNA. Comparisons of α/γ backbone torsional angles showed that torsional states changed marginally upon methylation for B-DNA, and Z-DNA. Methylation-induced conformational changes and lower energy differences may contribute to the transition to Z-DNA by methylated, over unmethylated, B-DNA and may be a contributing factor to biological function

Infrared Multiple Photon Dissociation Action Spectroscopy and Theoretical Studies of Diethyl Phosphate Complexes: Effects of Protonation and Sodium Cationization on Structure

The gas-phase structures of deprotonated, protonated, and sodium-cationized complexes of diethyl phosphate (DEP) including [DEP − H]−, [DEP + H]+, [DEP + Na]+, and [DEP − H + 2Na]+ are examined via infrared multiple photon dissociation (IRMPD) action spectroscopy using tunable IR radiation generated by a free electron laser, a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS) with an electrospray ionization (ESI) source, and theoretical electronic structure calculations. Measured IRMPD spectra are compared to linear IR spectra calculated at the B3LYP/6-31G(d,p) level of theory to identify the structures accessed in the experimental studies. For comparison, theoretical studies of neutral complexes are also performed. These experiments and calculations suggest that specific geometric changes occur upon the binding of protons and/or sodium cations, including changes correlating to nucleic acid backbone geometry, specifically P–O bond lengths and ∠OPO bond angles. Information from these observations may be used to gain insight into the structures of more complex systems, such as nucleotides and solvated nucleic acids

Dimensionality of Carbon Nanomaterials Determines the Binding and Dynamics of Amyloidogenic Peptides: Multiscale Theoretical Simulations

Experimental studies have demonstrated that nanoparticles can affect the rate of protein self-assembly, possibly interfering with the development of protein misfolding diseases such as Alzheimer's, Parkinson's and prion disease caused by aggregation and fibril formation of amyloid-prone proteins. We employ classical molecular dynamics simulations and large-scale density functional theory calculations to investigate the effects of nanomaterials on the structure, dynamics and binding of an amyloidogenic peptide apoC-II(60-70). We show that the binding affinity of this peptide to carbonaceous nanomaterials such as C60, nanotubes and graphene decreases with increasing nanoparticle curvature. Strong binding is facilitated by the large contact area available for π-stacking between the aromatic residues of the peptide and the extended surfaces of graphene and the nanotube. The highly curved fullerene surface exhibits reduced efficiency for π-stacking but promotes increased peptide dynamics. We postulate that the increase in conformational dynamics of the amyloid peptide can be unfavorable for the formation of fibril competent structures. In contrast, extended fibril forming peptide conformations are promoted by the nanotube and graphene surfaces which can provide a template for fibril-growth

The Formation and Stabilization of a Novel G-Quadruplex in the 5′-Flanking Region of the Relaxin Gene

It has been reported that binding of STAT3 protein to the 5′-flanking region of the relaxin gene may result in downregulation of the relaxin expression. There is a Guanine(G)-rich segment located in about 3.8 Kb upstream of the relaxin gene and very close to the STAT3's binding site. In our study, NMR spectroscopy revealed the formation of G-quadruplex by this G-rich strand, and the result was confirmed by ESI mass spectrometry and CD spectroscopy. The theoretical structure of RLX G-quadruplex was constructed and refined by molecular modeling. When this relaxin G-quadruplex was stabilized by berberine(ΔTm = 10°C), a natural alkaloid from a Chinese herb, the gene expression could be up-regulated in a dose-dependent manner which was proved by luciferase assay. This result is different from the general G-quadruplex function that inhibiting the telomere replication or down-regulating many oncogenes expression. Therefore, our study reported a novel G-quadruplex in the relaxin gene and complemented the regulation mechanism about gene expression by G-quadruplexes

Functional Analysis and Molecular Dynamics Simulation of LOX-1 K167N Polymorphism Reveal Alteration of Receptor Activity

The human lectin-like oxidized low density lipoprotein receptor 1 LOX-1, encoded by the ORL1 gene, is the major scavenger receptor for oxidized low density lipoprotein in endothelial cells. Here we report on the functional effects of a coding SNP, c.501G>C, which produces a single amino acid change (K>N at codon 167). Our study was aimed at elucidating whether the c.501G>C polymorphism changes the binding affinity of LOX-1 receptor altering its function. The presence of p.K167N mutation reduces ox-LDL binding and uptake. Ox-LDL activated extracellular signal-regulated kinases 1 and 2 (ERK 1/2) is inhibited. Furthermore, ox-LDL induced biosynthesis of LOX-1 receptors is dependent on the p.K167N variation. In human macrophages, derived from c.501G>C heterozygous individuals, the ox-LDL induced LOX-1 46 kDa band is markedly lower than in induced macrophages derived from c.501G>C controls. Investigation of p.K167N mutation through molecular dynamics simulation and electrostatic analysis suggests that the ox-LDL binding may be attributed to the coupling between the electrostatic potential distribution and the asymmetric flexibility of the basic spine residues. The N/N-LOX-1 mutant has either interrupted electrostatic potential and asymmetric fluctuations of the basic spine arginines

The non-octarepeat copper binding site of the prion protein is a key regulator of prion conversion

The conversion of the prion protein (PrP(C)) into prions plays a key role in transmissible spongiform encephalopathies. Despite the importance for pathogenesis, the mechanism of prion formation has escaped detailed characterization due to the insoluble nature of prions. PrP(C) interacts with copper through octarepeat and non-octarepeat binding sites. Copper coordination to the non-octarepeat region has garnered interest due to the possibility that this interaction may impact prion conversion. We used X-ray absorption spectroscopy to study copper coordination at pH 5.5 and 7.0 in human PrP(C) constructs, either wild-type (WT) or carrying pathological mutations. We show that mutations and pH cause modifications of copper coordination in the non-octarepeat region. In the WT at pH 5.5, copper is anchored to His96 and His111, while at pH 7 it is coordinated by His111. Pathological point mutations alter the copper coordination at acidic conditions where the metal is anchored to His111. By using in vitro approaches, cell-based and computational techniques, we propose a model whereby PrP(C) coordinating copper with one His in the non-octarepeat region converts to prions at acidic condition. Thus, the non-octarepeat region may act as the long-sought-after prion switch, critical for disease onset and propagation