Search
Close this search box.

by Leung SS, Mijalkovic J, Borrelli K, Jacobson MP

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.

June 25, 2012
Journal of Chemical Information and Modeling
Read online article

by Rafi SB, Hearn BR, Vedantham P, Jacobson MP, Renslo AR

We evaluate experimentally and computationally the membrane permeability of matched sets of peptidic small molecules bearing natural or bioisosteric unnatural amino acids. We find that the intentional introduction of hydrogen bond acceptor-donor pairs in such molecules can improve membrane permeability while retaining or improving other favorable drug-like properties. We employ an all-atom force field based method to calculate changes in free energy associated with the transfer of the peptidic molecules from water to membrane. This computational method correctly predicts rank order experimental permeability trends within congeneric series and is much more predictive than calculations (e.g., clogP) that do not consider three-dimensional conformation.

April 12, 2012
Journal of Medicinal Chemistry
Read online article

by White TR, Renzelman CM, Rand AC, Rezai T, McEwen CM, Gelev VM, Turner RA, Linington RG, Leung SS, Kalgutkar AS, Bauman JN, Zhang Y, Liras S, Price DA, Mathiowetz AM, Jacobson MP, Lokey RS

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.

September 25, 2011
Nature Chemical Biology
Read online article

by Elena Dolghih, Clifford Bryant, Adan R. Renslo, Matthew P. Jacobson

P-glycoprotein (P-gp) is an ATP-dependent transport protein that is selectively expressed at entry points of xenobiotics where, acting as an efflux pump, it prevents their entering sensitive organs. The protein also plays a key role in the absorption and blood-brain barrier penetration of many drugs, while its overexpression in cancer cells has been linked to multidrug resistance in tumors. The recent publication of the mouse P-gp crystal structure revealed a large and hydrophobic binding cavity with no clearly defined sub-sites that supports an “induced-fit” ligand binding model. We employed flexible receptor docking to develop a new prediction algorithm for P-gp binding specificity. We tested the ability of this method to differentiate between binders and nonbinders of P-gp using consistently measured experimental data from P-gp efflux and calcein-inhibition assays. We also subjected the model to a blind test on a series of peptidic cysteine protease inhibitors, confirming the ability to predict compounds more likely to be P-gp substrates. Finally, we used the method to predict cellular metabolites that may be P-gp substrates. Overall, our results suggest that many P-gp substrates bind deeper in the cavity than the cyclic peptide in the crystal structure and that specificity in P-gp is better understood in terms of physicochemical properties of the ligands (and the binding site), rather than being defined by specific sub-sites.

June 23, 2011
PLOS – Computational Biology
Read online article

by Combs DJ, Lokey RS

Peptoids (N-substituted polyglycines) represent a class of bioinspired oligomers that have unique physical and structural properties. Here we report the construction of “extended peptoids” based on aromatic building blocks, in which the N-alkylaminoacetyl group of the peptoid backbone has been replaced by an N-alkylaminomethylbenzoyl spacer. Both meta- and para-bromomethylbenzoic acids were synthesized, providing access to a new class of peptoids. Further, inclusion of hydrophilic side chains confers water solubility to these compounds, showing that, like simple peptoids, extended peptoids add an extra dimension to synthetic polyamide oligomers with potential application in a variety of biological contexts.

April 9, 2007
Tetrahedron Letters
Read online article

by Rezai T, Yu B, Millhauser GL, Jacobson MP, Lokey RS

Little is known about the effect of conformation on passive membrane diffusion rates in small molecules. Evidence suggests that intramolecular hydrogen bonding may play a role by reducing the energetic cost of desolvating hydrogen bond donors, especially amide N-H groups. We set out to test this hypothesis by investigating the passive membrane diffusion characteristics of a series of cyclic peptide diastereomers based on the sequence cyclo[Leu-Leu-Leu-Leu-Pro-Tyr]. We identified two cyclic hexapeptide diastereomers based on this sequence, whose membrane diffusion rates differed by nearly two log units. Results of solution NMR studies and hydrogen/deuterium (H/D) exchange experiments showed that membrane diffusion rates correlated with the degree of intramolecular hydrogen bonding and H/D exchange rates. The most permeable diastereomer, cyclo[d-Leu-d-Leu-Leu-d-Leu-Pro-Tyr] (1), exhibited a passive membrane diffusion rate comparable to that of the orally available drug cyclosporine A.

March 1, 2006
Journal of the American Chemical Society
Read online article