Events in the kinetic folding pathway of a small, all beta-sheet protein

Early Folding Events, Local Interactions, and Conservation of Protein Backbone Rigidity

Protein folding is in its early stages largely determined by the protein sequence and complex local interactions between amino acids, resulting in lower energy conformations that provide the context for further folding into the native state. We compiled a comprehensive data set of early folding residues based on pulsed labeling hydrogen deuterium exchange experiments. These early folding residues have corresponding higher backbone rigidity as predicted by DynaMine from sequence, an effect also present when accounting for the secondary structures in the folded protein. We then show that the amino acids involved in early folding events are not more conserved than others, but rather, early folding fragments and the secondary structure elements they are part of show a clear trend toward conserving a rigid backbone. We therefore propose that backbone rigidity is a fundamental physical feature conserved by proteins that can provide important insights into their folding mechanisms and stability.

Towards universal approach for bacterial production of three-finger Ly6/uPAR proteins: Case study of cytotoxin I from cobra N. oxiana

Cytotoxins or cardiotoxins is a group of polycationic toxins from cobra venom belonging to the 'three-finger' protein superfamily (Ly6/uPAR family) which includes small β-structural proteins (60-90 residues) with high disulfide bond content (4-5 disulfides). Due to a high cytotoxic activity for cancer cells, cytotoxins are considered as potential anticancer agents. Development of the high-throughput production methods is required for the prospective applications of cytotoxins. Here, efficient approach for bacterial production of recombinant analogue of cytotoxin I from N. oxiana containing additional N-terminal Met-residue (rCTX1) was developed. rCTX1 was produced in the form of E. coli inclusion bodies. Refolding in optimized conditions provided ∼6 mg of correctly folded protein from 1 L of bacterial culture. Cytotoxicity of rCTX1 for C6 rat glioma cells was found to be similar to the activity of wild type CTX1. The milligram quantities of ¹³C,¹⁵N-labeled rCTX1 were obtained. NMR study confirmed the similarity of the spatial structures of recombinant and wild-type toxins. Additional Met residue does not perturb the overall structure of the three-finger core. The analysis of available data for different Ly6/uPAR proteins of snake and human origin revealed that efficiency of their folding in vitro is correlated with the number of proline residues in the third loop and the surface area of hydrophobic residues buried within the protein interior. The obtained data indicate that hydrophobic core is important for the folding of proteins with high disulfide bond content. Developed expression method opens new possibilities for structure-function studies of CTX1 and other related three-finger proteins.

Amyloid-β peptide (1-42) aggregation induced by copper ions under acidic conditions

It is well known that the aggregation of amyloid-β peptide (Aβ) induced by Cu²⁺ is related to incubation time, solution pH, and temperature. In this work, the aggregation of Aβ₁₋₄₂ in the presence of Cu²⁺ under acidic conditions was studied at different incubation time and temperature (e.g. 25 and 37°C). Incubation temperature, pH, and the presence of Cu²⁺ in Aβ solution were confirmed to alter the morphology of aggregation (fibrils or amorphous aggregates), and the morphology is pivotal for Aβ neurotoxicity and Alzheimer disease (AD) development. The results of atomic force microscopy (AFM) indicated that the formation of Aβ fibrous morphology is preferred at lower pH, but Cu²⁺ induced the formation of amorphous aggregates. The aggregation rate of Aβ was increased with the elevation of temperature. These results were further confirmed by fluorescence spectroscopy and circular dichroism spectroscopy and it was found that the formation of β-sheet structure was inhibited by Cu²⁺ binding to Aβ. The result was consistent with AFM observation and the fibrillation process was restrained. We believe that the local charge state in hydrophilic domain of Aβ may play a dominant role in the aggregate morphology due to the strong steric hindrance. This research will be valuable for understanding of Aβ toxicity in AD.

Aggregation pathways of the amyloid β(1-42) peptide depend on its colloidal stability and ordered β-sheet stacking

Amyloid β (Aβ) fibrils are present as a major component in senile plaques, the hallmark of Alzheimer's disease (AD). Diffuse plaques (nonfibrous, loosely packed Aβ aggregates) containing amorphous Aβ aggregates are also formed in brain. This work examines the influence of Cu(2+) complexation by Aβ on the aggregation process in the context of charge and structural variations. Changes in the surface charges of Aβ molecules due to Cu(2+) binding, measured with a ζ-potential measurement device, were correlated with the aggregate morphologies examined by atomic force microscopy. As a result of the charge variation, the "colloid-like" stability of the aggregation intermediates, which is essential to the fibrillation process, is affected. Consequently, Cu(2+) enhances the amorphous aggregate formation. By monitoring variations in the secondary structures with circular dichroism spectroscopy, a direct transformation from the unstructured conformation to the β-sheet structure was observed for all types of aggregates observed (oligomers, fibrils, and/or amorphous aggregates). Compared to the Aβ aggregation pathway in the absence of Cu(2+) and taking other factors affecting Aβ aggregation (i.e., pH and temperature) into account, our investigation indicates that formations of amorphous and fibrous aggregates diverge from the same β-sheet-containing partially folded intermediate. This study suggests that the hydrophilic domain of Aβ also plays a role in the Aβ aggregation process. A kinetic model was proposed to account for the effects of the Cu(2+) binding on these two aggregation pathways in terms of charge and structural variations.

OneG: a computational tool for predicting cryptic intermediates in the unfolding kinetics of proteins under native conditions

Understanding the relationships between conformations of proteins and their stabilities is one key to address the protein folding paradigm. The free energy change (ΔG) of unfolding reactions of proteins is measured by traditional denaturation methods and native hydrogen-deuterium (H/D) exchange methods. However, the free energy of unfolding (ΔG(U)) and the free energy of exchange (ΔG(HX)) of proteins are not in good agreement, though the experimental conditions of both methods are well matching to each other. The anomaly is due to any one or combinations of the following reasons: (i) effects of cis-trans proline isomerisation under equilibrium unfolding reactions of proteins (ii) inappropriateness in accounting the baselines of melting curves (iii) presence of cryptic intermediates, which may elude the melting curve analysis and (iv) existence of higher energy metastable states in the H/D exchange reactions of proteins. Herein, we have developed a novel computational tool, OneG, which accounts the discrepancy between ΔG(U) and ΔG(HX) of proteins by systematically accounting all the four factors mentioned above. The program is fully automated and requires four inputs: three-dimensional structures of proteins, ΔG(U), ΔG(U)(*) and residue-specific ΔG(HX) determined under EX2-exchange conditions in the absence of denaturants. The robustness of the program has been validated using experimental data available for proteins such as cytochrome c and apocytochrome b(562) and the data analyses revealed that cryptic intermediates of the proteins detected by the experimental methods and the cryptic intermediates predicted by the OneG for those proteins were in good agreement. Furthermore, using OneG, we have shown possible existence of cryptic intermediates and metastable states in the unfolding pathways of cardiotoxin III and cobrotoxin, respectively, which are homologous proteins. The unique application of the program to map the unfolding pathways of proteins under native conditions have been brought into fore and the program is publicly available at http://sblab.sastra.edu/oneg.html.

Slow formation of aggregation-resistant beta-sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding-related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large beta-helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all-beta-sheet protein allows detailed analysis of the formation of beta-sheet structure in larger proteins. Using a combination of fluorescence and far-UV circular dichroism spectroscopy, we show that the pertactin beta-helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and beta-sheet-rich topology, pertactin refolding is reversible and not complicated by off-pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate-limiting step. Furthermore, site-specific labeling experiments indicate that the beta-helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, beta-sheet-rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, beta-sheet-rich refolding intermediates. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

INN-toxin, a highly lethal peptide from the venom of Indian cobra (Naja naja) venom-Isolation, characterization and pharmacological actions

A novel toxic polypeptide, INN-toxin, is purified from the venom of Naja naja using combination of gel-permeation and ion-exchange chromatography. It has a molecular mass of 6951.6Da as determined by MALDI-TOF/MS and the N-terminal sequence of LKXNKLVPLF. It showed both neurotoxic as well as cytotoxic activities. INN-toxin is lethal to mice with a LD(50) of 1.2mg/kg body weight. IgY raised in chicks against basic peptide pool neutralized the toxicity of INN-toxin. INN-toxin did not inhibit cholinesterase activity. It is toxic to Ehrlich ascites tumor (EAT) cells, but it is not toxic to leukocyte culture. The toxin appears to be specific in its mode of action. Interaction of N-bromosuccinamide (NBS) with the peptide resulted in the modification of tryptophan residues and loss of lethal toxicity of INN-toxin.

Predicting protein folding pathways at the mesoscopic level based on native interactions between secondary structure elements

Background

Since experimental determination of protein folding pathways remains difficult, computational techniques are often used to simulate protein folding. Most current techniques to predict protein folding pathways are computationally intensive and are suitable only for small proteins.

Results

By assuming that the native structure of a protein is known and representing each intermediate conformation as a collection of fully folded structures in which each of them contains a set of interacting secondary structure elements, we show that it is possible to significantly reduce the conformation space while still being able to predict the most energetically favorable folding pathway of large proteins with hundreds of residues at the mesoscopic level, including the pig muscle phosphoglycerate kinase with 416 residues. The model is detailed enough to distinguish between different folding pathways of structurally very similar proteins, including the streptococcal protein G and the peptostreptococcal protein L. The model is also able to recognize the differences between the folding pathways of protein G and its two structurally similar variants NuG1 and NuG2, which are even harder to distinguish. We show that this strategy can produce accurate predictions on many other proteins with experimentally determined intermediate folding states.

Conclusion

Our technique is efficient enough to predict folding pathways for both large and small proteins at the mesoscopic level. Such a strategy is often the only feasible choice for large proteins. A software program implementing this strategy (SSFold) is available at .

Investigation of the four cooperative unfolding units existing in the MICAL-1 CH domain

The solution structure of human MICAL-1 calpolnin homology (CH) domain is composed of six alpha helices and one 3(10) helix. To study the unfolding of this domain, we carry out native-state hydrogen exchange, intrinsic fluorescence and far-UV circular dichroism experiments. The free energy of unfolding, DeltaG(H2O), is calculated to be 7.11+/-0.58 kcal mol(-1) from GuHCl denaturation at pH 6.5. Four cooperative unfolding units are found using native-state hydrogen exchange experiment. Forty-seven slow-exchange residues can be studied by native-state hydrogen exchange experiments. From the concentration dependence of exchange rates, free energy of amide hydrogen with solvent, DeltaG(HX) and m-value (sensitivity of exposure to denaturant) are obtained, which reveal four cooperative unfolding units. The slowest exchanging protons are distributed throughout the whole hydrophobic core of the protein, which might be the folding core. These results will help us understand the structure of MICAL-1 CH domain more deeply.

The structural biology of protein aggregation diseases: Fundamental questions and some answers

Amyloid fibrils are found in association with at least two dozen fatal diseases. The tendency of numerous proteins to convert into amyloid-like fibrils poses fundamental questions for structural biology and for protein science in general. Among these are the following: What is the structure of the cross-beta spine, common to amyloid-like fibrils? Is there a sequence signature for proteins that form amyloid-like fibrils? What is the nature of the structural conversion from native to amyloid states, and do fibril-forming proteins have two distinct stable states, the native state and the amyloid state? What is the basis of protein complementarity, in which a protein chain can bind to itself? We offer tentative answers here, based on our own recent structural studies.

Molecular dynamics simulation of the enhancement of cobra cardiotoxin and E6 protein binding on mixed self-assembled monolayer molecules

Molecular dynamics simulations are performed on n-alkinethiol self-assembled monolayers (SAMs) and their mixture on a gold surface so that the orientations of the binding of cobra cardiotoxin and E6 protein molecules can be selected using the mixing ratio of CH3-terminated SAMs with different chain lengths. The simulations suggest that a SAM surface with different mixing ratios may provide a possible platform for aligning protein molecules with a desired orientation and for enhancing the binding energy of the protein on the designed surface.

Refolding of a small all beta-sheet protein proceeds with accumulation of kinetic intermediates

The refolding kinetics of Cobrotoxin (CBTX), a small all beta-sheet protein is investigated using a variety of biophysical techniques including quenched-flow hydrogen-deuterium (H/D) exchange in conjunction with two-dimensional NMR spectroscopy. Urea-induced equilibrium unfolding of CBTX follows a two-state mechanism with no distinct intermediates. The protein is observed to fold very rapidly within 250 ms. Both the refolding and the unfolding limbs of the chevron plot of CBTX show a prominent curvature suggesting the accumulation of kinetic intermediates. Quenched-flow H/D exchange data suggest the presence of a broad continuum of kinetic intermediates between the unfolded and native states of the protein. Comparison of the native state hydrogen exchange data and the results of the quenched-flow H/D exchange experiments, reveals that the residues constituting the folding core of CBTX are not a subset of the slow exchange core. To our knowledge, this is the first report wherein the refolding of a small all beta-sheet protein is shown to be a multi-step process involving the accumulation of kinetic intermediates.

Structure of the cross-beta spine of amyloid-like fibrils

Numerous soluble proteins convert to insoluble amyloid-like fibrils that have common properties. Amyloid fibrils are associated with fatal diseases such as Alzheimer's, and amyloid-like fibrils can be formed in vitro. For the yeast protein Sup35, conversion to amyloid-like fibrils is associated with a transmissible infection akin to that caused by mammalian prions. A seven-residue peptide segment from Sup35 forms amyloid-like fibrils and closely related microcrystals, from which we have determined the atomic structure of the cross-beta spine. It is a double beta-sheet, with each sheet formed from parallel segments stacked in register. Side chains protruding from the two sheets form a dry, tightly self-complementing steric zipper, bonding the sheets. Within each sheet, every segment is bound to its two neighbouring segments through stacks of both backbone and side-chain hydrogen bonds. The structure illuminates the stability of amyloid fibrils, their self-seeding characteristic and their tendency to form polymorphic structures.

Computation of conformational entropy from protein sequences using the machine-learning method-Application to the study of the relationship between structural conservation and local structural stability

A complete protein sequence can usually determine a unique conformation; however, the situation is different for shorter subsequences--some of them are able to adopt unique conformations, independent of context; while others assume diverse conformations in different contexts. The conformations of subsequences are determined by the interplay between local and nonlocal interactions. A quantitative measure of such structural conservation or variability will be useful in the understanding of the sequence-structure relationship. In this report, we developed an approach using the support vector machine method to compute the conformational variability directly from sequences, which is referred to as the sequence structural entropy. As a practical application, we studied the relationship between sequence structural entropy and the hydrogen exchange for a set of well-studied proteins. We found that the slowest exchange cores usually comprise amino acids of the lowest sequence structural entropy. Our results indicate that structural conservation is closely related to the local structural stability. This relationship may have interesting implications in the protein folding processes, and may be useful in the study of the sequence-structure relationship.

Monitoring protein folding at atomic resolution

Elucidation of the mechanisms by which proteins fold from disordered conformations to their unique native conformations is one of the most challenging tasks facing structural biologists. Understanding the mechanism(s) of protein folding involves the characterization of all structural species that occur in the protein-folding reaction. Nuclear magnetic resonance (NMR) spectroscopy is a powerful and versatile technique that provides an avenue to investigate the structures of the various conformational states at the residue level along the protein-folding reaction coordinate. In this Account, we provide a comprehensive review of the recent progress on the applications of NMR to monitor equilibrium and kinetic conformational states of the protein-folding reaction.

Investigation of the structural stability of the human acidic fibroblast growth factor by hydrogen-deuterium exchange

The conformational stability of the human acidic fibroblast growth factor (hFGF-1) is investigated using amide proton exchange and temperature-dependent chemical shifts, monitored by two-dimensional NMR spectroscopy. The change in free energy of unfolding (DeltaG(u)) of hFGF-1 is estimated to be 5.00 +/- 0.09 kcal.mol(-)(1). Amide proton-exchange rates of 74 residues (in hFGF-1) have been unambiguously measured, and the exchange process occurs predominately according to the conditions of the EX2 limit. The exchange rates of the fast-exchanging amide protons exposed to the solvent have been measured using the clean SEA-HSQC technique. The amide proton protection factor and temperature coefficient estimates show reasonably good correlation. Residues in beta-strands II and VI appear to constitute the stability core of the protein. Among the 12 beta-strands constituting the beta-barrel architecture of hFGF-1, beta-strand XI, located in the heparin binding domain, exhibits the lowest average protection factor value. Amide protons involved in the putative folding nucleation site in hFGF-1, identified by quench-flow NMR studies, do not represent the slow-exchanging core. Residues in portions of hFGF-1 experiencing high conformational flexibility mostly correspond to those involved in receptor recognition and binding.

Identification of rare partially unfolded states in equilibrium with the native conformation in an all beta-barrel protein

Human acidic fibroblast growth factor 1 (hFGF-1) is an all beta-barrel protein, and the secondary structural elements in the protein include 12 antiparallel beta-strands arranged into a beta-trefoil fold. In the present study, we investigate the stability of hFGF-1 by hydrogen-deuterium exchange as a function of urea concentration. Urea-induced equilibrium unfolding of hFGF-1 monitored by fluorescence and CD spectroscopy suggests that the protein unfolds by a two-state (native to denatured) mechanism. Hydrogen exchange in hFGF-1, under the experimental conditions used, occurs by the EX2 mechanism. In contrast to the equilibrium unfolding events monitored by optical probes, native state hydrogen exchange data show that the beta-trefoil architecture of hFGF-1 does not behave as a single cooperative unit. There are at least two structurally independent units with differing stabilities in hFGF-1. Beta-strands I, II, III, VI, VII, X, XI, and XII fit into the global unfolding isotherm. By contrast, residues in beta-strands IV, V, VIII, and IX exchange by the subfolding isotherm and could be responsible for the occurrence of high-energy partially unfolded state(s) in hFGF-1. There appears to be a broad continuum of stabilities among the four beta-strands (beta-strands IV, V, VIII, and IX) constituting the subglobal folding unit. The slow exchanging residues in hFGF-1 do not represent the folding nucleus of the protein.

Cytotoxic potency of cardiotoxin from Naja sputatrix: development of a new cytolytic assay

The possible involvement of specific regions/loops of cardiotoxin from Naja sputatrix venom in mediating its cytolytic activity is evaluated using a new cytolytic assay. In this assay, the amount of chloramphenicol acetyltransferase (CAT) that is released upon lysis of the cellular membranes by the cytotoxin has been measured as an index of cytolysis. This newly developed CAT system is more sensitive than the traditional haemolysis method utilizing red blood cells or the lactate dehydrogenase assay for cytolysis. Series of chimaeric toxin molecules have been constructed by swapping the loops between highly hydrophilic neurotoxin and highly hydrophobic cardiotoxin molecules from Naja sputatrix, which are known to exhibit structural similarity (three-finger conformation) but to have different functional properties. Comparison of the cytolytic activities of the recombinant chimaeric toxins demonstrated the possible involvement of all three loops of cardiotoxin in its cytolytic potency. However, the first two loops of the protein appear to make the major contribution to its lytic activity. cDNAs encoding cardiotoxin and the chimaeric toxins, when expressed in transfected cultured Chinese hamster ovary cells, resulted in cell lysis, indicating that these cDNAs can be developed as useful cytolytic agents.

Slow folding of three-fingered toxins is associated with the accumulation of native disulfide-bonded intermediates

Snake neurotoxins are short all-beta proteins that display a complex organization of the disulfide bonds: two bonds connect consecutive cysteine residues (C43-C54, C55-C60), and two bonds intersect when bridging (C3-C24, C17-C41) to form a particular structure classified as "disulfide beta-cross". We investigated the oxidative folding of a neurotoxin variant, named alpha62, to define the chemical nature of the three-disulfide intermediates that accumulate during the process in order to describe in detail its folding pathway. These folding intermediates were separated by reverse-phase HPLC, and their disulfide bonds were identified using a combination of tryptic hydrolysis, manual Edman degradation, and mass spectrometry. Two dominant intermediates containing three native disulfide bonds were identified, lacking the C43-C54 and C17-C41 pairing and therefore named des-[43-54] and des-[17-41], respectively. Both species were individually allowed to reoxidize under folding conditions, showing that des-[17-41] was a fast-forming nonproductive intermediate that had to interconvert into the des-[43-54] isomer before forming the native protein. Conversely, the des-[43-54] intermediate appeared to be the immediate precursor of the oxidized neurotoxin. A kinetic model for the folding of neurotoxin alpha62 which fits with the observed time-course accumulation of des-[17-41] and des-[43-54] is proposed. The effect of turn 2, located between residues 17 and 24, on the overall kinetics is discussed in view of this model.

Structural events during the refolding of an all beta-sheet protein

The refolding kinetics of the 140-residue, all beta-sheet, human fibroblast growth factor (hFGF-1) is studied using a variety of biophysical techniques such as stopped-flow fluorescence, stopped-flow circular dichroism, and quenched-flow hydrogen exchange in conjunction with multidimensional NMR spectroscopy. Urea-induced unfolding of hFGF-1 under equilibrium conditions reveals that the protein folds via a two-state (native <--> unfolded) mechanism without the accumulation of stable intermediates. However, measurement of the unfolding and refolding rates in various concentrations of urea shows that the refolding of hFGF-1 proceeds through accumulation of kinetic intermediates. Results of the quenched-flow hydrogen exchange experiments reveal that the hydrogen bonds linking the N- and C-terminal ends are the first to form during the refolding of hFGF-1. The basic beta-trefoil framework is provided by the simultaneous formation of beta-strands I, IV, IX, and X. The other beta-strands comprising the beta-barrel structure of hFGF-1 are formed relatively slowly with time constants ranging from 4 to 13 s.

Comparison of the structural stability of two homologous toxins isolated from the Taiwan cobra (Naja naja atra) venom

Cardiotoxin analogue III (CTX III) and cobrotoxin (CBTX) isolated from the Taiwan cobra venom (Naja naja atra) are structurally homologous, small molecular weight, all-beta-sheet proteins, cross-linked by four disulfide bonds at identical positions. The conformational stabilities of these toxins are compared based on temperature-dependent chemical shifts and amide proton exchange kinetics using two-dimensional NMR spectroscopy. The structure of CTX III is found to be significantly more stable than that of CBTX. In both the toxins, beta-strand III appears to constitute the stability core. In CTX III, the stability of the triple-stranded beta-sheet domain is observed to be markedly higher than the double-stranded beta-sheet segment. In contrast, in CBTX, both structural domains (double- and triple-stranded beta-sheet domains) appear to contribute equally to the stability of the protein. Estimation of the free energy of exchange (Delta G(ex)) of residues in CBTX and CTX III reveals that the enhanced stability of the structure of CTX III stems from the strong interactions among the beta-strands constituting the triple-stranded beta-sheet domain and also the molecular forces bridging the residues at the N- and C-terminal ends of the molecule.

Elucidation of the solution structure of cardiotoxin analogue V from the Taiwan cobra (Naja naja atra)--identification of structural features important for the lethal action of snake venom cardiotoxins

The aim of the present study is to understand the structural features responsible for the lethal activity of snake venom cardiotoxins. Comparison of the lethal potency of the five cardiotoxin isoforms isolated from the venom of Taiwan cobra (Naja naja atra) reveals that the lethal potency of CTX I and CTX V are about twice of that exhibited by CTX II, CTX III, and CTX IV. In the present study, the solution structure of CTX V has been determined at high resolution using multidimensional proton NMR spectroscopy and dynamical simulated annealing techniques. Comparison of the high resolution solution structures of CTX V with that of CTX IV reveals that the secondary structural elements in both the toxin isoforms consist of a triple and double-stranded antiparallel beta-sheet domains. Critical examination of the three-dimensional structure of CTX V shows that the residues at the tip of Loop III form a distinct "finger-shaped" projection comprising of nonpolar residues. The occurrence of the nonpolar "finger-shaped" projection leads to the formation of a prominent cleft between the residues located at the tip of Loops II and III. Interestingly, the occurrence of a backbone hydrogen bonding (Val27CO to Leu48NH) in CTX IV is found to distort the "finger-shaped" projection and consequently diminish the cleft formation at the tip of Loops II and III. Comparison of the solution structures and lethal potencies of other cardiotoxin isoforms isolated from the Taiwan cobra (Naja naja atra) venom shows that a strong correlation exists between the lethal potency and occurrence of the nonpolar "finger-shaped" projection at the tip of Loop III. Critical analysis of the structures of the various CTX isoforms from the Taiwan cobra suggest that the degree of exposure of the cationic charge (to the solvent) contributed by the invariant lysine residue at position 44 on the convex side of the CTX molecules could be another crucial factor governing their lethal potency.

Investigation of the structural stability of cardiotoxin analogue III from the Taiwan cobra by hydrogen-deuterium exchange kinetics

The conformational stability of a small ( approximately 7 kDa), all beta-sheet protein, cardiotoxin analogue III (CTX III), from the venom of the Taiwan cobra has been investigated by hydrogen-deuterium (H/D) exchange using two-dimensional NMR spectroscopy. The H/D exchange kinetics of backbone amide protons in CTX III has been monitored at pD 3.6 and 6.6 (at 25 degrees C), for over 5000 h. Examination of H/D exchange kinetics in the protein showed that a number of slowly exchanging residues are in the hydrophobic core of the protein. The average protection factor of the amide protons of residues belonging to the triple-stranded beta-sheet domain is about 20 times greater than that of those in the double-stranded beta-sheet segment. The residues in the C-terminal tail of the molecule, though structureless, have been found to exhibit significant protection against H/D exchange. Comparison of the quenched-flow H/D exchange data on CTX III with those obtained in the present study reveals that the most slowly exchanging portion constitutes the folding core of the protein.

The hydrogen exchange core and protein folding

A database of hydrogen-deuterium exchange results has been compiled for proteins for which there are published rates of out-exchange in the native state, protection against exchange during folding, and out-exchange in partially folded forms. The question of whether the slow exchange core is the folding core (Woodward C, 1993, Trends Biochem Sci 18:359-360) is reexamined in a detailed comparison of the specific amide protons (NHs) and the elements of secondary structure on which they are located. For each pulsed exchange or competition experiment, probe NHs are shown explicitly; the large number and broad distribution of probe NHs support the validity of comparing out-exchange with pulsed-exchange/competition experiments. There is a strong tendency for the same elements of secondary structure to carry NHs most protected in the native state, NHs first protected during folding, and NHs most protected in partially folded species. There is not a one-to-one correspondence of individual NHs. Proteins for which there are published data for native state out-exchange and theta values are also reviewed. The elements of secondary structure containing the slowest exchanging NHs in native proteins tend to contain side chains with high theta values or be connected to a turn/loop with high theta values. A definition for a protein core is proposed, and the implications for protein folding are discussed. Apparently, during folding and in the native state, nonlocal interactions between core sequences are favored more than other possible nonlocal interactions. Other studies of partially folded bovine pancreatic trypsin inhibitor (Barbar E, Barany G, Woodward C, 1995, Biochemistry 34:11423-11434; Barber E, Hare M, Daragan V, Barany G, Woodward C, 1998, Biochemistry 37:7822-7833), suggest that developing cores have site-specific energy barriers between microstates, one disordered, and the other(s) more ordered.

Secondary structure formation is the earliest structural event in the refolding of an all beta-sheet protein

The refolding kinetics of cobrotoxin (CBTX), a small-molecular-weight ( approximately 7 kDa) all beta-sheet protein, has been monitored using a variety of biophysical techniques. The secondary structure formation and hydrophobic collapse occur as distinct events during the refolding of the protein. Complete secondary structure formation occurs prior to the clustering of the hydrophobic residues. The late stage(s) of the refolding pathway of CBTX is characterized by change(s) in the local environment and optical asymmetry of the indole ring of the sole tryptophan residue. The results obtained in the present study, to our knowledge, represent the first unambiguous experimental support for the framework model of protein folding.

Binding of nucleotide triphosphates to cardiotoxin analogue II from the Taiwan cobra venom (Naja naja atra). Elucidation of the structural interactions in the dATP-cardiotoxin analogue ii complex

Snake venom cardiotoxins have been recently shown to block the enzymatic activity of phospholipid protein kinase and Na+,K+-ATPase. To understand the molecular basis for the inhibitory effects of cardiotoxin on the action of these enzymes, the nucleotide triphosphate binding ability of cardiotoxin analogue II (CTX II) from the Taiwan cobra (Naja naja atra) venom is investigated using a variety of spectroscopic techniques such as fluorescence, circular dichroism, and two-dimensional NMR. CTX II is found to bind to all the four nucleotide triphosphates (ATP, UTP, GTP, and CTP) with similar affinity. Detailed studies of the binding of dATP to CTX II indicated that the toxin molecule is significantly stabilized in the presence of the nucleotide. Molecular modeling, based on the NOEs observed for the dATP.CTX II complex, reveals that dATP binds to the CTX II molecule at the groove enclosed between the N- and C-terminal ends of the toxin molecule. Based on the results obtained in the present study, a molecular mechanism to account for the inhibition of the enzymatic activity of the phospholipid-sensitive protein kinase and Na+,K+-ATPase is also proposed.

Influence of disulfide bonds on the induction of helical conformation in proteins

The effect(s) of TFE (2,2,2-trifluoroethanol) on three different conformational states (native, denatured, and carboxymethylated) of CTX III and RNase A has been examined. Contrary to the general belief, the results of the present study reveal that TFE can induce helical conformation in a protein which has no sequence propensity to form a helix. It is found that the helix induction in TFE is intricately related to the destabilization of the tertiary structural conformation in proteins. More importantly, the disulfide bonds in proteins are found to have significant influence on the TFE-mediated helix induction. The results obtained in this study strongly suggest that information pertaining to the influence of disulfide bonds on helix induction need to be considered to improve the accuracy of secondary structure prediction algorithms.

Structurally homologous toxins isolated from the Taiwan cobra (Naja naja atra) differ significantly in their structural stability

Cardiotoxin and neurotoxin analogues isolated from snake venom sources are highly homologous proteins (>50% homology) with similar three-dimensional structures but exhibit drastically different biological properties. In the present study, we compare the conformational stability of cardiotoxin analogue III (CTX III) and cobrotoxin (CBTX), a neurotoxin analogue, from the Taiwan cobra (Naja naja atra), using circular dichroism spectroscopy and hydrogen-deuterium (H/D) exchange techniques in conjunction with two-dimensional NMR methods. Contrary to expectations, it is found that CTX III and CBTX differ significantly in their structural stabilities. The three-dimensional structure of CBTX is less stable than that of CTX III. The amide protons of residues at the N- and C-terminal ends of the CTX III molecule are strongly protected against H/D exchange, implying that the terminal ends of the molecule are bridged together by significant numbers of hydrogen bonds. However, in CBTX, amide protons at the terminal ends of the molecule do not exhibit an significant protection against H/D exchange. Comparison of the protection factors of the various amide protons in CTX III and CBTX reveals that the extraordinary stability of CTX III stems from the strong network of interactions among the residues at the N- and C-terminal ends and also due to the tight and ordered packing of the nonpolar residues involved in the triple-stranded, anti-parallel, beta-sheet segment of the molecule.

The length of a single turn controls the overall folding rate of "three-fingered" snake toxins

Snake curaremimetic toxins are short all-beta proteins, containing several disulfide bonds which largely contribute to their stability. The four disulfides present in snake toxins make a "disulfide beta-cross"-fold that was suggested to be a good protein folding template. Previous studies on the refolding of snake toxins (Ménez, A. et al. (1980) Biochemistry 19, 4166-4172) showed that this set of natural homologous proteins displays different rates of refolding. These studies suggested that the observed different rates could be correlated to the length of turn 2, one out of five turns present in the toxins structure and close to the "disulfide beta-cross". To demonstrate this hypothesis, we studied the refolding pathways and kinetics of two natural isotoxins, toxin alpha (Naja nigricollis) and erabutoxin b (Laticauda semifasciata), and two synthetic homologues, the alpha mutants, alpha60 and alpha62. These mutants were designed to probe the peculiar role of the turn 2 on the refolding process by deletion or insertion of one residue in the turn length that reproduced the natural heterogeneity at that locus. The refolding was studied by electrospray mass spectrometry (ESMS) time-course analysis. This analysis permitted both the identification and quantitation of the population of intermediates present during the process. All toxins were shown to share the same sequential scheme for disulfide bond formation despite large differences in their refolding rates. The results presented here demonstrate definitely that no residues except those forming turn 2 accounted for the observed differences in the refolding rate of toxins.