This subproject is one of many research subprojects utilizing theresources provided by a Center grant funded by NIH/NCRR. The subproject andinvestigator (PI) may have received primary funding from another NIH source,and thus could be represented in other CRISP entries. The institution listed isfor the Center, which is not necessarily the institution for the investigator.The purpose of this study is to determine the pKa of sphingosine 1-phosphate (S1P) when bound to a G Protein Coupled Receptor (GPCR), in order to better understand S1P-receptor interactions. Such understanding may contribute to development of models of agonists for the GPCR, which may have potential as future drugs. We previously derived models for the interaction of S1P with the S1P1 receptor, a GPCR family member receptor, using homology modeling and molecular mechanics. However, the charge and pKa of the S1P phosphate are unknown in the binding environment of S1P1. The binding environment contains two arginine, a lysine, an isoleucine, a tryptophan and glutamic acid residues. We will compute the pKa of the S1P phosphate in S1P1, by using the program GAMESS. The pKa determination is based on the free energy G of the proton exchange reaction: HA(aq) + R-(aq) A-(aq) + RH(aq) by the equation pKa = 5.66 + G/1.36 = 5.66 + {[G(A-)  G(HA)]  [G(R-)  G(RH)]}/1.36 Here, 5.66 is the experimental pKa value of O-phosphoethanolamine (R) at 298K, and A is S1P in the S1P1 receptor. The value 1.36 is RTln10 at the same temperature in kilocalories per mol. Each free energy term is a sum of solvation free energy and electronic energy. The electronic energy is composed of the potential energy of the electrons and nuclei and the kinetic energy of the electrons, and will be calculated using quantum mechanical technique [MP2/6-3+G(2d,p)]. The solvation energy will be calculated by the polarized continuum model (PCM)*, which represents the solvent as a continuum surrounding the solute. We will optimize the structure, by minimizing the total energy with a recently developed algorithm**, which includes solvation energy calculation, and therefore optimize the free energy at the level PCM/RHF/6-31G(d). Diffuse functions will be added to all heavy atoms in negatively charged groups, and the UAHF radii will define the molecular cavities. Reference: * Cances E, Mennucci B, Tomasi J. A new integral equation formalism for the polarizable continuum model: theoretical background and applications to isotropic and anisotropic dielectrics. J Chem Phys 1997; 107(8):3032-3041 ** Li H and Jensen J H. Improving the efficiency and convergence of geometry optimization with the polarizable continuum model: new energy gradients and molecular surface tessellation. J Comp Chem 2004; 25(12):1449-62

National Institute of Health (NIH)
National Center for Research Resources (NCRR)
Biotechnology Resource Grants (P41)
Project #
Application #
Study Section
Special Emphasis Panel (ZRG1-BCMB-Q (40))
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
Carnegie-Mellon University
Biostatistics & Other Math Sci
Schools of Arts and Sciences
United States
Zip Code
Simakov, Nikolay A; Kurnikova, Maria G (2018) Membrane Position Dependency of the pKa and Conductivity of the Protein Ion Channel. J Membr Biol 251:393-404
Yonkunas, Michael; Buddhadev, Maiti; Flores Canales, Jose C et al. (2017) Configurational Preference of the Glutamate Receptor Ligand Binding Domain Dimers. Biophys J 112:2291-2300
Hwang, Wonmuk; Lang, Matthew J; Karplus, Martin (2017) Kinesin motility is driven by subdomain dynamics. Elife 6:
Earley, Lauriel F; Powers, John M; Adachi, Kei et al. (2017) Adeno-associated Virus (AAV) Assembly-Activating Protein Is Not an Essential Requirement for Capsid Assembly of AAV Serotypes 4, 5, and 11. J Virol 91:
Murty, Vishnu P; Calabro, Finnegan; Luna, Beatriz (2016) The role of experience in adolescent cognitive development: Integration of executive, memory, and mesolimbic systems. Neurosci Biobehav Rev 70:46-58
Subramanian, Sandeep; Chaparala, Srilakshmi; Avali, Viji et al. (2016) A pilot study on the prevalence of DNA palindromes in breast cancer genomes. BMC Med Genomics 9:73
Ramakrishnan, N; Tourdot, Richard W; Radhakrishnan, Ravi (2016) Thermodynamic free energy methods to investigate shape transitions in bilayer membranes. Int J Adv Eng Sci Appl Math 8:88-100
Zhang, Yimeng; Li, Xiong; Samonds, Jason M et al. (2016) Relating functional connectivity in V1 neural circuits and 3D natural scenes using Boltzmann machines. Vision Res 120:121-31
Lee, Wei-Chung Allen; Bonin, Vincent; Reed, Michael et al. (2016) Anatomy and function of an excitatory network in the visual cortex. Nature 532:370-4
Jurkowitz, Marianne S; Patel, Aalapi; Wu, Lai-Chu et al. (2015) The YhhN protein of Legionella pneumophila is a Lysoplasmalogenase. Biochim Biophys Acta 1848:742-51

Showing the most recent 10 out of 292 publications