1
|
Casalino L, Dommer AC, Gaieb Z, Barros EP, Sztain T, Ahn SH, Trifan A, Brace A, Bogetti AT, Clyde A, Ma H, Lee H, Turilli M, Khalid S, Chong LT, Simmerling C, Hardy DJ, Maia JD, Phillips JC, Kurth T, Stern AC, Huang L, McCalpin JD, Tatineni M, Gibbs T, Stone JE, Jha S, Ramanathan A, Amaro RE. AI-driven multiscale simulations illuminate mechanisms of SARS-CoV-2 spike dynamics. Int J High Perform Comput Appl 2021; 35:432-451. [PMID: 38603008 PMCID: PMC8064023 DOI: 10.1177/10943420211006452] [Citation(s) in RCA: 58] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
We develop a generalizable AI-driven workflow that leverages heterogeneous HPC resources to explore the time-dependent dynamics of molecular systems. We use this workflow to investigate the mechanisms of infectivity of the SARS-CoV-2 spike protein, the main viral infection machinery. Our workflow enables more efficient investigation of spike dynamics in a variety of complex environments, including within a complete SARS-CoV-2 viral envelope simulation, which contains 305 million atoms and shows strong scaling on ORNL Summit using NAMD. We present several novel scientific discoveries, including the elucidation of the spike's full glycan shield, the role of spike glycans in modulating the infectivity of the virus, and the characterization of the flexible interactions between the spike and the human ACE2 receptor. We also demonstrate how AI can accelerate conformational sampling across different systems and pave the way for the future application of such methods to additional studies in SARS-CoV-2 and other molecular systems.
Collapse
Affiliation(s)
- Lorenzo Casalino
- University of California San Diego, La Jolla, CA, USA
- Authors with symbol indicate equal contribution
| | - Abigail C Dommer
- University of California San Diego, La Jolla, CA, USA
- Authors with symbol indicate equal contribution
| | - Zied Gaieb
- University of California San Diego, La Jolla, CA, USA
- Authors with symbol indicate equal contribution
| | | | - Terra Sztain
- University of California San Diego, La Jolla, CA, USA
| | - Surl-Hee Ahn
- University of California San Diego, La Jolla, CA, USA
| | - Anda Trifan
- Argonne National Lab, Lemont, IL, USA
- University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | | | | | - Austin Clyde
- Argonne National Lab, Lemont, IL, USA
- University of Chicago, Chicago, IL, USA
| | - Heng Ma
- Argonne National Lab, Lemont, IL, USA
| | | | | | | | | | | | - David J Hardy
- University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Julio Dc Maia
- University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | | | | | | | - Lei Huang
- Texas Advanced Computing Center, Austin, TX, USA
| | | | | | - Tom Gibbs
- NVIDIA Corporation, Santa Clara, CA, USA
| | - John E Stone
- University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Shantenu Jha
- Rutgers University, Piscataway, NJ, USA
- Brookhaven National Lab, Upton, NY, USA
| | | | | |
Collapse
|
2
|
Muthengi A, Wimalasena VK, Yosief HO, Bikowitz MJ, Sigua LH, Wang T, Li D, Gaieb Z, Dhawan G, Liu S, Erickson J, Amaro RE, Schönbrunn E, Qi J, Zhang W. Development of Dimethylisoxazole-Attached Imidazo[1,2- a]pyridines as Potent and Selective CBP/P300 Inhibitors. J Med Chem 2021; 64:5787-5801. [PMID: 33872011 DOI: 10.1021/acs.jmedchem.0c02232] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
The use of epigenetic bromodomain inhibitors as anticancer therapeutics has transitioned from targeting bromodomain extraterminal domain (BET) proteins into targeting non-BET bromodomains. The two most relevant non-BET bromodomain oncology targets are cyclic AMP response element-binding protein (CBP) and E1A binding protein P300 (EP300). To explore the growing CBP/EP300 interest, we developed a highly efficient two-step synthetic route for dimethylisoxazole-attached imidazo[1,2-a]pyridine scaffold-containing inhibitors. Our efficient two-step reactions enabled high-throughput synthesis of compounds designed by molecular modeling, which together with structure-activity relationship (SAR) studies facilitated an overarching understanding of selective targeting of CBP/EP300 over non-BET bromodomains. This led to the identification of a new potent and selective CBP/EP300 bromodomain inhibitor, UMB298 (compound 23, CBP IC50 72 nM and bromodomain 4, BRD4 IC50 5193 nM). The SAR we established is in good agreement with literature-reported CBP inhibitors, such as CBP30, and demonstrates the advantage of utilizing our two-step approach for inhibitor development of other bromodomains.
Collapse
Affiliation(s)
- Alex Muthengi
- Center for Green Chemistry and Department of Chemistry, University of Massachusetts Boston, 100 Morrissey Boulevard, Boston, Massachusetts 02125, United States
| | - Virangika K Wimalasena
- Department of Cancer Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, Massachusetts 02215, United States
| | - Hailemichael O Yosief
- Center for Green Chemistry and Department of Chemistry, University of Massachusetts Boston, 100 Morrissey Boulevard, Boston, Massachusetts 02125, United States
| | - Melissa J Bikowitz
- Drug Discovery Department, Moffit Cancer Center, Tampa, Florida 33612, United States.,Department of Molecular Medicine, USF Morsani College of Medicine, University of South Florida, Tampa, Florida 33612, United States
| | - Logan H Sigua
- Department of Cancer Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, Massachusetts 02215, United States
| | - Tingjian Wang
- Department of Cancer Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, Massachusetts 02215, United States
| | - Deyao Li
- Department of Cancer Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, Massachusetts 02215, United States
| | - Zied Gaieb
- Department of Chemistry & Biochemistry, University of California, San Diego, 9500 Gilman Dr, LA Jolla, California 92093, United States
| | - Gagan Dhawan
- Department of Biomedical Science, Acharya Narendra Dev College, University of Delhi, Delhi, New Delhi 110019, India
| | - Shuai Liu
- Center for Green Chemistry and Department of Chemistry, University of Massachusetts Boston, 100 Morrissey Boulevard, Boston, Massachusetts 02125, United States
| | - Jon Erickson
- Center for Green Chemistry and Department of Chemistry, University of Massachusetts Boston, 100 Morrissey Boulevard, Boston, Massachusetts 02125, United States
| | - Rommie E Amaro
- Department of Chemistry & Biochemistry, University of California, San Diego, 9500 Gilman Dr, LA Jolla, California 92093, United States
| | - Ernst Schönbrunn
- Drug Discovery Department, Moffit Cancer Center, Tampa, Florida 33612, United States.,Department of Molecular Medicine, USF Morsani College of Medicine, University of South Florida, Tampa, Florida 33612, United States
| | - Jun Qi
- Department of Cancer Biology, Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, Massachusetts 02215, United States.,Department of Medicine, Harvard Medical School, 220 Longwood Avenue, Boston, Massachusetts 02115, United States
| | - Wei Zhang
- Center for Green Chemistry and Department of Chemistry, University of Massachusetts Boston, 100 Morrissey Boulevard, Boston, Massachusetts 02125, United States
| |
Collapse
|
3
|
Yu A, Pak AJ, He P, Monje-Galvan V, Casalino L, Gaieb Z, Dommer AC, Amaro RE, Voth GA. A multiscale coarse-grained model of the SARS-CoV-2 virion. Biophys J 2021; 120:1097-1104. [PMID: 33253634 PMCID: PMC7695975 DOI: 10.1016/j.bpj.2020.10.048] [Citation(s) in RCA: 98] [Impact Index Per Article: 32.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2020] [Revised: 10/27/2020] [Accepted: 10/30/2020] [Indexed: 01/01/2023] Open
Abstract
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and ongoing development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data become publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.
Collapse
Affiliation(s)
- Alvin Yu
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois
| | - Alexander J Pak
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois
| | - Peng He
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois
| | - Viviana Monje-Galvan
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois
| | - Lorenzo Casalino
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Zied Gaieb
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Abigail C Dommer
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics and James Franck Institute, The University of Chicago, Chicago, Illinois.
| |
Collapse
|
4
|
Narkhede YB, Gautam AK, Hsu RV, Rodriguez W, Zewde NT, Harrison RES, Arantes PR, Gaieb Z, Gorham RD, Kieslich C, Morikis D, Sahu A, Palermo G. Role of Electrostatic Hotspots in the Selectivity of Complement Control Proteins Toward Human and Bovine Complement Inhibition. Front Mol Biosci 2021; 8:618068. [PMID: 33829039 PMCID: PMC8020814 DOI: 10.3389/fmolb.2021.618068] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Accepted: 01/08/2021] [Indexed: 11/13/2022] Open
Abstract
Poxviruses are dangerous pathogens, which can cause fatal infection in unvaccinated individuals. The causative agent of smallpox in humans, variola virus, is closely related to the bovine vaccinia virus, yet the molecular basis of their selectivity is currently incompletely understood. Here, we examine the role of the electrostatics in the selectivity of the smallpox protein SPICE and vaccinia protein VCP toward the human and bovine complement protein C3b, a key component of the complement immune response. Electrostatic calculations, in-silico alanine-scan and electrostatic hotspot analysis, as introduced by Kieslich and Morikis (PLoS Comput. Biol. 2012), are used to assess the electrostatic complementarity and to identify sites resistant to local perturbation where the electrostatic potential is likely to be evolutionary conserved. The calculations suggest that the bovine C3b is electrostatically prone to selectively bind its VCP ligand. On the other hand, the human isoform of C3b exhibits a lower electrostatic complementarity toward its SPICE ligand. Yet, the human C3b displays a highly preserved electrostatic core, which suggests that this isoform could be less selective in binding different ligands like SPICE and the human Factor H. This is supported by experimental cofactor activity assays revealing that the human C3b is prone to bind both SPICE and Factor H, which exhibit diverse electrostatic properties. Additional investigations considering mutants of SPICE and VCP that revert their selectivity reveal an “electrostatic switch” into the central modules of the ligands, supporting the critical role of the electrostatics in the selectivity. Taken together, these evidences provide insights into the selectivity mechanism of the complement regulator proteins encoded by the variola and vaccinia viruses to circumvent the complement immunity and exert their pathogenic action. These fundamental aspects are valuable for the development of novel vaccines and therapeutic strategies.
Collapse
Affiliation(s)
- Yogesh B Narkhede
- Department of Bioengineering, University of California, Riverside, CA
| | - Avneesh K Gautam
- National Centre for Cell Science, Pune University Campus, Ganeshkhind, India
| | - Rohaine V Hsu
- Department of Bioengineering, University of California, Riverside, CA
| | - Wilson Rodriguez
- Department of Bioengineering, University of California, Riverside, CA
| | - Nehemiah T Zewde
- Department of Bioengineering, University of California, Riverside, CA
| | - Reed E S Harrison
- Department of Bioengineering, University of California, Riverside, CA
| | - Pablo R Arantes
- Department of Bioengineering, University of California, Riverside, CA
| | - Zied Gaieb
- Department of Bioengineering, University of California, Riverside, CA
| | - Ronald D Gorham
- Department of Bioengineering, University of California, Riverside, CA
| | - Chris Kieslich
- Department of Bioengineering, University of California, Riverside, CA.,Department of Chemical Engineering, Auburn University, Auburn, AL
| | - Dimitrios Morikis
- Department of Bioengineering, University of California, Riverside, CA
| | - Arvind Sahu
- National Centre for Cell Science, Pune University Campus, Ganeshkhind, India
| | - Giulia Palermo
- Department of Bioengineering, University of California, Riverside, CA.,Department of Chemistry, University of California, Riverside, CA
| |
Collapse
|
5
|
Oliveira ASF, Ibarra AA, Bermudez I, Casalino L, Gaieb Z, Shoemark DK, Gallagher T, Sessions RB, Amaro RE, Mulholland AJ. A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors. Biophys J 2021; 120:983-993. [PMID: 33609494 PMCID: PMC7889469 DOI: 10.1016/j.bpj.2021.01.037] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Revised: 01/11/2021] [Accepted: 01/13/2021] [Indexed: 01/08/2023] Open
Abstract
Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4β2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom; Bristol Synthetic Biology Centre, BrisSynBio, Bristol, United Kingdom
| | - Amaurys Avila Ibarra
- Research Software Engineering, Advanced Computing Research Centre, University of Bristol, Bristol, United Kingdom
| | - Isabel Bermudez
- Department of Biological and Medical Sciences, Oxford Brookes University, Oxford, United Kingdom
| | - Lorenzo Casalino
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Zied Gaieb
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Deborah K Shoemark
- School of Biochemistry, University of Bristol, Bristol, United Kingdom; Bristol Synthetic Biology Centre, BrisSynBio, Bristol, United Kingdom
| | - Timothy Gallagher
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
| | | | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom.
| |
Collapse
|
6
|
Rathinaswamy MK, Gaieb Z, Fleming KD, Borsari C, Harris NJ, Moeller BE, Wymann MP, Amaro RE, Burke JE. Disease-related mutations in PI3Kγ disrupt regulatory C-terminal dynamics and reveal a path to selective inhibitors. eLife 2021; 10:e64691. [PMID: 33661099 PMCID: PMC7955810 DOI: 10.7554/elife.64691] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2020] [Accepted: 03/03/2021] [Indexed: 12/11/2022] Open
Abstract
Class I Phosphoinositide 3-kinases (PI3Ks) are master regulators of cellular functions, with the class IB PI3K catalytic subunit (p110γ) playing key roles in immune signalling. p110γ is a key factor in inflammatory diseases and has been identified as a therapeutic target for cancers due to its immunomodulatory role. Using a combined biochemical/biophysical approach, we have revealed insight into regulation of kinase activity, specifically defining how immunodeficiency and oncogenic mutations of R1021 in the C-terminus can inactivate or activate enzyme activity. Screening of inhibitors using HDX-MS revealed that activation loop-binding inhibitors induce allosteric conformational changes that mimic those in the R1021C mutant. Structural analysis of advanced PI3K inhibitors in clinical development revealed novel binding pockets that can be exploited for further therapeutic development. Overall, this work provides unique insights into regulatory mechanisms that control PI3Kγ kinase activity and shows a framework for the design of PI3K isoform and mutant selective inhibitors.
Collapse
Affiliation(s)
- Manoj K Rathinaswamy
- Department of Biochemistry and Microbiology, University of VictoriaVictoriaCanada
| | - Zied Gaieb
- Department of Chemistry and Biochemistry, University of California San DiegoSan DiegoUnited States
| | - Kaelin D Fleming
- Department of Biochemistry and Microbiology, University of VictoriaVictoriaCanada
| | - Chiara Borsari
- University of Basel, Department of BiomedicineBaselSwitzerland
| | - Noah J Harris
- Department of Biochemistry and Microbiology, University of VictoriaVictoriaCanada
| | - Brandon E Moeller
- Department of Biochemistry and Microbiology, University of VictoriaVictoriaCanada
| | | | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San DiegoSan DiegoUnited States
| | - John E Burke
- Department of Biochemistry and Microbiology, University of VictoriaVictoriaCanada
- Department of Biochemistry and Molecular Biology, The University of British ColumbiaVancouverCanada
| |
Collapse
|
7
|
Casalino L, Dommer A, Gaieb Z, Barros EP, Sztain T, Ahn SH, Trifan A, Brace A, Bogetti A, Ma H, Lee H, Turilli M, Khalid S, Chong L, Simmerling C, Hardy DJ, Maia JDC, Phillips JC, Kurth T, Stern A, Huang L, McCalpin J, Tatineni M, Gibbs T, Stone JE, Jha S, Ramanathan A, Amaro RE. AI-Driven Multiscale Simulations Illuminate Mechanisms of SARS-CoV-2 Spike Dynamics. bioRxiv 2020:2020.11.19.390187. [PMID: 33236007 PMCID: PMC7685317 DOI: 10.1101/2020.11.19.390187] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/17/2023]
Abstract
We develop a generalizable AI-driven workflow that leverages heterogeneous HPC resources to explore the time-dependent dynamics of molecular systems. We use this workflow to investigate the mechanisms of infectivity of the SARS-CoV-2 spike protein, the main viral infection machinery. Our workflow enables more efficient investigation of spike dynamics in a variety of complex environments, including within a complete SARS-CoV-2 viral envelope simulation, which contains 305 million atoms and shows strong scaling on ORNL Summit using NAMD. We present several novel scientific discoveries, including the elucidation of the spike's full glycan shield, the role of spike glycans in modulating the infectivity of the virus, and the characterization of the flexible interactions between the spike and the human ACE2 receptor. We also demonstrate how AI can accelerate conformational sampling across different systems and pave the way for the future application of such methods to additional studies in SARS-CoV-2 and other molecular systems.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Anda Trifan
- Argonne National Lab
- University of Illinois at Urbana-Champaign
| | | | | | | | - Hyungro Lee
- Rutgers University & Brookhaven National Lab
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
8
|
Barros EP, Casalino L, Gaieb Z, Dommer AC, Wang Y, Fallon L, Raguette L, Belfon K, Simmerling C, Amaro RE. The flexibility of ACE2 in the context of SARS-CoV-2 infection. Biophys J 2020; 120:1072-1084. [PMID: 33189680 PMCID: PMC7661960 DOI: 10.1016/j.bpj.2020.10.036] [Citation(s) in RCA: 78] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Revised: 10/22/2020] [Accepted: 10/27/2020] [Indexed: 12/13/2022] Open
Abstract
The coronavirus disease 2019 (COVID-19) pandemic has swept over the world in the past months, causing significant loss of life and consequences to human health. Although numerous drug and vaccine development efforts are underway, there are many outstanding questions on the mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral association to angiotensin-converting enzyme 2 (ACE2), its main host receptor, and host cell entry. Structural and biophysical studies indicate some degree of flexibility in the viral extracellular spike glycoprotein and at the receptor-binding domain (RBD)-receptor interface, suggesting a role in infection. Here, we perform explicitly solvated, all-atom, molecular dynamics simulations of the glycosylated, full-length, membrane-bound ACE2 receptor in both an apo and spike RBD-bound state to probe the intrinsic dynamics of the ACE2 receptor in the context of the cell surface. A large degree of fluctuation in the full-length structure is observed, indicating hinge bending motions at the linker region connecting the head to the transmembrane helix while still not disrupting the ACE2 homodimer or ACE2-RBD interfaces. This flexibility translates into an ensemble of ACE2 homodimer conformations that could sterically accommodate binding of the spike trimer to more than one ACE2 homodimer and suggests a mechanical contribution of the host receptor toward the large spike conformational changes required for cell fusion. This work presents further structural and functional insights into the role of ACE2 in viral infection that can potentially be exploited for the rational design of effective SARS-CoV-2 therapeutics.
Collapse
Affiliation(s)
- Emilia P Barros
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California
| | - Lorenzo Casalino
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California
| | - Zied Gaieb
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California
| | - Abigail C Dommer
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California
| | - Yuzhang Wang
- Department of Chemistry, Stony Brook University, Stony Brook, New York
| | - Lucy Fallon
- Department of Chemistry, Stony Brook University, Stony Brook, New York
| | - Lauren Raguette
- Department of Chemistry, Stony Brook University, Stony Brook, New York
| | - Kellon Belfon
- Department of Chemistry, Stony Brook University, Stony Brook, New York
| | - Carlos Simmerling
- Department of Chemistry, Stony Brook University, Stony Brook, New York; Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California.
| |
Collapse
|
9
|
Casalino L, Gaieb Z, Goldsmith JA, Hjorth CK, Dommer AC, Harbison AM, Fogarty CA, Barros EP, Taylor BC, McLellan JS, Fadda E, Amaro RE. Beyond Shielding: The Roles of Glycans in the SARS-CoV-2 Spike Protein. ACS Cent Sci 2020; 6:1722-1734. [PMID: 33140034 PMCID: PMC7523240 DOI: 10.1021/acscentsci.0c01056] [Citation(s) in RCA: 534] [Impact Index Per Article: 133.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Indexed: 05/04/2023]
Abstract
The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 28,000,000 infections and 900,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viral fusion proteins, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of the glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans and on the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift toward the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of the SARS-CoV-2 S protein, which may be exploited in the therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.
Collapse
Affiliation(s)
- Lorenzo Casalino
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
| | - Zied Gaieb
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
| | - Jory A. Goldsmith
- Department
of Molecular Biosciences, The University
of Texas at Austin, Austin, Texas 78712, United States
| | - Christy K. Hjorth
- Department
of Molecular Biosciences, The University
of Texas at Austin, Austin, Texas 78712, United States
| | - Abigail C. Dommer
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
| | - Aoife M. Harbison
- Department
of Chemistry and Hamilton Institute, Maynooth
University, Dublin, Ireland
| | - Carl A. Fogarty
- Department
of Chemistry and Hamilton Institute, Maynooth
University, Dublin, Ireland
| | - Emilia P. Barros
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
| | - Bryn C. Taylor
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
- Biomedical
Sciences Graduate Program, University of
California San Diego, La Jolla, California 92093, United States
| | - Jason S. McLellan
- Department
of Molecular Biosciences, The University
of Texas at Austin, Austin, Texas 78712, United States
| | - Elisa Fadda
- Department
of Chemistry and Hamilton Institute, Maynooth
University, Dublin, Ireland
| | - Rommie E. Amaro
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
- Biomedical
Sciences Graduate Program, University of
California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
10
|
Yu A, Pak AJ, He P, Monje-Galvan V, Casalino L, Gaieb Z, Dommer AC, Amaro RE, Voth GA. A Multiscale Coarse-grained Model of the SARS-CoV-2 Virion. bioRxiv 2020. [PMID: 33024966 DOI: 10.1101/2020.10.02.323915] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations, however, are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and on-going development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Structural data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data becomes publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion. Significance Statement This study reports the construction of a molecular model for the SARS-CoV-2 virion and details our multiscale approach towards model refinement. The resulting model and methods can be applied to and enable the simulation of SARS-CoV-2 virions.
Collapse
|
11
|
Barros EP, Casalino L, Gaieb Z, Dommer AC, Wang Y, Fallon L, Raguette L, Belfon K, Simmerling C, Amaro RE. The flexibility of ACE2 in the context of SARS-CoV-2 infection. bioRxiv 2020:2020.09.16.300459. [PMID: 32995769 PMCID: PMC7523095 DOI: 10.1101/2020.09.16.300459] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The COVID-19 pandemic has swept over the world in the past months, causing significant loss of life and consequences to human health. Although numerous drug and vaccine developments efforts are underway, many questions remain outstanding on the mechanism of SARS-CoV-2 viral association to angiotensin-converting enzyme 2 (ACE2), its main host receptor, and entry in the cell. Structural and biophysical studies indicate some degree of flexibility in the viral extracellular Spike glycoprotein and at the receptor binding domain-receptor interface, suggesting a role in infection. Here, we perform all-atom molecular dynamics simulations of the glycosylated, full-length membrane-bound ACE2 receptor, in both an apo and spike receptor binding domain (RBD) bound state, in order to probe the intrinsic dynamics of the ACE2 receptor in the context of the cell surface. A large degree of fluctuation in the full length structure is observed, indicating hinge bending motions at the linker region connecting the head to the transmembrane helix, while still not disrupting the ACE2 homodimer or ACE2-RBD interfaces. This flexibility translates into an ensemble of ACE2 homodimer conformations that could sterically accommodate binding of the spike trimer to more than one ACE2 homodimer, and suggests a mechanical contribution of the host receptor towards the large spike conformational changes required for cell fusion. This work presents further structural and functional insights into the role of ACE2 in viral infection that can be exploited for the rational design of effective SARS-CoV-2 therapeutics. STATEMENT OF SIGNIFICANCE As the host receptor of SARS-CoV-2, ACE2 has been the subject of extensive structural and antibody design efforts in aims to curtail COVID-19 spread. Here, we perform molecular dynamics simulations of the homodimer ACE2 full-length structure to study the dynamics of this protein in the context of the cellular membrane. The simulations evidence exceptional plasticity in the protein structure due to flexible hinge motions in the head-transmembrane domain linker region and helix mobility in the membrane, resulting in a varied ensemble of conformations distinct from the experimental structures. Our findings suggest a dynamical contribution of ACE2 to the spike glycoprotein shedding required for infection, and contribute to the question of stoichiometry of the Spike-ACE2 complex.
Collapse
|
12
|
Oliveira ASF, Ibarra AA, Bermudez I, Casalino L, Gaieb Z, Shoemark DK, Gallagher T, Sessions RB, Amaro RE, Mulholland AJ. Simulations support the interaction of the SARS-CoV-2 spike protein with nicotinic acetylcholine receptors. bioRxiv 2020:2020.07.16.206680. [PMID: 32743575 PMCID: PMC7386492 DOI: 10.1101/2020.07.16.206680] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Changeux et al. recently suggested that the SARS-CoV-2 spike (S) protein may interact with nicotinic acetylcholine receptors (nAChRs). Such interactions may be involved in pathology and infectivity. Here, we use molecular simulations of validated atomically detailed structures of nAChRs, and of the S protein, to investigate this 'nicotinic hypothesis'. We examine the binding of the Y674-R685 loop of the S protein to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results indicate that Y674-R685 has affinity for nAChRs and the region responsible for binding contains the PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. In particular, R682 has a key role in the stabilisation of the complexes as it forms interactions with loops A, B and C in the receptor's binding pocket. The conformational behaviour of the bound Y674-R685 region is highly dependent on the receptor subtype, adopting extended conformations in the α4β2 and α7 complexes and more compact ones when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket where it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1 and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of the simulations of the full-length S protein show that the Y674-R685 region is accessible for binding, and suggest a potential binding orientation of the S protein with nAChRs.
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | - Amaurys Avila Ibarra
- Research Software Engineering, Advanced Computing Research Centre, University of Bristol, Bristol BS1 5QD, UK
| | - Isabel Bermudez
- Department of Biological and Medical Sciences, Oxford Brookes University, Oxford OX30BP, UK
| | - Lorenzo Casalino
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, 92093 USA
| | - Zied Gaieb
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, 92093 USA
| | | | - Timothy Gallagher
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| | | | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA, 92093 USA
| | - Adrian J. Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| |
Collapse
|
13
|
Casalino L, Gaieb Z, Goldsmith JA, Hjorth CK, Dommer AC, Harbison AM, Fogarty CA, Barros EP, Taylor BC, McLellan JS, Fadda E, Amaro RE. Beyond Shielding: The Roles of Glycans in SARS-CoV-2 Spike Protein. bioRxiv 2020:2020.06.11.146522. [PMID: 32577644 PMCID: PMC7302197 DOI: 10.1101/2020.06.11.146522] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/11/2023]
Abstract
The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 15,000,000 infections and 600,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates the host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viruses, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans, and the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift towards the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of SARS-CoV-2 S protein, which may be exploited by therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.
Collapse
Affiliation(s)
- Lorenzo Casalino
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
| | - Zied Gaieb
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
| | - Jory A. Goldsmith
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, TX 78712, USA
| | - Christy K. Hjorth
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, TX 78712, USA
| | - Abigail C. Dommer
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
| | - Aoife M. Harbison
- Department of Chemistry and Hamilton Institute, Maynooth University, Dublin, Ireland
| | - Carl A. Fogarty
- Department of Chemistry and Hamilton Institute, Maynooth University, Dublin, Ireland
| | - Emilia P. Barros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
| | - Bryn C. Taylor
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
- Biomedical Sciences Graduate Program, University of California San Diego, La Jolla, CA, 92093, USA
| | | | - Elisa Fadda
- Department of Molecular Biosciences, The University of Texas at Austin, Austin, TX 78712, USA
- Department of Chemistry and Hamilton Institute, Maynooth University, Dublin, Ireland
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
- Biomedical Sciences Graduate Program, University of California San Diego, La Jolla, CA, 92093, USA
| |
Collapse
|
14
|
Parks C, Gaieb Z, Amaro RE. An Analysis of Proteochemometric and Conformal Prediction Machine Learning Protein-Ligand Binding Affinity Models. Front Mol Biosci 2020; 7:93. [PMID: 32671093 PMCID: PMC7328444 DOI: 10.3389/fmolb.2020.00093] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2020] [Accepted: 04/22/2020] [Indexed: 11/13/2022] Open
Abstract
Protein-ligand binding affinity is a key pharmacodynamic endpoint in drug discovery. Sole reliance on experimental design, make, and test cycles is costly and time consuming, providing an opportunity for computational methods to assist. Herein, we present results comparing random forest and feed-forward neural network proteochemometric models for their ability to predict pIC50 measurements for held out generic Bemis-Murcko scaffolds. In addition, we assess the ability of conformal prediction to provide calibrated prediction intervals in both a retrospective and semi-prospective test using the recently released Grand Challenge 4 data set as an external test set. In total, random forest and deep neural network proteochemometric models show quality retrospective performance but suffer in the semi-prospective setting. However, the conformal predictor prediction intervals prove to be well-calibrated both retrospectively and semi-prospectively showing that they can be used to guide hit discovery and lead optimization campaigns.
Collapse
Affiliation(s)
| | | | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA, United States
| |
Collapse
|
15
|
Bevins N, Sun S, Gaieb Z, Thorson JA, Murray SS. Comparison of commonly used solid tumor targeted gene sequencing panels for estimating tumor mutation burden shows analytical and prognostic concordance within the cancer genome atlas cohort. J Immunother Cancer 2020; 8:e000613. [PMID: 32217764 PMCID: PMC7174068 DOI: 10.1136/jitc-2020-000613] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/27/2020] [Indexed: 12/22/2022] Open
Abstract
BACKGROUND Tumor mutation burden (TMB) is a biomarker frequently reported by clinical laboratories, which is derived by quantifying of the number of single nucleotide or indel variants (mutations) identified by next-generation sequencing of tumors. TMB values can inform prognosis or predict the response of a patient's tumor to immune checkpoint inhibitor therapy. Methods for the calculation of TMB are not standardized between laboratories, with significant variables being the gene content of the panels sequenced and the inclusion or exclusion of synonymous variants in the calculations. The impact of these methodological differences has not been investigated and the concordance of reported TMB values between laboratories is unknown. METHODS Sequence variant lists from more than 9000 tumors of various types were downloaded from The Cancer Genome Atlas. Variant lists were filtered to include only appropriate variant types (ie, non-synonymous only or synonymous and non-synonymous variants) within the genes found in five commonly used targeted solid tumor gene panels as well as an in-house gene panel. Calculated TMB was paired with corresponding overall survival (OS) data of each patient. RESULTS Regression analysis indicates high concordance of TMB as derived from the examined panels. TMB derived from panels was consistently and significantly lower than that derived from a whole exome. TMB, as derived from whole exome or the examined panels, showed a significant correlation with OS in the examined data. CONCLUSIONS TMB derived from the examined gene panels was analytically equivalent between panels, but not between panels and whole-exome sequencing. Correlation between TMB and OS is significant if TMB method-specific cut-offs are used. These results suggest that TMB values, as derived from the gene panels examined, are analytically and prognostically equivalent.
Collapse
Affiliation(s)
| | - Shulei Sun
- Pathology, UC San Diego, La Jolla, California, USA
| | - Zied Gaieb
- Chemistry and Biochemistry, UC San Diego, La Jolla, California, USA
| | | | | |
Collapse
|
16
|
Parks CD, Gaieb Z, Chiu M, Yang H, Shao C, Walters WP, Jansen JM, McGaughey G, Lewis RA, Bembenek SD, Ameriks MK, Mirzadegan T, Burley SK, Amaro RE, Gilson MK. D3R grand challenge 4: blind prediction of protein-ligand poses, affinity rankings, and relative binding free energies. J Comput Aided Mol Des 2020; 34:99-119. [PMID: 31974851 PMCID: PMC7261493 DOI: 10.1007/s10822-020-00289-y] [Citation(s) in RCA: 60] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 01/13/2020] [Indexed: 12/11/2022]
Abstract
The Drug Design Data Resource (D3R) aims to identify best practice methods for computer aided drug design through blinded ligand pose prediction and affinity challenges. Herein, we report on the results of Grand Challenge 4 (GC4). GC4 focused on proteins beta secretase 1 and Cathepsin S, and was run in an analogous manner to prior challenges. In Stage 1, participant ability to predict the pose and affinity of BACE1 ligands were assessed. Following the completion of Stage 1, all BACE1 co-crystal structures were released, and Stage 2 tested affinity rankings with co-crystal structures. We provide an analysis of the results and discuss insights into determined best practice methods.
Collapse
Affiliation(s)
- Conor D Parks
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Zied Gaieb
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Michael Chiu
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Huanwang Yang
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Chenghua Shao
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | | | - Johanna M Jansen
- Novartis Institutes for BioMedical Research, Emeryville, CA, 94608, USA
| | | | - Richard A Lewis
- Novartis Institutes for BioMedical Research, Novartis Pharma AG, 4002, Basel, Switzerland
| | | | | | | | - Stephen K Burley
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Rommie E Amaro
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA.
- Department of Chemistry and Biochemistry, UC San Diego, La Jolla, CA, 92093-0340, USA.
| | - Michael K Gilson
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA.
- Skaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, 9500 Gilman Drive, MC0751, La Jolla, CA, 92093, USA.
| |
Collapse
|
17
|
Hwang SH, Gaieb Z, Morisseau C, Yang J, Wagner K, Gilson MK, Hammock BD. Dual Inhibitors of Cyclooxygenase‐2 and Soluble Epoxide Hydrolase: Studies of Binding Modes at the Active Sites and Time‐dependency of Inhibition, and Development of Water‐soluble Prodrugs. FASEB J 2018. [DOI: 10.1096/fasebj.2018.32.1_supplement.558.2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Sung Hee Hwang
- Entomology & NematologyUniversity of CaliforniaDavisDavisCA
| | - Zied Gaieb
- Skaggs school of pharmacy and pharmaceutical sciencesUniversity of CaliforniaSan DiegoLa JollaCA
| | | | - Jun Yang
- Entomology & NematologyUniversity of CaliforniaDavisDavisCA
| | - Karen Wagner
- Entomology & NematologyUniversity of CaliforniaDavisDavisCA
| | - Michael K. Gilson
- Skaggs school of pharmacy and pharmaceutical sciencesUniversity of CaliforniaSan DiegoLa JollaCA
| | | |
Collapse
|
18
|
Gaieb Z, Liu S, Gathiaka S, Chiu M, Yang H, Shao C, Feher VA, Walters P, Kuhn B, Rudolph MG, Burley SK, Gilson MK, Amaro RE. D3R Grand Challenge 2: blind prediction of protein-ligand poses, affinity rankings, and relative binding free energies. J Comput Aided Mol Des 2018; 32:1-20. [PMID: 29204945 PMCID: PMC5767524 DOI: 10.1007/s10822-017-0088-4] [Citation(s) in RCA: 129] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2017] [Accepted: 11/20/2017] [Indexed: 01/07/2023]
Abstract
The Drug Design Data Resource (D3R) ran Grand Challenge 2 (GC2) from September 2016 through February 2017. This challenge was based on a dataset of structures and affinities for the nuclear receptor farnesoid X receptor (FXR), contributed by F. Hoffmann-La Roche. The dataset contained 102 IC50 values, spanning six orders of magnitude, and 36 high-resolution co-crystal structures with representatives of four major ligand classes. Strong global participation was evident, with 49 participants submitting 262 prediction submission packages in total. Procedurally, GC2 mimicked Grand Challenge 2015 (GC2015), with a Stage 1 subchallenge testing ligand pose prediction methods and ranking and scoring methods, and a Stage 2 subchallenge testing only ligand ranking and scoring methods after the release of all blinded co-crystal structures. Two smaller curated sets of 18 and 15 ligands were developed to test alchemical free energy methods. This overview summarizes all aspects of GC2, including the dataset details, challenge procedures, and participant results. We also consider implications for progress in the field, while highlighting methodological areas that merit continued development. Similar to GC2015, the outcome of GC2 underscores the pressing need for methods development in pose prediction, particularly for ligand scaffolds not currently represented in the Protein Data Bank ( http://www.pdb.org ), and in affinity ranking and scoring of bound ligands.
Collapse
Affiliation(s)
- Zied Gaieb
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA 92093
| | - Shuai Liu
- Silicon Therapeutics, Boston MA 02210
| | | | - Michael Chiu
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA 92093
| | - Huanwang Yang
- RCSB Protein Data Bank Rutgers University, New Brunswick, NJ 08901
| | - Chenghua Shao
- RCSB Protein Data Bank Rutgers University, New Brunswick, NJ 08901
| | - Victoria A. Feher
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA 92093
| | | | - Bernd Kuhn
- Roche Pharmaceutical Research and Early Development (pRED), Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, 4070 Basel, Switzerland
| | - Markus G. Rudolph
- Roche Pharmaceutical Research and Early Development (pRED), Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Grenzacherstrasse 124, 4070 Basel, Switzerland
| | | | - Michael K. Gilson
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA 92093,Correspondence to:
| | - Rommie E. Amaro
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA 92093,Correspondence to:
| |
Collapse
|
19
|
Gaieb Z, Morikis D. Detection of Side Chain Rearrangements Mediating the Motions of Transmembrane Helices in Molecular Dynamics Simulations of G Protein-Coupled Receptors. Comput Struct Biotechnol J 2017; 15:131-137. [PMID: 28149485 PMCID: PMC5271675 DOI: 10.1016/j.csbj.2017.01.001] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2016] [Revised: 01/03/2017] [Accepted: 01/10/2017] [Indexed: 12/02/2022] Open
Abstract
Structure and dynamics are essential elements of protein function. Protein structure is constantly fluctuating and undergoing conformational changes, which are captured by molecular dynamics (MD) simulations. We introduce a computational framework that provides a compact representation of the dynamic conformational space of biomolecular simulations. This method presents a systematic approach designed to reduce the large MD simulation spatiotemporal datasets into a manageable set in order to guide our understanding of how protein mechanics emerge from side chain organization and dynamic reorganization. We focus on the detection of side chain interactions that undergo rearrangements mediating global domain motions and vice versa. Side chain rearrangements are extracted from side chain interactions that undergo well-defined abrupt and persistent changes in distance time series using Gaussian mixture models, whereas global domain motions are detected using dynamic cross-correlation. Both side chain rearrangements and global domain motions represent the dynamic components of the protein MD simulation, and are both mapped into a network where they are connected based on their degree of coupling. This method allows for the study of allosteric communication in proteins by mapping out the protein dynamics into an intramolecular network to reduce the large simulation data into a manageable set of communities composed of coupled side chain rearrangements and global domain motions. This computational framework is suitable for the study of tightly packed proteins, such as G protein-coupled receptors, and we present an application on a seven microseconds MD trajectory of CC chemokine receptor 7 (CCR7) bound to its ligand CCL21.
Collapse
Affiliation(s)
- Zied Gaieb
- Department of Bioengineering, University of California, Riverside 92521, USA
| | - Dimitrios Morikis
- Department of Bioengineering, University of California, Riverside 92521, USA
| |
Collapse
|
20
|
Abstract
Biased ligand binding to G protein-coupled receptors enables functional selectivity of intracellular effectors to mediate cellular function. Despite the significant advances made in characterizing the conformational states (transmembrane helical arrangements) capable of discriminating between G protein and arrestin binding, the role of the ligand in stabilizing such conformations remains unclear. To address this issue, we simulate microsecond dynamics of CC chemokine receptor 7 (CCR7) bound to its native biased ligands, CCL19 and CCL21, and detect a series of molecular switches that are mediated by various ligand-induced allosteric events. These molecular switches involve three tyrosine residues (Y112(3.32), Y255(6.51), and Y288(7.39)), three phenylalanine residues (F116(3.36), F208(5.47), and F248(6.44)), and a polar interaction between Q252(6.48) and R294(7.45) in the transmembrane domain of CCR7. Conformational changes within these switches, particularly hydrogen bond formation between Y112(3.32) and Y255(6.51), lead to global helical movements in the receptor's transmembrane helices and contribute to the transitioning of the receptor to distinct states. Ligand-induced helical movements in the receptor highlight the ability of biased ligands to stabilize the receptor in different states through a dynamic network of allosteric events.
Collapse
Affiliation(s)
- Zied Gaieb
- Department of Bioengineering, ‡Division of Biomedical Sciences, School of Medicine, University of California , Riverside, California 92521, United States
| | - David D Lo
- Department of Bioengineering, ‡Division of Biomedical Sciences, School of Medicine, University of California , Riverside, California 92521, United States
| | - Dimitrios Morikis
- Department of Bioengineering, ‡Division of Biomedical Sciences, School of Medicine, University of California , Riverside, California 92521, United States
| |
Collapse
|
21
|
Gaieb Z, Lo DD, Morikis D. Ligand-Specific Conformational Changes in CCR7 Coupled to Selecting Different Signaling Pathways upon CCL19 And CCL21 Ligand Binding. Biophys J 2016. [DOI: 10.1016/j.bpj.2015.11.1947] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
|