Dan S. Bolintineanu

Computational studies of protegrin antimicrobial peptides: A review

Antimicrobial peptides (AMPs) are small, naturally occurring peptides that exhibit strong antibacterial properties generally believed to be a result of selective bacterial membrane disruption. As a result, there has been significant interest in the development of therapeutic antibiotics based on AMPs; however, the poor understanding of the fundamental mechanism of action of these peptides has largely hampered such efforts. We present a summary of computational and theoretical investigations of protegrin, a particularly potent peptide that is both an excellent model for the mechanism of action of AMPs and a promising therapeutic candidate. Experimental investigations have shed light on many of the key steps in the action of protegrin: protegrin monomers are known to dimerize in various lipid environments; protegrin peptides interact strongly with lipid bilayer membranes, particularly anionic lipids; protegrins have been shown to form pores in lipid bilayers, which results in uncontrolled ion transport and may be a key factor in bacterial death. In this work, we present a comprehensive review of the computational and theoretical studies that have complemented and extended the information obtained from experimental work with protegrins, as well as a brief survey of the experimental biophysical studies that are most pertinent to such computational work. We show that a consistent, mechanistic description of the bactericidal mechanism of action of protegrins is emerging, and briefly outline areas where the current understanding is deficient. We hope that the research reviewed herein offers compelling evidence of the benefits of computational investigations of protegrins and other AMPs, as well as providing a useful guide to future work in this area.

Poisson-Nernst-Planck Models of Nonequilibrium Ion Electrodiffusion through a Protegrin Transmembrane Pore

Protegrin peptides are potent antimicrobial agents believed to act against a variety of pathogens by forming nonselective transmembrane pores in the bacterial cell membrane. We have employed 3D Poisson-Nernst-Planck (PNP) calculations to determine the steady-state ion conduction characteristics of such pores at applied voltages in the range of −100 to +100 mV in 0.1 M KCl bath solutions. We have tested a variety of pore structures extracted from molecular dynamics (MD) simulations based on an experimentally proposed octomeric pore structure. The computed single-channel conductance values were in the range of 290–680 pS. Better agreement with the experimental range of 40–360 pS was obtained using structures from the last 40 ns of the MD simulation, where conductance values range from 280 to 430 pS. We observed no significant variation of the conductance with applied voltage in any of the structures that we tested, suggesting that the voltage dependence observed experimentally is a result of voltage-dependent channel formation rather than an inherent feature of the open pore structure. We have found the pore to be highly selective for anions, with anionic to cationic current ratios (ICl−/IK+) on the order of 103. This is consistent with the highly cationic nature of the pore but surprisingly in disagreement with the experimental finding of only slight anionic selectivity. We have additionally tested the sensitivity of our PNP model to several parameters and found the ion diffusion coefficients to have a significant influence on conductance characteristics. The best agreement with experimental data was obtained using a diffusion coefficient for each ion set to 10% of the bulk literature value everywhere inside the channel, a scaling used by several other studies employing PNP calculations. Overall, this work presents a useful link between previous work focused on the structure of protegrin pores and experimental efforts aimed at investigating their conductance characteristics.

Effects of functional groups and ionization on the structure of alkanethiol coated gold nanoparticles

We report classical atomistic molecular dynamics simulations of alkanethiol-coated gold nanoparticles solvated in water and decane, as well as at water/vapor interfaces. The structure of the coatings is analyzed as a function of various functional end groups, including amine and carboxyl groups in various ionization states. We study both neutral and charged end groups for two different chain lengths (9 and 17 carbons). For the charged end groups, we simulated both mono- and divalent counterions. For the longer alkanes, we find significant local bundling of chains on the nanoparticle surface, which results in highly asymmetric coatings. In general, the charged end groups attenuate this effect by enhancing the water solubility of the nanoparticles. On the basis of the coating structures and density profiles, we can qualitatively infer the overall solubility of the nanoparticles. This asymmetry in the alkanethiol coatings is likely to have a significant effect on aggregation behavior. Our simulations elucidate the mechanism by which modulating the end group charge state can be used to control coating structure and therefore nanoparticle solubility and aggregation behavior.

Hydrogen-bonded aggregates in precise acid copolymers

We perform atomistic molecular dynamics simulations of melts of four precise acid copolymers, two poly(ethylene-co-acrylic acid) (PEAA) copolymers, and two poly(ethylene-co-sulfonic acid) (PESA) copolymers. The acid groups are spaced by either 9 or 21 carbons along the polymer backbones. Hydrogen bonding causes the acid groups to form aggregates. These aggregates give rise to a low wavevector peak in the structure factors, in agreement with X-ray scattering data for the PEAA materials. The structure factors for the PESA copolymers are very similar to those for the PEAA copolymers, indicating a similar distance between aggregates which depends on the spacer length but not on the nature of the acid group. The PEAA copolymers are found to form more dimers and other small aggregates than do the PESA copolymers, while the PESA copolymers have both more free acid groups and more large aggregates.

Heterogeneous Coordination Environments in Lithium-Neutralized Ionomers Identified Using 1H and 7Li MAS NMR

The carboxylic acid proton and the lithium coordination environments for precise and random Li-neutralized polyethylene acrylic acid P(E-AA) ionomers were explored using high speed solid-state 1H and 7Li MAS NMR. While the 7Li NMR revealed only a single Li coordination environment, the chemical shift temperature variation was dependent on the precise or random nature of the P(E-AA) ionomer. The 1H MAS NMR revealed two different carboxylic acid proton environments in these materials. By utilizing 1H-7Li rotational echo double resonance (REDOR) MAS NMR experiments, it was demonstrated that the proton environments correspond to different average 1H-7Li distances, with the majority of the protonated carboxylic acids having a close through space contact with the Li. Molecular dynamics simulations suggest that the shortest 1H-7Li distance corresponds to un-neutralized carboxylic acids directly involved in the coordination environment of Li clusters. These solid-state NMR results show that heterogeneous structural motifs need to be included when developing descriptions of these ionomer materials.

Molecular dynamics simulations of three protegrin-type antimicrobial peptides: interplay between charges at the termini, β-sheet structure and amphiphilic interactions.

We have carried out molecular dynamics simulations of the naturally occurring protegrin PG-1 peptide and two of its mutants, PC-9 and PC-13 in the presence of a dodecyl-phosphocholine (DPC) micelle. The effects of mutations that disrupt the β-sheet structure in the case of PC-9 and reduce the charge at the C-terminus in the case of PC-13 are analyzed. It is found that the surface-bound conformations of the peptides are severely affected by both mutations. PG-1 exhibits a conformation in which the C-terminus and the β-hairpin turn interact strongly with the micelle lipid head groups, while its N-terminal strand bends away from the micelle and resides in the aqueous region; PC-13 exhibits strong interactions with the micelle at its N-terminus as well as the β-hairpin turn region, while retaining a much more compact conformation than PG-1; PC-9 achieves a highly distorted conformation relative to the homologous PG-1 structure, which allows both its termini and the β-hairpin region to interact with the micelle. These significant differences observed as a result of seemingly minor mutations to the sequences of the three peptides are explained in terms of the interplay between residue charges, structural rigidity and amphiphilic interactions. Conservative inferences are made bridging these biophysical interactions and the pharmacological profiles of the peptides.

Multiscale Models of the Antimicrobial Peptide Protegrin-1 on Gram-Negative Bacteria Membranes

Antimicrobial peptides (AMPs) are naturally-occurring molecules that exhibit strong antibiotic properties against numerous infectious bacterial strains. Because of their unique mechanism of action, they have been touted as a potential source for novel antibiotic drugs. We present a summary of computational investigations in our lab aimed at understanding this unique mechanism of action, in particular the development of models that provide a quantitative connection between molecular-level biophysical phenomena and relevant biological effects. Our work is focused on protegrins, a potent class of AMPs that attack bacteria by associating with the bacterial membrane and forming transmembrane pores that facilitate the unrestricted transport of ions. Using fully atomistic molecular dynamics simulations, we have computed the thermodynamics of peptide-membrane association and insertion, as well as peptide aggregation. We also present a multi-scale analysis of the ion transport properties of protegrin pores, ranging from atomistic molecular dynamics simulations to mesoscale continuum models of single-pore electrodiffusion to models of transient ion transport from bacterial cells. Overall, this work provides a quantitative mechanistic description of the mechanism of action of protegrin antimicrobial peptides across multiple length and time scales.

Investigation of Changes in Tetracycline Repressor Binding upon Mutations in the Tetracycline Operator.

The tetracycline operon is an important gene network component, commonly used in synthetic biology applications because of its switch-like character. At the heart of this system is the highly specific interaction of the tet repressor protein (TetR) with its cognate DNA sequence (tetO). TetR binding on tetO practically stops expression of genes downstream of tetO by excluding RNA polymerase from binding the promoter and initiating transcription. Mutating the tetO sequence alters the strength of TetR–tetO binding and thus provides a tool to synthetic biologists to manipulate gene expression levels. We employ molecular dynamics (MD) simulations coupled with the free energy perturbation method to investigate the binding affinity of TetR to different tetO mutants. We also carry out in vivo tests in Escherichia coli for a series of promoters based on these mutants. We obtain reasonable agreement between experimental green fluorescent protein (GFP) repression levels and binding free energy differences computed from molecular simulations. In all cases, the wild-type tetO sequence yields the strongest TetR binding, which is observed both experimentally, in terms of GFP levels, and in simulation, in terms of free energy changes. Two of the four tetO mutants we tested yield relatively strong binding, whereas the other two mutants tend to be significantly weaker. The clustering and relative ranking of this subset of tetO mutants is generally consistent between our own experimental data, previous experiments with different systems and the free energy changes computed from our simulations. Overall, this work offers insights into an important synthetic biological system and demonstrates the potential, as well as limitations of molecular simulations to quantitatively explain biologically relevant behavior.

Effects of Morphology on Ion Transport in Ionomers for Energy Storage

Polymer electrolytes are essential elements of current and next-generation energy storage applications. An important class of polymer electrolytes is ionomers, in which one of the ions is covalently bound to the polymer backbone. Ionomers are currently used in fuel cells, and show extraordinary promise as solid electrolytes in batteries for transportation and portablepower applications. Solid electrolytes are desirable for a variety of reasons. A primary one is safety: the lack of solvent leads to fewer electrochemical reactions (e.g. with the electrodes) and the absence of flammable liquids. Solvent-free electrolytes allow for less packaging (and hence higher energy density batteries) and easier manufacture. Single-ion conductors such as ionomers also have the advantage of higher efficiency (high lithium transference numbers), since the anions are bound to the polymer backbone and the current is primarily due to the cations that actively participate in the electrochemical reactions. However, to date ionomeric materials do not have sufficiently high conductivities. Ion transport mechanisms in ionomers and their relation to molecular structure are poorly understood, although it is known that ion transport is coupled to polymeric motion and to the nanoscale morphology of ionic aggregates that often self-assemble in the polymer matrix. In this project we investigate the morphology and ion dynamics in ionomer melts using ab initio calculations, NMR experiments, and both atomistic and coarse-grained molecular dynamics simulations. We find that the aggregate morphology depends strongly on polymer architecture, neutralizing cation, and degree of neutralization. This morphology in turn affects the counterion dynamics; we find that counterions diffuse more quickly in systems with percolated aggregates than in systems with discrete aggregates.