Dual-resolution molecular dynamics simulation of antimicrobials in biomembranes

Free Energy Barriers for Passive Drug Transport through the Mycobacterium tuberculosis Outer Membrane: A Molecular Dynamics Study

The emergence of multi-drug-resistant tuberculosis strains poses a significant challenge to modern medicine. The development of new antituberculosis drugs is hindered by the low permeability of many active compounds through the extremely strong bacterial cell wall of mycobacteria. In order to estimate the ability of potential antimycobacterial agents to diffuse through the outer mycolate membrane, the free energy profiles, the corresponding activation barriers, and possible permeability modes of passive transport for a series of known antibiotics, modern antituberculosis drugs, and prospective active drug-like molecules were determined using molecular dynamics simulations with the all-atom force field and potential of mean-force calculations. The membranes of different chemical and conformational compositions, density, thickness, and ionization states were examined. The typical activation barriers for the low-mass molecules penetrating through the most realistic membrane model were 6-13 kcal/mol for isoniazid, pyrazinamide, and etambutol, and 19 and 25 kcal/mol for bedaquilin and rifampicin. The barriers for the ionized molecules are usually in the range of 37-63 kcal/mol. The linear regression models were derived from the obtained data, allowing one to estimate the permeability barriers from simple physicochemical parameters of the diffusing molecules, notably lipophilicity and molecular polarizability.

The potential of antifungal peptide Sesquin as natural food preservative

Sesquin is a wide spectrum antimicrobial peptide displaying a remarkable activity on fungi. Contrarily to most antimicrobial peptides, it presents an overall negative charge. In the present study, we elucidate the molecular basis of its mode of action towards biomimetic membranes by NMR and MD experiments. While a specific recognition of phosphatidylethanolamine (PE) might explain its activity in a variety of different organisms (including bacteria), a further interaction with ergosterol accounts for its strong antifungal activity. NMR data reveal a charge gradient along its amide protons allowing the peptide to reach the membrane phosphate groups despite its negative charge. Subsequently, the peptide gets structured inside the bilayer, reducing its order. MD simulations predict that its activity is retained in conditions commonly used for food preservation: low temperatures, high pressure, or the presence of electric field pulses, making Sesquin a good candidate as food preservative.

Molecular Dynamics Simulations of Curved Lipid Membranes

Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and thereby be trafficked to specific locations. The description of curvature sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow researcher to closely match the conditions of experimental studies of membrane curvatures.

Antimicrobial Bombinin-like Peptide 3 Selectively Recognizes and Inserts into Bacterial Biomimetic Bilayers in Multiple Steps

Bombinins are a wide family of antimicrobial peptides from Xenopus skin. By sequence clustering, we highlighted at least three families named A, B, and H, which might exert antibacterial activity by different modes of action. In this work, we study bombinin-like peptide 3 (BLP-3) as a nonhemolytic representative of the quite unexplored class A due to its appealing activity toward WHO-priority-list bacteria such as Neisseria, Pseudomonas aeruginosa, and Staphylococcus aureus. A marked preference for cardiolipin and phosphatidylglycerol head groups, typically found in bacteria, is proven with biomimetic membranes studied by liquid and solid NMR and MD simulations. BLP-3 gets structured upon interaction and penetrates deeply into the bilayer in two steps involving a superficial insertion of key side chains and subsequent internalization. All along the pathway, a fundamental role is played by lysine residues in the conserved region 11-19, which act in synergy with other key residues.

Molecular Basis of the Anticancer and Antibacterial Properties of CecropinXJ Peptide: An In Silico Study

Esophageal cancer is an aggressive lethal malignancy causing thousands of deaths every year. While current treatments have poor outcomes, cecropinXJ (CXJ) is one of the very few peptides with demonstrated in vivo activity. The great interest in CXJ stems from its low toxicity and additional activity against most ESKAPE bacteria and fungi. Here, we present the first study of its mechanism of action based on molecular dynamics (MD) simulations and sequence-property alignment. Although unstructured in solution, predictions highlight the presence of two helices separated by a flexible hinge containing P24 and stabilized by the interaction of W2 with target biomembranes: an amphipathic helix-I and a poorly structured helix-II. Both MD and sequence-property alignment point to the important role of helix I in both the activity and the interaction with biomembranes. MD reveals that CXJ interacts mainly with phosphatidylserine (PS) but also with phosphatidylethanolamine (PE) headgroups, both found in the outer leaflet of cancer cells, while salt bridges with phosphate moieties are prevalent in bacterial biomimetic membranes composed of PE, phosphatidylglycerol (PG) and cardiolipin (CL). The antibacterial activity of CXJ might also explain its interaction with mitochondria, whose phospholipid composition recalls that of bacteria and its capability to induce apoptosis in cancer cells.

Antimicrobial Peptide K11 Selectively Recognizes Bacterial Biomimetic Membranes and Acts by Twisting Their Bilayers

K11 is a synthetic peptide originating from the introduction of a lysine residue in position 11 within the sequence of a rationally designed antibacterial scaffold. Despite its remarkable antibacterial properties towards many ESKAPE bacteria and its optimal therapeutic index (320), a detailed description of its mechanism of action is missing. As most antimicrobial peptides act by destabilizing the membranes of the target organisms, we investigated the interaction of K11 with biomimetic membranes of various phospholipid compositions by liquid and solid-state NMR. Our data show that K11 can selectively destabilize bacterial biomimetic membranes and torque the surface of their bilayers. The same is observed for membranes containing other negatively charged phospholipids which might suggest additional biological activities. Molecular dynamic simulations reveal that K11 can penetrate the membrane in four steps: after binding to phosphate groups by means of the lysine residue at the N-terminus (anchoring), three couples of lysine residues act subsequently to exert a torque in the membrane (twisting) which allows the insertion of aromatic side chains at both termini (insertion) eventually leading to the flip of the amphipathic helix inside the bilayer core (helix flip and internalization).

Comparative Target Analysis of Chlorinated Biphenyl Antimicrobials Highlights MenG as a Molecular Target of Triclocarban

Triclocarban (TCC), a formerly used disinfectant, kills bacteria via an unknown mechanism of action. A structural hallmark is its N,N'-diaryl urea motif, which is also present in other antibiotics, including the recently reported small molecule PK150. We show here that, like PK150, TCC exhibits an inhibitory effect on Staphylococcus aureus menaquinone metabolism via inhibition of the biosynthesis protein demethylmenaquinone methyltransferase (MenG). However, the activity spectrum (MIC₉₀) of TCC across a broad range of multidrug-resistant staphylococcus and enterococcus strains was much narrower than that of PK150. Accordingly, TCC did not cause an overactivation of signal peptidase SpsB, a hallmark of the PK150 mode of action. Furthermore, we were able to rule out inhibition of FabI, a confirmed target of the diaryl ether antibiotic triclosan (TCS). Differences in the target profiles of TCC and TCS were further investigated by proteomic analysis, showing complex but rather distinct changes in the protein expression profile of S. aureus Downregulation of the arginine deiminase pathway provided additional evidence for an effect on bacterial energy metabolism by TCC.IMPORTANCE TCC's widespread use as an antimicrobial agent has made it a ubiquitous environmental pollutant despite its withdrawal due to ecological and toxicological concerns. With its antibacterial mechanism of action still being unknown, we undertook a comparative target analysis between TCC, PK150 (a recently discovered antibacterial compound with structural resemblance to TCC), and TCS (another widely employed chlorinated biphenyl antimicrobial) in the bacterium Staphylococcus aureus We show that there are distinct differences in each compound's mode of action, but also identify a shared target between TCC and PK150, the interference with menaquinone metabolism by inhibition of MenG. The prevailing differences, however, which also manifest in a remarkably better broad-spectrum activity of PK150, suggest that even high levels of TCC or TCS resistance observed by continuous environmental exposure may not affect the potential of PK150 or related N,N'-diaryl urea compounds as new antibiotic drug candidates against multidrug-resistant infections.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

Interaction of Antimicrobial Lipopeptides with Bacterial Lipid Bilayers

The resistance of pathogens to traditional antibiotics is currently a global issue of enormous concern. As the discovery and development of new antibiotics become increasingly challenging, synthetic antimicrobial lipopeptides (AMLPs) are now receiving renewed attention as a new class of antimicrobial agents. In contrast to traditional antibiotics, AMLPs act by physically disrupting the cell membrane (rather than targeting specific proteins), thus reducing the risk of inducing bacterial resistance. In this study, we use microsecond-timescale atomistic molecular dynamics simulations to quantify the interaction of a short AMLP (C16-KKK) with model bacterial lipid bilayers. In particular, we investigate how fundamental transmembrane properties change in relation to a range of lipopeptide concentrations. A number of structural, mechanical, and dynamical features are found to be significantly altered in a non-linear fashion. At 10 mol% concentration, lipopeptides have a condensing effect on bacterial bilayers, characterized by a decrease in the area per lipid and an increase in the bilayer order. Higher AMLP concentrations of 25 and 40 mol% destabilize the membrane by disrupting the bilayer core structure, inducing membrane thinning and water leakage. Important transmembrane properties such as the lateral pressure and dipole potential profiles are also affected. Potential implications on membrane function and associated proteins are discussed.

Effect of Triclosan and Chloroxylenol on Bacterial Membranes

Triclosan and chloroxylenol are broad-spectrum biocides used extensively in healthcare and consumer products. They have been suggested to perturb the structure of bacterial membranes, but studies so far have not considered that most bacterial membranes contain large amounts of branched-chain lipids. Here, molecular dynamics simulation is used to examine the effect of the two biocides on membranes consisting of lipids with methyl-branched chains, cyclopropanated chains, and nonbranched chains. It is shown that triclosan and chloroxylenol induced a phase transition in membranes from a liquid-crystalline to a liquid-ordered phase irrespective of the presence and nature of branching groups. At high concentration, chloroxylenol promoted chain interdigitation. Our results suggest that triclosan and chloroxylenol decrease the degree of fluidity of membranes and that this effect is more pronounced in bacterial membranes. As a result, their biocidal activity could be associated with a change in the function of membrane proteins.

From quantum to subcellular scales: multi-scale simulation approaches and the SIRAH force field

Modern molecular and cellular biology profits from astonishing resolution structural methods, currently even reaching the whole cell level. This is encompassed by the development of computational methods providing a deep view into the structure and dynamics of molecular processes happening at very different scales in time and space. Linking such scales is of paramount importance when aiming at far-reaching biological questions. Computational methods at the interface between classical and coarse-grained resolutions are gaining momentum with several research groups dedicating important efforts to their development and tuning. An overview of such methods is addressed herein, with special emphasis on the SIRAH force field for coarse-grained and multi-scale simulations. Moreover, we provide proof of concept calculations on the implementation of a multi-scale simulation scheme including quantum calculations on a classical fine-grained/coarse-grained representation of double-stranded DNA. This opens the possibility to include the effect of large conformational fluctuations in chromatin segments on, for instance, the reactivity of particular base pairs within the same simulation framework.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Effects of lipid composition on membrane permeation

Passive permeation through lipid membranes is an essential process in biology. In vivo membranes typically consist of mixtures of lamellar and nonlamellar lipids. Lamellar lipids are characterized by their tendency to form lamellar sheet-like structures, which are predominant in nature. Nonlamellar lipids, when isolated, instead form more geometrically complex nonlamellar phases. While mixed lamellar/nonlamellar lipid membranes tend to adopt the ubiquitous lamellar bilayer structure, the presence of nonlamellar lipids is known to have profound effects on key membrane properties, such as internal distributions of stress and elastic properties, which in turn may alter related biological processes. This work focuses on one such process, i.e., permeation, by utilising atomistic molecular dynamics simulations in order to obtain transfer free energy profiles, diffusion profiles and permeation coefficients for a series of thirteen small molecules and drugs. Each permeant is tested on two bilayer membranes of different lipid composition, i.e., purely lamellar and mixed lamellar/nonlamellar. Our results indicate that the presence of nonlamellar lipids reduces permeation for smaller molecules (molecular weight < 100) but facilitates it for the largest ones (molecular weight > 100). This work represents an advancement towards the development of more realistic in silico permeability assays, which may have a substantial future impact in the area of rational drug design.

Triclosan Is an Aminoglycoside Adjuvant for Eradication of Pseudomonas aeruginosa Biofilms

One of the most important clinical obstacles in cystic fibrosis (CF) treatment is antibiotic treatment failure due to biofilms produced by Pseudomonas aeruginosa The ability of this pathogen to survive eradication by tobramycin and pathoadapt into a hyperbiofilm state leading to chronic infections is key to its success. Retrospective studies have demonstrated that preventing this pathoadaptation by improving eradication is essential to extend the lives of CF patients. To identify adjuvants that enhance tobramycin eradication of P. aeruginosa, we performed a high-throughput screen of 6,080 compounds from four drug-repurposing libraries. We identified that the Food and Drug Administration (FDA)-approved compound triclosan, in combination with tobramycin, resulted in a 100-fold reduction of viable cells within biofilms at 6 h, but neither compound alone had significant antimicrobial activity against biofilms. This synergistic treatment significantly accelerated the killing of biofilms compared to that with tobramycin treatment alone, and the combination was effective against 6/7 CF clinical isolates compared to tobramycin treatment alone, including a tobramycin-resistant strain. Further, triclosan and tobramycin killed persister cells, causing a 100-fold reduction by 8 h and complete eradication by 24 h. Triclosan also enhances tobramycin killing of multiple Burkholderia cenocepacia and Staphylococcus aureus clinical isolates grown as biofilms. Additionally, triclosan showed synergy with other aminoglycosides, such as gentamicin or streptomycin. Triclosan is a well-tolerated aminoglycoside adjuvant shown to be safe for human use that could improve the treatment of biofilm-based infections.

Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics

The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells' membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane's mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.

Multiscale Methods in Drug Design Bridge Chemical and Biological Complexity in the Search for Cures

Drug action is inherently multiscale: it connects molecular interactions to emergent properties at cellular and larger scales. Simulation techniques at each of these different scales are already central to drug design and development, but methods capable of connecting across these scales will extend understanding of complex mechanisms and the ability to predict biological effects. Improved algorithms, ever-more-powerful computing architectures and the accelerating growth of rich datasets are driving advances in multiscale modeling methods capable of bridging chemical and biological complexity from the atom to the cell. Particularly exciting is the development of highly detailed, structure-based, physical simulations of biochemical systems, which are now able to access experimentally relevant timescales for large systems and, at the same time, achieve unprecedented accuracy. In this Perspective, we discuss how emerging data-rich, physics-based multiscale approaches are of the cusp of realizing long-promised impact in the discovery, design and development of novel therapeutics. We highlight emerging methods and applications in this growing field, and outline how different scales can be combined in practical modelling and simulation strategies.

Effects of High Pressure on Phospholipid Bilayers

The response of lipid membranes to changes in external pressure is important for many biological processes, and it can also be exploited for technological applications. In this work, we employ all-atom molecular dynamics simulations to characterize the changes in the physical properties of phospholipid bilayers brought about by high pressure (1000 bar). In particular, we study how the response differs, in relation to different chain unsaturation levels, by comparing monounsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and biunsaturated dioleoyl-phosphatidylcholine (DOPC) bilayers. Various structural, mechanical, and dynamical features are found to be altered by the pressure increase in both bilayers. Notably, for most properties, including bilayer area and thickness, lipid order parameters, lateral pressure profile, and curvature frustration energy, we observe significantly more pronounced effects for monounsaturated POPC than biunsaturated DOPC. Possible biological implications of the results obtained are discussed, especially in relation to how different lipids can control the structure and function of membrane proteins.

Hybrid All-Atom/Coarse-Grained Simulations of Proteins by Direct Coupling of CHARMM and PRIMO Force Fields

Hybrid all-atom/coarse-grained (AA/CG) simulations of proteins offer a computationally efficient compromise where atomistic details are only applied to biologically relevant regions while benefiting from the speedup of treating the remaining parts of a given system at the CG level. The recently developed CG model, PRIMO, allows a direct coupling with an atomistic force field with no additional modifications or coupling terms and the ability to carry out dynamic simulations without any restraints on secondary or tertiary structures. A hybrid AA/CG scheme based on combining all-atom CHARMM and coarse-grained PRIMO representations was validated via molecular dynamics and replica exchange simulations of soluble and membrane proteins. The AA/CG scheme was also tested in the calculation of the free energy profile for the transition from the closed to the open state of adenylate kinase via umbrella sampling molecular dynamics method. The overall finding is that the AA/CG scheme generates dynamics and energetics that are qualitatively and quantitatively comparable to AA simulations while offering the computational advantages of coarse-graining. This model opens the door to challenging applications where high accuracy is required only in parts of large biomolecular complexes.

Functional cardiotoxicity assessment of cosmetic compounds using human-induced pluripotent stem cell-derived cardiomyocytes

There is a large demand of a human relevant in vitro test system suitable for assessing the cardiotoxic potential of cosmetic ingredients and other chemicals. Using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), we have already established an in vitro cardiotoxicity assay and identified genomic biomarkers of anthracycline-induced cardiotoxicity in our previous work. Here, five cosmetic ingredients were studied by the new hiPSC-CMs test; kojic acid (KJA), triclosan (TS), triclocarban (TCC), 2,7-naphthalenediol (NPT), and basic red 51 (BR51) based on cytotoxicity as well as ATP assays, beating rate, and genomic biomarkers to determine the lowest observed effect concentration (LOEC) and no observed effect concentration (NOEC). The LOEC for beating rate were 400, 10, 3, >400, and 3 µM for KJA, TS, TCC, NPT, and BR51, respectively. The corresponding concentrations for cytotoxicity or ATP depletion were similar, with the exception of TS and TCC, where the cardiomyocyte-beating assay showed positive results at non-cytotoxic concentrations. Functional analysis also showed that the individual compounds caused different effects on hiPSC-CMs. While exposure to KJA, TS, TCC, and BR51 induced significant arrhythmic beating, NPT slightly decreased cell viability, but did not influence beating. Gene expression studies showed that TS and NPT caused down-regulation of cytoskeletal and cardiac ion homeostasis genes. Moreover, TS and NPT deregulated genomic biomarkers known to be affected also by anthracyclines. The present study demonstrates that hiPSC-CMs can be used to determine LOECs and NOECs in vitro, which can be compared to human blood concentrations to determine margins of exposure. Our in vitro assay, which so far has been tested with several anthracyclines and cosmetics, still requires validation by larger numbers of positive and negative controls, before it can be recommended for routine analysis.

Thermodynamic and kinetic stability of the Josephin Domain closed arrangement: evidences from replica exchange molecular dynamics

BACKGROUND

Molecular phenomena driving pathological aggregation in neurodegenerative diseases are not completely understood yet. Peculiar is the case of Spinocerebellar Ataxia 3 (SCA3) where the conformational properties of the AT-3 N-terminal region, also called Josephin Domain (JD), play a key role in the first step of aggregation, having the JD an amyloidogenic propensity itself. For this reason, unraveling the intimate relationship between JD structural features and aggregation tendency may lead to a step forward in understanding the pathology and rationally design a cure. In this connection, computational modeling has demonstrated to be helpful in exploring the protein molecular dynamics and mechanism of action.

RESULTS

Conformational dynamics of the JD is here finely investigated by replica exchange molecular dynamics simulations able to sample the microsecond time scale and to provide both a thermodynamic and kinetic description of the protein conformational changes. Accessible structural conformations of the JD have been identified in: open, intermediate and closed like arrangement. Data indicated the closed JD arrangement as the most likely protein arrangement. The protein transition from closed toward intermediate/open states was characterized by a rate constant higher than 700 ns. This result also explains the inability of classical molecular dynamics to explore transitions from closed to open JD configuration on a time scale of hundreds of nanoseconds.

CONCLUSION

This work provides the first kinetic estimation of the JD transition pathway from open-like to closed-like arrangement and vice-versa, indicating the closed-like arrangement as the most likely configuration for a JD in water environment. More widely, the importance of our results is also underscored considering that the ability to provide a kinetic description of the protein conformational changes is a scientific challenge for both experimental and theoretical approaches to date.

REVIEWERS

This article was reviewed by Oliviero Carugo, Bojan Zagrovic.

Free Energy Calculations for the Peripheral Binding of Proteins/Peptides to an Anionic Membrane. 1. Implicit Membrane Models

The binding of peptides and proteins to the surface of complex lipid membranes is important in many biological processes such as cell signaling and membrane remodeling. Computational studies can aid experiments by identifying physical interactions and structural motifs that determine the binding affinity and specificity. However, previous studies focused on either qualitative behaviors of protein/membrane interactions or the binding affinity of small peptides. Motivated by this observation, we set out to develop computational protocols for bimolecular binding to charged membranes that are applicable to both peptides and large proteins. In this work, we explore a method based on an implicit membrane/solvent model (generalized Born with a simple switching in combination with the Gouy-Chapman-Stern model for a charged interface), which we expect to lead to useful results when the binding does not implicate significant membrane deformation and local demixing of lipids. We show that the binding free energy can be efficiently computed following a thermodynamic cycle similar to protein-ligand binding calculations, especially when a Bennett acceptance ratio based protocol is used to consider both the membrane bound and solution conformational ensembles. Test calculations on a series of peptides show that our computational approach leads to binding affinities in encouraging agreement with experimental data, including for the challenging example of the bringing of flexible MARCKS-ED peptides to membranes. The calculations highlight that for a membrane with a significant fraction of anionic lipids, it is essential to include the effect of ion adsorption using the Stern model, which significantly modifies the effective surface charge. This implicit membrane model based computational protocol helps lay the groundwork for more systematic analysis of protein/peptide binding to membranes of complex shape and composition.

Effects of Lipid Composition on Bilayer Membranes Quantified by All-Atom Molecular Dynamics

Biological bilayer membranes typically contain varying amounts of lamellar and nonlamellar lipids. Lamellar lipids, such as dioleoylphosphatidylcholine (DOPC), are defined by their tendency to form the lamellar phase, ubiquitous in biology. Nonlamellar lipids, such as dioleoylphosphatidylethanolamine (DOPE), prefer instead to form nonlamellar phases, which are mostly nonbiological. However, nonlamellar lipids mix with lamellar lipids in biomembrane structures that remain overall lamellar. Importantly, changes in the lamellar vs nonlamellar lipid composition are believed to affect membrane function and modulate membrane proteins. In this work, we employ atomistic molecular dynamics simulations to quantify how a range of bilayer properties are altered by variations in the lamellar vs nonlamellar lipid composition. Specifically, we simulate five DOPC/DOPE bilayers at mixing ratios of 1/0, 3/1, 1/1, 1/3, and 0/1. We examine properties including lipid area and bilayer thickness, as well as the transmembrane profiles of electron density, lateral pressure, electric field, and dipole potential. While the bilayer structure is only marginally altered by lipid composition changes, dramatic effects are observed for the lateral pressure, electric field, and dipole potential profiles. Possible implications for membrane function are discussed.

A simple and transferable all-atom/coarse-grained hybrid model to study membrane processes

We present an efficient all-atom/coarse-grained hybrid model and apply it to membrane processes. This model is an extension of the all-atom/ELBA model applied previously to processes in water. Here, we improve the efficiency of the model by implementing a multiple-time step integrator that allows the atoms and the coarse-grained beads to be propagated at different timesteps. Furthermore, we fine-tune the interaction between the atoms and the coarse-grained beads by computing the potential of mean force of amino acid side chain analogs along the membrane normal and comparing to atomistic simulations. The model was independently validated on the calculation of small-molecule partition coefficients. Finally, we apply the model to membrane peptides. We studied the tilt angle of the Walp23 and Kalp23 helices in two different model membranes and the stability of the glycophorin A dimer. The model is efficient, accurate, and straightforward to use, as it does not require any extra interaction particles, layers of atomistic solvent molecules or tabulated potentials, thus offering a novel, simple approach to study membrane processes.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Triclosan Computational Conformational Chemistry Analysis for Antimicrobial Properties in Polymers

Triclosan is a diphenyl ether antimicrobial that has been analyzed by computational conformational chemistry for an understanding of Mechanomolecular Theory. Subsequent energy profile analysis combined with easily seen three-dimensional chemistry structure models for the nonpolar molecule Triclosan show how single bond rotations can alternate rapidly at a polar and nonpolar interface. Bond rotations for the center ether oxygen atom of the two aromatic rings then expose or hide nonbonding lone-pair electrons for the oxygen atom depending on the polar nature of the immediate local molecular environment. Rapid bond movements can subsequently produce fluctuations as vibration energy. Consequently, related mechanical molecular movements calculated as energy relationships by forces acting through different bond positions can help improve on current Mechanomolecular Theory. A previous controversy reported as a discrepancy in literature contends for a possible bacterial resistance from Triclosan antimicrobial. However, findings in clinical settings have not reported a single case for Triclosan bacterial resistance in over 40 years that has been documented carefully in government reports. As a result, Triclosan is recommended whenever there is a health benefit consistent with a number of approvals for use of Triclosan in healthcare devices. Since Triclosan is the most researched antimicrobial ever, literature meta analysis with computational chemistry can best describe new molecular conditions that were previously impossible by conventional chemistry methods. Triclosan vibrational energy can now explain the molecular disruption of bacterial membranes. Further, Triclosan mechanomolecular movements help illustrate use in polymer matrix composites as an antimicrobial with two new additive properties as a toughening agent to improve matrix fracture toughness from microcracking and a hydrophobic wetting agent to help incorporate strengthening fibers. Interrelated Mechanomolecular Theory by oxygen atom bond rotations or a nitrogen-type pyramidal inversion can be shown to produce energy at a polar and nonpolar boundary condition to better make clear membrane transport of other molecules, cell recognition/signaling/defense and enzyme molecular "mixing" action.

Signaling factor interactions with polysaccharide aggregates of bacterial biofilms

Biofilms are surface-attached colonies of bacteria embedded in an extracellular polymeric substance (EPS). Inside the eukaryotic hosts, bacterial biofilms interact with the host cells through signaling factors (SFs). These signaling processes play important roles in the interaction between bacteria and host cells and the outcome of infections and symbiosis. However, how host immune factors diffuse through biofilms is not well understood. Here, we describe synergistic molecular dynamics and experimental approaches for studying the translocation of signaling factors through polysaccharide chain aggregates present in the extracellular matrix of bacterial biofilms. The effect of polysaccharide chain degradation on the energetics of SF-EPS interactions was examined by simulating an EPS consisting of various polysaccharide chain lengths. It is shown that the SF stabilization energy, defined as the average potential of mean force difference between the environments outside and within the matrix, increases linearly with decreasing chain length. This effect has been explained based on the changes in the polysaccharide configurations around the SF. Specifically, shorter chains are packed tightly around the SF, promoting favorable SF-EPS interactions, while longer chains are packed loosely resulting in screening of interactions with neighboring chains. We further investigated the translocation of SFs through the host cell membrane using molecular dynamics simulations. Further, simulations predict the existence of energy barriers greater than 1000 kJ mol(-1) associated with the translocation of the signaling factors necrosis factor-alpha (TNF-α) and granulocyte macrophage colony stimulating factor (GM-CSF) across the lipid bilayer. The agreement of computational and experimental findings motivates future computational studies using a more detailed description of the EPS aimed at understanding the role of the extracellular matrix on biofilm drug resistance.

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) due to difficulties associated with expression and crystallization. As such, it is one area where computational studies, particularly Molecular Dynamics (MD) simulations, can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, though not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

Direct Mixing of Atomistic Solutes and Coarse-Grained Water

We present a new dual-resolution approach for coupling atomistic and coarse-grained models in molecular dynamics simulations of hydrated systems. In particular, a coarse-grained point dipolar water model is used to solvate molecules represented with standard all-atom force fields. A unique characteristic of our methodology is that the mixing of resolutions is direct, meaning that no additional or ad hoc scaling factors, intermediate regions, or extra sites are required. To validate the methodology, we compute the hydration free energy of 14 atomistic small molecules (analogs of amino acid side chains) solvated by the coarse-grained water. Remarkably, our predictions reproduce the experimental data as accurately as the predictions from state-of-the-art fully atomistic simulations. We also show that the hydration free energy of the coarse-grained water itself is in comparable or better agreement with the experimental value than the predictions from all but one of the most common multisite atomistic models. The coarse-grained water is then applied to solvate a typical atomistic protein containing both α-helix and β-strand elements. Moreover, parallel tempering simulations are performed to investigate the folding free energy landscape of a representative α helical and a β hairpin structure. For the simulations considered in this work, our dual-resolution method is found to be 3 to 6 times more computationally efficient than corresponding fully atomistic approaches.

A subset of membrane-altering agents and γ-secretase modulators provoke nonsubstrate cleavage by rhomboid proteases

Rhomboid proteases are integral membrane enzymes that regulate cell signaling, adhesion, and organelle homeostasis pathways, making substrate specificity a key feature of their function. Interestingly, we found that perturbing the membrane pharmacologically in living cells had little effect on substrate processing but induced inappropriate cleavage of nonsubstrates by rhomboid proteases. A subclass of drugs known to modulate γ-secretase activity acted on the membrane directly and induced nonsubstrate cleavage by rhomboid proteases but left true substrate cleavage sites unaltered. These observations highlight an active role for the membrane in guiding rhomboid selectivity and caution that membrane-targeted drugs should be evaluated for cross-activity against membrane-resident enzymes that are otherwise unrelated to the intended drug target. Furthermore, some γ-secretase-modulating activity or toxicity could partly result from global membrane effects.

Computational conformational antimicrobial analysis developing mechanomolecular theory for polymer biomaterials in materials science and engineering

Single-bond rotations or pyramidal inversions tend to either hide or expose relative energies that exist for atoms with nonbonding lone-pair electrons. Availability of lone-pair electrons depends on overall molecular electron distributions and differences in the immediate polarity of the surrounding pico/nanoenvironment. Stereochemistry three-dimensional aspects of molecules provide insight into conformations through single-bond rotations with associated lone-pair electrons on oxygen atoms in addition to pyramidal inversions with nitrogen atoms. When electrons are protected, potential energy is sheltered toward an energy minimum value to compatibilize molecularly with nonpolar environments. When electrons are exposed, maximum energy is available toward polar environment interactions. Computational conformational analysis software calculated energy profiles that exist during specific oxygen ether single-bond rotations with easy-to-visualize three-dimensional models for the trichlorinated bisaromatic ether triclosan antimicrobial polymer additive. As shown, fluctuating alternating bond rotations can produce complex interactions between molecules to provide entanglement strength for polymer toughness or alternatively disrupt weak secondary bonds of attraction to lower resin viscosity for new additive properties with nonpolar triclosan as a hydrophobic toughening/wetting agent. Further, bond rotations involving lone-pair electrons by a molecule at a nonpolar-hydrocarbon-membrane/polar-biologic-fluid interface might become sufficiently unstable to provide free mechanomolecular energies to disrupt weaker microbial membranes, for membrane transport of molecules into cells, provide cell signaling/recognition/defense and also generate enzyme mixing to speed reactions.

Theoretical frameworks for multiscale modeling and simulation

Biomolecular systems have been modeled at a variety of scales, ranging from explicit treatment of electrons and nuclei to continuum description of bulk deformation or velocity. Many challenges of interfacing between scales have been overcome. Multiple models at different scales have been used to study the same system or calculate the same property (e.g., channel conductance). Accurate modeling of biochemical processes under in vivo conditions and the bridging of molecular and subcellular scales will likely soon become reality.

The simulation approach to lipid-protein interactions

The interactions between lipids and proteins are crucial for a range of biological processes, from the folding and stability of membrane proteins to signaling and metabolism facilitated by lipid-binding proteins. However, high-resolution structural details concerning functional lipid/protein interactions are scarce due to barriers in both experimental isolation of native lipid-bound complexes and subsequent biophysical characterization. The molecular dynamics (MD) simulation approach provides a means to complement available structural data, yielding dynamic, structural, and thermodynamic data for a protein embedded within a physiologically realistic, modelled lipid environment. In this chapter, we provide a guide to current methods for setting up and running simulations of membrane proteins and soluble, lipid-binding proteins, using standard atomistically detailed representations, as well as simplified, coarse-grained models. In addition, we outline recent studies that illustrate the power of the simulation approach in the context of biologically relevant lipid/protein interactions.

Multiscale molecular dynamics simulations of membrane proteins

The time and length scales accessible by biomolecular simulations continue to increase. This is in part due to improvements in algorithms and computing performance, but is also the result of the emergence of coarse-grained (CG) potentials, which complement and extend the information obtainable from fully detailed models. CG methods have already proven successful for a range of applications that benefit from the ability to rapidly simulate spontaneous self-assembly within a lipid membrane environment, including the insertion and/or oligomerization of a range of "toy models," transmembrane peptides, and single- and multi-domain proteins. While these simplified approaches sacrifice atomistic level detail, it is now straightforward to "reverse map" from CG to atomistic descriptions, providing a strategy to assemble membrane proteins within a lipid environment, prior to all-atom simulation. Moreover, recent developments have been made in "dual resolution" techniques, allowing different molecules in the system to be modeled with atomistic or CG resolution simultaneously.

The ELBA force field for coarse-grain modeling of lipid membranes

A new coarse-grain model for molecular dynamics simulation of lipid membranes is presented. Following a simple and conventional approach, lipid molecules are modeled by spherical sites, each representing a group of several atoms. In contrast to common coarse-grain methods, two original (interdependent) features are here adopted. First, the main electrostatics are modeled explicitly by charges and dipoles, which interact realistically through a relative dielectric constant of unity (ε(r) = 1). Second, water molecules are represented individually through a new parametrization of the simple Stockmayer potential for polar fluids; each water molecule is therefore described by a single spherical site embedded with a point dipole. The force field is shown to accurately reproduce the main physical properties of single-species phospholipid bilayers comprising dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) in the liquid crystal phase, as well as distearoylphosphatidylcholine (DSPC) in the liquid crystal and gel phases. Insights are presented into fundamental properties and phenomena that can be difficult or impossible to study with alternative computational or experimental methods. For example, we investigate the internal pressure distribution, dipole potential, lipid diffusion, and spontaneous self-assembly. Simulations lasting up to 1.5 microseconds were conducted for systems of different sizes (128, 512 and 1058 lipids); this also allowed us to identify size-dependent artifacts that are expected to affect membrane simulations in general. Future extensions and applications are discussed, particularly in relation to the methodology's inherent multiscale capabilities.