1
|
Jin H, Merz KM. Modeling Zinc Complexes Using Neural Networks. J Chem Inf Model 2024; 64:3140-3148. [PMID: 38587510 PMCID: PMC11040731 DOI: 10.1021/acs.jcim.4c00095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Revised: 03/04/2024] [Accepted: 03/28/2024] [Indexed: 04/09/2024]
Abstract
Understanding the energetic landscapes of large molecules is necessary for the study of chemical and biological systems. Recently, deep learning has greatly accelerated the development of models based on quantum chemistry, making it possible to build potential energy surfaces and explore chemical space. However, most of this work has focused on organic molecules due to the simplicity of their electronic structures as well as the availability of data sets. In this work, we build a deep learning architecture to model the energetics of zinc organometallic complexes. To achieve this, we have compiled a configurationally and conformationally diverse data set of zinc complexes using metadynamics to overcome the limitations of traditional sampling methods. In terms of the neural network potentials, our results indicate that for zinc complexes, partial charges play an important role in modeling the long-range interactions with a neural network. Our developed model outperforms semiempirical methods in predicting the relative energy of zinc conformers, yielding a mean absolute error (MAE) of 1.32 kcal/mol with reference to the double-hybrid PWPB95 method.
Collapse
Affiliation(s)
- Hongni Jin
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Kenneth M. Merz
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department
of Biochemistry and Molecular Biology, Michigan
State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
2
|
Mozumder S, Bej A, Sengupta J. Ligand-Dependent Modulation of the Dynamics of Intracellular Loops Dictates Functional Selectivity of 5-HT 2AR. J Chem Inf Model 2022; 62:2522-2537. [PMID: 35324173 DOI: 10.1021/acs.jcim.2c00118] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The serotonin 2A receptor (5-HT2AR) subtype of the G protein-coupled receptor (GPCR) family is involved in a plethora of neuromodulatory functions (e.g., neurogenesis, sleep, and cognitive processes). 5-HT2AR is the target of pharmacologically distinct classes of ligands, binding of which either activate or inactivate the receptor. Although high-resolution structures of 5-HT2AR as well as several other 5-HT GPCRs provided snapshots of both active and inactive conformational states, these structures, representing a truncated form of the receptor, cannot fully explain the mechanism of conformational transitions during their function. Importantly, biochemical studies have suggested the importance of intracellular loops in receptor functions. In our previous study, a model of the ligand-free form of 5-HT2AR with the third intracellular loop (ICL3) has been meticulously built. Here, we have investigated the functional regulation of 5-HT2AR with intact intracellular loops in ligand-free and five distinct ligand-bound configurations using unbiased atomistic molecular dynamics (MD) simulations. The selected ligands belong to either of the full, partial, or inverse agonist classes, which exert distinct pharmacological responses. We have observed significant structural, dynamic, and thermodynamic differences within ligand-bound complexes. Our results revealed, for the first time, that either activation or inactivation of the receptor upon specific ligand binding is primarily achieved through conformational transitions of its second and third intracellular loops (ICL2 and ICL3). A remarkable allosteric cross-talk between the ligand-binding site and the distal intracellular parts of the receptor, where binding of a specific ligand thermodynamically controls (either stabilizes or destabilizes) the intracellular region, consisting of crucial dynamic elements ICL2 and ICL3, and differential conformational transitions of these loops determine ligand-dependent functional selectivity.
Collapse
Affiliation(s)
- Sukanya Mozumder
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, 4, Raja S.C. Mullick Road, Jadavpur, Kolkata 700032, India.,Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| | - Aritra Bej
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, 4, Raja S.C. Mullick Road, Jadavpur, Kolkata 700032, India
| | - Jayati Sengupta
- Structural Biology and Bioinformatics Division, CSIR-Indian Institute of Chemical Biology, 4, Raja S.C. Mullick Road, Jadavpur, Kolkata 700032, India.,Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| |
Collapse
|
3
|
Almeida BC, Kaczmarek JA, Figueiredo PR, Prather KLJ, Carvalho ATP. Transcription factor allosteric regulation through substrate coordination to zinc. NAR Genom Bioinform 2021; 3:lqab033. [PMID: 33987533 PMCID: PMC8092373 DOI: 10.1093/nargab/lqab033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2020] [Revised: 03/30/2021] [Accepted: 04/08/2021] [Indexed: 11/14/2022] Open
Abstract
The development of new synthetic biology circuits for biotechnology and medicine requires deeper mechanistic insight into allosteric transcription factors (aTFs). Here we studied the aTF UxuR, a homodimer of two domains connected by a highly flexible linker region. To explore how ligand binding to UxuR affects protein dynamics we performed molecular dynamics simulations in the free protein, the aTF bound to the inducer D-fructuronate or the structural isomer D-glucuronate. We then validated our results by constructing a sensor plasmid for D-fructuronate in Escherichia coli and performed site-directed mutagenesis. Our results show that zinc coordination is necessary for UxuR function since mutation to alanines prevents expression de-repression by D-fructuronate. Analyzing the different complexes, we found that the disordered linker regions allow the N-terminal domains to display fast and large movements. When the inducer is bound, UxuR can sample an open conformation with a more pronounced negative charge at the surface of the N-terminal DNA binding domains. In opposition, in the free and D-glucuronate bond forms the protein samples closed conformations, with a more positive character at the surface of the DNA binding regions. These molecular insights provide a new basis to harness these systems for biological systems engineering.
Collapse
Affiliation(s)
- Beatriz C Almeida
- CNC-Center for Neuroscience and Cell Biology, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, 3004-504 Coimbra, Portugal
| | - Jennifer A Kaczmarek
- MIT-Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Pedro R Figueiredo
- CNC-Center for Neuroscience and Cell Biology, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, 3004-504 Coimbra, Portugal
| | - Kristala L J Prather
- MIT-Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Alexandra T P Carvalho
- CNC-Center for Neuroscience and Cell Biology, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, 3004-504 Coimbra, Portugal
| |
Collapse
|
4
|
Tian W, Liu X, Wang L, Zheng B, Jiang K, Fu G, Feng W. Deciphering the selective binding mechanisms of anaplastic lymphoma kinase-derived neuroblastoma tumor neoepitopes to human leukocyte antigen. J Mol Model 2021; 27:134. [PMID: 33899124 DOI: 10.1007/s00894-021-04754-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2020] [Accepted: 04/15/2021] [Indexed: 12/14/2022]
Abstract
Neuroblastoma (NB), as a metastatic form of solid tumor, has a high fatality rate found in early childhood. The two anaplastic lymphoma kinase (ALK) neoepitopes nonamer and decamer used in cancer immunotherapy against NB cancer can selectively bind to the human leukocyte antigen (HLA-B*15:01) groove with high affinities, whereas the native self-peptide is unable to interact with the HLA-B*15:01. Here, we performed molecular dynamics (MD) simulations and subsequent molecular mechanics-generalized born surface area (MM-GBSA) binding free energy calculations to explore the selective binding mechanisms of nonamer and decamer to the HLA-B*15:01 against the self-peptide. MD simulations revealed the significant conformational dynamics of the self-peptide in the HLA-B*15:01 groove against the nonamer and decamer. Binding free energy calculations showed that the binding affinities of HLA-B*15:01-neoepitope complexes were followed in the order decamer > nonamer > self-peptide. Detailed analysis of HLA-B*15:01-neoepitope structural complexes showed that compared to the nonamer, the self-peptide tended to move outward to the solvent, whereas the decamer moved deep to the HLA-B*15:01 groove. These different dynamic observations of the ALK neoepitopes can explain the distinct binding affinities of self-peptide, nonamer, and decamer to the HLA-B*15:01. The results may be useful for the design of more selective ALK neoepitopes.
Collapse
Affiliation(s)
- Wenchao Tian
- Pediatric Surgery, Binzhou Medical University Hospital, Binzhou, 256603, Shandong, China
| | - Xianxian Liu
- Department of Infectious Diseases, Binzhou Medical University Hospital, Binzhou, 256603, Shandong, China
| | - Lulu Wang
- Department of Critical Care Medicine, Binzhou Medical University Hospital, Binzhou, 256603, Shandong, China
| | - Bufeng Zheng
- Pediatric Surgery, Binzhou Medical University Hospital, Binzhou, 256603, Shandong, China
| | - Kun Jiang
- Pediatric Surgery, Binzhou Medical University Hospital, Binzhou, 256603, Shandong, China
| | - Guoyong Fu
- Pediatric Surgery, Binzhou Medical University Hospital, Binzhou, 256603, Shandong, China
| | - Wenyu Feng
- Pediatric Surgery, Binzhou Medical University Hospital, Binzhou, 256603, Shandong, China.
| |
Collapse
|
5
|
Li X, Ye M, Wang Y, Qiu M, Fu T, Zhang J, Zhou B, Lu S. How Parkinson's disease-related mutations disrupt the dimerization of WD40 domain in LRRK2: a comparative molecular dynamics simulation study. Phys Chem Chem Phys 2021; 22:20421-20433. [PMID: 32914822 DOI: 10.1039/d0cp03171b] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The multidomain kinase enzyme leucine-rich-repeat kinase 2 (LRRK2), activated through a homodimerization manner, has been identified as an important pathogenic factor in Parkinson's disease (PD), the second most common neurodegenerative disease wordwide. The Trp-Asp-40 (WD40) domain, located in the C-terminal LRRK2, harbours one of the most frequent PD-related variants, G2385R. However, the detailed dynamics of WD40 during LRRK2 dimerization and the underlying mechanism through which the pathogenic mutations disrupt the formation of the WD40 dimer have remained elusive. Here, microsecond-scale molecular dynamics simulations were employed to provide a mechanistic view underlying the WD40 dimerization and unveil the structural basis by which the interface-based mutations G2385R, H2391D and R2394E compromise the corresponding process. The simulation results identified important residues, D2351, R2394, E2395, R2413, and R2443, involved in establishing the complex binding network along the dimerization interface, which was significantly weakened in the presence of interfacial mutations. A "sag-bulge" model was proposed to explain the unfavorable dimer formation in the mutant systems. In addition, mutations altered the community configuration in the wild-type system in which inter-monomeric interplay is prominent, leading to the destabilization of the WD40 dimer under mutation.
Collapse
Affiliation(s)
- Xinyi Li
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai, 200025, China.
| | - Mingyu Ye
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai, 200025, China.
| | - Yue Wang
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai, 200025, China.
| | - Ming Qiu
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai, 200025, China.
| | - Tingting Fu
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai, 200025, China.
| | - Jian Zhang
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai, 200025, China.
| | - Bin Zhou
- Department of Emergency, Changhai Hospital, Affiliated to Navy Military Medical University, Shanghai, 200433, China.
| | - Shaoyong Lu
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University, School of Medicine, Shanghai, 200025, China.
| |
Collapse
|