SciELO - Scientific Electronic Library Online

 
vol.78A Novel Biopolymer-based Nanomagnetic Catalyst for the Synthesis of 4H-pyran and Tetrahydro-4H-chromene DerivativesSurface-modified Water Hyacinth (Eichhornia crassipes) over Activated Carbon for Wastewater Treatment: A Comparative Account índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Artigo

Indicadores

Links relacionados

  • Em processo de indexaçãoCitado por Google
  • Em processo de indexaçãoSimilares em Google

Compartilhar


South African Journal of Chemistry

versão On-line ISSN 1996-840X
versão impressa ISSN 0379-4350

S.Afr.j.chem. (Online) vol.78  Durban  2024

http://dx.doi.org/10.17159/0379-4350/2024/v78a10 

RESEARCH ARTICLE

 

Application of Derandomisation to Uperin 3.x Peptides Circular Dichroism Spectra to Determine their secondary structure contents

 

 

Adewale OlamoyesanI, II, III, ; Alison RodgerI

IDepartment of Molecular Sciences, Faculty of Science and Engineering, Macquarie University, Macquarie Park, Australia
IIDepartment of Chemistry, University of Lagos, Akoka, Lagos, Nigeria
IIICollege of Science and Computing, Wigwe University, Isiokpo, Rivers, Nigeria

 

 


ABSTRACT

Survival of frogs and toads in hostile environments depends on their ability to secrete antimicrobial or host-defence peptides. The Uperin, U3.x, peptides investigated were originally isolated from the skin secretions of a frog endemic to Australia. In this paper, we extend the automated ad hoc approach for derandomisation of proteins to U3.x, whose spectra are collected in water and after the addition of buffer to eventually understand how this affects the structure's goodness of fit and spectral prediction. Systematically with the "CD app' disordered unit' varying fractions of random coil, RC spectrum (10-90% in increments of 10%) were removed from U3.x's spectra baseline corrected (in molar extinction units) to produce derandomised spectra. Self-organising map spectroscopy, SOMSpec, analysis gives the secondary structure of the derandomised U3.x, afterward, the RC component was added back to determine the percentage of the secondary structure motifs in the original protein or peptide. Most of the applications to U3.x wild types, WTs, and their mutants were straightforward, though, a few of them were not. These exceptions required visual inspection of the spectra. We also found that as the percentage of derandomisation increases, the observed behaviour is like proteins, where folded structures decrease.

Keywords: uperin, derandomised, normalised root mean square deviation, structure, circular dichorism


 

 

INTRODUCTION

A vital component of the innate immune systems of most plants and animals is the antimicrobial peptide, AMP, or host defence peptide.1 Suitable stress stimulus in the hostile environment where these organisms live, triggers the activity of these peptides against any invading microbe(s).2 For the past three decades, the skin secretions of Australian frogs and toads were the source of several host defence peptides. These isolated peptides are rich in AMPs and, or they have one or more neuropeptide activities. Their prevalence in nature demonstrates AMP's effectiveness as a defence compound.3-5 They are attractive as starting points in developing new peptide-based antimicrobial agents due to their high activity and selectivity for microbes. The interaction of AMPs with microorganisms' lipid bilayer is always fatal; and subsequently leads to the disruption of their membranes in a specific, but not receptor-mediated, process.2

Uperin 3.x (U3.x) family is a class of peptides originally isolated from the skin secretions of Uperoleia mjobergii (a frog endemic to Australia). This type of peptide is known to have antimicrobial activity against Gram-positive bacteria and mammalian cell lines.6,7 All three U3.x used herein have 17 amino acid residues, as shown in Table 1.8 Among them, U3.4 WT and U3.5 WT have identical peptide sequences, with U3.4 WT being slightly more nonpolar. More than half of the individual positions of the amino acids of U3.6 WT differ from equivalent positions on the other two U3.x WT peptides. Except for U3.5 wt that self-aggregates in a saline buffer the rest are unstructured either in water or buffer.9 Ray and her group members10 observed that U3.x WT aggregates less than their corresponding seventh position variants. Based on this behaviour, Ray et al. rationalised that substituting nonpolar alanine residue for the positively charged amino acid at the seventh position increases fibril formation. In addition, this group conducted mechanistic studies of the initial stages of self-aggregation of these U3.x WT peptides and their variants using molecular dynamics, MD, simulations and circular dichroism, CD, measurements. Their results show that the U3.x WT peptides demonstrated a lower propensity to β-sheet than the corresponding alanine variants. Likewise, the addition of salt increased hydrophobic interactions between peptides, which resulted in enhanced peptide aggregation.

A well-known assumption in structural biology is that the secondary structure of proteins is a prerequisite to their activity. To the best of our knowledge, there are quite a few methods to estimate the peptide secondary structure, 11-13 therefore, there is a need for new or existing methods never used for peptides to be adapted for this purpose. Though the secondary structure determination for proteins is well-established,14-18 that for peptides is less successful, e.g., as demonstrated by San Miguel and coworkers.13 The authors estimated helical structure by assuming that 100% helical peptides have 208 nm and 222 nm CD signals of -12 mol-1 dm3 cm-1 and the average 208 nm and 222 nm signal divided by -12 gave the fractional helix content. In general, a sample that looks partly folded on inspection of the CD spectrum will include a population of unfolded and folded peptides.

This work explores whether the ad hoc approach, systematically automated when complemented with self-organizing map spectroscopy, SOMSpec, is useful for the analysis of heterogeneous peptide populations. 19-21 It further discusses how the extent of derandomisation on the spectra of three U3.x peptides and their seventh position alanine variants in water and after the addition of buffer over time affects the overall goodness of fit of the structure and spectral output from SOMSpec.

 

MATERIALS AND METHODS

Lisa Martin of Monash University supplied us with the CD spectra for this work. The materials and methods used are stated in Ray et al.10

SOMSpec Input Data Preparation

We ensured all the input data for SOMSpec analysis was saved in text format. Firstly, the CD spectra were saved in a comma-separated variable, CSV, format with each spectrum placed in a column, with the reference spectra's structures forming five additional entries in each of its columns while that of unknown peptides (experimental spectra) lack structure entries. Secondly, the inspection of the worksheet named "SOMSpec" in the SP175 reference set file was carried out to ensure that all spectra have a wavelength range of 240-190 nm (the last five of the data points are for annotated structures).20 From this original reference set, two additional reference sets were created with the wavelength ranges of 240-195 nm and 240-200 nm, respectively having 46 and 41 data points for each spectrum with the same number of structure entries as that in 240 -190 nm annotated to it.

Thirdly, the raw CD spectra were created into test files truncated to a wavelength range of 240-190 nm, 240-195 nm, and 240-200 nm with a step size of 1 nm (so the data points are 51, 46, and 41 per spectrum for the wavelength ranges respectively). The test sets were converted to ∆ε per molar residue by a scale factor of 3298, then baseline corrected. Lastly, derandomised spectra were automatically created by removing different fractions of the RC spectrum from each ∆ε peptide spectrum.20,21

SOMSpec Methodology

A square map with 2 500 nodes is created using map size (50 x 50), iterations (50 000), best matching units, BMUs (5). Before training, the weight vectors to populate the nodes within the map are initialised. Next, as the iteration begins, the input vectors are randomly presented to a cluster of the nodes in the map, and the weight vectors associated with each adapt them. At certain iterations associating weight vectors become more identical with the input one, and the BMU input vector is determined using Euclidean distance. Vectors near the BMU with remarkably identical coordinates in the weight map (neighbouring nodes) are brought closer to BMU as training is repeated by a factor L (learning rate, 0.1). The learning rate and neighbouring function used in SOMSpec are based on power series and Gaussian function respectively. Then the BMU, neighbourhood nodes, and other vectors will be assigned to structure labels that will show in the output SOM. Expectedly the five structural labels in the reference spectra must add up to 1, as well as the predicted structure types. Afterwards, the trained map with the secondary structure, SS, property is used to test the unknown spectrum by considering the user-defined number of BMUs (set as 5) in the map. With this value, the nodes very similar to BMUs in the test spectrum are iteratively identified by a distance-dependent model known as NRMSD 22. The predicted spectrum is estimated by interpolating across the BMUs and is visualised in the plot accompanying the 2-D representation of SOM. Each trained map was used to predict the spectrum and structure of each peptide (test spectrum) in water only, and for which readings were taken at t = 0, t = 1, and t = 24 h after adding buffer with corresponding wavelength range.

 

RESULTS AND DISCUSSION

CD Spectra Tested Against Reference Set with a Wavelength Range of 240-190

The essence of taking measurements after adding buffer is to investigate the effect of buffer over time on peptides' structure

evolution or stability. Figure 1 shows the far-UV CD spectra for the aqueous solutions of U3.4 WT, U3.5 WT, U3.6 WT and their seventh position variants. All the peptides in water show a negative maximum at 198 nm, indicating they are mostly RC or unstructured peptides in solution. However, only U3.5 WT gained other spectral features after adding a buffer to the peptide solutions. The other two U3.x WT still maintain the negative maximum of the RC structure at 24 h (Figure 1a, c), while that for U3.5 WT shifted to 204 nm at 1 h. This is accompanied by the appearance of a broad negative maximum at 218 nm, which indicates a conformation transition from RC to a-helix and, or β-sheet. Finally, a positive maximum peak at 196 nm appeared together with a rise in the magnitude of the 218 nm peak, which seems to suggest a further increase in a-helix and, or β-sheet structure (from ~10% a-helix and ~26% β-sheet at t = 0 to ~15% a-helix and ~38% β-sheet at t =24 (Figures 1b, d, f and S21)) at 24 h. The transition from a few percentages of folded structures to several percentages of folded structures is more for the seventh position variants than their respective U3.x wt.

The direct structure estimates and spectra fitting for the experimental spectra with SOMSpec using the augmented SP175 reference set were unsatisfactory for all the U3.x peptides in water and after the addition ofthe buffer (Figure S1-S18). Figure 2 shows the spectral NRMSD plots for U3.4 WT, U3.5 WT, U3.6 WT, and U3.4 R7A recorded in water and at t = 0, 1, and 24 h after adding buffer before derandomisation of the spectra. See Table 2 for the predicted a-helix and β-sheet at 24 h after adding the buffer to these U3.x peptides.

The spectra were then systematically derandomised by the "CD app' disordered unit" dedicated for the task as described by Olamoyesan et al. 21,22 Equation 1, as stated in our previous work, 21 mathematically expresses how varying fractions of RC are removed from the baseline corrected spectra in molar extinction units. The thirty-six spectra resulting from this process for each peptide were fitted with SOMSpec to generate thirty-six spectral predictions with NRMSDs and associate structure estimates. The difficulty for SOMSpec to find perfect places on the trained map are indicated by the NRMSDs, which mostly are above the nominated target value (0.03), except for U3.5 R7A -30%RC at t =1 (0.025). We assessed multiple good fits by weighing the maxima, minima, and where the spectrum crosses zero while overlaying the predicted and the experimental spectra to determine the best of the fits.

Best Fit for U3.4 WTand U3.4 R7A Spectra Either Derandomised or Not

Assessment of the best options for these peptides was carried out using mainly their NRMSD values, except for two instances. Table 3 shows the best spectral NRMSDs and their associate structure estimates for an aqueous solution of U3.4 WT and U3.4 R7A only, after adding buffer, and subtracting (0.1-0.9) RC from the spectra. Almost all of the best choices of these peptides are obvious except for U3.4 R7A in water only with -60% RC and -70% RC, and at t = 0 h (immediately after adding buffer) with -20% RC and -40% RC. Though the NRMSDs were almost the same for the two options at each condition, the best choices were found to be -70% RC and -40% RC respectively before the addition of buffer and afterward, at t = 0 h (Figure 3b,c,e, and f) based on spectrum shapes match. Alternatively, the best choice can be decided without necessarily inspecting the spectra by implementing a refinement with the percentage RC between the two observed options (lowest and very similar NRMSDs). There is a clear indication that at t = 1 and 24 h the same extent of derandomisation (-30%RC) gives best fits with similar structure estimates.

There is a gradual decrease in β-sheet, bend and turn contents as derandomisation increases in the mutant of U3.4 WT in water only and after adding buffer as well as in the WT. The mutation causes a significant increase in a-helix with, say 120.51%, 176.82%, and 282.53% respectively at 0, 1, and 24 h of adding buffer in comparison with its WT. Moreover, the percentage of other folded structures aside helix increased with a mean average of 26.92%, except β-sheet at t =1, and 24 h after adding buffer which decreased respectively by 15.79% and 24.96% from the WT (Figure S19a and S20a). Although all the right answers for the WT were clear, the extent of derandomisation associated with the best fit increased from 0 to 0.9 after dropping by 0.1 and plummeting to 0 at t = 24 h.

Best Fit for U3.5 WT and U3.5 R7A Spectra Either Derandomised or Not

In general, there were no hidden challenges when assessing the best fit for U3.5 WT and U3.5 R7A spectral predictions and structure estimates for either derandomised or otherwise. The best fit for U3.5 WT in water only is with -90%RC (which results in 3.6% a-helix and 2.1% β-sheet (regenerated)). Immediately after adding the buffer (at t = 0 h), the derandomisation for best fit dropped to -80% RC, which gives rise to 4.3% a-helix and 6.9% β-sheet. This represents a 19.44% and 223.80% increment from the estimates of the initial condition. At 1 h the percentage derandomisation associated with the best fit fell to- 50% RC and remained the same at 24 h, which resulted in a further percentage increase of the most notable folded structures to 8.9% and 18.2% a-helix and 15.8% and 10.1% β-sheet, respectively from the initial conditions. Likewise, the structure evolution of the mutant of U3.5 WT follows a similar trend as the U3.5 WT respective to the change of extent of derandomisation associated with the best fit. We observed that the best choice for U3.5 R7A in water only has -80% RC removed from it with estimates of 4.7% a-helix, and 5.8% β-sheet. At t = 0 h after adding buffer the extent of derandomisation associated with the best fit is the same as in water only, and this resulted in a 20% increase of the a-helix while the β-sheet content dropped by 40%. However, at 1 h of adding the buffer, the best fit is with -40% RC which yielded an increment of 56.16% and 435% for each of the notable structures (11.4% a-helix and 21.4% β-sheet). Furthermore, the extent of deandomisation associated with best fit dropped to -10% RC at 24 h of placing the sample in water and buffer mixture; this gave rise to 16.6% a-helix and 31.2% β-sheet. Helical content of measurements taken in water only, at t = 0, and 1 h compared well, but moderately increased by 14.55 - 16.56% for both the WT and its mutant at t = 24 h after buffer addition (without no RC removed). Structure estimates of a-helix, β-sheet, bend and turn increase at 24 h after the addition of buffer. At the same time, unstructured content varies inversely to the β-sheet before derandomisation of the spectra for both U3.5 WT and U3.5 R7A. However, derandomisation of the spectra with varying fractions of RC causes a decrease in β-sheet, bend, and turn contents as the unstructured content increases. Overall, the drop or rise in the extent of derandomisation associated with the best answers is expected to cause a decrease or an increase in the unfolded structure contents (Table 4 and Figure S21-22).

Best Fit for U3.6 WT and U3.6 K7A Spectra Either Derandomised or Not

The results indicative of the best fits for both U3.6 WT and its mutant are summarised in Table 5. As can be seen, the extent of derandomisation that can be credited for the best fit increased to -90% RC from -80% RC for the WT after the addition of buffer at 1 h and subsequently plummeted to -0%RC at 24 h. This drastic drop resulted in 10.3% a-helix and 25.4% β-sheet (Figure S23). For its mutant, almost all cases with no derandomisation resulted primarily as best fits, except at 1 h after adding the buffer to the peptide solution rose to -90% RC. The structure estimates for helical content are identical in three instances (in water only, at t =0 and t =1), while at 24 h of adding buffer, the a-helix increased to 25.05% (Figure S24). It appears that as we sequentially increase the fractions of RC subtracted from these peptides (e.g., from -50%-60% RC), this directly leads to an 8% average increase of unstructured content at each step. At later times than 0 h, the mutation causes an increase of 3-151% in a-helix, and for all the conditions bend increases gradually from 0.24-23.32% at t = 24 h (Figure S23 and S24) and yet β-sheet and turn increase in three of the measured conditions.

CD Spectra Tested Against Reference Set with the Same Wavelength Range of 240-195 nm

There is a feeling that the reduction of wavelength range down to 195 or 200 nm affects the quality of the structure information obtainable from the spectra. Thus, the aim in this section is to explore whether this was the case by creating a reference set and test file with spectra down to 195 nm, and then putting the spectra through SOMSpec. This research question will be answered by comparing the SOMSpec NRMSDs for spectra down to 190, 195, and 200 nm, with or without derandomisation.

Best Fit for U3.4 WT and U3.4 R7A Spectra Either Derandomised or Not from 240-195 nm

Table 6 displays the best NRMSDs and their structure estimates and percentage derandomisation added to U3.4 WT and its mutant spectra truncated to 195 nm that effect best fits. The best fit for U3.4 WT in water only has a lower extent of derandomisation (-30%RC) and its NRMSD value (less 0.05), suggests the structure estimates and fitted spectral are of reasonable fit. This resulted in 8.2% a-helix and 18% β-sheet. Immediately after adding the buffer to the peptide solution, the extent of derandomisation associated with the best fit dropped to -10%RC, to give 10.2% a-helix and 21.5% β-sheet. At t = 1 h, the extent of derandomisation associated with the best fit rose to 80%RC, expectedly with a drop of about 68% for both a-helix and a-sheet compared to previous condition's estimates. Two similar lowest NRMSDs were observed at t = 24 h; and the best choice was determined to be with -10%RC following the approach used in the previous occurrence. Therefore, the estimates for both a-helix and β-sheet appreciate by 203% and 232% respectively. We recognised for the mutant that its extent of derandomisation that caused the best fit, dropped immediately after adding the buffer and later rose moderately at 1 h and at the end of the reading rose back to the initial value (-60%RC). It appears that spectra whose best fits have no or lower derandomisation contain more folded structures than those with a higher degree of derandomisation. For example, at t =1 h after adding buffer the estimates for a-helix rose by 171.26%.

Almost all the U3.4 WT spectra reduced down to 195 nm gave a better fit when compared to those down to 190 nm, except for three instances when dissolved in water only (with -0%RC and -80%RC) and immediately after the addition of buffer (with -90%RC). Expectedly, data down to 200 nm are worst, when compared to that down to 190 and 195 nm. In twenty-nine of the forty instances, data down to 195 nm gave the best fit when compared to the other wavelength ranges (190 nm and 200 nm) (Table S1). We hypothesised that these twenty-nine spectra are more similar to their BMUs in the training map than that of wavelength ranges down to 190 nm and 200 nm. As the percentage RC removed increases a gradual decrease in the four structure types (namely helix, sheet, bend, and turn) was observed, while unstructured content increases (Figure S25 and S26).

Best Fit for U3.5 WT and U3.5 R7A Spectra Either Derandomised or Not from 240-195 nm

The best fits for U3.5 WT in water only and at 24 h of adding buffer have the same extent of derandomisation, however, the structure estimates for a-helix and β-sheet differ by 10.6% and 56%, respectively. There is a significant drop in the extent of derandomisation required for best fit immediately after the addition of a buffer but a slight fall to -50%RC was observed at 1 h. Though the spectral trends for U3.5 R7A were similar to that of U3.5 WT, its extent of derandomisation for best fit immediately after the addition of a buffer dropped and further dropped at 1 h and 24 h, respectively to -30%RC and 0%RC. We observed its β-sheet content increased over time after adding a buffer to 34.5% at t = 24 h (Table 7). Twenty-six of the forty spectra of U3.5 WT down to 195 nm showed improvement in their quality and estimate of secondary structure when compared to original data, whereas overall data down to 200 nm were worst. The behaviour of U3.5 R7A spectra down to 195 nm was similar to the case of its U3.5 WT when compared to original data and that down to 200 nm (Table S2).

Best Fit for U3.6 WT and U3.6 K7A Spectra Either Derandomised or Not from 240-195 nm

Other measurements taking aside that of U3.6 WT in water only have the same percentage of derandomisation for their best fits and nearly exact structure estimates for all types (7.5% a-helix and 14.3% β-sheet). There is a significant drop in the extent of derandomisation associated with the mutant best fit immediately after the addition of buffer to 0%RC, as a result, the estimates appreciated by 234.78% and 64.04% respectively for a-helix and β-sheet. At 1 h, this dropped to -60%RC, and subsequently at 24 h further dropped to -30%RC resulting in 38.2% a-helix and 8.1% β-sheet. (Table 8). The exact numbers of structure elements that decrease as the unstructured content increases for both U3.6 wt and its mutant are found to be the same. (Figure S29 and S30). Of the three data sets those down to 195 nm were of the best quality while those down to 200 nm were of the poorest quality. However, in nine instances the predictions for spectrum down to 195 nm were of the worst quality of which four have the highest extent of derandomisation associated with them (Table S3).

CD Spectra tested Against reference Set with the Same Wavelength Range of 240-200 nm

To enable us to answer the question of whether the wavelength range 240-200 nm affects the quality of the spectral and structure predictions. We created spectra with a wavelength down to 200 nm by either removing 10 nm from the original data or 5 nm from the previous reference and test files (with a wavelength down to 195 nm). The reference set was trained with SOMSpec using the same input parameter as specified in section 2.2.5, the test file was tested against the trained map generated by SOMSpec to give the structure estimates and spectral fits.

Best Fit for U3.4 WT and U3.4 R7A Spectra Either Derandomised or Not from 240-200 nm

The best fits of the two initial conditions at which U3.4 WT spectra were collected have the same extent of derandomisation and their structure estimates for all structure types. Subsequently, the extent of derandomisation of the best fit increased to -80%RC at 1 h but later at 24 h dropped to -30%RC after the addition of a buffer. This results in an increase of folded structures (13.5% a-helix and 18.1% β-sheet) within the formulated peptide at 24 h. U3.4 R7A best fits are expected to behave differently regarding the extent of derandomisation associated with them to that of WT under the same conditions. At first, this significantly dropped to a lower fraction (-20%RC) and then remained the same 1 h later and subsequently slightly dropped to -10%RC at 24 h. We observed irrespective of the percentage derandomisation associated with the best fit(s), a-helix content increases over time after the addition of the buffer (e.g., at t = 24 h resulted in 22% a-helix). For cases with -80%RC ether WT or otherwise, their structure estimates are nearly the same (Table 9).

Best Fit for U3.5 wt and U3.5 R7A Spectra Either Derandomised or Not from 240-200 nm

Two of the four conditions for U3.5 wt considered have a higher extent of derandomisation associated with the best fits. As expected, the two notable structures' predictions for these are within a few percentages. There is a hidden challenge in determining the best fit at t = 0 h (immediately after the addition of the buffer to the peptide solution). Thus, two similar lowest NRMSDs were observed, and the best choice was determined to be with -10%RC by visually inspecting the spectra. For the mutant of U3.5 WT (U3.5 R7A), the extent of derandomisation associated with the best fit is the same in water only and immediately after the addition of buffer (-80%RC). Like with U3.5 WT, the notable structures' predictions for these are within a few percentages. Subsequently, as this dropped to -60%RC and further dropped to -40%RC over time the two notable structures' estimates increased to 12.5% and 16.2% respectively, as this corresponds to 135.85% and 82.02% increase from estimates at t = 1 h (Table 10). We found that β-sheet, bend, and turn decrease as the percentage of RC removed from the U3.6K7A spectra increases, while a-helix with -20%RC to 50%RC deviates from this behaviour (Figure S33 and S34).

Best Fit for U3.6 WT and U3.6 R7A Spectra Either Derandomised or Not from 240-200 nm

All the best fits for U3.6 WT explored are associated with the same extent of derandomisation (which is of a higher fraction) and the same structure estimates for each element. For 3.6 K7A, two of the four conditions explored behaviour is similar to that of the WT. They have the same higher extent of derandomisation associated with the best fits and the same structure estimates for each element. Though in water only two similar lowest NRMSDs were observed at first, and then the best choice was determined to be with -0%RC. Later at t = 1 h the associated extent of derandomisation to the best fit markedly increased to -80%RC and then at t = 24 h significantly dropped to -30%RC. The latter case gives rise to 13.5% a-helix and 18.1% β-sheet which are respectively higher than the corresponding estimates for the former at t = 1 h by 221.43% and 260.78% (Table 11). Overall, as the percentage of derandomisation increases there is a gradual decrease in the four structure elements, viz., helix, sheet, bend and turn (Figure S35 and S36).

 

CONCLUSIONS

In this study, we have been able to show that the derandomisation approach previously reported by our group can also be used for peptides, and evaluate how this affects the goodness of fits of the spectral and structure outputs as a function of interaction time with the buffer. The U3.x peptides derandomised are discovered on the skin of toadlet endemic to Australia and in their pristine state they exhibit antimicrobial properties. Almost all the U3.x peptides were straightforward to analyse, except in five instances where visual inspection of the spectra was required to determine the right answers. We found that the extent of derandomisation associated with the best fit could drop or jump or remain the same all through the conditions assessed. In general, it is recognised that an increase or a decrease in the two notable folded structures could have been caused by a drop or a rise of the percentage RC responsible for the right answer, as well as a gradual decrease in four to three structure elements as the percentage RC removed from peptide spectra increases. Our findings on whether the wavelength range affects the quality of spectral and structure predictions extend those of Spencer et al.23, confirming that spectra down to 195 nm are of the best quality while that down to 200 nm are of the poorest quality. We could infer that the structural information is less in the spectra whose wavelength is down to 200 nm, and the quality of the buffer data is only satisfactory down to 195 nm due to the chloride absorbance. There is a need for experiments to prove the underlying approach of derandomisation and its effect on the structure of the proteins or peptides in some processes is real.

 

ACKNOWLEDGEMENTS

The support from the international Macquarie University Research Excellence Scholarship (iMQRES) is gracefully acknowledged.

 

ORCID IDS

Adewale Olamoyesan: https://orcid.org/0000-0002-9879-6726

Alison Rodger: https://orcid.org/0000-0001-8817-4651

 

SUPPLEMENTARY MATERIAL

Supplementary tables and figures are available with this article.

 

REFERENCES

1. Nylén F, Miraglia E, Cederlund A, Ottosson H, Stromberg R, Gudmundsson GH, Agerberth B. Boosting innate immunity: development and validation of a cell-based screening assay to identify LL-37 inducers. Innate Immun. 2014;20(4):364-376. https://doi.org/10.1177/1753425913493338.         [ Links ]

2. Brogden KA. Antimicrobial peptides: pore formers or metabolic inhibitors in bacteria? Nat Rev Microbiol. 2005;3(3):238-250. https://doi.org/10.1038/nrmicro1098.         [ Links ]

3. Apponyi MA, Pukala TL, Brinkworth CS, Maselli VM, Bowie JH, Tyler MJ, Booker GW, Wallace JC, Carver JA, Separovic F, et al. Host-defence peptides of Australian anurans: structure, mechanism of action and evolutionary significance. Peptides. 2004;25(6):1035-1054. https://doi.org/10.1016/j.peptides.2004.03.006.         [ Links ]

4. Bowie JH, Separovic F, Tyler MJ. Host-defense peptides of Australian anurans. Part 2. Structure, activity, mechanism of action, and evolutionary significance. Peptides. 2012;37(1):174-188. https://doi.org/10.1016/j.peptides.2012.06.017.         [ Links ]

5. Pukala TL, Bowie JH, Maselli VM, Musgrave IF, Tyler MJ. Host-defence peptides from the glandular secretions of amphibians: structure and activity. Nat Prod Rep. 2006;23(3):368-393. https://doi.org/10.1039/b512118n.         [ Links ]

6. Bowie JH, Wegner KL, Wabnitz PA, Chia BCS, Carver JA, Wallace JC, Tyler MJ. Host defense antibacterial peptides from the skin secretions of Australian amphibians. Protein Pept Lett. 1999;6:259-269. https://doi.org/10.2174/092986650605221117113359.         [ Links ]

7. Chia BCS, Bowie JH, Carver JA, Mulhern TD. The solution structure of uperin 3.6, an antibiotic peptide from the granular dorsal glands of the Australian toadlet, Uperoleia mjobergii. J Pept Res. 1999;54(2):137-145. https://doi.org/10.1034/j.1399-3011.1999.00095.x.         [ Links ]

8. Bradford AM, Bowie JH, Tyler MJ, Wallace JC. New antibiotic uperin peptides from the dorsal glands of the Australian toadlet Uperoleia mjobergii. Aust J Chem. 1996;49(12):1325-1331. https://doi.org/10.1071/CH9961325.         [ Links ]

9. Calabrese AN, Liu Y, Wang T, Musgrave IF, Pukala TL, Tabor RF, Martin LL, Carver JA, Bowie JH. The amyloid fibril-forming properties of the amphibian antimicrobial peptide uperin 3.5. ChemBioChem. 2016;17(3):239-246. https://doi.org/10.1002/cbic.201500518.         [ Links ]

10. Ray S, Holden S, Prasad AK, Martin LL, Panwar AS. Exploring the role of peptide helical stability in the propensity of uperin 3. x peptides toward beta-aggregation. J Phys Chem B. 2020;124(51):11659-11670. https://doi.org/10.1021/acs.jpcb.0c10000.         [ Links ]

11. Juban MM, Javadpour MM, Barkley MD. Circular dichroism studies of secondary structure of peptides. In: Shafer WM, editor. Antibacterial Peptide Protocols. New Jersey: Humana Totowa; 1997. p. 73-78. https://doi.org/10.1385/0-89603-408-9:73.         [ Links ]

12. Chakrabartty A, Kortemme T, Padmanabhan S, Baldwin RLJB. Aromatic side-chain contribution to far-ultraviolet circular dichroism of helical peptides and its effect on measurement of helix propensities. Biochemistry. 1993;32(21):5560-5565. https://doi.org/10.1021/bi00072a010.         [ Links ]

13. San Miguel M, Marrington R, Rodger PM, Rodger A, Robinson C. An Escherichia coli twin-arginine signal peptide switches between helical and unstructured conformations depending on the hydrophobicity of the environment. Eur J Biochem. 2003;270(16):3345-3352. https://doi.org/10.1046/j.1432-1033.2003.03710.x.         [ Links ]

14. Johnson WC. Analysing protein circular dichroism spectra for accurate secondary structures. Proteins: Struct Funct Bioinf. 1999;35:307-312.         [ Links ]

15. Andrade M, Chacon P, Merelo JJ, Moran F. Evaluation of secondary structure of proteins from UV circular dichroism spectra using an unsupervised learning neural network. Protein Eng Des Sel. 1993;6(4):383-390. https://doi.org/10.1093/protein/6.4.383.         [ Links ]

16. Bulheller BM, Hirst JD. DichroCalc-circular and linear dichroism online. Bioinformatics. 2009;25(4):539-540. https://doi.org/10.1093/bioinformatics/btp016.         [ Links ]

17. Hall V, Nash A, Rodger A. SSNN, a method for neural network protein secondary structure fitting using circular dichroism data. Anal Methods. 2014;6(17):6721-6726. https://doi.org/10.1039/C3AY41831E        [ Links ]

18. Hall V, Sklepari M, Rodger A. Protein secondary structure prediction from circular dichroism spectra using a self-organising map with concentration correction. Chirality. 2014;26(9):471-482. https://doi.org/10.1002/chir.22338.         [ Links ]

19. Olamoyesan A. Protein Structure Fitting by a Self-Organising Map Spectroscopy (SOMSpec): Circular Dichroism and Infrared Spectroscopy; 2022. PhD thesis. Macquarie University, Sydney, Australia. 2022        [ Links ]

20. Olamoyesan A, Ang D, Rodger A. Circular dichroism for secondary structure determination of proteins with unfolded domains using a self-organising map algorithm SOMSpec. RSC Adv. 2021;11(39):23985-23991. https://doi.org/10.1039/D1RA02898G.         [ Links ]

21. Olamoyesan A. Automation of the ad hoc approach for derandomization of proteins: a tutorial for undergraduates in molecular sciences. J Chem Educ. Forthcoming.         [ Links ]

22. Natita W, Wiboonsak W, Dusadee S. Appropriate learning rate and neighborhood function of self-organising map (SOM) for specific humidity pattern classification over Southern Thailand. Int. J. Model. Optim. 2016;6:61-65.         [ Links ]

23. Spencer SE, Rodger A. Bayesian inference assessment of protein secondary structure analysis using circular dichroism data-how much structural information is contained in protein circular dichroism spectra? Anal Methods. 2021;13(3):359-368. https://doi.org/10.1039/D0AY01645D.         [ Links ]

 

 

Received 03 September 2022
Revised 29 December 2022
Accepted 12 February 2023

 

 

* To whom correspondence should be addressed. Email: adewale.olamoyesan@hdr.mq.edu.au; adewale.ola-moyesan@gmail.com

Creative Commons License Todo o conteúdo deste periódico, exceto onde está identificado, está licenciado sob uma Licença Creative Commons