In the presence of the other: How glyphosate and peptide molecules alter the dynamics of sorption on goethite

The interactions with soil mineral surfaces are among the factors that determine the mobility and bioavailability of organic contaminants and of nutrients present in dissolved organic matter (DOM) in soil and aquatic environments. While most studies focus on high molar mass organic matter fractions (e.g., humic and fulvic acids), very few studies investigate the impact of DOM constituents in competitive sorption. Here we assess the sorption behavior of a heavily used herbicide (i.e., glyphosate) and a component of DOM (i.e., a peptide) at the water/goethite interface, inclusive of potential glyphosate-peptide interactions. We used in-situ ATR-FTIR (attenuated total reflectance Fourier-transform infrared) spectroscopy to study sorption kinetics and mechanisms of interaction as well as conformational changes to the secondary structure of the peptide. NMR (nuclear magnetic resonance) spectroscopy was used to assess the level of interaction between glyphosate and the peptide and changes to the peptide' secondary structure in solution. For the first time, we illustrate competition for sorption sites results in co-sorption of glyphosate and peptide molecules that affects the extent, kinetics, and mechanism of interaction of each with the surface. In the presence of the peptide, the formation of outer-sphere glyphosate-goethite complexes is favored albeit inner-sphere glyphosate-goethite bonds (i.e., POFe) are still formed. The presence of glyphosate induces secondary structural shifts of the sorbed peptide that maximizes the formation of H-bonds with the goethite surface. However, glyphosate and the peptide do not seem to interact with one another in solution nor at the goethite surface upon sorption. The results of this work highlight potential consequences of competition for sorption sites, for example the transport of organic contaminants and nutrient-rich (i.e., nitrogen) DOM components in relevant environmental systems. Predicting the rate and extent with which organic pollutants are removed from solution by a given solid is also one of the most critical factors for the design of effective sorption systems in engineering applications.

Crystal structures of human neuropeptide Y (NPY) and peptide YY (PYY)

Neuropeptide Y (NPY), peptide YY (PYY) and pancreatic polypeptide (PP) form the evolutionarily conserved pancreatic polypeptide family. While the fold is widely utilized in nature, crystal structures remain elusive, particularly for the human forms, with only the structure of a distant avian form of PP reported. Here we utilize a crystallization chaperone (antibody Fab fragment), specifically recognizing the amidated peptide termini, to solve the structures of human NPY and human PYY. Intriguingly, and despite limited sequence identity (~50%), the structure of human PYY closely resembles that of avian PP, highlighting the broad structural conservation of the fold throughout evolution. Specifically, the PYY structure is characterized by a C-terminal amidated α-helix, preceded by a backfolded poly-proline N-terminus, with the termini in close proximity to each other. In contrast, in the structure of human NPY the N-terminal component is disordered, while the helical component of the peptide is observed in a four-helix bundle type arrangement, consistent with a propensity for multimerization suggested by NMR studies.

Schistocins: Novel antimicrobial peptides encrypted in the Schistosoma mansoni Kunitz Inhibitor SmKI-1

BACKGROUND

Here we describe a new class of cryptides (peptides encrypted within a larger protein) with antimicrobial properties, named schistocins, derived from SmKI-1, a key protein in Shistosoma mansoni survival. This is a multi-functional protein with biotechnological potential usage as a therapeutic molecule in inflammatory diseases and to control schistosomiasis.

METHODS

We used our algorithm enCrypted, to perform an in silico proteolysis of SmKI-1 and a screening for potential antimicrobial activity. The selected peptides were chemically synthesized, tested in vitro and evaluated by both structural (CD, NMR) and biophysical (ITC) studies to access their structure-function relationship.

RESULTS

EnCrypted was capable of predicting AMPs in SmKI-1. Our biophysical analyses described a membrane-induced conformational change from random coil-to-α-helix and a peptide-membrane equilibrium for all schistocins. Our structural data allowed us to suggest a well-known mode of peptide-membrane interaction in which electrostatic attraction between the cationic peptides and anionic membranes results in the bilayer disordering. Moreover, the NMR exchange H/D data with the higher entropic contribution observed for the peptide-membrane interaction showed that shistocins have different orientations upon the membrane.

CONCLUSIONS

This work demonstrate the robustness for using the physicochemical features of predicted peptides in the identification of new bioactive cryptides besides the relevance of combining these analyses with biophysical methods to understand the peptide-membrane affinity and improve further algorithms.

GENERAL SIGNIFICANCE

Bioprospecting cryptides can be conducted through data mining of protein databases demonstrating the success of our strategy. The peptides-based agents derived from SmKI-1 might have high impact for system-biology and biotechnology.

Influence of the conformations of αA-crystallin peptides on the isomerization rates of aspartic acid residues

The isomerization rate of aspartic acid (Asp) residue is known to be affected by the three-dimensional structures of peptides and proteins. Although the isomerized Asp residues were experimentally observed, structural features which affect the isomerization cannot be elucidated sufficiently because of protein denaturation and aggregation. In this study, molecular dynamics (MD) simulations were conducted on three αA-crystallin peptides (T6, T10, and T18), each containing a single Asp residue with different isomerization rate (T18 > T6 > T10) to clarify the structural factors of Asp isomerization tendency. For MD trajectories, distances between side-chain carboxyl carbon of Asp and main-chain amide nitrogen of (n + 1) residue (C_γ-N distances), root mean square fluctuations (RMSFs), and polar surface areas for main-chain amide nitrogen of (n + 1) residues (PSA_N) were calculated, because these structural features are considered to relate to the formations of cyclic imide intermediates. RMSFs and PSA_N are indexes of peptide backbone flexibilities and solvent exposure of the amide nitrogen, respectively. The average C_γ-N distances of T10 was longer than those of the other two peptides. In addition, the peptide containing Asp residue with a higher isomerization rate showed higher flexibility of the peptide backbone around the Asp residue. PSA_N for amide nitrogen in T18 were much larger than those of other two peptides. The computational results suggest that Asp-residue isomerization rates are affected by these factors.

Substance P in the Gas Phase: Conformational Changes and Dissociations Induced by Collisional Activation in a Drift Tube

The work presented below is related to our companion paper in this issue, entitled: Substance P in solution: trans-to-cis configurational changes of penultimate prolines initiate non-enzymatic peptide bond cleavages. Two-dimensional ion mobility spectrometry (IMS-IMS) and mass spectrometry techniques are used to investigate structural transitions for [M+3H]³⁺ ions of substance P (subP) upon collisional activation (CA) in the gas phase. In this approach, different conformations of ions having a specified mobility are selected after an initial IMS separation, collisionally activated to produce new conformers, and these product structures are separated again using a second IMS region. In this way, it is possible to follow folding and unfolding transitions of different conformations. The analysis shows evidence for five conformations. Unlike other systems, every transition is irreversible. Studies as a function of activation voltage are used to discern pathways of structural changes prior to reaching the energy required for dissociation. Thresholds associated with the onsets of transitions are calibrated to obtain estimates of the energetic barriers between different structures and semi-quantitative potential energy diagrams are presented. Overall, barriers associated with structural transitions of [subP+3H]³⁺ in the absence of solvent are on the order of ~ 40 kJ mol^-1, substantially lower than the ~ 90 kJ mol^-1 required for some similar structural transitions in solutions of ethanol. Comparisons of the transition energies in the gas phase with thermochemistry for similar transitions in solution provide clues about why reverse transitions are prohibited. Graphical Abstract.

Substance P in Solution: Trans-to-Cis Configurational Changes of Penultimate Prolines Initiate Non-enzymatic Peptide Bond Cleavages

We report ion mobility spectrometry and mass spectrometry studies of the non-enzymatic step-by-step degradation of substance P (subP), an 11-residue neuropeptide, with the sequence Arg¹-Pro²-Lys³-Pro⁴-Gln⁵-Gln⁶-Phe⁷-Phe⁸-Gly⁹-Leu¹⁰-Met¹¹-NH₂, in ethanol. At elevated solution temperatures (55 to 75 °C), several reactions are observed, including a protonation event, i.e., [subP+2H]²⁺ + H⁺ → [subP+3H]³⁺, that appears to be regulated by a configurational change and two sequential bond cleavages (the Pro²-Lys³ peptide bond is cleaved to form the smaller nonapeptide Lys³-Met¹¹-NH₂ [subP_(3-11)], and subsequently, subP_(3-11) is cleaved at the Pro⁴-Gln⁵ peptide bond to yield the heptapeptide Gln⁵-Met¹¹-NH₂ [subP_(5-11)]). Each of the product peptides [subP_(3-11) and subP_(5-11)] is accompanied by a complementary diketopiperazine (DKP): cyclo-Arg¹-Pro² (cRP) for the first cleavage, and cyclo-Lys³-Pro⁴ (cKP) for the second. Insight about the mechanism of degradation is obtained by comparing kinetics calculations of trial model mechanisms with experimental data. The best model of our experimental data indicates that the initial cleavage of subP is regulated by a conformational change, likely a trans→cis isomerization of the Arg¹-Pro² peptide bond. The subP_(3-11) product has a long lifetime (t_1/2 ~ 30 h at 55 °C) and appears to transition through several structural intermediates prior to dissociation, suggesting that subP_(3-11) is initially formed with a Lys³-trans-Pro⁴ peptide bond configuration and that slow trans→cis isomerization regulates the second bond cleavage event as well. From these data and our model mechanisms, we obtain transition state thermochemistry ranging from ΔH^‡ = 41 to 85 kJ mol^-1 and ΔS^‡ = - 43 to - 157 J mol^-1 K^-1 for each step in the reaction. Graphical Abstract.

Hydrogen-Deuterium Exchange and Electron Capture Dissociation to Interrogate the Conformation of Gaseous Melittin Ions

There is a need in the field of biological mass spectrometry for structural tools which can report on regional, rather than solely global, structure of gaseous protein ions. Site-specific hydrogen-deuterium (H/D) exchange has shown promise in fulfilling this need, but requires additional method development to prove its utility. In this study, we use H/D exchange and electron capture dissociation (ECD) to probe the gaseous structure of two peptides which are α-helical in solution and which differ by a single point mutation. Global H/D exchange levels, ECD fragmentation profiles, and region specific H/D exchange profiles are compared between wild type (WT) melittin, which adopts a hinged helix conformation in solution, and a mutant P14A melittin which folds into a single helix in solution. High protection from H/D exchange by both peptides is consistent with retention of secondary structure in the gas phase (or refolding into some other compact structure). The P14A mutant melittin exhibits lower ECD fragmentation efficiency than WT melittin, suggesting that it contains more secondary structure in the gas phase, which may indicate that these peptides retain some memory of their solution-phase structures. Examination of the isotopic distributions of fragment ions derived from H/D exchange with subsequent ECD reveals that the C-terminus of these peptides adopts multiple conformations. The results reported here offer insight into the stability of alpha helices in the gas phase, and also highlight the value of combining gas-phase H/D exchange with electron capture dissociation to interrogate gaseous peptide conformation.

Conformational perturbation of peptides in presence of polar organic solvents

The critical role played by solvent environment in maintaining the conformational integrity of peptides and proteins is accepted without question. Numerous experiments have suggested that perturbing the solvent environment of peptides and proteins by the addition of polar organic solvents have important consequences for the conformation of these molecules. However, experimental studies of such perturbations often report different kinds of effects depending on the solvent used and/or the sequence/structure of the molecule under study. In this work we report a simulation based comparative study on the effects of adding two common organic solvents viz. Dimethyl sulfoxide (DMSO) and Acetonitrile (MeCN) on the dynamical conformation of a test peptide Ace-Gly-X-Gly-Nme where X is any amino acid. Our studies identify important differences in peptide solvation by these two solvents, which we attempt to correlate with the kinetic stability of the conformation, as well as the identity of the central 'X' residue in the test peptide.

Structure-based design for binding peptides in anti-cancer therapy

The conventional anticancer therapeutics usually lack cancer specificity, leading to damage of normal tissues that patients find hard to tolerate. Ideally, anticancer therapeutics carrying payloads of drugs equipped with cancer targeting peptides can act like "guided missiles" with the capacity of targeted delivery toward many types of cancers. Peptides are amenable for conjugation to nano drugs for functionalization, thereby improving drug delivery and cellular uptake in cancer-targeting therapies. Peptide drugs are often more difficult to design through molecular docking and in silico analysis than small molecules, because peptide structures are more flexible, possess intricate molecular conformations, and undergo complex interactions. In this review, the development and application of strategies for structure-based design of cancer-targeting peptides against GRP78 are discussed. This Review also covers topics related to peptide pharmacokinetics and targeting delivery, including molecular docking studies, features that provide advantages for in vivo use, and properties that influence the cancer-targeting ability. Some advanced technologies and special peptides that can overcome the pharmacokinetic challenges have also been included.

In Situ Investigation of Peptide-Lipid Interaction Between PAP248-286 and Model Cell Membranes

Sum frequency generation vibrational spectroscopy (SFG) was utilized to investigate the interaction between PAP248-286 and the two lipid bilayer systems. The present study also provides spectroscopic evidence to confirm that, although PAP248-286 is unable to penetrate into the hydrophobic core of the lipid bilayers, it is capable of interacting more intimately with the fluid-phase POPG/POPC than with the gel-phase DPPG/DPPC lipid bilayer. The helical structure content of lipid-bound PAP248-286 was also observed to be high, in contrast to the results previously reported using nuclear magnetic resonance (NMR). Collectively, our SFG data suggest that lipid-bound PAP248-286 actually resembles its structure in 50 % 2,2,2-trifluoroethanol better than the structure when the peptide binds to SDS micelles. This present study questions the use of SDS micelles as the model membrane for NMR studies of PAP248-286 due to its protein denaturing activity.

In silico analysis and experimental verification of OSR1 kinase - Peptide interaction

The oxidative-stress-responsive kinase 1 (OSR1) and the STE20/SPS1-related proline/alanine-rich kinase (SPAK) are key enzymes in a signaling cascade regulating the activity of Na(+)-K(+)-2Cl(-) cotransporters (NKCC1-2) and Na(+)-Cl(-) cotransporter (NCC). Both kinases have a conserved carboxyl-terminal (CCT) domain, which recognizes a unique peptide motif present in OSR1- and SPAK-activating kinases (with-no-lysine kinase 1 (WNK1) and WNK4) as well as their substrates (NKCC1, NKCC2, and NCC). Utilizing various modalities of the Rosetta Molecular Modeling Software Suite including flexible peptide docking and protein design, we comprehensively explored the sequence space recognized by the CCT domain. Specifically, we studied single residue mutations as well as complete unbiased designs of a hexapeptide substrate. The computational study started from a crystal structure of the CCT domain of OSR1 in complex with a hexapeptide derived from WNK4. Point mutations predicted to be favorable include Arg to His or Trp substitutions at position 2 and a Phe to Tyr substitution at position 3 of the hexapeptide. In addition, de novo design yielded two peptides predicted to bind to the CCT domain: FRFQVT and TRFDVT. These results, which indicate a little bit more freedom in the composition of the peptide, were confirmed through the use of yeast two-hybrid screening.

Proton-detected C,H correlation NMR techniques for the complete assignment of all proton and carbon resonances of a cyclic peptide

The assignment of the proton and carbon spectra of the cyclic peptide cyclo(-Phe¹ -Pro² -Thr³ -Lys(Z)⁴ -Trp⁵ -Phe⁶ -) was accomplished by the application of multiple quantum proton-detected heteronuclear correlation spectroscopy. Since the proton spectrum shows severe overlap, the carbon chemical shifts were used to disentangle the proton resonances. The methodology described is useful even in cases where only limited quantities of materials are available. The combination of (i) a proton decoupled C,H correlation spectrum and (ii) a (not decoupled) H-relayed C,H correlation gave the assignments of all CH_n groups. The non-protonated carbons, i.e. the carbonyl carbons of the peptide bond, were assigned with a C,H correlation optimized for long-range couplings, an experiment that also gave helpful information about the conformational features of the hexapeptide. The cyclic peptide contains a Phe-Pro cis-peptide bond forming a βVI-like bend and a β-turn about the amino acids Thr-Lys-Trp-Phe. Although the conventional discussion of NMR parameters indicates a strong preference for one conformation, the quantitative evaluation of NOE-derived distances in restrained MD calculations proves that the type of β-turn in the last-mentioned region is not unique. Whereas the MD calculations converge to a βII' turn, the vicinal proton coupling constants are in better agreement with type βI. Thus a dynamic equilibrium of the backbone is proposed.