Dataset Information

Matrix metalloproteinase 2 inhibition: combined quantum mechanics and molecular mechanics studies of the inhibition mechanism of (4-phenoxyphenylsulfonyl)methylthiirane and its oxirane analogue.

ABSTRACT: The inhibition mechanism of matrix metalloproteinase 2 (MMP2) by the selective inhibitor (4-phenoxyphenylsulfonyl)methylthiirane (SB-3CT) and its oxirane analogue is investigated computationally. The inhibition mechanism involves C-H deprotonation with concomitant opening of the three-membered heterocycle. SB-3CT was docked into the active site of MMP2, followed by molecular dynamics simulation to prepare the complex for combined quantum mechanics and molecular mechanics (QM/MM) calculations. QM/MM calculations with B3LYP/6-311+G(d,p) for the QM part and the AMBER force field for the MM part were used to examine the reaction of these two inhibitors in the active site of MMP2. The calculations show that the reaction barrier for transformation of SB-3CT is 1.6 kcal/mol lower than its oxirane analogue, and the ring-opening reaction energy of SB-3CT is 8.0 kcal/mol more exothermic than that of its oxirane analogue. Calculations also show that protonation of the ring-opened product by water is thermodynamically much more favorable for the alkoxide obtained from the oxirane than for the thiolate obtained from the thiirane. A six-step partial charge fitting procedure is introduced for the QM/MM calculations to update atomic partial charges of the quantum mechanics region and to ensure consistent electrostatic energies for reactants, transition states, and products.

SUBMITTER: Tao P

PROVIDER: S-EPMC2818682 | biostudies-literature | 2009 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Matrix metalloproteinase 2 inhibition: combined quantum mechanics and molecular mechanics studies of the inhibition mechanism of (4-phenoxyphenylsulfonyl)methylthiirane and its oxirane analogue.

Tao Peng P Fisher Jed F JF Shi Qicun Q Vreven Thom T Mobashery Shahriar S Schlegel H Bernhard HB

Biochemistry 20091001 41

The inhibition mechanism of matrix metalloproteinase 2 (MMP2) by the selective inhibitor (4-phenoxyphenylsulfonyl)methylthiirane (SB-3CT) and its oxirane analogue is investigated computationally. The inhibition mechanism involves C-H deprotonation with concomitant opening of the three-membered heterocycle. SB-3CT was docked into the active site of MMP2, followed by molecular dynamics simulation to prepare the complex for combined quantum mechanics and molecular mechanics (QM/MM) calculations. QM ...[more]

PMID: 19754151

Similar Datasets

Project description:SB-3CT, (4-phenoxyphenylsulfonyl)methylthiirane, is a potent, mechanism-based inhibitor of the gelatinase sub-class of the matrix metalloproteinase (MMP) family of zinc proteases. The gelatinase MMPs are unusual in that there are several examples where both enantiomers of a racemic inhibitor have comparable inhibitory abilities. SB-3CT is one such example. Here, the inhibition mechanism of the MMP2 gelatinase by the (S)-SB-3CT enantiomer and its oxirane analogue is examined computationally, and compared to the mechanism of (R)-SB-3CT. Inhibition of MMP2 by (R)-SB-3CT was shown previously to involve enzyme-catalyzed C-H deprotonation adjacent to the sulfone, with concomitant opening by ?-elimination of the sulfur of the three-membered thiirane ring. Similarly to the R enantiomer, (S)-SB-3CT was docked into the active site of MMP2, followed by molecular dynamics simulation to prepare the complex for combined quantum mechanics and molecular mechanics (QM/MM) calculations. QM/MM calculations with B3LYP/6-311+G(d,p) for the QM part (46 atoms) and the AMBER force field for the MM part were used to compare the reaction of (S)-SB-3CT and its oxirane analogue in the active site of MMP2 (9208 atoms). These calculations show that the barrier for the proton abstraction coupled ring opening reaction of (S)-SB-3CT in the MMP2 active site is 4.4 kcal/mol lower than its oxirane analogue, and the ring opening reaction energy of (S)-SB-3CT is only 1.6 kcal/mol less exothermic than its oxirane analogue. Calculations also show that the protonation of the ring-opened products by water is thermodynamically much more favorable for the alkoxide obtained from the oxirane, than for the thiolate obtained from the thiirane. In contrast to (R)-SB-3CT and the R-oxirane analogue, the double bonds of the ring-opened products of (S)-SB-3CT and its S-oxirane analogue have the cis-configuration. Vibrational frequency and intrinsic reaction path calculations on a reduced size QM/MM model (2747 atoms) provide additional insight into the mechanism. These calculations yield 5.9 and 6.7 for the deuterium kinetic isotope effect for C-H bond cleavage in the transition state for the R and S enantiomers of SB-3CT, in good agreement with the experimental results.

Project description:(4-Phenoxyphenylsulfonyl)methylthiirane (SB-3CT) is the selective inhibitor of matrix metalloproteinase 2 (MMP2). The inhibition mechanism of MMP2 by SB-3CT involves C-H deprotonation with concomitant opening of the three-membered heterocycle. In this study, the energetics of the deprotonation-induced ring-opening of (4-phenoxyphenylsulfinyl)methylthiirane, the sulfoxide analogue of SB-3CT, are examined computationally using DFT and QM/MM calculations. A model system, 2-(methylsulfinylmethyl)thiirane, is used to study the stereoelectronic and conformational effects of reaction barriers in methanol. For the model system in methanol solution (using the polarizable continuum model), the reaction barriers range from 17 to 23 kcal/mol with significant stereoelectronic effects. However, the lowest barriers of the (R,R) and (S,R) diastereomers are similar. Two diastereomers of the sulfoxide analogue of SB-3CT are studied in the active site of MMP2 by QM/MM methods with an accurate partial charge fitting procedure. The ring-opening reactions of these two diastereomers have similar reaction energetics. Both are exothermic from the reactant to the ring-opening product (thiolate). The protonation of the thiolate by a water molecule is endothermic in both cases. However, the deprotonation/ring-opening barriers in the MMP2 active site using QM/MM methods for the (R,R) and (S,R) inhibitions are quite different (23.3 and 28.5 kcal/mol, respectively). The TSs identified in QM/MM calculations were confirmed by vibrational frequency analysis and following the reaction path. The (R,R) diastereomer has a hydrogen bond between the sulfoxide oxygen and the backbone NH of Leu191, while the (S,R) has a hydrogen bond between the sulfoxide oxygen and a water molecule. The dissimilar strengths of these hydrogen bonds as well as minor differences in the TS structures contribute to the difference between the barriers. Compared to SB-3CT, both diastereomers of the sulfoxide analogue have higher reaction barriers and have less exothermic reaction energies. This agrees well with the experiments, where SB-3CT is a more effective inhibitor of MMP2 than its sulfoxide analogue.

Project description:The significance of conical intersections in photophysics, photochemistry, and photodissociation of polyatomic molecules in gas phase has been demonstrated by numerous experimental and theoretical studies. Optimization of conical intersections of small- and medium-size molecules in gas phase has currently become a routine optimization process, as it has been implemented in many electronic structure packages. However, optimization of conical intersections of small- and medium-size molecules in solution or macromolecules remains inefficient, even poorly defined, due to large number of degrees of freedom and costly evaluations of gradient difference and nonadiabatic coupling vectors. In this work, based on the sequential quantum mechanics and molecular mechanics (QM/MM) and QM/MM-minimum free energy path methods, we have designed two conical intersection optimization methods for small- and medium-size molecules in solution or macromolecules. The first one is sequential QM conical intersection optimization and MM minimization for potential energy surfaces; the second one is sequential QM conical intersection optimization and MM sampling for potential of mean force surfaces, i.e., free energy surfaces. In such methods, the region where electronic structures change remarkably is placed into the QM subsystem, while the rest of the system is placed into the MM subsystem; thus, dimensionalities of gradient difference and nonadiabatic coupling vectors are decreased due to the relatively small QM subsystem. Furthermore, in comparison with the concurrent optimization scheme, sequential QM conical intersection optimization and MM minimization or sampling reduce the number of evaluations of gradient difference and nonadiabatic coupling vectors because these vectors need to be calculated only when the QM subsystem moves, independent of the MM minimization or sampling. Taken together, costly evaluations of gradient difference and nonadiabatic coupling vectors in solution or macromolecules can be reduced significantly. Test optimizations of conical intersections of cyclopropanone and acetaldehyde in aqueous solution have been carried out successfully.

Dataset Information

Matrix metalloproteinase 2 inhibition: combined quantum mechanics and molecular mechanics studies of the inhibition mechanism of (4-phenoxyphenylsulfonyl)methylthiirane and its oxirane analogue.

Publications

Matrix metalloproteinase 2 inhibition: combined quantum mechanics and molecular mechanics studies of the inhibition mechanism of (4-phenoxyphenylsulfonyl)methylthiirane and its oxirane analogue.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets