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.

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.

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.

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.

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

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.

Characterization of domain-selective inhibitor binding in angiotensin-converting enzyme using a novel derivative of lisinopril.

Human ACE (angiotensin-converting enzyme) (EC 3.4.15.1) is an important drug target because of its role in the regulation of blood pressure via the renin-angiotensin-aldosterone system. Somatic ACE comprises two homologous domains, the differing substrate preferences of which present a new avenue for domain-selective inhibitor design. We have co-crystallized lisW-S, a C-domain-selective derivative of the drug lisinopril, with human testis ACE and determined a structure using X-ray crystallography to a resolution of 2.30 A (1 A=0.1 nm). In this structure, lisW-S is seen to have a similar binding mode to its parent compound lisinopril, but the P2' tryptophan moiety takes a different conformation to that seen in other inhibitors having a tryptophan residue in this position. We have examined further the domain-specific interactions of this inhibitor by mutating C-domain-specific active-site residues to their N domain equivalents, then assessing the effect of the mutation on inhibition by lisW-S using a fluorescence-based assay. Kinetics analysis shows a 258-fold domain-selectivity that is largely due to the co-operative effect of C-domain-specific residues in the S2' subsite. The high affinity and selectivity of this inhibitor make it a good lead candidate for cardiovascular drug development.

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.

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 mutant chaperonin with rearranged inter-ring electrostatic contacts and temperature-sensitive dissociation.

The chaperonin GroEL assists protein folding through ATP-dependent, cooperative movements that alternately create folding chambers in its two rings. The substitution E461K at the interface between these two rings causes temperature-sensitive, defective protein folding in Escherichia coli. To understand the molecular defect, we have examined the mutant chaperonin by cryo-EM. The normal out-of-register alignment of contacts between subunits of opposing wild-type rings is changed in E461K to an in-register one. This is associated with loss of cooperativity in ATP binding and hydrolysis. Consistent with the loss of negative cooperativity between rings, the cochaperonin GroES binds simultaneously to both E461K rings. These GroES-bound structures were unstable at higher temperature, dissociating into complexes of single E461K rings associated with GroES. Lacking the allosteric signal from the opposite ring, these complexes cannot release their GroES and become trapped, dead-end states.

The cyanide hydratase from Neurospora crassa forms a helix which has a dimeric repeat.

The fungal cyanide hydratases form a functionally specialized subset of the nitrilases which catalyze the hydrolysis of cyanide to formamide with high specificity. These hold great promise for the bioremediation of cyanide wastes. The low resolution (3.0 nm) three-dimensional reconstruction of negatively stained recombinant cyanide hydratase fibers from the saprophytic fungus Neurospora crassa by iterative helical real space reconstruction reveals that enzyme fibers display left-handed D(1) S(5.4) symmetry with a helical rise of 1.36 nm. This arrangement differs from previously characterized microbial nitrilases which demonstrate a structure built along similar principles but with a reduced helical twist. The cyanide hydratase assembly is stabilized by two dyadic interactions between dimers across the one-start helical groove. Docking of a homology-derived atomic model into the experimentally determined negative stain envelope suggests the location of charged residues which may form salt bridges and stabilize the helix.

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.

Probing the basis of domain-dependent inhibition using novel ketone inhibitors of Angiotensin-converting enzyme.

Human angiotensin-converting enzyme (ACE) has two homologous domains, the N and C domains, with differing substrate preferences. X-ray crystal structures of the C and N domains complexed with various inhibitors have allowed identification of active site residues that might be important for the molecular basis of this selectivity. However, it is unclear to what extent the different residues contribute to substrate domain selectivity. Here, cocrystal structures of human testis ACE, equivalent to the C domain, have been determined with two novel C domain-selective ketomethylene inhibitors, (5 S)-5-[( N-benzoyl)amino]-4-oxo-6-phenylhexanoyl- l-tryptophan (kAW) and (5 S)-5-[( N-benzoyl)amino]-4-oxo-6-phenylhexanoyl- l-phenylalanine (kAF). The ketone groups of both inhibitors bind to the zinc ion as a hydrated geminal diolate, demonstrating the ability of the active site to catalyze the formation of the transition state. Moreover, active site residues involved in inhibitor binding have been mutated to their N domain counterparts, and the effect of the mutations on inhibitor binding has been determined. The C domain selectivity of these inhibitors was found to result from interactions between bulky hydrophobic side chain moieties and C domain-specific residues F391, V518, E376, and V380 (numbering of testis ACE). Mutation of these residues decreased the affinity for the inhibitors 4-20-fold. T282, V379, E403, D453, and S516 did not contribute individually to C domain-selective inhibitor binding. Further domain-selective inhibitor design should focus on increasing both the affinity and selectivity of the side chain moieties.

Genome mining of cyanide-degrading nitrilases from filamentous fungi.

A variety of fungal species are known to degrade cyanide through the action of cyanide hydratases, a specialized subset of nitrilases which hydrolyze cyanide to formamide. In this paper, we report on two previously unknown and uncharacterized cyanide hydratases from Neurospora crassa and Aspergillus nidulans. Recombinant forms of four cyanide hydratases from N. crassa, A. nidulans, Gibberella zeae, and Gloeocercospora sorghi were prepared after their genes were cloned with N-terminal hexahistidine purification tags, expressed in Escherichia coli, and purified using immobilized metal affinity chromatography. These enzymes were compared according to their relative specific activity, pH activity profiles, thermal stability, and ability to remediate cyanide contaminated waste water from silver and copper electroplating baths. Although all four were similar, the N. crassa cyanide hydratase (CHT) has the greatest thermal stability and widest pH range of >50% activity. N. crassa also demonstrated the highest rate of cyanide degradation in the presence of both heavy metals. The CHT of A. nidulans has the highest reaction rate of the four fungal nitrilases evaluated in this work. These data will help determine optimization procedures for the possible use of these enzymes in the bioremediation of cyanide-containing waste. Similar to known plant pathogenic fungi, both N. crassa and A. nidulans were induced to express CHT by growth in the presence of KCN.

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.

Post-translational cleavage of recombinantly expressed nitrilase from Rhodococcus rhodochrous J1 yields a stable, active helical form.

Nitrilases convert nitriles to the corresponding carboxylic acids and ammonia. The nitrilase from Rhodococcus rhodochrous J1 is known to be inactive as a dimer but to become active on oligomerization. The recombinant enzyme undergoes post-translational cleavage at approximately residue 327, resulting in the formation of active, helical homo-oligomers. Determining the 3D structure of these helices using electron microscopy, followed by fitting the stain envelope with a model based on homology with other members of the nitrilase superfamily, enables the interacting surfaces to be identified. This also suggests that the reason for formation of the helices is related to the removal of steric hindrance arising from the 39 C-terminal amino acids from the wild-type protein. The helical form can be generated by expressing only residues 1-327.

Three-dimensional structure of a type III glutamine synthetase by single-particle reconstruction.

GlnN, the type III glutamine synthetase (GSIII) from the medically important, anaerobic, opportunistic pathogen Bacteroides fragilis, has 82.8 kDa subunits that share only 9% sequence identity with the type I glutamine synthetases (GSI), the only family for which a structure is known. Active GlnN was found predominantly in a single peak that eluted from a calibrated gel-filtration chromatography column at a position equaivalent to 0.86(+/-0.08) MDa. Negative-stain electron microscopy enabled the identification of double-ringed particles and single hexameric rings ("pinwheels") resulting from partial staining. A 2D average of these pinwheels showed marked similarity to the corresponding structures found in preparations of GSI, except that the arms of the subunits were 40% longer. Reconstructions from particles embedded in vitreous ice showed that GlnN has a double-ringed, dodecameric structure with a 6-fold dihedral space group (D6) symmetry and dimensions of 17.0 nm parallel with the 6-fold axis and 18.3 nm parallel with the 2-fold axes. The structures, combined with a sequence alignment based on structural principles, showed how many aspects of the structure of GSI, and most notably the alpha/beta barrel fold active site were preserved. There was evidence for the presence of this structure in the reconstructed volume, thus, identifying the indentations between the pinwheel spokes as putative active sites and suggesting conservation of the overall molecular geometry found in GSI despite their low level of global homology. Furthermore, docking of GSI into the reconstruction left sufficient plausibly located unoccupied density to account for the additional residues in GSIII, thus validating the structure.