Search

Marine organisms produce a diverse range of toxins and bioactive peptides to support predation, competition, and defense. The peptide repertoires of stony corals (order Scleractinia) remain relatively understudied despite the presence of tentacles used for predation and defense that are likely to contain a range of bioactive compounds. Here, we show that a tentacle extract from the mushroom coral, Heliofungia actiniformis, contains numerous peptides with a range of molecular weights analogous to venom profiles from species such as cone snails. Using NMR spectroscopy and mass spectrometry we characterized a 12-residue peptide (Hact-1) with a new sequence (GCHYTPFGLICF) and well-defined β-hairpin structure stabilized by a single disulfide bond. The sequence is encoded within the genome of the coral and expressed in the polyp body tissue. The structure present is common among toxins and venom peptides, but Hact-1 does not show activity against select examples of Gram-positive and Gram-negative bacteria or a range of ion channels, common properties of such peptides. Instead, it appears to have a limited effect on human peripheral blood mononuclear cells, but the ecological function of the peptide remains unknown. The discovery of this peptide from H. actiniformis is likely to be the first of many from this and related species.

Delayed-onset muscle soreness (DOMS) is a very common condition in athletes and individuals not accustomed to physical activity that occurs after moderate/high-intensity exercise sessions. The activation of microglial Toll-like receptor 4 (TLR4) in the spinal cord has been described to be important for the induction and maintenance of persistent pain. Based on that, we hypothesize that 70 kilodalton heat-shock protein (Hsp70), a mediator released by exercise, could activate microglial TLR4 in the spinal cord, releasing proinflammatory cytokines and contributing to the start of DOMS. In fact, we found that the knockout of TLR4, myeloid differentiation primary response 88 (MyD88), interleukin-6 (IL-6), or both tumor necrosis factor-α (TNF-α) receptor 1 and TNF-α receptor 2 in mice prevented the development of DOMS following acute aerobic exercise in contrast to the findings in male C57BL/6 wild-type mice. Furthermore, DOMS in exercised wild-type mice was also prevented after pre-treatment with microglia inhibitor, TLR4 antagonist, and anti-Hsp70 antibody. During exercise-induced DOMS, Hsp70 mRNA, TLR4 mRNA, and protein levels, as well as Iba-1 (a microglial marker), IL-6, and TNF-α protein levels, were increased in the muscle and/or spinal cord. Together, these findings suggest that Hsp70 released during exercise-induced DOMS activates the microglial TLR4/IL-6/TNF-α pathway in the spinal cord. Thus, the blockade of TLR4 activation may be a new strategy to prevent the development of DOMS before intense exercise.

The majority of current drugs against diseases, such as cancer, can bind to one or more sites in a protein and inhibit its activity. There are, however, well known limits on the number of druggable proteins and complementing current drugs with compounds that could selectively target DNA or RNA would greatly enhance therapeutic progress and options. We are focusing on the design of sequence-specific DNA minor groove binders that, for example, target the promoter sites of transcription factors involved in a disease. We have started with AT specific minor groove binders that are known to enter human cells and have entered clinical trials. To broaden the sequence-specific recognition of these compounds, we have identified several modules that have H-bond acceptors that strongly and specifically recognize G•C base pairs. A lead module is a thiophene-N-alkyl-benzimidazole σ-hole based system with terminal phenyl-amidines where the optimum compounds have excellent affinity and selectivity for a G•C base pair in the minor groove. Efforts are now focused on optimizing this module. We have previously optimized the alkyl group. In the work described here we are evaluating modifications to the compound aromatic system with the goal of improving GC selectivity and affinity as with the N-alkyl modifications. The lead compounds from these studies retain the thiophene-N-alkyl-BI module but have halogen substituents adjacent to an amidine group on the terminal phenyl-amidine. Other improved compounds in this set have modified amidines and conversion of the amidine to an imidazoline, for example, resulted in a strong binding compound with good specificity.

[Image: see text] Interactions are reported of three representative naphthalenediimide derivatives with three quadruplex targets, from the promoter region of the telomerase (hTERT) gene, a human telomeric DNA quadruplex, and a telomeric RNA quadruplex (TERRA). Thermal melting studies showed that these compounds strongly stabilize the quadruplexes, with weak stabilization of a duplex DNA. Binding studies by surface plasmon resonance and fluorescence spectroscopy found that the compounds bind to the quadruplexes with nanomolar equilibrium dissociation constants. Plausible topologies for the quadruplex complexes were deduced from CD spectra, which together with the surface plasmon resonance data indicate that the quadruplexes with parallel quadruplex folds are preferred by two compounds, which was confirmed by qualitative molecular modeling.

Hairpin polyamides are synthetic small molecules that bind DNA minor groove sequence-selectively and, in many sequences, induce widening of the minor groove and compression of the major groove. The structural distortion of DNA caused by polyamides has enhanced our understanding of the regulation of DNA-binding proteins via polyamides. Polyamides have DNA binding affinities that are comparable to those proteins, therefore, can potentially be used as therapeutic agents to treat diseases caused by aberrant gene expression. In fact, many diseases are characterized by over- or under-expressed genes. PU.1 is a transcription factor that regulates many immune system genes. Aberrant expression of PU.1 has been associated with the development of acute myeloid leukemia (AML). We have, therefore, designed and synthesized ten hairpin polyamides to investigate their capacity in controlling the PU.1-DNA interaction. Our results showed that nine of the polyamides disrupt PU.1-DNA binding and the inhibition capacity strongly correlates with binding affinity. One molecule, FH1024, was observed forming a FH1024-PU.1-DNA ternary complex instead of inhibiting PU.1-DNA binding. This is the first report of a small molecule that is potentially a weak agonist that recruits PU.1 to DNA. This finding sheds light on the design of polyamides that exhibit novel regulatory mechanisms on protein-DNA binding.

The understanding of sequence-specific DNA minor groove interactions has recently made major steps forward and as a result, the goal of development of compounds that target the minor groove is an active research area. In an effort to develop biologically active minor groove agents, we are preparing and exploring the DNA interactions of diverse diamidine derivatives with a 5’-GAATTC-3’ binding site using a powerful array of methods including, biosensor-SPR methods, and X-ray crystallography. The benzimidazole-thiophene module provides an excellent minor groove recognition component. A central thiophene in a benzimidazole-thiophene-phenyl aromatic system provides essentially optimum curvature for matching the shape of the minor groove. Comparison of that structure to one with the benzimidazole replaced with an indole shows that the two structures are very similar, but have some interesting and important differences in electrostatic potential maps, the DNA minor groove binding structure based on x-ray crystallographic analysis, and inhibition of the major groove binding PU.1 transcription factor complex. The binding K(D) for both compounds is under 10 nM and both form amidine H-bonds to DNA bases. They both have bifurcated H-bonds from the benzimidazole or indole groups to bases at the center of the -AATT- binding site. Analysis of the comparative results provides an excellent understanding of how thiophene compounds recognize the minor groove and can act as transcription factor inhibitors.

Cyclic peptides are widespread throughout the plant kingdom, and display diverse sequences, structures and bioactivities. The potential applications attributed to these peptides and their unusual biosynthesis has captivated the attention of researchers for many years. Several gene sequences for plant cyclic peptides have been discovered over the last two decades but it is only recently that we are beginning to understand the intricacies associated with their biosynthesis. Recent studies have focussed on three main classes of plant derived cyclic peptides, namely orbitides, SFTI related peptides and cyclotides. In this mini-review, we discuss the expansion of the known sequence and structural diversity in these families, insights into the enzymes involved in the biosynthesis, the exciting applications which includes a cyclotide currently in clinical trials for the treatment of multiple sclerosis, and new production methods that are being developed to realise the potential of plant cyclic peptides as pharmaceutical or agricultural agents.

[Image: see text] Targeting of genomic quadruplexes is an approach to treating complex human cancers. We describe a series of tetra-substituted naphthalene diimide (ND) derivatives with a phenyl substituent directly attached to the ND core. The lead compound (SOP1812) has 10 times superior cellular and in vivo activity compared with previous ND compounds and nanomolar binding to human quadruplexes. The pharmacological properties of SOP1812 indicate good bioavailability, which is consistent with the in vivo activity in xenograft and genetic models for pancreatic cancer. Transcriptome analysis shows that it down-regulates several cancer gene pathways, including Wnt/β-catenin signaling.

Cone snails produce a fast-acting and often paralyzing venom, largely dominated by disulfide-rich conotoxins targeting ion channels. Although disulfide-poor conopeptides are usually minor components of cone snail venoms, their ability to target key membrane receptors such as GPCRs make them highly valuable as drug lead compounds. From the venom gland transcriptome of Conus miliaris, we report here on the discovery and characterization of two conopressins, which are nonapeptide ligands of the vasopressin/oxytocin receptor family. These novel sequence variants show unusual features, including a charge inversion at the critical position 8, with an aspartate instead of a highly conserved lysine or arginine residue. Both the amidated and acid C-terminal analogues were synthesized, followed by pharmacological characterization on human and zebrafish receptors and structural investigation by NMR. Whereas conopressin-M1 showed weak and only partial agonist activity at hV1bR (amidated form only) and ZFV1a1R (both amidated and acid form), both conopressin-M2 analogues acted as full agonists at the ZFV2 receptor with low micromolar affinity. Together with the NMR structures of amidated conopressins-M1, -M2 and -G, this study provides novel structure-activity relationship information that may help in the design of more selective ligands.

DB2277, a heterocyclic diamidine, is a successful design for mixed base pair (bp) DNA sequence recognition. The compound has a central aza-benzimidazole group that forms two H-bonds with a GC bp that has flanking AT bps. The nuclear magnetic resonance structure of the DB2277-DNA complex with an AAGATA recognition site sequence was determined, and here we report extended molecular dynamics (MD) simulations of the structure. DB2277 has two terminal phenyl-amidine groups, one of which is directly linked to the DB2277 heterocyclic core and the other through a flexible -OCH2- group. The flexibly linked phenyl is too far from the minor groove floor to make direct H-bonds but is linked to an AT bp through water-mediated H-bonds. The flexibly linked phenyl-amidine with water-mediated H-bonds to the bases at the floor of the minor groove suggested that it might rotate in time spans accessible in MD. To test this idea, we conducted multimicrosecond MD simulations to determine if these phenyl rotations could be observed for a bound compound. In a 3 μs simulation, highly dynamic torsional motions were observed for the -OCH2-linked phenyl but not for the other phenyl. The dynamics periodically reached a level to allow 180° rotation of the phenyl while it was still bound in the minor groove. This is the first observation of rotation of a phenyl bound to DNA, and the results provide mechanistic details about how a rotation can occur as well as how mixed bp recognition can occur for monomer compounds bound to the minor groove.

To extend the value of biosensor-SPR in the characterization of DNA recognition by nucleoproteins, we report a comparative analysis of DNA-facilitated target search by two ETS-family transcription factors: Elk1 and ETV6. ETS domains represent an attractive system for developing biosensor-based techniques due to a broad range of physicochemical properties encoded within a highly conserved DNA-binding motif. Building on a biosensor approach in which the protein is quantitatively sequestered and presented to immobilized cognate DNA as nonspecific complexes, we assessed the impact of intrinsic cognate and nonspecific affinities on long-range (intersegmental) target search. The equilibrium constants of DNA-facilitated binding were sensitive to the intrinsic binding properties of the proteins such that their relative specificity for cognate DNA were reinforced when binding occurred by transfer vs. without nonspecific DNA. Direct measurement of association and dissociation kinetics revealed ionic features of the activated complex that evidenced DNA-facilitated dissociation, even though Elk1 and ETV6 harbor only a single DNA-binding surface. At salt concentrations that masked the effects of nonspecific pre-binding at equilibrium, the dissociation kinetics of cognate binding were nevertheless distinct from conditions under which nonspecific DNA was absent. These results further strengthen the significance of long-range DNA-facilitated translocation in the physiologic environment.

Inflammatory bowel diseases (IBDs) are a set of complex and debilitating diseases for which there is no satisfactory treatment. Recent studies have shown that small peptides show promise for reducing inflammation in models of IBD. However, these small peptides are likely to be unstable and rapidly cleared from the circulation, and therefore, if not modified for better stability, represent non-viable drug leads. We hypothesized that improving the stability of these peptides by grafting them into a stable cyclic peptide scaffold may enhance their therapeutic potential. Using this approach, we have designed a novel cyclic peptide that comprises a small bioactive peptide from the annexin A1 protein grafted into a sunflower trypsin inhibitor cyclic scaffold. We used native chemical ligation to synthesize the grafted cyclic peptide. This engineered cyclic peptide maintained the overall fold of the naturally occurring cyclic peptide, was more effective at reducing inflammation in a mouse model of acute colitis than the bioactive peptide alone, and showed enhanced stability in human serum. Our findings suggest that the use of cyclic peptides as structural backbones offers a promising approach for the treatment of IBD and potentially other chronic inflammatory conditions.

Granulins are a family of protein growth factors that are involved in cell proliferation. An orthologue of granulin from the human parasitic liver fluke Opisthorchis viverrini, known as Ov-GRN-1, induces angiogenesis and accelerates wound repair. Recombinant Ov-GRN-1 production is complex and poses an obstacle for clinical development. To identify the bioactive region(s) of Ov-GRN-1, four truncated N-terminal analogues were synthesized and characterized structurally using NMR spectroscopy. Peptides that contained only two native disulfide bonds lack the characteristic granulin β-hairpin structure. Remarkably, the introduction of a non-native disulfide bond was critical for formation of β-hairpin structure. Despite this structural difference, both two and three disulfide-bonded peptides drove proliferation of a human cholangiocyte cell line and demonstrated potent wound healing in mice. Peptides derived from Ov-GRN-1 are leads for novel wound healing therapeutics, as they are likely less immunogenic than the full-length protein and more convenient to produce.

Minor-groove binding hairpin polyamides (PAs) bind specific DNA sequences. Synthetic modifications can improve PA-DNA binding affinity and include flexible modules, such as β-alanine (β) motifs to replace pyrroles (Py), and increasing compound charge using N-terminal cationic substituents. To better understand the variations in kinetics and affinities caused by these modifications on PA-DNA interactions, a comprehensive set of PAs with different numbers and positions of β and different types of N-cationic groups was systematically designed and synthesized to bind their cognate sequence, the λB motif. The λB motif is also a strong binding promoter site of the major groove targeting transcription factor PU.1. The PA binding affinities and kinetics were evaluated using a spectrum of powerful biophysical methods: thermal melting, biosensor surface plasmon resonance and circular dichroism. The results show that β inserts affect PA-DNA interactions in a number and position dependent manner. Specifically, a β replacement between two imidazole heterocycles (ImβIm) generally strengthens binding. In addition, N-terminal cationic groups can accelerate the association between PA and DNA, but the bulky size of TMG can cause steric hindrance and unfavourable repulsive electrostatic interactions in some PAs. The future design of stronger binding PA requires careful combination of βs and cationic substituents.

Specific targeting of protein–nucleic acid interactions is an area of current interest, for example, in the regulation of gene-expression. Most transcription factor proteins bind in the DNA major groove; however, we are interested in an approach using small molecules to target the minor groove to control expression by an allosteric mechanism. In an effort to broaden sequence recognition of DNA-targeted-small-molecules to include both A·T and G·C base pairs, we recently discovered that the heterocyclic diamidine, DB2277, forms a strong monomer complex with a DNA sequence containing 5΄-AAAGTTT-3΄. Competition mass spectrometry and surface plasmon resonance identified new monomer complexes, as well as unexpected binding of two DB2277 with certain sequences. Inherent microstructural differences within the experimental DNAs were identified through computational analyses to understand the molecular basis for recognition. These findings emphasize the critical nature of the DNA minor groove microstructure for sequence-specific recognition and offer new avenues to design synthetic small molecules for effective regulation of gene-expression.

The ends of nucleic acids oligomers alter the statistics of interior nonspecific ligand binding and act as binding sites with altered properties. While the former aspect of “end effects” has received much theoretical attention in the literature, the physical nature of end-binding, and hence its potential impact on a wide range of studies with oligomers, remains poorly known. Here, we report for the first time end-binding to DNA using a model helix-turn-helix motif, the DNA-binding domain of ETV6, as a function of DNA sequence length. Spectral analysis of ETV6 intrinsic tryptophan fluorescence by singular value decomposition showed that end-binding to nonspecific fragments was negligible at >0.2 kbp and accumulated to 8% of total binding to 23-bp oligomers. The affinity for end-binding was insensitive to salt but tracked the affinity of interior binding, suggesting translocation from interior sites rather than free solution as its mechanism. As the presence of a cognate site in the 23-bp oligomer suppressed end-binding, neglect of end-binding to the short cognate DNA does not introduce significant error. However, the same applies to nonspecific DNA only if longer fragments (>0.2 kbp) are used.

Sequence-specific binding to DNA is crucial for targeting transcription factor-DNA complexes to modulate gene expression. The heterocyclic diamidine, DB2277, specifically recognizes a single G•C base pair in the minor groove of mixed base pair sequences of the type AAAGTTT. NMR spectroscopy reveals the presence of major and minor species of the bound compound. To understand the principles that determine the binding affinity and orientation in mixed sequences of DNA, over thirty DNA hairpin substrates were examined by NMR and thermal melting. The NMR exchange dynamics between major and minor species shows that the exchange is much faster than compound dissociation determined from biosensor–surface plasmon resonance. Extensive modifications of DNA sequences resulted in a unique DNA sequence with binding site AAGATA that binds DB2277 in a single orientation. A molecular docking result agrees with the model representing rapid flipping of DB2277 between major and minor species. Imino spectral analysis of a 15N-labeled central G clearly shows the crucial role of the exocyclic amino group of G in sequence-specific recognition. Our results suggest that this approach can be expanded to additional modules for recognition of more sequence-specific DNA complexes. This approach provides substantial information about the sequence-specific, highly efficient, dynamic nature of minor groove binding agents.

Granulins are a family of unique protein growth factors which are found in a range of species and have several bioactivities that include cell proliferation and wound healing. They typically contain six disulfide bonds, but the sequences, structures and bioactivities vary significantly. We have previously shown that an N-terminally truncated version of a granulin from the human liver fluke, Opisthorchis viverrini, can fold independently into a “mini-granulin” structure and has potent wound healing properties in vivo. The incorporation of a non-native third disulfide bond, with respect to the full-length granulin module, was critical for the formation of regular secondary structure in the liver fluke derived peptide. By contrast, this third disulfide bond is not required for a carp granulin-1 truncated peptide to fold independently. This distinction led us to explore granulins from the zebrafish model organism. Here we show that the mini-granulin fold occurs in a naturally occurring paragranulin (half-domain) from zebrafish, and is also present in a truncated form of a full-length zebrafish granulin, suggesting this structure might be a common property in either naturally occurring or engineered N-terminally truncated granulins and the carp granulin-1 folding is an anomaly. The in vitro folding yield is significantly higher in the naturally occurring paragranulin, but only the truncated zebrafish granulin peptide promoted the proliferation of fibroblasts consistent with a growth factor function, and therefore the function of the paragranulin remains unknown. These findings provide insight into the folding and evolution of granulin domains and might be useful in the elucidation of the structural features important for bioactivity to aid the design of more potent and stable analogues for the development of novel wound healing agents.

Scorpion venoms are a rich source of bioactive molecules, but characterisation of toxin peptides affecting cytosolic Ca2+, central to cell signalling and cell death, is limited. We undertook a functional screening of the venom of the Australian scorpion Hormurus waigiensis to determine the breadth of Ca2+ mobilisation. A human embryonic kidney (HEK293) cell line stably expressing the genetically encoded Ca2+ reporter GCaMP5G and the rabbit type 1 ryanodine&nbsp, receptor (RyR1) was developed as a biosensor. Size-exclusion Fast Protein Liquid Chromatography separated the venom into 53 fractions, constituting 12 chromatographic peaks. Liquid chromatography mass spectroscopy identified 182 distinct molecules with 3 to 63 components per peak. The molecular weights varied from 258 Da&mdash, 13.6 kDa, with 53% under 1 kDa. The majority of the venom chromatographic peaks (tested as six venom pools) were found to reversibly modulate cell monolayer bioimpedance, detected using the xCELLigence platform (ACEA Biosciences). Confocal Ca2+ imaging showed 9/14 peak samples, with molecules spanning the molecular size range, increased cytosolic Ca2+ mobilization. H. waigiensis venom Ca2+ activity was correlated with changes in bio-impedance, reflecting multi-modal toxin actions on cell physiology across the venom proteome.

Cyanine compounds have previously shown excellent in vitro and promising in vivo antileishmanial efficacy, but the potential toxicity of these agents is a concern. A series of 22 analogs of thiazole orange ((Z)-1-methyl-4-((3-methylbenzo[d]thiazol-2(3H)-ylidene)methyl)quinolin-1-ium salt), a commercial cyanine dye with antileishmanial activity, were synthesized in an effort to increase the selectivity of such compounds while maintaining efficacy. Cyanines possessing substitutions on the quinolinium ring system displayed potency against Leishmania donovani axenic amastigotes that differed little from the parent compound (IC50 12–42 nM), while ring disjunction analogs were both less potent and less toxic. Changes in DNA melting temperature were modest when synthetic oligonucleotides were incubated with selected analogs (ΔTm ≤ 5 °C), with ring disjunction analogs showing the least effect on this parameter. Despite the high antileishmanial potency of the target compounds, their toxicity and relatively flat SAR suggests that further information regarding the target(s) of these molecules is needed to aid their development as antileishmanials.

AT specific heterocyclic cations that bind in the DNA duplex minor groove have had major successes as cell and nuclear stains and as therapeutic agents which can effectively enter human cells. Expanding the DNA sequence recognition capability of the minor groove compounds could also expand their therapeutic targets and have an impact in many areas, such as modulation of transcription factor biological activity. Success in the design of mixed sequence binding compounds has been achieved with N-methylbenzimidazole (N-MeBI) thiophenes which are preorganized to fit the shape of the DNA minor groove and H-bond to the –NH of G·C base pairs that projects into the minor groove. Initial compounds bind strongly to a single G·C base pair in an AT context with a specificity ratio of 50 (K(D) AT-GC/K(D) AT) or less and this is somewhat low for biological use. We felt that modifications of compound shape could be used to probe local DNA microstructure in target mixed base pair sequences of DNA and potentially improve the compound binding selectivity. Modifications were made by increasing the size of the benzimidazole N-substituent, for example, by using N-isobutyl instead of N-Me, and by changing the molecular twist by introducing substitutions at specific positions on the aromatic core of the compounds. In both cases, we have been able to achieve a dramatic increase in binding specificity, including no detectible binding to pure AT sequences, without a significant loss in affinity to mixed base pair target sequences.

Granulins are a family of growth factors involved in cell proliferation. The liver-fluke granulin, Ov-GRN-1, isolated from a carcinogenic liver fluke Opisthorchis viverrini, can significantly accelerate wound repair in vivo and in vitro. However, it is difficult to express Ov-GRN-1 in recombinant form at high yield, impeding its utility as a drug lead. Previously we reported that a truncated analogue (Ov-GRN12–35_3s) promotes healing of cutaneous wounds in mice. NMR analysis of this analogue indicates the presence of multiple conformations, most likely as a result of proline cis/trans isomerization. To further investigate whether the proline residues are involved in adopting the multiple confirmations, we have synthesized analogues involving mutation of the proline residues. We have shown that the proline residues have a significant influence on the structure, activity, and folding of Ov-GRN12–35_3s. These results provide insight into improving the oxidative folding yield and bioactivity of Ov-GRN12–35_3s...

Cone snails produce a fast-acting and often paralyzing venom that is usually injected into their prey or predator through a hypodermic needle-like modified radula tooth. Many diverse compounds are found in their venom including small molecules, peptides and enzymes. However, peptidic toxins called conotoxins (10&ndash, 40 residues and 2&ndash, 4 disulfide bonds) largely dominate these cocktails. These disulfide rich toxins are very valuable pharmacological tools for investigating the function of ions channels, G-protein coupled receptors, transporters and enzymes. Here, we report on the synthesis, structure determination and biological activities of two &alpha, conotoxins, CIA and CIB, found in the predatory venom of the piscivorous species Conus catus. CIA is a typical 3/5 &alpha, conotoxin that blocks the rat muscle type nAChR with an IC50 of 5.7 nM. Interestingly, CIA also inhibits the neuronal rat nAChR subtype &alpha, 3&beta, 2 with an IC50 of 2.06 &mu, M. CIB is a 4/7 &alpha, conotoxin that blocks rat neuronal nAChR subtypes, including &alpha, 2 (IC50 = 128.9 nM) and &alpha, 7 (IC50 = 1.51 &mu, M). High resolution NMR structures revealed typical &alpha, conotoxin folds for both peptides. We also investigated the in vivo effects of these toxins on fish, since both peptides were identified in the predatory venom of C. catus. Consistent with their pharmacology, CIA was highly paralytic to zebrafish (ED50 = 110 &mu, g/kg), whereas CIB did not affect the mobility of the fish. In conclusion, CIA likely participates in prey capture through muscle paralysis, while the putative ecological role of CIB remains to be elucidated.

CAZypedia was initiated in 2007 to create a comprehensive, living encyclopedia of the carbohydrate-active enzymes (CAZymes) and associated carbohydrate-binding modules involved in the synthesis, modification, and degradation of complex carbohydrates. CAZypedia is closely connected with the actively-curated CAZy database, which provides a sequence-based foundation for the biochemical, mechanistic, and structural characterization of these diverse proteins. Now celebrating its 10th anniversary online, CAZypedia is a successful example of dynamic, community-driven, and expert-based biocuration. CAZypedia is an open-access resource available at URL http://www.cazypedia.org.

DNA minor-groove-binding compounds have limited biological applications, in part due to problems with sequence specificity that cause off-target effects. A model to enhance specificity has been developed with the goal of preparing compounds that bind to two AT sites separated by G•C base pairs. Compounds of interest were probed using thermal melting, circular dichroism, mass spectrometry, biosensor-SPR, and molecular modeling methods. A new minor groove binder that can strongly and specifically recognize a single G•C base pair with flanking AT sequences has been prepared. This multi-site DNA recognition mode offers novel design principles to recognize entirely new DNA motifs.

The design, synthesis, and DNA binding properties of azaHx-PI or p-anisyl-4-aza-benzimidazole-pyrrole-imidazole (5) are described. AzaHx, 2-(p-anisyl)-4-aza-benzimidazole-5-carboxamide, is a novel, fluorescent DNA recognition element, derived from Hoechst 33258 to recognize G·C base pairs. Supported by theoretical data, the results from DNase I footprinting, CD, ΔTM, and SPR studies provided evidence that an azaHx/IP pairing, formed from antiparallel stacking of two azaHx-PI molecules in a side-by-side manner in the minor groove, selectively recognized a C–G doublet. AzaHx-PI was found to target 5′-ACGCGT-3′, the Mlu1 Cell Cycle Box (MCB) promoter sequence with specificity and significant affinity (Keq 4.0 ± 0.2 × 107 M−1).

Small-molecule targeting of the DNA minor groove is a promising approach to modulate genomic processes necessary for normal cellular function. For instance, dicationic diamindines, a well-known class of minor groove binding compounds, have been shown to inhibit interactions of transcription factors binding to genomic DNA. The applications of these compounds could be significantly expanded if we understand sequence-specific recognition of DNA better and could use the information to design more sequence-specific compounds. Aside from polyamides, minor groove binders typically recognize DNA at A-tract or alternating AT base pair sites. Targeting sites with GC base pairs, referred to here as mixed base pair sequences, is much more difficult than those rich in AT base pairs. Compound 1 is the first dicationic diamidine reported to recognize a mixed base pair site. It binds in the minor groove of ATGA sequences as a dimer with positive cooperativity. Due to the well-characterized behavior of 1 with ATGA and AT rich sequences, it provides a paradigm for understanding the elements that are key for recognition of mixed sequence sites. Electrospray ionization mass spectrometry (ESI-MS) is a powerful method to screen DNA complexes formed by analogues of 1 for specific recognition. We also report a novel approach to determine patterns of recognition by 1 for cognate ATGA and ATGA-mutant sequences. We found that functional group modifications and mutating the DNA target site significantly affect binding and stacking, respectively. Both compound conformation and DNA sequence directionality are crucial for recognition.

Sequence-specific recognition of DNA by small organic molecules offers a potentially effective approach for the external regulation of gene expression and is an important goal in cell biochemistry. Rational design of compounds from established modules can potentially yield compounds that bind strongly and selectively with specific DNA sequences. An initial approach is to start with common A·T bp recognition molecules and build in G·C recognition units. Here we report on the DNA interaction of a synthetic compound that specifically binds to a G·C bp in the minor groove of DNA by using an azabenzimidazole moiety. The detailed interactions were evaluated with biosensor-surface plasmon resonance (SPR), isothermal calorimetric (ITC), and mass spectrometry (ESI-MS) methods. The compound, DB2277, binds with single G·C bp containing sequences with sub-nanomolar potency and displays slow dissociation kinetics and high selectivity. A detailed thermodynamic and kinetic study at different experimental salt concentrations and temperatures shows that the binding free energy is salt concentration dependent but essentially temperature independent under our experimental conditions, and binding enthalpy is temperature dependent but salt concentration independent. The results show that in the proper compound structural context novel heterocyclic cations can be designed to strongly recognize complex DNA sequences.

Glycoside hydrolase (GH) family 48 is an understudied and increasingly important exoglucanase family found in the majority of bacterial cellulase systems. Moreover, many thermophilic enzyme systems contain GH48 enzymes. Deletion of GH48 enzymes in these microorganisms results in drastic reduction in biomass deconstruction. Surprisingly, given their importance for these microorganisms, GH48s have intrinsically low cellulolytic activity but even in low ratios synergize greatly with GH9 endoglucanases. In this study, we explore the structural and enzymatic diversity of these enzymes across a wide range of temperature optima. We have crystallized one new GH48 module from Bacillus pumilus in a complex with cellobiose and cellohexaose (BpumGH48). We compare this structure to other known GH48 enzymes in an attempt to understand GH48 structure/function relationships and draw general rules correlating amino acid sequences and secondary structures to thermostability in this GH family. Electronic supplementary material The online version of this article (10.1186/s13068-017-0951-5) contains supplementary material, which is available to authorized users.

Background Auxiliary activity (AA) enzymes are produced by numerous bacterial and fungal species to assist in the degradation of biomass. These enzymes are abundant but have yet to be fully characterized. Here, we report the X-ray structure of Thermobifida fusca AA10A (TfAA10A), investigate mutational characterization of key surface residues near its active site, and explore the importance of the various domains of Thermobifida fusca AA10B (TfAA10B). The structure of TfAA10A is similar to other bacterial LPMOs (lytic polysaccharide monooxygenases), including signs of photo-reduction and a distorted active site, with mixed features showing both type I and II copper coordination. The point mutation experiments of TfAA10A show that Trp82 and Asn83 are needed for binding, but only Trp82 affects activity. The TfAA10B domain truncation mutants reveal that CBM2 is crucial for the binding of substrate, but that the X1 module does not affect binding or activity. Results In TfAA10A, Trp82 and Asn83 are needed for binding, but only Trp82 affects activity. The TfAA10B domain truncation mutants reveal that CBM2 is crucial for substrate binding, but that the X1 module does not affect binding or activity. The structure of TfAA10A is similar to other bacterial lytic polysaccharide monooxygenases with mixed features showing both type I and II copper coordination. Conclusions The role of LPMOs and the variability of abundance in genomes are not fully explored. LPMOs likely perform initial attacks into crystalline cellulose to allow larger processive cellulases to bind and attack, but the precise nature of their synergistic behavior remains to be definitively characterized. Electronic supplementary material The online version of this article (doi:10.1186/s13068-017-0925-7) contains supplementary material, which is available to authorized users.

Spider venoms constitute incredibly diverse libraries of compounds, many of which are involved in prey capture and defence. Polyamines are often prevalent in the venom and target ionotropic glutamate receptors. Here we show that a novel spider polyamine, PA366, containing a hydroxyphenyl-based structure is present in the venom of several species of tarantula, and has selective toxicity against MCF-7 breast cancer cells. By contrast, a polyamine from an Australian funnel-web spider venom, which contains an identical polyamine tail to PA366 but an indole-based head-group, is only cytotoxic at high concentrations. Our results suggest that the ring structure plays a role in the cytotoxicity and that modification to the polyamine head group might lead to more potent and selective compounds with potential as novel cancer treatments.

To find their DNA target sites in complex solution environments containing excess heterogeneous DNA, sequence-specific DNA-binding proteins execute various translocation mechanisms known collectively as facilitated diffusion. For proteins harboring a single DNA contact surface, long-range translocation occurs by jumping between widely spaced DNA segments. We have configured biosensor-based surface plasmon resonance to directly measure the affinity and kinetics of this intersegmental jumping by the ETS-family transcription factor ETS variant 6 (ETV6). To isolate intersegmental target binding in a functionally defined manner, we pre-equilibrated ETV6 with excess salmon sperm DNA, a heterogeneous polymer, before exposing the nonspecifically bound protein to immobilized oligomeric DNA harboring a high-affinity ETV6 site. In this way, the mechanism of ETV6-target association could be toggled electrostatically through varying NaCl concentration in the bulk solution. Direct measurements of association and dissociation kinetics of the site-specific complex indicated that 1) freely diffusive binding by ETV6 proceeds through a nonspecific-like intermediate, 2) intersegmental jumping is rate-limited by dissociation from the nonspecific polymer, and 3) dissociation of the specific complex is independent of the history of complex formation. These results show that target searches by proteins with an ETS domain, such as ETV6, whose single DNA-binding domain cannot contact both source and destination sites simultaneously, are nonetheless strongly modulated by intersegmental jumping in heterogeneous site environments. Our findings establish biosensors as a general technique for directly and specifically measuring target site search by DNA-binding proteins via intersegmental translocation.

Thermostable α-amylases have many industrial applications and are therefore continuously explored from novel sources. We present the characterization of a novel putative α-amylase gene product (Tp-AmyS) cloned from Thermotoga petrophila. The purified recombinant enzyme is highly thermostable and able to hydrolyze starch into dextrin between 90 and 100°C, with optimum activity at 98°C and pH8.5. The activity increased in the presence of Rb1+, K1+ and Ca2+ ions, whereas other ions inhibited activity. The crystal structure of Tp-AmyS at 1.7A resolution showed common features of the GH-13 family, however was apparently found to be a dimer. Several residues from one monomer interacted with a docked acarbose, an inhibitor of Tp-AmyS, in the other monomer, suggesting catalytic cooperativity within the dimer. The most striking feature of the dimer was that it resembled the dimerization of salivary amylase from a previous crystal structure, and thus could be a functional feature of some amylases.

International audience; Conotoxins are a large family of disulfide-rich peptides that contain unique cysteine frameworks that target a broad range of ion channels and receptors. We recently discovered the 33-residue conotoxin Φ-MiXXVIIA from Conus miles with a novel cysteine framework comprising three consecutive cysteine residues and four disulfide bonds. Regioselective chemical synthesis helped decipher the disulfide bond connectivity and the structure of Φ-MiXXVIIA was determined by NMR spectroscopy. The 3D structure displays a unique topology containing two β-hairpins that resemble the N-terminal domain of granulin. Similar to granulin, Φ-MiXXVIIA promotes cell proliferation (EC50 17.85 μm) while inhibiting apoptosis (EC50 2.2 μm). Additional framework XXVII sequences were discovered with homologous signal peptides that define the new conotoxin superfamily G2. The novel structure and biological activity of Φ-MiXXVIIA expands the repertoire of disulfide-rich conotoxins that recognize mammalian receptors.

Bleomycins A5 and B2 were used to study the structural features in hairpin DNAs conducive to strong BLM–DNA interaction. Two members of a 10-hairpin DNA library previously found to bind most tightly to these BLMs were subsequently noted to share the sequence 5′-ACGC (complementary strand sequence 5′-GCGT). Each underwent double-strand cleavage at five sites within, or near, an eight base pair region of the DNA duplex which had been randomized to create the original library. A new hairpin DNA library was selected based on affinity for immobilized Fe(III)·BLM A5. Two of the 30 newly identified DNAs also contained the sequence 5′-ACGC/5′-GCGT. These DNAs bound to the Fe(II)·BLMs more tightly than any DNA characterized previously. Surface plasmon resonance confirmed tight Fe(III)·BLM B2 binding and gave an excellent fit for a 1:1 binding model, implying the absence of significant secondary binding sites. Fe(II)·BLM A5 was used to assess sites of double-strand DNA cleavage. Both hairpin DNAs underwent double-strand cleavage at five sites within or near the original randomized eight base region. For DNA 12, four of the five double-strand cleavages involved independent single-strand cleavage reactions; DNA 13 underwent double-strand DNA cleavage by independent single-strand cleavages at all five sites. DNA 14, which bound Fe·BLM poorly, was converted to a strong binder (DNA 15) by insertion of the sequence 5′-ACGC/5′-GCGT. These findings reinforce the idea that tighter DNA binding by Fe·BLM leads to increased double-strand cleavage by a novel mechanism and identify a specific DNA motif conducive to strong BLM binding and cleavage.

ETS family transcription factors regulate diverse genes through binding at cognate DNA sites that overlap substantially in sequence. The DNA-binding domains of ETS proteins (ETS domains) are highly conserved structurally yet share limited amino acid homology. To define the mechanistic implications of sequence diversity within the ETS family, we characterized the thermodynamics and kinetics of DNA site recognition by the ETS domains of Ets-1 and PU.1, which represent the extremes in amino acid divergence among ETS proteins. Even though the two ETS domains bind their optimal sites with similar affinities under physiologic conditions, their nature of site recognition differs strikingly in terms of the role of hydration and counter ion release. The data suggest two distinct mechanisms wherein Ets-1 follows a "dry" mechanism that rapidly parses sites through electrostatic interactions and direct protein-DNA contacts, whereas PU.1 utilizes hydration to interrogate sequence-specific sites and form a long-lived complex relative to the Ets-1 counterpart. The kinetic persistence of the high affinity PU.1 · DNA complex may be relevant to an emerging role of PU.1, but not Ets-1, as a pioneer transcription factor in vivo. In addition, PU.1 activity is critical to the development and function of macrophages and lymphocytes, which present osmotically variable environments, and hydration-dependent specificity may represent an important regulatory mechanism in vivo, a hypothesis that finds support in gene expression profiles of primary murine macrophages.

Minor groove binding compounds have been shown to induce changes in global DNA conformation, allosterically inhibiting DNA-protein interactions necessary for transcriptional processes. Many minor groove binders are specific for AT-base pairs but have little preference over alternating AT or A-tract sequences. Few compounds, other than polyamides, show selectivity for mixed sequences with AT and GC base pairs. Electrospray ionization mass spectrometry (ESI-MS) can provide insight on the stoichiometry and relative affinities in minor groove recognition of different DNA sequences with a library of minor groove binders. A goal in our current research is to develop new compounds that recognize mixed sequences of DNA. In an effort to optimize screening for compounds that target mixed AT and GC base pair sequences of DNA, ESI-MS was used to study the competitive binding of compounds with a mixed set of DNA sequences. The method identified preferred binding sites, relative affinities, and concentration-dependent binding stoichiometry for the minor groove binding compounds netropsin and DB75 with AT-rich sequences, and DB293 with ATGA and AT-sites.

Molecular dynamics simulations were used to study the possible catalytic role of an unusual conserved water-filled pore structure in the family 48 cellulase enzyme Cel48A from Thermobifida fusca. It was hypothesized that this pore serves as the pathway for the water molecules consumed in the hydrolysis catalyzed by the enzyme to reach the active site in a continuous stream to participate in the processive reactions. Theoretical mutants of this enzyme were created in which all of the residues lining the pore were made hydrophobic, which had the effect in molecular dynamics simulations of emptying the pore of water molecules and preventing any from passing through the pore on the simulation time scale. Mutants with smaller numbers of substitutions of this nature, which could be created experimentally by site-directed mutagenesis, were also identified from simulations, and these proteins were subsequently produced in Escherichia coli, expressed and purified, but were found to not fold in a manner similar to the wild type protein, preventing the determination of the importance of the water pore for activity. It is possible that the presence of a small vacuum in the pore was responsible for the instability of the mutants. In addition, alternate pathways were observed in the simulations that would allow water molecules to reach the active site of the enzyme, suggesting that the hypothesis that the pore has functional significance might be incorrect.

DB1255 is a symmetrical diamidinophenyl-dithiophene that exhibits cellular activity by binding to DNA and inhibiting binding of ERG, an ETS family transcription factor that is commonly overexpressed or translocated in leukemia and prostate cancer [Nhili, R., Peixoto, P., Depauw, S., Flajollet, S., Dezitter, X., Munde, M. M., Ismail, M. A., Kumar, A., Farahat, A. A., Stephens, C. E., Duterque-Coquillaud, M., Wilson, W. D., Boykin, D. W., and David-Cordonnier, M. H. (2013) Nucleic Acids Res. 41, 125–138]. Because transcription factor inhibition is complex but is an attractive area for anticancer and antiparasitic drug development, we have evaluated the DNA interactions of additional derivatives of DB1255 to gain an improved understanding of the biophysical chemistry of complex function and inhibition. DNase I footprinting, biosensor surface plasmon resonance, and circular dichroism experiments show that DB1255 has an unusual and strong monomer binding mode in minor groove sites that contain a single GC base pair flanked by AT base pairs, for example, 5′-ATGAT-3′. Closely related derivatives, such as compounds with the thiophene replaced with furan or selenophane, bind very weakly to GC-containing sequences and do not have biological activity. DB1255 is selective for the ATGAT site; however, a similar sequence, 5′-ATGAC-3′, binds DB1255 more weakly and does not produce a footprint. Molecular docking studies show that the two thiophene sulfur atoms form strong, bifurcated hydrogen bond-type interactions with the G-N-H sequence that extends into the minor groove while the amidines form hydrogen bonds to the flanking AT base pairs. The central dithiophene unit of DB1255 thus forms an excellent, but unexpected, single-GC base pair recognition module in a monomer minor groove complex.

Hairpin polyamides (PAs) are an important class of sequence-specific DNA minor groove binders, and frequently employ a flexible motif, β-alanine (β), to reduce the molecular rigidity to maintain the DNA recognition register. To better understand the diverse effects that β can have on DNA-PA binding affinity, selectivity, and especially kinetics, which have rarely been reported, we have initiated a detailed study for an eight-heterocyclic hairpin PA and its β derivatives with their cognate and mutant sequences. With these derivatives, all internal pyrroles of the parent PA are systematically substituted with single or double βs. A set of complementary experiments have been conducted to evaluate the molecular interactions in detail: UV-melting, biosensor-surface plasmon resonance, circular dichroism and isothermal titration calorimetry. The β substitutions generally weaken the binding affinities of these PAs with cognate DNA, and have large and diverse influences on PA binding kinetics in a position- and number-dependent manner. The DNA base mutations have also shown positional effects on the binding of a single PA. Besides the β substitutions, the monocationic Dp group [3-(dimethylamino)propylamine] in parent PA has been modified into a dicationic Ta group (3,3'-diamino-N-methyldipropylamine) to minimize the frequently observed PA aggregation with ITC experiments. The results clearly show that the Ta modification not only maintains the DNA binding mode and affinity of PA, but also significantly reduces PA aggregation and allows the complete thermodynamic signature of eight-ring hairpin PA to be determined for the first time. This combined set of results significantly extends our understanding of the energetic basis of specific DNA recognition by PAs.

It is now well established that, although only about 5% of the human genome codes for protein, most of the DNA has some function, such as synthesis of specific, functional RNAs and/or control of gene expression. These functional sequences open immense possibilities in both biotechnology and therapeutics for the use of cell-permeable, small molecules that can bind mixed-base pair sequences of DNA for regulation of genomic functions. Unfortunately very few types of modules have been designed to recognize mixed DNA sequences and for progress in targeting specific genes, it is essential to have additional classes of compounds. Compounds that can be rationally designed from established modules and which can bind strongly to mixed base pair DNA sequences are especially attractive. Based on extensive experience in design of minor-groove agents for AT recognition, a small library of compounds with two AT specific binding modules, connected through linkers which can recognize the G·C base pairs, were prepared. The compound-DNA interactions were evaluated with a powerful array of biophysical methods and the results show that some pyridyl-linked compounds bind with the target sequence with sub-nanomolar KD, with very slow dissociation kinetics and 200 times selectivity over the related sequence without a G·C base pair. Interestingly, a set of compounds with AT module connected by different linkers shows cooperative dimer recognition of related sequences. This type of design approach can be expanded to additional modules for recognition of a wide variety of sequences.

Heterocyclic diamidines are strong DNA minor-groove binders and have excellent antiparasitic activity. To extend the biological activity of these compounds, a series of arylimidamides (AIAs) analogues, which have better uptake properties in Leishmania and Trypanosoma cruizi than diamidines, was prepared. The binding of the AIAs to DNA was investigated by Tm , fluorescence displacement titration, circular dichroism, DNase I footprinting, biosensor surface plasmon resonance, X-ray crystallography and molecular modeling. These compounds form 1:1 complexes with AT sequences in the DNA minor groove, and the binding strength varies with substituent size, charge and polarity. These substituent-dependent structure and properties provide a SAR that can be used to estimate K values for binding to DNA in this series. The structural results and molecular modeling studies provide an explanation for the differences in binding affinities for AIAs.

Lignocellulosic biomass is a potential source of renewable, low-carbon-footprint liquid fuels. Biomass recalcitrance and enzyme cost are key challenges associated with the large-scale production of cellulosic fuel. Kinetic modeling of enzymatic cellulose digestion has been complicated by the heterogeneous nature of the substrate and by the fact that a true steady state cannot be attained. We present a two-parameter kinetic model based on the Michaelis-Menten scheme ( Michaelis, L., and Menten, M. L. ( 1913 ) Biochem. Z. , 49 , 333 - 369 ) with a time-dependent activity coefficient analogous to fractal-like kinetics formulated by Kopelman ( Kopelman, R. ( 1988 ) Science 241 , 1620 - 1626 ). We provide a mathematical derivation and experimental support to show that one of the parameters is a total activity coefficient and the other is an intrinsic constant that reflects the ability of the cellulases to overcome substrate recalcitrance. The model is applicable to individual cellulases and their mixtures at low-to-medium enzyme loads. Using biomass degrading enzymes from cellulolytic bacterium Thermobifida fusca , we show that the model can be used for mechanistic studies of enzymatic cellulose digestion. We also demonstrate that it applies to the crude supernatant of the widely studied cellulolytic fungus Trichoderma reesei ; thus it can be used to compare cellulases from different organisms. The two parameters may serve a similar role to Vmax, KM, and kcat in classical kinetics. A similar approach may be applicable to other enzymes with heterogeneous substrates and where a steady state is not achievable.

Orthogonally positioned diamino/dicationic polyamides (PAs) have good water solubility and enhanced binding affinity, whilst retaining DNA minor groove and sequence specificity compared to their monoamino/monocationic counterparts. The synthesis and DNA binding properties of the following diamino PAs: f-IPI (3a), f-IPP (4), f-PIP (5), and f-PPP (6) are described. P denotes the site where a 1-propylamino group is attached to the N1-position of the heterocycle. Binding of the diamino PAs to DNA was assessed by DNase I footprinting, thermal denaturation, circular dichroism titration, biosensor surface plasmon resonance (SPR), and isothermal titration calorimetry (ITC) studies. According to SPR studies, f-IPI (3a) bound more strongly (K(eq)=2.4×10(8) M(-1)) and with comparable sequence selectivity to its cognate sequence 5'-ACGCGT-3' when compared to its monoamino analog f-IPI (1). The binding of f-IPI (3a) to 5'-ACGCGT-3' via the stacked dimer motif was balanced between enthalpy and entropy, and that was quite different from the enthalpy-driven binding of its monoamino parent f-IPI (1). f-IPP (4) also bound more strongly to its cognate sequence 5'-ATGCAT-3' (K(eq)=7.4×10(6) M(-1)) via the side-by-side stacked motif than its monoamino analog f-IPP (2a). Although f-PPP (6) bound via a 1:1 motif, it bound strongly to its cognate sequence 5'-AAATTT-3' (K(eq)=4.8×10(7) M(-1)), 15-times higher than the binding of its monoamino analog f-PPP (2c), albeit f-PPP bound via the stacked motif. Finally, f-PIP (5) bound to its target sequence 5'-ATCGAT-3' as a stacked dimer and it has the lowest affinity among the diamino PAs tested (Keq

Synthetic N-methyl imidazole and N-pyrrole containing polyamides (PAs) that can form “stacked” dimers can be programmed to target and bind to specific DNA sequences and control gene expression. To accomplish this goal, the development of PAs with lower molecular mass which allows for the molecules to rapidly penetrate cells and localize in the nucleus, along with increased water solubility, while maintaining DNA binding sequence specificity and high binding affinity is key. To meet these challenges, six novel f-ImPy*Im PA derivatives that contain different orthogonally positioned moieties were designed to target 5′-ACGCGT-3′. The synthesis and biophysical characterization of six f-ImPy*Im were determined by CD, ΔTM, DNase I footprinting, SPR, and ITC studies, and were compared with those of their parent compound, f-ImPyIm. The results gave evidence for the minor groove binding and selectivity of PAs 1 and 6 for the cognate sequence 5′-ACGCGT-3′, and with strong affinity, Keq = 2.8 × 108 M−1 and Keq = 6.2 × 107 M−1, respectively. The six novel PAs presented in this study demonstrated increased water solubility, while maintaining low molecular mass, sequence specificity, and binding affinity, addressing key issues in therapeutic development. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 497–507, 2013.