Boophilin D1, a Kunitz type protease inhibitor, as a source of inhibitors for the ZIKA virus NS2B-NS3 protease.

Arboviruses are a global concern for a multitude of reasons, including their increased incidence and human mortality. Vectors associated with arboviruses include the mosquito Aedes sp., which is responsible for transmitting the Zika virus. Flaviviruses, like the Zika virus, present only one chymotrypsin-like serine protease (NS3) in their genome. Together with host enzymes, the NS2B co-factor NS3 protease complex are essential for the viral replication cycle by virus polyprotein processing. To search for Zika virus NS2B-NS3 protease (ZIKVPro) inhibitors, a phage display library was constructed using the Boophilin domain 1 (BoophD1), a thrombin inhibitor from the Kunitz family. A BoophilinD1 library mutated at positions P1-P4' was constructed, presenting a titer of 2.9x106 (cfu), and screened utilizing purified ZIKVPro. The results demonstrated at the P1-P4' positions the occurrence of 47% RALHA sequence (mut 12) and 11.8% RASWA sequence (mut14), SMRPT, or KALIP (wt) sequence. BoophD1-wt and mutants 12 and 14 were expressed and purified. The purified BoophD1 wt, mut 12 and 14, presented Ki values for ZIKVPro of 0.103, 0.116, and 0.101 µM, respectively. The BoophD1 mutant inhibitors inhibit the Dengue virus 2 protease (DENV2) with Ki values of 0.298, 0.271, and 0.379 µM, respectively. In conclusion, BoophD1 mut 12 and 14 selected for ZIKVPro demonstrated inhibitory activity like BoophD1-wt, suggesting that these are the strongest Zika inhibitors present in the BoophD1 mutated phage display library. Furthermore, BoophD1 mutants selected for ZIKVPro inhibit both Zika and Dengue 2 proteases making them potential pan-flavivirus inhibitors.

Modeling the role of charged residues in thermophilic proteins by rotamer and dynamic cross correlation analysis.

Discerning the determinants of protein thermostability is very important both from the theoretical and applied perspective. Different lines of evidence seem to indicate that a dynamical network of salt bridges/charged residues plays a fundamental role in the thermostability of enzymes. In this work, we applied measures of dynamic variance, like the Gini coefficients, Kullback-Leibler (KL) divergence and dynamic cross correlation (DCC) coefficients to compare the behavior of 3 pairs of homologous proteins from the thermophilic bacterium Thermus thermophilus and mesophilic Escherichia coli. Molecular dynamic (MD) simulations of these proteins were performed at 303 K and 363 K. In the characterization of their side chain rotamer distributions, the corresponding Gini coefficients and KL-divergence both revealed significant correlations with temperature. Similarly, a DCC analysis revealed a higher trend to de-correlate the movement of charged residues at higher temperatures in the thermophilic proteins, when compared with their mesophilic homologues. These results highlight the importance of dynamic electrostatic network interactions for the thermostability of enzymes.

A proteomics-MM/PBSA dual approach for the analysis of SARS-CoV-2 main protease substrate peptide specificity.

The main protease Mpro of SARS-CoV-2 is a well-studied major drug target. Additionally, it has been linked to this virus' pathogenicity, possibly through off-target effects. It is also an interesting diagnostic target. To obtain more data on possible substrates as well as to assess the enzyme's primary specificity a two-step approach was introduced. First, Terminal Amine Isobaric Labeling of Substrates (TAILS) was employed to identify novel Mpro cleavage sites in a mouse lung proteome library. In a second step, using a structural homology model, the MM/PBSA variant MM/GBSA (Molecular Mechanics Poisson-Boltzmann/Generalized Born Surface Area) free binding energy calculations were carried out to determine relevant interacting amino acids. As a result, 58 unique cleavage sites were detected, including six that displayed glutamine at the P1 position. Furthermore, modeling results indicated that Mpro has a far higher potential promiscuity towards substrates than expected. The combination of proteomics and MM/PBSA modeling analysis can thus be useful for elucidating the specificity of Mpro, and thus open novel perspectives for the development of future peptidomimetic drugs against COVID-19, as well as diagnostic tools.

Biochemical Screening of Potent Zika Virus Protease Inhibitors.

As the Zika virus protease is an essential and well-established target for the development of antiviral agents, we biochemically screened for inhibitors using a purified recombinantly expressed form of this enzyme. As a result, we were able to identify 10 new Zika virus protease inhibitors. These compounds are natural products and showed strong inhibition in the biochemical assays. Inhibitory constants values for the compounds ranged from 5 nM to 8 µM. Among the most potent inhibitors are flavonoids like irigenol hexa-acetate (Ki =0.28 µM), katacine (Ki =0.26 µM), theaflavin gallate (Ki =0.40 µM) and hematein (Ki =0.33 µM). Inhibitors from other groups of natural products include sennoside A (Ki =0.19 µM) and gossypol (Ki =0.70 µM). Several of the obtained compounds are known for their beneficial health effects and have acceptable pharmacokinetic characteristics. Thus, they could be of interest as lead compounds for the development of important and essential Zika antiviral drugs.

A proteomics-MM/PBSA dual approach for the analysis of SARS-CoV-2 main protease substrate peptide specificity

The main protease Mpro of SARS-CoV-2 is a well-studied major drug target. Additionally, it has been linked to this virus’ pathogenicity, possibly through off-target effects. It is also an interesting diagnostic target. To obtain more data on possible substrates as well as to assess the enzyme’s primary specificity a two-step approach was introduced. First, Terminal Amine Isobaric Labeling of Substrates (TAILS) was employed to identify novel Mpro cleavage sites in a mouse lung proteome library. In a second step, using a structural homology model, the MM/PBSA variant MM/GBSA (Molecular Mechanics Poisson-Boltzmann/Generalized Born Surface Area) free binding energy calculations were carried out to determine relevant interacting amino acids. As a result, 58 unique cleavage sites were detected, including six that displayed glutamine at the P1 position. Furthermore, modeling results indicated that Mpro has a far higher potential promiscuity towards substrates than expected. The combination of proteomics and MM/PBSA modeling analysis can thus be useful for elucidating the specificity of Mpro, and thus open novel perspectives for the development of future peptidomimetic drugs against COVID-19, as well as diagnostic tools.

Dynamic cross correlation analysis of Thermus thermophilus alkaline phosphatase and determinants of thermostability.

BACKGROUND: Understanding the determinants of protein thermostability is very important both from the theoretical and applied perspective. One emerging view in thermostable enzymes seems to indicate that a salt bridge/charged residue network plays a fundamental role in their thermostability. METHODS: The structure of alkaline phosphatase (AP) from Thermus thermophilus HB8 was solved by X-ray crystallography at 2.1 Å resolution. The obtained structure was further analyzed by molecular dynamics studies at different temperatures (303 K, 333 K and 363 K) and compared to homologous proteins from the cold-adapted organisms Shewanella sp. and Vibrio strain G15-21. To analyze differences in measures of dynamic variation, several data reduction techniques like principal component analysis (PCA), residue interaction network (RIN) analysis and rotamer analysis were used. Using hierarchical clustering, the obtained results were combined to determine residues showing high degree dynamical variations due to temperature jumps. Furthermore, dynamic cross correlation (DCC) analysis was carried out to characterize networks of charged residues. RESULTS: Top clustered residues showed a higher propensity for thermostabilizing mutations, indicating evolutionary pressure acting on thermophilic organisms. The description of rotamer distributions by Gini coefficients and Kullback-Leibler (KL) divergence both revealed significant correlations with temperature. DCC analysis revealed a significant trend to de-correlation of the movement of charged residues at higher temperatures. SIGNIFICANCE: The de-correlation of charged residues detected in Thermus thermophilus AP, highlights the importance of dynamic electrostatic network interactions for the thermostability of this enzyme.

Biochemical screening for SARS-CoV-2 main protease inhibitors.

The Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2) pandemic represents a global challenge. SARS-CoV-2's ability to replicate in host cells relies on the action of its non-structural proteins, like its main protease (Mpro). This cysteine protease acts by processing the viruses' precursor polyproteins. As proteases, together with polymerases, are main targets of antiviral drug design, we here have performed biochemical high throughput screening (HTS) with recombinantly expressed SARS-CoV-2 Mpro. A fluorescent assay was used to identify inhibitors in a compound library containing known drugs, bioactive molecules and natural products. These screens led to the identification of 13 inhibitors with IC50 values ranging from 0.2 µM to 23 µM. The screens confirmed several known SARS-CoV Mpro inhibitors as inhibitors of SARS-CoV-2 Mpro, such as the organo-mercuric compounds thimerosal and phenylmercuric acetate. Benzophenone derivatives could also be identified among the most potent screening hits. Additionally, Evans blue, a sulfonic acid-containing dye, could be identified as an Mpro inhibitor. The obtained compounds could be of interest as lead compounds for the development of future SARS-CoV-2 drugs.

Crystal structure of Thermus thermophilus methylenetetrahydrofolate dehydrogenase and determinants of thermostability.

The elucidation of mechanisms behind the thermostability of proteins is extremely important both from the theoretical and applied perspective. Here we report the crystal structure of methylenetetrahydrofolate dehydrogenase (MTHFD) from Thermus thermophilus HB8, a thermophilic model organism. Molecular dynamics trajectory analysis of this protein at different temperatures (303 K, 333 K and 363 K) was compared with homologous proteins from the less temperature resistant organism Thermoplasma acidophilum and the mesophilic organism Acinetobacter baumannii using several data reduction techniques like principal component analysis (PCA), residue interaction network (RIN) analysis and rotamer analysis. These methods enabled the determination of important residues for the thermostability of this enzyme. The description of rotamer distributions by Gini coefficients and Kullback-Leibler (KL) divergence both revealed significant correlations with temperature. The emerging view seems to indicate that a static salt bridge/charged residue network plays a fundamental role in the temperature resistance of Thermus thermophilus MTHFD by enhancing both electrostatic interactions and entropic energy dispersion. Furthermore, this analysis uncovered a relationship between residue mutations and evolutionary pressure acting on thermophilic organisms and thus could be of use for the design of future thermostable enzymes.