1
|
Shawky E, Nassra RA, El-Alkamy AMT, Sallam SM, El Sohafy SM. Unraveling the mechanisms of Fenugreek seed for managing different gynecological disorders: steroidal saponins and isoflavones revealed as key bioactive metabolites. J Pharm Biomed Anal 2024; 238:115865. [PMID: 38000191 DOI: 10.1016/j.jpba.2023.115865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Revised: 11/05/2023] [Accepted: 11/16/2023] [Indexed: 11/26/2023]
Abstract
This study employed network pharmacology-based analysis and reverse molecular docking to investigate the molecular targets and pathways associated with gynecological disorders, particularly those related to steroidal hormones and their receptors, and the potential therapeutic effects of fenugreek (Trigonella foenum-graecum L.) constituents. The STITCH 5.0 database was utilized to identify potential molecular targets, and a compound-target network was constructed. The main targets associated with gynecological disorders included estrogen receptor beta (ESR2), estrogen-related receptor gamma (GPER1), oxytocin receptor (OXTR), progesterone receptor (PGR), prolactin receptor (PRLR), and several enzymes involved in sex hormone biosynthesis. Additionally, network topological analysis revealed that specific compounds, such as quercetin, luteolin, genistein, and vitexin, had significant interactions with the identified targets. Reverse molecular docking analysis confirmed the interactions between the identified compounds and target proteins where quercetin, luteolin, genistein, 4'-methylgenistein, trigoneoside IIB, diosgenin, and vitexin possessed the highest combined docking scores, indicating their multi-target nature. The results highlighted the potential of steroidal saponins, isoflavones, and flavones as active constituents of fenugreek with implications for lactation, reproductive processes, and estrogenic activity. The chemical profiling of saponin-enriched and flavonoid-enriched fractions using UPLC/MS/MS further supported the presence of these bioactive compounds. In an animal model study, the steroidal saponins-enriched fraction of fenugreek seed exhibited a significant increase in the body weight of lactating female rats and serum prolactin levels while the flavonoids-enriched fraction showed an increase in serum estradiol levels and improved the histological structure of ovaries.
Collapse
Affiliation(s)
- Eman Shawky
- Department of Pharmacognosy, Faculty of Pharmacy, Alexandria University, Egypt.
| | - Rasha A Nassra
- Medical Biochemistry department, faculty of medicine, Alexandria University, Egypt
| | - Aliaa M T El-Alkamy
- Human Anatomy and Embryology Department, faculty of medicine, Alexandria University, Egypt
| | - Shaimaa M Sallam
- Department of Pharmacognosy, Faculty of Pharmacy, Alexandria University, Egypt
| | - Samah M El Sohafy
- Department of Pharmacognosy, Faculty of Pharmacy, Alexandria University, Egypt
| |
Collapse
|
2
|
The pH Effects on SARS-CoV and SARS-CoV-2 Spike Proteins in the Process of Binding to hACE2. Pathogens 2022; 11:pathogens11020238. [PMID: 35215181 PMCID: PMC8879864 DOI: 10.3390/pathogens11020238] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2021] [Revised: 01/28/2022] [Accepted: 02/07/2022] [Indexed: 02/01/2023] Open
Abstract
COVID-19 has been threatening human health since the late 2019, and has a significant impact on human health and economy. Understanding SARS-CoV-2 and other coronaviruses is important to develop effective treatments for COVID-19 and other coronavirus-caused diseases. In this work, we applied multi-scale computational approaches to study the electrostatic features of spike (S) proteins for SARS-CoV and SARS-CoV-2. From our results, we found that SARS-CoV and SARS-CoV-2 have similar charge distributions and electrostatic features when binding with the human angiotensin-converting enzyme 2 (hACE2). Energy pH-dependence calculations revealed that the complex structures of hACE2 and the S proteins of SARS-CoV/SARS-CoV-2 are stable at pH values ranging from 7.5 to 9. Three independent 100 ns molecular dynamics (MD) simulations were performed using NAMD to investigate the hydrogen bonds between S proteins RBD and hACE2 RBD. From MD simulations, we found that SARS-CoV-2 forms 19 pairs (average of three simulations) of hydrogen bonds with high occupancy (>50%) to hACE2, compared to 16 pairs between SARS-CoV and hACE2. Additionally, SARS-CoV viruses prefer sticking to the same hydrogen bond pairs, while SARS-CoV-2 tends to have a larger range of selections on hydrogen bonds acceptors. We also labelled key residues involved in forming the top five hydrogen bonds that were found in all three independent 100 ns simulations. This identification is important to potential drug designs for COVID-19 treatments. Our work will shed the light on current and future coronavirus-caused diseases.
Collapse
|
3
|
Guo W, Sun S, Sanchez JE, Lopez-Hernandez AE, Ale TA, Chen J, Afrin T, Qiu W, Xie Y, Li L. Using a comprehensive approach to investigate the interaction between Kinesin-5/Eg5 and the microtubule. Comput Struct Biotechnol J 2022; 20:4305-4314. [PMID: 36051882 PMCID: PMC9396395 DOI: 10.1016/j.csbj.2022.08.020] [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: 04/12/2022] [Revised: 08/01/2022] [Accepted: 08/08/2022] [Indexed: 01/02/2023] Open
Abstract
Kinesins are microtubule-based motor proteins that play important roles ranging from intracellular transport to cell division. Human Kinesin-5 (Eg5) is essential for mitotic spindle assembly during cell division. By combining molecular dynamics (MD) simulations with other multi-scale computational approaches, we systematically studied the interaction between Eg5 and the microtubule. We find the electrostatic feature on the motor domains of Eg5 provides attractive interactions to the microtubule. Additionally, the folding and binding energy analysis reveals that the Eg5 motor domain performs its functions best when in a weak acidic environment. Molecular dynamics analyses of hydrogen bonds and salt bridges demonstrate that, on the binding interfaces of Eg5 and the tubulin heterodimer, salt bridges play the most significant role in holding the complex. The salt bridge residues on the binding interface of Eg5 are mostly positive, while salt bridge residues on the binding interface of tubulin heterodimer are mostly negative. Such salt bridge residue distribution is consistent with electrostatic potential calculations. In contrast, the interface between α and β-tubulins is dominated by hydrogen bonds rather than salt bridges. Compared to the Eg5/α-tubulin interface, the Eg5/β-tubulin interface has a greater number of salt bridges and higher occupancy for salt bridges. This asymmetric salt bridge distribution may play a significant role in Eg5′s directionality. The residues involved in hydrogen bonds and salt bridges are identified in this work and may be helpful for anticancer drug design.
Collapse
Affiliation(s)
- Wenhan Guo
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | - Shengjie Sun
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | - Jason E. Sanchez
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | | | - Tolulope A. Ale
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | - Jiawei Chen
- Department of Physics, University of Texas at El Paso, El Paso, TX, USA
| | - Tanjina Afrin
- Department of Physics, University of Texas at El Paso, El Paso, TX, USA
- Department of Physics, Oregon State University, Corvallis, OR, USA
| | - Weihong Qiu
- Department of Physics, Oregon State University, Corvallis, OR, USA
| | - Yixin Xie
- Department of Information Technology, Kennesaw State University, Kennesaw, GA, USA
| | - Lin Li
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
- Department of Physics, University of Texas at El Paso, El Paso, TX, USA
- Corresponding author.
| |
Collapse
|
4
|
Xie Y, Guo W, Lopez-Hernadez A, Teng S, Li L. The pH Effects on SARS-CoV and SARS-CoV-2 Spike Proteins in the Process of Binding to hACE2. RESEARCH SQUARE 2021:rs.3.rs-871118. [PMID: 34518836 PMCID: PMC8437318 DOI: 10.21203/rs.3.rs-871118/v1] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
COVID-19 has been threatening human health since the late 2019, which has significant impact on human health and economy. Understanding the SARS-CoV-2 and other coronaviruses is important to develop effective treatments for COVID-19 and other coronaviruses-caused diseases. In this work, we applied multi-scale computational approaches to study the electrostatic features of spike (S) proteins for SARS-CoV and SARS-CoV-2. From our results, we found thatSARS-CoV and SARS-CoV-2 have similar charge distributions and electrostatic features when binding with the human angiotensin-converting enzyme 2 (hACE2). The energy pH-dependence calculation srevealed that the complex structures of hACE2 and the S proteins of SARS-CoV/SARS-CoV-2 are stable at pH values ranging from 7.5 to 9. Molecular dynamics simulations were performed using NAMD to investigate the hydrogen bonds between S proteins and hACE2. From the MD simulations it was found that SARS-CoV-2 has four pairsof essential hydrogenbonds (high occupancy, >80%), while SARS-CoV has three pairs, which indicates the SARS-CoV-2 S protein has relatively more robust binding strategy than SARS-CoVS protein.Four key residues forming essential hydrogen bonds from SARS-CoV-2 are identified, which are potential drug targets for COVID-19 treatments. The findings in this study shed lights on the current and future treatments for COVID-19 and other coronaviruses-caused diseases.
Collapse
Affiliation(s)
- Yixin Xie
- Computational Science Program, University of Texas at El Paso, El Paso, TX
| | - Wenhan Guo
- Computational Science Program, University of Texas at El Paso, El Paso, TX
| | | | - Shaolei Teng
- Department of Biology, Howard University, Washington, D.C
| | - Lin Li
- Computational Science Program, University of Texas at El Paso, El Paso, TX
- Department of Physics, University of Texas at El Paso, El Paso, TX
| |
Collapse
|
5
|
Guo W, Xie Y, Lopez-Hernandez AE, Sun S, Li L. Electrostatic features for nucleocapsid proteins of SARS-CoV and SARS-CoV-2. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2021; 18:2372-2383. [PMID: 33892550 PMCID: PMC8279046 DOI: 10.3934/mbe.2021120] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/24/2023]
Abstract
COVID-19 is increasingly affecting human health and global economy. Understanding the fundamental mechanisms of Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) is highly demanded to develop treatments for COVID-19. SARS-CoV and SARS-CoV-2 share 92.06% identity in their N protein RBDs' sequences, which results in very similar structures. However, the SARS-CoV-2 is more easily to spread. Utilizing multi-scale computational approaches, this work studied the fundamental mechanisms of the nucleocapsid (N) proteins of SARS-CoV and SARS-CoV-2, including their stabilities and binding strengths with RNAs at different pH values. Electrostatic potential on the surfaces of N proteins show that both the N proteins of SARS-CoV and SARS-CoV-2 have dominantly positive potential to attract RNAs. The binding forces between SARS-CoV N protein and RNAs at different distances are similar to that of SARS-CoV-2, both in directions and magnitudes. The electric filed lines between N proteins and RNAs are also similar for both SARS-CoV and SARS-CoV-2. The folding energy and binding energy dependence on pH revealed that the best environment for N proteins to perform their functions with RNAs is the weak acidic environment.
Collapse
Affiliation(s)
- Wenhan Guo
- Computational Science Program, University of Texas at El Paso, El Paso, TX 79968, USA
| | - Yixin Xie
- Computational Science Program, University of Texas at El Paso, El Paso, TX 79968, USA
| | | | - Shengjie Sun
- Computational Science Program, University of Texas at El Paso, El Paso, TX 79968, USA
| | - Lin Li
- Computational Science Program, University of Texas at El Paso, El Paso, TX 79968, USA
- Department of Physics, University of Texas at El Paso, El Paso, TX 79968, USA
| |
Collapse
|
6
|
Lopez-Hernandez AE, Xie Y, Guo W, Li L. The Electrostatic Features of Dengue Virus Capsid Assembly. JOURNAL OF COMPUTATIONAL BIOPHYSICS AND CHEMISTRY 2021. [DOI: 10.1142/s2737416520420089] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Dengue virus causes serious diseases and deaths in the world. Understanding the fundamental mechanisms of dengue virus is highly demanded to develop treatments for dengue virus caused diseases. Here, we present a computational work which focused on the stability of dengue viral capsid. The interactions among E proteins on the dengue viral capsid were studied using several computational approaches. It was found that the electrostatic distribution on a single E protein monomer is highly inhomogeneous, which makes an E protein strongly binding with another E protein. This is the reason why all the E proteins form homodimers as the basic units on the whole dengue viral capsids. The pKa calculations of E proteins demonstrated that the folding energy of an E protein is low and stable in the range of pH 6–10, which is different from many other proteins that are stable at certain pH. The pH dependence of binding energy of E protein homodimer shows that the binding energy is low and independent from pH when the pH is also in the range of 6–10. This finding implies that the dengue virus can survive in a wide range of pH, which can explain why the dengue virus is so widely distributed in the world and spreads fast.
Collapse
Affiliation(s)
| | - Yixin Xie
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | - Wenhan Guo
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | - Lin Li
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
- Department of Physics, University of Texas at El Paso, El Paso, TX, USA
| |
Collapse
|
7
|
Shashikala HBM, Chakravorty A, Panday SK, Alexov E. BION-2: Predicting Positions of Non-Specifically Bound Ions on Protein Surface by a Gaussian-Based Treatment of Electrostatics. Int J Mol Sci 2020; 22:ijms22010272. [PMID: 33383946 PMCID: PMC7794834 DOI: 10.3390/ijms22010272] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Revised: 12/25/2020] [Accepted: 12/26/2020] [Indexed: 01/16/2023] Open
Abstract
Ions play significant roles in biological processes—they may specifically bind to a protein site or bind non-specifically on its surface. Although the role of specifically bound ions ranges from actively providing structural compactness via coordination of charge–charge interactions to numerous enzymatic activities, non-specifically surface-bound ions are also crucial to maintaining a protein’s stability, responding to pH and ion concentration changes, and contributing to other biological processes. However, the experimental determination of the positions of non-specifically bound ions is not trivial, since they may have a low residential time and experience significant thermal fluctuation of their positions. Here, we report a new release of a computational method, the BION-2 method, that predicts the positions of non-specifically surface-bound ions. The BION-2 utilizes the Gaussian-based treatment of ions within the framework of the modified Poisson–Boltzmann equation, which does not require a sharp boundary between the protein and water phase. Thus, the predictions are done by the balance of the energy of interaction between the protein charges and the corresponding ions and the de-solvation penalty of the ions as they approach the protein. The BION-2 is tested against experimentally determined ion’s positions and it is demonstrated that it outperforms the old BION and other available tools.
Collapse
Affiliation(s)
- H. B. Mihiri Shashikala
- Department of Physics, Clemson University, Clemson, SC 29634, USA; (H.B.M.S.); (A.C.); (S.K.P.)
| | - Arghya Chakravorty
- Department of Physics, Clemson University, Clemson, SC 29634, USA; (H.B.M.S.); (A.C.); (S.K.P.)
- Department of Chemistry, University of Michigan, Ann Arbor, MI 48109, USA
| | - Shailesh Kumar Panday
- Department of Physics, Clemson University, Clemson, SC 29634, USA; (H.B.M.S.); (A.C.); (S.K.P.)
| | - Emil Alexov
- Department of Physics, Clemson University, Clemson, SC 29634, USA; (H.B.M.S.); (A.C.); (S.K.P.)
- Correspondence:
| |
Collapse
|
8
|
Baker S, Petukh M. Effect of pH on the Ability of N-Terminal Domain of Human NPC1 to Recognize, Bind, and Transfer Cholesterol. ACS OMEGA 2020; 5:29222-29230. [PMID: 33225153 PMCID: PMC7676335 DOI: 10.1021/acsomega.0c03983] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/18/2020] [Accepted: 10/21/2020] [Indexed: 06/11/2023]
Abstract
Niemann-Pick type C1 (NPC1) is a large multidomain transmembrane protein essential for transporting cholesterol (CLR) from late endosomes and lysosomes to the endoplasmic reticulum and other cellular compartments. The lumen-facing N-terminal domain (NTD), involved in direct binding of CLR, is expected to have an optimum activity at acidic pH = 4.5. Here, we show that acidic pH is vital for the functionality of NPC1(NTD) and should be taken into account when studying the protein activity. We applied evolutionary, structural, and physicochemical analyses to decipher the consequences of a change in pH from acidic (pH = 4.5) to neutral (pH = 7.2) on the structural integrity of the NTD and its ability to bind CLR. We revealed that the change in pH from 4.5 to 7.2 increases the potential energy of the protein in both apo- and holo-states making the system less energetically favorable. At neutral pH, the flexibility of the protein in the apo-state is decreased caused by the alteration of specific interactions, which in turn might have a high impact on ligand recognition and binding. In contrast, neutral pH significantly exaggerates the flexibility of the protein with bound CLR that causes a partial exposure of the ligand to the water phase and its mislocation inside the ligand-binding pocket, which might obstruct CLR translocation through the membrane.
Collapse
|
9
|
Flinck M, Kramer SH, Pedersen SF. Roles of pH in control of cell proliferation. Acta Physiol (Oxf) 2018; 223:e13068. [PMID: 29575508 DOI: 10.1111/apha.13068] [Citation(s) in RCA: 84] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2017] [Revised: 02/17/2018] [Accepted: 03/19/2018] [Indexed: 02/06/2023]
Abstract
Precise spatiotemporal regulation of intracellular pH (pHi ) is a prerequisite for normal cell function, and changes in pHi or pericellular pH (pHe ) exert important signalling functions. It is well established that proliferation of mammalian cells is dependent on a permissive pHi in the slightly alkaline range (7.0-7.2). It is also clear that mitogen signalling in nominal absence of HCO3- is associated with an intracellular alkalinization (~0.3 pH unit above steady-state pHi ), which is secondary to activation of Na+ /H+ exchange. However, it remains controversial whether this increase in pHi is part of the mitogenic signal cascade leading to cell cycle entry and progression, and whether it is relevant under physiological conditions. Furthermore, essentially all studies of pHi in mammalian cell proliferation have focused on the mitogen-induced G0-G1 transition, and the regulation and roles of pHi during the cell cycle remain poorly understood. The aim of this review is to summarize and critically discuss the possible roles of pHi and pHe in cell cycle progression. While the focus is on the mammalian cell cycle, important insights from studies in lower eukaryotes are also discussed. We summarize current evidence of links between cell cycle progression and pHi and discuss possible pHi - and pHe sensors and signalling pathways relevant to mammalian proliferation control. The possibility that changes in pHi during cell cycle progression may be an integral part of the checkpoint control machinery is explored. Finally, we discuss the relevance of links between pH and proliferation in the context of the perturbed pH homoeostasis and acidic microenvironment of solid tumours.
Collapse
Affiliation(s)
- M. Flinck
- Section for Cell Biology and Physiology; Department of Biology; Faculty of Science; University of Copenhagen; Copenhagen Denmark
| | - S. H. Kramer
- Section for Cell Biology and Physiology; Department of Biology; Faculty of Science; University of Copenhagen; Copenhagen Denmark
| | - S. F. Pedersen
- Section for Cell Biology and Physiology; Department of Biology; Faculty of Science; University of Copenhagen; Copenhagen Denmark
| |
Collapse
|
10
|
Harris RC, Tsai CC, Ellis CR, Shen J. Proton-Coupled Conformational Allostery Modulates the Inhibitor Selectivity for β-Secretase. J Phys Chem Lett 2017; 8:4832-4837. [PMID: 28927275 PMCID: PMC5713904 DOI: 10.1021/acs.jpclett.7b02309] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
Many important pharmaceutical targets, such as aspartyl proteases and kinases, exhibit pH-dependent dynamics, functions and inhibition. Accurate prediction of their binding free energies is challenging because current computational techniques neglect the effects of pH. Here we combine free energy perturbation calculations with continuous constant pH molecular dynamics to explore the selectivity of a small-molecule inhibitor for β-secretase (BACE1), an important drug target for Alzheimer's disease. The calculations predicted identical affinity for BACE1 and the closely related cathepsin D at high pH; however, at pH 4.6 the inhibitor is selective for BACE1 by 1.3 kcal/mol, in excellent agreement with experiment. Surprisingly, the pH-dependent selectivity can be attributed to the protonation of His45, which allosterically modulates a loop-inhibitor interaction. Allosteric regulation induced by proton binding is likely common in biology; considering such allosteric sites could lead to exciting new opportunities in drug design.
Collapse
Affiliation(s)
- Robert C Harris
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy , Baltimore, Maryland 21201, United States
| | - Cheng-Chieh Tsai
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy , Baltimore, Maryland 21201, United States
| | - Christopher R Ellis
- 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
|
11
|
Wang L, Zhang M, Alexov E. DelPhiPKa web server: predicting pKa of proteins, RNAs and DNAs. Bioinformatics 2016; 32:614-5. [PMID: 26515825 PMCID: PMC5963359 DOI: 10.1093/bioinformatics/btv607] [Citation(s) in RCA: 80] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2015] [Accepted: 10/14/2015] [Indexed: 12/25/2022] Open
Abstract
UNLABELLED A new pKa prediction web server is released, which implements DelPhi Gaussian dielectric function to calculate electrostatic potentials generated by charges of biomolecules. Topology parameters are extended to include atomic information of nucleotides of RNA and DNA, which extends the capability of pKa calculations beyond proteins. The web server allows the end-user to protonate the biomolecule at particular pH based on calculated pKa values and provides the downloadable file in PQR format. Several tests are performed to benchmark the accuracy and speed of the protocol. IMPLEMENTATION The web server follows a client-server architecture built on PHP and HTML and utilizes DelPhiPKa program. The computation is performed on the Palmetto supercomputer cluster and results/download links are given back to the end-user via http protocol. The web server takes advantage of MPI parallel implementation in DelPhiPKa and can run a single job on up to 24 CPUs. AVAILABILITY AND IMPLEMENTATION The DelPhiPKa web server is available at http://compbio.clemson.edu/pka_webserver.
Collapse
Affiliation(s)
- Lin Wang
- Computational Biophysics and Bioinformatics, Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA and,*To whom correspondence should be addressed
| | - Min Zhang
- Department of Computer Sciences, East Carolina University, Greenville, NC 27858, USA
| | - Emil Alexov
- Computational Biophysics and Bioinformatics, Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA and,*To whom correspondence should be addressed
| |
Collapse
|
12
|
Kim MO, McCammon JA. Computation of pH-dependent binding free energies. Biopolymers 2016; 105:43-9. [PMID: 26202905 PMCID: PMC4623928 DOI: 10.1002/bip.22702] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2015] [Accepted: 07/20/2015] [Indexed: 01/21/2023]
Abstract
Protein-ligand binding accompanies changes in the surrounding electrostatic environments of the two binding partners and may lead to changes in protonation upon binding. In cases where the complex formation results in a net transfer of protons, the binding process is pH-dependent. However, conventional free energy computations or molecular docking protocols typically employ fixed protonation states for the titratable groups in both binding partners set a priori, which are identical for the free and bound states. In this review, we draw attention to these important yet largely ignored binding-induced protonation changes in protein-ligand association by outlining physical origins and prevalence of the protonation changes upon binding. Following a summary of various theoretical methods for pKa prediction, we discuss the theoretical framework to examine the pH dependence of protein-ligand binding processes.
Collapse
Affiliation(s)
- M. Olivia Kim
- Department of Pharmacology, University of California San Diego, La Jolla, CA 92093, USA
| | - J. Andrew McCammon
- Department of Pharmacology, University of California San Diego, La Jolla, CA 92093, USA
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
- Howard Hughes Medical Institute, University of California San Diego, La Jolla, CA 92093, USA
- National Biomedical Computation Resource, University of California San Diego, La Jolla, CA 92093, USA
| |
Collapse
|
13
|
Conformational Dynamics and Binding Free Energies of Inhibitors of BACE-1: From the Perspective of Protonation Equilibria. PLoS Comput Biol 2015; 11:e1004341. [PMID: 26506513 PMCID: PMC4623973 DOI: 10.1371/journal.pcbi.1004341] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2015] [Accepted: 05/17/2015] [Indexed: 11/19/2022] Open
Abstract
BACE-1 is the β-secretase responsible for the initial amyloidogenesis in Alzheimer’s disease, catalyzing hydrolytic cleavage of substrate in a pH-sensitive manner. The catalytic mechanism of BACE-1 requires water-mediated proton transfer from aspartyl dyad to the substrate, as well as structural flexibility in the flap region. Thus, the coupling of protonation and conformational equilibria is essential to a full in silico characterization of BACE-1. In this work, we perform constant pH replica exchange molecular dynamics simulations on both apo BACE-1 and five BACE-1-inhibitor complexes to examine the effect of pH on dynamics and inhibitor binding properties of BACE-1. In our simulations, we find that solution pH controls the conformational flexibility of apo BACE-1, whereas bound inhibitors largely limit the motions of the holo enzyme at all levels of pH. The microscopic pKa values of titratable residues in BACE-1 including its aspartyl dyad are computed and compared between apo and inhibitor-bound states. Changes in protonation between the apo and holo forms suggest a thermodynamic linkage between binding of inhibitors and protons localized at the dyad. Utilizing our recently developed computational protocol applying the binding polynomial formalism to the constant pH molecular dynamics (CpHMD) framework, we are able to obtain the pH-dependent binding free energy profiles for various BACE-1-inhibitor complexes. Our results highlight the importance of correctly addressing the binding-induced protonation changes in protein-ligand systems where binding accompanies a net proton transfer. This work comprises the first application of our CpHMD-based free energy computational method to protein-ligand complexes and illustrates the value of CpHMD as an all-purpose tool for obtaining pH-dependent dynamics and binding free energies of biological systems. Formation of insoluble amyloid plaques in the vascular and hippocampal areas of the brain characterizes Alzheimer’s disease, a devastating neurodegenerative disorder causing dementia. Site-specific hydrolytic catalysis of β-secretase, or BACE-1, is responsible for production of oligomerative amyloid β-peptide. As the catalytic activity of BACE-1 is pH-dependent and its structural dynamics are intrinsic to the catalysis, we examine the dependence of dynamics of BACE-1 on solution pH and its implications on the catalytic mechanism of BACE-1. Also, we highlight the importance of accurate description of protonation states of the titratable groups in computer-aided drug discovery targeting BACE-1. We hope the understanding of pH dependence of the dynamics and inhibitor binding properties of BACE-1 will aid the structure-based inhibitor design efforts against Alzheimer’s disease.
Collapse
|
14
|
Sazanavets I, Warwicker J. Computational Tools for Interpreting Ion Channel pH-Dependence. PLoS One 2015; 10:e0125293. [PMID: 25915903 PMCID: PMC4411139 DOI: 10.1371/journal.pone.0125293] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2014] [Accepted: 03/21/2015] [Indexed: 01/23/2023] Open
Abstract
Activity in many biological systems is mediated by pH, involving proton titratable groups with pKas in the relevant pH range. Experimental analysis of pH-dependence in proteins focusses on particular sidechains, often with mutagenesis of histidine, due to its pKa near to neutral pH. The key question for algorithms that predict pKas is whether they are sufficiently accurate to effectively narrow the search for molecular determinants of pH-dependence. Through analysis of inwardly rectifying potassium (Kir) channels and acid-sensing ion channels (ASICs), mutational effects on pH-dependence are probed, distinguishing between groups described as pH-coupled or pH-sensor. Whereas mutation can lead to a shift in transition pH between open and closed forms for either type of group, only for pH-sensor groups does mutation modulate the amplitude of the transition. It is shown that a hybrid Finite Difference Poisson-Boltzmann (FDPB) – Debye-Hückel continuum electrostatic model can filter mutation candidates, providing enrichment for key pH-coupled and pH-sensor residues in both ASICs and Kir channels, in comparison with application of FDPB alone.
Collapse
Affiliation(s)
- Ivan Sazanavets
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Jim Warwicker
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- * E-mail:
| |
Collapse
|
15
|
Li L, Wang L, Alexov E. On the energy components governing molecular recognition in the framework of continuum approaches. Front Mol Biosci 2015; 2:5. [PMID: 25988173 PMCID: PMC4429657 DOI: 10.3389/fmolb.2015.00005] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Accepted: 02/04/2015] [Indexed: 01/14/2023] Open
Abstract
Molecular recognition is a process that brings together several biological macromolecules to form a complex and one of the most important characteristics of the process is the binding free energy. Various approaches exist to model the binding free energy, provided the knowledge of the 3D structures of bound and unbound molecules. Among them, continuum approaches are quite appealing due to their computational efficiency while at the same time providing predictions with reasonable accuracy. Here we review recent developments in the field emphasizing on the importance of adopting adequate description of physical processes taking place upon the binding. In particular, we focus on the efforts aiming at capturing some of the atomistic details of the binding phenomena into the continuum framework. When possible, the energy components are reviewed independently of each other. However, it is pointed out that rigorous approaches should consider all energy contributions on the same footage. The two major schemes for utilizing the individual energy components to predict binding affinity are outlined as well.
Collapse
Affiliation(s)
- Lin Li
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University Clemson, SC, USA
| | - Lin Wang
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University Clemson, SC, USA
| | - Emil Alexov
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University Clemson, SC, USA
| |
Collapse
|
16
|
Stefl S, Nishi H, Petukh M, Panchenko AR, Alexov E. Molecular mechanisms of disease-causing missense mutations. J Mol Biol 2013; 425:3919-36. [PMID: 23871686 DOI: 10.1016/j.jmb.2013.07.014] [Citation(s) in RCA: 187] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2013] [Revised: 07/04/2013] [Accepted: 07/10/2013] [Indexed: 12/23/2022]
Abstract
Genetic variations resulting in a change of amino acid sequence can have a dramatic effect on stability, hydrogen bond network, conformational dynamics, activity and many other physiologically important properties of proteins. The substitutions of only one residue in a protein sequence, so-called missense mutations, can be related to many pathological conditions and may influence susceptibility to disease and drug treatment. The plausible effects of missense mutations range from affecting the macromolecular stability to perturbing macromolecular interactions and cellular localization. Here we review the individual cases and genome-wide studies that illustrate the association between missense mutations and diseases. In addition, we emphasize that the molecular mechanisms of effects of mutations should be revealed in order to understand the disease origin. Finally, we report the current state-of-the-art methodologies that predict the effects of mutations on protein stability, the hydrogen bond network, pH dependence, conformational dynamics and protein function.
Collapse
Affiliation(s)
- Shannon Stefl
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University, Clemson, SC 29634, USA
| | | | | | | | | |
Collapse
|
17
|
Abstract
Formation of protein-ligand complexes causes various changes in both the receptor and the ligand. This review focuses on changes in pK and protonation states of ionizable groups that accompany protein-ligand binding. Physical origins of these effects are outlined, followed by a brief overview of the computational methods to predict them and the associated corrections to receptor-ligand binding affinities. Statistical prevalence, magnitude and spatial distribution of the pK and protonation state changes in protein-ligand binding are discussed in detail, based on both experimental and theoretical studies. While there is no doubt that these changes occur, they do not occur all the time; the estimated prevalence varies, both between individual complexes and by method. The changes occur not only in the immediate vicinity of the interface but also sometimes far away. When receptor-ligand binding is associated with protonation state change at particular pH, the binding becomes pH dependent: we review the interplay between sub-cellular characteristic pH and optimum pH of receptor-ligand binding. It is pointed out that there is a tendency for protonation state changes upon binding to be minimal at physiologically relevant pH for each complex (no net proton uptake/release), suggesting that native receptor-ligand interactions have evolved to reduce the energy cost associated with ionization changes. As a result, previously reported statistical prevalence of these changes - typically computed at the same pH for all complexes - may be higher than what may be expected at optimum pH specific to each complex. We also discuss whether proper account of protonation state changes appears to improve practical docking and scoring outcomes relevant to structure-based drug design. An overview of some of the existing challenges in the field is provided in conclusion.
Collapse
Affiliation(s)
- Alexey V Onufriev
- Department of Computer Science and Physics, 2050 Torgersen Hall, Virginia Tech, Blacksburg, VA 24061, USA.
| | | |
Collapse
|