Oxidative folding of the cystine knot motif in cyclotide proteins

The cyclotides are a large family of plant proteins that have a cyclic backbone and a knotted arrangement of three conserved disulfide bonds. Despite the apparent complexity of their cystine knot motif it is possible to efficiently fold these proteins, as exemplified by oxidative folding studies on the prototypic cyclotide, kalata B1. This mini-review reports on the current understanding of the folding process in cyclotides. The synthesis and folding of these molecules paves the way for their application as stable molecular templates.

Identification of initiation sites for T4 lysozyme folding using CD and NMR spectroscopy of peptide fragments

Using CD and 2D H-1 NMR spectroscopy, we have identified potential initiation sites for the folding of T4 lysozyme by examining the conformational preferences of peptide fragments corresponding to regions of secondary structure. CD spectropolarimetry showed most peptides were unstructured in water, but adopted partial helical conformations in TFE and SDS solution. This was also consistent with the H-1 NMR data which showed that the peptides were predominantly disordered in water, although in some cases, nascent or small populations of partially folded conformations could be detected. NOE patterns, coupling constants, and deviations from random coil Her chemical shift values complemented the CD data and confirmed that many of the peptides were helical in TFE and SDS micelles. In particular, the peptide corresponding to helix E in the native enzyme formed a well-defined helix in both TFE and SDS, indicating that helix E potentially forms an initiation site for T4 lysozyme folding. The data for the other peptides indicated that helices D, F, G, and H are dependent on tertiary interactions for their folding and/or stability. Overall, the results from this study, and those of our earlier studies, are in agreement with modeling and IID-deuterium exchange experiments, and support an hierarchical model of folding for T4 lysozyme.

Conformational analysis of LYS(11-36), a peptide derived from the beta-sheet region of T4 lysozyme, in TFE and SDS

The solution conformation of a peptide LYS(11-36), which corresponds to the beta-sheet region in T4 lysozyme, has been examined in aqueous solution, TFE, and SDS micelles by CD and H-1 NMR spectroscopy. Secondary structure predictions suggest some beta-sheet and turn character in aqueous solution but predict a helical conformation in a more hydrophobic environment. The predictions were supported by the CD and NMR studies which showed the peptide to be relatively unstructured in aqueous solution, although there was some evidence of a beta-turn conformer which was maintained in 200 mM SDS and, to a lesser extent, in 50% TFE. The peptide was significantly helical in the presence of either 50% TFE or 200 mM SDS. TFE and SDS titrations showed that the peptide could form helical, sheet, or extended structure depending on the TFE or SDS concentration. The studies indicate that peptide environment is the determining factor in secondary structure adopted by LYS(11-36).

Synthesis, NMR studies and conformational analysis of oxazolidine derivatives of the beta-adrenoreceptor antagonists metoprolol, atenolol and timolol

Formaldehyde-derived oxazolidine derivatives 4-7 of the beta-adrenoreceptor antagonists metoprolol 1, atenolol 2 and timolol 3 have been synthesised. Conformational analysis of 1-3 and the oxazolidine derivatives 4-7 has been performed using H-1 NMR spectroscopy and computational methods. The H-1 NMR studies show that for the aryloxypropanolamine beta-adrenoreceptor antagonists there is a predominance of the conformer in which the amine group is approximately antiperiplanar or trans to the aryloxymethylene group. Both H-1 NMR data and theoretical studies indicate that the oxazolidine derivatives 4-7 and the aryloxypropanolamine beta-adrenoreceptor antagonists 1-3 adopt similar conformations around the beta-amino alcohol moiety. Thus, oxazolidine ring formation does not dramatically alter the preferred conformation adopted by the beta-amino alcohol moiety of 1-3. Oxazolidine derivatives of aryloxypropanolamine beta-adrenoreceptor antagonists may therefore be appropriate as prodrugs, or semi-rigid analogues, when greater lipophilicity is required for drug delivery.

Structure determination of the three disulfide bond isomers of alpha-conotoxin GI: A model for the role of disulfide bonds in structural stability

The three possible disulfide bonded isomers of alpha-conotoxin GI have been selectively synthesised and their structures determined by H-1 NMR spectroscopy. alpha-Conotoxin GI derives from the venom of Conus geographus and is a useful neuropharmacological tool as it selectively binds to the nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel involved in nerve signal transmission. The peptide has the sequence ECCNPACGRHYSC-NH2, and the three disulfide bonded isomers are referred to as GI(2-7;3-13), GI(2-13;3-7) and GI(2-3;7-13). The NMR structure for the native isomer GI(2-7;3-13) is of excellent quality, with a backbone pairwise RMSD of 0.16 Angstrom for a family of 35 structures, and comprises primarily a distorted 3(10),, helix between residues 5 to 11. The two non-native isomers exhibit multiple conformers in solution, with the major populated forms being different in structure both from each other and from the native form. Structure-activity relationships for the native GI(2-7;3-13) as well as the role of the disulfide bonds on folding and stability of the three isomers are examined. It is concluded that the disulfide bonds in alpha-conotoxin GI play a crucial part in determining both the structure and stability of the peptide. A trend for increased conformational heterogeneity was observed in the order of GI(2-7;3-13) < GI(2-13;3-7) < GI(2-3;7-13). It was found that the peptide bond joining Cys2 to Cys3 in GI(2-3;7-13) is predominantly trans, rather than cis as theoretically predicted. These structural data are used to interpret the varying nAChR binding of the non-native forms. A model for the binding of native GI(2-7;3-13) to the mammalian nAChR is proposed, with an alpha-subunit binding face made up of Cys2, Asn4, Pro5, Ala6 and Cys7 and a selectivity face, comprised of Arg9 and His10. These two faces orient the molecule between the alpha and delta subunits of the receptor. The structure of the CCNPAC sequence of the native GI(2-7;3-13) is compared to the structure of the identical sequence from the toxic domain of heat-stable enterotoxins, which forms part of the receptor binding region of the enterotoxins, but which has a different disulfide connectivity. (C) 1998 Academic Press Limited.

Crystal structures of reduced and oxidized DsbA: investigation of domain motion and thiolate stabilization

Background: The redox proteins that incorporate a thioredoxin fold have diverse properties and functions. The bacterial protein-folding factor DsbA is the most oxidizing of the thioredoxin family. DsbA catalyzes disulfide-bond formation during the folding of secreted proteins, The extremely oxidizing nature of DsbA has been proposed to result from either domain motion or stabilizing active-site interactions in the reduced form. In the domain motion model, hinge bending between the two domains of DsbA occurs as a result of redox-related conformational changes. Results: We have determined the crystal structures of reduced and oxidized DsbA in the same crystal form and at the same pH (5.6). The crystal structure of a lower pH form of oxidized DsbA has also been determined (pH 5.0). These new crystal structures of DsbA, and the previously determined structure of oxidized DsbA at pH 6.5, provide the foundation for analysis of structural changes that occur upon reduction of the active-site disulfide bond. Conclusions: The structures of reduced and oxidized DsbA reveal that hinge bending motions do occur between the two domains. These motions are independent of redox state, however, and therefore do not contribute to the energetic differences between the two redox states, instead, the observed domain motion is proposed to be a consequence of substrate binding. Furthermore, DsbA's highly oxidizing nature is a result of hydrogen bond, electrostatic and helix-dipole interactions that favour the thiolate over the disulfide at the active site.

Folding, calcium binding, and structural characterization of a concatemer of the first and second ligand-binding modules of the low-density lipoprotein receptor

The ligand-binding domain of the low-density lipoprotein (LDL) receptor is comprised of seven tandemly repeated ligand-binding modules, each being approximately 40 amino acids long and containing six conserved cysteine residues. We have expressed and characterized a concatemer of the first two modules (LB1 and LB2) of the human LDL receptor. Oxidative folding of the recombinant concatemer (rLB(1-2)), in the presence of calcium ions, gave a single dominant isomer with six disulfide bonds. Peptic cleavage of the short Linker region that connects the last cysteine residue of LB1 and the first cysteine residue of LB2 yielded two discrete fragments, thus excluding the presence of intermodule disulfide bonds. The N-terminal module, LB1, reacted with a conformation-specific monoclonal antibody (IgG-C7) made to LB1 in the native LDL receptor. From this, we concluded that the first module was correctly folded, with the same set of disulfide bonds as LB1 of the LDL receptor. The disulfide bond connections of LB2 were identified from mass spectral analysis of fragments formed by digestion of the C-terminal peptic fragment with elastase. These data showed that the disulfide bonds of LB2 connected Cys(I) and Cys(III), Cys(II) and Cys(V), and Cys(IV) and Cys(VI). This pattern is identical to that found for recombinant LB1 and LB2. The concatemer has two high-affinity calcium-binding sites, one per module. An analysis of the secondary chemical shifts of C alpha protons shows that the conformations of LB1 and LB2 in the concatemer are very similar to those of the individual modules, with no evidence for strong interactions between the two modules.

Solution structure of methionine-oxidized amyloid beta-peptide (1-40). Does oxidation affect conformational switching?

The solution structure of A beta(1-40)Met(O), the methionine-oxidized form of amyloid beta-peptide A beta(1-40), has been investigated by CD and NMR spectroscopy. Oxidation of Met35 may have implications in the aetiology of Alzheimer's disease. Circular dichroism experiments showed that whereas A beta(1-40) and A beta(1-40)Met(O) both adopt essentially random coil structures in water (pH 4) at micromolar concentrations, the former aggregates within several days while the latter is stable for at least 7 days under these conditions. This remarkable difference led us to determine the solution structure of A beta(1-40)Met(O) using H-1 NMR spectroscopy. In a water-SDS micelle medium needed to solubilize both peptides at the millimolar concentrations required to measure NMR spectra, chemical shift and NOE data for A beta(1-40)Met(O) strongly suggest the presence of a helical region between residues 16 and 24. This is supported by slow H-D exchange of amide protons in this region and by structure calculations using simulated annealing with the program XPLOR. The remainder of the structure is relatively disordered. Our previously reported NMR data for A beta(1-40) in the same solvent shows that helices are present over residues 15-24 (helix 1) and 28-36 (helix 2), Oxidation of Met35 thus causes a local and selective disruption of helix 2. In addition to this helix-coil rearrangement in aqueous micelles, the CD data show that oxidation inhibits a coil-to-beta-sheet transition in water. These significant structural rearrangements in the C-terminal region of A beta may be important clues to the chemistry and biology of A beta(1-40) and A beta(1-42).

Chemical synthesis and folding pathways of large cyclic polypeptide: Studies of the cystine knot polypeptide kalata B1

Kalata B1 is a member of a new family of polypeptides, isolated from. plants, which have a cystine knot structure embedded within an amide-cyclized backbone. This family of molecules are the largest known cyclic peptides, and thus, the mechanism of synthesis and folding is of great interest. To provide information about both these phenomena, we have synthesized kalata B1 using two distinct strategies. In the first, oxidation of the cysteine residues of a linear precursor peptide to form the correct disulfide bonds results in folding of the three-dimensional structure and preorganization of the termini in close proximity for subsequent cyclization. The second approach involved cyclization prior to oxidation. In the first method, the correctly folded peptide was produced only in the presence of partially hydrophobic solvent conditions. These conditions are presumably required to stabilize the surface-exposed hydrophobic residues. However,; in the synthesis,involving cyclization prior to oxidation, the cyclic reduced peptide folded to a significant degree in the absence of hydrophobic solvents and even more efficiently in the presence of hydrophobic solvents. Cyclization clearly has a major effect on the folding pathway and facilitates formation of the correctly disulfide-bonded form in aqueous solution; In addition to facilitating folding to a compact stable structure cyclization has an important effect on biological activity as assessed by hemolytic activity.

Conformational selection of inhibitors and substrates by proteolytic enzymes: Implications for drug design and polypeptide processing

Inhibitors of proteolytic enzymes (proteases) are emerging as prospective treatments for diseases such as AIDS and viral infections, cancers, inflammatory disorders, and Alzheimer's disease. Generic approaches to the design of protease inhibitors are limited by the unpredictability of interactions between, and structural changes to, inhibitor and protease during binding. A computer analysis of superimposed crystal structures for 266 small molecule inhibitors bound to 48 proteases (16 aspartic, 17 serine, 8 cysteine, and 7 metallo) provides the first conclusive proof that inhibitors, including substrate analogues, commonly bind in an extended beta-strand conformation at the active sites of all these proteases. Representative superimposed structures are shown for (a) multiple inhibitors bound to a protease of each class, (b) single inhibitors each bound to multiple proteases, and (c) conformationally constrained inhibitors bound to proteases. Thus inhibitor/substrate conformation, rather than sequence/composition alone, influences protease recognition, and this has profound implications for inhibitor design. This conclusion is supported by NMR, CD, and binding studies for HIV-1 protease inhibitors/ substrates which, when preorganized in an extended conformation, have significantly higher protease affinity. Recognition is dependent upon conformational equilibria since helical and turn peptide conformations are not processed by proteases. Conformational selection explains the resistance of folded/structured regions of proteins to proteolytic degradation, the susceptibility of denatured proteins to processing, and the higher affinity of conformationally constrained 'extended' inhibitors/substrates for proteases. Other approaches to extended inhibitor conformations should similarly lead to high-affinity binding to a protease.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Multiple conformational epitopes are recognized by natural and induced immunity to the E7 protein of human papilloma virus type 16 in man

The reactivity of sera from patients with cervical cancer with the E7 protein of human papilloma virus type 16 (HPV16) was estimated using a novel non-radioactive immunoprecipitation assay and four established protein-and peptide-based immunoassays. Six of 14 sera from patients with cervical cancer and 1 of 10 sera from healthy laboratory staff showed repeated reactivity with E7 in at least one assay. Four of the 7 reactive sera were consistently reactive in more than one assay, but only one was reactive in all four assays. Following immunization with E7, 2 of 5 patients with cervical cancer had increased E7-specific reactivity, measurable in one or more assays. No single assay was particularly sensitive for E7 reactivity, or predictive of cervical cancer. Mapping of E7 reactivity to specific E7 peptides was unsuccessful, suggesting that natural or induced E7 reactivity in human serum is commonly directed to conformational epitopes of E7, These results suggest that each assay employed with is study measures a different aspect of E7 reactivity, and that various reactivities to E7 may manifest following HPV infection or immunization. This finding is of significance for monitoring of E7 immunotherapy and for serological screening for cervical cancer. Copyright (C) 2000 S.Karger, AG. Basel.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.