Molecular dynamics, crystallography and mutagenesis studies on the substrate gating mechanism of prolyl oligopeptidase

Altered prolyl oligopeptidase (PREP) activity is found in many common neurological and other genetic disorders, and in some cases PREP inhibition may be a promising treatment. The active site of PREP resides in an internal cavity; in addition to the direct interaction between active site and substrate or inhibitor, the pathway to reach the active site (the gating mechanism) must be understood for more rational inhibitor design and understanding PREP function. The gating mechanism of PREP has been investigated through molecular dynamics (MD) simulation combined with crystallographic and mutagenesis studies. The MD results indicate the inter-domain loop structure, comprised of 3 loops at residues, 189-209 (loop A), 577-608 (loop B), and 636-646 (loop C) (porcine PREP numbering), are important components of the gating mechanism. The results from enzyme kinetics of PREP variants also support this hypothesis: When loop A is (1) locked to loop B through a disulphide bridge, all enzyme activity is halted, (2) nicked, enzyme activity is increased, and (3) removed, enzyme activity is only reduced. Limited proteolysis study also supports the hypothesis of a loop A driven gating mechanism. The MD results show a stable network of H-bonds that hold the two protein domains together. Crystallographic study indicates that a set of known PREP inhibitors inhabit a common binding conformation, and this H-bond network is not significantly altered. Thus the domain separation, seen to occur in lower taxa, is not involved in the gating mechanism for mammalian PREP. In two of the MD simulations we observed a conformational change that involved the breaking of the H-bond network holding loops A and B together. We also found that this network was more stable when the active site was occupied, thus decreasing the likelihood of this transition.     

Steered molecular dynamics simulations reveal important mechanisms in reversible monoamine oxidase B inhibition

The monotopic membrane protein monoamine oxidase B (MAO B) is an important drug target for Parkinson's disease. In order to design more specific, and thereby more effective, inhibitors for this enzyme, it is necessary to determine what factors govern inhibitor specificity and the inhibitor binding process, including the roles of the lipid bilayer, the active site loop, and several key residues within the binding pocket. Atomistic molecular dynamics simulations of MAO B either embedded in a lipid bilayer or free in solution have been performed. The simulations suggest that the bilayer controls the availability of the active site cavity by regulating the degree of fluctuation in two key loops that form the greater part of the active site entrance (residues 85-110 and 155-165). In turn, the enzyme itself causes local thinning and a decrease in area per lipid of the surrounding bilayer environment. Additional MD simulations of MAO B in complex with seven different reversible inhibitors followed by nonequilibrium steered MD simulations of the inhibitor unbinding have also been performed. The simulations demonstrate that the average energy of interaction between inhibitor and MAO B residues during inhibitor egress is an effective indicator of inhibitor strength and is also useful for identifying key residues that govern inhibitor specificity. These data provide researchers with valuable tools for designing effective MAO B inhibitors as well as outline a method that can be translated to the study of other enzyme-inhibitor complexes.     

Molecular insight into pseudolysin inhibition using the MM-PBSA and LIE methods

Pseudolysin, the extracellullar elastase of Pseudomonas aeruginosa (EC: 3.4.24.26) plays an important role in the pathogenesis of P. aeruginosa infections. In the present study, molecular dynamics simulations and theoretical affinity predictions were used to gain molecular insight into pseudolysin inhibition. Four low molecular weight inhibitors were docked at their putative binding sites and molecular dynamics (MD) simulations were performed for 5.0 ns, and the free energy of binding was calculated by the linear interaction energy method. The number and the contact surface area of stabilizing hydrophobic, aromatic, and hydrogen bonding interactions appears to reflect the affinity differences between the inhibitors. The proteinaceous inhibitor, Streptomyces metalloproteinase inhibitor (SMPI) was docked in three different binding positions and MD simulations were performed for 3.0 ns. The MD trajectories were used for molecular mechanics-Poisson-Boltzmann surface area analysis of the three binding positions. Computational alanine scanning of the average pseudolysin-SMPI complexes after MD revealed residues at the pseudolysin-SMPI interface giving the main contribution to the free energy of binding. The calculations indicated that SMPI interacts with pseudolysin via the rigid active site loop, but that also contact sites outside this loop contribute significantly to the free energy of association.     

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

The tunnel region at triosephosphate isomerase (TIM)’s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM’s catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand’s effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop’s dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.     

Protonation state and substrate binding to B2 metallo-beta-lactamase CphA from Aeromonas hydrofila

The zinc enzymes metallo beta-lactamases counteract the beneficial action of beta-lactam antibiotics against bacterial infections, by hydrolyzing their beta-lactam rings. To understand structure/function relationships on a representative member of this class, the B2 M beta L CphA, we have investigated the H-bond pattern at the Zn enzymatic active site and substrate binding mode by molecular simulation methods. Extensive QM calculations at the DFT-BLYP level on eleven models of the protein active site, along with MD simulations of the protein in aqueous solution, allow us to propose two plausible protonation states for the unbound enzyme, which are probably in equilibrium. Docking procedures along with MD simulations and QM calculations suggest that in the complex between the enzyme and its substrate (biapenem), the latter is stable in only one of the two protonation states, in addition it exhibits two different binding modes, of which only one agrees with previous proposals. In both cases, the substrate is polarized as in aqueous solution. We conclude that addressing mechanistic issues on this class of enzymes requires a careful procedure to assign protonation states and substrate docking modes.     

Structural determinants of ligand imprinting: a molecular dynamics simulation study of subtilisin in aqueous and apolar solvents

The phenomenon known as "ligand imprinting" or "ligand-induced enzyme memory" was first reported in 1988, when Russell and Klibanov observed that lyophilizing subtilisin in the presence of competitive inhibitors (that were subsequently removed) could significantly enhance its activity in an apolar solvent. (Russell and Klibanov, J Biol Chem 1988;263:11624-11626). They further observed that this enhancement did not occur when similar assays were carried out in water. Herein, we shed light on the molecular determinants of ligand imprinting using a molecular dynamics (MD) approach. To simulate the effect of placing an enzyme in the presence of a ligand before its lyophilization, an inhibitor was docked in the active site of subtilisin and 20 ns MD simulations in water were performed. The ligand was then removed and the resulting structure was used for subsequent MD runs using hexane and water as solvents. As a control, the same simulation setup was applied using the structure of subtilisin in the absence of the inhibitor. We observed that the ligand maintains the active site in an open conformation and that this configuration is retained after the removal of the inhibitor, when the simulations are carried out in hexane. In agreement with experimental findings, the structural configuration induced by the ligand is lost when the simulations take place in water. Our analysis of fluctuations indicates that this behavior is a result of the decreased flexibility displayed by enzymes in an apolar solvent, relatively to the aqueous situation.     

Conformational mobility of active and E-64-inhibited actinidin

Background

Actinidin, a protease from kiwifruit, belongs to the C1 family of cysteine proteases. Cysteine proteases were found to be involved in many disease states and are valid therapeutic targets. Actinidin has a wide pH activity range and wide substrate specificity, which makes it a good model system for studying enzyme–substrate interactions.

Methods

The influence of inhibitor (E-64) binding on the conformation of actinidin was examined by 2D PAGE, circular dichroism (CD) spectroscopy, hydrophobic ligand binding assay, and molecular dynamics simulations.

Results

Significant differences were observed in electrophoretic mobility of proteolytically active and E-64-inhibited actinidin. CD spectrometry and hydrophobic ligand binding assay revealed a difference in conformation between active and inhibited actinidin. Molecular dynamics simulations showed that a loop defined by amino-acid residues 88–104 had greater conformational mobility in the inhibited enzyme than in the active one. During MD simulations, the covalently bound inhibitor was found to change its conformation from extended to folded, with the guanidino moiety approaching the carboxylate.

Conclusions

Conformational mobility of actinidin changes upon binding of the inhibitor, leading to a sequence of events that enables water and ions to protrude into a newly formed cavity of the inhibited enzyme. Drastic conformational mobility of E-64, a common inhibitor of cysteine proteases found in many crystal structures stored in PDB, was also observed.

General significance

The analysis of structural changes which occur upon binding of an inhibitor to a cysteine protease provides a valuable starting point for the future design of therapeutic agents.     

Successful molecular dynamics simulation of the zinc-bound farnesyltransferase using the cationic dummy atom approach

Farnesyltransferase (FT) inhibitors can suppress tumor cell proliferation without substantially interfering with normal cell growth, thus holding promise for cancer treatment. A structure-based approach to the design of improved FT inhibitors relies on knowledge of the conformational flexibility of the zinc-containing active site of FT. Although several X-ray structures of FT have been reported, detailed information regarding the active site conformational flexibility of the enzyme is still not available. Molecular dynamics (MD) simulations of FT can offer the requisite information, but have not been applied due to a lack of effective methods for simulating the four-ligand coordination of zinc in proteins. Here, we report in detail the problems that occurred in the conventional MD simulations of the zinc-bound FT and a solution to these problems by employing a simple method that uses cationic dummy atoms to impose orientational requirement for zinc ligands. A successful 1.0 ns (1.0 fs time step) MD simulation of zinc-bound FT suggests that nine conserved residues (Asn127alpha, Gln162alpha, Asn165alpha, Gln195alpha, His248beta, Lys294beta, Leu295beta, Lys353beta, and Ser357beta) in the active site of mammalian FT are relatively mobile. Some of these residues might be involved in the ligand-induced active site conformational rearrangement upon binding and deserve attention in screening and design of improved FT inhibitors for cancer chemotherapy.     

Brownian and Essential Dynamics Studies of the HIV-1 Integrase Catalytic Domain

Abstract

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computional approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a “gating”-like motion with respect to the active site.     

Brownian and essential dynamics studies of the HIV-1 integrase catalytic domain

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computational approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.     

Loop dynamics and ligand binding kinetics in the reaction catalyzed by the Yersinia protein tyrosine phosphatase

The Yersinia protein tyrosine phosphatase (YopH) contains a loop of ten amino acids (the WPD loop) that covers the entrance of the active site of the enzyme during substrate binding. In this work the substrate mimicking competitive inhibitor p-nitrocatechol sulfate (PNC) is used as a probe of the active site. The dynamics of the WPD loop was determined by subjecting an equilibrated system containing YopH, PNC, and YopH bound to PNC to a laser induced temperature jump, and subsequently following the change in equilibrium due to the perturbation. Using this methodology the dynamics associated with substrate binding in YopH have been determined. These results indicate that substrate binding is coupled to the WPD loop motion, and WPD loop dynamics occur in the sub-millisecond time scale. The significance of these dynamic results is interpreted in terms of the catalytic cycle of the enzyme.     

MD simulations reveal alternate conformations of the oxyanion hole in the Zika virus NS2B/NS3 protease

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple 100-ns conventional MD simulations in both the presence and absence of inhibitors. The S1, S2, and S3 pockets are stable as well. However, in two of eight simulations, the A132-G133 peptide bond in the binding pocket of S1' spontaneously flips to form a 3₁₀-helix that corresponds to the inactive conformation of the oxyanion hole, and then maintains this conformation until the end of the 100-ns conventional MD simulations without inversion of the flip. This conformational change affects the S1' pocket in ZIKV NS2B/NS3 protease active site, which is important for small molecule binding. The simulation results provide evidence at the atomic level that the inactive conformation of the oxyanion hole is more favored energetically when no specific interactions are formed between substrate/inhibitor and oxyanion hole residues. Interestingly, however, transition between the active and inactive conformation of the oxyanion hole can be observed by boosting the valley potential in accelerated MD simulations. This supports a proposed induced-fit mechanism of ZIKV NS2B/NS3 protease from computational methods and provides useful direction to enhance inhibitor binding predictions in structure-based drug design.     

Nanosecond dynamics of acetylcholinesterase near the active center gorge

To delineate the role of peptide backbone flexibility and rapid molecular motion in acetylcholinesterase catalysis and inhibitor association, we investigated the decay of fluorescence anisotropy at three sites of fluorescein conjugation to cysteine-substitution mutants of the enzyme. One cysteine was placed in a loop at the peripheral site near the rim of the active center gorge (H287C); a second was in a helical region outside of the active center gorge (T249C); a third was at the tip of a small, flexible omega loop well separated from the gorge (A262C). Mutation and fluorophore conjugation did not appreciably alter catalytic or inhibitor binding parameters of the enzyme. The results show that each site examined was associated with a high degree of segmental motion; however, the A262C and H287C sites were significantly more flexible than the T249C site. Association of the active center inhibitor, tacrine, and the peripheral site peptide inhibitor, fasciculin, had no effect on the anisotropy decay of fluorophores at positions 249 and 262. Fasciculin, but not tacrine, on the other hand, dramatically altered the decay profile of the fluorophore at the 287 position, in a manner consistent with fasciculin reducing the segmental motion of the peptide chain in this local region. The results suggest that the motions of residues near the active center gorge and across from the Cys(69)-Cys(96) omega loop are uncoupled and that ligand binding at the active center or the peripheral site does not influence acetylcholinesterase conformational dynamics globally, but induces primarily domain localized decreases in flexibility proximal to the bound ligand.     

Investigation for antimicrobial resistance-modulating activity of diethyl malate and 1-methyl malate against beta-lactamase class A from Bacillus licheniformis by molecular dynamics,in vitro and in vivo studies

Resistance to antibiotics in bacteria, is one of the major problems of mankind. Each year, a large number of patients due to infection, lose their lives. One of the main mechanisms of antibiotic resistance is beta-lactamase secretion. This enzyme hydrolyzes the amide bond of a lactam ring in beta-lactam antibiotics. Bacillus licheniformis is a mesophilic gram-positive bacterium, which has a high potential to produce beta-lactamase class A. In this study, the inhibitory effects of some malate analogous were studied by in vitro and in vivo studies. In addition, the effects of inhibitor binding on beta-lactamase were studied using MD simulations. Our results showed that diethyl malate and 1-methyl malate can decrease the MIC value of benzyl penicillin by sixteen and eight-fold, respectively. Data derived from in vitro studies revealed that decrease in MIC values is correlated with beta-lactamase inhibition. Molecular docking studies predicted the binding mode of inhibitors with the beta-lactamase active site. The structural analysis from MD simulations exhibits that binding of citrate and diethyl malate causes earlier equilibrium of beta-lactamase. After binding, the fluctuation of Ser 70 is also decreased. Based on our data, diethyl malate can be used to design the potent inhibitor against beta-lactamase class A.     

Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

Binding affinity of hydroxamate inhibitors of matrix metalloproteinase-2

Here we report molecular dynamics (MD) and free energy perturbation (FEP) simulations applied to hydroxamate-matrix metalloproteinase-2 (MMP-2) complex systems. We have developed some new force field parameters for the hydroxamate functional group that were not included in the AMBER94 force field but were necessary in our simulations. For the representation of the active zinc center, a bonded model was adopted in which restrained electrostatic potential fitting (RESP) charges were used as the electrostatic representation of this model. Using the resulted bonded model, FEP simulations predict the relative binding free energy in good agreement with the experimental value. By analyzing the molecular dynamics (MD) trajectories of the two complex systems, we can provide an explanation of why one of the two inhibitors is favored over the other. The results provide a chemical insight into the interactions between inhibitor and enzyme, and can indicate changes in the inhibitor that would enhance inhibitor–enzyme interactions.Figure The scheme of the binding site     

Serotype-selective, small-molecule inhibitors of the zinc endopeptidase of botulinum neurotoxin serotype A

Botulinum neurotoxin serotype A (BoNTA) is one of the most toxic substances known. Currently, there is no antidote to BoNTA. Small molecules identified from high-throughput screening reportedly inhibit the endopeptidase--the zinc-bound, catalytic domain of BoNTA--at a drug concentration of 20 microM. However, optimization of these inhibitors is hampered by challenges including the computational evaluation of the ability of a zinc ligand to compete for coordination with nearby residues in the active site of BoNTA. No improved inhibitor of the endopeptidase has been reported. This article reports the development of a serotype-selective, small-molecule inhibitor of BoNTA with a K(i) of 12 microM. This inhibitor was designed to coordinate the zinc ion embedded in the active site of the enzyme for affinity and to interact with a species-specific residue in the active site for selectivity. It is the most potent small-molecule inhibitor of BoNTA reported to date. The results suggest that multiple molecular dynamics simulations using the cationic dummy atom approach are useful to structure-based design of zinc protease inhibitors.     

Structural and dynamical correlations in PfHGXPRT oligomers: A molecular dynamics simulation study

PfHGXPRT is a key enzyme involved in purine nucleotide salvage pathway of the malarial parasite, Plasmodium falciparum. Atomistic molecular dynamics simulations have been performed on two types of PfHGXPRT dimers (D1 and D3) and its tetramer in their apo and ligand-bound states. A significant event in the catalytic cycle is the dynamics of a gate that provides access for the ligand molecules to the reaction center. The gate is formed by loops II and IV, the former being the most flexible. Large amplitude conformational changes have been observed in active site loop II. Upon complete occupancy of the active site, loop II gets stabilized due to specific interactions between its residues and the ligand molecules. Remote loop, X, is seen to be less fluxional in the D3 dimer than in D1 which is rationalized as due to the greater number of inter-subunit contacts in the former. The presence of ligand molecules in subunits of the tetramer further reduces the flexibility of loop X epitomizing a communication between this region and the active sites in the tetramer. These observations are in accordance with the outcomes of several experimental investigations. Participation of loop X in the oligomerization process has also been discerned. Between the two types of dimers in solution, D1 tetramerizes readily and thus would not be present as free dimers. We conjecture an equilibrium to exist between D3 and the tetramer in solution; upon binding of the ligand molecules to the D3 dimer, this equilibrium shifts toward the tetramer.     

Variations in the Binding Pocket of an Inhibitor of the Bacterial Division Protein FtsZ across Genotypes and Species

The recent increase in antibiotic resistance in pathogenic bacteria calls for new approaches to drug-target selection and drug development. Targeting the mechanisms of action of proteins involved in bacterial cell division bypasses problems associated with increasingly ineffective variants of older antibiotics; to this end, the essential bacterial cytoskeletal protein FtsZ is a promising target. Recent work on its allosteric inhibitor, PC190723, revealed in vitro activity on Staphylococcus aureus FtsZ and in vivo antimicrobial activities. However, the mechanism of drug action and its effect on FtsZ in other bacterial species are unclear. Here, we examine the structural environment of the PC190723 binding pocket using PocketFEATURE, a statistical method that scores the similarity between pairs of small-molecule binding sites based on 3D structure information about the local microenvironment, and molecular dynamics (MD) simulations. We observed that species and nucleotide-binding state have significant impacts on the structural properties of the binding site, with substantially disparate microenvironments for bacterial species not from the Staphylococcus genus. Based on PocketFEATURE analysis of MD simulations of S. aureus FtsZ bound to GTP or with mutations that are known to confer PC190723 resistance, we predict that PC190723 strongly prefers to bind Staphylococcus FtsZ in the nucleotide-bound state. Furthermore, MD simulations of an FtsZ dimer indicated that polymerization may enhance PC190723 binding. Taken together, our results demonstrate that a drug-binding pocket can vary significantly across species, genetic perturbations, and in different polymerization states, yielding important information for the further development of FtsZ inhibitors.