Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mutation”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Conformational heterogeneity in the dGsw purine riboswitch: role of Mg²⁺ and 2’-dG in aptamer folding

Recent advancements in RNA structural biology have focused on unraveling the complexities of non-coding mRNA elements like riboswitches. These cis-acting regulatory regions undergo structural changes in response to specific cellular metabolites, leading to up or downregulation of downstream genes. The purine riboswitch family regulates many prokaryotic genes involved in purine degradation and biosynthesis. They feature an aptamer domain organized around a 3-way helical junction, where ligand encapsulation occurs at the junctional core. In our study, we chemically probed the aptamer domain of the 2’-dG-sensing purine riboswitch from Mesoplasma florum (dGsw) under various solution conditions to understand how Mg²⁺ and 2’-dG influence riboswitch folding. Here, we find that efficient 2’-dG binding strongly depends on Mg²⁺, indicating that Mg²⁺ is essential for priming dGsw for ligand interactions. We identified a previously undescribed sequence in the 5’ tail of dGsw that is complementary to a conserved helix. The inclusion of this region in a construct led to intramolecular competition between the alternate helix, Palt, and P1. Mutational analysis confirmed that 5’ flanking end of the aptamer domain forms an alternate helix in the absence of ligand. Molecular dynamics simulations revealed that this alternative conformation is stable. This helix may, therefore, facilitate the formation of an anti-terminator helix by opening the 3-way junction surrounding the 2’-dG binding site. Our study further establishes the importance of a closed terminal P1 helix conformation for metabolite binding and suggests that the delicate interplay between P1 and Palt may fine-tune downstream gene regulation. These insights offer a new perspective on riboswitch structure and enhance our understanding of the role that a conformational ensemble plays in riboswitch activity and regulation.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Data for "Enhancing Lipid Production in Plant Cells through Automated High-Throughput Genome Engineering and Phenotyping"

Plant bioengineering is a time-consuming and labor-intensive process with no guarantee of achieving desired traits. Here, we present a fast, automated, scalable, high-throughput pipeline for plant bioengineering (FAST-PB) in maize (Zea mays) and Nicotiana benthamiana. FAST-PB enables genome editing and product characterization by integrating automated biofoundry engineering of callus and protoplast cells with single-cell matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS). We first demonstrated that FAST-PB could streamline Golden Gate cloning, with the capacity to construct 96 vectors in parallel. Using FAST-PB in protoplasts, we found that PEG2050 increased transfection efficiency by over 45%. For proof-of-concept, we established a reporter-gene-free method for CRISPR editing and phenotyping via mutation of high chlorophyll fluorescence 136. We show that diverse lipids were enhanced up to 6-fold using CRISPR activation of lipid controlling genes. In callus cells, an automated transformation platform was employed to regenerate plants with enhanced lipid traits through introducing multigene cassettes. Lastly, FAST-PB enabled high-throughput single-cell lipid profiling by integrating MALDI-MS with the biofoundry, protoplast, and callus cells, differentiating engineered and unengineered cells using single-cell lipidomics. These innovations massively increase the throughput of synthetic biology, genome editing, and metabolic engineering and change what is possible using single-cell metabolomics in plants.

AI/ML↗
DOE OSTI2020

Data from A Mass Spectrometry-Based High-Throughput Screening Method for Engineering Fatty Acid Synthases with Improved Production of Medium Chain Fatty Acids

Microbial cell factories have been extensively engineered to produce free fatty acids (FFAs) as key components of crucial nutrients, soaps, industrial chemicals, and fuels. However, our ability to control the composition of microbially synthesized FFAs is still limited, particularly, for producing medium‐chain fatty acids (MCFAs). This is mainly due to the lack of high‐throughput approaches for FFA analysis to engineer enzymes with desirable product specificity. Here we report a mass spectrometry (MS)‐based method for rapid profiling of MCFAs in Saccharomyces cerevisiae by using membrane lipids as a proxy. In particular, matrix‐assisted laser desorption/ionization time‐of‐flight (MALDI‐ToF) MS was used to detect shorter acyl chain phosphatidylcholines from membrane lipids and a higher m/z peak ratio at 730 and 758 was used as an indication for improved MCFA production. This colony‐based method can be performed at a rate of ~2 s per sample, representing a substantial improvement over gas chromatography‐MS (typically >30 min per sample) as the gold standard method for FFA detection. To demonstrate the power of this method, we performed site‐saturation mutagenesis of the yeast fatty acid synthase and identified nine missense mutations that resulted in improved MCFA production relative to the wild‐type strain. Colony‐based MALDI‐ToF MS screening provides an effective approach for engineering microbial fatty acid compositions in a high‐throughput manner.

Conversion↗
DOE OSTI2024

Data for Rapid and Efficient in planta Genome Editing in Sorghum Using Foxtail Mosaic Virus-mediated sgRNA Delivery

The requirement of in vitro tissue culture for the delivery of gene editing reagents limits the application of gene editing to commercially relevant varieties of many crop species. To overcome this bottleneck, plant RNA viruses have been deployed as versatile tools for in planta delivery of recombinant RNA. Viral delivery of single-guide RNAs (sgRNAs) to transgenic plants that stably express CRISPR-associated (Cas) endonuclease has been successfully used for targeted mutagenesis in several dicotyledonous and few monocotyledonous plants. Progress with this approach in monocotyledonous plants is limited so far by the availability of effective viral vectors. We engineered a set of foxtail mosaic virus (FoMV) and barley stripe mosaic virus (BSMV) vectors to deliver the fluorescent protein AmCyan to track viral infection and movement in Sorghum bicolor . We further used these viruses to deliver and express sgRNAs to Cas9 and Green Fluorescent Protein (GFP) expressing transgenic sorghum lines, targeting Phytoene desaturase (PDS), Magnesium-chelatase subunit I (MgCh), 4-hydroxy-3-methylbut-2-enyl diphosphate reductase, orthologs of maize Lemon white1 (Lw1) or GFP. The recombinant BSMV did neither infect sorghum nor deliver or express AmCyan and sgRNAs. In contrast, the recombinant FoMV systemically spread throughout sorghum plants and induced somatic mutations with frequencies reaching up to 60%. This mutagenesis led to visible phenotypic changes, demonstrating the potential of FoMV for in planta gene editing and functional genomics studies in sorghum.

Feedstock Production↗
DOE OSTI2020

Structural analyses of the PKA RIIβ holoenzyme containing the oncogenic DnaJB1-PKAc fusion protein reveal protomer asymmetry and fusion-induced allosteric perturbations in fibrolamellar hepatocellular carcinoma

When the J-domain of the heat shock protein DnaJB1 is fused to the catalytic (C) subunit of cAMP-dependent protein kinase (PKA), replacing exon 1, this fusion protein, J-C subunit (J-C), becomes the driver of fibrolamellar hepatocellular carcinoma (FL-HCC). Here, we use cryo-electron microscopy (cryo-EM) to characterize J-C bound to RIIβ, the major PKA regulatory (R) subunit in liver, thus reporting the first cryo-EM structure of any PKA holoenzyme. We report several differences in both structure and dynamics that could not be captured by the conventional crystallography approaches used to obtain prior structures. Most striking is the asymmetry caused by the absence of the second cyclic nucleotide binding (CNB) domain and the J-domain in one of the RIIβ:J-C protomers. Using molecular dynamics (MD) simulations, we discovered that this asymmetry is already present in the wild-type (WT) RIIβ 2 C 2 but had been masked in the previous crystal structure. This asymmetry may link to the intrinsic allosteric regulation of all PKA holoenzymes and could also explain why most disease mutations in PKA regulatory subunits are dominant negative. The cryo-EM structure, combined with small-angle X-ray scattering (SAXS), also allowed us to predict the general position of the Dimerization/Docking (D/D) domain, which is essential for localization and interacting with membrane-anchored A-Kinase-Anchoring Proteins (AKAPs). This position provides a multivalent mechanism for interaction of the RIIβ holoenzyme with membranes and would be perturbed in the oncogenic fusion protein. The J-domain also alters several biochemical properties of the RIIβ holoenzyme: It is easier to activate with cAMP, and the cooperativity is reduced. These results provide new insights into how the finely tuned allosteric PKA signaling network is disrupted by the oncogenic J-C subunit, ultimately leading to the development of FL-HCC.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The novel anti-CRISPR AcrIIA22 relieves DNA torsion in target plasmids and impairs SpyCas9 activity

To overcome CRISPR-Cas defense systems, many phages and mobile genetic elements (MGEs) encode CRISPR-Cas inhibitors called anti-CRISPRs (Acrs). Nearly all characterized Acrs directly bind Cas proteins to inactivate CRISPR immunity. Here, using functional metagenomic selection, we describe AcrIIA22, an unconventional Acr found in hypervariable genomic regions of clostridial bacteria and their prophages from human gut microbiomes. AcrIIA22 does not bind strongly to SpyCas9 but nonetheless potently inhibits its activity against plasmids. To gain insight into its mechanism, we obtained an X-ray crystal structure of AcrIIA22, which revealed homology to PC4-like nucleic acid–binding proteins. Based on mutational analyses and functional assays, we deduced that acrIIA22 encodes a DNA nickase that relieves torsional stress in supercoiled plasmids. This may render them less susceptible to SpyCas9, which uses free energy from negative supercoils to form stable R-loops. Modifying DNA topology may provide an additional route to CRISPR-Cas resistance in phages and MGEs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A new class of antibodies that overcomes a steric barrier to cross-group neutralization of influenza viruses

Antibody titers that inhibit the influenza virus hemagglutinin (HA) from engaging its receptor are the accepted correlate of protection from infection. Many potent antibodies with broad, intra-subtype specificity bind HA at the receptor binding site (RBS). One barrier to broad H1-H3 cross-subtype neutralization is an insertion (133a) between positions 133 and 134 on the rim of the H1 HA RBS. We describe here a class of antibodies that overcomes this barrier. These genetically unrestricted antibodies are abundant in the human B cell memory compartment. Analysis of the affinities of selected members of this class for historical H1 and H3 isolates suggest that they were elicited by H3 exposure and broadened or diverted by later exposure(s) to H1 HA. RBS mutations in egg-adapted vaccine strains cause the new H1 specificity of these antibodies to depend on the egg adaptation. The results suggest that suitable immunogens might elicit 133a-independent, H1-H3 cross neutralization by RBS-directed antibodies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Early prediction of antigenic transitions for influenza A/H3N2

Influenza A/H3N2 is a rapidly evolving virus which experiences major antigenic transitions every two to eight years. Anticipating the timing and outcome of transitions is critical to developing effective seasonal influenza vaccines. Using a published phylodynamic model of influenza transmission, we identified indicators of future evolutionary success for an emerging antigenic cluster and quantified fundamental trade-offs in our ability to make such predictions. The eventual fate of a new cluster depends on its initial epidemiological growth rate––which is a function of mutational load and population susceptibility to the cluster––along with the variance in growth rate across co-circulating viruses. Logistic regression can predict whether a cluster at 5% relative frequency will eventually succeed with ~80% sensitivity, providing up to eight months advance warning. As a cluster expands, the predictions improve while the lead-time for vaccine development and other interventions decreases. However, attempts to make comparable predictions from 12 years of empirical influenza surveillance data, which are far sparser and more coarse-grained, achieve only 56% sensitivity. By expanding influenza surveillance to obtain more granular estimates of the frequencies of and population-wide susceptibility to emerging viruses, we can better anticipate major antigenic transitions. This provides added incentives for accelerating the vaccine production cycle to reduce the lead time required for strain selection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Dose optimization of an adjuvanted peptide-based personalized neoantigen melanoma vaccine

The advancements in next-generation sequencing have made it possible to effectively detect somatic mutations, which has led to the development of personalized neoantigen cancer vaccines that are tailored to the unique variants found in a patient’s cancer. These vaccines can provide significant clinical benefit by leveraging the patient’s immune response to eliminate malignant cells. However, determining the optimal vaccine dose for each patient is a challenge due to the heterogeneity of tumors. To address this challenge, we formulate a mathematical dose optimization problem based on a previous mathematical model that encompasses the immune response cascade produced by the vaccine in a patient. We propose an optimization approach to identify the optimal personalized vaccine doses, considering a fixed vaccination schedule, while simultaneously minimizing the overall number of tumor and activated T cells. To validate our approach, we perform in silico experiments on six real-world clinical trial patients with advanced melanoma. We compare the results of applying an optimal vaccine dose to those of a suboptimal dose (the dose used in the clinical trial and its deviations). Our simulations reveal that an optimal vaccine regimen of higher initial doses and lower final doses may lead to a reduction in tumor size for certain patients. Our mathematical dose optimization offers a promising approach to determining an optimal vaccine dose for each patient and improving clinical outcomes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Chromoanagenesis from radiation-induced genome damage in Populus

Chromoanagenesis is a genomic catastrophe that results in chromosomal shattering and reassembly. These extreme single chromosome events were first identified in cancer, and have since been observed in other systems, but have so far only been formally documented in plants in the context of haploid induction crosses. The frequency, origins, consequences, and evolutionary impact of such major chromosomal remodeling in other situations remain obscure. Here, we demonstrate the occurrence of chromoanagenesis in poplar ( Populus sp . ) trees produced from gamma-irradiated pollen. Specifically, in this population of siblings carrying indel mutations, two individuals exhibited highly frequent copy number variation (CNV) clustered on a single chromosome, one of the hallmarks of chromoanagenesis. Using short-read sequencing, we confirmed the presence of clustered segmental rearrangement. Independently, we identified and validated novel DNA junctions and confirmed that they were clustered and corresponded to these rearrangements. Our reconstruction of the novel sequences suggests that the chromosomal segments have reorganized randomly to produce a novel rearranged chromosome but that two different mechanisms might be at play. Our results indicate that gamma irradiation can trigger chromoanagenesis, suggesting that this may also occur when natural or induced mutagens cause DNA breaks. We further demonstrate that such events can be tolerated in poplar, and even replicated clonally, providing an attractive system for more in-depth investigations of their consequences.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

PHYTOCHROME C regulation of photoperiodic flowering via PHOTOPERIOD1 is mediated by EARLY FLOWERING 3 in Brachypodium distachyon

Daylength sensing in many plants is critical for coordinating the timing of flowering with the appropriate season. Temperate climate-adapted grasses such as Brachypodium distachyon flower during the spring when days are becoming longer. The photoreceptor PHYTOCHROME C is essential for long-day (LD) flowering in B. distachyon. PHYC is required for the LD activation of a suite of genes in the photoperiod pathway including PHOTOPERIOD1 (PPD1) that, in turn, result in the activation of FLOWERING LOCUS T (FT1)/FLORIGEN, which causes flowering. Thus, B. distachyon phyC mutants are extremely delayed in flowering. Here we show that PHYC-mediated activation of PPD1 occurs via EARLY FLOWERING 3 (ELF3), a component of the evening complex in the circadian clock. The extreme delay of flowering of the phyC mutant disappears when combined with an elf3 loss-of-function mutation. Moreover, the dampened PPD1 expression in phyC mutant plants is elevated in phyC/elf3 mutant plants consistent with the rapid flowering of the double mutant. We show that loss of PPD1 function also results in reduced FT1 expression and extremely delayed flowering consistent with results from wheat and barley. Additionally, elf3 mutant plants have elevated expression levels of PPD1, and we show that overexpression of ELF3 results in delayed flowering associated with a reduction of PPD1 and FT1 expression, indicating that ELF3 represses PPD1 transcription consistent with previous studies showing that ELF3 binds to the PPD1 promoter. Indeed, PPD1 is the main target of ELF3-mediated flowering as elf3/ppd1 double mutant plants are delayed flowering. Our results indicate that ELF3 operates downstream from PHYC and acts as a repressor of PPD1 in the photoperiod flowering pathway of B. distachyon.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Machine learning identifies novel signatures of antifungal drug resistance in Saccharomycotina yeasts

Antifungal drug resistance is a major challenge in fungal infection management. Numerous genomic changes are known to contribute to acquired drug resistance in clinical isolates of specific pathogens, but whether they broadly explain natural resistance across entire lineages is unknown. We leveraged genomic, ecological, and phenotypic trait data from naturally sampled strains from nearly all known species in subphylum Saccharomycotina to examine the evolution of resistance to eight antifungal drugs. The phylogenetic distribution of drug resistance varied by drug; fluconazole resistance was widespread, while 5-fluorocytosine resistance was rare, except in Lipomycetales. A random forest algorithm trained on genomic data predicted drug-resistant yeasts with 54–75% accuracy. Fluconazole resistance was consistently predicted with the highest accuracy (75.2%). Furthermore, fluconazole resistance prediction accuracy was similar between models trained on genome-wide variation in the presence and number of InterPro protein annotations across Saccharomycotina (75.2%) and those trained on amino acid sequence alignment data of Erg11, a protein known to be involved in fluconazole resistance (74.3-74.9%). Interestingly, the top Erg11 residues for predicting fluconazole resistance across Saccharomycotina do not overlap with, are not spatially close to, and are less conserved than those previously linked to resistance in clinical isolates of Candida albicans. In silico deep mutational scanning of the C. albicans Erg11 protein reveals that amino acid variants implicated in clinical cases of resistance are almost universally destabilizing while variants in our most informative residues are energetically more neutral, explaining why the latter are much more common than the former in natural populations. Importantly, previous experimental analyses of C. albicans Erg11 have shown that amino acid variation in our most informative residues, despite having never been directly implicated in clinical cases, can directly contribute to resistance. Our results suggest that studies of natural resistance in yeast species never encountered in the clinic will yield a fuller understanding of antifungal drug resistance.

Harrison, Marie-Claire [Vanderbilt Univ., Nashvill↗
DOE OSTI2020

The effect of enrofloxacin on enteric Escherichia coli : Fitting a mathematical model to in vivo data

Antimicrobial drugs administered systemically may cause the emergence and dissemination of antimicrobial resistance among enteric bacteria. To develop logical, research-based recommendations for food animal veterinarians, we must understand how to maximize antimicrobial drug efficacy while minimizing risk of antimicrobial resistance. Our objective is to evaluate the effect of two approved dosing regimens of enrofloxacin (a single high dose or three low doses) on Escherichia coli in cattle. We look specifically at bacteria above and below the epidemiological cutoff (ECOFF), above which the bacteria are likely to have an acquired or mutational resistance to enrofloxacin. We developed a differential equation model for the antimicrobial drug concentrations in plasma and colon, and bacteria populations in the feces. The model was fit to animal data of drug concentrations in the plasma and colon obtained using ultrafiltration probes. Fecal E. coli counts and minimum inhibitory concentrations were measured for the week after receiving the antimicrobial drug. We predict that the antimicrobial susceptibility of the bacteria above the ECOFF pre-treatment strongly affects the composition of the bacteria following treatment. Faster removal of the antimicrobial drugs from the colon throughout the study leads to improved clearance of bacteria above the ECOFF in the low dose regimen. If we assume a fitness cost is associated with bacteria above the ECOFF, the increased fitness costs leads to reduction of bacteria above the ECOFF in the low dose study. These results suggest the initial E. coli susceptibility is a strong indicator of how steers respond to antimicrobial drug treatment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Predicting the viability of beta-lactamase: How folding and binding free energies correlate with beta-lactamase fitness

One of the long-standing holy grails of molecular evolution has been the ability to predict an organism’s fitness directly from its genotype. With such predictive abilities in hand, researchers would be able to more accurately forecast how organisms will evolve and how proteins with novel functions could be engineered, leading to revolutionary advances in medicine and biotechnology. In this work, we assemble the largest reported set of experimental TEM-1 β-lactamase folding free energies and use this data in conjunction with previously acquired fitness data and computational free energy predictions to determine how much of the fitness of β-lactamase can be directly predicted by thermodynamic folding and binding free energies. We focus upon β-lactamase because of its long history as a model enzyme and its central role in antibiotic resistance. Based upon a set of 21 β-lactamase single and double mutants expressly designed to influence protein folding, we first demonstrate that modeling software designed to compute folding free energies such as FoldX and PyRosetta can meaningfully, although not perfectly, predict the experimental folding free energies of single mutants. Interestingly, while these techniques also yield sensible double mutant free energies, we show that they do so for the wrong physical reasons. We then go on to assess how well both experimental and computational folding free energies explain single mutant fitness. We find that folding free energies account for, at most, 24% of the variance in β-lactamase fitness values according to linear models and, somewhat surprisingly, complementing folding free energies with computationally-predicted binding free energies of residues near the active site only increases the folding-only figure by a few percent. This strongly suggests that the majority of β-lactamase’s fitness is controlled by factors other than free energies. Overall, our results shed a bright light on to what extent the community is justified in using thermodynamic measures to infer protein fitness as well as how applicable modern computational techniques for predicting free energies will be to the large data sets of multiply-mutated proteins forthcoming.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Identifying key determinants and dynamics of SARS-CoV-2/ACE2 tight interaction

SARS-CoV-2 virus, the causative agent of Covid-19, has fired up a global pandemic. The virus interacts with the human receptor angiotensin-converting enzyme 2 (ACE2) for an invasion via receptor binding domain (RBD) on its spike protein. To provide a deeper understanding of this interaction, we performed microsecond simulations of the RBD-ACE2 complex for SARS-CoV-2 and compared it with the closely related SARS-CoV discovered in 2003. We show residues in the RBD of SARS-CoV-2 that were mutated from SARS-CoV, collectively help make the RBD anchor much stronger to the N-terminal part of ACE2 than the corresponding residues on RBD of SARS-CoV. This would result in a reduced dissociation rate of SARS-CoV-2 from human receptor protein compared to SARS-CoV. The phenomenon was consistently observed in simulations beyond 500 ns and was reproducible across different force fields. Altogether, our study adds more insight into the critical dynamics of the key residues at the virus spike and human receptor binding interface and potentially aids the development of diagnostics and therapeutics to combat the pandemic efficiently.

60 APPLIED LIFE SCIENCES↗