Structure of the Bacterial Cell Envelope and Interactions with Antimicrobials: Insights from Molecular Dynamics Simulations.

Bacteria have evolved over 3 billion years, shaping our intrinsic and symbiotic coexistence with these single-celled organisms. With rising populations of drug-resistant strains, the search for novel antimicrobials is an ongoing area of research. Advances in high-performance computing platforms have led to a variety of molecular dynamics simulation strategies to study the interactions of antimicrobial molecules with different compartments of the bacterial cell envelope of both Gram-positive and Gram-negative species. In this review, we begin with a detailed description of the structural aspects of the bacterial cell envelope. Simulations concerned with the transport and associated free energy of small molecules and ions through the outer membrane, peptidoglycan, inner membrane and outer membrane porins are discussed. Since surfactants are widely used as antimicrobials, a section is devoted to the interactions of surfactants with the cell wall and inner membranes. The review ends with a discussion on antimicrobial peptides and the insights gained from the molecular simulations on the free energy of translocation. Challenges involved in developing accurate molecular models and coarse-grained strategies that provide a trade-off between atomic details with a gain in sampling time are highlighted. The need for efficient sampling strategies to obtain accurate free energies of translocation is also discussed. Molecular dynamics simulations have evolved as a powerful tool that can potentially be used to design and develop novel antimicrobials and strategies to effectively treat bacterial infections.

Elucidating collective translocation of nanoparticles across the skin lipid matrix: a molecular dynamics study.

The top layer of skin, the stratum corneum, provides a formidable barrier to the skin. Nanoparticles are utilized and further explored for personal and health care applications related to the skin. In the past few years, several researchers have studied the translocation and permeation of nanoparticles of various shapes, sizes, and surface chemistry through cell membranes. Most of these studies focused on a single nanoparticle and a simple bilayer system, whereas skin has a highly complex lipid membrane architecture. Moreover, it is highly unlikely that a nanoparticle formulation applied on the skin will not have multiple nanoparticle-nanoparticle and skin-nanoparticle interactions. In this study, we have utilized coarse-grained MARTINI molecular dynamics simulations to assess the interactions of two types (bare and dodecane-thiol coated) of nanoparticles with two models (single bilayer and double bilayer) of skin lipid membranes. The nanoparticles were found to be partitioned from the water layer to the lipid membrane as an individual entity as well as in the cluster form. It was discovered that each nanoparticle reached the interior of both single bilayer and double bilayer membranes irrespective of the nanoparticle type and concentration, though coated particles were observed to efficiently traverse across the bilayer when compared with bare particles. The coated nanoparticles also created a single large cluster inside the membrane, whereas the bare nanoparticles were found in small clusters. Both the nanoparticles exhibited preferential interactions with cholesterol molecules present in the lipid membrane as compared to other lipid components of the membrane. We have also observed that the single membrane model exhibited unrealistic instability at moderate to higher concentrations of nanoparticles, and hence for translocation study, a minimum double bilayer model should be employed.

Interactions of Surfactants with the Bacterial Cell Wall and Inner Membrane: Revealing the Link between Aggregation and Antimicrobial Activity.

Surfactants with their intrinsic ability to solubilize lipid membranes are widely used as antibacterial agents, and their interactions with the bacterial cell envelope are complicated by their differential aggregation tendencies. We present a combined experimental and molecular dynamics investigation to unravel the molecular basis for the superior antimicrobial activity and faster kill kinetics of shorter-chain fatty acid surfactant, laurate, when compared with the longer-chain surfactants studied in contact time assays with live Escherichia coli (E. coli). From all-atom molecular dynamics simulations, translocation events across peptidoglycan were the highest for laurate followed by sodium dodecyl sulfate, myristate, palmitate, oleate, and stearate. The translocation kinetics were positively correlated with the critical micellar concentration, which determined the free monomer surfactant concentration available for translocation across peptidoglycan. Interestingly, aggregates showed a lower propensity to translocate across the peptidoglycan layer and longer translocation times were observed for oleate, thereby revealing an intrinsic sieving property of the bacterial cell wall. Molecular dynamics simulations with surfactant-incorporated bacterial inner membranes revealed the greatest hydrophobic mismatch and membrane thinning in the presence of laurate when compared with the other surfactants. The enhanced antimicrobial efficacy of laurate over oleate was further verified by experiments with giant unilamellar vesicles, and electroporation molecular dynamics simulations revealed greater inner membrane poration tendency in the presence of laurate when compared with the longer-chain surfactants. Our study provides molecular insights into surfactant translocation across peptidoglycan and chain length-induced structural disruption of the inner membrane, which correlate with contact time kill efficacies observed as a function of chain length with E. coli. The insights gained from our study uncover unexplored barrier properties of the bacterial cell envelope to rationalize the development of antimicrobial formulations and therapeutics.

Differentiating interactions of antimicrobials with Gram-negative and Gram-positive bacterial cell walls using molecular dynamics simulations.

Developing molecular models to capture the complex physicochemical architecture of the bacterial cell wall and to study the interaction with antibacterial molecules is an important aspect of assessing and developing novel antimicrobial molecules. We carried out molecular dynamics simulations using an atomistic model of peptidoglycan to represent the architecture for Gram-positive S. aureus. The model is developed to capture various structural features of the Staphylococcal cell wall, such as the peptide orientation, area per disaccharide, glycan length distribution, cross-linking, and pore size. A comparison of the cell wall density and electrostatic potentials is made with a previously developed cell wall model of Gram-negative bacteria, E. coli, and properties for both single and multilayered structures of the Staphylococcal cell wall are studied. We investigated the interactions of the antimicrobial peptide melittin with peptidoglycan structures. The depth of melittin binding to peptidoglycan is more pronounced in E. coli than in S. aureus, and consequently, melittin has greater contacts with glycan units of E. coli. Contacts of melittin with the amino acids of peptidoglycan are comparable across both the strains, and the D-Ala residues, which are sites for transpeptidation, show enhanced interactions with melittin. A low energetic barrier is observed for translocation of a naturally occurring antimicrobial thymol with the four-layered peptidoglycan model. The molecular model developed for Gram-positive peptidoglycan allows us to compare and contrast the cell wall penetrating properties with Gram-negative strains and assess for the first time binding and translocation of antimicrobial molecules for Gram-positive cell walls.

A Molecular Dynamics Study of Antimicrobial Peptide Interactions with the Lipopolysaccharides of the Outer Bacterial Membrane.

With rising bacterial resistance, antimicrobial peptides (AMPs) have been widely investigated as potential antibacterial molecules to replace conventional antibiotics. Our understanding of the molecular mechanisms for membrane disruption are largely based on AMP interactions with the inner phospholipid bilayers of both Gram-negative and Gram-positive bacteria. Mechanisms for AMP translocation across the outer membrane of Gram-negative bacteria composed of lipopolysaccharides and the asymmetric lipid bilayer are complicated by the secondary structure adopted by the peptide in the different membrane environments. We have employed atomistic molecular dynamics and umbrella-sampling simulations with an aggregate duration of [Formula: see text] 6 microseconds to obtain the free energy landscape of CM15 peptide translocating through the lipopolysaccharide region of Gram-negative bacteria, E. coli. The peptide has a favorable binding-free energy (- 130 kJ mol[Formula: see text]) in the O-antigen region with a large barrier (150 kJ mol[Formula: see text]) at the interface between the anionic core saccharides and upper bilayer leaflet made up of lipid-A molecules. Restraint-free molecular dynamics simulations show that the random coil structure is favored over the helix in both the extracellular aqueous region and the cation-rich core-saccharide regions of the outer membrane. The peptide and membrane properties are analyzed at each of the 100 ns duration of the umbrella-sampling windows to illustrate changes in peptide length, orientation, and hydration. Our study provides insights into the free energy landscape for the insertion of the AMP CM15 in the outer membrane of Gram-negative bacteria, and we discuss the implications of our findings with the broader question of how AMPs overcome this barrier during antimicrobial activity.

Evaluating Coarse-Grained MARTINI Force-Fields for Capturing the Ripple Phase of Lipid Membranes.

Phospholipids, which are an integral component of cell membranes, exhibit a rich variety of lamellar phases modulated by temperature and composition. Molecular dynamics (MD) simulations have greatly enhanced our understanding of phospholipid membranes by capturing experimentally observed phases and phase transitions at molecular resolution. However, the ripple (Pß') membrane phase, observed as an intermediate phase below the main gel-to-liquid crystalline transition with some lipids, has been challenging to capture with MD simulations, both at all-atom and coarse-grained (CG) resolutions. Here, with an aggregate ∼2.5 µs all-atom and ∼122 µs CGMD simulations, we systematically assess the ability of six CG MARTINI 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid and water force-field (FF) variants, parametrized to capture the DPPC gel and fluid phases, for their ability to capture the Pß' phase, and compared observations with those from an all-atom FF. Upon cooling from the fluid phase to below the phase transition temperature with smaller (380-lipid) and larger (>2200-lipid) MARTINI and all-atom (CHARMM36 FF) DPPC lipid bilayers, we observed that smaller bilayers with both all-atom and MARTINI FFs sampled interdigitated Pß' and ripple-like states, respectively. However, while all-atom simulations of the larger DPPC membranes exhibited the formation of the Pß' phase, MARTINI membranes did not sample interdigitated ripple-like states at larger system sizes. We then demonstrated that the ripple-like states in smaller MARTINI membranes were kinetically trapped structures caused by finite size effects rather than being representative of true Pß' phases. We showed that a MARTINI FF variant that could capture the tilted Lß' gel phase, a prerequisite for stabilizing the Pß' phase, was unable to capture the rippled phase upon cooling. Our study reveals that the current MARTINI FFs (including MARTINI3) may require specific reparametrization of the interaction potentials to stabilize lipid interdigitation, a characteristic of the ripple phase.

Assessing Barriers for Antimicrobial Penetration in Complex Asymmetric Bacterial Membranes: A Case Study with Thymol.

The bacterial cell envelope is a complex multilayered structure evolved to protect bacteria in hostile environments. An understanding of the molecular basis for the interaction and transport of antibacterial therapeutics with the bacterial cell envelope will enable the development of drug molecules to combat bacterial infections and suppress the emergence of drug-resistant strains. Here we report the successful creation of an in vitro supported lipid bilayer (SLB) platform of the outer membrane (OM) of E. coli, an archetypical Gram-negative bacterium, containing the full smooth lipopolysaccharide (S-LPS) architecture of the membrane. Using this platform, we performed fluorescence correlation spectroscopy (FCS) in combination with molecular dynamics (MD) simulations to measure lipid diffusivities and provide molecular insights into the transport of natural antimicrobial agent thymol. Lipid diffusivities measured on symmetric supported lipid bilayers made up of inner membrane lipids show a distinct increase in the presence of thymol as also corroborated by MD simulations. However, lipid diffusivities in the asymmetric OM consisting of only S-LPS are invariant upon exposure to thymol. Increasing the phospholipid content in the LPS-containing outer leaflet improved the penetration toward thymol as reflected in slightly higher relative diffusivity changes in the inner leaflet when compared with the outer leaflet. Free-energy computations reveal the presence of a barrier (∼6 kT) only in the core-saccharide region of the OM for the translocation of thymol while the external O-antigen part is easily traversed. In contrast, thymol spontaneously inserts into the inner membrane. In addition to providing leaflet-resolved penetration barriers in bacterial membranes, we also assess the ability of small molecules to penetrate various membrane components. With rising bacterial resistance, our study opens up the possibility of screening potential antimicrobial drug candidates using these realistic model platforms for Gram-negative bacteria.

Developing a Coarse-Grained Model for Bacterial Cell Walls: Evaluating Mechanical Properties and Free Energy Barriers.

The bacterial cell envelope of Gram-negative bacteria is a complex biological barrier with multiple layers consisting of the inner membrane, periplasm of peptidoglycan, and the outer membrane with lipopolysaccharides (LPS). With rising antimicrobial resistance there is increasing interest in understanding interactions of small molecules with the cell membrane to aid in the development of novel drug molecules. Hence suitable representations of the bacterial membrane are required to carry out meaningful molecular dynamics simulations. Given the complexity of the cell envelope, fully atomistic descriptions of the cell membrane with explicit solvent are computationally prohibitive, allowing limited sampling with small system sizes. However, coarse-grained (CG) models such as MARTINI allow one to study phenomena at physiologically relevant length and time scales. Although MARTINI models for lipids and the LPS are available in literature, a suitable CG model of peptidoglycan is lacking. Using an all-atom model described by Gumbart et al. [PLoS Comput. Biol. 2014, 10, e1003475], we develop a CG model of the peptidoglycan network within the MARTINI framework. The model is parametrized to reproduce the end-to-end distance of glycan strands. The structural properties such as the equilibrium angle between adjacent peptides along the strands, area per disaccharide, and cavity size distributions agree well with the atomistic simulation results. Mechanical properties such as the area compressibility and the bending modulus are accurately reproduced. While developing novel antibiotics it is important to assess barrier properties of the peptidogylcan network. We evaluate and compare the free energy of insertion for a thymol molecule using umbrella sampling on both the MARTINI and all-atom peptidoglycan models. The insertion free energy was found to be less than kBT for both the MARTINI and all-atom models. Additional restraint free simulations reveal rapid translocation of thymol across peptidogylcan. We expect that the proposed MARTINI model for peptidoglycan will be useful in understanding phenomena associated with bacterial cell walls at larger length and time scales, thereby overcoming the current limitations of all-atom models.