Configuration Flipping in Distal Pocket of Multidrug Transporter MexB Impacts the Efflux Inhibitory Mechanism

MexAB-OprM efflux pumps, found in Pseudomonas aeruginosa, play a major role in drug resistance by extruding out drugs and antibiotic molecules from cells. Inhibitors are used to cease the potency of the efflux pumps. In this study, in-silico models are used to investigate the nature of the binding pocket of the MexAB-OprM efflux pump. First, we have performed classical molecular dynamics (MD) simulations to shed light on different aspects of protein–inhibitor interaction in the binding pocket of the pump. Using classical MD simulations, quantum mechanics/molecular mechanics (QM/MM), and various types of analyses, it is found that D13-9001 has a higher binding affinity towards the binding pocket compared to D1 and D2; the results are in sync with the experimental dat. Two stable configurations of D13-9001 are discovered inside the distal pocket which could be one of the primary reasons for the greater efficacy of D13-9001. The free energy barrier upon changing one state to another is calculated by employing umbrella sampling method. Finally, F178 is mutated to have the complete picture as it contributes significantly to the binding energy irrespective of the three inhibitors. Our results may help to design a new generation of inhibitors for such an efflux pump.     

Quantum Mechanics/Molecular Mechanics Studies of the Mechanism of Cysteine Proteases Inhibition by Dipeptidyl Nitroalkenes

In this work a computational study of the mechanism of inhibition of cruzain, rhodesain, and cathepsin L cysteine proteases by the dipeptidyl nitroalkene Cbz-Phe-Ala-CH=CH-NO₂ has been carried out by means of molecular dynamics simulations with hybrid QM/MM potentials. The free-energy surfaces confirmed that the inhibition takes place by the formation of a covalent bond between the protein and the β-carbon atom of the inhibitor. According to the results, the tested inhibitor should be a much more efficient inhibitor of cruzain than of rhodesain, and little activity would be expected against cathepsin L, in total correspondence with the available experimental data. The origin of these differences may lie in the different stabilizing electrostatic interactions established between the inhibitor and the residues of the active site and S2 pocket of these enzymes. These results may be useful for the rational design of new dipeptidyl nitroalkenes with higher and more selective inhibitory activity against cysteine proteases.     

The Role of Hydrogen Bond in Catalytic Triad of Serine Proteases

In order to investigate the origin of catalytic power for serine proteases, the role of the hydrogen bond in the catalytic triad was studied in the proteolysis process of the peptides chymotrypsin inhibitor 2 (CI2), MCTI-A, and a hexapeptide (SUB), respectively. We first calculated the free energy profile of the proton transfer between His and Asp residues of the catalytic triad in the enzyme-substrate state and transition state by employing QM/MM molecular dynamics simulations. The results show that a low-barrier hydrogen bond (LBHB) only forms in the transition state of the acylation of CI2, while it is a normal hydrogen bond in the acylation of MCTI-A or SUB. In addition, the change of the hydrogen bond strength is much larger in CI2 and SUB systems than in MCTI-A system, which decreases the acylation energy barrier significantly for CI2 and SUB. Clearly, a LBHB formed in the transition state region helps accelerate the acylation reaction. But to our surprise, a normal hydrogen bond can also help to decrease the energy barrier. The key to reducing the reaction barrier is the increment of hydrogen bond strength in the transition state state, whether it is a LBHB or not. Our studies cast new light on the role of the hydrogen bond in the catalytic triad, and help to understand the catalytic triad of serine proteases.     

Implementation of an adaptive umbrella sampling method for the calculation of multidimensional potential of mean force of chemical reactions in solution

We describe the implementation of an adaptive umbrella sampling method, making use of the weighted histogram analysis method, for computing multidimensional potential of mean force for chemical reaction in solution. The approach is illustrated by investigating the effect of aqueous solution on the free energy surface for the proton transfer reaction of [H(3)N-H-NH(3)](+) using a combined quantum mechanical and molecular mechanical AM1/TIP3P potential.     

Mechanism of chemical reactions in the active site of aspartate N-acetyltransferase NAT8L revealed by molecular modeling

The results of a computational study of the synthesis of a key brain metabolite, N-acetyl-l-aspartate, catalyzed by aspartate N-acetyltransferase, encoded by the NAT8L gene, are reported. The reaction Gibbs energy profiles were computed using molecular dynamics simulations with interaction potentials estimated on-the-fly by the quantum mechanics/molecular mechanics QM(PBE0/6-31G**)/MM(CHARMM) approach. The revealed reaction mechanism includes four elementary steps with corresponding activation energies not exceeding 14 kcal mol^?1     

Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations

Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

The O to S substitution in urea brings inhibition activity against thiocyanate dehydrogenase

1. Download : Download high-res image (150KB)
2. Download : Download full-size image

     

Is inhibition process better described with MD(QM/MM) simulations? The case of urokinase type plasminogen activator inhibitors

Urokinase plasminogen activator (uPA) is an enzyme involved in cancer growth and metastasis. Therefore, the design of inhibitors of uPA is of high therapeutic value, and several chemical families have been explored, even if none has still emerged, emphasizing the need of a rationalized approach. This work represents a complete computational study of uPA complexed with five inhibitors, which present weak similarities. Molecular dynamics simulations in explicit solvent were conducted, and structural analyses, along with molecular mechanics (MM)/Poisson-Boltzmann surface area free energies estimations, yield precious structure-activity relationships of these inhibitors. Besides, we realized supplemental QM/MM computations that improved drastically the quality of our models providing original information on the hydrogen bonds and charge transfer effects, which are, most often, neglected in other studies. We suggest that these simulations and analyses could be reproduced for other systems involving protein/ligand molecular recognitions.     

Two Tags in One Probe: Combining Fluorescence- and Biotin-based Detection of the Trypanosomal Cysteine Protease Rhodesain

Rhodesain is the major cysteine protease of the protozoan parasite Trypanosoma brucei and a therapeutic target for sleeping sickness, a fatal neglected tropical disease. We designed, synthesized and characterized a bimodal activity-based probe that binds to and inactivates rhodesain. This probe exhibited an irreversible mode of action and extraordinary potency for the target protease with a k_inac/K_i value of 37,000 M⁻¹s⁻¹. Two reporter tags, a fluorescent coumarin moiety and a biotin affinity label, were incorporated into the probe and enabled highly sensitive detection of rhodesain in a complex proteome by in-gel fluorescence and on-blot chemiluminescence. Furthermore, the probe was employed for microseparation and quantification of rhodesain and for inhibitor screening using a competition assay. The developed bimodal rhodesain probe represents a new proteomic tool for studying Trypanosoma pathobiochemistry and antitrypanosomal drug discovery.     

Minimum MD simulation length required to achieve reliable results in free energy perturbation calculations: Case study of relative binding free energies of fructose‐1,6‐bisphosphatase inhibitors

In an attempt to establish the criteria for the length of simulation to achieve the desired convergence of free energy calculations, two studies were carried out on chosen complexes of FBPase‐AMP mimics. Calculations were performed for varied length of simulations and for different starting configurations using both conventional‐ and QM/MM‐FEP methods. The results demonstrate that for small perturbations, 1248 ps simulation time could be regarded a reasonable yardstick to achieve convergence of the results. As the simulation time is extended, the errors associated with free energy calculations also gradually tapers off. Moreover, when starting the simulation from different initial configurations of the systems, the results are not changed significantly, when performed for 1248 ps. This study carried on FBPase‐AMP mimics corroborates well with our previous successful demonstration of requirement of simulation time for solvation studies, both by conventional and ab initio FEP. The establishment of aforementioned criteria of simulation length serves a useful benchmark in drug design efforts using FEP methodologies, to draw a meaningful and unequivocal conclusion. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Molecular dynamics simulations of the detoxification of paraoxon catalyzed by phosphotriesterase

Combined QM(PM3)/MM molecular dynamics simulations together with QM(DFT)/MM optimizations for key configurations have been performed to elucidate the enzymatic catalysis mechanism on the detoxification of paraoxon by phosphotriesterase (PTE). In the simulations, the PM3 parameters for the phosphorous atom were reoptimized. The equilibrated configuration of the enzyme/substrate complex showed that paraoxon can strongly bind to the more solvent‐exposed metal ion Zn_β, but the free energy profile along the binding path demonstrated that the binding is thermodynamically unfavorable. This explains why the crystal structures of PTE with substrate analogues often exhibit long distances between the phosphoral oxygen and Zn_β. The subsequent S_N2 reaction plays the key role in the whole process, but controversies exist over the identity of the nucleophilic species, which could be either a hydroxide ion terminally coordinated to Zn_α or the μ‐hydroxo bridge between the α‐ and β‐metals. Our simulations supported the latter and showed that the rate‐limiting step is the distortion of the bound paraoxon to approach the bridging hydroxide. After this preparation step, the bridging hydroxide ion attacks the phosphorous center and replaces the diethyl phosphate with a low barrier. Thus, a plausible way to engineer PTE with enhanced catalytic activity is to stabilize the deformed paraoxon. Conformational analyses indicate that Trp131 is the closest residue to the phosphoryl oxygen, and mutations to Arg or Gln or even Lys, which can shorten the hydrogen bond distance with the phosphoryl oxygen, could potentially lead to a mutant with enhanced activity for the detoxification of organophosphates. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD)

Extensive combined quantum mechanical (B3LYP/6‐31G*) and molecular mechanical (QM/MM) molecular dynamics simulations have been performed to elucidate the hydrolytic deamination mechanism of cytosine to uracil catalyzed by the yeast cytosine deaminase (yCD). Though cytosine has no direct binding to the zinc center, it reacts with the water molecule coordinated to zinc, and the adjacent conserved Glu64 serves as a general acid/base to shuttle protons from water to cytosine. The overall reaction consists of several proton‐transfer processes and nucleophilic attacks. A tetrahedral intermediate adduct of cytosine and water binding to zinc is identified and similar to the crystal structure of yCD with the inhibitor 2‐pyrimidinone. The rate‐determining step with the barrier of 18.0 kcal/mol in the whole catalytic cycle occurs in the process of uracil departure where the proton transfer from water to Glu64 and nucleophilic attack of the resulting hydroxide anion to C2 of the uracil ring occurs synchronously. © 2016 Wiley Periodicals, Inc.     

Theoretical analysis of the kinetic isotope effect on carboxylation in RubisCO

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO) fixes atmospheric carbon dioxide into bioavailable sugar molecules. It is also well known that a kinetic isotope effect (KIE; CO₂ carbon atoms) accompanies the carboxylation process. To describe the reaction and the KIE α, two different types of molecular dynamics (MD) simulations (ab initio MD and classical MD) have been performed with an Own N-layered Integrated molecular Orbitals and molecular Mechanics (ONIOM)-hybrid model. A channel structure for CO₂ transport has been observed during the MD simulation in RubisCO, and assuming the reaction path from the inlet to the product through the coordinate complex with Mg²⁺, simulations have been performed on several molecular configuration models fixing several distances between CO₂ and ribulose-1,5-bisphosphate along the channel. Free energy analysis and diffusion coefficient analysis have been evaluated for different phases of the process. It is confirmed that the isotopic fractionation effect for CO₂ containing either ¹³C or ¹²C would appear through the transiting path in the channel structure identified in RubisCO. The estimated isotope fractionation constant was quite close to the experimental value.     

Two Sides of Quantum-Based Modeling of Enzyme-Catalyzed Reactions: Mechanistic and Electronic Structure Aspects of the Hydrolysis by Glutamate Carboxypeptidase

We report the results of a computational study of the hydrolysis reaction mechanism of N-acetyl-l-aspartyl-l-glutamate (NAAG) catalyzed by glutamate carboxypeptidase II. Analysis of both mechanistic and electronic structure aspects of this multistep reaction is in the focus of this work. In these simulations, model systems are constructed using the relevant crystal structure of the mutated inactive enzyme. After selection of reaction coordinates, the Gibbs energy profiles of elementary steps of the reaction are computed using molecular dynamics simulations with ab initio type QM/MM potentials (QM/MM MD). Energies and forces in the large QM subsystem are estimated in the DFT(PBE0-D3/6-31G**) approximation. The established mechanism includes four elementary steps with the activation energy barriers not exceeding 7 kcal/mol. The models explain the role of point mutations in the enzyme observed in the experimental kinetic studies; namely, the Tyr552Ile substitution disturbs the “oxyanion hole”, and the Glu424Gln replacement increases the distance of the nucleophilic attack. Both issues diminish the substrate activation in the enzyme active site. To quantify the substrate activation, we apply the QTAIM-based approaches and the NBO analysis of dynamic features of the corresponding enzyme-substrate complexes. Analysis of the 2D Laplacian of electron density maps allows one to define structures with the electron density deconcentration on the substrate carbon atom, i.e., at the electrophilic site of reactants. The similar electronic structure element in the NBO approach is a lone vacancy on the carbonyl carbon atom in the reactive species. The electronic structure patterns revealed in the NBO and QTAIM-based analyses consistently clarify the reactivity issues in this system.     

Preorganization and reorganization as related factors in enzyme catalysis: the chorismate mutase case

In this paper a deeper insight into the chorismate-to prephenate-rearrangement, catalyzed by Bacillus subtilis chorismate mutase, is provided by means of a combination of statistical quantum mechanics/molecular mechanics simulation methods and hybrid potential energy surface exploration techniques. The main aim of this work is to present an estimation of the preorganization and reorganization terms of the enzyme catalytic rate enhancement. To analyze the first of these, we have studied different conformational equilibria of chorismate in aqueous solution and in the enzyme active site. Our conclusion is that chorismate mutase preferentially binds the reactive conformer of the substrate--that presenting a structure similar to the transition state of the reaction to be catalyzed--with shorter distances between the carbon atoms to be bonded and more diaxial character. With respect to the reorganization effect, an energy decomposition analysis of the potential energies of the reactive reactant and of the reaction transition state in aqueous solution and in the enzyme shows that the enzyme structure is better adapted to the transition structure. This means not only a more negative electrostatic interaction energy with the transition state but also a low enzyme deformation contribution to the energy barrier. Our calculations reveal that the structure of the enzyme is responsible for stabilizing the transition state structure of the reaction, with concomitant selection of the reactive form of the reactants. This is, the same enzymatic pattern that stabilizes the transition structure also promotes those reactant structures closer to the transition structure (i.e., the reactive reactants). In fact, both reorganization and preorganization effects have to be considered as the two faces of the same coin, having a common origin in the effect of the enzyme structure on the energy surface of the substrate.     

Interactions of Solvated Electrons with Nucleobases: The Effect of Base Pairing

We have investigated the effect of base pairing on the electron attachment to nucleobases in bulk water, taking the guanine−cytosine (GC) base pair as a test case. The presence of the complementary base reinforces the stabilization effect provided by water and preferentially stabilizes the anion by hydrogen bonding. The electron attachment in bulk-solvated GC happens through a doorway mechanism, where the initial electron attached state is water bound, and it subsequently gets converted to a GC bound state. The additional electron in the final GC bound state is localized on the cytosine, similar to that in the gas phase. The transfer of the electron from the initial water-bound state to the final GC bound state happens due to the mixing of electronic and nuclear degrees of freedom and takes place at a picosecond time scale.     

锌酶的计算模拟:挑战与最新进展

锌酶在人体中分布非常广泛,种类繁多,是当前最受关注的金属酶之一。由于在锌配位结构上的多样性以及Zn²⁺饱和的d轨道带来的“光谱寂静”性,导致许多实验研究手段受限。计算模拟在锌酶的研究中发挥着越来越重要的作用,已经成为不可或缺的研究工具。现代量子化学计算模拟方法,特别是被视为研究生物大分子体系非常有效的QM/MM组合方法,目前已经被广泛应用于探讨复杂多变的锌配位结构以及锌酶催化反应机理。通过在QM/MM水平下开展的分子动力学模拟,可以揭示锌酶体系中结构与功能间的相互关系。此外,分子力场方法在锌酶研究中同样发挥了不可替代的作用,由于传统力场普遍无法正确描述锌配位结构,因此,锌酶分子力场的开发具有迫切的现实意义。本文总结了近年来锌酶计算模拟领域的最新进展,提出了锌酶计算研究中还有待解决的一些问题。     

On the effect of a variation of the force field, spatial boundary condition and size of the QM region in QM/MM MD simulations

During the past years, the use of combined quantum-classical, QM/MM, methods for the study of complex biomolecular processes, such as enzymatic reactions and photocycles, has increased considerably. The quality of the results obtained from QM/MM calculations is largely dependent on five aspects to be considered when setting up a molecular model: the QM Hamiltonian, the MM Hamiltonian or force field, the boundary and coupling between the QM and MM regions, the size of the QM region and the boundary condition for the MM region. In this study, we systematically investigate the influence of a variation of the molecular mechanics force field and the size of the QM region in QM/MM MD simulations on properties of the photoactive part of the blue light photoreceptor protein AppA. For comparison, we additionally performed classical MD simulations and studied the effect of a variation of the type of spatial boundary condition. The classical boundary conditions and the force field used in a QM/MM MD simulation are shown to have non-neglegible effects upon the structural and energetic properties of the protein which makes it advisable to minimize computational artifacts in QM/MM MD simulations by application of periodic boundary conditions and a thermodynamically calibrated force field. A comparison of the structural and energetic properties of MD simulations starting from two alternative, different X-ray structures for the blue light utilizing flavin protein in its dark state indicates a slight preference of the two force fields used for the so-called Anderson structure over the Jung structure.