Structure–function discrepancy inClostridium botulinumC3 toxin for its rational prioritization as a subunit vaccine

Low SARS-CoV-2 viral load among vaccinated individuals infected with Delta B.1.617.2 and Omicron BA.1.1.529 but not with Omicron BA.1.1 and BA.2 variants

The rapid spread of SARS-CoV-2 variants in the global population is indicative of the development of selective advantages in emerging virus strains. Here, we performed a case-control investigation of the clinical and demographic characteristics, clinical history, and virological markers to predict disease progression in hospitalized adults for COVID-19 between December 2021 and January 2022 in Chennai, India. COVID-19 diagnosis was made by a commercial TaqPath COVID-19 RT-PCR, and WGS was performed with the Ion Torrent Next Generation Sequencing System. High-quality (<5% of N) complete sequences of 73 Omicron B.1.1.529 variants were randomly selected for phylogenetic analysis. SARS-CoV-2 viral load, number of comorbidities, and severe disease presentation were independently associated with a shorter time-to-death. Strikingly, this was observed among individuals infected with Omicron BA.2 but not among those with the BA.1.1.529, BA.1.1, or the Delta B.1.617.2 variants. Phylogenetic analysis revealed severe cases predominantly clustering under the BA.2 lineage. Sequence analyses showed 30 mutation sites in BA.1.1.529 and 33 in BA.1.1. The mutations unique to BA.2 were T19I, L24S, P25del, P26del, A27S, V213G, T376A, D405N and R408S. Low SARS-CoV-2 viral load among vaccinated individuals infected with Delta B.1.617.2 and the Omicron BA.1.1.529 variant but not with Omicron BA.1.1 or BA.2 suggests that the newer strains are largely immune escape variants. The number of vaccine doses received was independently associated with increased odds of developing asymptomatic disease or recovery. We propose that the novel mutations reported herein could likely bear a significant impact on the clinical characteristics, disease progression, and epidemiological aspects of COVID-19. Surging rates of mutations and the emergence of eclectic variants of SARS-CoV-2 appear to impact disease dynamics.

Identification of genotypic variants and its proteomic mutations of Brazilian SARS-CoV-2 isolates

The second wave of COVID-19 caused by severe acute respiratory syndrome virus (SARS-CoV-2) is rapidly spreading over the world. Mechanisms behind the flee from current antivirals are still unclear due to the continuous occurrence of SARS-CoV-2 genetic variants. Brazil is the world's second-most COVID-19 affected country. In the present study, we identified the genomic and proteomic variants of Brazilian SARS-CoV-2 isolates. We identified 16 different genotypic variants were found among the 27 isolates. The genotypes of three isolates such as Bra/1236/2021 (G15), Bra/MASP2C844R2/2020 (G11), and Bra/RJ-DCVN5/2020 (G9) have a unique mutant in NSP4 (S184N), 2′O-Mutase (R216N), membrane protein (A2V) and Envelope protein (V5A). A mutation in RdRp of SARS-CoV-2, particularly the change of Pro-to Leu-at 323 resulted in the stabilization of the structure in BRA/CD1739-P4/2020. NSP4, NSP5 protein mutants are more virulent in genotype 15 and 16. A fast protein folding rate changes the structural stability and leads to escape for current antivirals. Thus, our findings help researchers to develop the best potent antivirals based on the new mutant of Brazilian isolates.

A hijack mechanism of Indian SARS-CoV-2 isolates for relapsing contemporary antiviral therapeutics

Coronavirus disease (COVID-19) rapidly expands to a global pandemic and its impact on public health varies from country to country. It is caused by a new virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It is imperative for relapsing current antiviral therapeutics owing to randomized genetic drift in global SARS-CoV-2 isolates. A molecular mechanism behind the emerging genomic variants is not yet understood for the prioritization of selective antivirals. The present computational study was aimed to repurpose existing antivirals for Indian SARS-CoV-2 isolates by uncovering a hijack mechanism based on structural and functional characteristics of protein variants. Forty-one protein mutations were identified in 12 Indian SARS-CoV-2 isolates by analysis of genome variations across 460 genome sequences obtained from 30 geographic sites in India. Two unique mutations such as W6152R and N5928H found in exonuclease of Surat (GBRC275b) and Gandhinagar (GBRC239) isolates. We report for the first time the impact of folding rate on stabilizing/retaining a sequence-structure-function-virulence link of emerging protein variants leading to accommodate hijack ability from current antivirals. Binding affinity analysis revealed the effect of point mutations on virus infectivity and the drug-escaping efficiency of Indian isolates. Emodin and artinemol suggested herein as repurposable antivirals for the treatment of COVID-19 patients infected with Indian isolates. Our study concludes that a protein folding rate is a key structural and evolutionary determinant to enhance the receptor-binding specificity and ensure hijack ability from the prevalent antiviral therapeutics.

Deciphering Molecular Virulence Mechanism of Mycobacterium tuberculosis Dop isopeptidase Based on Its Sequence-Structure-Function Linkage

The pupylation pathway marks proteins for prokaryotic ubiquitin-like protein (Pup)-proteasomal degradation and survival strategy of mycobacteria inside of the host macrophages. Deamidase of Pup (Dop) plays a central role in the pupylation pathway. It is still a matter of investigation to know the function of Dop in virulence of mycobacterial lineage. Hence, the present study was intended to describe the sequence-structure-function-virulence link of Dop for understanding the molecular virulence mechanism of Mycobacterium tuberculosis H37Rv (Mtb). Phylogenetic analysis of this study indicated that Dop has extensively diverged across the proteasome-harboring bacteria. The functional part of Dop was converged across the pathogenic mycobacterial lineage. The genome-wide analysis pointed out that the pupylation gene locus was identical to each other, but its genome neighborhood differed from species to species. Molecular modeling and dynamic studies proved that the predicted structure of Mtb Dop was energetically stable and low conformational freedom. Moreover, evolutionary constraints in Mtb Dop were intensively analyzed for inferring its sequence-structure-function relationships for the full virulence of Mtb. It indicated that evolutionary optimization was extensively required to stabilize its local structural environment at the side chains of mutable residues. The sequence-structure-function-virulence link of Dop might have retained in Mtb by reordering hydrophobic and hydrogen bonding patterns in the local structural environment. Thus, the results of our study provide a quest to understand the molecular virulence and pathogenesis mechanisms of Mtb during the infection process.

Substrate-imprinted docking of Agrobacterium tumefaciens uronate dehydrogenase for increased substrate selectivity

Agrobacterium tumefaciens uronate dehydrogenase (AtuUdh) belongs to the short-chain dehydrogenase superfamily, specifically those acting on the CH-OH group of donor with NAD⁺ or NADP⁺ as acceptor. It is apparently required for the production of D-glucaric acid. AtuUdh-catalyzed reaction is reversible with dual substrate-specific activity (D-galacturonic acid and D-glucuronic acid) in nature. In our study, 34 mutants were pre-screened from 155 mutants generated from AtuUdh (wild-type) and selected 10 structurally stable mutants with increased substrate selectivity. The specificity, efficiency, and selectivity of these mutants for different substrates and cofactors were predicted from 121 docked models using a substrate-imprinted docking approach. Q14F, S36L, and S75T mutants have shown a high binding affinity to D-glucuronic acid and its substrate intermediates such as D-glucaro-1,4-lactone and D-glucaro-1,5-lactone. These mutants exhibited a low binding affinity to the substrate and cofactor required for D-galactaric acid. D34S, N112E and S165E mutants found to show a high selectivity of D-galacturonic acid and its substrate intermediates for D-galactaric acid production. Ser75, Ser165, and Arg174 are active residues playing an imperative role in the substrate selectivity and also contributed in the conjecture the mechanism of transition state stabilization catalyzed by AtuUdh mutants. The present approach was used to reveal the substrate binding mechanism of AtuUdh mutants for a better understanding of the structural basis for selectivity and function.

Structural constraints-based evaluation of immunogenic avirulent toxins from Clostridium botulinum C2 and C3 toxins as subunit vaccines

Clostridium botulinum (group-III) is an anaerobic bacterium producing C2 and C3 toxins in addition to botulinum neurotoxins in avian and mammalian cells. C2 and C3 toxins are members of bacterial ADP-ribosyltransferase superfamily, which modify the eukaryotic cell surface proteins by ADP-ribosylation reaction. Herein, the mutant proteins with lack of catalytic and pore forming function derived from C2 (C2I and C2II) and C3 toxins were computationally evaluated to understand their structure-function integrity. We have chosen many structural constraints including local structural environment, folding process, backbone conformation, conformational dynamic sub-space, NAD-binding specificity and antigenic determinants for screening of suitable avirulent toxins. A total of 20 avirulent mutants were identified out of 23 mutants, which were experimentally produced by site-directed mutagenesis. No changes in secondary structural elements in particular to α-helices and β-sheets and also in fold rate of all-β classes. Structural stability was maintained by reordered hydrophobic and hydrogen bonding patterns. Molecular dynamic studies suggested that coupled mutations may restrain the binding affinity to NAD(+) or protein substrate upon structural destabilization. Avirulent toxins of this study have stable energetic backbone conformation with a common blue print of folding process. Molecular docking studies revealed that avirulent mutants formed more favorable hydrogen bonding with the side-chain of amino acids near to conserved NAD-binding core, despite of restraining NAD-binding specificity. Thus, structural constraints in the avirulent toxins would determine their immunogenic nature for the prioritization of protein-based subunit vaccine/immunogens to avian and veterinary animals infected with C. botulinum.