Structural basis for binding of Smac/DIABLO to the XIAP BIR3 domain

The inhibitor-of-apoptosis proteins (IAPs) regulate programmed cell death by inhibiting members of the caspase family of enzymes. Recently, a mammalian protein called Smac (also named DIABLO) was identified that binds to the IAPs and promotes caspase activation. Although undefined in the X-ray structure, the amino-terminal residues of Smac are critical for its function. To understand the structural basis for molecular recognition between Smac and the IAPs, we determined the solution structure of the BIR3 domain of X-linked IAP (XIAP) complexed with a functionally active nine-residue peptide derived from the N terminus of Smac. The peptide binds across the third beta-strand of the BIR3 domain in an extended conformation with only the first four residues contacting the protein. The complex is stabilized by four intermolecular hydrogen bonds, an electrostatic interaction involving the N terminus of the peptide, and several hydrophobic interactions. This structural information, along with the binding data from BIR3 and Smac peptide mutants reported here, should aid in the design of small molecules that may be used for the treatment of cancers that overexpress IAPs.

Rationale for Bcl-xL/Bad peptide complex formation from structure, mutagenesis, and biophysical studies

The three-dimensional structure of the anti-apoptotic protein Bcl-xL complexed to a 25-residue peptide from the death promoting region of Bad was determined using NMR spectroscopy. Although the overall structure is similar to Bcl-xL bound to a 16-residue peptide from the Bak protein (Sattler et al., 1997), the Bad peptide forms additional interactions with Bcl-xL. However, based upon site-directed mutagenesis experiments, these additional contacts do not account for the increased affinity of the Bad 25-mer for Bcl-xL compared to the Bad 16-mer. Rather, the increased helix propensity of the Bad 25-mer is primarily responsible for its greater affinity for Bcl-xL. Based on this observation, a pair of 16-residue peptides were designed and synthesized that were predicted to have a high helix propensity while maintaining the interactions important for complexation with Bcl-xL. Both peptides showed an increase in helix propensity compared to the wild-type and exhibited an enhanced affinity for Bcl-xL.

An approach for high-throughput structure determination of proteins by NMR spectroscopy

An approach is described for rapidly determining protein structures by NMR that utilizes proteins containing 13C-methyl labeled Val, Leu, and Ile (delta1) and protonated Phe and Tyr in a deuterated background. Using this strategy, the key NOEs that define the hydrophobic core and overall fold of the protein are easily obtained. NMR data are acquired using cryogenic probe technology which markedly reduces the spectrometer time needed for data acquisition. The approach is demonstrated by determining the overall fold of the antiapoptotic protein, Bcl-xL, from data collected in only 4 days. Refinement of the Bcl-xL structure to a backbone rmsd of 0.95 A was accomplished with data collected in an additional 3 days. A distance analysis of 180 different proteins and structure calculations using simulated data suggests that our method will allow the global folds of a wide variety of proteins to be determined.

NMR structure and mutagenesis of the third Bir domain of the inhibitor of apoptosis protein XIAP

The inhibitor of apoptosis proteins (IAPs) regulate the caspase family of cysteine proteases, which play an important role in the execution of programmed cell death. Human X-linked inhibitor of apoptosis protein (XIAP) is a potent inhibitor of caspases-3, -7, and -9. Here we show that the Bir3 domain is the minimal region of XIAP that is needed for potent caspase-9 inhibition. The three-dimensional structure of the Bir3 domain of XIAP, determined by NMR spectroscopy, resembles a classical zinc finger and consists of five alpha-helices, a three-stranded beta-sheet, and a zinc atom chelated to three cysteines and one histidine. The structure of the Bir3 domain is similar to that of the Bir2 domain of XIAP but differs from the previously determined structure of the Bir3 domain of MIHB. Based on site-directed mutagenesis, we have identified the regions of the Bir3 domain of XIAP that are important for inhibiting caspase-9. Despite the structural similarities of the Bir2 and Bir3 domain of XIAP, a different set of residues were found to be critical for inhibiting the individual caspases. These results suggest that XIAP inhibits caspase-3 and caspase-9 in a different manner.

NMR structure and mutagenesis of the inhibitor-of-apoptosis protein XIAP

The inhibitor-of-apoptosis (IAP) family of proteins, originally identified in baculoviruses, regulate programmed cell death in a variety of organisms. IAPs inhibit specific enzymes (caspases) in the death cascade and contain one to three modules of a common 70-amino-acid motif called the BIR domain. Here we describe the nuclear magnetic resonance structure of a region encompassing the second BIR domain (BIR2) of a human IAP family member, XIAP (also called hILP or MIHA). The structure of the BIR domain consists of a three-stranded antiparallel beta-sheet and four alpha-helices and resembles a classical zinc finger. Unexpectedly, conserved amino acids within the linker region between the BIR1 and BIR2 domains were found to be critical for inhibiting caspase-3. The absence or presence of these residues may explain the differences in caspase inhibition observed for different truncated and full-length IAPs. Our data further indicate that these residues may bind to the active site and that the BIR domain may interact with an adjacent site on the enzyme.

NMR-based screening in drug discovery

1. Introduction 211

2. Screening methods 213

2.1 Chemical shifts 213

2.2 Diffusion 214

2.3 Transverse relaxation 218

2.4 Nuclear Overhauser effects 218

3. Strategies for drug discovery and design 221

3.1 Fragment-based methods 221

3.1.1 Linked-fragment approach 221

3.1.2 Directed combinatorial libraries 222

3.1.3 Modification of high-affinity ligands 223

3.1.4 Solvent mapping techniques 223

3.2 High-throughput NMR-based screening 224

3.3 Enzymatic assays 226

4. Discovery of novel ligands 227

4.1 High-affinity ligands for FKBP 227

4.2 Potent inhibitors of stromelysin 229

4.3 Ligands for the DNA-binding domain of the E2 protein 233

4.4 Discovery of Erm methyltransferase inhibitors 233

4.5 Phosphotyrosine mimetics for SH2 domains 236

5. Conclusions 237

6. References 237

A critical step in the drug discovery process is the identification of high-affinity ligands for macromolecular targets. Traditionally, the identification of such lead compounds has been accomplished through the high-throughout screening (HTS) of corporate compound repositories. Conventional HTS methodology has enjoyed widespread application and success in the pharmaceutical industry and, through recent technological advances in screening (Fernandes, 1998; Oldenburg et al. 1998; Silverman et al. 1998) and combinatorial chemistry (Fauchere et al. 1998; Fecik et al. 1998), these programs will continue to have a prominent role in drug discovery. However, suitable leads cannot always be found using conventional methods. This is not surprising since typical corporate libraries contain fewer than 106 compounds compared with the estimated 1050–1080 universe of compounds (Martin, 1997). In addition, most conventional assays are limited to screening libraries of compounds against proteins with known function, excluding the large number of targets becoming available from genomics research.

Recently, a number of NMR-based screening methods have been employed to identify and design lead ligands for protein targets (see Table 1). These NMR-based strategies can augment ongoing conventional HTS for identifying leads and can be used to aid in lead optimization. All of these techniques take advantage of the fact that upon complex formation between a target molecule and a ligand, significant perturbations can be observed in NMR-sensitive parameters of either the target or the ligand. These perturbations can be used qualitatively to detect ligand binding or quantitatively to assess the strength of the binding interaction. In addition, some of the techniques allow the identification of the ligand binding site or which part of the ligand is responsible for interacting with the target. In this article, the current state of NMR-based screening is reviewed.

Accurate measurement of H(N)-H(alpha) residual dipolar couplings in proteins

A method for accurately measuring H(N)-H(alpha) residual dipolar couplings is described. Using this technique, both the sign and magnitude of the coupling can be determined easily. Residual dipolar coupling between H(N)(i)-H(alpha)(i) and H(N)(i)-H(alpha)(i-1) were measured for the FK506 binding protein complexed to FK506. The experimental values were in excellent agreement with predictions based on an X-ray crystal structure of the protein/ligand complex, suggesting that these residual dipolar couplings will provide accurate structural constraints for the refinement of protein structures determined by NMR.

NMR structure and mutagenesis of the N-terminal Dbl homology domain of the nucleotide exchange factor Trio

Guanine nucleotide exchange factors for the Rho family of GTPases contain a Dbl homology (DH) domain responsible for catalysis and a pleckstrin homology (PH) domain whose function is unknown. Here we describe the solution structure of the N-terminal DH domain of Trio that catalyzes nucleotide exchange for Rac1. The all-alpha-helical protein has a very different structure compared to other exchange factors. Based on site-directed mutagenesis, functionally important residues of the DH domain were identified. They are all highly conserved and reside in close proximity on two a helices. In addition, we have discovered a unique capability of the PH domain to enhance nucleotide exchange in DH domain-containing proteins.

Solution structure of the cytohesin-1 (B2-1) Sec7 domain and its interaction with the GTPase ADP ribosylation factor 1

Cytohesin-1 (B2-1) is a guanine nucleotide exchange factor for human ADP ribosylation factor (Arf) GTPases, which are important for vesicular protein trafficking and coatamer assembly in the cell. Cytohesin-1 also has been reported to promote cellular adhesion via binding to the beta2 integrin cytoplasmic domain. The solution structure of the Sec7 domain of cytohesin-1, which is responsible for both the protein's guanine nucleotide exchange factor function and beta2 integrin binding, was determined by NMR spectroscopy. The structure consists of 10 alpha-helices that form a unique tertiary fold. The binding between the Sec7 domain and a soluble, truncated version of human Arf-1 was investigated by examining 1H-15N and 1H-13C chemical shift changes between the native protein and the Sec7/Arf-1 complex. We show that the binding to Arf-1 occurs through a large surface on the C-terminal subdomain that is composed of both hydrophobic and polar residues. Structure-based mutational analysis of the cytohesin-1 Sec7 domain has been used to identify residues important for binding to Arf and for mediating nucleotide exchange. Investigations into the interaction between the Sec7 domain and the beta2 integrin cytoplasmic domain suggest that the two proteins do not interact in the solution phase.

NMR structure and mutagenesis of the FADD (Mort1) death-effector domain

When activated, membrane-bound receptors for Fas and tumour-necrosis factor initiate programmed cell death by recruiting the death domain of the adaptor protein FADD to the membrane. FADD then activates caspase 8 (also known as FLICE or MACH) through an interaction between the death-effector domains of FADD and caspase 8. This ultimately leads to the apoptotic response. Death-effector domains and homologous protein modules known as caspase-recruitment domains have been found in several proteins and are important regulators of caspase (FLICE) activity and of apoptosis. Here we describe the solution structure of a soluble, biologically active mutant of the FADD death-effector domain. The structure consists of six antiparallel, amphipathic alpha-helices and resembles the overall fold of the death domains of Fas and p75. Despite this structural similarity, mutations that inhibit protein-protein interactions involving the Fas death domain have no effect when introduced into the FADD death-effector domain. Instead, a hydrophobic region of the FADD death-effector domain that is not present in the death domains is vital for binding to FLICE and for apoptotic activity.

Structure of Bcl-xL-Bak peptide complex: recognition between regulators of apoptosis

Heterodimerization between members of the Bcl-2 family of proteins is a key event in the regulation of programmed cell death. The molecular basis for heterodimer formation was investigated by determination of the solution structure of a complex between the survival protein Bcl-xL and the death-promoting region of the Bcl-2-related protein Bak. The structure and binding affinities of mutant Bak peptides indicate that the Bak peptide adopts an amphipathic alpha helix that interacts with Bcl-xL through hydrophobic and electrostatic interactions. Mutations in full-length Bak that disrupt either type of interaction inhibit the ability of Bak to heterodimerize with Bcl-xL.

NMR structure and mutagenesis of the Fas (APO-1/CD95) death domain

Programmed cell death (apoptosis) mediated by the cytokine receptor Fas is critical for the removal of autoreactive T cells, the mechanism of immune privilege, and for maintenance of immune-system homeostasis. Signalling of programmed cell death involves the self-association of a conserved cytoplasmic region of Fas called the death domain and interaction with another death-domain-containing protein, FADD (also known as MORT1). Although death domains are found in several proteins, their three-dimensional structure and the manner in which they interact is unknown. Here we describe the solution structure of the Fas death domain, as determined by NMR spectroscopy. The structure consists of six antiparallel, amphipathic alpha-helices arranged in a novel fold. From the structure and from site-directed mutagenesis, we have identified the region of the death domain involved in self-association and binding to the downstream signalling partner FADD.

Discovering high-affinity ligands for proteins: SAR by NMR

A nuclear magnetic resonance (NMR)-based method is described in which small organic molecules that bind to proximal subsites of a protein are identified, optimized, and linked together to produce high-affinity ligands. The approach is called "SAR by NMR" because structure-activity relationships (SAR) are obtained from NMR. With this technique, compounds with nanomolar affinities for the FK506 binding protein were rapidly discovered by tethering two ligands with micromolar affinities. The method reduces the amount of chemical synthesis and time required for the discovery of high-affinity ligands and appears particularly useful in target-directed drug research.

X-ray and NMR structure of human Bcl-xL, an inhibitor of programmed cell death

THE Bcl-2 family of proteins regulate programmed cell death by an unknown mechanism. Here we describe the crystal and solution structures of a Bcl-2 family member, Bcl-xL (ref. 2). The structures consist of two central, primarily hydrophobic alpha-helices, which are surrounded by amphipathic helices. A 60-residue loop connecting helices alpha1 and alpha2 was found to be flexible and non-essential for anti-apoptotic activity. The three functionally important Bcl-2 homology regions (BH1, BH2 and BH3) are in close spatial proximity and form an elongated hydrophobic cleft that may represent the binding site for other Bcl-2 family members. The arrangement of the alpha-helices in Bcl-xL is reminiscent of the membrane translocation domain of bacterial toxins, in particular diphtheria toxin and the colicins. The structural similarity may provide a clue to the mechanism of action of the Bcl-2 family of proteins.

Structural basis for IL-4 receptor phosphopeptide recognition by the IRS-1 PTB domain

We present the NMR structure of the PTB domain of insulin receptor substrate-1 (IRS-1) complexed to a tyrosine-phosphorylated peptide derived from the IL-4 receptor. Despite the lack of sequence homology and different binding specificity, the overall fold of the protein is similar to that of the Shc PTB domain and closely resembles that of PH domains. However, the PTB domain of IRS-1 is smaller than that of Shc (110 versus 170 residues) and binds to phosphopeptides in a distinct manner. We explain the phosphopeptide binding specificity based on the structure of the complex and results of site-directed mutagenesis experiments.

Solution structure of the DNA-binding domain of a human papillomavirus E2 protein: evidence for flexible DNA-binding regions

The three-dimensional structure of the DNA-binding domain of the E2 protein from human papillomavirus-31 was determined by using multidimensional heteronuclear nuclear magnetic resonance (NMR) spectroscopy. A total of 1429 NMR-derived distance and dihedral angle restraints were obtained for each of the 83-residue subunits of this symmetric dimer. The average root mean square deviations of 20 structures calculated using a distance geometry-simulated annealing protocol are 0.59 and 0.90 angstroms for the backbone and all heavy atoms, respectively, for residues 2-83. The structure of the human virus protein free in solution consists of an eight-stranded beta-barrel and two pairs of alpha-helices. Although the overall fold of the protein is similar to the crystal structure of the bovine papillomavirus-1 E2 protein when complexed to DNA, several small but interesting differences were observed between these two structures at the subunit interface. In addition, a beta-hairpin that contacts DNA in the crystal structure of the protein-DNA complex is disordered in the NMR structures, and steady-state 1H-15N heteronuclear NOE measurements indicate that this region is highly mobile in the absence of DNA. The recognition helix also appears to be flexible, as evidenced by fast amide exchange rates. This phenomenon has also been observed for a number of other DNA-binding proteins and may constitute a common theme in protein/DNA recognition.

Structure and ligand recognition of the phosphotyrosine binding domain of Shc

The nuclear magnetic resonance structure of the phosphotyrosine binding (PTB) domain of Shc complexed to a phosphopeptide reveals an alternative means of recognizing tyrosine-phosphorylated proteins. Unlike in SH2 domains, the phosphopeptide forms an antiparallel beta-strand with a beta-sheet of the protein, interacts with a hydrophobic pocket through the (pY-5) residue, and adopts a beta-turn. The PTB domain is structurally similar to pleckstrin homology domains (a beta-sandwich capped by an alpha-helix) and binds to acidic phospholipids, suggesting a possible role in membrane localization.

Solution structure of the Shc SH2 domain complexed with a tyrosine-phosphorylated peptide from the T-cell receptor

She is a widely expressed adapter protein that plays an important role in signaling via a variety of cell surface receptors and has been implicated in coupling the stimulation of growth factor, cytokine, and antigen receptors to the Ras signaling pathway. She interacts with several tyrosine-phosphorylated receptors through its C-terminal SH2 domain, and one of the mechanisms of T-cell receptor-mediated Ras activation involves the interaction of the Shc SH2 domain with the tyrosine-phosphorylated zeta chain of the T-cell receptor. Here we describe a high-resolution NMR structure of the Shc SH2 domain complexed to a phosphopeptide (GHDGLpYQGLSTATK) corresponding to a portion of the zeta chain of the T-cell receptor. Although the overall architecture of the protein is similar to other SH2 domains, distinct structural differences were observed in the smaller beta-sheet, BG loop, (pY + 3) phosphopeptide-binding site, and relative position of the bound phosphopeptide.

A genetic algorithm-based protocol for docking ensembles of small ligands using experimental restraints

A genetic algorithm (GA) based method for docking ensembles of small, flexible ligands to receptor proteins using NMR-derived constraints is described. In this method, three translations and rotations of the ligand and the dihedral angles of the ligand are represented by binary strings and evolve under the genetic operators of cross-over, mutation, migration and selection. The fitness function for the selection process includes distance and dihedral restraints and a repulsive van der Waals term. The GA was applied to a three-atom model system as well as to the streptavidin-biotin complex using simulated intermolecular distance restraints. In both systems, the GA was able to obtain low-energy conformations when only a single binding site was simulated. Calculations were also performed using distance restraints from two distinct binding sites. In these simulations, the GA was able to obtain low-energy conformations corresponding to ligand molecules in each of the two sites. The inclusion of additional ligands in the ensemble did not result in an energetic benefit, confirming that only two ligand conformations were necessary to fulfill the distance restraints. This method allows for a direct investigation of the minimum number of ligand orientations necessary to fulfill experimental distance restraints, and simultaneously yields detailed structural information about each site.

Backbone dynamics of a two-domain protein: 15N relaxation studies of the amino-terminal fragment of urokinase-type plasminogen activator

The amino-terminal fragment (ATF) of urokinase-type plasminogen activator (u-PA) is a two-domain protein which consists of a kringle and a growth factor domain (GFD). The dynamics of uniformly 15N-labeled ATF was examined by measuring the longitudinal (T1) and transverse (T2) 15N relaxation times and heteronuclear NOEs. The data were interpreted in terms of the model-independent spectral density function. The GFD was found to exhibit a high degree of anisotropy, whereas the kringle domain of ATF undergoes isotropic reorientation. This difference in anisotropy is best explained by the two domains moving independently such as differently shaped beads on a string. With the exception of the N- and C-terminal regions of the protein, the most flexible region of ATF was the seven-residue omega loop (N22-I28) of the GFD which has been implicated in the binding of u-PA to its receptor. The amides of the linker region between the domains displayed high values of the order parameter, indicating restricted motion on the picosecond time scale. This is in contrast to the flexible linker of calmodulin [Barbato et al. (1992) Biochemistry 31, 5269-5278], which displayed low values of S2 and unrestricted motion in the linker region.

Solution structure of the ets domain of Fli-1 when bound to DNA

Members of the ets family of transcription factors share a conserved DNA-binding domain, the ets domain. By using multidimensional NMR, we have determined the structure of the ets domain of human Fli-1 in the DNA-bound form. It consists of three alpha-helices and a four-stranded beta-sheet, similar to structures of the class of helix-turn-helix DNA binding proteins first found in the catabolite activator protein of Escherichia coli. NMR and mutagenesis experiments suggest that in comparison to structurally related proteins, the ets domain uses a new variation of the helix-turn-helix motif for binding to DNA.

Solution structure of a low molecular weight protein tyrosine phosphatase

Protein tyrosine phosphatases (PTPs) are important enzymes involved in signal transduction, cell cycle regulation, and the control of differentiation. Despite the importance of this class of enzymes in the control of critical cell processes, very little structural information is available for this family of proteins. In this paper, we present the first solution structure of a protein tyrosine phosphatase. This protein is a low molecular weight cytosolic PTP that was initially isolated from bovine heart. The structure that was determined from 1747 NMR-derived restraints consists of a central four-stranded parallel beta-sheet surrounded by four alpha-helices and a short 3(10) helix. The phosphate binding site, identified by chemical shift changes upon the addition of the competitive inhibitors phosphate and vanadate, is in a loop region connecting the C-terminal end of the first beta-strand with the first alpha-helix. Residues in the second, fourth, and fifth alpha-helices and in some of the loop regions connecting the elements of regular secondary structure also contribute to the binding site. The structure determined here is consistent with previous mutagenesis and chemical modification studies conducted on this protein.

Solution structure of a pleckstrin-homology domain

Pleckstrin, the major protein kinase C substrate of platelets, contains domains of about 100 amino acids at the amino and carboxy termini that have been found in a number of proteins, including serine/threonine kinases, GTPase-activating proteins, phospholipases and cytoskeletal proteins. These conserved sequences, termed pleckstrin-homology (PH) domains, are thought to be involved in signal transduction. But the details of the function and binding partners of the PH domains have not been characterized. Here we report the solution structure of the N-terminal pleckstrin-homology domain of pleckstrin determined using heteronuclear three-dimensional nuclear magnetic resonance spectroscopy. The structure consists of an up-and-down beta-barrel of seven antiparallel beta-strands and a C-terminal amphiphilic alpha-helix that caps one end of the barrel. The overall topology of the domain is similar to that of the retinol-binding protein family of structures.

Solution structure of the amino-terminal fragment of urokinase-type plasminogen activator

The amino-terminal fragment (ATF) of urokinase-type plasminogen activator is a two domain protein which consists of a growth factor and a kringle domain. The 1H, 13C, and 15N chemical shifts of this protein have been assigned using heteronuclear two- and three-dimensional NMR experiments on selective and uniformly 15N- and 15N/13C-labeled protein isolated from mammalian cells that overexpress the protein. The chemical shift assignments were used to interpret the NOE data which resulted in a total of 1299 NOE restraints. The NOE restraints were used along with 27 phi angle restraints and 21 hydrogen-bonding restraints to produce 15 low energy structures. The individual domains in the structures are highly converged, but the two domains are structurally independent. The root mean square deviations (rmsd) between residues 11-46 in the growth factor domain and the mean atomic coordinates were 0.99 +/- 0.2 for backbone heavy atoms and 1.65 +/- 0.2 for all non-hydrogen atoms. For residues 55-130 in the kringle domain, the rmsd was 0.84 +/- 0.2 for backbone heavy atoms and 1.42 +/- 0.2 for all non-hydrogen atoms. The overall structures of the individual domains are very similar to the structures of homologous proteins. However, important structural differences between the growth factor and other homologous proteins were observed in the region which has been implicated in binding the urokinase receptor which may explain, in part, why other growth factors show no appreciable affinity for the urokinase receptor.

A computer-based protocol for semiautomated assignments and 3D structure determination of proteins

A strategy is presented for the semiautomated assignment and 3D structure determination of proteins from heteronuclear multidimensional Nuclear Magnetic Resonance (NMR) data. This approach involves the computer-based assignment of the NMR signals, identification of distance restraints from nuclear Overhauser effects, and generation of 3D structures by using the NMR-derived restraints. The protocol is described in detail and illustrated on a resonance assignment and structure determination of the FK506 binding protein (FKBP, 107 amino acids) complexed to the immunosuppressant, ascomycin. The 3D structures produced from this automated protocol attained backbone and heavy atom rmsd of 1.17 and 1.69 A, respectively. Although more highly resolved structures of the complex have been obtained by standard interpretation of NMR data (Meadows et al. (1993) Biochemistry, 32, 754-765), the structures generated with this automated protocol required minimal manual intervention during the spectral assignment and 3D structure calculations stages. Thus, the protocol may yield an approximate order of magnitude reduction in the time required for the generation of 3D structures of proteins from NMR data.

Heteronuclear 3D NMR studies of water bound to an FK506 binding protein/immunosuppressant complex

From a series of 15N-resolved 3D ROESY-HMQC and 13C-resolved 3D NOESY-HMQC spectra of the FK506 binding protein (FKBP)/ascomycin complex in H2O, the locations of three tightly bound water molecules were identified. These waters are all buried within the interior of the complex and form an integral part of its structure via a network of hydrogen bonds. Water molecules in identical locations exhibiting a similar hydrogen bonding pattern were also observed in the X-ray crystal structures of FKBP/FK506 [Van Duyne, G. D., Standaert, R. F., Karplus, P. A., Schreiber, S. L., & Clardy, J. (1991) Science 252, 839-842] and FKBP/rapamycin [Van Duyne, G. D., Standaert, R. F., Schreiber, S. L., & Clardy, J. (1991) J. Am. Chem. Soc. 113, 7433-7434]. However, none of the surface waters observed in the X-ray structures were detected in the NMR experiments due to their fast exchange with bulk water. In order to examine the effects of the three internal water molecules on NMR structure determinations of the FKBP/ascomycin complex, two sets of NMR structures were calculated either with or without the waters. By including the three internal waters in the structure calculations, a decrease in the root mean square deviation and improved angular order parameters was observed for FKBP residues in the vicinity of the water molecules. In addition, subtle conformational differences were observed between NMR structures generated either with or without the waters.(ABSTRACT TRUNCATED AT 250 WORDS)

Three-dimensional structure of the FK506 binding protein/ascomycin complex in solution by heteronuclear three- and four-dimensional NMR

A high-resolution three-dimensional solution structure of the FKBP/ascomycin complex has been determined using heteronuclear multidimensional nuclear magnetic resonance spectroscopy (NMR) and a distance geometry/simulated annealing protocol. A total of 43 structures of the complex, including 3 tightly bound water molecules, were obtained using 1958 experimental restraints consisting of 1724 nuclear Overhauser effect (NOE) derived distances, 66 chi 1 and 46 phi angular restraints, and 122 hydrogen bond restraints. The root mean square (rms) deviations between the 43 FKBP/ascomycin solution structures and the mean atomic coordinates were 0.43 +/- 0.08 A for the backbone heavy atoms and 0.80 +/- 0.08 A for all non-hydrogen atoms. Angular order parameters for the family of 43 conformations were calculated to determine dihedral convergence. Order parameters for phi, psi, and chi 1 angles exhibited mean values of 0.98, 0.97, and 0.95, respectively, while the mean of the chi 2 order parameter was 0.63. Comparisons were made between the FKBP/ascomycin complex and two NMR-derived solution structures of unbound FKBP and the X-ray crystal structure of an FKBP/FK506 complex. Differences were observed between the FKBP/ascomycin complex and uncomplexed FKBP for residues 33-45 and 78-92. In contrast, the NMR-derived solution structure of the FKBP/ascomycin complex and the X-ray crystal structure of the FKBP/FK506 complex were very similar. Differences between the two complexes were mainly observed in the conformations of some highly solvent exposed side chains.

Effect of Elicitation and Changes in Extracellular pH on the Cytoplasmic and Vacuolar pH of Suspension-Cultured Soybean Cells

We have employed both (31)P nuclear magnetic resonance spectroscopy and two intracellular fluorescent pH indicator dyes to monitor the pH of the vacuole and cytoplasm of suspension-cultured soybean cells (Glycine max Merr cv Kent). For the (31)P nuclear magnetic resonance studies, a flow cell was constructed that allowed perfusion of the cells in oxygenated growth medium throughout the experiment. When the perfusion medium was transiently adjusted to a pH higher than that of the ambient growth medium, a rapid elevation of vacuolar pH was observed followed by a slow (approximately 30 minute) return to near resting pH. In contrast, the concurrent pH changes in the cytoplasm were usually fourfold smaller. These data indicate that extracellular pH changes are rapidly communicated to the vacuole in soybean cells without significantly perturbing cytoplasmic pH. When elicitors were dissolved in a medium of altered pH and introduced into the cell suspension, the pH of the vacuole, as above, quickly reflected the pH of the added elicitor solution. In contrast, when the pH of either a polygalacturonic acid or Verticillium dahliae elicitor preparation was adjusted to the same pH as the ambient medium, no significant change in either vacuolar or cytoplasmic pH was observed during the 35 minute experiment. These results were confirmed in experiments with pH-sensitive fluorescent dyes. We conclude that suspension-cultured soybean cells do not respond to elicitation by significantly changing the pH of their vacuolar or cytoplasmic compartments.

How accurately can oligonucleotide structures be determined from the hybrid relaxation rate matrix/NOESY distance restrained molecular dynamics approach?

The accuracy and precision of structures derived from a combined hybrid relaxation rate matrix/NOESY distance restrained molecular dynamics methodology were examined with simulations that included typical experimental errors. NOESY data were simulated for a DNA dodecamer duplex, d-(CGCGAATTCGCG)2, with added volume error of approximately 20% and low-level thermal noise. Distances derived from a hybrid relaxation matrix analysis of the NOE data were used as constraints in molecular dynamics driven structural refinements of several initial model geometries. The final structures were compared against results obtained from the traditional isolated two-spin approximation treatment of these NOESY volumes and also against refined structures that employed error-free data. Results show that the structures derived from the relaxation rate matrix analysis of the NOESY data are more accurate than those derived from a simple two-spin approximation analysis and it is possible to achieve refinement to the level of simulated experimental error. Results may be significantly improved with the use of either more accurately measured NOESY volumes or additional matrix-derived constraints. Many of the helical parameters and backbone torsional angles may be accurately reproduced by the hybrid matrix methodology.

Two-dimensional NMR and structure determination of salmon calcitonin in methanol

The structure of the 32-residue peptide salmon calcitonin (sCT) in 90% MeOH-10% H2O has been investigated by two-dimensional NMR techniques and molecular modeling. Sequential assignments for nearly all of the 32 spin systems have been obtained, and results indicate that the heptaresidue loop formed by the disulfide bond between Cys-1 and Cys-7 is followed by an alpha-helical segment from Val-8 through Tyr-22. A region of conformational heterogeneity is observed for residues 20-25, resulting from the slow isomerism of the cis and trans forms of Pro-23. The C-terminal segment is found to exist in an extended conformation.

NMR structural refinement of a tandem G.A mismatched decamer d(CCAAGATTGG)2 via the hybrid matrix procedure

A complete relaxation matrix approach employing a matrix eigenvalue/eigenvector solution to the Bloch equations is used to evaluate the NMR solution structure of a tandemly positioned G.A double mismatch decamer oligodeoxyribonucleotide duplex, d(CCAAGATTGG)2. An iterative refinement method using a hybrid relaxation matrix combined with restrained molecular dynamics calculations is shown to provide structures having good agreement with the experimentally derived structures. Distances incorporated into the MD simulations have been calculated from the relaxation rate matrix evaluated from a hybrid NOESY volume matrix whose elements are obtained from the merging of experimental and calculated NOESY intensities. Starting from both A- and B-DNA and mismatch syn and anti models, it is possible to calculate structures that are in good atomic RMS agreement with each other (less than 1.6 A RMS) but differ from the reported crystal structure (greater than 3.6 A). Importantly, the hybrid matrix derived structures are in excellent agreement with the experimental solution conformation as determined by comparison of the 200-ms simulated and experimental NOESY spectra, while the crystallographic data provide spectra that are grossly different.

NMR structural refinement of an extrahelical adenosine tridecamer d(CGCAGAATTCGCG)2 via a hybrid relaxation matrix procedure

Until very recently interproton distances from NOESY experiments have been derived solely from the two-spin approximation method. Unfortunately, even at short mixing times, there is a significant error in many of these distances. A complete relaxation matrix approach employing a matrix eigenvalue/eigenvector solution to the Bloch equations avoids the approximation of the two-spin method. We have calculated the structure of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex, d(CGCAGAATTCGCG)2, by an iterative refinement approach using a hybrid relaxation matrix method combined with restrained molecular dynamics calculations. Distances from the 2D NOESY spectra have been calculated from the relaxation rate matrix which has been evaluated from a hybrid NOESY volume matrix comprising elements from the experiment and those calculated from an initial structure. The hybrid matrix derived distances have then been used in a restrained molecular dynamics procedure to obtain a new structure that better approximates the NOESY spectra. The resulting partially refined structure is then used to calculate an improved theoretical NOESY volume matrix which is once again merged with the experimental matrix until refinement is complete. Although the crystal structure of the tridecamer clearly shows the extrahelical adenosine looped out way from the duplex, the NOESY distance restrained hybrid matrix/molecular dynamics structural refinement establishes that the extrahelical adenosine stacks into the duplex.