Pharmacological profile of opicapone, a third-generation nitrocatechol catechol-O-methyl transferase inhibitor, in the rat.

BACKGROUND AND PURPOSE: Catechol-O-methyltransferase (COMT) is an important target in the levodopa treatment of Parkinson's disease; however, the inhibitors available have problems, and not all patients benefit from their efficacy. Opicapone was developed to overcome those limitations. In this study, opicapone's pharmacological properties were evaluated as well as its potential cytotoxic effects. EXPERIMENTAL APPROACH: The pharmacodynamic effects of opicapone were explored by evaluating rat COMT activity and levodopa pharmacokinetics, in the periphery through microdialysis and in whole brain. The potential cytotoxicity risk of opicapone was explored in human hepatocytes by assessing cellular ATP content and mitochondrial membrane potential. KEY RESULTS: Opicapone inhibited rat peripheral COMT with ED50 values below 1.4 mg⋅kg(-1) up to 6 h post-administration. The effect was sustained over the first 8 h and by 24 h COMT had not returned to control values. A single administration of opicapone resulted in increased and sustained plasma levodopa levels with a concomitant reduction in 3-O-methyldopa from 2 h up to 24 h post-administration, while tolcapone produced significant effects only at 2 h post-administration. The effects of opicapone on brain catecholamines after levodopa administration were sustained up to 24 h post-administration. Opicapone was also the least potent compound in decreasing both the mitochondrial membrane potential and the ATP content in human primary hepatocytes after a 24 h incubation period. CONCLUSIONS AND IMPLICATIONS: Opicapone has a prolonged inhibitory effect on peripheral COMT, which extends the bioavailability of levodopa, without inducing toxicity. Thus, it exhibits some improved properties compared to the currently available COMT inhibitors.

Comparative study of ortho- and meta-nitrated inhibitors of catechol-O-methyltransferase: interactions with the active site and regioselectivity of O-methylation.

In this work, we present a comparative case study of "ortho-" and "meta-nitrated" catecholic inhibitors of catechol-O-methyltransferase (COMT), with regard to their interaction with the catalytic site of the enzyme and the in vitro regioselective formation of their mono-O-methyl ether metabolites. In particular, the effects of altering the attachment position of the inhibitors' side-chain substituent, within the classic nitrocatechol pharmacophore, were investigated. For this purpose, we compared two simple regioisomeric nitrocatechol-type inhibitors of COMT, BIA 3-228 and BIA 8-176, which contain the benzoyl substituent attached at the meta and ortho positions, respectively, relative to the nitro group. The two compounds were slowly O-methylated by COMT in vitro, but the particular substitution pattern of each compound was shown to have a profound impact on the regioselectivity of their O-methylation. To provide a plausible interpretation of these results, a comprehensive analysis of the protein-inhibitor interactions and of the relative chemical susceptibility to O-methylation of the catechol hydroxyl groups was performed by means of docking simulations and ab initio molecular orbital calculations. The major structural and chemical factors that determine the enzyme regioselectivity of O-methylation were identified, and the X-ray structure of the complex of COMT with S-adenosyl-l-methionine and BIA 8-176 is herein disclosed. This is the first reported structure of the soluble form of COMT complexed with a nitrocatecholic inhibitor having a bulky substituent group in adjacent position (ortho) to the nitro group. Structural and dynamic aspects of this complex are analyzed and discussed, in the context of the present study.

Molecular modeling and metabolic studies of the interaction of catechol-O-methyltransferase and a new nitrocatechol inhibitor.

Catechol-O-methyltransferase (COMT, EC 2.1.1.6) plays a central role in the metabolic inactivation of neurotransmitters and neuroactive xenobiotics possessing a catechol motif. 1-(3,4-Dihydroxy-5-nitrophenyl)-2-phenyl-ethanone (BIA 3-202) is a novel nitrocatechol-type inhibitor of COMT, the potential clinical benefit of which is currently being evaluated in the treatment of Parkinson's disease. In the present work we characterize the molecular interactions of BIA 3-202 within the active site of COMT and discuss their implication on the regioselectivity of metabolic O-methylation. Unrestrained flexible-docking simulations suggest that the solution structure of this complex is better described as an ensemble of alternative binding modes, in contrast to the well defined bound configuration revealed by the X-ray structures of related nitrocatechol inhibitors, co-crystallized with COMT. The docking results wherein presented are well supported by experimental evidence, where the pattern of in vitro enzymatic O-methylation and O-demethylation reactions are analyzed. We propose a plausible explanation for the paradoxical in vivo regioselectivity of O-methylation of BIA 3-202, as well as of its related COMT inhibitor tolcapone. Both compounds undergo in vivo O-methylation by COMT at either meta or para catechol hydroxyl groups. However, results herein presented suggest that, in a subsequent step, the p-O-methyl derivatives are selectively demethylated by a microsomal enzyme system. The overall balance is the accumulation of the m-O-methylated metabolites over the para-regioisomers. The implications for the general recognition of nitrocatechol-type inhibitors by COMT and the regioselectivity of their metabolic O-methylation are discussed.

A novel approach for assessing macromolecular complexes combining soft-docking calculations with NMR data.

We present a novel and efficient approach for assessing protein-protein complex formation, which combines ab initio docking calculations performed with the protein docking algorithm BiGGER and chemical shift perturbation data collected with heteronuclear single quantum coherence (HSQC) or TROSY nuclear magnetic resonance (NMR) spectroscopy. This method, termed "restrained soft-docking," is validated for several known protein complexes. These data demonstrate that restrained soft-docking extends the size limitations of NMR spectroscopy and provides an alternative method for investigating macromolecular protein complexes that requires less experimental time, effort, and resources. The potential utility of this novel NMR and simulated docking approach in current structural genomic initiatives is discussed.

BiGGER: a new (soft) docking algorithm for predicting protein interactions.

A new computationally efficient and automated "soft docking" algorithm is described to assist the prediction of the mode of binding between two proteins, using the three-dimensional structures of the unbound molecules. The method is implemented in a software package called BiGGER (Bimolecular Complex Generation with Global Evaluation and Ranking) and works in two sequential steps: first, the complete 6-dimensional binding spaces of both molecules is systematically searched. A population of candidate protein-protein docked geometries is thus generated and selected on the basis of the geometric complementarity and amino acid pairwise affinities between the two molecular surfaces. Most of the conformational changes observed during protein association are treated in an implicit way and test results are equally satisfactory, regardless of starting from the bound or the unbound forms of known structures of the interacting proteins. In contrast to other methods, the entire molecular surfaces are searched during the simulation, using absolutely no additional information regarding the binding sites. In a second step, an interaction scoring function is used to rank the putative docked structures. The function incorporates interaction terms that are thought to be relevant to the stabilization of protein complexes. These include: geometric complementarity of the surfaces, explicit electrostatic interactions, desolvation energy, and pairwise propensities of the amino acid side chains to contact across the molecular interface. The relative functional contribution of each of these interaction terms to the global scoring function has been empirically adjusted through a neural network optimizer using a learning set of 25 protein-protein complexes of known crystallographic structures. In 22 out of 25 protein-protein complexes tested, near-native docked geometries were found with C(alpha) RMS deviations < or =4.0 A from the experimental structures, of which 14 were found within the 20 top ranking solutions. The program works on widely available personal computers and takes 2 to 8 hours of CPU time to run any of the docking tests herein presented. Finally, the value and limitations of the method for the study of macromolecular interactions, not yet revealed by experimental techniques, are discussed.

Heteronuclear NMR and soft docking: an experimental approach for a structural model of the cytochrome c553-ferredoxin complex.

The combination of docking algorithms with NMR data has been developed extensively for the studies of protein-ligand interactions. However, to extend this development for the studies of protein-protein interactions, the intermolecular NOE constraints, which are needed, are more difficult to access. In the present work, we describe a new approach that combines an ab initio docking calculation and the mapping of an interaction site using chemical shift variation analysis. The cytochrome c553-ferredoxin complex is used as a model of numerous electron-transfer complexes. The 15N-labeling of both molecules has been obtained, and the mapping of the interacting site on each partner, respectively, has been done using HSQC experiments. 1H and 15N chemical shift analysis defines the area of both molecules involved in the recognition interface. Models of the complex were generated by an ab initio docking software, the BiGGER program (bimolecular complex generation with global evaluation and ranking). This program generates a population of protein-protein docked geometries ranked by a scoring function, combining relevant stabilization parameters such as geometric complementarity surfaces, electrostatic interactions, desolvation energy, and pairwise affinities of amino acid side chains. We have implemented a new module that includes experimental input (here, NMR mapping of the interacting site) as a filter to select the accurate models. Final structures were energy minimized using the X-PLOR software and then analyzed. The best solution has an interface area (1037.4 A2) falling close to the range of generally observed recognition interfaces, with a distance of 10.0 A between the redox centers.

The surface-charge asymmetry and dimerisation of cytochrome c550 from Paracoccus denitrificans--implications for the interaction with cytochrome c peroxidase.

The implications of the dimeric state of cytochrome c550 for its binding to Paracoccus cytochrome c peroxidase and its delivery of the two electrons required to restore the active enzyme during catalysis have been investigated. The amino acid sequence of cytochrome c550 of Paracoccus denitrificans strain LMD 52.44 was determined and showed 21 differences from that of strain LMD 22.21. Based on the X-ray structure of the latter, a structure for the cytochrome c550 monomer from strain 52.44 is proposed and a dipole moment of 945 debye was calculated with an orientation close to the exposed haem edge. The behaviour of the cytochrome on molecular-exclusion chromatography is indicative of an ionic strength-dependent monomer (15 kDa)/dimer (30 kDa) equilibrium that can also be detected by 1H-NMR spectroscopy. The apparent mass of 50 kDa observed at very low ionic strength was consistent with the presence of a strongly asymmetric dimer. This was confirmed by cross-linking studies, which showed that a cross-linked species of mass 30 kDa on SDS behaved with an apparent mass of 50 kDa on molecular-exclusion chromatography. A programme which carried out and evaluated molecular docking of two monomers to give a dimer generated a most probable dimer in which the monomer dipoles lay almost antiparallel to each other. The resultant dipole moment of the dimer is therefore small. Although this finding calls into question the possibility of preorientation of a strongly asymmetrically charged cytochrome as it collides with a redox partner, the stoichiometry of complex formation with cytochrome c peroxidase as studied by 1H-NMR spectroscopy shows that it is the monomer that binds.

Structure of the tetraheme cytochrome from Desulfovibrio desulfuricans ATCC 27774: X-ray diffraction and electron paramagnetic resonance studies.

The three-dimensional X-ray structure of cytochrome c3 from a sulfate reducing bacterium, Desulfovibrio desulfuricans ATCC 27774 (107 residues, 4 heme groups), has been determined by the method of molecular replacement [Frazão et al. (1994) Acta Crystallogr. D50, 233-236] and refined at 1.75 A to an R-factor of 17.8%. When compared with the homologous proteins isolated from Desulfovibrio gigas, Desulfovibrio vulgaris Hildenborough, Desulfovibrio vulgaris Miyazaki F, and Desulfomicrobium baculatus, the general outlines of the structure are essentialy kept [heme-heme distances, heme-heme angles, His-His (axial heme ligands) dihedral angles, and the geometry of the conserved aromatic residues]. The three-dimensional structure of D. desulfuricans ATCC 27774 cytochrome c3Dd was modeled on the basis of the crystal structures available and amino acid sequence comparisons within this homologous family of multiheme cytochromes [Palma et al. (1994) Biochemistry 33, 6394-6407]. This model is compared with the refined crystal structure now reported, in order to discuss the validity of structure prediction methods and critically evaluate the steps used to predict protein structures by homology modeling. The four heme midpoint redox potentials were determined by using deconvoluted electron paramagnetic resonance (EPR) redox titrations. Structural criteria (electrostatic potentials, heme ligand orientation, EPR g values, heme exposure, data from protein-protein interaction studies) are invoked to assign the redox potentials corresponding to each specific heme in the three-dimensional structure.

Evidence for a ternary complex formed between flavodoxin and cytochrome c3: 1H-NMR and molecular modeling studies.

Small electron-transfer proteins such as flavodoxin (16 kDa) and the tetraheme cytochrome c3 (13 kDa) have been used to mimic, in vitro, part of the complex electron-transfer chain operating between substrate electron donors and respiratory electron acceptors, in sulfate-reducing bacteria (Desulfovibrio species). The nature and properties of the complex formed between these proteins are revealed by 1H-NMR and molecular modeling approaches. Our previous study with the Desulfovibrio vulgaris proteins [Moura, I., Moura, J.J. G., Santos, M.H., & Xavier, A. V. (1980) Cienc. Biol. (Portugal) 5, 195-197; Stewart, D.E. LeGall, J., Moura, I., Moura, J. J. G., Peck, H.D. Jr., Xavier, A. V., Weiner, P. K., & Wampler, J.E. (1988) Biochemistry 27, 2444-2450] indicated that the complex between cytochrome c3 and flavodoxin could be monitored by changes in the NMR signals of the heme methyl groups of the cytochrome and that the electrostatic surface charge (Coulomb's law) on the two proteins favored interaction between one unique heme of the cytochrome with flavodoxin. If the interaction is indeed driven by the electrostatic complementarity between the acidic flavodoxin and a unique positive region of the cytochrome c3, other homologous proteins from these two families of proteins might be expected to interact similarly. In this study, three homologous Desulfovibrio cytochromes c3 were used, which show a remarkable variation in their individual isoelectric points (ranging from 5.5 to 9.5). On the basis of data obtained from protein-protein titrations followed at specific proton NMR signals (i.e., heme methyl resonances), a binding model for this complex has been developed with evaluation of stoichiometry and binding constants. This binding model involves one site on the cytochromes c3 and two sites on the flavodoxin, with formation of a ternary complex at saturation. In order to understand the potential chemical form of the binding model, a structural model for the hypothetical ternary complex, formed between one molecule of Desulfovibrio salexigens flavodoxin and two molecules of cytochrome c3, is proposed. These molecular models of the complexes were constructed on the basis of complementarity of Coulombic electrostatic surface potentials, using the available X-ray structures of the isolated proteins and, when required, model structures (D. salexigens flavodoxin and Desulfovibrio desulfuricans ATCC 27774 cytochrome c3) predicted by homology modeling.

Primary sequence, oxidation-reduction potentials and tertiary-structure prediction of Desulfovibrio desulfuricans ATCC 27774 flavodoxin.

Flavodoxin was isolated and purified from Desulfovibrio desulfuricans ATCC 27774, a sulfate-reducing organism that can also utilize nitrate as an alternative electron acceptor. Mid-point oxidation-reduction potentials of this flavodoxin were determined by ultraviolet/visible and EPR methods coupled to potentiometric measurements and their pH dependence studied in detail. The redox potential E2, for the couple oxidized/semiquinone forms at pH 6.7 and 25 degrees C is -40 mV, while the value for the semiquinone/hydroquinone forms (E1), at the same pH, -387 mV. E2 varies linearly with pH, while E1 is independent of pH at high values. However, at low pH (< 7.0), this value is less negative, compatible with a redox-linked protonation of the flavodoxin hydroquinone. A comparative study is presented for Desulfovibrio salexigens NCIB 8403 flavodoxin [Moura, I., Moura, J.J.G., Bruschi, M. & LeGall, J. (1980) Biochim. Biophys. Acta 591, 1-8]. The complete primary amino acid sequence was obtained by automated Edman degradation from peptides obtained by chemical and enzymic procedures. The amino acid sequence was confirmed by FAB/MS. Using the previously determined tridimensional structure of Desulfovibrio vulgaris flavodoxin as a model [similarity, 48.6%; Watenpaugh, K.D., Sieker, L.C., Jensen, L.H., LeGall, J. & Dubourdieu M. (1972) Proc. Natl Acad. Sci. USA 69, 3185-3188], the tridimensional structure of D. desulfuricans ATCC 27774 flavodoxin was predicted using AMBER force-field calculations.