If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
The octapeptins are lipopeptide antibiotics that are structurally similar to polymyxins yet retain activity against polymyxin-resistant Gram-negative pathogens, suggesting they might be used to treat recalcitrant infections. However, the basis of their unique activity is unclear because of the difficulty in generating high-resolution experimental data of the interaction of antimicrobial peptides with lipid membranes. To elucidate these structure–activity relationships, we employed all-atom molecular dynamics simulations with umbrella sampling to investigate the conformational and energetic landscape of octapeptins interacting with bacterial outer membrane (OM). Specifically, we examined the interaction of octapeptin C4 and FADDI-115, lacking a single hydroxyl group compared with octapeptin C4, with the lipid A–phosphoethanolamine modified OM of Acinetobacter baumannii. Octapeptin C4 and FADDI-115 both penetrated into the OM hydrophobic center but experienced different conformational transitions from an unfolded to a folded state that was highly dependent on the structural flexibility of their respective N-terminal fatty acyl groups. The additional hydroxyl group present in the fatty acyl group of octapeptin C4 resulted in the molecule becoming trapped in a semifolded state, leading to a higher free energy barrier for OM penetration. The free energy barrier for the translocation through the OM hydrophobic layer was ∼72 kcal/mol for octapeptin C4 and 62 kcal/mol for FADDI-115. Our results help to explain the lower antimicrobial activity previously observed for octapeptin C4 compared with FADDI-115 and more broadly improve our understanding of the structure–function relationships of octapeptins. These findings may facilitate the discovery of next-generation octapeptins against polymyxin-resistant Gram-negative ‘superbugs.‘
). Of particular concern are the Gram-negative ‘superbugs,‘ including Acinetobacter baumannii, that show resistance to almost all currently available antibiotics (
). Given that the development of new antibiotics has stagnated in recent years, polymyxins (i.e. polymyxin B and colistin) are increasingly used as a last-line therapy against these multidrug-resistant Gram-negative pathogens (
). Unfortunately, polymyxin resistance is increasing, such as the emergence and spread of plasmid-borne mcr genes that encode phosphoethanolamine transferases for the modification of lipid A (
Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study.
). Although structurally similar to the polymyxins (Table S1), there are major differences in the antibacterial activity between the two classes of lipopeptides. Whereas the polymyxins are only active against Gram-negative bacteria (
). Importantly, the octapeptins do not exhibit cross-resistance with polymyxins and retain activity against polymyxin-resistant Gram-negative bacteria (
Given the structural similarities to the polymyxins, the octapeptins are believed to have a similar mode of action that involves insertion into the OM of Gram-negative bacteria followed by membrane disorganization, lysis, and eventual cell death (
). Resistance to the polymyxins generally develops through modification of the lipid A component of lipopolysaccharide in the OM with cationic moieties such as ethanolamine and 4′-aminoarabinose (
). Interestingly, lipid A–based NMR models revealed that the octapeptins can specifically bind to the ethanolamine- or 4′-aminoarabinose–modified lipid A molecule through strong hydrophobic contacts (
). Minor structural modifications to the octapeptin core scaffold not involving the α-γ-diaminobutyric acid (Dab) residues have also been shown to dramatically hinder their bactericidal activity (
). These findings suggest that the activity of the octapeptins depends upon their unique interaction pattern with the OM rather than simply acting as a cationic detergent. However, the interaction between the octapeptins and the bacterial OM has never been investigated at the atomic level, significantly limiting our understanding of their mechanism of activity and hindering the discovery of novel octapeptins.
In the present study, an asymmetric lipid A–phosphoethanolamine (lipid A–pEtN; Fig. S1)–modified OM model was constructed based on our quantitative lipidomics results of polymyxin-resistant A. baumannii (Fig. 1B) (
). For the first time, we employed all-atom molecular dynamics (MD) simulations and umbrella sampling technique to investigate the conformation and energetics governing the penetration of two octapeptins (octapeptin C4 and FADDI-115) into the lipid A–pEtN modified OM. Despite their close structural similarity, significant changes in their free energy profiles and conformational transitions during OM penetration were observed, revealing new insight into the structure–interaction relationship of octapeptins.
Figure 1Chemical structures and description of molecular dynamics simulation systems.A, chemical structures of octapeptin C4 and the octapeptin analog FADDI-115. The additional hydroxyl group in octapeptin C4 is shown in red. B, the equilibrated simulation box used for molecular dynamics simulations. The water shell is shown as transparent surface model, the octapeptin molecule is displayed as blue spheres, and the bacterial outer membrane is shown as gray spheres with the phosphate groups in pink. The lipid composition of the outer membrane is shown for each membrane leaflet. Lipid A–pEtN, phosphoethanolamine-modified lipid A; PE, phosphatidylethanolamine; PG, phosphatidylglycerol; CL, cardiolipin.
Interaction of octapeptin C4 and FADDI-115 with the bacterial OM
The motion trajectories of the octapeptin molecules in the steered simulations were tracked relative to the center of the OM (Fig. 2A). After ∼5 ns, both the octapeptin C4 and FADDI-115 molecules approached the surface of the outer leaflet of the OM, traversing the polar headgroup region after ∼12 ns. After a short plateau phase (∼2 ns for octapeptin C4 and 4 ns for FADDI-115), the octapeptin molecule completely penetrated into the hydrophobic region of the OM after roughly 36 ns. We investigated the interactions between the octapeptin molecules and the OM at different stages of penetration. At the surface-binding stage, the positively charged side chains of d-Dab1 and Dab6 of octapeptin C4 closely interacted with the negatively charged phosphate groups of adjacent lipid A–pEtN molecules (Fig. 2B, panel i). Similarly, the side chains of d-Dab1 and Dab3 in FADDI-115 bound to the phosphate groups of lipid A–pEtN molecules (Fig. 2C, panel i). These results indicate that the initial binding of the octapeptins to the bacterial OM was mainly driven by the long-range electrostatic interactions. Subsequently, both octapeptin C4 and FADDI-115 inserted their fatty acyl groups into the hydrocarbon region of the OM, although their cyclic heptapeptide segments formed multiple hydrogen bonds with the headgroups of lipid A–pEtN molecules (Fig. 2, B and C, panels ii). After this, both octapeptins formed a folded conformation within the hydrophobic center of the OM. For octapeptin C4, its N-terminal fatty acyl group and the d-Phe4 oriented toward the OM center and formed hydrophobic contacts with the hydrocarbon chains of lipid A–pEtN and phospholipid molecules, whereas the polar heptapeptide ring made intensive interactions with the phosphate groups of lipid A–pEtN molecules (Fig. 2B, panel iii). The interaction was somewhat different for FADDI-115; the fatty acyl chain and Leu8 were embedded in the hydrocarbon region of the OM, whereas the cyclic head interacted with the phosphate groups of lipid A–pEtN molecules (Fig. 2C, panel iii). During the OM penetration process, the positively charged Dab residues of both octapeptins formed strong interactions with the headgroups of lipid A–pEtN molecules, which triggered ∼150° rotation of the headgroups of the adjacent lipid A–pEtN molecules (Fig. 2D) and the local structural disorganization of the OM (Fig. S2).
Figure 2Octapeptin interaction with the bacterial outer membrane.A, distance of the center of mass of the octapeptin molecule from the center of the outer membrane. The polar headgroup region of the outer membrane is labeled by the two orange dashed lines. B and C show the structural snapshots of the penetration of octapeptin C4 and FADDI-115 in different regions of the outer membrane, respectively. The octapeptin molecule is shown as blue sticks, the lipid A–pEtN molecules are shown as teal lines, and the phospholipid molecules are shown as gray lines. The hydrogen bonds between the octapeptin molecule and the lipid A–pEtN molecules in the outer membrane are shown as purple dashed lines. D, the orientation of lipid A–pEtN molecule relative to the normal of the outer membrane.
Before the penetration of the octapeptins, water molecules could only reach the polar headgroup region of the OM and were largely excluded from the hydrophobic region. However, following the penetration of the octapeptin molecule into the center of the OM, a large number of water molecules entered into the hydrophobic layer of the OM through the penetration site of the octapeptin molecules; the number of water molecules within the membrane–water interface increased from ∼350 to 650 in the octapeptin C4 system (Fig. 3A) and to 550 in the FADDI-115 system (Fig. 3B). This increase in the water permeability of the bacterial OM promoted by the octapeptins may lead to further destabilization of the bacterial OM and help explain the antimicrobial activity of the octapeptins.
Figure 3Water permeation of the outer membrane caused by octapeptin C4 (A) and FADDI-115 (B).Z = 0 indicates the membrane center. The structural snapshots show the different stages of the penetration of octapeptins into the outer membrane. The octapeptin molecule is shown as blue sticks. The water molecules are shown as spheres with red oxygen atoms and white hydrogen atoms. The reaction coordinates along z axis represent the positions relative to the center of the outer membrane.
). Through a series of umbrella sampling simulations, we characterized the free energy profiles for the octapeptins moving from the bulk water to the hydrophobic center of the OM (Fig. 4). The whole free energy profile was divided into three components. The ΔGbind describes the binding of the octapeptin to the surface of the OM outer leaflet (2.5 nm < Z < 3.0 nm), whereas the ΔGtranslocate-1 and ΔGtranslocate-2 values describe the translocations of the octapeptin molecule through the headgroup region (1.8 nm < Z < 2.5 nm) and the hydrocarbon region (0 nm < Z < 1.8 nm) of the OM outer leaflet, respectively. The ΔGbind was approximately −3.5 kcal/mol (favorable) for both octapeptin C4 and FADDI-115. The ΔGtranslocate-1 was comparatively unfavorable (∼8 kcal/mol) because the octapeptin C4 and FADDI-115 passed through the headgroup region of the OM. Curiously, the free energy profile diverged when traversing the hydrocarbon region of the OM for each octapeptin, especially in the range Z < 1.2 nm. To pass through this hydrophobic region, the ΔGtranslocate-2 was ∼72 kcal/mol for octapeptin C4 and 62 kcal/mol for FADDI-115 (energy-unfavorable). Interestingly, the free energy profiles of octapeptin C4 and FADDI-115 significantly diverged when they penetrated through the inner leaflet of the OM. The free energy curve reached the plateau at Z = −1 nm in the system with FADDI-115, and the free energy for FADDI-115 was ∼50 kcal/mol lower than that for octapeptin C4 at Z = −3 nm (Fig. S3). Because the only structural difference between these two octapeptins is an additional hydroxyl group on the fatty acyl chain of octapeptin C4 (Fig. 1A), these results show that this structural variation did not affect their binding affinity with the bacterial OM and penetrating ability through the polar headgroup region but significantly changed their ability to traverse the hydrophobic region of the OM.
Figure 4Free energy profiles for the penetration of octapeptin C4 and FADDI-115 into the outer membrane. The free energy barriers for octapeptin binding to the outer membrane, passing through the headgroup and hydrocarbon regions of the outer leaflet are depicted by ΔGbind, ΔGtranslocate-1, and ΔGtranslocate-2, respectively. The reaction coordinates along z axis represent the positions relative to the center of the outer membrane, with Z = 0 indicating the hydrophobic center of the outer membrane.
Conformational transition of octapeptin C4 and FADDI-11 in the bacterial OM
As shown in Fig. 2 (B and C, panels iii), both octapeptin C4 and FADDI-115 formed a folded conformation in the hydrophobic center of the OM, consistent with the previous NMR findings (
). To examine the mechanism of the conformational transition, we characterized the conformational ensembles of the octapeptin C4 and FADDI-115 molecules by measuring their radius of gyration (Rg) around the z axis during the penetration process. The octapeptin C4 molecule displayed three major conformational clusters during the entire penetration process (Fig. 5A). It adopted an unfolded state (Rg > 0.6 nm) in the polar headgroup region of the OM (Z > 1.8 nm), transformed into an intermediate, semifolded state (Rg = ∼0.54 nm) in the headgroup–hydrocarbon interface region of the OM (1.0 nm < Z < 1.5 nm), and finally stayed in a more compact, folded state (Rg = ∼0.45 nm) near the OM hydrophobic center (Z < 0.5 nm). Similarly, the FADDI-115 molecule also maintained an unfolded state in the polar headgroup region of the OM and formed a folded state near the OM center (Fig. 5B). However, the transition between these two states was achieved through an intermediate aisle (0.5 nm < Z < 1 nm), rather than the intermediate conformational cluster observed in the octapeptin C4 system.
Figure 5Conformational dynamics of octapeptin C4 (A) and FADDI-115 (B) in the outer membrane, with the superimposition of the octapeptin molecule at different conformational states shown for octapeptin C4 (C) and FADDI-115 (D). The conformational landscape is defined by sampling the radius of gyration of the octapeptin molecule around the z axis and the reaction coordinates of the penetration. The conformations of the octapeptin molecule are represented by points, and the color spectrum represents the density of points. The pathway of conformational transition from the unfolded state (US) to intermediate state (IM) and folded state (FS) is shown by a black line.
Interestingly, structural superimposition showed that the cyclic heptapeptide component of the octapeptin C4 and FADDI-115 molecules was well-aligned among their unfolded, intermediate, and folded states. The variation in the radius of gyration observed was due to conformational swing of their respective fatty acyl groups (Fig. 5, C and D). This may suggest that the hydroxyl group on the fatty acyl group of octapeptin C4 is a key factor for their observed different conformational transitions in the OM.
Structure–interaction relationships of octapeptins
The all-atom simulation results were utilized to examine the interactions between the OM and the octapeptin molecule in its folded state. A single octapeptin C4 molecule was observed to bind to three lipid A–pEtN molecules, and the complex was stabilized by a combination of electrostatic and hydrophobic interactions (Fig. 6A). The positively charged side chains of Dab7 interacted with the negatively charged phosphate group of the lipid A–pEtN molecule; the side chain of d-Dab1 made a hydrogen bond with the sugar of the lipid A–pEtN molecule; the main chain of the residues on the cyclic segment contributed five hydrogen bonds with the lipid A–pEtN molecules, and the N-terminal fatty acyl group and d-Phe4 formed hydrophobic contacts with the hydrocarbon chains of the adjacent lipid A–pEtN molecules. Notably, the side chains of Dab3 and Dab6 made few interactions with the lipid A–pEtN molecules, likely because of the shielding from the hydrophobic d-Phe4, Leu5 and Leu8. Similarly, the FADDI-115 molecule bound to three lipid A–pEtN molecules (Fig. 6B). The side chains of d-Dab1, Dab3, and Dab7 interacted with the phosphate groups of the lipid A–pEtN molecules; the main chains of Dab3 and Dab6 made three hydrogen bonds with the lipid A–pEtN molecules; and the N-terminal fatty acyl group and Leu8 formed hydrophobic contacts with the hydrocarbon tails of the adjacent lipid A–pEtN molecules.
Figure 6Interactions between the octapeptin and lipid A–pEtN molecules in the outer membrane.A, octapeptin C4. B, FADDI-115. C, the hydrogen bonds formed between the hydroxyl group on the fatty acyl chain of octapeptin C4 and the outer membrane lipids.
Additionally, we examined the interactions associated with the hydroxyl group on the fatty acyl group of octapeptin C4 (Fig. 6C). The hydroxyl group formed up to two hydrogen bonds with the OM lipids during the entire penetration process and likely contributed two hydrogen bonds when the octapeptin C4 molecule reached the headgroup–hydrocarbon interface region of the OM (1.0 nm < Z < 1.5 nm). These intermolecular hydrogen bonds with octapeptin C4 appear to change the conformational dynamics of the fatty acyl group, which stabilized a large intermediate, semifolded conformational cluster at this interface region (Fig. 5A).
Discussion
The octapeptins are structurally similar to the polymyxins, with both possessing an N-terminal fatty acyl group and three cationic Dab residues located at the same positions in a heptapeptide ring. However, despite these similarities, the octapeptins retain activity against many polymyxin-resistant Gram-negative bacteria (
). As a class of membrane-targeting antibiotics, the interaction between the octapeptins and the bacterial OM has not been investigated at the atomic level to date. In this study, all-atom MD simulations were employed to investigate how the octapeptins interact with the lipid A–pEtN modified OM that renders many Gram-negative bacteria polymyxin-resistant. Unlike previous simulation studies in which the bacterial OM was simplified into a symmetric or asymmetric lipopolysaccharide bilayer with model phospholipid compositions (
Polymyxin binding to the bacterial outer membrane reveals cation displacement and increasing membrane curvature in susceptible but not in resistant lipopolysaccharide chemotypes.
J. Chem. Inf. Model.2017; 57 (28805387): 2181-2193
), we utilized the quantitative membrane lipidomics results from polymyxin-resistant A. baumannii 5075R to construct a more realistic OM model that allowed for better precision when characterizing interactions with the octapeptins (
). Analysis of the simulations revealed that a minor structural variation between octapeptin C4 and FADDI-115, namely, the presence of a hydroxyl group on the fatty acyl group of octapeptin C4 (Fig. 1A), was sufficient to impact their conformational transitions and resulted in different free energy profiles for OM penetration. These results provided novel atomic-scale insights into the structure–activity relationship of the octapeptins against polymyxin-resistant Gram-negative bacteria.
Free energy profiles quantitatively evaluated the energetics governing the penetration of the octapeptin molecules into the bacterial OM (Fig. 4). The free energy barrier for the translocation through the hydrophobic layer of the OM (ΔGtranslocate-2 = ∼72 and 62 kcal/mol for octapeptin C4 and FADDI-115, respectively; Fig. 4) was much lower than that of colistin A (∼82 kcal/mol) (
Outer membranes of polymyxin-resistant Acinetobacter baumannii with phosphoethanolamine-modified lipid a and lipopolysaccharide loss display different atomic-scale interactions with polymyxins.
). The high free energy barrier in this region suggests that the hydrophobic layer of the OM constitutes the major defense to the action of both octapeptins and polymyxins. Lacking the hydroxyl group found on the fatty acyl tail of octapeptin C4 (Fig. 1A), FADDI-115 is slightly more hydrophobic than octapeptin C4. The lower free energy required for FADDI-115 to penetrate the membrane indicates that the additional hydroxyl group of octapeptin C4 negatively impacted OM penetration. This may explain the greater bacterial killing of FADDI-115 compared with octapeptin C4 against polymyxin-resistant A. baumannii strains (
). Importantly, the lower free energy barrier for the octapeptins compared with the polymyxins illustrates that their superior antimicrobial activity against polymyxin-resistant Gram-negative bacteria is very likely caused by their enhanced ability to traverse the hydrophobic region of the OM.
Our conformational analysis suggests that free energy differences between octapeptin C4 and FADDI-115 were linked to their different conformational dynamics within the OM. Specifically, a large conformational cluster populating an intermediate state was observed for the octapeptin C4 at the headgroup–hydrocarbon interface region of the OM and not for FADDI-115 (Fig. 5). This may indicate that during penetration into the membrane, octapeptin C4 is readily trapped in its intermediate state, making the transition to the folded state more difficult to achieve compared with FADDI-115. That octapeptin C4 formed the intermediate conformation close to the headgroup–hydrocarbon interface region of the OM (the region where the free energy profiles of octapeptin C4 and FADDI-115 diverged) suggests that the formation of the folded conformation plays a key role in regulating their penetration into the bacterial OM. With the folded conformation, the hydrophobic residues and charged Dab residues were separated into two domains, conferring structural amphipathicity (Fig. 6). This special structural arrangement suits the interactions with both the hydrocarbon tails and phosphate groups of lipid A–pEtN molecules.
A similar folded conformation was also found in an NMR-based interaction model of octapeptin C4 with a single lipid A molecule (
). Importantly, the specific interaction pairs derived from the NMR-based model and our simulation-based model are quite similar. In both models, the fatty acyl group and d-Phe4 formed hydrophobic contacts with the carbon tails of lipid A; the d-Dab1 residue interacted with the sugar of lipid A; and Dab3, Dab6, and Dab8 residues interacted with the phosphate groups of the lipid A (Fig. 6A). However, the octapeptin C4 molecule adopted a more compact conformation in the NMR-based model (
) compared with our simulation-based model, very likely because the octapeptin C4 interacted with only a single lipid A molecule in the NMR model but with three lipid A molecules in the membrane condition. Moreover, our recent studies discovered that the polymyxins also adopted a similar folded conformation in WT lipid A OM but not in lipid A–pEtN modified OM (
Outer membranes of polymyxin-resistant Acinetobacter baumannii with phosphoethanolamine-modified lipid a and lipopolysaccharide loss display different atomic-scale interactions with polymyxins.
). Collectively, the reported NMR study and our simulation results reveal that the folding is critical for the OM penetration of both octapeptins and polymyxins. Once in the folded conformation, octapeptin and polymyxin molecules are able to fully penetrate the membrane, which leads to membrane destabilization. Given the bactericidal activity of the polymyxins involves the disorganization of the OM (
), the ability of the octapeptins (but not the polymyxins) to form the folded conformation in the OM of polymyxin-resistant bacteria may explain their activity against these polymyxin-resistant organisms. Furthermore, we also found that the octapeptins formed the folded conformation within the WT, nonmodified lipid A OM of A. baumannii and interacted with the surrounding lipid A molecules via both electrostatic interaction (between Dab residues and phosphate groups of lipid A) and hydrophobic interaction (between hydrophobic regions of octapeptin and hydrocarbon tails of lipid A) (Fig. S4). These results support the similar activity of octapeptins against polymyxin-susceptible and polymyxin-resistant Gram-negative bacteria (
In the interaction model of the two octapeptins (Fig. 6), both the hydrophobic interactions (involving the fatty acyl group and d-Phe4 in octapeptin C4; or fatty acyl group and Leu8 in FADDI-115) and polar interactions (involving the Dab residues) stabilized the interaction complex with multiple lipid A–pEtN molecules. Our previous NMR results emphasized the importance of the hydrophobic interactions for octapeptin activity (
). In our simulations, we found that the main chain of the heptapeptide ring of octapeptins formed three to five hydrogen bonds with the phosphate groups of lipid A–pEtN molecules (Fig. 6). In contrast, although similar hydrophobic interactions were still present with polymyxin B, only the side chains of its Dab residues interacted with the phosphate groups of lipid A (
Molecular dynamics simulations informed by membrane lipidomics reveal structure–interaction relationship of polymyxins with lipid A–based outer membrane of Acinetobacter baumannii.
J. Antimicrob. Chemother.2020; (32911540): dkaa376
). These results indicate that the interactions with the bacterial OM of polymyxins, but not octapeptins, are more dependent on the electrostatic interactions between the positively charged side chains of Dab residues and the negatively charged phosphate groups of lipid A (
). This explains why the addition of positively charged moieties (e.g. ethanolamine or 4′-aminoarabinose) to the phosphate of lipid A significantly attenuates the interaction with polymyxins, but not octapeptins. In the lipid A–pEtN modified OM, the octapeptins strongly interacted with the lipid A–pEtN molecules through a combination of main chain–mediated polar interactions and hydrophobic contacts (Fig. 6), thus retaining the strong interaction with the outer membrane and subsequent membrane penetration. As noted above, this may explain why the octapeptins maintain activity against polymyxin-resistant bacteria (
). That the interaction of octapeptins with the bacterial OM is not highly dependent on electrostatic interactions with the lipid A phosphate groups may also explain their broad antimicrobial spectrum against Gram-positive bacteria, yeast, fungi, and protozoa that do not contain lipid A in their membranes (
Another interesting finding of our study was that the hydroxyl group on the fatty acyl group of octapeptin C4 formed multiple hydrogen bonds with the OM lipids (Fig. 6C), inhibiting the conformational transition of the molecule and increasing the free energy barrier for its penetration into the OM. A minor structural variation causing such a dramatic impact as observed in the present study (i.e. trapping the molecule in a semifolded state and creating a higher free energy barrier for OM penetration) indicates that the fatty acyl group of octapeptins is highly sensitive to structural modifications. This has important implications for antimicrobial activity. Indeed, previous studies with the octapeptins have shown that changes in the length of the fatty acyl group, as well as certain other structural modifications, significantly affect activity against Gram-negative bacteria (
). Although the present study is the first to examine the effect on membrane penetration of altering the hydroxyl group on the fatty acyl group of the octapeptins, the effects of other structural modifications at the remaining positions on membrane penetration and antimicrobial activity are warranted, with the aim of better understanding the structure–activity relationship of the octapeptins.
To the best of our knowledge, this is the first mechanistic study on the interaction between the octapeptins and the bacterial OM at the atomic level. The thermodynamics governing the penetration of the octapeptins into the OM and their conformational transitions within the OM environment were investigated. Importantly, our results help to explain why the octapeptins retain antimicrobial activity against polymyxin-resistant bacteria despite sharing a similar structural scaffold with the polymyxins. These novel findings enhance our understanding of the mechanism underlying the antimicrobial activity of the octapeptins and will facilitate the development of novel antibiotics for the treatment of infections caused by life-threatening Gram-negative ‘superbugs.‘
Experimental procedures
System preparation
A representative natural octapeptin (octapeptin C4) and a novel lead (FADDI-115) from our in-house library of octapeptin analogs were employed as the models (
). The only difference between these two octapeptins is that the hydroxyl group on the fatty acyl group of octapeptin C4 is absent in FADDI-115 (Fig. 1A). Importantly, both octapeptins show excellent activity against polymyxin-resistant Gram-negative bacteria (
). The structural models of the octapeptins were constructed in Chem3D. Energy minimization was performed based on the MM2 force field to relieve any intramolecular steric clashes (
), the octapeptin molecule formed a folded conformation in the simulations, which verified the topologies of the octapeptins. An asymmetric OM model was constructed using the CHARMM-GUI membrane builder (
). Notably, the accurate lipid composition of the OM (Fig. 1B) was based on our quantitative lipidomics results of polymyxin-resistant A. baumannii 5075R (
). Initially, the octapeptin molecules were placed above the outer leaflet of the OM using the gmx insert-molecules tool by replacing the overlapping water molecules (
). TIP3P water molecules and counter ions (CaCl2) were employed to hydrate and neutralize the simulation systems, respectively.
Steered MD
Because of the limitation of conventional MD simulations, it is not feasible to observe the spontaneous penetration of antimicrobials into the bacterial membrane within the time scale of hundreds to thousands of nanoseconds (
Polymyxin binding to the bacterial outer membrane reveals cation displacement and increasing membrane curvature in susceptible but not in resistant lipopolysaccharide chemotypes.
J. Chem. Inf. Model.2017; 57 (28805387): 2181-2193
) were conducted to construct the reaction coordinates of the penetration of the octapeptins into the lipid A–pEtN modified OM. A harmonic potential with a force constant of 500 KJ mol−1 nm−2 was applied on the centers of mass of both the octapeptin molecule and the OM along the z axis and changed adaptively according to the interaction between the octapeptin molecule and the OM. For each system, three independent steered simulations were performed.
Umbrella sampling
To characterize the equilibrated thermodynamics governing the penetration of the octapeptins into the OM, umbrella sampling simulations were performed (
). The system configurations from the steered simulation trajectory were used to construct a series of umbrella sampling windows in which the octapeptin molecule was harmonically restrained at various fixed positions at the z axis. The position interval between neighboring windows was 0.2 nm, and 17 simulation windows were constructed for each of the octapeptin C4 and FADDI-115 systems. Each window was simulated for 50 ns to enhance the conformational sampling; hence, a total of 850 ns of simulations were performed for each octapeptin system to obtain the free energy profile that quantitatively described the penetration of the octapeptin from bulk water to the hydrophobic center of the OM. The free energy profiles calculated using 10-ns time blocks and the overlap between position distributions of octapeptins in simulation windows were used to check the convergence of free energy calculations (Fig. S5). The WHAM integration method was used to calculate the free energy profile (
). The reaction coordinates were set based on the position of the octapeptin molecule relative to the OM center at z axis. Z > 0 indicated the bulk water shell and the outer leaflet of the OM; Z = 0 corresponded to the hydrophobic center of the OM.
MD simulations
GROMACS 5.1.2 was used to conduct all molecular simulations with the CHARMM36 all-atom force field (
). Energy minimizations were performed using the steepest descent method with the maximum force tolerance of 1,000 KJ mol−1 nm−1. To fully equilibrate the simulation system, six equilibration cycles were carried out by gradually turning off the position restraints on the lipids. Periodic boundary conditions were considered. All production simulations were conducted at constant temperature and pressure (NPT ensemble). The temperature and pressure were maintained at 313 K using the Nose–Hoover algorithm and 1 bar using the semi-isotropic pressure coupling method with Parrinello–Rhaman barostat (
). The electrostatic and Van der Waals interactions were calculated using the particle mesh Ewald method and Lennard–Jones potential algorithm, respectively (
Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study.
Polymyxin binding to the bacterial outer membrane reveals cation displacement and increasing membrane curvature in susceptible but not in resistant lipopolysaccharide chemotypes.
J. Chem. Inf. Model.2017; 57 (28805387): 2181-2193
Outer membranes of polymyxin-resistant Acinetobacter baumannii with phosphoethanolamine-modified lipid a and lipopolysaccharide loss display different atomic-scale interactions with polymyxins.
Molecular dynamics simulations informed by membrane lipidomics reveal structure–interaction relationship of polymyxins with lipid A–based outer membrane of Acinetobacter baumannii.
J. Antimicrob. Chemother.2020; (32911540): dkaa376
Author contributions—X. J. data curation; X. J. formal analysis; K. Y., B. Y., and F. S. methodology; B. G., L. Wan, N. A. P., K. D. R., and L. Wang resources; J. D. S., K. D. R., T. V., and J. L. writing-review and editing; T. V. and J. L. funding acquisition; J. L. conceptualization; J. L. supervision.
Funding and additional information—This work was supported by NIAID, National Institutes of Health Grant R01 AI132154 (to J. L. and T. V.). X. J. is recipient of the 2019 Faculty Bridging Fellowship from Monash University. J. L. is an Australia National Health Medical Research Council Principal Research Fellow. The simulations were performed on the supercomputer at Monash University (Australia) and the HPC Cloud Platform (National Key Research and Development Project 2016YFB0201702) at Shandong University (China). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Conflict of interest—The authors declare that they have no conflicts of interest with the contents of this article.
50419-6&id=gr1.jpg){kind=link}
50419-6&id=gr2.jpg){kind=link}
50419-6&id=gr3.jpg){kind=link}
50419-6&id=gr4.jpg){kind=link}
50419-6&id=gr5.jpg){kind=link}
50419-6&id=gr6.jpg){kind=link}