Search

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.Significance StatementThe SARS-CoV-2 proteins have constantly been evolving. These variants assist the virus to survive, adapt and evade the host immune responses. While the main focus has been on structural proteins like Spike, there is very limited structural and functional information on the effects of emerging mutations on other essential non-structural viral proteins. One such protein is NSP10, an essential cofactor for NSP14 and NSP16. This study demonstrates that NSP10 is more resistant to genetic variations than other SARS-CoV-2 non-structural proteins and that the presence of mutations conserve structural and dynamic changes in NSP10. The effects of naturally occurring mutations reflect the evolutionary relationship between structurally conserved essential cofactors, their function and the role they play in the survival of the virus.

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 (DFT) calculations indicate that both active site states are stable. The hydroperoxo formed is potentially an intermediate in the mechanism of hydrogen peroxide formation in the absence of substrate. We establish that the N-terminal amino group of the copper coordinating His1 remains doubly protonated directly following molecular oxygen reduction by copper. Aided by mining minima free energy calculations we establish His157 conformational flexibility in solution that is abolished by steric hindrance in the crystal. A neutron crystal structure of NcLPMO9D at low pH supports occlusion of the active site which prevents protonation of His157 at acidic conditions.

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 world-wide pandemic in a matter of weeks.1 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.

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.

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.

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.

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 homologs, 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 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 analogs. 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 B(2) subsite ligand preferences, our results show lack of typical B(2) subsite interactions normally observed in homologous ligand-bound RNases. Comparison of the dynamical properties of RNase 6 in its apo, substrate-, and product-bound states highlight 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.

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.

Understanding protein motions and their role in enzymatic reactions is an important and timely topic in enzymology. Protein motions that are involved in the chemical step of catalysis are particularly intriguing but difficult to identify. A global network of coupled residues in Escherichia coli dihydrofolate reductase (E. coli DHFR), which assists in catalyzing the chemical step, has previously been demonstrated through quantum mechanical/molecular mechanical and molecular dynamics simulations, as well as bioinformatic analyses. A few specific residues — M42, G121, F125 and I14 — were shown to function synergistically with measurements of single turnover rates and the temperature dependence of intrinsic kinetic isotope effects (KIEs(int)) of site directed mutants. The current study hypothesizes that the global network of residues involved in the chemical step is evolutionarily conserved and probes homologous residues of the potential global network in human DHFR through measurements of the temperature dependence of KIEs(int) and computer simulations based on empirical valence bond (EVB) method. We study the mutants M53W and S145V. Both of these remote residues are homologous to network residues in E. coli DHFR. Non-additive isotope effects on activation energy are observed between M53 and S145, indicating their synergistic effect on the chemical step in human DHFR, which suggests that both of these residues are part of a network affecting the chemical step in enzyme catalysis. This finding supports the hypothesis that human and E. coli DHFR share similar networks, consistent with evolutionary preservation of such networks.

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.

Dihydrofolate reductase (DHFR) reduces dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. Due to its role in one carbon metabolism, chromosomal DHFR is the target of the antibacterial drug, trimethoprim. Resistance to trimethoprim has resulted in a type II DHFR that is not structurally related to the chromosomal enzyme target. Because of its metabolic significance, understanding DHFR kinetics and ligand binding behavior in more cell-like conditions, where the total macromolecule concentration can be as great as 300 mg/mL, is important. The progress-curve kinetics and ligand binding properties of the drug target (chromosomal E. coli DHFR) and the drug resistant (R67 DHFR) enzymes were studied in the presence of macromolecular cosolutes. There were varied effects on NADPH oxidation and binding to the two DHFRs, with some cosolutes increasing affinity and others weakening binding. However, DHF binding and reduction in both DHFRs decreased in the presence of all cosolutes. The decreased binding of ligands is mostly attributed to weak associations with the macromolecules, as opposed to crowder effects on the DHFRs. Computer simulations found weak, transient interactions for both ligands with several proteins. The net charge of protein cosolutes correlated with effects on NADP(+) binding, with near neutral and positively charged proteins having more detrimental effects on binding. For DHF binding, effects correlated more with the size of binding pockets on the protein crowders. These non-specific interactions between DHFR ligands and proteins predict that the in vivo efficiency of DHFRs may be much lower than expected from their in vitro rates.

Even after a century of investigation, our understanding of how enzymes work remains far from complete. In particular, several factors that enable enzymes to achieve high catalytic efficiencies remain only poorly understood. A number of theories have been developed, which propose or reaffirm that enzymes work as structural scaffolds, serving to bring together and properly orient the participants so that the reaction can proceed; therefore, leading to enzymes being viewed as only passive participants in the catalyzed reaction. A growing body of evidence shows that enzymes are not rigid structures but are constantly undergoing a wide range of internal motions and conformational fluctuations. In this Perspective, on the basis of studies from our group, we discuss the emerging biophysical model of enzyme catalysis that provides a detailed understanding of the interconnection among internal protein motions, conformational substates, enzyme mechanisms, and the catalytic efficiency of enzymes. For a number of enzymes, networks of conserved residues that extend from the surface of the enzyme all the way to the active site have been discovered. These networks are hypothesized to serve as pathways of energy transfer that enables thermodynamical coupling of the surrounding solvent with enzyme catalysis and play a role in promoting enzyme function. Additionally, the role of enzyme structure and electrostatic effects has been well acknowledged for quite some time. Collectively, the recent knowledge gained about enzyme mechanisms suggests that the conventional paradigm of enzyme structure encoding function is incomplete and needs to be extended to structure encodes dynamics, and together these enzyme features encode function including catalytic rate acceleration.

The biliverdin reductase B (BLVRB) class of enzymes catalyze the NADPH-dependent reduction of multiple flavin substrates and are emerging as critical players in cellular redox regulation. However, the role of dynamics and allostery have not been addressed, prompting studies here that have revealed a position 15 Å away from the active site within human BLVRB (T164) that is inherently dynamic and can be mutated to control global micro-millisecond motions and function. By comparing the inherent dynamics through nuclear magnetic resonance (NMR) relaxation approaches of evolutionarily distinct BLVRB homologues and by applying our previously developed Relaxation And Single Site Multiple Mutations (RASSMM) approach that monitors both the functional and dynamic effects of multiple mutations to the single T164 site, we have discovered that the most dramatic mutagenic effects coincide with evolutionary changes and these modulate coenzyme binding. Thus, evolutionarily changing sites distal to the active site serve as dynamic “dials” to globally modulate motions and function. Despite the distal dynamic and functional coupling modulated by this site, micro-millisecond motions span an order of magnitude in their apparent kinetic rates of motions. Thus, global dynamics within BLVRB are a collection of partially coupled motions tied to catalytic function.

Graphical abstract, 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.

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., eLife digest 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.

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.

Biliverdin reductase B (BLVRB) family members are general flavin reductases critical in maintaining cellular redox with recent findings revealing that BLVRB alone can dictate cellular fate. However, as opposed to most enzymes, the BLVRB family remains enigmatic with an evolutionarily changing active site and unknown structural and functional consequences. Here, we applied a multi-faceted approach that combines X-ray crystallography, NMR and kinetics methods to elucidate the structural and functional basis of the evolutionarily changing BLVRB active site. Using a panel of three BLVRB isoforms (human, lemur and hyrax) and multiple human BLVRB mutants, our studies reveal a novel evolutionary mechanism where coenzyme ‘clamps’ formed by arginine side chains at two co-evolving positions within the active site serve to slow coenzyme release (Positions 14 and 78). We find that coenzyme release is further slowed by the weaker binding substrate, resulting in relatively slow turnover numbers. However, different BLVRB active sites imposed by either evolution or mutagenesis exhibit a surprising inverse relationship between coenzyme release and substrate turnover that is independent of the faster chemical step of hydride transfer also measured here. Collectively, our studies have elucidated the role of the evolutionarily changing BLVRB active site that serves to modulate coenzyme release and has revealed that coenzyme release is coupled to substrate turnover.

The genome of the hyperthermophile Thermotoga maritima contains three isoforms of maltose binding protein (MBP) that are high-affinity receptors for di-, tri-, and tetrasaccharides. Two of these proteins (tmMBP1 and tmMBP2) share significant sequence identity, approximately 90%, while the third (tmMBP3) shares less than 40% identity. MBP from Escherichia coli (ecMBP) shares 35% sequence identity with the tmMBPs. This subset of MBP isoforms offers an interesting opportunity to investigate the mechanisms underlying the evolution of substrate specificity and affinity profiles in a genome where redundant MBP genes are present. In this study, the X-ray crystal structures of tmMBP1, tmMBP2, and tmMBP3 are reported in the absence and presence of oligosaccharides. tmMBP1 and tmMBP2 have binding pockets that are larger than that of tmMBP3, enabling them to bind to larger substrates, while tmMBP1 and tmMBP2 also undergo substrate-induced hinge bending motions (∼52°) that are larger than that of tmMBP3 (∼35°). Small-angle X-ray scattering was used to compare protein behavior in solution, and computer simulations provided insights into dynamics of these proteins. Comparing quantitative protein-substrate interactions and dynamical properties of tmMBPs with those of the promiscuous ecMBP and disaccharide selective Thermococcus litoralis MBP provides insights into the features that enable selective binding. Collectively, the results provide insights into how the structure and dynamics of tmMBP homologues enable them to differentiate between a myriad of chemical entities while maintaining their common fold.

Optimal enzyme activity depends on a number of factors, including structure and dynamics. The role of enzyme structure is well recognized, however, the linkage between protein dynamics and enzyme activity has given rise to a contentious debate. We have developed an approach that uses an aqueous mixture of organic solvent to control the functionally relevant enzyme dynamics (without changing the structure), which in turn modulates the enzyme activity. Using this approach, we predicted that the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) from Escherichia coli in aqueous mixtures of isopropanol (IPA) with water will decrease by ~3 fold at 20% (v/v) IPA concentration. Stop flow kinetic measurements find that the pH-independent k(hydride) rate decreases by 2.2 fold. X-ray crystallographic enzyme structures shows no noticeable differences, while computational studies indicate that the transition-state and electrostatic effects were identical for water and mixed solvent conditions; and quasi-elastic neutron scattering studies show that the dynamical enzyme motions are suppressed. Our approach provides a unique avenue to modulating enzyme activity through changes in enzyme dynamics. Further it provides vital insights that show the altered motions of DHFR cause significant changes in the enzyme’s ability to access its functionally relevant conformational sub-states, explaining the decreased k(hydride) rate. This approach has important implications for obtaining fundamental insights into the role of rate-limiting dynamics in catalysis and as well as for enzyme engineering.

R67 dihydrofolate reductase (DHFR) is a homotetramer with a single active site pore and no sequence or structural homology with chromosomal DHFRs. The R67 enzyme provides resistance to trimethoprim, an active site-directed inhibitor of Escherichia coli DHFR. Sixteen to twenty N-terminal amino acids are intrinsically disordered in the R67 dimer crystal structure. Chymotrypsin cleavage of 16 N-terminal residues results in an active enzyme with a decreased stability. The space sampled by the disordered N-termini of R67 DHFR was investigated using small angle neutron scattering. From a combined analysis using molecular dynamics and the program SASSIE (http://www.smallangles.net/sassie/SASSIE_HOME.html), the apoenzyme displays a radius of gyration (Rg) of 21.46 ± 0.50 Å. Addition of glycine betaine, an osmolyte, does not result in folding of the termini as the Rg increases slightly to 22.78 ± 0.87 Å. SASSIE fits of the latter SANS data indicate that the disordered N-termini sample larger regions of space and remain disordered, suggesting they might function as entropic bristles. Pressure perturbation calorimetry also indicated that the volume of R67 DHFR increases upon addition of 10% betaine and decreased at 20% betaine because of the dehydration of the protein. Studies of the hydration of full-length R67 DHFR in the presence of the osmolytes betaine and dimethyl sulfoxide find around 1250 water molecules hydrating the protein. Similar studies with truncated R67 DHFR yield around 400 water molecules hydrating the protein in the presence of betaine. The difference of ∼900 waters indicates the N-termini are well-hydrated.

The ligand-induced conformational changes of periplasmic binding proteins (PBP) play a key role in the acquisition of metabolites in ATP binding cassette (ABC) transport systems. This conformational change allows for differential recognition of the ligand occupancy of the PBP by the ABC transporter. This minimizes futile ATP hydrolysis in the transporter, a phenomenon in which ATP hydrolysis is not coupled to metabolite transport. In many systems, the PBP conformational change is insufficient at eliminating futile ATP hydrolysis. Here we identify an additional state of the PBP that is also allosterically regulated by the ligand. Ligand binding to the homodimeric apo PBP leads to a tightening of the interface α-helices so that the hydrogen bonding pattern shifts to that of a 310 helix, in-turn altering the contacts and the dynamics of the protein interface so that the monomer exists in the presence of ligand.

Lytic polysaccharide monooxygenases have attracted vast attention owing to their abilities to disrupt glycosidic bonds via oxidation instead of hydrolysis and to enhance enzymatic digestion of recalcitrant substrates including chitin and cellulose. We have determined high-resolution X-ray crystal structures of an enzyme from Neurospora crassa in the resting state and of a copper(II) dioxo intermediate complex formed in the absence of substrate. X-ray crystal structures also revealed "pre-bound" molecular oxygen adjacent to the active site. An examination of protonation states enabled by neutron crystallography and density functional theory calculations identified a role for a conserved histidine in promoting oxygen activation. These results provide a new structural description of oxygen activation by substrate free lytic polysaccharide monooxygenases and provide insights that can be extended to reactivity in the enzyme-substrate complex.

The single residue mutation of butyrylcholinesterase (BChE(G117H)) hydrolyzes a number of organophosphosphorus (OP) anticholinesterases. While other active site/proximal mutations have been investigated, none are sufficiently active to be prophylactically useful. In a fundamentally different, computer simulation-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 K(M) and catalytic efficiency were noted in others. Simulations suggested weaker interactions between OP versus choline substrates and the active site of all enzymes. The results indicate that 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.

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.

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-a-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.

Folate, or vitamin B9, is an important compound in one-carbon metabolism. Previous studies have found weaker binding of dihydrofolate to dihydrofolate reductase in the presence of osmolytes. In other words, osmolytes are more difficult to remove from the dihydrofolate solvation shell than water; this shifts the equilibrium toward the free ligand and protein species. This study uses vapor-pressure osmometry to explore the interaction of folate with the model osmolyte, glycine betaine. This method yields a preferential interaction potential (μ23/RT value). This value is concentration-dependent as folate dimerizes. The μ23/RT value also tracks the deprotonation of folate’s N3–O4 keto–enol group, yielding a pKa of 8.1. To determine which folate atoms interact most strongly with betaine, the interaction of heterocyclic aromatic compounds (as well as other small molecules) with betaine was monitored. Using an accessible surface area approach coupled with osmometry measurements, deconvolution of the μ23/RT values into α values for atom types was achieved. This allows prediction of μ23/RT values for larger molecules such as folate. Molecular dynamics simulations of folate show a variety of structures from extended to L-shaped. These conformers possess μ23/RT values from −0.18 to 0.09 m–1, where a negative value indicates a preference for solvation by betaine and a positive value indicates a preference for water. This range of values is consistent with values observed in osmometry and solubility experiments. As the average predicted folate μ23/RT value is near zero, this indicates folate interacts almost equally well with betaine and water. Specifically, the glutamate tail prefers to interact with water, while the aromatic rings prefer betaine. In general, the more protonated species in our small molecule survey interact better with betaine as they provide a source of hydrogens (betaine is not a hydrogen bond donor). Upon deprotonation of the small molecule, the preference swings toward water interaction because of its hydrogen bond donating capacities.

International audience; After publication of this article, the authors noticed that a 15N-13C dimension error was unwillingly coded in the 3D NMR spectrum "fid.com" processing script used to perform backbone assignments for this enzyme. The authors noticed that some OBS, CAR and LAB values in the "fid.com" had been switched in the y and z dimensions, probably resulting from a wrong NMRPipe selection when reading the Varian NMR experimental parameters. They have carefully re-processed, re-analyzed, re-assigned, in addition to checking all scripts to evaluate the extent of this processing error on the published assignments. Authors determined that the "fid.com" error resulted in a significant number of incorrect backbone resonance assignments, requiring us to issue corrections in Figs. 2, 3 and 4 of this published manuscript, in addition to Table S1. New versions of these figures and table are provided below. The corresponding BMRB entry has also been revised. The authors note that these modifications do not change the global message, conclusions, and molecular dynamics simulations presented in this article. The authors are grateful to David N. Bernard (INRS) for help with finding and correcting these errors.

Enzyme catalysis is a complex process involving several steps along the reaction coordinates, including substrate recognition and binding, chemical transformation, and product release. Evidence continues to emerge linking the functional and evolutionary role of conformational exchange processes in optimal catalytic activity. Ligand binding changes the conformational landscape of enzymes, inducing long-range conformational rearrangements. Using functionally distinct members of the pancreatic ribonuclease superfamily as a model system, we characterized the structural and conformational changes associated with the binding of two mononucleotide ligands. By combining NMR chemical shift titration experiments with the chemical shift projection analysis (CHESPA) and relaxation dispersion experiments, we show that biologically distinct members of the RNase superfamily display discrete chemical shift perturbations upon ligand binding that are not conserved even in structurally related members. Amino acid networks exhibiting coordinated chemical shift displacements upon binding of the two ligands are unique to each of the RNases analyzed. Our results reveal the contribution of conformational rearrangements to the observed chemical shift perturbations. These observations provide important insights into the contribution of the different ligand binding specificities and effects of conformational exchange on the observed perturbations associated with ligand binding for functionally diverse members of the pancreatic RNase superfamily.