Rapid simulation of protein motion: merging flexibility, rigidity and normal mode analyses.

Protein function frequently involves conformational changes with large amplitude on timescales which are difficult and computationally expensive to access using molecular dynamics. In this paper, we report on the combination of three computationally inexpensive simulation methods--normal mode analysis using the elastic network model, rigidity analysis using the pebble game algorithm, and geometric simulation of protein motion--to explore conformational change along normal mode eigenvectors. Using a combination of ElNemo and First/Froda software, large-amplitude motions in proteins with hundreds or thousands of residues can be rapidly explored within minutes using desktop computing resources. We apply the method to a representative set of six proteins covering a range of sizes and structural characteristics and show that the method identifies specific types of motion in each case and determines their amplitude limits.

Inhibition of HIV-1 protease: the rigidity perspective.

MOTIVATION: HIV-1 protease is a key drug target due to its role in the life cycle of the HIV-1 virus. Rigidity analysis using the software First is a computationally inexpensive method for inferring functional information from protein crystal structures. We evaluate the rigidity of 206 high-resolution (2 Å or better) X-ray crystal structures of HIV-1 protease and compare the effects of different inhibitors binding to the enzyme. RESULTS: Inhibitor binding has little effect on the overall rigidity of the protein homodimer, including the rigidity of the active site. The principal effect of inhibitor binding on rigidity is to constrain the flexibility of the ß-hairpin flaps, which move to allow access to the active site of the enzyme. We show that commercially available antiviral drugs which target HIV-1 protease can be divided into two classes, those which significantly affect flap rigidity and those which do not. The non-peptidic inhibitor tipranavir is distinctive in its consistently strong effect on flap rigidity. CONTACT: jack.heal@warwick.ac.uk; r.roemer@warwick.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The refolding and reassembly of Escherichia coli heat-labile enterotoxin B-subunit: analysis of reassembly-competent and reassembly-incompetent unfolded states.

The B-subunit pentamer of Escherichia coli heat-labile enterotoxin (EtxB) is an exceptionally stable protein maintaining its quaternary structure over the pH value range 2.0-11.0. Up to 80% yields of reassembled pentamer can be obtained in vitro from material disassembled for very short incubation periods in KCl-HCl, pH 1.0. However, when the incubation period in acid is extended, the reassembly yield decreases to no more than 20% (Ruddock et al. (1996) J. Biol. Chem. 271 19118-19123). Here we demonstrate that the ion species present in the disassembly conditions strongly influence the reassembly competence of EtxB showing that 60% reassembly yields can be achieved, even after prolonged incubations, by the use of a phosphate buffer for acid disassembly. Using this system, we have fully characterized the disassembly and reassembly behavior of EtxB by electrophoretic, immunochemical, and spectroscopic techniques and compared it with that previously observed. Depending on the denaturation system used, the acid-denatured monomer is either in a predominantly reassembly-competent state (H(3)PO(4) system) or in a predominantly reassembly-incompetent conformation (KCl-HCl system). Interconversion between these two conformations in the denatured state is possible by the addition of salts to the denatured protein. The results are consistent with the previous hypothesis that the conversion between reassembly-competent and -incompetent states corresponds to a cis/trans isomerization of a peptide bond, presumably that to Pro93.

The disassembly and reassembly of mutants of Escherichia coli heat-labile enterotoxin: replacement of proline 93 does not abolish the reassembly-competent and reassembly-incompetent states.

The carrier moiety of heat-labile enterotoxin of Escherichia coli (EtxB) is formed by the noncovalent association of identical monomeric subunits, which assemble, in vivo and in vitro, into exceptionally stable pentameric complexes. In vitro, acid disassembly followed by neutralization results in reassembly yields of between 20% and 60% depending on the identity of the salts present during the acid denaturation process. Loss of reassembly competence has been attributed to isomerization of the native cis-proline residue at position 93. To characterize this phenomenon further, two mutants of EtxB at proline 93 (P93G and P93A) were generated and purified. The proline variants reveal only minor differences in their biophysical and biochemical properties relative to wild-type protein, but major changes were observed in the kinetics of pentamer disassembly and reassembly. Additionally, a loss of assembly competence was observed following longer term acid treatment, which was even more marked than that of the wild-type protein. We present evidence that the loss of assembly competence of these mutants is best explained by a cis/trans peptidyl isomerization of the unfolded mutant subunits in acid conditions; this limited reassembly competence and the biophysical properties of the native P93 mutant pentamers imply the retention of the native cis conformation in the nonproline peptide bond between residues 92 and 93 in the mutated proteins.

The influence of His94 and Pro149 in modulating the activity of V. cholerae DsbA.

DsbA is the primary catalyst of disulfide bond formation in the periplasm of gram-negative bacteria. Numerous theoretical and experimental studies have been undertaken to determine the molecular mechanisms by which DsbA acts as a potent oxidant, whereas the homologous cytoplasmic protein, thioredoxin, acts as a reductant. Many of these studies have focused on the nature of the two residues that lie between the active-site cysteines. Although these are clearly important, they are not solely responsible for the differences in activity between these thiol-disulfide oxidoreductases. Q97 in the helical domain of E. coli DsbA has been implicated in influencing the redox potential of E. coli DsbA. In V. cholerae DsbA, the analogous residue is H94. In this study, the effect of H94 on the oxidase activity of DsbA is examined, along with the role of the conserved cis-proline residue P149. The DsbA mutant H94L shows a nearly fourfold increase in activity over the wild-type enzyme. To our knowledge, this is the first time an increase in the normal activity of a thiol-disulfide oxidoreductase has been reported. Potential reasons for this increase in activity are discussed.

The pancreas-specific protein disulphide-isomerase PDIp interacts with a hydroxyaryl group in ligands.

Using a cross-linking approach, we have recently demonstrated that radiolabelled model peptides or misfolded proteins specifically interact in vitro with two members of the protein disulphide- isomerase family, namely PDI and PDIp, in a crude extract from sheep pancreas microsomes. In addition, we have shown that tyrosine and tryptophan residues within a peptide are the recognition motifs for the binding to PDIp. Here we examine non-peptide ligands and present evidence that a hydroxyaryl group is a structural motif for the binding to PDIp; simple constructs containing this group and certain xenobiotics and phytoestrogens, which contain an unmodified hydroxyaryl group, can all efficiently inhibit peptide binding to PDIp. To our knowledge this is the first time that the recognition motif of a molecular chaperone or folding catalyst has been specified as a simple chemical structure.

Metabolic control of recombinant protein N-glycan processing in NS0 and CHO cells.

Chinese hamster ovary and murine myeloma NS0 cells are currently favored host cell types for the production of therapeutic recombinant proteins. In this study, we compared N-glycan processing in GS-NS0 and GS-CHO cells producing the same model recombinant glycoprotein, tissue inhibitor of metalloproteinases 1. By manipulation of intracellular nucleotide-sugar content, we examined the feasibility of implementing metabolic control strategies aimed at reducing the occurrence of murine-specific glycan motifs on NS0-derived recombinant proteins, such as Galalpha1,3Galbeta1,4GlcNAc. Although both CHO and NS0-derived oligosaccharides were predominantly of the standard complex type with variable sialylation, 30% of N-glycan antennae associated with NS0-derived TIMP-1 terminated in alpha1,3-linked galactose residues. Furthermore, NS0 cells conferred a greater proportion of terminal N-glycolylneuraminic (sialic) acid residues as compared with the N-acetylneuraminic acid variant. Inclusion of the nucleotide-sugar precursors, glucosamine (10 mM, plus 2 mM uridine) and N-acetylmannosamine (20 mM), in culture media were shown to significantly increase the intracellular pools of UDP-N-acetylhexosamine and CMP-sialic acid, respectively, in both NS0 and CHO cells. The elevated UDP-N-acetylhexosamine content induced by the glucosamine/uridine treatment was associated with an increase in the antennarity of N-glycans associated with TIMP-1 produced in CHO cells but not N-glycans associated with TIMP-1 from NS0 cells. In addition, elevated UDP-N-acetylhexosamine content was associated with a slight decrease in sialylation in both cell lines. The elevated CMP-sialic acid content induced by N-acetylmannosamine had no effect on the overall level of sialylation of TIMP-1 produced by both CHO and NS0 cells, although the ratio of N-glycolylneuraminic acid:N-acetylneuraminic acid associated with NS0-derived TIMP-1 changed from 1:1 to 1:2. These data suggest that manipulation of nucleotide-sugar metabolism can promote changes in N-glycan processing that are either conserved between NS0 and CHO cells or specific to either NS0 cells or CHO cells.

Domains b' and a' of protein disulfide isomerase fulfill the minimum requirement for function as a subunit of prolyl 4-hydroxylase. The N-terminal domains a and b enhances this function and can be substituted in part by those of ERp57.

Protein disulfide isomerase (PDI) is a modular polypeptide consisting of four domains, a, b, b', and a', plus an acidic C-terminal extension, c. PDI carries out multiple functions, acting as the beta subunit in the animal prolyl 4-hydroxylases and in the microsomal triglyceride transfer protein and independently acting as a protein folding catalyst. We report here that the minimum sequence requirement for the assembly of an active prolyl 4-hydroxylase alpha(2)beta(2) tetramer in insect cell coexpression experiments is fulfilled by the PDI domain construct b'a' but that the sequential addition of the b and a domains greatly increases the level of enzyme activity obtained. In the assembly of active prolyl 4-hydroxylase tetramers, the a and b domains of PDI, but not b' and a', can in part be substituted by the corresponding domains of ERp57, a PDI isoform that functions naturally in association with the lectins calnexin and calreticulin. The a' domain of PDI could not be substituted by the PDI a domain, suggesting that both b' and a' domains contain regions critical for prolyl 4-hydroxylase assembly. All PDI domain constructs and PDI/ERp57 hybrids that contain the b' domain can bind the 14-amino acid peptide Delta-somatostatin, as measured by cross-linking; however, binding of the misfolded protein "scrambled" RNase required the addition of domains ab or a' of PDI. The human prolyl 4-hydroxylase alpha subunit has at least two isoforms, alpha(I) and alpha(II), which form with the PDI polypeptide the (alpha(I))(2)beta(2) and (alpha(II))(2)beta(2) tetramers. We report here that all the PDI domain constructs and PDI/ERp57 hybrid polypeptides tested were more effectively associated with the alpha(II) subunit than the alpha(I) subunit.

Mutations that destabilize the a' domain of human protein-disulfide isomerase indirectly affect peptide binding.

Protein-disulfide isomerase (PDI) is a catalyst of folding of disulfide-bonded proteins and also a multifunctional polypeptide that acts as the beta-subunit in the prolyl 4-hydroxylase alpha(2)beta(2)-tetramer (P4H) and the microsomal triglyceride transfer protein alphabeta-dimer. The principal peptide-binding site of PDI is located in the b' domain, but all domains contribute to the binding of misfolded proteins. Mutations in the C-terminal part of the a' domain have significant effects on the assembly of the P4H tetramer and other functions of PDI. In this study we have addressed the question of whether these mutations in the C-terminal part of the a' domain, which affect P4H assembly, also affect peptide binding to PDI. We observed a strong correlation between P4H assembly competence and peptide binding; mutants of PDI that failed to form a functional P4H tetramer were also inactive in peptide binding. However, there was also a correlation between inactivity in these assays and indicators of conformational disruption, such as protease sensitivity. Peptide binding activity could be restored in inactive, protease-sensitive mutants by selective proteolytic removal of the mutated a' domain. Hence we propose that structural changes in the a' domain indirectly affect peptide binding to the b' domain.

Specificity in substrate binding by protein folding catalysts: tyrosine and tryptophan residues are the recognition motifs for the binding of peptides to the pancreas-specific protein disulfide isomerase PDIp.

Using a cross-linking approach, we recently demonstrated that radiolabeled peptides or misfolded proteins specifically interact in vitro with two luminal proteins in crude extracts from pancreas microsomes. The proteins were the folding catalysts protein disulfide isomerase (PDI) and PDIp, a glycosylated, PDI-related protein, expressed exclusively in the pancreas. In this study, we explore the specificity of these proteins in binding peptides and related ligands and show that tyrosine and tryptophan residues in peptides are the recognition motifs for their binding by PDIp. This peptide-binding specificity may reflect the selectivity of PDIp in binding regions of unfolded polypeptide during catalysis of protein folding.

Early intermediates in the PDI-assisted folding of ribonuclease A.

The oxidative refolding of ribonuclease A has been investigated in several experimental conditions using a variety of redox systems. All these studies agree that the formation of disulfide bonds during the process occurs through a nonrandom mechanism with a preferential coupling of certain cysteine residues. We have previously demonstrated that in the presence of glutathione the refolding process occurs through the reiteration of two sequential reactions: a mixed disulfide with glutathione is produced first which evolves to form an intramolecular S-S bond. In the same experimental conditions, protein disulfide isomerase (PDI) was shown to catalyze formation and reduction of mixed disulfides with glutathione as well as formation of intramolecular S-S bonds. This paper reports the structural characterization of the one-disulfide intermediate population during the oxidative refolding of Ribonuclease A under the presence of PDI and glutathione with the aim of defining the role of the enzyme at the early stages of the reaction. The one-disulfide intermediate population occurring at the early stages of both the uncatalyzed and the PDI-catalyzed refolding was purified and structurally characterized by proteolytic digestion followed by MALDI-MS and LC/ESIMS analyses. In the uncatalyzed refolding, a total of 12 disulfide bonds out of the 28 theoretical possible cysteine couplings was observed, confirming a nonrandom distribution of native and nonnative disulfide bonds. Under the presence of PDI, only two additional nonnative disulfides were detected. Semiquantitative LC/ESIMS analysis of the distribution of the S-S bridged peptides showed that the most abundant species were equally populated in both the uncatalyzed and the catalyzed process. This paper shows the first structural characterization of the one-disulfide intermediate population formed transiently during the refolding of ribonuclease A in quasi-physiological conditions that mimic those present in the ER lumen. At the early stages of the process, three of the four native disulfides are detected, whereas the Cys26-Cys84 pairing is absent. Most of the nonnative disulfide bonds identified are formed by nearest-neighboring cysteines. The presence of PDI does not significantly alter the distribution of S-S bonds, suggesting that the ensemble of single-disulfide species is formed under thermodynamic control.

The effect of matrix metalloproteinase complex formation on the conformational mobility of tissue inhibitor of metalloproteinases-2 (TIMP-2).

The backbone mobility of the N-terminal domain of tissue inhibitor of metalloproteinases-2 (N-TIMP-2) was determined both for the free protein and when bound to the catalytic domain of matrix metalloproteinase-3 (N-MMP-3). Regions of the protein with internal motion were identified by comparison of the T(1) and T(2) relaxation times and (1)H-(15)N nuclear Overhauser effect values for the backbone amide (15)N signals for each residue in the sequence. This analysis revealed rapid internal motion on the picosecond to nanosecond time scale for several regions of free N-TIMP-2, including the extended beta-hairpin between beta-strands A and B, which forms part of the MMP binding site. Evidence of relatively slow motion indicative of exchange between two or more local conformations on a microsecond to millisecond time scale was also found in the free protein, including two other regions of the MMP binding site (the CD and EF loops). On formation of a tight N-TIMP-2. N-MMP-3 complex, the rapid internal motion of the AB beta-hairpin was largely abolished, a change consistent with tight binding of this region to the MMP-3 catalytic domain. The extended AB beta-hairpin is not a feature of all members of the TIMP family; therefore, the binding of this highly mobile region to a site distant from the catalytic cleft of the MMPs suggests a key role in TIMP-2 binding specificity.

High resolution structure of the N-terminal domain of tissue inhibitor of metalloproteinases-2 and characterization of its interaction site with matrix metalloproteinase-3.

The high resolution structure of the N-terminal domain of tissue inhibitor of metalloproteinases-2 (N-TIMP-2) in solution has been determined using multidimensional heteronuclear NMR spectroscopy, with the structural calculations based on an extensive set of constraints, including 3132 nuclear Overhauser effect-based distance constraints, 56 hydrogen bond constraints, and 220 torsion angle constraints (an average of 26.9 constraints/residue). The core of the protein consists of a five-stranded beta-barrel that is homologous to the beta-barrel found in the oligosaccharide/oligonucleotide binding protein fold. The binding site for the catalytic domain of matrix metalloproteinases-3 (N-MMP-3) on N-TIMP-2 has been mapped by determining the changes in chemical shifts on complex formation for signals from the protein backbone (15N, 13C, and 1H). This approach identified a discrete N-MMP-3 binding site on N-TIMP-2 composed of the N terminus of the protein and the loops between beta-strands AB, CD, and EF. The beta-hairpin formed from strands A and B in N-TIMP-2 is significantly longer than the equivalent structure in TIMP-1, allowing it to make more extensive binding interactions with the MMP catalytic domain. A detailed comparison of the N-TIMP-2 structure with that of TIMP-1 bound to N-MMP-3 (Gomis-Ruth, F.-X., Maskos, K., Betz, M., Bergner, A., Huber, R., Suzuki, K., Yoshida, N., Nagase, H. , Brew, K., Bourne, G. P., Bartunik, H. & Bode, W. (1997) Nature 389, 77-80) revealed that the core beta-barrels are very similar in topology but that the loop connecting beta-strands CD (P67-C72) would need to undergo a large conformational change for TIMP-2 to bind in a similar manner to TIMP-1.

A pancreas-specific glycosylated protein disulphide-isomerase binds to misfolded proteins and peptides with an interaction inhibited by oestrogens.

Using a cross-linking approach, we have demonstrated that radiolabeled model peptides or misfolded proteins specifically interact in vitro with two different luminal proteins in a crude extract from sheep pancreas microsomes. One of the proteins was identified as protein disulphide-isomerase (PDI), the other one was a related protein (PDIp). We have shown that PDIp was expressed exclusively in the pancreas. Interspecies conservation of PDIp was confirmed and, unlike other members of the PDI family, PDIp from various sources was found to be a glycoprotein. PDIp interacted with peptides and also a misfolded protein, but not with native proteins, suggesting that it might act as a molecular chaperone. The initial binding process was independent of the presence of Cys residues in the probed peptides. Certain oestrogens strongly inhibited the interaction between peptides and PDIp, with 17beta-oestradiol being the most potent inhibitor.

Protein folding: a missing redox link in the endoplasmic reticulum.

Native disulphide-bond formation during protein folding in the endoplasmic reticulum requires oxidative machinery, the components and mechanism of which are not yet fully understood. Two recent papers have identified a novel protein component that appears to play a key role in this important redox pathway.

Novel disulfide oxidoreductase in search of a function.

The high-resolution structure of an oxidoreductase from the archaeon Pyrococcus furiosus shows that it comprises two adjacent domains each with the thioredoxin (trx/grx) fold. Its functional properties are not yet fully defined but it may be related to the multi-domain eukaryotic protein disulfide-isomerases.

Experimental and theoretical analyses of the domain architecture of mammalian protein disulphide-isomerase.

The high resolution structure of full-length protein disulphide-isomerase (PDI) has not been determined, but the polypeptide is generally assumed to comprise a series of consecutive domains. Models of its domain organisation have been proposed on the basis of various sequence-based criteria and, more recently, from structural studies on recombinant fragments corresponding to putative domains. We here describe direct studies of the domain architecture of full-length mammalian PDI based on limited proteolysis of the native enzyme. The results are consistent with an emerging model based on the existence of 4 consecutive domains each with the thioredoxin fold. The model was further tested by expressing recombinant fragments corresponding to alternative domain models and to truncated domains; the observed properties of these purified fragments supported the 4-domain model. A multiple alignment of many PDI-like sequences was generated to test whether domain boundaries could be predicted from any features of the alignment, such as sequence variability or hydrophilicity; neither of these parameters reliably predicted the domain boundaries determined by experiment.

The b' domain provides the principal peptide-binding site of protein disulfide isomerase but all domains contribute to binding of misfolded proteins.

Protein disulfide isomerase (PDI) is a very efficient catalyst of folding of many disulfide-bonded proteins. A great deal is known about the catalytic functions of PDI, while little is known about its substrate binding. We recently demonstrated by cross-linking that PDI binds peptides and misfolded proteins, with high affinity but broad specificity. To characterize the substrate-binding site of PDI, we investigated the interactions of various recombinant fragments of human PDI, expressed in Escherichia coli, with different radiolabelled model peptides. We observed that the b' domain of human PDI is essential and sufficient for the binding of small peptides. In the case of larger peptides, specifically a 28 amino acid fragment derived from bovine pancreatic trypsin inhibitor, or misfolded proteins, the b' domain is essential but not sufficient for efficient binding, indicating that contributions from additional domains are required. Hence we propose that the different domains of PDI all contribute to the binding site, with the b' domain forming the essential core.

Mapping the binding site for matrix metalloproteinase on the N-terminal domain of the tissue inhibitor of metalloproteinases-2 by NMR chemical shift perturbation.

Changes in the NMR chemical shift of backbone amide nuclei (1H and 15N) have been used to map the matrix metalloproteinase (MMP) binding site on the N-terminal domain of the tissue inhibitor of metalloproteinase-2 (N-TIMP-2). Amide chemical shift changes were measured on formation of a stable complex with the catalytic domain of stromelysin-1 (N-MMP-3). Residues with significantly shifted amide signals mapped specifically to a broad site covering one face of the molecule. This site (the MMP binding site) consists primarily of residues 1-11, 27-41, 68-73, 87-90, and 97-104. The site overlaps with the OB-fold binding site seen in other proteins that share the same five-stranded beta-barrel topology. Sequence conservation data and recent site-directed mutagenesis studies are discussed in relation to the MMP binding site identified in this work.

Interactions between protein disulphide isomerase and peptides.

There is growing evidence that protein disulphide isomerase (PDI) has a common chaperone function in the endoplasmic reticulum. To characterise this function, we investigated the interaction of purified PDI with radiolabelled model peptides, somatostatin and mastoparan, by cross-linking. The interaction between the peptides and PDI was specific, for it showed saturation and was abolished by denaturation of PDI. The interaction between a hydrophobic peptide without cysteine residues was much more sensitive to Triton X-100 than the interaction between PDI and a more hydrophilic peptide with or without cysteine residues. We therefore propose that hydrophobic interactions between protein disulphide isomerase and peptides play an important role in the binding process. The interaction between PDI and the bound peptide therefore is enhanced by the formation of mixed disulphide bonds.