Self-derived peptides from the SARS-CoV-2 spike glycoprotein disrupting shaping and stability of the homotrimer unit

The structural spike (S) protein from the SARS-CoV-2 β-coronavirus is shown to make different pre- and post-fusion conformations within its homotrimer unit. To support the ongoing novel vaccine design and development strategies, we report the structure-based design approach to develop self-derived S peptides. A dataset of crucial regions from the S protein were transformed into linear motifs that could act as the blockers or stabilizers for the S protein homotrimer unit. Among these distinct S peptides, the pep02 (537-QQFGRDIAD-545) and pep07 (821-RDLICAQKFNGLTVLPPLLTDE-842) were found making stable folded binding with the S protein (550–750 and 950–1050 regions). Upon inserting SARS-CoV-2 S variants in the peptide destabilized the complexed S protein structure, resulting an allosteric effect in different functional regions of the protein. Particularly, the molecular dynamics revealed that A544D mutation in the pep02 peptide induced instability for the complexed S protein, whereas the N943K variant from pep09 exhibited an opposite behavior. An increased protein-peptide binding affinity and the stable structural folding were observed in mutated systems, compared to that of the wild type systems. The presence of mutation has induced an “up” active conformation of the spike (RBD) domain, responsible for interacting the host cell receptor. Among the lower affinity peptide datasets (e.g., pep01), the S1 and S2 subunit in the protein formed an “open” conformation, whereas with higher affinity peptides (e.g., pep07) these domains gained a “closed” conformation. These findings propose that our designed self-derived S peptides could replace a single S protein monomer, blocking the homotrimer formation or inducing stability

Recognition dynamics of cancer mutations on the ERp57-Tapasin interface

Down regulation of the major histocompatibility class (MHC) I pathway plays an important role in tumour development, and can be achieved by suppression of HLA expression or mutations in the MHC peptide-binding pocket. The peptide-loading complex (PLC) loads peptides on the MHC-I molecule in a dynamic multi-step assembly process. The effects of cancer variants on ERp57 and tapasin components from the MHC-I pathway is less known, and they could have an impact on antigen presentation. Applying computational approaches, we analysed whether the ERp57-tapasin binding might be altered by missense mutations. The variants H408R(ERp57) and P96L, D100A, G183R(tapasin) at the protein–protein interface improved protein stability (ΔΔG) during the initial screen of 14 different variants. The H408R(ERp57) and P96L(tapasin) variants, located close to disulphide bonds, were further studied by molecular dynamics (MD). Identifying intramolecular a-a’ domain interactions, MD revealed open and closed conformations of ERp57 in the presence and absence of tapasin. In wild-type and mutant ERp57-tapasin complexes, residues Val97, Ser98, Tyr100, Trp405, Gly407(ERp57) and Asn94, Cys95, Arg97, Asp100(tapasin) formed common H-bond interactions. Moreover, comparing the H-bond networks for P96L and H408R with each other, suggests that P96L(tapasin) improved ERp57-tapasin binding more than the H408R(ERp57) mutant. During MD, the C-terminus domain (that binds MHC-I) in tapasin from the ERp57(H408R)-tapasin complex moved away from the PLC, whereas in the ERp57-tapasin(P96L) system was oppositely displaced. These findings can have implications for the function of PLC and, ultimately, for the presentation of MHC-I peptide complex on the tumour cell surface

The Elephant Evolved p53 Isoforms that Escape MDM2-Mediated Repression and Cancer

The p53 tumor suppressor is a transcription factor with roles in cell development, apoptosis, oncogenesis, aging, and homeostasis in response to stresses and infections. p53 is tightly regulated by the MDM2 E3 ubiquitin ligase. The p53-MDM2 pathway has coevolved, with MDM2 remaining largely conserved, whereas the TP53 gene morphed into various isoforms. Studies on prevertebrate ancestral homologs revealed the transition from an environmentally induced mechanism activating p53 to a tightly regulated system involving cell signaling. The evolution of this mechanism depends on structural changes in the interacting protein motifs. Elephants such as Loxodonta africana constitute ideal models to investigate this coevolution as they are large and long-living as well as having 20 copies of TP53 isoformic sequences expressing a variety of BOX-I MDM2-binding motifs. Collectively, these isoforms would enhance sensitivity to cellular stresses, such as DNA damage, presumably accounting for strong cancer defenses and other adaptations favoring healthy aging. Here we investigate the molecular evolution of the p53-MDM2 system by combining in silico modeling and in vitro assays to explore structural and functional aspects of p53 isoforms retaining the MDM2 interaction, whereas forming distinct pools of cell signaling. The methodology used demonstrates, for the first time that in silico docking simulations can be used to explore functional aspects of elephant p53 isoforms. Our observations elucidate structural and mechanistic aspects of p53 regulation, facilitate understanding of complex cell signaling, and suggest testable hypotheses of p53 evolution referencing Peto's Paradox

Structural determinants of peptide-dependent TAP1-TAP2 transit passage targeted by viral proteins and altered by cancer-associated mutations

The TAP1-TAP2 complex transports antigenic peptide substrates into the endoplasmic reticulum (ER). In ER, the peptides are further processed and loaded on the major histocompatibility class (MHC) I molecules by the peptide loading complex (PLC). The TAP transporters are linked with the PLC; a target for cancers and viral immune evasion. But the mechanisms whereby the cancer-derived mutations in TAP1-TAP2 or viral factors targeting the PLC, interfere peptide transport are only emerging. This study describes that transit of peptides through TAP can take place via two different channels (4 or 8 helices) depending on peptide length and sequence. Molecular dynamics and binding affinity predictions of peptide-transporters demonstrated that smaller peptides (8–10 mers; e.g. AAGIGILTV, SIINFEKL) can transport quickly through the transport tunnel compared to longer peptides (15-mer; e.g. ENPVVHFFKNIVTPR). In line with a regulated and selective peptide transport by TAPs, the immunopeptidome upon IFN-γ treatment in melanoma cells induced the shorter length (9-mer) peptide presentation over MHC-I that exhibit a relatively weak binding affinity with TAP. A conserved distance between N and C terminus residues of the studied peptides in the transport tunnel were reported. Furthermore, by adversely interacting with the TAP transport passage or affecting TAPNBD domains tilt movement, the viral proteins and cancer-derived mutations in TAP1-TAP2 may induce allosteric effects in TAP that block conformation of the tunnel (closed towards ER lumen). Interestingly, some cancer-associated mutations (e.g. TAP1R372Q and TAP2R373H) can specifically interfere with selective transport channels (i.e. for longer-peptides). These results provide a model for how viruses and cancer-associated mutations targeting TAP interfaces can affect MHC-I antigen presentation, and how the IFN-γ pathway alters MHC-I antigen presentation via the kinetics of peptide transport

Insights into the effects of cancer associated mutations at the UPF2 and ATP binding sites of NMD master regulator:UPF1

10.3390/ijms20225644International Journal of Molecular Sciences2022564

Highly Conserved Homotrimer Cavity Formed by the SARS-CoV-2 Spike Glycoprotein: A Novel Binding Site

An important stage in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) life cycle is the binding of the spike (S) protein to the angiotensin converting enzyme-2 (ACE2) host cell receptor. Therefore, to explore conserved features in spike protein dynamics and to identify potentially novel regions for drugging, we measured spike protein variability derived from 791 viral genomes and studied its properties by molecular dynamics (MD) simulation. The findings indicated that S2 subunit (heptad-repeat 1 (HR1), central helix (CH), and connector domain (CD) domains) showed low variability, low fluctuations in MD, and displayed a trimer cavity. By contrast, the receptor binding domain (RBD) domain, which is typically targeted in drug discovery programs, exhibits more sequence variability and flexibility. Interpretations from MD simulations suggest that the monomer form of spike protein is in constant motion showing transitions between an “up” and “down” state. In addition, the trimer cavity may function as a “bouncing spring” that may facilitate the homotrimer spike protein interactions with the ACE2 receptor. The feasibility of the trimer cavity as a potential drug target was examined by structure based virtual screening. Several hits were identified that have already been validated or suggested to inhibit the SARS-CoV-2 virus in published cell models. In particular, the data suggest an action mechanism for molecules including Chitosan and macrolides such as the mTOR (mammalian target of Rapamycin) pathway inhibitor Rapamycin. These findings identify a novel small molecule binding-site formed by the spike protein oligomer, that might assist in future drug discovery programs aimed at targeting the coronavirus (CoV) family of viruses

Identification of novel interferon responsive protein partners of human leukocyte antigen A (HLA-A) using cross-linking mass spectrometry (CLMS) approach

The interferon signalling system elicits a robust cytokine response against a wide range of environmental pathogenic and internal pathological signals, leading to induction of a subset of interferon-induced proteins. We applied DSS (disuccinimidyl suberate) mediated cross-linking mass spectrometry (CLMS) to capture novel protein–protein interactions within the realm of interferon induced proteins. In addition to the expected interferon-induced proteins, we identified novel inter- and intra-molecular cross-linked adducts for the canonical interferon induced proteins, such as MX1, USP18, OAS3, and STAT1. We focused on orthogonal validation of a cohort of novel interferon-induced protein networks formed by the HLA-A protein (H2BFS-HLA-A-HMGA1) using co-immunoprecipitation assay, and further investigated them by molecular dynamics simulation. Conformational dynamics of the simulated protein complexes revealed several interaction sites that mirrored the interactions identified in the CLMS findings. Together, we showcase a proof-of-principle CLMS study to identify novel interferon-induced signaling complexes and anticipate broader use of CLMS to identify novel protein interaction dynamics within the tumour microenvironment

The Binding Specificity of PAB1 with Poly(A) mRNA, Regulated by Its Structural Folding

The poly(A)-binding protein cytoplasmic 1 (PAB1 or PABPC1) protein is associated with the long poly(A) mRNA tails, inducing stability. Herein, we investigated the dynamics of the PABPC1 protein, along with tracing its mRNA binding specificity. During molecular dynamics simulations (MDS), the R176-Y408 amino acids (RRM3–4 domains; RNA recognition motifs) initiated a folded structure that resulted in the formation of different conformations. The RRM4 domain formed high-frequency intramolecular interactions, despite such induced flexibility. Residues D45, Y54, Y56, N58, Q88, and N100 formed long-lasting interactions, and specifically, aromatic residues (Y14, Y54, Y56, W86, and Y140) gained a unique binding pattern with the poly(A) mRNA. In addition, the poly(A) mRNA motif assembled a PABPC1-specific conformation, by inducing movement of the center three nucleotides to face towards RRM1–2 domains. The majority of the high-frequency cancer mutations in PAB1 reside within the RRM4 domain and amino acids engaging in high-frequency interactions with poly(A) mRNA were found to be preserved in different cancer types. Except for the G123C variant, other studied cancer-derived mutants hindered the stability of the protein. Molecular details from this study will provide a detailed understanding of the PABPC1 structure, which can be used to modulate the activity of this gene, resulting in production of mutant peptide or neoantigens in cancer

Molecular determinants and specificity of mrna with alternatively-spliced upf1 isoforms, influenced by an insertion in the 'regulatory loop'

The nonsense-mediated mRNA decay (NMD) pathway rapidly detects and degrades mRNA containing premature termination codons (PTCs). UP-frameshift 1 (UPF1), the master regulator of the NMD process, has two alternatively-spliced isoforms; one carries 353-GNEDLVIIWLR-363 insertion in the ‘regulatory loop (involved in mRNA binding)’. Such insertion can induce catalytic and/or ATPase activity, as determined experimentally; however, the kinetics and molecular level information are not fully understood. Herein, applying all-atom molecular dynamics, we probe the binding specificity of UPF1 with different GC-and AU-rich mRNA motifs and the influence of insertion to the viable control over UPF1 catalytic activity. Our results indicate two distinct conformations between 1B and RecA2 domains of UPF1: ‘open (isoform_2; without insertion)’ and ‘closed (isoform_1; with insertion)’. These structural movements correspond to an important stacking pattern in mRNA motifs, i.e., absence of stack formation in mRNA, with UPF1 isoform_2 results in the ‘open conformation’. Particularly, for UPF1 isoform_1, the increased distance between 1B and RecA2 domains has resulted in reducing the mRNA–UPF1 interactions. Lower fluctuating GC-rich mRNA motifs have better binding with UPF1, compared with AU-rich sequences. Except CCUGGGG, all other GC-rich motifs formed a 4-stack pattern with UPF1. High occupancy R363, D364, T627, and G862 residues were common binding GC-rich motifs, as were R363, N535, and T627 for the AU-rich motifs. The GC-rich motifs behave distinctly when bound to either of the isoforms; lower stability was observed with UPF1 isoform_2. The cancer-associated UPF1 variants (P533L/T and A839T) resulted in decreased protein–mRNA binding efficiency. Lack of mRNA stacking poses in the UPF1P533T system significantly decreased UPF1-mRNA binding efficiency and increased distance between 1B-RecA2. These novel findings can serve to further inform NMD-associated mechanistic and kinetic studies

Multivalent Display of SARS-CoV-2 Spike (RBD Domain) of COVID-19 to Nanomaterial, Protein Ferritin Nanocages

SARS-CoV-2, or COVID-19, has a devastating effect on our society, both in terms of quality of life and death rates; hence, there is an urgent need for developing safe and effective therapeutics against SARS-CoV-2. The most promising strategy to fight against this deadly virus is to develop an effective vaccine. Internalization of SARS-CoV-2 into the human host cell mainly occurs through the binding of the coronavirus spike protein (a trimeric surface glycoprotein) to the human angiotensin-converting enzyme 2 (ACE2) receptor. The spike-ACE2 protein–protein interaction is mediated through the receptor-binding domain (RBD) of the spike protein. Mutations in the spike RBD can significantly alter interactions with the ACE2 host receptor. Due to its important role in virus transmission, the spike RBD is considered to be one of the key molecular targets for vaccine development. In this study, a spike RBD-based subunit vaccine was designed by utilizing a ferritin protein nanocage as a scaffold. Several fusion protein constructs were designed in silico by connecting the spike RBD via a synthetic linker (different sizes) to different ferritin subunits (H-ferritin and L-ferritin). The stability and the dynamics of the engineered nanocage constructs were tested by extensive molecular dynamics simulation (MDS). Based on our MDS analysis, a five amino acid-based short linker (S-Linker) was the most effective for displaying the spike RBD over the surface of ferritin. The behavior of the spike RBD binding regions from the designed chimeric nanocages with the ACE2 receptor was highlighted. These data propose an effective multivalent synthetic nanocage, which might form the basis for new vaccine therapeutics designed against viruses such as SARS-CoV-2