Search PubMed⌕ Search

Biomedical subjects

J Aqvist

Publications and source records attributed to J Aqvist.

At least 37 records · Page 2Linked to original sources
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↗
PubMed1994

A new method for predicting binding affinity in computer-aided drug design.

A new semi-empirical method for calculating free energies of binding from molecular dynamics (MD) simulations is presented. It is based on standard thermodynamic cycles and on a linear approximation of polar and non-polar free energy contributions from the corresponding MD averages. The method is tested on a set of endothiapepsin inhibitors and found to give accurate results both for absolute as well as relative free energies.

Aspartic Acid Endopeptidases↗
PubMed1992

Computer simulation of the initial proton transfer step in human carbonic anhydrase I.

The initial water proteolysis step in the proton transfer "half-reaction" of human carbonic anhydrase I is simulated using the empirical valence bond method in combination with free energy perturbation molecular dynamics calculations. A free energy profile for the enzyme catalysed reaction and the corresponding pKa associated with ionization of the zinc-bound water is calculated. The obtained pKa value of 7 to 8 appears to be in good agreement with experimental observations and the calculated rate constant for this step is also compatible with kinetic data. The simulations clearly emphasize the important electrostatic effect associated with the catalytic zinc ion.

Carbonic Anhydrases↗
PubMed1992

A molecular model for the retinol binding protein-transthyretin complex.

A three-dimensional model for the complex between human serum retinol binding protein and transthyretin (formerly named prealbumin) is presented. The model was obtained by interactive rigid-body computer graphics docking and the characterization of the molecular surfaces in terms of fractal dimension. Available experimental data, as well as results from molecular dynamics calculations, support the proposed model.

Computer Graphics↗
PubMed1992

Computer simulations of enzymatic reactions: examination of linear free-energy relationships and quantum-mechanical corrections in the initial proton-transfer step of carbonic anhydrase.

Computer simulation approaches can provide a powerful tool for correlating the structure of enzymes with their catalytic activity. One of the most effective ways of simulating enzymatic reactions is provided by the empirical valence bond method. The general applicability of this method has been demonstrated in several enzymatic reactions and it is reexamined here in a study of the initial proton-transfer step in the catalytic reaction of carbonic anhydrase. The simulations produce a rate constant which is in agreement with the observed kinetic data and emphasizes the importance of the electrostatic effect associated with the catalytic zinc ion. The calculations are also used to examine the validity of linear free-energy relationships (LFERs) in enzyme catalysis and to evaluate quantum-mechanical corrections of the calculated rate constant. It is found that LFERs are valid in the present case and it is argued that this reflects the fact that the protein responds linearly to the development of electrostatic forces during the reaction. It is concluded that the present approach can be used to augment experimental studies in establishing the general validity of LFERs. It is noted, however, that such relationships are much more valid for transitions between different resonance structures than for transitions between reactants and product states.

Carbonic Anhydrases↗
PubMed1991

Dipoles localized at helix termini of proteins stabilize charges.

The charge-stabilization effect associated with alpha-helices in proteins has been reexamined by microscopic calculations without any a priori assumptions about the dielectric constant of the protein. The calculations reproduce the observed charge stabilization effect of a helix in two well-defined test cases: a histidine residue situated at the C-terminal end of a helix in barnase and the sulfate ligand located near the N-terminal end of a helix of the sulfate-binding protein. They also show that the effective dielectric constant for helix-charge interactions is much larger than previously assumed and that the stabilizing effect of the helix is not associated with the helix macrodipole but rather with a few localized dipoles confined mostly to the first turn of the helix. It is predicted that mutations at one end of the helix should have very small effects on the stabilization of charges at the opposite terminus. It is pointed out that the relatively short-ranged effect of the helix is essentially similar to other cases in which localized dipoles play key roles in electrostatic stabilization.

Bacterial Proteins↗
PubMed1989

Calculations of free energy profiles for the staphylococcal nuclease catalyzed reaction.

Calculations of the free energy profile for the first two (rate-limiting) steps of the staphylococcal nuclease catalyzed reaction are reported. The calculations are based on the empirical valence bond method in combination with free energy perturbation molecular dynamics simulations. The calculated activation free energy is in good agreement with experimental kinetic data, and the catalytic effect of the enzyme is reproduced without any arbitrary adjustment of parameters. The enormous reduction of the activation barrier (relative to the reference reaction in water) appears to be largely associated with the strong electrostatic effect of the Ca2+ ion and the two arginine residues in the active site. This favorable electrostatic environment reduces the cost of the general-base catalysis step by almost 15 kcal/mol (by stabilizing the OH- nucleophile) and then stabilizes the developing negative charge on the 5'-phosphate group in the second step of the reaction by about 19 kcal/mol. The basic features of the originally postulated enzyme mechanism (Cotton et al., 1979) are found to be compatible with the observed activation free energy. However, the proposed modification of the mechanism (Sepersu et al., 1987), in which Arg 87 interacts only with the pentacoordinated transition state, is supported by the simulations. Further calculations on the D21E mutant also give results in good agreement with kinetic data.

Binding Sites↗