Tying down the arm in Bacillus dUTPase: structure and mechanism

Structural Outlier Detection and Zernike-Canterakis Moments for Molecular Surface Meshes-Fast Implementation in Python

Object retrieval systems measure the degree of similarity of the shape of 3D models. They search for the elements of the 3D model databases that resemble the query model. In structural bioinformatics, the query model is a protein tertiary/quaternary structure and the objective is to find similarly shaped molecules in the Protein Data Bank. With the ever-growing size of the PDB, a direct atomic coordinate comparison with all its members is impractical. To overcome this problem, the shape of the molecules can be encoded by fixed-length feature vectors. The distance of a protein to the entire PDB can be measured in this low-dimensional domain in linear time. The state-of-the-art approaches utilize Zernike-Canterakis moments for the shape encoding and supply the retrieval process with geometric data of the input structures. The BioZernike descriptors are a standard utility of the PDB since 2020. However, when trying to calculate the ZC moments locally, the issue of the deficiency of libraries readily available for use in custom programs (i.e., without relying on external binaries) is encountered, in particular programs written in Python. Here, a fast and well-documented Python implementation of the Pozo-Koehl algorithm is presented. In contrast to the more popular algorithm by Novotni and Klein, which is based on the voxelized volume, the PK algorithm produces ZC moments directly from the triangular surface meshes of 3D models. In particular, it can accept the molecular surfaces of proteins as its input. In the presented PK-Zernike library, owing to Numba's just-in-time compilation, a mesh with 50,000 facets is processed by a single thread in a second at the moment order 20. Since this is the first time the PK algorithm is used in structural bioinformatics, it is employed in a novel, simple, but efficient protein structure retrieval pipeline. The elimination of the outlying chain fragments via a fast PCA-based subroutine improves the discrimination ability, allowing for this pipeline to achieve an 0.961 area under the ROC curve in the BioZernike validation suite (0.997 for the assemblies). The correlation between the results of the proposed approach and of the 3D Surfer program attains values up to 0.99.

Improved Assessment of Globularity of Protein Structures and the Ellipsoid Profile of the Biological Assemblies from the PDB

In this paper, we present an update to the ellipsoid profile algorithm (EP), a simple technique for the measurement of the globularity of protein structures without the calculation of molecular surfaces. The globularity property is understood in this context as the ability of the molecule to fill a minimum volume enclosing ellipsoid (MVEE) that approximates its assumed globular shape. The more of the interior of this ellipsoid is occupied by the atoms of the protein, the better are its globularity metrics. These metrics are derived from the comparison of the volume of the voxelized representation of the atoms and the volume of all voxels that can fit inside that ellipsoid (a uniform unit Å cube lattice). The so-called ellipsoid profile shows how the globularity changes with the distance from the center. Two of its values, the so-called ellipsoid indexes, are used to classify the structure as globular, semi-globular or non-globular. Here, we enhance the workflow of the EP algorithm via an improved outlier detection subroutine based on principal component analysis. It is capable of robust distinguishing between the dense parts of the molecules and, for example, disordered chain fragments fully exposed to the solvent. The PCA-based method replaces the current approach based on kernel density estimation. The improved EP algorithm was tested on 2124 representatives of domain superfamilies from SCOP 2.08. The second part of this work is dedicated to the survey of globularity of 3594 representatives of biological assemblies from molecules currently deposited in the PDB and analyzed by the 3DComplex database (monomers and complexes up to 60 chains).

Symmetrization in the Calculation Pipeline of Gauss Function-Based Modeling of Hydrophobicity in Protein Structures

In this paper, we show, discuss, and compare the effects of symmetrization in two calculation subroutines of the Fuzzy Oil Drop model, a coarse-grained model of density of hydrophobicity in proteins. In the FOD model, an input structure is enclosed in an axis-aligned ellipsoid called a drop. Two profiles of hydrophobicity are then calculated for its residues: theoretical (based on the 3D Gauss function) and observed (based on pairwise hydrophobic interactions). Condition of the hydrophobic core is revealed by comparing those profiles through relative entropy, while analysis of their local differences allows, in particular, determination of the starting location for the search for protein–protein and protein–ligand interaction areas. Here, we improve the baseline workflow of the FOD model by introducing symmetry to the hydrophobicity profile comparison and ellipsoid bounding procedures. In the first modification (FOD–JS), Kullback–Leibler divergence is enhanced with its Jensen–Shannon variant. In the second modification (FOD-PCA), the molecule is optimally aligned with the axes of the coordinate system via principal component analysis, and the size of its drop is determined by the standard deviation of all its effective atoms, making it less susceptible to structural outliers. Tests on several molecules with various shapes and functions confirm that the proposed modifications improve the accuracy, robustness, speed, and usability of Gauss function-based modeling of the density of hydrophobicity in protein structures.

Wilavidin* - a novel member of the avidin family that forms unique biotin-binding hexamers

Nature's optimization of protein functions is a highly intricate evolutionary process. In addition to optimal tertiary folding, the intramolecular recognition among the monomers that generate higher-order quaternary arrangements is driven by stabilizing interactions that have a pivotal role for ideal activity. Homotetrameric avidin and streptavidin are regularly utilized in many applications, whereby their ultra-high affinity toward biotin is dependent on their quaternary arrangements. In recent years, a new subfamily of avidins was discovered that comprises homodimers rather than tetramers, in which the high affinity toward biotin is maintained. Intriguingly, several of the respective dimers have been shown to assemble into higher-order cylindrical hexamers or octamers that dissociate into dimers upon biotin binding. Here, we present wilavidin, a newly discovered member of the dimeric subfamily, forming hexamers in the apo form, which are uniquely maintained upon biotin binding with six high-affinity binding sites. Removal of the short C-terminal segment of wilavidin resulted in the presence of the dimer only, thus emphasizing the role of this segment in stabilizing the hexamer. Utilization of a hexavalent biotin-binding form of avidin would be beneficial for expanding the biotechnological toolbox. Additionally, this unique family of dimeric avidins and their propensity to oligomerize to hexamers or octamers can serve as a basis for protein oligomerization and intermonomeric recognition as well as cumulative interactions that determine molecular assemblies.

The Bacillus subtilis dCMP deaminase ComEB acts as a dynamic polar localization factor for ComGA within the competence machinery

Bacillus subtilis can import DNA from the environment by an uptake machinery that localizes to a single cell pole. We investigated the roles of ComEB and of the ATPase ComGA during the state of competence. We show that ComEB plays an important role during competence, possibly because it is necessary for the recruitment of GomGA to the cell pole. ComEB localizes to the cell poles even upon expression during exponential phase, indicating that it can serve as polar marker. ComEB is also a deoxycytidylate monophosphate (dCMP) deaminase, for the function of which a conserved cysteine residue is important. However, cysteine-mutant ComEB is still capable of natural transformation, while a comEB deletion strain is highly impaired in competence, indicating that ComEB confers two independent functions. Single-molecule tracking (SMT) reveals that both proteins exchange at the cell poles between bound and unbound in a time scale of a few milliseconds, but turnover of ComGA increases during DNA uptake, whereas the mobility of ComEB is not affected. Our data reveal a highly dynamic role of ComGA during DNA uptake and an unusual role for ComEB as a mediator of polar localization, localizing by diffusion-capture on an extremely rapid time scale and functioning as a moonlighting enzyme.

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.

Crystal structures of the complex of a kallikrein inhibitor from Bauhinia bauhinioides with trypsin and modeling of kallikrein complexes

Structures of a recombinant Kunitz-type serine protease inhibitor from Bauhinia bauhinioides (BbKI) complexed with bovine trypsin were determined in two crystal forms. The crystal structure with the L55R mutant of BbKI was determined in space group P64 at 1.94 Å resolution and that with native BbKI in the monoclinic space group P21 at 3.95 Å resolution. The asymmetric unit of the latter crystals contained 44 independent complexes, thus representing one of the largest numbers of independent objects deposited in the Protein Data Bank. Additionally, the structure of the complex with native BbKI was determined at 2.0 Å resolution from P64 crystals isomorphous to those of the mutant. Since BbKI has previously been found to be a potent inhibitor of the trypsin-like plasma kallikrein, it was also tested against several tissue kallikreins. It was found that BbKI is a potent inhibitor of human tissue kallikrein 4 (KLK4) and the chymotrypsin-like human tissue kallikrein 7 (KLK7). Structures of BbKI complexed with the catalytic domain of human plasma kallikrein were modeled, as well as those with KLK4 and KLK7, and the structures were analyzed in order to identify the interactions that are responsible for inhibitory potency.

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.

Preventive DNA repair by sanitizing the cellular (deoxy)nucleoside triphosphate pool

The occurrence of modified bases in DNA is attributed to some major factors: incorporation of altered nucleotide building blocks and chemical reactions or radiation effects on bases within the DNA structure. Several enzyme families are involved in preventing the incorporation of noncanonical bases playing a 'sanitizing' role. The catalytic mechanism of action of these enzymes has been revealed for a number of representatives in clear structural and kinetic detail. In this review, we focus in detail on those examples where clear evidence has been produced using high-resolution structural studies. Comparing the protein fold and architecture of the enzyme active sites, two main classes of sanitizing deoxyribonucleoside triphosphate pyrophosphatases can be assigned that are distinguished by the site of nucleophilic attack. In enzymes associated with attack at the α-phosphorus, it is shown that coordination of the γ-phosphate group is also ensured by multiple interactions. By contrast, enzymes catalyzing attack at the β-phosphorus atom mainly coordinate the α- and the β-phosphate only. Characteristic differences are also observed with respect to the role of the metal ion cofactor (Mg(2+) ) and the coordination of nucleophilic water. Using different catalytic mechanisms embedded in different protein folds, these enzymes present a clear example of convergent evolution.

Catalytic mechanism of α-phosphate attack in dUTPase is revealed by X-ray crystallographic snapshots of distinct intermediates, 31P-NMR spectroscopy and reaction path modelling

Enzymatic synthesis and hydrolysis of nucleoside phosphate compounds play a key role in various biological pathways, like signal transduction, DNA synthesis and metabolism. Although these processes have been studied extensively, numerous key issues regarding the chemical pathway and atomic movements remain open for many enzymatic reactions. Here, using the Mason-Pfizer monkey retrovirus dUTPase, we study the dUTPase-catalyzed hydrolysis of dUTP, an incorrect DNA building block, to elaborate the mechanistic details at high resolution. Combining mass spectrometry analysis of the dUTPase-catalyzed reaction carried out in and quantum mechanics/molecular mechanics (QM/MM) simulation, we show that the nucleophilic attack occurs at the α-phosphate site. Phosphorus-31 NMR spectroscopy ((31)P-NMR) analysis confirms the site of attack and shows the capability of dUTPase to cleave the dUTP analogue α,β-imido-dUTP, containing the imido linkage usually regarded to be non-hydrolyzable. We present numerous X-ray crystal structures of distinct dUTPase and nucleoside phosphate complexes, which report on the progress of the chemical reaction along the reaction coordinate. The presently used combination of diverse structural methods reveals details of the nucleophilic attack and identifies a novel enzyme-product complex structure.