Structural evidence for ligand specificity in the binding domain of the human androgen receptor. Implications for pathogenic gene mutations

The crystal structures of the human androgen receptor (hAR) and human progesterone receptor ligand-binding domains in complex with the same ligand metribolone (R1881) have been determined. Both three-dimensional structures show the typical nuclear receptor fold. The change of two residues in the ligand-binding pocket between the human progesterone receptor and hAR is most likely the source for the specificity of R1881 to the hAR. The structural implications of the 14 known mutations in the ligand-binding pocket of the hAR ligand-binding domains associated with either prostate cancer or the partial or complete androgen receptor insensitivity syndrome were analyzed. The effects of most of these mutants could be explained on the basis of the crystal structure.

Crystallization and preliminary x-ray diffraction analysis of beta-cinnamomin, an elicitin secreted by the phytopathogenic fungus Phytophthora cinnamomin

Cinnamomin (CIN) belongs to a family of 10 kDa proteins designated as elicitins. Some of these proteins induce a hypersensitive response in diverse plant species, leading to resistance against fungal and bacterial plant pathogens. CIN was crystallized by the vapour-diffusion method using either ammonium sulfate or polyethyleneglycol (PEG) as precipitants in solutions buffered at around pH 7. These crystals are isomorphous and belong to the triclinic space group, with unit-cell parameters a = 31.69, b = 36. 99, c = 44.09 A, alpha = 76.86, beta = 84.41, gamma = 80.26 degrees. A frozen crystal diffracted X-rays beyond 1.45 A resolution on a synchrotron-radiation source.

Improved modeling of side-chains in proteins with rotamer-based methods: a flexible rotamer model

Side-chain modeling has a widespread application in many current methods for protein tertiary structure determination, prediction, and design. Of the existing side-chain modeling methods, rotamer-based methods are the fastest and most efficient. Classically, a rotamer is conceived as a single, rigid conformation of an amino acid sidechain. Here, we present a flexible rotamer model in which a rotamer is a continuous ensemble of conformations that cluster around the classic rigid rotamer. We have developed a thermodynamically based method for calculating effective energies for the flexible rotamer. These energies have a one-to-one correspondence with the potential energies of the rigid rotamer. Therefore, the flexible rotamer model is completely general and may be used with any rotamer-based method in substitution of the rigid rotamer model. We have compared the performance of the flexible and rigid rotamer models with one side-chain modeling method in particular (the self-consistent mean field theory method) on a set of 20 high quality crystallographic protein structures. For the flexible rotamer model, we obtained average predictions of 85.8% for chi1, 76.5% for chi1+2 and 1.34 A for root-mean-square deviation (RMSD); the corresponding values for core residues were 93.0%, 87.7% and 0.70 A, respectively. These values represent improvements of 7.3% for chi1, 8.1% for chi1+2 and 0.23 A for RMSD over the predictions obtained with the rigid rotamer model under otherwise identical conditions; the corresponding improvements for core residues were 6.9%, 10.5% and 0.43 A, respectively. We found that the predictions obtained with the flexible rotamer model were also significantly better than those obtained for the same set of proteins with another state-of-the-art side-chain placement method in the literature, especially for core residues. The flexible rotamer model represents a considerable improvement over the classic rigid rotamer model. It can, therefore, be used with considerable advantage in all rotamer-based methods commonly applied to protein tertiary structure determination, prediction, and design and also in predictions of free energies in mutational studies.

An iterative structure-assisted approach to sequence alignment and comparative modeling

Correct alignment of the sequence of a target protein with those of homologues of known three-dimensional structure is a key step in comparative modeling. Usually an iterative approach that takes account of the local and overall structural features is required. We describe such an approach that exploits databases of structural alignments of homologous proteins (HOMSTRAD, http:/(/)www-cryst.bioc.cam.ac.uk/ approximately homstrad) and protein superfamilies (CAMPASS, http:/(/)www-cryst.bioc.cam.ac.uk/ approximately campass), in which structure-based alignments are analyzed and formatted with the program JOY (http:/(/)www-cryst.bioc.cam.ac.uk/ approximately joy) to reveal conserved local structural features. The databases facilitate the recognition of a family or superfamily, they assist in the selection of useful parent structures, they are helpful in alignment of the target sequences with the parent set, and are useful for deriving relationships that can be used in validating models. In the iterative approach, a model is constructed on the basis of the proposed sequence alignment and this is then reexpressed in the JOY format and realigned with the parent set. This is repeated until the model and sequence alignment is optimized. We examine the case for comparison and use of multiple structures of family members, rather than a single parent structure. We use the targets attempted by our group in CASP3 to assess the value of such procedures.

Crystal structure of cardosin A, a glycosylated and Arg-Gly-Asp-containing aspartic proteinase from the flowers of Cynara cardunculus L

Aspartic proteinases (AP) have been widely studied within the living world, but so far no plant AP have been structurally characterized. The refined cardosin A crystallographic structure includes two molecules, built up by two glycosylated peptide chains (31 and 15 kDa each). The fold of cardosin A is typical within the AP family. The glycosyl content is described by 19 sugar rings attached to Asn-67 and Asn-257. They are localized on the molecular surface away from the conserved active site and show a new glycan of the plant complex type. A hydrogen bond between Gln-126 and Manbeta4 renders the monosaccharide oxygen O-2 sterically inaccessible to accept a xylosyl residue, therefore explaining the new type of the identified plant glycan. The Arg-Gly-Asp sequence, which has been shown to be involved in recognition of a putative cardosin A receptor, was found in a loop between two beta-strands on the molecular surface opposite the active site cleft. Based on the crystal structure, a possible mechanism whereby cardosin A might be orientated at the cell surface of the style to interact with its putative receptor from pollen is proposed. The biological implications of these findings are also discussed.

The caa3 terminal oxidase of the thermohalophilic bacterium Rhodothermus marinus: a HiPIP:oxygen oxidoreductase lacking the key glutamate of the D-channel

The respiratory chain of the thermohalophilic bacterium Rhodothermus marinus contains a novel complex III and a high potential iron-sulfur protein (HiPIP) as the main electron shuttle (Pereira et al., Biochemistry 38 (1999) 1268-1275 and 1276-1283). In this paper, one of the terminal oxidases expressed in this bacterium is extensively characterised. It is a caa3-type oxidase, isolated with four subunits (apparent molecular masses of 42, 19 and 15 kDa and a C-haem containing subunit of 35 kDa), which has haems of the A(s) type. This oxidase is capable of using TMPD and horse heart cytochrome c as substrates, but has a higher turnover with HiPIP, being the first example of a HiPIP:oxygen oxidoreductase. The oxidase has unusually low reduction potentials of 260 (haem C), 255 (haem A) and 180 mV (haem A3). Subunit I of R. marinus caa3 oxidase has an overall significant homology with the subunits I of the COX type oxidases, namely the metal binding sites and most residues considered to be functionally important for proton uptake and pumping (K- and D-channels). However, a major difference is present: the putative essential glutamate (E278 in Paraccocus denitrificans) of the D-channel is missing in the R. marinus oxidase. Homology modelling of the R. marinus oxidase shows that the phenol group of a tyrosine residue may occupy a similar spatial position as the glutamate carboxyl, in relation to the binuclear centre. Moreover, sequence comparisons reveal that several enzymes lacking that glutamate have a conserved substitution pattern in helix VI: -YSHPXV- instead of -XGHPEV-. These observations are discussed in terms of the mechanisms for proton uptake and it is suggested that, in these enzymes, tyrosine may play the role of the glutamate in the proton channel.

Improvement of side-chain modeling in proteins with the self-consistent mean field theory method based on an analysis of the factors influencing prediction

With the objective of improving side-chain conformation prediction, we have analyzed the influence of various factors on prediction by the Self-Consistent Mean Field Theory method, applied to a set of high resolution x-ray protein structure models. These factors may be classed as variations in the mean field optimization protocol, variations in the potential energy function, and variations in rotamer library completeness. We have developed an optimization protocol that consistently reached lower mean field conformational free energies than two other protocols. This protocol led to an important improvement in prediction. We observed a major improvement in prediction with two more detailed van der Waals parameter sets, which we found to be due mainly to the introduction of scaling of 1-4 interactions. In a comparison of two knowledge-based rotamer libraries of considerably different size, we observed an unexpected decrease in prediction with an increase in library completeness. However, when we introduced a torsion potential term in the potential energy function, we found an important increase in average prediction and in the prediction of almost all residue types with a more complete rotamer set. The two knowledge-based rotamer libraries now became equivalent in terms of average prediction. The results we obtained in an analysis of the effect of the introduction of an additional electrostatic term in the potential energy function were largely inconclusive. However, we found a small increase in average prediction for an electrostatic potential term with a fixed dielectric constant of 15. The combined effect of all the factors we analyzed in this study resulted in average prediction accuracies of 79.9% for X1, 68.1% for X1 + 2, and 1.590 A for global rms deviation (RMSD); the corresponding values for core residues were 88.2%, 78.6%, and 1.171 A. These values represent improvements in average prediction of 6.5% for X1, 9.1% for X1 + 2, and 0.163 A for global RMSD over the original conditions; the corresponding improvements in the core were 5.9%, 9.0%, and 0.180 A, respectively.

Crystallization and preliminary diffraction data analysis of both single and pseudo-merohedrally twinned crystals of rubredoxin oxygen oxidoreductase from Desulfovibrio gigas

Crystals of rubredoxin oxygen oxidoreductase have been obtained and characterized. They belong to space group P2(1)2(1)2, with unit-cell dimensions a = 88.24 (15), b = 101.25 (7), c = 90.80 (3) A. The homodimer (86 kDa) in the asymmetric unit is related by a non-crystallographic twofold rotation axis parallel to the ab 'diagonal' direction, as shown by the self-rotation maximum in the section with chi = 180 degrees. This pseudo-crystallographic symmetry element was also found to be the twinning axis of pseudo-merohedrally twinned crystals, leading to apparent pseudo-tetragonal P42(1)2 crystal symmetry.

Nine-haem cytochrome c from Desulfovibrio desulfuricans ATCC 27774:primary sequence determination, crystallographic refinement at 1.8 and modelling studies of its interaction with the tetrahaem cytochrome c3

A monomeric nine-haem cytochrome c (9Hcc) with 292 amino acid residues was isolated from cells of the sulfate- and nitrate-reducing bacterium Desulfovibrio desulfuricans ATCC 27774 grown under both nitrate- and sulfate-respiring conditions. The nucleotide sequence encoding the 292 residues was determined, allowing the correction of about 10% of the previous primary structure, determined from 1.8 A electron density maps. The refinement at 1.8 A resolution of the structural model was completed, giving an R-value of 16.5%. The nine haem groups are arranged into two tetrahaem clusters, located at both ends of the molecule, with Fe-Fe distances and local protein fold very similar to tetrahaem cytochromes c3, and the extra haem is located asymmetrically between the two regions. The new primary sequence determination confirmed the 39% sequence homology found between this cytochrome and the C-terminal region (residues 229-514) of the high-molecular-weight cytochrome c (Hmc) from D. vulgaris Hildenborough, providing strong evidence of structural similarity between 9Hcc and the C-terminal region of Hmc. The interaction between 9Hcc and the tetrahaem cytochrome c3 from the same organism was studied by modelling methods, and the results suggest that a specific interaction is possible between haem 4 of tetrahaem cytochrome c3 and haem 1 or haem 2 of 9Hcc, in agreement with previous kinetic experiments which showed the catalytic effect of the tetrahaem cytochrome c3 upon the reduction of 9Hcc by the [NiFe] hydrogenase from D. desulfuricans ATCC 27774. These studies suggest a role for 9Hcc as part of the assembly of redox proteins involved in recycling the molecular hydrogen released by the cell as a result of substrate oxidation.

Ab initio structure solution of a dimeric cytochrome c3 from Desulfovibrio gigas containing disulfide bridges

The 1.2 A resolution crystal structure of the 29 kDa di-tetrahaem cytochrome c3 from the sulfate reducing bacterium Desulfovibrio gigas was solved by ab initio methods, making this the largest molecule to be solved by this procedure. The actual refined model of the cysteine-linked dimeric molecule reveals that this molecule is very similar to the non-covalently linked symmetrical dimer of the di-tetrahaem cytochrome c3 from Desulfomicrobium norvegicum. Each monomer has the typical polypeptide fold, haem arrangement and iron coordination found for the tetrahaem cytochrome c3 molecules. The interface between the covalently linked monomers in the asymmetric unit of the crystal shows a pseudo two-fold arrangement, disturbed from symmetry by crystal packing forces. The fact that D. gigas contains a dimeric tetrahaem cytochrome with solvent accessible disulfide bridges and that this cytochrome specifically couples hydrogen oxidation to thiosulfate reduction in bacterial extracts provides an interesting aspect related to disulfide exchange reactions in this microorganism.

The primary and three-dimensional structures of a nine-haem cytochrome c from Desulfovibrio desulfuricans ATCC 27774 reveal a new member of the Hmc family

BACKGROUND

Haem-containing proteins are directly involved in electron transfer as well as in enzymatic functions. The nine-haem cytochrome c (9Hcc), previously described as having 12 haem groups, was isolated from cells of Desulfovibrio desulfuricans ATCC 27774, grown under both nitrate- and sulphate-respiring conditions.

RESULTS

Models for the primary and three-dimensional structures of this cytochrome, containing 292 amino acid residues and nine haem groups, were derived using the multiple wavelength anomalous dispersion phasing method and refined using 1.8 A diffraction data to an R value of 17.0%. The nine haem groups are arranged into two tetrahaem clusters, with Fe-Fe distances and local protein fold similar to tetrahaem cytochromes c3, while the extra haem is located asymmetrically between the two clusters.

CONCLUSIONS

This is the first known three-dimensional structure in which multiple copies of a tetrahaem cytochrome c3-like fold are present in the same polypeptide chain. Sequence homology was found between this cytochrome and the C-terminal region (residues 229-514) of the high molecular weight cytochrome c from Desulfovibrio vulgaris Hildenborough (DvH Hmc). A new haem arrangement in domains III and IV of DvH Hmc is proposed. Kinetic experiments showed that 9Hcc can be reduced by the [NiFe] hydrogenase from D. desulfuricans ATCC 27774, but that this reduction is faster in the presence of tetrahaem cytochrome c3. As Hmc has never been found in D. desulfuricans ATCC 27774, we propose that 9Hcc replaces it in this organism and is therefore probably involved in electron transfer across the membrane.

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.

Crystallization and preliminary X-ray crystallographic studies of the plant aspartic proteinase cardosin A

The plant aspartic proteinase cardosin A was crystallized using vapour diffusion. Crystals belong to the monoclinic space group C2, cell dimensions a = 116.9 (2), b = 87.2 (8), c = 81.3 (1) A, beta = 104.4 (4) degrees, and contain two molecules in the asymmetric unit related by a non-crystallographic twofold axis. Diffraction data were collected at room temperature with radiation from a synchrotron source up to 2.85 A resolution. When the crystals were flash cooled to 110 K in a nitrogen stream the same resolution limit could also be obtained on a rotating-anode source. Recently, synchrotron radiation together with flash cooling led to an improvement of the diffraction data to 1.72 A resolution.

Molecular dynamics simulation of cytochrome c3: studying the reduction processes using free energy calculations

The tetraheme cytochrome c3 from Desulfovibrio vulgaris Hildenborough is studied using molecular dynamics simulation studies in explicit solvent. The high heme content of the protein, which has its core almost entirely made up of c-type heme, presents specific problems in the simulation. Instability in the structure is observed in long simulations above 1 ns, something that does not occur in a monoheme cytochrome, suggesting problems in heme parametrization. Given these stability problems, a partially restrained model, which avoids destruction of the structure, was created with the objective of performing free energy calculations of heme reduction, studies that require long simulations. With this model, the free energy of reduction of each individual heme was calculated. A correction in the long-range electrostatic interactions of charge groups belonging to the redox centers had to be made in order to make the system physically meaningful. Correlation is obtained between the calculated free energies and the experimental data for three of four hemes. However, the relative scale of the calculated energies is different from the scale of the experimental free energies. Reasons for this are discussed. In addition to the free energy calculations, this model allows the study of conformational changes upon reduction. Even if the precise details of the structural changes that take place in this system upon individual heme reduction are probably out of the reach of this study, it appears that these structural changes are small, similarly to what is observed for other redox proteins. This does not mean that their effect is minor, and one example is the conformational change observed in propionate D from heme I when heme II becomes reduced. A motion of this kind could be the basis of the experimentally observed cooperativity effects between heme reduction, namely positive cooperativity.

The glycosylation of the aspartic proteinases from barley (Hordeum vulgare L.) and cardoon (Cynara cardunculus L.)

Plant aspartic proteinases characterised at the molecular level contain one or more consensus N-glycosylation sites [Runeberg-Roos, P., Tŏrmäkangas, K. & Ostman, A. (1991) Eur. J. Biochem. 202, 1021-1027; Asakura, T., Watanabe, H., Abe, K. & Arai, S. (1995) Eur. J. Biochem, 232, 77-83; Veríssimo, P., Faro, C., Moir, A. J. G., Lin, Y., Tang, J. & Pires, E. (1996) Eur. J. Biochem. 235, 762-768]. We found that the glycosylation sites are occupied for the barley (Hordeum vulgare L.) aspartic proteinase (Asn333) and the cardoon (Cynara cardunculus L.) aspartic proteinase, cardosin A (Asn70 and Asn363). The oligosaccharides from each site were released from peptide pools by enzymatic hydrolysis with peptide-N-glycanase A or by hydrazinolysis and their structures were determined by exoglycosidase sequencing combined with matrix-assisted laser desorption/ionization time of flight mass spectrometry. It was observed that 6% of the oligosaccharides from the first glycosylation site of cardosin A are of the oligomannose type. Modified type glycans with proximal Fuc and without Xyl account for about 82%, 14% and 3% of the total oligosaccharides from the first and the second glycosylation sites of cardosin A and from H. vulgare aspartic proteinase, respectively. Oligosaccharides with Xyl but without proximal Fuc were only detected in the latter proteinase (4%). Glycans with proximal Fuc and Xyl account for 6%, 86% and 92% of total oligosaccharides from the first and second glycosylation sites of cardosin A and from H. vulgare aspartic proteinase, respectively.

Preliminary crystallographic analysis and further characterization of a dodecaheme cytochrome c from Desulfovibrio desulfuricans ATCC 27774

Dodecaheme cytochrome c has been purified from Desulfovibrio (D.) desulfuricans ATCC 27774 cells grown under both nitrate and sulfate-respiring conditions. Therefore, it is likely to play a role in the electron-transfer system of both respiratory chains. Its molecular mass (37768 kDa) was determined by electrospray mass spectrometry. Its first 39 amino acids were sequenced and a motif was found between amino acids 32 and 37 that seems to exist in all the cytochromes of the c(3) type from sulfate-reducing bacteria sequenced at present. The midpoint redox potentials of this cytochrome were estimated to be -68, -120, -248 and -310 mV. Electron paramagnetic resonance spectroscopy of the oxidized cytochrome shows several low-spin components with a g(max) spreading from 3.254 to 2.983. Two crystalline forms were obtained by vapour diffusion from a solution containing 2% PEG 6000 and 0.25-0.75 M acetate buffer pH = 5.5. Both crystals belong to monoclinic space groups: one is P2(1), with a = 61.00, b = 106.19, c = 82.05 A, beta = 103.61 degrees, and the other is C2 with a = 152.17, b = 98.45, c = 89.24 A, beta = 119.18 degrees. Density measurements of the P2(1) crystals suggest that there are two independent molecules in the asymmetric unit. Self-rotation function calculations indicate, in both crystal forms, the presence of a non-crystallographic axis perpendicular to the crystallographic twofold axis. This result and the calculated values for the volume per unit molecular weight of the C2 crystals suggest the presence of two or four molecules in the asymmetric unit.