Backbone N-methylation is common among peptide natural products and has a significant impact on both the physical properties and the conformational states of cyclic peptides. However, the specific impact of N-methylation on passive membrane diffusion in cyclic peptides has not been investigated systematically. Here we report a method for the selective, on-resin N-methylation of cyclic peptides to generate compounds with drug-like membrane permeability and oral bioavailability. The selectivity and degree of N-methylation of the cyclic peptide was determined by backbone stereochemistry, suggesting that conformation dictates the regiochemistry of the N-methylation reaction. The permeabilities of the N-methyl variants were corroborated by computational studies on a 1024-member virtual library of N-methyl cyclic peptides. One of the most permeable compounds, a cyclic hexapeptide (MW = 755) with three N-methyl groups, showed an oral bioavailability of 28% in rat.

The effects of addition of a methyl group to a lead compound on biological activity are examined. A literature analysis of >2000 cases reveals that an activity boost of a factor of ten or more is found with an 8% frequency, and a 100-fold boost is a 1 in 200 event. Four cases in the latter category are analyzed in depth to elucidate any unusual aspects of the protein-ligand binding, distribution of water molecules, and changes in conformational energetics. The analyses include Monte Carlo/Free Energy Perturbation (MC/FEP) calculations for methyl replacements in inhibitor series for p38α MAP kinase, ACK1, PTP1B, and thrombin. Methyl substitutions ortho to an aryl ring can be particularly effective at improving activity by inducing a propitious conformational change. The greatest improvements in activity arise from coupling the conformational gain with the burial of the methyl group in a hydrophobic region of the protein.

Drug design efforts are turning to a new generation of therapeutic targets, such as protein-protein interactions (PPIs), that had previously been considered "undruggable" by typical small molecules. There is an emerging view that accessing these targets will require molecules that are larger and more complex than typical small molecule drugs. Here, we present a methodology for the discovery of geometrically diverse, membrane permeable cyclic peptide scaffolds based on the synthesis and permeability screening of a combinatorial library, followed by deconvolution of membrane-permeable scaffolds to identify cyclic peptides with good to excellent passive cell permeabilities. We use a combination of experimental and computational approaches to investigate structure-permeability relationships in one of these scaffolds, and uncover structural and conformational factors that govern passive membrane diffusion in a related set of cyclic peptide diastereomers. Further, we investigate the dependency of permeability on side-chain identity of one of these scaffolds through single-point diversifications to show the adaptability of these scaffolds toward development of permeability-biased libraries suitable for bioactivity screens. Overall, our results demonstrate that many novel, cell permeable scaffolds exist beyond those found in extant natural products, and that such scaffolds can be rapidly identified using a combination of synthesis and deconvolution which can, in principle, be applied to any type of macrocyclic template.

The biophysical basis of passive membrane permeability is well understood, but most methods for predicting membrane permeability in the context of drug design are based on statistical relationships that indirectly capture the key physical aspects. Here, we investigate molecular mechanics-based models of passive membrane permeability and evaluate their performance against different types of experimental data, including parallel artificial membrane permeability assays (PAMPA), cell-based assays, in vivo measurements, and other in silico predictions. The experimental data sets we use in these tests are diverse, including peptidomimetics, congeneric series, and diverse FDA approved drugs. The physical models are not specifically trained for any of these data sets; rather, input parameters are based on standard molecular mechanics force fields, such as partial charges, and an implicit solvent model. A systematic approach is taken to analyze the contribution from each component in the physics-based permeability model. A primary factor in determining rates of passive membrane permeation is the conformation-dependent free energy of desolvating the molecule, and this measure alone provides good agreement with experimental permeability measurements in many cases. Other factors that improve agreement with experimental data include deionization and estimates of entropy losses of the ligand and the membrane, which lead to size-dependence of the permeation rate.

A series of cyclic peptides were designed and prepared to investigate the physicochemical properties that affect oral bioavailabilty of this chemotype in rats. In particular, the ionization state of the peptide was examined by the incorporation of naturally occurring amino acid residues that are charged in differing regions of the gut. In addition, data was generated in a variety of in vitro assays and the usefulness of this data in predicting the subsequent oral bioavailability observed in the rat is discussed.

BackgroundChagas Disease is the leading cause of heart failure in Latin America. Current drug therapy is limited by issues of both efficacy and severe side effects. Trypansoma cruzi, the protozoan agent of Chagas Disease, is closely related to two other major global pathogens, Leishmania spp., responsible for leishmaniasis, and Trypansoma brucei, the causative agent of African Sleeping Sickness. Both T. cruzi and Leishmania parasites have an essential requirement for ergosterol, and are thus vulnerable to inhibitors of sterol 14α-demethylase (CYP51), which catalyzes the conversion of lanosterol to ergosterol. Clinically employed anti-fungal azoles inhibit ergosterol biosynthesis in fungi, and specific azoles are also effective against both Trypanosoma and Leishmania parasites. However, modification of azoles to enhance efficacy and circumvent potential drug resistance has been problematic for both parasitic and fungal infections due to the lack of structural insights into drug binding.Methodology/Principal FindingsWe have determined the crystal structures for CYP51 from T. cruzi (resolutions of 2.35 Å and 2.27 Å), and from the related pathogen T. brucei (resolutions of 2.7 Å and 2.6 Å), co-crystallized with the antifungal drugs fluconazole and posaconazole. Remarkably, both drugs adopt multiple conformations when binding the target. The fluconazole 2,4-difluorophenyl ring flips 180° depending on the H-bonding interactions with the BC-loop. The terminus of the long functional tail group of posaconazole is bound loosely in the mouth of the hydrophobic substrate binding tunnel, suggesting that the major contribution of the tail to drug efficacy is for pharmacokinetics rather than in interactions with the target.Conclusions/SignificanceThe structures provide new insights into binding of azoles to CYP51 and mechanisms of potential drug resistance. Our studies define in structural detail the CYP51 therapeutic target in T. cruzi, and offer a starting point for rationally designed anti-Chagasic drugs with improved efficacy and reduced toxicity.

This paper introduces a self-consistent computational protocol for modeling protein electrostatic potentials according to static point-charge model distributions. The protocol involves a simple space-domain decomposition scheme where individual molecular domains are modeled as Quantum-Mechanical (QM) layers embedded in the otherwise classical Molecular-Mechanics (MM) protein environment. ElectroStatic-Potential (ESP) atomic charges of the constituent molecular domains are computed, to account for mutual polarization effects, and iterated until obtaining a self-consistent point-charge model of the protein electrostatic potential. The novel protocol achieves quantitative agreement with full QM calculations in the description of electrostatic potentials of small polypeptides where polarization effects are significant, showing a remarkable improvement relative to the corresponding electrostatic potentials obtained with popular MM force fields. The capabilities of the method are demonstrated in several applications, including calculations of the electrostatic potential in the potassium channel protein and the description of protein-protein electrostatic interactions.

The effect of peptide-to-peptoid substitutions on the passive membrane permeability of an N-methylated cyclic hexapeptide is examined. In general, substitutions maintained permeability but increased conformational heterogeneity. Diversification with nonproteinogenic side chains increased permeability up to 3-fold. Additionally, the conformational impact of peptoid substitutions within a β-turn are explored. Based on these results, the strategic incorporation of peptoid residues into cyclic peptides can maintain or improve cell permeability, while increasing access to diverse side-chain functionality.