Construction of a 3D model of oligopeptidase B, a potential processing enzyme in prokaryotes.

A three dimensional structural model of oligopeptidase B (OpB) was constructed by homology modeling. High resolution X-ray structure of prolyl oligopeptidase (PEP), the only protein with sequential and functional homology was used as a template. Initial models of OpB were built by the MODELLER and were analysed by the PROCHECK programs. The best quality model was chosen for further refinement by two different techniques--either constrained molecular dynamics simulations or simulated annealing calculations starting from 500 K. The overall quality of each of the refined models was evaluated and the simulated annealing procedure found to be more effective. The refined model was analysed by different protein analysis programs including PROCHECK for the evaluation of the Ramachandran plot quality, PROSA for testing interaction energies and WHATIF for the calculation of packing quality. This structure was found to be satisfactory and also stable at room temperature as demonstrated by a 300 ps long unconstrained molecular dynamics simulation. Calculation of molecular electrostatic potentials revealed that the binding site of OpB is more negative than that of PEP, in accordance with the experimentally observed selectivity of OpB towards proteolysis at dibasic sites. A recently developed Monte Carlo docking method was used provide a structural rationale for the affinity differences measured between Z-Arg and Z-Arg-Arg substrates.

Chemical fragmentation in quantum mechanical methods.

We give a survey on the application of the chemical fragmentation concept in computer modelling of extended covalent systems. It will be stressed that information on molecular topology, as well as location and composition of the reaction centre allows the construction of a reasonable initial guess for the wave function and thus facilitates the solution of the Schrödinger equation. For systems, where the chemical changes are localised to a few atoms, while others play the role of essentially electrostatic perturbation, a partition into active site and environment is possible providing a background to hybrid quantum mechanical/molecular mechanical (QM/MM) methods. Full molecular orbital treatment of large covalent systems at the minimal basis, semiempirical level becomes possible in the frame of the fragment self-consistent field (FSCF) method which was developed in the past two decades in our laboratory. As an application, we discuss the hydride shift reaction step in xylose isomerase catalysis.

A new potent calmodulin antagonist with arylalkylamine structure: crystallographic, spectroscopic and functional studies.

An arylalkylamine-type calmodulin antagonist, N-(3, 3-diphenylpropyl)-N'-[1-R-(3, 4-bis-butoxyphenyl)ethyl]-propylene-diamine (AAA) is presented and its complexes with calmodulin are characterized in solution and in the crystal. Near-UV circular dichroism spectra show that AAA binds to calmodulin with 2:1 stoichiometry in a Ca(2+)-dependent manner. The crystal structure with 2:1 stoichiometry is determined to 2.64 A resolution. The binding of AAA causes domain closure of calmodulin similar to that obtained with trifluoperazine. Solution and crystal data indicate that each of the two AAA molecules anchors in the hydrophobic pockets of calmodulin, overlapping with two trifluoperazine sites, i.e. at a hydrophobic pocket and an interdomain site. The two AAA molecules also interact with each other by hydrophobic forces. A competition enzymatic assay has revealed that AAA inhibits calmodulin-activated phosphodiesterase activity at two orders of magnitude lower concentration than trifluoperazine. The apparent dissociation constant of AAA to calmodulin is 18 nM, which is commensurable with that of target peptides. On the basis of the crystal structure, we propose that the high-affinity binding is mainly due to a favorable entropy term, as the AAA molecule makes multiple contacts in its complex with calmodulin.

Mechanism of action of D-xylose isomerase.

The present knowledge on the stereochemical mechanism of action of glucose (or xylose) isomerase, one of the highest tonnage industrial enzymes, is summarized. First we deal shortly with experimental methods applied to study the structure and function of this enzyme: enzyme kinetics, protein engineering, X-ray crystallography, nuclear magnetic and electron paramagnetic resonance spectroscopy. Computational methods like homology modeling, molecular orbital, molecular dynamics and continuum electrostatic methods are also shortly treated. We discuss mostly those results and their contribution to the elucidation of the mechanism of action that have been published in the last decade. Structural characteristics of free xylose isomerase as well as its complexes with various ligands are depicted. This information provides a tool for the study of structural details of the enzyme mechanism. We present a general mechanism where the first step is ring opening, which is followed by the extension of the substrate to an open-chain conformation, a proton shuttle with the participation of a structural water molecule and the rate-determining hydride shift. The role of metal ions in the catalytic process is discussed in detail. Finally we present main trends in efforts of engineering the enzyme and delineate the prospective future lines. The review is completed by an extended bibliography with over 100 citations.

Comparative redox and pKa calculations on cytochrome c3 from several Desulfovibrio species using continuum electrostatic methods.

A comparative study of the pH-dependent redox mechanisms of several members of the cytochrome c3 family has been carried out. In a previous work, the molecular determinants of this dependency (the so-called redox-Bohr effect) were investigated for one species using continuum electrostatic methods to find groups with a titrating range and strength of interaction compatible with a mediating role in the redox-Bohr effect. Here we clarify these aspects in the light of new and improved pKa calculations, our findings supporting the hypothesis of propionate D from heme I being the main effector in the pH-dependent modulation of the cytochrome c3 redox potentials in all the c3 molecules studied here. However, the weaker (but significant) role of other titrating groups cannot be excluded, their importance and identity changing with the particular molecule under study. We also calculate the relative redox potentials of the four heme centers among the selected members of the c3 family, using a continuum electrostatic method that takes into account both solvation and interaction effects. Comparison of the calculated values with available data for the microscopic redox potentials was undertaken, the quality of the agreement being dependent upon the choice of the dielectric constant for the protein interior. We find that high dielectric constants give best correlations, while low values result in better magnitudes for the calculated potentials. The possibility that the crystallographic calcium ion in c3 from Desulfovibrio gigas may be present in the solution structure was tested, and found to be likely.

Structure determination and refinement of the Al3+ complex of the D254,256E mutant of Arthrobacter D-xylose isomerase at 2.40 A resolution. Further evidence for inhibitor-induced metal ion movement.

The structure of the D254.256E double mutant of Arthrobacter xylose isomerase with Al3+ at both metal-binding sites was determined by the molecular replacement method at a conventional R-factor of 0.179. Binding of the two Al3+ does not alter the overall structure significantly. However, there are local rearrangements in the octahedral co-ordination sphere of the Al3+. The inhibitor molecule moves somewhat away from the active site. Furthermore, evidence was revealed for metal ion movement from site 2(1) to site 2(2) upon double mutation. Xylose isomerase requires two divalent metal cations for activation. The catalytic metal ion is translocated 1.8 A away from its initial position during the catalytic reaction. The fact that both activating and inactivating metals (including Al3+) were found exclusively at a single location in the double mutant was an indication that the consequently missing shuttle may account for the crippled catalytic efficiency.

The three-dimensional structure of Asp189Ser trypsin provides evidence for an inherent structural plasticity of the protease.

Trypsin mutant Asp189Ser, first described by Gráf et al. [Gráf, L., Jancsó, A., Szilágyi, L., Hegyi, G., Pintér, K., Náray-Szabó, G., Hepp, J., Medzihradszky, K. & Rutter, W.J. (1988) Proc. Natl Acad. Sci. USA 85, 4961-4965] has played an important role in recent studies on the structural basis of substrate-specific catalysis by serine proteases. The present work reports the three-dimensional structure of this mutant crystallized in unliganded form: the first unliganded rat trypsin structure reported. The X-ray structure of the Asp189Ser trypsin mutant in complex with bovine pancreatic trypsin inhibitor is already known. The X-ray structure of free Asp189Ser rat trypsin revealed that the single amino acid mutation at the bottom of the substrate binding pocket of trypsin resulted in extensive structural changes around the mutated site and in dimerization of the mutant, in contrast with the complexed enzyme the structure of which is practically the same as that of wild-type trypsin. The structural rearrangement in the mutant was shown to be restricted to the activation domain region providing further evidence for the allosteric property of this structural-functional unit of the enzyme. This study supports our view that the plasticity of the activation domain may play an important role in the mechanism of substrate-specific serine protease action.

Application of neural networks in structure-activity relationships.

Methodology and application of artificial neural networks in structure-activity relationships are reviewed focusing on the most frequently used three-layer feedforward back-propagation procedure. Two applications of neural networks are presented and a comparison of the performance with those of CoMFA and a classical QSAR analysis is also discussed.

Simultaneous binding of drugs with different chemical structures to Ca2+-calmodulin: crystallographic and spectroscopic studies.

The modulatory action of Ca2+-calmodulin on multiple targets is inhibited by trifluoperazine, which competes with target proteins for calmodulin binding. The structure of calmodulin crystallized with two trifluoperazine molecules is determined by X-ray crystallography at 2.74 A resolution. The X-ray data together with the characteristic and distinct signals obtained by circular dichroism in solution allowed us to identify the binding domains as well as the order of the binding of two trifluoperazine molecules to calmodulin. Accordingly, the binding of trifluperazine to the C-terminal hydrophobic pocket is followed by the interaction of the second drug molecule with an interdomain site. Recently, we demonstrated that the two bisindole derivatives, vinblastine and KAR-2 [3"-(beta-chloroethyl)-2",4"-dioxo-3, 5"-spirooxazolidino-4-deacetoxyvinblastine], interact with calmodulin with comparable affinity; however, they display different functional effects [Orosz et al. (1997) British J. Pharmacol. 121, 955-962]. The structural basis responsible for these effects were investigated by circular dichroism and fluorescence spectroscopy. The data provide evidence that calmodulin can simultaneously accommodate trifluoperazine and KAR-2 as well as vinblastine and KAR-2, but not trifluoperazine and vinblastine. The combination of the binding and structural data suggests that distinct binding sites exist on calmodulin for vinblastine and KAR-2 which correspond, at least partly, to that of trifluoperazine at the C-terminal hydrophobic pocket and at an interdomain site, respectively. This structural arrangement can explain why these drugs display different anticalmodulin activities. Calmodulin complexed with melittin is also able to bind two trifluoperazine molecules, the binding of which appears to be cooperative. Results obtained with intact and proteolytically cleaved calmodulin reveal that the central linker region of the protein is indispensable for simultanous interactions with two molecules of either identical or different ligands.

[Methods of computer-assisted molecular modeling]. / Számítógépes molekulamodellezési módszerek.

We present an overview on modern computer methods of molecular modelling. After treating three main steps of drug evaluation, namely target identification, lead identification and lead optimisation, we shortly discuss molecular graphics, molecular mechanics, molecular orbital and molecular dynamics methods. These are suitable for the more-or-less adequate modelling of real molecular processes both at the microscopic and the macroscopic levels. Molecular graphics provides beautiful pictures for the specialist that allow inspection and manipulation of detailed molecular models. An especially useful tool for the visualisation of molecular entities is the display of various properties on the molecular surface that allows rapid recognition of important relationships. Molecular mechanics is able to predict properties (e.g. geometric parameters, conformer stability) of several classes of molecules with an accuracy close to the experimental one, therefore it plays an important role in complementing molecular graphics. The performance of molecular orbital methods increased considerably in the last decade thus we can calculate parameters for isolated or interacting molecules that are not easily amenable to experiment (e.g. structure and energetics of unstable species or activation energies of elementary processes). Computer simulation methods provide a link between gas-phase models of microscopic structures or processes and macroscopic properties or events that may be derived from the former. Thus, it became possible to apply computerised methods for an adequate simulation of important events, like chemical and biochemical reactions, drug-target interactions, drug delivery and the similar that determine drug action. It is stressed that the hardware and software for computer-aided molecular modelling may not be absent from the arsenal of a drug designer.

Crystallization and preliminary X-ray analysis of porcine muscle prolyl oligopeptidase.

Prolyl oligopeptidase from pig muscle has been crystallized in complex with an inhibitor, using PEG 8000 and calcium acetate as precipitants. The crystals are orthorombic and the space group is P212121 with cell dimensions a = 111.8, b = 101.8, c = 72.4 A. The asymmetric unit contains a single chain of prolyl oligopeptidase, corresponding to a specific volume of 2.55 A3 Da-1 and a solvent content of 52%. The observed diffraction pattern extends to 2.3 A resolution and the native crystals are well suited for structural analysis by X-ray diffraction methods.

Role of electrostatics at the catalytic metal binding site in xylose isomerase action: Ca(2+)-inhibition and metal competence in the double mutant D254E/D256E.

The catalytic metal binding site of xylose isomerase from Arthrobacter B3728 was modified by protein engineering to diminish the inhibitory effect of Ca2+ and to study the competence of metals on catalysis. To exclude Ca2+ from Site 2 a double mutant D254E/D256E was designed with reduced space available for binding. In order to elucidate structural consequences of the mutation the binary complex of the mutant with Mg2+ as well as ternary complexes with bivalent metal ions and the open-chain inhibitor xylitol were crystallized for x-ray studies. We determined the crystal structures of the ternary complexes containing Mg2+, Mn2+, and Ca2+ at 2.2 to 2.5 A resolutions, and refined them to R factors of 16.3, 16.6, and 19.1, respectively. We found that all metals are liganded by both engineered glutamates as well as by atoms O1 and O2 of the inhibitor. The similarity of the coordination of Ca2+ to that of the cofactors as well as results with Be2+ weaken the assumption that geometry differences should account for the catalytic noncompetence of this ion. Kinetic results of the D254E/D256E mutant enzyme showed that the significant decrease in Ca2+ inhibition was accompanied by a similar reduction in the enzymatic activity. Qualitative argumentation, based on the protein electrostatic potential, indicates that the proximity of the negative side chains to the substrate significantly reduces the electrostatic stabilization of the transition state. Furthermore, due to the smaller size of the catalytic metal site, no water molecule, coordinating the metal, could be observed in ternary complexes of the double mutant. Consequently, the proton shuttle step in the overall mechanism should differ from that in the wild type. These effects can account for the observed decrease in catalytic efficiency of the D254E/D256E mutant enzyme.

Crystallization and preliminary diffraction analysis of Ca(2+)-calmodulin-drug and apocalmodulin-drug complexes.

Ca(2+)-calmodulin is crystallized with two new and potent drugs: a bisindol derivative (KAR-2, 3"-(beta-chloroethyl)-2",4"-dioxo-3,5"- spiro-oxazolidino-4-deacetoxy-vinblastine) with antitumor activity and an arylalkylamine fendiline analogue (N-(3,3-diphenylpropyl)-N'-[1-(3,4- di-n-butoxy-phenyl)-ethyl]-1,3-diaminopropane) with anticalmodulin activity. The crystals diffract beyond 2.8 A and differ in unit cell parameters from each other as well as from crystals of Ca(2+)-calmodulin or Ca(2+)-calmodulin-ligand complexes, as reported thus far. Attempts to crystallize Ca(2+)-free calmodulin without drugs failed, in consonance with earlier results; however, single Ca(2+)-free calmodulin crystals diffracting-beyond 2.5 A resolution were grown in the presence of KAR-2. Results indicate that binding of the two drugs to apocalmodulin or Ca(2+)-calmodulin may induce unique novel protein conformers, targets of further detailed X-ray studies.

Molecular modelling of xylose isomerase catalysis: the role of electrostatics and charge transfer to metals.

The two main steps of the mechanism of xylose-xylulose conversion catalysed by D-xylose isomerase, the ring opening of xylose and the isomerization of the opened product by hydride transfer, were investigated by molecular mechanical and molecular orbital techniques. The activation energies calculated for these reactions clearly showed that hydrogen transfer is the rate-determining step of the enzymatic isomerization and that Mg2+ ions activate whereas Zn2+ ions inhibit the reaction, in agreement with the experiments. The remarkable differences between the net charges of these ions found by molecular orbital calculations and the inspection of the protein electrostatic potential around the reaction intermediates indicate that the main role of bivalent metal ions should be the electrostatic stabilization of the substrate transition states. In order to propose a more detailed mechanism, an attempt was made to clarify the effects of nearby residues (e.g. His54, Asp57, Lys183, Asp257) in the reaction. Different isomerization mechanisms, such as through an enediol intermediate, were examined and could be excluded, in addition to the charge-relay mechanism during the ring opening.

Analysis of molecular recognition: steric electrostatic and hydrophobic complementarity.

We discuss three important aspects of molecular recognition: steric, electrostatic and hydrophobic. Steric fit means that interacting atoms may not approach each other beyond their van der Waals radii and, simultaneously, crevices should be filled as densely as possible. Electrostatic fit requires the maximum ionic and polar (hydrogen bond or other) interaction between host and guest atoms while the hydrophobic fit corresponds to the association trend between apolar groups in an aqueous medium. Space-filling models, obtained by molecular graphics, illustrate steric complementarity while we use molecular electrostatic potentials (MEPs) and fields (MEFs) to investigate electrostatic and hydrophobic matching. Molecular regions with negative and positive MEPs attract and repel a positive probe charge, respectively, so we consider them as attracting each other. Furthermore we postulate that regions with MEFs of similar magnitude tend to associate more strongly than those with very different fields (similis simili gaudet principle). We apply the above rules to the study of complementarity in the trypsin-BPTI complex and in a crystalline association between styrene epoxide as guest and a camphor-based anthracene derivative as host. We discuss molecular similarity on the same footing as complementarity and give some examples on the application of the concept to the rationalization of relative strengths of trypsin inhibitors.

Electrostatic complementarity in molecular associations.

In this paper, I attempt to summarize the main qualitative features of electrostatic complementarity and similarity, important determinants of molecular recognition. The two aspects, Coulombic and hydrophobic matching, can be formulated in terms of molecular electrostatic potentials and fields. The Coulombic aspect is equivalent to the requirement to produce a potential pattern in the host cavity that is opposite in sign to that emerging from a guest. Hydrophobic complementarity is best described by the similis simili gaudet principle. This means that field patterns near the interacting molecular surfaces must be of similar magnitude. The above rules, which may find useful application in molecular graphics, were studied for different cases of enzyme-ligand interactions in trypsin. A further example, a noncovalent structural model of the catalytic diad in Streptomyces Griseus protease A, supports the observation that the same molecular entities form similar associations even in different environments, as is the case in the complex of small species in a crystal and amino acid residues with structural water molecules in a protein.

How do serine proteases really work?

Recent advances in genetic engineering have led to a growing acceptance of the fact that enzymes work like other catalysts by reducing the activation barriers of the corresponding reactions. However, the key question about the action of enzymes is not related to the fact that they stabilize transition states but to the question to how they accomplish this task. This work considers the catalytic reaction of serine proteases and demonstrates how one can use a combination of calculations and experimental information to elucidate the key contributions to the catalytic free energy. Recent reports about genetic modifications of the buried aspartic group in serine proteases, which established the large effect of this group (but could not determine its origin), are analyzed. Two independent methods indicate that the buried aspartic group in serine proteases stabilizes the transition state by electrostatic interactions rather than by alternative mechanisms. Simple free energy considerations are used to eliminate the double proton-transfer mechanism (which is depicted in many textbooks as the key catalytic factor in serine proteases). The electrostatic stabilization of the oxyanion side of the transition state is also considered. It is argued that serine proteases and other enzymes work by providing electrostatic complementarity to the changes in charge distribution occurring during the reactions they catalyze.

Electrostatic complementarity within the substrate-binding pocket of trypsin.

The aspartic residue (Asp-189) at the base of the substrate-binding pocket of trypsin was replaced by serine (present in a similar position in chymotrypsin) through site-directed mutagenesis. The wild-type (with Asp-189 in the mature trypsin sequence) and mutant (Ser-189) trypsinogens were expressed in Escherichia coli, purified to homogeneity, activated by enterokinase, and tested with a series of fluorogenic tetrapeptide substrates with the general formula succinyl-Ala-Ala-Pro-Xaa-AMC, where AMC is 7-amino-4-methyl-coumarin and Xaa is Lys, Arg, Tyr, Phe, Leu, or Trp. As compared to [Asp189]trypsin, the activity of [Ser189]trypsin on lysyl and arginyl substrates decreased by about 5 orders of magnitude while its Km values increased only 2- to 6-fold. In contrast, [Ser189]trypsin was 10-50 times more active on the less preferred, chymotrypsin-type substrates (tyrosyl, phenylalanyl, leucyl, and tryptophanyl). The activity of [Ser189]trypsin on lysyl substrate was about 100-fold greater at pH 10.5 than at pH 7.0, indicating that the unprotonated lysine is preferred. Assuming the reaction mechanisms of the wild-type and mutant enzymes to be the same, we calculated the changes in the transition-state energies for various enzyme-substrate pairs to reflect electrostatic and hydrogen-bond interactions. The relative binding energies (E) in the transition state are as follows: EII greater than EPP greater than EPA greater than EIP approximately equal to EIA, where I = ionic, P = nonionic but polar, and A = apolar residues in the binding pocket. These side-chain interactions become prominent during the transition of the Michaelis complex to the tetrahedral transition-state complex.

Role of counter ions in trypsin acylation. NaCl effect.

The Asp102-carboxylate negative charge of the trypsin catalytic-triad has been substituted in part by Cl- counter-ions. A His57-imidazolium cation and Ser195-tetrahedral oxyanion ionpair generated in the acylation step of catalysis is stabilized by the negative charge of Asp102. The importance of correct location of this negative charge has been investigated by experimental analysis of the NaCl influence on the acylation-rate of trypsin, as well as by electrostatic, potential calculations. The acylation-rate was determined with stopped-flow technique under pseudo-first order conditions, by using 4-nitrophenyl-4'-guanidinium benzoate active site titrant at pH 4.25 +/- 0.04, in an unbuffered solution of I = O or 0.5 M NaCl. The acylation-rate constants, kappa 2, are: kappa H2O = 0.32 +/- 0.02 s-1 and kappa NaCl = 3.5 +/- 0.5 s-1, and they correlate to beta-trypsin (the most rapid single-chained form of the enzyme). The rate increasing effect of NaCl, together with the calculated electrostatic energies, indicate that the negative charge contribution of the Asp102-carboxylate to the stabilization of the imidazolium cation and tetrahedral oxyanion intermediate is larger of more orders, than that of the Cl- counter-ion, located in a less favourable position.