Search PubMed⌕ Search

Biomedical subjects

J Aqvist

Publications and source records attributed to J Aqvist.

At least 19 recordsLinked to original sources
PubMed2000

Computational predictions of binding affinities to dihydrofolate reductase: synthesis and biological evaluation of methotrexate analogues.

The relative binding affinities to human dihydrofolate reductase of four new potential antifolates, containing ester linkages between the two aromatic systems, were estimated by free energy perturbation simulations. The ester analogue, predicted to exhibit the highest binding affinity to human dihydrofolate reductase, and a reference ester (more structurally related to methotrexate) were synthesized. As deduced from the measured IC(50) values, the calculated ranking of the ligands was correct although a greater difference in affinity was indicated by the experimental measurements. Among the new antifolates the most potent inhibitor exhibited a similar pharmacokinetic profile to methotrexate but lacked activity in a complex antiarthritic model in rat in vivo.

Animals↗
PubMed2000

A computational study of ion binding and protonation states in the KcsA potassium channel.

We report results from microscopic molecular dynamics and free energy perturbation simulations of the KcsA potassium channel based on its experimental atomic structure. Conformational properties of selected amino acid residues as well as equilibrium positions of K(+) ions inside the selectivity filter and the internal water cavity are examined. Positions three and four (counting from the extracellular site) in the experimental structure correspond to distinctly separate binding sites for K(+) ions inside the selectivity filter. The protonation states of Glu71 and Asp80, which are close to each other and to the selectivity filter, as well as K(+) binding energies are determined using free energy perturbation calculations. The Glu71 residue which is buried inside a protein cavity is found to be most stable in the neutral form while the solvent exposed Asp80 is ionized. The channel altogether exothermically binds up to three ions, where two of them are located inside the selectivity filter and one in the internal water cavity. Ion permeation mechanisms are discussed in relation to these results.

Bacterial Proteins↗
PubMed2000

Computer simulation of primary kinetic isotope effects in the proposed rate-limiting step of the glyoxalase I catalyzed reaction.

The proposed rate-limiting step of the glyoxalase I catalyzed reaction is the proton abstraction from the C1 carbon of the substrate by Glu(172). Here we examine primary kinetic isotope effects and the influence of quantum dynamics on this process by computer simulations. The calculations utilize the empirical valence bond method in combination with the molecular dynamics free energy perturbation technique and path integral simulations. For the enzyme-catalyzed reaction a H/D kinetic isotope effect of 5.0 +/- 1. 3 is predicted in reasonable agreement with the experimental result of about 3. Furthermore, the magnitude of quantum mechanical effects is found to be very similar for the enzyme reaction and the corresponding uncatalyzed process in solution, in agreement with other studies. The problems associated with attaining the required accuracy in order for the present approach to be useful as a diagnostic tool for the study of enzyme reactions are also discussed.

Catalysis↗
PubMed2000

Ion permeation mechanism of the potassium channel.

Ion-selective channels enable the specific permeation of ions through cell membranes and provide the basis of several important biological functions; for example, electric signalling in the nervous system. Although a large amount of electrophysiological data is available, the molecular mechanisms by which these channels can mediate ion transport remain a significant unsolved problem. With the recently determined crystal structure of the representative K+ channel (KcsA) from Streptomyces lividans, it becomes possible to examine ion conduction pathways on a microscopic level. K+ channels utilize multi-ion conduction mechanisms, and the three-dimensional structure also shows several ions present in the channel. Here we report results from molecular dynamics free energy perturbation calculations that both establish the nature of the multiple ion conduction mechanism and yield the correct ion selectivity of the channel. By evaluating the energetics of all relevant occupancy states of the selectivity filter, we find that the favoured conduction pathway involves transitions only between two main states with a free difference of about 5 kcal mol(-1). Other putative permeation pathways can be excluded because they would involve states that are too high in energy.

Bacterial Proteins↗
PubMed2000

Prediction of a ligand-induced conformational change in the catalytic core of Cdc25A.

The cell cycle control phosphatases Cdc25 are dual specificity phosphatases that dephosphorylate both phosphothreonine and phosphotyrosine residues on their substrate proteins. The determination of the apo-protein structure of Cdc25A revealed that this enzyme has a completely different fold compared to all other phosphatases crystallised to date. The conformation of the active site residues does not seem very suitable for catalysis in this unliganded structure. We have studied some structural features of the Cdc25A apo-structure and a modelled Cdc25A-ligand complex by molecular dynamics simulations. The simulations predict a conformational change in the peptide backbone of the complex, which is not observed in the apo-structure. This ligand-induced conformational change yields a structure that is similar to other protein tyrosine phosphatase-ligand complexes that have been crystallised. The change in conformation takes place in the position between a serine and a glutamic acid residue in the phosphate binding loop. We suggest that this type of conformational change is an important molecular switch in the catalytic process.

Binding Sites↗
PubMed1999

Long-range electrostatic effects on peptide folding.

The reversible folding/unfolding of a short peptide in solution is studied by molecular dynamics simulations. The effects of long-range electrostatic interactions are examined and found to be important both for the equilibrium between folded and unfolded states and the dynamics of the folding process. The neglect of long-range electrostatics leads to an increased population of unfolded states and increased structural fluctuations. When such interactions are taken into account, the peptide unfolds and folds to the experimentally determined structure several times during a 25 ns trajectory, with approximately equal populations of folded and unfolded states in the neighborhood of its proposed melting temperature. The effect of using spherical boundary conditions rather than periodic ones does not appear to have any major effect on the folding dynamics.

Computer Simulation↗
PubMed1999

Mechanism of substrate dephosphorylation in low Mr protein tyrosine phosphatase.

Substrate dephosphorylation by the low molecular weight protein tyrosine phosphatases proceeds via nucleophilic substitution at the phosphorous atom yielding a cysteinyl phosphate intermediate. However, several questions regarding the exact reaction mechanism remain unanswered. Starting from the crystal structure of the enzyme we study the energetics of this reaction, using the empirical valence bond method in combination with molecular dynamics and free energy perturbation simulations. The free energy profiles of two mechanisms corresponding to different protonation states of the reacting groups are examined along stepwise and concerted pathways. The activation barriers calculated relative to the enzyme-substrate complex are very similar for both monoanionic and dianionic substrates, but taking the substrate binding step into account shows that hydrolysis of monoanionic substrates is strongly favored by the enzyme, because a dianionic substrate will not bind when the reacting cysteine is ionized. The calculated activation barrier for dephosphorylation of monoanionic phenyl phosphate according to this novel mechanism is 14 kcal mol(-1), which is in good agreement with experimental data. Proteins 1999;36:370-379.

Catalytic Domain↗
PubMed1999

Computational modeling of the rate limiting step in low molecular weight protein tyrosine phosphatase.

Hydrolysis of the phosphoenzyme intermediate is the second and rate limiting step of the reaction catalyzed by the protein tyrosine phosphatases (PTPs). The cysteinyl phosphate thioester bond is cleaved by nucleophilic displacement where an active site water molecule attacks the phosphorus atom. Starting from the crystal structure of the low molecular weight PTP, we study the energetics of this reaction utilizing the empirical valence bond method in combination with molecular dynamics and free energy perturbation simulations. The reactions of the wild-type as well as the D129A and C17S mutants are modeled. For the D129A mutant, which lacks the general acid/base residue Asp-129, an alternative reaction mechanism is proposed. The calculated activation barriers are in all cases in good agreement with experimental reaction rates. The present results together with earlier computational and experimental work now provide a detailed picture of the complete reaction mechanism in many PTPs. The key role played by the structurally invariant signature motif in stabilizing a double negative charge is reflected by its control of the energetics of both transition states and the reaction intermediate.

Animals↗
PubMed1999

Energetics of the proposed rate-determining step of the glyoxalase I reaction.

The proposed rate-limiting step of the reaction catalyzed by glyoxalase I is the proton abstraction from the C1 carbon atom of the substrate by a glutamate residue, resulting in a high-energy enolate intermediate. This proton transfer reaction was modelled using molecular dynamics and free energy perturbation simulations, with the empirical valence bond method describing the potential energy surface of the system. The calculated rate constant for the reaction is approximately 300-1500 s(-1) with Zn2+, Mg2+ or Ca2+ bound to the active site, which agrees well with observed kinetics of the enzyme. Furthermore, the results imply that the origin of the catalytic rate enhancement is mainly associated with enolate stabilization by the metal ion.

Catalytic Domain↗
PubMed1999

Mechanistic alternatives in phosphate monoester hydrolysis: what conclusions can be drawn from available experimental data?

Phosphate monoester hydrolysis reactions in enzymes and solution are often discussed in terms of whether the reaction pathway is associative or dissociative. Although experimental results for solution reactions have usually been considered as evidence for the second alternative, a closer thermodynamic analysis of observed linear free energy relationships shows that experimental information is consistent with the associative, concerted and dissociative alternatives.

Animals↗
PubMed1998

Computation of affinity and selectivity: binding of 2,4-diaminopteridine and 2,4-diaminoquinazoline inhibitors to dihydrofolate reductases.

Binding energy calculations for complexes of mutant and wild-type human dihydrofolate reductases with 2,4-diaminopteridine and 2,4-diaminoquinazoline inhibitors are reported. Quantitative insight into binding energetics of these molecules is obtained from calculations based on force field energy evaluation and thermal sampling by molecular dynamics simulations. The calculated affinity of methotrexate for wild-type and mutant enzymes is reasonably well reproduced. Truncation of the methotrexate glutamate tail results in a loss of affinity by several orders of magnitude. No major difference in binding strength is predicted between the pteridines and the quinazolìnes, while the N-methyl group present in methotrexate appears to confer significantly stronger binding. The recent improvement, which is used here, of our linear interaction energy method for binding affinity prediction, as well as problems with treating charged and flexible ligands are discussed. This approach should be suitable in a drug discovery context for prediction of binding energies of new inhibitors prior to their synthesis, when some information about the binding mode is available.

Computer Simulation↗
PubMed1998

Ligand binding affinity prediction by linear interaction energy methods.

A recent method for estimating ligand binding affinities is extended. This method employs averages of interaction potential energy terms from molecular dynamics simulations or other thermal conformational sampling techniques. Incorporation of systematic deviations from electrostatic linear response, derived from free energy perturbation studies, into the absolute binding free energy expression significantly enhances the accuracy of the approach. This type of method may be useful for computational prediction of ligand binding strengths, e.g., in drug design applications.

Binding Sites↗
PubMed1997

Cyclic HIV-1 protease inhibitors derived from mannitol: synthesis, inhibitory potencies, and computational predictions of binding affinities.

Ten C2-symmetric cyclic urea and sulfamide derivatives have been synthesized from L-mannonic gamma-lactone and D-mannitol. The results of experimental measurement of their inhibitory potencies against HIV-1 protease were compared to calculated free energies of binding derived from molecular dynamics (MD) simulations. The compounds were selected, firstly, to enable elucidation of the role of stereochemistry for binding affinity (1a-d) and, secondly, to allow evaluation of the effects of variation in the link to the P1 and P1' phenyl groups on affinity (1a and 2-5). Thirdly, compounds with hydrogen bond-accepting or-donating groups attached to the phenyl groups in the P2 and P2' side chains (6 and 7) were selected. Binding free energies were estimated by a linear response method, whose predictive power for estimating binding affinities from MD simulations was demonstrated.

Anti-HIV Agents↗
PubMed1997

Energetics of nucleophile activation in a protein tyrosine phosphatase.

The nucleophilic attack by cysteine 12 in the low-molecular-weight protein tyrosine phosphatase is believed to be carried out by the thiolate anion form of this residue. We here study the energetics of proton transfer between the thiol group of cysteine 12 and a substrate phosphate oxygen atom, to examine the effects of the enzymic environment on the stability of the thiolate nucleophile. This is done by molecular dynamics and free energy perturbation simulations, utilizing the empirical valence bond method to describe the potential surface of the system. The calculations show that the protein environment significantly stabilizes the thiolate ion, thereby setting the stage for the nucleophilic attack. We compare these results with those from further simulations of a mutant enzyme, and demonstrate the importance of serine 19 in thiolate stabilization.

Computer Simulation↗
PubMed1996

Computer simulation of the triosephosphate isomerase catalyzed reaction.

A major challenge for theoretical simulation methods is the calculation of enzymic reaction rates directly from the three-dimensional protein structure together with some idea of the chemical reaction mechanism. Here, we report the evaluation of a complete free energy profile for all the elementary steps of the triosephosphate isomerase catalyzed reaction using such an approach. The results are compatible with available experimental data and also suggest which of the possible reaction intermediates is kinetically observable. In addition to previously identified catalytic residues, the simulations show that a crystallographically observed active site water molecule plays an important role during catalysis and an intersubunit interaction that could explain the low activity of the monomeric enzyme is also observed. The calculations clearly demonstrate the important catalytic effects associated with stabilization of charged high energy intermediates and reduction of reorganization energy, which are likely to be general principles of enzyme catalyzed charge transfer and separation reactions.

Catalysis↗
PubMed1995

Sugar recognition by a glucose/galactose receptor. Evaluation of binding energetics from molecular dynamics simulations.

A new theoretical method for free energy calculations is used to compute the absolute binding constants for beta-D-glucose and methyl-beta-D-galactoside to the periplasmic glucose/galactose receptor from Salmonella typhimurium. The computer simulation results agree well with available experimental data and make it possible to assess the sources of both the high affinity as well as the specificity for glucose. It was found that the major contribution to the binding energy comes from electrostatic interactions and particularly hydrogen bonds of the charge-dipole type. We also predict the structure of the complex with methyl-galactoside as this has not yet been experimentally determined.

Computer Simulation↗
PubMed1995

Estimation of binding free energies for HIV proteinase inhibitors by molecular dynamics simulations.

Absolute binding free energies for three inhibitors of HIV-1 proteinase were estimated from molecular dynamics simulations by a recently reported linear approximation procedure. The results were in fairly good agreement with experimental binding data. Two of the inhibitors were very similar and, for comparison, their relative free energies of binding were also calculated by free energy perturbation methods, giving virtually the same result. Effects of cut-off radii and charge states of the protein model were examined. The effects of pH on binding of one of the inhibitors were predicted.

Computer Simulation↗