Structure and folding of disulfide-rich miniproteins: insights from molecular dynamics simulations and MM-PBSA free energy calculations

Graph-Directed Approach for Downselecting Toxins for Experimental Structure Determination

Conotoxins are short, cysteine-rich peptides of great interest as novel therapeutic leads and of great concern as lethal biological agents due to their high affinity and specificity for various receptors involved in neuromuscular transmission. Currently, of the approximately 6000 known conotoxin sequences, only about 3% have associated structural characterization, which leads to a bottleneck in rapid high-throughput screening (HTS) for identification of potential leads or threats. In this work, we combine a graph-based approach with homology modeling to expand the library of conotoxin structures and to identify those conotoxin sequences that are of the greatest value for experimental structural characterization. The latter would allow for the rapid expansion of the known structural space for generating high quality template-based models. Our approach generalizes to other evolutionarily-related, short, cysteine-rich venoms of interest. Overall, we present and validate an approach for venom structure modeling and experimental guidance and employ it to produce a 290%-larger library of approximate conotoxin structures for HTS. We also provide a set of ranked conotoxin sequences for experimental structure determination to further expand this library.

Cyclotides: From Structure to Function

This Review explores the class of plant-derived macrocyclic peptides called cyclotides. We include an account of their discovery, characterization, and distribution in the plant kingdom as well as a detailed analysis of their sequences and structures, biosynthesis and chemical synthesis, biological functions, and applications. These macrocyclic peptides are around 30 amino acids in size and are characterized by their head-to-tail cyclic backbone and cystine knot motif, which render them to be exceptionally stable, with resistance to thermal or enzymatic degradation. Routes to their chemical synthesis have been developed over the past two decades, and this capability has facilitated a wide range of mutagenesis and structure-activity relationship studies. In turn, these studies have both led to an increased understanding of their mechanisms of action as well as facilitated a range of applications in agriculture and medicine, as ecofriendly crop protection agents, and as drug leads or scaffolds for pharmaceutical design. Our overall objective in this Review is to provide readers with a comprehensive overview of cyclotides that we hope will stimulate further work on this fascinating family of peptides.

Interaction of rat α9α10 nicotinic acetylcholine receptor with α-conotoxin RgIA and Vc1.1: Insights from docking, molecular dynamics and binding free energy contributions

The α9α10 nicotinic acetylcholine receptor (nAChR) is an effective therapeutic target for neuropathic pain. α-Conotoxin RgIA and Vc1.1 are two well-known peptides blocking α9α10 nAChR potently and selectively, which have been extensively investigated as drug candidates. Several key residues were established in previous experimental research. However, the mechanism of the specific interaction still needs to be elucidated in more detail. In this work, we explored the interaction mechanism between RgIA/Vc1.1 and rat α9α10 nAChR using docking and molecular dynamics (MD) simulations. Energy and network analysis programs were used to reveal key residues responsible for their interaction. Our results indicated that the most critical residues were in accord with previous studies. Importantly, several novel residues, including Tyr95, Trp151 in α9 (+)α10 (-) interface as well as Tyr196, Arg59in α10 (+)α9 (-) interface, were found in our models. Furthermore, we analyzed noncovalent interaction energies between RgIA/Vc1.1 and rat α9α10 nAChR. The results showed that three negatively charged residues (Glu197 in α10 subunit, Asp168 in α9 subunit and Asp205 in α10 subunit) were involved in the interaction between RgIA and rat α9α10 nAChR. In contrast, the interaction between Vc1.1 and rat α9α10 nAChR was mediated by the positively charged residues Arg59, Arg81 in α9 (-) subunit. These findings provided further insights into the molecular mechanisms of interaction between RgIA and Vc1.1 and rat α9α10 nAChR.

Cytosolic expression, solution structures, and molecular dynamics simulation of genetically encodable disulfide-rich de novo designed peptides

Disulfide-rich peptides represent an important protein family with broad pharmacological potential. Recent advances in computational methods have made it possible to design new peptides which adopt a stable conformation de novo. Here, we describe a system to produce disulfide-rich de novo peptides using Escherichia coli as the expression host. The advantage of this system is that it enables production of uniformly 13 C- and 15 N-labeled peptides for solution nuclear magnetic resonance (NMR) studies. This expression system was used to isotopically label two previously reported de novo designed peptides, and to determine their solution structures using NMR. The ensemble of NMR structures calculated for both peptides agreed well with the design models, further confirming the accuracy of the design protocol. Collection of NMR data on the peptides under reducing conditions revealed a dependency on disulfide bonds to maintain stability. Furthermore, we performed long-time molecular dynamics (MD) simulations with tempering to assess the stability of two families of de novo designed peptides. Initial designs which exhibited a stable structure during simulations were more likely to adopt a stable structure in vitro, but attempts to utilize this method to redesign unstable peptides to fold into a stable state were unsuccessful. Further work is therefore needed to assess the utility of MD simulation techniques for de novo protein design.

Role of individual disulfide bridges in the conformation and activity of spinoxin (α-KTx6.13), a potassium channel toxin from Heterometrus spinifer scorpion venom

Spinoxin (SPX; α-KTx6.13), isolated from venom of the scorpion Heterometrus spinifer, is a K⁺ channel-specific peptide toxin (KTx), which adopts a cysteine-stabilized α/β scaffold that is cross-linked by four disulfide bridges (Cys1-Cys5, Cys2-Cys6, Cys3-Cys7, and Cys4-Cys8). To investigate the role of the individual disulfide bonds in the structure-activity relationship of SPX, we synthesized four SPX analogs in which each pair of cysteine residues was replaced by alanine residues. The analysis of circular dichroism spectra and inhibitory activity against Kv1.3 channels showed that the SPX analogs lacking any of three specific disulfide bonds (Cys1-Cys5, Cys2-Cys6, and Cys3-Cys7) were unable to form the native secondary structure and completely lost inhibitory activities. Thus, we conclude that Cys1-Cys5, Cys2-Cys6, and Cys3-Cys7 are required for the inhibition of the Kv1.3 channel by SPX. In contrast, the analog lacking Cys4-Cys8 retained both native secondary structure and inhibitory activity. Interestingly, one of the isomers of the analog lacking Cys1-Cys5 also showed inhibitory activities, although its inhibition was ∼18-fold weaker than native SPX. This isomer had an atypical disulfide bond pairing (Cys3-Cys4 and Cys7-Cys8) that corresponds to that of maurotoxin (MTX), another α-KTx6 family member. These results indicate that the Cys1-Cys5 and Cys2-Cys6 bonds are important for restricting the toxin from forming an atypical (MTX-type) disulfide bond pairing among the remaining four cysteine residues (Cys3, Cys4, Cys7, and Cys8) in native SPX.

Probing pattern and dynamics of disulfide bridges using synthesis and NMR of an ion channel blocker peptide toxin with multiple diselenide bonds

Anuroctoxin (AnTx), a 35-amino-acid scorpion toxin containing four disulfide bridges, is a high affinity blocker of the voltage-gated potassium channel Kv1.3, but also blocks Kv1.2. To improve potential therapeutic use of the toxin, we have designed a double substituted analog, [N17A/F32T]-AnTx, which showed comparable Kv1.3 affinity to the wild-type peptide, but also a 2500-fold increase in the selectivity for Kv1.3 over Kv1.2. In the present study we have achieved the chemical synthesis of a Sec-analog in which all cysteine (Cys) residues have been replaced by selenocysteine (Sec) forming four diselenide bonds. To the best of our knowledge this is the first time to replace, by chemical synthesis, all disulfide bonds with isosteric diselenides in a peptide/protein. Gratifyingly, the key pharmacological properties of the Sec-[N17A/F32T]-AnTx are retained since the peptide is functionally active. We also propose here a combined experimental and theoretical approach including NOE- and 77Se-based NMR supplemented by MD simulations for conformational and dynamic characterization of the Sec-[N17A/F32T]-AnTx. Using this combined approach allowed us to attain unequivocal assignment of all four diselenide bonds and supplemental MD simulations allowed characterization of the conformational dynamics around each disulfide/diselenide bridge.

The knottin-like Blufensin family regulates genes involved in nuclear import and the secretory pathway in barley-powdery mildew interactions

Plants have evolved complex regulatory mechanisms to control a multi-layered defense response to microbial attack. Both temporal and spatial gene expression are tightly regulated in response to pathogen ingress, modulating both positive and negative control of defense. BLUFENSINs, small knottin-like peptides in barley, wheat, and rice, are highly induced by attack from fungal pathogens, in particular, the obligate biotrophic fungus, Blumeria graminis f. sp. hordei (Bgh), causal agent of barley powdery mildew. Previous research indicated that Blufensin1 (Bln1) functions as a negative regulator of basal defense mechanisms. In the current report, we show that BLN1 and BLN2 can both be secreted to the apoplast and Barley stripe mosaic virus (BSMV)-mediated overexpression of Bln2 increases susceptibility of barley to Bgh. Bimolecular fluorescence complementation (BiFC) assays signify that BLN1 and BLN2 can interact with each other, and with calmodulin. We then used BSMV-induced gene silencing to knock down Bln1, followed by Barley1 GeneChip transcriptome analysis, to identify additional host genes influenced by Bln1. Analysis of differential expression revealed a gene set enriched for those encoding proteins annotated to nuclear import and the secretory pathway, particularly Importin α1-b and Sec61 γ subunits. Further functional analysis of these two affected genes showed that when silenced, they also reduced susceptibility to Bgh. Taken together, we postulate that Bln1 is co-opted by Bgh to facilitate transport of disease-related host proteins or effectors, influencing the establishment of Bgh compatibility on its barley host.

Evolution of a horizontally acquired legume gene, albumin 1, in the parasitic plant Phelipanche aegyptiaca and related species

Background

Parasitic plants, represented by several thousand species of angiosperms, use modified structures known as haustoria to tap into photosynthetic host plants and extract nutrients and water. As a result of their direct plant-plant connections with their host plant, parasitic plants have special opportunities for horizontal gene transfer, the nonsexual transmission of genetic material across species boundaries. There is increasing evidence that parasitic plants have served as recipients and donors of horizontal gene transfer (HGT), but the long-term impacts of eukaryotic HGT in parasitic plants are largely unknown.

Results

Here we show that a gene encoding albumin 1 KNOTTIN-like protein, closely related to the albumin 1 genes only known from papilionoid legumes, where they serve dual roles as food storage and insect toxin, was found in Phelipanche aegyptiaca and related parasitic species of family Orobanchaceae, and was likely acquired by a Phelipanche ancestor via HGT from a legume host based on phylogenetic analyses. The KNOTTINs are well known for their unique “disulfide through disulfide knot” structure and have been extensively studied in various contexts, including drug design. Genomic sequences from nine related parasite species were obtained, and 3D protein structure simulation tests and evolutionary constraint analyses were performed. The parasite gene we identified here retains the intron structure, six highly conserved cysteine residues necessary to form a KNOTTIN protein, and displays levels of purifying selection like those seen in legumes. The albumin 1 xenogene has evolved through >150 speciation events over ca. 16 million years, forming a small family of differentially expressed genes that may confer novel functions in the parasites. Moreover, further data show that a distantly related parasitic plant, Cuscuta, obtained two copies of albumin 1 KNOTTIN-like genes from legumes through a separate HGT event, suggesting that legume KNOTTIN structures have been repeatedly co-opted by parasitic plants.

Conclusions

The HGT-derived albumins in Phelipanche represent a novel example of how plants can acquire genes from other plants via HGT that then go on to duplicate, evolve, and retain the specialized features required to perform a unique host-derived function.

Engineered pathways for correct disulfide bond oxidation

Correct formation of disulfide bonds is critical for protein folding. We find that cells lacking protein disulfide isomerases (PDIs) can use alternative mechanisms for correct disulfide bond formation. By linking correct disulfide bond formation to antibiotic resistance, we selected mutants that catalyze correct disulfide formation in the absence of DsbC, Escherichia coli's PDI. Most of our mutants massively overproduce the disulfide oxidase DsbA and change its redox status. They enhance DsbA's ability to directly form the correct disulfides by increasing the level of mixed disulfides between DsbA and substrate proteins. One mutant operates via a different mechanism; it contains mutations in DsbB and CpxR that alter the redox environment of the periplasm and increases the level of the chaperone/protease DegP, allowing DsbA to gain disulfide isomerase ability in vivo. Thus, given the proper expression level, redox status, and chaperone assistance, the oxidase DsbA can readily function in vivo to catalyze the folding of proteins with complex disulfide bond connectivities. Our selection reveals versatile strategies for correct disulfide formation in vivo. Remarkably, our evolution of new pathways for correct disulfide bond formation in E. coli mimics eukaryotic PDI, a highly abundant partially reduced protein with chaperone activity.

Molecular jamming--the cystine slipknot mechanical clamp in all-atom simulations

A recent survey of 17 134 proteins has identified a new class of proteins which are expected to yield stretching induced force peaks in the range of 1 nN. Such high force peaks should be due to forcing of a slip-loop through a cystine ring, i.e., by generating a cystine slipknot. The survey has been performed in a simple coarse grained model. Here, we perform all-atom steered molecular dynamics simulations on 15 cystine knot proteins and determine their resistance to stretching. In agreement with previous studies within a coarse grained structure based model, the level of resistance is found to be substantially higher than in proteins in which the mechanical clamp operates through shear. The large stretching forces arise through formation of the cystine slipknot mechanical clamp and the resulting steric jamming. We elucidate the workings of such a clamp in an atomic detail. We also study the behavior of five top strength proteins with the shear-based mechanostability in which no jamming is involved. We show that in the atomic model, the jamming state is relieved by moving one amino acid at a time and there is a choice in the selection of the amino acid that advances the first. In contrast, the coarse grained model also allows for a simultaneous passage of two amino acids.

Design and therapeutic applications of cyclotides

Cyclotides are plant-derived peptides with a cyclic backbone and knotted topology of disulfide bonds. Their extreme stability and natural sequence variation has led to the suggestion that they might be useful as scaffolds to stabilize bioactive sequences. Recent studies have shown that anti-angiogenic activity and protease inhibitory activity against a foot and mouth disease protease can be grafted onto the cyclotide framework. There has also been significant progress made in determining the mechanism of cyclization of cyclotides and in producing cyclotides using bacterial expression and plant cell culture. There is a wide range of disease states that can be targeted using the cyclotide framework and the advances that have been made in the production of cyclotides will facilitate their development as pharmaceutical templates.

Knottin cyclization: impact on structure and dynamics

Background

Present in various species, the knottins (also referred to as inhibitor cystine knots) constitute a group of extremely stable miniproteins with a plethora of biological activities. Owing to their small size and their high stability, knottins are considered as excellent leads or scaffolds in drug design. Two knottin families contain macrocyclic compounds, namely the cyclotides and the squash inhibitors. The cyclotide family nearly exclusively contains head-to-tail cyclized members. On the other hand, the squash family predominantly contains linear members. Head-to-tail cyclization is intuitively expected to improve bioactivities by increasing stability and lowering flexibility as well as sensitivity to proteolytic attack.

Results

In this paper, we report data on solution structure, thermal stability, and flexibility as inferred from NMR experiments and molecular dynamics simulations of a linear squash inhibitor EETI-II, a circular squash inhibitor MCoTI-II, and a linear analog lin-MCoTI. Strikingly, the head-to-tail linker in cyclic MCoTI-II is by far the most flexible region of all three compounds. Moreover, we show that cyclic and linear squash inhibitors do not display large differences in structure or flexibility in standard conditions, raising the question as to why few squash inhibitors have evolved into cyclic compounds. The simulations revealed however that the cyclization increases resistance to high temperatures by limiting structure unfolding.

Conclusion

In this work, we show that, in contrast to what could have been intuitively expected, cyclization of squash inhibitors does not provide clear stability or flexibility modification. Overall, our results suggest that, for squash inhibitors in standard conditions, the circularization impact might come from incorporation of an additional loop sequence, that can contribute to the miniprotein specificity and affinity, rather than from an increase in conformational rigidity or protein stability. Unfolding simulations showed however that cyclization is a stabilizing factor in strongly denaturing conditions. This information should be useful if one wants to use the squash inhibitor scaffold in drug design.