Tripleurin XIIc: Peptide Folding Dynamics in Aqueous and Hydrophobic Environment Mimic Using Accelerated Molecular Dynamics

Modeling the Effect on a Novel Fungal Peptaibol Placed in an All-Atom Bacterial Membrane Mimicking System via Accelerated Molecular Dynamics Simulations

We previously reported on a novel peptaibol, named Tripleurin XIIc (TPN), an 18-residue long sequence produced by the fungus Trichoderma pleuroti. We elucidated its 3D structure via classical and accelerated molecular dynamics simulation (aMD) methods and reported the folding dynamics of TPN in water and chloroform solvents. Peptaibols, in general, are insoluble in water, as they are amphipathic and may prefer hydrophobic environments like transmembrane regions. In this study, we attempted to use aMD simulations to model an all-atom bacterial membrane system while placing a TPN molecule in its vicinity. The results highlighted that TPN was able to introduce some disorder into the membrane and caused lipid clustering. It could also enter the transmembrane region from the water-bilayer interface. The structural dynamics of TPN in the transmembrane region revealed a single energetically stable conformation similar to the one obtained from water and chloroform solvent simulations reported by us previously. However, this linear structure was found to be at the local energy minimum (stable) in water but at a metastable intermediate state (higher energy) in chloroform. Therefore, it could be said that the water solvent can be successfully used for folding simulations of peptaibols.

Unraveling the Mechanism of Action of Myricetin in the Inhibition of hUba1∼Ubiquitin Thioester Bond Formation via In Silico Molecular Modeling Techniques

Ubiquitination is a crucial type of protein modification which helps to control substrate degradation and maintain cell homeostasis. Recent studies suggest that ubiquitination and deubiquitination are involved in regulating metabolic reprogramming in cancer cells and maintaining cancer stem cells. Uba1, a crucial protein in the ubiquitination cascade, can be targeted to develop effective inhibitors for cancer treatment. In previous work, we showed that myricetin (Myr) acts as a potential human Uba1 (hUba1) inhibitor. In this study, we have utilized computational modeling techniques to attempt to illustrate the mechanism of action of Myr. Through extra-precision docking, we confirmed that Myr binds to the adenosine triphosphate (ATP)-binding site of hUba1 (referred to as hotspot 1) with the highest binding affinity. The dynamics of this interaction revealed that hUba1 undergoes a conformational shift from open to closed upon binding of Myr. Myr also migrates outward to interact with the crossover loop simultaneously as the rotational shift of the ubiquitin fold domain (UFD) takes place, thereby blocking access to the ubiquitin binding interface of hUba1 and the crossover loop. The outward migration also explains the reversible nature of Myr binding to hUba1 in previous experiments. We hypothesize that Myr acts as an inhibitor of Uba1∼Ub thioester bond formation by causing a large domain shift toward a closed conformation. Few other analogues of Myr containing the same flavone skeleton showed promising docking scores against hUba1 and could be considered for further validation. We propose that Myr and some of its analogues reported in this study may be promising candidates for developing effective Uba1 inhibitors for cancer treatment.

Structure-activity correlations for peptaibols obtained from clade Longibrachiatum of Trichoderma: A combined experimental and computational approach

Integrated disease management and plant protection have been discussed with much fervor in the past decade due to the rising environmental concerns of using industrially produced pesticides. Members of the genus Trichoderma are a subject of considerable research today due to their several properties as biocontrol agents. In our study, the peptaibol production of Trichoderma longibrachiatum SZMC 1775, T. longibrachiatum f. bissettii SZMC 12546, T. reesei SZMC 22616, T. reesei SZMC 22614, T. saturnisporum SZMC 22606 and T. effusum SZMC 22611 were investigated to elucidate structure-activity relationships (SARs) between the properties of peptaibols and their 3D structures. The effects of peptaibol mixtures obtained from every Trichoderma strain were examined against nine commonly known bacteria. The lowest minimum inhibitory concentrations (MIC, mg ml^-1) were exerted by T. longibrachiatum f. bissettii SZMC 12546 against Gram-positive bacteria, which was also able to inhibit the plant pathogenic Gram-negative Rhizobium radiobacter. Accelerated molecular dynamics (aMD) simulations were performed in aqueous solvent to explore the folding dynamics of 12 selected peptaibol sequences. The most characteristic difference between the peptaibols from group A and B relies in the 'Gly-Leu-Aib-Pro' and 'Gly-Aib-Aib-Pro' motifs ('Aib' stands for α-aminoisobutyric acid), which imparted a significant effect on the folding dynamics in water and might be correlated with their expressed bioactivity. In our aMD simulation experiments, Group A peptaibols showed more restricted folding dynamics with well-folded helical conformations as the most stable representative structures. This structural stability and dynamics may contribute to their bioactivity against the selected bacterial species.

Mode of inhibitory binding of epigallocatechin gallate to the ubiquitin-activating enzyme Uba1 via accelerated molecular dynamics

The green tea polyphenol (−)-epigallocatechin-3-gallate (EGCG) and some of its analogs potently inhibit the ubiquitin-activating enzyme Uba1. In an effort to understand the possible molecular basis of inhibitory activity of EGCG, we conducted a molecular docking and molecular dynamics simulation study. We found that EGCG and its two selected analogs, (−)-epicatechin-3-gallate (ECG) and (−)-epigallocatechin (EGC), bind favorably at two likely hot spots for small-molecule ligand binding on human Uba1. The compounds bind with energetics that mirror their experimental potency for inhibition of Uba1∼ubiquitin thioester formation. The binding of EGCG, ECG, and EGC at one of the hot spots, in particular, recapitulates the rank order of potency determined experimentally and suggests a possible mechanism for inhibition. A hinge-like conformational change of the second catalytic cysteine domain and the opposing ubiquitin-fold domain observed during accelerated molecular dynamics simulations of the EGCG-bound Uba1 complex that results in disruption of the ubiquitin-binding interfaces could explain the compounds' inhibitory activity. These results shed light on the possible molecular mechanism of EGCG and related catechins in the inhibition of Uba1.

The hinge-like movement of the SCCH domain upon ligand binding closes the ubiquitin binding site and disrupts the interfaces crucial for thioester bond formation.

Analysis of mutations of defensin protein using accelerated molecular dynamics simulations

Plant defensins possess diverse biological functions that include antifungal and antibacterial activities and α-amylase and trypsin inhibitory properties. Two mutations, G9R and V39R, were confirmed to increase the antifungal activity of Raphanus sativus antifungal protein 2 (RsAFP2). Accelerated Molecular Dynamics (aMD) were carried out to examine the conformational changes present in these RsAFP2 mutants, and its two closest homologs compared to the wild-type protein. Specifically, the root mean square fluctuation values for the eight cysteine amino acids involved in the four disulfide bonds were low in the V39R mutant compared to the wild-type. Additionally, analysis of the free energy change revealed that G9R and V39R mutations exert a neutral and stabilizing effect on RsAFP2 conformation, and this is supported by the observed lower total energy of mutants compared to the wild-type, suggesting that enhanced stability of the mutants. However, MD simulations to a longer time scale would aid in capturing more conformational state of the wild-type and mutants defensin protein. Furthermore, the aMD simulations on fungal mimic membranes with RsAFP2 and its mutants and homologs showed that the mutant proteins caused higher deformation and water diffusion than the native RsAFP2, especially the V39R mutant. The mutant variants seem to interact by specifically targeting the POPC and POPI lipids amongst others. This work highlights the stabilizing effect of mutations at the 9^th and 39^th positions of RsAFP2 and their increased membrane deformation activity.

Molecular Modeling of Chemosensory Protein 3 from Spodoptera litura and Its Binding Property with Plant Defensive Metabolites

Chemosensory perception in insects involves a broad set of chemosensory proteins (CSPs) that identify the bouquet of chemical compounds present in the external environment and regulate specific behaviors. The current study is focused on the Spodoptera litura (Fabricius) chemosensory-related protein, SlitCSP3, a midgut-expressed CSP, which demonstrates differential gene expression upon different diet intake. There is an intriguing possibility that SlitCSP3 can perceive food-derived chemical signals and modulate insect feeding behavior. We predicted the three-dimensional structure of SlitCSP3 and subsequently performed an accelerated molecular dynamics (aMD) simulation of the best-modeled structure. SlitCSP3 structure has six α-helices arranged as a prism and a hydrophobic binding pocket predominated by leucine and isoleucine. We analyzed the interaction of selected host plant metabolites with the modeled structure of SlitCSP3. Out of two predicted binding pockets in SlitCSP3, the plant-derived defensive metabolites 2-b-D-glucopyranosyloxy-4-hydroxy-7-methoxy-1, 4-benzoxazin-3-one (DIMBOA), 6-Methoxy-2–benzoxazolinone (MBOA), and nicotine were found to interact preferably to the hydrophobic site 1, compared to site 2. The current study provides the potential role of CSPs in recognizing food-derived chemical signals, host-plant specialization, and adaptation to the varied ecosystem. Our work opens new perspectives in designing novel pest-management strategies. It can be further used in the development of CSP-based advanced biosensors.

Changes in Peptaibol Production of Trichoderma Species during In Vitro Antagonistic Interactions with Fungal Plant Pathogens

Trichoderma species are widely used as biofungicides for the control of fungal plant pathogens. Several studies have been performed to identify the main genes and compounds involved in Trichoderma–plant–microbial pathogen cross-talks. However, there is not much information about the exact mechanism of this profitable interaction. Peptaibols secreted mainly by Trichoderma species are linear, 5–20 amino acid residue long, non-ribosomally synthesized peptides rich in α-amino isobutyric acid, which seem to be effective in Trichoderma–plant pathogenic fungus interactions. In the present study, reversed phase (RP) high-performance liquid chromatography (HPLC) coupled with electrospray ionization (ESI) mass spectrometry (MS) was used to detect peptaibol profiles of Trichoderma strains during interactions with fungal plant pathogens. MS investigations of the crude extracts deriving from in vitro confrontations of Trichoderma asperellum and T. longibrachiatum with different plant pathogenic fungi (Fusarium moniliforme, F. culmorum, F. graminearum, F. oxysporum species complex, Alternaria solani and Rhizoctonia solani) were performed to get a better insight into the role of these non-ribosomal antimicrobial peptides. The results revealed an increase in the total amount of peptaibols produced during the interactions, as well as some differences in the peptaibol profiles between the confrontational and control tests. Detection of the expression level of the peptaibol synthetase tex1 by qRT-PCR showed a significant increase in T. asperellum/R. solani interaction in comparison to the control. In conclusion, the interaction with plant pathogens highly influenced the peptaibol production of the examined Trichoderma strains.

Accelerated Molecular Dynamics Applied to the Peptaibol Folding Problem

The use of enhanced sampling molecular dynamics simulations to facilitate the folding of proteins is a relatively new approach which has quickly gained momentum in recent years. Accelerated molecular dynamics (aMD) can elucidate the dynamic path from the unfolded state to the near-native state, “flattened” by introducing a non-negative boost to the potential. Alamethicin F30/3 (Alm F30/3), chosen in this study, belongs to the class of peptaibols that are 7–20 residue long, non-ribosomally synthesized, amphipathic molecules that show interesting membrane perturbing activity. The recent studies undertaken on the Alm molecules and their transmembrane channels have been reviewed. Three consecutive simulations of ~900 ns each were carried out where N-terminal folding could be observed within the first 100 ns, while C-terminal folding could only be achieved almost after 800 ns. It took ~1 μs to attain the near-native conformation with stronger potential boost which may take several μs worth of classical MD to produce the same results. The Alm F30/3 hexamer channel was also simulated in an E. coli mimicking membrane under an external electric field that correlates with previous experiments. It can be concluded that aMD simulation techniques are suited to elucidate peptaibol structures and to understand their folding dynamics.

Structural Diversity and Bioactivities of Peptaibol Compounds From the Longibrachiatum Clade of the Filamentous Fungal Genus Trichoderma

This study examined the structural diversity and bioactivity of peptaibol compounds produced by species from the phylogenetically separated Longibrachiatum Clade of the filamentous fungal genus Trichoderma, which contains several biotechnologically, agriculturally and clinically important species. HPLC-ESI-MS investigations of crude extracts from 17 species of the Longibrachiatum Clade (T. aethiopicum, T. andinense, T. capillare, T. citrinoviride, T. effusum, T. flagellatum, T. ghanense, T. konilangbra, T. longibrachiatum, T. novae-zelandiae, T. pinnatum, T. parareesei, T. pseudokoningii, T. reesei, T. saturnisporum, T. sinensis, and T. orientale) revealed several new and recurrent 20-residue peptaibols related to trichobrachins, paracelsins, suzukacillins, saturnisporins, trichoaureocins, trichocellins, longibrachins, hyporientalins, trichokonins, trilongins, metanicins, trichosporins, gliodeliquescins, alamethicins and hypophellins, as well as eight 19-residue sequences from a new subfamily of peptaibols named brevicelsins. Non-ribosomal peptide synthetase genes were mined from the available genome sequences of the Longibrachiatum Clade. Their annotation and product prediction were performed in silico and revealed full agreement in 11 out of 20 positions regarding the amino acids predicted based on the signature sequences and the detected amino acids incorporated. Molecular dynamics simulations were performed for structural characterization of four selected peptaibol sequences: paracelsins B, H and their 19-residue counterparts brevicelsins I and IV. Loss of position R6 in brevicelsins resulted in smaller helical structures with higher atomic fluctuation for every residue than the structures formed by paracelsins. We observed the formation of highly bent, almost hairpin-like, helical structures throughout the trajectory, along with linear conformation. Bioactivity tests were performed on the purified peptaibol extract of T. reesei on clinically and phytopathologically important filamentous fungi, mammalian cells, and Arabidopsis thaliana seedlings. Porcine kidney cells and boar spermatozoa proved to be sensitive to the purified peptaibol extract. Peptaibol concentrations ≥0.3 mg ml⁻¹ deterred the growth of A. thaliana. However, negative effects to plants were not detected at concentrations below 0.1 mg ml⁻¹, which could still inhibit plant pathogenic filamentous fungi, suggesting that those peptaibols reported here may have applications for plant protection.