Determinants of backbone dynamics in native BPTI: cooperative influence of the 14-38 disulfide and the Tyr35 side-chain

Conformational dynamics in peptide toxins: Implications for receptor interactions and molecular design

Peptide toxins are potent and often exquisitely selective probes of the structure and function of ion channels and receptors, and are therefore of significant interest to the pharmaceutical and biotech industries as both pharmacological tools and therapeutic leads. The three-dimensional structures of peptide toxins are essential as a basis for understanding their structure-activity relationships and their binding to target receptors, as well as in guiding the design of analogues with modified potency and/or selectivity for key targets. NMR spectroscopy has played a key role in elucidating the structures of peptide toxins and probing their structure-function relationships. In this article, we highlight the additional important contribution of NMR to characterising the dynamics of peptide toxins. We also compare the information available from NMR measurements with that afforded by molecular dynamics simulations. We describe several examples of the importance of dynamics measurements over a range of timescales for understanding the structure-function relationships of peptide toxins and their receptor engagement. Peptide toxins that inhibit the voltage-gated potassium channel K_V1.3 with pM affinities display different degrees of conformational flexibility, even though they contain multiple disulfide bonds, and this flexibility can affect the relative orientation of residues that have been shown to be critical for channel binding. Information on the dynamic properties of peptide toxins is important in the design of analogues or mimetics where receptor-bound structures are not available.

BPTI folding revisited: switching a disulfide into methylene thioacetal reveals a previously hidden path

The folding mechanism of the model protein bovine pancreatic trypsin inhibitor was revisited. By switching the solvent exposed disulfide bond with methylene thioacetal we uncovered a hidden pathway in its folding mechanism. In addition, this moiety enhanced protein stability while fully maintaining the protein structure and biological function.

Bovine pancreatic trypsin inhibitor (BPTI) is a 58-residue protein that is stabilized by three disulfide bonds at positions 5–55, 14–38 and 30–51. Widely studied for about 50 years, BPTI represents a folding model for many disulfide-rich proteins. In the study described below, we replaced the solvent exposed 14–38 disulfide bond with a methylene thioacetal bridge in an attempt to arrest the folding pathway of the protein at its two well-known intermediates, N′ and N*. The modified protein was expected to be unable to undergo the rate-determining step in the widely accepted BPTI folding mechanism: the opening of the 14–38 disulfide bond followed by rearrangements that leads to the native state, N. Surprisingly, instead of halting BPTI folding at N′ and N*, we uncovered a hidden pathway involving a direct reaction between the N* intermediate and the oxidizing reagent glutathione (GSSG) to form the disulfide-mixed intermediate N*–SG, which spontaneously folds into N. On the other hand, N′ was unable to fold into N. In addition, we found that the methylene thioacetal bridge enhances BPTI stability while fully maintaining its structure and biological function. These findings suggest a general strategy for enhancing protein stability without compromising on function or structure, suggesting potential applications for future therapeutic protein production.

The influence of disulfide bonds on the mechanical stability of proteins is context dependent

Disulfide bonds play a crucial role in proteins, modulating their stability and constraining their conformational dynamics. A particularly important case is that of proteins that need to withstand forces arising from their normal biological function and that are often disulfide bonded. However, the influence of disulfides on the overall mechanical stability of proteins is poorly understood. Here, we used single-molecule force spectroscopy (smFS) to study the role of disulfide bonds in different mechanical proteins in terms of their unfolding forces. For this purpose, we chose the pilus protein FimG from Gram-negative bacteria and a disulfide-bonded variant of the I91 human cardiac titin polyprotein. Our results show that disulfide bonds can alter the mechanical stability of proteins in different ways depending on the properties of the system. Specifically, disulfide-bonded FimG undergoes a 30% increase in its mechanical stability compared with its reduced counterpart, whereas the unfolding force of I91 domains experiences a decrease of 15% relative to the WT form. Using a coarse-grained simulation model, we rationalized that the increase in mechanical stability of FimG is due to a shift in the mechanical unfolding pathway. The simple topology-based explanation suggests a neutral effect in the case of titin. In summary, our results indicate that disulfide bonds in proteins act in a context-dependent manner rather than simply as mechanical lockers, underscoring the importance of considering disulfide bonds both computationally and experimentally when studying the mechanical properties of proteins.

FF12MC: A revised AMBER forcefield and new protein simulation protocol

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened CH bonds, (ii) removal of torsions involving a nonperipheral sp(3) atom, and (iii) reduced 1-4 interaction scaling factors of torsions ϕ and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric-isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left- and right-handed configurations for C14-C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp-cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general-purpose AMBER forcefield FF14SB locks the C14-C38 bond to the right-handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric-isothermal MD simulations 2-4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490-1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Functional and structural roles of the Cys14-Cys38 disulfide of bovine pancreatic trypsin inhibitor

The disulfide bond between Cys14 and Cys38 of bovine pancreatic trypsin inhibitor lies on the surface of the inhibitor and forms part of the protease-binding region. The functional properties of three variants lacking this disulfide, with one or both of the Cys residues replaced with Ser, were examined, and X-ray crystal structures of the complexes with bovine trypsin were determined and refined to the 1.58-A resolution limit. The crystal structure of the complex formed with the mutant with both Cys residues replaced was nearly identical with that of the complex containing the wild-type protein, with the Ser oxygen atoms positioned to replace the disulfide bond with a hydrogen bond. The two structures of the complexes with single replacements displayed small local perturbations with alternate conformations of the Ser side chains. Despite the absence of the disulfide bond, the crystallographic temperature factors show no evidence of increased flexibility in the complexes with the mutant inhibitors. All three of the variants were cleaved by trypsin more rapidly than the wild-type inhibitor, by as much as 10,000-fold, indicating that the covalent constraint normally imposed by the disulfide contributes to the remarkable resistance to hydrolysis displayed by the wild-type protein. The rates of hydrolysis display an unusual dependence on pH over the range of 3.5-8.0, decreasing at the more alkaline values, as compared with the increased hydrolysis rates for normal substrates under these conditions. These observations can be accounted for by a model for inhibition in which an acyl-enzyme intermediate forms at a significant rate but is rapidly converted back to the enzyme-inhibitor complex by nucleophilic attack by the newly created amino group. The model suggests that a lack of flexibility in the acyl-enzyme intermediate, rather than the enzyme-inhibitor complex, may be a key factor in the ability of bovine pancreatic trypsin inhibitor and similar inhibitors to resist hydrolysis.

Three-dimensional solution structure and conformational plasticity of the N-terminal scavenger receptor cysteine-rich domain of human CD5

The lymphocyte receptor CD5 influences cell activation by modifying the strength of the intracellular response initiated by antigen engagement. Regulation through CD5 involves the interaction of one or more of its three scavenger receptor cysteine-rich domains present in the extracellular region. Here, we present the 3D solution structure of a non-glycosylated double mutant of the N-terminal domain of human CD5 expressed in Escherichia coli (eCD5d1m), which has enhanced solubility compared to the non-glycosylated wild-type (eCD5d1). In common with a glycosylated form expressed in Pichia pastoris, the [(15)N,(1)H]-correlation spectra of both eCD5d1 and eCD5d1m exhibit non-uniform temperature-dependent signal intensities, indicating extensive conformational fluctuations on the micro-millisecond timescale. Although approximately one half of the signals expected for the domain are absent at 298 K, essentially complete resonance assignments and a solution structure could be obtained at 318 K. Because of the sparse nature of the experimental restraint data and the potentially important contribution of conformational exchange to the nuclear Overhauser effect peak intensity, we applied inferential structure determination to calculate the eCD5d1m structure. The inferential structure determination ensemble has similar features to that obtained by traditional simulated annealing methods, but displays superior definition and structural quality. The eCD5d1m structure is similar to other members of the scavenger receptor cysteine-rich superfamily, but the position of the lone alpha helix differs due to interactions with the unique N-terminal region of the domain. The availability of an experimentally tractable form of CD5d1, together with its 3D structure, provides new tools for further investigation of its function within intact CD5.

Disulfide bond influence on protein structural dynamics probed with 2D-IR vibrational echo spectroscopy

Intramolecular disulfide bonds are understood to play a role in regulating protein stability and activity. Because disulfide bonds covalently link different components of a protein, they influence protein structure. However, the effects of disulfide bonds on fast (subpicosecond to approximately 100 ps) protein equilibrium structural fluctuations have not been characterized experimentally. Here, ultrafast 2D-IR vibrational echo spectroscopy is used to examine the constraints an intramolecular disulfide bond places on the structural fluctuations of the protein neuroglobin (Ngb). Ngb is a globin family protein found in vertebrate brains that binds oxygen reversibly. Like myoglobin (Mb), Ngb has the classical globin fold and key residues around the heme are conserved. Furthermore, the heme-ligated CO vibrational spectra of Mb (Mb-CO) and Ngb (Ngb-CO) are virtually identical. However, in contrast to Mb, human Ngb has an intramolecular disulfide bond that affects its oxygen affinity and protein stability. By using 2D-IR vibrational echo spectroscopy, we investigated the equilibrium protein dynamics of Ngb-CO by observing the CO spectral diffusion (time dependence of the 2D-IR line shapes) with and without the disulfide bond. Despite the similarity of the linear FTIR spectra of Ngb-CO with and without the disulfide bond, 2D-IR measurements reveal that the equilibrium sampling of different protein configurations is accelerated by disruption of the disulfide bond. The observations indicate that the intramolecular disulfide bond in Ngb acts as an inhibitor of fast protein dynamics even though eliminating it does not produce significant conformational change in the protein's structure.

Rigidification of a flexible protease inhibitor variant upon binding to trypsin

The Tyr35-->Gly replacement in bovine pancreatic trypsin inhibitor (BPTI) has previously been shown to dramatically enhance the flexibility of the trypsin-binding region of the free inhibitor and to destabilize the interaction with the protease by about 3 kcal/mol. The effects of this replacement on the enzyme-inhibitor interaction were further studied here by X-ray crystallography and isothermal titration calorimetry (ITC). The co-crystal structure of Y35G BPTI bound to trypsin was determined using 1.65 A resolution X-ray diffraction data collected from cryopreserved crystals, and a new structure of the complex with wild-type BPTI under the same conditions was determined using 1.62 A data. These structures reveal that, in contrast to the free protein, Y35G BPTI adopts a conformation nearly identical with that of the wild-type protein, with a water-filled cavity in place of the missing Tyr side-chain. The crystallographic temperature factors for the two complexes indicate that the mutant inhibitor is nearly as rigid as the wild-type protein when bound to trypsin. Calorimetric measurements show that the change in enthalpy upon dissociation of the complex is 2.5 kcal/mol less favorable for the complex containing Y35G BPTI than for the complex with the wild-type inhibitor. Thus, the destabilization of the complex resulting from the Y35G replacement is due to a more favorable change in entropy upon dissociation. The heat capacity changes for dissociation of the mutant and wild-type complexes were very similar, suggesting that the entropic effects probably do not arise from solvation effects, but are more likely due to an increase in protein conformational entropy upon dissociation of the mutant inhibitor. These results define the biophysical role of a highly conserved core residue located outside of a protein-binding interface, demonstrating that Tyr35 has little impact on the trypsin-bound BPTI structure and acts primarily to define the structure of the free protein so as to maximize binding affinity.

Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: calculation of structural and dynamic properties

Long dynamics simulations were carried out on the B1 immunoglobulin-binding domain of streptococcal protein G (ProtG) and bovine pancreatic trypsin inhibitor (BPTI) using atomistic descriptions of the proteins and a continuum representation of solvent effects. To mimic frictional and random collision effects, Langevin dynamics (LD) were used. The main goal of the calculations was to explore the stability of tens-of-nanosecond trajectories as generated by this molecular mechanics approximation and to analyze in detail structural and dynamical properties. Conformational fluctuations, order parameters, cross correlation matrices, residue solvent accessibilities, pKa values of titratable groups, and hydrogen-bonding (HB) patterns were calculated from all of the trajectories and compared with available experimental data. The simulations comprised over 40 ns per trajectory for ProtG and over 30 ns per trajectory for BPTI. For comparison, explicit water molecular dynamics simulations (EW/MD) of 3 ns and 4 ns, respectively, were also carried out. Two continuum simulations were performed on each protein using the CHARMM program, one with the all-atom PAR22 representation of the protein force field (here referred to as PAR22/LD simulations) and the other with the modifications introduced by the recently developed CMAP potential (CMAP/LD simulations). The explicit solvent simulations were performed with PAR22 only. Solvent effects are described by a continuum model based on screened Coulomb potentials (SCP) reported earlier, i.e., the SCP-based implicit solvent model (SCP-ISM). For ProtG, both the PAR22/LD and the CMAP/LD 40-ns trajectories were stable, yielding C(alpha) root mean square deviations (RMSD) of about 1.0 and 0.8 A respectively along the entire simulation time, compared to 0.8 A for the EW/MD simulation. For BPTI, only the CMAP/LD trajectory was stable for the entire 30-ns simulation, with a C(alpha) RMSD of approximately 1.4 A, while the PAR22/LD trajectory became unstable early in the simulation, reaching a C(alpha) RMSD of about 2.7 A and remaining at this value until the end of the simulation; the C(alpha) RMSD of the EW/MD simulation was about 1.5 A. The source of the instabilities of the BPTI trajectories in the PAR22/LD simulations was explored by an analysis of the backbone torsion angles. To further validate the findings from this analysis of BPTI, a 35-ns SCP-ISM simulation of Ubiquitin (Ubq) was carried out. For this protein, the CMAP/LD simulation was stable for the entire simulation time (C(alpha) RMSD of approximately 1.0 A), while the PAR22/LD trajectory showed a trend similar to that in BPTI, reaching a C(alpha) RMSD of approximately 1.5 A at 7 ns. All the calculated properties were found to be in agreement with the corresponding experimental values, although local deviations were also observed. HB patterns were also well reproduced by all the continuum solvent simulations with the exception of solvent-exposed side chain-side chain (sc-sc) HB in ProtG, where several of the HB interactions observed in the crystal structure and in the EW/MD simulation were lost. The overall analysis reported in this work suggests that the combination of an atomistic representation of a protein with a CMAP/CHARMM force field and a continuum representation of solvent effects such as the SCP-ISM provides a good description of structural and dynamic properties obtained from long computer simulations. Although the SCP-ISM simulations (CMAP/LD) reported here were shown to be stable and the properties well reproduced, further refinement is needed to attain a level of accuracy suitable for more challenging biological applications, particularly the study of protein-protein interactions.

Alteration of the disulfide-coupled folding pathway of BPTI by circular permutation

The kinetics of disulfide-coupled folding and unfolding of four circularly permuted forms of bovine pancreatic trypsin inhibitor (BPTI) were studied and compared with previously published results for both wild-type BPTI and a cyclized form. Each of the permuted proteins was found to be less stable than either the wild-type or circular proteins, by 3-8 kcal/mole. These stability differences were used to estimate effective concentrations of the chain termini in the native proteins, which were 1 mM for the wild-type protein and 2.5 to 4000 M for the permuted forms. The circular permutations increased the rates of unfolding and caused a variety of effects on the kinetics of refolding. For two of the proteins, the rates of a direct disulfide-formation pathway were dramatically increased, making this process as fast or faster than the competing disulfide rearrangement mechanism that predominates in the folding of the wild-type protein. These two permutations break the covalent connectivity among the beta-strands of the native protein, and removal of these constraints appears to facilitate direct formation and reduction of nearby disulfides that are buried in the folded structure. The effects on folding kinetics and mechanism do not appear to be correlated with relative contact order, a measure of overall topological complexity. These observations are consistent with the results of other recent experimental and computational studies suggesting that circular permutation may generally influence folding mechanisms by favoring or disfavoring specific interactions that promote alternative pathways, rather than through effects on the overall topology of the native protein.

Accelerated Exchange of a Buried Water Molecule in Selectively Disulfide-Reduced Bovine Pancreatic Trypsin Inhibitor

Using magnetic relaxation dispersion (MRD), we have previously shown that the four internal water molecules in bovine pancreatic trypsin inhibitor (BPTI) exchange with bulk water on time scales between 10(-8) and 10(-4) s at room temperature. Because this exchange is controlled by the protein structure, internal water molecules can be used to probe rare conformational fluctuations. Here, we report (2)H and (17)O MRD data at three temperatures for wild-type BPTI and two BPTI variants where the 14-38 disulfide bond has been cleaved by a double Cys --> Ser mutation or by disulfide reduction and carboxamidomethylation. The MRD data show that the internal water molecules are conserved on disulfide cleavage. However, the exchange rate of the water molecule buried near the disulfide bond is enhanced by 2-4 orders of magnitude. The relation of water exchange to other dynamic processes in BPTI is discussed.

Identification of a Residue Critical for Maintaining the Functional Conformation of BPTI

The effects of amino acid replacements on the backbone dynamics of bovine pancreatic trypsin inhibitor (BPTI) were examined using 15N NMR relaxation experiments. Previous studies have shown that backbone amide groups within the trypsin-binding region of the wild-type protein undergo conformational exchange processes on the micros time scale, and that replacement of Tyr35 with Gly greatly increases the number of backbone atoms involved in such motions. In order to determine whether these mutational effects are specific to the replacement of this residue with Gly, six additional replacements were examined in the present study. In two of these, Tyr35 was replaced with either Ala or Leu, and the other four were single replacements of Tyr23, Phe33, Asn43 or Asn44, all of which are highly buried in the native structure and conserved in homologous proteins. The Y35A and Y35L mutants displayed dynamic properties very similar to those of the Y35G mutant, with the backbone segments including residues 10-19 and 32-44 undergoing motions revealed by enhanced 15N transverse relaxation rates. On the other hand, the Y23L, N43G and N44A substitutions caused almost no detectable changes in backbone dynamics, on either the ns-ps or ms-micros time scales, even though each of these replacements significantly destabilizes the native conformation. Replacement of Phe33 with Leu caused intermediate effects, with several residues that have previously been implicated in motions in the wild-type protein displaying enhanced transverse relaxation rates. These results demonstrate that destabilizing amino acid replacements can be accommodated in a native protein with dramatically different effects on conformational dynamics and that Tyr35 plays a particularly important role in defining the conformation of the trypsin-binding site of BPTI.

Solution structure and backbone dynamics of the TGFbeta type II receptor extracellular domain

Isoforms of transforming growth factor beta (TGFbeta) are 25 kDa homodimeric polypeptides that signal by binding and bringing together two related, functionally distinct cell surface receptors designated as TbetaR1 and TbetaR2. Here, we report the solution structure of the 13.8 kDa extracellular domain of human TbetaR2 (ecTbetaR2) as calculated from N(N)-H(N), C(alpha)-H(alpha), and C(alpha)-C(O) residual dipolar coupling restraints in conjunction with NOE distance, dihedral angle, and scalar coupling restraints. Comparison of the free ecTbetaR2 solution structure with the TGFbeta3-bound ecTbetaR2 crystal structure reveals backbone conformations that superimpose with RMSDs of 1.0 A over the regions of regular secondary structure and 1.4 A overall. The differences in structure fall mainly in loop regions that are either poorly defined by the available NMR data or are involved in crystal contacts. The noted similarities between the NMR structure of the free form and the crystal structure of the TGFbeta-bound form are also consistent with the close correspondence, 0.16 A RMSD for regions of secondary structure and 0.51 A RMSD overall, for the crystal structure of free ecTbetaR2 as compared to the crystal structure of TGFbeta3-bound ecTbetaR2. Despite the apparent similarities between the free and the bound forms, there appears to be small but significant differences in structure involving the interfacial contact region of the receptor. Measurements of backbone (15)N relaxation times and interpretation of these by the model-free formalism with axial diffusional anisotropy further reveal significant ms to micros time scale motions centered about two of the conserved disulfide bonds and in several residues that comprise the TGFbeta binding surface. Together, these observations indicate that binding likely occurs through a mechanism with a small component of induced fit character, whereby flexibility within the receptor facilitates the transition to the TGFbeta-bound state.

Native-like interactions favored in the unfolded bovine pancreatic trypsin inhibitor have different roles in folding

Folding kinetics of a series of bovine pancreatic trypsin inhibitor (BPTI) variants with similar stabilities and structures have been measured. All are strongly destabilized relative to WT. In Y21A, F22A, Y23A, G37A, and F45A, the three native disulfide bonds are retained. In RM(14-38), Cys14 and Cys38 thiols are methylated while C30-C51 and C5-C55 disulfides remain intact. At pH 2 and 20 degrees C, relaxation rate constants of the major kinetic phase range from approximately 10 ms to 0.71 s in the absence of denaturant. All mutants except G37A exhibit standard two-state behavior. Y21A, F22A, and Y23A fold much more slowly than other mutants. The experiments were designed to test the hypothesis that native-like structure detected in the unfolded BPTI is important in folding. Two native-like contacts are implied by NOEs in reduced and unfolded BPTI, between residues Tyr23 and Ala25, and between Gly37 NH and the Tyr35 ring. The results support an earlier hypothesis that formation of the central beta-hairpin, monitored by a local native interaction between Tyr23 and Ala25, is crucial to initiation of BPTI folding. The second native-like contact is important, not in folding initiation, but in preventing a kinetic trap later in the process. Evidence for this comes from mutant G37A, which behaves very differently from the others in displaying a phenomenon called rollover. G37A is, to our knowledge, the first reported case in which a single-site replacement causes rollover, while the wild type and all other known mutants of the same protein show typical two-state chevron plots. The best explanation is that the G37A mutation introduces a kinetic trap of the type described by Chan and Dill [(1998) Proteins 30, 2-33]. In native BPTI, there is an unusual polar interaction between the ring of Tyr35 and the backbone NH of Gly37. Our results suggest that the NH-aromatic interaction between residues 37 and 35 is important throughout folding in stabilizing native-like loop conformations and in preventing the flexible loops from being trapped in nonfunctional conformations during later stages of folding.

A highly destabilizing mutation, G37A, of the bovine pancreatic trypsin inhibitor retains the average native conformation but greatly increases local flexibility

A point mutation, G37A, on the surface of bovine pancreatic trypsin inhibitor (BPTI) destabilizes the protein by approximately 5 kcal/mol, which is very high for addition of one methyl group. In wild-type (WT) BPTI, Gly 37 HN is in an unusual NH-aromatic-NH network of interactions with the ring of Tyr 35 and the side chain HN of Asn 44. G37A was designed to disrupt this interaction, since the phi and psi backbone angles of G37 are not favorable for an amino acid containing a beta-carbon. Investigations of the structure and dynamics by NMR methods show that G37A retains the average WT structure. The NH-aromatic-NH interactions remain intact, as indicated by NOEs and the large upfield ring current shift (approximately 4 ppm) of A37 HN. The NMR structure, confirmed by molecular modeling calculations, requires phi and psi backbone angles that are highly destabilizing when alanine is in position 37. Although the average structure is essentially unchanged, the dynamics are altered dramatically. Many residues in the region of the mutation have increased flexibility, as probed by aromatic ring flip rates and native state hydrogen exchange. We conclude that a large fraction of the destabilization arises from maintaining A37 in a high-energy conformation. This suggests that disruption of the NH-aromatic-NH network is energetically very costly, and may involve other cooperatively linked interactions. The results illustrate the importance of the Gly-Gly sequence at positions 36 and 37 and the 37 HN-35 aromatic interaction to the stability, folding, and dynamics of the BPTI.

Mutational analysis of hydrogen bonding residues in the BPTI folding pathway

Nine BPTI variants with replacements that remove one or more hydrogen bonds from the native protein were constructed, and the folding pathways of these proteins were determined by isolating and identifying the disulfide-bonded intermediates that accumulated during unfolding and refolding. The forward and reverse rate constants for the individual steps in the folding pathways for each protein were measured, providing a detailed description of the energetic effects of the substitutions. The native forms of eight of the nine variants were measurably destabilized, by 1-7 kcal/mol (1 cal=4.184 J), with an average effect of 1.6 kcal/mol per hydrogen bond removed. The folding pathways for the variants were found to be similar to that previously described for the wild-type protein, with the kinetically preferred mechanism involving intramolecular rearrangements of intermediates with two disulfide bonds. Some of the substitutions, however, significantly destabilized the major intermediates and broadened the distribution of species with one or two disulfide bonds, thus identifying residues that play important roles in stabilizing the normal intermediates and defining specificity in the folding process. The kinetic data also suggest that one residue, Asn43, may play a distinctive role in defining the BPTI folding mechanism. Replacement of this residue with either Gly or Ala appeared to stabilize the major transition states for folding and unfolding. In the native protein, the side-chain of Asn43 participates directly in the hydrogen bonding pattern of the central beta-sheet, and the kinetic behavior of the Asn43 variants suggests that the major energy barriers in folding and unfolding may be due in part to the steric constraints imposed by this structural element, together with those imposed by the chemical transition states for thiol-disulfide exchange.

NMR-detected order in core residues of denatured bovine pancreatic trypsin inhibitor

The NMR characteristics of [14-38]Abu, a synthetic variant of BPTI that is partially folded in aqueous buffer near neutral pH, support a model of early folding events which begin with stabilization of the nativelike, slow exchange core [Barbar, E., Hare, M., Daragan, V., Barany, G., and Woodward, C. (1998) Biochemistry 37, 7822-7833 (1)]. In partially folded [14-38]Abu, urea denaturation profiles for representative amide protons show that global unfolding is non-two-state and that core residues require a higher concentration of urea to unfold. Dynamic properties of pH-denatured [14-38]Abu and fully reduced and unfolded BPTI analogue were determined from heteronuclear NMR relaxation measurements at similar solution conditions. Differences at various sites in the polypeptide chain were evaluated from spectral density functions determined from T1, T2, and steady-state heteronuclear NOE data. Although denatured [14-38]Abu contains no persistent secondary structure, its most ordered residues are those that, in native BPTI, fold into the slow exchange core. The fully reduced analogue is significantly more mobile and shows less heterogeneous dynamics, but at 1 degree C, restricted motion is observed for residues in the central segments of the polypeptide chain. These observations indicate that there is a developing core or cores even in highly unfolded species. Apparently the effect of 14-38 disulfide on unfolded

Refolding of Taiwan cobra neurotoxin: intramolecular cross-link affects its refolding reaction

In order to explore the effect of intramolecular cross-linking in the folding reaction of cobrotoxin from Naja naja atra (Taiwan cobra) venom, the toxin molecule was modified with glutaraldehyde (GA). The monomeric GA-modified cobrotoxin (mGA-cobrotoxin) was separated from the dimeric and trimeric derivatives using gel filtration. The results of electrophoretic and chromatographic analyses revealed that mGA-cobrotoxin comprised two modified derivatives, which contained modified Lys residues at positions 26 and 27 and at positions 26, 27, and 47, respectively. Moreover, an intramolecular cross-linking of loops II and III by Lys residues was noted with the monomeric derivative containing three modified Lys residues. In sharp contrast to cobrotoxin observations, the folding rate of mGA-cobrotoxin decreased in the presence of GSH/ GSSG, but notably increased in the absence of thiol compounds. Particularly, the accelerated effect of GSH/GSSG on the refolding reaction was affected by the presence of the intramolecular cross-link. Comparative analyses on cobrotoxin and mGA-cobrotoxin CD spectra revealed that modification with the GA reagent caused a change in the gross conformation of cobrotoxin. Fluorescence measurement revealed that the stability of the microenvironment around the single Trp-29 in mGA-cobrotoxin and unfolded mGA-cobrotoxin was appreciably higher than in cobrotoxin and unfolded toxin. Moreover, the ordered structure formation around Trp-29 in refolded mGA-cobrotoxin was faster than in refolded cobrotoxin as evidenced by fluorescence quenching studies. Taken together, these results suggest that the structural flexibility of unfolded cobrotoxin should be favorable for the thiol catalyst to exert its action in the refolding reaction after modification with GA.

Solution structure and backbone dynamics of an omega-conotoxin precursor

Nuclear magnetic resonance spectroscopy was used to characterize the solution structure and backbone dynamics of a putative precursor form of omega-conotoxin MVIIA, a 25-amino-acid residue peptide antagonist of voltage-gated Ca(2+) channels. The mature peptide is found in the venom of a fish-hunting marine snail Conus magus and contains an amidated carboxyl terminus that is generated by oxidative cleavage of a Gly residue. The form examined in this study is identical to the mature peptide except for the presence of the unmodified carboxy-terminal Gly. This form, referred to as omega-MVIIA-Gly, has previously been shown to refold and form its disulfides more efficiently than the mature form, suggesting that the presence of the terminal Gly may favor folding in vivo. The nuclear magnetic resonance (NMR) structure determination indicated that the fold of omega-MVIIA-Gly is very similar to that previously determined for the mature form, but revealed that the terminal Gly residue participates in a network of hydrogen bonds involving both backbone and side chain atoms, very likely accounting for the enhanced stability and folding efficiency. (15)N relaxation experiments indicated that the backbone is well ordered on the nanosecond time scale but that residues 9-15 undergo a conformational exchange processes with a time constant of approximately 35 microseconds. Other studies have implicated this segment in the binding of the peptide to its physiological target, and the observed motions may play a role in allowing the peptide to enter the binding site

The binding of a glycoprotein 120 V3 loop peptide to HIV-1 neutralizing antibodies. Structural implications

The structural and antigenic properties of a peptide ("CRK") derived from the V3 loop of HIV-1 gp120 protein were studied using NMR and SPR techniques. The sequence of CRK corresponds to the central portion of the V3 loop containing the highly conserved "GPGR" residue sequence. Although the biological significance of this conserved sequence is unknown, the adoption of conserved secondary structure (type II beta-turn) in this region has been proposed. The tendency of CRK (while free or conjugated to protein), to adopt such structure and the influence of such structure upon CRK antigenicity were investigated by NMR and SPR, respectively. Regardless of conjugation, CRK is conformationally averaged in solution but a weak tendency of the CRK "GPGR" residues to adopt a beta-turn conformation was observed after conjugation. The influence of GPGR structure upon CRK antigenicity was investigated by measuring the affinities of two cognate antibodies: "5023A" and "5025A," for CRK, protein-conjugated CRK and gp120 protein. Each antibody bound to all the antigens with nearly the same affinity. From these data, it appears that: (a) antibody binding most likely involves an induced fit of the peptide and (b) the gp120 V3 loop is probably conformationally heterogeneous. Since 5023A and 5025A are HIV-1 neutralizing antibodies, neutralization in these cases appears to be independent of adopted GPGR beta-turn structure.

Structural and dynamic characterization of omega-conotoxin MVIIA: the binding loop exhibits slow conformational exchange

omega-Conotoxin MVIIA is a 25-residue, disulfide-bridged polypeptide from the venom of the sea snail Conus magus that binds to neuronal N-type calcium channels. It forms a compact folded structure, presenting a loop between Cys8 and Cys15 that contains a set of residues critical for its binding. The loop does not have a unique defined structure, nor is it intrinsically flexible. Broadening of a subset of resonances in the NMR spectrum at low temperature, anomalous temperature dependence of the chemical shifts of some resonances, and exchange contributions to J(0) from (13)C relaxation measurements reveal that conformational exchange affects the residues in this loop. The effects of this exchange on the calculated structure of omega-conotoxin MVIIA are discussed. The exchange appears to be associated with a change in the conformation of the disulfide bridge Cys8-Cys20. The implications for the use of the omega-conotoxins as a scaffold for carrying other functions is discussed.

Leap-dynamics: efficient sampling of conformational space of proteins and peptides in solution

A molecular simulation scheme, called Leap-dynamics, that provides efficient sampling of protein conformational space in solution is presented. The scheme is a combined approach using a fast sampling method, imposing conformational 'leaps' to force the system over energy barriers, and molecular dynamics (MD) for refinement. The presence of solvent is approximated by a potential of mean force depending on the solvent accessible surface area. The method has been successfully applied to N-acetyl-L-alanine-N-methylamide (alanine dipeptide), sampling experimentally observed conformations inaccessible to MD alone under the chosen conditions. The method predicts correctly the increased partial flexibility of the mutant Y35G compared to native bovine pancreatic trypsin inhibitor. In particular, the improvement over MD consists of the detection of conformational flexibility that corresponds closely to slow motions identified by nuclear magnetic resonance techniques.

Tyrosine, phenylalanine, and disulfide contributions to the circular dichroism of proteins: circular dichroism spectra of wild-type and mutant bovine pancreatic trypsin inhibitor

Improved descriptions of the lowest energy excited states of tyrosine and phenylalanine side chains have been developed in order to extend the capabilities of calculating the circular dichroism (CD) spectra of proteins. Four transitions (Lb, La, Bb, and Ba) for each of the side-chain chromophores were considered, and the transition monopole charges were obtained from a CNDO/S calculation on models representing the individual groups. Monopole charges at midpoints of the bonds, corresponding to the maximum transition charge densities in the Lb band, and monopole charges representing the vibronic coupling with the B transitions for the La transition were also included. The aromatic transitions were combined with the peptide transitions (npi, pi0pi n'pi, and pi+pi) and disulfide transitions (n1sigma and n4sigma) in the framework of the origin-independent matrix method to compute the CD spectra of different crystal forms and Y --> L and F --> L mutants of bovine pancreatic trypsin inhibitor (BPTI). The structures of the mutants were obtained by replacing the appropriate tyrosine or phenylalanine residue by leucine in the wild-type crystal structure. The CD calculations were performed on the energy-minimized structures. The CD spectrum calculated for the form II crystal structure of BPTI showed the best agreement with experiment. In the far UV, the calculated and experimental CD spectra agree to various extents for the wild-type and mutant BPTI. Among the mutants, the calculated CD spectra of Y4L, Y10L, Y23L, and F45L showed reasonable agreement with experiment, while those of Y21L and F22L, the two residues interacting with most aromatic groups, showed poor agreement. In the near UV, the negative bands predicted for the wild-type and mutant BPTI have much less intensity than observed experimentally.