Dimerization of Protegrin-1 in Different Environments

The dimerization of the cationic β-hairpin antimicrobial peptide protegrin-1 (PG1) is investigated in three different environments: water, the surface of a lipid bilayer membrane, and the core of the membrane. PG1 is known to kill bacteria by forming oligomeric membrane pores, which permeabilize the cells. PG1 dimers are found in two distinct, parallel and antiparallel, conformations, known as important intermediate structural units of the active pore oligomers. What is not clear is the sequence of events from PG1 monomers in solution to pores inside membranes. The step we focus on in this work is the dimerization of PG1. In particular, we are interested in determining where PG1 dimerization is most favorable. We use extensive molecular dynamics simulations to determine the potential of mean force as a function of distance between two PG1 monomers in the aqueous subphase, the surface of model lipid bilayers and the interior of these bilayers. We investigate the two known distinct modes of dimerization that result in either a parallel or an antiparallel β-sheet orientation. The model bilayer membranes are composed of anionic palmitoyl-oleoyl-phosphatidylglycerol (POPG) and palmitoyl-oleoyl-phosphatidylethanolamine (POPE) in a 1:3 ratio (POPG:POPE). We find the parallel PG1 dimer association to be more favorable than the antiparallel one in water and inside the membrane. However, we observe that the antiparallel PG1 β-sheet dimer conformation is somewhat more stable than the parallel dimer association at the surface of the membrane. We explore the role of hydrogen bonds and ionic bridges in peptide dimerization in the three environments. Detailed knowledge of how networks of ionic bridges and hydrogen bonds contribute to peptide stability is essential for the purpose of understanding the mechanism of action for membrane-active peptides as well as for designing peptides which can modulate membrane properties. The findings are suggestive of the dominant pathways leading from individual PG1 molecules in solution to functional pores in bacterial membranes

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

Assessment of Solvated Interaction Energy Function for Ranking Antibody–Antigen Binding Affinities

Affinity modulation of antibodies and antibody fragments of therapeutic value is often required in order to improve their clinical efficacies. Virtual affinity maturation has the potential to quickly focus on the critical hotspot residues without the combinatorial explosion problem of conventional display and library approaches. However, this requires a binding affinity scoring function that is capable of ranking single-point mutations of a starting antibody. We focus here on assessing the solvated interaction energy (SIE) function that was originally developed for and is widely applied to scoring of protein–ligand binding affinities. To this end, we assembled a structure–function data set called Single-Point Mutant Antibody Binding (SiPMAB) comprising several antibody–antigen systems suitable for this assessment, i.e., based on high-resolution crystal structures for the parent antibodies and coupled with high-quality binding affinity measurements for sets of single-point antibody mutants in each system. Using this data set, we tested the SIE function with several mutation protocols based on the popular methods SCWRL, Rosetta, and FoldX. We found that the SIE function coupled with a protocol limited to sampling only the mutated side chain can reasonably predict relative binding affinities with a Spearman rank-order correlation coefficient of about 0.6, outperforming more aggressive sampling protocols. Importantly, this performance is maintained for each of the seven system-specific component subsets as well as for other relevant subsets including non-alanine and charge-altering mutations. The transferability and enrichment in affinity-improving mutants can be further enhanced using consensus ranking over multiple methods, including the SIE, Talaris, and FOLDEF energy functions. The knowledge gained from this study can lead to successful prospective applications of virtual affinity maturation

Assessment of Solvated Interaction Energy Function for Ranking Antibody–Antigen Binding Affinities

Affinity modulation of antibodies and antibody fragments of therapeutic value is often required in order to improve their clinical efficacies. Virtual affinity maturation has the potential to quickly focus on the critical hotspot residues without the combinatorial explosion problem of conventional display and library approaches. However, this requires a binding affinity scoring function that is capable of ranking single-point mutations of a starting antibody. We focus here on assessing the solvated interaction energy (SIE) function that was originally developed for and is widely applied to scoring of protein–ligand binding affinities. To this end, we assembled a structure–function data set called Single-Point Mutant Antibody Binding (SiPMAB) comprising several antibody–antigen systems suitable for this assessment, i.e., based on high-resolution crystal structures for the parent antibodies and coupled with high-quality binding affinity measurements for sets of single-point antibody mutants in each system. Using this data set, we tested the SIE function with several mutation protocols based on the popular methods SCWRL, Rosetta, and FoldX. We found that the SIE function coupled with a protocol limited to sampling only the mutated side chain can reasonably predict relative binding affinities with a Spearman rank-order correlation coefficient of about 0.6, outperforming more aggressive sampling protocols. Importantly, this performance is maintained for each of the seven system-specific component subsets as well as for other relevant subsets including non-alanine and charge-altering mutations. The transferability and enrichment in affinity-improving mutants can be further enhanced using consensus ranking over multiple methods, including the SIE, Talaris, and FOLDEF energy functions. The knowledge gained from this study can lead to successful prospective applications of virtual affinity maturation

Assessment of Solvated Interaction Energy Function for Ranking Antibody–Antigen Binding Affinities

Affinity modulation of antibodies and antibody fragments of therapeutic value is often required in order to improve their clinical efficacies. Virtual affinity maturation has the potential to quickly focus on the critical hotspot residues without the combinatorial explosion problem of conventional display and library approaches. However, this requires a binding affinity scoring function that is capable of ranking single-point mutations of a starting antibody. We focus here on assessing the solvated interaction energy (SIE) function that was originally developed for and is widely applied to scoring of protein–ligand binding affinities. To this end, we assembled a structure–function data set called Single-Point Mutant Antibody Binding (SiPMAB) comprising several antibody–antigen systems suitable for this assessment, i.e., based on high-resolution crystal structures for the parent antibodies and coupled with high-quality binding affinity measurements for sets of single-point antibody mutants in each system. Using this data set, we tested the SIE function with several mutation protocols based on the popular methods SCWRL, Rosetta, and FoldX. We found that the SIE function coupled with a protocol limited to sampling only the mutated side chain can reasonably predict relative binding affinities with a Spearman rank-order correlation coefficient of about 0.6, outperforming more aggressive sampling protocols. Importantly, this performance is maintained for each of the seven system-specific component subsets as well as for other relevant subsets including non-alanine and charge-altering mutations. The transferability and enrichment in affinity-improving mutants can be further enhanced using consensus ranking over multiple methods, including the SIE, Talaris, and FOLDEF energy functions. The knowledge gained from this study can lead to successful prospective applications of virtual affinity maturation