X-ray crystallography and computational docking for the detection and development of protein-ligand interactions.

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disorder characterised by the selective dysfunction and death of the upper and lower motor neurons. Median survival rates are between 3 and 5 years after diagnosis. Mutations in the gene encoding Cu/Zn superoxide dismutase (SOD1) have been linked to a subset of familial forms of ALS (fALS). Herein, we describe a fragment- based drug discovery (FBDD) approach for the investigation of small molecule binding sites in SOD1. X-ray crystallography has been used as the primary screening method and has been shown to directly detect protein-ligand interactions which cannot be unambiguously identified using other biophysical methods. The structural requirements for effective binding at Trp32 are detailed for a series of quinazoline-containing compounds. The investigation of an additional site that binds a range of catecholamines and the use of computational modelling to assist fragment evolution is discussed. This study also highlights the importance of ligand solubility for successful Xray crystallographic campaigns in lead compound design.

Crystal structure of human prion protein bound to a therapeutic antibody.

Prion infection is characterized by the conversion of host cellular prion protein (PrP(C)) into disease-related conformers (PrP(Sc)) and can be arrested in vivo by passive immunization with anti-PrP monoclonal antibodies. Here, we show that the ability of an antibody to cure prion-infected cells correlates with its binding affinity for PrP(C) rather than PrP(Sc). We have visualized this interaction at the molecular level by determining the crystal structure of human PrP bound to the Fab fragment of monoclonal antibody ICSM 18, which has the highest affinity for PrP(C) and the highest therapeutic potency in vitro and in vivo. In this crystal structure, human PrP is observed in its native PrP(C) conformation. Interactions between neighboring PrP molecules in the crystal structure are mediated by close homotypic contacts between residues at position 129 that lead to the formation of a 4-strand intermolecular beta-sheet. The importance of this residue in mediating protein-protein contact could explain the genetic susceptibility and prion strain selection determined by polymorphic residue 129 in human prion disease, one of the strongest common susceptibility polymorphisms known in any human disease.

XAFS study of the high-affinity copper-binding site of human PrP(91-231) and its low-resolution structure in solution.

Here, we describe the structure of a C-terminal high-affinity copper-binding site within a truncated recombinant human PrP containing residues 91-231, which lacks the octapeptide repeat region. We show that at least two extra co-ordinating groups are involved in binding this copper(II) ion in conjunction with histidine residues 96 and 111 in a region of the molecule known to be critical in conferring strain type. In addition, using X-ray solution scattering, a low-resolution shape of PrP(91-231) is provided. The restored molecular envelope is consistent with the picture where the N-terminal segment, residues 91-120, extends out from the previously known globular domain containing residues 121-231.

X-ray structure of a blue copper nitrite reductase at high pH and in copper-free form at 1.9 A resolution.

Copper-containing nitrite reductases possess a trimeric structure where the catalytic Cu site, located at the monomer-monomer interface, resembles the catalytic sites of a number of Zn enzymes. Nitrite reductase from Alcaligenes xylosoxidans has optimum activity at pH 5.2 which decreases to a negligible level at pH 8. The structure of this nitrite reductase has previously been determined at pH 4.6. It has now been crystallized under new conditions at pH 8.5. Its crystallographic structure provides a structural explanation for the greatly reduced activity of the enzyme at high pH. Characterization of overexpressed protein in solution by EXAFS suggested that the protein lacked Cu in the catalytic type 2 Cu site and that the site was most probably occupied by Zn. Using the anomalous signals from Cu and Zn, the crystal structure revealed that the expressed protein was devoid of Cu in the catalytic site and that only a trace amount (<10%) of Zn was present at this site in the crystal. Despite the close structural similarity of the catalytic site to a number of Zn enzymes, these data suggest that Zn, if it binds at the catalytic copper site, binds weakly in nitrite reductase.

Structure of the M148Q mutant of rusticyanin at 1.5 A: a model for the copper site of stellacyanin.

The small blue copper protein rusticyanin from Thiobacillus ferrooxidans contains a type 1 Cu centre with a single axial ligand, Met148, which together with the His-Cys-His trigonal planar ligands produces a distorted trigonal pyramidal coordination geometry to copper. Type 1 Cu sites are found in cupredoxins and several multicopper proteins, including oxidases and nitrite reductases. The role of the axial ligand has been extensively debated in terms of its function in the fine tuning of the redox potential and spectroscopic properties of type 1 Cu sites. Numerous mutations of the Met ligand in azurins have been studied, but interpretation of the results has been complicated by the presence of the additional carbonyl oxygen ligand from Gly45, a neighbouring residue to the coordinating His46. The importance of the axial ligand has been further emphasized by the finding that the type 1 centre in Rhus vernicifera stellacyanin, with the lowest redox potential in a type 1 Cu site of 184 mV, has Gln as the axial ligand, whilst fungal laccase and ceruloplasmin, which have redox potentials of 550-800 mV, have a Leu in this position. Here, the crystal structure of the M148Q mutant of rusticyanin at 1.5 A resolution is presented. This is a significantly higher resolution than that of the structures of native rusticyanin. In addition, the M148Q structure is that of the oxidized protein while the native structures to date are of the reduced protein. The mutant protein crystallizes with two molecules per asymmetric unit, in contrast to the one present in the native crystal form. This mutant's redox potential (550 mV at pH 3.2) is lowered compared with that of the native protein ( approximately 670 mV at pH 3.2) by about 120 mV. The type 1 Cu site of M148Q closely mimics the structural characteristics of the equivalent site in non-glycosylated cucumber stellacyanin (redox potential approximately 260 mV) and, owing to the absence in rusticyanin of the fifth, carbonyl ligand present in azurin, may provide a better model for the R. vernicifera stellacyanin (redox potential approximately 184 mV) type 1 Cu site, which also lacks the fifth ligand. Furthermore, the presence of two molecules in the asymmetric unit cell indicates a potential binding region of the redox partners.

Electronic characterization of the oxidized state of the blue copper protein rusticyanin by 1H NMR: is the axial methionine the dominant influence for the high redox potential?

The oxidized state of rusticyanin, the blue copper protein with the highest redox potential in its class, has been investigated through (1)H nuclear magnetic resonance applied to its cobalt(II) derivative. The assignment of the protons belonging to the coordinated residues has been performed. Many other amino acids situated in the vicinity of the metal ion, including six hydrophobic residues (isoleucine140 and five phenylalanines) have also been identified. The orientation of the main axes of the magnetic susceptibility tensor for the cobalt(II)-rusticyanin as well as its axial, Deltachi(ax), and rhombic, Deltachi(rh), magnetic susceptibility anisotropy components have been determined. A comparison of the present results with those previously obtained for cobalt(II)azurin [Donaire, A., Salgado, J., Moratal, J. M. (1998) Biochemistry 37, 8659-8673] allows us to provide further insights into the reasons for the high redox potential of this protein. According to our results, the interaction between the metal ion and the thioether Sdelta of the axial methionine is not as influential as the strong destabilizing effect that the hydrophobic residues close to the metal ion undergo in the oxidized state.

Evidence for the selective population of FeMo cofactor sites in MoFe protein and its molecular recognition by the Fe protein in transition state complex analogues of nitrogenase.

We have collected synchrotron x-ray solution scattering data for the MoFe protein of Klebsiella pneumoniae nitrogenase and show that the molecular conformation of the protein that contains only one molybdenum per alpha(2)beta(2) tetramer is different from that of the protein that has full occupancy i.e. two molybdenums per molecule. This structural finding is consistent with the existence of MoFe protein molecules that contain only one FeMo cofactor site occupied and provides a rationale for the 50% loss of the specific activity of such preparations. A stable inactive transition state complex has been shown to form in the presence of MgADP and AlF(4)(-). Gel filtration chromatography data show that the MoFe protein lacking a full complement of the cofactor forms initially a 1:1 complex before forming a low affinity 1:2 complex. A similar behavior is found for the MoFe protein with both cofactors occupied, but the high affinity 1:2 complex is formed at a lower ratio of Fe protein/MoFe protein. The 1:1 complex, MoFe protein-Fe protein x (ADP x AlF(4)(-))(2), formed with MoFe protein that lacks one of the cofactors, is stable. X-ray scattering studies of this complex have enabled us to obtain its low resolution structure at approximately 20-A resolution, which confirms the gel filtration finding that only one molecule of the Fe protein binds the MoFe protein. By comparison with the low resolution structure of purified MoFe protein that contains only one molybdenum per tetramer, we deduce that the Fe protein interacts with the FeMo cofactor-binding alpha-subunit of the MoFe protein. This observation demonstrates that the conformation of the alpha-subunit or the alpha beta subunit pair that lacks the FeMo cofactor is altered and that the change is recognized by the Fe protein. The structure of the 1:1 complex reveals a similar change in the conformation of the Fe protein as has been observed in the low resolution scattering mask and the high resolution crystallographic study of the 1:2 complex where both cofactors are occupied and with the Fe protein bound to both subunits. This extensive conformational change observed for the Fe protein in the complexes is, however, not observed when MgATP or MgADP binds to the isolated Fe protein. Thus, the large scale conformational change of the Fe protein is associated with the complex formation of the two proteins.

Conformational variability of the Cu site in one subunit of bovine CuZn superoxide dismutase: the importance of mobility in the Glu119-Leu142 loop region for catalytic function.

The structure of the catalytic site in one subunit of bovine CuZn superoxide dismutase is shown to be highly variable. A series of crystal structures at approximately 1.7 A have been determined using data collected from different crystals. Several conformations are observed for the copper site from one of the subunits. These conformations lie between those expected for the Cu(II) and Cu(I) forms of the enzyme and may represent a slow positional rearrangement of the Cu site during the crystallisation process due to the presence of a trace reductant in the mother liquor. These states perhaps indicate some functionally relevant structural steps that ultimately result in the breakage of the imidazolate bridge between the two metal sites. This behaviour is not observed for the second subunit of the dimeric enzyme, which remains in the five-coordinate, distorted square planar geometry in all cases. We suggest that this asymmetric behaviour may be caused by the lack of mobility for the Glu119-Leu142 loop region in the second subunit caused by crystal contacts. This region forms one wall of the active-site cavity, and its mobility has been suggested, via molecular dynamics studies, to be important for the catalytic mechanism.

Implementation of cluster analysis for ab initio phasing using the molecular envelope from solution X-ray scattering.

Solution of the phase problem is central to crystallographic structure determination. The conventional methods of isomorphous replacement (MIR or SIR) and molecular replacement are ineffective in the absence of a suitable isomorphous heavy-atom derivative or knowledge of the structure of a homologous protein. A recent method utilizing the low-resolution molecular shape determined from solution X-ray scattering data has shown to be successful in locating the molecular shape within the crystallographic unit cell in the case of the trimer nitrite reductase (NiR, 105 kDa) [Hao et al. (1999), Acta Cryst. D55, 243-246]. This was achieved by performing a direct real-space search for orientation and translation using knowledge of the orientation of the polar angles of the non-crystallographic axis obtained by performing a self-rotation on crystallographic data. This effectively reduces the potential six-dimensional search to a four-dimensional one (Eulerian angle gamma and three translational parameters). In the case of NiR, the direct four-dimensional search produced a clear solution that was in good agreement with the known structure. The program FSEARCH incorporating this method has been generalized to handle molecules from all space groups and in particular those in possession of non-crystallographic symmetry. However, the method employed was initially unsuccessful when applied to the small dimeric molecule superoxide dismutase (SOD, 32 kDa) owing to the absence of strong reflections at low resolution caused by saturation at the detector. The determined solution deviated greatly from that of the known structure [Hough & Hasnain (1999), J. Mol. Biol. 287, 579-592]. It was found that once these absent reflections were replaced by a series of randomly generated intensity values and cluster analysis was performed on the output, the signal-to-noise ratio was improved and a most probable solution was found. The electron-density map of the stochastically determined solution agrees well with the known structure; the phase error calculated from this map was 67 degrees within 14 A resolution.

Structures of oxidized and reduced azurin II from Alcaligenes xylosoxidans at 1.75 A resolution.

Crystallographic structures of oxidized and reduced forms of azurin II are reported at 1.75 A resolution. Data were collected using one crystal in each case and by translating the crystal after each oscillation range to minimize photoreduction. Very small differences are observed at the Cu site upon reduction and these cannot be determined with confidence at current resolution. A comparison with the three-dimensional EXAFS reveals a good correspondence for all the ligand distances except for Cu-His46, where a larger deviation of approximately 0.12-0.18 A is observed, indicating that this ligand is more tightly restrained in the crystallographic refinement at the current resolution.

3D EXAFS refinement of the Cu site of azurin sheds light on the nature of structural change at the metal centre in an oxidation-reduction process: an integrated approach combining EXAFS and crystallography.

Three-dimensional information is obtained for the Cu site in azurin at very high resolution by combining high-resolution crystallographic structures and EXAFS data for the oxidized and reduced form of the protein. This combined approach has allowed us to define the subtle structural changes (<0.1 A) which take place at the Cu site during a single-electron redox process.

Low-resolution molecular structures of isolated functional units from arthropodan and molluscan hemocyanin.

Synchrotron x-ray scattering measurements were performed on dilute solutions of the purified hemocyanin subunit (Bsin1) from scorpion (Buthus sindicus) and the N-terminal functional unit (Rta) from a marine snail (Rapana thomasiana). The model-independent approach based on spherical harmonics was applied to calculate the molecular envelopes directly from the scattering profiles. Their molecular shapes in solution could be restored at 2-nm resolution. We show that these units represent stable, globular building blocks of the two hemocyanin families and emphasize their conformational differences on a subunit level. Because no crystallographic or electron microscopy data are available for isolated functional units, this study provides for the first time structural information for isolated, monomeric functional subunits from both hemocyanin families. This has been made possible through the use of low protein concentrations (< or = 1 mg/ml). The observed structural differences may offer advantages in building very different overall molecular architectures of hemocyanin by the two phyla.

Role of the axial ligand in type 1 Cu centers studied by point mutations of met148 in rusticyanin.

Type 1 Cu centers in cupredoxins, nitrite reductases, and multi-copper oxidases utilize the same trigonal core ligation to His-Cys-His, with a weak axial ligand generally provided by a Met sulfur. In azurin, an additional axial ligand, a carbonyl oxygen from a Gly, is present. The importance of these axial ligands and in particular the Met has been debated extensively in terms of their role in fine-tuning the redox potential, spectroscopic properties, and rack-induced or entatic state properties of the copper sites. Extensive site-directed mutagenesis of the Met ligand has been carried out in azurin, but the presence of an additional carbonyl oxygen axial ligand has made it difficult to interpret the effects of these substitutions. Here, the axial methionine ligand (Met148) in rusticyanin is replaced with Leu, Gln, Lys, and Glu to examine the effect on the redox potential, acid stability, and copper site geometry. The midpoint redox potential varies from 363 (Met148Lys) to 798 mV (Met148Leu). The acid stability of the oxidized proteins is reduced except for the Met148Gln mutant. The Gln mutant remains blue at all pH values between 2.8 and 8, and has a redox potential of 563 mV at pH 3.2. The optical and rhombic EPR properties of this mutant closely resemble those of stellacyanin, which has the lowest redox potential among single-type 1 copper proteins (185 mV). The Met148Lys mutant exhibits type 2 Cu EPR and optical spectra in this pH range. The Met148Glu mutant exhibits a type 2 Cu EPR spectrum above pH 3 and a mixture of type 1 and type 2 Cu spectra at lower pH. The Met148Leu mutant exhibits the highest redox potential ( approximately 800 mV at pH 3.2) which is similar to the values in fungal laccase and in the type 1 Cu site of ceruloplasmin where this axial ligand is also a Leu.

Is single-wavelength anomalous scattering sufficient for solving phases? A comparison of different methods for a 2.1 A structure solution.

The structure of rusticyanin is the largest unknown structure (M(r) = 16.8 kDa) which has been recently solved by the direct-methods approach using only single-wavelength anomalous scattering (SAS) data from the native protein [Harvey et al. (1998). Acta Cryst. D54, 629-635]. Here, the results of the Sim distribution approach [Hendrickson & Teeter (1981). Nature (London), 290, 107-113] and of the CCP4 procedure MLPHARE [Collaborative Computational Project, Number 4 (1994). Acta Cryst. D50, 760-763] are compared with those from direct methods. Analysis against the final refined model shows that direct methods produced significantly better phases (average phase error 56 degrees ) and therefore significantly better electron-density maps than the Sim distribution and MLPHARE approaches (average phase error was around 63 degrees in both cases).

Iron coordination structures of oxygen sensor FixL characterized by Fe K-edge extended x-ray absorption fine structure and resonance raman spectroscopy.

FixL is a heme-based O(2) sensor protein involved in a two-component system of a symbiotic bacterium. In the present study, the iron coordination structure in the heme domain of Rhizobium meliloti FixLT (RmFixLT, a soluble truncated FixL) was examined using Fe K-edge extended x-ray absorption fine structure (EXAFS) and resonance Raman spectroscopic techniques. In the EXAFS analyses, the interatomic distances and angles of the Fe-ligand bond and the iron displacement from the heme plane were obtained for RmFixLT in the Fe(2+), Fe(2+)O(2), Fe(2+)CO, Fe(3+), Fe(3+)F(-), and Fe(3+)CN(-) states. An apparent correlation was found between the heme-nitrogen (proximal His-194) distance in the heme domain and the phosphorylation activity of the histidine kinase domain. Comparison of the Fe-CO coordination geometry between RmFixLT and RmFixLH (heme domain of RmFixL), based on the EXAFS and Raman results, has suggested that the kinase domain directly or indirectly influences steric interaction between the iron-bound ligand and the heme pocket. Referring to the crystal structure of the heme domain of Bradyrhizobium japonicum FixL (Gong, W., Hao, B., Mansy, S. S., Gonzalez, G., Gilles-Gonzalez, M. A., and Chan, M. K. (1998) Proc. Natl. Acad. Sci. U. S. A. 95, 15177-15182), we discussed details of the iron coordination structure of RmFixLT and RmFixLH in relation to an intramolecular signal transduction mechanism in its O(2) sensing.