Jose M. Borreguero

Mantid—Data analysis and visualization package for neutron scattering and μ SR experiments

The Mantid framework is a software solution developed for the analysis and visualization of neutron scattering and muon spin measurements. The framework is jointly developed by software engineers and scientists at the ISIS Neutron and Muon Facility and the Oak Ridge National Laboratory. The objectives, functionality and novel design aspects of Mantid are described.

Redox-Promoting Protein Motions in Rubredoxin

Proteins are dynamic objects, constantly undergoing conformational fluctuations, yet the linkage between internal protein motion and function is widely debated. This study reports on the characterization of temperature-activated collective and individual atomic motions of oxidized rubredoxin, a small 53 residue protein from thermophilic Pyrococcus furiosus (RdPf). Computational modeling allows detailed investigations of protein motions as a function of temperature, and neutron scattering experiments are used to compare to computational results. Just above the dynamical transition temperature which marks the onset of significant anharmonic motions of the protein, the computational simulations show both a significant reorientation of the average electrostatic force experienced by the coordinated Fe(3+) ion and a dramatic rise in its strength. At higher temperatures, additional anharmonic modes become activated and dominate the electrostatic fluctuations experienced by the ion. At 360 K, close to the optimal growth temperature of P. furiosus, simulations show that three anharmonic modes including motions of two conserved residues located at the protein active site (Ile7 and Ile40) give rise to the majority of the electrostatic fluctuations experienced by the Fe(3+) ion. The motions of these residues undergo displacements which may facilitate solvent access to the ion.

BEAM: A Computational Workflow System for Managing and Modeling Material Characterization Data in HPC Environments

Abstract Improvements in scientific instrumentation allow imaging at mesoscopic to atomic length scales, many spectroscopic modes, and now—with the rise of multimodal acquisition systems and the associated processing capability—the era of multidimensional, informationally dense data sets has arrived. Technical issues in these combinatorial scientific fields are exacerbated by computational challenges best summarized as a necessity for drastic improvement in the capability to transfer, store, and analyze large volumes of data. The Bellerophon Environment for Analysis of Materials (BEAM) platform provides material scientists the capability to directly leverage the integrated computational and analytical power of High Performance Computing (HPC) to perform scalable data analysis and simulation via an intuitive, cross-platform client user interface. This framework delivers authenticated, “push-button” execution of complex user workflows that deploy data analysis algorithms and computational simulations utilizing the converged compute-and-data infrastructure at Oak Ridge National Laboratorys (ORNL) Compute and Data Environment for Science (CADES) and HPC environments like Titan at the Oak Ridge Leadership Computing Facility (OLCF). In this work we address the underlying HPC needs for characterization in the material science community, elaborate how BEAMs design and infrastructure tackle those needs, and present a small sub-set of user cases where scientists utilized BEAM across a broad range of analytical techniques and analysis modes.

Protein-Style Dynamical Transition in a Non-Biological Polymer and a Non-Aqueous Solvent

Temperature-dependent onset of apparent anharmonicity in the microscopic dynamics of hydrated proteins and other biomolecules has been known as protein dynamical transition for the last quarter of a century. Using neutron scattering and molecular dynamics simulation, techniques most often associated with protein dynamical transition studies, we have investigated the microscopic dynamics of one of the most common polymers, polystyrene, which was exposed to toluene vapor, mimicking the process of protein hydration from water vapor. Polystyrene with adsorbed toluene is an example of a solvent-solute system, which, unlike biopolymers, is anhydrous and lacks hydrogen bonding. Nevertheless, it exhibits the essential traits of the dynamical transition in biomolecules, such as a specific dependence of the microscopic dynamics of both solvent and host on the temperature and the amount of solvent adsorbed. We conclude that the protein dynamical transition is a manifestation of a universal solvent-solute dynamical relationship, which is not specific to either biomolecules as solute, or aqueous media as solvent, or even a particular type of interactions between solvent and solute.

Self-assembly and structural relaxation in a model ionomer melt

Molecular dynamics simulations are used to understand the self-assembly and structural relaxation in ionomer melts containing less than 10% degree of ionization on the backbone. The self-assembly of charged sites and counterions shows structural ordering and agglomeration with a range of structures that can be achieved by changing the dielectric constant of the medium. The intermediate scattering function shows a decoupling of charge and counterion relaxation at longer length scales for only high dielectric constant and at shorter length scales for all dielectric constants. Overall, the slow structural decay of counterions in the strongly correlated ionomer system closely resembles transport properties of semi-flexible polymers.

An automated analysis workflow for optimization of force-field parameters using neutron scattering data

Large-scale simulations and data analysis are often required to explain neutron scattering experiments to establish a connection between the fundamental physics at the nanoscale and data probed by neutrons. However, to perform simulations at experimental conditions it is critical to use correct force-field (FF) parameters which are unfortunately not available for most complex experimental systems. In this work, we have developed a workflow optimization technique to provide optimized FF parameters by comparing molecular dynamics (MD) to neutron scattering data. We describe the workflow in detail by using an example system consisting of tRNA and hydrophilic nanodiamonds in a deuterated water (D2O) environment. Quasi-elastic neutron scattering (QENS) data show a faster motion of the tRNA in the presence of nanodiamond than without the ND. To compare the QENS and MD results quantitatively, a proper choice of FF parameters is necessary. We use an efficient workflow to optimize the FF parameters between the hydrophilic nanodiamond and water by comparing to the QENS data. Our results show that we can obtain accurate FF parameters by using this technique. The workflow can be generalized to other types of neutron data for FF optimization, such as vibrational spectroscopy and spin echo.

Integrating advanced materials simulation techniques into an automated data analysis workflow at the Spallation Neutron Source

This presentation will review developments on the integration of advanced modeling and simulation techniques into the analysis step of experimental data obtained at the Spallation Neutron Source. A workflow framework for the purpose of refining molecular mechanics force-fields against quasi-elastic neutron scattering data is presented. The workflow combines software components to submit model simulations to remote high performance computers, a message broker interface for communications between the optimizer engine and the simulation production step, and tools to convolve the simulated data with the experimental resolution. A test application shows the correction to a popular fixed-charge water model in order to account polarization effects due to the presence of solvated ions. Future enhancements to the refinement workflow are discussed. This work is funded through the DOE Center for Accelerating Materials Modeling.

Precise determination of water exchanges on a mineral surface

Solvent exchanges on solid surfaces and dissolved ions are a fundamental property important for understanding chemical reactions, but the rates of fast exchanges are poorly constrained. We probed the diffusional motions of water adsorbed onto nanoparticles of the mineral barite (BaSO4) using quasi-elastic neutron scattering (QENS) and classical molecular dynamics (MD) to reveal the complex dynamics of water exchange along mineral surfaces. QENS data as a function of temperature and momentum transfer (Q) were fit using scattering functions derived from MD trajectories. The simulations reproduce the dynamics measured in the experiments at ambient temperatures, but as temperature is lowered the simulations overestimate slower motions. Decomposition of the MD-computed QENS intensity into contributions from adsorbed and unbound water shows that the majority of the signal arises from adsorbed species, although the dynamics of unbound water cannot be dismissed. The mean residence times of water on each of the four surface sites present on the barite {001} were calculated using MD: at room temperature the low barium site is 194 ps, whereas the high barium site contains two distributions of motions at 84 and 2.5 ps. These contrast to 13 ps residence time on both sulfate sites, with an additional surface diffusion exchange of 66 ps. Surface exchanges are similar to those of the aqueous ions calculated using the same force field: Baaq2+ is 208 ps and SO4aq2- is 5.8 ps. This work demonstrates how MD can be a reliable method to deconvolute solvent exchange reactions when quantitatively validated by QENS measurements.

Modulating Enzyme Activity by Altering Protein Dynamics with Solvent

Optimal enzyme activity depends on a number of factors, including structure and dynamics. The role of enzyme structure is well recognized; however, the linkage between protein dynamics and enzyme activity has given rise to a contentious debate. We have developed an approach that uses an aqueous mixture of organic solvent to control the functionally relevant enzyme dynamics (without changing the structure), which in turn modulates the enzyme activity. Using this approach, we predicted that the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) from Escherichia coli in aqueous mixtures of isopropanol (IPA) with water will decrease by ∼3 fold at 20% (v/v) IPA concentration. Stopped-flow kinetic measurements find that the pH-independent khydride rate decreases by 2.2 fold. X-ray crystallographic enzyme structures show no noticeable differences, while computational studies indicate that the transition state and electrostatic effects were identical for water and mixed solvent conditions; quasi-elastic neutron scattering studies show that the dynamical enzyme motions are suppressed. Our approach provides a unique avenue to modulating enzyme activity through changes in enzyme dynamics. Further it provides vital insights that show the altered motions of DHFR cause significant changes in the enzymes ability to access its functionally relevant conformational substates, explaining the decreased khydride rate. This approach has important implications for obtaining fundamental insights into the role of rate-limiting dynamics in catalysis and as well as for enzyme engineering.

Dynamics Of Charged Species in Ionic-Neutral Block Copolymer and Surfactant Complexes

Structure-property relationships of ionic block copolymer (BCP) surfactant complexes are critical toward the progress of favorable engineering design of efficient charge-transport materials. In this article, molecular dynamics simulations are used to understand the dynamics of charged-neutral BCP and surfactant complexes. The dynamics are examined for two different systems: charged-neutral double-hydrophilic and hydrophobic-hydrophilic block copolymers with oppositely charged surfactant moieties. The dynamics of the surfactant head, tails, and charges are studied for five different BCP volume fractions. We observe that the dynamics of the different species solely depend on the balance between electrostatic and entropic interactions between the charged species and the neutral monomers. The favorable hydrophobic-hydrophobic interactions and the unfavorable hydrophobic-hydrophilic interactions determine the mobilities of the monomers. The dynamical properties of the charge species influence complex formation. Structural relaxations exhibit length-scale dependent behavior, with slower relaxation at the radius of gyration length-scale and faster relaxation at the segmental length-scale, consistent with previous results. The dynamical analysis correlates ion-exchange kinetics to the self-assembly behavior of the complexes.