1
|
Ose NJ, Campitelli P, Modi T, Kazan IC, Kumar S, Ozkan SB. Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics. eLife 2024; 12:RP92063. [PMID: 38713502 PMCID: PMC11076047 DOI: 10.7554/elife.92063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2024] Open
Abstract
We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.
Collapse
Affiliation(s)
- Nicholas James Ose
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Paul Campitelli
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Tushar Modi
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - I Can Kazan
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| | - Sudhir Kumar
- Institute for Genomics and Evolutionary Medicine, Temple UniversityPhiladelphiaUnited States
- Department of Biology, Temple UniversityPhiladelphiaUnited States
- Center for Genomic Medicine Research, King Abdulaziz UniversityJeddahSaudi Arabia
| | - Sefika Banu Ozkan
- Department of Physics and Center for Biological Physics, Arizona State UniversityTempeUnited States
| |
Collapse
|
2
|
Chen M, Jiang X, Zhang L, Chen X, Wen Y, Gu Z, Li X, Zheng M. The emergence of machine learning force fields in drug design. Med Res Rev 2024; 44:1147-1182. [PMID: 38173298 DOI: 10.1002/med.22008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2023] [Revised: 11/29/2023] [Accepted: 12/05/2023] [Indexed: 01/05/2024]
Abstract
In the field of molecular simulation for drug design, traditional molecular mechanic force fields and quantum chemical theories have been instrumental but limited in terms of scalability and computational efficiency. To overcome these limitations, machine learning force fields (MLFFs) have emerged as a powerful tool capable of balancing accuracy with efficiency. MLFFs rely on the relationship between molecular structures and potential energy, bypassing the need for a preconceived notion of interaction representations. Their accuracy depends on the machine learning models used, and the quality and volume of training data sets. With recent advances in equivariant neural networks and high-quality datasets, MLFFs have significantly improved their performance. This review explores MLFFs, emphasizing their potential in drug design. It elucidates MLFF principles, provides development and validation guidelines, and highlights successful MLFF implementations. It also addresses potential challenges in developing and applying MLFFs. The review concludes by illuminating the path ahead for MLFFs, outlining the challenges to be overcome and the opportunities to be harnessed. This inspires researchers to embrace MLFFs in their investigations as a new tool to perform molecular simulations in drug design.
Collapse
Affiliation(s)
- Mingan Chen
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Physical Science and Technology, ShanghaiTech University, Shanghai, China
- Lingang Laboratory, Shanghai, China
| | - Xinyu Jiang
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
| | - Lehan Zhang
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
| | - Xiaoxu Chen
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
- School of Pharmaceutical Science and Technology, Hangzhou Institute for Advanced Study, UCAS, Hangzhou, China
| | - Yiming Wen
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
- School of Pharmaceutical Science and Technology, Hangzhou Institute for Advanced Study, UCAS, Hangzhou, China
| | - Zhiyong Gu
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
- School of Pharmaceutical Science and Technology, Hangzhou Institute for Advanced Study, UCAS, Hangzhou, China
| | - Xutong Li
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
| | - Mingyue Zheng
- Drug Discovery and Design Center, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
- School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
- School of Pharmaceutical Science and Technology, Hangzhou Institute for Advanced Study, UCAS, Hangzhou, China
| |
Collapse
|
3
|
Kožić M, Bertoša B. Trajectory maps: molecular dynamics visualization and analysis. NAR Genom Bioinform 2024; 6:lqad114. [PMID: 38226394 PMCID: PMC10789246 DOI: 10.1093/nargab/lqad114] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Revised: 12/01/2023] [Accepted: 01/08/2024] [Indexed: 01/17/2024] Open
Abstract
Molecular dynamics simulations generate trajectories that depict system's evolution in time and are analyzed visually and quantitatively. Commonly conducted analyses include RMSD, Rgyr, RMSF, and more. However, those methods are all limited by their strictly statistical nature. Here we present trajectory maps, a novel method to analyze and visualize protein simulation courses intuitively and conclusively. By plotting protein's backbone movements during the simulation as a heatmap, trajectory maps provide new tools to directly visualize protein behavior over time, compare multiple simulations, and complement established methods. A user-friendly Python application developed for this purpose is presented, alongside detailed documentation for easy usage and implementation. The method's validation is demonstrated on three case studies. Considering its benefits, trajectory maps are expected to adopt broad application in obtaining and communicating meaningful results of protein molecular dynamics simulations in many associated fields such as biochemistry, structural biology, pharmaceutical research etc.
Collapse
Affiliation(s)
- Matej Kožić
- Department of Chemistry, Faculty of Science, University of Zagreb, Horvatovac 102a, HR-10000 Zagreb, Croatia
| | - Branimir Bertoša
- Department of Chemistry, Faculty of Science, University of Zagreb, Horvatovac 102a, HR-10000 Zagreb, Croatia
| |
Collapse
|
4
|
Mushebenge AGA, Ugbaja SC, Mbatha NA, B. Khan R, Kumalo HM. Assessing the Potential Contribution of In Silico Studies in Discovering Drug Candidates That Interact with Various SARS-CoV-2 Receptors. Int J Mol Sci 2023; 24:15518. [PMID: 37958503 PMCID: PMC10647470 DOI: 10.3390/ijms242115518] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 10/18/2023] [Accepted: 10/23/2023] [Indexed: 11/15/2023] Open
Abstract
The COVID-19 pandemic has spurred intense research efforts to identify effective treatments for SARS-CoV-2. In silico studies have emerged as a powerful tool in the drug discovery process, particularly in the search for drug candidates that interact with various SARS-CoV-2 receptors. These studies involve the use of computer simulations and computational algorithms to predict the potential interaction of drug candidates with target receptors. The primary receptors targeted by drug candidates include the RNA polymerase, main protease, spike protein, ACE2 receptor, and transmembrane protease serine 2 (TMPRSS2). In silico studies have identified several promising drug candidates, including Remdesivir, Favipiravir, Ribavirin, Ivermectin, Lopinavir/Ritonavir, and Camostat Mesylate, among others. The use of in silico studies offers several advantages, including the ability to screen a large number of drug candidates in a relatively short amount of time, thereby reducing the time and cost involved in traditional drug discovery methods. Additionally, in silico studies allow for the prediction of the binding affinity of the drug candidates to target receptors, providing insight into their potential efficacy. This study is aimed at assessing the useful contributions of the application of computational instruments in the discovery of receptors targeted in SARS-CoV-2. It further highlights some identified advantages and limitations of these studies, thereby revealing some complementary experimental validation to ensure the efficacy and safety of identified drug candidates.
Collapse
Affiliation(s)
- Aganze Gloire-Aimé Mushebenge
- Discipline of Pharmaceutical Sciences, University of KwaZulu-Natal, Westville, Durban 4000, South Africa;
- Drug Research and Innovation Unit, Discipline of Medical Biochemistry, School of Laboratory Medicine and Medical Science, University of KwaZulu-Natal, Durban 4000, South Africa
- Faculty of Pharmaceutical Sciences, University of Lubumbashi, Lubumbashi 1825, Democratic Republic of the Congo
| | - Samuel Chima Ugbaja
- Drug Research and Innovation Unit, Discipline of Medical Biochemistry, School of Laboratory Medicine and Medical Science, University of KwaZulu-Natal, Durban 4000, South Africa
- Africa Health Research Institute, University of KwaZulu-Natal, Durban 4000, South Africa
| | - Nonkululeko Avril Mbatha
- KwaZulu-Natal Research Innovation and Sequencing Platform, School of Laboratory Medicine and Medical Science, University of KwaZulu-Natal, Durban 4000, South Africa
| | - Rene B. Khan
- Drug Research and Innovation Unit, Discipline of Medical Biochemistry, School of Laboratory Medicine and Medical Science, University of KwaZulu-Natal, Durban 4000, South Africa
| | - Hezekiel M. Kumalo
- Drug Research and Innovation Unit, Discipline of Medical Biochemistry, School of Laboratory Medicine and Medical Science, University of KwaZulu-Natal, Durban 4000, South Africa
| |
Collapse
|
5
|
Pipitò L, Reynolds CA, Deganutti G. Is the Stalk of the SARS-CoV-2 Spike Protein Druggable? Viruses 2022; 14:2789. [PMID: 36560795 PMCID: PMC9786045 DOI: 10.3390/v14122789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Revised: 11/28/2022] [Accepted: 12/12/2022] [Indexed: 12/23/2022] Open
Abstract
The spike protein is key to SARS-CoV-2 high infectivity because it facilitates the receptor binding domain (RBD) encounter with ACE2. As targeting subunit S1 has not yet delivered an ACE2-binding inhibitor, we have assessed the druggability of the conserved segment of the spike protein stalk within subunit S2 by means of an integrated computational approach that combines the molecular docking of an optimized library of fragments with high-throughput molecular dynamics simulations. The high propensity of the spike protein to mutate in key regions that are responsible for the recognition of the human angiotensin-converting enzyme 2 (hACE2) or for the recognition of antibodies, has made subunit S1 of the spike protein difficult to target. Despite the inherent flexibility of the stalk region, our results suggest two hidden interhelical binding sites, whose accessibility is only partially hampered by glycan residues.
Collapse
Affiliation(s)
| | | | - Giuseppe Deganutti
- Centre for Sport, Exercise and Life Sciences, Faculty of Health and Life Sciences, Coventry University, Coventry CV1 5FB, UK
| |
Collapse
|
6
|
Chan KC, Song Y, Xu Z, Shang C, Zhou R. SARS-CoV-2 Delta Variant: Interplay between Individual Mutations and Their Allosteric Synergy. Biomolecules 2022; 12:biom12121742. [PMID: 36551170 PMCID: PMC9775976 DOI: 10.3390/biom12121742] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 11/17/2022] [Accepted: 11/20/2022] [Indexed: 11/25/2022] Open
Abstract
Since its first appearance in April 2021, B.1.617.2, also termed variant Delta, catalyzed one major worldwide wave dominating the second year of coronavirus disease 2019 (COVID-19) pandemic. Despite its quick disappearance worldwide, the strong virulence caused by a few point mutations remains an unsolved problem largely. Along with the other two sublineages, the Delta variant harbors an accumulation of Spike protein mutations, including the previously identified L452R, E484Q, and the newly emerged T478K on its receptor binding domain (RBD). We used molecular dynamics (MD) simulations, in combination with free energy perturbation (FEP) calculations, to examine the effects of two combinative mutation sets, L452R + E484Q and L452R + T478K. Our dynamic trajectories reveal an enhancement in binding affinity between mutated RBD and the common receptor protein angiotensin converting enzyme 2 (ACE2) through a net increase in the buried molecular surface area of the binary complex. This enhanced binding, mediated through Gln493, sets the same stage for all three sublineages due to the presence of L452R mutation. The other mutation component, E484Q or T478K, was found to impact the RBD-ACE2 binding and help the variant to evade several monoclonal antibodies (mAbs) in a distinct manner. Especially for L452R + T478K, synergies between mutations are mediated through a complex residual and water interaction network and further enhance its binding to ACE2. Taking together, this study demonstrates that new variants of SARS-CoV-2 accomplish both "attack" (infection) and "defense" (antibody neutralization escape) with the same "polished sword" (mutated Spike RBD).
Collapse
Affiliation(s)
- Kevin C. Chan
- Institute of Quantitative Biology, College of Life Sciences, Zhejiang University, Hangzhou 310058, China
- Shanghai Institute for Advanced Study, Zhejiang University, 799 Dangui Road, Shanghai 201203, China
| | - Yi Song
- Institute of Quantitative Biology, College of Life Sciences, Zhejiang University, Hangzhou 310058, China
| | - Zheng Xu
- BirenTech Research, Shanghai 201112, China
| | - Chun Shang
- Institute of Quantitative Biology, College of Life Sciences, Zhejiang University, Hangzhou 310058, China
| | - Ruhong Zhou
- Institute of Quantitative Biology, College of Life Sciences, Zhejiang University, Hangzhou 310058, China
- Shanghai Institute for Advanced Study, Zhejiang University, 799 Dangui Road, Shanghai 201203, China
- Department of Chemistry, Columbia University, New York, NY 10027, USA
- The First Affiliated Hospital, School of Medicine, Zhejiang University, Hangzhou 310058, China
- Correspondence:
| |
Collapse
|
7
|
Aksenova AY, Likhachev IV, Grishin SY, Galzitskaya OV. The Increased Amyloidogenicity of Spike RBD and pH-Dependent Binding to ACE2 May Contribute to the Transmissibility and Pathogenic Properties of SARS-CoV-2 Omicron as Suggested by In Silico Study. Int J Mol Sci 2022; 23:13502. [PMID: 36362302 PMCID: PMC9655063 DOI: 10.3390/ijms232113502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Revised: 10/19/2022] [Accepted: 11/01/2022] [Indexed: 11/06/2022] Open
Abstract
SARS-CoV-2 is a rapidly evolving pathogen that has caused a global pandemic characterized by several consecutive waves. Based on epidemiological and NGS data, many different variants of SARS-CoV-2 were described and characterized since the original variant emerged in Wuhan in 2019. Notably, SARS-CoV-2 variants differ in transmissibility and pathogenicity in the human population, although the molecular basis for this difference is still debatable. A significant role is attributed to amino acid changes in the binding surface of the Spike protein to the ACE2 receptor, which may facilitate virus entry into the cell or contribute to immune evasion. We modeled in silico the interaction between Spike RBDs of Wuhan-Hu-1, Delta, and Omicron BA.1 variants and ACE2 at different pHs (pH 5 and pH 7) and showed that the strength of this interaction was higher for the Omicron BA.1 RBD compared to Wuhan-Hu-1 or Delta RBDs and that the effect was more profound at pH 5. This finding is strikingly related to the increased ability of Omicron variants to spread in the population. We also noted that during its spread in the population, SARS-CoV-2 evolved to a more charged, basic composition. We hypothesize that the more basic surface of the Omicron variant may facilitate its spread in the upper respiratory tract but not in the lower respiratory tract, where pH estimates are different. We calculated the amyloidogenic properties of Spike RBDs in different SARS-CoV-2 variants and found eight amyloidogenic regions in the Spike RBDs for each of the variants predicted by the FoldAmyloid program. Although all eight regions were almost identical in the Wuhan to Gamma variants, two of them were significantly longer in both Omicron variants, making the Omicron RBD more amyloidogenic. We discuss how the increased predicted amyloidogenicity of the Omicron variants RBDs may be important for protein stability, influence its interaction with ACE2 and contribute to immune evasion.
Collapse
Affiliation(s)
- Anna Y. Aksenova
- Laboratory of Amyloid Biology, St. Petersburg State University, 199034 St. Petersburg, Russia
| | - Ilya V. Likhachev
- Institute of Protein Research, Russian Academy of Sciences, 142290 Pushchino, Russia
- Institute of Mathematical Problems of Biology RAS, The Branch of Keldysh Institute of Applied Mathematics, Russian Academy of Sciences, 142290 Pushchino, Russia
| | - Sergei Y. Grishin
- Institute of Protein Research, Russian Academy of Sciences, 142290 Pushchino, Russia
- Institute of Environmental and Agricultural Biology (X-BIO), Tyumen State University, 625003 Tyumen, Russia
| | - Oxana V. Galzitskaya
- Institute of Protein Research, Russian Academy of Sciences, 142290 Pushchino, Russia
- Institute of Theoretical and Experimental Biophysics, Russian Academy of Sciences, 142290 Pushchino, Russia
| |
Collapse
|
8
|
Pipitò L, Reynolds CA, Mobarec JC, Vickery O, Deganutti G. A Pathway Model to Understand the Evolution of Spike Protein Binding to ACE2 in SARS-CoV-2 Variants. Biomolecules 2022; 12:1607. [PMID: 36358957 PMCID: PMC9687612 DOI: 10.3390/biom12111607] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 10/25/2022] [Accepted: 10/28/2022] [Indexed: 01/16/2024] Open
Abstract
After the SARS-CoV-2 Wuhan variant that gave rise to the pandemic, other variants named Delta, Omicron, and Omicron-2 sequentially became prevalent, with mutations spread around the viral genome, including on the spike (S) protein; in order to understand the resultant in gains in infectivity, we interrogated in silico both the equilibrium binding and the binding pathway of the virus' receptor-binding domain (RBD) to the angiotensin-converting enzyme 2 (ACE2) receptor. We interrogated the molecular recognition between the RBD of different variants and ACE2 through supervised molecular dynamics (SuMD) and classic molecular dynamics (MD) simulations to address the effect of mutations on the possible S protein binding pathways. Our results indicate that compensation between binding pathway efficiency and stability of the complex exists for the Omicron BA.1 receptor binding domain, while Omicron BA.2's mutations putatively improved the dynamic recognition of the ACE2 receptor, suggesting an evolutionary advantage over the previous strains.
Collapse
Affiliation(s)
- Ludovico Pipitò
- Centre for Sport, Exercise and Life Sciences, Faculty of Health and Life Sciences, Coventry University, Coventry CV1 5FB, UK
| | - Christopher A. Reynolds
- Centre for Sport, Exercise and Life Sciences, Faculty of Health and Life Sciences, Coventry University, Coventry CV1 5FB, UK
- School of Life Sciences, University of Essex, Wivenhoe Park, Colchester CO4 3SQ, UK
| | - Juan Carlos Mobarec
- Mechanistic and Structural Biology, Discovery Sciences, R&D, AstraZeneca, Cambridge, UK
| | - Owen Vickery
- Mechanistic and Structural Biology, Discovery Sciences, R&D, AstraZeneca, Cambridge, UK
| | - Giuseppe Deganutti
- Centre for Sport, Exercise and Life Sciences, Faculty of Health and Life Sciences, Coventry University, Coventry CV1 5FB, UK
| |
Collapse
|