The pulmonary surfactant protein A (SP-A) is a constitutively expressed immune-protective collagenous lectin (collectin) in the lung. It binds to the cell membrane of immune cells and opsonizes infectious agents such as bacteria, fungi, and viruses through glycoprotein binding. SARS-CoV-2 enters airway epithelial cells by ligating the Angiotensin Converting Enzyme 2 (ACE2) receptor on the cell surface using its Spike glycoprotein (S protein). We hypothesized that SP-A binds to the SARS-CoV-2 S protein and this binding interferes with ACE2 ligation. To study this hypothesis, we used a hybrid quantum and classical in silico modeling technique that utilized protein graph pruning. This graph pruning technique determines the best binding sites between amino acid chains by utilizing the Quantum Approximate Optimization Algorithm (QAOA)-based MaxCut (QAOA-MaxCut) program on a Near Intermediate Scale Quantum (NISQ) device. In this, the angles between every neighboring three atoms were Fourier-transformed into microwave frequencies and sent to a quantum chip that identified the chemically irrelevant atoms to eliminate based on their chemical topology. We confirmed that the remaining residues contained all the potential binding sites in the molecules by the Universal Protein Resource (UniProt) database. QAOA-MaxCut was compared with GROMACS with T-REMD using AMBER, OPLS, and CHARMM force fields to determine the differences in preparing a protein structure docking, as well as with Goemans-Williamson, the best classical algorithm for MaxCut. The relative binding affinity of potential interactions between the pruned protein chain residues of SP-A and SARS-CoV-2 S proteins was assessed by the ZDOCK program. Our data indicate that SP-A could ligate the S protein with a similar affinity to the ACE2-Spike binding. Interestingly, however, the results suggest that the most tightly-bound SP-A binding site is localized to the S2 chain, in the fusion region of the SARS-CoV-2 S protein, that is responsible for cell entry Based on these findings we speculate that SP-A may not directly compete with ACE2 for the binding site on the S protein, but interferes with viral entry to the cell by hindering necessary conformational changes or the fusion process.
Determining an optimal protein configuration for the employment of protein binding analysis as completed by Temperature based Replica Exchange Molecular Dynamics (T-REMD) is an important process in the accurate depiction of a protein’s behavior in different solvent environments, especially when determining a protein’s top binding sites for use in protein-ligand and protein-protein docking studies. However, the completion of this analysis, which pushes out top binding sites through configurational changes, is an polynomial-state computational problem that can take multiple hours to compute, even on the fastest supercomputers. In this study, we aim to determine if graph cutting provide approximated solutions the MaxCut problem can be used as a method to provide similar results to T-REMD in the determination of top binding sites of Surfactant Protein A (SP-A) for binding analysis. Additionally, we utilize a quantum-hybrid algorithm within Iff Technology’s Polar+ package using an actual quantum processor unit (QPU), an implementation of Polar+ using an emulated QPU, or Quantum Abstract Machine (QAM), on a large scale classical computing device, and an implementation of a classical MaxCut algorithm on a supercomputer in order to determine the types of advantages that can be gained through using a quantum computing device for this problem, or even using quantum algorithms on a classical device. This study found that Polar+ provides a dramatic speedup over a classical implementation of a MaxCut approximation algorithm or the use of GROMACS T-REMD, and produces viable results, in both its QPU and QAM implementations. However, the lack of direct configurational changes carried out onto the structure of SP-A after the use of graph cutting methods produces different final binding results than those produced by GROMACS T-REMD. Thus, further work needs to be completed into translating quantum-based probabilities into configurational changes based on a variety of noise conditions to better determine the accuracy advantage that quantum algorithms and quantum devices can provide in the near future.
Determining an optimal protein configuration for the employment of protein binding analysis as completed by Temperature based Replica Exchange Molecular Dynamics (T-REMD) is an important process in the accurate depiction of a protein’s behavior in different solvent environments, especially when determining a protein’s top binding sites for use in protein-ligand and protein-protein docking studies. However, the completion of this analysis, which pushes out top binding sites through configurational changes, is an polynomial-state computational problem that can take multiple hours to compute, even on the fastest supercomputers. In this study, we aim to determine if graph cutting provide approximated solutions the MaxCut problem can be used as a method to provide similar results to T-REMD in the determination of top binding sites of Surfactant Protein A (SP-A) for binding analysis. Additionally, we utilize a quantum-hybrid algorithm within Iff Technology’s Polar+ package using an actual quantum processor unit (QPU), an implementation of Polar+ using an emulated QPU, or Quantum Abstract Machine (QAM), on a large scale classical computing device, and an implementation of a classical MaxCut algorithm on a supercomputer in order to determine the types of advantages that can be gained through using a quantum computing device for this problem, or even using quantum algorithms on a classical device. This study found that Polar+ provides a dramatic speedup over a classical implementation of a MaxCut approximation algorithm or the use of GROMACS T-REMD, and produces viable results, in both its QPU and QAM implementations. However, the lack of direct configurational changes carried out onto the structure of SP-A after the use of graph cutting methods produces different final binding results than those produced by GROMACS T-REMD. Thus, further work needs to be completed into translating quantum-based probabilities into configurational changes based on a variety of noise conditions to better determine the accuracy advantage that quantum algorithms and quantum devices can provide in the near future.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.