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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

HIV influences microtubule associated protein-2: potential marker of HIV-associated neurocognitive disorders

Postmortem brains of patients diagnosed with HIV-1-associated neurocognitive disorders (HAND) exhibit loss of dendrites. However, the mechanisms by which synapses are damaged are not fully understood. Dendrite length and remodeling occurs via microtubules, the dynamics of which are regulated by microtubule-binding proteins, including microtubule-associated protein 2 (MAP2). The HIV protein gp120 is neurotoxic and interferes with neuronal microtubules. We measured MAP2 concentrations in human cerebrospinal fluid (CSF) and MAP2 immunoreactivity in rat cortical neurons exposed to HIV and gp120. First, we examined whether HIV affects MAP2 levels by analyzing the CSF of 27 persons living with HIV (PLH) whose neurocognitive performance had been characterized. We then used rat cortical neurons to study the mechanisms of HIV-mediated dendritic loss. PLH who had HAND had greater MAP2 concentrations within the CSF than cognitive normal PLH. In cortical neurons, the deleterious effect of HIV on MAP2-positive dendrites occurred through a gp120-mediated mechanism. The neurotoxic effect of HIV was blocked by a CCR5 antagonist and prevented by Helix-A, a peptide that displaces gp120 from binding to microtubules, conjugated to a nanolipoprotein particle delivery platform. Our findings support that HIV at least partially effects its neurotoxicity via neuronal cytoskeleton modifications and provide evidence of a new therapeutic compound that could be used to prevent the HIV-associated neuropathology.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Thermodynamic and Kinetic Modeling of Co-utilization of Glucose and Xylose for 2,3-BDO Production by Zymomonas mobilis

Prior engineering of the ethanologen Zymomonas mobilis has enabled it to metabolize xylose and to produce 2,3-butanediol (2,3-BDO) as a dominant fermentation product. When co-fermenting with xylose, glucose is preferentially utilized, even though xylose metabolism generates ATP more efficiently during 2,3-BDO production on a BDO-mol basis. To gain a deeper understanding of Z. mobilis metabolism, we first estimated the kinetic parameters of the glucose facilitator protein of Z. mobilis by fitting a kinetic uptake model, which shows that the maximum transport capacity of glucose is seven times higher than that of xylose, and glucose is six times more affinitive to the transporter than xylose. With these estimated kinetic parameters, we further compared the thermodynamic driving force and enzyme protein cost of glucose and xylose metabolism. It is found that, although 20% more ATP can be yielded stoichiometrically during xylose utilization, glucose metabolism is thermodynamically more favorable with 6% greater cumulative Gibbs free energy change, more economical with 37% less enzyme cost required at the initial stage and sustains the advantage of the thermodynamic driving force and protein cost through the fermentation process until glucose is exhausted. Glucose-6-phosphate dehydrogenase (g6pdh), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and phosphoglycerate mutase (pgm) are identified as thermodynamic bottlenecks in glucose utilization pathway, as well as two more enzymes of xylose isomerase and ribulose-5-phosphate epimerase in xylose metabolism. Acetolactate synthase is found as potential engineering target for optimized protein cost supporting unit metabolic flux. Pathway analysis was then extended to the core stoichiometric matrix of Z. mobilis metabolism. Growth was simulated by dynamic flux balance analysis and the model was validated showing good agreement with experimental data. Dynamic FBA simulations suggest that a high agitation is preferable to increase 2,3-BDO productivity while a moderate agitation will benefit the 2,3-BDO titer. Taken together, this work provides thermodynamic and kinetic insights of Z. mobilis metabolism on dual substrates, and guidance of bioengineering efforts to increase hydrocarbon fuel production.

23-butanediol↗
DOE OSTI2023

A previously uncharacterized divisome-associated lipoprotein, DalA, is needed for normal cell division in Rhodobacterales

ABSTRACT The bacterial cell envelope is a key subcellular compartment with important roles in antibiotic resistance, nutrient acquisition, and cell morphology. We seek to gain a better understanding of proteins that contribute to the function of the cell envelope in Alphaproteobacteria . Using Rhodobacter sphaeroides , we show that a previously uncharacterized protein, RSP_1200, is an outer membrane (OM) lipoprotein that non-covalently binds peptidoglycan (PG). Using a fluorescently tagged version of this protein, we find that RSP_1200 undergoes a dynamic repositioning during the cell cycle and is enriched at the septum during cell division. We show that the position of RSP_1200 mirrors the location of FtsZ rings, leading us to propose that RSP_1200 is a newly identified component of the R. sphaeroides ’ divisome. Additional support for this hypothesis includes the co-precipitation of RSP_1200 with FtsZ, the Pal protein, and several predicted PG L,D-transpeptidases. We also find that a ∆ RSP_1200 mutation leads to defects in cell division, sensitivity to PG-active antibiotics, and results in the formation of OM protrusions at the septum during cell division. Based on these results, we propose to name RSP_1200 DalA (for division-associated lipoprotein A) and postulate that DalA serves as a scaffold to position or modulate the activity of PG transpeptidases that are needed to form envelope invaginations during cell division. We find that DalA homologs are present in members of the Rhodobacterales order within Alphaproteobacteria . Therefore, we propose that further analysis of this and related proteins will increase our understanding of the macromolecular machinery and proteins that participate in cell division in Gram-negative bacteria. IMPORTANCE Multi-protein complexes of the bacterial cell envelope orchestrate key processes like growth, division, biofilm formation, antimicrobial resistance, and production of valuable compounds. The subunits of these protein complexes are well studied in some bacteria, and differences in their composition and function are linked to variations in cell envelope composition, shape, and proliferation. However, some envelope protein complex subunits have no known homologs across the bacterial phylogeny. We find that Rhodobacter sphaeroides RSP_1200 is a newly identified lipoprotein (DalA) and that loss of this protein causes defects in cell division and changes the sensitivity to compounds, affecting cell envelope synthesis and function. We find that DalA forms a complex with proteins needed for cell division, binds the cell envelope polymer peptidoglycan, and colocalizes with enzymes involved in the assembly of this macromolecule. The analysis of DalA provides new information on the cell division machinery in this and possibly other Alphaproteobacteria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Deuteron off-resonance rotating frame relaxation for the characterization of slow motions in rotating and static solid-state proteins

Here, we demonstrate the feasibility of deuterium solid-state NMR off-resonance rotating frame relaxation measurements for studies of slow motions in biomolecular solids. The pulse sequence, which includes adiabatic pulses for magnetization alignment, is illustrated for static and magic-angle spinning conditions away from rotary resonances. We apply the measurements for three systems with selective deuterium labels at methyl groups: a) a model compound, Fluorenylmethyloxycarbonyl methionine-D 3 amino acid, for which the principles of the measurements and corresponding motional modeling based on rotameric interconversions are demonstrated; b) amyloid-β 1-40 fibrils labeled at a single alanine methyl group located in the disordered N-terminal domain. This system has been extensively studied in prior work and here serves as a test of the method for complex biological systems. The essential features of the dynamics consist of large-scale rearrangements of the disordered N-terminal domain and the conformational exchange between the free and bound forms of the domain, the latter one due to transient interactions with the structured core of the fibrils. and c) a 15-residue helical peptide which belongs to the predicted α-helical domain near the N-terminus of apolipoprotein B. The peptide is solvated with triolein and incorporates a selectively labeled leucine methyl groups. The method permits model refinement, indicating rotameric interconversions with a distribution of rate constants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Expression, purification, and characterization of diacylated Lipo-YcjN from Escherichia coli

YcjN is a putative substrate binding protein expressed from a cluster of genes involved in carbohydrate import and metabolism in Escherichia coli. Here, we determine the crystal structure of YcjN to a resolution of 1.95 Å, revealing that its three-dimensional structure is similar to substrate binding proteins in subcluster D-I, which includes the well-characterized maltose binding protein. Furthermore, we found that recombinant overexpression of YcjN results in the formation of a lipidated form of YcjN that is posttranslationally diacylated at cysteine 21. Comparisons of size-exclusion chromatography profiles and dynamic light scattering measurements of lipidated and nonlipidated YcjN proteins suggest that lipidated YcjN aggregates in solution via its lipid moiety. Additionally, bioinformatic analysis indicates that YcjN-like proteins may exist in both Bacteria and Archaea, potentially in both lipidated and nonlipidated forms. Together, our results provide a better understanding of the aggregation properties of recombinantly expressed bacterial lipoproteins in solution and establish a foundation for future studies that aim to elucidate the role of these proteins in bacterial physiology.

Escherichia coli↗
DOE OSTI2025

ER stress and viral defense: Advances and future perspectives on plant unfolded protein response in pathogenesis

Viral infections pose significant threats to crop productivity and agricultural sustainability. The frequency and severity of these infections are increasing, and pathogens are evolving rapidly under the influence of climate change. This underscores the importance of exploring the fundamental mechanisms by which plants defend themselves against dynamic viral threats. One such mechanism is the unfolded protein response (UPR), which is activated when the protein folding demand exceeds the capacity of the endoplasmic reticulum, particularly under adverse environmental conditions. While the key regulators of the UPR in response to viral infections have been identified, our understanding of how they modulate the UPR to suppress plant viral infections at the molecular and genetic levels is still in its infancy. Recent findings have shown that, in response to plant viral infections, the UPR swiftly reprograms transcriptional changes to support cellular, metabolic, and physiological processes associated with cell viability. However, the underlying mechanisms and functional outcomes of these changes remain largely unexplored. Here, we highlight recent advances in plant UPR research and summarize key findings related to viral infection–induced UPR, focusing on the balance between prosurvival and prodeath strategies. We also discuss the potential of systems-level approaches to uncover the full extent of the functional link between the UPR and plant responses to viral infections.

ER stress↗
NASA NTRS2003

Physiological Maturation of Regenerating Hair Cells

The bullfrog saccule, a sensor of gravity and substrate-borne vibration, is a model system for hair cell transduction. Saccular hair cells also increase in number throughout adult life and rapidly recover after hair cell damage, making this organ an ideal system for studying hair cell development, repair, and regeneration. We have used of hair cell and supporting cell immunocytochemical markers to identify damaged hair cells and hair cell precursors in organotypic cultures of the bullfrog saccule. We then used an innovative combination of confocal, electron, and time-lapse microscopy to study the fate of damaged hair cells and the origin of new hair cells after gentamicin ototoxicity in normal and mitotically blocked saccular cultures. These studies have shown that gentamicin ototoxicity produces both lethal and sublethal hair cell damage. They have also shown that hair cell recovery in this organ takes place by both the repair of sublethally damaged hair cells and by the replacement of lost hair cells by mitotic regeneration. In parallel studies, we have used biophysical and molecular biological techniques to study the differentiation and innervation of developing, repairing, and regenerating hair cells. More specifically, we have used RT-PCR to obtain the bullfrog homologues of L-type voltage- gated calcium (L-VGCC) and large-conductance Ca(2+)-activated potassium (BK) channel genes. We have then obtained probes for these genes and, using in situ hybridization, begun to examine their expression in the bullfrog saccule and amphibian papilla. We have also used fluorescent-labeled channel toxins and channel toxin derivatives to determine the time of appearance of L-type voltage-gated calcium (L-VGCC) and Ca(2+)-activated potassium (BK) channels and to study dynamic changes in the number, distribution, and co-localization of these proteins in developing, repairing, and regenerating hair cells. Using time-lapse microscopy, we are also studying the dynamic relationship between ion channel clustering and synaptic formation in hair cells and afferent neurons. In future studies, we will determine when hair cell precursors acquire electrical tuning, and, using whole-cell patch-clamp techniques, identify and characterize their L-VGCC and BK currents. We will also use biophysical techniques to determine the number of L-VGCC and BK channels and the size and gating kinetics of their underlying L-VGCC and BK conductances, correlating these variables with the amplitude and frequency of membrane oscillations produced by intracellular current steps. We expect these studies to determine how hair cells regulate ion channel expression to achieve specific physiological responses.

Baird, Richard A.↗
DOE OSTI2023

Engineering activatable promoters for scalable and multi-input CRISPRa/i circuits

Dynamic, multi-input gene regulatory networks (GRNs) are ubiquitous in nature. Multilayer CRISPR-based genetic circuits hold great promise for building GRNs akin to those found in naturally occurring biological systems. We develop an approach for creating high-performing activatable promoters that can be assembled into deep, wide, and multi-input CRISPR-activation and -interference (CRISPRa/i) GRNs. By integrating sequence-based design and in vivo screening, we engineer activatable promoters that achieve up to 1,000-fold dynamic range in an Escherichia coli-based cell-free system. These components enable CRISPRa GRNs that are six layers deep and four branches wide. We show the generalizability of the promoter engineering workflow by improving the dynamic range of the light-dependent EL222 optogenetic system from 6-fold to 34-fold. Additionally, high dynamic range promoters enable CRISPRa systems mediated by small molecules and protein–protein interactions. We apply these tools to build input-responsive CRISPRa/i GRNs, including feedback loops, logic gates, multilayer cascades, and dynamic pulse modulators. Our work provides a generalizable approach for the design of high dynamic range activatable promoters and enables classes of gene regulatory functions in cell-free systems.

42 ENGINEERING↗
DOE OSTI2023

Systematic mapping of the conformational landscape and dynamism of soluble fibrinogen

Background: Fibrinogen is a soluble, multisubunit, and multidomain dimeric protein, which, upon its proteolytic cleavage by thrombin, is converted to insoluble fibrin, initiating polymerization that substantially contributes to clot growth. Fibrinogen contains numerous, transiently accessible “cryptic” epitopes for hemostatic and immunologic proteins, suggesting that fibrinogen exhibits conformational flexibility, which may play functional roles in its temporal and spatial interactions. Hitherto, there have been limited integrative approaches characterizing the solution structure and internal flexibility of fibrinogen. Methods: Here, utilizing a multipronged, biophysical approach involving 2 solution-based techniques, temperature-dependent hydrogen-deuterium exchange mass spectrometry and small angle X-ray scattering, corroborated by negative stain electron microscopy, we present a holistic, conformationally dynamic model of human fibrinogen in solution. Results: Our data reveal 4 major and distinct conformations of fibrinogen accommodated by a high degree of internal protein flexibility along its central scaffold. We propose that the fibrinogen structure in the solution consists of a complex, conformational landscape with multiple local minima. This is further supported by the location of numerous point mutations that are linked to dysfibrinogenemia and posttranslational modifications, residing near the identified fibrinogen flexions. Conclusion: Further, this work provides a molecular basis for the structural “dynamism” of fibrinogen that is expected to influence the broad swath of its functionally diverse macromolecular interactions and fine-tune the structural and mechanical properties of blood clots.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2010

Interactions of Ku70/80 with Double-Strand DNA: Energetic, Dynamics, and Functional Implications

Space radiation is a proficient inducer of DNA damage leading to mutation, aberrant cell signaling, and cancer formation. Ku is among the first responding proteins in nucleus to recognize and bind the DNA double strand breaks (DSBs) whenever they are introduced. Once loaded Ku works as a scaffold to recruit other repair factors of non-homologous end joining and facilitates the following repair processes. The crystallographic study of the Ku70/80 heterodimer indicate the core structure of this protein shows virtually no conformational change after binding with DNA. To investigate the dynamical features as well as the energetic characteristics of Ku-DNA binding, we conduct multi-nanosecond molecular dynamics simulations of a modeled Ku70/80 structure and several complexes with two 24-bp DNA duplexes. Free energy calculations show significant energy differences between the complexes with Ku bound at DSBs and those with Ku associated at an internal site of a chromosome. The results also reveal detailed interactions between different nucleotides and the amino acids along the DNA-binding cradle of Ku, indicating subtle binding preference of Ku at specific DNA sequences. The covariance matrix analyses along the trajectories demonstrate the protein is stimulated to undergo correlated motions of different domains once bound to DNA ends. Additionally, principle component analyses identify these low frequency collective motions suitable for binding with and translocation along duplex DNA. It is proposed that the modification of dynamical properties of Ku upon binding with DSBs may provide a signal for the further recruitment of other repair factors such as DNA-PKcs, XLF, and XRCC4.

Hu, Shaowen↗
DOE OSTI2023

Biochemical Process Modeling and Simulation (BPMS)

The Biochemical Process Modeling and Simulation project aims to reduce the cost and time of research by applying theory, modeling, and simulation to the most relevant bottlenecks in the biochemical process. We use molecular modeling, quantum mechanics, metabolic modeling, fluid dynamics, and reaction-diffusion methods in close collaboration with pretreatment, hydrolysis, upgrading, and TEA. The project's outcomes are increased yields and efficiency of the biochemical process, added value to products, and reduced price of fuels by specifically targeting catalytic efficiency, reactor design, enzyme efficiency, and microbial design. We work closely with experimental projects to identify problems and iterate with experiments to find and refine solutions. By working with experimentalists, we decide on problems that can be solved with simulation that could otherwise not be solved or would take too long with experiment alone to reach BETO's targets. Over the years, we have produced solutions that have resulted in determining the most likely fatty-acid derivative for passive transport out of bacteria that upgrade biomass, and we have also designed enzyme mutations for enhanced lignin upgrading. Metabolic models have been developed to tune the activity of 2,3 butanediol production for the 2030 target. A computational method to deliver understanding of how complex omics data can be interpreted in the metabolic pathways of organisms used in the Agile Biofoundry. We have found methods to overcome specific barriers and continue to develop those methods. Our reactor studies have guided the design of both the microbes and reactors for aerobic and micro-aerobic production at all scales and have been instrumental in improving the accuracy of techno-economic analysis models. This project is essential in the process of selecting the final processes for 2030 SAF production targets. More specifically, recently, we have: 1) Predicted the strength of the basic structural interactions in commodity plastics to provide guidance for plastics upcycling strategies. 2) Developed computational tool to improve the characterization of lignin-derived compounds 3) Developed new methodologies to enable Machine Learning-based Directed Evolution for protein engineering. 4) Developed Machine Learning methods to predict protein promiscuity and mutations to further improve microbial and enzymatic driven processes and demonstrated the utility of ML approaches to engineering proteins from sparse experimental datasets. 5) Developed new methods to enable high-fidelity simulation of aerobic fermentation at industrial scale and resolving mismatch of time scales through subcycling/operator splitting 7) Identified the difficulty in preventing local high-oxygen conditions in industrial bubble columns, which leads to less-desirable acetoin production, suggesting future research directions in alternative reactor configurations (e.g loop reactors, shallow-channel reactors).

BIOMASS FUELS↗
DOE OSTI2024

Mapping structural and dynamic divergence across the MBOAT family

Membrane-bound O-acyltransferases (MBOATs) are membrane-embedded enzymes that catalyze acyl chain transfer to a diverse group of substrates, including lipids, small molecules, and proteins. MBOATs share a conserved structural core, despite wide-ranging functional specificity across both prokaryotes and eukaryotes. The structural basis of catalytic specificity, regulation and interactions with the surrounding environment remain uncertain. Here, we combine comparative molecular dynamics (MD) simulations with bioinformatics to assess molecular and interactional divergence across the family. In simulations, MBOATs differentially distort the bilayer depending on their substrate type. Additionally, we identify lipid binding sites surrounding reactant gates in the surrounding membrane. Complementary bioinformatic analyses reveal a conserved role for re-entrant loop-2 in MBOAT fold stabilization and a key hydrogen bond bridging DGAT1 dimerization. Finally, we predict differences in MBOAT solvation and water gating properties. These data are pertinent to the design of MBOAT-specific inhibitors that encompass dynamic information within cellular mimetic environments.

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

Packaging “vegetable oils”: Insights into plant lipid droplet proteins

Abstract Plant neutral lipids, also known as “vegetable oils”, are synthesized within the endoplasmic reticulum (ER) membrane and packaged into subcellular compartments called lipid droplets (LDs) for stable storage in the cytoplasm. The biogenesis, modulation, and degradation of cytoplasmic LDs in plant cells are orchestrated by a variety of proteins localized to the ER, LDs, and peroxisomes. Recent studies of these LD-related proteins have greatly advanced our understanding of LDs not only as steady oil depots in seeds but also as dynamic cell organelles involved in numerous physiological processes in different tissues and developmental stages of plants. In the past 2 decades, technology advances in proteomics, transcriptomics, genome sequencing, cellular imaging and protein structural modeling have markedly expanded the inventory of LD-related proteins, provided unprecedented structural and functional insights into the protein machinery modulating LDs in plant cells, and shed new light on the functions of LDs in nonseed plant tissues as well as in unicellular algae. Here, we review critical advances in revealing new LD proteins in various plant tissues, point out structural and mechanistic insights into key proteins in LD biogenesis and dynamic modulation, and discuss future perspectives on bridging our knowledge gaps in plant LD biology.

Cai, Yingqi (ORCID:0000000203575809)↗
DOE OSTI2023

RNA structures and dynamics with Å resolution revealed by x-ray free-electron lasers

RNA macromolecules, like proteins, fold to assume shapes that are intimately connected to their broadly recognized biological functions; however, because of their high charge and dynamic nature, RNA structures are far more challenging to determine. We introduce an approach that exploits the high brilliance of x-ray free-electron laser sources to reveal the formation and ready identification of angstrom-scale features in structured and unstructured RNAs. Previously unrecognized structural signatures of RNA secondary and tertiary structures are identified through wide-angle solution scattering experiments. With millisecond time resolution, we observe an RNA fold from a dynamically varying single strand through a base-paired intermediate to assume a triple-helix conformation. While the backbone orchestrates the folding, the final structure is locked in by base stacking. This method may help to rapidly characterize and identify structural elements in nucleic acids in both equilibrium and time-resolved experiments.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1997

A Rationale for System-Dependent Advantages and Disadvantages of Solution Crystal Growth at Low Gravity

Protein crystallization experiments at reduced gravity have yielded crystals that, depending on the specific material, are either superior or inferior in their structural perfection compared to counterparts grown at normal gravity. A reduction of the crystals' quality due to their growth at low gravity cannot be understood from existing models. Our experimental investigations of the ground-based crystallization of the protein lysozyme have revealed pronounced unsteady growth layer dynamics and associated defect formation under steady external conditions. Through scaling analysis and numerical simulations we show that the observed fluctuations originate from the coupling of bulk transport with non-linear interface kinetics under mixed kinetics-transport control of the growth rate. The amplitude of the fluctuations is smallest when either transport or interfacial kinetics dominate the control of the crystallization process. Thus, depending on the specific system, crystal quality may be improved by either enhancing or suppressing the transport in the solution. These considerations provide, for the first time, a material-dependent rationale for the advantages, as well as the disadvantages, of reduced gravity for (protein) crystallization.

Rosenberger, Franz↗
DOE OSTI2021

Toho-1 β-lactamase: backbone chemical shift assignments and changes in dynamics upon binding with avibactam

Backbone chemical shift assignments for the Toho-1 β-lactamase (263 amino acids, 28.9 kDa) are reported based on triple resonance solution-state NMR experiments performed on a uniformly 2 H, 13 C, 15 N-labeled sample. These assignments allow for subsequent site-specific characterization at the chemical, structural, and dynamical levels. At the chemical level, titration with the non-β-lactam β-lactamase inhibitor avibactam is found to give chemical shift perturbations indicative of tight covalent binding that allow for mapping of the inhibitor binding site. At the structural level, protein secondary structure is predicted based on the backbone chemical shifts and protein residue sequence using TALOS-N and found to agree well with structural characterization from X-ray crystallography. At the dynamical level, model-free analysis of 15 N relaxation data at a single field of 16.4 T reveals well-ordered structures for the ligand-free and avibactam-bound enzymes with generalized order parameters of ~0.85. Complementary relaxation dispersion experiments indicate that there is an escalation in motions on the millisecond timescale in the vicinity of the active site upon substrate binding. The combination of high rigidity on short timescales and active site flexibility on longer timescales is consistent with hypotheses for achieving both high catalytic efficiency and broad substrate specificity: the induced active site dynamics allows variously sized substrates to be accommodated and increases the probability that the optimal conformation for catalysis will be sampled.

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

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

3DRISM↗