Discovery of Potent and Orally Bioavailable Macrocyclic Peptide-Peptoid Hybrid CXCR7 Modulators.

The chemokine receptor CXCR7 is an attractive target for a variety of diseases. While several small-molecule modulators of CXCR7 have been reported, peptidic macrocycles may provide advantages in terms of potency, selectivity, and reduced off-target activity. We produced a series of peptidic macrocycles that incorporate an N-linked peptoid functionality where the peptoid group enabled us to explore side-chain diversity well beyond that of natural amino acids. At the same time, theoretical calculations and experimental assays were used to track and reduce the polarity while closely monitoring the physicochemical properties. This strategy led to the discovery of macrocyclic peptide-peptoid hybrids with high CXCR7 binding affinities (Ki < 100 nM) and measurable passive permeability (Papp > 5 × 10-6 cm/s). Moreover, bioactive peptide 25 (Ki = 9 nM) achieved oral bioavailability of 18% in rats, which was commensurate with the observed plasma clearance values upon intravenous administration.

Simple Predictive Models of Passive Membrane Permeability Incorporating Size-Dependent Membrane-Water Partition.

We investigate the relationship between passive permeability and molecular size, in the context of solubility-diffusion theory, using a diverse compound set with molecular weights ranging from 151 to 828, which have all been characterized in a consistent manner using the RRCK cell monolayer assay. Computationally, each compound was subjected to extensive conformational search and physics-based permeability prediction, and multiple linear regression analyses were subsequently performed to determine, empirically, the relative contributions of hydrophobicity and molecular size to passive permeation in the RRCK assay. Additional analyses of Log D and PAMPA data suggest that these measurements are not size selective, a possible reason for their sometimes weak correlation with cell-based permeability.

Peptide to Peptoid Substitutions Increase Cell Permeability in Cyclic Hexapeptides.

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.

Cell-permeable cyclic peptides from synthetic libraries inspired by natural products.

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.

Chemical-biological characterization of a cruzain inhibitor reveals a second target and a mammalian off-target.

Inhibition of the Trypanosoma cruzi cysteine protease cruzain has been proposed as a therapeutic approach for the treatment of Chagas' disease. Among the best-studied cruzain inhibitors to date is the vinylsulfone K777 (1), which has proven effective in animal models of Chagas' disease. Recent structure-activity studies aimed at addressing potential liabilities of 1 have now produced analogues such as N-[(2S)-1-[[(E,3S)-1-(benzenesulfonyl)-5-phenylpent-1-en-3-yl]amino]-3-(4-methylphenyl)-1-oxopropan-2-yl]pyridine-4-carboxamide (4), which is trypanocidal at ten-fold lower concentrations than for 1. We now find that the trypanocidal activity of 4 derives primarily from the inhibition of T. cruzi 14-α-demethylase (TcCYP51), a cytochrome P450 enzyme involved in the biosynthesis of ergosterol in the parasite. Compound 4 also inhibits mammalian CYP isoforms but is trypanocidal at concentrations below those required to significantly inhibit mammalian CYPs in vitro. A chemical-proteomics approach employing an activity-based probe derived from 1 was used to identify mammalian cathepsin B as a potentially important off-target of 1 and 4. Computational docking studies and the evaluation of truncated analogues of 4 reveal structural determinants for TcCYP51 binding, information that will be useful in further optimization of this new class of inhibitors.

Optimizing PK properties of cyclic peptides: the effect of side chain substitutions on permeability and clearance().

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.

Testing physical models of passive membrane permeation.

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.

Methyl effects on protein-ligand binding.

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 10 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.

On-resin N-methylation of cyclic peptides for discovery of orally bioavailable scaffolds.

Backbone N-methylation is common among peptide natural products and has a substantial 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 dependent on 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 1,024-member virtual library of N-methyl cyclic peptides. One of the most permeable compounds, a cyclic hexapeptide (molecular mass = 755 Da) with three N-methyl groups, showed an oral bioavailability of 28% in rat.

Synthesis and evaluation of selected key methyl ether derivatives of vancomycin aglycon.

A select series of methyl ether derivatives of vancomcyin aglycon were prepared and examined for antimicrobial activity against vancomycin-sensitive Staphylococcus aureus and vancomycin-resistant Enterococci faecalis as well as their binding affinity for D-Ala-D-Ala and D-Ala-D-Lac. The intent of the study was to elucidate the role selected key methyl groups may play in the improvement of the in vitro antimicrobial profile of the tetra methyl ether derivative of vancomycin aglycon against vancomycin-resistant Enterococci faecalis previously reported. In these studies, methodology for selective derivatization of the A-, B-, and D-ring was developed that defines the relative reactivity of the four phenols of vancomycin aglycon, providing a foundation for future efforts for site-directed modification of the vancomycin aglycon core.

A nonazole CYP51 inhibitor cures Chagas' disease in a mouse model of acute infection.

Chagas' disease, the leading cause of heart failure in Latin America, is caused by the kinetoplastid protozoan Trypanosoma cruzi. The sterols of T. cruzi resemble those of fungi, both in composition and in biosynthesis. Azole inhibitors of sterol 14alpha-demethylase (CYP51) successfully treat fungal infections in humans, and efforts to adapt the success of antifungal azoles posaconazole and ravuconazole as second-use agents for Chagas' disease are under way. However, to address concerns about the use of azoles for Chagas' disease, including drug resistance and cost, the rational design of nonazole CYP51 inhibitors can provide promising alternative drug chemotypes. We report the curative effect of the nonazole CYP51 inhibitor LP10 in an acute mouse model of T. cruzi infection. Mice treated with an oral dose of 40 mg LP10/kg of body weight twice a day (BID) for 30 days, initiated 24 h postinfection, showed no signs of acute disease and had histologically normal tissues after 6 months. A very stringent test of cure showed that 4/5 mice had negative PCR results for T. cruzi, and parasites were amplified by hemoculture in only two treated mice. These results compare favorably with those reported for posaconazole. Electron microscopy and gas chromatography-mass spectrometry (GC-MS) analysis of sterol composition confirmed that treatment with LP10 blocked the 14alpha-demethylation step and induced breakdown of parasite cell membranes, culminating in severe ultrastructural and morphological alterations and death of the clinically relevant amastigote stage of the parasite.

Structural characterization of CYP51 from Trypanosoma cruzi and Trypanosoma brucei bound to the antifungal drugs posaconazole and fluconazole.

BACKGROUND: Chagas 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 14alpha-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 FINDINGS: We have determined the crystal structures for CYP51 from T. cruzi (resolutions of 2.35 A and 2.27 A), and from the related pathogen T. brucei (resolutions of 2.7 A and 2.6 A), 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 degrees 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/SIGNIFICANCE: The 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.

Vancomycin analogs: Seeking improved binding of d-Ala-d-Ala and d-Ala-d-Lac peptides by side-chain and backbone modifications.

In order to seek vancomycin analogs with improved performance against VanA and VanB resistant bacterial strains, extensive computational investigations have been performed to examine the effects of side-chain and backbone modifications. Changes in binding affinities for tripeptide cell-wall precursor mimics, Ac(2)-l-Lys-d-Ala-d-Ala (3) and Ac(2)-l-Lys-d-Ala-d-Lac (4), with vancomycin analogs were computed with Monte Carlo/free energy perturbation (MC/FEP) calculations. Replacements of the 3-hydroxyl group in residue 7 with small alkyl or alkoxy groups, which improve contacts with the methyl side chain of the ligands'd-Ala residue, are predicted to be the most promising to enhance binding for both ligands. The previously reported amine backbone modification as in 5 is shown to complement the hydrophobic modifications for binding monoacetylated tripeptides. In addition, replacement of the hydroxyl groups in residues 5 and 7 by fluorine is computed to have negligible impact on binding the tripeptides, though it may be pharmacologically advantageous.

Vancomycin resistance: modeling backbone variants with D-Ala-D-Ala and D-Ala-D-Lac peptides.

To seek vancomycin analogs with broader antibacterial activity, effects of backbone modifications for the agylcon 2 on binding with D-Ala-D-Ala- and D-Ala-D-Lac-containing peptides were investigated by Monte Carlo/free energy perturbation (MC/FEP) calculations. The experimental trend in binding affinities for 2 with three tripeptides was well reproduced. Possible modifications of the peptide bond between residues 4 and 5 were then considered, specifically for conversion of the OCNH linkage to CH(2)NH(2)(+) (6), FCCH (7), HCCH (8), and HNCO (9). The MC/FEP results did not yield binding improvements for 7, 8, and 9, though the fluorovinyl replacement is relatively benign. The previously reported analog 6 remains as the only variant that exhibits improved affinity for the D-Ala-D-Lac sequence and acceptable affinity for the D-Ala-D-Ala sequence.

A Self-Consistent Space-Domain Decomposition Method for QM/MM Computations of Protein Electrostatic Potentials.

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.