1
|
Clayton J, Romany A, Matenoglou E, Gavathiotis E, Poulikakos PI, Shen J. Mechanism of Dimer Selectivity and Binding Cooperativity of BRAF Inhibitors. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.12.12.571293. [PMID: 38168366 PMCID: PMC10760002 DOI: 10.1101/2023.12.12.571293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
Aberrant signaling of BRAF V600E is a major cancer driver. Current FDA-approved RAF inhibitors selectively inhibit the monomeric BRAF V600E and suffer from tumor resistance. Recently, dimer-selective and equipotent RAF inhibitors have been developed; however, the mechanism of dimer selectivity is poorly understood. Here, we report extensive molecular dynamics (MD) simulations of the monomeric and dimeric BRAF V600E in the apo form or in complex with one or two dimer-selective (PHI1) or equipotent (LY3009120) inhibitor(s). The simulations uncovered the unprecedented details of the remarkable allostery in BRAF V600E dimerization and inhibitor binding. Specifically, dimerization retrains and shifts the αC helix inward and increases the flexibility of the DFG motif; dimer compatibility is due to the promotion of the αC-in conformation, which is stabilized by a hydrogen bond formation between the inhibitor and the αC Glu501. A more stable hydrogen bond further restrains and shifts the αC helix inward, which incurs a larger entropic penalty that disfavors monomer binding. This mechanism led us to propose an empirical way based on the co-crystal structure to assess the dimer selectivity of a BRAF V600E inhibitor. Simulations also revealed that the positive cooperativity of PHI1 is due to its ability to preorganize the αC and DFG conformation in the opposite protomer, priming it for binding the second inhibitor. The atomically detailed view of the interplay between BRAF dimerization and inhibitor allostery as well as cooperativity has implications for understanding kinase signaling and contributes to the design of protomer selective RAF inhibitors.
Collapse
|
2
|
Peeples CA, Liu R, Shen J. Force Field Limitations of All-Atom Continuous Constant pH Molecular Dynamics. J Phys Chem B 2024; 128:11616-11624. [PMID: 39531617 DOI: 10.1021/acs.jpcb.4c05971] [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: 11/16/2024]
Abstract
All-atom constant pH molecular dynamics simulations offer a powerful tool for understanding pH-mediated and proton-coupled biological processes. As the protonation equilibria of protein side chains are shifted by electrostatic interactions and desolvation energies, pKa values calculated from the constant pH simulations may be sensitive to the underlying protein force field and water model. Here we investigated the force field dependence of the all-atom particle mesh Ewald (PME) continuous constant pH (PME-CpHMD) simulations of a mini-protein BBL. The replica-exchange titration simulations based on the Amber ff19sb and ff14sb force fields with the respective water models showed significantly overestimated pKa downshifts for a buried histidine (His166) and for two glutamic acids (Glu141 and Glu161) that are involved in salt-bridge interactions. These errors (due to undersolvation of neutral histidines and overstabilization of salt bridges) are consistent with the previously reported pKa's based on the CHARMM c22/CMAP force field, albeit in larger magnitudes. The pKa calculations also demonstrated that ff19sb with OPC water is significantly more accurate than ff14sb with TIP3P water, and the salt-bridge related pKa downshifts can be partially alleviated by the atom-pair specific Lennard-Jones corrections (NBFIX). Together, these data suggest that the accuracies of the protonation equilibria of proteins from constant pH simulations can significantly benefit from improvements of force fields.
Collapse
Affiliation(s)
- Craig A Peeples
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Ruibin Liu
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| |
Collapse
|
3
|
Li P, Shi M, Wang Y, Liu Q, Du X, Wang X. pH-Dependent Assembly and Stability of Toll-Like Receptor 3/dsRNA Signaling Complex: Insights from Constant pH Molecular Dynamics and Metadynamics Simulations. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2024:e2411445. [PMID: 39520076 DOI: 10.1002/advs.202411445] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/17/2024] [Revised: 10/26/2024] [Indexed: 11/16/2024]
Abstract
The pH-dependent assembly of Toll-like receptors (TLRs), which triggers a threshold-like response, is a key principle in immune signaling. While crystallography has revealed the intricate structure of these assembly complexes, the mechanisms underlying their pH dependency remain unclear. Herein, constant pH simulations and metadynamics are employed to investigate the pH-dependent assembly and stability of the TLR3/dsRNA signaling complex. The findings demonstrate that system pH regulates complex assembly and stability by modulating the protonation and charge states of histidines. Histidines in TLR3 act as pH-dependent, positively charged binding sites that capture negatively charged dsRNA. Additionally, these histidines form a [H682⁺]-[E626⁻] dipole, facilitating the assembly of two TLR3 molecules into an antisymmetric dimer through dipole-dipole interactions. Surprisingly, TLR3 can shift the pKa values of key histidines from their model pKa of 6.5, increasing protonation likelihood and enhancing ligand binding. Notably, the aromatic residue Phe84, located within the dsRNA binding site [His39⁺-His60⁺-Phe84-His108⁺], alters the pKa of His60 through cation-π interactions with its protonated state. This study offers new insights into the molecular mechanisms underlying pH-dependent immune signaling via higher-order assemblies and suggests potential applications for histidine in self-assembling biomaterials.
Collapse
Affiliation(s)
- Penghui Li
- Shenzhen Key Laboratory of Marine Biotechnology and Ecology, College of Life Sciences & Oceanography, Shenzhen University, Shenzhen, 518055, China
- Key Laboratory of Optoelectronic Devices and System of Ministry of Education and Guangdong Province, College Physics and Optoelectronic Engineering, Shenzhen University, Shenzhen, 518060, China
| | - Mingsong Shi
- NHC Key Laboratory of Nuclear Technology Medical Transformation, Mianyang Central Hospital, School of Medicine, University of Electronic Science and Technology of China, Mianyang, Sichuan, 621099, China
| | - Yibo Wang
- Laboratory of Chemical Biology, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, China
| | - Qiong Liu
- Shenzhen Key Laboratory of Marine Biotechnology and Ecology, College of Life Sciences & Oceanography, Shenzhen University, Shenzhen, 518055, China
- Key Laboratory of Optoelectronic Devices and System of Ministry of Education and Guangdong Province, College Physics and Optoelectronic Engineering, Shenzhen University, Shenzhen, 518060, China
- Shenzhen-Hong Kong Institute of Brain Science, Shenzhen Fundamental Research Institutions, Shenzhen, 518055, China
| | - Xiubo Du
- Shenzhen Key Laboratory of Marine Biotechnology and Ecology, College of Life Sciences & Oceanography, Shenzhen University, Shenzhen, 518055, China
| | - Xiaohui Wang
- Laboratory of Chemical Biology, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, China
- School of Applied Chemistry and Engineering, University of Science and Technology of China, Hefei, 230026, China
| |
Collapse
|
4
|
Peeples CA, Liu R, Shen J. Force Field Limitations of All-Atom Continuous Constant pH Molecular Dynamics. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.09.03.611076. [PMID: 39282392 PMCID: PMC11398383 DOI: 10.1101/2024.09.03.611076] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 09/21/2024]
Abstract
All-atom constant pH molecular dynamics simulations offer a powerful tool for understanding pH-mediated and proton-coupled biological processes. As the protonation equilibria of protein sidechains are shifted by electrostatic interactions and desolvation energies, pK a values calculated from the constant pH simulations may be sensitive to the underlying protein force field and water model. Here we investigated the force field dependence of the all-atom particle mesh Ewald (PME) continuous constant pH (PME-CpHMD) simulations of a mini-protein BBL. The replica-exchange titration simulations based on the Amber ff19sb and ff14sb force fields with the respective water models showed significantly overestimated pK a downshifts for a buried histidine (His166) and for two glutamic acids (Glu141 and Glu161) that are involved in salt-bridge interactions. These errors (due to undersolvation of neutral histidines and overstabilization of salt bridges) are consistent with the previously reported pK a's based on the CHARMM c22/CMAP force field, albeit in larger magnitudes. The pK a calculations also demonstrated that ff19sb with OPC water is significantly more accurate than ff14sb with TIP3P water, and the salt-bridge related pK a downshifts can be partially alleviated by the atom-pair specific Lennard-Jones corrections (NBFIX). Together, these data suggest that the accuracies of the protonation equilibria of proteins from constant pH simulations can significantly benefit from improvements of force fields.
Collapse
Affiliation(s)
- Craig A Peeples
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, MD 21201
| | - Ruibin Liu
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, MD 21201
| | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, MD 21201
| |
Collapse
|
5
|
Van Dyck PK, Piszkin L, Gorski EA, Nascimento ET, Abebe JA, Hoffmann LM, Peng JW, White KA. Ionizable networks mediate pH-dependent allostery in SH2 signaling proteins. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.08.21.608875. [PMID: 39229188 PMCID: PMC11370553 DOI: 10.1101/2024.08.21.608875] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
IntroductionTransient intracellular pH dynamics1regulate mammalian proliferation2,3, migration4, and differentiation5. However, for many pH-dependent cell processes, the molecular mediators are unknown6. Prior work identified histidine residues as molecular switches in pH-sensitive proteins, but how other ionizable residues contribute to pH-dependent protein allostery is understudied. Here, we develop anin silicocomputational pipeline to identify putative pH-sensitive proteins and their molecular mechanisms. We first apply this pipeline to SHP2, a known pH-sensitive signaling protein with an uncharacterized molecular mechanism. We show wild-type SHP2 phosphatase activity is pH-sensitivein vitroand in cells, and mutation of identified H116 and E252 to non-titratable alanine residues abolishes pH-sensitive function. We also show that c-Src is a previously unrecognized pH-dependent kinase, and mutation of the identified ionizable network again abolishes pH-sensitive activity. Constant pH molecular dynamics simulations support a conserved allosteric mechanism of pH-dependent binding of inhibitory SH2 domains to the functional catalytic domains of SHP2 and c-Src. We apply our computational pipeline across SH2 domain-containing signaling proteins and identify evolutionarily conserved putative pH-sensing networks. Our results reveal that pH is an allosteric regulator of SH2 domain-containing signaling proteins providing insight into normal pH-dependent cell biology and diseases where pHi is dysregulated, such as cancer.
Collapse
Affiliation(s)
- Papa Kobina Van Dyck
- Department of Chemistry and Biochemistry, University of Notre Dame 251 Nieuwland Science Hall Notre Dame, IN 46556 USA
- Harper Cancer Research Institute, University of Notre Dame 1234 N. Notre Dame Avenue South Bend, IN 46617 USA
| | - Luke Piszkin
- Department of Chemistry and Biochemistry, University of Notre Dame 251 Nieuwland Science Hall Notre Dame, IN 46556 USA
| | - Elijah A. Gorski
- Department of Chemistry and Biochemistry, University of Notre Dame 251 Nieuwland Science Hall Notre Dame, IN 46556 USA
- Harper Cancer Research Institute, University of Notre Dame 1234 N. Notre Dame Avenue South Bend, IN 46617 USA
| | - Eduarda Tartarella Nascimento
- Department of Chemistry and Biochemistry, University of Notre Dame 251 Nieuwland Science Hall Notre Dame, IN 46556 USA
- Harper Cancer Research Institute, University of Notre Dame 1234 N. Notre Dame Avenue South Bend, IN 46617 USA
- Saint Mary’s College, Notre Dame, IN 46556 USA
| | - Joshua A. Abebe
- Department of Chemistry and Biochemistry, University of Notre Dame 251 Nieuwland Science Hall Notre Dame, IN 46556 USA
- Harper Cancer Research Institute, University of Notre Dame 1234 N. Notre Dame Avenue South Bend, IN 46617 USA
| | - Logan M. Hoffmann
- Department of Chemistry and Biochemistry, University of Notre Dame 251 Nieuwland Science Hall Notre Dame, IN 46556 USA
- Harper Cancer Research Institute, University of Notre Dame 1234 N. Notre Dame Avenue South Bend, IN 46617 USA
| | - Jeffrey W. Peng
- Department of Chemistry and Biochemistry, University of Notre Dame 251 Nieuwland Science Hall Notre Dame, IN 46556 USA
| | - Katharine A. White
- Department of Chemistry and Biochemistry, University of Notre Dame 251 Nieuwland Science Hall Notre Dame, IN 46556 USA
- Harper Cancer Research Institute, University of Notre Dame 1234 N. Notre Dame Avenue South Bend, IN 46617 USA
| |
Collapse
|
6
|
Zlobin A, Smirnov I, Golovin A. Dynamic interchange between two protonation states is characteristic of active sites of cholinesterases. Protein Sci 2024; 33:e5100. [PMID: 39022909 PMCID: PMC11255601 DOI: 10.1002/pro.5100] [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: 02/22/2024] [Revised: 05/28/2024] [Accepted: 06/19/2024] [Indexed: 07/20/2024]
Abstract
Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.
Collapse
Affiliation(s)
- Alexander Zlobin
- Institute for Drug DiscoveryLeipzig University Medical SchoolLeipzigGermany
- Faculty of Bioengineering and BioinformaticsLomonosov Moscow State UniversityMoscowRussia
| | - Ivan Smirnov
- Shemyakin and Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of SciencesMoscowRussia
| | - Andrey Golovin
- Faculty of Bioengineering and BioinformaticsLomonosov Moscow State UniversityMoscowRussia
- Shemyakin and Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of SciencesMoscowRussia
- Belozersky Institute of Physico‐Chemical BiologyLomonosov Moscow State UniversityMoscowRussia
| |
Collapse
|
7
|
Zheng R, Zhao M, Du JS, Sudarshan TR, Zhou Y, Paravastu AK, De Yoreo JJ, Ferguson AL, Chen CL. Assembly of short amphiphilic peptoids into nanohelices with controllable supramolecular chirality. Nat Commun 2024; 15:3264. [PMID: 38627405 PMCID: PMC11021492 DOI: 10.1038/s41467-024-46839-y] [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/14/2023] [Accepted: 03/12/2024] [Indexed: 04/19/2024] Open
Abstract
A long-standing challenge in bioinspired materials is to design and synthesize synthetic materials that mimic the sophisticated structures and functions of natural biomaterials, such as helical protein assemblies that are important in biological systems. Herein, we report the formation of a series of nanohelices from a type of well-developed protein-mimetics called peptoids. We demonstrate that nanohelix structures and supramolecular chirality can be well-controlled through the side-chain chemistry. Specifically, the ionic effects on peptoids from varying the polar side-chain groups result in the formation of either single helical fiber or hierarchically stacked helical bundles. We also demonstrate that the supramolecular chirality of assembled peptoid helices can be controlled by modifying assembling peptoids with a single chiral amino acid side chain. Computational simulations and theoretical modeling predict that minimizing exposure of hydrophobic domains within a twisted helical form presents the most thermodynamically favorable packing of these amphiphilic peptoids and suggests a key role for both polar and hydrophobic domains on nanohelix formation. Our findings establish a platform to design and synthesize chiral functional materials using sequence-defined synthetic polymers.
Collapse
Affiliation(s)
- Renyu Zheng
- Department of Chemical Engineering, University of Washington, Seattle, WA, 98195, USA
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, WA, 99352, USA
| | - Mingfei Zhao
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, 60637, USA
| | - Jingshan S Du
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, WA, 99352, USA
| | - Tarunya Rao Sudarshan
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, GA, 30332, USA
| | - Yicheng Zhou
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, WA, 99352, USA
| | - Anant K Paravastu
- School of Chemical and Biomolecular Engineering, Georgia Institute of Technology, Atlanta, GA, 30332, USA
- Parker H. Petit Institute for Bioengineering and Biosciences, Georgia Institute of Technology, Atlanta, GA, 30332, USA
| | - James J De Yoreo
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, WA, 99352, USA
- Department of Materials Science, University of Washington, Seattle, WA, 98195, USA
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, 60637, USA
| | - Chun-Long Chen
- Department of Chemical Engineering, University of Washington, Seattle, WA, 98195, USA.
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, WA, 99352, USA.
| |
Collapse
|
8
|
Omar MN, Rahman RNZRA, Noor NDM, Latip W, Knight VF, Ali MSM. Exploring the Antarctic aminopeptidase P from Pseudomonas sp. strain AMS3 through structural analysis and molecular dynamics simulation. J Biomol Struct Dyn 2024:1-13. [PMID: 38555730 DOI: 10.1080/07391102.2024.2331093] [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: 11/30/2023] [Accepted: 03/11/2024] [Indexed: 04/02/2024]
Abstract
Aminopeptidase P (APPro) is a crucial metalloaminopeptidase involved in amino acid cleavage from peptide N-termini, playing essential roles as versatile biocatalysts with applications ranging from pharmaceuticals to industrial processes. Despite acknowledging its potential for catalysis in lower temperatures, detailed molecular basis and biotechnological implications in cold environments are yet to be explored. Therefore, this research aims to investigate the molecular mechanisms underlying the cold-adapted characteristics of APPro from Pseudomonas sp. strain AMS3 (AMS3-APPro) through a detailed analysis of its structure and dynamics. In this study, structure analysis and molecular dynamics (MD) simulation of a predicted model of AMS3-APPro has been performed at different temperatures to assess structural flexibility and thermostability across a temperature range of 0-60 °C over 100 ns. The MD simulation results revealed that the structure were able to remain stable at low temperatures. Increased temperatures present a potential threat to the overall stability of AMS3-APPro by disrupting the intricate hydrogen bond networks crucial for maintaining structural integrity, thereby increasing the likelihood of protein unfolding. While the metal binding site at the catalytic core exhibits resilience at higher temperatures, highlighting its local structural integrity, the overall enzyme structure undergoes fluctuations and potential denaturation. This extensive structural instability surpasses the localized stability observed at the metal binding site. Consequently, these assessments offer in-depth understanding of the cold-adapted characteristics of AMS3-APPro, highlighting its capability to uphold its native conformation and stability in low-temperature environments. In summary, this research provides valuable insights into the cold-adapted features of AMS3-APPro, suggesting its efficient operation in low thermal conditions, particularly relevant for potential biotechnological applications in cold environments.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Muhamad Nadzmi Omar
- Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
- Department of Biochemistry, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
| | - Raja Noor Zaliha Raja Abd Rahman
- Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
- Department of Biochemistry, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
| | - Noor Dina Muhd Noor
- Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
- Department of Biochemistry, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
| | - Wahhida Latip
- Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
| | - Victor Feizal Knight
- Research Centre for Chemical Defence, National Defence University of Malaysia, Kem Perdana Sungai Besi, Kuala Lumpur, Malaysia
| | - Mohd Shukuri Mohamad Ali
- Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
- Department of Biochemistry, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia
| |
Collapse
|
9
|
Lasham J, Djurabekova A, Zickermann V, Vonck J, Sharma V. Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures. J Phys Chem B 2024; 128:2304-2316. [PMID: 38430110 PMCID: PMC11389979 DOI: 10.1021/acs.jpcb.3c07421] [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: 03/03/2024]
Abstract
Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.
Collapse
Affiliation(s)
- Jonathan Lasham
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland
| | - Amina Djurabekova
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland
| | - Volker Zickermann
- Institute of Biochemistry II, University Hospital, Goethe University, 60590 Frankfurt am Main, Germany
- Centre for Biomolecular Magnetic Resonance, Institute for Biophysical Chemistry, Goethe University, 60438 Frankfurt am Main, Germany
| | - Janet Vonck
- Department of Structural Biology, Max Planck Institute of Biophysics, 60438 Frankfurt am Main, Germany
| | - Vivek Sharma
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland
- HiLIFE Institute of Biotechnology, University of Helsinki, 00014 Helsinki, Finland
| |
Collapse
|
10
|
Jansen A, Aho N, Groenhof G, Buslaev P, Hess B. phbuilder: A Tool for Efficiently Setting up Constant pH Molecular Dynamics Simulations in GROMACS. J Chem Inf Model 2024; 64:567-574. [PMID: 38215282 PMCID: PMC10865341 DOI: 10.1021/acs.jcim.3c01313] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 12/05/2023] [Accepted: 12/06/2023] [Indexed: 01/14/2024]
Abstract
Constant pH molecular dynamics (MD) is a powerful technique that allows the protonation state of residues to change dynamically, thereby enabling the study of pH dependence in a manner that has not been possible before. Recently, a constant pH implementation was incorporated into the GROMACS MD package. Although this implementation provides good accuracy and performance, manual modification and the preparation of simulation input files are required, which can be complicated, tedious, and prone to errors. To simplify and automate the setup process, we present phbuilder, a tool that automatically prepares constant pH MD simulations for GROMACS by modifying the input structure and topology as well as generating the necessary parameter files. phbuilder can prepare constant pH simulations from both initial structures and existing simulation systems, and it also provides functionality for performing titrations and single-site parametrizations of new titratable group types. The tool is freely available at www.gitlab.com/gromacs-constantph. We anticipate that phbuilder will make constant pH simulations easier to set up, thereby making them more accessible to the GROMACS user community.
Collapse
Affiliation(s)
- Anton Jansen
- Department
of Applied Physics and Swedish e-Science Research Center, Science
for Life Laboratory, KTH Royal Institute
of Technology, 100 44 Stockholm, Sweden
| | - Noora Aho
- Nanoscience
Center and Department of Chemistry, University
of Jyväskylä, 40014 Jyväskylä, Finland
| | - Gerrit Groenhof
- Nanoscience
Center and Department of Chemistry, University
of Jyväskylä, 40014 Jyväskylä, Finland
| | - Pavel Buslaev
- Nanoscience
Center and Department of Chemistry, University
of Jyväskylä, 40014 Jyväskylä, Finland
| | - Berk Hess
- Department
of Applied Physics and Swedish e-Science Research Center, Science
for Life Laboratory, KTH Royal Institute
of Technology, 100 44 Stockholm, Sweden
| |
Collapse
|
11
|
Clayton J, de Oliveira VM, Ibrahim MF, Sun X, Mahinthichaichan P, Shen M, Hilgenfeld R, Shen J. Integrative Approach to Dissect the Drug Resistance Mechanism of the H172Y Mutation of SARS-CoV-2 Main Protease. J Chem Inf Model 2023; 63:3521-3533. [PMID: 37199464 PMCID: PMC10237302 DOI: 10.1021/acs.jcim.3c00344] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Indexed: 05/19/2023]
Abstract
Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of Paxlovid, a drug approved by the U.S. Food and Drug Administration for high-risk COVID-19 patients. Recently, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir's inhibitory activity. As the COVID-19 cases skyrocket in China and the selective pressure of antiviral therapy builds in the US, there is an urgent need to characterize and understand how the H172Y mutation confers drug resistance. Here, we investigated the H172Y Mpro's conformational dynamics, folding stability, catalytic efficiency, and inhibitory activity using all-atom constant pH and fixed-charge molecular dynamics simulations, alchemical and empirical free energy calculations, artificial neural networks, and biochemical experiments. Our data suggest that the mutation significantly weakens the S1 pocket interactions with the N-terminus and perturbs the conformation of the oxyanion loop, leading to a decrease in the thermal stability and catalytic efficiency. Importantly, the perturbed S1 pocket dynamics weaken the nirmatrelvir binding in the P1 position, which explains the decreased inhibitory activity of nirmatrelvir. Our work demonstrates the predictive power of the combined simulation and artificial intelligence approaches, and together with biochemical experiments, they can be used to actively surveil continually emerging mutations of SARS-CoV-2 Mpro and assist the optimization of antiviral drugs. The presented approach, in general, can be applied to characterize mutation effects on any protein drug targets.
Collapse
Affiliation(s)
- Joseph Clayton
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, USA
| | - Vinicius Martins de Oliveira
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, USA
| | | | - Xinyuanyuan Sun
- Institute of Molecular Medicine, University of Lübeck, Lübeck 23562, Germany
| | - Paween Mahinthichaichan
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, USA
| | - Mingzhe Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, USA
| | - Rolf Hilgenfeld
- Institute for Molecular Medicine, University of Lübeck, Lübeck 23562, Germany
- German Center for Infection Research (DZIF), Hamburg – Lübeck – Borstel – Riems Site, University of Lübeck, Lübeck 23562, Germany
| | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, USA
| |
Collapse
|
12
|
Cai Z, Liu T, Lin Q, He J, Lei X, Luo F, Huang Y. Basis for Accurate Protein p Ka Prediction with Machine Learning. J Chem Inf Model 2023; 63:2936-2947. [PMID: 37146199 DOI: 10.1021/acs.jcim.3c00254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
pH regulates protein structures and the associated functions in many biological processes via protonation and deprotonation of ionizable side chains where the titration equilibria are determined by pKa's. To accelerate pH-dependent molecular mechanism research in the life sciences or industrial protein and drug designs, fast and accurate pKa prediction is crucial. Here we present a theoretical pKa data set PHMD549, which was successfully applied to four distinct machine learning methods, including DeepKa, which was proposed in our previous work. To reach a valid comparison, EXP67S was selected as the test set. Encouragingly, DeepKa was improved significantly and outperforms other state-of-the-art methods, except for the constant-pH molecular dynamics, which was utilized to create PHMD549. More importantly, DeepKa reproduced experimental pKa orders of acidic dyads in five enzyme catalytic sites. Apart from structural proteins, DeepKa was found applicable to intrinsically disordered peptides. Further, in combination with solvent exposures, it is revealed that DeepKa offers the most accurate prediction under the challenging circumstance that hydrogen bonding or salt bridge interaction is partly compensated by desolvation for a buried side chain. Finally, our benchmark data qualify PHMD549 and EXP67S as the basis for future developments of protein pKa prediction tools driven by artificial intelligence. In addition, DeepKa built on PHMD549 has been proven an efficient protein pKa predictor and thus can be applied immediately to, for example, pKa database construction, protein design, drug discovery, and so on.
Collapse
Affiliation(s)
- Zhitao Cai
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Tengzi Liu
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Qiaoling Lin
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Jiahao He
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Xiaowei Lei
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Fangfang Luo
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Yandong Huang
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| |
Collapse
|
13
|
Shahab M, Danial M, Khan T, Liang C, Duan X, Wang D, Gao H, Zheng G. In Silico Identification of Lead Compounds for Pseudomonas Aeruginosa PqsA Enzyme: Computational Study to Block Biofilm Formation. Biomedicines 2023; 11:961. [PMID: 36979940 PMCID: PMC10046026 DOI: 10.3390/biomedicines11030961] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 03/09/2023] [Accepted: 03/14/2023] [Indexed: 03/30/2023] Open
Abstract
Pseudomonas aeruginosa is an opportunistic Gram-negative bacterium implicated in acute and chronic nosocomial infections and a leading cause of patient mortality. Pseudomonas aeruginosa infections are frequently associated with the development of biofilms, which give the bacteria additional drug resistance and increase their virulence. The goal of this study was to find strong compounds that block the Anthranilate-CoA ligase enzyme made by the pqsA gene. This would stop the P. aeruginosa quorum signaling system. This enzyme plays a crucial role in the pathogenicity of P. aeruginosa by producing autoinducers for cell-to-cell communication that lead to the production of biofilms. Pharmacophore-based virtual screening was carried out utilizing a library of commercially accessible enzyme inhibitors. The most promising hits obtained during virtual screening were put through molecular docking with the help of MOE. The virtual screening yielded 7/160 and 10/249 hits (ZINC and Chembridge). Finally, 2/7 ZINC hits and 2/10 ChemBridge hits were selected as potent lead compounds employing diverse scaffolds due to their high pqsA enzyme binding affinity. The results of the pharmacophore-based virtual screening were subsequently verified using a molecular dynamic simulation-based study (MDS). Using MDS and post-MDS, the stability of the complexes was evaluated. The most promising lead compounds exhibited a high binding affinity towards protein-binding pocket and interacted with the catalytic dyad. At least one of the scaffolds selected will possibly prove useful for future research. However, further scientific confirmation in the form of preclinical and clinical research is required before implementation.
Collapse
Affiliation(s)
- Muhammad Shahab
- State Key Laboratories of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| | - Muhammad Danial
- Shenzhen Institute of Advanced Technology, University of Chinese Academy of Sciences, Shenzhen 518055, China
| | - Taimur Khan
- State Key Laboratories of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| | - Chaoqun Liang
- State Key Laboratories of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| | - Xiuyuan Duan
- State Key Laboratories of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| | - Daixi Wang
- State Key Laboratories of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| | - Hanzi Gao
- State Key Laboratories of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| | - Guojun Zheng
- State Key Laboratories of Chemical Resources Engineering, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
14
|
Harris JA, Liu R, Martins de Oliveira V, Vázquez-Montelongo EA, Henderson JA, Shen J. GPU-Accelerated All-Atom Particle-Mesh Ewald Continuous Constant pH Molecular Dynamics in Amber. J Chem Theory Comput 2022; 18:7510-7527. [PMID: 36377980 PMCID: PMC10130738 DOI: 10.1021/acs.jctc.2c00586] [Citation(s) in RCA: 26] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Constant pH molecular dynamics (MD) simulations sample protonation states on the fly according to the conformational environment and user specified pH conditions; however, the current accuracy is limited due to the use of implicit-solvent models or a hybrid solvent scheme. Here, we report the first GPU-accelerated implementation, parametrization, and validation of the all-atom continuous constant pH MD (CpHMD) method with particle-mesh Ewald (PME) electrostatics in the Amber22 pmemd.cuda engine. The titration parameters for Asp, Glu, His, Cys, and Lys were derived for the CHARMM c22 and Amber ff14sb and ff19sb force fields. We then evaluated the PME-CpHMD method using the asynchronous pH replica-exchange titration simulations with the c22 force field for six benchmark proteins, including BBL, hen egg white lysozyme (HEWL), staphylococcal nuclease (SNase), thioredoxin, ribonuclease A (RNaseA), and human muscle creatine kinase (HMCK). The root-mean-square deviation from the experimental pKa's of Asp, Glu, His, and Cys is 0.76 pH units, and the Pearson's correlation coefficient for the pKa shifts with respect to model values is 0.80. We demonstrated that a finite-size correction or much enlarged simulation box size can remove a systematic error of the calculated pKa's and improve agreement with experiment. Importantly, the simulations captured the relevant biology in several challenging cases, e.g., the titration order of the catalytic dyad Glu35/Asp52 in HEWL and the coupled residues Asp19/Asp21 in SNase, the large pKa upshift of the deeply buried catalytic Asp26 in thioredoxin, and the large pKa downshift of the deeply buried catalytic Cys283 in HMCK. We anticipate that PME-CpHMD will offer proper pH control to improve the accuracies of MD simulations and enable mechanistic studies of proton-coupled dynamical processes that are ubiquitous in biology but remain poorly understood due to the lack of experimental tools and limitation of current MD simulations.
Collapse
Affiliation(s)
- Julie A Harris
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland21201, United States
| | - Ruibin Liu
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland21201, United States
| | - Vinicius Martins de Oliveira
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland21201, United States.,Lilly Biotechnology Center, San Diego, California92121, United States
| | | | - Jack A Henderson
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland21201, United States
| | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland21201, United States
| |
Collapse
|