Influenza viral hemagglutinin complicated shape is advantageous to its binding affinity for sialosaccharide receptor

System truncation accelerates binding affinity calculations with the fragment molecular orbital method: A benchmark study

The fragment molecular orbital (FMO) method is a fast quantum-mechanical method that divides systems into pieces of fragments and performs ab initio calculations. The system truncation enables further speed improvement. In this article, we systematically study the effects of system truncations on binding affinity calculations obtained with FMO in combination with either the polarizable continuum model (FMO/PCM) or in combination with the Møller-Plesset method (FMO-MP2). We have used five protein complexes with ligands of several charged states. The calculated binding energies of the size variants of the truncated system, including only a restricted number of atoms around the ligand, are compared to the energy obtained from a full system. The result shows that the systems could be truncated to a radius of 8 Å from neutral ligands within an error of 0.7 kcal/mol, and 12 Å from charged ligands within an error of 1.1 kcal/mol for calculating the binding energy in solution.

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.

Systematic Interaction Analysis of Anti-Human Immunodeficiency Virus Type-1 Neutralizing Antibodies with High Mannose Glycans Using Fragment Molecular Orbital and Molecular Dynamics Methods

A series of broadly neutralizing antibodies called PGT have been shown to be bound directly to human immunodeficiency virus type-1 via high mannose glycans on glycoprotein gp120. Despite the sequence similarities of amino acids of the antibodies, their affinities to the glycan differ. Glycan-antibody interactions among these antibodies are systematically compared with quantum chemical fragment molecular orbital calculations and molecular dynamics simulations. The differences among structural stability of the glycan in the active site of the complexes and total interaction energies as well as binding free energies between the glycan and antibodies agree well with the experimentally shown affinities of the glycan to the antibodies. The terminal saccharide, Man D3, is structurally stable and responsible for the glycan-antibody binding through electrostatic and dispersion interactions. The structural stability of nonterminal saccharides such as Man 4 or Man C plays substantial roles in the interaction via direct hydrogen bonds. © 2019 Wiley Periodicals, Inc.

Energy decomposition analysis approaches and their evaluation on prototypical protein–drug interaction patterns

The partitioning of the interaction energy into chemical components such as electrostatics, polarization, and charge transfer is possible with energy decomposition analysis approaches. We review and evaluate these for biomolecular applications.

Comparison of the local structural stabilities of mammalian prion protein (PrP) by fragment molecular orbital calculations

Bovine spongiform encephalopathy (BSE), a member of the prion diseases, is a fatal neurodegenerative disorder suspected to be caused by a malfunction of prion protein (PrP). Although BSE prions have been reported to be transmitted to a wide range of animal species, dogs and hamsters are known to be BSE-resistant animals. Analysis of canine and hamster PrP could elucidate the molecular mechanisms supporting the species barriers to BSE prion transmission. The structural stability of 6 mammalian PrPs, including human, cattle, mouse, hamster, dog and cat, was analyzed. We then evaluated intramolecular interactions in PrP by fragment molecular orbital (FMO) calculations. Despite similar backbone structures, the PrP side-chain orientations differed among the animal species examined. The pair interaction energies between secondary structural elements in the PrPs varied considerably, indicating that the local structural stabilities of PrP varied among the different animal species. Principal component analysis (PCA) demonstrated that different local structural stability exists in bovine PrP compared with the PrP of other animal species examined. The results of the present study suggest that differences in local structural stabilities between canine and bovine PrP link diversity in susceptibility to BSE prion infection.

An Efficient Method to Evaluate Intermolecular Interaction Energies in Large Systems Using Overlapping Multicenter ONIOM and the Fragment Molecular Orbital Method

We propose an approach based on the overlapping multicenter ONIOM to evaluate intermolecular interaction energies in large systems and demonstrate its accuracy on several representative systems in the complete basis set limit at the MP2 and CCSD(T) level of theory. In the application to the intermolecular interaction energy between insulin dimer and 4'-hydroxyacetanilide at the MP2/CBS level, we use the fragment molecular orbital method for the calculation of the entire complex assigned to the lowest layer in three-layer ONIOM. The developed method is shown to be efficient and accurate in the evaluation of the protein-ligand interaction energies.

Rapid estimation of binding activity of influenza virus hemagglutinin to human and avian receptors

A critical step for avian influenza viruses to infect human hosts and cause epidemics or pandemics is acquisition of the ability of the viral hemagglutinin (HA) to bind to human receptors. However, current global influenza surveillance does not monitor HA binding specificity due to a lack of rapid and reliable assays. Here we report a computational method that uses an effective scoring function to quantify HA-receptor binding activities with high accuracy and speed. Application of this method reveals receptor specificity changes and its temporal relationship with antigenicity changes during the evolution of human H3N2 viruses. The method predicts that two amino acid differences at 222 and 225 between HAs of A/Fujian/411/02 and A/Panama/2007/99 viruses account for their differences in binding to both avian and human receptors; this prediction was verified experimentally. The new computational method could provide an urgently needed tool for rapid and large-scale analysis of HA receptor specificities for global influenza surveillance.

Influenza vaccines: the good, the bad, and the eggs

Outbreaks of influenza A viruses continue to cause morbidity and mortality worldwide. The global disease burden of influenza is substantial. While antiviral therapies are available, influenza vaccines are the mainstay of efforts to reduce the substantial health burden from seasonal influenza. Inactivated influenza vaccines have been available since the 1940s, with live attenuated, cold-adapted vaccines becoming available in the United States in 2003. In spite of the successes, more research is needed to develop more effective seasonal influenza vaccines that provide long-lasting immunity and broad protection against strains that differ antigenically from vaccine viruses. This review introduces the virus and its disease, the current state of seasonal and pandemic influenza vaccines, and the challenges we face in the future.

Extending the Power of Quantum Chemistry to Large Systems with the Fragment Molecular Orbital Method

Following the brief review of the modern fragment-based methods and other approaches to perform quantum-mechanical calculations of large systems, the theoretical development of the fragment molecular orbital method (FMO) is covered in detail, with the emphasis on the physical properties, which can be computed with FMO. The FMO-based polarizable continuum model (PCM) for treating the solvent effects in large systems and the pair interaction energy decomposition analysis (PIEDA) are described in some detail, and a range of applications of FMO to biological studies is introduced. The factors determining the relative stability of polypeptide conformers (alpha-helix, beta-turn, and extended form) are elucidated using FMO/PCM and PIEDA, and the interactions in the Trp-cage miniprotein construct (PDB: 1L2Y) are analyzed using PIEDA.