1
|
Kairys V, Baranauskiene L, Kazlauskiene M, Zubrienė A, Petrauskas V, Matulis D, Kazlauskas E. Recent advances in computational and experimental protein-ligand affinity determination techniques. Expert Opin Drug Discov 2024; 19:649-670. [PMID: 38715415 DOI: 10.1080/17460441.2024.2349169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2024] [Accepted: 04/25/2024] [Indexed: 05/22/2024]
Abstract
INTRODUCTION Modern drug discovery revolves around designing ligands that target the chosen biomolecule, typically proteins. For this, the evaluation of affinities of putative ligands is crucial. This has given rise to a multitude of dedicated computational and experimental methods that are constantly being developed and improved. AREAS COVERED In this review, the authors reassess both the industry mainstays and the newest trends among the methods for protein - small-molecule affinity determination. They discuss both computational affinity predictions and experimental techniques, describing their basic principles, main limitations, and advantages. Together, this serves as initial guide to the currently most popular and cutting-edge ligand-binding assays employed in rational drug design. EXPERT OPINION The affinity determination methods continue to develop toward miniaturization, high-throughput, and in-cell application. Moreover, the availability of data analysis tools has been constantly increasing. Nevertheless, cross-verification of data using at least two different techniques and careful result interpretation remain of utmost importance.
Collapse
Affiliation(s)
- Visvaldas Kairys
- Department of Bioinformatics, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Lina Baranauskiene
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | | | - Asta Zubrienė
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Vytautas Petrauskas
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Daumantas Matulis
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Egidijus Kazlauskas
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| |
Collapse
|
2
|
Avenali M, Cerri S, Palmieri I, Ongari G, Stiuso R, Buongarzone G, Tassorelli C, Biagini T, Valente M, Cereda C, Mazza T, Gana S, Pacchetti C, Valente EM. Functional Study of SNCA p.V15A Variant: Further Linking α-Synuclein and Glucocerebrosidase. Mov Disord 2024; 39:1060-1065. [PMID: 38436488 DOI: 10.1002/mds.29736] [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: 06/20/2023] [Revised: 12/16/2023] [Accepted: 01/22/2024] [Indexed: 03/05/2024] Open
Abstract
BACKGROUND SNCA p.V15A was reported in five families. In vitro models showed increased aggregation and seeding activity, mitochondrial damage, and apoptosis. Mutant flies had reduced flying ability and survival. OBJECTIVES To clinically and functionally evaluate SNCA p.V15A in a large Italian family with Parkinson's disease (PD). METHODS Genetic diagnosis was reached through next-generation sequencing. Pathogenicity was assessed by molecular dynamics simulation and biochemical studies on peripheral blood mononuclear cells (PBMCs). RESULTS Five siblings carried SNCA p.V15A; three developed bradykinetic-rigid PD in their 50s with rapid motor progression and variable cognitive impairment. A fourth sibling had isolated mood disturbance, whereas the fifth was still unaffected at age 47. The mutant protein showed decreased stability and an unstable folded structure. Proband's PBMCs showed elevated total and phosphorylated α-synuclein (α-syn) levels and significantly reduced glucocerebrosidase activity. CONCLUSION This study demonstrates accumulation of α-synV15A in PBMCs and strengthens the link between α-syn pathophysiology and glucocerebrosidase dysfunction. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.
Collapse
Affiliation(s)
- Micol Avenali
- Department of Brain and Behavioural Sciences, University of Pavia, Pavia, Italy
- Parkinson's Disease and Movement Disorders Unit, IRCCS Mondino Foundation, Pavia, Italy
| | - Silvia Cerri
- Cellular and Molecular Neurobiology Section, IRCCS Mondino Foundation, Pavia, Italy
| | | | - Gerardo Ongari
- Cellular and Molecular Neurobiology Section, IRCCS Mondino Foundation, Pavia, Italy
| | - Rita Stiuso
- Cellular and Molecular Neurobiology Section, IRCCS Mondino Foundation, Pavia, Italy
| | | | - Cristina Tassorelli
- Department of Brain and Behavioural Sciences, University of Pavia, Pavia, Italy
- Neurorehabilitation Unit, IRCCS Mondino Foundation, Pavia, Italy
| | - Tommaso Biagini
- Bioinformatics Unit, IRCCS Casa Sollievo della Sofferenza, San Giovanni Rotondo, Italy
| | - Marialuisa Valente
- Clinical Pathology Unit, Medical Genetics Section, SS. Annunziata Hospital, ASL Taranto, Taranto, Italy
| | - Cristina Cereda
- Center of Functional Genomics and Rare Diseases, Department of Pediatrics, Buzzi Children's Hospital, Milan, Italy
| | - Tommaso Mazza
- Bioinformatics Unit, IRCCS Casa Sollievo della Sofferenza, San Giovanni Rotondo, Italy
| | - Simone Gana
- Neurogenetics Unit, IRCCS Mondino Foundation, Pavia, Italy
| | - Claudio Pacchetti
- Parkinson's Disease and Movement Disorders Unit, IRCCS Mondino Foundation, Pavia, Italy
| | - Enza Maria Valente
- Neurogenetics Unit, IRCCS Mondino Foundation, Pavia, Italy
- Department of Molecular Medicine, University of Pavia, Pavia, Italy
| |
Collapse
|
3
|
Chen J, Song Y, Ma L, Jin Y, Yu J, Guo Y, Huang Y, Pu X. Computational insights into diverse binding modes of the allosteric modulator and their regulation on dopamine D1 receptor. Comput Biol Med 2024; 173:108283. [PMID: 38552278 DOI: 10.1016/j.compbiomed.2024.108283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2024] [Revised: 02/29/2024] [Accepted: 03/12/2024] [Indexed: 04/17/2024]
Abstract
Allosteric drugs hold the promise of addressing many challenges in the current drug development of GPCRs. However, the molecular mechanism underlying their allosteric modulations remain largely elusive. The dopamine D1 receptor (DRD1), a member of Class A GPCRs, is critical for treating psychiatric disorders, and LY3154207 serves as its promising positive allosteric modulator (PAM). In the work, we utilized extensive Gaussian-accelerated molecular dynamics simulations (a total of 41μs) for the first time probe the diverse binding modes of the allosteric modulator and their regulation effects, based on the DRD1 and LY3154207 as representative. Our simulations identify four binding modes of LY3154207 (one boat mode, two metastable vertical modes and a novel cleft-anchored mode), in which the boat mode is the most stable while there three modes are similar in the stability. However, it is interesting to observed that the most stable boat mode inversely exhibits the weakest positive allosteric effect on influencing the orthosteric ligand binding and maintaining the activity of the transducer binding site. It should result from its induced weaker correlation between the allosteric site and the orthosteric site, and between the orthosteric site and the transducer binding site than the other three binding modes, as well as its weakened interaction between a crucial activation-related residue (S2025.46) and the orthosteric ligand (dopamine). Overall, the work offers atomic-level information to advance our understanding of the complex allosteric regulation on GPCRs, which is beneficial to the allosteric modulator design and development.
Collapse
Affiliation(s)
- Jianfang Chen
- College of Chemistry, Sichuan University, Chengdu, 610064, China.
| | - Yuanpeng Song
- College of Chemistry, Sichuan University, Chengdu, 610064, China.
| | - Luhan Ma
- College of Chemistry, Sichuan University, Chengdu, 610064, China.
| | - Yizhou Jin
- College of Chemistry, Sichuan University, Chengdu, 610064, China.
| | - Jin Yu
- Department of Physics and Astronomy, University of California, Irvine, CA, 92697, USA.
| | - Yanzhi Guo
- College of Chemistry, Sichuan University, Chengdu, 610064, China.
| | - Yan Huang
- College of Chemistry, Sichuan University, Chengdu, 610064, China.
| | - Xuemei Pu
- College of Chemistry, Sichuan University, Chengdu, 610064, China.
| |
Collapse
|
4
|
Peng C, Lv X, Zhang Z, Lin J, Li D. The Recognition Pathway of the SARS-CoV-2 Spike Receptor-Binding Domain to Human Angiotensin-Converting Enzyme 2. Molecules 2024; 29:1875. [PMID: 38675695 PMCID: PMC11054751 DOI: 10.3390/molecules29081875] [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: 03/11/2024] [Revised: 04/16/2024] [Accepted: 04/17/2024] [Indexed: 04/28/2024] Open
Abstract
COVID-19 caused by SARS-CoV-2 has spread around the world. The receptor-binding domain (RBD) of the spike protein of SARS-CoV-2 is a critical component that directly interacts with host ACE2. Here, we simulate the ACE2 recognition processes of RBD of the WT, Delta, and OmicronBA.2 variants using our recently developed supervised Gaussian accelerated molecular dynamics (Su-GaMD) approach. We show that RBD recognizes ACE2 through three contact regions (regions I, II, and III), which aligns well with the anchor-locker mechanism. The higher binding free energy in State d of the RBDOmicronBA.2-ACE2 system correlates well with the increased infectivity of OmicronBA.2 in comparison with other variants. For RBDDelta, the T478K mutation affects the first step of recognition, while the L452R mutation, through its nearby Y449, affects the RBDDelta-ACE2 binding in the last step of recognition. For RBDOmicronBA.2, the E484A mutation affects the first step of recognition, the Q493R, N501Y, and Y505H mutations affect the binding free energy in the last step of recognition, mutations in the contact regions affect the recognition directly, and other mutations indirectly affect recognition through dynamic correlations with the contact regions. These results provide theoretical insights for RBD-ACE2 recognition and may facilitate drug design against SARS-CoV-2.
Collapse
Affiliation(s)
- Can Peng
- State Key Laboratory of Medicinal Chemical Biology, College of Pharmacy and Tianjin Key Laboratory of Molecular Drug Research, Nankai University, Haihe Education Park, 38 Tongyan Road, Tianjin 300350, China; (C.P.); (X.L.)
| | - Xinyue Lv
- State Key Laboratory of Medicinal Chemical Biology, College of Pharmacy and Tianjin Key Laboratory of Molecular Drug Research, Nankai University, Haihe Education Park, 38 Tongyan Road, Tianjin 300350, China; (C.P.); (X.L.)
| | - Zhiqiang Zhang
- Xiongan Institute of Innovation, Xiong’an New Area 070001, China;
| | - Jianping Lin
- State Key Laboratory of Medicinal Chemical Biology, College of Pharmacy and Tianjin Key Laboratory of Molecular Drug Research, Nankai University, Haihe Education Park, 38 Tongyan Road, Tianjin 300350, China; (C.P.); (X.L.)
| | - Dongmei Li
- State Key Laboratory of Medicinal Chemical Biology, College of Pharmacy and Tianjin Key Laboratory of Molecular Drug Research, Nankai University, Haihe Education Park, 38 Tongyan Road, Tianjin 300350, China; (C.P.); (X.L.)
| |
Collapse
|
5
|
Valdés-Albuernes JL, Díaz-Pico E, Alfaro S, Caballero J. Modeling of noncovalent inhibitors of the papain-like protease (PLpro) from SARS-CoV-2 considering the protein flexibility by using molecular dynamics and cross-docking. Front Mol Biosci 2024; 11:1374364. [PMID: 38601323 PMCID: PMC11004324 DOI: 10.3389/fmolb.2024.1374364] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Accepted: 03/11/2024] [Indexed: 04/12/2024] Open
Abstract
The papain-like protease (PLpro) found in coronaviruses that can be transmitted from animals to humans is a critical target in respiratory diseases linked to Severe Acute Respiratory Syndrome (SARS-CoV). Researchers have proposed designing PLpro inhibitors. In this study, a set of 89 compounds, including recently reported 2-phenylthiophenes with nanomolar inhibitory potency, were investigated as PLpro noncovalent inhibitors using advanced molecular modeling techniques. To develop the work with these inhibitors, multiple structures of the SARS-CoV-2 PLpro binding site were generated using a molecular sampling method. These structures were then clustered to select a group that represents the flexibility of the site. Subsequently, models of the protein-ligand complexes were created for the set of inhibitors within the chosen conformations. The quality of the complex models was assessed using LigRMSD software to verify similarities in the orientations of the congeneric series and interaction fingerprints to determine the recurrence of chemical interactions. With the multiple models constructed, a protocol was established to choose one per ligand, optimizing the correlation between the calculated docking energy values and the biological activities while incorporating the effect of the binding site's flexibility. A strong correlation (R2 = 0.922) was found when employing this flexible docking protocol.
Collapse
Affiliation(s)
| | | | | | - Julio Caballero
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, Talca, Chile
| |
Collapse
|
6
|
Wen K, Tao Y, Jiang W, Jiang L, Zhu J, Li Q. (De)carboxylation mechanisms of heteroaromatic substrates catalyzed by prenylated FMN-dependent UbiD decarboxylases: An in-silico study. Int J Biol Macromol 2024; 260:129294. [PMID: 38211929 DOI: 10.1016/j.ijbiomac.2024.129294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2023] [Revised: 12/24/2023] [Accepted: 01/05/2024] [Indexed: 01/13/2024]
Abstract
The UbiD enzymes are proposed to catalyze reversible (de)carboxylation reaction of unsaturated carboxylic acids using prenylated flavin mononucleotide (prFMN) as a cofactor. This positions UbiD enzymes as promising candidates for converting CO2 into valuable chemicals. However, their industrial-scale biotransformation is currently constrained by low conversion rates attributed to thermodynamic limitations. To enhance the carboxylation activity of UbiD enzymes, a molecular-level understanding of the (de)carboxylation mechanisms is necessary. In this study, we investigated the reaction mechanisms of heteroaromatic substrates catalyzed by PtHmfF, PaHudA, and AnlnD enzymes using molecular dynamics (MD) simulations and free energy calculations. Our extensive mechanistic study elucidates the mechanisms involved in the formation of the initial prFMN-substrate intermediate. Specifically, we observed nucleophilic attack during decarboxylation, while carboxylation reactions involving furoic acid, pyrrole, and indole tend to favor a 1,3-dipolar cycloaddition mechanism. Furthermore, we identified proton transfer as the rate-limiting step in the carboxylation reaction. In addition, we considered the perspectives of reaction energies and electron transfer to understand the distinct mechanisms underlying decarboxylation and carboxylation. Our calculated free energies are consistent with available experimental kinetics data. Finally, we explored how different rotamers of catalytic residues influence the efficiency of the initial intermediate formation.
Collapse
Affiliation(s)
- Kai Wen
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China
| | - Yu Tao
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China
| | - Wenyan Jiang
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China
| | - Liyan Jiang
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China
| | - Jingxuan Zhu
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China.
| | - Quanshun Li
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China; Center for Supramolecular Chemical Biology, Jilin University, Changchun 130012, China.
| |
Collapse
|
7
|
Liu X, Brooks Iii CL. Enhanced Sampling of Buried Charges in Free Energy Calculations Using Replica Exchange with Charge Tempering. J Chem Theory Comput 2024; 20:1051-1061. [PMID: 38232295 DOI: 10.1021/acs.jctc.3c00993] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
Buried ionizable groups in proteins often play important structural and functional roles. However, it is generally challenging to study the detailed molecular mechanisms solely based on experimental measurements. Free energy calculations using atomistic simulations, on the other hand, complement experimental studies and can provide high temporal and spatial resolution information that can lead to mechanistic insights. Nevertheless, it is also well recognized that sufficient sampling of such atomistic simulations can be challenging, considering that structural changes related to the buried charges may be very slow. In the present study, we describe a simple but effective enhanced sampling technique called replica exchange with charge tempering (REChgT) with a novel free energy method, multisite λ dynamics (MSλD), to study two systems containing buried charges, pKa prediction of a small molecule, orotate, in complex with the dihydroorotate dehydrogenase, and relative stability of a Glu-Lys pair buried in the hydrophobic core of two variants of Staphylococcal nuclease. Compared to the original MSλD simulations, the usage of REChgT dramatically increases sampling in both conformational and alchemical spaces, which directly translates into a significant reduction of wall time to converge the free energy calculations. This study highlights the importance of sufficient sampling toward developing improved free energy methods.
Collapse
Affiliation(s)
- Xiaorong Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Charles L Brooks Iii
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Biophysics Program, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
8
|
Sk MF, Samanta S, Poddar S, Kar P. Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey. J Comput Aided Mol Des 2024; 38:8. [PMID: 38324213 DOI: 10.1007/s10822-023-00548-8] [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: 10/05/2023] [Accepted: 12/24/2023] [Indexed: 02/08/2024]
Abstract
The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.
Collapse
Affiliation(s)
- Md Fulbabu Sk
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Simrol, MP, 453552, India
- Theoretical and Computational Biophysics Group, Beckman Institute for Advanced Science and Technology, NIH Resource for Macromolecular Modeling and Visualization, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA
| | - Sunanda Samanta
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Simrol, MP, 453552, India
| | - Sayan Poddar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Simrol, MP, 453552, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Simrol, MP, 453552, India.
| |
Collapse
|
9
|
Brooks CL, MacKerell AD, Post CB, Nilsson L. Biomolecular dynamics in the 21st century. Biochim Biophys Acta Gen Subj 2024; 1868:130534. [PMID: 38065235 PMCID: PMC10842176 DOI: 10.1016/j.bbagen.2023.130534] [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: 09/26/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 01/03/2024]
Abstract
The relevance of motions in biological macromolecules has been clear since the early structural analyses of proteins by X-ray crystallography. Computer simulations have been applied to provide a deeper understanding of the dynamics of biological macromolecules since 1976, and are now a standard tool in many labs working on the structure and function of biomolecules. In this mini-review we highlight some areas of current interest and active development for simulations, in particular all-atom molecular dynamics simulations.
Collapse
Affiliation(s)
- Charles L Brooks
- University of Michigan, Department of Chemistry, Ann Arbor, MI 48109, USA.
| | | | - Carol B Post
- Purdue University, Department of Medicinal Chemistry and Molecular Pharmacology, West Lafayette, IN 47907-2091, USA.
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition, SE-14183 Huddinge, Sweden.
| |
Collapse
|
10
|
Yu Z, Wang Z, Cui X, Cao Z, Zhang W, Sun K, Hu G. Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics. Molecules 2024; 29:645. [PMID: 38338389 PMCID: PMC10856033 DOI: 10.3390/molecules29030645] [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: 11/27/2023] [Revised: 01/01/2024] [Accepted: 01/25/2024] [Indexed: 02/12/2024] Open
Abstract
The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.
Collapse
Affiliation(s)
- Zhiping Yu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| | - Zhen Wang
- Pingyin People’s Hospital, Jinan 250400, China; (Z.W.); (X.C.)
| | - Xiuzhen Cui
- Pingyin People’s Hospital, Jinan 250400, China; (Z.W.); (X.C.)
| | - Zanxia Cao
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| | - Wanyunfei Zhang
- School of Science, Xi’an Polytechnic University, Xi’an 710048, China; (W.Z.); (K.S.)
| | - Kunxiao Sun
- School of Science, Xi’an Polytechnic University, Xi’an 710048, China; (W.Z.); (K.S.)
| | - Guodong Hu
- Shandong Key Laboratory of Biophysics, Dezhou University, Dezhou 253023, China; (Z.Y.); (Z.C.)
| |
Collapse
|
11
|
Madhu MK, Shewani K, Murarka RK. Biased Signaling in Mutated Variants of β 2-Adrenergic Receptor: Insights from Molecular Dynamics Simulations. J Chem Inf Model 2024; 64:449-469. [PMID: 38194225 DOI: 10.1021/acs.jcim.3c01481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
The molecular basis of receptor bias in G protein-coupled receptors (GPCRs) caused by mutations that preferentially activate specific intracellular transducers over others remains poorly understood. Two experimentally identified biased variants of β2-adrenergic receptors (β2AR), a prototypical GPCR, are a triple mutant (T68F, Y132A, and Y219A) and a single mutant (Y219A); the former bias the receptor toward the β-arrestin pathway by disfavoring G protein engagement, while the latter induces G protein signaling explicitly due to selection against GPCR kinases (GRKs) that phosphorylate the receptor as a prerequisite of β-arrestin binding. Though rigorous characterizations have revealed functional implications of these mutations, the atomistic origin of the observed transducer selectivity is not clear. In this study, we investigated the allosteric mechanism of receptor bias in β2AR using microseconds of all-atom Gaussian accelerated molecular dynamics (GaMD) simulations. Our observations reveal distinct rearrangements in transmembrane helices, intracellular loop 3, and critical residues R1313.50 and Y3267.53 in the conserved motifs D(E)RY and NPxxY for the mutant receptors, leading to their specific transducer interactions. Moreover, partial dissociation of G protein from the receptor core is observed in the simulations of the triple mutant in contrast to the single mutant and wild-type receptor. The reorganization of allosteric communications from the extracellular agonist BI-167107 to the intracellular receptor-transducer interfaces drives the conformational rearrangements responsible for receptor bias in the single and triple mutants. The molecular insights into receptor bias of β2AR presented here could improve the understanding of biased signaling in GPCRs, potentially opening new avenues for designing novel therapeutics with fewer side-effects and superior efficacy.
Collapse
Affiliation(s)
- Midhun K Madhu
- Department of Biological Sciences, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| | - Kunal Shewani
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| | - Rajesh K Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal, Madhya Pradesh 462066, India
| |
Collapse
|
12
|
Liu Y, Zhao W, Jiang Y, Xing S, Li W. Study on the Mechanism of Interaction between Dipeptidyl Peptidase 4 and Inhibitory Peptides Based on Gaussian Accelerated Molecular Dynamic Simulation. Int J Mol Sci 2024; 25:839. [PMID: 38255913 PMCID: PMC10815451 DOI: 10.3390/ijms25020839] [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: 11/24/2023] [Revised: 01/05/2024] [Accepted: 01/06/2024] [Indexed: 01/24/2024] Open
Abstract
Dipeptidyl peptidase 4 (DPP4) inhibitors can effectively inhibit the activity of DPP4, increasing the concentrations of glucagon-like peptide-1 (GLP-1) and glucose-dependent insulinotropic polypeptide (GIP), which allows for them to effectively contribute to the reduction of blood sugar levels. Leu-Pro-Ala-Val-Thr-Ile-Arg (LPAVTIR) and Leu-Pro-Pro-Glu-His-Asp-Trp-Arg (LPPEHDWR) were the two peptides with the strongest inhibitory activity against DPP4 selected from silkworm pupa proteins. In this study, four systems were established: Apo (ligand-free DPP4), IPI (IPI-bound DPP4), LPAVTIR (LPAVTIR-bound DPP4), LPPEHDWR (LPPEHDWR-bound DPP4), and Gaussian accelerated molecular dynamic (GaMD) simulation was conducted to investigate the mechanism of action of two inhibitory peptides binding to DPP4. Our study revealed that the LPAVTIR peptide possessed a more stable structure and exhibited a tighter binding to the Ser630 active site in DPP4, thus exhibiting a favorable competitive inhibition effect. In contrast, the LPPEHDWR peptide caused the horizontal α-helix (residues 201-215) composed of Glu205 and Glu206 residues in DPP4 to disappear. The spatial arrangement of active sites Ser630 relative to Glu205 and Glu206 was disrupted, resulting in enzyme inactivation. Moreover, the size of the substrate channel and cavity volume was significantly reduced after the binding of the inhibitory peptide to the protein, which was an important factor in the inhibition of the enzyme activity. A similar effect was also found from IPI (our positive control). By stabilizing the active site of DPP4, the IPI peptide induced the disappearance of the horizontal α-helix and a notable reduction in the active cavity volume. In conclusion, our study provided a solid theoretical foundation for the inhibitory mechanisms of IPI, LPAVTIR, and LPPEHDWR on DPP4, offering valuable insights for advancing the development of drug targets for type 2 diabetes.
Collapse
Affiliation(s)
- Yuyang Liu
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (Y.L.); (W.Z.)
| | - Wencheng Zhao
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (Y.L.); (W.Z.)
| | - Yongxin Jiang
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun 130012, China;
| | - Shu Xing
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (Y.L.); (W.Z.)
| | - Wannan Li
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Changchun 130012, China; (Y.L.); (W.Z.)
| |
Collapse
|
13
|
Pawnikar S, Magenheimer BS, Munoz EN, Haldane A, Maser RL, Miao Y. Activation of Polycystin-1 Signaling by Binding of Stalk-derived Peptide Agonists. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.06.574465. [PMID: 38260358 PMCID: PMC10802338 DOI: 10.1101/2024.01.06.574465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2024]
Abstract
Polycystin-1 (PC1) is the membrane protein product of the PKD1 gene whose mutation is responsible for 85% of the cases of autosomal dominant polycystic kidney disease (ADPKD). ADPKD is primarily characterized by the formation of renal cysts and potential kidney failure. PC1 is an atypical G protein-coupled receptor (GPCR) consisting of 11 transmembrane helices and an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal (NTF) and membrane-embedded C-terminal (CTF) fragments. Recently, signaling activation of the PC1 CTF was shown to be regulated by a stalk tethered agonist (TA), a distinct mechanism observed in the adhesion GPCR family. A novel allosteric activation pathway was elucidated for the PC1 CTF through a combination of Gaussian accelerated molecular dynamics (GaMD), mutagenesis and cellular signaling experiments. Here, we show that synthetic, soluble peptides with 7 to 21 residues derived from the stalk TA, in particular, peptides including the first 9 residues (p9), 17 residues (p17) and 21 residues (p21) exhibited the ability to re-activate signaling by a stalkless PC1 CTF mutant in cellular assays. To reveal molecular mechanisms of stalk peptide-mediated signaling activation, we have applied a novel Peptide GaMD (Pep-GaMD) algorithm to elucidate binding conformations of selected stalk peptide agonists p9, p17 and p21 to the stalkless PC1 CTF. The simulations revealed multiple specific binding regions of the stalk peptide agonists to the PC1 protein including an "intermediate" bound yet inactive state. Our Pep-GaMD simulation findings were consistent with the cellular assay experimental data. Binding of peptide agonists to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of the PC1 CTF signaling activation mechanism. Using sequence covariation analysis of PC1 homologs, we further showed that the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. Therefore, structural dynamic insights into the mechanisms of PC1 activation by stalk-derived peptide agonists have enabled an in-depth understanding of PC1 signaling. They will form a foundation for development of PC1 as a therapeutic target for the treatment of ADPKD.
Collapse
Affiliation(s)
- Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Brenda S. Magenheimer
- Departments of Biochemistry and Molecular Biology, Temple University, Philadelphia, PA 19122
- Clinical Laboratory Sciences, Temple University, Philadelphia, PA 19122
| | | | - Allan Haldane
- Dept of Physics, and Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA 19122
| | - Robin L. Maser
- Departments of Biochemistry and Molecular Biology, Temple University, Philadelphia, PA 19122
- Clinical Laboratory Sciences, Temple University, Philadelphia, PA 19122
- The Jared Grantham Kidney Institute, University of Kansas Medical Center, Kansas City, KS 66160
| | - Yinglong Miao
- Department of Pharmacology and Computational Medicine Program, University of North Carolina – Chapel Hill, Chapel Hill, NC 27599
| |
Collapse
|
14
|
Adediwura VA, Miao Y. Mechanistic Insights into Peptide Binding and Deactivation of an Adhesion G Protein-Coupled Receptor. Molecules 2023; 29:164. [PMID: 38202747 PMCID: PMC10780249 DOI: 10.3390/molecules29010164] [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: 11/27/2023] [Revised: 12/20/2023] [Accepted: 12/23/2023] [Indexed: 01/12/2024] Open
Abstract
Adhesion G protein-coupled receptors (ADGRGs) play critical roles in the reproductive, neurological, cardiovascular, and endocrine systems. In particular, ADGRG2 plays a significant role in Ewing sarcoma cell proliferation, parathyroid cell function, and male fertility. In 2022, a cryo-EM structure was reported for the active ADGRG2 bound by an optimized peptide agonist IP15 and the Gs protein. The IP15 peptide agonist was also modified to antagonists 4PH-E and 4PH-D with mutations of the 4PH residue to Glu and Asp, respectively. However, experimental structures of inactive antagonist-bound ADGRs remain to be resolved, and the activation mechanism of ADGRs such as ADGRG2 is poorly understood. Here, we applied Gaussian accelerated molecular dynamics (GaMD) simulations to probe conformational dynamics of the agonist- and antagonist-bound ADGRG2. By performing GaMD simulations, we were able to identify important low-energy conformations of ADGRG2 in the active, intermediate, and inactive states, as well as explore the binding conformations of each peptide. Moreover, our simulations revealed critical peptide-receptor residue interactions during the deactivation of ADGRG2. In conclusion, through GaMD simulations, we uncovered mechanistic insights into peptide (agonist and antagonist) binding and deactivation of the ADGRG2. These findings will potentially facilitate rational design of new peptide modulators of ADGRG2 and other ADGRs.
Collapse
Affiliation(s)
| | - Yinglong Miao
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA;
| |
Collapse
|
15
|
Aguilar-Toalá JE, Vidal-Limon A, Liceaga AM, Zambrano-Zaragoza ML, Quintanar-Guerrero D. Application of Molecular Dynamics Simulations to Determine Interactions between Canary Seed ( Phalaris canariensis L.) Bioactive Peptides and Skin-Aging Enzymes. Int J Mol Sci 2023; 24:13420. [PMID: 37686226 PMCID: PMC10487734 DOI: 10.3390/ijms241713420] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Revised: 08/22/2023] [Accepted: 08/28/2023] [Indexed: 09/10/2023] Open
Abstract
Food bioactive peptides are well recognized for their health benefits such as antimicrobial, antioxidant, and antihypertensive benefits, among others. Their drug-like behavior has led to their potential use in targeting skin-related aging factors like the inhibition of enzymes related with the skin-aging process. In this study, canary seed peptides (CSP) after simulated gastrointestinal digestion (<3 kDa) were fractioned by RP-HPLC and their enzyme-inhibition activity towards elastase and tyrosinase was evaluated in vitro. CSP inhibited elastase (IC50 = 6.2 mg/mL) and tyrosinase (IC50 = 6.1 mg/mL), while the hydrophobic fraction-VI (0.2 mg/mL) showed the highest inhibition towards elastase (93%) and tyrosinase (67%). The peptide fraction with the highest inhibition was further characterized by a multilevel in silico workflow, including physicochemical descriptor calculations, antioxidant activity predictions, and molecular dynamics-ensemble docking towards elastase and tyrosinase. To gain insights into the skin permeation process during molecular dynamics simulations, based on their docking scores, five peptides (GGWH, VPPH, EGLEPNHRVE, FLPH, and RPVNKYTPPQ) were identified to have favorable intermolecular interactions, such as hydrogen bonding of polar residues (W, H, and K) to lipid polar groups and 2-3 Å van der Waals close contact of hydrophobic aliphatic residues (P, V, and L). These interactions can play a critical role for the passive insertion of peptides into stratum corneum model skin-membranes, suggesting a promising application of CSP for skin-aging treatments.
Collapse
Affiliation(s)
- José E. Aguilar-Toalá
- Departamento de Ciencias de la Alimentación, División de Ciencias Biológicas y de la Salud, Universidad Autónoma Metropolitana, Unidad Lerma. Av. de las Garzas 10. Col. El Panteón, Lerma de Villada 52005, Estado de México, Mexico;
| | - Abraham Vidal-Limon
- Red de Estudios Moleculares Avanzados, Instituto de Ecología A.C. (INECOL), Carretera Antigua a Coatepec 351, Xalapa 91073, Veracruz, Mexico
| | - Andrea M. Liceaga
- Protein Chemistry and Bioactive Peptides Laboratory, Purdue University, 745 Agriculture Mall, West Lafayette, IN 47907, USA
| | - Maria L. Zambrano-Zaragoza
- Laboratorio de Procesos de Transformación y Tecnologías Emergentes de Alimentos-UIM, FES-Cuautitlán, Universidad Nacional Autónoma de México, Cuautitlán Izcalli 54714, Estado de México, Mexico;
| | - David Quintanar-Guerrero
- Laboratorio de Posgrado en Tecnología Farmacéutica, FES-Cuautitlán, Universidad Nacional Autónoma de México, Av. 1o de Mayo s/n, Cuautitlán Izcalli 54714, Estado de México, Mexico;
| |
Collapse
|
16
|
Vani BP, Aranganathan A, Wang D, Tiwary P. AlphaFold2-RAVE: From Sequence to Boltzmann Ranking. J Chem Theory Comput 2023; 19:4351-4354. [PMID: 37171364 PMCID: PMC10524496 DOI: 10.1021/acs.jctc.3c00290] [Citation(s) in RCA: 19] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
While AlphaFold2 is rapidly being adopted as a new standard in protein structure predictions, it is limited to single structures. This can be insufficient for the inherently dynamic world of biomolecules. In this Letter, we propose AlphaFold2-RAVE, an efficient protocol for obtaining Boltzmann-ranked ensembles from sequence. The method uses structural outputs from AlphaFold2 as initializations for artificial intelligence-augmented molecular dynamics. We release the method as an open-source code and demonstrate results on different proteins.
Collapse
Affiliation(s)
- Bodhi P. Vani
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, USA
| | - Akashnathan Aranganathan
- Biophysics Program and Institute for Physical Science and Technology, University of Maryland, College Park 20742, USA
| | - Dedi Wang
- Biophysics Program and Institute for Physical Science and Technology, University of Maryland, College Park 20742, USA
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park 20742, USA
| |
Collapse
|
17
|
Martínez-Rosas V, Hernández-Ochoa B, Morales-Luna L, Ortega-Cuellar D, González-Valdez A, Arreguin-Espinosa R, Rufino-González Y, Calderón-Jaimes E, Castillo-Rodríguez RA, Wong-Baeza C, Baeza-Ramírez I, Pérez de la Cruz V, Vidal-Limón A, Gómez-Manzo S. Nitazoxanide Inhibits the Bifunctional Enzyme GlG6PD::6PGL of Giardia lamblia: Biochemical and In Silico Characterization of a New Druggable Target. Int J Mol Sci 2023; 24:11516. [PMID: 37511272 PMCID: PMC10380810 DOI: 10.3390/ijms241411516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Revised: 07/08/2023] [Accepted: 07/12/2023] [Indexed: 07/30/2023] Open
Abstract
Giardiasis, which is caused by Giardia lamblia infection, is a relevant cause of morbidity and mortality worldwide. Because no vaccines are currently available to treat giardiasis, chemotherapeutic drugs are the main options for controlling infection. Evidence has shown that the nitro drug nitazoxanide (NTZ) is a commonly prescribed treatment for giardiasis; however, the mechanisms underlying NTZ's antigiardial activity are not well-understood. Herein, we identified the glucose-6-phosphate::6-phosphogluconate dehydrogenase (GlG6PD::6PGL) fused enzyme as a nitazoxanide target, as NTZ behaves as a GlG6PD::6PGL catalytic inhibitor. Furthermore, fluorescence assays suggest alterations in the stability of GlG6PD::6PGL protein, whereas the results indicate a loss of catalytic activity due to conformational and folding changes. Molecular docking and dynamic simulation studies suggest a model of NTZ binding on the active site of the G6PD domain and near the structural NADP+ binding site. The findings of this study provide a novel mechanistic basis and strategy for the antigiardial activity of the NTZ drug.
Collapse
Affiliation(s)
- Víctor Martínez-Rosas
- Laboratorio de Bioquímica Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City 04530, Mexico
- Programa de Posgrado en Biomedicina y Biotecnología Molecular, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Mexico City 11340, Mexico
| | - Beatriz Hernández-Ochoa
- Laboratorio de Inmunoquímica, Hospital Infantil de México Federico Gómez, Secretaría de Salud, Mexico City 06720, Mexico
| | - Laura Morales-Luna
- Laboratorio de Bioquímica Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City 04530, Mexico
- Posgrado en Ciencias Biológicas, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
| | - Daniel Ortega-Cuellar
- Laboratorio de Nutrición Experimental, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City 04530, Mexico
| | - Abigail González-Valdez
- Departamento de Biología Molecular y Biotecnología, Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
| | - Roberto Arreguin-Espinosa
- Departamento de Química de Biomacromoléculas, Instituto de Química, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
| | - Yadira Rufino-González
- Laboratorio de Parasitología Experimental, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City 04530, Mexico
| | - Ernesto Calderón-Jaimes
- Laboratorio de Inmunoquímica, Hospital Infantil de México Federico Gómez, Secretaría de Salud, Mexico City 06720, Mexico
| | | | - Carlos Wong-Baeza
- Laboratorio de Biomembranas, Departamento de Bioquímica, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Mexico City 11340, Mexico
| | - Isabel Baeza-Ramírez
- Laboratorio de Biomembranas, Departamento de Bioquímica, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Mexico City 11340, Mexico
| | - Verónica Pérez de la Cruz
- Neurobiochemistry and Behavior Laboratory, National Institute of Neurology and Neurosurgery "Manuel Velasco Suárez", Mexico City 14269, Mexico
| | - Abraham Vidal-Limón
- Red de Estudios Moleculares Avanzados, Clúster Científico y Tecnológico BioMimic®, Instituto de Ecología A.C. (INECOL), Carretera Antigua a Coatepec 351, El Haya, Xalapa 91073, Mexico
| | - Saúl Gómez-Manzo
- Laboratorio de Bioquímica Genética, Instituto Nacional de Pediatría, Secretaría de Salud, Mexico City 04530, Mexico
| |
Collapse
|
18
|
Castillo-Campos L, Velázquez-Libera JL, Caballero J. Computational study of the binding orientation and affinity of noncovalent inhibitors of the papain-like protease (PLpro) from SARS-CoV-1 considering the protein flexibility by using molecular dynamics and cross-docking. Front Mol Biosci 2023; 10:1215499. [PMID: 37426421 PMCID: PMC10326900 DOI: 10.3389/fmolb.2023.1215499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Accepted: 06/12/2023] [Indexed: 07/11/2023] Open
Abstract
The papain-like protease (PLpro) from zoonotic coronaviruses (CoVs) has been identified as a target with an essential role in viral respiratory diseases caused by Severe Acute Respiratory Syndrome-associated coronaviruses (SARS-CoVs). The design of PLpro inhibitors has been proposed as an alternative to developing potential drugs against this disease. In this work, 67 naphthalene-derived compounds as noncovalent PLpro inhibitors were studied using molecular modeling methods. Structural characteristics of the bioactive conformations of these inhibitors and their interactions at the SARS-CoV-1 PLpro binding site were reported here in detail, taking into account the flexibility of the protein residues. Firstly, a molecular docking protocol was used to obtain the orientations of the inhibitors. After this, the orientations were compared, and the recurrent interactions between the PLpro residues and ligand chemical groups were described (with LigRMSD and interaction fingerprints methods). In addition, efforts were made to find correlations between docking energy values and experimentally determined binding affinities. For this, the PLpro was sampled by using Gaussian Accelerated Molecular Dynamics (GaMD), generating multiple conformations of the binding site. Diverse protein conformations were selected and a cross-docking experiment was performed, yielding models of the 67 naphthalene-derived compounds adopting different binding modes. Representative complexes for each ligand were selected to obtain the highest correlation between docking energies and activities. A good correlation (R 2 = 0.948) was found when this flexible docking protocol was performed.
Collapse
|
19
|
Yang L, Liu T, Zhuo Y, Li D, Li D, Liu J, Gao H, Zhang L, Lin J, Wang X. Verbenalin alleviates acute lung injury induced by sepsis and IgG immune complex through GPR18 receptor. Cell Signal 2023:110768. [PMID: 37315751 DOI: 10.1016/j.cellsig.2023.110768] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Revised: 04/24/2023] [Accepted: 06/10/2023] [Indexed: 06/16/2023]
Abstract
Acute lung injury is significantly associated with the aberrant activation and pyroptosis of alveolar macrophages. Targeting the GPR18 receptor presents a potential therapeutic approach to mitigate inflammation. Verbenalin, a prominent component of Verbena in Xuanfeibaidu (XFBD) granules, is recommended for treating COVID-19. In this study, we demonstrate the therapeutic effect of verbenalin on lung injury through direct binding to the GPR18 receptor. Verbenalin inhibits the activation of inflammatory signaling pathways induced by lipopolysaccharide (LPS) and IgG immune complex (IgG IC) via GPR18 receptor activation. The structural basis for verbenalin's effect on GPR18 activation is elucidated through molecular docking and molecular dynamics simulations. Furthermore, we establish that IgG IC induces macrophage pyroptosis by upregulating the expression of GSDME and GSDMD through CEBP-δ activation, while verbenalin inhibits this process. Additionally, we provide the first evidence that IgG IC promotes the formation of neutrophil extracellular traps (NETs), and verbenalin suppresses NETs formation. Collectively, our findings indicate that verbenalin functions as a "phytoresolvin" to promote inflammation regression and suggests that targeting the C/EBP-δ/GSDMD/GSDME axis to inhibit macrophage pyroptosis may represent a novel strategy for treating acute lung injury and sepsis.
Collapse
Affiliation(s)
- Lei Yang
- Tianjin Key Laboratory of Acute Abdomen Disease Associated Organ Injury and ITCWM Repair, Tianjin Nankai Hospital, Tianjin, China
| | - Tianyu Liu
- Tianjin Key Laboratory of Acute Abdomen Disease Associated Organ Injury and ITCWM Repair, Tianjin Nankai Hospital, Tianjin, China; Graduate School, Tianjin Medical University, Tianjin, China
| | - Yuzhen Zhuo
- Tianjin Key Laboratory of Acute Abdomen Disease Associated Organ Injury and ITCWM Repair, Tianjin Nankai Hospital, Tianjin, China
| | - Dongmei Li
- College of Pharmacy, Nankai University, Tianjin, China
| | - Dihua Li
- Tianjin Key Laboratory of Acute Abdomen Disease Associated Organ Injury and ITCWM Repair, Tianjin Nankai Hospital, Tianjin, China
| | - Junhong Liu
- Tianjin Key Laboratory of Acute Abdomen Disease Associated Organ Injury and ITCWM Repair, Tianjin Nankai Hospital, Tianjin, China
| | - Hejun Gao
- Graduate School, Tianjin Medical University, Tianjin, China
| | - Lanqiu Zhang
- Tianjin Key Laboratory of Acute Abdomen Disease Associated Organ Injury and ITCWM Repair, Tianjin Nankai Hospital, Tianjin, China.
| | - Jianping Lin
- College of Pharmacy, Nankai University, Tianjin, China.
| | - Ximo Wang
- Tianjin Key Laboratory of Acute Abdomen Disease Associated Organ Injury and ITCWM Repair, Tianjin Nankai Hospital, Tianjin, China; Graduate School, Tianjin Medical University, Tianjin, China; Tianjin Key Laboratory of Acute Abdomen Disease Associated Organ Injury and ITCWM Repair, Institute of Integrative Medicine for Acute Abdominal Diseases, Tianjin University, Tianjin, China.
| |
Collapse
|
20
|
Li Y, Acharya A, Yang L, Liu J, Tajkhorshid E, Zgurskaya HI, Jackson M, Gumbart JC. Insights into substrate transport and water permeation in the mycobacterial transporter MmpL3. Biophys J 2023; 122:2342-2352. [PMID: 36926696 PMCID: PMC10257117 DOI: 10.1016/j.bpj.2023.03.018] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Revised: 02/04/2023] [Accepted: 03/10/2023] [Indexed: 03/17/2023] Open
Abstract
Mycobacteria, such as Mycobacterium tuberculosis, are characterized by a uniquely thick and waxy cell envelope that consists of two membranes, with a variety of mycolates comprising their outer membrane (OM). The protein Mycobacterial membrane protein Large 3 (MmpL3) is responsible for the transport of a primary OM component, trehalose monomycolate (TMM), from the inner (cytoplasmic) membrane (IM) to the periplasmic space, a process driven by the proton gradient. Although multiple structures of MmpL3 with bound substrates have been solved, the exact pathway(s) for TMM or proton transport remains elusive. Here, employing molecular dynamics simulations we investigate putative pathways for either transport species. We hypothesized that MmpL3 will cycle through similar conformational states as the related transporter AcrB, which we used as targets for modeling the conformation of MmpL3. A continuous water pathway through the transmembrane region was found in one of these states, illustrating a putative pathway for protons. Additional equilibrium simulations revealed that TMM can diffuse from the membrane into a binding pocket in MmpL3 spontaneously. We also found that acetylation of TMM, which is required for transport, makes it more stable within MmpL3's periplasmic cavity compared with the unacetylated form.
Collapse
Affiliation(s)
- Yupeng Li
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois Urbana-Champaign, Urbana, Illinois; Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois
| | - Atanu Acharya
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia
| | - Lixinhao Yang
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia
| | - Jinchan Liu
- Department of Molecular Biophysics and Biochemistry (MB&B), Yale University, New Haven, Connecticut
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois Urbana-Champaign, Urbana, Illinois; Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois; Department of Biochemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois
| | - Helen I Zgurskaya
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma
| | - Mary Jackson
- Mycobacteria Research Laboratories, Department of Microbiology, Immunology and Pathology, Colorado State University, Fort Collins, Colorado
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia; School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia.
| |
Collapse
|
21
|
Vuckovic Z, Wang J, Pham V, Mobbs JI, Belousoff MJ, Bhattarai A, Burger WAC, Thompson G, Yeasmin M, Nawaratne V, Leach K, van der Westhuizen ET, Khajehali E, Liang YL, Glukhova A, Wootten D, Lindsley CW, Tobin A, Sexton P, Danev R, Valant C, Miao Y, Christopoulos A, Thal DM. Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics. eLife 2023; 12:83477. [PMID: 37248726 DOI: 10.7554/elife.83477] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Accepted: 04/12/2023] [Indexed: 05/31/2023] Open
Abstract
Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M4 muscarinic acetylcholine receptor (M4 mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M4 mAChR. By determining the cryo-electron microscopy structures of the M4 mAChR, bound to a cognate Gi1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand-receptor-transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.
Collapse
Affiliation(s)
- Ziva Vuckovic
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
| | - Vi Pham
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Jesse I Mobbs
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Matthew J Belousoff
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
| | - Wessel A C Burger
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Geoff Thompson
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Mahmuda Yeasmin
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Vindhya Nawaratne
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Katie Leach
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Emma T van der Westhuizen
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Elham Khajehali
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Yi-Lynn Liang
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Alisa Glukhova
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Denise Wootten
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Craig W Lindsley
- Department of Pharmacology, Warren Center for Neuroscience Drug Discovery and Department of Chemistry, Warren Center for Neuroscience Drug Discovery, Vanderbilt University, Nashville, United States
| | - Andrew Tobin
- The Centre for Translational Pharmacology, Advanced Research Centre (ARC), College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, United Kingdom
| | - Patrick Sexton
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Radostin Danev
- Graduate School of Medicine, University of Tokyo, Tokyo, Japan
| | - Celine Valant
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, United States
| | - Arthur Christopoulos
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- Neuromedicines Discovery Centre, Monash University, Parkville, Australia
| | - David M Thal
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
- ARC Centre for Cryo-electron Microscopy of Membrane Proteins, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, Australia
| |
Collapse
|
22
|
Rodríguez-Hernández A, Bonilla-Landa I, Vidal-Limon A, Ibarra-Juárez A, Barrera-Méndez F, Medrano FJE, Díaz de León RE, Olivares-Romero JL. Synthesis, insecticidal activity, and ensembled docking of nitroguanidines bearing S- and R-proline. PEST MANAGEMENT SCIENCE 2023; 79:1912-1921. [PMID: 36656155 DOI: 10.1002/ps.7368] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 12/19/2022] [Accepted: 01/19/2023] [Indexed: 06/17/2023]
Abstract
BACKGROUND The amino acids R- and S-proline were used to synthesize novel neonicotinoid derivatives that, after being characterized by 1 H, DEPTQ 135, and HRMS-QTOF, were evaluated for use as insecticides against Galleria mellonella (caterpillar), Sitophilus zeamais, Xylosandrus morigerus, Xyleborus affinis, and Xyleborus ferrugineus. RESULTS Comparisons of biological activity and absolute configuration showed that the R enantiomer had excellent and outstanding insecticidal activity against the insects tested, with up to 100% mortality after 12 h compared with dinotefuran at the same concentration. CONCLUSIONS The results suggest that compound R6 is an excellent lead enantiopure insecticide for future development in the field of crop protection. Furthermore, intermolecular interactions between nicotinic acetylcholine receptors and the R enantiomer displays a lower score which mean a higher affinity to the nAChR receptor and the π-π interactions are more stable than the S derivative. © 2023 Society of Chemical Industry.
Collapse
Affiliation(s)
- Alfredo Rodríguez-Hernández
- Instituto de Ecología, A.C, Red de Estudios Moleculares Avanzados, Clúster Científico y Tecnológico BioMimic, Veracruz, Mexico
| | - Israel Bonilla-Landa
- Instituto de Ecología, A.C, Red de Estudios Moleculares Avanzados, Clúster Científico y Tecnológico BioMimic, Veracruz, Mexico
| | - Abraham Vidal-Limon
- Instituto de Ecología, A.C, Red de Estudios Moleculares Avanzados, Clúster Científico y Tecnológico BioMimic, Veracruz, Mexico
| | - Arturo Ibarra-Juárez
- Instituto de Ecología, A.C, Red de Estudios Moleculares Avanzados, Clúster Científico y Tecnológico BioMimic, Veracruz, Mexico
| | - Felipe Barrera-Méndez
- Instituto de Ecología, A.C, Red de Estudios Moleculares Avanzados, Clúster Científico y Tecnológico BioMimic, Veracruz, Mexico
- Cátedra CONACyT en el Instituto de Ecología, A. C, Veracruz, Mexico
| | | | | | - José Luis Olivares-Romero
- Instituto de Ecología, A.C, Red de Estudios Moleculares Avanzados, Clúster Científico y Tecnológico BioMimic, Veracruz, Mexico
| |
Collapse
|
23
|
Ojha AA, Srivastava A, Votapka LW, Amaro RE. Selectivity and Ranking of Tight-Binding JAK-STAT Inhibitors Using Markovian Milestoning with Voronoi Tessellations. J Chem Inf Model 2023; 63:2469-2482. [PMID: 37023323 PMCID: PMC10131228 DOI: 10.1021/acs.jcim.2c01589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
Abstract
Janus kinases (JAK), a group of proteins in the nonreceptor tyrosine kinase (NRTKs) family, play a crucial role in growth, survival, and angiogenesis. They are activated by cytokines through the Janus kinase-signal transducer and activator of a transcription (JAK-STAT) signaling pathway. JAK-STAT signaling pathways have significant roles in the regulation of cell division, apoptosis, and immunity. Identification of the V617F mutation in the Janus homology 2 (JH2) domain of JAK2 leading to myeloproliferative disorders has stimulated great interest in the drug discovery community to develop JAK2-specific inhibitors. However, such inhibitors should be selective toward JAK2 over other JAKs and display an extended residence time. Recently, novel JAK2/STAT5 axis inhibitors (N-(1H-pyrazol-3-yl)pyrimidin-2-amino derivatives) have displayed extended residence times (hours or longer) on target and adequate selectivity excluding JAK3. To facilitate a deeper understanding of the kinase-inhibitor interactions and advance the development of such inhibitors, we utilize a multiscale Markovian milestoning with Voronoi tessellations (MMVT) approach within the Simulation-Enabled Estimation of Kinetic Rates v.2 (SEEKR2) program to rank order these inhibitors based on their kinetic properties and further explain the selectivity of JAK2 inhibitors over JAK3. Our approach investigates the kinetic and thermodynamic properties of JAK-inhibitor complexes in a user-friendly, fast, efficient, and accurate manner compared to other brute force and hybrid-enhanced sampling approaches.
Collapse
Affiliation(s)
- Anupam Anand Ojha
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Ambuj Srivastava
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Lane William Votapka
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
24
|
Hasse T, Zhang Z, Huang YMM. Molecular dynamics study reveals key disruptors of MEIG1-PACRG interaction. Proteins 2023; 91:555-566. [PMID: 36444670 PMCID: PMC10374433 DOI: 10.1002/prot.26449] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Revised: 11/07/2022] [Accepted: 11/15/2022] [Indexed: 11/30/2022]
Abstract
Interactions between the meiosis-expressed gene 1 (MEIG1) and Parkin co-regulated gene (PACRG) protein are critical in the formation of mature sperm cells. Targeting either MEIG1 or PACRG protein could be a contraceptive strategy. The W50A and Y68A mutations on MEIG1 are known to interrupt the MEIG1-PACRG interactions resulting in defective sperm cells. However, the details about how the mutants disrupt the protein-protein binding are not clear. In this study, we reveal insights on MEIG1 and PACRG protein dynamics by applying Gaussian-accelerated molecular dynamics (GaMD) simulations and post-GaMD analysis. Our results show that the mutations destabilize the protein-protein interfacial interaction. The effect of the Y68A mutation is more significant than W50A as Y68 forms stronger polar interactions with PACRG. Because both human and mouse models demonstrate similar dynamic properties, the findings from mouse proteins can be applied to the human system. Moreover, we report a potential ligand binding pocket on the MEIG1 and PACRG interaction surface that could be a target for future drug design to inhibit the MEIG1-PACRG interaction. PACRG shows more qualified pockets along the protein-protein interface, implying that it is a better target than MEIG1. Our work provides a fundamental understanding of MEIG1 and PACRG protein dynamics, paving the way for drug discovery in male-based contraception.
Collapse
Affiliation(s)
- Timothy Hasse
- Department of Physics and Astronomy, Wayne State University, Detroit, Michigan, USA
| | - Zhibing Zhang
- Department of Physiology, Wayne State University School of Medicine, Detroit, Michigan, USA.,Department of Obstetrics and Gynecology, Wayne State University, Detroit, Michigan, USA
| | - Yu-Ming M Huang
- Department of Physics and Astronomy, Wayne State University, Detroit, Michigan, USA
| |
Collapse
|
25
|
Ojha AA, Thakur S, Ahn SH, Amaro RE. DeepWEST: Deep Learning of Kinetic Models with the Weighted Ensemble Simulation Toolkit for Enhanced Sampling. J Chem Theory Comput 2023; 19:1342-1359. [PMID: 36719802 DOI: 10.1021/acs.jctc.2c00282] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Recent advances in computational power and algorithms have enabled molecular dynamics (MD) simulations to reach greater time scales. However, for observing conformational transitions associated with biomolecular processes, MD simulations still have limitations. Several enhanced sampling techniques seek to address this challenge, including the weighted ensemble (WE) method, which samples transitions between metastable states using many weighted trajectories to estimate kinetic rate constants. However, initial sampling of the potential energy surface has a significant impact on the performance of WE, i.e., convergence and efficiency. We therefore introduce deep-learned kinetic modeling approaches that extract statistically relevant information from short MD trajectories to provide a well-sampled initial state distribution for WE simulations. This hybrid approach overcomes any statistical bias to the system, as it runs short unbiased MD trajectories and identifies meaningful metastable states of the system. It is shown to provide a more refined free energy landscape closer to the steady state that could efficiently sample kinetic properties such as rate constants.
Collapse
Affiliation(s)
- Anupam Anand Ojha
- Department of Chemistry, University of California San Diego, La Jolla, California92093, United States
| | - Saumya Thakur
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, Maharashtra400076, India
| | - Surl-Hee Ahn
- Department of Chemical Engineering, University of California Davis, Davis, California95616, United States
| | - Rommie E Amaro
- Department of Chemistry, University of California San Diego, La Jolla, California92093, United States
| |
Collapse
|
26
|
Sayed AM, Ibrahim AH, Tajuddeen N, Seibel J, Bodem J, Geiger N, Striffler K, Bringmann G, Abdelmohsen UR. Korupensamine A, but not its atropisomer, korupensamine B, inhibits SARS-CoV-2 in vitro by targeting its main protease (M pro). Eur J Med Chem 2023; 251:115226. [PMID: 36893625 PMCID: PMC9972725 DOI: 10.1016/j.ejmech.2023.115226] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2022] [Revised: 02/07/2023] [Accepted: 02/20/2023] [Indexed: 03/06/2023]
Abstract
By combining docking and molecular dynamics simulations, we explored a library of 65 mostly axially chiral naphthylisoquinoline alkaloids and their analogues, with most different molecular architectures and structural analogues, for their activity against SARS-CoV-2. Although natural biaryls are often regarded without consideration of their axial chirality, they can bind to protein targets in an atroposelective manner. By combining docking results with steered molecular dynamics simulations, we identified one alkaloid, korupensamine A, that atropisomer-specifically inhibited the main protease (Mpro) activity of SARS-CoV-2 significantly in comparison to the reference covalent inhibitor GC376 (IC50 = 2.52 ± 0.14 and 0.88 ± 0.15 μM, respectively) and reduced viral growth by five orders of magnitude in vitro (EC50 = 4.23 ± 1.31 μM). To investigate the binding pathway and mode of interaction of korupensamine A within the active site of the protease, we utilized Gaussian accelerated molecular dynamics simulations, which reproduced the docking pose of korupensamine A inside the active site of the enzyme. The study presents naphthylisoquinoline alkaloids as a new class of potential anti-COVID-19 agents.
Collapse
Affiliation(s)
- Ahmed M Sayed
- Department of Pharmacognosy, Faculty of Pharmacy, Nahda University, Beni-Suef, 62513, Egypt
| | - Alyaa Hatem Ibrahim
- Department of Pharmacognosy, Faculty of Pharmacy, Sohag University, Sohag, 82524, Egypt
| | - Nasir Tajuddeen
- Department of Chemistry, Ahmadu Bello University, 15 Sokoto Road Samaru, Zaria, 810107, Nigeria
| | - Jürgen Seibel
- Institute of Organic Chemistry, University of Würzburg, Am Hubland, 97074, Würzburg, Germany
| | - Jochen Bodem
- Institute of Virology and Immunobiology, University of Würzburg, Versbacher Str. 7, 97078, Würzburg, Germany
| | - Nina Geiger
- Institute of Virology and Immunobiology, University of Würzburg, Versbacher Str. 7, 97078, Würzburg, Germany
| | - Kathrin Striffler
- Institute of Virology and Immunobiology, University of Würzburg, Versbacher Str. 7, 97078, Würzburg, Germany
| | - Gerhard Bringmann
- Institute of Organic Chemistry, University of Würzburg, Am Hubland, 97074, Würzburg, Germany.
| | - Usama Ramadan Abdelmohsen
- Department of Pharmacognosy, Faculty of Pharmacy, Minia University, Minia, 61519, Egypt; Department of Pharmacognosy, Faculty of Pharmacy, Deraya University, Universities Zone, New Minia City, 61111, Egypt.
| |
Collapse
|
27
|
Huerta-Miranda GA, García-García WI, Vidal-Limon A, Miranda-Hernández M. Use of simplified models for theoretical prediction of the interactions between available antibodies and the receptor-binding domain of SARS-CoV-2 spike protein. J Biomol Struct Dyn 2023; 41:1018-1027. [PMID: 34935602 DOI: 10.1080/07391102.2021.2019123] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The negative impact of infectious diseases like COVID-19 on public health and the global economy is evident. This pandemic represents a significant challenge for the scientific community to develop new practical analytical methods for accurately diagnosing emerging cases. Due to their selectivity and sensitivity, new methodologies based on antigen/antibody interactions to detect COVID-19 biomarkers are necessary. In this context, the theoretical, computational modeling reduces experimental efforts and saves resources for rational biosensor design. This study proposes using molecular dynamics to predict the interactions between the Receptor Binding Domain (RBD) of the SARS-CoV-2 spike protein simplified model and a set of highly characterized antibodies. The binding free energy of the antigen/antibody complexes was calculated for the simplified models and compared against the complete SARS-CoV-2 ectodomain to validate the methodology. The structural data derived from our molecular dynamics and end-point free energy calculations showed a positive correlation between both approximations, with a 0.82 Pearson correlation coefficient; t = 3.661, df = 3, p-value = 0.03522, with a 95% confident interval. Furthermore, we identified the interfacial residues that could generate covalent bonds with a specific chemical surface without perturbing the binding dynamics to develop highly sensitive and specific diagnostic devices. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- G A Huerta-Miranda
- Instituto de Energías Renovables, Universidad Nacional Autónoma de México, Temixco, Morelos, México
| | - W I García-García
- Instituto de Energías Renovables, Universidad Nacional Autónoma de México, Temixco, Morelos, México
| | - A Vidal-Limon
- Instituto de Energías Renovables, Universidad Nacional Autónoma de México, Temixco, Morelos, México
| | - M Miranda-Hernández
- Instituto de Energías Renovables, Universidad Nacional Autónoma de México, Temixco, Morelos, México
| |
Collapse
|
28
|
Williams RV, Rogals MJ, Eletsky A, Huang C, Morris LC, Moremen KW, Prestegard JH. AssignSLP_GUI, a software tool exploiting AI for NMR resonance assignment of sparsely labeled proteins. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2022; 345:107336. [PMID: 36442299 PMCID: PMC9742323 DOI: 10.1016/j.jmr.2022.107336] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Revised: 11/12/2022] [Accepted: 11/15/2022] [Indexed: 05/06/2023]
Abstract
Not all proteins are amenable to uniform isotopic labeling with 13C and 15N, something needed for the widely used, and largely deductive, triple resonance assignment process. Among them are proteins expressed in mammalian cell culture where native glycosylation can be maintained, and proper formation of disulfide bonds facilitated. Uniform labeling in mammalian cells is prohibitively expensive, but sparse labeling with one or a few isotopically enriched amino acid types is an option for these proteins. However, assignment then relies on accessing the best match between a variety of measured NMR parameters and predictions based on 3D structure, often from X-ray crystallography. Finding this match is a challenging process that has benefitted from many computational tools, including trained neural nets for chemical shift prediction, genetic algorithms for searches through a myriad of assignment possibilities, and now AI-based prediction of high-quality structures for protein targets. AssignSLP_GUI, a new version of a software package for assignment of resonances from sparsely-labeled proteins, uses many of these tools. These tools and new additions to the package are highlighted in an application to a sparsely-labeled domain from a glycoprotein, CEACAM1.
Collapse
Affiliation(s)
- Robert V Williams
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602, USA
| | - Monique J Rogals
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602, USA
| | - Alexander Eletsky
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602, USA
| | - Chin Huang
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602, USA
| | - Laura C Morris
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602, USA
| | - Kelley W Moremen
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602, USA
| | - James H Prestegard
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602, USA.
| |
Collapse
|
29
|
Iida S, Tomoshi K. Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules. Biophys Rev 2022; 14:1303-1314. [PMID: 36659997 PMCID: PMC9842846 DOI: 10.1007/s12551-022-01036-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Accepted: 11/26/2022] [Indexed: 12/30/2022] Open
Abstract
Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time-mean-first passage time-from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation.
Collapse
Affiliation(s)
- Shinji Iida
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-Ku, Tokyo, 135-0064 Japan
| | - Kameda Tomoshi
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-Ku, Tokyo, 135-0064 Japan
| |
Collapse
|
30
|
Chen X, Yuan Y, Chen Y, Yu J, Wang J, Chen J, Guo Y, Pu X. Biased Activation Mechanism Induced by GPCR Heterodimerization: Observations from μOR/δOR Dimers. J Chem Inf Model 2022; 62:5581-5600. [PMID: 36377848 DOI: 10.1021/acs.jcim.2c00962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
GPCRs regulate multiple intracellular signaling cascades. Biasedly activating one signaling pathway over the others provides additional clinical utility to optimize GPCR-based therapies. GPCR heterodimers possess different functions from their monomeric states, including their selectivity to different transducers. However, the biased signaling mechanism induced by the heterodimerization remains unclear. Motivated by the issue, we select an important GPCR heterodimer (μOR/δOR heterodimer) as a case and use microsecond Gaussian accelerated molecular dynamics simulation coupled with potential of mean force and protein structure network (PSN) to probe mechanisms regarding the heterodimerization-induced constitutive β-arrestin activity and efficacy change of the agonist DAMGO. The results show that only the lowest energy state of the μOR/δOR heterodimer, which adopts a slightly outward shift of TM6 and an ICL2 conformation close to the receptor core, can selectively accommodate β-arrestins. PSN further reveals important roles of H8, ICL1, and ICL2 in regulating the constitutive β-arrestin-biased activity for the apo μOR/δOR heterodimer. In addition, the heterodimerization can allosterically alter the binding mode of DAMGO mainly by means of W7.35. Consequently, DAMGO transmits the structural signal mainly through TM6 and TM7 in the dimer, rather than TM3 similar to the μOR monomer, thus changing the efficacy of DAMGO from a balanced agonist to the β-arrestin-biased one. On the other side, the binding of DAMGO to the heterodimer can stabilize μOR/δOR heterodimers through a stronger interaction of TM1/TM1 and H8/H8, accordingly enhancing the interaction of μOR with δOR and the binding affinity of the dimer to the β-arrestin. The agonist DAMGO does not change main compositions of the regulation network from the dimer interface to the transducer binding pocket of the μOR protomer, but induces an increase in the structural communication of the network, which should contribute to the enhanced β-arrestin coupling. Our observations, for the first time, reveal the molecular mechanism of the biased signaling induced by the heterodimerization for GPCRs, which should be beneficial to more comprehensively understand the GPCR bias signaling.
Collapse
Affiliation(s)
- Xin Chen
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Yuan Yuan
- College of Management, Southwest University for Nationalities, Chengdu610041, China
| | - Yichi Chen
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Jin Yu
- Department of Physics and Astronomy, University of California, Irvine, California92697, United States
| | - Jingzhou Wang
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Jianfang Chen
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Yanzhi Guo
- College of Chemistry, Sichuan University, Chengdu610064, China
| | - Xuemei Pu
- College of Chemistry, Sichuan University, Chengdu610064, China
| |
Collapse
|
31
|
Wu F, Jin S, Jiang Y, Jin X, Tang B, Niu Z, Liu X, Zhang Q, Zeng X, Li SZ. Pre-Training of Equivariant Graph Matching Networks with Conformation Flexibility for Drug Binding. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2022; 9:e2203796. [PMID: 36202759 PMCID: PMC9685463 DOI: 10.1002/advs.202203796] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Revised: 09/07/2022] [Indexed: 05/16/2023]
Abstract
The latest biological findings observe that the motionless "lock-and-key" theory is not generally applicable and that changes in atomic sites and binding pose can provide important information for understanding drug binding. However, the computational expenditure limits the growth of protein trajectory-related studies, thus hindering the possibility of supervised learning. A spatial-temporal pre-training method based on the modified equivariant graph matching networks, dubbed ProtMD which has two specially designed self-supervised learning tasks: atom-level prompt-based denoising generative task and conformation-level snapshot ordering task to seize the flexibility information inside molecular dynamics (MD) trajectories with very fine temporal resolutions is presented. The ProtMD can grant the encoder network the capacity to capture the time-dependent geometric mobility of conformations along MD trajectories. Two downstream tasks are chosen to verify the effectiveness of ProtMD through linear detection and task-specific fine-tuning. A huge improvement from current state-of-the-art methods, with a decrease of 4.3% in root mean square error for the binding affinity problem and an average increase of 13.8% in the area under receiver operating characteristic curve and the area under the precision-recall curve for the ligand efficacy problem is observed. The results demonstrate a strong correlation between the magnitude of conformation's motion in the 3D space and the strength with which the ligand binds with its receptor.
Collapse
Affiliation(s)
- Fang Wu
- School of EngineeringWestlake UniversityHangzhou310024China
- MindRank AI Ltd.Hangzhou310000China
| | - Shuting Jin
- MindRank AI Ltd.Hangzhou310000China
- School of InformaticsXiamen UniversityXiamen361005China
| | | | | | | | | | - Xiangrong Liu
- School of InformaticsXiamen UniversityXiamen361005China
| | - Qiang Zhang
- ZJU‐Hangzhou Global Scientific and Technological Innovation CenterHangzhou311200China
- College of Computer Science and TechnologyZhejiang UniversityHangzhou310013China
| | - Xiangxiang Zeng
- School of Information Science and EngineeringHunan UniversityHunan410082China
| | - Stan Z. Li
- School of EngineeringWestlake UniversityHangzhou310024China
| |
Collapse
|
32
|
Hofer F, Fischer AL, Kamenik AS, Waibl F, Fernández-Quintero ML, Liedl KR. pH-dependent structural diversity of profilin allergens determines thermal stability. FRONTIERS IN ALLERGY 2022; 3:1007000. [PMID: 36324331 PMCID: PMC9618696 DOI: 10.3389/falgy.2022.1007000] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Accepted: 09/20/2022] [Indexed: 11/05/2022] Open
Abstract
The family of profilin allergens is a common class of proteins found in plants, viruses and various eukaryotes including mammals. Profilins are characterized by an evolutionary conserved structural fold, which is responsible for their cross-reactive nature of Immunoglobulin E (IgE) antibodies. Despite their high overall structural similarity, they exhibit substantial differences in their biophysical properties, such as thermal and pH stability. To understand the origin of these functional differences of Amb a 8, Art v 4 and Bet v 2, we performed constant pH molecular dynamics simulation in combination with Gaussian accelerated MD simulations. Depending on the respective protonation at different pH levels, we find distinct differences in conformational flexibility, which are consistent with experimentally determined melting temperatures. These variations in flexibility are accompanied by ensemble shifts in the conformational landscape and quantified and localized by residue-wise B-factors and dihedral entropies. These findings strengthen the link between flexibility of profilin allergens and their thermal stability. Thus, our results clearly show the importance of considering protonation dependent conformational ensembles in solution to elucidate biophysical differences between these structurally similar allergens.
Collapse
|
33
|
Geethakumari AM, Ahmed WS, Rasool S, Fatima A, Nasir Uddin SM, Aouida M, Biswas KH. A genetically encoded BRET-based SARS-CoV-2 M pro protease activity sensor. Commun Chem 2022; 5:117. [PMID: 36187754 PMCID: PMC9516532 DOI: 10.1038/s42004-022-00731-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2022] [Accepted: 09/07/2022] [Indexed: 01/27/2023] Open
Abstract
The main protease, Mpro, is critical for SARS-CoV-2 replication and an appealing target for designing anti-SARS-CoV-2 agents. Therefore, there is a demand for the development of improved sensors to monitor its activity. Here, we report a pair of genetically encoded, bioluminescence resonance energy transfer (BRET)-based sensors for detecting Mpro proteolytic activity in live cells as well as in vitro. The sensors were generated by sandwiching peptides containing the Mpro N-terminal autocleavage sites, either AVLQSGFR (short) or KTSAVLQSGFRKME (long), in between the mNeonGreen and NanoLuc proteins. Co-expression of the sensors with Mpro in live cells resulted in their cleavage while mutation of the critical C145 residue (C145A) in Mpro completely abrogated their cleavage. Additionally, the sensors recapitulated the inhibition of Mpro by the well-characterized pharmacological agent GC376. Further, in vitro assays with the BRET-based Mpro sensors revealed a molecular crowding-mediated increase in the rate of Mpro activity and a decrease in the inhibitory potential of GC376. The sensors developed here will find direct utility in studies related to drug discovery targeting the SARS-CoV-2 Mpro and functional genomics application to determine the effect of sequence variation in Mpro.
Collapse
Affiliation(s)
- Anupriya M. Geethakumari
- Division of Biological and Biomedical Sciences, College of Health & Life Sciences, Hamad Bin Khalifa University, Education City, Qatar Foundation, Doha, 34110 Qatar
| | - Wesam S. Ahmed
- Division of Biological and Biomedical Sciences, College of Health & Life Sciences, Hamad Bin Khalifa University, Education City, Qatar Foundation, Doha, 34110 Qatar
| | - Saad Rasool
- Division of Genomics and Precision Medicine, College of Health & Life Sciences, Hamad Bin Khalifa University, Education City, Qatar Foundation, Doha, 34110 Qatar
| | - Asma Fatima
- Division of Biological and Biomedical Sciences, College of Health & Life Sciences, Hamad Bin Khalifa University, Education City, Qatar Foundation, Doha, 34110 Qatar
| | - S. M. Nasir Uddin
- Division of Biological and Biomedical Sciences, College of Health & Life Sciences, Hamad Bin Khalifa University, Education City, Qatar Foundation, Doha, 34110 Qatar
| | - Mustapha Aouida
- Division of Biological and Biomedical Sciences, College of Health & Life Sciences, Hamad Bin Khalifa University, Education City, Qatar Foundation, Doha, 34110 Qatar
| | - Kabir H. Biswas
- Division of Biological and Biomedical Sciences, College of Health & Life Sciences, Hamad Bin Khalifa University, Education City, Qatar Foundation, Doha, 34110 Qatar
| |
Collapse
|
34
|
Allele-specific activation, enzyme kinetics, and inhibitor sensitivities of EGFR exon 19 deletion mutations in lung cancer. Proc Natl Acad Sci U S A 2022; 119:e2206588119. [PMID: 35867821 PMCID: PMC9335329 DOI: 10.1073/pnas.2206588119] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Epidermal growth factor receptor (EGFR) mutations are detected in approximately 30% of all lung adenocarcinomas, and the most common EGFR mutation occurring in ∼50% of patients is termed “exon 19 deletion” (ex19del). Despite the existence of dozens of different genomic variants comprising what is generically referred to clinically as ex19del, clinicians currently do not distinguish between ex19del variants in considering treatment options, and differences between ex19del variants are largely unstudied in the broader scientific community. Herein, we describe functional differences between distinct EGFR ex19del variants attributable to the structural features of each variant. These findings suggest a possible explanation for observed differences in patient outcomes stratified by ex19del subtype and reinforce the need for allele-specific considerations in clinical treatment decision-making. Oncogenic mutations within the epidermal growth factor receptor (EGFR) are found in 15 to 30% of all non–small-cell lung carcinomas. The term exon 19 deletion (ex19del) is collectively used to refer to more than 20 distinct genomic alterations within exon 19 that comprise the most common EGFR mutation subtype in lung cancer. Despite this heterogeneity, clinical treatment decisions are made irrespective of which EGFR ex19del variant is present within the tumor, and there is a paucity of information regarding how individual ex19del variants influence protein structure and function. Herein, we identified allele-specific functional differences among ex19del variants attributable to recurring sequence and structure motifs. We built all-atom structural models of 60 ex19del variants identified in patients and combined molecular dynamics simulations with biochemical and biophysical experiments to analyze three ex19del mutations (E746_A750, E746_S752 > V, and L747_A750 > P). We demonstrate that sequence variation in ex19del alters oncogenic cell growth, dimerization propensity, enzyme kinetics, and tyrosine kinase inhibitor (TKI) sensitivity. We show that in contrast to E746_A750 and E746_S752 > V, the L747_A750 > P variant forms highly active ligand-independent dimers. Enzyme kinetic analysis and TKI inhibition experiments suggest that E746_S752 > V and L747_A750 > P display reduced TKI sensitivity due to decreased adenosine 5′-triphosphate Km. Through these analyses, we propose an expanded framework for interpreting ex19del variants and considerations for therapeutic intervention.
Collapse
|
35
|
Talluri S. Engineering and Design of Programmable Genome Editors. J Phys Chem B 2022; 126:5140-5150. [PMID: 35819243 DOI: 10.1021/acs.jpcb.2c03761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Programmable genome editors are enzymes that can be targeted to a specific location in the genome for making site-specific alterations or deletions. The engineering, design, and development of sequence-specific editors has resulted in a dramatic increase in the precision of editing for nucleotide sequences. These editors can target specific locations in a genome, in vivo. The genome editors are being deployed for the development of genetically modified organisms for agriculture and industry, and for gene therapy of inherited human genetic disorders, cancer, and immunotherapy. Experimental and computational studies of structure, binding, activity, dynamics, and folding, reviewed here, have provided valuable insights that have the potential for increasing the functional efficiency of these gene/genome editors. Biochemical and biophysical studies of the specificities of natural and engineered genome editors reveal that increased binding affinity can be detrimental because of the increase of off-target effects and that the engineering and design of genome editors with higher specificity may require modulation and control of the conformational dynamics.
Collapse
Affiliation(s)
- Sekhar Talluri
- Department of Biotechnology, GITAM, Visakhapatnam, India 530045
| |
Collapse
|
36
|
Shekhar M, Smith Z, Seeliger MA, Tiwary P. Protein Flexibility and Dissociation Pathway Differentiation Can Explain Onset of Resistance Mutations in Kinases. Angew Chem Int Ed Engl 2022; 61:e202200983. [PMID: 35486370 DOI: 10.1002/anie.202200983] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Indexed: 12/14/2022]
Abstract
Understanding how mutations render a drug ineffective is a problem of immense relevance. Often the mechanism through which mutations cause drug resistance can be explained purely through thermodynamics. However, the more perplexing situation is when two proteins have the same drug binding affinities but different residence times. In this work, we demonstrate how all-atom molecular dynamics simulations using recent developments grounded in statistical mechanics can provide a detailed mechanistic rationale for such variances. We discover dissociation mechanisms for the anti-cancer drug Imatinib (Gleevec) against wild-type and the N368S mutant of Abl kinase. We show how this point mutation triggers far-reaching changes in the protein's flexibility and leads to a different, much faster, drug dissociation pathway. We believe that this work marks an efficient and scalable approach to obtain mechanistic insight into resistance mutations in biomolecular receptors that are hard to explain using a structural perspective.
Collapse
Affiliation(s)
- Mrinal Shekhar
- Center for Development of Therapeutics, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Zachary Smith
- Biophysics Program and Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
| | - Markus A Seeliger
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY 11794-8651, USA
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry and Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA
| |
Collapse
|
37
|
Multifunctional Analysis of Chia Seed (Salvia hispanica L.) Bioactive Peptides Using Peptidomics and Molecular Dynamics Simulations Approaches. Int J Mol Sci 2022; 23:ijms23137288. [PMID: 35806294 PMCID: PMC9266559 DOI: 10.3390/ijms23137288] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Revised: 06/24/2022] [Accepted: 06/27/2022] [Indexed: 02/01/2023] Open
Abstract
Chia seed peptides (CSP) can be a source of multifunctional biopeptides to treat non-communicable diseases. However, interactions and binding affinity involved in targeting specific receptors remains unexplored. In this study, molecular simulation techniques were used as virtual screening of CSP to determine drug-like candidates using a multi-target-directed ligand approach. CSP fraction with the best bioactivities in vitro was sequenced. Then, a prediction model was built using physicochemical descriptors (hydrophobicity, hydrophilicity, intestinal stability, antiangiogenic, antihypertensive, and anti-inflammatory) to calculate potential scores and rank possible biopeptides. Furthermore, molecular dynamics simulations (MDS) and ensemble molecular docking analysis were carried out using four human protein targets (ACE, angiotensin converting enzyme; VEGF, vascular endothelial growth factor; GLUC, glucocorticoid and MINC, mineralocorticoid receptors). Five known-sequence peptides (NNVFYPF, FNIVFPG, SRPWPIDY, QLQRWFR, GSRFDWTR) and five de novo peptides (DFKF, DLRF, FKAF, FRSF, QFRF) had the lowest energy score and higher affinity for ACE and VEGF. The therapeutic effects of these selected peptides can be related to the inhibition of the enzymes involved in angiogenesis and hypertension, due to formation of stable complexes with VEGF and ACE binding sites, respectively. The application of MDS is a good resource for identifying bioactive peptides for future experimental validation.
Collapse
|
38
|
Molecular dynamics simulation guided distal mutation of Thermotoga naphthophila β-glucosidase for significantly enhanced synthesis of galactooligosaccharides and expanded product scope. Int J Biol Macromol 2022; 210:21-32. [DOI: 10.1016/j.ijbiomac.2022.05.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2022] [Revised: 04/13/2022] [Accepted: 05/01/2022] [Indexed: 11/18/2022]
|
39
|
Calvó-Tusell C, Maria-Solano MA, Osuna S, Feixas F. Time Evolution of the Millisecond Allosteric Activation of Imidazole Glycerol Phosphate Synthase. J Am Chem Soc 2022; 144:7146-7159. [PMID: 35412310 PMCID: PMC9052757 DOI: 10.1021/jacs.1c12629] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Deciphering the molecular
mechanisms of enzymatic allosteric regulation
requires the structural characterization of functional states and
also their time evolution toward the formation of the allosterically
activated ternary complex. The transient nature and usually slow millisecond
time scale interconversion between these functional states hamper
their experimental and computational characterization. Here, we combine
extensive molecular dynamics simulations, enhanced sampling techniques,
and dynamical networks to describe the allosteric activation of imidazole
glycerol phosphate synthase (IGPS) from the substrate-free form to
the active ternary complex. IGPS is a heterodimeric bienzyme complex
whose HisH subunit is responsible for hydrolyzing glutamine and delivering
ammonia for the cyclase activity in HisF. Despite significant advances
in understanding the underlying allosteric mechanism, essential molecular
details of the long-range millisecond allosteric activation of IGPS
remain hidden. Without using a priori information
of the active state, our simulations uncover how IGPS, with the allosteric
effector bound in HisF, spontaneously captures glutamine in a catalytically
inactive HisH conformation, subsequently attains a closed HisF:HisH
interface, and finally forms the oxyanion hole in HisH for efficient
glutamine hydrolysis. We show that the combined effector and substrate
binding dramatically decreases the conformational barrier associated
with oxyanion hole formation, in line with the experimentally observed
4500-fold activity increase in glutamine hydrolysis. The allosteric
activation is controlled by correlated time-evolving dynamic networks
connecting the effector and substrate binding sites. This computational
strategy tailored to describe millisecond events can be used to rationalize
the effect of mutations on the allosteric regulation and guide IGPS
engineering efforts.
Collapse
Affiliation(s)
- Carla Calvó-Tusell
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain
| | - Miguel A Maria-Solano
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain.,Global AI Drug Discovery Center, College of Pharmacy and Graduate School of Pharmaceutical Science, Ewha Womans University, 03760 Seoul, Republic of Korea
| | - Sílvia Osuna
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain.,Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Catalonia, Spain
| | - Ferran Feixas
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Universitat de Girona, c/Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain
| |
Collapse
|
40
|
Human Vitamin K Epoxide Reductase as a Target of Its Redox Protein. Int J Mol Sci 2022; 23:ijms23073899. [PMID: 35409257 PMCID: PMC8998853 DOI: 10.3390/ijms23073899] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Revised: 03/29/2022] [Accepted: 03/30/2022] [Indexed: 02/05/2023] Open
Abstract
Human vitamin K epoxide reductase (hVKORC1) enzymatic activity requires an initial activation by a specific redox protein, a less studied step in the hVKORC1 vital cycle. Significant steric conditions must be met by enzymes, being that to adapt their configurations is mandatory for hVKORC1 activation. We studied, by molecular dynamics (MD) simulations, the folding and conformational plasticity of hVKORC1 in its inactive (fully oxidised) state using available structures, crystallographic and from de novo modelling. According to the obtained results, hVKORC1 is a modular protein composed of the stable transmembrane domain (TMD) and intrinsically disordered luminal (L) loop, possessing the great plasticity/adaptability required to perform various steps of the activation process. The docking (HADDOCK) of Protein Disulfide Isomerase (PDI) onto different hVKORC1 conformations clearly indicated that the most interpretable solutions were found on the target closed L-loop form, a prevalent conformation of hVKORC1’s oxidised state. We also suggest that the cleaved L-loop is an appropriate entity to study hVKORC1 recognition/activation by its redox protein. Additionally, the application of hVKORC1 (membrane protein) in aqueous solution is likely to prove to be very useful in practice in either in silico studies or in vitro experiments.
Collapse
|
41
|
Pawnikar S, Miao Y. Mechanism of Peptide Agonist Binding in CXCR4 Chemokine Receptor. Front Mol Biosci 2022; 9:821055. [PMID: 35359589 PMCID: PMC8963245 DOI: 10.3389/fmolb.2022.821055] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 02/14/2022] [Indexed: 01/07/2023] Open
Abstract
Chemokine receptors are key G-protein-coupled receptors (GPCRs) that control cell migration in immune system responses, development of cardiovascular and central nervous systems, and numerous diseases. In particular, the CXCR4 chemokine receptor promotes metastasis, tumor growth and angiogenesis in cancers. CXCR4 is also used as one of the two co-receptors for T-tropic HIV-1 entry into host cells. Therefore, CXCR4 serves as an important therapeutic target for treating cancers and HIV infection. Apart from the CXCL12 endogenous peptide agonist, previous studies suggested that the first 17 amino acids of CXCL12 are sufficient to activate CXCR4. Two 17-residue peptides with positions 1-4 mutated to RSVM and ASLW functioned as super and partial agonists of CXCR4, respectively. However, the mechanism of peptide agonist binding in CXCR4 remains unclear. Here, we have investigated this mechanism through all-atom simulations using a novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method. The Pep-GaMD simulations have allowed us to explore representative binding conformations of each peptide and identify critical low-energy states of CXCR4 activated by the super versus partial peptide agonists. Our simulations have provided important mechanistic insights into peptide agonist binding in CXCR4, which are expected to facilitate rational design of new peptide modulators of CXCR4 and other chemokine receptors.
Collapse
|
42
|
Wang J, Miao Y. Protein-Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD): Characterization of Protein Binding Thermodynamics and Kinetics. J Chem Theory Comput 2022; 18:1275-1285. [PMID: 35099970 DOI: 10.1021/acs.jctc.1c00974] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Protein-protein interactions (PPIs) play key roles in many fundamental biological processes such as cellular signaling and immune responses. However, it has proven challenging to simulate repetitive protein association and dissociation in order to calculate binding free energies and kinetics of PPIs due to long biological timescales and complex protein dynamics. To address this challenge, we have developed a new computational approach to all-atom simulations of PPIs based on a robust Gaussian accelerated molecular dynamics (GaMD) technique. The method, termed "PPI-GaMD", selectively boosts interaction potential energy between protein partners to facilitate their slow dissociation. Meanwhile, another boost potential is applied to the remaining potential energy of the entire system to effectively model the protein's flexibility and rebinding. PPI-GaMD has been demonstrated on a model system of the ribonuclease barnase interactions with its inhibitor barstar. Six independent 2 μs PPI-GaMD simulations have captured repetitive barstar dissociation and rebinding events, which enable calculations of the protein binding thermodynamics and kinetics simultaneously. The calculated binding free energies and kinetic rate constants agree well with the experimental data. Furthermore, PPI-GaMD simulations have provided mechanistic insights into barstar binding to barnase, which involves long-range electrostatic interactions and multiple binding pathways, being consistent with previous experimental and computational findings of this model system. In summary, PPI-GaMD provides a highly efficient and easy-to-use approach for binding free energy and kinetics calculations of PPIs.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
43
|
Do HN, Wang J, Bhattarai A, Miao Y. GLOW: A Workflow Integrating Gaussian-Accelerated Molecular Dynamics and Deep Learning for Free Energy Profiling. J Chem Theory Comput 2022; 18:1423-1436. [PMID: 35200019 PMCID: PMC9773012 DOI: 10.1021/acs.jctc.1c01055] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We introduce a Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and free energy profiling workflow (GLOW) to predict molecular determinants and map free energy landscapes of biomolecules. All-atom GaMD-enhanced sampling simulations are first performed on biomolecules of interest. Structural contact maps are then calculated from GaMD simulation frames and transformed into images for building DL models using a convolutional neural network. Important structural contacts are further determined from DL models of attention maps of the structural contact gradients, which allow us to identify the system reaction coordinates. Finally, free energy profiles are calculated for the selected reaction coordinates through energetic reweighting of the GaMD simulations. We have also successfully demonstrated GLOW for the characterization of activation and allosteric modulation of a G protein-coupled receptor, using the adenosine A1 receptor (A1AR) as a model system. GLOW findings are highly consistent with previous experimental and computational studies of the A1AR, while also providing further mechanistic insights into the receptor function. In summary, GLOW provides a systematic approach to mapping free energy landscapes of biomolecules. The GLOW workflow and its user manual can be downloaded at http://miaolab.org/GLOW.
Collapse
Affiliation(s)
- Hung N. Do
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Jinan Wang
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Apurba Bhattarai
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Yinglong Miao
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047,Corresponding author:
| |
Collapse
|
44
|
Ni D, Liu Y, Kong R, Yu Z, Lu S, Zhang J. Computational elucidation of allosteric communication in proteins for allosteric drug design. Drug Discov Today 2022; 27:2226-2234. [DOI: 10.1016/j.drudis.2022.03.012] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2021] [Revised: 01/22/2022] [Accepted: 03/17/2022] [Indexed: 02/07/2023]
|
45
|
ElHefnawi M, Jo E, Tolba MM, Fares M, Yang J, Shahbaaz M, Windisch MP. Drug repurposing through virtual screening and in vitro validation identifies tigecycline as a novel putative HCV polymerase inhibitor. Virology 2022; 570:9-17. [DOI: 10.1016/j.virol.2022.02.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Revised: 01/25/2022] [Accepted: 02/26/2022] [Indexed: 10/18/2022]
|
46
|
Madhu MK, Debroy A, Murarka RK. Molecular Insights into Phosphorylation-Induced Allosteric Conformational Changes in a β 2-Adrenergic Receptor. J Phys Chem B 2022; 126:1917-1932. [PMID: 35196859 DOI: 10.1021/acs.jpcb.1c08610] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The large conformational flexibility of G protein-coupled receptors (GPCRs) has been a puzzle in structural and pharmacological studies for the past few decades. Apart from structural rearrangements induced by ligands, enzymatic phosphorylations by GPCR kinases (GRKs) at the carboxy-terminal tail (C-tail) of a GPCR also make conformational alterations to the transmembrane helices and facilitates the binding of one of its transducer proteins named β-arrestin. The phosphorylation-induced conformational transition of the receptor that causes specific binding to β-arrestin but prevents the association of other transducers such as G proteins lacks atomistic understanding and is elusive to experimental studies. Using microseconds of all-atom conventional and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigate the allosteric mechanism of phosphorylation induced-conformational changes in β2-adrenergic receptor, a well-characterized GPCR model system. Free energy profiles reveal that the phosphorylated receptor samples a new conformational state in addition to the canonical active state corroborating with recent nuclear magnetic resonance experimental findings. The new state has a smaller intracellular cavity that is likely to accommodate β-arrestin better than G protein. Using contact map and inter-residue interaction energy calculations, we found the phosphorylated C-tail adheres to the cytosolic surface of the transmembrane domain of the receptor. Transfer entropy calculations show that the C-tail residues drive the correlated motions of TM residues, and the allosteric signal is relayed via several residues at the cytosolic surface. Our results also illustrate how the redistribution of inter-residue nonbonding interaction couples with the allosteric communication from the phosphorylated C-tail to the transmembrane. Atomistic insight into phosphorylation-induced β-arrestin specific conformation is therapeutically important to design drugs with higher efficacy and fewer side effects. Our results, therefore, open novel opportunities to fine-tune β-arrestin bias in GPCR signaling.
Collapse
Affiliation(s)
- Midhun K Madhu
- Department of Biological Sciences, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| | - Annesha Debroy
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| | - Rajesh K Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India
| |
Collapse
|
47
|
Vidal-Limon A, Aguilar-Toalá JE, Liceaga AM. Integration of Molecular Docking Analysis and Molecular Dynamics Simulations for Studying Food Proteins and Bioactive Peptides. JOURNAL OF AGRICULTURAL AND FOOD CHEMISTRY 2022; 70:934-943. [PMID: 34990125 DOI: 10.1021/acs.jafc.1c06110] [Citation(s) in RCA: 95] [Impact Index Per Article: 47.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
In silico tools, such as molecular docking, are widely applied to study interactions and binding affinity of biological activity of proteins and peptides. However, restricted sampling of both ligand and receptor conformations and use of approximated scoring functions can produce results that do not correlate with actual experimental binding affinities. Molecular dynamics simulations (MDS) can provide valuable information in deciphering functional mechanisms of proteins/peptides and other biomolecules, overcoming the rigid sampling limitations in docking analysis. This review will discuss the information related to the traditional use of in silico models, such as molecular docking, and its application for studying food proteins and bioactive peptides, followed by an in-depth introduction to the theory of MDS and description of why these molecular simulation techniques are important in the theoretical prediction of structural and functional dynamics of food proteins and bioactive peptides. Applications, limitations, and future prospects of MDS will also be discussed.
Collapse
Affiliation(s)
- Abraham Vidal-Limon
- Red de Estudios Moleculares Avanzados, Clúster Científico y Tecnológico BioMimic, Instituto de Ecología A.C. (INECOL), Carretera Antigua a Coatepec 351, El Haya, Xalapa, Veracruz 91073, Mexico
| | - José E Aguilar-Toalá
- Departamento de Ciencias de la Alimentación, División de Ciencias Biológicas y de la Salud, Universidad Autónoma Metropolitana Unidad Lerma, Avenida de las Garzas 10, Colonia El Panteón, Lerma de Villada, Estado de México 52005, Mexico
| | - Andrea M Liceaga
- Protein Chemistry and Bioactive Peptides Laboratory. Department of Food Science, Purdue University, 745 Agriculture Mall Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
48
|
Célerse F, Inizan TJ, Lagardère L, Adjoua O, Monmarché P, Miao Y, Derat E, Piquemal JP. An Efficient Gaussian-Accelerated Molecular Dynamics (GaMD) Multilevel Enhanced Sampling Strategy: Application to Polarizable Force Fields Simulations of Large Biological Systems. J Chem Theory Comput 2022; 18:968-977. [PMID: 35080892 DOI: 10.1021/acs.jctc.1c01024] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We introduce a novel multilevel enhanced sampling strategy grounded on Gaussian-accelerated Molecular Dynamics (GaMD). First, we propose a GaMD multi-GPUs-accelerated implementation within the Tinker-HP molecular dynamics package. We introduce the new "dual-water" mode and its use with the flexible AMOEBA polarizable force field. By adding harmonic boosts to the water stretching and bonding terms, it accelerates the solvent-solute interactions while enabling speedups, thanks to the use of fast multiple-time step integrators. To further reduce the time-to-solution, we couple GaMD to Umbrella Sampling (US). The GaMD─US/dual-water approach is tested on the 1D Potential of Mean Force (PMF) of the solvated CD2-CD58 system (168 000 atoms), allowing the AMOEBA PMF to converge within 1 kcal/mol of the experimental value. Finally, Adaptive Sampling (AS) is added, enabling AS-GaMD capabilities but also the introduction of the new Adaptive Sampling-US-GaMD (ASUS-GaMD) scheme. The highly parallel ASUS-GaMD setup decreases time to convergence by, respectively, 10 and 20 times, compared to GaMD-US and US. Overall, beside the acceleration of PMF computations, Tinker-HP now allows for the simultaneous use of Adaptive Sampling and GaMD-"dual water" enhanced sampling approaches increasing the applicability of polarizable force fields to large-scale simulations of biological systems.
Collapse
Affiliation(s)
- Frédéric Célerse
- LCT, UMR 7616 CNRS, Sorbonne Université, Paris 75005, France.,IPCM, UMR 8232 CNRS, Sorbonne Université, Paris 75005, France
| | | | - Louis Lagardère
- LCT, UMR 7616 CNRS, Sorbonne Université, Paris 75005, France.,IP2CT, FR 2622 CNRS, Sorbonne Université, Paris 75005, France
| | - Olivier Adjoua
- LCT, UMR 7616 CNRS, Sorbonne Université, Paris 75005, France
| | - Pierre Monmarché
- LCT, UMR 7616 CNRS, Sorbonne Université, Paris 75005, France.,LJLL, UMR 7598 CNRS, Sorbonne Université, Paris 75005, France
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66045, United States
| | - Etienne Derat
- IPCM, UMR 8232 CNRS, Sorbonne Université, Paris 75005, France
| | - Jean-Philip Piquemal
- LCT, UMR 7616 CNRS, Sorbonne Université, Paris 75005, France.,The University of Texas at Austin, Department of Biomedical Engineering, Austin, Texas 78705, United States.,Institut Universitaire de France, Paris 75005, France
| |
Collapse
|
49
|
Wang Y, Li M, Liang W, Shi X, Fan J, Kong R, Liu Y, Zhang J, Chen T, Lu S. Delineating the activation mechanism and conformational landscape of a class B G protein-coupled receptor glucagon receptor. Comput Struct Biotechnol J 2022; 20:628-639. [PMID: 35140883 PMCID: PMC8801358 DOI: 10.1016/j.csbj.2022.01.015] [Citation(s) in RCA: 33] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Revised: 01/17/2022] [Accepted: 01/17/2022] [Indexed: 02/09/2023] Open
Abstract
Class B G protein-coupled receptors (GPCRs) are important targets in the treatment of metabolic syndrome and diabetes. Although multiple structures of class B GPCRs-G protein complexes have been elucidated, the detailed activation mechanism of the receptors remains unclear. Here, we combine Gaussian accelerated molecular dynamics simulations and Markov state models (MSM) to investigate the activation mechanism of a canonical class B GPCR, human glucagon receptor-GCGR, including the negative allosteric modulator-bound inactive state, the agonist glucagon-bound active state, and both glucagon- and Gs-bound fully active state. The free-energy landscapes of GCGR show the conformational ensemble consisting of three activation-associated states: inactive, active, and fully active. The structural analysis indicates the high dynamics of GCGR upon glucagon binding with both active and inactive conformations in the ensemble. Significantly, the H8 and TM6 exhibits distinct features from the inactive to the active states. The additional simulations demonstrate the role of H8 in the recruitment of Gs. Gs binding presents a crucial function of stabilizing the glucagon binding site and MSM highlights the absolute requirement of Gs to help the GCGR reach the fully active state. Together, our results reveal the detailed activation mechanism of GCGR from the view of conformational dynamics.
Collapse
Affiliation(s)
- Ying 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
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China
| | - Mingyu 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
| | - Wenqi Liang
- Department of Emergency, Changhai Hospital, Naval Medical University, Shanghai 200433, China
| | - Xinchao Shi
- 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
| | - Jigang Fan
- 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
| | - Ren Kong
- Institute of Bioinformatics and Medical Engineering, School of Electrical and Information Engineering, Jiangsu University of Technology, Changzhou 213001, China
| | - Yaqin Liu
- Medicinal Chemistry and Bioinformatics Center, 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
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China
| | - Ting Chen
- Department of Cardiology, Changzheng Hospital, Naval Medical University, Shanghai 200023, 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
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University, School of Medicine, Shanghai 200025, China
| |
Collapse
|
50
|
Pawnikar S, Bhattarai A, Wang J, Miao Y. Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives. Adv Appl Bioinform Chem 2022; 15:1-19. [PMID: 35023931 PMCID: PMC8747661 DOI: 10.2147/aabc.s247950] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2020] [Accepted: 02/20/2021] [Indexed: 12/12/2022] Open
Abstract
Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.
Collapse
Affiliation(s)
- Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| |
Collapse
|