Search PubMed⌕ Search

Biomedical subjects

Johan Aqvist

Publications and source records attributed to Johan Aqvist.

27 records · Page 2Linked to original sources

Potent inhibitors of the Plasmodium falciparum enzymes plasmepsin I and II devoid of cathepsin D inhibitory activity.

The hemoglobin-degrading aspartic proteases plasmepsin I (Plm I) and plasmepsin II (Plm II) of the malaria parasite Plasmodium falciparum have lately emerged as putative drug targets. A series of C(2)-symmetric compounds encompassing the 1,2-dihydroxyethylene scaffold and a variety of elongated P1/P1' side chains were synthesized via microwave-assisted palladium-catalyzed coupling reactions. Binding affinity calculations with the linear interaction energy method and molecular dynamics simulations reproduced the experimental binding data obtained in a Plm II assay with very good accuracy. Bioactive conformations of the elongated P1/P1' chains were predicted and agreed essentially with a recent X-ray structure. The compounds exhibited picomolar to nanomolar inhibition constants for the plasmepsins and no measurable affinity to the human enzyme cathepsin D. Some of the compounds also demonstrated significant inhibition of parasite growth in cell culture. To the best of our knowledge, these plasmepsin inhibitors represent the most selective reported to date and constitute promising lead compounds for further optimization.

Amides↗

Design, synthesis, and computational affinity prediction of ester soft drugs as inhibitors of dihydrofolate reductase from Pneumocystis carinii.

A series of dihydrofolate reductase (DHFR) inhibitors, where the methylenamino-bridge of non-classical inhibitors was replaced with an ester function, have been prepared as potential soft drugs intended for inhalation against Pneumocystis carinii pneumonia (PCP). Several of the new ester-based inhibitors that should serve as good substrates for the ubiquitous esterases and possibly constitute safer alternatives to metabolically stable DHFR inhibitors administered orally, were found to be potent inhibitors of P. carinii DHFR (pcDHFR). Although the objectives of the present program is to achieve a favorable toxicity profile by applying the soft drug concept, a high preference for inhibition of the fungal DHFR versus the mammalian DHFR is still desirable to suppress host toxicity at the site of administration. Compounds with a slight preference for the fungal enzyme were identified. The selection of the target compounds for synthesis was partly guided by an automated docking and scoring procedure as well as molecular dynamics simulations. The modest selectivity of the synthesized inhibitors was reasonably well predicted, although a correct ranking of the relative affinities was not successful in all cases.

Drug Design↗

Structure-activity relationship for extracellular block of K+ channels by tetraalkylammonium ions.

External tetraalkylammonium ion binding to potassium channels is studied using microscopic molecular modelling methods and the experimental structure of the KcsA channel. Relative binding free energies of the KcsA complexes with Me4N+, Et4N+, and n-Pr4N+ are calculated with the molecular dynamics free energy perturbation approach together with automated ligand docking. The four-fold symmetry of the entrance cavity formed by the Tyr82 residues is found to provide stronger binding for the D2d than for the S4 conformation of the ligands. In agreement with experiment the Et4N+ blocker shows several kcal/mol better binding than the other tetraalkylammonium ions.

Bacterial Proteins↗

Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine.

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.

Amino Acid Sequence↗

C2-symmetric inhibitors of Plasmodium falciparum plasmepsin II: synthesis and theoretical predictions.

A series of C(2)-symmetric compounds with a mannitol-based scaffold has been investigated, both theoretically and experimentally, as Plm II inhibitors. Four different stereoisomers with either benzyloxy or allyloxy P1/P1' side chains were studied. Computational ranking of the binding affinities of the eight compounds was carried out using the linear interaction energy (LIE) method relying on a complex previously determined by crystallography. Within both series of isomers the theoretical binding energies were in agreement with the enzymatic measurements, illustrating the power of the LIE method for the prediction of ligand affinities prior to synthesis. The structural models of the enzyme-inhibitor complexes obtained from the MD simulations provided a basis for interpretation of further structure-activity relationships. Hence, the affinity of a structurally similar ligand, but with a different P2/P2' substituent was examined using the same procedure. The predicted improvement in binding constant agreed well with experimental results.

Animals↗

The catalytic power of ketosteroid isomerase investigated by computer simulation.

Ketosteroid isomerase (KSI) catalyzes the isomerization of Delta(5)-3-ketosteroids and Delta(4)-3-ketosteroids at very high rates. Here we examine the principles underlying the catalytic efficiency of KSI by computer simulations using the empirical valence bond method in combination with molecular dynamics free energy perturbation simulations. The simulations reproduce available kinetic and structural data very well and allow us to examine several features of the catalytic mechanism in detail. It is found that about 60% of the rate enhancement is due to stabilization of the negatively charged dienolate intermediate by hydrogen bonding. The critical H-bond between Tyr16 and the intermediate is found to be a normal ionic H-bond with the preferred proton location on the tyrosine residue. The remaining 40% of the catalytic effect originates from a reduction of the reorganization energy of the reaction. The possibility of an active site water molecule occupying the empty cavity adjacent to the catalytic base (Asp40) is also addressed. The existence of such a water molecule could explain how the enzyme manages to maintain a low pK(a) for the general base residue.

Binding Sites↗

Computational study of the influence of solvent on (16)O/(18)O equilibrium isotope effects in phosphate deprotonation reactions.

Results from theoretical calculations of (16)O/(18)O equilibrium isotope effects (EIEs) on deprotonation of phosphate and methyl phosphate monoanions as well as their deuterated counterparts are reported. The EIEs are calculated from the Bigeleisen equation using harmonic vibrational frequencies from several quantum mechanical methods (HF, DFT, MP2, and AM1). All methods correctly predict the qualitative trends in the EIEs related to the different isotope substitutions. However, the calculated gas-phase values are found to be systematically higher than those experimentally observed in aqueous solution. On the other hand, the addition of explicit solvent molecules (up to 24 waters) in the first solvation shells of the phosphate ion substantially improves the calculated EIE, which approaches the experimental value with increasing size of the water cluster. The large effects of surrounding water molecules on the phosphate deprotonation EIE can be explained by the strong solute-solvent interactions, which result in solvent coupled vibrational modes of the phosphate ions.

Deuterium↗

Ligand binding affinities from MD simulations.

Simplified free energy calculations based on force field energy estimates of ligand-receptor interactions and thermal conformational sampling have emerged as a useful tool in structure-based ligand design. Here we give an overview of the linear interaction energy (LIE) method for calculating ligand binding free energies from molecular dynamics simulations. A notable feature is that the binding energetics can be predicted by considering only the intermolecular interactions of the ligand in the associated and dissociated states. The approximations behind this approach are examined, and different parametrizations of the model are discussed. LIE-type methods appear particularly promising for computational "lead optimization". Recent applications to protein-protein interactions and ion channel blocking are also discussed.

Animals↗