Assessment of Covalently Binding Warhead Compounds in the Validation of the Cytomegalovirus Nuclear Egress Complex as an Antiviral Target

Herpesviral nuclear egress is a regulated process of viral capsid nucleocytoplasmic release. Due to the large capsid size, a regular transport via the nuclear pores is unfeasible, so that a multistageregulated export pathway through the nuclear lamina and both leaflets of the nuclear membrane has evolved. This process involves regulatory proteins, which support the local distortion of the nuclear envelope. For human cytomegalovirus (HCMV), the nuclear egress complex (NEC) is determined by the pUL50–pUL53 core that initiates multicomponent assembly with NEC-associated proteins and capsids. The transmembrane NEC protein pUL50 serves as a multi-interacting determinant that recruits regulatory proteins by direct and indirect contacts. The nucleoplasmic core NEC component pUL53 is strictly associated with pUL50 in a structurally defined hook-into-groove complex and is considered as the potential capsid-binding factor. Recently, we validated the concept of blocking the pUL50–pUL53 interaction by small molecules as well as cell-penetrating peptides or an overexpression of hook-like constructs, which can lead to a pronounced degree of antiviral activity. In this study, we extended this strategy by utilizing covalently binding warhead compounds, originally designed as binders of distinct cysteine residues in target proteins, such as regulatory kinases. Here, we addressed the possibility that warheads may likewise target viral NEC proteins, building on our previous crystallization-based structural analyses that revealed distinct cysteine residues in positions exposed from the hook-into-groove binding surface. To this end, the antiviral and NEC-binding properties of a selection of 21 warhead compounds were investigated. The combined findings are as follows: (i) warhead compounds exhibited a pronounced anti-HCMV potential in cell-culturebased infection models; (ii) computational analysis of NEC primary sequences and 3D structures revealed cysteine residues exposed to the hook-into-groove interaction surface; (iii) several of the active hit compounds exhibited NEC-blocking activity, as shown at the single-cell level by confocal imaging; (iv) the clinically approved warhead drug ibrutinib exerted a strong inhibitory impact on the pUL50–pUL53 core NEC interaction, as demonstrated by the NanoBiT assay system; and (v) the generation of recombinant HCMV ∆UL50-ΣUL53, allowing the assessment of viral replication under conditional expression of the viral core NEC proteins, was used for characterizing viral replication and a mechanistic evaluation of ibrutinib antiviral efficacy. Combined, the results point to a ratelimiting importance of the HCMV core NEC for viral replication and to the option of exploiting this determinant by the targeting of covalently NEC-binding warhead compounds

Mutations in the B.1.1.7 SARS-CoV-2 Spike Protein Reduce Receptor-Binding Affinity and Induce a Flexible Link to the Fusion Peptide

The B.1.1.7 variant of the SARS-CoV-2 virus shows enhanced infectiousness over the wild type virus, leading to increasing patient numbers in affected areas. Amino acid exchanges within the SARS-CoV-2 spike protein variant of B.1.1.7 affect inter-monomeric contact sites within the trimer (A570D and D614G) as well as the ACE2-receptor interface region (N501Y), which comprises the receptor-binding domain (RBD) of the spike protein. However, the molecular consequences of mutations within B.1.1.7 on spike protein dynamics and stability or ACE2 binding are largely unknown. Here, molecular dynamics simulations comparing SARS-CoV-2 wild type with the B.1.1.7 variant revealed inter-trimeric contact rearrangements, altering the structural flexibility within the spike protein trimer. Furthermore, we found increased flexibility in direct spatial proximity of the fusion peptide due to salt bridge rearrangements induced by the D614G mutation in B.1.1.7. This study also implies a reduced binding affinity for B.1.1.7 with ACE2, as the N501Y mutation restructures the RBD–ACE2 interface, significantly decreasing the linear interaction energy between the RBD and ACE2. Our results demonstrate how mutations found within B.1.1.7 enlarge the flexibility around the fusion peptide and change the RBD–ACE2 interface. We anticipate our findings to be starting points for in depth biochemical and cell biological analyses of B.1.1.7

Computational decomposition reveals reshaping of the SARS‐CoV‐2–ACE2 interface among viral variants expressing the N501Y mutation

Abstract Variants of concern of the SARS‐CoV‐2 virus with an asparagine‐to‐tyrosine substitution at position 501 (N501Y) in the receptor‐binding domain (RBD) show enhanced infectivity compared to wild‐type, resulting in an altered pandemic situation in affected areas. These SARS‐Cov‐2 variants comprise the two Alpha variants (B.1.1.7, United Kingdom and B.1.1.7 with the additional E484K mutation), the Beta variant (B.1.351, South Africa), and the Gamma variant (P.1, Brazil). Understanding the binding modalities between these viral variants and the host cell receptor ACE2 allows to depict changes, but also common motifs of virus–host cell interaction. The trimeric spike protein expressed at the viral surface contains the RBD that forms the molecular interface with ACE2. All the above‐mentioned variants carry between one and three amino acid exchanges within the interface‐forming region of the RBD, thereby altering the binding interface with ACE2. Using molecular dynamics (MD) simulations and decomposition of intermolecular contacts between the RBD and ACE2, we identified phenylalanine 486, glutamine 498, threonine 500, and tyrosine 505 as important interface‐forming residues across viral variants. However, especially the N501Y exchange increased contact formation for this residue and also induced some local conformational changes. Comparing here, the in silico generated B.1.1.7 RBD–ACE2 complex with the now available experimentally solved structure reveals very similar behavior during MD simulation. We demonstrate, how computational methods can help to identify differences in conformation as well as contact formation for newly emerging viral variants. Altogether, we provide extensive data on all N501Y expressing SARS‐CoV‐2 variants of concern with respect to their interaction with ACE2 and how this induces reshaping of the RBD–ACE2 interface

Characterization of SARS-CoV-2 Escape Mutants to a Pair of Neutralizing Antibodies Targeting the RBD and the NTD

Mutations in the spike protein of SARS-CoV-2 can lead to evasion from neutralizing antibodies and affect the efficacy of passive and active immunization strategies. Immunization of mice harboring an entire set of human immunoglobulin variable region gene segments allowed to identify nine neutralizing monoclonal antibodies, which either belong to a cluster of clonally related RBD or NTD binding antibodies. To better understand the genetic barrier to emergence of SARS-CoV-2 variants resistant to these antibodies, escape mutants were selected in cell culture to one antibody from each cluster and a combination of the two antibodies. Three independently derived escape mutants to the RBD antibody harbored mutations in the RBD at the position T478 or S477. These mutations impaired the binding of the RBD antibodies to the spike protein and conferred resistance in a pseudotype neutralization assay. Although the binding of the NTD cluster antibodies were not affected by the RBD mutations, the RBD mutations also reduced the neutralization efficacy of the NTD cluster antibodies. The mutations found in the escape variants to the NTD antibody conferred resistance to the NTD, but not to the RBD cluster antibodies. A variant resistant to both antibodies was more difficult to select and only emerged after longer passages and higher inoculation volumes. VOC carrying the same mutations as the ones identified in the escape variants were also resistant to neutralization. This study further underlines the rapid emergence of escape mutants to neutralizing monoclonal antibodies in cell culture and indicates the need for thorough investigation of escape mutations to select the most potent combination of monoclonal antibodies for clinical use

XCR1 expression distinguishes human conventional dendritic cell type 1 with full effector functions from their immediate precursors

Dendritic cells (DCs) are major regulators of innate and adaptive immune responses. DCs can be classified into plasmacytoid DCs and conventional DCs (cDCs) type 1 and 2. Murine and human cDC1 share the mRNA expression of XCR1. Murine studies indicated a specific role of the XCR1-XCL1 axis in the induction of immune responses. Here, we describe that human cDC1 can be distinguished into XCR1$^{-}$ and XCR1$^{+}$ cDC1 in lymphoid as well as nonlymphoid tissues. Steady-state XCR1$^{+}$ cDC1 display a preactivated phenotype compared to XCR1$^{-}$ cDC1. Upon stimulation, XCR1$^{+}$ cDC1, but not XCR1$^{-}$ cDC1, secreted high levels of inflammatory cytokines as well as chemokines. This was associated with enhanced activation of NK cells mediated by XCR1$^{+}$ cDC1. Moreover, XCR1$^{+}$ cDC1 excelled in inhibiting replication of Influenza A virus. Further, under DC differentiation conditions, XCR1$^{-}$ cDC1 developed into XCR1$^{+}$ cDC1. After acquisition of XCR1 expression, XCR1$^{-}$ cDC1 secreted comparable level of inflammatory cytokines. Thus, XCR1 is a marker of terminally differentiated cDC1 that licenses the antiviral effector functions of human cDC1, while XCR1$^{-}$ cDC1 seem to represent a late immediate precursor of cDC1

Mutations in the B.1.1.7 SARS-CoV-2 Spike Protein Reduce Receptor-Binding Affinity and Induce a Flexible Link to the Fusion Peptide

The B.1.1.7 variant of the SARS-CoV-2 virus shows enhanced infectiousness over the wild type virus, leading to increasing patient numbers in affected areas. Amino acid exchanges within the SARS-CoV-2 spike protein variant of B.1.1.7 affect inter-monomeric contact sites within the trimer (A570D and D614G) as well as the ACE2-receptor interface region (N501Y), which comprises the receptor-binding domain (RBD) of the spike protein. However, the molecular consequences of mutations within B.1.1.7 on spike protein dynamics and stability or ACE2 binding are largely unknown. Here, molecular dynamics simulations comparing SARS-CoV-2 wild type with the B.1.1.7 variant revealed inter-trimeric contact rearrangements, altering the structural flexibility within the spike protein trimer. Furthermore, we found increased flexibility in direct spatial proximity of the fusion peptide due to salt bridge rearrangements induced by the D614G mutation in B.1.1.7. This study also implies a reduced binding affinity for B.1.1.7 with ACE2, as the N501Y mutation restructures the RBD-ACE2 interface, significantly decreasing the linear interaction energy between the RBD and ACE2. Our results demonstrate how mutations found within B.1.1.7 enlarge the flexibility around the fusion peptide and change the RBD-ACE2 interface. We anticipate our findings to be starting points for in depth biochemical and cell biological analyses of B.1.1.7

Mimicking titration experiments with MD simulations: A protocol for the investigation of pH-dependent effects on proteins

Protein structure and function are highly dependent on the environmental pH. However, the temporal or spatial resolution of experimental approaches hampers direct observation of pH-induced conformational changes at the atomic level. Molecular dynamics (MD) simulation strategies (e.g. constant pH MD) have been developed to bridge this gap. However, one frequent problem is the sampling of unrealistic conformations, which may also lead to poor pKa predictions. To address this problem, we have developed and benchmarked the pH-titration MD (pHtMD) approach, which is inspired by wet-lab titration experiments. We give several examples how the pHtMD protocol can be applied for pKa calculation including peptide systems, Staphylococcus nuclease (SNase), and the chaperone HdeA. For HdeA, pHtMD is also capable of monitoring pH-dependent dimer dissociation in accordance with experiments. We conclude that pHtMD represents a versatile tool for pKa value calculation and simulation of pH-dependent effects in proteins

Probing the Structure of the Escherichia coli Periplasmic Proteins HdeA and YmgD by Molecular Dynamics Simulations

HdeA and YmgD are structurally homologous proteins in the periplasm of Escherichia coli. HdeA has been shown to represent an acid-activated chaperone, whereas the function of YmgD has not yet been characterized. We performed pH-titrating molecular dynamics simulations (pHtMD) to investigate the structural changes of both proteins and to assess whether YmgD may also exhibit an unfolding behavior similar to that of HdeA. The unfolding pathway of HdeA includes partially unfolded dimer structures, which represent a prerequisite for subsequent dissociation. In contrast to the coupled unfolding and dissociation of HdeA, YmgD displays dissociation of the folded subunits, and the subunits do not undergo significant unfolding even at low pH values. The differences in subunit stability between HdeA and YmgD may be explained by the structural features of helix D, which represents the starting point of unfolding in HdeA. In summary, the present study suggests that YmgD either is not an acid-activated chaperone or, at least, does not require unfolding for activation

Conformational Dynamics of Herpesviral NEC Proteins in Different Oligomerization States

All herpesviruses use a heterodimeric nuclear egress complex (NEC) to transport capsids out of host cell nuclei. Despite their overall similar structure, NECs may differ significantly in sequence between different viruses. Up to now, structural information is limited to isolated NEC heterodimers and to large hexagonal lattices made up of hexagonal ring-like structures (“Hexagons”). The present study aimed to expand the existing structural knowledge with information on the dynamics of NECs from different viruses and in different oligomerization states. For this task, comparative molecular dynamics simulations were performed of the free NEC heterodimers from three different viruses (HCMV (human cytomegalovirus), HSV-1 (herpes simplex virus 1), and PRV (pseudorabies virus)). In addition, higher oligomerization states comprising two or six NEC heterodimers were characterized for HCMV and HSV-1. The study revealed that the isolated NEC heterodimers from α- (HSV-1, PRV) and β-herpesviruses (HCMV) differ significantly in their dynamics, which can be attributed to a poorly conserved interface region between the NEC subdomains. These differences become smaller for higher oligomerization states, and both HCMV and HSV-1 individual Hexagons exhibit a common region of enhanced dynamics, which might be of functional relevance for the formation of curved vesicle structures or the recognition of hexameric capsid proteins