AI-driven multiscale simulations illuminate mechanisms of SARS-CoV-2 spike dynamics.

We develop a generalizable AI-driven workflow that leverages heterogeneous HPC resources to explore the time-dependent dynamics of molecular systems. We use this workflow to investigate the mechanisms of infectivity of the SARS-CoV-2 spike protein, the main viral infection machinery. Our workflow enables more efficient investigation of spike dynamics in a variety of complex environments, including within a complete SARS-CoV-2 viral envelope simulation, which contains 305 million atoms and shows strong scaling on ORNL Summit using NAMD. We present several novel scientific discoveries, including the elucidation of the spike's full glycan shield, the role of spike glycans in modulating the infectivity of the virus, and the characterization of the flexible interactions between the spike and the human ACE2 receptor. We also demonstrate how AI can accelerate conformational sampling across different systems and pave the way for the future application of such methods to additional studies in SARS-CoV-2 and other molecular systems.

Specific cation effects at aqueous solution-vapor interfaces: Surfactant-like behavior of Li+ revealed by experiments and simulations.

It is now well established by numerous experimental and computational studies that the adsorption propensities of inorganic anions conform to the Hofmeister series. The adsorption propensities of inorganic cations, such as the alkali metal cations, have received relatively little attention. Here we use a combination of liquid-jet X-ray photoelectron experiments and molecular dynamics simulations to investigate the behavior of K+ and Li+ ions near the interfaces of their aqueous solutions with halide ions. Both the experiments and the simulations show that Li+ adsorbs to the aqueous solution-vapor interface, while K+ does not. Thus, we provide experimental validation of the "surfactant-like" behavior of Li+ predicted by previous simulation studies. Furthermore, we use our simulations to trace the difference in the adsorption of K+ and Li+ ions to a difference in the resilience of their hydration shells.

Simulation and theory of ions at atmospherically relevant aqueous liquid-air interfaces.

Chemistry occurring at or near the surface of aqueous droplets and thin films in the atmosphere influences air quality and climate. Molecular dynamics simulations are becoming increasingly useful for gaining atomic-scale insight into the structure and reactivity of aqueous interfaces in the atmosphere. Here we review simulation studies of atmospherically relevant aqueous liquid-air interfaces, with an emphasis on ions that play important roles in the chemistry of atmospheric aerosols. In addition to surveying results from simulation studies, we discuss challenges to the refinement and experimental validation of the methodology for simulating ion adsorption to the air-water interface and recent advances in elucidating the driving forces for adsorption. We also review the recent development of a dielectric continuum theory capable of reproducing simulation and experimental data on ion behavior at aqueous interfaces.

Interactions of gaseous HNO3 and water with individual and mixed alkyl self-assembled monolayers at room temperature.

The major removal processes for gaseous nitric acid (HNO3) in the atmosphere are dry and wet deposition onto various surfaces. The surface in the boundary layer is often covered with organic films, but the interaction of gaseous HNO3 with them is not well understood. To better understand the factors controlling the uptake of gaseous nitric acid and its dissociation in organic films, studies were carried out using single component and mixtures of C8 and C18 alkyl self-assembled monolayers (SAMs) attached to a germanium (Ge) attenuated total reflectance (ATR) crystal upon which a thin layer of SiOx had been deposited. For comparison, diffuse reflectance infrared Fourier transform spectrometry (DRIFTS) studies were also carried out using a C18 SAM attached to the native oxide layer on the surface of silicon powder. These studies show that the alkyl chain length and order/disorder of the SAMs does not significantly affect the uptake or dissociation/recombination of molecular HNO3. Thus, independent of the nature of the SAM, molecular HNO3 is observed up to 70-90% relative humidity. After dissociation, molecular HNO3 is regenerated on all SAM surfaces when water is removed. Results of molecular dynamics simulations are consistent with experiments and show that defects and pores on the surfaces control the uptake, dissociation and recombination of molecular HNO3. Organic films on surfaces in the boundary layer will certainly be more irregular and less ordered than SAMs studied here, therefore undissociated HNO3 may be present on surfaces in the boundary layer to a greater extent than previously thought. The combination of this observation with the results of recent studies showing enhanced photolysis of nitric acid on surfaces suggests that renoxification of deposited nitric acid may need to be taken into account in atmospheric models.

Thermodynamics of iodide adsorption at the instantaneous air-water interface.

We performed molecular dynamics simulations using both polarizable and non-polarizable force fields to study the adsorption of iodide to the air-water interface. A novel aspect of our analysis is that the progress of ion adsorption is measured as the distance from the instantaneous interface, which is defined by a coarse-graining scheme proposed recently by Willard and Chandler ["Instantaneous liquid interfaces," J. Phys. Chem. B 114, 1954-1958 (2010)]. Referring structural and thermodynamic quantities to the instantaneous interface unmasks molecular-scale details that are obscured by thermal fluctuations when the same quantities are referred to an average measure of the position of the interface, such as the Gibbs dividing surface. Our results suggest that an ion adsorbed at the interface resides primarily in the topmost water layer, and the interfacial location of the ion is favored by enthalpy and opposed by entropy.

Understanding hydrogen sorption in a polar metal-organic framework with constricted channels.

A high fidelity molecular model is developed for a metal-organic framework (MOF) with narrow (approximately 7.3 Å) nearly square channels. MOF potential models, both with and neglecting explicit polarization, are constructed. Atomic partial point charges for simulation are derived from both fragment-based and fully periodic electronic structure calculations. The molecular models are designed to accurately predict and retrodict material gas sorption properties while assessing the role of induction for molecular packing in highly restricted spaces. Thus, the MOF is assayed via grand canonical Monte Carlo (GCMC) for its potential in hydrogen storage. The confining channels are found to typically accommodate between two to three hydrogen molecules in close proximity to the MOF framework at or near saturation pressures. Further, the net attractive potential energy interactions are dominated by van der Waals interactions in the highly polar MOF - induction changes the structure of the sorbed hydrogen but not the MOF storage capacity. Thus, narrow channels, while providing reasonably promising isosteric heat values, are not the best choice of topology for gas sorption applications from both a molecular and gravimetric perspective.

Atomic charges derived from electrostatic potentials for molecular and periodic systems.

We present a method for fitting atomic charges to the electrostatic potential (ESP) of periodic and nonperiodic systems. This method is similar to the method of Campa et al. [ J. Chem. Theory Comput. 2009, 5, 2866]. We compare the Wolf and Ewald long-range electrostatic summation methods in calculating the ESP for periodic systems. We find that the Wolf summation is computationally more efficient than the Ewald summation by about a factor of 5 with comparable accuracy. Our analysis shows that the choice of grid mesh size influences the fitted atomic charges, especially for systems with buried (highly coordinated) atoms. We find that a maximum grid spacing of 0.2−0.3 A is required to obtain reliable atomic charges. The effect of the exclusion radius for point selection is assessed; we find that the common choice of using the van der Waals (vdW) radius as the exclusion radius for each atom may result in large deviations between the ESP generated from the ab initio calculations and that computed from the fitted charges, especially for points closest to the exclusion radii. We find that a larger value of exclusion radius than commonly used, 1.3 times the vdW radius, provides more reliable results. We find that a penalty function approach for fitting charges for buried atoms, with the target charge taken from Bader charge analysis, gives physically reasonable results.

On the mechanism of hydrogen storage in a metal-organic framework material.

Monte Carlo simulations were performed modeling hydrogen sorption in a recently synthesized metal-organic framework material (MOF) that exhibits large molecular hydrogen uptake capacity. The MOF is remarkable because at 78 K and 1.0 atm it sorbs hydrogen at a density near that of liquid hydrogen (at 20 K and 1.0 atm) when considering H2 density in the pores. Unlike most other MOFs that have been investigated for hydrogen storage, it has a highly ionic framework and many relatively small channels. The simulations demonstrate that it is both of these physical characteristics that lead to relatively strong hydrogen interactions in the MOF and ultimately large hydrogen uptake. Microscopically, hydrogen interacts with the MOF via three principle attractive potential energy contributions: Van der Waals, charge-quadrupole, and induction. Previous simulations of hydrogen storage in MOFs and other materials have not focused on the role of polarization effects, but they are demonstrated here to be the dominant contribution to hydrogen physisorption. Indeed, polarization interactions in the MOF lead to two distinct populations of dipolar hydrogen that are identified from the simulations that should be experimentally discernible using, for example, Raman spectroscopy. Since polarization interactions are significantly enhanced by the presence of a charged framework with narrow pores, MOFs are excellent hydrogen storage candidates.

Electrochemical Surface Potential Due to Classical Point Charge Models Drives Anion Adsorption to the Air-Water Interface.

We demonstrate that the driving forces for ion adsorption to the air-water interface for point charge models result from both cavitation and a term that is of the form of a negative electrochemical surface potential. We carefully characterize the role of the free energy due to the electrochemical surface potential computed from simple empirical models and its role in ionic adsorption within the context of dielectric continuum theory. Our research suggests that the electrochemical surface potential due to point charge models provides anions with a significant driving force for adsoprtion to the air-water interface. This is contrary to the results of ab initio simulations that indicate that the average electrostatic surface potential should favor the desorption of anions at the air-water interface. The results have profound implications for the studies of ionic distributions in the vicinity of hydrophobic surfaces and proteins.

Dissociation of strong acid revisited: X-ray photoelectron spectroscopy and molecular dynamics simulations of HNO3 in water.

Molecular-level insight into the dissociation of nitric acid in water is obtained from X-ray photoelectron spectroscopy and first-principles molecular dynamics (MD) simulations. Our combined studies reveal surprisingly abrupt changes in solvation configurations of undissociated nitric acid at approximately 4 M concentration. Experimentally, this is inferred from shifts of the N1s binding energy of HNO(3)(aq) as a function of concentration and is associated with variations in the local electronic structure of the nitrogen atom. It also shows up as a discontinuity in the degree of dissociation as a function of concentration, determined here from the N1s photoelectron signal intensity, which can be separately quantified for undissociated HNO(3)(aq) and dissociated NO(3)(-)(aq). Intermolecular interactions within the nitric acid solution are discussed on the basis of MD simulations, which reveal that molecular HNO(3) interacts remarkably weakly with solvating water molecules at low concentration; around 4 M there is a turnover to a more structured solvation shell, accompanied by an increase in hydrogen bonding between HNO(3) and water. We suggest that the driving force behind the more structured solvent configuration of HNO(3) is the overlap of nitric acid solvent shells that sets in around 4 M concentration.

An Accurate and Transferable Intermolecular Diatomic Hydrogen Potential for Condensed Phase Simulation.

An anisotropic many-body H2 potential energy function has been developed for use in heterogeneous systems. The intermolecular potential has been derived from first principles and expressed in a form that is readily transferred to exogenous systems, e.g. in modeling H2 sorption in solid-state materials. Explicit many-body polarization effects, known to be important in simulating hydrogen at high density, are incorporated. The analytic form of the potential energy function is suitable for methods of statistical physics, such as Monte Carlo or Molecular Dynamics simulation. The model has been validated on dense supercritical hydrogen and demonstrated to reproduce the experimental data with high accuracy.