Cation-π interactions of methylated ammonium ions: a quantum mechanical study.

Cation-π interactions of methylated ammonium ions play a key role in a broad range of biochemical systems. These include methyl-lysine binding proteins which bind to methylated sites on histone proteins, lysine demethylase enzymes which demethylate these sites, and neurotransmitter receptor complexes which bind choline-derived ligands. Recognition in these systems is achieved through an 'aromatic cage' motif in the binding site.

The role of tyrosine sulfation in the dimerization of the CXCR4:SDF-1 complex.

Oligomerization of G protein-coupled receptors is a recognized mode of regulation of receptor activities, with alternate oligomeric states resulting in different signaling functions. The CXCR4 chemokine receptor is a G protein-coupled receptor that is post-translationally modified by tyrosine sulfation at three sites on its N-terminus (Y7, Y12, Y21), leading to enhanced affinity for its ligand, stromal cell derived factor (SDF-1, also called CXCL12). The complex has been implicated in cancer metastasis and is a therapeutic target in cancer treatment.

Hydrogen bond strengths in phosphorylated and sulfated amino acid residues.

Post-translational modification by the addition of an oxoanion functional group, usually a phosphate group and less commonly a sulfate group, leads to diverse structural and functional consequences in protein systems. Building upon previous studies of the phosphoserine residue (pSer), we address the distinct nature of hydrogen bonding interactions in phosphotyrosine (pTyr) and sulfotyrosine (sTyr) residues. We derive partial charges for these modified residues and then study them in the context of molecular dynamics simulation of model tripeptides and sulfated protein complexes, potentials of mean force for interacting residue pairs, and a survey of the interactions of modified residues among experimental protein structures.

A molecular mechanics approach to modeling protein-ligand interactions: relative binding affinities in congeneric series.

We introduce the "Prime-ligand" method for ranking ligands in congeneric series. The method employs a single scoring function, the OPLS-AA/GBSA molecular mechanics/implicit solvent model, for all stages of sampling and scoring. We evaluate the method using 12 test sets of congeneric series for which experimental binding data is available in the literature, as well as the structure of one member of the series bound to the protein.

Automated site preparation in physics-based rescoring of receptor ligand complexes.

Hydrogen atoms are not typically observable in X-ray crystal structures, but inferring their locations is often important in structure-based drug design. In addition, protonation states of the protein can change in response to ligand binding, as can the orientations of OH groups, a subtle form of "induced fit." We implement and evaluate an automated procedure for optimizing polar hydrogens in protein-binding sites in complex with ligands.

Prediction of protein loop geometries in solution.

The ability to determine the structure of a protein in solution is a critical tool for structural biology, as proteins in their native state are found in aqueous environments. Using a physical chemistry based prediction protocol, we demonstrate the ability to reproduce protein loop geometries in experimentally derived solution structures. Predictions were run on loops drawn from (1)NMR entries in the Protein Databank (PDB), and from (2) the RECOORD database in which NMR entries from the PDB have been standardized and re-refined in explicit solvent.

Strengths of hydrogen bonds involving phosphorylated amino acid side chains.

Post-translational phosphorylation plays a key role in regulating protein function. Here, we provide a quantitative assessment of the relative strengths of hydrogen bonds involving phosphorylated amino acid side chains (pSer, pAsp) with several common donors (Arg, Lys, and backbone amide groups). We utilize multiple levels of theory, consisting of explicit solvent molecular dynamics, implicit solvent molecular mechanics, and quantum mechanics with a self-consistent reaction field treatment of solvent.

Crystal packing effects on protein loops.

The effects of crystal packing on protein loop structures are examined by (1) a comparison of loops in proteins that have been crystallized in alternate packing arrangements, and (2) theoretical prediction of loops both with and without the inclusion of the crystal environment. Results show that in a minority of cases, loop geometries are dependent on crystal packing effects. Explicit representation of the crystal environment in a loop prediction algorithm can be used to model these effects and to reconstruct the structures, and relative energies, of a loop in alternative packing environments.

A hierarchical approach to all-atom protein loop prediction.

The application of all-atom force fields (and explicit or implicit solvent models) to protein homology-modeling tasks such as side-chain and loop prediction remains challenging both because of the expense of the individual energy calculations and because of the difficulty of sampling the rugged all-atom energy surface. Here we address this challenge for the problem of loop prediction through the development of numerous new algorithms, with an emphasis on multiscale and hierarchical techniques. As a first step in evaluating the performance of our loop prediction algorithm, we have applied it to the problem of reconstructing loops in native structures; we also explicitly include crystal packing to provide a fair comparison with crystal structures.