Biomedical subjects

Johan Aqvist

Publications and source records attributed to Johan Aqvist.

At least 19 recordsLinked to original sources

PubMed2007

Energetics of codon-anticodon recognition on the small ribosomal subunit.

Recent crystal structures of the small ribosomal subunit have made it possible to examine the detailed energetics of codon recognition on the ribosome by computational methods. The binding of cognate and near-cognate anticodon stem loops to the ribosome decoding center, with mRNA containing the Phe UUU and UUC codons, are analyzed here using explicit solvent molecular dynamics simulations together with the linear interaction energy (LIE) method. The calculated binding free energies are in excellent agreement with experimental binding constants and reproduce the relative effects of mismatches in the first and second codon position versus a mismatch at the wobble position. The simulations further predict that the Leu2 anticodon stem loop is about 10 times more stable than the Ser stem loop in complex with the Phe UUU codon. It is also found that the ribosome significantly enhances the intrinsic stability differences of codon-anticodon complexes in aqueous solution. Structural analysis of the simulations confirms the previously suggested importance of the universally conserved nucleotides A1492, A1493, and G530 in the decoding process.

Anticodon↗

PubMed2006

Computational analysis of plasmepsin IV bound to an allophenylnorstatine inhibitor.

The plasmepsin proteases from the malaria parasite Plasmodium falciparum are attracting attention as putative drug targets. A recently published crystal structure of Plasmodium malariae plasmepsin IV bound to an allophenylnorstatine inhibitor [Clemente, J.C. et al. (2006) Acta Crystallogr. D 62, 246-252] provides the first structural insights regarding interactions of this family of inhibitors with plasmepsins. The compounds in this class are potent inhibitors of HIV-1 protease, but also show nM binding affinities towards plasmepsin IV. Here, we utilize automated docking, molecular dynamics and binding free energy calculations with the linear interaction energy LIE method to investigate the binding of allophenylnorstatine inhibitors to plasmepsin IV from two different species. The calculations yield excellent agreement with experimental binding data and provide new information regarding protonation states of active site residues as well as conformational properties of the inhibitor complexes.

Animals↗

PubMed2006

Computational study of the binding affinity and selectivity of the bacterial ammonium transporter AmtB.

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.

Binding Sites↗

PubMed2006

Calculations of solute and solvent entropies from molecular dynamics simulations.

The translational, rotational and conformational (vibrational) entropy contributions to ligand-receptor binding free energies are analyzed within the standard formulation of statistical thermodynamics. It is shown that the partitioning of the binding entropy into different components is to some extent arbitrary, but an appropriate method to calculate both translational and rotational entropy contributions to noncovalent association is by estimating the configurational volumes of the ligand in the bound and free states. Different approaches to calculating solute entropies using free energy perturbation calculations, configurational volumes based on root-mean-square fluctuations and covariance matrix based quasiharmonic analysis are illustrated for some simple molecular systems. Numerical examples for the different contributions demonstrate that theoretically derived results are well reproduced by the approximations. Calculation of solvent entropies, either using total potential energy averages or van't Hoff plots, are carried out for the case of ion solvation in water. Although convergence problems will persist for large and complex simulation systems, good agreement with experiment is obtained here for relative and absolute ion hydration entropies. We also outline how solvent and solute entropic contributions are taken into account in empirical binding free energy calculations using the linear interaction energy method. In particular it is shown that empirical scaling of the nonpolar intermolecular ligand interaction energy effectively takes into account size dependent contributions to the binding free energy.

Chymotrypsin↗

PubMed2006

Inhibitor binding to the plasmepsin IV aspartic protease from Plasmodium falciparum.

Plasmepsin IV (Plm IV) is one of the aspartic proteases present in the food vacuole of the malaria parasite Plasmodium falciparum involved in host hemoglobin degradation by the parasite. Using a series of previously synthesized plasmepsin inhibitors [Ersmark, K., et al. (2005) J. Med. Chem. 48, 6090-106], we report here experimental data and theoretical analysis of their inhibitory activity toward Plm IV. All compounds share a 1,2-dihydroxyethylene unit as the transition state mimic. They possess symmetric P1 and P1' side chains and either a diacylhydrazine, a five-membered oxadiazole ring, or a retroamide at the P2 and P2' positions. Experimental binding affinities are compared to those predicted by the linear interaction energy (LIE) method and an empirical scoring function, using both a crystal structure and a homology model for the enzyme. Molecular dynamics (MD) simulations of the modeled complexes allow a rational interpretation of the structural determinants for inhibitor binding. A ligand bearing a P2 and P2' symmetric oxadiazole which is devoid of amide bonds is identified both experimentally and theoretically as the most potent inhibitor of Plm IV. For the P2 and P2' asymmetric compounds, the results are consistent with earlier predictions regarding the mode of binding of this class of inhibitors to Plm II. Theoretical estimation of selectivity for some compounds is also reported. Significant features of the Plm IV binding pocket are discussed in comparison to related enzymes, and the results obtained here should be helpful for further optimization of inhibitors.

Animals↗

PubMed2006

Free energy calculations show that acidic P1 variants undergo large pKa shifts upon binding to trypsin.

Serine proteinases and their protein inhibitors belong to one of the most comprehensively studied models of protein-protein interactions. It is well established that the narrow trypsin specificity is caused by the presence of a negatively charged aspartate at the specificity pocket. X-ray crystallography as well as association measurements revealed, surprisingly, that BPTI with glutamatic acid as the primary binding (P1) residue was able to bind to trypsin. Previous free energy calculations showed that there was a substantially unfavorable binding free energy associated with accommodation of ionized P1 Glu at the S1-site of trypsin. In this study, the binding of P1 Glu to trypsin has been systematically investigated in terms of the protonation states of P1 Glu and Asp189, the orientation of Gln192, as well as the possible presence of counterions using the linear interaction energy (LIE) approach and the free energy perturbation (FEP) method. Twenty-four conceivable binding arrangements were evaluated and quantitative agreement with experiments is obtained when the P1 Glu binds in its protonated from. The results suggest that P1 Glu is one of the variants of BPTI that inhibit trypsin strongest at low pH, contrary to the specificity profile of trypsin, suggesting a new regulation mechanism of trypsin-like enzymes.

Aprotinin↗

PubMed2006

Catalysis and linear free energy relationships in aspartic proteases.

Aspartic proteases are receiving considerable attention as potential drug targets in several serious diseases, such as AIDS, malaria, and Alzheimer's disease. These enzymes cleave polypeptide chains, often between specific amino acid residues, but despite the common reaction mechanism, they exhibit large structural differences. Here, the catalytic mechanism of aspartic proteases plasmepsin II, cathepsin D, and HIV-1 protease is examined by computer simulations utilizing the empirical valence bond approach in combination with molecular dynamics and free energy perturbation calculations. Free energy profiles are established for four different substrates, each six amino acids long and containing hydrophobic side chains in the P1 and P1' positions. Our simulations reproduce the catalytic effect of these enzymes, which accelerate the reaction rate by a factor of approximately 10(10) compared to that of the corresponding uncatalyzed reaction in water. The calculations elucidate the origin of the catalytic effect and allow a rationalization of the fact that, despite large structural differences between plasmepsin II/cathepsin D and HIV-1 protease, the magnitude of their rate enhancement is very similar. Amino acid residues surrounding the active site together with structurally conserved water molecules are found to play an important role in catalysis, mainly through dipolar (electrostatic) stabilization. A linear free energy relationship for the reactions in the different enzymes is established that also demonstrates the reduced reorganization energy in the enzymes compared to that in the uncatalyzed water reaction.

Aspartic Acid Endopeptidases↗

PubMed2006

Continuum solvation models in the linear interaction energy method.

The linear interaction energy (LIE) method in combination with two different continuum solvent models has been applied to calculate protein-ligand binding free energies for a set of inhibitors against the malarial aspartic protease plasmepsin II. Ligand-water interaction energies are calculated from both Poisson-Boltzmann (PB) and Generalized Born (GB) continuum models using snapshots from explicit solvent simulations of the ligand and protein-ligand complex. These are compared to explicit solvent calculations, and we find close agreement between the explicit water and PB solvation models. The GB model overestimates the change in solvation energy, and this is caused by consistent underestimation of the effective Born radii in the protein-ligand complex. The explicit solvent LIE calculations and LIE-PB, with our standard parametrization, reproduce absolute experimental binding free energies with an average unsigned error of 0.5 and 0.7 kcal/mol, respectively. The LIE-GB method, however, requires a constant offset to approach the same level of accuracy.

Models, Molecular↗

PubMed2006

Analysis of predictions for the catalytic mechanism of ribosomal peptidyl transfer.

The reaction mechanism of peptide bond formation on the ribosome is now becoming established by results from both experiments and computer simulations. Here, we analyze predictions from molecular dynamics simulations, as well as from new crystal structures, and examine their implications for the mechanisms of peptidyl transfer and peptidyl-tRNA hydrolysis. A number of computational predictions for the peptidyl transfer reaction, including quantitative energetics, stereochemistry, hydrogen bonding network, and role of solvent molecules, are found to be supported and confirmed by kinetic and structural data. The results show that this type of reaction calculations can provide important links between structure and function that cannot be obtained by experimental means.

Catalysis↗

PubMed2005

Probing the effect of point mutations at protein-protein interfaces with free energy calculations.

We have studied the effect of point mutations of the primary binding residue (P1) at the protein-protein interface in complexes of chymotrypsin and elastase with the third domain of the turkey ovomucoid inhibitor and in trypsin with the bovine pancreatic trypsin inhibitor, using molecular dynamics simulations combined with the linear interaction energy (LIE) approach. A total of 56 mutants have been constructed and docked into their host proteins. The free energy of binding could be reliably calculated for 52 of these mutants that could unambiguously be fitted into the binding sites. We find that the predicted binding free energies are in very good agreement with experimental data with mean unsigned errors between 0.50 and 1.03 kcal/mol. It is also evident that the standard LIE model used to study small drug-like ligand binding to proteins is not suitable for protein-protein interactions. Three different LIE models were therefore tested for each of the series of protein-protein complexes included, and the best models for each system turn out to be very similar. The difference in parameterization between small drug-like compounds and protein point mutations is attributed to the preorganization of the binding surface. Our results clearly demonstrate the potential of free energy calculations for probing the effect of point mutations at protein-protein interfaces and for exploring the principles of specificity of hot spots at the interface.

Animals↗

PubMed2005

Synthesis of malarial plasmepsin inhibitors and prediction of binding modes by molecular dynamics simulations.

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.

Animals↗

PubMed2005

Mechanism of peptide bond synthesis on the ribosome.

With the emergence of atomic-resolution crystal structures of bacterial ribosomal subunits, major advances in eliciting structure-function relationships of the translation process are underway. Nevertheless, the detailed mechanism of peptide bond synthesis that occurs on the large ribosomal subunit remains unknown. Separate x-ray structures of aminoacyl-tRNA and peptidyl-tRNA analogues bound to the ribosomal A- and P-sites, however, allow for structural modeling of the active complex in catalysis. Here, we combine available structural data to construct such a model of the peptidyl transfer reaction center with bound substrates. Molecular dynamics and free energy perturbation simulations then are used in combination with an empirical valence bond description of the reaction energy surface to examine possible catalytic mechanisms. Already, simulations of the reactant and tetrahedral intermediate states reveal a stable, preorganized H-bond network poised for catalysis. The most favorable mechanism is found not to involve any general acid-base catalysis by ribosomal groups but an intra-reactant proton shuttling via the P-site adenine O2' oxygen, which follows the attack of the A-site alpha-amino group on the P-site ester. The calculated rate enhancement for this mechanism is approximately 10(5), and the catalytic effect is found to be entirely of entropic origin, in accordance with recent experimental data, and is associated with the reduction of solvent reorganization energy rather than with substrate alignment or proximity. This mechanism also explains the inability of 2'-deoxyadenine P-site substrates to promote peptidyl transfer. The observed H-bond network suggests an important structural role of several universally conserved rRNA residues.

Binding Sites↗

PubMed2005

Resolving the energy paradox of chaperone/usher-mediated fibre assembly.

Periplasmic chaperone/usher machineries are used for assembly of filamentous adhesion organelles of Gram-negative pathogens in a process that has been suggested to be driven by folding energy. Structures of mutant chaperone-subunit complexes revealed a final folding transition (condensation of the subunit hydrophobic core) on the release of organelle subunit from the chaperone-subunit pre-assembly complex and incorporation into the final fibre structure. However, in view of the large interface between chaperone and subunit in the pre-assembly complex and the reported stability of this complex, it is difficult to understand how final folding could release sufficient energy to drive assembly. In the present paper, we show the X-ray structure for a native chaperone-fibre complex that, together with thermodynamic data, shows that the final folding step is indeed an essential component of the assembly process. We show that completion of the hydrophobic core and incorporation into the fibre results in an exceptionally stable module, whereas the chaperone-subunit pre-assembly complex is greatly destabilized by the high-energy conformation of the bound subunit. This difference in stabilities creates a free energy potential that drives fibre formation.

Bacterial Proteins↗

PubMed2005

Exploring blocker binding to a homology model of the open hERG K+ channel using docking and molecular dynamics methods.

Binding of blockers to the human voltage-gated hERG potassium channel is studied using a combination of homology modelling, automated docking calculations and molecular dynamics simulations, where binding affinities are evaluated using the linear interaction energy method. A homology model was constructed based on the available crystal structure of the bacterial KvAP channel and the affinities of a series of sertindole analogues predicted using this model. The calculations reproduce the relative binding affinities of these compounds very well and indicate that both polar interactions near the intracellular opening of the selectivity filter as well as hydrophobic complementarity in the region around F656 are important for blocker binding. These results are consistent with recent alanine scanning mutation experiments on the blocking of the hERG channel by other compounds.

Amino Acid Sequence↗

PubMed2005

Absolute and relative entropies from computer simulation with applications to ligand binding.

A comparison between two related methods, Schlitter's formula and quasiharmonic analysis, for calculating absolute entropies from the covariance matrix of atomic fluctuations using molecular dynamics (MD) simulations is presented. Calculations for a set of organic compounds in the gas phase are compared to the corresponding statistical thermodynamics results for translational and rotational entropies and to experimental data for vibrational entropies. Encouraging agreement is obtained for translational entropies, but for the rotational contribution, both methods fail to reproduce the theoretically calculated values. Absolute and relative vibrational entropies are found to be better reproduced using quasiharmonic analysis compared to Schlitter's formula. For rotational entropies, we propose a method based on the variances in Euler angles, which gives good agreement with theory. Alternative methods for estimating translational entropies based on principal root mean-square (rms) fluctuations of the center of mass are also presented, and these reproduce theoretically calculated values well. These methodologies are applied to the binding of benzene to T4-lysozyme, where close agreement with the literature is obtained for translational and rotational entropies.

Bacteriophage T4↗

PubMed2004

Computational prediction of structure, substrate binding mode, mechanism, and rate for a malaria protease with a novel type of active site.

The histo-aspartic protease (HAP) from the malaria parasite P. falciparum is one of several new promising targets for drug intervention. The enzyme possesses a novel type of active site, but its 3D structure and mechanism of action are still unknown. Here we use a combination of homology modeling, automated docking searches, and molecular dynamics/reaction free energy profile simulations to predict the enzyme structure, conformation of bound substrate, catalytic mechanism, and rate of the peptide cleavage reaction. We find that the computational tools are sufficiently reliable both for identifying substrate binding modes and for distinguishing between different possible reaction mechanisms. It is found that the favored pathway only involves direct participation by the catalytic aspartate, with the neighboring histidine providing critical stabilization (by a factor of approximately 10000) along the reaction. The calculated catalytic rate constant of about 0.1 s(-1) for a hexapeptide substrate derived from the alpha chain of human hemoglobin is in excellent agreement with experimental kinetic data for a similar peptide fragment.

Animals↗

PubMed2004

Binding affinity prediction with different force fields: examination of the linear interaction energy method.

A systematic study of the linear interaction energy (LIE) method and the possible dependence of its parameterization on the force field and system (receptor binding site) is reported. We have calculated the binding free energy for nine different ligands in complex with P450cam using three different force fields (Amber95, Gromos87, and OPLS-AA). The results from these LIE calculations using our earlier parameterization give relative free energies of binding that agree remarkably well with the experimental data. However, the absolute energies are too positive for all three force fields, and it is clear that an additional constant term (gamma) is required in this case. Out of five examined LIE models, the same one emerges as the best for all three force fields, and this, in fact, corresponds to our earlier one apart from the addition of the constant gamma, which is almost identical for the three force fields. Thus, the present free energy calculations clearly indicate that the coefficients of the LIE method are independent of the force field used. Their relation to solvation free energies is also demonstrated. The only free parameter of the best model is gamma, which is found to depend on the hydrophobicity of the binding site. We also attempt to quantify the binding site hydrophobicity of four different proteins which shows that the ordering of gamma's for these sites reflects the fraction of hydrophobic surface area.

Journal Article↗

PubMed2004

Comparative analysis of putative agonist-binding modes in the human A1 adenosine receptor.

A recent study reported a model of the human A(1) adenosine receptor and its agonist binding site, proposing two putative binding modes in the same binding site for the natural agonist, adenosine. The present work investigates the flexibility of this binding site by exhaustive exploration with the natural agonist and with three other adenosine derivatives: N6-cyclopentyladenosine (CPA), 2-chloro-N6-cyclopentyladenosine (CCPA), and 5'-N-ethylcarboxamidoadenosine (NECA). Our aim was to find a common binding mode for agonists that would explain the role in the binding process of the different substitutions allowed at the 2, N6, and 5' positions of adenosine. This problem was addressed through docking simulations, molecular dynamics studies, and estimations of the ligand-binding free energy with both the AUTODOCK scoring function and the linear interaction energy (LIE) approach. The results point to a single receptor-binding position that explains the effects of the different chemical modifications on the adenosine derivatives considered here.

Adenosine↗