Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Pushing and Pulling: A Dual pH Trigger Controlled by Varying the Alkyl Tail Length in Heme Coordinating Peptide Amphiphiles

Some organisms in nature that undergo anaerobic respiration utilize 1D nanoscale arrays of densely packed cytochromes containing the molecule heme. The assemblies can be mimicked with 1D nanoscale fibrils composed of peptide amphiphiles designed to coordinate heme in dense arrays. To create such materials and assemblies it is critical to understand the assembly process and what controls the various aspects of hierarchical assembly. MD simulations suggest that shorter alkyl chains on the peptide lead to more dynamic structures than the peptides with longer chains that yield kinetically trapped states. The hydration parameters manifest themselves experimentally through the observation of a dual pH trigger which controls the peptide assembly rate, the heme binding affinity, and heme organization kinetics. Great strides in understanding the relative complexity of the self-assembly process in relation to incorporating a functional moiety like heme opens up many possibilities in developing abiotic assemblies for bioelectronic devices and assemblies

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Concentration-Dependent Solvation Structure and Dynamics of Aqueous Sulfuric Acid Using Multinuclear NMR and DFT

Sulfuric acid is a ubiquitous compound for industrial processes and aqueous sulfate solutions also play a critical role in electrolytes for many prominent battery chemistries. While the thermodynamic literature for it is quite well-developed, comprehensive studies of the solvation structure, and particularly molecular-scale dynamical and transport properties, are less available. This study applies a multinuclear NMR approach to the elucidation of the solvation structure and dynamics over a wide temperature (-10 – 50?C) and concentration (0 – 18 M) range, combining 17O shift, linewidth, and T1 relaxation measurements, 33S shift and linewidth measurements, and 1H PFG-NMR measurements of the proton self-diffusivity. In conjunction, these results indicate a crossover between two regimes of solvation structure and dynamics, occurring above the concentration associated with the deep eutectic point (~ 4.5 M), with the high-concentration regime dominated by strong water-sulfate correlation. This description was borne out in detail by the activation energy trends with increasing concentration derived from the relaxation of both the H2O/H3O+ and H2SO4/HSO4-/SO42- 17O resonances and the 1H self-diffusivity. However, the 17O chemical shift difference between the H2O/H3O+ and H2SO4/HSO4-/SO42- resonances across the entire temperature range is nevertheless strikingly linear. A computational approach coupling MD simulations and DFT NMR shift calculations to reproduce this trend is presented, which will be the subject of further development. This combination of multinuclear, dynamical NMR and computational methods, and the results furnished by this study, will provide a platform for future studies on battery electrolytes where aqueous sulfate chemistry plays a central role in the solution structure.

Bazak, Jonathan D.↗
DOE OSTI2022

Modeling Interactions within and between Peptide Amphiphile Supramolecular Filaments

Many peptides are able to self-assemble into one-dimensional (1D) nanostructures, such as cylindrical fibers or ribbons of variable widths, but the relationship between the morphology of 1D objects and their molecular structure is not well understood. Here, we use coarse-grained molecular dynamics (CG-MD) simulations to study the nanostructures formed by self-assembly of different peptide amphiphiles (PAs). The results show that ribbons are hierarchical superstructures formed by laterally assembled cylindrical fibers. Simulations starting from bilayer structures demonstrate the formation of filaments, whereas other simulations starting from filaments indicate varying degrees of interaction among them depending on chemical structure. These interactions are verified by observations using atomic force microscopy of the various systems. The interfilament interactions are predicted to be strongest in supramolecular assemblies that display hydrophilic groups on their surfaces, while those with hydrophobic ones are predicted to interact more weakly as confirmed by viscosity measurements. The simulations also suggest that peptide amphiphiles with hydrophobic termini bend to reduce their interfacial energy with water, which may explain why these systems do not collapse into superstructures of bundled filaments. The simulations suggest that future experiments will need to address mechanistic questions about the self-assembly of these systems into hierarchical structures, namely, the preformation of interactive filaments vs equilibration of large assemblies into superstructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Hydrophilic and Apolar Hydration in Densely Grafted Cationic Brushes and Counterions with Large Mobilities

We employ all-atom molecular dynamics (MD) simulation framework to unravel water microstructure and ion properties for cationic [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] (PMETAC) brushes with chloride ions as counterions. First, we identify locally separate water domains (or first hydration shells) each around the {N(CH 3 ) 3 } + and the C=O functional groups of the PMETAC chain and one around the Cl - ion. These first hydration shells around the respective moieties overlap and the extent of the overlap depends on the nature of the species triggering it. Second, despite the overlap, the water molecules in these domains demonstrate disparate properties dictated by the properties of the atoms/groups around which they are located. For example, the presence of the methyl groups make the {N(CH 3 ) 3 } + group trigger apolar hydration as evidenced by the corresponding orientation of the dipole of the water molecules around the {N(CH 3 ) 3 } + moiety. These water molecules around the N(CH 3 ) 3 } + group also have enhanced tetrahedrality as compared to the water molecules constituting the hydration layer around the C=O group and the Cl - counterion. Our simulations also identify that there is an intervening water layer between the Cl - ion and {N(CH 3 ) 3 } + group: this layer prevents the Cl - ion from coming very close to the {N(CH 3 ) 3 } + group. As a consequence, there is a significantly large mobility of the Cl - ions inside the PMETAC brush layer. Furthermore, the C=O group of the PE chain, due to the partial negative charge on the oxygen atom and the specific structure of the PMETAC brush system, demonstrates strongly hydrophilic behavior and enforces a specific dipole response of water molecules analogous to that experienced by water around anionic species of high charge density. As a result, our findings confirm that PMETAC brushes undergo hydrophilic hydration at one site, apolar hydration at another site, and ensures large mobility of the supported Cl- counterions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Thermodynamic Analysis of Silk Fibroin–Graphite Hybrid Materials and Their Morphology

Silk fibroin (SF) is a β-sheet-rich protein that is responsible for the remarkable tensile strength of silk. In addition to its mechanical properties, SF is biocompatible and biodegradable, making it an attractive candidate for use in biotic/abiotic hybrid materials. A pairing of particular interest is the use of SF with graphene-based nanomaterials (GBNs). The properties of this interface drive the formation of well-ordered nanostructures and can improve the electronic properties of the resulting hybrid. It was previously demonstrated that SF can form lamellar nanostructures in the presence of graphite; however, the equilibrium morphology and associated driving interactions are not fully understood. Here, in this study, we characterize these interactions between SF and SF lamellar with graphite using molecular dynamics (MD) simulations and umbrella sampling (US). We find that SF lamellar nanostructures have strong orientational and spatial preferences on graphite that are driven by the hydrophobic effect, destabilizing solvent–protein interactions and stabilizing protein–protein and protein–graphite interactions. Finally, we show how careful consideration of these underlying interactions can be applied to rationally modify the nanostructure morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Probing the Solvation Shells of Lithium Ions in Glyme-Based Electrolytes

Glymes have been extensively studied as solvents for Li-battery electrolytes, most recently in equimolar mixtures with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI), due to their ability to form stable solvates. However, directly quantifying free and coordinated glyme molecules in the liquid state has been challenging due to several experimental limitations. Here, in this work, new vibrational probes are demonstrated for studying the solvation structures of diglyme and triglyme in LiTFSI electrolytes. These IR probes make use of an amine group to report the solvation state of glymes at salt-to-solvent molar ratios ranging from 1:5 to 1:10. Characterization of the thermodynamic properties of the solvent exchange occurring in the first solvation shell of lithium ions (Li + ) showed an equilibrium constant for these probes close to unity at room temperature. This result demonstrates that the probes exhibit a similar solvation behavior to their glyme analogue. Concentration dependence studies also revealed a lack of significant amounts of contact ion pairs at the studied concentrations. Moreover, the first solvation shell of Li + appears to be formed by two partially chelating glyme molecules, establishing that even triglyme with multiple chelation sites does not fully coordinate the cation. Complementary molecular dynamics (MD) simulations agree with the experimental results and suggest that at these concentrations, TFSI – predominantly forms solvent-separated ion pairs. However, the simulations do not properly capture the partial solvation structure of the glyme molecules in the solvation shell of Li + as derived from the experiments.

electrolytes↗
DOE OSTI2026

Ionic Liquid-Enhanced Interfaces to Boost Reactive C O2 Capture

The addition of ionic liquids (ILs) to a mixture containing a molecular solvent and other ionic species can induce the heterogeneous redistribution of cations and anions at the gas–liquid interface. This nonuniform redistribution of cations and anions driven by the differences in the solvophilicity of ions can improve the thermophysical and interfacial properties of such mixtures, creating a local chemical environment that is conducive to some reactions. In this work, ILs are added to a mixture of potassium hydroxide (KOH) and ethylene glycol (EG), used as a reactive absorbent and electrolyte in the migration-assisted moisture-gradient (MAMG) process for CO 2 capture. Molecular dynamics (MD) simulations are employed to probe into the effects of complex ion–ion and ion–solvent interactions and to examine the chemical composition at the gas–liquid interface. A total of 12 systems are investigated using molecular simulations to identify trends in the performance of IL additives based on the choice of cation, anion, and IL concentration. The cation effects are studied using IL additives based on 1-ethyl-3-methylimidazolium ([EMIM] + ) and 1-butyl-3-methylimidazolium ([BMIM] + ), while the impact of anions is examined using additives based on dicyanamide [DCA] − , triflate [TfO] − , bistriflimide [NTf 2 ] − , and hexafluorophosphate [PF 6 ] − anions, respectively. The influence of the IL concentration is also evaluated at molar concentrations between 1% and 4%. The simulation results indicate that the use of IL additives can affect the physical CO 2 solubility, surface tension, and the localization of CO 2 around the [OH] − ions at the gas–liquid interface. It is also evident that the choice of cations, anions, and IL concentration determines the extent to which the IL additives impact the local physicochemical properties. Physical dissolution, diffusive transport, and interaction with [OH] − are critical intermediate steps toward reactive CO 2 capture using a liquid absorbent. Hence, the improvement in one or more of these properties, aided by IL additives, is expected to improve the overall CO 2 capture performance. Experiments reaffirmed the impact of IL additives on CO 2 capture performance and the sensitivity to the choice of the cation, anion, and concentration of the IL additive.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Phosphoserine Charge State Drives Ion Condensation and Spatial Polyamine Presentation in Multirepeat Silaffin

Diatom silaffins direct silica biomineralization through heavily post-translationally modified repeat domains, yet how these modifications reshape the multirepeat conformational ensemble remains unknown. We report all-atom MD simulations of a 195-residue construct spanning repeats R1–R7 of Sil1p from Cylindrotheca fusiformis, carrying the full complement of native PTMs: phosphoserine (pSer), long-chain polyamines (LCPA), dimethyllysine (MLY), and trimethylhydroxylysine phosphate (TPL). We simulated three variants (Native, P1/singly deprotonated phosphate, and P2/doubly deprotonated phosphate) at two NaCl concentrations in triplicate for 500 ns each. All systems disorder from the AlphaFold 3 starting structure. Phosphate charge state, not ionic strength, is the dominant control of ensemble compaction and ion organization. Doubly deprotonated phosphate organizes an extensive Na + condensation shell (∼100 ions, 20% of box Na + ) and a heterogeneous bridging network that integrates both pSer and TPL phosphate groups. The resulting ensemble is compact with LCPA side chains exhibiting above-median solvent accessibility in 84% of simulation frames in P2 at 300 mM. This is higher than any other condition we simulated and indicates that polyamine groups are preferentially surface-presented in the most compact, ion-organized state. A charge-neutralization control confirms that this compact state is a structured intermediate maintained by the bridging network, not a simple collapsed globule. Here, this repeat-scale spatial organization is not captured by single-repeat peptide studies. Understanding the dynamics, mechanism, and spatial organization of PTM-rich silaffin at the repeat scale is a step closer to hierarchical biomimetic materials beyond simple silica morphologies.

Amines↗
DOE OSTI2021

The Role of Phase Separation on Rayleigh-Plateau Type Instabilities in Alloys

Classical molecular dynamics (MD) simulations are used to investigate the role of phase separation (PS) on the Rayleigh-Plateau (RP) instability. Ni–Ag bulk structures are created at temperatures (2000 K and 1400 K) that generate different PS length scales, λ PS , relative to the RP instability length scale, λRP. Rectanguloids are then cut from the bulk structures and patterned with a perturbation of certain amplitude and wavelength, λ RP . It is found that when λ PS << λ RP (2000 K), the patterned rectanguloids break up into nanoparticles in a manner consistent with classical RP theory, whereas when λ PS << λ RP (1400 K), soluto-capillarity affects the RP instability significantly. Specifically, since Ag has a lower surface energy than Ni, Ag migrates to cover neighboring Ni regions, therefore modifying the RP instability. Thus, we demonstrate that the phase separation length scale of an immiscible alloy can be exploited to direct the assembly of functional bimetallic alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Excitons in Hematite Fe 2 O 3 : Short-Time Dynamics from TD-DFT and Non-Adiabatic Dynamics Theories

We present a first-principles study of the short-time dynamics of excitons in hematite Fe 2 O 3 . We used time-dependent density functional theory (TD-DFT) with an underlying DFT+U treatment of electron interactions to characterize the electronic structure of excitons and nonadiabatic molecular dynamics theory (NA-MD) to determine their recombination (electronic ground-state recovery) and relaxation dynamics. Decoherence-corrected trajectory surface hopping approaches in NA-MD simulations yielded recovery times of ~1.1 to 1.8 ns and “higher-lying” exciton relaxation times of ~60 to 70 fs, in accord with experimentally derived lifetimes. With hematite phonons in the range of ~100 to 700 cm –1 , higher-lying excitons relax within one or two oscillations of the phonons before getting trapped into an electron–hole pair Exc-3 structure on the first excited state potential energy surface. This structure resembles already a pair of polarons (electron polaron plus hole polaron) with associated lattice distortions three (3) basal planes away. On longer time scales, the electron–hole bipolaronic pair hops to structures Exc-5, then Exc-7, then Exc-9, ... with the electron polaron and hole polaron separated by 5, 7, 9, ... basal planes in a process of charge separation. The largest frequency phonon ~672 cm –1 for the Exc-3 exciton structure is associated with the electron polaron moiety of the exciton. This phonon is a good candidate for giving rise to the recently observed and reported postexcitation transient IR absorption peak.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atomistic Insights into the Reactive Diffusion of CO 2 in Guanidine-Based Facilitated Transport Membranes

The pressing need to address climate change has led to significant advancements in carbon dioxide (CO 2 ) capture technologies. Notably, facilitated transport membranes (FTMs) are distinguished by their exceptional selectivity and permeance, attributed to their reversible chemical reactions with CO 2 . This study, for the first time, sheds light on the reactive diffusion mechanism of CO 2 in FTMs, utilizing 1,1,3,3-tetramethylguanidine (TMG) as a mobile carrier. Specifically, state-ofthe- art molecular dynamics (MD) simulations, augmented by a reparameterized reactive force field (ReaxFF) capable of describing atomistic interactions and reaction pathways, are conducted to investigate the transport of CO 2 in TMG. The analysis of mean squared displacement (MSD) and diffusion coefficients reveals a clear hierarchy in the mobility of reaction components. Our findings highlight a unique hopping diffusion mechanism between bicarbonate ions and TMG molecules, increasing the diffusivity of reacted CO 2 by 1.4 times. The hopping events observed not only enhance our understanding of molecular mobility but also offer a means to boost the performance of FTMs in CO 2 capture applications. Overall, this research lays the groundwork for the future design of FTMs with optimal carrier properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Modeling Framework to Predict Melting Dynamics at Microstructural Defects in TNT-HMX High Explosive Composites

Many high explosive (HE) formulations are composite materials whose microstructure is understood to impact functional characteristics. Interfaces are known to mediate the formation of hot spots that control their safety and initiation. Here, to study such processes at molecular scales, we developed all-atom force fields (FFs) for Octol, a prototypical HE formulation comprised of TNT (2,4,6-trinitrotoluene) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We extended a FF for TNT and recasted it in a form that can be readily combined with a well-established FF for HMX. The resulting FF was extensively validated against experimental results and density functional theory calculations. We applied the new combined TNT-HMX FF to predict and rank surface and interface energies, which indicate that there is an energetic driver for coarsening of microstructural grains in TNT-HMX composites. Finally, we assess the impact of several microstructural environments on the dynamic melting of TNT crystal under ultrafast thermal loading. We find that both free surfaces and planar material interfaces are effective nucleation points for TNT melting. However, MD simulations show that TNT crystal is prone to superheating by at least 50 K on subnanosecond time scales and that the degree of superheating is inversely correlated with surface and interface energy. The modeling framework presented here will enable future studies on hot spot formation processes in accident scenarios that are governed by strong coupling between microstructural interfaces, material mechanics, momentum and energy transport, phase transitions, and chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemically Mediated Au–C(sp 2 ) Anchors for Molecular Electronics

Terminal anchor groups play a key role in the stability and electronic properties of molecular junctions. Single molecule junctions typically consist of two preinstalled terminal anchors linking organic molecules to metal electrodes. Here, in this work, we show that p -terphenyl derivatives containing only a single terminal anchor show conductance features similar to junctions with two preinstalled terminal anchors. A set of p -terphenyl derivatives with one terminal anchor was prepared using automated chemical synthesis and characterized using single molecule electronics experiments, molecular dynamics (MD) simulations, bulk electrochemistry and spectroscopy, and nonequilibrium Green’s function-density functional theory (NEGF-DFT) calculations. Our results show that 4-amino- p -terphenyl (PPP) and related analogs exhibit a well-defined high conductance state that is diminished or absent in other p-terphenyl derivatives lacking a preinstalled amine terminal anchor or fluorine or methyl substitutions at the terminal para position. However, a low conductance state is observed in all amino- p -terphenyl derivatives with one preinstalled anchor due to molecular junctions formed by noncovalent dimeric π–π stacking interactions. The observed high conductance state diminishes upon the addition of reducing agents and is restored upon the addition of an oxidizing agent. Our results suggest that the high conductance state arises due to Au–C(sp 2 ) bond formation facilitated by a single electron oxidation event at the electrode surface. A series of control experiments with different anchor groups shows that primary amines play a key role in forming Au–C bonds for molecular junctions. Overall, these results suggest that Au–C bond formation gives rise to high conductance pathways in organic molecules containing only one preinstalled terminal anchor. Insights from this work can be leveraged in the design of molecular electronic devices, particularly in understanding the mechanisms of molecular binding and junction formation.

charge transport↗
DOE OSTI2022

Squeezing Out Interfacial Solvation: The Role of Hydrogen-Bonding in the Structural and Orientational Freedom of Molecular Self-Assembly

We report bioinspired membrane molecules with improved physical properties and enhanced stability can serve as functional models for conventional lipid or amphiphilic species. Importantly, these molecules can also provide new insights into emergent phenomena that manifest during self-assembly at interfaces. Here, we elucidate the structural response and mechanistic steps underlying the self-assembly of the amphiphilic, charged oligodimethylsiloxane imidazolium cation (ODMS-MIM+) at the air–aqueous interface using Langmuir trough methods with coincident surface-specific vibrational sum-frequency generation (SFG) spectroscopy. We find evidence for a new compression-induced desolvation step that precedes commonly known disordered-to-ordered phase transitions to form nanoscopic assemblies. The experimental data was supported by atomistic molecular dynamics (MD) simulations to provide a detailed mechanistic picture underlying the assembly and the role of water in these phase transitions. The sensitivity of the hydrophobic ODMS tail conformations to compression-owing to distinct water–ODMS interactions and tail–tail solvation properties-offers new strategies for the design of interfaces that can be further used to develop soft-matter electronics and low-dimensional materials using physical and chemical controls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Binding Interactions between Receptor-Binding Domain of Spike Protein and Human Angiotensin Converting Enzyme-2 in Omicron Variant

The emergence of new SARS-CoV-2 Omicron variant of concern (OV) has exacerbated the COVID-19 pandemic because of a large number of mutations in the spike protein, particularly in the receptor-binding domain (RBD), resulting in highly contagious and/or vaccine-resistant strains. Herein, we present a systematic analysis based on detailed molecular dynamics (MD) simulations in order to understand how the OV RBD mutations affect the ACE2 binding. We show that the OV RBD binds to ACE2 more efficiently and tightly predominantly because of strong electrostatic interactions, thereby promoting increased infectivity and transmissibility compared to other strains. Some of the OV RBD mutations are predicted to affect the antibody neutralization either through their role in the S-protein conformational changes, such as S371L, S373P, and S375F, or through changing its surface charge distribution, such as G339D, N440K, T478K, and E484A. Other mutations, such as K417N, G446S, and Y505H, decrease the ACE2 binding, whereas S447N, Q493R, G496S, Q498R, and N501Y tend to increase it.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Guanosine Dianions Hydrated by One to Four Water Molecules

Intermolecular interactions as those present in molecule…water complexes may profoundly influence the physico-chemical properties of molecules. Here, we carried out an experimentally-computational study on the doubly-deprotonated guanosine monophosphate…water clusters, [dGMP-2H] 2- •nH 2 O (n=1–4), using a combination of negative anion photoelectron spectroscopy (NIPES) with molecular dynamics (MD) and quantum chemical (QM) calculations. Successive addition of water molecules to [dGMP-2H] 2- increases the experimental adiabatic detachment (ADE) and vertical detachment energy (VDE) by 0.5 to 0.1 eV, depending on the cluster size. In order to choose the representative conformations, we combined MD simulations with a clustering procedure to identify low energy geometries for which ADEs and VDEs were computed at the CAM-B3LYP/6-31++G(d,p) level. Our results demonstrate that the assumed approach lead to the sound geometries and energetics of the studied micro-solvates since the calculated ADEs and VDEs are in pretty good agreement with the experimental characteristics. The evolution of hydrogen bonding with cluster size indicates possibility of the occurrence of proton transfer for clusters comprising a larger number of water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗