key: cord-0715731-xed3w47t
authors: Suzuki, Yoshiyuki
title: Predicting receptor functionality of signaling lymphocyte activation molecule for measles virus hemagglutinin by docking simulation
date: 2017-05-17
journal: Microbiol Immunol
DOI: 10.1111/1348-0421.12484
sha: 5eeca2a28e7dc3e0a4a65f80dc2d331e6c78c1ba
doc_id: 715731
cord_uid: xed3w47t

Predicting susceptibility of various species to a virus assists assessment of risk of interspecies transmission. Evaluation of receptor functionality may be useful in screening for susceptibility. In this study, docking simulation was conducted for measles virus hemagglutinin (MV‐H) and immunoglobulin‐like variable domain of signaling lymphocyte activation molecule (SLAM‐V). It was observed that the docking scores for MV‐H and SLAM‐V correlated with the activity of SLAM as an MV receptor. These results suggest that the receptor functionality may be predicted from the docking scores of virion surface proteins and cellular receptor molecules.

for establishment of susceptibility, evaluation of the former may facilitate screening for the latter.

Receptor functionality may reflect the binding affinity of virion surface proteins and cellular receptor molecules. Binding affinity may be evaluated by using docking simulation analyzing 3-D structures, binding potential being measured as a docking score (9) . Docking scores for DPP4 and MERS-CoV S protein reportedly correlate with the activity of DPP4 as a MERS-CoV receptor in mammals (6) . It is therefore of interest to determine whether this relationship is also applicable to other viruses.

Measles virus is a member of the genus Morbillivirus in the family Paramyxoviridae in the order Mononegavirales (10) . SLAM is known to be the principal cellular entry receptor for MV (11) . The 3-D structure of the binding complex of MV-H and SLAM has been resolved (12) . In addition, SLAM's activity as an MV receptor has been investigated by transfection experiments (11) (12) (13) (14) (15) (16) (17) . The purpose of the present study was to investigate the relationship between docking scores and activity of SLAM as an MV receptor by conducting docking simulation of MV-H and SLAM.

A type I transmembrane protein, SLAM belongs to the Ig superfamily (18) . The Ig-like variable domain in the ectodomain of SLAM (SLAM-V) is responsible for binding to MV-H (14) . The activity of SLAM as MV receptor has been investigated by transfection experiments for human (Homo sapiens), dog (Canis familaris), cow (Bos taurus), cotton rat (Sigmodon hispidus), and mouse (Mus musculus). Human is the single natural host of MV; indeed, human SLAM has been shown to act as an MV receptor (11) (12) (13) (14) (15) (16) (17) . Nevertheless, it has also been demonstrated that SLAMs from dog, cow, and cotton rat act as MV receptors, although less efficiently than human SLAM (13, 17) . In contrast, mouse SLAM reportedly does not exhibit MV receptor functionality (14, 15) . SLAM sequences of these species have been retrieved from the International Nucleotide Sequence Database; their accession numbers are listed in Table 1 .

Transfection experiments have also been performed for chimeras of human and mouse SLAMs (15, 16) . in the binding interface as described above, these positions may be critical for binding of SLAM-V and MV-H.

For each of human, dog, cow, cotton rat, and mouse SLAM-Vs, and for human SLAM-V the amino acid positions 60, 61, and 63 of which were mutated to those of mouse SLAM-V (mutated human SLAM-V), and for mouse SLAM-V the corresponding positions of which were mutated to those of human SLAM-V (mutated mouse SLAM-V), pairwise alignment of amino acid sequences was achieved with cotton-top tamarin SLAM-V (International Nucleotide Sequence Database accession number: AF257239) using computer program MAFFT (version 7.305b) (19) . Thence, the 3-D structure was constructed by homology modelling using MOD-ELLER (version 9.17) (20) and using the 3-D structure of cotton-top tamarin SLAM-V as the template (PDB ID: 3ALZ, B chain).

The 3-D structure of SLAM-V obtained as described above and that of MV-H (PDB ID: 3ALZ, A chain) were employed for docking simulation with ClusPro (version 2.0), which is one of the best, fully automated web servers available for protein-protein docking (21) . In ClusPro, missing atoms and polar hydrogens are automatically added to the proteins before docking simulation. Amino acid positions 61 and 63 of SLAM-V and positions 552 and 554 of MV-H were assumed to attract one another. Several candidate structures of binding complexes together with their docking scores were generated for each SLAM-V. The candidate structures that appeared to be the most closely related to the original structure of the binding complex of MV-H and cotton-top tamarin SLAM-V (PDB ID: 3ALZ) (12) was selected as follows. In each candidate structure, physical distances were computed for 55 pairs of amino acid positions that are considered to form the binding interface in the original structure, as described above. On the assumption that interacting amino acid positions are located within 10 Å (22, 23) , the structure in which the greatest number of the 55 pairs of amino acid positions are located within 10 Aẘ as considered to be the most closely related to the original structure. In fact, more than 90% (50 of the 55 pairs) were located within 10 Åin the original structure.

From 44 to 52 of the 55 pairs of amino acid positions were found to be located within 10 Åin the selected structures; these were largely comparable to the original structure (Table 1) . Transfection experiments have shown that SLAM is functional as an MV receptor in human, functional but with less efficiency in dog, cow and cotton rat, and non-functional in mouse (11) (12) (13) (14) (15) (16) (17) . The docking score for MV-H and SLAM-V is consistently lowest (indicating the strongest interaction) in human, intermediate in dog, cow and cotton rat, and highest (indicating the weakest interaction) in mouse ( Table 1 ). The probability of obtaining this order of docking scores ([human] < [dog, cow and cotton rat] < [mouse]) by chance is 0.05, which is of marginal significance.

In additional transfection experiments, recombinant human SLAM containing the amino acid sequence of mutated human SLAM-V did not act as an MV receptor, whereas recombinant mouse SLAM containing the sequence of mutated mouse SLAM-V did (15, 16) . Accordingly, the docking scores for mutated human SLAM-V were found to be higher than those for SLAM-Vs functional as MV receptor, whereas the docking scores for mutated mouse SLAM-V were lower than those for SLAM-Vs non-functional as MV receptor ( Table 1 ). The probability of obtaining this order of docking scores ([human, dog, cow, cotton rat, and mutated mouse] < [mouse and mutated human]) by chance is 0.0476, which is statistically significant. Thus, there is a correlation between docking scores for MV-H and SLAM-V and the activity of SLAM as an MV receptor.

In the above analysis, interacting amino acid positions were assumed to be located within 10 Å (22, 23) . To examine the robustness of the results against violation of this assumption, the above analysis was repeated assuming that interacting amino acid positions were located within 5 Åor 15 Å. In the original structure of the binding complex of MV-H and cotton-top tamarin SLAM-V (PDB ID: 3ALZ) (12), 27 and 55 pairs of the 55 pairs of interacting amino acid positions are located within 5 Åand 15 Å, respectively. When 5 Åwas used as the threshold value, the selected structures sometimes differed from those in Table 1 (dog, mouse and mutated human SLAM-Vs) (Table S1 ). However, the correlation between the docking score of MV-H and SLAM-V and the activity of SLAM as an MV receptor was unaffected. Similarly, different structures were sometimes selected when 15 Aẘ as used as the threshold value (dog, cow and mutated human SLAM-Vs) ( Table S2 ). In this case, however, multiple structures were selected for two SLAM-Vs (cow and mutated human), and in one case (cow) the docking score of the additional structure did not support the correlation. This result was apparently attributable to the threshold value being too large to distinguish different structures. Indeed, all 55 pairs of interacting amino acid positions were found to be located within 15 Åin multiple selected structures. These observations suggest that the results obtained in the present study are robust provided the threshold value assumed for the physical distance between interacting amino acid positions is appropriate (22, 23) .

Previous docking simulation studies of MV-H and SLAM-V have examined the relationship between docking scores and susceptibility to MV in mammals (24, 25) . However, the relationship has remained unclear, probably because susceptibility is determined by multiple infection-related processes. In fact, although dog, cow and cotton rat SLAMs have been shown to act as MV receptors in transfection experiments (13, 17) , these species are not considered to be susceptible to MV (24, 25) . It is possible that infection processes of MV that are irrelevant to receptor functionality are inhibited in these species. Cross-immunity to MV elicited by infection of other morbilliviruses, such as canine distemper virus and peste-des-petits-ruminants virus, may also suppress infection of MV in these species (26) .

In the present study, it was observed that the docking scores of MV-H and SLAM-V correlated with the activity of SLAM as an MV receptor, suggesting that receptor functionality may be predicted from the docking scores of virion surface proteins and cellular receptor molecules. This approach may be useful for screening for susceptibility of various species to a virus when assessing the risk of interspecies transmission (1, 2) . It will be of interest to investigate whether this relationship is further applicable to other viruses (6) . nonhuman primates provides partial protection against infection with canine distemper virus. J Virol 88: 4423-33.

Additional supporting information may be found in the online version of this article at the publisher's web-site. Table S1 . Results from docking simulation of MV-H and immunoglobulin-like SLAM-V assuming that interacting amino acid positions are located within 5 AT able S2. Results from docking simulation of MV-H and immunoglobulin-like SLAM-V assuming that interacting amino acid positions are located within 15 A©

The glycoproteins of all filovirus species use the same host factors for entry into bat and human cells but entry efficiency is species dependent

Dynamics of a morbillivirus at the domestic-wildlife interface: Canine distemper virus in domestic dogs and lions

Virus entry and uncoating

Functional receptor molecules CD300lf and CD300ld within the CD300 family enable murine noroviruses to infect cells

Discovery of a proteinaceous cellular receptor for a norovirus

Host species restriction of Middle East respiratory syndrome coronavirus through its receptor, dipeptidyl peptidase 4

Host range of poliovirus is restricted to simians because of a rapid sequence change of the poliovirus receptor gene during evolution

Ancient positive selection on CD155 as a possible cause for susceptibility to poliovirus infection in simians

How good is automated protein docking?

Measles virus

SLAM (CDw150) is a cellular receptor for measles virus

Structure of the measles virus hemagglutinin bound to its cellular receptor SLAM

Morbilliviruses use signaling lymphocyte activation molecules (CD150) as cellular receptors

V domain of human SLAM (CDw150) is essential for its function as a measles virus receptor

Histidine at position 61 and its adjacent amino acid residues are critical for the ability of SLAM (CD150) to act as a cellular receptor for measles virus

Identification of amino acid residues involved in the interaction between measles virus haemagglutin (MVH) and its human cell receptor

Cotton rat (Sigmodon hispidus) signaling lymphocyte activation molecule (CD150) is an entry receptor for measles virus

The morbillivirus receptor SLAM (CD150)

MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform

Comparative protein modelling by satisfaction of spatial restraints

The ClusPro web server for protein-protein docking

Four distances between pairs of amino acids provide a precise description of their interaction

Advances in human-protein interactioninteractive and immersive molecular simulations

Analysis of global ecosystem ecology by fragment molecular orbital (FMO) method -analyses of the interactions between virus hemagglutinins and their receptors

Analysis of global ecosystem ecology by fragment molecular orbital (FMO) method -analyses of the interactions between virus hemagglutinins and their receptors

Measles vaccination of

The author thanks two anonymous reviewers for valuable comments. This study was supported by a Grant-in-Aid for Research in Nagoya City University to Y.S.

The author declares no conflict of interest.