In Silico Screening for Solid Dispersions: The Trouble with Solubility Parameters and χFH.

The problem of predicting small molecule-polymer compatibility is relevant to many areas of chemistry and pharmaceutical science but particularly drug delivery. Computational methods based on Hildebrand and Hansen solubility parameters, and the estimation of the Flory-Huggins parameter, χ, have proliferated across the literature. Focusing on the need to develop amorphous solid dispersions to improve the bioavailability of poorly soluble drug candidates, an innovative, high-throughput 2D printing method has been employed to rapidly assess the compatibility of 54 drug-polymer pairings (nine drug compounds in six polymers). In this study, the first systematic assessment of the in silico methods for this application, neither the solubility parameter approach nor the calculated χ, correctly predicted drug-polymer compatibility. The theoretical limitations of the solubility parameter approach are discussed and used to explain why this approach is fundamentally unsuitable for predicting polymer-drug interactions. Examination of the original sources describing the method for calculating χ shows that only the enthalpic contributions to the term have been included, and the corrective entropic term is absent. The development and application of new in silico techniques, that consider all parts of the free energy of mixing, are needed in order to usefully predict small molecule-polymer compatibility and to realize the ambition of a drug-polymer screening method.

Cooperativity and site selectivity in the ileal lipid binding protein.

The ileal lipid binding protein (ILBP or I-BABP) binds bile salts with positive cooperativity and has unusual site selectivity, whereby cholic acid binds preferentially in one site and chenodeoxycholic in another, despite both sites having an affinity for both ligands and the ligands only differing by a single hydroxyl group. Previous studies of the human variant have assumed that the ligand/protein binding ratio is 2:1, but we show, using electrospray ionization mass spectroscopy, that human ILBP binds bile acids with a 3:1 ratio, even at low protein and ligand concentrations. Docking calculations and molecular dynamics (MD) simulations identify an allosterically active binding site on the protein exterior that induces a change from a closed conformation to an open one, characterized by a movement of one of the α-helices by ~10° with respect to the ß-clam shell. Additional independent MD simulations of several hundred nanoseconds implicate the change between conformations in the mechanisms of both cooperativity and ligand site selectivity.

Stereoselective disulfide formation stabilizes the local peptide conformation in nisin mimics.

Nisin is a polymacrocyclic peptide antimicrobial with high activity against Gram-positive bacteria. Lanthionine and methyllanthionine bridges, closing the macrocycles, are stabilized by thioether bonds, formed between cysteines and dehydrated serine or threonine. The role of polypeptide backbone conformation in the formation of macrocycles A and B within cysteine mutants of nisin residues 1−12 is investigated here by molecular dynamics simulations. Enantiomeric combinational space of Cys3 and Cys7 and of Cys8 and Cys11 is examined for the preference of disulfide bond formation over helical turn formation within this region. A clear preference for spontaneous disulfide formation and closure of rings 3,7 and 8,11 is demonstrated for the D-Cys3, D-Cys7, L-Cys8, L-Cys11 nisin homologue, while interlinked rings A and B are obtained through disulfide bridges between L-Cys3 and D-Cys8 and between D-Cys7 and D-Cys11. This study offers a simple designer approach to solid phase synthesis of macrocyclic peptides and lantibiotic analogues.

Docking and molecular dynamics simulations of the ternary complex nisin2:lipid II.

Lanthionine antibiotics are an important class of naturally-occurring antimicrobial peptides. The best-known, nisin, is a commercial food preservative. However, structural and mechanistic details on nisin-lipid II membrane complexes are currently lacking. Recently, we have developed empirical force-field parameters to model lantibiotics. Docking and molecular dynamics (MD) simulations have been used to study the nisin2:lipid II complex in bacterial membranes, which has been put forward as the building block of nisin/lipid II binary membrane pores. An Ile1Trp mutation of the N-terminus of nisin has been modelled and docked onto lipid II models; the computed binding affinity increased compared to wild-type. Wild-type nisin was also docked onto three different lipid II structures and a stable 2:1 nisin:lipid II complex formed. This complex was inserted into a membrane. Six independent MD simulations revealed key interactions in the complex, specifically the N-terminal engagement of nisin with lipid II at the pyrophosphate and C-terminus of the pentapeptide chain. Nisin(2) inserts into the membrane and we propose this as the first step in pore formation, mediated by the nisin N-terminus-lipid II pentapeptide hydrogen bond. The lipid II undecaprenyl chain adopted different conformations in the presence of nisin, which may also have implications for pore formation.