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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2020

How cholesterol stiffens unsaturated lipid membranes

Cholesterol is an integral component of eukaryotic cell membranes and a key molecule in controlling membrane fluidity, organization, and other physicochemical parameters. It also plays a regulatory function in antibiotic drug resistance and the immune response of cells against viruses, by stabilizing the membrane against structural damage. While it is well understood that, structurally, cholesterol exhibits a densification effect on fluid lipid membranes, its effects on membrane bending rigidity are assumed to be nonuniversal; i.e., cholesterol stiffens saturated lipid membranes, but has no stiffening effect on membranes populated by unsaturated lipids, such as 1,2-dioleoyl- sn -glycero-3-phosphocholine (DOPC). This observation presents a clear challenge to structure–property relationships and to our understanding of cholesterol-mediated biological functions. Here, using a comprehensive approach—combining neutron spin-echo (NSE) spectroscopy, solid-state deuterium NMR ( 2 H NMR) spectroscopy, and molecular dynamics (MD) simulations—we report that cholesterol locally increases the bending rigidity of DOPC membranes, similar to saturated membranes, by increasing the bilayer’s packing density. All three techniques, inherently sensitive to mesoscale bending fluctuations, show up to a threefold increase in effective bending rigidity with increasing cholesterol content approaching a mole fraction of 50%. Our observations are in good agreement with the known effects of cholesterol on the area-compressibility modulus and membrane structure, reaffirming membrane structure–property relationships. The current findings point to a scale-dependent manifestation of membrane properties, highlighting the need to reassess cholesterol’s role in controlling membrane bending rigidity over mesoscopic length and time scales of important biological functions, such as viral budding and lipid–protein interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Molecular-level design of alternative media for energy-saving pilot-scale fibrillation of nanocellulose

The outstanding mechanical properties, light weight, and biodegradability of cellulose nanofibrils (CNFs) make them promising components of renewable and sustainable next-generation reinforced composite biomaterials and bioplastics. Manufacturing CNFs at a pilot scale requires disc-refining fibrillation of dilute cellulose fibers in aqueous pulp suspensions to shear the fibers apart into their nanodimensional forms, which is, however, an energy-intensive process. Here, in this study, we used atomistic molecular dynamics (MD) simulation to examine media that might facilitate the reduction of interactions between cellulose fibers, thereby reducing energy consumption in fibrillation. The most suitable medium found by the simulations was an aqueous solution with 0.007:0.012 wt.% NaOH:urea, and indeed this was found in pilot-scale experiments to reduce the fibrillation energy by ~21% on average relative to water alone. The NaOH:urea-mediated CNFs have similar crystallinity, morphology, and mechanical strength to those formed in water. The NaOH and urea act synergistically on CNFs to aid fibrillation but at different length scales. NaOH deprotonates hydroxyl groups leading to mesoscale electrostatic repulsion between fibrils, whereas urea forms hydrogen bonds with protonated hydroxyl groups thus disrupting interfibril hydrogen bonds. This suggests a general mechanism in which an aqueous medium that contains a strong base and a small organic molecule acting as a hydrogen-bond acceptor and/or donor may be effectively employed in materials processes where dispersion of deprotonable polymers is required. The study demonstrates how atomic-detail computer simulation can be integrated with pilot-scale experiments in the rational design of materials processes for the circular bioeconomy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling prebiotic chemistries with quantum accuracy at classical costs

Molecular Dynamics (MD) simulations using classical force-fields are commonly employed in numerous scientific investigations. However, many natural processes involve bond breaking and quantum forces. This complexity is compounded by the presence of multiple competing length and timescales. For example, accurately modeling the thermodynamics and dynamics of a chemical reaction requires accounting for the concerted movements of numerous solvent molecules and ions with their own fast or slow timescales. While widely used static Density Functional Theory (DFT) calculations at 0 temperature can be beneficial for such investigations, they do not account for dynamics, and lack precision in describing the molecular environments. They particularly fail at correct, rigorous treatments of finite-temperature fluctuations, and thus generalization to experimentally relevant conditions. In PNAS Benayad et al develop a scalable, generalizable approach for designing Neural Network Potentials (NNPs) that can handle chemical reactivity in solvated systems with quantum accuracy at classical costs. Specifically, they study phosphoester bond formation and rupture, which is fundamentally relevant to the Phosphorus-Oxygen bond formation central to life, and especially for the RNA world hypothesis. The framework developed here has the potential to generalize to different chemical reactions of energy and biological relevance.

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

Effects of interlayer spacing and applied pressure on the lanthanide transport in MoS 2 -based two-dimensional channels

Rare-earth elements (REEs) are critical to modern industry but difficult to separate due to their subtle and monotonic changes in physicochemical properties. MoS 2 -based two-dimensional (2D) materials offer novel opportunities for enhancing REE separation, exhibiting a distinct volcano-shaped transport performance distribution that peaks at Sm 3+ . However, the specific contributions of thermodynamic and kinetic factors to ion transport within 2D confinement remain unclear. In this study, we conducted a series of non-equilibrium all-atom molecular dynamics (MD) simulations to explore the effects of interlayer spacing and external pressure on the transport of lanthanide ions in Å-scale acetate functionalized 2D MoS 2 (MoS 2 -COOH) channels. We examined ion entry and permeation rates, water flux, dehydration, and binding modes. The simulation results reveal that the transport trends of lanthanide ions are jointly driven by the dehydration degree and the relative-binding strengths of ions to water and to the acetate within the 2D channels. Notably, the dehydration pattern of lanthanide ions during permeation is closely linked to kinetic factors. Overall, this study provides a detailed atomistic understanding of the mechanisms underlying lanthanide ion transport under confinement. These findings point to the significant potential for tuning confinement and chemical functionalization within Å-scale channels for more efficient REE separation.

2D channels↗
DOE OSTI2020

Quasi-equilibrium predictions of water desorption kinetics from rapidly-heated metal oxide surfaces

Controlling sub-microsecond desorption of water and other impurities from electrode surfaces at high heating rates is crucial for pulsed power applications. Despite the short time scales involved, quasi-equilibrium ideas based on transition state theory (TST) and Arrhenius temperature dependence have been widely applied to fit desorption activation free energies. Here, we apply molecular dynamics (MD) simulations in conjunction with equilibrium potential-of-mean-force (PMF) techniques to directly compute the activation free energies (ΔG*) associated with desorption of intact water molecules from Fe 2 O 3 and Cr 2 O 3 (0001) surfaces. The desorption free energy profiles are diffuse, without maxima, and have substantial dependences on temperature and surface water coverage. Incorporating the predicted ΔG* into an analytical form gives rate equations that are in reasonable agreement with non-equilibrium molecular dynamics desorption simulations. We also show that different ΔG* analytical functional forms which give similar predictions at a particular heating rate can yield desorption times that differ by up to a factor of four or more when the ramp rate is extrapolated by 8 orders of magnitude. This highlights the importance of constructing a physically-motivated ΔG* functional form to predict fast desorption kinetics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Interplay of folded domains and the disordered low-complexity domain in mediating hnRNPA1 phase separation

Abstract Liquid–liquid phase separation underlies the membrane-less compartmentalization of cells. Intrinsically disordered low-complexity domains (LCDs) often mediate phase separation, but how their phase behavior is modulated by folded domains is incompletely understood. Here, we interrogate the interplay between folded and disordered domains of the RNA-binding protein hnRNPA1. The LCD of hnRNPA1 is sufficient for mediating phase separation in vitro. However, we show that the folded RRM domains and a folded solubility-tag modify the phase behavior, even in the absence of RNA. Notably, the presence of the folded domains reverses the salt dependence of the driving force for phase separation relative to the LCD alone. Small-angle X-ray scattering experiments and coarse-grained MD simulations show that the LCD interacts transiently with the RRMs and/or the solubility-tag in a salt-sensitive manner, providing a mechanistic explanation for the observed salt-dependent phase separation. These data point to two effects from the folded domains: (i) electrostatically-mediated interactions that compact hnRNPA1 and contribute to phase separation and (ii) increased solubility at higher ionic strengths mediated by the folded domains. The interplay between disordered and folded domains can modify the dependence of phase behavior on solution conditions and can obscure signatures of physicochemical interactions underlying phase separation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Cryo neutron crystallography demonstrates influence of RNA 2'-OH orientation on conformation, sugar pucker and water structure

The ribose 2'-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2'-OH orientation (H2'-C2'-O2'-HO2' torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2'-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3', 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2'-OD backbone-orientation. GAGA-tetraloop riboses display a 2'-OD base-orientation. An atypical C2'-endo sugar pucker is strictly correlated with a 2'-OD sugar-orientation. Neutrons reveal the strong preference of the 2'-OH to donate in H-bonds and that 2'-OH orientation affects both backbone geometry and ribose pucker. We discuss 2'-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Regioselectivity mechanism of the Thunbergia alata Δ 6 -16:0-acyl carrier protein desaturase

Plant plastidial acyl–acyl carrier protein (ACP) desaturases are a soluble class of diiron-containing enzymes that are distinct from the diiron-containing integral membrane desaturases found in plants and other organisms. The archetype of this class is the stearoyl-ACP desaturase which converts stearoyl-ACP into oleoyl (18:1Δ 9 cis )-ACP. Several variants expressing distinct regioselectivity have been described including a Δ 6 -16:0-ACP desaturase from black-eyed Susan vine ( Thunbergia alata ). We solved a crystal structure of the T. alata desaturase at 2.05 Å resolution. Using molecular dynamics (MD) simulations, we identified a low-energy complex between 16:0-ACP and the desaturase that would position C6 and C7 of the acyl chain adjacent to the diiron active site. The model complex was used to identify mutant variants that could convert the T. alata Δ 6 desaturase to Δ 9 regioselectivity. Additional modeling between ACP and the mutant variants confirmed the predicted regioselectivity. Furthermore, to validate the in-silico predictions, we synthesized two variants of the T. alata desaturase and analyzed their reaction products using gas chromatography-coupled mass spectrometry. Assay results confirmed that mutants designed to convert T. alata Δ 6 to Δ 9 selectivity exhibited the predicted changes. In complementary experiments, variants of the castor desaturase designed to convert Δ 9 to Δ 6 selectivity lost some of their Δ9 desaturation ability and gained the ability to desaturate at the Δ 6 position. The computational workflow for revealing the mechanistic understanding of regioselectivity presented herein lays a foundation for designing acyl-ACP desaturases with novel selectivities to increase the diversity of monoenes available for bioproduct applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

First-principles investigation of structure-property relationships in stable and metastable MXenes

Understanding the structure–property relationships in layered transition-metal carbides or nitrides, known as MXenes, is of critical importance for their rational design, synthesis, and application. However, the vast chemical and structural diversity of MXenes, stemming from their wide range of M and X elements, surface terminations, and different atomic coordination environments, makes it challenging to clearly understand these structure–property relationships. Here, in this work, we perform first-principles density functional theory (DFT) calculations and molecular dynamics (MD) simulations to comprehensively investigate the stability and a variety of physical properties of MXenes with different coordination environments. Using Ti- and Mo-based carbide MXenes as model systems, energetic calculations reveal that Ti-based MXenes are most stable in octahedral coordination, whereas Mo-based MXenes preferentially adopt prismatic coordination. This fundamental difference in preferred atomic coordination gives rise to markedly distinct properties between these two systems as a function of the fraction of octahedral and prismatic sites. For instance, the in-plane stiffness of Ti-based MXenes increases as octahedral coordination becomes dominant, but it decreases in the Mo-based MXenes under the same conditions. Additional stability analyses based on mechanical, lattice-dynamical, and temperature-dependent thermodynamic properties demonstrate that many metastable MXenes not only satisfy the strict stability criteria but can also undergo phase transitions among different structures and even become stabilized at elevated temperatures. Although surface terminations, such as F and O atoms, do not alter the energetic ordering or the overall stiffness trends among stable and metastable MXenes, they influence other material properties. For instance, O termination can induce semiconducting behavior in both stable and metastable Ti 2 ⁢CO 2 MXenes. This study significantly advances the fundamental understanding of structure–property relationships in MXenes and provides valuable guidance for developing coordination-based design principles to precisely engineer MXenes with improved properties.

Oyeniran, Noah [University of Alabama, Tuscaloosa,↗
DOE OSTI2021

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
DOE OSTI2025

Generalized approach for rapid entropy calculation of liquids and solids

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity autocorrelation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unravelling the Glycan Code: Molecular Dynamics and Quantum Chemistry Reveal How O-Glycan Functional Groups Govern OgpA Selectivity in Mucin Degradation by Akkermansia muciniphila

Mucins, heavily O-glycosylated glycoproteins, are a key component of mucus, and certain gut microbiota, including Akkermansia muciniphila , can utilise mucin glycans as a carbon source. Akkermansia muciniphila produces the O-glycopeptidase enzyme OgpA, which cleaves peptide bonds at the N-terminus of serine (Ser) or threonine (Thr) residues carrying O-glycan substitutions, with selectivity influenced by the O-glycan functional groups. Using molecular dynamics (MD) simulations and quantum chemistry calculations, we explored how different O-glycan groups affect OgpA's selectivity. Our results show that peptides bind to the enzyme via hydrogen bonds, π–π interactions, van der Waals forces and electrostatic interactions, with key residues, including Tyr90, Val138, Gly176, Tyr210 and Glu91, playing important roles. The primary determinant of selectivity is the interaction between the peptide's functional group and the enzyme's binding cavity, while peptide–enzyme interface interactions are secondary. Quantum chemistry calculations reveal that OgpA prefers peptides with a lower electrophilic character. This study provides new insights into mucin degradation by gut microbiota enzymes, advancing our understanding of this critical biological process.

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

Enhancing Thermal Transport in Polymeric Composites Via Engineered Noncovalent Filler–Polymer Interactions

Understanding thermal transport mechanisms in polymeric composites allows us to expand the boundaries of thermal conductivity in them, either increasing it for more efficient heat dissipation or decreasing it for better thermal insulation. But, these mechanisms are not fully understood. Systematic experimental investigations remain limited. Practical strategies to tune the interfacial thermal resistance (ITR) between fillers and polymers and the thermal conductivity of composites remain elusive. Here, we studied the thermal transport in representative polymer composites, using polyethylene (PE) or polyaniline (PANI) as matrices and graphite as fillers. PANI, with aromatic rings in its backbone, interacts with graphite through strong noncovalent π–π stacking interactions, whereas PE lacks such interactions. We can then quantify how π–π stacking interactions between graphite and polymers enhance thermal transport in composites. PE/graphite and PANI/graphite composites with the same 1.5% filler volume fractions show a ∼22.82% and ∼34.85% enhancement in thermal conductivity compared to pure polymers, respectively. Calculated ITRs in PE/graphite and PANI/graphite are ~6 x 10 -8 m 2 KW -1 and ~1 x 10 -8 m 2 KW -1 , respectively, highlighting how π–π stacking interactions reduce ITR. Molecular dynamics (MD) simulations suggest that π–π stacking interactions between PANI chains and graphite surfaces enhance alignment of PANI's aromatic rings with graphite surfaces. This allows more carbon atoms from PANI chains to interact with graphite surfaces at a shorter distance compared to PE chains. Finally, our work indicates that tuning the π–π stacking interactions between polymers and fillers is an effective approach to reduce the ITR and enhance the thermal conductivity of composites.

atoms↗
DOE OSTI2025

A molecular conveyor belt-associated protein controls the rotational direction of the bacterial type 9 secretion system

ABSTRACT Many bacteria utilize the type 9 secretion system (T9SS) for gliding motility, surface colonization, and pathogenesis. This dual-function motor supports both gliding motility and protein secretion, where rotation of the T9SS plays a central role. Fueled by the energy of the stored proton motive force and transmitted through the torque of membrane-anchored stator units, the rotary T9SS propels an adhesin-coated conveyor belt along the bacterial outer membrane like a molecular snowmobile, thereby enabling gliding motion. However, the mechanisms controlling the rotational direction and gliding motility of T9SS remain elusive. Shedding light on this mechanism, we find that in the gliding bacteriumFlavobacterium johnsoniae, deletion of the C-terminus of the conveyor belt-associated protein GldJ controls and, in fact, reverses the rotational direction of T9SS from counterclockwise (CCW) to clockwise (CW). This suggests that the interface between the conveyor belt-associated protein GldJ and the T9SS ring protein GldK plays an important role in controlling the directionality of T9SS, potentially by modulating its interaction with the stator complex GldLM, which drives motor rotation. Combined with MD simulation of the T9SS stator units GldLM, we suggest a “tri-component gearset” model where GldJ controls the rotational direction of its driver, the T9SS, thus providing adaptive sensory feedback to influence the motility of the gliding bacterium. IMPORTANCE The type 9 secretion system (T9SS) is fundamental to bacterial gliding motility, pathogenesis, and surface colonization. Our findings reveal that the C-terminal region of the conveyor belt-associated protein GldJ functions as a molecular switch which is capable of reversing the rotational direction of T9SS. Through the coordinated actions of the T9SS stator units (akin to a driving motor), the GldK ring (the gear that converts rotational energy into linear movement), and GldJ, this machinery forms a smart conveyor belt system reminiscent of flexible or cognitive mechanical conveyors. Such advanced conveyors can alter their direction to adapt to shifting demands. Here, we show that the bacterial T9SS similarly adjusts its rotational bias based on feedback from the conveyor belt-associated protein GldJ. This dual-role feedback mechanism underscores an evolved, controllable biological snowmobile, offering new avenues for studying how bacteria fine-tune motility in dynamic environments.

Microbiology↗
DOE OSTI2021

Transport Phenomena in Low Temperature Lithium-Ion Battery Electrolytes

Lithium-ion batteries face low temperature performance issues, limiting the adoption of technologies ranging from electric vehicles to stationary grid storage. This problem is thought to be exacerbated by slow transport within the electrolyte, which in turn may be influenced by ion association, solvent viscosity, and cation transference number. How these factors collectively impact low temperature transport phenomena, however, remains poorly understood. Here we show using all-atom classical molecular dynamics (MD) simulations that the dominant factor influencing low temperature transport in LP57 (1 M LiPF 6 in 3:7 ethylene carbonate (EC)/ethyl methyl carbonate (EMC)) is solvent viscosity, rather than ion aggregation or cation transference number. We find that ion association decreases with decreasing temperature, while the cation transference number is positive and roughly independent of temperature. In an effort to improve low temperature performance, we introduce γ -butyrolactone (GBL) as a low viscosity co-solvent to explore two alternative formulations: 1 M LiPF 6 in 15:15:70 EC/GBL/EMC and 3:7 GBL/EMC. While GBL reduces solution viscosity, its low dielectric constant results in increased ion pairing, yielding neither improved bulk ionic conductivity nor appreciably altered ion transport mechanisms. We expect that these results will enhance understanding of low temperature transport and inform the development of superior electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2020

Permeability of the HIV-1 capsid to metabolites modulates viral DNA synthesis

Reverse transcription, an essential event in the HIV-1 life cycle, requires deoxynucleotide triphosphates (dNTPs) to fuel DNA synthesis, thus requiring penetration of dNTPs into the viral capsid. The central cavity of the capsid protein (CA) hexamer reveals itself as a plausible channel that allows the passage of dNTPs into assembled capsids. Nevertheless, the molecular mechanism of nucleotide import into the capsid remains unknown. Employing all-atom molecular dynamics (MD) simulations, we established that cooperative binding between nucleotides inside a CA hexamer cavity results in energetically favorable conditions for passive translocation of dNTPs into the HIV-1 capsid. Furthermore, binding of the host cell metabolite inositol hexakisphosphate (IP 6 ) enhances dNTP import, while binding of synthesized molecules like benzenehexacarboxylic acid (BHC) inhibits it. The enhancing effect on reverse transcription by IP 6 and the consequences of interactions between CA and nucleotides were corroborated using atomic force microscopy, transmission electron microscopy, and virological assays. Collectively, our results provide an atomistic description of the permeability of the HIV-1 capsid to small molecules and reveal a novel mechanism for the involvement of metabolites in HIV-1 capsid stabilization, nucleotide import, and reverse transcription.

60 APPLIED LIFE SCIENCES↗