Computational Analysis of the Interaction Energies between Amino Acid Residues of the Measles Virus Hemagglutinin and Its Receptors

Protein-Protein Interaction Modelling with the Fragment Molecular Orbital Method

Fragment molecular orbital (FMO) method enables ab initio quantum-chemical calculations for biomolecular systems with high accuracy and moderate computational cost. Through this analysis we can evaluate the inter-fragment interaction energies (IFIEs) that provide useful measures for effective interactions between the fragments representing amino-acid residues and ligand molecules. Here I describe how to prepare the input structures and perform the FMO calculations for protein-protein complex system. In addition to the pre-processing, some useful tools for the post-processing analysis are also illustrated.

Novel and classical morbilliviruses: Current knowledge of three divergent morbillivirus groups

Currently, seven species of morbillivirus have been classified. Six of these species (Measles morbillivirus, Rinderpest morbillivirus, Small ruminant morbillivirus, Canine morbillivirus, Phocine morbillivirus, and Cetacean morbillivirus) are highly infectious and cause serious systemic diseases in humans, livestock, domestic dogs, and wild animals. These species commonly use the host proteins signaling lymphocytic activation molecule (SLAM) and nectin-4 as receptors, and this usage contributes to their virulence. The seventh species (Feline morbillivirus: FeMV) is phylogenetically divergent from the six SLAM-using species. FeMV differs from the SLAM-using morbillivirus group in pathogenicity and infectivity, and is speculated to use non-SLAM receptors. Recently, novel species of morbilliviruses have been discovered in bats, rodents, and domestic pigs. Because the ability to use SLAM and nectin-4 is closely related to the infectivity and pathogenicity of morbilliviruses, investigation of the potential usage of these receptors is useful for estimating infectivity and pathogenicity. The SLAM-binding sites in the receptor-binding protein show high similarity among the SLAM-using morbilliviruses. This feature may help to estimate whether novel morbillivirus species can use SLAM as a receptor. A novel morbillivirus species isolated from wild mice diverged from the classified morbilliviruses in the phylogenetic tree, forming a third group separate from the SLAM-using morbillivirus group and FeMV. This suggests that the novel rodent morbillivirus may exhibit a different risk from the SLAM-using morbillivirus group, and analyses of its viral pathogenicity and infectivity toward humans are warranted.

The Importance of Charge Transfer and Solvent Screening in the Interactions of Backbones and Functional Groups in Amino Acid Residues and Nucleotides

Quantum mechanical (QM) calculations at the level of density-functional tight-binding are applied to a protein–DNA complex (PDB: 2o8b) consisting of 3763 atoms, averaging 100 snapshots from molecular dynamics simulations. A detailed comparison of QM and force field (Amber) results is presented. It is shown that, when solvent screening is taken into account, the contributions of the backbones are small, and the binding of nucleotides in the double helix is governed by the base–base interactions. On the other hand, the backbones can make a substantial contribution to the binding of amino acid residues to nucleotides and other residues. The effect of charge transfer on the interactions is also analyzed, revealing that the actual charge of nucleotides and amino acid residues can differ by as much as 6 and 8% from the formal integer charge, respectively. The effect of interactions on topological models (protein -residue networks) is elucidated.

FeMV is a cathepsin-dependent unique morbillivirus infecting the kidneys of domestic cats

Feline morbillivirus (FeMV) is a recently discovered pathogen of domestic cats and has been classified as a morbillivirus in the Paramyxovirus family. We determined the complete sequence of FeMV^US5 directly from an FeMV-positive urine sample without virus isolation or cell passage. Sequence analysis of the viral genome revealed potential divergence from characteristics of archetypal morbilliviruses. First, the virus lacks the canonical polybasic furin cleavage signal in the fusion (F) glycoprotein. Second, conserved amino acids in the hemagglutinin (H) glycoprotein used by all other morbilliviruses for binding and/or fusion activation with the cellular receptor CD150 (signaling lymphocyte activation molecule [SLAM]/F1) are absent. We show that, despite this sequence divergence, FeMV H glycoprotein uses feline CD150 as a receptor and cannot use human CD150. We demonstrate that the protease responsible for cleaving the FeMV F glycoprotein is a cathepsin, making FeMV a unique morbillivirus and more similar to the closely related zoonotic Nipah and Hendra viruses. We developed a reverse genetics system for FeMV^US5 and generated recombinant viruses expressing Venus fluorescent protein from an additional transcription unit located either between the phospho-protein (P) and matrix (M) genes or the H and large (L) genes of the genome. We used these recombinant FeMVs to establish a natural infection and demonstrate that FeMV causes an acute morbillivirus-like disease in the cat. Virus was shed in the urine and detectable in the kidneys at later time points. This opens the door for long-term studies to address the postulated role of this morbillivirus in the development of chronic kidney disease.

Molecular signatures in cetacean morbillivirus and host species proteomes: Unveiling the evolutionary dynamics of an enigmatic pathogen?

Cetacean morbillivirus (CeMV) infects marine mammals often causing a fatal respiratory and neurological disease. Recently, CeMV has expanded its geographic and host species range, with cases being reported worldwide among dolphins, whales, seals, and other aquatic mammalian species, and therefore has emerged as the most threatening nonanthropogenic factor affecting marine mammal's health and conservation. Extensive research efforts have aimed to understand CeMV epidemiology and ecology, however, the molecular mechanisms underlying its transmission and pathogenesis are still poorly understood. In particular, the field suffers from a knowledge gap on the structural and functional properties of CeMV proteins and their host interactors. Nevertheless, the body of scientific literature produced in recent years has inaugurated new investigational trends, driving future directions in CeMV molecular research. In this mini-review, the most recent literature has been summarized in the context of such research trends, and categorized into four priority research topics, such as (1) the interaction between CeMV glycoprotein and its host cell receptors across several species; (2) the CeMV molecular determinants responsible for different disease phenotype; (3) the host molecular determinants responsible for differential susceptibility to CeMV infection; (4) the CeMV molecular determinants responsible for difference virulence among circulating CeMV strains. Arguably, these are the most urgent topics that need to be investigated and that most promisingly will help to shed light on the details of CeMV evolutionary dynamics in the immediate future.

Intermolecular Interaction Analyses on SARS-CoV-2 Spike Protein Receptor Binding Domain and Human Angiotensin-Converting Enzyme 2 Receptor-Blocking Antibody/Peptide Using Fragment Molecular Orbital Calculation

The spike glycoprotein (S-protein) mediates SARS-CoV-2 entry via intermolecular interaction with human angiotensin-converting enzyme 2. The receptor binding domain (RBD) of the S-protein has been considered critical for this interaction and acts as the target of numerous neutralizing antibodies and antiviral peptides. This study used the fragment molecular orbital method to analyze the interactions between the RBD and antibodies/peptides and extracted crucial residues that can be used as epitopes. The interactions evaluated as interfragment interaction energy values between the RBD and 12 antibodies/peptides showed a fairly good correlation with the experimental activity pIC50 (R2 = 0.540). Nine residues (T415, K417, Y421, F456, A475, F486, N487, N501, and Y505) were confirmed as being crucial. Pair interaction energy decomposition analyses showed that hydrogen bonds, electrostatic interactions, and π-orbital interactions are important. Our results provide essential information for understanding SARS-CoV-2-antibody/peptide binding and may play roles in future antibody/antiviral drug design.

Porcine Complement Regulatory Protein CD46 Is a Major Receptor for Atypical Porcine Pestivirus but Not for Classical Swine Fever Virus

Pestiviruses comprise animal pathogens such as classical swine fever virus (CSFV) and bovine viral diarrhea virus (BVDV) that cause notifiable diseases with great economic impact. Several additional pestivirus species affecting animal health were recently identified, including atypical porcine pestivirus (APPV).

Pestiviruses such as bovine viral diarrhea virus (BVDV) and classical swine fever virus (CSFV) belong to the family Flaviviridae and represent pathogens of outstanding veterinary relevance. Pestiviruses enter cells via receptor-mediated endocytosis. For entry in bovine cells, complement regulatory protein CD46_bov serves as a cellular receptor for BVDV. In this study, the role of porcine CD46_pig in cellular entry was investigated for the recently discovered atypical porcine pestivirus (APPV), CSFV, and Bungowannah virus (BuPV) in order to elucidate the observed differences in host cell tropism. A cell culture-adapted APPV variant, which shows enhanced viral replication in vitro, was generated and demonstrated a strict tropism of APPV for porcine cells. One of the porcine cell lines displayed areas of CD46_pig-expressing cells and areas of nonexpressing cells, and one single cell line revealed not to express any CD46_pig. The CD46_pig-deficient porcine lymphoma cell line, known to facilitate CSFV replication, was the only porcine cell line nonpermissive to APPV, indicating a significant difference in the entry mechanism of APPV and CSFV. Infection experiments with a set of genetically engineered CD46_pig knockout cells confirmed that CD46_pig is a major receptor of APPV as CD46_bov is for BVDV. In contrast, it is apparently not an essential determinant in host cell entry of other porcine pestiviruses such as CSFV and BuPV. Existence of a CD46_pig-independent entry mechanism illustrates that the pestiviral entry process is more diverse than previously recognized.

IMPORTANCE Pestiviruses comprise animal pathogens such as classical swine fever virus (CSFV) and bovine viral diarrhea virus (BVDV) that cause notifiable diseases with great economic impact. Several additional pestivirus species affecting animal health were recently identified, including atypical porcine pestivirus (APPV). APPV is associated with health problems in piglets and is highly abundant in pig populations worldwide. Complement control protein CD46 serves as a receptor for diverse bacterial and viral pathogens, including particular adenoviruses, herpesviruses, measles virus (MeV), and BVDV. Porcine CD46 (CD46_pig) was suggested to be a major receptor for CSFV. Here, we identified remarkable differences in relevance of CD46_pig during entry of porcine pestiviruses. Resembling BVDV, efficient APPV infection in cell culture depends on CD46_pig, while other porcine pestiviruses can efficiently enter and infect cells in the absence of CD46_pig. Thus, the study provides insights into the entry process of these pathogens and may help to understand differences in their biology.

Molecular recognition of SARS-CoV-2 spike glycoprotein: quantum chemical hot spot and epitope analyses

Due to the COVID-19 pandemic, researchers have attempted to identify complex structures of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S-protein) with angiotensin-converting enzyme 2 (ACE2) or a blocking antibody. However, the molecular recognition mechanism-critical information for drug and antibody design-has not been fully clarified at the amino acid residue level. Elucidating such a microscopic mechanism in detail requires a more accurate molecular interpretation that includes quantum mechanics to quantitatively evaluate hydrogen bonds, XH/π interactions (X = N, O, and C), and salt bridges. In this study, we applied the fragment molecular orbital (FMO) method to characterize the SARS-CoV-2 S-protein binding interactions with not only ACE2 but also the B38 Fab antibody involved in ACE2-inhibitory binding. By analyzing FMO-based interaction energies along a wide range of binding interfaces carefully, we identified amino acid residues critical for molecular recognition between S-protein and ACE2 or B38 Fab antibody. Importantly, hydrophobic residues that are involved in weak interactions such as CH-O hydrogen bond and XH/π interactions, as well as polar residues that construct conspicuous hydrogen bonds, play important roles in molecular recognition and binding ability. Moreover, through these FMO-based analyses, we also clarified novel hot spots and epitopes that had been overlooked in previous studies by structural and molecular mechanical approaches. Altogether, these hot spots/epitopes identified between S-protein and ACE2/B38 Fab antibody may provide useful information for future antibody design, evaluation of the binding property of the SARS-CoV-2 variants including its N501Y, and small or medium drug design against the SARS-CoV-2.

Interaction analyses of SARS-CoV-2 spike protein based on fragment molecular orbital calculations

At the stage of SARS-CoV-2 infection in human cells, the spike protein consisting of three chains, A, B, and C, with a total of 3300 residues plays a key role, and thus its structural properties and the binding nature of receptor proteins to host human cells or neutralizing antibodies has attracted considerable interest. Here, we report on interaction analyses of the spike protein in both closed (PDB-ID: 6VXX) and open (6VYB) structures, based on large-scale fragment molecular orbital (FMO) calculations at the level of up to the fourth-order Møller–Plesset perturbation with singles, doubles, and quadruples (MP4(SDQ)). Inter-chain interaction energies were evaluated for both structures, and a mutual comparison indicated considerable losses of stabilization energies in the open structure, especially in the receptor binding domain (RBD) of chain-B. The role of charged residues in inter-chain interactions was illuminated as well. By two separate calculations for the RBD complexes with angiotensin-converting enzyme 2 (ACE2) (6M0J) and B38 Fab antibody (7BZ5), it was found that the binding with ACE2 or antibody partially compensated for this stabilization loss of RBD.

Visualized IFIE results seen from chain-B of spike protein.

Identification of correlated inter-residue interactions in protein complex based on the fragment molecular orbital method

A theoretical scheme to systematically describe correlated (network-like) interactions between molecular fragments is proposed within the framework of the fragment molecular orbital (FMO) method. The method is mathematically based on the singular value decomposition (SVD) of the inter-fragment interaction energy (IFIE) matrix obtained by the FMO calculation, and can be applied to a comprehensive description of protein-protein interactions in the context of molecular recognition. In the present study we apply the proposed method to a complex of measles virus hemagglutinin and human SLAM receptor, thus finding a usefulness for efficiently eliciting the correlated interactions among the amino-acid residues involved in the two proteins. Additionally, collective interaction networks by amino-acid residues important for mutation experiments can be clearly visualized.

Fragment Molecular Orbital Based Interaction Analyses on COVID-19 Main Protease - Inhibitor N3 Complex (PDB ID: 6LU7)

The worldwide spread of COVID-19 (new coronavirus found in 2019) is an emergent issue to be tackled. In fact, a great amount of works in various fields have been made in a rather short period. Here, we report a fragment molecular orbital (FMO) based interaction analysis on a complex between the SARS-CoV-2 main protease (Mpro) and its peptide-like inhibitor N3 (PDB ID: 6LU7). The target inhibitor molecule was segmented into five fragments in order to capture site specific interactions with amino acid residues of the protease. The interaction energies were decomposed into several contributions, and then the characteristics of hydrogen bonding and dispersion stabilization were made clear. Furthermore, the hydration effect was incorporated by the Poisson-Boltzmann (PB) scheme. From the present FMO study, His41, His163, His164, and Glu166 were found to be the most important amino acid residues of Mpro in interacting with the inhibitor, mainly due to hydrogen bonding. A guideline for optimizations of the inhibitor molecule was suggested as well based on the FMO analysis.

Identification of amino acid residues involved in the interaction between peste-des-petits-ruminants virus haemagglutinin protein and cellular receptors

Peste-des-petits-ruminants virus (PPRV) haemagglutinin (H) protein mediates binding to cellular receptors and then initiates virus entry. To identify the key residues of PPRV H (Hv) protein of the Nigeria 75/1 strain involved in binding to receptors, interaction of the Hv and mutated Hv (mHv) proteins with receptors (SLAM and Nectin 4) and their mutants (mSLAM1, mSLAM2, mSLAM3 and mNectin 4) was investigated using surface plasmon resonance imaging (SPRi) and coimmunoprecipitation (co-IP) assays. The results showed that the Hv protein failed to interact with mSLAM3, but interacted at a strong or medium intensity with SLAM, mSLAM2, Nectin 4 and mNectin 4, and at a low level with mSLAM1. The mHv protein was unable to interact with SLAM and its mutants, but bound to Nectin 4 and mNectin 4 with medium and weak intensity, respectively. Further analysis showed that the Hv protein could precipitate mSLAM1, mSLAM2 and mNectin 4, but not mSLAM3. The mHv protein failed to coprecipitate with SLAM and its mutants. The binding activities of mNectin 4 and Nectin 4 to mHv were less than 30.36 and 51.94 % of the wild-type levels, respectively. Based on the results obtained, amino acids at positions R389, L464, I498, R503, R533, Y541, Y543, F552 and Y553 of H protein and I61, H62, L64, K76, K78, E123, H130, I210, A211, S226 and R227 in SLAM were identified to be essential for the specificity of H-SLAM interaction, while the critical residues of H-Nectin 4 interaction require further study. These findings would improve our understanding of the invasive mechanisms of PPRV.

Canine and Phocine Distemper Viruses: Global Spread and Genetic Basis of Jumping Species Barriers

Canine distemper virus (CDV) and phocine distemper (PDV) are closely-related members of the Paramyxoviridae family, genus morbillivirus, in the order Mononegavirales. CDV has a broad host range among carnivores. PDV is thought to be derived from CDV through contact between terrestrial carnivores and seals. PDV has caused extensive mortality in Atlantic seals and other marine mammals, and more recently has spread to the North Pacific Ocean. CDV also infects marine carnivores, and there is evidence of morbillivirus infection of seals and other species in Antarctica. Recently, CDV has spread to felines and other wildlife species in the Serengeti and South Africa. Some CDV vaccines may also have caused wildlife disease. Changes in the virus haemagglutinin (H) protein, particularly the signaling lymphocyte activation molecule (SLAM) receptor binding site, correlate with adaptation to non-canine hosts. Differences in the phosphoprotein (P) gene sequences between disease and non-disease causing CDV strains may relate to pathogenicity in domestic dogs and wildlife. Of most concern are reports of CDV infection and disease in non-human primates raising the possibility of zoonosis. In this article we review the global occurrence of CDV and PDV, and present both historical and genetic information relating to these viruses crossing species barriers.

Solvent Screening in Zwitterions Analyzed with the Fragment Molecular Orbital Method

Based on induced solvent charges, a new model of solvent screening is developed in the framework of the fragment molecular orbital combined with the polarizable continuum model. The developed model is applied to analyze interactions in a prototypical zwitterionic system, sodium chloride in water, and it is shown that the large underestimation of the interaction in the original solvent screening based on local charges is successfully corrected. The model is also applied to a complex of the Trp-cage (PDB: 1L2Y ) miniprotein with an anionic ligand, and the physical factors determined protein-ligand binding in solution are unraveled.

Host Cellular Receptors for the Peste des Petits Ruminant Virus

Peste des Petits Ruminant (PPR) is an important transboundary, OIE-listed contagious viral disease of primarily sheep and goats caused by the PPR virus (PPRV), which belongs to the genus Morbillivirus of the family Paramyxoviridae. The mortality rate is 90–100%, and the morbidity rate may reach up to 100%. PPR is considered economically important as it decreases the production and productivity of livestock. In many endemic poor countries, it has remained an obstacle to the development of sustainable agriculture. Hence, proper control measures have become a necessity to prevent its rapid spread across the world. For this, detailed information on the pathogenesis of the virus and the virus host interaction through cellular receptors needs to be understood clearly. Presently, two cellular receptors; signaling lymphocyte activation molecule (SLAM) and Nectin-4 are known for PPRV. However, extensive information on virus interactions with these receptors and their impact on host immune response is still required. Hence, a thorough understanding of PPRV receptors and the mechanism involved in the induction of immunosuppression is crucial for controlling PPR. In this review, we discuss PPRV cellular receptors, viral host interaction with cellular receptors, and immunosuppression induced by the virus with reference to other Morbilliviruses.

Marine Morbilliviruses: Diversity and Interaction with Signaling Lymphocyte Activation Molecules

Epidemiological reports of phocine distemper virus (PDV) and cetacean morbillivirus (CeMV) have accumulated since their discovery nearly 30 years ago. In this review, we focus on the interaction between these marine morbilliviruses and their major cellular receptor, the signaling lymphocyte activation molecule (SLAM). The three-dimensional crystal structure and homology models of SLAMs have demonstrated that 35 residues are important for binding to the morbillivirus hemagglutinin (H) protein and contribute to viral tropism. These 35 residues are essentially conserved among pinnipeds and highly conserved among the Caniformia, suggesting that PDV can infect these animals, but are less conserved among cetaceans. Because CeMV can infect various cetacean species, including toothed and baleen whales, the CeMV-H protein is postulated to have broader specificity to accommodate more divergent SLAM interfaces and may enable the virus to infect seals. In silico analysis of viral H protein and SLAM indicates that each residue of the H protein interacts with multiple residues of SLAM and vice versa. The integration of epidemiological, virological, structural, and computational studies should provide deeper insight into host specificity and switching of marine morbilliviruses.