Biomedical subjects

Johannes C Hermann

Publications and source records attributed to Johannes C Hermann.

8 recordsLinked to original sources

PubMed2006

Predicting substrates by docking high-energy intermediates to enzyme structures.

With the emergence of sequences and even structures for proteins of unknown function, structure-based prediction of enzyme activity has become a pragmatic as well as an interesting question. Here we investigate a method to predict substrates for enzymes of known structure by docking high-energy intermediate forms of the potential substrates. A database of such high-energy transition-state analogues was created from the KEGG metabolites. To reduce the number of possible reactions to consider, we restricted ourselves to enzymes of the amidohydrolase superfamily. We docked each metabolite into seven different amidohydrolases in both the ground-state and the high-energy intermediate forms. Docking the high-energy intermediates improved the discrimination between decoys and substrates significantly over the corresponding standard ground-state database, both by enrichment of the true substrates and by geometric fidelity. To test this method prospectively, we attempted to predict the enantioselectivity of a set of chiral substrates for phosphotriesterase, for both wild-type and mutant forms of this enzyme. The stereoselectivity ratios of the six enzymes considered for those four substrate enantiomer pairs differed over a range of 10- to 10,000-fold and underwent 20 switches in stereoselectivities for favored enantiomers, compared to the wild type. The docking of the high-energy intermediates correctly predicted the stereoselectivities for 18 of the 20 substrate/enzyme combinations when compared to subsequent experimental synthesis and testing. The possible applications of this approach to other enzymes are considered.

Algorithms↗

PubMed2006

Resolution of chiral phosphate, phosphonate, and phosphinate esters by an enantioselective enzyme library.

An array of 16 enantiomeric pairs of chiral phosphate, phosphonate, and phosphinate esters was used to establish the breadth of the stereoselective discrimination inherent within the bacterial phosphotriesterase and 15 mutant enzymes. For each substrate, the leaving group was 4-hydroxyacetophenone while the other two groups attached to the phosphorus core consisted of an asymmetric mixture of methyl, methoxy, ethyl, ethoxy, isopropoxy, phenyl, phenoxy, cyclohexyl, and cyclohexoxy substituents. For the wild-type enzyme, the relative rates of hydrolysis for the two enantiomers ranged from 3 to 5.4 x 10(5). Various combinations of site-specific mutations within the active site were used to create modified enzymes with alterations in their enantioselective properties. For the single-site mutant enzyme, G60A, the stereoselectivity is enhanced relative to that of the wild-type enzyme by 1-3 orders of magnitude. Additional mutants were obtained where the stereoselectivity is inverted relative to the wild-type enzyme for 13 of the 16 pairs of enantiomers tested for this investigation. The most dramatic example was obtained for the hydrolysis of 4-acetylphenyl methyl phenyl phosphate. The G60A mutant preferentially hydrolyzes the SP-enantiomer by a factor of 3.7 x 10(5). The I106G/F132G/H257Y mutant preferentially hydrolyzes the RP-enantiomer by a factor of 9.7 x 10(2). This represents an enantioselective discrimination of 3.6 x 10(8) between these two mutants, with a total of only four amino acid changes. The rate differential between the two enantiomers for any given mutant enzyme is postulated to be governed by the degree of nonproductive binding within the enzyme active site and stabilization of the transition state. This hypothesis is supported by computational docking of the high-energy, pentavalent form of the substrates to modeled structures of the mutant enzyme; the energies of the docked transition-state analogues qualitatively capture the enantiomeric preferences of the various mutants for the different substrates. These results demonstrate that the catalytic properties of the wild-type phosphotriesterase can be exploited for the kinetic resolution of a wide range of phosphate, phosphonate, and phosphinate esters and that the active site of this enzyme is remarkably amenable to structural perturbations via amino acid substitution.

Amino Acid Substitution↗

PubMed2005

Molecular mechanisms of antibiotic resistance: QM/MM modelling of deacylation in a class A beta-lactamase.

Modelling of the first step of the deacylation reaction of benzylpenicillin in the E. coli TEM1 beta-lactamase (with B3LYP/6-31G + (d)//AM1-CHARMM22 quantum mechanics/molecular mechanics methods) shows that a mechanism in which Glu166 acts as the base to deprotonate a conserved water molecule is both energetically and structurally consistent with experimental data; the results may assist the design of new antibiotics and beta-lactamase inhibitors.

Acylation↗

PubMed2005

QM/MM modelling of oleamide hydrolysis in fatty acid amide hydrolase (FAAH) reveals a new mechanism of nucleophile activation.

Fatty acid amide hydrolase (FAAH), a promising target for the treatment of several central and peripheral nervous system disorders, such as anxiety, pain and hypertension, has an unusual catalytic site, and its mechanism has been uncertain; hybrid quantum mechanics/molecular mechanics (QM/MM) calculations reveal a new mechanism of nucleophile activation (involving a Lys-Ser-Ser catalytic triad), with potentially crucial insights for the design of potent and selective inhibitors.

Amidohydrolases↗

PubMed2005

Mechanisms of antibiotic resistance: QM/MM modeling of the acylation reaction of a class A beta-lactamase with benzylpenicillin.

Understanding the mechanisms by which beta-lactamases destroy beta-lactam antibiotics is potentially vital in developing effective therapies to overcome bacterial antibiotic resistance. Class A beta-lactamases are the most important and common type of these enzymes. A key process in the reaction mechanism of class A beta-lactamases is the acylation of the active site serine by the antibiotic. We have modeled the complete mechanism of acylation with benzylpenicillin, using a combined quantum mechanical and molecular mechanical (QM/MM) method (B3LYP/6-31G+(d)//AM1-CHARMM22). All active site residues directly involved in the reaction, and the substrate, were treated at the QM level, with reaction energies calculated at the hybrid density functional (B3LYP/6-31+Gd) level. Structures and interactions with the protein were modeled by the AM1-CHARMM22 QM/MM approach. Alternative reaction coordinates and mechanisms have been tested by calculating a number of potential energy surfaces for each step of the acylation mechanism. The results support a mechanism in which Glu166 acts as the general base. Glu166 deprotonates an intervening conserved water molecule, which in turn activates Ser70 for nucleophilic attack on the antibiotic. This formation of the tetrahedral intermediate is calculated to have the highest barrier of the chemical steps in acylation. Subsequently, the acylenzyme is formed with Ser130 as the proton donor to the antibiotic thiazolidine ring, and Lys73 as a proton shuttle residue. The presented mechanism is both structurally and energetically consistent with experimental data. The QM/MM energy barrier (B3LYP/ 6-31G+(d)//AM1-CHARMM22) for the enzymatic reaction of 9 kcal mol(-1) is consistent with the experimental activation energy of about 12 kcal mol(-1). The effects of essential catalytic residues have been investigated by decomposition analysis. The results demonstrate the importance of the "oxyanion hole" in stabilizing the transition state and the tetrahedral intermediate. In addition, Asn132 and a number of charged residues in the active site have been identified as being central to the stabilizing effect of the enzyme. These results will be potentially useful in the development of stable beta-lactam antibiotics and for the design of new inhibitors.

Acylation↗

PubMed2004

A combined QM/MM approach to protein--ligand interactions: polarization effects of the HIV-1 protease on selected high affinity inhibitors.

HIV-1 protease inhibitors are one of the two widely used therapeutic agents for the treatment of HIV-infected patients. The investigation of HIV-1 protease-inhibitor interactions can provide further insight for developing new compounds that are still required due to the growing problem of drug resistance. To this end, a combined QM/MM approach was used to determine electrostatic and polarization interactions on three high affinity inhibitors, nelfinavir, mozenavir, and tipranavir. The present computational results show that explicit treatment of the polarization effect is particularly important since it can contribute as much as one-third of the total electrostatic interaction energy. Further, an amino acid decomposition analysis was applied to determine contributions of individual residues to the enzyme--inhibitor interactions. It was found that the 4-hydroxy-dihydropyrone substructure of tipranavir is especially suited for extended charge delocalization by interacting with the catalytic aspartates and isoleucines of the HIV-1 protease. The calculated electron density difference maps reaffirm and provide a means of visualizing these results.

HIV Protease↗

PubMed2004

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction.

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.

Catalysis↗

PubMed2003

Identification of Glu166 as the general base in the acylation reaction of class A beta-lactamases through QM/MM modeling.

Bacterial class A beta-lactamases are responsible for the most known resistance against beta-lactam antibiotics. With the continuing rise in antibiotic resistance, improved knowledge of the mechanisms of action of these enzymes is needed in the development of effective therapeutic agents and strategies. The mechanism of the deacylation step in class A beta-lactamases is well accepted. In contrast, the mechanism of the acylation step has been uncertain, with several conflicting proposals put forward. We have modeled the acylation step in a class A beta-lactamase, using a combined quantum mechanics/molecular mechanics approach. The results provide an atomic level description of the reaction and show that Glu166 acts as the general base in the reaction, deprotonating Ser70 via an intervening water molecule. Ser70 acts as the nucleophile for attack on the lactam ring in a concerted reaction. The results do not rule out the importance of Lys73 in catalysis, in agreement with experimental data.

Acylation↗