Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein dynamics”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Mass spectrometry-based technologies for probing the 3D world of plant proteins

Abstract Over the past two decades, mass spectrometric (MS)-based proteomics technologies have facilitated the study of signaling pathways throughout biology. Nowhere is this needed more than in plants, where an evolutionary history of genome duplications has resulted in large gene families involved in posttranslational modifications and regulatory pathways. For example, at least 5% of the Arabidopsis thaliana genome (ca. 1,200 genes) encodes protein kinases and protein phosphatases that regulate nearly all aspects of plant growth and development. MS-based technologies that quantify covalent changes in the side-chain of amino acids are critically important, but they only address one piece of the puzzle. A more crucially important mechanistic question is how noncovalent interactions—which are more difficult to study—dynamically regulate the proteome’s 3D structure. The advent of improvements in protein 3D technologies such as cryo-electron microscopy, nuclear magnetic resonance, and X-ray crystallography has allowed considerable progress to be made at this level, but these methods are typically limited to analyzing proteins, which can be expressed and purified in milligram quantities. Newly emerging MS-based technologies have recently been developed for studying the 3D structure of proteins. Importantly, these methods do not require protein samples to be purified and require smaller amounts of sample, opening the wider proteome for structural analysis in complex mixtures, crude lysates, and even in intact cells. These MS-based methods include covalent labeling, crosslinking, thermal proteome profiling, and limited proteolysis, all of which can be leveraged by established MS workflows, as well as newly emerging methods capable of analyzing intact macromolecules and the complexes they form. In this review, we discuss these recent innovations in MS-based “structural” proteomics to provide readers with an understanding of the opportunities they offer and the remaining challenges for understanding the molecular underpinnings of plant structure and function.

Plant Sciences↗
DOE OSTI2022

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Beneath the surface: Unsolved questions in soil virus ecology

Soil virus ecology is an exciting but still nascent field of research in soil microbiology. While there has been a recent surge in soil virus research studies, many fundamental questions remain unanswered, and a range of technical and bioinformatic challenges need to be overcome. In this perspective article, we present a series of key questions that highlight fruitful research areas for ongoing and future efforts. These include describing the challenges involved in understanding soil viral abundance and activity, spatiotemporal dynamics, life strategy prevalence, virus-mediated biogeochemical impacts, viral protein function, host prediction, and soil RNA virus discovery. In the near term, combining approaches (e.g., cultivation-based, meta-omics, biogeochemical, experimental, and bioinformatic) will be key to assessing the ecological and biogeochemical impacts of soil viruses from the microscopic to the field and global scales. Still, we stress that results must be tempered by current methodological limitations and highlight knowledge gaps that are most pressing to fill via new methods or measurements, such as the prevalence of different viral replication strategies across soils, the fate of microbial necromass carbon after viral lysis, the frequency of virus-host encounters that do not lead to successful infections yet could be bioinformatically mistaken as infections, and the diversity and ecological impacts of RNA viruses in soil.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Enzymes in 3D: Synthesis, remodelling, and hydrolysis of cell wall (1,3;1,4)-β-glucans

Abstract Recent breakthroughs in structural biology have provided valuable new insights into enzymes involved in plant cell wall metabolism. More specifically, the molecular mechanism of synthesis of (1,3;1,4)-β-glucans, which are widespread in cell walls of commercially important cereals and grasses, has been the topic of debate and intense research activity for decades. However, an inability to purify these integral membrane enzymes or apply transgenic approaches without interpretative problems associated with pleiotropic effects has presented barriers to attempts to define their synthetic mechanisms. Following the demonstration that some members of the CslF sub-family of GT2 family enzymes mediate (1,3;1,4)-β-glucan synthesis, the expression of the corresponding genes in a heterologous system that is free of background complications has now been achieved. Biochemical analyses of the (1,3;1,4)-β-glucan synthesized in vitro, combined with 3-dimensional (3D) cryogenic-electron microscopy and AlphaFold protein structure predictions, have demonstrated how a single CslF6 enzyme, without exogenous primers, can incorporate both (1,3)- and (1,4)-β-linkages into the nascent polysaccharide chain. Similarly, 3D structures of xyloglucan endo-transglycosylases and (1,3;1,4)-β-glucan endo- and exohydrolases have allowed the mechanisms of (1,3;1,4)-β-glucan modification and degradation to be defined. X-ray crystallography and multi-scale modeling of a broad specificity GH3 β-glucan exohydrolase recently revealed a previously unknown and remarkable molecular mechanism with reactant trajectories through which a polysaccharide exohydrolase can act with a processive action pattern. The availability of high-quality protein 3D structural predictions should prove invaluable for defining structures, dynamics, and functions of other enzymes involved in plant cell wall metabolism in the immediate future.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural analysis of wild-type and Val120Thr mutant Candida boidinii formate dehydrogenase by X-ray crystallography

Candida boidinii NAD + -dependent formate dehydrogenase (CbFDH) has gained significant attention for its potential application in the production of biofuels and various industrial chemicals from inorganic carbon dioxide. The present study reports the atomic X-ray crystal structures of wild-type CbFDH at cryogenic and ambient temperatures, as well as that of the Val120Thr mutant at cryogenic temperature, determined at the Turkish Light Source ‘Turkish DeLight’. The structures reveal new hydrogen bonds between Thr120 and water molecules in the active site of the mutant CbFDH, suggesting increased stability of the active site and more efficient electron transfer during the reaction. Further experimental data is needed to test these hypotheses. Collectively, these findings provide invaluable insights into future protein-engineering efforts that could potentially enhance the efficiency and effectiveness of CbFDH.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Molecular recognition of a host protein by NS1 of pandemic and seasonal influenza A viruses

The 1918 influenza A virus (IAV) caused the most severe flu pandemic in recorded human history. Nonstructural protein 1 (NS1) is an important virulence factor of the 1918 IAV. NS1 antagonizes host defense mechanisms through interactions with multiple host factors. One pathway by which NS1 increases virulence is through the activation of phosphoinositide 3-kinase (PI3K) by binding to its p85β subunit. Here we present the mechanism underlying the molecular recognition of the p85β subunit by 1918 NS1. Using X-ray crystallography, we determine the structure of 1918 NS1 complexed with p85β of human PI3K. We find that the 1918 NS1 effector domain (1918 NS1 ED ) undergoes a conformational change to bind p85β. Using NMR relaxation dispersion and molecular dynamics simulation, we identify that free 1918 NS1 ED exists in a dynamic equilibrium between p85β-binding–competent and –incompetent conformations in the submillisecond timescale. Moreover, we discover that NS1 ED proteins of 1918 (H1N1) and Udorn (H3N2) strains exhibit drastically different conformational dynamics and binding kinetics to p85β. These results provide evidence of strain-dependent conformational dynamics of NS1. Using kinetic modeling based on the experimental data, we demonstrate that 1918 NS1 ED can result in the faster hijacking of p85β compared to Ud NS1 ED , although the former has a lower affinity to p85β than the latter. Our results suggest that the difference in binding kinetics may impact the competition with cellular antiviral responses for the activation of PI3K. We anticipate that our findings will increase the understanding of the strain-dependent behaviors of influenza NS1 proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Macroscale structural changes of thylakoid architecture during high light acclimation in Chlamydomonas reinhardtii

Abstract Photoprotection mechanisms are ubiquitous among photosynthetic organisms. The photoprotection capacity of the green alga Chlamydomonas reinhardtii is correlated with protein levels of stress-related light-harvesting complex (LHCSR) proteins, which are strongly induced by high light (HL). However, the dynamic response of overall thylakoid structure during acclimation to growth in HL has not been fully understood. Here, we combined live-cell super-resolution microscopy and analytical membrane subfractionation to investigate macroscale structural changes of thylakoid membranes during HL acclimation in Chlamydomonas . Subdiffraction-resolution live-cell imaging revealed that the overall thylakoid structures became thinned and shrunken during HL acclimation. The stromal space around the pyrenoid also became enlarged. Analytical density-dependent membrane fractionation indicated that the structural changes were partly a consequence of membrane unstacking. The analysis of both an LHCSR loss-of-function mutant , npq4 lhcsr1 , and a regulatory mutant that over-expresses LHCSR, spa1-1 , showed that structural changes occurred independently of LHCSR protein levels, demonstrating that LHCSR was neither necessary nor sufficient to induce the thylakoid structural changes associated with HL acclimation. In contrast, stt7-9 , a mutant lacking a kinase of major light-harvesting antenna proteins, had a slower thylakoid structural response to HL relative to all other lines tested but still showed membrane unstacking. These results indicate that neither LHCSR- nor antenna-phosphorylation-dependent HL acclimation are required for the observed macroscale structural changes of thylakoid membranes in HL conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A dynamic charge-charge interaction modulates PP2A:B56 substrate recruitment

The recruitment of substrates by the ser/thr protein phosphatase 2A (PP2A) is poorly understood, limiting our understanding of PP2A-regulated signaling. Recently, the first PP2A:B56 consensus binding motif, LxxIxE, was identified. However, most validated LxxIxE motifs bind PP2A:B56 with micromolar affinities, suggesting that additional motifs exist to enhance PP2A:B56 binding. Here, we report the requirement of a positively charged motif in a subset of PP2A:B56 interactors, including KIF4A, to facilitate B56 binding via dynamic, electrostatic interactions. Using molecular and cellular experiments, we show that a conserved, negatively charged groove on B56 mediates dynamic binding. We also discovered that this positively charged motif, in addition to facilitating KIF4A dephosphorylation, is essential for condensin I binding, a function distinct and exclusive from PP2A-B56 binding. Together, these results reveal how dynamic, charge-charge interactions fine-tune the interactions mediated by specific motifs, providing a new framework for understanding how PP2A regulation drives cellular signaling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Staphylococcus aureus FtsZ and PBP4 bind to the conformationally dynamic N-terminal domain of GpsB

In the Firmicutes phylum, GpsB is a membrane associated protein that coordinates peptidoglycan synthesis with cell growth and division. Although GpsB has been studied in several bacteria, the structure, function, and interactome of Staphylococcus aureus GpsB is largely uncharacterized. To address this knowledge gap, we solved the crystal structure of the N-terminal domain of S. aureus GpsB, which adopts an atypical, asymmetric dimer, and demonstrates major conformational flexibility that can be mapped to a hinge region formed by a three-residue insertion exclusive to Staphylococci. When this three-residue insertion is excised, its thermal stability increases, and the mutant no longer produces a previously reported lethal phenotype when overexpressed in Bacillus subtilis. In S. aureus, we show that these hinge mutants are less functional and speculate that the conformational flexibility imparted by the hinge region may serve as a dynamic switch to fine-tune the function of the GpsB complex and/or to promote interaction with its various partners. Furthermore, we provide the first biochemical, biophysical, and crystallographic evidence that the N-terminal domain of GpsB binds not only PBP4, but also FtsZ, through a conserved recognition motif located on their C-termini, thus coupling peptidoglycan synthesis to cell division. Taken together, the unique structure of S. aureus GpsB and its direct interaction with FtsZ/PBP4 provide deeper insight into the central role of GpsB in S. aureus cell division.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Distinct binding conformations of epinephrine with α- and β-adrenergic receptors

Abstract Agonists targeting α 2 -adrenergic receptors (ARs) are used to treat diverse conditions, including hypertension, attention-deficit/hyperactivity disorder, pain, panic disorders, opioid and alcohol withdrawal symptoms, and cigarette cravings. These receptors transduce signals through heterotrimeric Gi proteins. Here, we elucidated cryo-EM structures that depict α 2A -AR in complex with Gi proteins, along with the endogenous agonist epinephrine or the synthetic agonist dexmedetomidine. Molecular dynamics simulations and functional studies reinforce the results of the structural revelations. Our investigation revealed that epinephrine exhibits different conformations when engaging with α-ARs and β-ARs. Furthermore, α 2A -AR and β 1 -AR (primarily coupled to Gs, with secondary associations to Gi) were compared and found to exhibit different interactions with Gi proteins. Notably, the stability of the epinephrine–α 2A -AR–Gi complex is greater than that of the dexmedetomidine–α 2A -AR–Gi complex. These findings substantiate and improve our knowledge on the intricate signaling mechanisms orchestrated by ARs and concurrently shed light on the regulation of α-ARs and β-ARs by epinephrine.

Biochemistry & Molecular Biology↗
DOE OSTI2023

Energy dissipation efficiency in the CP43 assembly intermediate complex of photosystem II

Photosystem II in oxygenic organisms is a large membrane bound rapidly turning over pigment protein complex. During its biogenesis, multiple assembly intermediates are formed, including the CP43-preassembly complex (pCP43). Here, to understand the energy transfer dynamics in pCP43, we first engineered a His-tagged version of the CP43 in a CP47-less strain of the cyanobacterium Synechocystis 6803. Isolated pCP43 from this engineered strain was subjected to advanced spectroscopic analysis to evaluate its excitation energy dissipation characteristics. These included measurements of steady-state absorption and fluorescence emission spectra for which correlation was tested with Stepanov relation. Comparison of fluorescence excitation and absorptance spectra determined that efficiency of energy transfer from β-carotene to chlorophyll a is 39 %. Time-resolved fluorescence images of pCP43-bound Chl a were recorded on streak camera, and fluorescence decay dynamics were evaluated with global fitting. These demonstrated that the decay kinetics strongly depends on temperature and buffer used to disperse the protein sample and fluorescence decay lifetime was estimated in 3.2–5.7 ns time range, depending on conditions. The pCP43 complex was also investigated with femtosecond and nanosecond time-resolved absorption spectroscopy upon excitation of Chl a and β-carotene to reveal pathways of singlet excitation relaxation/decay, Chl a triplet dynamics and Chl a → β-carotene triplet state sensitization process. The latter demonstrated that Chl a triplet in the pCP43 complex is not efficiently quenched by carotenoids. Finally, detailed kinetic analysis of the rise of the population of β-carotene triplets determined that the time constant of the carotenoid triplet sensitization is 40 ns.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

HURP regulates Kif18A recruitment and activity to synergistically control microtubule dynamics

Abstract During mitosis, microtubule dynamics are regulated to ensure proper alignment and segregation of chromosomes. The dynamics of kinetochore-attached microtubules are regulated by hepatoma-upregulated protein (HURP) and the mitotic kinesin-8 Kif18A, but the underlying mechanism remains elusive. Using single-molecule imaging in vitro, we demonstrate that Kif18A motility is regulated by HURP. While sparse decoration of HURP activates the motor, higher concentrations hinder processive motility. To shed light on this behavior, we determine the binding mode of HURP to microtubules using cryo-EM. The structure helps rationalize why HURP functions as a microtubule stabilizer. Additionally, HURP partially overlaps with the microtubule-binding site of the Kif18A motor domain, indicating that excess HURP inhibits Kif18A motility by steric exclusion. We also observe that HURP and Kif18A function together to suppress dynamics of the microtubule plus-end, providing a mechanistic basis for how they collectively serve in microtubule length control.

Science & Technology - Other Topics↗
DOE OSTI2023

PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid

Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β -sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Finally, fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗
DOE OSTI2023

Nuclear–Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics

Simulating the nuclear–electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear–electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear–electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron–proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. Furthermore, these examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Using enzymes to tame nitrogen-centred radicals for enantioselective hydroamination

The formation of C–N bonds—of great importance to the pharmaceutical industry—can be facilitated enzymatically using nucleophilic and nitrene transfer mechanisms. However, neither natural nor engineered enzymes are known to generate and control nitrogen-centred radicals, which serve as valuable species for C–N bond formation. Here we use flavin-dependent ‘ene’-reductases with an exogenous photoredox catalyst to selectively generate amidyl radicals within the protein active site. These enzymes are engineered through directed evolution to catalyse 5-exo, 6-endo, 7-endo, 8-endo, and intermolecular hydroamination reactions with high levels of enantioselectivity. Mechanistic studies suggest that radical initiation occurs via an enzyme-gated mechanism, where the protein thermodynamically activates the substrate for reduction by the photocatalyst. Molecular dynamics studies indicate that the enzymes bind substrates using non-canonical binding interactions, which may serve as a handle to further manipulate reactivity. Furthermore, this approach demonstrates the versatility of these enzymes for controlling the reactivity of high-energy radical intermediates and highlights the opportunity for synergistic catalyst strategies to unlock previously inaccessible enzymatic functions.

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

Neutron spin echo spectroscopy

Neutron spin echo (NSE) spectroscopy is a powerful technique used to probe the internal dynamics and diffusion processes within matter, allowing researchers to quantify how materials respond to changes in external conditions. To fully understand the bulk properties of a material, it is often crucial to investigate dynamic processes occurring on nanometre length scales and nanosecond timescales. NSE spectroscopy provides an exceptional capability to measure such dynamics across a broad spectrum of condensed matter systems, ranging from proteins and quantum magnets to catalysts. In this Primer, we present an overview of NSE spectroscopy, highlighting the instruments used, methods of data collection and representative applications in both soft-matter and hard-matter sciences.

Czakkel, Orsolya [Institut Max von Laue-Paul Lange↗
DOE OSTI2021

Probing remdesivir nucleotide analogue insertion to SARS-CoV-2 RNA dependent RNA polymerase in viral replication

Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. Here in this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (~1.5 kcal mol -1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick base pairing with the template nt, and inserts into the active site with a medium low free energy barrier (~2.6 kcal mol -1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.

36 MATERIALS SCIENCE↗