Homology modeling and molecular dynamics simulations of MUC1-9/H-2Kb complex suggest novel binding interactions

International audienceHuman MUC1 is over-expressed in human adenocarcinomas and has been used as a target for immunotherapy studies. The 9-mer MUC1-9 peptide has been identified as one of the peptides which binds to murine MHC class I H-2K. The structure of MUC1-9 in complex with H-2K has been modeled and simulated with classical molecular dynamics, based on the x-ray structure of the SEV9 peptide/H-2K complex. Two independent trajectories with the solvated complex (10 ns in length) were produced. Approximately 12 hydrogen bonds were identified during both trajectories to contribute to peptide/MHC complex, as well as 1-2 water mediated hydrogen bonds. Stability of the complex was also confirmed by buried surface area analysis, although the corresponding values were about 20% lower than those of the original x-ray structure. Interestingly, a bulged conformation of the peptide's central region, partially characterized as a -turn, was found exposed form the binding groove. In addition, P1 and P9 residues remained bound in the A and F binding pockets, even though there was a suggestion that P9 was more flexible. The complex lacked numerous water mediated hydrogen bonds that were present in the reference peptide x-ray structure. Moreover, local displacements of residues Asp4, Thr5 and Pro9 resulted in loss of some key interactions with the MHC molecule. This might explain the reduced affinity of the MUC1-9 peptide, relatively to SEV9, for the MHC class I H-2K

Three-dimensional structure of the histidine-containing phosphocarrier protein (HPr) from Enterococcus faecalis in solution

The histidine-containing phosphocarrier protein (HPr) transfers a phosphate group between components of the prokaryotic phosphoenolpyruvate-dependent phosphotransferase system (PTS), which is finally used to phosphorylate the carbohydrate transported by the PTS through the cell membrane. Recently it has also been found to act as an intermediate in the signaling cascade that regulates transcription of genes related to the carbohydrate-response system. Both functions involve phosphorylation/dephosphorylation reactions, but at different sites. Using multidimensional (1)H-NMR spectroscopy and angular space simulated annealing calculations, we determined the structure of HPr from Enterococcus faecalis in aqueous solution using 1469 distance and 44 angle constraints derived from homonuclear NMR data. It has a similar overall fold to that found in HPrs from other organisms. Four beta strands, A, B, C, D, encompassing residues 2-7, 32-37, 40-42 and 60-66, form an antiparallel beta sheet lying opposite the two antiparallel alpha helices, a and c (residues 16-26 and 70-83). A short alpha helix, b, from residues 47-53 is also observed. The pairwise root mean square displacement for the backbone heavy atoms of the mean of the 16 NMR structures to the crystal structure is 0.164 nm. In contrast with the crystalline state, in which a torsion angle strain in the active-center loop has been described [Jia, Z., Vandonselaar, M., Quail, J.W. & Delbaere, L.T.J. (1993) Nature (London) 361, 94-97], in the solution structure, the active-site His15 rests on top of helix a, and the phosphorylation site N(delta 1) of the histidine ring is oriented towards the surface, making it easily accessible to the solvent. Back calculation of the 2D NOESY NMR spectra from both the NMR and X-ray structures shows that the active-center structure derived by X-ray crystallography is not compatible with experimental data recorded in solution. The observed torsional strain must either be a crystallization artefact or represents a conformational state that exists only to a small extent in solution

Determination of mean and standard deviation of dihedral angles

Backbone torsional angles are a characteristic and useful parameter for the description and characterisation of protein structures determined by x-ray crystallography or NMR spectroscopy. For the comparison of an ensemble of three-dimensional structures the calculation of the statistical parameters mean and standard deviation would be very useful. However, they are not defined unambiguously for periodic quantities such as the dihedral angles. In this paper a plausible and unique definition of these parameters is introduced and a straightforward method for their calculation is given

Pressure-induced local unfolding of the Ras binding domain of RalGDS

The reliable prediction of the precise three-dimensional structure of proteins from their amino acid sequence is a major, still unresolved problem in biochemistry. Pressure is a parameter that controls folding/unfolding transitions of proteins through the volume change DeltaV of the protein-solvent system. By varying the pressure from 30 to 2,000 bar we detected using 15N/ 1H 2D NMR spectroscopy a unique equilibrium unfolding intermediate I in the Ras binding domain of the Ral guanine nucleotide dissociation stimulator (Ral GDS). It is characterized by a local melting of specific structural elements near hydrophobic cavities while the overall folded structure is maintained