Molecular recognition in the sphingosine 1-phosphate receptor family.

Computational modeling and its application in ligand screening and ligand receptor interaction studies play important roles in structure-based drug design. A series of sphingosine 1-phosphate (S1P) receptor ligands with varying potencies and receptor selectivities were docked into homology models of the S1P(1-5) receptors. These studies provided molecular insights into pharmacological trends both across the receptor family as well as at single receptors. This study identifies ligand recognition features that generalize across the S1P receptor family, features unique to the S1P(4) and S1P(5) receptors, and suggests significant structural differences of the S1P(2) receptor. Docking results reveal a previously unknown sulfur-aromatic interaction between the S1P(4) C5.44 sulfur atom and the phenyl ring of benzimidazole as well as pi-pi interaction between F3.33 of S1P(1,4,5) and aromatic ligands. The findings not only confirm the importance of a cation-pi interaction between W4.64 and the ammonium of S1P at S1P(4) but also predict the same interaction at S1P(5). S1P receptor models are validated for pharmacophore development including database mining and new ligand discovery and serve as tools for ligand optimization to improve potency and selectivity.

Sphingosine 1-phosphate pKa and binding constants: intramolecular and intermolecular influences.

The dissociation constant for an ionizable ligand binding to a receptor is dependent on its charge and therefore on its environmentally-influenced pKa value. The pKa values of sphingosine 1-phosphate (S1P) were studied computationally in the context of the wild type S1P1 receptor and the following mutants: E3.29Q, E3.29A, and K5.38A. Calculated pKa values indicate that S1P binds to S1P1 and its site mutants with a total charge of -1, including a +1 charge on the ammonium group and a -2 charge on the phosphate group. The dissociation constant of S1P binding to these receptors was studied as well. The models of wild type and mutant proteins originated from an active receptor model that was developed previously. We used ab initio RHF/6-31+G(d) to optimize our models in aqueous solution, where the solvation energy derivatives are represented by conductor-like polarizable continuum model (C-PCM) and integral equation formalism polarizable continuum model (IEF-PCM). Calculation of the dissociation constant for each mutant was determined by reference to the experimental dissociation constant of the wild type receptor. The computed dissociation constants of the E3.29Q and E3.29A mutants are three to five orders of magnitude higher than those for the wild type receptor and K5.38A mutant, indicating vital contacts between the S1P phosphate group and the carboxylate group of E3.29. Computational dissociation constants for K5.38A, E3.29A, and E3.29Q mutants were compared with experimentally determined binding and activation data. No measurable binding of S1P to the E3.29A and E3.29Q mutants was observed, supporting the critical contacts observed computationally. These results validate the quantitative accuracy of the model.

Determinants of cysteine pKa values in creatine kinase and alpha1-antitrypsin.

The structural determinants of the unusually low pK(a) values of Cys282 in human creatine kinase and Cys232 in alpha1-antitrypsin were studied computationally. We have demonstrated that hydrogen bonding to the cysteine residue is the prime determinant for both proteins. In the case of creatine kinase, the hydrogen bond donors are a serine side chain and an amide NH-group, while in alpha1-antitrypsin the donor is an amide NH. Each hydrogen bond lowers the pK(a) by between 0.8 and 1.5 pH units. The 1.1-unit lowering due to the Ser284-Cys282 hydrogen bond is in good agreement with the 1.2-unit difference between the Cys282 pK(a) value of wild-type and the S284A mutant of creatine kinase.

Exploring the role of the active site cysteine in human muscle creatine kinase.

All known guanidino kinases contain a conserved cysteine residue that interacts with the non-nucleophilic eta1-nitrogen of the guanidino substrate. Site-directed mutagenesis studies have shown that this cysteine is important, but not essential for activity. In human muscle creatine kinase (HMCK) this residue, Cys283, forms part of a conserved cysteine-proline-serine (CPS) motif and has a pKa about 3 pH units below that of a regular cysteine residue. Here we employ a computational approach to predict the contribution of residues in this motif to the unusually low cysteine pKa. We calculate that hydrogen bonds to the hydroxyl and to the backbone amide of Ser285 would both contribute approximately 1 pH unit, while the presence of Pro284 in the motif lowers the pKa of Cys283 by a further 1.2 pH units. Using UV difference spectroscopy the pKa of the active site cysteine in WT HMCK and in the P284A, S285A, and C283S/S285C mutants was determined experimentally. The pKa values, although consistently about 0.5 pH unit lower, were in broad agreement with those predicted. The effect of each of these mutations on the pH-rate profile was also examined. The results show conclusively that, contrary to a previous report (Wang et al. (2001) Biochemistry 40, 11698-11705), Cys283 is not responsible for the pKa of 5.4 observed in the WT V/K(creatine) pH profile. Finally we use molecular dynamics simulations to demonstrate that, in order to maintain the linear alignment necessary for associative inline transfer of a phosphoryl group, Cys283 needs to be ionized.

Proton transport through water-filled carbon nanotubes.

Proton transfer along 1D chains of water molecules inside carbon nanotubes is studied by simulations. Ab initio molecular dynamics and an empirical valence bond model yield similar structures and time scales. The proton mobility along 1D water chains exceeds that in bulk water by a factor of 40, but is reduced if orientational defects are present. Excess protons interact with hydrogen-bonding defects through long-range electrostatics, resulting in coupled motion of protons and defects.