Ensembles Are Required to Handle Aleatoric and Parametric Uncertainty in Molecular Dynamics Simulation

Classical molecular dynamics is a computer simulation technique that is in widespread use across many areas of science, from physics and chemistry to materials, biology, and medicine. The method continues to attract criticism due its oft-reported lack of reproducibility which is in part due to a failure to submit it to reliable uncertainty quantification (UQ). Here we show that the uncertainty arises from a combination of (i) the input parameters and (ii) the intrinsic stochasticity of the method controlled by the random seeds. To illustrate the situation, we make a systematic UQ analysis of a widely used molecular dynamics code (NAMD), applied to estimate binding free energy of a ligand-bound to a protein. In particular, we replace the usually fixed input parameters with random variables, systematically distributed about their mean values, and study the resulting distribution of the simulation output. We also perform a sensitivity analysis, which reveals that, out of a total of 175 parameters, just six dominate the variance in the code output. Furthermore, we show that binding energy calculations dampen the input uncertainty, in the sense that the variation around the mean output free energy is less than the variation around the mean of the assumed input distributions, if the output is ensemble-averaged over the random seeds. Without such ensemble averaging, the predicted free energy is five times more uncertain. The distribution of the predicted properties is thus strongly dependent upon the random seed. Owing to this substantial uncertainty, robust statistical measures of uncertainty in molecular dynamics simulation require the use of ensembles in all contexts

Automated Variance-Based Sensitivity Analysis of a Heterogeneous Atomistic-Continuum System

A fully automated computational tool for the study of the uncertainty in a mathematical-computational model of a heterogeneous multi-scale atomistic-continuum coupling system is implemented and tested in this project. This tool can facilitate quantitative assessments of the model’s overall uncertainty for a given specific range of variables. The computational approach here is based on the polynomial chaos expansion using projection variance, a pseudo-spectral method. It also supports regression variance, a point collocation method with nested quadrature point where the random sampling method takes a dictionary of the names of the parameters which are manually defined to vary with corresponding distributions. The tool in conjunction with an existing platform for verification, validation, and uncertainty quantification offers a scientific simulation environment and data processing workflows that enables the execution of simulation and analysis tasks on a cluster or supercomputing platform with remote submission capabilities

Crystallization and preliminary diffraction studies of native and selenomethionine CcmG (CycY, DsbE)

t Disulfide-bond (Dsb) proteins are a family of redox proteins containing a Cys-X-X-Cys motif. They are essential for disulfide-bond exchange in the bacterial periplasm and are necessary for the correct folding and function of many secreted proteins. CcmG (DsbE) is a reducing Dsb protein required for cytochrome c maturation. Crystals of Bradyrhizobium japonicum CcmG have been obtained that diffract X-rays to 1.14 Angstrom resolution. The crystals are orthorhombic, space group P2(1)2(1)2(1), with unit-cell parameters a = 35.1, b = 48.2, c = 90.2 Angstrom. Selenomethionine CcmG was expressed without using a methionine auxotroph or methionine-pathway inhibition and was purified without reducing agents

Ensembles are required to handle aleatoric and parametric uncertainty in molecular dynamics simulation

Classical molecular dynamics is a computer simulation technique that is in widespread use across many areas of science, from physics and chemistry to materials, biology, and medicine. The method continues to attract criticism due its oft-reported lac

Structural requirements for PACSIN/Syndapin operation during zebrafish embryonic notochord development.

PACSIN/Syndapin proteins are membrane-active scaffolds that participate in endocytosis. The structure of the Drosophila Syndapin N-terminal EFC domain reveals a crescent shaped antiparallel dimer with a high affinity for phosphoinositides and a unique membrane-inserting prong upon the concave surface. Combined structural, biochemical and reverse genetic approaches in zebrafish define an important role for Syndapin orthologue, Pacsin3, in the early formation of the notochord during embryonic development. In pacsin3-morphant embryos, midline convergence of notochord precursors is defective as axial mesodermal cells fail to polarize, migrate and differentiate properly. The pacsin3 morphant phenotype of a stunted body axis and contorted trunk is rescued by ectopic expression of Drosophila Syndapin, and depends critically on both the prong that protrudes from the surface of the bowed Syndapin EFC domain and the ability of the antiparallel dimer to bind tightly to phosphoinositides. Our data confirm linkage between directional migration, endocytosis and cell specification during embryonic morphogenesis and highlight a key role for Pacsin3 in this coupling in the notochord

Tutorial applications for Verification, Validation and Uncertainty Quantification using VECMA toolkit

The VECMA toolkit enables automated Verification, Validation and Uncertainty Quantification (VVUQ) for complex applications that can be deployed on emerging exascale platforms and provides support for software applications for any domain of interest. The toolkit has four main components including EasyVVUQ for VVUQ workflows, FabSim3 for automation and tool integration, MUSCLE3 for coupling multiscale models and QCG tools to execute application workflows on high performance computing (HPC). A more recent addition to the VECMAtk is EasySurrogate for various types of surrogate methods. In this paper, we present five tutorials from different application domains that apply these VECMAtk components to perform uncertainty quantification analysis, use surrogate models, couple multiscale models and execute sensitivity analysis on HPC. This paper aims to provide hands-on experience for practitioners aiming to test and contrast with their own applications

The Structure, Dynamics and Electronic Structure of Liquid Ag-Se Alloys Investigated by Ab Initio Simulation

Ab initio molecular-dynamics simulations have been used to investigate the structure, dynamics and electronic properties of the liquid alloy Ag(1-x)Se(x) at 1350 K and at the three compositions x=0.33, 0.42 and 0.65. The calculations are based on density-functional theory in the local density approximation and on the pseudopotential plane-wave method. The reliability of the simulations is confirmed by detailed comparisons with very recent neutron diffraction results for the partial structure factors and radial distribution functions (RDF) of the stoichiometric liquid Ag2Se. The simulations show a dramatic change of the Se-Se RDF with increasing Se content. This change is due to the formation of Se clusters bound by covalent bonds, the Se-Se bond length being almost the same as in pure c-Se and l-Se. The clusters are predominantly chain-like, but for higher x a large fraction of 3-fold coordinated Se atoms is also found. It is shown that the equilibrium fractions of Se present as isolated atoms and in clusters can be understood on a simple charge-balance model based on an ionic interpretation. The Ag and Se diffusion coefficients both increase with Se content, in spite of the Se clustering. An analysis of the Se-Se bond dynamics reveals surprisingly short bond lifetimes of less than 1 ps. The changes in the density of states with composition arise directly from the formation of Se-Se covalent bonds. Results for the electronic conductivity obtained using the Kubo-Greenwood approximation are in adequate agreement with experiment for l-Ag2Se, but not for the high Se contents. Possible reasons for this are discussed.Comment: 14 pages, Revtex, 14 Postscript figures embedded in the tex

Structure and relaxations in liquid and amorphous Selenium

We report a molecular dynamics simulation of selenium, described by a three-body interaction. The temperatures T_g and T_c and the structural properties are in agreement with experiment. The mean nearest neighbor coordination number is 2.1. A small pre-peak at about 1 AA^-1 can be explained in terms of void correlations. In the intermediate self-scattering function, i.e. the density fluctuation correlation, classical behavior, alpha- and beta-regimes, is found. We also observe the plateau in the beta-regime below T_g. In a second step, we investigated the heterogeneous and/or homogeneous behavior of the relaxations. At both short and long times the relaxations are homogeneous (or weakly heterogeneous). In the intermediate time scale, lowering the temperature increases the heterogeneity. We connect these different domains to the vibrational (ballistic), beta- and alpha-regimes. We have also shown that the increase in heterogeneity can be understood in terms of relaxations