Disulfide Bond Assignments by Mass Spectrometry of Native Natural Peptides: Cysteine Pairing in Disulfide Bonded Conotoxins

Contryphan sequence diversity: Messy N-terminus processing, effects on chromatographic behaviour and mass spectrometric fragmentation

Contryphans, peptides containing a single disulfide bond, are found abundantly in cone snail venom. The analysis of a large dataset of available contryphan sequences permits a classification based on the occurrence of proline residues at positions 2 and 5 within the macrocyclic 23-membered disulfide loop. Further sequence diversity is generated by variable proteolytic processing of the contryphan precursor proteins. In the majority of contryphans, presence of Pro at position 2 and a D-residue at position 3 leads to a slow conformational dynamics, manifesting as anomalous chromatographic profiles during LC analysis. LC-MS analysis of diverse contryphans suggests that elution profiles may be used as a rapid diagnostic for the presence of the Pro2-^DXxx3 motif. Natural sequences from C.inscriptus and C.frigidus together with synthetic analogs permit the delineation of the features necessary for abnormal chromatographic behaviour. A diagnostic for the presence of Pro at position 5 is obtained by the observation of non-canonical fragment ions, generated by N-C^α bond cleavage at the dehydroalanine residue formed by disulfide cleavage. Anomalous LC profiles supports Pro at position 2, while non-canonical mass spectral fragments established Pro at position 5, providing a rapid method for contryphan analysis from LC-ESI-MS/MS profiles of crude Conus venom. SIGNIFICANCE: Contryphans are peptides, widely distributed in cone snail venom, which display extensive sequence diversity. Heterogeneity of proteolytic processing of contryphan precursor proteins, together with post-translational modifications contributes to contryphan diversity. Contryphans, identified by a combination of mass spectrometry and transcriptomic analysis, are classified on the basis of sequence features, primarily the number of proline residues within the disulfide loop. Conformational diversity arises in contryphans by cis-trans isomerization of Cys-Pro bonds, resulting in characteristic chromatographic profiles, permitting identification even in crude venom mixtures. Rapid identification of contryphans in cone snail peptide libraries is also facilitated by diagnostic mass spectral fragments arising by non-canonical cleavage of the N-C^α bond at Cys(7).

Mass spectral and theoretical investigations of N-Cα bond cleavages in the disulfide-containing peptide TTCPYCKK and its analogues

RATIONALE

For disulfide-containing peptides, mass spectrometric analyses are rarely comparably studied between their dithiol and disulfide forms. Persulfide ions afforded from peptides with a disulfide ring are from either an unusual N-C^α bond cleavage or a canonical peptide bond cleavage; their isomeric structures are, however, not identified just from peaks of mass spectra.

METHODS

Isomeric structures of [C³ P⁴ X⁵ |C⁶ _M ], [C³ _MA P⁴ X⁵ |C⁶ _MB ] and [P⁴ X⁵ C⁶ |C³ _M ] were identified from a series of the X⁵ substituted dicysteine octapeptides using electrospray ionization tandem mass spectrometry for both their dithiol and disulfide forms. Formation mechanisms of different persulfide ions were investigated systematically by theoretical methods. Moreover, electrostatic potential-mapped molecular van der Waals surfaces were used to determine the stabilities of the intermediates, which gave a further evaluation of favored bond cleavage.

RESULTS

Mass spectral analyses indicated that the fragmented ions changed largely when an intramolecular disulfide bond was formed. New types of disulfide-containing fragmented ions [C³ P⁴ X⁵ |C⁶ _M ] or [C³ _MA P⁴ X⁵ |C⁶ _MB ] were thus proposed. Energy analysis showed that the N-C^α cleavage was not competitive energetically with that of the amide bond for Y⁵ and its phosphorylated analogue. However, the N-C^α cleavage products dominated for the S⁵ - and T⁵ -containing peptides. Stabilities of the intermediates were found to be related with the electrostatic potential-mapped molecular van der Waals surfaces.

CONCLUSIONS

Persulfide ions containing more residues than previously found were proposed not only from b₇ ions but also from y₆ ions. In addition, a new kind of phosphorylated analogue, [C³ P⁴ _p Y⁵ |C⁶ _M ], is reported in this work. Our study provides convincing results for separating isomeric structures in the cases of N-C^α cleavages, which greatly assists in the structural identification of disulfide-containing peptides.

Mass Spectral Signatures of Complex Post-Translational Modifications in Proteins: A Proof-of-Principle Based on X-ray Irradiated Vancomycin

Characterizing post-translational modifications (PTM) of proteins is of key relevance for the understanding of many biological processes, as these covalent modifications strongly influence or even determine protein function. Among the different analytical techniques available, mass spectrometry is attracting growing attention because recent instrumental and computational improvements have led to a massive rise of the number of PTM sites that can be identified and quantified. However, multiple PTM occurring at adjacent amino acid residues can lead to complex and dense chemical patterns that are a challenge to characterize. By means of X-ray synchrotron radiation coupled to mass spectrometry, and through the test-case of the glycopeptide antibiotic vancomycin, we show that such a pattern has a unique and robust signature in terms of photon energy and molecular environment. This highlights the potential of this technique in proteomics and its value as a tool to understand the biological roles of PTM.

Mechanistic Insights into an Unusual Side-Chain-Mediated N-Cα Bond Cleavage under Collision-Induced Dissociation Conditions in the Disulfide-Containing Peptide Conopressin

Conopressin, a nonapeptide disulfide CFIRNCPKG amide present in cone snail venom, undergoes a facile cleavage at the Cys6-Pro7 peptide bond to yield a disulfide bridged b₆ ion. Analysis of the mass spectral fragmentation pattern reveals the presence of a major fragment ion, which is unambiguously assigned as the tripeptide sequence IRN amide. The sequence dependence of this unusual fragmentation process has been investigated by comparing it with the fragmentation patterns of related peptides, oxytocin (CYIQNCPLG amide), Lys-vasopressin (CYFQNCPKG amide), and a series of synthetic analogues. The results establish the role of the Arg4 residue in facilitating the unusual N-C^α bond cleavage at Cys6. Structures are proposed for a modified disulfide bridged fragment containing the Cys1 and Cys6 residues. Gas-phase molecular dynamics simulations provide evidence for the occurrence of conformational states that permit close approach of the Arg4 side chain to the Cys6 C^β methylene protons.

Structural plasticity of mini-M conotoxins - expression of all mini-M subtypes by Conus regius

The mini-M conotoxins are peptidic scaffolds found in the venom of cones snails. These scaffolds are tightly folded structures held together by three disulfide bonds with a CC-C-C-CC arrangement (conotoxin framework III) and belong to the M Superfamily of conotoxins. Here, we describe mini-M conotoxins from the venom of Conus regius, a Western Atlantic worm-hunting cone snail species using transcriptomic and peptidomic analyses. These C. regius conotoxins belong to three different subtypes: M1, M2, and M3. The subtypes show little sequence homology, and their loop sizes (intercysteine amino acid chains) vary significantly. The mini-Ms isolated from dissected venom contains preferentially hydroxylated proline residues, thus augmenting the structural reach of this conotoxin class. Using 2D-NMR methods, we have determined the 3D structure of reg3b, an M2 subtype conotoxin, which shows a constrained multi-turn scaffold. The structural diversity found within mini-M conotoxin scaffolds of C. regius is indicative of structural hypervariability of the conotoxin M superfamily that is not seen in other superfamilies. These stable minimalistic scaffolds may be investigated for the development of engineered peptides for therapeutic applications.

DATABASES

Sequences are available in GenBank under accession numbers MF588935-MF588952. Structural data are available in the RCSB protein database under the accession code 6BX9.

Aqueous ionic liquids influence the disulfide bond isoform equilibrium in conotoxin AuIB: a consequence of the Hofmeister effect?

The appearance of several disulfide bond isoforms in multiple cysteine containing venom peptides poses a significant challenge in their synthesis and purification under laboratory conditions. Recent experiments suggest that careful tuning of solvent and temperature conditions can propel the disulfide bond isoform equilibrium in favor of the most potent, native form. Certain aqueous ionic liquids (ILs) have proven significantly useful as solvents for this purpose, while exceptions have also been noted. To elucidate the molecular level origin behind such a preference, we report a detailed explicit solvent replica exchange molecular dynamics study of a conotoxin, AuIB, in pure water and four different aqueous IL solutions (~45-60% v/v). The ILs studied here are comprised of cations like 1-ethyl-3-methyl-imidazolium (Im₂₁⁺) or 1-butyl-3-methyl-imidazolium (Im₄₁⁺) coupled with either acetate (OAc^-) or chloride (Cl^-) as the counter anion. Our simulations unfold interesting features of the conformational spaces sampled by the peptide and its solvation in pure water and aqueous IL solutions. Detailed investigation into populations of the globular disulfide bond isoform of AuIB in aqueous IL solutions reveal distinct trends which might be related to the Hofmeister effect of the cation and anion of the IL and of specific interactions of the aqueous IL solutions with the peptide. In accordance with experimental observations, the aqueous [Im₂₁][OAc] solution is found to promote the highest globular isoform population in AuIB.

Im10A, a short conopeptide isolated from Conus imperialis and possesses two highly concentrated disulfide bridges and analgesic activity

In the present study, we isolated, synthesized and NMR structurally characterized a novel conopeptide Im10A consisting of 11 amino acids (NTICCEGCMCY-NH2) from Conus imperialis. Unlike other conopeptides with four cysteine residues, Im10A had only two residues in loop 1 and one residue in loop 2 (CC-loop1-C-loop2-C), which formed a stable disulfide connectivity "I-IV, II- III" (framework X) with a type I β-turn. Interestingly, Im10A exhibited 50.7% analgesic activity on rat partial sciatic nerve ligation (PNL) at 2h after Im10A administration. However, 10μM Im10A exhibited no apparent effect on neuronal nicotinic acetylcholine receptor, and it did not target DRG voltage-dependent sodium, potassium and calcium ion channels and opioid receptor. To our knowledge, Im10A had the most concentrated disulfide bridges among conopeptides with four cysteine residues. This finding provided a new motif for the future development of biomimetic compounds.

A Disulfide Stabilized β-Sandwich Defines the Structure of a New Cysteine Framework M-Superfamily Conotoxin

The structure of a new cysteine framework (-C-CC-C-C-C-) "M"-superfamily conotoxin, Mo3964, shows it to have a β-sandwich structure that is stabilized by inter-sheet cross disulfide bonds. Mo3964 decreases outward K(+) currents in rat dorsal root ganglion neurons and increases the reversal potential of the NaV1.2 channels. The structure of Mo3964 (PDB ID: 2MW7 ) is constructed from the disulfide connectivity pattern, i.e., 1-3, 2-5, and 4-6, that is hitherto undescribed for the "M"-superfamily conotoxins. The tertiary structural fold has not been described for any of the known conus peptides. NOE (549), dihedral angle (84), and hydrogen bond (28) restraints, obtained by measurement of (h3)JNC' scalar couplings, were used as input for structure calculation. The ensemble of structures showed a backbone root mean square deviation of 0.68 ± 0.18 Å, with 87% and 13% of the backbone dihedral (ϕ, ψ) angles lying in the most favored and additional allowed regions of the Ramachandran map. The conotoxin Mo3964 represents a new bioactive peptide fold that is stabilized by disulfide bonds and adds to the existing repertoire of scaffolds that can be used to design stable bioactive peptide molecules.

A sleep-inducing peptide from the venom of the Indian cone snail Conus araneosus

The marine snail Conus araneosus has unusual significance due to its confined distribution to coastal regions of southeast India and Sri Lanka. Due to its relative scarceness, this species has been poorly studied. In this work, we characterized the venom of C. araneosus to identify new venom peptides. We identified 14 novel compounds. We determined amino acid sequences from chemically-modified and unmodified crude venom using liquid chromatography-electrospray ionization mass spectrometry and matrix assisted laser desorption ionization time-of-flight mass spectrometry. Ten sequences showed six Cys residues arranged in a pattern that is most commonly associated with the M-superfamily of conotoxins. Four other sequences had four Cys residues in a pattern that is most commonly associated with the T-superfamily of conotoxins. The post-translationally modified residue (pyroglutamate) was determined at the N-terminus of two sequences, ar3h and ar3i respectively. In addition, two sequences, ar3g and ar3h were C-terminally amidated. At a dose of 2 nmol, peptide ar3j elicited sleep when injected intraperitoneally into mice. To our knowledge, this is the first report of a peptide from a molluscivorous cone snail with sleep-inducing effects in mice. The novel peptides characterized herein extend the repertoire of unique peptides derived from cone snails and may add value to the therapeutic promise of conotoxins.

Fragmentation follows structure: top-down mass spectrometry elucidates the topology of engineered cystine-knot miniproteins

Over the last decades the field of pharmaceutically relevant peptides has enormously expanded. Among them, several peptide families exist that contain three or more disulfide bonds. In this context, elucidation of the disulfide patterns is extremely important as these motifs are often prerequisites for folding, stability, and activity. An example of this structure-determining pattern is a cystine knot which comprises three constrained disulfide bonds and represents a core element in a vast number of mechanically interlocked peptidic structures possessing different biological activities. Herein, we present our studies on disulfide pattern determination and structure elucidation of cystine-knot miniproteins derived from Momordica cochinchinensis peptide MCoTI-II, which act as potent inhibitors of human matriptase-1. A top-down mass spectrometric analysis of the oxidised and bioactive peptides is described. Following the detailed sequencing of the peptide backbone, interpretation of the MS(3) spectra allowed for the verification of the knotted topology of the examined miniproteins. Moreover, we found that the fragmentation pattern depends on the knottin's folding state, hence, tertiary structure, which to our knowledge has not been described for a top-down MS approach before.

Characterization of a dual-active enzyme, DcpA, involved in cyclic diguanosine monophosphate turnover in Mycobacterium smegmatis

We have reported previously that the long-term survival of Mycobacterium smegmatis is facilitated by a dual-active enzyme MSDGC-1 (renamed DcpA), which controls the cellular turnover of cyclic diguanosine monophosphate (c-di-GMP). Most mycobacterial species possess at least a single copy of a DcpA orthologue that is highly conserved in terms of sequence similarity and domain architecture. Here, we show that DcpA exists in monomeric and dimeric forms. The dimerization of DcpA is due to non-covalent interactions between two protomers that are arranged in a parallel orientation. The dimer shows both synthesis and hydrolysis activities, whereas the monomer shows only hydrolysis activity. In addition, we have shown that DcpA is associated with the cytoplasmic membrane and exhibits heterogeneous cellular localization with a predominance at the cell poles. Finally, we have also shown that DcpA is involved in the change in cell length and colony morphology of M. smegmatis. Taken together, our study provides additional evidence about the role of the bifunctional protein involved in c-di-GMP signalling in M. smegmatis.

Strategies in mass spectrometry for the assignment of Cys-Cys disulfide connectivities in proteins

Elucidating disulfide linkage patterns is a crucial part of protein characterization, for which mass spectrometry (MS) is now an indispensable analytical tool. In many cases, MS-based disulfide connectivity assignment is straightforwardly achieved using one-step protein fragmentation in the unreduced form followed by mass measurement of bridged fragments. By contrast, venom proteins, which are receiving increasing interest as potential therapeutics, are a challenge for MS-based disulfide assignment due to their numerous closely spaced cysteines and knotted disulfide structure, requiring creative strategies to determine their connectivity. Today, these include the use of an array of reagents for enzymatic and/or chemical cleavage, partial reduction, differential cysteine labeling and tandem MS. This review aims to describe the toolkit of techniques available to MS users approaching both straightforward and complex disulfide bridge assignments, with a particular focus on strategies utilizing standard instrumentation found in a well-equipped analytical or proteomics laboratory.

Collision-induced dissociation fragmentation inside disulfide C-terminal loops of natural non-tryptic peptides

Collision-induced dissociation (CID) spectra of long non-tryptic peptides are usually quite complicated and rather difficult to interpret. Disulfide bond formed by two cysteine residues at C-terminus of frog skin peptides precludes one to determine sequence inside the forming loop. Thereby, chemical modification of S-S bonds is often used in "bottom up" sequencing approach. However, low-energy CID spectra of natural non-tryptic peptides with C-terminal disulfide cycle demonstrate an unusual fragmentation route, which may be used to elucidate the "hidden" C-terminal sequence. Low charge state protonated molecules experience peptide bond cleavage at the N-terminus of C-terminal cysteine. The forming isomeric acyclic ions serve as precursors for a series of b-type ions revealing sequence inside former disulfide cycle. The reaction is preferable for peptides with basic lysine residues inside the cycle. It may also be activated by acidic protons of Asp and Glu residues neighboring the loop. The observed cleavages may be quite competitive, revealing the sequence inside disulfide cycle, although S-S bond rupture does not occur in this case.

Resolving double disulfide bond patterns in SNAP25B using liquid chromatography-ion trap mass spectrometry

Complex disulfide bond patterns in synaptosomal-associated protein of 25 kD B (SNAP25B) are thought to regulate neurotransmitter release in response to oxidative stress. However, the steric feasibility of each possible disulfide pattern in SNAP25B has not been assessed. To assess the steric feasibility of hypothesized closely spaced complex disulfide patterning in SNAP25B and also the feasibility of identifying complex disulfide bond patterns with MS, we have developed a novel probabilistic analysis to unambiguously resolve complex double disulfide bond patterns by using an ion trap mass spectrometer. We analyzed fragmentation patterns of singly linked peptides to determine likely fragmentation events in an ion trap mass spectrometer and observed double and single backbone cleavage along with heterolytic cleavage of the disulfide bond. We modeled these same events in the doubly disulfide linked SNAP25B peptide and used a cumulative hypergeometric distribution with top-down scoring to both identify and differentiate these bonding patterns. Because of the presence of unique MS/MS peaks, two of the bonding patterns were directly identified. The third was assigned on the basis of full chromatographic separation and confirmed by modeling triple breakage fragments. In total, this work demonstrates the feasibility--and also limitations--of identification of complex intradisulfide patterns by using ion trap-based collision-induced dissociation-based fragmentation methods.

Mammalian neuronal sodium channel blocker μ-conotoxin BuIIIB has a structured N-terminus that influences potency

Among the μ-conotoxins that block vertebrate voltage-gated sodium channels (VGSCs), some have been shown to be potent analgesics following systemic administration in mice. We have determined the solution structure of a new representative of this family, μ-BuIIIB, and established its disulfide connectivities by direct mass spectrometric collision induced dissociation fragmentation of the peptide with disulfides intact. The major oxidative folding product adopts a 1-4/2-5/3-6 pattern with the following disulfide bridges: Cys5-Cys17, Cys6-Cys23, and Cys13-Cys24. The solution structure reveals that the unique N-terminal extension in μ-BuIIIB, which is also present in μ-BuIIIA and μ-BuIIIC but absent in other μ-conotoxins, forms part of a short α-helix encompassing Glu3 to Asn8. This helix is packed against the rest of the toxin and stabilized by the Cys5-Cys17 and Cys6-Cys23 disulfide bonds. As such, the side chain of Val1 is located close to the aromatic rings of Trp16 and His20, which are located on the canonical helix that displays several residues found to be essential for VGSC blockade in related μ-conotoxins. Mutations of residues 2 and 3 in the N-terminal extension enhanced the potency of μ-BuIIIB for NaV1.3. One analogue, [d-Ala2]BuIIIB, showed a 40-fold increase, making it the most potent peptide blocker of this channel characterized to date and thus a useful new tool with which to characterize this channel. On the basis of previous results for related μ-conotoxins, the dramatic effects of mutations at the N-terminus were unanticipated and suggest that further gains in potency might be achieved by additional modifications of this region.

Structural characterization of a recombinant fusion protein by instrumental analysis and molecular modeling

Conbercept is a genetically engineered homodimeric protein for the treatment of wet age-related macular degeneration (wet AMD) that functions by blocking VEGF-family proteins. Its huge, highly variable architecture makes characterization and development of a functional assay difficult. In this study, the primary structure, number of disulfide linkages and glycosylation state of conbercept were characterized by high-performance liquid chromatography, mass spectrometry, and capillary electrophoresis. Molecular modeling was then applied to obtain the spatial structural model of the conbercept–VEGF-A complex, and to study its inter-atomic interactions and dynamic behavior. This work was incorporated into a platform useful for studying the structure of conbercept and its ligand binding functions.

Simple approach to assign disulfide connectivity using extracted ion chromatograms of electron transfer dissociation spectra

Increasing interest in production of protein-based pharmaceuticals (biotherapeutics) is accompanied by an increased need for verification of protein folding and correct disulfide bonding. Recombinant protein expression may produce aberrant disulfide bonds and could result in safety concerns or decreased efficacy. Thus, the thorough analysis of disulfide bonding is a necessity for protein therapeutics. The use of electron transfer dissociation (ETD) facilitates this analysis because disulfide bonds are preferentially cleaved when subjected to ETD. Here, we make use of this well-characterized reaction to assign disulfide bonding networks by coupling the use of extracted ion chromatograms (XICs) of cysteine-containing peptides with ETD analysis to produce an efficient assignment approach for disulfide bonding. This method can be used to assign a disulfide pattern in a de novo fashion, to detect disulfide shuffling, and to provide information on heterogeneity, when more than one disulfide bonding pattern is present. The method was applied for assigning the disulfide-bonding network of a recombinant monomer of the HIV envelope protein gp120. It was found that one region of the protein, the V1/V2 loops, had significant heterogeneity in the disulfide bonds.

Distinct disulfide isomers of μ-conotoxins KIIIA and KIIIB block voltage-gated sodium channels

In the preparation of synthetic conotoxins containing multiple disulfide bonds, oxidative folding can produce numerous permutations of disulfide bond connectivities. Establishing the native disulfide connectivities thus presents a significant challenge when the venom-derived peptide is not available, as is increasingly the case when conotoxins are identified from cDNA sequences. Here, we investigate the disulfide connectivity of μ-conotoxin KIIIA, which was predicted originally to have a [C1-C9,C2-C15,C4-C16] disulfide pattern based on homology with closely related μ-conotoxins. The two major isomers of synthetic μ-KIIIA formed during oxidative folding were purified and their disulfide connectivities mapped by direct mass spectrometric collision-induced dissociation fragmentation of the disulfide-bonded polypeptides. Our results show that the major oxidative folding product adopts a [C1-C15,C2-C9,C4-C16] disulfide connectivity, while the minor product adopts a [C1-C16,C2-C9,C4-C15] connectivity. Both of these peptides were potent blockers of Na(V)1.2 (K(d) values of 5 and 230 nM, respectively). The solution structure for μ-KIIIA based on nuclear magnetic resonance data was recalculated with the [C1-C15,C2-C9,C4-C16] disulfide pattern; its structure was very similar to the μ-KIIIA structure calculated with the incorrect [C1-C9,C2-C15,C4-C16] disulfide pattern, with an α-helix spanning residues 7-12. In addition, the major folding isomers of μ-KIIIB, an N-terminally extended isoform of μ-KIIIA identified from its cDNA sequence, were isolated. These folding products had the same disulfide connectivities as μ-KIIIA, and both blocked Na(V)1.2 (K(d) values of 470 and 26 nM, respectively). Our results establish that the preferred disulfide pattern of synthetic μ-KIIIA and μ-KIIIB folded in vitro is 1-5/2-4/3-6 but that other disulfide isomers are also potent sodium channel blockers. These findings raise questions about the disulfide pattern(s) of μ-KIIIA in the venom of Conus kinoshitai; indeed, the presence of multiple disulfide isomers in the venom could provide a means of further expanding the snail's repertoire of active peptides.

Mass spectrometry in India

This review emphasizes the mass spectrometry research being performed at academic and established research institutions in India. It consists of three main parts covering the work done in organic, atomic and biological mass spectrometry. The review reveals that the use of mass spectrometry techniques started in the middle of the 20th century and was applied to research in the fields of organic, nuclear, geographical and atomic chemistry. Later, with the advent of soft and atmospheric ionization techniques it has been applied to pharmaceutical and biological research. In due course, several research centers with advanced mass spectrometry facilities have been established for specific areas of research such as gas-phase ion chemistry, ion-molecule reactions, proscribed chemicals, pesticide residues, pharmacokinetics, protein/peptide chemistry, nuclear chemistry, geochronological studies, archeology, petroleum industry, proteomics, lipidomics and metabolomics. Day-by-day the mass spectrometry centers/facilities in India have attracted young students for their doctoral research and other advanced research applications.

PADLOC: a powerful tool to assign disulfide bond connectivities in peptides and proteins by NMR spectroscopy

The determination of the disulfide bond connectivity in a peptide or protein represents a significant challenge. It is notoriously difficult to use NMR spectroscopy to assign disulfide connectivities because NMR spectra lack direct evidence for disulfide bonds. These bonds are typically inferred from three-dimensional structure calculations, which can result in ambiguous disulfide assignment. Here, we present a new NMR based methodology, in which the disulfide connectivity is obtained by applying Bayesian rules of inference to the local topology of cysteine residues. We illustrate how this approach successfully predicts the disulfide connectivity in proteins for which crystal structures are available in the protein data bank (PDB). We also demonstrate how this methodology is used with experimental NMR data for peptides with complex disulfide topologies, including hepcidin, Kalata-B1, and μ-Conotoxin KIIIA. In the case of μ-Conotoxin KIIIA, the PADLOC connectivity (1-15,2-9,4-16) differs from previously published results; additional evidence is presented demonstrating unequivocally that this newly proposed connectivity is correct.