Interdiction of Protein Folding for Therapeutic Drug Development in SARS CoV-2

Accurate predictions of SARS-CoV-2 infectivity from comprehensive analysis

An unprecedented amount of SARS-CoV-2 data has been accumulated compared with previous infectious diseases, enabling insights into its evolutionary process and more thorough analyses. This study investigates SARS-CoV-2 features as it evolved to evaluate its infectivity. We examined viral sequences and identified the polarity of amino acids in the receptor binding motif (RBM) region. We detected an increased frequency of amino acid substitutions to lysine (K) and arginine (R) in variants of concern (VOCs). As the virus evolved to Omicron, commonly occurring mutations became fixed components of the new viral sequence. Furthermore, at specific positions of VOCs, only one type of amino acid substitution and a notable absence of mutations at D467 were detected. We found that the binding affinity of SARS-CoV-2 lineages to the ACE2 receptor was impacted by amino acid substitutions. Based on our discoveries, we developed APESS, an evaluation model evaluating infectivity from biochemical and mutational properties. In silico evaluation using real-world sequences and in vitro viral entry assays validated the accuracy of APESS and our discoveries. Using Machine Learning, we predicted mutations that had the potential to become more prominent. We created AIVE, a web-based system, accessible at https://ai-ve.org to provide infectivity measurements of mutations entered by users. Ultimately, we established a clear link between specific viral properties and increased infectivity, enhancing our understanding of SARS-CoV-2 and enabling more accurate predictions of the virus.

A Perspective on Interdicting in Protein Misfolding for Therapeutic Drug Design: Modulating the Formation of Nonlocal Contacts in α-Synuclein as a Strategy against Parkinson's Disease

In recent work we proposed that interdiction in the earliest contact-formation events along the folding pathway of key viral proteins could provide a novel avenue for therapeutic drug design. In this Perspective we explore the potential applicability of the protein folding interdiction strategy in the realm of neurodegenerative diseases with a specific focus on synucleinopathies. In order to fulfill this goal we review the interdiction proposal and its practical challenges, and we present new results concerning design strategies for possible peptide drugs that could be useful in preventing α-synuclein aggregation.

Potential Self-Peptide Inhibitors of the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) plays an essential role in viral replication, cleaving viral polyproteins into functional proteins. This makes Mpro an important drug target. Mpro consists of an N-terminal catalytic domain and a C-terminal α-helical domain (MproC). Previous studies have shown that peptides derived from a given protein sequence (self-peptides) can affect the folding and, in turn, the function of that protein. Since the SARS-CoV-1 MproC is known to stabilize its Mpro and regulate its function, we hypothesized that SARS-CoV-2 MproC-derived self-peptides may modulate the folding and the function of SARS-CoV-2 Mpro. To test this, we studied the folding of MproC in the presence of various self-peptides using coarse-grained structure-based models and molecular dynamics simulations. In these simulations of MproC and one self-peptide, we found that two self-peptides, the α1-helix and the loop between α4 and α5 (loop4), could replace the equivalent native sequences in the MproC structure. Replacement of either sequence in full-length Mpro should, in principle, be able to perturb Mpro function albeit through different mechanisms. Some general principles for the rational design of self-peptide inhibitors emerge: The simulations show that prefolded self-peptides are more likely to replace native sequences than those which do not possess structure. Additionally, the α1-helix self-peptide is kinetically stable and once inserted rarely exchanges with the native α1-helix, while the loop4 self-peptide is easily replaced by the native loop4, making it less useful for modulating function. In summary, a prefolded α1-derived peptide should be able to inhibit SARS-CoV-2 Mpro function.

Interdiction in the Early Folding of the p53 DNA-Binding Domain Leads to Its Amyloid-Like Misfolding

In this article, we investigate two issues: (a) the initial contact formation events along the folding pathway of the DNA-binding domain of the tumor suppressor protein p53 (core p53); and (b) the intermolecular events leading to its conversion into a prion-like form upon incubation with peptide P8(250-257). In the case of (a), the calculations employ the sequential collapse model (SCM) to identify the segments involved in the initial contact formation events that nucleate the folding pathway. The model predicts that there are several possible initial non-local contacts of comparative stability. The most stable of these possible initial contacts involve the protein segments 159AMAIY163 and 251ILTII255, and it is the only native-like contact. Thus, it is predicted to constitute “Nature’s shortcut” to the native structure of the core domain of p53. In the case of issue (b), these findings are then combined with experimental evidence showing that the incubation of the core domain of p53 with peptide P8(250-257), which is equivalent to the native protein segment 250PILTIITL257, leads to an amyloid conformational transition. It is explained how the SCM predicts that P8(250-257) effectively interdicts in the formation of the most stable possible initial contact and, thereby, disrupts the subsequent normal folding. Interdiction by polymeric P8(250-257) seeds is also studied. It is then hypothesized that enhanced folding through one or several of the less stable contacts could play a role in P8(250-257)-promoted core p53 amyloid misfolding. These findings are compared to previous results obtained for the prion protein. Experiments are proposed to test the hypothesis presented regarding core p53 amyloid misfolding.

Protein Folding Interdiction Strategy for Therapeutic Drug Development in Viral Diseases: Ebola VP40 and Influenza A M1

In a recent paper, we proposed the folding interdiction target region (FITR) strategy for therapeutic drug design in SARS-CoV-2. This paper expands the application of the FITR strategy by proposing therapeutic drug design approaches against Ebola virus disease and influenza A. We predict target regions for folding interdicting drugs on correspondingly relevant structural proteins of both pathogenic viruses: VP40 of Ebola, and matrix protein M1 of influenza A. Identification of the protein targets employs the sequential collapse model (SCM) for protein folding. It is explained that the model predicts natural peptide candidates in each case from which to start the search for therapeutic drugs. The paper also discusses how these predictions could be tested, as well as some challenges likely to be found when designing effective therapeutic drugs from the proposed peptide candidates. The FITR strategy opens a potential new avenue for the design of therapeutic drugs that promises to be effective against infectious diseases.

Comparative study of the interaction of ivermectin with proteins of interest associated with SARS-CoV-2: A computational and biophysical approach

The SARS-CoV-2 pandemic has accelerated the study of existing drugs. The mixture of homologs called ivermectin (avermectin-B1a [HB1a] + avermectin-B1b [HB1b]) has shown antiviral activity against SARS-CoV-2 in vitro. However, there are few reports on the behavior of each homolog. We investigated the interaction of each homolog with promising targets of interest associated with SARS-CoV-2 infection from a biophysical and computational-chemistry perspective using docking and molecular dynamics. We observed a differential behavior for each homolog, with an affinity of HB1b for viral structures, and of HB1a for host structures considered. The induced disturbances were differential and influenced by the hydrophobicity of each homolog and of the binding pockets. We present the first comparative analysis of the potential theoretical inhibitory effect of both avermectins on biomolecules associated with COVID-19, and suggest that ivermectin through its homologs, has a multiobjective behavior.

Identification of Two Early Folding Stage Prion Non-Local Contacts Suggested to Serve as Key Steps in Directing the Final Fold to Be Either Native or Pathogenic

The initial steps of the folding pathway of the C-terminal domain of the murine prion protein mPrP(90–231) are predicted based on the sequential collapse model (SCM). A non-local dominant contact is found to form between the connecting region between helix 1 and β-sheet 1 and the C-terminal region of helix 3. This non-local contact nucleates the most populated molten globule-like intermediate along the folding pathway. A less stable early non-local contact between segments 120–124 and 179–183, located in the middle of helix 2, promotes the formation of a less populated molten globule-like intermediate. The formation of the dominant non-local contact constitutes an example of the postulated Nature’s Shortcut to the prion protein collapse into the native structure. The possible role of the less populated molten globule-like intermediate is explored as the potential initiation point for the folding for three pathogenic mutants (T182A, I214V, and Q211P in mouse prion numbering) of the prion protein.

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.

Nano-size dependence in the adsorption by the SARS-CoV-2 spike protein over gold colloid

Gold nano-particles were coated with the spike protein (S protein) of SARS-CoV-2 and exposed to increasingly acidic conditions. Their responses were investigated by monitoring the surface plasmon resonance (SPR) band shift. As the external pH was gradually changed from neutral pH to pH ∼2 the peak of the SPR band showed a significant red-shift, with a sigmoidal feature implying the formation of the gold-protein aggregates. The coating of S protein changed the surface property of the gold enough to extract the coverage fraction of protein over nano particles, Θ, which did not exhibit clear nano-size dependence. The geometrical simulation to explain Θ showed the average axial length to be a = 7. 25 nm and b =8.00 nm when the S-protein was hypothesized as a prolate shape with spiking-out orientation. As the pH value externally hopped between pH∼3 and pH∼10, a behavior of reversible protein folding was observed for particles with diameters >30 nm. It was concluded that S protein adsorption conformation was impacted by the size (diameter, d) of a core nano-gold, where head-to-head dimerized S protein was estimated for d ≤ 80 nm and a parallel in opposite directions formation for d = 100 nm.