SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62Linked to original sources

PubMed1997

QXP: powerful, rapid computer algorithms for structure-based drug design.

New methods for docking, template fitting and building pseudo-receptors are described. Full conformational searches are carried out for flexible cyclic and acyclic molecules. QXP (quick explore) search algorithms are derived from the method of Monte Carlo perturbation with energy minimization in Cartesian space. An additional fast search step is introduced between the initial perturbation and energy minimization. The fast search produces approximate low-energy structures, which are likely to minimize to a low energy. For template fitting, QXP uses a superposition force field which automatically assigns short-range attractive forces to similar atoms in different molecules. The docking algorithms were evaluated using X-ray data for 12 protein-ligand complexes. The ligands had up to 24 rotatable bonds and ranged from highly polar to mostly nonpolar. Docking searches of the randomly disordered ligands gave rms differences between the lowest energy docked structure and the energy-minimized X-ray structure, of less than 0.76 A for 10 of the ligands. For all the ligands, the rms difference between the energy-minimized X-ray structure and the closest docked structure was less than 0.4 A, when parts of one of the molecules which are in the solvent were excluded from the rms calculation. Template fitting was tested using four ACE inhibitors. Three ACE templates have been previously published. A single run using QXP generated a series of templates which contained examples of each of the three. A pseudo-receptor, complementary to an ACE template, was built out of small molecules, such as pyrrole, cyclopentanone and propane. When individually energy minimized in the pseudo-receptor, each of the four ACE inhibitors moved with an rms of less than 0.25 A. After random perturbation, the inhibitors were docked into the pseudo-receptor. Each lowest energy docked structure matched the energy-minimized geometry with an rms of less than 0.08 A. Thus, the pseudo-receptor shows steric and chemical complementarity to all four molecules. The QXP program is reliable, easy to use and sufficiently rapid for routine application in structure-based drug design.

Algorithms↗

PubMed2006

Function and structure of the molybdenum cofactor carrier protein from Chlamydomonas reinhardtii.

The molybdenum cofactor (Moco) forms the catalytic site in all eukaryotic molybdenum enzymes and is synthesized by a multistep biosynthetic pathway. The mechanism of transfer, storage, and insertion of Moco into the appropriate apo-enzyme is poorly understood. In Chlamydomonas reinhardtii, a Moco carrier protein (MCP) has been identified and characterized recently. Here we show biochemical evidence that MCP binds Moco as well as the tungstate-substituted form of the cofactor (Wco) with high affinity, whereas molybdopterin, the ultimate cofactor precursor, is not bound. This binding selectivity points to a specific metal-mediated interaction with MCP, which protects Moco and Wco from oxidation with t((1/2)) of 24 and 96 h, respectively. UV-visible spectroscopy showed defined absorption bands at 393, 470, and 570 nm pointing to ene-diothiolate and protein side-chain charge transfer bonds with molybdenum. We have determined the crystal structure of MCP at 1.6 Angstrom resolution using seleno-methionated and native protein. The monomer constitutes a Rossmann fold with two homodimers forming a symmetrical tetramer in solution. Based on conserved surface residues, charge distribution, shape, in silico docking studies, structural comparisons, and identification of an anionbinding site, a prominent surface depression was proposed as a Moco-binding site, which was confirmed by structure-guided mutagenesis coupled to substrate binding studies.

Amino Acid Sequence↗

PubMed1998

Molecular modeling studies of the DNA-topoisomerase I ternary cleavable complex with camptothecin.

The present studies provide a three-dimensional model for the postulated ternary cleavable complex of topoisomerase I (top1), DNA, and camptothecin (CPT). Molecular simulations were done using the AMBER force field. The results suggest that a ternary cleavable complex might be stabilized by several hydrogen bonds in the binding site. In this proposed "drug-stacking" model, CPT is pseudointercalated in the top1-linked DNA cleavage site and interacts with the protein near its catalytic tyrosine through hydrogen bonding and stacking. The structural model is consistent with the following experimental observations: (i) the N3 position of the 5' terminal purine of the cleaved DNA strand is readily alkylated by 7-chloromethyl 10,11-methylenedioxy CPT; (ii) CPT generally tolerates substituents at positions 7, 9, and 10 but is inactivated by additions at position 12; (iii) 10,11-methylenedioxy (MDO) CPT is much more potent than 10,11-dimethoxy (DMO) CPT; (iv) the lactone portion of CPT is essential for top1 inhibitory activity; (v) 20S derivatives of CPT are much more potent than the 20R analogues; (vi) a catalytic tyrosine hydroxyl in top1 covalently links to the 3' terminal base, T, of the cleaved DNA strand; and (vii) top1 mutation Asn722Ser leads to CPT resistance. A total of 18 camptothecin derivatives with different DNA cleavage potencies were docked into the hypothetical cleavable complex binding site to test and refine the model. These studies provide insight into a possible mechanism of top1 inhibition by CPT derivatives and suggest rational approaches for the design of new CPT derivatives.

Amino Acid Substitution↗

PubMed2004

Type I collagen N-telopeptides adopt an ordered structure when docked to their helix receptor during fibrillogenesis.

The in vitro rate and specificity of fibrillogenesis in type I collagen depends on the integrity of the amino (N)-telopeptide domain. In vivo an intact N-telopeptide domain is also required for normal fibril assembly. Although Chou-Fasman predictions and NMR studies suggested that a type I beta-turn could be induced in alpha1(I) N-telopeptide chains, computer modeling did not identify ordered structures. Nevertheless, X-ray analysis and electron tomography studies have shown that the N-telopeptide is in one of the most highly ordered fibril domains. This study was undertaken to determine if the docking of the N-telopeptide to its helix receptor domain could induce the telopeptides to take up a specific conformation. With use of molecular modeling suite of programs, a (Gly-Pro-Pro)(n) triple-helical structure was built on the basis of high-resolution X-ray crystallographic coordinates and then replaced with the actual bovine collagen residues 924-938, the triple-helical alpha1(I)-N-telopeptide-receptor sequences. Energy minimization produced a modified triple-helical conformation. The bovine alpha1(I) N-telopeptide sequence was similarly minimized and docked to this receptor. The docking induced an ordered conformation with a stabilizing hydrogen bond in the N-telopeptide and, importantly, a reciprocal reordering of the triple-helical conformation in the binding domain. This docked structure placed Lys residues in both telopeptide and helix in the correct locations for cross-link formation. The modeling has been extended to the three-chain N-telopeptide domain and finally to the construction of the Hulmes-Miller quasi-hexagonal packing structure. Each N-telopeptide domain can form linkages with two adjacent, aligned helix receptor domains. The telopeptides and the order of staggering of the three chains in the helix play crucial roles in the packing and intrafibrillar cross-linking patterns and the relative azimuthal orientations of adjacent molecules in the fibril. The models confirm the high order in the N-telopeptide 4D overlap zone.

Amino Acid Sequence↗

PubMed2005

Receptor flexibility in de novo ligand design and docking.

One of the major problems in computational drug design is incorporation of the intrinsic flexibility of protein binding sites. This is particularly crucial in ligand binding events, when induced fit can lead to protein structure rearrangements. As a consequence of the huge conformational space available to protein structures, receptor flexibility is rarely considered in ligand design procedures. In this work, we present an algorithm for integrating protein binding-site flexibility into de novo ligand design and docking processes. The approach allows dynamic rearrangement of amino acid side chains during the docking and design simulations. The impact of protein conformational flexibility is investigated in the docking of highly active inhibitors in the binding sites of acetylcholinesterase and human collagenase (matrix metalloproteinase-1) and in the design of ligands in the S1' pocket of MMP-1. The results of corresponding simulations for both rigid and flexible binding sites are compared in order to gauge the influence of receptor flexibility in drug discovery protocols.

Acetylcholinesterase↗

PubMed2005

Searching for protein-protein interaction sites and docking by the methods of molecular dynamics, grid scoring, and the pairwise interaction potential of amino acid residues.

In CAPRI Rounds 1 and 2, we assumed that because there are many ionic charges that weaken electrostatic interaction forces in living cells, the hydrophobic interaction force might be important entropically. As a result of Rounds 1 and 2, the predictions for binding sites and geometric centers were acceptable, but those of the binding axes were poor, because only the largest benzene cluster was used for generating the initial docking structures. These were generated by fitting of benzene clusters formed on the surface of receptor and ligand. In CAPRI Rounds 3-5, the grid-scoring sum on the protein-protein interaction surface and the pairwise potential of the amino acid residues, which were indicated as coming easily into the protein-protein interaction regions, were used as the calculation methods, along with the smaller benzene clusters that participated in benzene cluster fitting. Good predicted models were obtained for Targets 11 and 12. When the modeled receptor proteins were superimposed on the experimental structures, the smallest ligand root-mean-square deviation (RMSD) values corresponding to the RMSD between the model and experimental structures were 6.2 A and 7.3 A, respectively.

Algorithms↗

PubMed2002

Dynamic docking and electron transfer between myoglobin and cytochrome b(5).

The interaction of trypsin-digested bovine cytochrome b(5) (cyt b(5)) with horse heart myoglobin (Mb) and the interprotein electron transfer (ET) between these redox partners have been studied to gain better understanding of ET processes between weakly bound protein partners. The bimolecular rate constant ( k(2)) for photo-induced ET between zinc-substituted Mb (ZnMb) and cyt b(5) decreases with increasing ionic strength, consistent with the predominantly electrostatic character of this complex. The formation of a protein-protein complex has been confirmed and the binding affinities of metMb and ZnMb for cyt b(5) have been measured by two techniques: (1)H NMR titrations at pH 6.0 give binding constants of K(a) approximately (1.0+/-0.1)x10(3) M(-1) for metMb and K(a) approximately (0.75+/-0.1)x10(3) M(-1) for ZnMb; isothermal calorimetry gives K(a) approximately (0.35+/-0.1)x10(3) M(-1) for ZnMb. Brownian dynamic (BD) simulations show that cyt b(5) binds over a broad surface of Mb that includes its heme edge. The experimental results are described in terms of a dynamic docking model which proposes that Mb binds cyt b(5) in a large ensemble of protein binding conformations, not one or a few dominant ones, but that only a small subset are ET reactive. Aided by the BD simulations, this model explains why k(2) decreases with increasing pH: increasing pH not only weakens the binding affinity but also reduces the number of binding conformations with high ET reactivity.

Animals↗

PubMed1998

The structure of the two amino-terminal domains of human ICAM-1 suggests how it functions as a rhinovirus receptor and as an LFA-1 integrin ligand.

The normal function of human intercellular adhesion molecule-1 (ICAM-1) is to provide adhesion between endothelial cells and leukocytes after injury or stress. ICAM-1 binds to leukocyte function-associated antigen (LFA-1) or macrophage-1 antigen (Mac-1). However, ICAM-1 is also used as a receptor by the major group of human rhinoviruses and is a catalyst for the subsequent viral uncoating during cell entry. The three-dimensional atomic structure of the two amino-terminal domains (D1 and D2) of ICAM-1 has been determined to 2.2-A resolution and fitted into a cryoelectron microscopy reconstruction of a rhinovirus-ICAM-1 complex. Rhinovirus attachment is confined to the BC, CD, DE, and FG loops of the amino-terminal Ig-like domain (D1) at the end distal to the cellular membrane. The loops are considerably different in structure to those of human ICAM-2 or murine ICAM-1, which do not bind rhinoviruses. There are extensive charge interactions between ICAM-1 and human rhinoviruses, which are mostly conserved in both major and minor receptor groups of rhinoviruses. The interaction of ICAMs with LFA-1 is known to be mediated by a divalent cation bound to the insertion (I)-domain on the alpha chain of LFA-1 and the carboxyl group of a conserved glutamic acid residue on ICAMs. Domain D1 has been docked with the known structure of the I-domain. The resultant model is consistent with mutational data and provides a structural framework for the adhesion between these molecules.

Amino Acid Sequence↗

PubMed1999

Inhibition of dihydrofolate reductases from Toxoplasma gondii, Pneumocystis carinii, and rat liver by rotationally restricted analogues of pyrimethamine and metoprine.

Twenty-one conformationally restricted tricyclic pyrimethamine and metoprine analogues with one or two chlorine atoms, or other substituents, at different positions of the phenyl ring were tested for potency and species selectivity against dihydrofolate reductase (DHFR) from Toxoplasma gondii, Pneumocystis carinii, and rat liver. Heterocyclic systems studied included indeno[2,1-d]pyrimidines, benzo[f]quinazolines, and benzo[3,4]cyclohepta[1,2-d]pyrimidines. All but one of the analogues were more potent against T. gondii and rat liver DHFR than against P. carinii DHFR, and those with a one-carbon (CH2) bridge were generally less potent than those with a two-carbon (CH2CH2, CH=CH) or three-carbon (CH2CH2CH2) bridge. Although a number of compounds with a two- and three-carbon bridge were more potent than pyrimethamine against P. carinii DHFR, and especially T. gondii DHFR, none of them were selective for the P. carinii versus the mammalian enzyme, and only those with a one-carbon bridge showed selectivity approaching that of pyrimethamine for the T. gondii enzyme. Computer-simulated docking into the active site pocket of P. carinii and human DHFR suggested that, as a group, the rotationally restricted tricyclic structures are at a disadvantage relative to pyrimethamine and metoprine, in that torsional relief of unfavorable steric interactions between the chlorine atoms and two critical serine and threonine residues in the active site is prevented by the bridge.

Animals↗

PubMed2004

Brownian adhesive dynamics (BRAD) for simulating the receptor-mediated binding of viruses.

Current viral docking models have relied upon the assumption that bond formation and breakage are independent of viral and docking surface geometry, as well as the forces exerted on the bonds. This assumption, known as the equivalent site hypothesis (ESH), is examined in detail using a newly developed simulation technique-Brownian adhesive dynamics (BRAD). The simulation couples the thermal motion of viral particles with adhesive dynamics models to characterize the effect of bonding on viral motion. We use the binding of HIV-like particles to CD4 expressing cells as a model system to illustrate the utility of BRAD. Comparison of the transition rates between bound states predicted by ESH and the rates resulting from BRAD simulations show dramatic differences; at values of the equilibrium crosslinking constant, K(x)R(T), where ESH suggests all virus adhesion proteins will be bound (K(x)R(T) = 10(6)), BRAD predicts not all virus adhesion proteins will be bound. At values of the equilibrium crosslinking constant used in typical ESH calculations of virus docking (K(x)R(T) = 1) we find BRAD simulations predict no binding. The mean bond density from BRAD models is often much lower than that predicted by ESH for equivalent parameter values. BRAD suggests that the viruses are much less well bound than ESH predicts. The differences suggest that binding models for viruses need to be reexamined closely. BRAD is a simulation technique that will be useful for quantifying the receptor-mediated binding of a wide variety of viruses to cells.

Algorithms↗

PubMed1994

Molecular evolution and domain structure of plasminogen-related growth factors (HGF/SF and HGF1/MSP).

Plasminogen-related growth factors, a new family of polypeptide growth factors with the basic domain organization and mechanism of activation of the blood proteinase plasminogen, include hepatocyte growth factor/scatter factor (HGF/SF), a potent effector of the growth, movement, and differentiation of epithelia and endothelia, and hepatocyte growth factor-like/macrophage stimulating protein (HGF1/MSP), an effector of macrophage chemotaxis and phagocytosis. Phylogeny of the serine proteinase domains and analysis of intron-exon boundaries and kringle sequences indicate that HGF/SF, HGF1/MSP, plasminogen, and apolipoprotein (a) have evolved from a common ancestral gene that consisted of an N-terminal domain corresponding to plasminogen activation peptide (PAP), 3 copies of the kringle domain, and a serine proteinase domain. Models of the N domains of HGF/SF, HGF1/MSP, and plasminogen, characterized by the presence of 4 conserved Cys residues forming a loop in a loop, have been modeled based on disulfide-bond constraints. There is a distinct pattern of charged and hydrophobic residues in the helix-strand-helix motif proposed for the PAP domain of HGF/SF; these may be important for receptor interaction. Three-dimensional structures of the 4 kringle and the serine proteinase domains of HGF/SF were constructed by comparative modeling using the suite of programs COMPOSER and were energy minimized. Docking of a lysine analogue indicates a putative lysine-binding pocket within kringle 2 (and possibly another in kringle 4). The models suggest a mechanism for the formation of a noncovalent HGF/SF homodimer that may be responsible for the activation of the Met receptor. These data provide evidence for the divergent evolution and structural similarity of plasminogen, HGF/SF, and HGF1/MSP, and highlight a new strategy for growth factor evolution, namely the adaptation of a proteolytic enzyme to a role in receptor activation.

Amino Acid Sequence↗

PubMed1998

Light-activated rhodopsin induces structural binding motif in G protein alpha subunit.

A large superfamily of transmembrane receptors control cellular responses to diverse extracellular signals by catalyzing activation of specific types of heterotrimeric GTP-binding proteins. How these receptors recognize and promote nucleotide exchange on G protein alpha subunits to initiate signal amplification is unknown. The three-dimensional structure of the transducin (Gt) alpha subunit C-terminal undecapeptide Gtalpha(340-350) IKENLKDCGLF was determined by transferred nuclear Overhauser effect spectroscopy while it was bound to photoexcited rhodopsin. Light activation of rhodopsin causes a dramatic shift from a disordered conformation of Gtalpha(340-350) to a binding motif with a helical turn followed by an open reverse turn centered at Gly-348, a helix-terminating C capping motif of an alphaL type. Docking of the NMR structure to the GDP-bound x-ray structure of Gt reveals that photoexcited rhodopsin promotes the formation of a continuous helix over residues 325-346 terminated by the C-terminal helical cap with a unique cluster of crucial hydrophobic side chains. A molecular mechanism by which activated receptors can control G proteins through reversible conformational changes at the receptor-G protein interface is demonstrated.

Amino Acid Sequence↗

PubMed2005

Computer-assisted docking of flavodoxin with the ATP:Co(I)rrinoid adenosyltransferase (CobA) enzyme reveals residues critical for protein-protein interactions but not for catalysis.

The activity of the housekeeping ATP:co(I)rrinoid adenosyltransferase (CobA) enzyme of Salmonella enterica sv. Typhimurium is required to adenosylate de novo biosynthetic intermediates of adenosylcobalamin and to salvage incomplete and complete corrinoids from the environment of this bacterium. In vitro, reduced flavodoxin (FldA) provides an electron to generate the co(I)rrinoid substrate in the CobA active site. To understand how CobA and FldA interact, a computer model of a CobA.FldA complex was generated. This model was used to guide the introduction of mutations into CobA using site-directed mutagenesis and the synthesis of a peptide mimic of FldA. Residues Arg-9 and Arg-165 of CobA were critical for FldA-dependent adenosylation but were catalytically as competent as the wild-type protein when cob(I)alamin was provided as substrate. These results indicate that Arg-9 and Arg-165 are important for CobA.FldA docking but not to catalysis. A truncation of the 9-amino acid N-terminal helix of CobA reduced its FldA-dependent cobalamin adenosyltransferase activity by 97.4%. The same protein, however, had a 4-fold higher specific activity than the native enzyme when cob(I)alamin was generated chemically in situ.

Alkyl and Aryl Transferases↗

PubMed2004

Design and synthesis of novel dihydroxyindole-2-carboxylic acids as HIV-1 integrase inhibitors.

In a search for new HIV-1 integrase (IN) inhibitors, we synthesized and evaluated the biological activity of 5,6-dihydroxyindole-2-carboxylic acid (DHICA) and a series of its derivatives. These compounds were designed as conformationally constrained analogues of the acrylate moiety of caffeic acid phenethyl ester (CAPE). DHICA, an intermediate in the biosynthesis of melanins, was prepared as a monomeric unit by a novel synthetic route. In order to perform coherent SAR studies, two series of DHICA amides were synthesized. First, to validate the utility of a previously identified three-point pharmacophore based on CAPE in inhibitor design, we prepared a series of benzyl- or phenylethylamine substituted derivatives lacking and containing hydroxyl groups. Second, dimers of DHICA containing various aminoalkylamine linkers were also prepared with a goal to increase potency. All compounds were tested against purified IN and the C65S mutant in enzyme-based assays. They were also tested for cytotoxicity in an ovarian carcinoma cell line and antiviral activity in HIV-1-infected CEM cells. Seven compounds inhibited catalytic activities of purified IN with IC50 values below 10 microM. Further computational docking studies were performed to determine the title compounds' mode of interaction with the IN active site. The residues K156, K159 and D64 were the most important for potency against purified IN.

Catalysis↗

PubMed2006

Mutagenic potency of food-derived heterocyclic amines.

The understanding of mutagenic potency has been primarily approached using "quantitative structure-activity relationships" (QSAR). Often this method allows the prediction of mutagenic potency of the compound based on its structure. But it does not give the underlying reason why the mutagenic activities differ. We have taken a set of heterocyclic amine structures and used molecular dynamic calculations to dock these molecules into the active site of a computational model of the cytochrome P4501A2 enzyme. The calculated binding strength using Boltzman distribution constants was then compared to the QSAR value (HF/6-31G* optimized structures) and the Ames/Salmonella mutagenic potency. Further understanding will only come from knowing the complete set of mutagenic determinants. These include the nitrenium ion half-life, DNA adduct half-life, efficiency of repair of the adduct, and ultimately fixation of the mutation through cellular processes. For two isomers, PhIP and 3-Me-PhIP, we showed that for the 100-fold difference in the mutagenic potency a 5-fold difference can be accounted for by differences in the P450 oxidation. The other factor of 20 is not clearly understood but is downstream from the oxidation step. The application of QSAR (chemical characteristics) to biological principles related to mutagenesis is explored in this report.

Amines↗

PubMed2006

In silico prediction of drug binding to CYP2D6: identification of a new metabolite of metoclopramide.

Patients with cancer often take many different classes of drugs to treat the effects of their malignancy and the side effects of treatment, as well as their comorbidities. The potential for drug-drug interactions that may affect the efficacy of anticancer treatment is high, and a major source of such interactions is competition for the drug-metabolizing enzymes, cytochromes P450 (P450s). We have examined a series of 20 drugs commonly prescribed to cancer patients to look for potential interactions via CYP2D6. We used a homology model of CYP2D6, together with molecular docking techniques, to perform an in silico screen for binding to CYP2D6. Experimental IC50 values were determined for these compounds and compared with the model predictions to reveal a correlation with a regression coefficient of r2= 0.61. Importantly, the docked conformation of the commonly prescribed antiemetic metoclopramide predicted a new site of metabolism that was further investigated through in vitro analysis with recombinant CYP2D6. An aromatic N-hydroxy metabolite of metoclopramide, consistent with predictions from our modeling studies, was identified by high-performance liquid chromatography/mass spectrometry. This metabolite was found to represent a major product of metabolism in human liver microsomes, and CYP2D6 was identified as the main P450 isoform responsible for catalyzing its formation. In view of the prevalence of interindividual variation in the CYP2D6 genotype and phenotype, we suggest that those experiencing adverse reactions with metoclopramide, e.g., extrapyramidal syndrome, are likely to have a particular CYP2D6 genotype/phenotype. This warrants further investigation.

Antiemetics↗

PubMed1999

Modelling of factor Xa-inhibitor complexes: a computational flexible docking approach.

In order to understand the structural basis of Factor Xa (FXa) specificity, structural complexes of FXa with its synthetic inhibitors are determined using a computational docking approach. The AutoDock suite of programs is used to determine the binding modes of the synthetic inhibitors such as 3- and 4-amidinobenzylphenyl ether (ABP), amidinophenyl pyruvic acid (APPA), diamidinobenzofuranyl ethene (DABE), and DX-9065a 2-(5'-amidino-2'-benzofuranyl)-3-(7'amidino-2'-napthyl)-propionic acid (ABAP) to FXa. The synthetic inhibitors docked in the present study are different in size, nature of linkage, and properties. Two sets of simulations were carried out for synthetic inhibitors docking to FXa. In the first set of simulations, no explicit water molecules were included. In the second set of simulations two explicit solvent molecules were considered. In all the computationally predicted synthetic inhibitor complexes of FXa, the specificity pocket residue Asp-189 is involved in hydrogen bonding with the bound inhibitor. The active site water molecule WAT522 is involved in hydrogen bonding with all the bound inhibitors. The computed energies clearly discriminate the high affinity from low affinity binders.

Amidines↗

PubMed1992

Molecular dynamics characterization of the active cavity of carboxypeptidase A and some of its inhibitor adducts.

Molecular dynamics (MD) calculations have been performed on carboxypeptidase A and on its adducts with inhibitors, such as d-phenylalanine (dPhe) and acetate. The catalytically essential zinc ion present in the protein was explicitly included in all the simulations. The simulation was carried out over a sphere of 15 A centered on the zinc ion. The crystallographic water molecules were explicitly taken into account; then the protein was solvated with a 18 A sphere of water molecules. MD calculations were carried out for 45-60 ps. There is no large deviation from the available X-ray structures of native and the dPhe adduct for the MD structures. Average MD structures were calculated starting from the X-ray structure of the dPhe adduct, and, from a structure obtained by docking the inhibitor in the native structure. Comparison between these two structures and with that of the native protein shows that some of the key variations produced by inhibitor binding are reproduced by MD calculations. Addition of acetate induces structural changes relevant for the understanding of the interaction network in the active cavity. The structural variations induced by different inhibitors are examined. The effects of these interactions on the catalytic mechanism and on the binding of substrate are discussed.

Acetates↗