Accelerated Molecular Dynamics Applied to the Peptaibol Folding Problem

The Functionality of Membrane-Inserting Proteins and Peptides: Curvature Sensing, Generation, and Pore Formation

Proteins and peptides with hydrophobic and amphiphilic segments are responsible for many biological functions. The sensing and generation of membrane curvature are the functions of several protein domains or motifs. While some specific membrane proteins play an essential role in controlling the curvature of distinct intracellular membranes, others participate in various cellular processes such as clathrin-mediated endocytosis, where several proteins sort themselves at the neck of the membrane bud. A few membrane-inserting proteins form nanopores that permeate selective ions and water to cross the membrane. In addition, many natural and synthetic small peptides and protein toxins disrupt the membrane by inducing nonspecific pores in the membrane. The pore formation causes cell death through the uncontrolled exchange between interior and exterior cellular contents. In this article, we discuss the insertion depth and orientation of protein/peptide helices, and their role as a sensor and inducer of membrane curvature as well as a pore former in the membrane. We anticipate that this extensive review will assist biophysicists to gain insight into curvature sensing, generation, and pore formation by membrane insertion.

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.

Peptaibiotics: Harnessing the potential of microbial secondary metabolites for mitigation of plant pathogens

Agricultural systems are in need of low-cost, safe antibiotics to protect crops from pests and diseases. Peptaibiotics, a family of linear, membrane-active, amphipathic polypeptides, have been shown to exhibit antibacterial, antifungal, and antiviral activity, and to be inducers of plant resistance against a wide range of phytopathogens. Peptaibiotics belong to the new generation of alternatives to agrochemicals, aligned with the United Nations Sustainable Development Goals and the One Health approach toward ensuring global food security and safety. Despite that, these fungi-derived, non-ribosomal peptides remain surprisingly understudied, especially in agriculture, where only a small number has been tested against a reduced number of phytopathogens. This lack of adoption stems from peptaibiotics' poor water solubility and the difficulty to synthesize and purify them in vitro, which compromises their delivery and inclusion in formulations. In this review, we offer a comprehensive analysis of peptaibiotics' classification, biosynthesis, relevance to plant protection, and mode of action against phytopathogens, along with the techniques enabling researchers to extract, purify, and elucidate their structure, and the databases holding such valuable data. It is also discussed how chemical synthesis and ionic liquids could increase their solubility, how genetic engineering and epigenetics could boost in vitro production, and how omics can reduce screenings' workload through in silico selection of the best candidates. These strategies could turn peptaibiotics into effective, ultra-specific, biodegradable tools for phytopathogen control.

Flexible Gaussian Accelerated Molecular Dynamics to Enhance Biological Sampling

Molecular dynamics simulations often struggle to obtain sufficient sampling to study complex molecular events due to high energy barriers separating the minima of interest. Multiple enhanced sampling techniques have been developed and improved over the years to tackle this issue. Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that works by adding a biasing potential, lifting the energy landscape up, and decreasing the height of its barriers. GaMD allows one to increase the sampling of events of interest without the need of a priori knowledge of the system or the relevant coordinates. All required acceleration parameters can be obtained from a previous search run. Upon its development, several improvements for the methodology have been proposed, among them selective GaMD in which the boosting potential is selectively applied to the region of interest. There are currently four selective GaMD methods that have shown promising results. However, all of these methods are constrained on the number, location, and scenarios in which this selective boosting potential can be applied to ligands, peptides, or protein-protein interactions. In this work, we showcase a GROMOS implementation of the GaMD methodology with a fully flexible selective GaMD approach that allows the user to define, in a straightforward way, multiple boosting potentials for as many regions as desired. We show and analyze the advantages of this flexible selective approach on two previously used test systems, the alanine dipeptide and the chignolin peptide, and extend these examples to study its applicability and potential to study conformational changes of glycans and glycosylated proteins.

Accelerated Molecular Dynamics for Peptide Folding: Benchmarking Different Combinations of Force Fields and Explicit Solvent Models

Accelerated molecular dynamics (aMD) protocols were assessed on predicting the secondary structure of eight peptides, of which two are helical, three are β-hairpins, and three are disordered. Protocols consisted of combinations of three force fields (ff99SB, ff14SB, ff19SB) and two explicit solvation models (TIP3P and OPC), and were evaluated in two independent aMD simulations, one starting from an extended conformation, the other starting from a misfolded conformation. The results of these analyses indicate that all three combinations performed well on helical peptides. As for β-hairpins, ff19SB performed well with both solvation methods, with a slight preference for the TIP3P solvation model, even though performance was dependent on both peptide sequence and initial conformation. The ff19SB/OPC combination had the best performance on intrinsically disordered peptides. In general, ff14SB/TIP3P suffered the strongest helical bias.

Sigmoid Accelerated Molecular Dynamics: An Efficient Enhanced Sampling Method for Biosystems

Gaussian accelerated molecular dynamics (GaMD) is recognized as a popular enhanced sampling method for tackling long-standing challenges in biomolecular simulations. Inspired by GaMD, Sigmoid accelerated molecular dynamics (SaMD) is proposed in this work by adding a Sigmoid boost potential to improve the balance between the highest acceleration and accurate reweighting. Compared with GaMD, SaMD extends the accessible time scale and improves the computational efficiency as tested in three tasks. In the alanine dipeptide task, SaMD can produce the free energy landscape with better accuracy and efficiency. In the chignolin folding task, the estimated Gibbs free energy difference can converge to the experimental value ∼30% faster. In the protein-ligand binding task, the bound conformations are closer to the crystal structure with a minimal ligand root-mean-square deviation of 1.7 Å. The binding of the ligand XK263 to the HIV protease is reproduced by SaMD in ∼60% less simulation time.

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.

Dynorphin A induces membrane permeabilization by formation of proteolipidic pores. Insights from electrophysiology and computational simulations

Dynorphins are endogenous neuropeptides that function as ligands for the κ-opioid receptor. In addition to opioid activity, dynorphins can induce several pathological effects such as neurological dysfunctions and cell death. Previous studies have suggested that Dynorphin A (DynA) mediates some pathogenic actions through formation of transient pores in lipid domains of the plasma membrane. Here, we use planar bilayer electrophysiology to show that DynA induces pore formation in negatively charged membranes. We find a large variability in pore conformations showing equilibrium conductance fluctuations, what disregards electroporation as the dominant mechanism of pore formation. Ion selectivity measurements showing cationic selectivity indicate that positive protein charges of DynA are stabilized by phosphatidyl serine negative charges in the formation of combined structures. We complement our study with computational simulations that assess the stability of diverse peptide arrangements in the hydrophobic core of the bilayer. We show that DynA is capable of assembling in charged membranes to form water-filled pores that conduct ions.

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives

Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Plant antimicrobial peptides: structures, functions, and applications

Antimicrobial peptides (AMPs) are a class of short, usually positively charged polypeptides that exist in humans, animals, and plants. Considering the increasing number of drug-resistant pathogens, the antimicrobial activity of AMPs has attracted much attention. AMPs with broad-spectrum antimicrobial activity against many gram-positive bacteria, gram-negative bacteria, and fungi are an important defensive barrier against pathogens for many organisms. With continuing research, many other physiological functions of plant AMPs have been found in addition to their antimicrobial roles, such as regulating plant growth and development and treating many diseases with high efficacy. The potential applicability of plant AMPs in agricultural production, as food additives and disease treatments, has garnered much interest. This review focuses on the types of plant AMPs, their mechanisms of action, the parameters affecting the antimicrobial activities of AMPs, and their potential applications in agricultural production, the food industry, breeding industry, and medical field.

Fibril fragments from the amyloid core of lysozyme: An accelerated molecular dynamics study

Protein aggregation and formation of amyloid fibrils are associated with many diseases and present a ubiquitous problem in protein science. Hen egg white lysozyme (HEWL) can form fibrils both from the full length protein and from its fragments. In the present study, we simulated unfolding of the amyloidogenic fragment of HEWL encompassing residues 49-101 to study the conformational aspects of amyloidogenesis. The accelerated molecular dynamics approach was used to speed up the sampling of the fragment conformers under enhanced temperature. Analysis of conformational transformation and intermediate structures was performed. During the unfolding, the novel short-living and long-living β-structures are formed along with the unstructured random coils. Such β-structure enriched monomers can interact with each other and propagate into fibril-like forms. The stability of oligomers assembled from these monomers was evaluated in the course of MD simulations with explicit water. The residues playing a key role in fibril stabilization were determined. The work provides new insights into the processes occurring at the early stages of amyloid fibril assembly.

Estimation of pore dimensions in lipid membranes induced by peptides and other biomolecules: A review

The cytoplasmic membrane is one of the most frequent cell targets of antimicrobial peptides (AMPs) and other biomolecules. Understanding the mechanism of action of AMPs at the molecular level is of utmost importance for designing of new membrane-specific molecules. In particular, the formation of pores, the structure and size of these pores are of great interest and require nanoscale resolution approaches, therefore, biophysical strategies are essential to achieve an understanding of these processes at this scale. In the case of membrane active peptides, pore formation or general membrane disruption is usually the last step before cell death, and so, pore size is generally directly associated to pore structure and stability and loss of cellular homeostasis, implicated in overall peptide activity. Up to date, there has not been a critical review discussing the methods that can be used specifically for estimating the pore dimensions induced by membrane active peptides. In this review we discuss the scope, relevance and popularity of the different biophysical techniques such as liposome leakage experiments, advanced microscopy, neutron or X-ray scattering, electrophysiological techniques and molecular dynamics studies, all of them useful for determining pore structure and dimension.