Next Article in Journal
Study of the Molecular Mechanisms of the Therapeutic Properties of Extracellular Vesicles
Next Article in Special Issue
Gibbs Energy and Gene Expression Combined as a New Technique for Selecting Drug Targets for Inhibiting Specific Protein–Protein Interactions
Previous Article in Journal
Dynamics of the Equine Placental DNA Methylome and Transcriptome from Mid- to Late Gestation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

HSV-1 Glycoprotein D and Its Surface Receptors: Evaluation of Protein–Protein Interaction and Targeting by Triazole-Based Compounds through In Silico Approaches

1
Dipartimento di Scienze e Tecnologie Biologiche Chimiche e Farmaceutiche (STEBICEF), Università degli Studi di Palermo, Via Archirafi 32, 90123 Palermo, Italy
2
Dipartimento di Scienze della Salute, Università degli Studi “Magna Græcia” di Catanzaro, Campus “S. Venuta”, Viale Europa, 88100 Catanzaro, Italy
3
Net4Science Academic Spin-Off, Università degli Studi “Magna Græcia” di Catanzaro, Campus “S. Venuta”, Viale Europa, 88100 Catanzaro, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(8), 7092; https://doi.org/10.3390/ijms24087092
Submission received: 3 March 2023 / Revised: 27 March 2023 / Accepted: 9 April 2023 / Published: 11 April 2023
(This article belongs to the Special Issue Protein-Protein Interactions as Druggable Targets: Recent Advances)

Abstract

:
Protein–protein interactions (PPI) represent attractive targets for drug design. Thus, aiming at a deeper insight into the HSV-1 envelope glycoprotein D (gD), protein–protein docking and dynamic simulations of gD-HVEM and gD-Nectin-1 complexes were performed. The most stable complexes and the pivotal key residues useful for gD to anchor human receptors were identified and used as starting points for a structure-based virtual screening on a library of both synthetic and designed 1,2,3-triazole-based compounds. Their binding properties versus gD interface with HVEM and Nectin-1 along with their structure-activity relationships (SARs) were evaluated. Four [1,2,3]triazolo[4,5-b]pyridines were identified as potential HSV-1 gD inhibitors, for their good theoretical affinity towards all conformations of HSV-1 gD. Overall, this study suggests promising basis for the design of new antiviral agents targeting gD as a valuable strategy to prevent viral attachment and penetration into the host cell.

1. Introduction

Herpesviridae is a large family of enveloped double-stranded DNA viruses that includes eight human herpesviruses (HHV) divided into three subfamilies: alphaherpesvirinae (herpes simplex virus [HSV] and varicella-zoster virus [VZV]); betaherpesvirinae (cytomegalovirus [CMV], HHV-6, and HHV-7); and gammaherpesvirinae (Epstein–Barr virus [EBV] and HHV-8) [1]. All herpesviruses, especially those belonging to the alphaherpesvirinae subfamily, are able to establish latent infections that can be reactivated by endogenous or exogenous stimuli, causing clinical symptoms [2]. HSV infections are widely spread around the world, thus representing a considerable public health issue. There are two main types of HSV that infect humans, HSV-1 and HSV-2. Both are responsible for genital herpes, although HSV-1 mostly causes oral infections and is acquired during childhood [3]. HSV-1 infection can also lead to complications, such as keratitis, encephalitis, meningitis, and systemic disease in neonates and immune-compromised patients [4,5]. Currently, antiviral drugs available for the treatment of herpesvirus infections include acyclovir (ACV), penciclovir, valacyclovir, famciclovir, foscarnet, and cidofovir, which terminate viral DNA synthesis by inhibiting the viral DNA polymerase [6]. However, these drugs are not effective in the complete eradication of the infection, but only in reducing the frequency and duration of the episodes. In addition, both the emergence of drug-resistant HSV-1 strains, in particular for immune-compromised patients, and drug toxicity issues, increase the need for new antiviral agents acting with different mechanisms of action from those currently in use [1,7].
An interesting strategy for the treatment of viral infections is to prevent viral attachment and entry into the host cell, a complex process mediated by the interaction of different viral envelope glycoproteins with specific host surface receptors [8]. The HSV-1 envelope contains 11 glycoproteins involved in the early stages of viral attachment and penetration [2]. Among them, gB, gC, gD, gH, and gL are considered essential for cell entry. Attachment to the host membrane begins with the interaction of gB and gC with cell surface heparan sulfate (HS) proteoglycans, which brings the virus closer to the cell [9,10]. The binding of gD to one of the specific cellular receptors and the consequent conformational change that occurs in the structure of gD trigger a cascade of events that promote the formation of the core fusion machinery with gB and the gH/gL complex that allows virus entry [11]. HSV-1 gD can bind three classes of receptors, depending on the cell type:
  • Herpesvirus entry mediator (HVEM), a member of the tumor necrosis factor receptor superfamily (TNFR) expressed on activate lymphocytes and in other human tissues including kidney, lung, and liver;
  • Nectin-1 and 2, immunoglobulin (Ig)-like cell adhesion molecules (CAMs) expressed on the surface of neuronal and epithelial cells [12];
  • 3-O-sulfated heparan sulfate (3-OS HS), whose biological role has not yet been well clarified [13].
Due to the crucial role that gD plays in host cell fusion, a thorough knowledge of its structural features and of the molecular mechanisms regulating its activity may be useful to develop new inhibitors.
The current study aims at a deeper insight into gD and its cellular receptor interactions by means of computational methods, taking advantage of the HSV-1 gD experimentally solved structures in complex to HVEM and Nectin-1. The complexes resulting from protein–protein docking experiments were submitted to molecular dynamic simulations (MDs). GBPM and free energy calculation analyses revealed the most stable complexes and their pivotal points of interaction. Considering that protein–protein interactions (PPI) represent attractive targets for drug design, a structure-based virtual screening (SBVS) was performed taking into account the key residues useful for gD to anchor human receptors.
Much effort has been devoted by our research group in the identification of small molecules as bioactive compounds [14,15,16,17] and among them, we have recently reported a small set of [1,2,3]triazolo[4,5-h][1,6]naphthyridines as promising photosensitizer agents [18]. Both 1,2,3-triazole based compounds [19,20,21,22] and naphthyridine derivatives have been reported in the literature for their promising antiviral activity, refs. [23,24,25,26,27] more specifically, 1,2,3 triazole ring and 1,6-naphthyridine core appear as recurring moieties in heterocyclic compounds exhibiting anti-HSV-1 properties at the micromolar level [28,29,30,31,32]. Hence, we used a small set of 12 synthetic [1,2,3]triazolo[4,5-h][1,6]naphthyridines (I, Figure 1 and 1–12, Table 1) to evaluate their ability to target gD binding interface with HVEM and Nectin-1. Moreover, 63 novel [1,2,3]triazolo[4,5-h][1,6]naphthyridines were designed (19–81, Table 1) in order to evaluate their binding properties into gD pockets and to explore their structure-activity relationships (SARs). Furthermore, since triazolopyridine systems have been extensively studied as interesting scaffolds exhibiting antiviral activities, [33,34,35] we also investigated the binding affinity of [1,2,3]triazolo[4,5-b]pyridines (II, Figure 1 and 13–18 Table 1), synthetic precursors of the aforementioned triazolo-naphthyridines, versus gD interface with HVEM and Nectin-1. The workflow is reported in Figure 1.

2. Results and Discussion

In order to investigate the binding affinity of gD for its surface receptors and to gain useful structural insights regarding the pivotal interactions occurring at the gD interface, we performed a knowledge-based protein–protein docking employing HADDOCK 2.4 web-server [36]. Due to the involvement of gD both in cell adhesion function and viral entry mechanism depending on its cellular receptors, we considered the X-ray structure of gD both in complex to HVEM (PDB code: 1JMA) and Nectin-1 (PDB code: 3U82), respectively at 2.65 and 3.16 Å resolution. The complexes were prepared using Protein Preparation Wizard and uploaded in PDB format to the HADDOCK server. The input docking parameters were discussed in the Section 3.1. All the generated complexes were clustered according to their HADDOCK score, calculated as a weighted sum of a variety of energy terms (including van der Waals, electrostatic, desolvation, and restraint violation energies) and buried surface area (BSA); the Z-score value was also calculated to select the best cluster with respect to all obtained clusters: the most negative Z-score is indicative of the top ranked cluster. For the gD-HVEM and gD-Nectin-1 complexes, we obtained 188 structures gathered in 9 clusters (Figure S1), and 198 structures in 2 clusters (Figure S2), respectively. For both complexes, the docking results of each cluster were reported in Table 2.
The best clusters (the lowest in energy) were analysed in terms of interactions using Maestro graphical user interface (Figure 2) and were aligned to the experimental structures through the Protein Structure Alignment tool.
By overlapping (Figure S3), it resulted that for gD-HVEM and gD-Nectin-1 complexes, the RMSD value between the best docked pose and the X-ray structure was 1.54 and 0.80 Å, respectively. All the residues involved in the PPIs are summarized in Table 3.
The best obtained docked structures for each complex were refined using Protein Preparation Wizard and energy minimized with OPLS_2005 force field [37]. In order to explore any potential conformational changes of the gD interface, we submitted 100 ns of MDs for both complexes. The stability of MDs trajectories was monitored by the RMSD trend of the protein’s backbone atoms from its initial conformation. The average RMSD values of gD-HVEM and gD-Nectin-1 complexes were 3.05 Å and 2.50 Å, respectively. By monitoring the distances between the interface’s residues for gD-HVEM, we observed that several interactions previously detailed in Table 3 were maintained during MDs run, except for Asn15, Gly19, and Lys122 of gD, which were subjected to greater fluctuations, thus preventing the contacts with HVEM residues. Instead, Ala7 gained contacts with Ser20 during the whole MDs (Figure 3).
Regarding gD in complex to Nectin-1, the interactions engaged by Gln27 and Gln132 of gD at the interface with Nectin-1 were lost, but it was found that Tyr38 was located in proximity to Gly86 of Nectin-1 during MDs (Figure 4).
To better characterize the gD interface, a deeper analysis using GBPM was performed on all frames of MDs for each system. This method helps to map the key hotspots responsible for PPI by combining GRID molecular interaction fields (MIFs) according to the GRAB tool algorithm [38]. We considered gD as guest and HVEM and Nectin-1 as hosts. Three GRID probes, such as DRY, N1, and O, were chosen to mimic the hydrophobic, H-bond donor, and acceptor areas, respectively. Taking into account an energy threshold above the 30% from the global energy minimum GRID points, we summarized the pivotal residues up to 3 Å from GBPM points. The contribution of each residue was derived by the summa of its GBPM points energy in the matching frames (see Tables S1 and S2). After calculating the average score based on the total number of frames, the key hotspots were split into quartiles to increase the clarity for the reader: quartile 1 (Q1) includes the residues with the major contribution to PPI until quartile 4 (Q4), which contains residues with the weakest interactions, during the entire trajectories (Table 4).
For both protein–protein complexes, 100 frames extracted by MDs were adopted to calculate the relative binding free energy (ΔGbind) using Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) methodology, as applied in a recent study [39]. The results of the calculated ΔGbind trend for gD-HVEM and gD-Nectin-1 are depicted in Figure 5. Even showing a similar energy profile, the average values of −134.36 and −138.26 kcal/mol were turned out for HVEM and Nectin-1 in complex to gD, respectively, thus resulting in a more stabilizing effect of gD towards Nectin-1.
After investigating the PPIs’ structural details, we focused on the key residues of gD liable for the interaction with the analyzed human receptors to find potential ligands able to prevent the connection with them. According to the literature [40], several gD N-terminus residues at position 7–15 and 26–29 are responsible for the bind with HVEM, whereas two residues, such as Asp26 and Tyr38, also belonging to gD N-terminus, are relevant for the connection with Nectin-1 [41]. Apart from the shared Asp26 and Gln27 residues, both human receptors interact with different portions of gD. It was underlined that HVEM and Nectin-1 show non-reciprocal competition for binding to gD [42]. Indeed, the interaction between HVEM and gD induces a conformational change in gD that results in the formation of the N-terminal hairpin structure that masks the binding site of Nectin-1 [42]. Similarly, the interaction with Nectin-1 induces a new conformation of the gD N-terminus that prevents the HVEM binding [4]. Considering therefore that gD is not able to bind both receptors at the same time [41], we evaluated the possible conformational changes of gD during MDs. As previously reported [43], MDs trajectories were clustered based on RMSD matrix using the average hierarchical clustering linkage method, obtaining three representative structures of gD. Taking into account that the pivotal interactions of gD with HVEM and Nectin-1 involve distinct residues, we used all three representative structures generated from each complex (Figure S4) as starting point for SBVS.
Thus, we performed molecular docking simulations of an in-house small library of 12 [1,2,3] triazolo [4,5-h] [1,6] naphthyridines (1–12, Table 1) with their precursors [1,2,3] triazolo [4,5-b] pyridines (13–18, Table 1) on each binding site, separately, using the docking program Glide in SP mode. Additionally, with the aim of exploring the influence of electron-withdrawing and/or electron-donating substituents at N-3 and C-8 positions of the triazolo-naphthyridine core on gD-ligand interaction, a focused library of 63 triazolo [4,5-h] [1,6] naphthyridines (19–81, Table 1) was also designed, along with the corresponding [1,2,3]triazolo [4,5-b] pyridines (82–85, Table 1).
We calculated the average G-score value for each cluster, aiming to focus on the compounds able to better recognize both Pocket 1 and 2 (Table S3). From the docking results, it clearly emerged that the tricyclic [1,2,3] triazolo [4,5-h] [1,6] naphthyridine moiety showed a lower ability to recognize both the binding pockets compared to the bicyclic triazolo [4,5-b] pyridine derivatives. Furthermore, regardless of the electronic nature of the substituents at the N-3 phenyl ring of the triazolo-naphthyridine core, no significant differences could be observed.
Concerning triazolo [4,5-b] pyridine derivatives, the presence of the 3,4,5-trimethoxylphenyl ring at N-3 reduced the theoretical binding affinity to gD. Conversely, better results were obtained with p-phenyl substituted derivatives bearing the hydroxymethyl group at C-6 (14, 16, 83, 85).
Overall, among the 85 investigated compounds, we focused on the binding mode of the compounds able to recognize all the representative structures of both pockets with average G-score values lower than −5.00 kcal/mol. Accordingly, four triazole-pyridines derivatives, 14 (Figure S5), 16 (Figure S6), 83 (Figure S7), and 85 (Figure S8) showed a favored energetic profile in complex with all gD conformations. In detail, considering the most populated cluster resulted for gD-Pocket 1, the hydroxyl group of 14 (Figure 6A), 16 (Figure 6B), 83 (Figure 6C), and 85 (Figure 6D) was anchored to Leu25. Moreover, 16 and 83 formed an additional H-bond between the amino group and Leu25.
In gD-Pocket 2, 14, 16, 83, and 85 engaged several H-bonds with Leu28, Asp30, Asn227, and a π-π stacking interaction with Phe223, probably due to the absence of the N-terminal extension (Figure 7).
The best docked poses of 14, 16, 83, and 85 in complex with the most representative cluster of gD-Pocket 1 and gD-Pocket 2 were submitted to 500 ns of MDs using Desmond [44]. The results of MDs were investigated in terms of stability and conformational flexibility in the presence of the selected compounds. The stability of the complexes was evaluated by calculating the RMSD of the protein’s backbone atoms from its initial to final conformation over the whole simulation. By RMSD analysis, we observed that the most promising compounds maintained overall stability throughout MDs in both Pockets, as shown in Figure 8. In particular, for Pocket 1, the average RMSD values of 2.27, 3.09, 2.58, and 2.59 Å (Figure 8A) were computed for 14, 16, 83, and 85, respectively. On the other hand, for Pocket 2, we observed average RMSD values in the range of 2.17–2.51 Å (Figure 8B).
Furthermore, for each system, the binding free energy ΔGbind was calculated using MM/GBSA methodology, extracting 100 snapshots from 500 ns of MDs. MM/GBSA analysis showed that the average calculated ΔGbind of 14, 16, 83, and 85 complexed with gD-Pocket 1 were −38.88, −36.52, −40.22, and −31.14 kcal/mol, respectively, during the entire trajectories. The average ΔGbind of 14, 16, 83, and 85 in complex to gD-Pocket 2 were −31.78, −29.71, −36.28, and −34.63 kcal/mol, respectively. ADME parameters of the most promising compounds were predicted using the SwissADME server, and the obtained data are reported in Table 5.
All the investigated compounds are achiral, fit Lipinski’s rule of five, and were not found to be pan assay interference compounds (PAINS), thus resulting suitable for further HSV-1 in vitro investigation.

3. Materials and Methods

3.1. Computational Methods

3.1.1. Protein–protein Preparation of gD-HVEM and gD-Nectin-1 and Docking Simulations

All computational studies were carried out using Schrödinger Suite 2018-1 [45]. The X-ray crystallographic structures of gD in complex with HVEM (PDB code: 1JMA) [46] and with Nectin-1 (PDB code: 3U82) [41] were used. Both gD structures were prepared and optimized using the Maestro Protein Preparation Wizard tool [47] with OPLS_2005 [37] as force field at pH 7.4. The structures were also optimized by the addition of missing loops using Prime software [44,48] and the determination of the protonation state of the ionizable amino acid residues by means of the Epik program [49]. Knowledge-based protein–protein docking of gD with HVEM and Nectin-1 was performed with HADDOCK 2.4 (High Ambiguity-Driven biomolecular DOCKing) web-server [36]. For each complex, the sampling parameters were as follows: 1000 structures for rigid-body docking, 200 structures for the final refinement, and a cut-off equal to 5.0 to define neighboring flexible regions. For the complex gD-HVEM from 1JMA model, gD amino acids Ala7, Ser8, Leu9, Lys10, Met11, Ala12, Asp13, Pro14, Asn15, Val24, Leu25, Asp26, Gln27, Leu28, Thr29, Asp30, Pro31, Pro32 and HVEM amino acids Pro17, Lys18, Cys19, Ser20, Pro21, Gly22, Tyr23, Arg24, Val25, Lys26, Gly30, Glu31, Leu32 Thr33, Gly34, Thr35, Val36, Cys37, Glu38, Pro39, Ser74, Arg75, and Thr76 were considered as active residues, whereas passive residues were automatically identified as residues surrounding the active ones before submitting the docking job. For the complex gD-Nectin-1 from 3U82 model, gD amino acids Pro23, Leu25, Gln27, Arg36, Val37, Tyr38, His39, Gln132, Val214, Asp215, Ser216, Ile217, Gly218, Met219, Leu220, Pro221, Arg222, Phe223, Thr230, Val231, Tyr234 and Nectin-1 amino acids Ser59, Lys61, Thr63, Gln64, Thr66, Gln68, Lys75, Gln76, Asn77, Ile80, Tyr81, Asn82, Met85, Gly86, Val87, Ser88, Leu90, Glu125, Ala127, Thr128, Phe129, Pro130, Thr131, Gly132, and Asn133 were considered as active residues, whereas passive residues were again automatically identified. For each complex, the docked structures were loaded on Maestro interface of the Schrödinger software in PDB format for visual inspection and for the following in silico analysis.

3.1.2. MDs, GBPM, and MM/GBSA Calculations of gD-HVEM and gD-Nectin-1 Complexes

The best docked poses of each complex were submitted to 100 ns of MDs using Desmond ver. 4.2 [44]. To perform simulations in an aqueous biological environment, an appropriate system was built using OPLS_2005 as force field and an orthorhombic box with TIP4P water model extending of 10 Å outside the complex in all sides. The systems were maintained at a salt concentration of 0.15 M by adding appropriate Cl counter ions to neutralize them to maintain the physiological condition. After optimization of the solvated models, we relaxed the systems with the Martyna−Tobias−Klein isobaric−isothermal ensemble (MTK_NPT). Finally, 100 ns unconstrained MDs were carried out using the following conditions: the NPT ensemble, a constant temperature of 300 K, a pressure of 1 bar, and a recording interval equal to 100 ps both for energy and for trajectory collecting 1000 frames for each simulation.
For both complexes, all frames were considered for GBPM analysis [38]. As previously reported [50], in order to evidence hydrophobic and hydrogen bond donors and acceptors spots, we used DRY, N1, and O GRID probes, respectively. For each complex, gD was seen as guest and HVEM and Nectin-1 as hosts. The selected residues at the interface covered a maximum distance of 3 Å from the most relevant interaction energy points (GBPM features) of the computed molecular interaction fields (MIFs). After selecting an energy cutoff 30% above the global minimum, the pivotal hotpots were resulted by the summa of the related GBPM features interaction energy.
One thousand snapshots from 100 ns of MDs were applied for the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) free energy calculations [51] based on the following equation:
ΔGbind = Gcomp − Gpro −Glig = ΔEele + ΔEvdw + ΔEint + ΔEGB + ΔEsurf
where Gcomp, Gpro and Glig denotes the free energy of the complex, protein, and the ligand; by splitting the energy contribution, it referred to ΔEele, ΔEvdw and ΔEint as the gas-phase interaction energy between protein and ligand, thus including the electrostatic energy term, the van der Waals energy term, and the bond, angle, and dihedral terms, respectively. On the other hand, ΔEGB and ΔEsurf indicate the polar and nonpolar desolvation free energy, respectively. The implicit solvation was calculated using the GB model [52], and the non-polar solvation energy was calculated using the solvent accessible surface area algorithm. The ΔGbind reported in this study omitted the entropy contribution due to its relatively high computational demand and the lack of information of the conformational entropy that could lead to the introduction of additional error into the results [53].

3.2. Structure-Based Virtual Screening of [1,2,3]triazolo[4,5-h][1,6]naphthyridines, and [1,2,3]triazolo[4,5-b]Pyridines on gD-Pocket 1 and gD-Pocket 2

The MDs trajectories were clustered based on the RMSD matrix of backbone atoms, and we obtained three representative structures for each complex. After removing the HVEM and Nectin-1 structures, we used a total of six gD conformations for the SBVS. The target binding sites were defined by a regular grid of about 20 Å centered on the residues responsible for binding with the cell receptor [8,40,41,54]. The residues that defined the gD binding site at the interface with HVEM, called for clarity “gD-Pocket 1”, were as follows: Ala7, Ser8, Leu9, Lys10, Met11, Ala12, Asp13, Pro14, Asn15, Val24, Leu25, Asp26, Gln27, Leu28, Thr29, Asp30, Pro31, and Pro32. The gD binding site at the interface with Nectin-1, called “gD-Pocket 2”, was characterized by the following residues: Pro23, Leu25, Gln27, Arg36, Val37, Tyr38, His39, Gln132, Val214, Asp215, Ser216, Ile217, Gly218, Met219, Leu220, Pro221, Arg222, Phe223, Thr230, Val231, and Tyr234.
The selected ligands were taken from an in-house small library of [1,2,3]triazolo[4,5-h][1,6]naphthyridines 1–12 including also their synthetic precursor [1,2,3]triazolo[4,5-b]pyridines 13–18 [18]. Moreover, a focused library of 63 triazolo [4,5-h][1,6]naphthyridines (19–81, Table 1) was also designed, along with the corresponding [1,2,3]triazolo[4,5-b]pyridines (82–85, Table 1). The 2D structures, reported in Table 1, were drawn using the ChemDraw Ultra 7.0 software and converted into 3D form using the import structures panel from Schrodinger maestro interface.
All compounds were optimized via the Ligprep module, [55] considering their ionization state at pH 7.4, and energy minimized using OPLS_2005 as force field [37].
As reported in other studies [56,57,58,59], the docking simulations of our focused library were computed using the Glide [60] ligand flexible algorithm at the standard-precision (SP) level, generating 10 possible poses for each site. The best docked poses for gD-Pocket 1 and gD-Pocket 2 were submitted to 500 ns of MDs in order to define the structural and energy profile of the best ligands in complex with both gD-Pockets. The simulations were carried out under the above-mentioned conditions. All simulations were performed by Desmond package [44] and “Simulation Interactions Diagram” panel was used as a post-MD analysis tool for exploring protein–ligand interactions. MM/GBSA free energy calculations of the best generated complexes were conducted along 100 frames of 500 ns of MDs. ADME descriptors and pharmacokinetic properties of the promising compounds were predicted by means of SwissADME tool [61].

4. Conclusions

In the present computational investigation, the pivotal interactions occurring at the interface of gD and its human surface receptors were carefully explored. Since the availability of X-ray structures of gD in complex with HVEM and Nectin-1, this study aimed at highlighting the key residues involved in the binding interface during MDs. Our reasoning was based on the idea that molecular recognition entails a two-way influence between the interacting partners, whereby the flexibility of gD and its human receptors determined the optimal conformation for the complexes formation. This aspect was better appreciated analyzing the behavior of the complexes after MDs, and the resulting dynamic adaptation was used to define the principal residues responsible for forming a stable complex. MM/GBSA analysis revealed a greater binding affinity of gD towards Nectin-1 with respect to HVEM, given that the rearrangement of N-terminal hairpin is not necessary for gD-Nectin-1 interaction, which instead ensured a more rigid gD-binding pocket. By analyzing HSV-1 PPIs, through the GBPM method, gD N-terminus residues at position 10, 12, 13, 14, 15, 26, 27, 29, 33, and 35 were found pivotal for gD-HVEM binding, whereas nine residues, at position 26, 27, 38, 215, 218, 219, 220, 221, and 222 were relevant for the connection with Nectin-1. Afterwards, an SBVS targeting different conformations of HSV-1 gD binding interfaces was applied. Developing PPI inhibitors is challenging [62], owing to issues such as the general lack of small-molecule starting points for drug design, the typical flatness of the interface, the difficulty of distinguishing real from artefact binding, and the size and character of traditional libraries. In silico approaches are a consolidated strategy to speed up the drug discovery process, as demonstrated during the pandemic emergency, increasing our understanding of how biological systems work and translating this knowledge into new molecules with interesting therapeutic potential. In this context, the inhibition of HSV-1 through targeting viral gD protein represents a good way to act with virus adsorption and membrane fusion. Herein we identified four triazole-pyridines as potential HSV-1 gD inhibitors, with a favourable pharmacokinetic profile and a good theoretical affinity towards all conformations of HSV-1 gD. The surface plasmon resonance technique could provide further information into the mechanism of action at a molecular level and could move forward novel antiviral compounds. Our study suggests a promising basis for the design of a new generation of anti-HSV-1 drugs targeting gD-receptor interfaces.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24087092/s1.

Author Contributions

Conceptualization, I.R., A.M. and R.B.; methodology, I.R. and R.B.; software, S.A.; validation, I.R., R.B., M.B.; formal analysis, R.B.; investigation, I.R. and A.M.; resources, S.A.; data curation, M.B.; writing—original draft preparation, I.R., R.B.; writing—review and editing, M.B., A.M.; visualization, S.A., P.B.; supervision, A.M.; project administration, S.A., P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Additional results found in this work are presented in Supplementary Materials, which will be attached to this document.

Acknowledgments

The authors acknowledge MIUR-PON R&I 2014-2020, Asse I “Investimenti in Capitale Umano”—Azione I “Dottorati innovativi con caratterizzazione industriale”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gilbert, C.; Bestman-Smith, J.; Boivin, G. Resistance of herpesviruses to antiviral drugs: Clinical impacts and molecular mechanisms. Drug Resist. Updat. 2002, 5, 88–114. [Google Scholar] [CrossRef] [PubMed]
  2. De Mello, C.P.P.; Bloom, D.C.; Paixao, I.C. Herpes simplex virus type-1: Replication, latency, reactivation and its antiviral targets. Antivir. Ther. 2016, 21, 277–286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gupta, R.; Warren, T.; Wald, A. Genital herpes. Lancet 2007, 370, 2127–2137. [Google Scholar] [CrossRef]
  4. Wu, J.; Power, H.; Miranda-Saksena, M.; Valtchev, P.; Schindeler, A.; Cunningham, A.L.; Dehghani, F. Identifying HSV-1 Inhibitors from Natural Compounds via Virtual Screening Targeting Surface Glycoprotein D. Pharmaceuticals 2022, 15, 361. [Google Scholar] [CrossRef] [PubMed]
  5. Whitley, R.; Baines, J. Clinical management of herpes simplex virus infections: Past, present, and future. F1000Research 2018, 7, F1000 Faculty Rev-1726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sadowski, L.A.; Upadhyay, R.; Greeley, Z.W.; Margulies, B.J. Current drugs to treat infections with herpes simplex viruses-1 and-2. Viruses 2021, 13, 1228. [Google Scholar] [CrossRef] [PubMed]
  7. Jiang, Y.-C.; Feng, H.; Lin, Y.-C.; Guo, X.-R. New strategies against drug resistance to herpes simplex virus. Int. J. Oral Sci. 2016, 8, 1–6. [Google Scholar] [CrossRef]
  8. Connolly, S.A.; Jackson, J.O.; Jardetzky, T.S.; Longnecker, R. Fusing structure and function: A structural view of the herpesvirus entry machinery. Nat. Rev. Microbiol. 2011, 9, 369–381. [Google Scholar] [CrossRef]
  9. Spear, P.G. Herpes simplex virus: Receptors and ligands for cell entry. Cell. Microbiol. 2004, 6, 401–410. [Google Scholar] [CrossRef]
  10. Heldwein, E.; Krummenacher, C. Entry of herpesviruses into mammalian cells. Cell. Mol. Life Sci. 2008, 65, 1653–1668. [Google Scholar] [CrossRef]
  11. Connolly, S.A.; Jardetzky, T.S.; Longnecker, R. The structural basis of herpesvirus entry. Nat. Rev. Microbiol. 2021, 19, 110–121. [Google Scholar] [CrossRef] [PubMed]
  12. Petermann, P.; Thier, K.; Rahn, E.; Rixon, F.J.; Bloch, W.; Özcelik, S.; Krummenacher, C.; Barron, M.J.; Dixon, M.J.; Scheu, S. Entry mechanisms of herpes simplex virus 1 into murine epidermis: Involvement of nectin-1 and herpesvirus entry mediator as cellular receptors. J. Virol. 2015, 89, 262–274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tiwari, V.; O’Donnell, C.D.; Oh, M.-J.; Valyi-Nagy, T.; Shukla, D. A role for 3-O-sulfotransferase isoform-4 in assisting HSV-1 entry and spread. Biochem. Biophys. Res. Commun. 2005, 338, 930–937. [Google Scholar] [CrossRef] [PubMed]
  14. Labbozzetta, M.; Barreca, M.; Spanò, V.; Raimondi, M.V.; Poma, P.; Notarbartolo, M.; Barraja, P.; Montalbano, A. Novel insights on [1, 2] oxazolo [5, 4-e] isoindoles on multidrug resistant acute myeloid leukemia cell line. Drug Dev. Res. 2022, 83, 1331–1341. [Google Scholar] [CrossRef] [PubMed]
  15. Barreca, M.; Spanò, V.; Rocca, R.; Bivacqua, R.; Abel, A.-C.; Maruca, A.; Montalbano, A.; Raimondi, M.V.; Tarantelli, C.; Gaudio, E. Development of [1, 2] oxazoloisoindoles tubulin polymerization inhibitors: Further chemical modifications and potential therapeutic effects against lymphomas. Eur. J. Med. Chem. 2022, 243, 114744. [Google Scholar] [CrossRef]
  16. Barreca, M.; Spanò, V.; Raimondi, M.V.; Tarantelli, C.; Spriano, F.; Bertoni, F.; Barraja, P.; Montalbano, A. Recurrence of the oxazole motif in tubulin colchicine site inhibitors with anti-tumor activity. Eur. J. Med. Chem. Rep. 2021, 1, 100004. [Google Scholar] [CrossRef]
  17. Barreca, M.; Ingarra, A.M.; Raimondi, M.V.; Spanò, V.; De Franco, M.; Menilli, L.; Gandin, V.; Miolo, G.; Barraja, P.; Montalbano, A. Insight on pyrimido [5, 4-g] indolizine and pyrimido [4, 5-c] pyrrolo [1, 2-a] azepine systems as promising photosensitizers on malignant cells. Eur. J. Med. Chem. 2022, 237, 114399. [Google Scholar] [CrossRef]
  18. Frasson, I.; Spano, V.; Di Martino, S.; Nadai, M.; Doria, F.; Parrino, B.; Carbone, A.; Cascioferro, S.M.; Diana, P.; Cirrincione, G. Synthesis and photocytotoxic activity of [1, 2, 3] triazolo [4, 5-h][1, 6] naphthyridines and [1, 3] oxazolo [5, 4-h][1, 6] naphthyridines. Eur. J. Med. Chem. 2019, 162, 176–193. [Google Scholar] [CrossRef]
  19. Seck, I.; Nguemo, F. Triazole, imidazole, and thiazole-based compounds as potential agents against coronavirus. Results Chem. 2021, 3, 100132. [Google Scholar] [CrossRef]
  20. Dheer, D.; Singh, V.; Shankar, R. Medicinal attributes of 1, 2, 3-triazoles: Current developments. Bioorganic Chem. 2017, 71, 30–54. [Google Scholar] [CrossRef]
  21. Bozorov, K.; Zhao, J.; Aisa, H.A. 1, 2, 3-Triazole-containing hybrids as leads in medicinal chemistry: A recent overview. Bioorganic Med. Chem. 2019, 27, 3511–3531. [Google Scholar] [CrossRef] [PubMed]
  22. Bivacqua, R.; Barreca, M.; Spanò, V.; Raimondi, M.V.; Romeo, I.; Alcaro, S.; Andrei, G.; Barraja, P.; Montalbano, A. Insight into non-nucleoside triazole-based systems as viral polymerases inhibitors. Eur. J. Med. Chem. 2023, 115136. [Google Scholar] [CrossRef] [PubMed]
  23. Zhuang, L.; Wai, J.S.; Embrey, M.W.; Fisher, T.E.; Egbertson, M.S.; Payne, L.S.; Guare, J.P.; Vacca, J.P.; Hazuda, D.J.; Felock, P.J. Design and synthesis of 8-hydroxy-[1, 6] naphthyridines as novel inhibitors of HIV-1 integrase in vitro and in infected cells. J. Med. Chem. 2003, 46, 453–456. [Google Scholar] [CrossRef] [PubMed]
  24. Falardeau, G.; Lachance, H.; St-Pierre, A.; Yannopoulos, C.G.; Drouin, M.; Bédard, J.; Chan, L. Design and synthesis of a potent macrocyclic 1, 6-napthyridine anti-human cytomegalovirus (HCMV) inhibitors. Bioorganic Med. Chem. Lett. 2005, 15, 1693–1695. [Google Scholar] [CrossRef] [PubMed]
  25. Embrey, M.W.; Wai, J.S.; Funk, T.W.; Homnick, C.F.; Perlow, D.S.; Young, S.D.; Vacca, J.P.; Hazuda, D.J.; Felock, P.J.; Stillmock, K.A. A series of 5-(5, 6)-dihydrouracil substituted 8-hydroxy-[1, 6] naphthyridine-7-carboxylic acid 4-fluorobenzylamide inhibitors of HIV-1 integrase and viral replication in cells. Bioorganic Med. Chem. Lett. 2005, 15, 4550–4554. [Google Scholar] [CrossRef]
  26. Chan, L.; Stefanac, T.; Lavallée, J.-F.; Jin, H.; Bédard, J.; May, S.; Falardeau, G. Design and synthesis of new potent human cytomegalovirus (HCMV) inhibitors based on internally hydrogen-bonded 1, 6-naphthyridines. Bioorganic Med. Chem. Lett. 2001, 11, 103–105. [Google Scholar] [CrossRef]
  27. Chan, L.; Jin, H.; Stefanac, T.; Lavallée, J.-F.; Falardeau, G.; Wang, W. Discovery of 1, 6-naphthyridines as a novel class of potent and selective human cytomegalovirus inhibitors. J. Med. Chem. 1999, 42, 3023–3025. [Google Scholar] [CrossRef]
  28. Jordão, A.K.; Ferreira, V.F.; Souza, T.M.; de Souza Faria, G.G.; Machado, V.; Abrantes, J.L.; de Souza, M.C.; Cunha, A.C. Synthesis and anti-HSV-1 activity of new 1, 2, 3-triazole derivatives. Bioorganic Med. Chem. 2011, 19, 1860–1865. [Google Scholar] [CrossRef] [Green Version]
  29. Głowacka, I.E.; Balzarini, J.; Wróblewski, A.E. The synthesis, antiviral, cytostatic and cytotoxic evaluation of a new series of acyclonucleotide analogues with a 1, 2, 3-triazole linker. Eur. J. Med. Chem. 2013, 70, 703–722. [Google Scholar] [CrossRef]
  30. Głowacka, I.E.; Andrei, G.; Schols, D.; Snoeck, R.; Gawron, K. Design, Synthesis, and the Biological Evaluation of a New Series of Acyclic 1,2,3-Triazole Nucleosides. Arch. Der Pharm. 2017, 350, 1700166. [Google Scholar] [CrossRef] [Green Version]
  31. Cunha, A.C.; Ferreira, V.F.; Vaz, M.G.; Cassaro, R.A.A.; Resende, J.A.; Sacramento, C.Q.; Costa, J.; Abrantes, J.L.; Souza, T.M.L.; Jordão, A.K. Chemistry and anti-herpes simplex virus type 1 evaluation of 4-substituted-1 H-1, 2, 3-triazole-nitroxyl-linked hybrids. Mol. Divers. 2021, 25, 2035–2043. [Google Scholar] [CrossRef] [PubMed]
  32. Bernardino, A.M.; Castro, H.C.; Frugulhetti, I.C.; Loureiro, N.I.; Azevedo, A.R.; Pinheiro, L.C.; Souza, T.M.; Giongo, V.; Passamani, F.; Magalhaes, U.O. SAR of a series of anti-HSV-1 acridone derivatives, and a rational acridone-based design of a new anti-HSV-1 3H-benzo [b] pyrazolo [3, 4-h]-1, 6-naphthyridine series. Bioorganic Med. Chem. 2008, 16, 313–321. [Google Scholar] [CrossRef] [PubMed]
  33. Parcella, K.; Patel, M.; Tu, Y.; Eastman, K.; Peese, K.; Gillis, E.; Belema, M.; Dicker, I.B.; McAuliffe, B.; Ding, B. Scaffold modifications to the 4-(4, 4-dimethylpiperidinyl) 2, 6-dimethylpyridinyl class of HIV-1 allosteric integrase inhibitors. Bioorganic Med. Chem. 2022, 116833. [Google Scholar] [CrossRef]
  34. Karypidou, K.; Ribone, S.R.; Quevedo, M.A.; Persoons, L.; Pannecouque, C.; Helsen, C.; Claessens, F.; Dehaen, W. Synthesis, biological evaluation and molecular modeling of a novel series of fused 1, 2, 3-triazoles as potential anti-coronavirus agents. Bioorganic Med. Chem. Lett. 2018, 28, 3472–3476. [Google Scholar] [CrossRef] [PubMed]
  35. Hartwich, A.; Zdzienicka, N.; Schols, D.; Andrei, G.; Snoeck, R.; Głowacka, I.E. Design, synthesis and antiviral evaluation of novel acyclic phosphonate nucleotide analogs with triazolo [4, 5-b] pyridine, imidazo [4, 5-b] pyridine and imidazo [4, 5-b] pyridin-2 (3 H)-one systems. Nucl. Nucl. Nucleic Acids 2020, 39, 542–591. [Google Scholar] [CrossRef] [Green Version]
  36. Van Zundert, G.; Rodrigues, J.; Trellet, M.; Schmitz, C.; Kastritis, P.; Karaca, E.; Melquiond, A.; van Dijk, M.; De Vries, S.; Bonvin, A. The HADDOCK2. 2 web server: User-friendly integrative modeling of biomolecular complexes. J. Mol. Biol. 2016, 428, 720–725. [Google Scholar] [CrossRef] [Green Version]
  37. Jorgensen, W.L.; Maxwell, D.S.; Tirado-Rives, J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. [Google Scholar] [CrossRef]
  38. Ortuso, F.; Langer, T.; Alcaro, S. GBPM: GRID-based pharmacophore model: Concept and application studies to protein–protein recognition. Bioinformatics 2006, 22, 1449–1455. [Google Scholar] [CrossRef] [Green Version]
  39. Ongaro, A.; Oselladore, E.; Memo, M.; Ribaudo, G.; Gianoncelli, A. Insight into the LFA-1/SARS-CoV-2 Orf7a complex by protein–protein docking, molecular dynamics, and MM-GBSA calculations. J. Chem. Inf. Model. 2021, 61, 2780–2787. [Google Scholar] [CrossRef]
  40. Eisenberg, R.; Atanasiu, D.; Cairns, T.; Gallagher, J.; Krummenacher, C.; Cohen, G. Herpes virus fusion and entry: A story with many characters. Viruses 2012, 4, 800–832. [Google Scholar] [CrossRef]
  41. Zhang, N.; Yan, J.; Lu, G.; Guo, Z.; Fan, Z.; Wang, J.; Shi, Y.; Qi, J.; Gao, G.F. Binding of herpes simplex virus glycoprotein D to nectin-1 exploits host cell adhesion. Nat. Commun. 2011, 2, 577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Lazear, E.; Whitbeck, J.C.; Zuo, Y.; Carfí, A.; Cohen, G.H.; Eisenberg, R.J.; Krummenacher, C. Induction of conformational changes at the N-terminus of herpes simplex virus glycoprotein D upon binding to HVEM and nectin-1. Virology 2014, 448, 185–195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Prejanò, M.; Romeo, I.; La Serra, M.A.; Russo, N.; Marino, T. Computational Study Reveals the Role of Water Molecules in the Inhibition Mechanism of LAT1 by 1, 2, 3-Dithiazoles. J. Chem. Inf. Model. 2021, 61, 5883–5892. [Google Scholar] [CrossRef] [PubMed]
  44. Schrödinger Release 2018-1: Desmond Molecular Dynamics System; D.E. Shaw Research, Maestro-Desmond Interoperability Tools; Schrödinger LLC: New York, NY, USA, 2018.
  45. Schrödinger Release 2018-1: Maestro; Schrödinger LLC: New York, NY, USA, 2018.
  46. Carfí, A.; Willis, S.H.; Whitbeck, J.C.; Krummenacher, C.; Cohen, G.H.; Eisenberg, R.J.; Wiley, D.C. Herpes simplex virus glycoprotein D bound to the human receptor HveA. Mol. Cell 2001, 8, 169–179. [Google Scholar] [CrossRef]
  47. Schrödinger Release 2018-1: Protein Preparation Wizard; Epik, 2019; Schrödinger LLC: New York, NY, USA, 2018.
  48. Jacobson, M.P.; Pincus, D.L.; Rapp, C.S.; Day, T.J.; Honig, B.; Shaw, D.E.; Friesner, R.A. A hierarchical approach to all-atom protein loop prediction. Proteins: Struct. Funct. Bioinform. 2004, 55, 351–367. [Google Scholar] [CrossRef] [Green Version]
  49. Greenwood, J.R.; Calkins, D.; Sullivan, A.P.; Shelley, J.C. Towards the comprehensive, rapid, and accurate prediction of the favorable tautomeric states of drug-like molecules in aqueous solution. J. Comput.-Aided Mol. Des. 2010, 24, 591–604. [Google Scholar] [CrossRef]
  50. Ortuso, F.; Mercatelli, D.; Guzzi, P.H.; Giorgi, F.M. Structural genetics of circulating variants affecting the SARS-CoV-2 spike/human ACE2 complex. J. Biomol. Struct. Dyn. 2022, 40, 6545–6555. [Google Scholar] [CrossRef]
  51. Hou, T.; Wang, J.; Li, Y.; Wang, W. Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking. J. Comput. Chem. 2011, 32, 866–877. [Google Scholar] [CrossRef] [Green Version]
  52. Kopitz, H.; Cashman, D.A.; Pfeiffer-Marek, S.; Gohlke, H. Influence of the solvent representation on vibrational entropy calculations: Generalized born versus distance-dependent dielectric model. J. Comput. Chem. 2012, 33, 1004–1013. [Google Scholar] [CrossRef]
  53. Genheden, S.; Kuhn, O.; Mikulskis, P.; Hoffmann, D.; Ryde, U. The normal-mode entropy in the MM/GBSA method: Effect of system truncation, buffer region, and dielectric constant. J. Chem. Inf. Model. 2012, 52, 2079–2088. [Google Scholar] [CrossRef] [Green Version]
  54. Krummenacher, C.; Supekar, V.M.; Whitbeck, J.C.; Lazear, E.; Connolly, S.A.; Eisenberg, R.J.; Cohen, G.H.; Wiley, D.C.; Carfí, A. Structure of unliganded HSV gD reveals a mechanism for receptor-mediated activation of virus entry. EMBO J. 2005, 24, 4144–4153. [Google Scholar] [CrossRef]
  55. Schrödinger Release 2018-1: LigPrep; Schrödinger LLC: New York, NY, USA, 2018.
  56. Prejanò, M.; Romeo, I.; Sgrizzi, L.; Russo, N.; Marino, T. Why hydroxy-proline improves the catalytic power of the peptidoglycan N-deacetylase enzyme: Insight from theory. Phys. Chem. Chem. Phys. 2019, 21, 23338–23345. [Google Scholar] [CrossRef] [PubMed]
  57. Maruca, A.; Ambrosio, F.A.; Lupia, A.; Romeo, I.; Rocca, R.; Moraca, F.; Talarico, C.; Bagetta, D.; Catalano, R.; Costa, G. Computer-based techniques for lead identification and optimization I: Basics. Phys. Sci. Rev. 2019, 4. [Google Scholar] [CrossRef]
  58. Lupia, A.; Moraca, F.; Bagetta, D.; Maruca, A.; Ambrosio, F.A.; Rocca, R.; Catalano, R.; Romeo, I.; Talarico, C.; Ortuso, F. Computer-based techniques for lead identification and optimization II: Advanced search methods. Phys. Sci. Rev. 2019, 5, 20180114. [Google Scholar] [CrossRef]
  59. Costa, G.; Rocca, R.; Moraca, F.; Talarico, C.; Romeo, I.; Ortuso, F.; Alcaro, S.; Artese, A. A Comparative Docking Strategy to Identify Polyphenolic Derivatives as Promising Antineoplastic Binders of G-quadruplex DNA c-myc and bcl-2 Sequences. Mol. Inform. 2016, 35, 391–402. [Google Scholar] [CrossRef]
  60. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein-Ligand Complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [Green Version]
  61. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [Green Version]
  62. Ferraro, M.; Colombo, G. Targeting difficult protein-protein interactions with plain and general computational approaches. Molecules 2018, 23, 2256. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Graphical workflow of the applied in silico approaches on HSV-1 gD Pockets.
Figure 1. Graphical workflow of the applied in silico approaches on HSV-1 gD Pockets.
Ijms 24 07092 g001
Figure 2. Surface representation of (A) gD-HVEM and (B) gD-Nectin-1 complexes and focus on the interface with key interactions labelled and displayed in carbon sticks.
Figure 2. Surface representation of (A) gD-HVEM and (B) gD-Nectin-1 complexes and focus on the interface with key interactions labelled and displayed in carbon sticks.
Ijms 24 07092 g002
Figure 3. Plots of the distances between the residues of gD involved in the interaction at the interface with HVEM after 100 ns of MDs.
Figure 3. Plots of the distances between the residues of gD involved in the interaction at the interface with HVEM after 100 ns of MDs.
Ijms 24 07092 g003
Figure 4. Plots of the distances between the residues of gD involved in the interaction at the interface with Nectin-1 after 100 ns of MDs.
Figure 4. Plots of the distances between the residues of gD involved in the interaction at the interface with Nectin-1 after 100 ns of MDs.
Ijms 24 07092 g004
Figure 5. Plot of MM/GBSA trend for HVEM and Nectin-1 in complex to gD during 100 ns of MDs.
Figure 5. Plot of MM/GBSA trend for HVEM and Nectin-1 in complex to gD during 100 ns of MDs.
Ijms 24 07092 g005
Figure 6. 3D representation of (A) 14, (B) 16, (C) 83, and (D) 85 in complex to the most populated cluster of gD-Pocket 1. gD-Pocket 1 is illustrated in gray, with the residues involved in pivotal contacts shown as carbon sticks. 14, 16, 83, and 85 are depicted as pink, green, violet, and cyan carbon sticks, and H-bonds are indicated as yellow dashed lines.
Figure 6. 3D representation of (A) 14, (B) 16, (C) 83, and (D) 85 in complex to the most populated cluster of gD-Pocket 1. gD-Pocket 1 is illustrated in gray, with the residues involved in pivotal contacts shown as carbon sticks. 14, 16, 83, and 85 are depicted as pink, green, violet, and cyan carbon sticks, and H-bonds are indicated as yellow dashed lines.
Ijms 24 07092 g006
Figure 7. 3D representation of (A) 14, (B) 16, (C) 83, and (D) 85 in complex to the most populated cluster of gD-Pocket 2. gD-Pocket 2 is illustrated in salmon, with the residues involved in pivotal contacts shown as carbon sticks. Compounds 14, 16, 83, and 85 are depicted as pink, green, violet, and cyan carbon sticks, whereas H-bonds and π-π interactions are indicated as yellow and cyan dashed lines, respectively.
Figure 7. 3D representation of (A) 14, (B) 16, (C) 83, and (D) 85 in complex to the most populated cluster of gD-Pocket 2. gD-Pocket 2 is illustrated in salmon, with the residues involved in pivotal contacts shown as carbon sticks. Compounds 14, 16, 83, and 85 are depicted as pink, green, violet, and cyan carbon sticks, whereas H-bonds and π-π interactions are indicated as yellow and cyan dashed lines, respectively.
Ijms 24 07092 g007
Figure 8. RMSD plots of 14, 16, 83, and 85 compounds in complex with (A) gD-Pocket 1 and (B) gD-Pocket 2, calculated on protein’s backbone atoms during 500 ns of MDs.
Figure 8. RMSD plots of 14, 16, 83, and 85 compounds in complex with (A) gD-Pocket 1 and (B) gD-Pocket 2, calculated on protein’s backbone atoms during 500 ns of MDs.
Ijms 24 07092 g008
Table 1. Synthetic (1–12) and designed (19–81) triazolo [4,5-h][1,6]naphthyridines and synthetic (13–18) and designed (82–85) triazolo [4,5-b] pyridines.
Table 1. Synthetic (1–12) and designed (19–81) triazolo [4,5-h][1,6]naphthyridines and synthetic (13–18) and designed (82–85) triazolo [4,5-b] pyridines.
Ijms 24 07092 i001
CpdRR1CpdRR1
1HMe103,4,5-(OMe)3Ph
2HPh113,4,5-(OMe)32-OMe-Ph
3H2-OMe-Ph123,4,5-(OMe)34-OMe-Ph
4H4-OMe-Ph13HCOOEt
54-OMeMe14HCH2OH
64-OMePh154-OMeCOOEt
74-OMe2-OMe-Ph164-OMeCH2OH
84-OMe4-OMe-Ph173,4,5-(OMe)3COOEt
93,4,5-(OMe)3Me183,4,5-(OMe)3CH2OH
19H2-F-Ph554-Cl2-OMe-Ph
20H2-Br-Ph564-Cl2-F-Ph
21H2-Cl-Ph574-Cl2-Br-Ph
22H2-NO2-Ph584-Cl2-Cl-Ph
23H2-NH2-Ph594-Cl2-NO2-Ph
24H2-CF3-Ph604-Cl2-NH2-Ph
25H2,3-(OMe)2-Ph614-Cl2-CF3-Ph
26H2,4-(OMe)2-Ph624-Cl2,3-(OMe)2-Ph
27H2,5-(OMe)2-Ph634-Cl2,4-(OMe)2-Ph
28H2-OMe-5-Cl-Ph644-Cl2,5-(OMe)2-Ph
29H2-OMe-5-Br-Ph654-Cl2-OMe-5-Cl-Ph
304-OMe2-F-Ph664-Cl2-OMe-5-Br-Ph
314-OMe2-Br-Ph674-F4-OMe-Ph
324-OMe2-Cl-Ph684-FCH3
334-OMe2-NO2-Ph694-FPh
344-OMe2-NH2-Ph704-F2-OMe-Ph
354-OMe2-CF3-Ph714-F2-F-Ph
364-OMe2,3-(OMe)2-Ph724-F2-Br-Ph
374-OMe2,4-(OMe)2-Ph734-F2-Cl-Ph
384-OMe2,5-(OMe)2-Ph744-F2-NO2-Ph
394-OMe2-OMe-5-Cl-Ph754-F2-NH2-Ph
404-OMe2-OMe-5-Br-Ph764-F2-CF3-Ph
413,4,5-(OMe)32-F-Ph774-F2,3-(OMe)2-Ph
423,4,5-(OMe)32-Br-Ph784-F2,4-(OMe)2-Ph
433,4,5-(OMe)32-Cl-Ph794-F2,5-(OMe)2-Ph
443,4,5-(OMe)32-NO2-Ph804-F2-OMe-5-Cl-Ph
453,4,5-(OMe)32-NH2-Ph814-F2-OMe-5-Br-Ph
463,4,5-(OMe)32-CF3-Ph824-ClCOOEt
473,4,5-(OMe)32,3-(OMe)2-Ph834-ClCH2OH
483,4,5-(OMe)32,4-(OMe)2-Ph844-FCOOEt
493,4,5-(OMe)32,5-(OMe)2-Ph854-FCH2OH
503,4,5-(OMe)32-OMe-5-Cl-Ph
513,4,5-(OMe)32-OMe-5-Br-Ph
524-Cl4-OMe-Ph
534-ClCH3
544-ClPh
Table 2. Number of generated clusters, HADDOCK scores, cluster size, Z-score, and BSA values of gD-HVEM and gD-Nectin-1 complexes.
Table 2. Number of generated clusters, HADDOCK scores, cluster size, Z-score, and BSA values of gD-HVEM and gD-Nectin-1 complexes.
Complex ClusterHADDOCK ScoreCluster SizeZ-ScoreBSA
gD-HVEM1−124.3 ± 5.696−2.12144.2 ± 59.6
2−89.1± 1.924−0.21939.7 ± 69.3
3−97.5 ± 7.717−0.71925.6 ± 72.1
4−93.7 ± 2.115−0.51965.3 ± 124.9
5−54.1 ± 3.3101.61401.2 ± 37.0
6−82.2 ± 4.490.11543.9 ± 50.6
7−76.3 ± 3.770.51864.7 ± 63.3
8−70.1 ± 7.660.81652.0 ± 72.3
9−75.7 ± 20.740.51701.5 ± 193.4
gD-Nectin-11−146.0 ± 2.2178−1.02107.6 ± 20.8
2−101.4 ± 3.0201.01887 ± 108.6
Table 3. Residues involved in the interaction of gD with HVEM and Nectin-1.
Table 3. Residues involved in the interaction of gD with HVEM and Nectin-1.
gD-HVEMgD-Nectin-1
gDHVEMgDNectin-1
Hydrogen bondsLys10Asp7Asp26Lys61
Ala12Tyr23Gln27Lys61
Pro14Arg75Gln27Thr63
Asn15Ser74Gln132Gln64
Gly19Arg75Asp215Asn77
Gln27Cys37Asp215Asn77
Thr29Thr35Gly218Gln68
Lys122Glu31Arg222Gly132
Salt bridgesLys10Asp7Arg222Glu125
Asp26Lys26Asp26Lys61
Lys122Glu31Arg222Glu125
Table 4. GBPM average scores and quartile distribution of the pivotal hotspots of gD in complex to HVEM and Nectin-1, for all frames of MDs.
Table 4. GBPM average scores and quartile distribution of the pivotal hotspots of gD in complex to HVEM and Nectin-1, for all frames of MDs.
gD-HVEMgD-Nectin-1
ResidueAverage ScoreQuartileResidueAverage ScoreQuartile
Ala12−6.75Q1Arg222−1.80Q1
Asp13−2.08Q1Asp215−3.74Q1
Asn15−3.13Q1Asp26−1.01Q1
Asp26−1.97Q1Gln27−1.13Q1
Gln27−5.48Q1Gly218−0.97Q1
Thr29−10.00Q1Met219−1.28Q1
Gly33−1.21Q1Pro221−1.83Q1
Arg35−14.64Q1Ser200−1.85Q1
Lys10−0.21Q2Tyr38−1.09Q1
Lys122−0.47Q2Arg196−0.20Q2
Pro14−0.30Q2Arg64−0.83Q2
Lys186−0.20Q2Gln132−0.40Q2
Gly19−0.07Q2Leu25−0.44Q2
Leu257−0.12Q2Lys186−0.66Q2
Arg64−0.20Q2Pro23−0.23Q2
Ala7−0.67Q2Ser235−0.64Q2
Phe17−0.01Q3Tyr137−0.17Q2
Lys20−0.02Q3Arg184−0.04Q3
Leu25−0.03Q3Asn136−0.03Q3
Glu259−0.02Q3Asp139−0.16Q3
Asp30−0.06Q3Pro199−0.03Q3
Pro32−0.04Q3Ser140−0.04Q3
Ala5−0.04Q3Ser216−0.08Q3
Glu63−0.02Q3Tyr234−0.15Q3
Met11−0.01Q4Val24−0.05Q3
Arg18−0.01Q4Ala185−0.01Q4
Val240.00Q4Arg36−0.01Q4
Lys245−0.01Q4Asn2270.00Q4
Ser258−0.01Q4Ile217−0.02Q4
Pro310.00Q4Lys190−0.02Q4
Val34−0.01Q4Phe2230.00Q4
Leu40.00Q4Thr230−0.01Q4
Val231−0.02Q4
Val37−0.01Q4
Table 5. Predicted ADME parameters: Molecular Weight (MW); Number of H-bond Acceptors and Donors; Octanol/Water Partition Coefficient (logP); Topological Polar Surface Area in A2 (TPSA); and Water Solubility (LogS).
Table 5. Predicted ADME parameters: Molecular Weight (MW); Number of H-bond Acceptors and Donors; Octanol/Water Partition Coefficient (logP); Topological Polar Surface Area in A2 (TPSA); and Water Solubility (LogS).
CpdMWH-Bond
Acceptors
H-Bond
Donors
LogPTPSALogS
14255.28421.0589.85−2.77
16285.30521.0699.08−2.82
83289.72421.7189.85−3.35
85273.27521.6189.85−2.91
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bivacqua, R.; Romeo, I.; Barreca, M.; Barraja, P.; Alcaro, S.; Montalbano, A. HSV-1 Glycoprotein D and Its Surface Receptors: Evaluation of Protein–Protein Interaction and Targeting by Triazole-Based Compounds through In Silico Approaches. Int. J. Mol. Sci. 2023, 24, 7092. https://doi.org/10.3390/ijms24087092

AMA Style

Bivacqua R, Romeo I, Barreca M, Barraja P, Alcaro S, Montalbano A. HSV-1 Glycoprotein D and Its Surface Receptors: Evaluation of Protein–Protein Interaction and Targeting by Triazole-Based Compounds through In Silico Approaches. International Journal of Molecular Sciences. 2023; 24(8):7092. https://doi.org/10.3390/ijms24087092

Chicago/Turabian Style

Bivacqua, Roberta, Isabella Romeo, Marilia Barreca, Paola Barraja, Stefano Alcaro, and Alessandra Montalbano. 2023. "HSV-1 Glycoprotein D and Its Surface Receptors: Evaluation of Protein–Protein Interaction and Targeting by Triazole-Based Compounds through In Silico Approaches" International Journal of Molecular Sciences 24, no. 8: 7092. https://doi.org/10.3390/ijms24087092

APA Style

Bivacqua, R., Romeo, I., Barreca, M., Barraja, P., Alcaro, S., & Montalbano, A. (2023). HSV-1 Glycoprotein D and Its Surface Receptors: Evaluation of Protein–Protein Interaction and Targeting by Triazole-Based Compounds through In Silico Approaches. International Journal of Molecular Sciences, 24(8), 7092. https://doi.org/10.3390/ijms24087092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop