Ancestral sequence reconstruction dissects structural and functional differences among eosinophil ribonucleases.

Evolutionarily conserved structural folds can give rise to diverse biological functions, yet predicting atomic-scale interactions that contribute to the emergence of novel activities within such folds remains challenging. Pancreatic-type ribonucleases illustrate this complexity, sharing a core structure that has evolved to accommodate varied functions. In this study, we used ancestral sequence reconstruction to probe evolutionary and molecular determinants that distinguish biological activities within eosinophil members of the RNase 2/3 subfamily. Our investigation unveils functional, structural, and dynamical behaviors that differentiate the evolved ancestral ribonuclease (AncRNase) from its contemporary eosinophil RNase orthologs. Leveraging the potential of ancestral reconstruction for protein engineering, we used AncRNase predictions to design a minimal 4-residue variant that transforms human RNase 2 into a chimeric enzyme endowed with the antimicrobial and cytotoxic activities of RNase 3 members. This work provides unique insights into mutational and evolutionary pathways governing structure, function, and conformational states within the eosinophil RNase subfamily, offering potential for targeted modulation of RNase-associated functions.

Emerging variants of SARS-CoV-2 NSP10 highlight strong functional conservation of its binding to two non-structural proteins, NSP14 and NSP16.

The coronavirus SARS-CoV-2 protects its RNA from being recognized by host immune responses by methylation of its 5' end, also known as capping. This process is carried out by two enzymes, non-structural protein 16 (NSP16) containing 2'-O-methyltransferase and NSP14 through its N7 methyltransferase activity, which are essential for the replication of the viral genome as well as evading the host's innate immunity. NSP10 acts as a crucial cofactor and stimulator of NSP14 and NSP16. To further understand the role of NSP10, we carried out a comprehensive analysis of >13 million globally collected whole-genome sequences (WGS) of SARS-CoV-2 obtained from the Global Initiative Sharing All Influenza Data (GISAID) and compared it with the reference genome Wuhan/WIV04/2019 to identify all currently known variants in NSP10. T12I, T102I, and A104V in NSP10 have been identified as the three most frequent variants and characterized using X-ray crystallography, biophysical assays, and enhanced sampling simulations. In contrast to other proteins such as spike and NSP6, NSP10 is significantly less prone to mutation due to its crucial role in replication. The functional effects of the variants were examined for their impact on the binding affinity and stability of both NSP14-NSP10 and NSP16-NSP10 complexes. These results highlight the limited changes induced by variant evolution in NSP10 and reflect on the critical roles NSP10 plays during the SARS-CoV-2 life cycle. These results also indicate that there is limited capacity for the virus to overcome inhibitors targeting NSP10 via the generation of variants in inhibitor binding pockets.

Altered Nucleotide Insertion Mechanisms of Disease-Associated TERT Variants.

Telomere biology disorders (TBDs) are a spectrum of diseases that arise from mutations in genes responsible for maintaining telomere integrity. Human telomerase reverse transcriptase (hTERT) adds nucleotides to chromosome ends and is frequently mutated in individuals with TBDs. Previous studies have provided insight into how relative changes in hTERT activity can lead to pathological outcomes. However, the underlying mechanisms describing how disease-associated variants alter the physicochemical steps of nucleotide insertion remain poorly understood. To address this, we applied single-turnover kinetics and computer simulations to the Tribolium castaneum TERT (tcTERT) model system and characterized the nucleotide insertion mechanisms of six disease-associated variants. Each variant had distinct consequences on tcTERT's nucleotide insertion mechanism, including changes in nucleotide binding affinity, rates of catalysis, or ribonucleotide selectivity. Our computer simulations provide insight into how each variant disrupts active site organization, such as suboptimal positioning of active site residues, destabilization of the DNA 3' terminus, or changes in nucleotide sugar pucker. Collectively, this work provides a holistic characterization of the nucleotide insertion mechanisms for multiple disease-associated TERT variants and identifies additional functions of key active site residues during nucleotide insertion.

Future COVID19 surges prediction based on SARS-CoV-2 mutations surveillance.

COVID19 has aptly revealed that airborne viruses such as SARS-CoV-2 with the ability to rapidly mutate combined with high rates of transmission and fatality can cause a deadly worldwide pandemic in a matter of weeks (Plato et al., 2021). Apart from vaccines and post-infection treatment options, strategies for preparedness will be vital in responding to the current and future pandemics. Therefore, there is wide interest in approaches that allow predictions of increase in infections ('surges') before they occur. We describe here real-time genomic surveillance particularly based on mutation analysis, of viral proteins as a methodology for a priori determination of surge in number of infection cases. The full results are available for SARS-CoV-2 at http://pandemics.okstate.edu/covid19/, and are updated daily as new virus sequences become available. This approach is generic and will also be applicable to other pathogens.

Conformational exchange divergence along the evolutionary pathway of eosinophil-associated ribonucleases.

The evolutionary role of conformational exchange in the emergence and preservation of function within structural homologs remains elusive. While protein engineering has revealed the importance of flexibility in function, productive modulation of atomic-scale dynamics has only been achieved on a finite number of distinct folds. Allosteric control of unique members within dynamically diverse structural families requires a better appreciation of exchange phenomena. Here, we examined the functional and structural role of conformational exchange within eosinophil-associated ribonucleases. Biological and catalytic activity of various EARs was performed in parallel to mapping their conformational behavior on multiple timescales using NMR and computational analyses. Despite functional conservation and conformational seclusion to a specific domain, we show that EARs can display similar or distinct motional profiles, implying divergence rather than conservation of flexibility. Comparing progressively more distant enzymes should unravel how this subfamily has evolved new functions and/or altered their behavior at the molecular level.

Capture of activated dioxygen intermediates at the copper-active site of a lytic polysaccharide monooxygenase.

Metalloproteins perform a diverse array of redox-related reactions facilitated by the increased chemical functionality afforded by their metallocofactors. Lytic polysaccharide monooxygenases (LPMOs) are a class of copper-dependent enzymes that are responsible for the breakdown of recalcitrant polysaccharides via oxidative cleavage at the glycosidic bond. The activated copper-oxygen intermediates and their mechanism of formation remains to be established. Neutron protein crystallography which permits direct visualization of protonation states was used to investigate the initial steps of oxygen activation directly following active site copper reduction in Neurospora crassa LPMO9D. Herein, we cryo-trap an activated dioxygen intermediate in a mixture of superoxo and hydroperoxo states, and we identify the conserved second coordination shell residue His157 as the proton donor. Density functional theory calculations indicate that both superoxo and hydroperoxo active site states are stable. The hydroperoxo formed is potentially an early LPMO catalytic reaction intermediate or the first step in the mechanism of hydrogen peroxide formation in the absence of substrate. We observe that the N-terminal amino group of the copper coordinating His1 remains doubly protonated directly following molecular oxygen reduction by copper. Aided by molecular dynamics and mining minima free energy calculations we establish that the conserved second-shell His161 in MtPMO3* displays conformational flexibility in solution and that this flexibility is also observed, though to a lesser extent, in His157 of NcLPMO9D. The imidazolate form of His157 observed in our structure following oxygen intermediate protonation can be attributed to abolished His157 flexibility due steric hindrance in the crystal as well as the solvent-occluded active site environment due to crystal packing. A neutron crystal structure of NcLPMO9D at low pH further supports occlusion of the active site since His157 remains singly protonated even at acidic conditions.

Mechanism of nucleotide discrimination by the translesion synthesis polymerase Rev1.

Rev1 is a translesion DNA synthesis (TLS) polymerase involved in the bypass of adducted-guanine bases and abasic sites during DNA replication. During damage bypass, Rev1 utilizes a protein-template mechanism of DNA synthesis, where the templating DNA base is evicted from the Rev1 active site and replaced by an arginine side chain that preferentially binds incoming dCTP. Here, we utilize X-ray crystallography and molecular dynamics simulations to obtain structural insight into the dCTP specificity of Rev1. We show the Rev1 R324 protein-template forms sub-optimal hydrogen bonds with incoming dTTP, dGTP, and dATP that prevents Rev1 from adopting a catalytically competent conformation. Additionally, we show the Rev1 R324 protein-template forms optimal hydrogen bonds with incoming rCTP. However, the incoming rCTP adopts an altered sugar pucker, which prevents the formation of a catalytically competent Rev1 active site. This work provides novel insight into the mechanisms for nucleotide discrimination by the TLS polymerase Rev1.

Altered APE1 activity on abasic ribonucleotides is mediated by changes in the nucleoside sugar pucker.

Ribonucleotides (rNTPs) are predicted to be incorporated into the genome at a rate of up to 3 million times per cell division, making rNTPs the most common non-standard nucleotide in the human genome. Typically, misinserted ribonucleotides are repaired by the ribonucleotide excision repair (RER) pathway, which is initiated by RNase H2 cleavage. However, rNTPs are susceptible to spontaneous depurination generating abasic ribonucleotides (rAPs), which are unable to be processed by RNase H2. Additionally, rAPs have been found in nascent RNA and coupled to R-loops. Recent work identified that base excision repair (BER) protein AP-Endonuclease 1 (APE1) is responsible for the initial processing of rAPs embedded in DNA and in R-loops. APE1 is a well characterized AP endonuclease that cleaves 5' of abasic sites, but its ability to cleave at rAPs remains poorly understood. Here, we utilize enzyme kinetics, X-ray crystallography, and molecular dynamics simulations to provide insight into rAP processing by APE1. Enzyme kinetics were used to determine pre-steady-state rates of APE1 cleavage on DNA substrates containing rAP, revealing a decrease in activity compared to cleavage at a canonical deoxy-AP substrate. Using X-ray crystallography, we identified novel contacts between the rAP and the APE1 active site. We demonstrate that the rAP sugar pucker is accommodated in the active site in a C3'-endo conformation, influencing its position and contributing to a decrease in activity compared to the deoxy-AP site. Together, this work provides molecular level insights into rAP processing by APE1 and advances our understanding of ribonucleotide processing within genomic DNA.

Allosteric communication in class A ß-lactamases occurs via cooperative coupling of loop dynamics.

Understanding allostery in enzymes and tools to identify it offer promising alternative strategies to inhibitor development. Through a combination of equilibrium and nonequilibrium molecular dynamics simulations, we identify allosteric effects and communication pathways in two prototypical class A ß-lactamases, TEM-1 and KPC-2, which are important determinants of antibiotic resistance. The nonequilibrium simulations reveal pathways of communication operating over distances of 30 Å or more. Propagation of the signal occurs through cooperative coupling of loop dynamics. Notably, 50% or more of clinically relevant amino acid substitutions map onto the identified signal transduction pathways. This suggests that clinically important variation may affect, or be driven by, differences in allosteric behavior, providing a mechanism by which amino acid substitutions may affect the relationship between spectrum of activity, catalytic turnover, and potential allosteric behavior in this clinically important enzyme family. Simulations of the type presented here will help in identifying and analyzing such differences.

Antibiotics are crucial drugs for treating and preventing bacterial infections, but some bacteria are evolving ways to resist their effects. This 'antibiotic resistance' threatens lives and livelihoods worldwide. ß-lactam antibiotics, like penicillin, are some of the most commonly used, but some bacteria can now make enzymes called ß-lactamases, which destroy these antibiotics. Dozens of different types of ß-lactamases now exist, each with different properties. Two of the most medically important are TEM-1 and KPC-2. One way to counteract ß-lactamases is with drugs called inhibitors that stop the activity of these enzymes. The approved ß-lactamase inhibitors work by blocking the part of the enzyme that binds and destroys antibiotics, known as the 'active site'. The ß-lactamases have evolved, some of which have the ability to resist the effects of known inhibitors. It is possible that targeting parts of ß-lactamases far from the active site, known as 'allosteric sites', might get around these new bacterial defences. A molecule that binds to an allosteric site might alter the enzyme's shape, or restrict its movement, making it unable to do its job. Galdadas, Qu et al. used simulations to understand how molecules binding at allosteric sites affect enzyme movement. The experiments examined the structures of both TEM-1 and KPC-2, looking at how their shapes changed as molecules were removed from the allosteric site. This revealed how the allosteric sites and the active site are linked together. When molecules were taken out of the allosteric sites, they triggered ripples of shape change that travelled via loop-like structures across the surface of the enzyme. These loops contain over half of the known differences between the different types of ß-lactamases, suggesting mutations here may be responsible for changing which antibiotics each enzyme can destroy. In other words, changes in the 'ripples' may be related to the ability of the enzymes to resist particular antibiotics. Understanding how changes in one part of a ß-lactamase enzyme reach the active site could help in the design of new inhibitors. It might also help to explain how ß-lactamases evolve new properties. Further work could show why different enzymes are more or less active against different antibiotics.

On the Case of the Misplaced Hydrogens.

Few other elements play a more central role in biology than hydrogen. The interactions, bonding and movement of hydrogen atoms are central to biological catalysis, structure and function. Yet owing to the elusive nature of a single hydrogen atom few experimental and computational techniques can precisely determine its location. This is exemplified in short hydrogen bonds (SHBs) where the location of the hydrogen atom is indicative of the underlying strength of the bonds, which can vary from 1-5 kcal/mol in canonical hydrogen bonds, to an almost covalent nature in single-well hydrogen bonds. Owing to the often-times inferred position of hydrogen, the role of SHBs in biology has remained highly contested and debated. This has also led to discrepancies in computational, biochemical and structural studies of proteins thought to use SHBs in performing chemistry and stabilizing interactions. Herein, we discuss in detail two distinct examples, namely the conserved catalytic triad and the photoreceptor, photoactive yellow protein, where studies of these SHB-containing systems have permitted contextualization of the role these unique hydrogen bonds play in biology.

Visualizing Rev1 catalyze protein-template DNA synthesis.

During DNA replication, replicative DNA polymerases may encounter DNA lesions, which can stall replication forks. One way to prevent replication fork stalling is through the recruitment of specialized translesion synthesis (TLS) polymerases that have evolved to incorporate nucleotides opposite DNA lesions. Rev1 is a specialized TLS polymerase that bypasses abasic sites, as well as minor-groove and exocyclic guanine adducts. Lesion bypass is accomplished using a unique protein-template mechanism in which the templating base is evicted from the DNA helix and the incoming dCTP hydrogen bonds with an arginine side chain of Rev1. To understand the protein-template mechanism at an atomic level, we employed a combination of time-lapse X-ray crystallography, molecular dynamics simulations, and DNA enzymology on the Saccharomyces cerevisiae Rev1 protein. We find that Rev1 evicts the templating base from the DNA helix prior to binding the incoming nucleotide. Binding the incoming nucleotide changes the conformation of the DNA substrate to orient it for nucleotidyl transfer, although this is not coupled to large structural changes in Rev1 like those observed with other DNA polymerases. Moreover, we found that following nucleotide incorporation, Rev1 converts the pyrophosphate product to two monophosphates, which drives the reaction in the forward direction and prevents pyrophosphorolysis. Following nucleotide incorporation, the hydrogen bonds between the incorporated nucleotide and the arginine side chain are broken, but the templating base remains extrahelical. These postcatalytic changes prevent potentially mutagenic processive synthesis by Rev1 and facilitate dissociation of the DNA product from the enzyme.

Myelin Basic Protein Phospholipid Complexation Likely Competes with Deimination in Experimental Autoimmune Encephalomyelitis Mouse Model.

Multiple sclerosis has complex pathogenesis encompassing a variety of components (immunologic, genetic, and environmental). The autoimmunogenicity against the host's myelin basic protein is a major contributor. An increase in myelin basic protein deimination (a post-translational modification) and a change in phospholipid composition have been associated with multiple sclerosis. The interaction of myelin basic protein with phospholipids in the myelin membrane is an important contributor to the stability and maintenance of proper myelin sheath function. The study of this aspect of multiple sclerosis is an area that has yet to be fully explored and that the present study seeks to understand. Several biochemical methods, a capillary electrophoresis coupled system and mass spectrometry, were used in this study. These methods identified four specific phospholipids complexing with myelin basic protein. We show that lysophosphatidylcholine 18:1 provides a robust competitive effect against hyper-deimination. Our data suggest that lysophosphatidylcholine 18:1 has a different biochemical behavior when compared to other phospholipids and lysophosphatidylcholines 14:0, 16:0, and 18:0.

AP-endonuclease 1 sculpts DNA through an anchoring tyrosine residue on the DNA intercalating loop.

Base excision repair (BER) maintains genomic stability through the repair of DNA damage. Within BER, AP-endonuclease 1 (APE1) is a multifunctional enzyme that processes DNA intermediates through its backbone cleavage activity. To accomplish these repair activities, APE1 must recognize and accommodate several diverse DNA substrates. This is hypothesized to occur through a DNA sculpting mechanism where structural adjustments of the DNA substrate are imposed by the protein; however, how APE1 uniquely sculpts each substrate within a single rigid active site remains unclear. Here, we utilize structural and biochemical approaches to probe the DNA sculpting mechanism of APE1, specifically by characterizing a protein loop that intercalates the minor groove of the DNA (termed the intercalating loop). Pre-steady-state kinetics reveal a tyrosine residue within the intercalating loop (Y269) that is critical for AP-endonuclease activity. Using X-ray crystallography and molecular dynamics simulations, we determined the Y269 residue acts to anchor the intercalating loop on abasic DNA. Atomic force microscopy reveals the Y269 residue is required for proper DNA bending by APE1, providing evidence for the importance of this mechanism. We conclude that this previously unappreciated tyrosine residue is key to anchoring the intercalating loop and stabilizing the DNA in the APE1 active site.

Electrostatic Switching Controls Channel Dynamics of the Sensor Protein VirB10 in A. tumefaciens Type IV Secretion System.

Type IV secretion systems are large nanomachines assembled across the bacterial cell envelope for effector translocation and conjugation. VirB10 traverses the inner and outer membranes, sensing cellular signals for coordinating the conformational switch for pilus biogenesis and/or secretion. Mutations uncoupling secretion from pilus biogenesis were identified in Agrobacterium tumefaciens VirB10 including a gating defect mutation G272R that made VirB10 unresponsive to intracellular ATP, causing unregulated secretion of VirE2 in a contact-independent manner. Comparative long-timescale molecular dynamics of the wild type and G272R mutant of the A. tumefaciens VirB10CTD tetradecamer reveals how the G272R mutation locks the oligomer in a rigid conformation by swapping the ionic interactions between the loops from the ß-barrel close to the inner leaflet of the outer membrane. This electrostatic switching changes the allosteric communication pathway from the extracellular loop to the base of the barrel, suggesting that the local conformational dynamics in the loops can gate information across VirB10.

Insights into Structural and Dynamical Changes Experienced by Human RNase 6 upon Ligand Binding.

Ribonuclease 6 (RNase 6) is one of eight catalytically active human pancreatic-type RNases that belong to a superfamily of rapidly evolving enzymes. Like some of its human homologues, RNase 6 exhibits host defense properties such as antiviral and antibacterial activities. Recently solved crystal structures of this enzyme in its nucleotide-free form show the conservation of the prototypical kidney-shaped fold preserved among vertebrate RNases, in addition to revealing the presence of a unique secondary active site. In this study, we determine the structural and conformational properties experienced by RNase 6 upon binding to substrate and product analogues. We present the first crystal structures of RNase 6 bound to a nucleotide ligand (adenosine 5'-monophosphate), in addition to RNase 6 bound to phosphate ions. While the enzyme preserves B2 subsite ligand preferences, our results show a lack of typical B2 subsite interactions normally observed in homologous ligand-bound RNases. A comparison of the dynamical properties of RNase 6 in its apo-, substrate-, and product-bound states highlight the unique dynamical properties experienced on time scales ranging from nano- to milliseconds. Overall, our results confirm the specific evolutionary adaptation of RNase 6 relative to its unique catalytic and biological activities.

Application Health Monitoring for Extreme-scale Resiliency using Cooperative Fault Management.

Resiliency is and will be a critical factor in determining scientific productivity on current and exascale supercomputers, and beyond. Applications oblivious to and incapable of handling transient soft and hard errors could waste supercomputing resources or, worse, yield misleading scientific insights. We introduce a novel application-driven silent error detection and recovery strategy based on application health monitoring. Our methodology uses application output that follows known patterns as indicators of an application's health, and knowledge that violation of these patterns could be indication of faults. Information from system monitors that report hardware and software health status is used to corroborate faults. Collectively, this information is used by a fault coordinator agent to take preventive and corrective measures by applying computational steering to an application between checkpoints. This cooperative fault management system uses the Fault Tolerance Backplane as a communication channel. The benefits of this framework are demonstrated with two real application case studies, molecular dynamics and quantum chemistry simulations, on scalable clusters with simulated memory and I/O corruptions. The developed approach is general and can be easily applied to other applications.

Enzyme dynamics: Looking beyond a single structure.

Conventional understanding of how enzymes function strongly emphasizes the role of structure. However, increasing evidence clearly indicates that enzymes do not remain fixed or operate exclusively in or close to their native structure. Different parts of the enzyme (from individual residues to full domains) undergo concerted motions on a wide range of time-scales, including that of the catalyzed reaction. Information obtained on these internal motions and conformational fluctuations has so far uncovered and explained many aspects of enzyme mechanisms, which could not have been understood from a single structure alone. Although there is wide interest in understanding enzyme dynamics and its role in catalysis, several challenges remain. In addition to technical difficulties, the vast majority of investigations are performed in dilute aqueous solutions, where conditions are significantly different than the cellular milieu where a large number of enzymes operate. In this review, we discuss recent developments, several challenges as well as opportunities related to this topic. The benefits of considering dynamics as an integral part of the enzyme function can also enable new means of biocatalysis, engineering enzymes for industrial and medicinal applications.

SARS-CoV-2 Genomes From Oklahoma, United States.

Genomic sequencing has played a major role in understanding the pathogenicity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). With the current pandemic, it is essential that SARS-CoV-2 viruses are sequenced regularly to determine mutations and genomic modifications in different geographical locations. In this study, we sequenced SARS-CoV-2 from five clinical samples obtained in Oklahoma, United States during different time points of pandemic presence in the state. One sample from the initial days of the pandemic in the state and four during the peak in Oklahoma were sequenced. Previously reported mutations including D614G in S gene, P4715L in ORF1ab, S194L, R203K, and G204R in N gene were identified in the genomes sequenced in this study. Possible novel mutations were also detected in the S gene (G1167V), ORF1ab (A6269S and P3371S), ORF7b (T28I), and ORF8 (G96R). Phylogenetic analysis of the genomes showed similarity to other SARS-CoV-2 viruses reported from across the globe. Structural characterization indicates that the mutations in S gene possibly influences conformational flexibility and motion of the spike protein, and the mutations in N gene are associated with disordered linker region within the nucleocapsid protein.

Low-Barrier and Canonical Hydrogen Bonds Modulate Activity and Specificity of a Catalytic Triad.

The position, bonding and dynamics of hydrogen atoms in the catalytic centers of proteins are essential for catalysis. The role of short hydrogen bonds in catalysis has remained highly debated and led to establishment of several distinctive geometrical arrangements of hydrogen atoms vis-à-vis the heavier donor and acceptor counterparts, that is, low-barrier, single-well or short canonical hydrogen bonds. Here we demonstrate how the position of a hydrogen atom in the catalytic triad of an aminoglycoside inactivating enzyme leads to a thirty-fold increase in catalytic turnover. A low-barrier hydrogen bond is present in the enzyme active site for the substrates that are turned over the best, whereas a canonical hydrogen bond is found with the least preferred substrate. This is the first comparison of these hydrogen bonds involving an identical catalytic network, while directly demonstrating how active site electrostatics adapt to the electronic nature of substrates to tune catalysis.

Engineering Dynamic Surface Peptide Networks on ButyrylcholinesteraseG117H for Enhanced Organophosphosphorus Anticholinesterase Catalysis.

The single residue mutation of butyrylcholinesterase (BChEG117H) hydrolyzes a number of organophosphosphorus (OP) anticholinesterases. Whereas other BChE active site/proximal mutations have been investigated, none are sufficiently active to be prophylactically useful. In a fundamentally different computer simulations driven strategy, we identified a surface peptide loop (residues 278-285) exhibiting dynamic motions during catalysis and modified it via residue insertions. We evaluated these loop mutants using computer simulations, substrate kinetics, resistance to inhibition, and enzyme reactivation assays using both the choline ester and OP substrates. A slight but significant increase in reactivation was noted with paraoxon with one of the mutants, and changes in KM and catalytic efficiency were noted in others. Simulations suggested weaker interactions between OP versus choline substrates and the active site of all engineered versions of the enzyme. The results indicate that an improvement of OP anticholinesterase hydrolysis through surface loop engineering may be a more effective strategy in an enzyme with higher intrinsic OP compound hydrolase activity.