Structural Characterization of Arginine Fingers: Identification of an Arginine Finger for the Pyrophosphatase dUTPases

Emerging Mechanisms of Metal-Catalyzed RNA and DNA Modifications

Metal ions play a critical role in various chemical, biological, and environmental processes. This review reports on emerging chemical mechanisms in the catalysis of DNA and RNA. We provide an overview of the metal-dependent mechanisms of DNA cleavage in CRISPR (clustered regularly interspaced short palindromic repeats)-Cas systems that are transforming life sciences through genome editing technologies, and showcase intriguing metal-dependent mechanisms of RNA cleavages. We show that newly discovered CRISPR-Cas complexes operate as protein-assisted ribozymes, highlighting RNA's versatility and the enhancement of CRISPR-Cas functions through strategic metal ion use. We demonstrate the power of computer simulations in observing chemical processes as they unfold and in advancing structural biology through innovative approaches for refining cryo-electron microscopy maps. Understanding metal ion involvement in nucleic acid catalysis is crucial for advancing genome editing, aiding therapeutic interventions for genetic disorders, and improving the editing tools' specificity and efficiency.

The homodimerization domain of the Stl repressor is crucial for efficient inhibition of mycobacterial dUTPase

The dUTPase is a key DNA repair enzyme in Mycobacterium tuberculosis, and it may serve as a novel promising anti-tuberculosis target. Stl repressor from Staphylococcus aureus was shown to bind to and inhibit dUTPases from various sources, and its expression in mycobacterial cells interfered with cell growth. To fine-tune and optimize Stl-induced inhibition of mycobacterial dUTPase, we aimed to decipher the molecular details of this interaction. Structural background of the complex between dUTPase and a truncated Stl lacking the repressor C-terminal homodimerization domain has been described, however, the effects of this truncation of Stl on enzyme binding and inhibition are still not known. Using several independent biophysical, structural and enzyme kinetic methods, here we show that lack of the repressor homodimerization domain strongly perturbs both enzyme binding and inhibition. We also investigated the role of a mycobacteria-specific loop in the Stl-interaction. Our results show that removal of this loop leads to a ten-fold increase in the apparent inhibition constant of Stl. We present a high-resolution three-dimensional structure of mycobacterial dUTPase lacking the genus-specific loop for structural insight. Our present data suggest that potent inhibition of mycobacterial dUTPase by Stl requires the wild-type full-length protein context.

Systematic QM/MM Study for Predicting 31P NMR Chemical Shifts of Adenosine Nucleotides in Solution and Stages of ATP Hydrolysis in a Protein Environment

NMR (nuclear magnetic resonance) spectroscopy allows for important atomistic insights into the structure and dynamics of biological macromolecules; however, reliable assignments of experimental spectra are often difficult. Herein, quantum mechanical/molecular mechanical (QM/MM) calculations can provide crucial support. A major problem for the simulations is that experimental NMR signals are time-averaged over much longer time scales, and since computed chemical shifts are highly sensitive to local changes in the electronic and structural environment, sufficiently large averages over representative structural ensembles are essential. This entails high computational demands for reliable simulations. For NMR measurements in biological systems, a nucleus of major interest is 31P since it is both highly present (e.g., in nucleic acids) and easily observable. The focus of our present study is to develop a robust and computationally cost-efficient framework for simulating 31P NMR chemical shifts of nucleotides. We apply this scheme to study the different stages of the ATP hydrolysis reaction catalyzed by p97. Our methodology is based on MM molecular dynamics (MM-MD) sampling, followed by QM/MM structure optimizations and NMR calculations. Overall, our study is one of the most comprehensive QM-based 31P studies in a protein environment and the first to provide computed NMR chemical shifts for multiple nucleotide states in a protein environment. This study sheds light on a process that is challenging to probe experimentally and aims to bridge the gap between measured and calculated NMR spectroscopic properties.

Structures and Functions of the Human GATOR1 Complex

Eukaryotic cells coordinate available nutrients with their growth through the mechanistic target of rapamycin complex 1 (mTORC1) pathway, in which numerous evolutionarily conserved protein complexes survey and transmit nutrient inputs toward mTORC1. mTORC1 integrates these inputs and activates downstream anabolic or catabolic programs that are in tune with cellular needs, effectively maintaining metabolic homeostasis. The GAP activity toward Rags-1 (GATOR1) protein complex is a critical negative regulator of the mTORC1 pathway and, in the absence of amino acid inputs, is activated to turn off mTORC1 signaling. GATOR1-mediated inhibition of mTORC1 signaling is tightly regulated by an ensemble of protein complexes that antagonize or promote its activity in response to the cellular nutrient environment. Structural, biochemical, and biophysical studies of the GATOR1 complex and its interactors have advanced our understanding of how it regulates cellular metabolism when amino acids are limited. Here, we review the current research with a focus on GATOR1 structure, its enzymatic mechanism, and the growing group of proteins that regulate its activity. Finally, we discuss the implication of GATOR1 dysregulation in physiology and human diseases.

Mechanism-Based Redesign of GAP to Activate Oncogenic Ras

Ras GTPases play a crucial role in cell signaling pathways. Mutations of the Ras gene occur in about one third of cancerous cell lines and are often associated with detrimental clinical prognosis. Hot spot residues Gly12, Gly13, and Gln61 cover 97% of oncogenic mutations, which impair the enzymatic activity in Ras. Using QM/MM free energy calculations, we present a two-step mechanism for the GTP hydrolysis catalyzed by the wild-type Ras.GAP complex. We found that the deprotonation of the catalytic water takes place via the Gln61 as a transient Brønsted base. We also determined the reaction profiles for key oncogenic Ras mutants G12D and G12C using QM/MM minimizations, matching the experimentally observed loss of catalytic activity, thereby validating our reaction mechanism. Using the optimized reaction paths, we devised a fast and accurate procedure to design GAP mutants that activate G12D Ras. We replaced GAP residues near the active site and determined the activation barrier for 190 single mutants. We furthermore built a machine learning for ultrafast screening, by fast prediction of the barrier heights, tested both on the single and double mutations. This work demonstrates that fast and accurate screening can be accomplished via QM/MM reaction path optimizations to design protein sequences with increased catalytic activity. Several GAP mutations are predicted to re-enable catalysis in oncogenic G12D, offering a promising avenue to overcome aberrant Ras-driven signal transduction by activating enzymatic activity instead of inhibition. The outlined computational screening protocol is readily applicable for designing ligands and cofactors analogously.

Molecular detection and genetic characteristics of a novel porcine circovirus (porcine circovirus 4) and porcine reproductive and respiratory syndrome virus in Shaanxi and Henan Provinces of China

Porcine circovirus 4 (PCV4) is a recently discovered circovirus that was first reported in 2019 in several pigs with severe clinical disease in Hunan province of China, and also identified in pigs infected with porcine reproductive and respiratory syndrome virus (PRRSV). To further investigate the epidemic profile and genetic characteristics of the two viruses, 150 clinical samples were collected from 9 swine farms in Shaanxi and Henan provinces of China, and a SYBR Green I-based duplex quantitative real-time polymerase chain reaction (qPCR) was developed for detecting PCV4 and PRRSV simultaneously. The results showed the limits of detection were 41.1 copies/μL and 81.5 copies/μL for PCV4 and PRRSV, respectively. The detection rates of PCV4 and PRRSV were 8.00% (12/150) and 12.00% (18/150) respectively, and a case of co-infection with PCV4 and PRRSV was found in the lung tissue of a suckling pig with respiratory symptom. Subsequently, the complete genomic sequences of five PCV4 strains were obtained, of which one PCV4 strain (SX-ZX) was from Shaanxi province, and these strains were 1770 nucleotides in length and had 97.7%-99.4% genomic identity with 59 PCV4 reference strains. The genome characteristic of the SX-ZX strain was evaluated from three aspects, a "stem-loop" structure, ORF1 and ORF2. As essential elements for the replication, the 17-bp iterative sequence was predicted as the stem structure, in which three non-tandem hexamers were found at downstream with H1/H2 (12-CGGCACACTTCGGCAC-27) as the minimal binding site. Three of the five PCV4 strains were clustered into PCV4b, which was composed of Suidae, fox, dairy cow, dog and raccoon dog. Phylogenetic analysis revealed that seven PRRSV strains from the present study were clustered into the PRRSV-2 genotype. Collectively, these data extend our understanding of the genome characteristic of PCV4 as well as the molecular epidemiology and the genetic profile of PCV4 and PRRSV.

Continuous millisecond conformational cycle of a DEAH box helicase reveals control of domain motions by atomic-scale transitions

Helicases are motor enzymes found in every living organism and viruses, where they maintain the stability of the genome and control against false recombination. The DEAH-box helicase Prp43 plays a crucial role in pre-mRNA splicing in unicellular organisms by translocating single-stranded RNA. The molecular mechanisms and conformational transitions of helicases are not understood at the atomic level. We present a complete conformational cycle of RNA translocation by Prp43 in atomic detail based on molecular dynamics simulations. To enable the sampling of such complex transition on the millisecond timescale, we combined two enhanced sampling techniques, namely simulated tempering and adaptive sampling guided by crystallographic data. During RNA translocation, the center-of-mass motions of the RecA-like domains followed the established inchworm model, whereas the domains crawled along the RNA in a caterpillar-like movement, suggesting an inchworm/caterpillar model. However, this crawling required a complex sequence of atomic-scale transitions involving the release of an arginine finger from the ATP pocket, stepping of the hook-loop and hook-turn motifs along the RNA backbone, and several others. These findings highlight that large-scale domain dynamics may be controlled by complex sequences of atomic-scale transitions.

The Electronic Structure of Genome Editors from the First Principles

Genome editing based on the CRISPR-Cas9 system has paved new avenues for medicine, pharmaceutics, biotechnology, and beyond. This article reports the role of first-principles (ab-initio) molecular dynamics (MD) in the CRISPR-Cas9 revolution, achieving a profound understanding of the enzymatic function and offering valuable insights for enzyme engineering. We introduce the methodologies and explain the use of ab-initio MD simulations to characterize the two-metal dependent mechanism of DNA cleavage in the RuvC domain of the Cas9 enzyme, and how a second catalytic domain, HNH, cleaves the target DNA with the aid of a single metal ion. A detailed description of how ab-initio MD is combined with free-energy methods - i.e., thermodynamic integration and metadynamics - to break and form chemical bonds is given, explaining the use of these methods to determine the chemical landscape and establish the catalytic mechanism in CRISPR-Cas9. The critical role of classical methods is also discussed, explaining theory and application of constant pH MD simulations, used to accurately predict the catalytic residues' protonation states. Overall, first-principles methods are shown to unravel the electronic structure of the Cas9 enzyme, providing valuable insights that can serve for the design of genome editing tools with improved catalytic efficiency or controllable activity.

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

Structural Dynamics and Catalytic Mechanism of ATP13A2 (PARK9) from Simulations

ATP13A2 is a gene encoding a protein of the P5B subfamily of ATPases and is a PARK gene. Molecular defects of the gene are mainly associated with variations of Parkinson's disease (PD). Despite the established importance of the protein in regulating neuronal integrity, the three-dimensional structure of the protein currently remains unresolved crystallographically. We have modeled the structure and reactivity of the full-length protein in its E1-ATP state. Using molecular dynamics (MD), quantum cluster, and quantum mechanical/molecular mechanical (QM/MM) methods, we aimed at describing the main catalytic reaction, leading to the phosphorylation of Asp513. Our MD simulations suggest that two positively charged Mg²⁺ cations are present at the active site during the catalytic reaction, stabilizing a specific triphosphate binding mode. Using QM/MM calculations, we subsequently calculated the reaction profiles for the phosphoryl transfer step in the presence of one and two Mg²⁺ cations. The calculated barrier heights in both cases are found to be ∼12.5 and 7.5 kcal mol^-1, respectively. We elucidated details of the catalytically competent ATP conformation and the binding mode of the second Mg²⁺ cofactor. We also examined the role of the conserved Arg686 and Lys859 catalytic residues. We observed that by significantly lowering the barrier height of the ATP cleavage reaction, Arg686 had major effect on the reaction. The removal of Arg686 increased the barrier height for the ATP cleavage by more than 5.0 kcal mol^-1 while the removal of key electrostatic interactions created by Lys859 to the γ-phosphate and Asp513 destabilizes the reactant state. When missense mutations occur in close proximity to an active site residue, they can interfere with the barrier height of the reaction, which can halt the normal enzymatic rate of the protein. We also found large binding pockets in the full-length structure, including a transmembrane domain pocket, which is likely where the ATP13A2 cargo binds.

Mechanism of Guanosine Triphosphate Hydrolysis by the Visual Proteins Arl3-RP2: Free Energy Reaction Profiles Computed with Ab Initio Type QM/MM Potentials

We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme-substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.

Catalytic mechanism of the colistin resistance protein MCR-1

The mcr-1 gene encodes a membrane-bound Zn2+-metalloenzyme, MCR-1, which catalyses phosphoethanolamine transfer onto bacterial lipid A, making bacteria resistant to colistin, a last-resort antibiotic. Mechanistic understanding of this process remains incomplete. Here, we investigate possible catalytic pathways using DFT and ab initio calculations on cluster models and identify a complete two-step reaction mechanism. The first step, formation of a covalent phosphointermediate via transfer of phosphoethanolamine from a membrane phospholipid donor to the acceptor Thr285, is rate-limiting and proceeds with a single Zn2+ ion. The second step, transfer of the phosphoethanolamine group to lipid A, requires an additional Zn2+. The calculations suggest the involvement of the Zn2+ orbitals directly in the reaction is limited, with the second Zn2+ acting to bind incoming lipid A and direct phosphoethanolamine addition. The new level of mechanistic detail obtained here, which distinguishes these enzymes from other phosphotransferases, will aid in the development of inhibitors specific to MCR-1 and related bacterial phosphoethanolamine transferases.

Light-Induced Change of Arginine Conformation Modulates the Rate of Adenosine Triphosphate to Cyclic Adenosine Monophosphate Conversion in the Optogenetic System Containing Photoactivated Adenylyl Cyclase

We report the first computational characterization of an optogenetic system composed of two photosensing BLUF (blue light sensor using flavin adenine dinucleotide) domains and two catalytic adenylyl cyclase (AC) domains. Conversion of adenosine triphosphate (ATP) to the reaction products, cyclic adenosine monophosphate (cAMP) and pyrophosphate (PPi), catalyzed by ACs initiated by excitation in photosensing domains has emerged in the focus of modern optogenetic applications because of the request in photoregulated enzymes that modulate cellular concentrations of signaling messengers. The photoactivated AC from the soil bacterium Beggiatoa sp. (bPAC) is an important model showing a considerable increase in the ATP to cAMP conversion rate in the catalytic domain after the illumination of the BLUF domain. The 1 μs classical molecular dynamics simulations reveal that the activation of the BLUF domain leading to tautomerization of Gln49 in the chromophore-binding pocket results in switching of the position of the side chain of Arg278 in the active site of AC. Allosteric signal transmission pathways between Gln49 from BLUF and Arg278 from AC were revealed by the dynamical network analysis. The Gibbs energy profiles of the ATP → cAMP + PPi reaction computed using QM(DFT(ωB97X-D3/6-31G**))/MM(CHARMM) molecular dynamics simulations for both Arg278 conformations in AC clarify the reaction mechanism. In the light-activated system, the corresponding arginine conformation stabilizes the pentacoordinated phosphorus of the α-phosphate group in the transition state, thus lowering the activation energy. Simulations of the bPAC system with the Tyr7Phe replacement in the BLUF demonstrate occurrence of both arginine conformations in an equal ratio, explaining the experimentally observed intermediate catalytic activity of the bPAC-Y7F variant as compared with the dark and light states of the wild-type bPAC.

Structural comparisons of host and African swine fever virus dUTPases reveal new clues for inhibitor development

African swine fever, caused by the African swine fever virus (ASFV), is among the most significant swine diseases. There are currently no effective treatments against ASFV. ASFV contains a gene encoding a dUTPase (E165R), which is required for viral replication in swine macrophages, making it an attractive target for inhibitor development. However, the full structural details of the ASFV dUTPase and those of the comparable swine enzyme are not available, limiting further insights. Herein, we determine the crystal structures of ASFV dUTPase and swine dUTPase in both their ligand-free and ligand-bound forms. We observe that the swine enzyme employs a classical dUTPase architecture made up of three-subunit active sites, whereas the ASFV enzyme employs a novel two-subunit active site. We then performed a comparative analysis of all dUTPase structures uploaded in the Protein Data Bank (PDB), which showed classical and non-classical types were mainly determined by the C-terminal β-strand orientation, and the difference was mainly related to the four amino acids behind motif IV. Thus, our study not only explains the reason for the structural diversity of dUTPase but also reveals how to predict dUTPase type, which may have implications for the dUTPase family. Finally, we tested two dUTPase inhibitors developed for the Plasmodium falciparum dUTPase against the swine and ASFV enzymes. One of these compounds inhibited the ASFV dUTPase at low micromolar concentrations (Kd = 15.6 μM) and with some selectivity (∼2x) over swine dUTPase. In conclusion, our study expands our understanding of the dUTPase family and may aid in the development of specific ASFV inhibitors.

Long-range dynamic correlations regulate the catalytic activity of the bacterial tyrosine kinase Wzc

BY-kinases represent a highly conserved family of protein tyrosine kinases unique to bacteria without eukaryotic orthologs. BY-kinases are regulated by oligomerization-enabled transphosphorylation on a C-terminal tyrosine cluster through a process with sparse mechanistic detail. Using the catalytic domain (CD) of the archetypal BY-kinase, Escherichia coli Wzc, and enhanced-sampling molecular dynamics simulations, isothermal titration calorimetry and nuclear magnetic resonance measurements, we propose a mechanism for its activation and nucleotide exchange. We find that the monomeric Wzc CD preferentially populates states characterized by distortions at its oligomerization interfaces and by catalytic element conformations that allow high-affinity interactions with ADP but not with ATP·Mg2+ We propose that oligomer formation stabilizes the intermonomer interfaces and results in catalytic element conformations suitable for optimally engaging ATP·Mg2+, facilitating exchange with bound ADP. This sequence of events, oligomerization, i.e., substrate binding, before engaging ATP·Mg2+, facilitates optimal autophosphorylation by preventing a futile cycle of ATP hydrolysis.

Tetramerisation of the CRISPR ring nuclease Crn3/Csx3 facilitates cyclic oligoadenylate cleavage

Type III CRISPR systems detect foreign RNA and activate the cyclase domain of the Cas10 subunit, generating cyclic oligoadenylate (cOA) molecules that act as a second messenger to signal infection, activating nucleases that degrade the nucleic acid of both invader and host. This can lead to dormancy or cell death; to avoid this, cells need a way to remove cOA from the cell once a viral infection has been defeated. Enzymes specialised for this task are known as ring nucleases, but are limited in their distribution. Here, we demonstrate that the widespread CRISPR associated protein Csx3, previously described as an RNA deadenylase, is a ring nuclease that rapidly degrades cyclic tetra-adenylate (cA4). The enzyme has an unusual cooperative reaction mechanism involving an active site that spans the interface between two dimers, sandwiching the cA4 substrate. We propose the name Crn3 (CRISPR associated ring nuclease 3) for the Csx3 family.

Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H+ and Mg2+ ions in enzyme reactions

Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg²⁺ ions, and recent advances in their transient involvement in the catalytic mechanisms.

Structural Insight into African Swine Fever Virus dUTPase Reveals a Novel Folding Pattern in the dUTPase Family

The African swine fever virus (ASFV) is the deadly pathogen of African swine fever (ASF) that induces high mortality, approaching 100% in domestic pigs, causes enormous losses to the global pig industry, and threatens food security. Currently, there is no effective treatment or preventive countermeasure. dUTPases (deoxyuridine 5'-triphosphate pyrophosphatases) are ubiquitous enzymes that are essential for the hydrolysis of dUTP and prevent the misincorporation of dUTP into newly synthesized DNA. Here, we present the crystal structures of the ASFV dUTPase in complex with the product dUMP and cofactor Mg²⁺ at a resolution of 2.2 Å. We observed that a unique "turning point" at G125 plays an unexpected critical role in the swapping region of the C-terminal segment, which is further stabilized by the interactions of the last C-terminal β strand with the β1 and β2 strands, thereby positioning the catalytic motif 5 into the active site of its own subunit instead of into a third subunit. Therefore, the ASFV dUTPase employs a novel two-subunit active site that is different than the classic trimeric dUTPase active site, which is composed of all three subunits. Meanwhile, further results confirmed that the configuration of motifs 1 to 5 has high structural homology with and a catalytic mechanism similar to that of the known trimeric dUTPases. In general, our study expands the information not only on the structural diversity of the conserved dUTPase family but also on the details needed to utilize this dUTPase as a novel target in the treatment of ASF.IMPORTANCE African swine fever virus (AFSV), a large enveloped double-stranded DNA virus, causes a deadly infection in domestic pigs. In addition to Africa, Europe, and South America, countries in Asia, such as China, Vietnam, and Mongolia, have suffered the hazards posed by ASFV outbreaks in recent years. Until now, there has been no vaccine for protection from ASFV infection or effective treatments to cure ASF. Here, we solved the crystal structure of the ASFV dUTPase-dUMP-Mg²⁺ complex. The ASFV dUTPase displays a noncanonical folding pattern that differs from that of the classic homotrimeric dUTPase, in which the active site is composed of two subunits. In addition, several nonconserved residues within the 3-fold axis channel play a vital role in ASFV dUTPase homotrimer stability. Our finding on these unique structural features of the ASFV dUTPase could be explored for the design of potential specific inhibitors that target this unique enzyme.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

Controlling the Revolving and Rotating Motion Direction of Asymmetric Hexameric Nanomotor by Arginine Finger and Channel Chirality

Nanomotors in nanotechnology are as important as engines in daily life. Many ATPases are nanoscale biomotors classified into three categories based on the motion mechanisms in transporting substrates: linear, rotating, and the recently discovered revolving motion. Most biomotors adopt a multisubunit ring-shaped structure that hydrolyzes ATP to generate force. How these biomotors control the motion direction and regulate the sequential action of their multiple subunits is intriguing. Many ATPases are hexameric with each monomer containing a conserved arginine finger. This review focuses on recent findings on how the arginine finger controls motion direction and coordinates adjacent subunit interactions in both revolving and rotating biomotors. Mechanisms of intersubunit interactions and sequential movements of individual subunits are evidenced by the asymmetrical appearance of one dimer and four monomers in high-resolution structural complexes. The arginine finger is situated at the interface of two subunits and extends into the ATP binding pocket of the downstream subunit. An arginine finger mutation results in deficiency in ATP binding/hydrolysis, substrate binding, and transport, highlighting the importance of the arginine finger in regulating energy transduction and motor function. Additionally, the roles of channel chirality and channel size are discussed as related to controlling one-way trafficking and differentiating the revolving and rotating mechanisms. Finally, the review concludes by discussing the conformational changes and entropy conversion triggered by ATP binding/hydrolysis, offering a view different from the traditional concept of ATP-mediated mechanochemical energy coupling. The elucidation of the motion mechanism and direction control in ATPases could facilitate nanomotor fabrication in nanotechnology.

Structure and Dynamics of the CRISPR-Cas9 Catalytic Complex

CRISPR-Cas9 is a bacterial immune system with exciting applications for genome editing. In spite of extensive experimental characterization, the active site chemistry of the RuvC domain-which performs DNA cleavages-has remained elusive. Its knowledge is key for structure-based engineering aimed at improving DNA cleavages. Here, we deliver an in-depth characterization by using quantum-classical (QM/MM) molecular dynamics (MD) simulations and a Gaussian accelerated MD method, coupled with bioinformatics analysis. We disclose a two-metal aided architecture in the RuvC active site, which is poised to operate DNA cleavages, in analogy with other DNA/RNA processing enzymes. The conformational dynamics of the RuvC domain further reveals that an "arginine finger" stably contacts the scissile phosphate, with the function of stabilizing the active complex. Remarkably, the formation of a catalytically competent state of the RuvC domain is only observed upon the conformational activation of the other nuclease domain of CRISPR-Cas9-i.e., the HNH domain-such allowing concerted cleavages of double stranded DNA. This structure is in agreement with the available experimental data and remarkably differs from previous models based on classical mechanics, demonstrating also that only quantum mechanical simulations can accurately describe the metal-aided active site in CRISPR-Cas9. This fully catalytic structure-in which both the HNH and RuvC domains are prone to perform DNA cleavages-constitutes a stepping-stone for understanding DNA cleavage and specificity. It calls for novel experimental verifications and offers the structural foundations for engineering efforts aimed at improving the genome editing capability of CRISPR-Cas9.

Eukaryotic Circular Rep-Encoding Single-Stranded DNA (CRESS DNA) Viruses: Ubiquitous Viruses With Small Genomes and a Diverse Host Range

While single-stranded DNA (ssDNA) was once thought to be a relatively rare genomic architecture for viruses, modern metagenomics sequencing has revealed circular ssDNA viruses in most environments and in association with diverse hosts. In particular, circular ssDNA viruses encoding a homologous replication-associated protein (Rep) have been identified in the majority of eukaryotic supergroups, generating interest in the ecological effects and evolutionary history of circular Rep-encoding ssDNA viruses (CRESS DNA) viruses. This review surveys the explosion of sequence diversity and expansion of eukaryotic CRESS DNA taxonomic groups over the last decade, highlights similarities between the well-studied geminiviruses and circoviruses with newly identified groups known only through their genome sequences, discusses the ecology and evolution of eukaryotic CRESS DNA viruses, and speculates on future research horizons.

Structural model of human dUTPase in complex with a novel proteinaceous inhibitor

Human deoxyuridine 5'-triphosphate nucleotidohydrolase (dUTPase), essential for DNA integrity, acts as a survival factor for tumor cells and is a target for cancer chemotherapy. Here we report that the Staphylococcal repressor protein Stl_SaPIBov1 (Stl) forms strong complex with human dUTPase. Functional analysis reveals that this interaction results in significant reduction of both dUTPase enzymatic activity and DNA binding capability of Stl. We conducted structural studies to understand the mechanism of this mutual inhibition. Small-angle X-ray scattering (SAXS) complemented with hydrogen-deuterium exchange mass spectrometry (HDX-MS) data allowed us to obtain 3D structural models comprising a trimeric dUTPase complexed with separate Stl monomers. These models thus reveal that upon dUTPase-Stl complex formation the functional homodimer of Stl repressor dissociates, which abolishes the DNA binding ability of the protein. Active site forming dUTPase segments were directly identified to be involved in the dUTPase-Stl interaction by HDX-MS, explaining the loss of dUTPase activity upon complexation. Our results provide key novel structural insights that pave the way for further applications of the first potent proteinaceous inhibitor of human dUTPase.

The dUTPase of white spot syndrome virus assembles its active sites in a noncanonical manner

dUTPases are essential enzymes for maintaining genome integrity and have recently been shown to play moonlighting roles when containing extra sequences. Interestingly, the trimeric dUTPase of white spot syndrome virus (wDUT) harbors a sequence insert at the position preceding the C-terminal catalytic motif V (pre-V insert), rarely seen in other dUTPases. However, whether this extra sequence endows wDUT with additional properties is unknown. Herein, we present the crystal structures of wDUT in both ligand-free and ligand-bound forms. We observed that the pre-V insert in wDUT forms an unusual β-hairpin structure in the domain-swapping region and thereby facilitates a unique orientation of the adjacent C-terminal segment, positioning the catalytic motif V onto the active site of its own subunit instead of a third subunit. Consequently, wDUT employs two-subunit active sites, unlike the widely accepted paradigm that the active site of trimeric dUTPase is contributed by all three subunits. According to results from local structural comparisons, the active-site configuration of wDUT is similar to that of known dUTPases. However, we also found that residues in the second-shell region of the active site are reconfigured in wDUT as an adaption to its unique C-terminal orientation. We also show that deletion of the pre-V insert significantly reduces wDUT's enzymatic activity and thermal stability. We hypothesize that this rare structural arrangement confers additional functionality to wDUT. In conclusion, our study expands the structural diversity in the conserved dUTPase family and illustrates how sequence insertion and amino acid substitution drive protein evolution cooperatively.

Divalent Metal Ion Activation of a Guanine General Base in the Hammerhead Ribozyme: Insights from Molecular Simulations

The hammerhead ribozyme is a well-studied nucleolytic ribozyme that catalyzes the self-cleavage of the RNA phosphodiester backbone. Despite experimental and theoretical efforts, key questions remain about details of the mechanism with regard to the activation of the nucleophile by the putative general base guanine (G12). Straightforward interpretation of the measured activity-pH data implies the pK_a value of the N1 position in the G12 nucleobase is significantly shifted by the ribozyme environment. Recent crystallographic and biochemical work has identified pH-dependent divalent metal ion binding at the N7/O6 position of G12, leading to the hypothesis that this binding mode could induce a pK_a shift of G12 toward neutrality. We present computational results that support this hypothesis and provide a model that unifies the interpretation of available structural and biochemical data and paints a detailed mechanistic picture of the general base step of the reaction. Experimentally testable predictions are made for mutational and rescue effects on G12, which will give further insights into the catalytic mechanism. These results contribute to our growing knowledge of the potential roles of divalent metal ions in RNA catalysis.