1
|
Wujieti B, Feng X, Liu E, Li D, Hao M, Zhou L, Cui W. A theoretical study on the activity and selectivity of IDO/TDO inhibitors. Phys Chem Chem Phys 2024; 26:16747-16764. [PMID: 38818624 DOI: 10.1039/d3cp06036e] [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: 06/01/2024]
Abstract
Indoleamine 2,3-dioxygenase 1 (IDO) is a tryptophan (Trp) metabolic enzyme along the kynurenine (NFK) pathway. Under pathological conditions, IDO overexpressed by tumor cells causes depletion of tryptophan and the accumulation of metabolic products, which inhibit the local immune response and form immune escape. Therefore, the suppression of IDO activity is one of the strategies for tumor immunotherapy, and drug design for this target has been the focus of research for more than two decades. Apart from IDO, tryptophan dioxygenase (TDO) of the same family can also catalyze the same biochemical reaction in the human body, but it has different tissue distribution and substrate selectivity from IDO. Based on the principle of drug design with high potency and low cross-reactivity to specific targets, in this subject, the activity and selectivity of IDO and TDO toward small molecular inhibitors were studied from the perspective of thermodynamics and kinetics. The aim was to elucidate the structural requirements for achieving favorable biological activity and selectivity of IDO and TDO inhibitors. Specifically, the interactions of inhibitors from eight families with IDO and TDO were initially investigated through molecular docking and molecular dynamics simulations, and the thermodynamic data for binding of inhibitors were predicted by the molecular mechanics/generalized Born surface area (MM/GBSA) method. Secondly, we explored the free energy landscape of JKloops, the kinetic control element of IDO/TDO, using temperature replica exchange molecular dynamics (T-REMD) simulations and elucidated the connection between the rules of IDO/TDO conformational changes and the inhibitor selectivity mechanism. Furthermore, the binding and dissociation processes of the C1 inhibitor (NLG919) were simulated by the adaptive steering molecular dynamics (ASMD) method, which not only addressed the possible stable, metastable, and transition states for C1 inhibitor-IDO/TDO interactions, but also accurately predicted kinetic data for C1 inhibitor binding and dissociation. In conclusion, we have constructed a complete process from enzyme (IDO/TDO) conformational activation to inhibitor binding/dissociation and used the thermodynamic and kinetic data of each link as clues to verify the control mechanism of IDO/TDO on inhibitor selectivity. This is of great significance for us to understand the design principles of tumor immunotherapy drugs and to avoid drug resistance of immunotherapy drugs.
Collapse
Affiliation(s)
- Baerlike Wujieti
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China.
| | - Xinping Feng
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China.
| | - Erxia Liu
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China.
| | - Deqing Li
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China.
| | - Mingtian Hao
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China.
| | - Luqi Zhou
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China.
| | - Wei Cui
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China.
| |
Collapse
|
2
|
Gentile R, Modric M, Thiele B, Jaeger KE, Kovacic F, Schott-Verdugo S, Gohlke H. Molecular Mechanisms Underlying Medium-Chain Free Fatty Acid-Regulated Activity of the Phospholipase PlaF from Pseudomonas aeruginosa. JACS AU 2024; 4:958-973. [PMID: 38559719 PMCID: PMC10976570 DOI: 10.1021/jacsau.3c00725] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Revised: 01/08/2024] [Accepted: 01/09/2024] [Indexed: 04/04/2024]
Abstract
PlaF is a membrane-bound phospholipase A1 from Pseudomonas aeruginosa that is involved in remodeling membrane glycerophospholipids (GPLs) and modulating virulence-associated signaling and metabolic pathways. Previously, we identified the role of medium-chain free fatty acids (FFAs) in inhibiting PlaF activity and promoting homodimerization, yet the underlying molecular mechanism remained elusive. Here, we used unbiased and biased molecular dynamics simulations and free energy computations to assess how PlaF interacts with FFAs localized in the water milieu surrounding the bilayer or within the bilayer and how these interactions regulate PlaF activity. Medium-chain FFAs localized in the upper bilayer leaflet can stabilize inactive dimeric PlaF, likely through interactions with charged surface residues, as has been experimentally validated. Potential of mean force (PMF) computations indicate that membrane-bound FFAs may facilitate the activation of monomeric PlaF by lowering the activation barrier for changing into a tilted, active configuration. We estimated that the coupled equilibria of PlaF monomerization-dimerization and tilting at the physiological concentration of PlaF lead to the majority of PlaF forming inactive dimers when in a cell membrane loaded with decanoic acid (C10). This is in agreement with a suggested in vivo product feedback loop and gas chromatography-mass spectrometry profiling results, indicating that PlaF catalyzes the release of C10 from P. aeruginosa membranes. Additionally, we found that C10 in the water milieu can access the catalytic site of active monomeric PlaF, contributing to the competitive component of C10-mediated PlaF inhibition. Our study provides mechanistic insights into how medium-chain FFAs may regulate the activity of PlaF, a potential bacterial drug target.
Collapse
Affiliation(s)
- Rocco Gentile
- Institute
for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Matea Modric
- Institute
of Molecular Enzyme Technology, Heinrich
Heine University Düsseldorf, Forschungszentrum Jülich
GmbH, 52425 Jülich, Germany
| | - Björn Thiele
- Institute
of Bio- and Geosciences (IBG-2: Plant Sciences and IBG-3: Agrosphere), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Karl-Erich Jaeger
- Institute
of Molecular Enzyme Technology, Heinrich
Heine University Düsseldorf, Forschungszentrum Jülich
GmbH, 52425 Jülich, Germany
- Institute
of Bio- and Geosciences (IBG-1: Biotechnology), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Filip Kovacic
- Institute
of Molecular Enzyme Technology, Heinrich
Heine University Düsseldorf, Forschungszentrum Jülich
GmbH, 52425 Jülich, Germany
| | - Stephan Schott-Verdugo
- Institute
of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Holger Gohlke
- Institute
for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
- Institute
of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| |
Collapse
|
3
|
Zhuang Y, Quirk S, Stover ER, Bureau HR, Allen CR, Hernandez R. Tertiary Plasticity Drives the Efficiency of Enterocin 7B Interactions with Lipid Membranes. J Phys Chem B 2024; 128:2100-2113. [PMID: 38412510 PMCID: PMC10926100 DOI: 10.1021/acs.jpcb.3c08199] [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: 02/29/2024]
Abstract
The ability of antimicrobial peptides to efficiently kill their bacterial targets depends on the efficiency of their binding to the microbial membrane. In the case of enterocins, there is a three-part interaction: initial binding, unpacking of helices on the membrane surface, and permeation of the lipid bilayer. Helical unpacking is driven by disruption of the peptide hydrophobic core when in contact with membranes. Enterocin 7B is a leaderless enterocin antimicrobial peptide produced from Enterococcus faecalis that functions alone, or with its cognate partner enterocin 7A, to efficiently kill a wide variety of Gram-stain positive bacteria. To better characterize the role that tertiary structural plasticity plays in the ability of enterocin 7B to interact with the membranes, a series of arginine single-site mutants were constructed that destabilize the hydrophobic core to varying degrees. A series of experimental measures of structure, stability, and function, including CD spectra, far UV CD melting profiles, minimal inhibitory concentrations analysis, and release kinetics of calcein, show that decreased stabilization of the hydrophobic core is correlated with increased efficiency of a peptide to permeate membranes and in killing bacteria. Finally, using the computational technique of adaptive steered molecular dynamics, we found that the atomistic/energetic landscape of peptide mechanical unfolding leads to free energy differences between the wild type and its mutants, whose trends correlate well with our experiment.
Collapse
Affiliation(s)
- Yi Zhuang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Stephen Quirk
- Kimberly-Clark Corporation, Atlanta, Georgia 30076-2199, United States
| | - Erica R Stover
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hailey R Bureau
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Caley R Allen
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Rigoberto Hernandez
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Department of Chemical & Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Department of Materials Science & Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
4
|
Lemcke S, Appeldorn JH, Wand M, Speck T. Toward a structural identification of metastable molecular conformations. J Chem Phys 2023; 159:114105. [PMID: 37712784 DOI: 10.1063/5.0164145] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Accepted: 08/21/2023] [Indexed: 09/16/2023] Open
Abstract
Interpreting high-dimensional data from molecular dynamics simulations is a persistent challenge. In this paper, we show that for a small peptide, deca-alanine, metastable states can be identified through a neural net based on structural information alone. While processing molecular dynamics data, dimensionality reduction is a necessary step that projects high-dimensional data onto a low-dimensional representation that, ideally, captures the conformational changes in the underlying data. Conventional methods make use of the temporal information contained in trajectories generated through integrating the equations of motion, which forgoes more efficient sampling schemes. We demonstrate that EncoderMap, an autoencoder architecture with an additional distance metric, can find a suitable low-dimensional representation to identify long-lived molecular conformations using exclusively structural information. For deca-alanine, which exhibits several helix-forming pathways, we show that this approach allows us to combine simulations with different biasing forces and yields representations comparable in quality to other established methods. Our results contribute to computational strategies for the rapid automatic exploration of the configuration space of peptides and proteins.
Collapse
Affiliation(s)
- Simon Lemcke
- Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudingerweg 7-9, 55128 Mainz, Germany
| | - Jörn H Appeldorn
- Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudingerweg 7-9, 55128 Mainz, Germany
| | - Michael Wand
- Institut für Informatik, Johannes Gutenberg-Universität Mainz, Staudingerweg 9, 55128 Mainz, Germany
| | - Thomas Speck
- Institut für Theoretische Physik IV, Universität Stuttgart, Heisenbergstr. 3, 70569 Stuttgart, Germany
| |
Collapse
|
5
|
Liu X, Liang L, Wu B, Zhang X, Zeng X, Deng Y, Peng B, Zhang X, Zheng L. Effect of the R126C mutation on the structure and function of the glucose transporter GLUT1: A molecular dynamics simulation study. J Mol Graph Model 2022; 116:108227. [PMID: 35671570 DOI: 10.1016/j.jmgm.2022.108227] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2021] [Revised: 05/16/2022] [Accepted: 05/17/2022] [Indexed: 12/14/2022]
Abstract
Glucose transporter 1 (GLUT1) is responsible for basal glucose uptake and is expressed in most tissues under normal conditions. GLUT1 mutations can cause early-onset absence epilepsy and myoclonus dystonia syndrome (MDS), with MDS potentially lethal. In this study, the effect of the R126C mutation, which is associated with MDS, on structural stability and substrate transport of GLUT1 was investigated. Various bioinformatics tools were used to predict the stability of GLUT1, revealing that the R126C mutation reduces the structural stability of GLUT1. Molecular dynamics (MD) simulations were used to further characterize the effect of the R126C mutation on GLUT1 structural stability. Based on the MD simulations, specific conformational changes and dominant motions of the GLUT1 mutant were characterized by Principal component analysis (PCA). The mutation disrupts hydrogen bonds between substrate-binding residues and glucose, thus likely reducing substrate affinity. The R126C mutation reduces the conformational stability of the protein, and fewer intramolecular hydrogen bonds were present in the mutated GLUT1 when compared with that of wild-type GLUT1. The mutation increased the free energy of glucose transport through GLUT1 significantly, especially at the mutation site, indicating that passage of glucose through the channel is hindered, and this mutant may even release cytoplasmic glucose. This study provides a detailed atomic-level explanation for the reduced structural stability and substrate transport capacity of a GLUT1 mutant. The results aid our understanding of the structure of GLUT1 and provide a framework for developing drugs to treat GLUT1-related diseases, such as MDS.
Collapse
Affiliation(s)
- Xiaoliu Liu
- Department of Laboratory Medicine, Nanfang Hospital, Southern Medical University, Guangzhou, 510515, China; Medical Laboratory of Shenzhen Luohu People's Hospital, 518001, China
| | - Luguang Liang
- School of Laboratory Medicine, Guangdong Medical University, Dongguan, China
| | - Bodeng Wu
- Department of Laboratory Medicine, Nanfang Hospital, Southern Medical University, Guangzhou, 510515, China
| | - Xin Zhang
- Department of Laboratory Medicine, Nanfang Hospital, Southern Medical University, Guangzhou, 510515, China
| | | | - Yurong Deng
- Medical Laboratory of Shenzhen Luohu People's Hospital, 518001, China
| | - Bin Peng
- Department of Laboratory Medicine, Nanfang Hospital, Southern Medical University, Guangzhou, 510515, China
| | - Xiuming Zhang
- Medical Laboratory of Shenzhen Luohu People's Hospital, 518001, China.
| | - Lei Zheng
- Department of Laboratory Medicine, Nanfang Hospital, Southern Medical University, Guangzhou, 510515, China.
| |
Collapse
|
6
|
Zhuang Y, Thota N, Quirk S, Hernandez R. Implementation of Telescoping Boxes in Adaptive Steered Molecular Dynamics. J Chem Theory Comput 2022; 18:4649-4659. [PMID: 35830368 DOI: 10.1021/acs.jctc.2c00498] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Long-time dynamical processes, such as those involving protein unfolding and ligand interactions, can be accelerated and realized through steered molecular dynamics (SMD). The challenge has been the extraction of information from such simulations that generalize for complex nonequilibrium processes. The use of Jarzynski's equality opened the possibility of determining the free energy along the steered coordinate, but sampling over the nonequilibrium trajectories is slow to converge. Adaptive steered molecular dynamics (ASMD) and other related techniques have been introduced to overcome this challenge through the use of stages. Here, we take advantage of these stages to address the numerical cost that arises from the required use of very large solvent boxes. We introduce telescoping box schemes within adaptive steered molecular dynamics (ASMD) in which we adjust the solvent box between stages and thereby vary (and optimize) the required number of solvent molecules. We have benchmarked the method on a relatively long α-helical peptide, Ala30, with respect to the potential of mean force and hydrogen bonds. We show that the use of telescoping boxes introduces little numerical error while significantly reducing the computational cost.
Collapse
Affiliation(s)
- Yi Zhuang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Nikhil Thota
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Stephen Quirk
- Kimberly-Clark Corporation, Atlanta, Georgia 30076-2199, United States
| | - Rigoberto Hernandez
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States.,Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States.,Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
7
|
Allen C, Bureau HR, McGee TD, Quirk S, Hernandez R. Benchmarking Adaptive Steered Molecular Dynamics (ASMD) on CHARMM Force Fields. Chemphyschem 2022; 23:e202200175. [PMID: 35594194 PMCID: PMC9543079 DOI: 10.1002/cphc.202200175] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Revised: 05/19/2022] [Indexed: 11/07/2022]
Abstract
The potentials of mean force (PMFs) along the end‐to‐end distance of two different helical peptides have been obtained and benchmarked using the adaptive steered molecular dynamics (ASMD) method. The results depend strongly on the choice of force field driving the underlying all‐atom molecular dynamics, and are reported with respect to the three most popular CHARMM force field versions: c22, c27 and c36. Two small peptides, ALA10
and 1PEF, serve as the particular case studies. The comparisons between the versions of the CHARMM force fields provides both a qualitative and quantitative look at their performance in forced unfolding simulations in which peptides undergo large changes in structural conformations. We find that ASMD with the underlying c36 force field provides the most robust results for the selected benchmark peptides.
Collapse
Affiliation(s)
- Caley Allen
- Johns Hopkins University Zanvyl Krieger School of Arts and Sciences, Chemistry, UNITED STATES
| | - Hailey R Bureau
- Johns Hopkins University Zanvyl Krieger School of Arts and Sciences, Department of Chemistry, UNITED STATES
| | | | | | - Rigoberto Hernandez
- Johns Hopkins University - Homewood Campus: Johns Hopkins University, Chemistry, 3400 N Charles St, 21218, Baltimore, UNITED STATES
| |
Collapse
|
8
|
Sun Z, He Q. Seeding the multi-dimensional nonequilibrium pulling for Hamiltonian variation: indirect nonequilibrium free energy simulations at QM levels. Phys Chem Chem Phys 2022; 24:8800-8819. [PMID: 35352744 DOI: 10.1039/d2cp00355d] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The combination of free energy simulations in the alchemical and configurational spaces provides a feasible route to access the thermodynamic profiles under a computationally demanding target Hamiltonian. Normally, due to the significant differences between the computational cost of ab initio quantum mechanics (QM) calculations and those of semi-empirical quantum mechanics (SQM) and molecular mechanics (MM), this indirect method could be used to obtain the QM thermodynamics by combining the SQM or MM results and the SQM-to-QM or MM-to-QM corrections. In our previous work, a multi-dimensional nonequilibrium pulling framework for Hamiltonian variations was introduced based on bidirectional pulling and bidirectional reweighting. The method performs nonequilibrium free energy simulations in the configurational space to obtain the thermodynamic profile along the conformational change pathway under a selected computationally efficient Hamiltonian, and uses the nonequilibrium alchemical method to correct or perturb the thermodynamic profile to that under the target Hamiltonian. The BAR-based method is designed to achieve the best generality and transferability and thus leads to modest (∼20 fold) speedup. In this work, we explore the possibility of further accelerating the nonequilibrium free energy simulation by employing unidirectional pulling and using the selection criterion to obtain the initial configurations used to initiate nonequilibrium trajectories following the idea of adaptive steered molecular dynamics (ASMD). A single initial condition is used to seed the whole multi-dimensional nonequilibrium free energy simulation and the sampling is performed fully in the nonequilibrium ensemble. Introducing very short ps-length equilibrium sampling to grab more initial seeds could also be helpful. The ASMD scheme estimates the free energy difference with the unidirectional exponential average (EXP), but it does not follow exactly the requirements of the EXP estimator. Another deficiency of the seeding simulation is the inherently sequential or serial pulling due to the inter-segment dependency, which triggers some problems in the parallelizability of the simulation. Numerical tests are performed to grasp some insights and guidelines for using this selection-criterion-based ASMD scheme. The presented selection-criterion-based multi-dimensional ASMD scheme follows the same perturbation network of the BAR-based method, and thus could be used in various Hamiltonian-variation cases.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China.
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio-Engineering Co., Ltd, North Century Avenue 333, 315100 Ningbo, China
| |
Collapse
|
9
|
Tomić A, Tomić S. Demystifying DPP III Catalyzed Peptide Hydrolysis—Computational Study of the Complete Catalytic Cycle of Human DPP III Catalyzed Tynorphin Hydrolysis. Int J Mol Sci 2022; 23:ijms23031858. [PMID: 35163780 PMCID: PMC8836397 DOI: 10.3390/ijms23031858] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Revised: 01/28/2022] [Accepted: 02/03/2022] [Indexed: 12/04/2022] Open
Abstract
Dipeptidyl peptides III (DPP III) is a dual-domain zinc exopeptidase that hydrolyzes peptides of varying sequence and size. Despite attempts to elucidate and narrow down the broad substrate-specificity of DPP III, there is no explanation as to why some of them, such as tynorphin (VVYPW), the truncated form of the endogenous heptapeptide spinorphin, are the slow-reacting substrates of DPP III compared to others, such as Leu-enkephalin. Using quantum molecular mechanics calculations followed by various molecular dynamics techniques, we describe for the first time the entire catalytic cycle of human DPP III, providing theoretical insight into the inhibitory mechanism of tynorphin. The chemical step of peptide bond hydrolysis and the substrate binding to the active site of the enzyme and release of the product were described for DPP III in complex with tynorphin and Leu-enkephalin and their products. We found that tynorphin is cleaved by the same reaction mechanism determined for Leu-enkephalin. More importantly, we showed that the product stabilization and regeneration of the enzyme, but not the nucleophilic attack of the catalytic water molecule and inversion at the nitrogen atom of the cleavable peptide bond, correspond to the rate-determining steps of the overall catalytic cycle of the enzyme.
Collapse
|
10
|
Gimenez-Dejoz J, Numata K. Molecular dynamics study of the internalization of cell-penetrating peptides containing unnatural amino acids across membranes. NANOSCALE ADVANCES 2022; 4:397-407. [PMID: 36132688 PMCID: PMC9419563 DOI: 10.1039/d1na00674f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 11/10/2021] [Indexed: 06/16/2023]
Abstract
Peptide-based delivery systems that deliver target molecules into cells have been gaining traction. These systems need cell-penetrating peptides (CPPs), which have the remarkable ability to penetrate into biological membranes and help internalize different cargoes into cells through the cell membranes. The molecular internalization mechanism and structure-function relationships of CPPs are not clear, although the incorporation of nonproteinogenic amino acids such as α-aminoisobutyric acid (Aib) has been reported to increase their helicity, biostability and penetration efficiencies. Here, we used molecular dynamics to study two Aib-containing CPPs, poly(LysAibAla)3 (KAibA) and poly(LysAibGly)3 (KAibG), that previously showed high cell internalization efficiency. KAibA and KAibG displayed the lowest internalization energies among the studied CPPs, showing distinct internalization mechanisms depending on the lipid composition of the model membranes. The presence of Aib residues allows these CPPs to adopt amphipathic folding to efficiently penetrate through the membranes. Elucidating how Aib incorporation affects CPP-membrane binding and interactions is beneficial for the design of CPPs for efficient intracellular delivery.
Collapse
Affiliation(s)
- Joan Gimenez-Dejoz
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science Saitama Japan
| | - Keiji Numata
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science Saitama Japan
- Department of Material Chemistry, Graduate School of Engineering, Kyoto University Kyoto Japan
| |
Collapse
|
11
|
Qin J, Wilson KA, Sarkar S, Heras B, O'Mara ML, Totsika M. Conserved FimH mutations in the global Escherichia coli ST131 multi-drug resistant lineage weaken interdomain interactions and alter adhesin function. Comput Struct Biotechnol J 2022; 20:4532-4541. [PMID: 36090810 PMCID: PMC9428848 DOI: 10.1016/j.csbj.2022.08.040] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Revised: 08/16/2022] [Accepted: 08/16/2022] [Indexed: 11/27/2022] Open
Abstract
The binding of the type 1 fimbrial adhesin FimH to mannosylated receptors is allosterically regulated to enhance the fitness of uropathogenic Escherichia coli (UPEC) during urinary tract infection (UTI). Mutations in the two FimH domains (pilin and lectin) located outside the mannose binding pocket have been shown to influence mannose binding affinity, yet the details of the allostery mechanism are not fully elucidated. Here we characterised different FimH conformational states (termed low-affinity tense and high-affinity relaxed conformations) of natural FimH variants using molecular dynamics (MD) simulation techniques and report key structural dynamics differences between them. The clinically dominant FimH30 variant from the pandemic multidrug resistant E. coli ST131 lineage contains an R166H mutation that weakens FimH interdomain interactions and allows enhanced mannose interactions with pre-existing high-affinity relaxed conformations. When expressed in an isogenic ST131 strain background, FimH30 mediated high human cell adhesion and invasion, and enhanced biofilm formation over other variants. Collectively, our computational and experimental findings support a model of FimH protein allostery that is mediated by shifts in the pre-existing conformational equilibrium of FimH, additional to the sequential step-wise process of structural perturbations transmitted from one site to another within the protein. Importantly, it is the first study to shed light into how natural mutations in a clinically dominant FimH variant influence the protein’s conformational landscape optimising its function for ST131 fitness at intestinal and extraintestinal niches.
Collapse
|
12
|
Ahmad S, Strunk CH, Schott-Verdugo SN, Jaeger KE, Kovacic F, Gohlke H. Substrate Access Mechanism in a Novel Membrane-Bound Phospholipase A of Pseudomonas aeruginosa Concordant with Specificity and Regioselectivity. J Chem Inf Model 2021; 61:5626-5643. [PMID: 34748335 DOI: 10.1021/acs.jcim.1c00973] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
PlaF is a cytoplasmic membrane-bound phospholipase A1 from Pseudomonas aeruginosa that alters the membrane glycerophospholipid (GPL) composition and fosters the virulence of this human pathogen. PlaF activity is regulated by a dimer-to-monomer transition followed by tilting of the monomer in the membrane. However, how substrates reach the active site and how the characteristics of the active site tunnels determine the activity, specificity, and regioselectivity of PlaF for natural GPL substrates have remained elusive. Here, we combined unbiased and biased all-atom molecular dynamics (MD) simulations and configurational free-energy computations to identify access pathways of GPL substrates to the catalytic center of PlaF. Our results map out a distinct tunnel through which substrates access the catalytic center. PlaF variants with bulky tryptophan residues in this tunnel revealed decreased catalysis rates due to tunnel blockage. The MD simulations suggest that GPLs preferably enter the active site with the sn-1 acyl chain first, which agrees with the experimentally demonstrated PLA1 activity of PlaF. We propose that the acyl chain-length specificity of PlaF is determined by the structural features of the access tunnel, which results in favorable free energy of binding of medium-chain GPLs. The suggested egress route conveys fatty acid (FA) products to the dimerization interface and, thus, contributes to understanding the product feedback regulation of PlaF by FA-triggered dimerization. These findings open up opportunities for developing potential PlaF inhibitors, which may act as antibiotics against P. aeruginosa.
Collapse
Affiliation(s)
- Sabahuddin Ahmad
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Christoph Heinrich Strunk
- Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Stephan N Schott-Verdugo
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany.,Centro de Bioinformática y Simulación Molecular (CBSM), Faculty of Engineering, University of Talca, 3460000 Talca, Chile.,John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry) & Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Karl-Erich Jaeger
- Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany.,Institute of Bio- and Geosciences (IBG-1: Biotechnology), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Filip Kovacic
- Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany.,John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry) & Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| |
Collapse
|
13
|
Chowdhury UD, Bhargava BL. Helix-coil transition and conformational deformity in A β42-monomer: a case study using the Zn 2+ cation. J Biomol Struct Dyn 2021; 40:8949-8960. [PMID: 34018465 DOI: 10.1080/07391102.2021.1927190] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The metal ions (like Fe2+, Zn2+, Cu2+) are known to influence the amyloid beta (Aβ) aggregation. In this study, we have examined the conformational and dynamical changes during the coordination of Aβ-monomer with the Zn2+ ion using all-atom molecular dynamics (MD) simulations using explicit solvent models. We have probed the unfolding of the full-length Aβ42 monomer both inclusive and exclusive of the Zn2+ cation, with 1:1 ratio of the peptide and the Zn2+ cation. The inclusion of the Zn2+ cation shows differential intra-peptide interactions which has been probed using various analyses. The Helix - Coil transition of the wild type Aβ42 monomer is studied using the steered molecular dynamics simulations by taking the end-to-end C-α distance across the peptide. This gives an idea of the unequal intra - peptide and peptide - water interactions being found across the length of the Aβ monomer. The transition of an α-helix dominated wild-type (WT) Aβ structure to the unfolded coil structure gives significant evidence of the intra-peptide hydrogen bonding shifts in the presence of the Zn2+ cation. This accounts for the structural and the dynamical variations that take place in the Aβ monomer in the presence of the Zn2+ cation to mimic the conditions/environment at the onset of fibrillation.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Unmesh D Chowdhury
- School of Chemical Sciences, National Institute of Science Education and Research - Bhubaneswar, HBNI, Khurda, Odisha, India
| | - B L Bhargava
- School of Chemical Sciences, National Institute of Science Education and Research - Bhubaneswar, HBNI, Khurda, Odisha, India
| |
Collapse
|
14
|
Zhuang Y, Bureau HR, Lopez C, Bucher R, Quirk S, Hernandez R. Energetics and structure of alanine-rich α-helices via adaptive steered molecular dynamics. Biophys J 2021; 120:2009-2018. [PMID: 33775636 PMCID: PMC8204395 DOI: 10.1016/j.bpj.2021.03.017] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Revised: 03/03/2021] [Accepted: 03/18/2021] [Indexed: 12/12/2022] Open
Abstract
The energetics and hydrogen bonding profiles of the helix-to-coil transition were found to be an additive property and to increase linearly with chain length, respectively, in alanine-rich α-helical peptides. A model system of polyalanine repeats was used to establish this hypothesis for the energetic trends and hydrogen bonding profiles. Numerical measurements of a synthesized polypeptide Ac-Y(AEAAKA)kF-NH2 and a natural α-helical peptide a2N (1-17) provide evidence of the hypothesis's generality. Adaptive steered molecular dynamics was employed to investigate the mechanical unfolding of all of these alanine-rich polypeptides. We found that the helix-to-coil transition is primarily dependent on the breaking of the intramolecular backbone hydrogen bonds and independent of specific side-chain interactions and chain length. The mechanical unfolding of the α-helical peptides results in a turnover mechanism in which a 310-helical structure forms during the unfolding, remaining at a near constant population and thereby maintaining additivity in the free energy. The intermediate partially unfolded structures exhibited polyproline II helical structure as previously seen by others. In summary, we found that the average force required to pull alanine-rich α-helical peptides in between the endpoints-namely the native structure and free coil-is nearly independent of the length or the specific primary structure.
Collapse
Affiliation(s)
- Yi Zhuang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland
| | - Hailey R Bureau
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland
| | - Christine Lopez
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland
| | - Ryan Bucher
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland
| | | | - Rigoberto Hernandez
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland; Departments of Chemical and Biomolecular Engineering, and Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland.
| |
Collapse
|
15
|
Han D, Wang H, Wujieti B, Zhang B, Cui W, Chen BZ. Insight into the drug resistance mechanisms of GS-9669 caused by mutations of HCV NS5B polymerase via molecular simulation. Comput Struct Biotechnol J 2021; 19:2761-2774. [PMID: 34093991 PMCID: PMC8134009 DOI: 10.1016/j.csbj.2021.04.026] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2021] [Revised: 04/08/2021] [Accepted: 04/09/2021] [Indexed: 12/17/2022] Open
Abstract
GS-9669 is a non-nucleos(t)ide inhibitor (NNI) binding to the thumb site II of the Hepatitis C virus (HCV) NS5B polymerase and has advanced into phase II trials. To clarify the drug resistance mechanisms of GS-9669 caused by M423T/I/V, L419M, R422K, and I482L mutations of NS5B polymerase (GT1b) and the receptor-ligand interactions during the binding process, a series of molecular simulation methods including molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed for the wild-type (WT) and six mutant NS5B/GS-9669 complexes. The calculated results indicate that the binding free energies of the mutant systems are less negative than that of the WT system, indicating that these mutations will indeed cause NS5B to produce different degrees of resistance to GS-9669. The mutation-induced drug resistances are mainly caused by the loss of binding affinities of Leu419 and Trp528 with GS-9669 or the formation of multiple solvent bridges. Moreover, the ASMD calculations show that GS-9669 binds to the thumb II sites of the seven NS5B polymerases in distinct pathways without any obvious energy barriers. Although the recognition methods and binding pathways are distinct, the binding processes of GS-9669 with the WT and mutant NS5B polymerases are basically controlled thermodynamically. This study clearly reveals the resistance mechanisms of GS-9669 caused by M423T/I/V, L419M, R422K, and I482L mutations of HCV NS5B polymerase and provides some valuable clues for further optimization and design of novel NS5B inhibitors.
Collapse
Affiliation(s)
- Di Han
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China.,School of Medical Engineering, Xinxiang Medical University, Xinxiang 453003, China
| | - Huiqun Wang
- Department of Medicinal Chemistry, School of Pharmacy, Virginia Commonwealth University, 800 E Leigh Street, Richmond, VA 23298, USA
| | - Baerlike Wujieti
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China
| | - Beibei Zhang
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China
| | - Wei Cui
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China
| | - Bo-Zhen Chen
- School of Chemical Sciences, University of Chinese Academy of Sciences, No. 19A, YuQuan Road, Beijing 100049, China
| |
Collapse
|
16
|
Ladds MJGW, Popova G, Pastor-Fernández A, Kannan S, van Leeuwen IMM, Håkansson M, Walse B, Tholander F, Bhatia R, Verma CS, Lane DP, Laín S. Exploitation of dihydroorotate dehydrogenase (DHODH) and p53 activation as therapeutic targets: A case study in polypharmacology. J Biol Chem 2020; 295:17935-17949. [PMID: 32900849 PMCID: PMC7939445 DOI: 10.1074/jbc.ra119.012056] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2019] [Revised: 08/18/2020] [Indexed: 01/17/2023] Open
Abstract
The tenovins are a frequently studied class of compounds capable of inhibiting sirtuin activity, which is thought to result in increased acetylation and protection of the tumor suppressor p53 from degradation. However, as we and other laboratories have shown previously, certain tenovins are also capable of inhibiting autophagic flux, demonstrating the ability of these compounds to engage with more than one target. In this study, we present two additional mechanisms by which tenovins are able to activate p53 and kill tumor cells in culture. These mechanisms are the inhibition of a key enzyme of the de novo pyrimidine synthesis pathway, dihydroorotate dehydrogenase (DHODH), and the blockage of uridine transport into cells. These findings hold a 3-fold significance: first, we demonstrate that tenovins, and perhaps other compounds that activate p53, may activate p53 by more than one mechanism; second, that work previously conducted with certain tenovins as SirT1 inhibitors should additionally be viewed through the lens of DHODH inhibition as this is a major contributor to the mechanism of action of the most widely used tenovins; and finally, that small changes in the structure of a small molecule can lead to a dramatic change in the target profile of the molecule even when the phenotypic readout remains static.
Collapse
Affiliation(s)
- Marcus J G W Ladds
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden; SciLifeLab, Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden.
| | - Gergana Popova
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
| | - Andrés Pastor-Fernández
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden; SciLifeLab, Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
| | | | | | | | | | - Fredrik Tholander
- Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
| | - Ravi Bhatia
- Division of Hematology and Oncology, O'Neal Comprehensive Cancer Center, University of Alabama, Birmingham, Alabama, USA
| | - Chandra S Verma
- Bioinformatics Institute (BII), A*STAR, Singapore; Department of Biological Sciences, National University of Singapore, Singapore; School of Biological Sciences, Nanyang Technological University, Singapore
| | - David P Lane
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden; SciLifeLab, Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden
| | - Sonia Laín
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden; SciLifeLab, Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden.
| |
Collapse
|
17
|
Moghadam SA, Preto J, Klobukowski M, Tuszynski JA. Testing amino acid-codon affinity hypothesis using molecular docking. Biosystems 2020; 198:104251. [PMID: 32966852 DOI: 10.1016/j.biosystems.2020.104251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2020] [Revised: 09/08/2020] [Accepted: 09/08/2020] [Indexed: 10/23/2022]
Abstract
Genetic code refers to a set of rules that assign trinucleotides called codons to amino acids in the process of protein synthesis. Investigating the genetic code's logic and its evolutionary origin has always been both intriguing and challenging. While the correspondence rules between codons and amino acids in the genetic code are well-known, it is still unclear whether those assignments can be explained based on energetic or/and entropic arguments. As an attempt at deciphering basic thermodynamic rules governing DNA translation, we used molecular docking to investigate the ability of amino acids to bind to their corresponding anticodon compared to other codons. The total number of 1280 direct docking interactions have been performed for each amino acid-codon/anti-codon case to find whether the amino acids have a preference to bind to their cognate anticodons or codons. Based on docking scores which are expected to correlate with binding affinity, no correlation with genetic correspondence rules was observed suggesting a more subtle process, other than direct binding, to explain codon-amino-acid specificity.
Collapse
Affiliation(s)
- S Arbabi Moghadam
- Department of Physics, University of Alberta, Edmonton, AB, T6G 2E1, Canada
| | - J Preto
- INSERM 1052, CNRS 5286, Centre de Recherche en Cancérologie de Lyon, Université Claude Bernard Lyon 1, Lyon, France
| | - M Klobukowski
- Department of Chemistry, University of Alberta, Edmonton, AB, T6G 2G2, Canada
| | - J A Tuszynski
- Department of Physics, University of Alberta, Edmonton, AB, T6G 2E1, Canada; Department of Oncology, University of Alberta, Edmonton, Alberta, T6G 1Z2, Canada; DIMEAS, Politecnico di Torino, Turin, Italy.
| |
Collapse
|
18
|
Zhuang Y, Bureau HR, Quirk S, Hernandez R. Adaptive steered molecular dynamics of biomolecules. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1807542] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Yi Zhuang
- Department of Chemistry, Johns Hopkins University, Baltimore, MD, USA
| | - Hailey R. Bureau
- Department of Chemistry, Johns Hopkins University, Baltimore, MD, USA
| | | | | |
Collapse
|
19
|
Matić S, Kekez I, Tomin M, Bogár F, Šupljika F, Kazazić S, Hanić M, Jha S, Brkić H, Bourgeois B, Madl T, Gruber K, Macheroux P, Matković-Čalogović D, Matovina M, Tomić S. Binding of dipeptidyl peptidase III to the oxidative stress cell sensor Kelch-like ECH-associated protein 1 is a two-step process. J Biomol Struct Dyn 2020; 39:6870-6881. [PMID: 32811353 DOI: 10.1080/07391102.2020.1804455] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
This work is about synergy of theory and experiment in revealing mechanism of binding of dipeptidyl peptidase III (DPP III) and Kelch-like ECH-associated protein 1 (KEAP1), the main cellular sensor of oxidative stress. The NRF2 ̶ KEAP1 signaling pathway is important for cell protection, but it is also impaired in many cancer cells where NRF2 target gene expression leads to resistance to chemotherapeutic drugs. DPP III competitively binds to KEAP1 in the conditions of oxidative stress and induces release of NRF2 and its translocation into nucleus. The binding is established mainly through the ETGE motif of DPP III and the Kelch domain of KEAP1. However, although part of a flexible loop, ETGE itself is firmly attached to the DPP III surface by strong hydrogen bonds. Using combined computational and experimental study, we found that DPP III ̶ Kelch binding is a two-step process comprising the endergonic loop detachment and exergonic DPP III ̶ Kelch interaction. Substitution of arginines, which keep the ETGE motif attached, decreases the work needed for its release and increases DPP III ̶ Kelch binding affinity. Interestingly, mutations of one of these arginine residues have been reported in cBioPortal for cancer genomics, implicating its possible involvement in cancer development. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sara Matić
- Divison of Organic Chemistry and Biochemistry, Ruđer Bošković Institute, Zagreb, Croatia
| | - Ivana Kekez
- Department of Chemistry, Faculty of Science, University of Zagreb, Zagreb, Croatia
| | - Marko Tomin
- Divison of Organic Chemistry and Biochemistry, Ruđer Bošković Institute, Zagreb, Croatia
| | - Ferenc Bogár
- Department of Medical Chemistry, MTA-SZTE Biomimetic Systems Research Group, University of Szeged, Szeged, Hungary
| | - Filip Šupljika
- Faculty of Food Technology and Biotechnology, University of Zagreb, Zagreb, Croatia
| | - Saša Kazazić
- Divison of Physical Chemistry, Ruđer Bošković Institute, Zagreb, Croatia
| | - Maja Hanić
- Divison of Organic Chemistry and Biochemistry, Ruđer Bošković Institute, Zagreb, Croatia
| | - Shalinee Jha
- Institute of Biochemistry, Graz University of Technology, Graz, Austria
| | - Hrvoje Brkić
- Faculty of Medicine, Josip Juraj Strossmayer University of Osijek, Osijek, Croatia
| | - Benjamin Bourgeois
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging, Molecular Biology & BiocBioTechMed-Graz, Graz, Austria.,BioTechMed-Graz, Graz, Austria
| | - Tobias Madl
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging, Molecular Biology & BiocBioTechMed-Graz, Graz, Austria.,BioTechMed-Graz, Graz, Austria
| | - Karl Gruber
- Institute of Molecular Biosciences, University of Graz, Graz, Austria
| | - Peter Macheroux
- Institute of Biochemistry, Graz University of Technology, Graz, Austria
| | | | - Mihaela Matovina
- Divison of Organic Chemistry and Biochemistry, Ruđer Bošković Institute, Zagreb, Croatia
| | - Sanja Tomić
- Divison of Organic Chemistry and Biochemistry, Ruđer Bošković Institute, Zagreb, Croatia
| |
Collapse
|
20
|
Moesgaard L, Reinholdt P, Wüstner D, Kongsted J. Modeling the Sterol-Binding Domain of Aster-A Provides Insight into Its Multiligand Specificity. J Chem Inf Model 2020; 60:2268-2281. [PMID: 32233488 DOI: 10.1021/acs.jcim.0c00086] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Intracellular transport of cholesterol and related sterols relies to a large degree on nonvesicular mechanisms, which are only partly understood at the molecular level. Aster proteins belonging to the Lam family of sterol transfer proteins have recently been identified as important catalysts of nonvesicular sterol exchange between the plasma membrane (PM) and endoplasmic reticulum (ER). Here, we used a range of computational tools to study the molecular mechanisms underlying sterol binding as well as multisterol ligand specificity of Aster-A. Our study focused primarily on gaining atomistic insight into the bound ligand-protein complex and was, on this basis, performed in the absence of any membrane. Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations provide a rationale for the experimentally found ranking of binding affinities of various sterols to Aster-A. In particular, the polarity of the sterols and the length of their alkyl chain could be identified as being critical determinants of ligand affinity. A Gibbs free energy decomposition identified a charged residue, Glu444, at the base of the binding pose as an important control point for sterol binding. Removing its net charge via protonation was found to cause significant changes to the environment surrounding this residue. In addition, the protonation of Glu444 was found to be paralleled by a large redistribution of molecular flexibility in the Aster domain. This finding was supplemented by multiple branched adaptive steered molecular dynamics (MB-ASMD) simulations by which we defined a possible molecular path for sterol release and demonstrated the importance of Glu444 in this process.
Collapse
Affiliation(s)
- Laust Moesgaard
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Daniel Wüstner
- Department of Biochemistry and Molecular Biology, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| |
Collapse
|
21
|
Bureau HR, Quirk S, Hernandez R. The relative stability of trpzip1 and its mutants determined by computation and experiment. RSC Adv 2020; 10:6520-6535. [PMID: 35495997 PMCID: PMC9049704 DOI: 10.1039/d0ra00920b] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2019] [Accepted: 02/04/2020] [Indexed: 11/21/2022] Open
Abstract
The single-point mutations of tprzip1 are indicated at left, and their relative energetics are compared at right.
Collapse
|
22
|
Zhang Y, Zheng QC. What are the effects of the serine triad on proton conduction of an influenza B M2 channel? An investigation by molecular dynamics simulations. Phys Chem Chem Phys 2019; 21:8820-8826. [PMID: 30968902 DOI: 10.1039/c9cp00612e] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The tetrameric influenza B M2 channel (BM2), an acid activated proton channel, is important in the influenza virus B lifecycle. A conserved HxxxW motif is responsible for proton conduction and channel gating. In this study, to explore the effects of the serine triad (S9, S12 and S16) on proton conduction, we performed classical molecular dynamics (CMD) simulations and adaptive steered molecular dynamics (ASMD) simulations at different protonation states of the H19 tetrad. The results of the pore radius and the C-terminal tilt angle show that the electrostatic repulsion induced by protonated H19 is the key driving force for opening the BM2 channel. The open states could be stabilized by the hydrogen bonds between S16 and protonated H19. The solvent accessible surface area and water density indicate that the polar hydrophilic environment provided by the serine triad facilitates the formation of a water wire, and then exhibits favourable effects on proton conduction. The mutant research verifies and supports these views. Our work clarifies the effects of the serine triad on proton conduction in the BM2 channel, which would help us deeply understand the proton conduction mechanism in BM2 and provides a new perspective for antiviral drug design against BM2.
Collapse
Affiliation(s)
- Yue Zhang
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun 130023, People's Republic of China.
| | | |
Collapse
|
23
|
Wang X, Tu X, Deng B, Zhang JZH, Sun Z. BAR-based optimum adaptive steered MD for configurational sampling. J Comput Chem 2019; 40:1270-1289. [PMID: 30762879 DOI: 10.1002/jcc.25784] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2018] [Revised: 11/05/2018] [Accepted: 01/06/2019] [Indexed: 11/08/2022]
Abstract
The equilibrium and nonequilibrium adaptive alchemical free energy simulation methods optimum Bennett's acceptance ratio and optimum crooks' equation (OCE), based on the statistically optimal bidirectional reweighting estimator named Bennett's Acceptance Ratio or Crooks' equation, perform initial sampling in the staging alchemical transformation and then determine the importance rank of different states via the time-derivative of the variance. The method is proven to give speedups compared with the equal time rule. In the current work, we extend the time derivative of variance guided adaptive sampling method to the configurational space, falling in the term of steered MD (SMD). The SMD approach biasing physically meaningful collective variable (CV) such as one dihedral or one distance to pulling the system from one conformational state to another. By minimizing the variance of the free energy differences along the pathway in an optimized way, a new type of adaptive SMD (ASMD) is introduced. As exhibits in the alchemical case, this adaptive sampling method outperforms the traditional equal-time SMD in nonequilibrium stratification. Also, the method gives much more efficient calculation of potential of mean force than the selection criterion-based ASMD scheme, which is proven to be more efficient than traditional SMD. The OCE workflow is periodicity-of-CV dependent while ASMD is not. The performance is demonstrated in a dihedral flipping case and two distance pulling cases, accounting for periodic and nonperiodic CVs, respectively. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Xiaohui Wang
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Institute of Computational Science, Università della Svizzera italiana (USI), CH-6900, Lugano, Ticino, Switzerland
| | - Xingzhao Tu
- Institute of Organic Chemistry, University of Leipzig, Leipzig 04103, Germany
| | - Boming Deng
- Laboratory of Oil Analysis, Beijing Hangfengkewei Equipment Technology Co., Ltd., Beijing 100141, China
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University, New York, New York, 10003
| | - Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200062, China.,Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Juelich, Jülich 52425, Germany
| |
Collapse
|
24
|
Quirk S, Hopkins MM, Bureau H, Lusk RJ, Allen C, Hernandez R, Bain DL. Mutational Analysis of Neuropeptide Y Reveals Unusual Thermal Stability Linked to Higher-Order Self-Association. ACS OMEGA 2018; 3:2141-2154. [PMID: 29619413 PMCID: PMC5876621 DOI: 10.1021/acsomega.7b01949] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Accepted: 02/08/2018] [Indexed: 06/08/2023]
Abstract
Neuropeptide Y (NPY) is a 36-residue peptide, abundant in the central and peripheral nervous system. The peptide interacts with membrane-bound receptors to control processes such as food intake, vasoconstriction, and memory retention. The N-terminal polyproline sequence of NPY folds back onto a C-terminal α-helix to form a hairpin structure. The hairpin undergoes transient unfolding to allow the monomer to interact with its target membranes and receptors and to form reversible dimers in solution. Using computational, functional, and biophysical approaches, we characterized the role of two conserved tyrosines (Y20 and Y27) located within the hydrophobic core of the hairpin fold. Successive mutation of the tyrosines to more hydrophobic phenylalanines increased the thermal stability of NPY and reduced functional activity, consistent with computational studies predicting a more stable hairpin structure. However, mutant stability was high relative to wild-type: melting temperatures increased by approximately 20 °C for the single mutants (Y20F and Y27F) and by 30 °C for the double mutant (Y20F + Y27F). These findings suggested that the mutations were not just simply enhancing hairpin structure stability, but might also be driving self-association to dimer. Using analytical ultracentrifugation, we determined that the mutations indeed increased self-association, but shifted the equilibrium toward hexamer-like species. Notably, these latter species were not unique to the NPY mutants, but were found to preexist at low levels in the wild-type population. Collectively, the findings indicate that NPY self-association is more complex than previously recognized and that the ensemble of NPY quaternary states is tunable by modulating hairpin hydrophobicity.
Collapse
Affiliation(s)
- Stephen Quirk
- Archeus
Bioscience, 7094 Peachtree
Industrial Blvd., Norcross, Georgia 30071, United
States
| | - Mandi M. Hopkins
- Department
of Pharmaceutical Sciences, University of
Colorado Anschutz Medical Campus, 12850 E Montview Blvd., Aurora, Colorado 80045, United
States
| | - Hailey Bureau
- Center
for Computational and Molecular Science and Technology, School of
Chemistry and Biochemistry, Georgia Institute
of Technology, 901 Atlantic
Dr, Atlanta, Georgia 30332, United States
| | - Ryan J. Lusk
- Department
of Pharmaceutical Sciences, University of
Colorado Anschutz Medical Campus, 12850 E Montview Blvd., Aurora, Colorado 80045, United
States
| | - Caley Allen
- Department
of Chemistry, Johns Hopkins University, 3400 N Charles Street, Baltimore, Maryland 21218, United States
| | - Rigoberto Hernandez
- Center
for Computational and Molecular Science and Technology, School of
Chemistry and Biochemistry, Georgia Institute
of Technology, 901 Atlantic
Dr, Atlanta, Georgia 30332, United States
- Department
of Chemistry, Johns Hopkins University, 3400 N Charles Street, Baltimore, Maryland 21218, United States
| | - David L. Bain
- Department
of Pharmaceutical Sciences, University of
Colorado Anschutz Medical Campus, 12850 E Montview Blvd., Aurora, Colorado 80045, United
States
| |
Collapse
|
25
|
Fan JR, Li H, Zhang HX, Zheng QC. Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine. Biopolymers 2018; 109:e23108. [PMID: 29484634 DOI: 10.1002/bip.23108] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2017] [Revised: 02/06/2018] [Accepted: 02/07/2018] [Indexed: 12/25/2022]
Abstract
The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure-function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.
Collapse
Affiliation(s)
- Jing-Rong Fan
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| | - Heng Li
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, College of Life Science, Jilin University, Changchun, 130012, People's Republic of China
| | - Hong-Xing Zhang
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, 130023, People's Republic of China
| | - Qing-Chuan Zheng
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, 130023, People's Republic of China.,Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, College of Life Science, Jilin University, Changchun, 130012, People's Republic of China
| |
Collapse
|
26
|
Sun DR, Zheng QC, Zhang HX. Molecular dynamics investigation of stereoselective inhibition mechanism of HIF-2α/ARNT heterodimer. J Mol Recognit 2017; 31. [PMID: 28990233 DOI: 10.1002/jmr.2675] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2017] [Revised: 08/27/2017] [Accepted: 09/05/2017] [Indexed: 12/19/2022]
Abstract
Hypoxia-inducible factors (HIFs) are heterodimeric transcription factors related with the onset and progression of solid tumors. Studies demonstrated a class of tetrazole containing chiral inhibitors could stereoselectively disrupt the HIF-2 dimerization and reduce the target gene expression. However, the dynamical features and structural motifs of the HIF-2 heterodimer caused by the binding of enantiomers have not been rationalized at the atomistic level. In this work, molecular dynamics (MD) simulations combined with adaptive steered MD (ASMD) simulations were used to investigate stereoselective interrupting mechanism of HIF-2. Our results decipher that the binding of ligand A (S, R)-24 begets the significant conformation changes of β-sheets and interrupts the HIF-2α/ARNT heterodimerization, which may be attributed to the disruption of the hydrogen bond and salt bridge interactions formed by the 4 foremost residues (Asp240, Arg247, Glu362, and Arg366) and the destruction of hydrophobic interactions on the binding interface. By contrast, the binding of ligand B (R, S)-24 does not disrupt protein dimerization and causes the motion of Fα helix in HIF-2α PAS-B domain to further change the major tunnel for ligand ingress and engress. The present work provides important molecular-level insight into the effect of the binding enantiomers on HIF-2 heterodimerization and bridges the gap between theory and the experimental results, which may conduce to develop highly potent antagonists for intervening the HIF-2-driven tumors.
Collapse
Affiliation(s)
- Dong-Ru Sun
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, China
| | - Qing-Chuan Zheng
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, Jilin University, Changchun, China.,Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, China
| | - Hong-Xing Zhang
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, International Joint Research Laboratory of Nano-Micro Architecture Chemistry, Jilin University, Changchun, China
| |
Collapse
|
27
|
Liu H, Chen F, Sun H, Li D, Hou T. Improving the Efficiency of Non-equilibrium Sampling in the Aqueous Environment via Implicit-Solvent Simulations. J Chem Theory Comput 2017; 13:1827-1836. [DOI: 10.1021/acs.jctc.6b01139] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Hui Liu
- College
of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Fu Chen
- College
of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Huiyong Sun
- College
of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Dan Li
- College
of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Tingjun Hou
- College
of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
- State Key Lab of CAD & CG, Zhejiang University, Hangzhou, Zhejiang 310058, China
| |
Collapse
|
28
|
Thurston BA, Tovar JD, Ferguson AL. Thermodynamics, morphology, and kinetics of early-stage self-assembly of π-conjugated oligopeptides. MOLECULAR SIMULATION 2016. [DOI: 10.1080/08927022.2015.1125997] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Bryce A. Thurston
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - John D. Tovar
- Department of Chemistry, Johns Hopkins University, Baltimore, MD, USA
- Institute of NanoBioTechnology, Johns Hopkins University, Baltimore, MD, USA
- Department of Materials Science and Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Andrew L. Ferguson
- Department of Materials Science and Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| |
Collapse
|
29
|
Bureau HR, Hershkovits E, Quirk S, Hernandez R. Determining the Energetics of Small β-Sheet Peptides using Adaptive Steered Molecular Dynamics. J Chem Theory Comput 2016; 12:2028-37. [PMID: 26930270 DOI: 10.1021/acs.jctc.5b01110] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Mechanically driven unfolding is a useful computational tool for extracting the energetics and stretching pathway of peptides. In this work, two representative β-hairpin peptides, chignolin (PDB: 1UAO ) and trpzip1 (PDB: 1LE0 ), were investigated using an adaptive variant of the original steered molecular dynamics method called adaptive steered molecular dynamics (ASMD). The ASMD method makes it possible to perform energetic calculations on increasingly complex biological systems. Although the two peptides are similar in length and have similar secondary structures, their unfolding energetics are quite different. The hydrogen bonding profile and specific residue pair interaction energies provide insight into the differing stabilities of these peptides and reveal which of the pairs provides the most significant stabilization.
Collapse
Affiliation(s)
- Hailey R Bureau
- Center for Computational and Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Eli Hershkovits
- Center for Computational and Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Stephen Quirk
- Kimberly-Clark Corporation , Atlanta, Georgia 30076-2199, United States
| | - Rigoberto Hernandez
- Center for Computational and Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| |
Collapse
|