Publications

Export 41 results:
Sort by: Author Title Type [ Year  (Desc)]
2020
Designed Metal-ATCUN Derivatives: Redox and Non-redox-Based Applications Relevant for Chemistry, Biology, and Medicine, Maiti, B. K., Govil N., Kundu T., and J.J.G. Moura , iScience, Volume 23, p.101792, (2020)
2017
Molybdenum and tungsten-containing enzymes: an overview, Maia, L. B., Moura I., and Moura J. J. G. , Molybdenum and Tungsten Enzymes: Biochemistry, RSC Metallobiology Series No. 5 (ISBN: 978-1-78262-089-1). , p.1-80, (2017) mo_w_enzymes-rsc_book_biochemistry-chap_1.pdf
2016
Predicting Protein-Protein Interactions Using BiGGER: Case Studies, Almeida, R. M., Dell'Acqua S., Krippahl L., Moura J. J. G., and Pauleta S. R. , Molecules, Volume 21, p.1037, (2016) Website
2011
Crystal structure of the zinc-, cobalt-, and iron-containing adenylate kinase from Desulfovibrio gigas: a novel metal-containing adenylate kinase from Gram-negative bacteria, Mukhopadhyay, A., Kladova A. V., Bursakov S. A., Gavel O. Y., Calvete J. J., Shnyrov V. L., Moura I., Moura J. J., Romao M. J., and Trincao J. , J Biol Inorg Chem, Jan, Volume 16, Number 1, p.51-61, (2011) AbstractWebsite

Adenylate kinases (AK) from Gram-negative bacteria are generally devoid of metal ions in their LID domain. However, three metal ions, zinc, cobalt, and iron, have been found in AK from Gram-negative bacteria. Crystal structures of substrate-free AK from Desulfovibrio gigas with three different metal ions (Zn(2+), Zn-AK; Co(2+), Co-AK; and Fe(2+), Fe-AK) bound in its LID domain have been determined by X-ray crystallography to resolutions 1.8, 2.0, and 3.0 A, respectively. The zinc and iron forms of the enzyme were crystallized in space group I222, whereas the cobalt-form crystals were C2. The presence of the metals was confirmed by calculation of anomalous difference maps and by X-ray fluorescence scans. The work presented here is the first report of a structure of a metal-containing AK from a Gram-negative bacterium. The native enzyme was crystallized, and only zinc was detected in the LID domain. Co-AK and Fe-AK were obtained by overexpressing the protein in Escherichia coli. Zn-AK and Fe-AK crystallized as monomers in the asymmetric unit, whereas Co-AK crystallized as a dimer. Nevertheless, all three crystal structures are very similar to each other, with the same LID domain topology, the only change being the presence of the different metal atoms. In the absence of any substrate, the LID domain of all holoforms of AK was present in a fully open conformational state. Normal mode analysis was performed to predict fluctuations of the LID domain along the catalytic pathway.

2009
Cobalt-, zinc- and iron-bound forms of adenylate kinase (AK) from the sulfate-reducing bacterium Desulfovibrio gigas: purification, crystallization and preliminary X-ray diffraction analysis, Kladova, A. V., Gavel O. Y., Mukhopaadhyay A., Boer D. R., Teixeira S., Shnyrov V. L., Moura I., Moura J. J., Romao M. J., Trincao J., and Bursakov S. A. , Acta Crystallogr Sect F Struct Biol Cryst Commun, Sep 1, Volume 65, Number Pt 9, p.926-9, (2009) AbstractWebsite

Adenylate kinase (AK; ATP:AMP phosphotransferase; EC 2.7.4.3) is involved in the reversible transfer of the terminal phosphate group from ATP to AMP. AKs contribute to the maintenance of a constant level of cellular adenine nucleotides, which is necessary for the energetic metabolism of the cell. Three metal ions, cobalt, zinc and iron(II), have been reported to be present in AKs from some Gram-negative bacteria. Native zinc-containing AK from Desulfovibrio gigas was purified to homogeneity and crystallized. The crystals diffracted to beyond 1.8 A resolution. Furthermore, cobalt- and iron-containing crystal forms of recombinant AK were also obtained and diffracted to 2.0 and 3.0 A resolution, respectively. Zn(2+)-AK and Fe(2+)-AK crystallized in space group I222 with similar unit-cell parameters, whereas Co(2+)-AK crystallized in space group C2; a monomer was present in the asymmetric unit for both the Zn(2+)-AK and Fe(2+)-AK forms and a dimer was present for the Co(2+)-AK form. The structures of the three metal-bound forms of AK will provide new insights into the role and selectivity of the metal in these enzymes.

Camelid nanobodies raised against an integral membrane enzyme, nitric oxide reductase, Conrath, K., Pereira A. S., Martins C. E., Timoteo C. G., Tavares P., Spinelli S., Kinne J., Flaudrops C., Cambillau C., Muyldermans S., Moura I., Moura J. J., Tegoni M., and Desmyter A. , Protein Sci, Mar, Volume 18, Number 3, p.619-28, (2009) AbstractWebsite

Nitric Oxide Reductase (NOR) is an integral membrane protein performing the reduction of NO to N(2)O. NOR is composed of two subunits: the large one (NorB) is a bundle of 12 transmembrane helices (TMH). It contains a b type heme and a binuclear iron site, which is believed to be the catalytic site, comprising a heme b and a non-hemic iron. The small subunit (NorC) harbors a cytochrome c and is attached to the membrane through a unique TMH. With the aim to perform structural and functional studies of NOR, we have immunized dromedaries with NOR and produced several antibody fragments of the heavy chain (VHHs, also known as nanobodies). These fragments have been used to develop a faster NOR purification procedure, to proceed to crystallization assays and to analyze the electron transfer of electron donors. BIAcore experiments have revealed that up to three VHHs can bind concomitantly to NOR with affinities in the nanomolar range. This is the first example of the use of VHHs with an integral membrane protein. Our results indicate that VHHs are able to recognize with high affinity distinct epitopes on this class of proteins, and can be used as versatile and valuable tool for purification, functional study and crystallization of integral membrane proteins.

Molybdenum induces the expression of a protein containing a new heterometallic Mo-Fe cluster in Desulfovibrio alaskensis, Rivas, M. G., Carepo M. S., Mota C. S., Korbas M., Durand M. C., Lopes A. T., Brondino C. D., Pereira A. S., George G. N., Dolla A., Moura J. J., and Moura I. , Biochemistry, Feb 10, Volume 48, Number 5, p.873-82, (2009) AbstractWebsite

The characterization of a novel Mo-Fe protein (MorP) associated with a system that responds to Mo in Desulfovibrio alaskensis is reported. Biochemical characterization shows that MorP is a periplasmic homomultimer of high molecular weight (260 +/- 13 kDa) consisting of 16-18 monomers of 15321.1 +/- 0.5 Da. The UV/visible absorption spectrum of the as-isolated protein shows absorption peaks around 280, 320, and 570 nm with extinction coefficients of 18700, 12800, and 5000 M(-1) cm(-1), respectively. Metal content, EXAFS data and DFT calculations support the presence of a Mo-2S-[2Fe-2S]-2S-Mo cluster never reported before. Analysis of the available genomes from Desulfovibrio species shows that the MorP encoding gene is located downstream of a sensor and a regulator gene. This type of gene arrangement, called two component system, is used by the cell to regulate diverse physiological processes in response to changes in environmental conditions. Increase of both gene expression and protein production was observed when cells were cultured in the presence of 45 microM molybdenum. Involvement of this system in Mo tolerance of sulfate reducing bacteria is proposed.

2008
Purification, crystallization and preliminary X-ray diffraction analysis of adenosine triphosphate sulfurylase (ATPS) from the sulfate-reducing bacterium Desulfovibrio desulfuricans ATCC 27774, Gavel, O. Y., Kladova A. V., Bursakov S. A., Dias J. M., Texeira S., Shnyrov V. L., Moura J. J., Moura I., Romao M. J., and Trincao J. , Acta Crystallogr Sect F Struct Biol Cryst Commun, Jul 1, Volume 64, Number Pt 7, p.593-5, (2008) AbstractWebsite

Native zinc/cobalt-containing ATP sulfurylase (ATPS; EC 2.7.7.4; MgATP:sulfate adenylyltransferase) from Desulfovibrio desulfuricans ATCC 27774 was purified to homogeneity and crystallized. The orthorhombic crystals diffracted to beyond 2.5 A resolution and the X-ray data collected should allow the determination of the structure of the zinc-bound form of this ATPS. Although previous biochemical studies of this protein indicated the presence of a homotrimer in solution, a dimer was found in the asymmetric unit. Elucidation of this structure will permit a better understanding of the role of the metal in the activity and stability of this family of enzymes.

2006
Modelling the electron-transfer complex between aldehyde oxidoreductase and flavodoxin, Krippahl, Ludwig, Palma Nuno P., Moura Isabel, and Moura Jose J. G. , European Journal of Inorganic Chemistry, Oct 2, Number 19, p.3835-3840, (2006) AbstractWebsite

Three-dimensional protein structures of the xanthine oxidase family show different solutions for the problem of transferring electrons between the flavin adenine dinucleotide (FAD) group and the molybdenum cofactor. In xanthine oxidase all the cofactors he within domains of the same protein chain, whereas in CO dehydrogenase the Fe-S centres, FAD and Mo cofactors are enclosed in separate chains and the enzyme exists as a stable complex of all three. In aldehyde oxidore-ductase, only Fe-S and Mo co-factors are present in a single protein chain. Flavodoxin is docked to aldehyde oxidoreductase to mimic the flavin component on the intramolecular electron transfer chain of aanthine oxidase and CO dehydrogenase and, remarkably, the main features of the electron-transfer pathway are observed.

2005
Interactions of vanadium(V)-citrate complexes with the sarcoplasmic reticulum calcium pump, Aureliano, M., Tiago T., Gandara R. M., Sousa A., Moderno A., Kaliva M., Salifoglou A., Duarte R. O., and Moura J. J. , J Inorg Biochem, Dec, Volume 99, Number 12, p.2355-61, (2005) AbstractWebsite

Among the biotargets interacting with vanadium is the calcium pump from the sarcoplasmic reticulum (SR). To this end, initial research efforts were launched with two vanadium(V)-citrate complexes, namely (NH(4))(6)[V(2)O(4)(C(6)H(4)O(7))(2)].6H(2)O and (NH(4))(6)[V(2)O(2)(O(2))(2)(C(6)H(4)O(7))(2)].4H(2)O, potentially capable of interacting with the SR calcium pump by combining kinetic studies with (51)V NMR spectroscopy. Upon dissolution in the reaction medium (concentration range: 4-0.5mM), both vanadium(V):citrate (VC) and peroxovanadium(V):citrate (PVC) complexes are partially converted into vanadate oligomers. A 1mM solution of the PVC complex, containing 184microM of the PVC complex, 94microM oxoperoxovanadium(V) (PV) species, 222microM monomeric (V1), 43microM dimeric (V2) and 53microM tetrameric (V4) species, inhibits Ca(2+) accumulation by 75 %, whereas a solution of the VC complex of the same vanadium concentration, containing 98microM of the VC complex, 263microM monomeric (V1), 64microM dimeric (V2) and 92microM tetrameric (V4) species inhibits the calcium pump activity by 33 %. In contrast, a 1 mM metavanadate solution, containing 460microM monomeric (V1), 90.2microM dimeric (V2) and 80microM tetrameric (V4) species, has no effect on Ca(2+) accumulation. The NMR signals from the VC complex (-548.0ppm), PVC complex (-551.5ppm) and PV (-611.1ppm) are broadened upon SR vesicle addition (2.5mg/ml total protein). The relative order for the half width line broadening of the NMR signals, which reflect the interaction with the protein, was found to be V4>PVC>VC>PV>V2=V1=1, with no effect observed for the V1 and V2 signals. Putting it all together the effects of two vanadium(V)-citrate complexes on the modulation of calcium accumulation and ATP hydrolysis by the SR calcium pump reflected the observed variable reactivity into the nature of key species forming upon dissolution of the title complexes in the reaction media.

Synechocystis ferredoxin/ferredoxin-NADP(+)-reductase/NADP+ complex: Structural model obtained by NMR-restrained docking, Palma, P. N., Lagoutte B., Krippahl L., Moura J. J., and Guerlesquin F. , FEBS Lett, Aug 29, Volume 579, Number 21, p.4585-90, (2005) AbstractWebsite

Ferredoxin (Fd) and ferredoxin-NADP(+)-reductase (FNR) are two terminal physiological partners of the photosynthetic electron transport chain. Based on a nuclear magnetic resonance (NMR)-restrained-docking approach, two alternative structural models of the Fd-FNR complex in the presence of NADP+ are proposed. The protein docking simulations were performed with the software BiGGER. NMR titration revealed a 1:1 stoichiometry for the complex and allowed the mapping of the interacting residues at the surface of Fd. The NMR chemical shifts were encoded into distance constraints and used with theoretically calculated electronic coupling between the redox cofactors to propose experimentally validated docked complexes.

2003
A further investigation of the cytochrome b5-cytochrome c complex, Banci, L., Bertini I., Felli I. C., Krippahl L., Kubicek K., Moura J. J., and Rosato A. , J Biol Inorg Chem, Sep, Volume 8, Number 7, p.777-86, (2003) AbstractWebsite

The interaction of reduced rabbit cytochrome b(5) with reduced yeast iso-1 cytochrome c has been studied through the analysis of (1)H-(15)N HSQC spectra, of (15)N longitudinal ( R(1)) and transverse ( R(2)) relaxation rates, and of the solvent exchange rates of protein backbone amides. For the first time, the adduct has been investigated also from the cytochrome c side. The analysis of the NMR data was integrated with docking calculations. The result is that cytochrome b(5) has two negative patches capable of interacting with a single positive surface area of cytochrome c. At low protein concentrations and in equimolar mixture, two different 1:1 adducts are formed. At high concentration and/or with excess cytochrome c, a 2:1 adduct is formed. All the species are in fast exchange on the scale of differences in chemical shift. By comparison with literature data, it appears that the structure of one 1:1 adduct changes with the origin or primary sequence of cytochrome b(5).

Electron transfer complexes of cytochrome c peroxidase from Paracoccus denitrificans containing more than one cytochrome, Pettigrew, G. W., Pauleta S. R., Goodhew C. F., Cooper A., Nutley M., Jumel K., Harding S. E., Costa C., Krippahl L., Moura I., and Moura J. , Biochemistry, Oct 21, Volume 42, Number 41, p.11968-81, (2003) AbstractWebsite

According to the model proposed in previous papers [Pettigrew, G. W., Prazeres, S., Costa, C., Palma, N., Krippahl, L., and Moura, J. J. (1999) The structure of an electron-transfer complex containing a cytochrome c and a peroxidase, J. Biol. Chem. 274, 11383-11389; Pettigrew, G. W., Goodhew, C. F., Cooper, A., Nutley, M., Jumel, K., and Harding, S. E. (2003) Electron transfer complexes of cytochrome c peroxidase from Paracoccus denitrificans, Biochemistry 42, 2046-2055], cytochrome c peroxidase of Paracoccus denitrificans can accommodate horse cytochrome c and Paracoccus cytochrome c(550) at different sites on its molecular surface. Here we use (1)H NMR spectroscopy, analytical ultracentrifugation, molecular docking simulation, and microcalorimetry to investigate whether these small cytochromes can be accommodated simultaneously in the formation of a ternary complex. The pattern of perturbation of heme methyl and methionine methyl resonances in binary and ternary solutions shows that a ternary complex can be formed, and this is confirmed by the increase in the sedimentation coefficient upon addition of horse cytochrome c to a solution in which cytochrome c(550) fully occupies its binding site on cytochrome c peroxidase. Docking experiments in which favored binary solutions of cytochrome c(550) bound to cytochrome c peroxidase act as targets for horse cytochrome c and the reciprocal experiments in which favored binary solutions of horse cytochrome c bound to cytochrome c peroxidase act as targets for cytochrome c(550) show that the enzyme can accommodate both cytochromes at the same time on adjacent sites. Microcalorimetric titrations are difficult to interpret but are consistent with a weakened binding of horse cytochrome c to a binary complex of cytochrome c peroxidase and cytochrome c(550) and binding of cytochrome c(550) to the cytochrome c peroxidase that is affected little by the presence of horse cytochrome c in the other site. The presence of a substantial capture surface for small cytochromes on the cytochrome c peroxidase has implications for rate enhancement mechanisms which ensure that the two electrons required for re-reduction of the enzyme after reaction with hydrogen peroxide are delivered efficiently.

Modeling protein complexes with BiGGER, Krippahl, L., Moura J. J., and Palma P. N. , Proteins, Jul 1, Volume 52, Number 1, p.19-23, (2003) AbstractWebsite

This article describes the method and results of our participation in the Critical Assessment of PRediction of Interactions (CAPRI) experiment, using the protein docking program BiGGER (Bimolecular complex Generation with Global Evaluation and Ranking) (Palma et al., Proteins 2000;39:372-384). Of five target complexes (CAPRI targets 2, 4, 5, 6, and 7), only one was successfully predicted (target 6), but BiGGER generated reasonable models for targets 4, 5, and 7, which could have been identified if additional biochemical information had been available.

Electrochemical studies on small electron transfer proteins using membrane electrodes, dos Santos, M. M. C., de Sousa P. M. P., Goncalves M. L. S., Krippahl L., Moura J. J. G., Lojou E., and Bianco P. , Journal of Electroanalytical Chemistry, Jan 16, Volume 541, p.153-162, (2003) AbstractWebsite

Membrane electrodes (ME) were constructed using gold, glassy carbon and pyrolytic graphite supports and a dialysis membrane, and used to study the electrochemical behavior of small size electron transfer proteins: monohemic cytochrome c(522) from Pseudomonas nautica and cytochrome c(533) as well as rubredoxin from Desulfovibrio vulgaris. Different electrochemical techniques were used including cyclic voltammetry (CV), square wave voltammetry (SW) and differential pulse voltammetry (DP). A direct electrochemical response was obtained in all cases except with rubredoxin where a facilitator was added to the protein solution entrapped between the membrane and the electrode surface. Formal potentials and heterogeneous charge transfer rate constants were determined from the voltammetric data. The influence of the ionic strength and the pH of the medium on the electrochemical response at the ME were analyzed. The benefits from the use of the ME in protein electrochemistry and its role in modulating the redox behavior are analyzed. A critical comparison is presented with data obtained at non-MEs. Finally, the interactions that must be established between the proteins and the electrode surfaces are discussed, thereby modeling molecular interactions that occur in biological systems. (C) 2002 Elsevier Science B.V. All rights reserved.

Kinetic behavior of Desulfovibrio gigas aldehyde oxidoreductase encapsulated in reverse micelles, Andrade, S. L., Brondino C. D., Kamenskaya E. O., Levashov A. V., and Moura J. J. , Biochem Biophys Res Commun, Aug 15, Volume 308, Number 1, p.73-8, (2003) AbstractWebsite

We report the kinetic behavior of the enzyme aldehyde oxidoreductase (AOR) from the sulfate reducing bacterium Desulfovibrio gigas (Dg) encapsulated in reverse micelles of sodium bis-(2-ethylhexyl) sulfosuccinate in isooctane using benzaldehyde, octaldehyde, and decylaldehyde as substrates. Dg AOR is a 200-kDa homodimeric protein that catalyzes the conversion of aldehydes to carboxylic acids. Ultrasedimentation analysis of Dg AOR-containing micelles showed the presence of 100-kDa molecular weight species, confirming that the Dg AOR subunits can be dissociated. UV-visible spectra of encapsulated Dg AOR are indistinguishable from the enzyme spectrum in solution, suggesting that both protein fold and metal cofactor are kept intact upon encapsulation. The catalytic constant (k(cat)) profile as a function of the micelle size W(0) (W(0)=[H(2)O]/[AOT]) using benzaldehyde as substrate showed two bell-shaped activity peaks at W(0)=20 and 26. Furthermore, enzymatic activity for octaldehyde and decylaldehyde was detected only in reverse micelles. Like for the benzaldehyde kinetics, two peaks with both similar k(cat) values and W(0) positions were obtained. EPR studies using spin-labeled reverse micelles indicated that octaldehyde and benzaldehyde are intercalated in the micelle membrane. This suggests that, though Dg AOR is found in the cytoplasm of bacterial cells, the enzyme may catalyze the reaction of substrates incorporated into a cell membrane.

2000
Redox potential measurements of the Mycobacterium tuberculosis heme protein KatG and the isoniazid-resistant enzyme KatG(S315T): insights into isoniazid activation, Wengenack, N. L., Lopes H., Kennedy M. J., Tavares P., Pereira A. S., Moura I., Moura J. J., and Rusnak F. , Biochemistry, Sep 19, Volume 39, Number 37, p.11508-13, (2000) AbstractWebsite

Mycobacterium tuberculosis KatG is a multifunctional heme enzyme responsible for activation of the antibiotic isoniazid. A KatG(S315T) point mutation is found in >50% of isoniazid-resistant clinical isolates. Since isoniazid activation is thought to involve an oxidation reaction, the redox potential of KatG was determined using cyclic voltammetry, square wave voltammetry, and spectroelectrochemical titrations. Isoniazid activation may proceed via a cytochrome P450-like mechanism. Therefore, the possibility that substrate binding by KatG leads to an increase in the heme redox potential and the possibility that KatG(S315T) confers isoniazid resistance by altering the redox potential were examined. Effects of the heme spin state on the reduction potentials of KatG and KatG(S315T) were also determined. Assessment of the Fe(3+)/Fe(2+) couple gave a midpoint potential of ca. -50 mV for both KatG and KatG(S315T). In contrast to cytochrome P450s, addition of substrate had no significant effect on either the KatG or KatG(S315T) redox potential. Conversion of the heme to a low-spin configuration resulted in a -150 to -200 mV shift of the KatG and KatG(S315T) redox potentials. These results suggest that isoniazid resistance conferred by KatG(S315T) is not mediated through changes in the heme redox potential. The redox potentials of isoniazid were also determined using cyclic and square wave voltammetry, and the results provide evidence that the ferric KatG and KatG(S315T) midpoint potentials are too low to promote isoniazid oxidation without formation of a high-valent enzyme intermediate such as compounds I and II or oxyferrous KatG.

Neelaredoxin, an iron-binding protein from the syphilis spirochete, Treponema pallidum, is a superoxide reductase, Jovanovic, T., Ascenso C., Hazlett K. R., Sikkink R., Krebs C., Litwiller R., Benson L. M., Moura I., Moura J. J., Radolf J. D., Huynh B. H., Naylor S., and Rusnak F. , J Biol Chem, Sep 15, Volume 275, Number 37, p.28439-48, (2000) AbstractWebsite

Treponema pallidum, the causative agent of venereal syphilis, is a microaerophilic obligate pathogen of humans. As it disseminates hematogenously and invades a wide range of tissues, T. pallidum presumably must tolerate substantial oxidative stress. Analysis of the T. pallidum genome indicates that the syphilis spirochete lacks most of the iron-binding proteins present in many other bacterial pathogens, including the oxidative defense enzymes superoxide dismutase, catalase, and peroxidase, but does possess an orthologue (TP0823) for neelaredoxin, an enzyme of hyperthermophilic and sulfate-reducing anaerobes shown to possess superoxide reductase activity. To analyze the potential role of neelaredoxin in treponemal oxidative defense, we examined the biochemical, spectroscopic, and antioxidant properties of recombinant T. pallidum neelaredoxin. Neelaredoxin was shown to be expressed in T. pallidum by reverse transcriptase-polymerase chain reaction and Western blot analysis. Recombinant neelaredoxin is a 26-kDa alpha(2) homodimer containing, on average, 0.7 iron atoms/subunit. Mossbauer and EPR analysis of the purified protein indicates that the iron atom exists as a mononuclear center in a mixture of high spin ferrous and ferric oxidation states. The fully oxidized form, obtained by the addition of K(3)(Fe(CN)(6)), exhibits an optical spectrum with absorbances at 280, 320, and 656 nm; the last feature is responsible for the protein's blue color, which disappears upon ascorbate reduction. The fully oxidized protein has a A(280)/A(656) ratio of 10.3. Enzymatic studies revealed that T. pallidum neelaredoxin is able to catalyze a redox equilibrium between superoxide and hydrogen peroxide, a result consistent with it being a superoxide reductase. This finding, the first description of a T. pallidum iron-binding protein, indicates that the syphilis spirochete copes with oxidative stress via a primitive mechanism, which, thus far, has not been described in pathogenic bacteria.

Heteronuclear NMR and soft docking: an experimental approach for a structural model of the cytochrome c553-ferredoxin complex, Morelli, X., Dolla A., Czjzek M., Palma P. N., Blasco F., Krippahl L., Moura J. J., and Guerlesquin F. , Biochemistry, Mar 14, Volume 39, Number 10, p.2530-7, (2000) AbstractWebsite

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.

BiGGER: a new (soft) docking algorithm for predicting protein interactions, Palma, P. N., Krippahl L., Wampler J. E., and Moura J. J. , Proteins, Jun 1, Volume 39, Number 4, p.372-84, (2000) AbstractWebsite

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.

1999
Biochemical and spectroscopic characterization of overexpressed fuscoredoxin from Escherichia coli, Pereira, A. S., Tavares P., Krebs C., Huynh B. H., Rusnak F., Moura I., and Moura J. J. , Biochem Biophys Res Commun, Jun 24, Volume 260, Number 1, p.209-15, (1999) AbstractWebsite

Fuscoredoxin is a unique iron containing protein of yet unknown function originally discovered in the sulfate reducers of the genus Desulfovibrio. It contains two iron-sulfur clusters: a cubane [4Fe-4S] and a mixed oxo- and sulfido-bridged 4Fe cluster of unprecedented structure. The recent determination of the genomic sequence of Escherichia coli (E. coli) has revealed a homologue of fuscoredoxin in this facultative microbe. The presence of this gene in E. coli raises interesting questions regarding the function of fuscoredoxin and whether this gene represents a structural homologue of the better-characterized Desulfovibrio proteins. In order to explore the latter, an overexpression system for the E. coli fuscoredoxin gene was devised. The gene was cloned from genomic DNA by use of the polymerase chain reaction into the expression vector pT7-7 and overexpressed in E. coli BL21(DE3) cells. After two chromatographic steps a good yield of recombinant protein was obtained (approximately 4 mg of pure protein per liter of culture). The purified protein exhibits an optical spectrum characteristic of the homologue from D. desulfuricans, indicating that cofactor assembly was accomplished. Iron analysis indicated that the protein contains circa 8 iron atoms/molecule which were shown by EPR and Mossbauer spectroscopies to be present as two multinuclear clusters, albeit with slightly altered spectroscopic features. A comparison of the primary sequences of fuscoredoxins is presented and differences on cluster coordination modes are discussed on the light of the spectroscopic data.

The structure of an electron transfer complex containing a cytochrome c and a peroxidase, Pettigrew, G. W., Prazeres S., Costa C., Palma N., Krippahl L., Moura I., and Moura J. J. , J Biol Chem, Apr 16, Volume 274, Number 16, p.11383-9, (1999) AbstractWebsite

Efficient biological electron transfer may require a fluid association of redox partners. Two noncrystallographic methods (a new molecular docking program and 1H NMR spectroscopy) have been used to study the electron transfer complex formed between the cytochrome c peroxidase (CCP) of Paracoccus denitrificans and cytochromes c. For the natural redox partner, cytochrome c550, the results are consistent with a complex in which the heme of a single cytochrome lies above the exposed electron-transferring heme of the peroxidase. In contrast, two molecules of the nonphysiological but kinetically competent horse cytochrome bind between the two hemes of the peroxidase. These dramatically different patterns are consistent with a redox active surface on the peroxidase that may accommodate more than one cytochrome and allow lateral mobility.

1998
The structural origin of nonplanar heme distortions in tetraheme ferricytochromes c3, Ma, J. G., Zhang J., Franco R., Jia S. L., Moura I., Moura J. J., Kroneck P. M., and Shelnutt J. A. , Biochemistry, Sep 8, Volume 37, Number 36, p.12431-42, (1998) AbstractWebsite

Resonance Raman (RR) spectroscopy, molecular mechanics (MM) calculations, and normal-coordinate structural decomposition (NSD) have been used to investigate the conformational differences in the hemes in ferricytochromes c3. NSD analyses of heme structures obtained from X-ray crystallography and MM calculations of heme-peptide fragments of the cytochromes c3 indicate that the nonplanarity of the hemes is largely controlled by a fingerprint peptide segment consisting of two heme-linked cysteines, the amino acids between the cysteines, and the proximal histidine ligand. Additional interactions between the heme and the distal histidine ligand and between the heme propionates and the protein also influence the heme conformation, but to a lesser extent than the fingerprint peptide segment. In addition, factors that influence the folding pattern of the fingerprint peptide segment may have an effect on the heme conformation. Large heme structural differences between the baculatum cytochromes c3 and the other proteins are uncovered by the NSD procedure [Jentzen, W., Ma, J.-G., and Shelnutt, J. A. (1998) Biophys. J. 74, 753-763]. These heme differences are mainly associated with the deletion of two residues in the covalently linked segment of hemes 4 for the baculatum proteins. Furthermore, some of these structural differences are reflected in the RR spectra. For example, the frequencies of the structure-sensitive lines (nu4, nu3, and nu2) in the high-frequency region of the RR spectra are lower for the Desulfomicrobium baculatum cytochromes c3 (Norway 4 and 9974) than for the Desulfovibrio (D.) gigas, D. vulgaris, and D. desulfuricans strains, consistent with a more ruffled heme. Spectral decompositions of the nu3 and nu10 lines allow the assignment of the sublines to individual hemes and show that ruffling, not saddling, is the dominant factor influencing the frequencies of the structure-sensitive Raman lines. The distinctive spectra of the baculatum strains investigated are a consequence of hemes 2 and 4 being more ruffled than is typical of the other proteins.

Oxovanadium(IV) complexes of the dipeptides glycyl-L-aspartic acid, L-aspartylglycine and related ligands; a spectroscopic and potentiometric study, Pessoa, J. C., Gajda T., Gillard R. D., Kiss T., Luz S. M., Moura J. J. G., Tomaz I., Telo J. P., and Torok I. , Journal of the Chemical Society-Dalton Transactions, Nov 7, Number 21, p.3587-3600, (1998) AbstractWebsite

The equilibria in the systems VO2+ + L (L = Gly-L-Asp, L-Asp-Gly, N-acetyl-L-aspartic acid or succinic acid) have been studied at 25 degrees C and 0.2 mol dm(3) K(CI) medium by a combination of potentiometric and spectroscopic methods (ESR, circular dichroism and visible absorption). Formation constants were calculated from pH-metric data with total oxovanadium(Iv) concentrations of(0.6-4) x 10(-3) mol dm(-3) and ligand-to-metal (L:M) ratios of 2-8 (AspGly) or 4-15: 1 (other systems). The position of the Asp residue in the peptide chain affects the co-ordination mode of the ligands: while in the GlyAsp system bis complexes start to form at pH less than 2, for AspGly only 1 : 1 complexes form, with relatively high CD signal. The co-ordination behaviour of N-acetyl-L-aspartic and succinic acids is similar. The results of potentiometric and spectroscopic methods are self consistent. Isomeric structures are discussed for each stoichiometry proposed and the results compared with those for L-aspartic acid and dipeptides with non-coordinating side chains.

Spectroscopic characterization of a novel tetranuclear Fe cluster in an iron-sulfur protein isolated from Desulfovibrio desulfuricans, Tavares, P., Pereira A. S., Krebs C., Ravi N., Moura J. J., Moura I., and Huynh B. H. , Biochemistry, Mar 3, Volume 37, Number 9, p.2830-42, (1998) AbstractWebsite

Mossbauer and EPR spectroscopies were used to characterize the Fe clusters in an Fe-S protein isolated from Desulfovibrio desulfuricans (ATCC 27774). This protein was previously thought to contain hexanuclear Fe clusters, but a recent X-ray crystallographic measurement on a similar protein isolated from Desulfovibrio vulgaris showed that the protein contains two tetranuclear clusters, a cubane-type [4Fe-4S] cluster and a mixed-ligand cluster of novel structure [Lindley et al. (1997) Abstract, Chemistry of Metals in Biological Systems, European Research Conference, Tomar, Portugal]. Three protein samples poised at different redox potentials (as-purified, 40 and 320 mV) were investigated. In all three samples, the [4Fe-4S] cluster was found to be present in the diamagnetic 2+ oxidation state and exhibited typical Mossbauer spectra. The novel-structure cluster was found to be redox active. In the 320-mV and as-purified samples, the cluster is at a redox equilibrium between its fully oxidized and one-electron reduced states. In the 40-mV sample, the cluster is in a two-electron reduced state. Distinct spectral components associated with the four Fe sites of cluster 2 in the three oxidation states were identified. The spectroscopic parameters obtained for the Fe sites reflect different ligand environments, making it possible to assign the spectral components to individual Fe sites. In the fully oxidized state, all four iron ions are high-spin ferric and antiferromagnetically coupled to form a diamagnetic S = 0 state. In the one-electron and two-electron reduced states, the reducing electrons were found to localize, consecutively, onto two Fe sites that are rich in oxygen/nitrogen ligands. Based on the X-ray structure and the Mossbauer parameters, attempts could be made to identify the reduced Fe sites. For the two-electron reduced cluster, EPR and Mossbauer data indicate that the cluster is paramagnetic with a nonzero interger spin. For the one-electron reduced cluster, the data suggest a half-integer spin of 9/2. Characteristic fine and hyperfine parameters for all four Fe sites were obtained. Structural implications and the nature of the spin-coupling interactions are discussed.