Engineering enhanced thermostability into the Geobacillus pallidus nitrile hydratase.

Nitrile hydratases (NHases) are important biocatalysts for the enzymatic conversion of nitriles to industrially-important amides such as acrylamide and nicotinamide. Although thermostability in this enzyme class is generally low, there is not sufficient understanding of its basis for rational enzyme design. The gene expressing the Co-type NHase from the moderate thermophile, Geobacillus pallidus RAPc8 (NRRL B-59396), was subjected to random mutagenesis. Four mutants were selected that were 3 to 15-fold more thermostable than the wild-type NHase, resulting in a 3.4-7.6 â€‹kJ/mol increase in the activation energy of thermal inactivation at 63 â€‹°C. High resolution X-ray crystal structures (1.15-1.80 â€‹Å) were obtained of the wild-type and four mutant enzymes. Mutant 9E, with a resolution of 1.15 â€‹Å, is the highest resolution crystal structure obtained for a nitrile hydratase to date. Structural comparisons between the wild-type and mutant enzymes illustrated the importance of salt bridges and hydrogen bonds in enhancing NHase thermostability. These additional interactions variously improved thermostability by increased intra- and inter-subunit interactions, preventing cooperative unfolding of α-helices and stabilising loop regions. Some hydrogen bonds were mediated via a water molecule, specifically highlighting the significance of structured water molecules in protein thermostability. Although knowledge of the mutant structures makes it possible to rationalize their behaviour, it would have been challenging to predict in advance that these mutants would be stabilising.

Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function.

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Selection analysis identifies unusual clustered mutational changes in Omicron lineage BA.1 that likely impact Spike function.

Among the 30 non-synonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (i) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (ii) interactions of Spike with ACE2 receptors, and (iii) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any genomes within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron over all previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Biochemical and structural insights into the cytochrome P450 reductase from Candida tropicalis.

Cytochrome P450 reductases (CPRs) are diflavin oxidoreductases that supply electrons to type II cytochrome P450 monooxygenases (CYPs). In addition, it can also reduce other proteins and molecules, including cytochrome c, ferricyanide, and different drugs. Although various CPRs have been functionally and structurally characterized, the overall mechanism and its interaction with different redox acceptors remain elusive. One of the main problems regarding electron transfer between CPRs and CYPs is the so-called "uncoupling", whereby NAD(P)H derived electrons are lost due to the reduced intermediates' (FAD and FMN of CPR) interaction with molecular oxygen. Additionally, the decay of the iron-oxygen complex of the CYP can also contribute to loss of reducing equivalents during an unproductive reaction cycle. This phenomenon generates reactive oxygen species (ROS), leading to an inefficient reaction. Here, we present the study of the CPR from Candida tropicalis (CtCPR) lacking the hydrophobic N-terminal part (Δ2-22). The enzyme supports the reduction of cytochrome c and ferricyanide, with an estimated 30% uncoupling during the reactions with cytochrome c. The ROS produced was not influenced by different physicochemical conditions (ionic strength, pH, temperature). The X-ray structures of the enzyme were solved with and without its cofactor, NADPH. Both CtCPR structures exhibited the closed conformation. Comparison with the different solved structures revealed an intricate ionic network responsible for the regulation of the open/closed movement of CtCPR.

Substrate specificity of plant nitrilase complexes is affected by their helical twist.

Nitrilases are oligomeric, helix-forming enzymes from plants, fungi and bacteria that are involved in the metabolism of various natural and artificial nitriles. These biotechnologically important enzymes are often specific for certain substrates, but directed attempts at modifying their substrate specificities by exchanging binding pocket residues have been largely unsuccessful. Thus, the basis for their selectivity is still unknown. Here we show, based on work with two highly similar nitrilases from the plant Capsella rubella, that modifying nitrilase helical twist, either by exchanging an interface residue or by imposing a different twist, without altering any binding pocket residues, changes substrate preference. We reveal that helical twist and substrate size correlate and when binding pocket residues are exchanged between two nitrilases that show the same twist but different specificities, their specificities change. Based on these findings we propose that helical twist influences the overall size of the binding pocket.

Residue Y70 of the Nitrilase Cyanide Dihydratase from Bacillus pumilus Is Critical for Formation and Activity of the Spiral Oligomer.

Nitrilases pose attractive alternatives to the chemical hydrolysis of nitrile compounds. The activity of bacterial nitrilases towards substrate is intimately tied to the formation of large spiral-shaped oligomers. In the nitrilase CynD (cyanide dihydratase) from Bacillus pumilus, mutations in a predicted oligomeric surface region altered its oligomerization and reduced its activity. One mutant, CynD Y70C, retained uniform oligomer formation however it was inactive, unlike all other inactive mutants throughout that region all of which significantly perturbed oligomer formation. It was hypothesized that Y70 is playing an additional role necessary for CynD activity beyond influencing oligomerization. Here, we performed saturation mutagenesis at residue 70 and demonstrated that only tyrosine or phenylalanine is permissible for CynD activity. Furthermore, we show that other residues at this position are not only inactive, but have altered or disrupted oligomer conformations. These results suggest that Y70's essential role in activity is independent of its role in the formation of the spiral oligomer.

Bacillus pumilus Cyanide Dihydratase Mutants with Higher Catalytic Activity.

Cyanide degrading nitrilases are noted for their potential to detoxify industrial wastewater contaminated with cyanide. However, such application would benefit from an improvement to characteristics such as their catalytic activity and stability. Following error-prone PCR for random mutagenesis, several cyanide dihydratase mutants from Bacillus pumilus were isolated based on improved catalysis. Four point mutations, K93R, D172N, A202T, and E327K were characterized and their effects on kinetics, thermostability and pH tolerance were studied. K93R and D172N increased the enzyme's thermostability whereas E327K mutation had a less pronounced effect on stability. The D172N mutation also increased the affinity of the enzyme for its substrate at pH 7.7 but lowered its k cat. However, the A202T mutation, located in the dimerization or the A surface, destabilized the protein and abolished its activity. No significant effect on activity at alkaline pH was observed for any of the purified mutants. These mutations help confirm the model of CynD and are discussed in the context of the protein-protein interfaces leading to the protein quaternary structure.

Probing C-terminal interactions of the Pseudomonas stutzeri cyanide-degrading CynD protein.

The cyanide dihydratases from Bacillus pumilus and Pseudomonas stutzeri share high amino acid sequence similarity throughout except for their highly divergent C-termini. However, deletion or exchange of the C-termini had different effects upon each enzyme. Here we extended previous studies and investigated how the C-terminus affects the activity and stability of three nitrilases, the cyanide dihydratases from B. pumilus (CynDpum) and P. stutzeri (CynDstut) and the cyanide hydratase from Neurospora crassa. Enzymes in which the C-terminal residues were deleted decreased in both activity and thermostability with increasing deletion lengths. However, CynDstut was more sensitive to such truncation than the other two enzymes. A domain of the P. stutzeri CynDstut C-terminus not found in the other enzymes, 306GERDST311, was shown to be necessary for functionality and explains the inactivity of the previously described CynDstut-pum hybrid. This suggests that the B. pumilus C-terminus, which lacks this motif, may have specific interactions elsewhere in the protein, preventing it from acting in trans on a heterologous CynD protein. We identify the dimerization interface A-surface region 195-206 (A2) from CynDpum as this interaction site. However, this A2 region did not rescue activity in C-terminally truncated CynDstutΔ302 or enhance the activity of full-length CynDstut and therefore does not act as a general stability motif.

Probing an Interfacial Surface in the Cyanide Dihydratase from Bacillus pumilus, A Spiral Forming Nitrilase.

Nitrilases are of significant interest both due to their potential for industrial production of valuable products as well as degradation of hazardous nitrile-containing wastes. All known functional members of the nitrilase superfamily have an underlying dimer structure. The true nitrilases expand upon this basic dimer and form large spiral or helical homo-oligomers. The formation of this larger structure is linked to both the activity and substrate specificity of these nitrilases. The sequences of the spiral nitrilases differ from the non-spiral forming homologs by the presence of two insertion regions. Homology modeling suggests that these regions are responsible for associating the nitrilase dimers into the oligomer. Here we used cysteine scanning across these two regions, in the spiral forming nitrilase cyanide dihydratase from Bacillus pumilus (CynD), to identify residues altering the oligomeric state or activity of the nitrilase. Several mutations were found to cause changes to the size of the oligomer as well as reduction in activity. Additionally one mutation, R67C, caused a partial defect in oligomerization with the accumulation of smaller oligomer variants. These results support the hypothesis that these insertion regions contribute to the unique quaternary structure of the spiral microbial nitrilases.

Synthetic Hemozoin (ß-Hematin) Crystals Nucleate at the Surface of Neutral Lipid Droplets that Control Their Sizes.

Emulsions of monopalmitoylglycerol (MPG) and of a neutral lipid blend (NLB), consisting of MPG, monostearoylglycerol, dipalmitoylglycerol, dioleoylglycerol and dilineoylglycerol (4:2:1:1:1), the composition associated with hemozoin from the malaria parasite Plasmodium falciparum, have been used to mediate the formation of ß-hematin microcrystals. Transmission electron microscopy (TEM), electron diffraction and electron spectroscopic imaging/electron energy loss spectroscopy (ESI/EELS) have been used to characterize both the lipid emulsion and ß-hematin crystals. The latter have been compared with ß-hematin formed at a pentanol/aqueous interface and with hemozoin both within P. falciparum parasites and extracted from the parasites. When lipid and ferriprotoporphyrin IX solutions in 1:9 v/v acetone/methanol were thoroughly pre-mixed either using an extruder or ultrasound, ß-hematin crystals were found formed in intimate association with the lipid droplets. These crystals resembled hemozoin crystals, with prominent {100} faces. Lattice fringes in TEM indicated that these faces made contact with the lipid surface. The average length of these crystals was 0.62 times the average diameter of NLB droplets and their size distributions were statistically equivalent after 10 min incubation, suggesting that the lipid droplets also controlled the sizes of the crystals. This most closely resembles hemozoin formation in the helminth worm Schistosoma mansoni, while in P. falciparum, crystal formation appears to be associated with the much more gently curved digestive vacuole membrane which apparently leads to formation of much larger hemozoin crystals, similar to those formed at the flat pentanol-water interface.

The mechanism of the amidases: mutating the glutamate adjacent to the catalytic triad inactivates the enzyme due to substrate mispositioning.

All known nitrilase superfamily amidase and carbamoylase structures have an additional glutamate that is hydrogen bonded to the catalytic lysine in addition to the Glu, Lys, Cys "catalytic triad." In the amidase from Geobacillus pallidus, mutating this glutamate (Glu-142) to a leucine or aspartate renders the enzyme inactive. X-ray crystal structure determination shows that the structural integrity of the enzyme is maintained despite the mutation with the catalytic cysteine (Cys-166), lysine (Lys-134), and glutamate (Glu-59) in positions similar to those of the wild-type enzyme. In the case of the E142L mutant, a chloride ion is located in the position occupied by Glu-142 O(ε1) in the wild-type enzyme and interacts with the active site lysine. In the case of the E142D mutant, this site is occupied by Asp-142 O(δ1.) In neither case is an atom located at the position of Glu-142 O(ε2) in the wild-type enzyme. The active site cysteine of the E142L mutant was found to form a Michael adduct with acrylamide, which is a substrate of the wild-type enzyme, due to an interaction that places the double bond of the acrylamide rather than the amide carbonyl carbon adjacent to the active site cysteine. Our results demonstrate that in the wild-type active site a crucial role is played by the hydrogen bond between Glu-142 O(ε2) and the substrate amino group in positioning the substrate with the correct stereoelectronic alignment to enable the nucleophilic attack on the carbonyl carbon by the catalytic cysteine.

Structural insight into African horsesickness virus infection.

African horsesickness (AHS) is a devastating disease of horses. The disease is caused by the double-stranded RNA-containing African horsesickness virus (AHSV). Using electron cryomicroscopy and three-dimensional image reconstruction, we determined the architecture of an AHSV serotype 4 (AHSV-4) reference strain. The structure revealed triple-layered AHS virions enclosing the segmented genome and transcriptase complex. The innermost protein layer contains 120 copies of VP3, with the viral polymerase, capping enzyme, and helicase attached to the inner surface of the VP3 layer on the 5-fold axis, surrounded by double-stranded RNA. VP7 trimers form a second, T=13 layer on top of VP3. Comparative analyses of the structures of bluetongue virus and AHSV-4 confirmed that VP5 trimers form globular domains and VP2 trimers form triskelions, on the virion surface. We also identified an AHSV-7 strain with a truncated VP2 protein (AHSV-7 tVP2) which outgrows AHSV-4 in culture. Comparison of AHSV-7 tVP2 to bluetongue virus and AHSV-4 allowed mapping of two domains in AHSV-4 VP2, and one in bluetongue virus VP2, that are important in infection. We also revealed a protein plugging the 5-fold vertices in AHSV-4. These results shed light on virus-host interactions in an economically important orbivirus to help the informed design of new vaccines.

Immobilisation and characterisation of biocatalytic co-factor recycling enzymes, glucose dehydrogenase and NADH oxidase, on aldehyde functional ReSyn™ polymer microspheres.

The use of enzymes in industrial applications is limited by their instability, cost and difficulty in their recovery and re-use. Immobilisation is a technique which has been shown to alleviate these limitations in biocatalysis. Here we describe the immobilisation of two biocatalytically relevant co-factor recycling enzymes, glucose dehydrogenase (GDH) and NADH oxidase (NOD) on aldehyde functional ReSyn™ polymer microspheres with varying functional group densities. The successful immobilisation of the enzymes on this new high capacity microsphere technology resulted in the maintenance of activity of ∼40% for GDH and a maximum of 15.4% for NOD. The microsphere variant with highest functional group density of ∼3500 µmol g⁻¹ displayed the highest specific activity for the immobilisation of both enzymes at 33.22 U mg⁻¹ and 6.75 U mg⁻¹ for GDH and NOD with respective loading capacities of 51% (0.51 mg mg⁻¹) and 129% (1.29 mg mg⁻¹). The immobilised GDH further displayed improved activity in the acidic pH range. Both enzymes displayed improved pH and thermal stability with the most pronounced thermal stability for GDH displayed on ReSyn™ A during temperature incubation at 65 °C with a 13.59 fold increase, and NOD with a 2.25-fold improvement at 45 °C on the same microsphere variant. An important finding is the suitability of the microspheres for stabilisation of the multimeric protein GDH.

Engineering pH-tolerant mutants of a cyanide dihydratase.

Cyanide dihydratase is an enzyme in the nitrilase family capable of transforming cyanide to formate and ammonia. This reaction has been exploited for the bioremediation of cyanide in wastewater streams, but extending the pH operating range of the enzyme would improve its utility. In this work, we describe mutants of Bacillus pumilus C1 cyanide dihydratase (CynD(pum)) with improved activity at higher pH. Error-prone PCR was used to construct a library of CynD(pum) mutants, and a high-throughput screening system was developed to screen the library for improved activity at pH 10. Two mutant alleles were identified that allowed cells to degrade cyanide in solutions at pH 10, whereas the wild-type was inactive above pH 9. The mutant alleles each encoded three different amino acid substitutions, but for one of those, a single change, E327G, accounted for the phenotype. The purified proteins containing multiple mutations were five times more active than the wild-type enzyme at pH 9, but all purified enzymes lost activity at pH 10. The mutation Q86R resulted in the formation of significantly longer fibers at low pH, and both E327G and Q86R contributed to the persistence of active oligomeric assemblies at pH 9. In addition, the mutant enzymes proved to be more thermostable than the wild type, suggesting improved physical stability rather than any change in chemistry accounts for their increased pH tolerance.

Symmetry-restrained flexible fitting for symmetric EM maps.

Many large biological macromolecules have inherent structural symmetry, being composed of a few distinct subunits, repeated in a symmetric array. These complexes are often not amenable to traditional high-resolution structural determination methods, but can be imaged in functionally relevant states using cryo-electron microscopy (cryo-EM). A number of methods for fitting atomic-scale structures into cryo-EM maps have been developed, including the molecular dynamics flexible fitting (MDFF) method. However, quality and resolution of the cryo-EM map are the major determinants of a method's success. In order to incorporate knowledge of structural symmetry into the fitting procedure, we developed the symmetry-restrained MDFF method. The new method adds to the cryo-EM map-derived potential further restraints on the allowed conformations of a complex during fitting, thereby improving the quality of the resultant structure. The benefit of using symmetry-based restraints during fitting, particularly for medium to low-resolution data, is demonstrated for three different systems.

Proteolysis of the type III glutamine synthetase from Bacteroides fragilis causes expedient crystal-packing rearrangements.

This work details the intentional modifications that led to the first structure of a type III glutamine synthetase enzyme (GSIII). This approach followed the serendipitous discovery of digestion caused by an extracellular protease from a contaminating bacterium, Pseudomonas fluorescens. The protease only cleaves the GSIII protein at a single site, leaving the oligomer intact but allowing the protein to crystallize in a different space group. This transition from space group P1 to space group C222(1) is accompanied by improved growth characteristics, more reproducible diffraction and enhanced mechanical stability. The crystallographic analyses presented here provide the structural basis of the altered molecular packing in the full-length and digested crystal forms and suggest modifications for future structural studies.

The role of a topologically conserved isoleucine in glutathione transferase structure, stability and function.

The common fold shared by members of the glutathione-transferase (GST) family has a topologically conserved isoleucine residue at the N-terminus of helix 3 which is involved in the packing of helix 3 against two beta-strands in domain 1. The role of the isoleucine residue in the structure, function and stability of GST was investigated by replacing the Ile71 residue in human GSTA1-1 by alanine or valine. The X-ray structures of the I71A and I71V mutants resolved at 1.75 and 2.51 A, respectively, revealed that the mutations do not alter the overall structure of the protein compared with the wild type. Urea-induced equilibrium unfolding studies using circular dichroism and tryptophan fluorescence suggest that the mutation of Ile71 to alanine or valine reduces the stability of the protein. A functional assay with 1-chloro-2,4-dinitrobenzene shows that the mutation does not significantly alter the function of the protein relative to the wild type. Overall, the results suggest that conservation of the topologically conserved Ile71 maintains the structural stability of the protein but does not play a significant role in catalysis and substrate binding.

Structural and biochemical characterization of a nitrilase from the thermophilic bacterium, Geobacillus pallidus RAPc8.

Geobacillus pallidus RAPc8 (NRRL: B-59396) is a moderately thermophilic gram-positive bacterium, originally isolated from Australian lake sediment. The G. pallidus RAPc8 gene encoding an inducible nitrilase was located and cloned using degenerate primers coding for well-conserved nitrilase sequences, coupled with inverse PCR. The nitrilase open reading frame was cloned into an expression plasmid and the expressed recombinant enzyme purified and characterized. The protein had a monomer molecular weight of 35,790 Da, and the purified functional enzyme had an apparent molecular weight of approximately 600 kDa by size exclusion chromatography. Similar to several plant nitrilases and some bacterial nitrilases, the recombinant G. pallidus RAPc8 enzyme produced both acid and amide products from nitrile substrates. The ratios of acid to amide produced from the substrates we tested are significantly different to those reported for other enzymes, and this has implications for our understanding of the mechanism of the nitrilases which may assist with rational design of these enzymes. Electron microscopy and image classification showed complexes having crescent-like, "c-shaped", circular and "figure-8" shapes. Protein models suggested that the various complexes were composed of 6, 8, 10 and 20 subunits, respectively.

A modeled structure of an aptamer-gp120 complex provides insight into the mechanism of HIV-1 neutralization.

The HIV-1 envelope glycoprotein, gp120, is a key target for a class of drugs called entry inhibitors. Here we used molecular modeling to construct a three-dimensional model of an anti-gp120 RNA aptamer, B40t77, alone and in complex with gp120. An initial model of B40t77 was built from the predicted secondary structure and then subjected to a combination of energy minimization and molecular dynamics. To model the B40t77-gp120 complex, we docked the B40t77 predicted structure onto the CD4-induced epitope of the gp120 crystal structure. A series of gp120 point mutations in the predicted B40t77-gp120 interface were measured for their binding affinity for B40t77 by surface plasmon resonance. According to the model, of the 10 gp120 amino acids that showed a reduction in the level of binding when mutated to alanine, all of them are modeled as making direct contact with B40t77 as part of a hydrogen bonding network. Comparison by electron microscopy of the B40t77-gp120 complex with gp120 alone revealed that only the longest dimension of the complex significantly increased in length, in a manner consistent with the predicted model. Binding assays revealed that B40t77 can weaken the binding of gp120 to the monoclonal antibodies B6, B12, and 2G12, none of which have binding sites that overlap with B40t77, as well as strengthen the binding to the antibody 19b. Thus, B40t77 may induce distant conformational changes in gp120 that disrupt its association with host cells and may suggest a mechanism for aptamer neutralization of HIV-1.