Molecular Dynamics Simulations in Protein-Protein Docking.

Concerted interactions between all the cell components form the basis of biological processes. Protein-protein interactions (PPIs) constitute a tremendous part of this interaction network. Deeper insight into PPIs can help us better understand numerous diseases and lead to the development of new diagnostic and therapeutic strategies. PPI interfaces, until recently, were considered undruggable. However, it is now believed that the interfaces contain "hot spots," which could be targeted by small molecules. Such a strategy would require high-quality structural data of PPIs, which are difficult to obtain experimentally. Therefore, in silico modeling can complement or be an alternative to in vitro approaches. There are several computational methods for analyzing the structural data of the binding partners and modeling of the protein-protein dimer/oligomer structure. The major problem with in silico structure prediction of protein assemblies is obtaining sufficient sampling of protein dynamics. One of the methods that can take protein flexibility and the effects of the environment into account is Molecular Dynamics (MD). While sampling of the whole protein-protein association process with plain MD would be computationally expensive, there are several strategies to harness the method to PPI studies while maintaining reasonable use of resources. This chapter reviews known applications of MD in the PPI investigation workflows.

Magainin 2 and PGLa in bacterial membrane mimics IV: Membrane curvature and partitioning.

We previously reported that the synergistically enhanced antimicrobial activity of magainin 2 (MG2a) and PGLa is related to membrane adhesion and fusion. Here, we demonstrate that equimolar mixtures of MG2a and L18W-PGLa induce positive monolayer curvature stress and sense, at the same time, positive mean and Gaussian bilayer curvatures already at low amounts of bound peptide. The combination of both abilities-membrane curvature sensing and inducing-is most likely the base for the synergistically enhanced peptide activity. In addition, our coarse-grained simulations suggest that fusion stalks are promoted by decreasing the free-energy barrier for their formation rather than by stabilizing their shape. We also interrogated peptide partitioning as a function of lipid and peptide concentration using tryptophan fluorescence spectroscopy and peptide-induced leakage of dyes from lipid vesicles. In agreement with a previous report, we find increased membrane partitioning of L18W-PGLa in the presence of MG2a. However, this effect does not prevail to lipid concentrations higher than 1 mM, above which all peptides associate with the lipid bilayers. This implies that synergistic effects of MG2a and L18W-PGLa in previously reported experiments with lipid concentrations >1 mM are due to peptide-induced membrane remodeling and not their specific membrane partitioning.

Magainin 2 and PGLa in bacterial membrane mimics III: Membrane fusion and disruption.

We previously speculated that the synergistically enhanced antimicrobial activity of Magainin 2 and PGLa is related to membrane adhesion, fusion, and further membrane remodeling. Here we combined computer simulations with time-resolved in vitro fluorescence microscopy, cryoelectron microscopy, and small-angle X-ray scattering to interrogate such morphological and topological changes of vesicles at nanoscopic and microscopic length scales in real time. Coarse-grained simulations revealed formation of an elongated and bent fusion zone between vesicles in the presence of equimolar peptide mixtures. Vesicle adhesion and fusion were observed to occur within a few seconds by cryoelectron microscopy and corroborated by small-angle X-ray scattering measurements. The latter experiments indicated continued and time-extended structural remodeling for individual peptides or chemically linked peptide heterodimers but with different kinetics. Fluorescence microscopy further captured peptide-dependent adhesion, fusion, and occasional bursting of giant unilamellar vesicles a few seconds after peptide addition. The synergistic interactions between the peptides shorten the time response of vesicles and enhance membrane fusogenic and disruption properties of the equimolar mixture compared with the individual peptides.

Enhanced translocation of amphiphilic peptides across membranes by transmembrane proteins.

Cell membranes are phospholipid bilayers with a large number of embedded transmembrane proteins. Some of these proteins, such as scramblases, have properties that facilitate lipid flip-flop from one membrane leaflet to another. Scramblases and similar transmembrane proteins could also affect the translocation of other amphiphilic molecules, including cell-penetrating or antimicrobial peptides. We studied the effect of transmembrane proteins on the translocation of amphiphilic peptides through the membrane. Using two very different models, we consistently demonstrate that transmembrane proteins with a hydrophilic patch enhance the translocation of amphiphilic peptides by stabilizing the peptide in the membrane. Moreover, there is an optimum amphiphilicity because the peptide could become overstabilized in the transmembrane state, in which the peptide-protein dissociation is hampered, limiting the peptide translocation. The presence of scramblases and other proteins with similar properties could be exploited for more efficient transport into cells. The described principles could also be utilized in the design of a drug-delivery system by the addition of a translocation-enhancing peptide that would integrate into the membrane.

Advances in Molecular Understanding of α-Helical Membrane-Active Peptides.

Biological membranes separate the interior of cells or cellular compartments from their outer environments. This barrier function of membranes can be disrupted by membrane-active peptides, some of which can spontaneously penetrate through the membranes or open leaky transmembrane pores. However, the origin of their activity/toxicity is not sufficiently understood for the development of more potent peptides. To this day, there are no design rules that would be generally valid, and the role of individual amino acids tends to be sequence-specific.In this Account, we describe recent progress in understanding the design principles that govern the activity of membrane-active peptides. We focus on α-helical amphiphilic peptides and their ability to (1) translocate across phospholipid bilayers, (2) form transmembrane pores, or (3) act synergistically, i.e., to produce a significantly more potent effect in a mixture than the individual components.We refined the description of peptide translocation using computer simulations and demonstrated the effect of selected residues. Our simulations showed the necessity to explicitly include charged residues in the translocation description to correctly sample the membrane perturbations they can cause. Using this description, we calculated the translocation of helical peptides with and without the kink induced by the proline/glycine residue. The presence of the kink had no effect on the translocation barrier, but it decreased the peptide affinity to the membrane and reduced the peptide stability inside the membrane. Interestingly, the effects were mainly caused by the peptide's increased polarity, not the higher flexibility of the kink.Flexibility plays a crucial role in pore formation and affects distinct pore structures in different ways. The presence of a kink destabilizes barrel-stave pores, because the kink prevents the tight packing of peptides in the bundle, which is characteristic of the barrel-stave structure. In contrast, the kink facilitates the formation of toroidal pores, where the peptides are only loosely arranged and do not need to closely assemble. The exact position of the kink in the sequence further determines the preferred arrangement of peptides in the pore, i.e., an hourglass or U-shaped structure. In addition, we demonstrated that two self-associated (via termini) helical peptides could mimic the behavior of peptides with a helix-kink-helix motif.Finally, we review the recent findings on the peptide synergism of the archetypal mixture of Magainin 2 and PGLa peptides. We focused on a bacterial plasma membrane mimic that contains negatively charged lipids and lipids with negative intrinsic curvature. We showed that the synergistic action of peptides was highly dependent on the lipid composition. When the lipid composition and peptide/lipid ratios were changed, the systems exhibited more complex behavior than just the previously reported pore formation. We observed membrane adhesion, fusion, and even the formation of the sponge phase in this regime. Furthermore, enhanced adhesion/partitioning to the membrane was reported to be caused by lipid-induced peptide aggregation.In conclusion, the provided molecular insight into the complex behavior of membrane-active peptides provides clues for the design and modification of antimicrobial peptides or toxins.

Selecting Collective Variables and Free-Energy Methods for Peptide Translocation across Membranes.

The selective permeability of cellular membranes is a crucial property for controlled transport into and out of cells. Molecules that can bypass the cellular machinery and spontaneously translocate across membranes could be used as therapeutics or drug carriers. Peptides are a prominent class of such molecules, which include natural and man-developed antimicrobial and cell-penetrating peptides. However, the necessary peptide properties for translocation remain elusive. Computer simulations could uncover these properties once we have a good collective variable (CV) that accurately describes the translocation process. Here, we developed a new CV, which includes a description of peptide insertion, local membrane deformation, and peptide internal degrees of freedom related to its charged groups. By comparison of CVs, we demonstrated that all these components are necessary for an accurate description of peptide translocation. Moreover, the advantages and disadvantages of three common methods for free-energy calculations with our CV were evaluated using the MARTINI coarse-grained model: umbrella sampling, umbrella sampling with replica exchange, and metadynamics. The developed CV leads to the reliable and effective calculation of the free energy of peptide translocation, and thus, it could be useful in the design of spontaneously translocating peptides.

Effect of Helical Kink on Peptide Translocation across Phospholipid Membranes.

Biological membranes present a major obstacle for the delivery of therapeutic agents into cells. Some peptides have been shown to translocate across the membrane spontaneously, and they could be thus used as drug-carriers. However, the advantageous peptide properties for the translocation remain unclear. Of particular interest is the effect of a proline-induced kink in α-helical peptides, because the kink was previously reported to both increase and decrease the antimicrobial activity. The antimicrobial activity of peptides could be related to their translocation across the membrane as is the case of the buforin 2 peptide investigated here. Using computer simulations with two independent models, we consistently showed that the presence of the kink has (1) no effect on the translocation barrier, (2) reduces the peptide affinity to the membrane, and (3) disfavors the transmembrane state. Moreover, we were able to determine that these effects are mainly caused by the peptide increased polarity, not the increased flexibility of the kink. The provided molecular understanding can be utilized for the design of cell-penetrating and drug-carrying peptides.

Effect of helical kink in antimicrobial peptides on membrane pore formation.

Every cell is protected by a semipermeable membrane. Peptides with the right properties, for example Antimicrobial peptides (AMPs), can disrupt this protective barrier by formation of leaky pores. Unfortunately, matching peptide properties with their ability to selectively form pores in bacterial membranes remains elusive. In particular, the proline/glycine kink in helical peptides was reported to both increase and decrease antimicrobial activity. We used computer simulations and fluorescence experiments to show that a kink in helices affects the formation of membrane pores by stabilizing toroidal pores but disrupting barrel-stave pores. The position of the proline/glycine kink in the sequence further controls the specific structure of toroidal pore. Moreover, we demonstrate that two helical peptides can form a kink-like connection with similar behavior as one long helical peptide with a kink. The provided molecular-level insight can be utilized for design and modification of pore-forming antibacterial peptides or toxins.

Magainin 2 and PGLa in Bacterial Membrane Mimics II: Membrane Fusion and Sponge Phase Formation.

We studied the synergistic mechanism of equimolar mixtures of magainin 2 (MG2a) and PGLa in phosphatidylethanolamine/phosphatidylglycerol mimics of Gram-negative cytoplasmic membranes. In a preceding article of this series, we reported on the early onset of parallel heterodimer formation of the two antimicrobial peptides already at low concentrations and the resulting defect formation in the membranes. Here, we focus on the structures of the peptide-lipid aggregates occurring in the synergistic regime at elevated peptide concentrations. Using a combination of calorimetric, scattering, electron microscopic, and in silico techniques, we demonstrate that the two peptides, even if applied individually, transform originally large unilamellar vesicles into multilamellar vesicles with a collapsed interbilayer spacing resulting from peptide-induced adhesion. Interestingly, the adhesion does not lead to a peptide-induced lipid separation of charged and charge-neutral species. In addition to this behavior, equimolar mixtures of MG2a and PGLa formed surface-aligned fibril-like structures, which induced adhesion zones between the membranes and the formation of transient fusion stalks in molecular dynamics simulations and a coexisting sponge phase observed by small-angle x-ray scattering. The previously reported increased leakage of lipid vesicles of identical composition in the presence of MG2a/PGLa mixtures is therefore related to a peptide-induced cross-linking of bilayers.

Magainin 2 and PGLa in Bacterial Membrane Mimics I: Peptide-Peptide and Lipid-Peptide Interactions.

We addressed the onset of synergistic activity of the two well-studied antimicrobial peptides magainin 2 (MG2a) and PGLa using lipid-only mimics of Gram-negative cytoplasmic membranes. Specifically, we coupled a joint analysis of small-angle x-ray and neutron scattering experiments on fully hydrated lipid vesicles in the presence of MG2a and L18W-PGLa to all-atom and coarse-grained molecular dynamics simulations. In agreement with previous studies, both peptides, as well as their equimolar mixture, were found to remain upon adsorption in a surface-aligned topology and to induce significant membrane perturbation, as evidenced by membrane thinning and hydrocarbon order parameter changes in the vicinity of the inserted peptide. These effects were particularly pronounced for the so-called synergistic mixture of 1:1 (mol/mol) L18W-PGLa/MG2a and cannot be accounted for by a linear combination of the membrane perturbations of two peptides individually. Our data are consistent with the formation of parallel heterodimers at concentrations below a synergistic increase of dye leakage from vesicles. Our simulations further show that the heterodimers interact via salt bridges and hydrophobic forces, which apparently makes them more stable than putatively formed antiparallel L18W-PGLa and MG2a homodimers. Moreover, dimerization of L18W-PGLa and MG2a leads to a relocation of the peptides within the lipid headgroup region as compared to the individual peptides. The early onset of dimerization of L18W-PGLa and MG2a at low peptide concentrations consequently appears to be key to their synergistic dye-releasing activity from lipid vesicles at high concentrations.

Optimal Hydrophobicity and Reorientation of Amphiphilic Peptides Translocating through Membrane.

Cell-penetrating and some antimicrobial peptides can translocate across lipid bilayers without disrupting the membrane structure. However, the molecular properties required for efficient translocation are not fully understood. We employed the Metropolis Monte Carlo method together with coarse-grained models to systematically investigate free-energy landscapes associated with the translocation of secondary amphiphilic peptides. We studied α-helical peptides with different length, amphiphilicity, and distribution of hydrophobic content and found a common translocation path consisting of adsorption, tilting, and insertion. In the adsorbed state, the peptides are parallel to the membrane plane, whereas, in the inserted state, the peptides are perpendicular to the membrane. Our simulations demonstrate that, for all tested peptides, there is an optimal ratio of hydrophilic/hydrophobic content at which the peptides cross the membrane the easiest. Moreover, we show that the hydrophobicity of peptide termini has an important effect on the translocation barrier. These results provide general guidance to optimize peptides for use as carriers of molecular cargos or as therapeutics themselves.

Synergism of Antimicrobial Frog Peptides Couples to Membrane Intrinsic Curvature Strain.

Mixtures of the frog peptides magainin 2 and PGLa are well-known for their pronounced synergistic killing of Gram-negative bacteria. We aimed to gain insight into the underlying biophysical mechanism by interrogating the permeabilizing efficacies of the peptides as a function of stored membrane curvature strain. For Gram-negative bacterial-inner-membrane mimics, synergism was only observed when the anionic bilayers exhibited significant negative intrinsic curvatures imposed by monounsaturated phosphatidylethanolamine. In contrast, the peptides and their mixtures did not exhibit significant activities in charge-neutral mammalian mimics, including those with negative curvature, which is consistent with the requirement of charge-mediated peptide binding to the membrane. Our experimental findings are supported by computer simulations showing a significant decrease of the peptide-insertion free energy in membranes upon shifting intrinsic curvatures toward more positive values. The physiological relevance of our model studies is corroborated by a remarkable agreement with the peptide's synergistic activity in Escherichia coli. We propose that synergism is related to a lowering of a membrane-curvature-strain-mediated free-energy barrier by PGLa that assists membrane insertion of magainin 2, and not by strict pairwise interactions of the two peptides as suggested previously.

Behavior of BsoBI endonuclease in the presence and absence of DNA.

BsoBI is a type II restriction endonuclease belonging to the EcoRI family. There is only one previously published X-ray structure for this endonuclease: it shows a homodimer of BsoBI completely encircling DNA in a tunnel. In this work, molecular dynamics simulations were employed to elucidate possible ways in which DNA is loaded into this complex prior to its cleavage. We found that the dimer does not open spontaneously when DNA is removed from the complex on the timescale of our simulations (~ 0.5 µs). A biased simulation had to be used to facilitate the opening, which revealed a possible way for the two catalytic domains to separate. The α-helices connecting the catalytic and helical domains were found to act as a hinge during the separation. In addition, we found that the opening of the BsoBI dimer was influenced by the type of counterions present in the environment. A reference simulation of the BsoBI/DNA complex further showed spontaneous reorganization of the active sites due to the binding of solvent ions, which led to an active-site structure consistent with other experimental structures of type II restriction endonucleases determined in the presence of metal ions.

Optimal conditions for opening of membrane pore by amphiphilic peptides.

Amphiphilic peptides can interact with biological membranes and severely affect their barrier and signaling functions. These peptides, including antimicrobial peptides, can self-assemble into transmembrane pores that cause cell death. Despite their medical importance, the conditions required for pore formation remain elusive. Monte Carlo simulations with coarse-grained models enabled us to calculate the free energies of pore opening under various conditions. In agreement with oriented circular dichroism experiments, a high peptide-to-lipid ratio was found to be necessary for spontaneous pore assembly. The peptide length has a non-monotonic impact on pore formation, and the optimal length matches with the membrane thickness. Furthermore, the hydrophobicity of the peptide ends and the mutual positions of peptides on the membrane play a role.