We report a consistent set of AMBER force-field parameters for the most common phosphorylated amino acids, phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine in different protonation states. The calculation of atomic charges followed the original restrained electrostatic potential fitting procedure used to determine the charges for the parm94/99 parameter set, taking alpha-helical and beta-strand conformations of the corresponding ACE-/NME-capped model peptide backbone into account. Missing force-field parameters were taken directly from the general AMBER force field (gaff) and the parm99 data set with minor modifications, or were newly generated based on ab initio calculations for model systems. Final parameters were validated by geometry optimizations and molecular-dynamics simulations. Template libraries for the phosphorylated amino acids in Leap format and corresponding frcmod parameter files are made available. [Figure: see text].

Assembling phenylpiperazines with 7a-azaindole via different spacer elements, we developed subtype selective dopamine receptor ligands of types 1a,c, 2a, and 3a preferentially interacting with D4, D2, and D3, respectively. To complete this set, the methylthio analogues 2b and 3b exceeding the affinity of 2a and 3a by one order of magnitude and the structural intermediate 1b were synthesized. These chemically similar but biologically divergent target compounds served as molecular probes for radioligand displacement experiments, mutagenesis, and docking studies on homology models based on the recent crystal structure of the beta2-adrenergic receptor. Specific interactions with the highly conserved amino acids Asp3.32 and His6.55 and less conserved residues at positions 2.61, 2.64, 3.28, and 3.29 were identified. Inclusion of a carefully modeled extracellular loop 2 displayed two nonconserved residues in EL2 that differently contribute to ligand binding. Obviously, subtype selectivity is caused by nonconserved but frequently mediated by conserved amino acids.

We present a computational model study designed to simulate the results of time-resolved fluorescence spectra of tryptophan in proteins. In such measurements, the occurrence of more than one fluorescence lifetime is generally attributed to the existence of several tryptophan rotamers and/or structural conformations of the protein structure. The protein system we chose for this initial study is the tetracycline repressor (TetR), an interesting model system for the investigation of the mechanisms of transcriptional regulation. Fluorescence resonance energy transfer (FRET) from tryptophan to tetracycline is frequently observed in complexes of the TetR with the antibiotic tetracycline. We use a combined classical/quantum mechanical approach to model the structure and the spectroscopic properties of the TetR-tetracycline complex. A classical molecular dynamics simulation provides input geometries for semiempirical quantum mechanical/molecular mechanical (QM/MM) single-point configuration interaction (CI) calculations, which are used to calculate tryptophan vertical absorption and fluorescence energies and intensities as well as relative FRET rate constants. These rate constants together with the Einstein coefficients for spontaneous emission and an assumed rate for nonradiative deactivation allow us to simulate fluorescence decay curves with and without FRET and for the entire ensemble as well as for individual rotamers. Our results indicate that the classical "rotamer model", used to explain the multiexponential fluorescence-decay curves of time-resolved tryptophan emission spectra, can be extended to systems with FRET acceptors present in the protein matrix but that the interpretation of the fitted lifetimes is different to that usually used.

Based on the recently described crystal structure of the β2 adrenergic receptor - Gs-protein complex, we report the first molecular-dynamics simulations of ternary GPCR complexes designed to identify the selectivity determinants for receptor-G-protein binding. Long-term molecular dynamics simulations of agonist-bound β2AR-Gαs and D2R-Gαi complexes embedded in a hydrated bilayer environment and computational alanine-scanning mutagenesis identified distinct residues of the N-terminal region of intracellular loop 3 to be crucial for coupling selectivity. Within the G-protein, specific amino acids of the α5-helix, the C-terminus of the Gα-subunit and the regions around αN-β1 and α4-β6 were found to determine receptor recognition. Knowledge of these determinants of receptor-G-protein binding selectivity is essential for designing drugs that target specific receptor/G-protein combinations.

Focusing on the similarity and divergence of GPCR subtypes and their ligand interactions, we generated dopamine D2, D3, and D4 receptor models based on the rhodopsin crystal structure and refined these with an extensive MM/MD protocol. After validation by diagnostic experimental data, subtype-specific relative positions of TM1, 2, 6, and 7 and bending angles of TM7 were found. To sample the conformational space of the complex, we performed simulated-annealing runs of the receptor protein with the sub-nanomolar antagonist spiperone. Docking a representative set of ligands, we were able to identify one superior model for each subtype when excellent correlations between predicted energies of binding and experimental affinities (r2 = 0.72 for D2, 0.91 for D3 and 0.77 for D4) could be observed. Further analysis revealed general ligand interactions with ASP3.32 and aromatic residues in TM6/7 and individual key interactions with TM1 and TM2 residues of the D3 and D4 receptor models, respectively.