No abstract

We report results from microscopic molecular dynamics and free energy perturbation simulations of substrate binding and selectivity for the Escherichia coli high-affinity ammonium transporter AmtB. The simulation system consists of the protein embedded in a model membrane/water surrounding. The calculated absolute binding free energies for the external NH(4)(+) ions are between -5.8 and -7.3 kcal/mol and are in close agreement with experimental data. The apparent pK(a) of the bound NH(4)(+) increases by more than 4 units, indicating a preference for binding ammonium ion and not neutral ammonia. The external binding site is also selective for NH(4)(+) toward monovalent metal cations by 2.4-4.4 kcal/mol. The externally bound NH(4)(+) shows strong electrostatic interactions with the proximal buried Asp160, stabilized in the anionic form, whereas the interactions with the aromatic rings of Phe107 and Trp148, lining the binding cavity, are less pronounced. Simulated mutation of the highly conserved Asp160 to Asn reduces the pK(a) of the bound ammonium ion by approximately 7 units and causes loss of its binding. The calculations further predict that the substrate affinity of E. coli AmtB depends on the ionization state of external histidines. The computed free energies of hypothetical intermediate states related to transfer of NH(3), NH(4)(+), or H(2)O from the external binding site to the first position inside the internal channel pore favor permeation of the neutral species through the channel interior. However, the predicted change in the apparent pK(a) of NH(4)(+) upon translocation from the external site, Am1, to the first internal site, Am2, indicates that ammonium ion becomes deprotonated only when it enters the channel interior.

A series of inhibitors of the malarial aspartic proteases Plm I and II have been synthesized with L-mannitol as precursor. These inhibitors are characterized by either a diacylhydrazine or a five-membered oxadiazole ring replacing backbone amide functionalities. Molecular dynamics simulations were applied in the design process. The computationally predicted Plm II Ki values were generally in excellent agreement with the biological results. The diacylhydrazine was found to be superior over the oxadiazole as an amide bond replacement in the Plm I and II inhibitors studied. An extensive flexibility of the S2' pocket was captured by the simulations predicting the binding mode of the unsymmetrical inhibitors. Plm I and II inhibitors with single digit nanomolar Ki values devoid of inhibitory activity toward human Cat D were identified. One compound, lacking amide bonds, was found to be Plm IV selective and very potent, with a Ki value of 35 nM.