Intraprotein electrostatics derived from first principles: Divide-and-conquer approaches for QM/MM calculations

Two divide-and-conquer (DAQ) approaches for building multipole-based molecular electrostatic potentials of proteins are presented and evaluated for use in QM/MM calculations. One approach is a further development of the neutralization method of Bellido and Rullmann (J Comput Chem 1989, 10, 479-487) while the other is based on removing part of the electron density before performing the multipole expansion. Both methods create systems with integer charges without using charge renormalization. To determine their performance in terms of location of cuts and distance to QM region, the new DAQ approaches are tested in calculations of the proton affinity of N(zeta) of Lys55 in the inhibitor turkey ovomucoid third domain. Finally, the two methods are used to build a variety of MM regions, applied to calculations of the pK(a) of Lys55, and compared to other computational methodologies in which force field charges are employed.

The solution structure of the Ga(III)-bleomycin A2 complex resolved by NMR and molecular modeling; interaction with d(CCAGGCCTGG)

The solution structure of the Ga(III)-bleomycin A2 complex (GaBLM) has been determined using 2D NMR methods in combination with molecular dynamics calculations. Complete assignment of the amide and amine protons, observation of 80 NOEs and measurement of 15 (3)JH(-H) coupling constants provided us with a well-defined structure using a restrained simulated annealing protocol. On the basis of distance and dihedral angle constraints agreement, along with potential energy considerations, the favored model is a five-coordinate complex with the primary amine of beta-aminoalanine holding the axial position of a distorted tetragonal pyramid. The disaccharide moiety of GaBLM is not a ligand, sharing the same side of the equatorial plane with the axial amine ligand. Titration of the self-complementary oligonucleotide d(CCAGGCCTGG) with GaBLM results in the formation of only one 1:1 complex in slow exchange on the NMR time scale. Our data indicate that the bithiazole moiety intercalates between the C6*G15 and C7*G14 base pairs, in a similar mode to that reported by earlier studies. Structural implications and comparisons to other metallo-bleomycins are discussed.

Conformational analysis and automated receptor docking of selective arylacetamide-based kappa-opioid agonists

The three-dimensional structure, dynamics, and binding modes of representative kappa-opioid agonists of the arylacetamide class (U50, 488; U69,593; U62,066; CI-977; ICI199,441; ICI197,067; BRL52,537; and BRL52,656) have been investigated using molecular modeling techniques. Systematic exploration of the conformational space of the ligand combined with molecular dynamics (MD) simulations in water revealed consistent conformational preferences for all the kappa-agonists in this series. The results were further compared with available X-ray and 1D- and 2D-NMR data to identify potential "lead" conformers for molecular docking. Ligand binding modes were initially determined using automated docking of two of the ligands (U50,488 and BRL52,537) to the kappa-opioid receptor. Extrapolation of the predicted binding mode to other members in this ligand series revealed similar docking preferences, with each ligand docked along the receptor helical axis. The binding modes were further refined using MD simulations of the receptor-ligand complexes. The results show a that salt bridge is formed between the amino proton of the ligands and the carboxylate group of Asp138 in TM3. This interaction most likely serves as a key anchoring point for the agonist association. Additional ligand contacts were noted with kappa-specific residues Ile294, Leu295, and Ala298, which may, in part, explain the kappa-selectivity in this series. In comparing the arylacetamides with opiate-based ligands, no evidence was found to link these classes through a common binding motif (except for the ion pair). The binding site model was also applied to explain the enantiomeric preference of U50,488 and to provide insight to the mu/kappa-selectivity of representative ligands in this series. Overall, the results provide a structure-based rationale for ligand recognition that is consistent both with site-directed mutagenesis experiments and structure-function relationship data.