In silico mechanistic analysis of IRF3 inactivation and high-risk HPV E6 species-dependent drug response

Unveiling the multifaceted realm of human papillomavirus: a comprehensive exploration of biology, interactions, and advances in cancer management

Human Papillomavirus (HPV), an extensive family of DNA viruses, manifests as a persistent global health challenge. Persistent HPV infection is now firmly established as a significant aetiological factor for a spectrum of malignancies. In this review, we examine the latest insights into HPV biology and its intricate relationship with the host. We delve into the complex dynamics of co-infections involving HPV alongside other viruses, such as HIV, EBV, and HSV, as well as the burgeoning role of the microbiome in cancer development. We also explore recent advancements in understanding the specific contributions of HPV in the development of various cancers, encompassing cancers of the anogenital region, head and neck, as well as breast, lung, and prostate. Moreover, we focus on the current preventive strategies, including vaccination and screening methods, and therapeutic interventions that range from traditional approaches like surgery and chemotherapy to emerging modalities such as targeted therapies and immunotherapies. Additionally, we provide a forward-looking view on the future directions of HPV research, highlighting potential areas of exploration to further our understanding and management of HPV and its associated cancers. Collectively, this review is positioned to deepen readers' understanding of HPV biology and its complex interplay with cancer biology. It presents innovative strategies for the prevention, management, and therapeutic intervention of HPV-associated malignancies.

Loxapine inhibits replication of hepatitis A virus in vitro and in vivo by targeting viral protein 2C

No antiviral drugs currently are available for treatment of infection by hepatitis A virus (HAV), a causative agent of acute hepatitis, a potentially life-threatening disease. Chemical screening of a small-compound library using nanoluciferase-expressing HAV identified loxapine succinate, a selective dopamine receptor D2 antagonist, as a potent inhibitor of HAV propagation in vitro. Loxapine succinate did not inhibit viral entry nor internal ribosome entry site (IRES)-dependent translation, but exhibited strong inhibition of viral RNA replication. Blind passage of HAV in the presence of loxapine succinate resulted in the accumulation of viruses containing mutations in the 2C-encoding region, which contributed to resistance to loxapine succinate. Analysis of molecular dynamics simulations of the interaction between 2C and loxapine suggested that loxapine binds to the N-terminal region of 2C, and that resistant mutations impede these interactions. We further demonstrated that administration of loxapine succinate to HAV-infected Ifnar1-/- mice (which lack the type I interferon receptor) results in decreases in the levels of fecal HAV RNA and of intrahepatic HAV RNA at an early stage of infection. These findings suggest that HAV protein 2C is a potential target for antivirals, and provide novel insights into the development of drugs for the treatment of hepatitis A.

Structure of the p53 degradation complex from HPV16

Human papillomavirus (HPV) is a significant contributor to the global cancer burden, and its carcinogenic activity is facilitated in part by the HPV early protein 6 (E6), which interacts with the E3-ligase E6AP, also known as UBE3A, to promote degradation of the tumor suppressor, p53. In this study, we present a single-particle cryoEM structure of the full-length E6AP protein in complex with HPV16 E6 (16E6) and p53, determined at a resolution of ~3.3 Å. Our structure reveals extensive protein-protein interactions between 16E6 and E6AP, explaining their picomolar binding affinity. These findings shed light on the molecular basis of the ternary complex, which has been pursued as a potential therapeutic target for HPV-driven cervical, anal, and oropharyngeal cancers over the last two decades. Understanding the structural and mechanistic underpinnings of this complex is crucial for developing effective therapies to combat HPV-induced cancers. Our findings may help to explain why previous attempts to disrupt this complex have failed to generate therapeutic modalities and suggest that current strategies should be reevaluated.

Implications of oral dysbiosis and HPV infection in head and neck cancer: from molecular and cellular mechanisms to early diagnosis and therapy

Head and neck cancer (HNC) is the sixth most common type of cancer, with more than half a million new cases annually. This review focuses on the role of oral dysbiosis and HPV infection in HNCs, presenting the involved taxons, molecular effectors and pathways, as well as the HPV-associated particularities of genetic and epigenetic changes and of the tumor microenvironment occurred in different stages of tumor development. Oral dysbiosis is associated with the evolution of HNCs, through multiple mechanisms such as inflammation, genotoxins release, modulation of the innate and acquired immune response, carcinogens and anticarcinogens production, generation of oxidative stress, induction of mutations. Thus, novel microbiome-derived biomarkers and interventions could significantly contribute to achieving the desideratum of personalized management of oncologic patients, regarding both early diagnosis and treatment. The results reported by different studies are not always congruent regarding the variations in the abundance of different taxons in HNCs. However, there is a consistent reporting of a higher abundance of Gram-negative species such as Fusobacterium, Leptotrichia, Treponema, Porphyromonas gingivalis, Prevotella, Bacteroidetes, Haemophilus, Veillonella, Pseudomonas, Enterobacterales, which are probably responsible of chronic inflammation and modulation of tumor microenvironment. Candida albicans is the dominant fungi found in oral carcinoma being also associated with shorter survival rate. Specific microbial signatures (e.g., F. nucleatum, Bacteroidetes and Peptostreptococcus) have been associated with later stages and larger tumor, suggesting their potential to be used as biomarkers for tumor stratification and prognosis. On the other hand, increased abundance of Corynebacterium, Kingella, Abiotrophia is associated with a reduced risk of HNC. Microbiome could also provide biomarkers for differentiating between oropharyngeal and hypopharyngeal cancers as well as between HPV-positive and HPV-negative tumors. Ongoing clinical trials aim to validate non-invasive tests for microbiome-derived biomarkers detection in oral and throat cancers, especially within high-risk populations. Oro-pharyngeal dysbiosis could also impact the HNCs therapy and associated side-effects of radiotherapy, chemotherapy, and immunotherapy. HPV-positive tumors harbor fewer mutations, as well as different DNA methylation pattern and tumor microenvironment. Therefore, elucidation of the molecular mechanisms by which oral microbiota and HPV infection influence the HNC initiation and progression, screening for HPV infection and vaccination against HPV, adopting a good oral hygiene, and preventing oral dysbiosis are important tools for advancing in the battle with this public health global challenge.

Molecular dynamics of the ERRγ ligand-binding domain bound with agonist and inverse agonist

Estrogen-related receptor gamma (ERRγ), the latest member of the ERR family, does not have any known reported natural ligands. Although the crystal structures of the apo, agonist-bound, and inverse agonist-bound ligand-binding domain (LBD) of ERRγ have been solved previously, their dynamic behavior has not been studied. Hence, to explore the intrinsic dynamics of the apo and ligand-bound forms of ERRγ, we applied long-range molecular dynamics (MD) simulations to the crystal structures of the apo and ligand-bound forms of the LBD of ERRγ. Using the MD trajectories, we performed hydrogen bond and binding free energy analysis, which suggested that the agonist displayed more hydrogen bonds with ERRγ than the inverse agonist 4-OHT. However, the binding energy of 4-OHT was higher than that of the agonist GSK4716, indicating that hydrophobic interactions are crucial for the binding of the inverse agonist. From principal component analysis, we observed that the AF-2 helix conformation at the C-terminal domain was similar to the initial structures during simulations, indicating that the AF-2 helix conformation is crucial with respect to the agonist or inverse agonist for further functional activity of ERRγ. In addition, we performed residue network analysis to understand intramolecular signal transduction within the protein. The betweenness centrality suggested that few of the amino acids are important for residue signal transduction in apo and ligand-bound forms. The results from this study may assist in designing better therapeutic compounds against ERRγ associated diseases.

Rational strategies for enhancing mAb binding to SARS-CoV-2 variants through CDR diversification and antibody-escape prediction

Since the emergence of SARS-CoV-2, dozens of variants of interest and half a dozen variants of concern (VOCs) have been documented by the World Health Organization. The emergence of these VOCs due to the continuous evolution of the virus is a major concern for COVID-19 therapeutic antibodies and vaccines because they are designed to target prototype/previous strains and lose effectiveness against new VOCs. Therefore, there is a need for time- and cost-effective strategies to estimate the immune escape and redirect therapeutic antibodies against newly emerging variants. Here, we computationally predicted the neutralization escape of the SARS-CoV-2 Delta and Omicron variants against the mutational space of RBD-mAbs interfaces. Leveraging knowledge of the existing RBD-mAb interfaces and mutational space, we fine-tuned and redirected CT-p59 (Regdanvimab) and Etesevimab against the escaped variants through complementarity-determining regions (CDRs) diversification. We identified antibodies against the Omicron lineage BA.1 and BA.2 and Delta variants with comparable or better binding affinities to that of prototype Spike. This suggests that CDRs diversification by hotspot grafting, given an existing insight into the Ag-Abs interface, is an exquisite strategy to redirect antibodies against preselected epitopes and combat the neutralization escape of emerging SARS-CoV-2 variants.

Human Papillomavirus-Associated Head and Neck Cancers. Where are We Now? A Systematic Review

BACKGROUND

Human papillomavirus targets the skin and mucous membranes, producing benign hyperplastic lesions and precancerous and cancerous lesions. An increasing number of head and neck cancersin particular, oropharyngeal squamous cell carcinoma, laryngeal squamous cell carcinoma, and oral squamous cell carcinoma, are attributable to HPV infection. HPV-induced HNCs typically affect younger, nonsmoking patients with no prior history of heavy alcohol use, more extensive sexual history, and higher socioeconomic status.

AIM

The purpose of the review is to present the most recent and well-established findings concerning HPV-induced head and neck cancers and consequently to provide medical specialists with essential information regarding the epidemiology, the role of HPV in HNC cancerogenesis, prevention, diagnosis, and treatment.

MATERIAL AND METHODS

All authors independently have searched The EMbase, Medline/Pubmed, and Cochrane databases by using the following keywords "head and neck cancer", "human papillomavirus", "HPV", "HPV biology", "oropharyngeal squamous cell carcinoma", "carcinogenesis", "transoral surgery", "robotic surgery". The last search was conducted in March 2022. The references of the publications of interest were also screened for relevant papers. There were no limitations in regard to the publication date.

CONCLUSION

Aiming to avoid the epidemic of HPV-induced HNC, it is paramount to improve the access to vaccination as well as resolve parental concerns regarding vaccine safety. Physicians should rely on reduced-dose radiation and aim to reduce the overall treatment time. Thanks to a more elaborate understanding of the genomic background of HPV-induced HNC, precision medicine could become a relevant part of patients' management. In comparison to traditional techniques and non-operative treatment, transoral robotic surgery (TORS) offers similar oncologic and functional outcomes, with a possible benefit on long-term quality of life. However, more research is needed to establish clear guidelines indicating when TORS resections should be supported with adjuvant therapy.

High-Risk Mucosal Human Papillomavirus 16 (HPV16) E6 Protein and Cutaneous HPV5 and HPV8 E6 Proteins Employ Distinct Strategies To Interfere with Interferon Regulatory Factor 3-Mediated Beta Interferon Expression

Persistent infection with some mucosal α-genus human papillomaviruses (HPVs; the most prevalent one being HPV16) can induce cervical carcinoma, anogenital cancers, and a subset of head and neck squamous cell carcinoma (HNSCC). Cutaneous β-genus HPVs (such as HPV5 and HPV8) associate with skin lesions that can progress into squamous cell carcinoma with sun exposure in Epidermodysplasia verruciformis patients and immunosuppressed patients. Here, we analyzed mechanisms used by E6 proteins from the α- and β-genus to inhibit the interferon-β (IFNB1) response. HPV16 E6 mediates this effect by a strong direct interaction with interferon regulatory factor 3 (IRF3). The binding site of E6 was localized within a flexible linker between the DNA-binding domain and the IRF-activation domain of IRF3 containing an LxxLL motif. The crystallographic structure of the complex between HPV16 E6 and the LxxLL motif of IRF3 was solved and compared with the structure of HPV16 E6 interacting with the LxxLL motif of the ubiquitin ligase E6AP. In contrast, cutaneous HPV5 and HPV8 E6 proteins bind to the IRF3-binding domain (IBiD) of the CREB-binding protein (CBP), a key transcriptional coactivator in IRF3-mediated IFN-β expression. IMPORTANCE Persistent HPV infections can be associated with the development of several cancers. The ability to persist depends on the ability of the virus to escape the host immune system. The type I interferon (IFN) system is the first-line antiviral defense strategy. HPVs carry early proteins that can block the activation of IFN-I. Among mucosal α-genus HPV types, the HPV16 E6 protein has a remarkable property to strongly interact with the transcription factor IRF3. Instead, cutaneous HPV5 and HPV8 E6 proteins bind to the IRF3 cofactor CBP. These results highlight the versatility of E6 proteins to interact with different cellular targets. The interaction between the HPV16 E6 protein and IRF3 might contribute to the higher prevalence of HPV16 than that of other high-risk mucosal HPV types in HPV-associated cancers.

SARS-CoV-2 pan-variant inhibitory peptides deter S1-ACE2 interaction and neutralize delta and omicron pseudoviruses

Approved neutralizing antibodies that target the prototype Spike are losing their potency against the emerging variants of concern (VOCs) of SARS-CoV-2, particularly Omicron. Although SARS-CoV-2 is continuously adapting the host environment, emerging variants recognize the same ACE2 receptor for cell entry. Protein and peptide decoys derived from ACE2 or Spike proteins may hold the pan-variant inhibitory potential. Here, we deployed interactive structure- and pharmacophore-based approaches to design short and stable peptides -Coronavirus Spike Neutralizing Peptides (CSNPs)- capable of neutralizing all SARS-CoV-2 VOCs. After in silico structural stability investigation and free energies perturbation of the isolated and target-bound peptides, nine candidate peptides were evaluated for the biophysical interaction through SPR assay. CSNP1, CSNP2, and Pep1 dose-dependently bind the S1 domain of the prototype Spike, whereas CSNP4 binds both S1 and ACE2. After safety and immunocytochemistry evaluation, peptides were probed for their pan-variant inhibitory effects. CSNP1, CSNP2, and CSNP4 inhibited all VOCs dose-dependently, whereas Pep1 had a moderate effect. CSNP2 and CSNP4 could neutralize the wild-type pseudovirus up to 80 % when treated at 0.5 µM. Furthermore, CSNP4 synergize the neutralization effect of monoclonal antibody and CSNP1 in Delta variant pseudovirus assay as they target different regions on the RBD. Thus, we suggest that CSNPs are SARS-CoV-2 pan-variant inhibitory candidates for COVID-19 therapy, which may pave the way for combating the emerging immune-escaping variants. We also propose that CSNP1/2-CSNP4 peptide cocktail or CSNP1/4 mAbs cocktail with no overlapping epitopes could be effective therapeutic strategies against COVID-19.

Unique Mode of Antiviral Action of a Marine Alkaloid against Ebola Virus and SARS-CoV-2

Lamellarin α 20-sulfate is a cell-impenetrable marine alkaloid that can suppress infection that is mediated by the envelope glycoprotein of human immunodeficiency virus type 1. We explored the antiviral action and mechanisms of this alkaloid against emerging enveloped RNA viruses that use endocytosis for infection. The alkaloid inhibited the infection of retroviral vectors that had been pseudotyped with the envelope glycoprotein of Ebola virus and SARS-CoV-2. The antiviral effects of lamellarin were independent of the retrovirus Gag-Pol proteins. Interestingly, although heparin and dextran sulfate suppressed the cell attachment of vector particles, lamellarin did not. In silico structural analyses of the trimeric glycoprotein of the Ebola virus disclosed that the principal lamellarin-binding site is confined to a previously unappreciated cavity near the NPC1-binding site and fusion loop, whereas those for heparin and dextran sulfate were dispersed across the attachment and fusion subunits of the glycoproteins. Notably, lamellarin binding to this cavity was augmented under conditions where the pH was 5.0. These results suggest that the final action of the alkaloid against Ebola virus is specific to events following endocytosis, possibly during conformational glycoprotein changes in the acidic environment of endosomes. Our findings highlight the unique biological and physicochemical features of lamellarin α 20-sulfate and should lead to the further use of broadly reactive antivirals to explore the structural mechanisms of virus replication.

Human Papillomaviruses as Infectious Agents in Gynecological Cancers. Oncogenic Properties of Viral Proteins

Human papillomaviruses (HPVs), which belong to the Papillomaviridae family, constitute a group of small nonenveloped double-stranded DNA viruses. HPV has a small genome that only encodes a few proteins, and it is also responsible for 5% of all human cancers, including cervical, vaginal, vulvar, penile, anal, and oropharyngeal cancers. HPV types may be classified as high- and low-risk genotypes (HR-HPVs and LR-HPVs, respectively) according to their oncogenic potential. HR-HPV 16 and 18 are the most common types worldwide and are the primary types that are responsible for most HPV-related cancers. The activity of the viral E6 and E7 oncoproteins, which interfere with critical cell cycle points such as suppressive tumor protein p53 (p53) and retinoblastoma protein (pRB), is the major contributor to HPV-induced neoplastic initiation and progression of carcinogenesis. In addition, the E5 protein might also play a significant role in tumorigenesis. The role of HPV in the pathogenesis of gynecological cancers is still not fully understood, which indicates a wide spectrum of potential research areas. This review focuses on HPV biology, the distribution of HPVs in gynecological cancers, the properties of viral oncoproteins, and the molecular mechanisms of carcinogenesis.

Structure-Guided Creation of an Anti-HA Stalk Antibody F11 Derivative That Neutralizes Both F11-Sensitive and -Resistant Influenza A(H1N1)pdm09 Viruses

The stalk domain of influenza virus envelope glycoprotein hemagglutinin (HA) constitutes the axis connecting the head and transmembrane domains, and plays pivotal roles in conformational rearrangements of HA for virus infection. Here we characterized molecular interactions between the anti-HA stalk neutralization antibody F11 and influenza A(H1N1)pdm09 HA to understand the structural basis of the actions and modifications of this antibody. In silico structural analyses using a model of the trimeric HA ectodomain indicated that the F11 Fab fragment has physicochemical properties, allowing it to crosslink two HA monomers by binding to a region near the proteolytic cleavage site of the stalk domain. Interestingly, the F11 binding allosterically caused a marked suppression of the structural dynamics of the HA cleavage loop and flanking regions. Structure-guided mutagenesis of the F11 antibody revealed a critical residue in the F11 light chain for the F11-mediated neutralization. Finally, the mutagenesis led to identification of a unique F11 derivative that can neutralize both F11-sensitive and F11-resistant A(H1N1)pdm09 viruses. These results raise the possibility that F11 sterically and physically disturbs proteolytic cleavage of HA for the ordered conformational rearrangements and suggest that in silico guiding experiments can be useful to create anti-HA stalk antibodies with new phenotypes.

Elucidating Molecular Interactions of Ten Natural Compounds Targeting E6 HPV High Risk Oncoproteins Using Microsecond Molecular Dynamics Simulations

BACKGROUND

Cervical cancer is a major public health issue worldwide, occurring in the vast majority of cases (85%) in low-income countries. Human papillomavirus (HPV) mainly infects the mucosal epithelium, and a small portion causes over 600,000 cases every year worldwide at various anatomical spots, mainly leading to anogenital and head and neck.

INTRODUCTION

The E6 oncoprotein encoded by cancer-associated alpha HPV can transform epithelial cells into tumorigenic tissue. Therapy for this infection and blocking of the HPV E6 oncoprotein could be provided with cost-effective and abundant natural products which are an exponentially growing topic in the literature. Finding an active natural compound that readily blocks HPV E6 oncoprotein which could be available for developing countries without expensive extraction processes or costly synthetic pathways is of major interest.

METHODS

Molecular dynamics simulation was performed using the most up-to-date AMBER protein force field ff14SB and a GPU enabled high performance computing cluster.

RESULTS

In this research, we present a study of the binding properties between 10 selected natural compounds that are readily available with two variants of the E6 oncoprotein types (HPV-16 and HPV-18) using 10+ microsecond molecular dynamics simulations.

CONCLUSION

Our results suggest that crocetin, ergosterol peroxide and κ-carrageenan natural products bind strongly to both HPV-16 and HPV-18 and could potentially serve as a scaffolding for further drug development.

In Silico Approaches: A Way to Unveil Novel Therapeutic Drugs for Cervical Cancer Management

Cervical cancer (CC) is the fourth most common pathology in women worldwide and presents a high impact in developing countries due to limited financial resources as well as difficulties in monitoring and access to health services. Human papillomavirus (HPV) is the leading cause of CC, and despite the approval of prophylactic vaccines, there is no effective treatment for patients with pre-existing infections or HPV-induced carcinomas. High-risk (HR) HPV E6 and E7 oncoproteins are considered biomarkers in CC progression. Since the E6 structure was resolved, it has been one of the most studied targets to develop novel and specific therapeutics to treat/manage CC. Therefore, several small molecules (plant-derived or synthetic compounds) have been reported as blockers/inhibitors of E6 oncoprotein action, and computational-aided methods have been of high relevance in their discovery and development. In silico approaches have become a powerful tool for reducing the time and cost of the drug development process. Thus, this review will depict small molecules that are already being explored as HR HPV E6 protein blockers and in silico approaches to the design of novel therapeutics for managing CC. Besides, future perspectives in CC therapy will be briefly discussed.

Interpretations on the Interaction between Protein Tyrosine Phosphatase and E7 Oncoproteins of High and Low-Risk HPV: A Computational Perception

The most prevalent and common sexually transmitted infection is caused by human papillomavirus (HPV) among sexually active women. Numerous genotypes of HPV are available, among which the major oncoproteins E6 and E7 lead to the progression of cervical cancer. The E7 oncoprotein interacts with cytoplasmic tumor suppressor protein PTPN14, which is the key regulator of cellular growth control pathways effecting the reduction of steady-state level. Disrupting the interaction between the tumor suppressor and the oncoprotein is vital to cease the development of cancer. Hence, the mechanism of interaction between E7 and tumor suppressor is explored through protein-protein and protein-ligand binding along with the conformational stability studies. The obtained results state that the LXCXE domain of HPV E7 of high and low risks binds with the tumor suppressor protein. Also, the small molecules bind in the interface of E7-PTPN14 that disrupts the interaction between the tumor suppressor and oncoprotein. These results were further supported by the dynamics simulation stating the stability over the bounded complex and the energy maintained during postdocking as well as postdynamics calculations. These observations possess an avenue in the drug discovery that leads to further validation and also proposes a potent drug candidate to treat cervical cancer caused by HPV.

PPI Modulators of E6 as Potential Targeted Therapeutics for Cervical Cancer: Progress and Challenges in Targeting E6

Advanced cervical cancer is primarily managed using cytotoxic therapies, despite evidence of limited efficacy and known toxicity. There is a current lack of alternative therapeutics to treat the disease more effectively. As such, there have been more research endeavors to develop targeted therapies directed at oncogenic host cellular targets over the past 4 decades, but thus far, only marginal gains in survival have been realized. The E6 oncoprotein, a protein of human papillomavirus origin that functionally inactivates various cellular antitumor proteins through protein–protein interactions (PPIs), represents an alternative target and intriguing opportunity to identify novel and potentially effective therapies to treat cervical cancer. Published research has reported a number of peptide and small-molecule modulators targeting the PPIs of E6 in various cell-based models. However, the reported compounds have rarely been well characterized in animal or human subjects. This indicates that while notable progress has been made in targeting E6, more extensive research is needed to accelerate the optimization of leads. In this review, we summarize the current knowledge and understanding of specific E6 PPI inhibition, the progress and challenges being faced, and potential approaches that can be utilized to identify novel and potent PPI inhibitors for cervical cancer treatment.

Structure of High-Risk Papillomavirus 31 E6 Oncogenic Protein and Characterization of E6/E6AP/p53 Complex Formation

The degradation of p53 is a hallmark of high-risk human papillomaviruses (HPVs) of the alpha genus and HPV-related carcinogenicity. The oncoprotein E6 forms a ternary complex with the E3 ubiquitin ligase E6-associated protein (E6AP) and tumor suppressor protein p53 targeting p53 for ubiquitination. The extent of p53 degradation by different E6 proteins varies greatly, even for the closely related HPV16 and HPV31. HPV16 E6 and HPV31 E6 display high sequence identity (∼67%). We report here, for the first time, the structure of HPV31 E6 bound to the LxxLL motif of E6AP. HPV16 E6 and HPV31 E6 are structurally very similar, in agreement with the high sequence conservation. Both E6 proteins bind E6AP and degrade p53. However, the binding affinities of 31 E6 to the LxxLL motif of E6AP and p53, respectively, are reduced 2-fold and 5.4-fold compared to 16 E6. The affinity of E6-E6AP-p53 ternary complex formation parallels the efficacy of the subsequent reaction, namely, degradation of p53. Therefore, closely related E6 proteins addressing the same cellular targets may still diverge in their binding efficiencies, possibly explaining their different phenotypic or pathological impacts.IMPORTANCE Variations of carcinogenicity of human papillomaviruses are related to variations of the E6 and E7 interactome. While different HPV species and genera are known to target distinct host proteins, the fine differences between E6 and E7 of closely related HPVs, supposed to target the same cellular protein pools, remain to be addressed. We compare the oncogenic E6 proteins of the closely related high-risk HPV31 and HPV16 with regard to their structure and their efficiency of ternary complex formation with their cellular targets p53 and E6AP, which results in p53 degradation. We solved the crystal structure of 31 E6 bound to the E6AP LxxLL motif. HPV16 E6 and 31 E6 structures are highly similar, but a few sequence variations lead to different protein contacts within the ternary complex and, as quantified here, an overall lower binding affinity of 31 E6 than 16 E6. These results align with the observed lower p53 degradation potential of 31 E6.

Mutations in the SARS-CoV-2 spike RBD are responsible for stronger ACE2 binding and poor anti-SARS-CoV mAbs cross-neutralization

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes coronavirus disease 2019 (COVID-19), is a novel beta coronavirus. SARS-CoV-2 uses spike glycoprotein to interact with host angiotensin-converting enzyme 2 (ACE2) and ensure cell recognition. High infectivity of SARS-CoV-2 raises questions on spike-ACE2 binding affinity and its neutralization by anti-SARS-CoV monoclonal antibodies (mAbs). Here, we observed Val-to-Lys417 mutation in the receptor-binding domains (RBD) of SARS-CoV-2, which established a Lys-Asp electrostatic interaction enhancing its ACE2-binding. Pro-to-Ala475 substitution and Gly482 insertion in the AGSTPCNGV-loop of RBD possibly hinders neutralization of SARS-CoV-2 by anti-SARS-CoV mAbs. In addition, we identified unique and structurally conserved conformational-epitopes on RBDs, which can be potential therapeutic targets. Collectively, we provide new insights into the mechanisms underlying the high infectivity of SARS-CoV-2 and development of effective neutralizing agents.

Subversion of Host Innate Immunity by Human Papillomavirus Oncoproteins

The growth of human papillomavirus (HPV)-transformed cells depends on the ability of the viral oncoproteins E6 and E7, especially those from high-risk HPV16/18, to manipulate the signaling pathways involved in cell proliferation, cell death, and innate immunity. Emerging evidence indicates that E6/E7 inhibition reactivates the host innate immune response, reversing what until then was an unresponsive cellular state suitable for viral persistence and tumorigenesis. Given that the disruption of distinct mechanisms of immune evasion is an attractive strategy for cancer therapy, the race is on to gain a better understanding of E6/E7-induced immune escape and cancer progression. Here, we review recent literature on the interplay between E6/E7 and the innate immune signaling pathways cGAS/STING/TBK1, RIG-I/MAVS/TBK1, and Toll-like receptors (TLRs). The overall emerging picture is that E6 and E7 have evolved broad-spectrum mechanisms allowing for the simultaneous depletion of multiple rather than single innate immunity effectors. The cGAS/STING/TBK1 pathway appears to be the most heavily impacted, whereas the RIG-I/MAVS/TBK1, still partially functional in HPV-transformed cells, can be activated by the powerful RIG-I agonist M8, triggering the massive production of type I and III interferons (IFNs), which potentiates chemotherapy-mediated cell killing. Overall, the identification of novel therapeutic targets to restore the innate immune response in HPV-transformed cells could transform the way HPV-associated cancers are treated.

Allosteric Regulation of HIV-1 Capsid Structure for Gag Assembly, Virion Production, and Viral Infectivity by a Disordered Interdomain Linker

The retroviral Gag capsid (Gag-CA) interdomain linker is an unstructured peptide segment connecting structured N-terminal and C-terminal domains. Although the region is reported to play roles in virion morphogenesis and infectivity, underlying molecular mechanisms remain unexplored. To address this issue, we determined biological and molecular phenotypes of HIV-1 CA linker mutants by experimental and in silico approaches. Among the nine linker mutants tested, eight exhibited attenuation of viral particle production to various extents mostly in parallel with a reduction in viral infectivity. Sucrose density gradient, confocal microscopy, and live-cell protein interaction analyses indicated that the defect is accompanied by attenuation of Gag-Gag interactions following Gag plasma membrane targeting in the cells. In silico analyses revealed distinct distributions of interaction-prone hydrophobic patches between immature and mature CA proteins. Molecular dynamics simulations predicted that the linker mutations can allosterically alter structural fluctuations, including the interaction surfaces apart from the mutation sites in both the immature and mature CA proteins. These results suggest that the HIV-1 CA interdomain linker is a cis-modulator of the CA interaction surfaces to optimize efficiency of Gag assembly, virion production, and viral infectivity.IMPORTANCE HIV-1 particle production and infection are highly ordered processes. Viral Gag proteins play a central role in the assembly and disassembly of viral molecules. Of these, capsid protein (CA) is a major contributor to the Gag-Gag interactions. CA consists of two structured domains, i.e., N-terminal (NTD) and C-terminal (CTD) domains, connected by an unstructured domain named the interdomain linker. While multiple regions in the NTD and CTD are reported to play roles in virion morphogenesis and infectivity, the roles of the linker region in Gag assembly and virus particle formation remain elusive. In this study, we showed by biological and molecular analyses that the linker region functions as an intramolecular modulator to tune Gag assembly, virion production, and viral infectivity. Our study thus illustrates a hitherto-unrecognized mechanism, an allosteric regulation of CA structure by the disordered protein element, for HIV-1 replication.

Linear and Rationally Designed Stapled Peptides Abrogate TLR4 Pathway and Relieve Inflammatory Symptoms in Rheumatoid Arthritis Rat Model

A mounting evidence exists for the despicable role of the aberrant immune response in the pathogenesis of rheumatoid arthritis (RA), where toll-like receptor 4 (TLR4) can activate synovial fibroblasts that lead to the chronic inflammation and joint destruction, thus making TLR4 a potent drug target in RA. We report that novel TLR4-antagonizing peptide, PIP2, inhibits the induction of inflammatory biomarkers in vitro as well as in vivo. Systemically, PIP2 inhibits the lipopolysaccharide (LPS)-elicited TNF-α, IL-6, and IL-12p40 in a mouse model. The rationally designed cyclic derivative, cPIP2, is capable of inhibiting LPS-induced proinflammatory cytokines at significantly lower concentration as compared to PIP2 (PIP2 IC₅₀ = 20 μM, cPIP2 IC₅₀ = 5 μM). Finally, cPIP2 was able to relieve the inflammatory symptoms and synovial tissue destruction in the RA rat model. Cumulatively, these data suggest that PIP2 and cPIP2 hold strong promise for the development of peptide-based immunotherapeutics that could be of great value in curbing TLR-related immune complications including RA.

Extraction of molecular features for the drug discovery targeting protein-protein interaction of Helicobacter pylori CagA and tumor suppressor protein ASSP2

Half of the world population is infected by the Gram-negative bacterium Helicobacter pylori (H. pylori). It colonizes in the stomach and is associated with severe gastric pathologies including gastric cancer and peptic ulceration. The most virulent factor of H. pylori is the cytotoxin-associated gene A (CagA) that is injected into the host cell. CagA interacts with several host proteins and alters their function, thereby causing several diseases. The most well-known target of CagA is the tumor suppressor protein ASPP2. The subdomain I at the N-terminus of CagA interacts with the proline-rich motif of ASPP2. Here, in this study, we carried out alanine scanning mutagenesis and an extensive molecular dynamics simulation summing up to 3.8 μs to find out hot spot residues and discovered some new protein-protein interaction (PPI)-modulating molecules. Our findings are in line with previous biochemical studies and further suggested new residues that are crucial for binding. The alanine scanning showed that mutation of Y207 and T211 residues to alanine decreased the binding affinity. Likewise, dynamics simulation and molecular mechanics with generalized Born surface area (MMGBSA) analysis also showed the importance of these two residues at the interface. A four-feature pharmacophore model was developed based on these two residues, and top 10 molecules were filtered from ZINC, NCI, and ChEMBL databases. The good binding affinity of the CHEMBL17319 and CHEMBL1183979 molecules shows the reliability of our adopted protocol for binding hot spot residues. We believe that our study provides a new insight for using CagA as the therapeutic target for gastric cancer treatment and provides a platform for a future experimental study.

Molecular modeling simulation studies reveal new potential inhibitors against HPV E6 protein

High-risk strains of human papillomavirus (HPV) have been identified as the etiologic agent of some anogenital tract, head, and neck cancers. Although prophylactic HPV vaccines have been approved; it is still necessary a drug-based treatment against the infection and its oncogenic effects. The E6 oncoprotein is one of the most studied therapeutic targets of HPV, it has been identified as a key factor in cell immortalization and tumor progression in HPV-positive cells. E6 can promote the degradation of p53, a tumor suppressor protein, through the interaction with the cellular ubiquitin ligase E6AP. Therefore, preventing the formation of the E6-E6AP complex is one of the main strategies to inhibit the viability and proliferation of infected cells. Herein, we propose an in silico pipeline to identify small-molecule inhibitors of the E6-E6AP interaction. Virtual screening was carried out by predicting the ADME properties of the molecules and performing ensemble-based docking simulations to E6 protein followed by binding free energy estimation through MM/PB(GB)SA methods. Finally, the top-three compounds were selected, and their stability in the E6 docked complex and their effect in the inhibition of the E6-E6AP interaction was corroborated by molecular dynamics simulation. Therefore, this pipeline and the identified molecules represent a new starting point in the development of anti-HPV drugs.

Structural-dynamic insights into the H. pylori cytotoxin-associated gene A (CagA) and its abrogation to interact with the tumor suppressor protein ASPP2 using decoy peptides

Helicobacter pylori (H. pylori) is one of the most extensively studied Gram-negative bacteria due to its implication in gastric cancer. The oncogenicity of H. pylori is associated with cytotoxin-associated gene A (CagA), which is injected into epithelial cells lining the stomach. Both the C- and N-termini of CagA are involved in the interaction with several host proteins, thereby disrupting vital cellular functions, such as cell adhesion, cell cycle, intracellular signal transduction, and cytoskeletal structure. The N-terminus of CagA interacts with the tumor-suppressing protein, apoptosis-stimulating protein of p53 (ASPP2), subsequently disrupting the apoptotic function of tumor suppressor gene p53. Here, we present the in-depth molecular dynamic mechanism of the CagA-ASPP2 interaction and highlight hot-spot residues through in silico mutagenesis. Our findings are in agreement with previous studies and further suggest other residues that are crucial for the CagA-ASPP2 interaction. Furthermore, the ASPP2-binding pocket possesses potential druggability and could be engaged by decoy peptides, identified through a machine-learning system and suggested in this study. The binding affinities of these peptides with CagA were monitored through extensive computational procedures and reported herein. While CagA is crucial for the oncogenicity of H. pylori, our designed peptides possess the potential to inhibit CagA and restore the tumor suppressor function of ASPP2.

Structural insights into the Middle East respiratory syndrome coronavirus 4a protein and its dsRNA binding mechanism

Middle East respiratory syndrome coronavirus (MERS-CoV) has evolved to navigate through the sophisticated network of a host's immune system. The immune evasion mechanism including type 1 interferon and protein kinase R-mediated antiviral stress responses has been recently attributed to the involvement of MERS-CoV protein 4a (p4a) that masks the viral dsRNA. However, the structural mechanism of how p4a recognizes and establishes contacts with dsRNA is not well explained. In this study, we report a dynamic mechanism deployed by p4a to engage the viral dsRNA and make it unavailable to the host immune system. Multiple variants of p4a-dsRNA were created and investigated through extensive molecular dynamics procedures to highlight crucial interfacial residues that may be used as potential pharmacophores for future drug development. The structural analysis revealed that p4a exhibits a typical αβββα fold structure, as found in other dsRNA-binding proteins. The α1 helix and the β1-β2 loop play a crucial role in recognizing and establishing contacts with the minor grooves of dsRNA. Further, mutational and binding free energy analyses suggested that in addition to K63 and K67, two other residues, K27 and W45, might also be crucial for p4a-dsRNA stability.

Structure-Activity Relationship in TLR4 Mutations: Atomistic Molecular Dynamics Simulations and Residue Interaction Network Analysis

Toll-like receptor 4 (TLR4), a vital innate immune receptor present on cell surfaces, initiates a signaling cascade during danger and bacterial intrusion. TLR4 needs to form a stable hexamer complex, which is necessary to dimerize the cytoplasmic domain. However, D299G and T399I polymorphism may abrogate the stability of the complex, leading to compromised TLR4 signaling. Crystallography provides valuable insights into the structural aspects of the TLR4 ectodomain; however, the dynamic behavior of polymorphic TLR4 is still unclear. Here, we employed molecular dynamics simulations (MDS), as well as principal component and residue network analyses, to decipher the structural aspects and signaling propagation associated with mutations in TLR4. The mutated complexes were less cohesive, displayed local and global variation in the secondary structure, and anomalous decay in rotational correlation function. Principal component analysis indicated that the mutated complexes also exhibited distinct low-frequency motions, which may be correlated to the differential behaviors of these TLR4 variants. Moreover, residue interaction networks (RIN) revealed that the mutated TLR4/myeloid differentiation factor (MD) 2 complex may perpetuate abnormal signaling pathways. Cumulatively, the MDS and RIN analyses elucidated the mutant-specific conformational alterations, which may help in deciphering the mechanism of loss-of-function mutations.

A structural insight into the negative effects of opioids in analgesia by modulating the TLR4 signaling: An in silico approach

Opioids are considered the gold standard therapy for pain. However, TLR-dependent negative effects in analgesia have highlighted the complexities in the pharmacodynamics of opioids. While successive studies have reported that morphine and Morphine-3-glucuronide (M3G) activate the TLR4 pathway, the structural details of this mechanism are lacking. Here, we have utilized various computational tools to reveal the structural dynamics of the opioid-bound TLR4/MD2 complex, and have proposed a potential TLR4 activation mechanism. Our results support previous findings, and include the novel insight that the stable binding of morphine and naloxone, but not M3G, in the MD2 cavity, is TLR4 dependent. Morphine interacts with MD2 near its Phe126 loop to induce the active conformation (MD2^C); however, this binding is likely reversible, and the complex gains stability upon interaction with TLR4. M3G also induces the MD2^C state, with both the Phe126 loop and the H1 loop being involved in MD2-M3G complex stability. Remarkably, naloxone, which requires TLR4 interaction for complex stability, switches the conformation of the gating loop to the inactive state (MD2°). Cumulatively, our findings suggest that ligand binding and receptor clustering occur successively in opioid-induced TLR4 signaling, and that MD2 plasticity and pocket hydrophobicity are crucial for the recognition and accommodation of ligands.

Structural and conformational insights into SOX2/OCT4-bound enhancer DNA: a computational perspective

The roles of SOX2 and OCT4 are critical in stem cell maintenance either in the context of iPSCs generation or cancer stem cell growth; therefore, it is imperative to study their cooperative binding and SOX2/OCT4-induced DNA conformational switching.