An atomistic model of passive membrane permeability: application to a series of FDA approved drugs

We apply an atomistic model of passive membrane permeability to a series of weakly basic drugs. The computational model uses conformational sampling in combination with an all-atom force field and implicit solvent model to estimate relative passive membrane permeabilities. The model does not require the use of training data for rank-ordering compounds, and as such represents a different approach from the more commonly employed QSPR models. We compare the computational results to previously published experimental PAMPA and Caco-2 permeabilities.     

Improving on Nature: Making a Cyclic Heptapeptide Orally Bioavailable

The use of peptides in medicine is limited by low membrane permeability, metabolic instability, high clearance, and negligible oral bioavailability. The prediction of oral bioavailability of drugs relies on physicochemical properties that favor passive permeability and oxidative metabolic stability, but these may not be useful for peptides. Here we investigate effects of heterocyclic constraints, intramolecular hydrogen bonds, and side chains on the oral bioavailability of cyclic heptapeptides. NMR‐derived structures, amide H–D exchange rates, and temperature‐dependent chemical shifts showed that the combination of rigidification, stronger hydrogen bonds, and solvent shielding by branched side chains enhances the oral bioavailability of cyclic heptapeptides in rats without the need for N‐methylation.     

cis‐Peptide Bonds: A Key for Intestinal Permeability of Peptides?

Recent structural studies on libraries of cyclic hexapeptides led to the identification of common backbone conformations that may be instrumental to the oral availability of peptides. Furthermore, the observation of differential Caco‐2 permeabilities of enantiomeric pairs of some of these peptides strongly supports the concept of conformational specificity driven uptake and also suggests a pivotal role of carrier‐mediated pathways for peptide transport, especially for scaffolds of polar nature. This work presents investigations on the Caco‐2 and PAMPA permeability profiles of 13 selected N‐methylated cyclic pentaalanine peptides derived from the basic cyclo(‐D ‐Ala‐Ala₄‐) template. These molecules generally showed moderate to low transport in intestinal epithelia with a few of them exhibiting a Caco‐2 permeability equal to or slightly higher than that of mannitol, a marker for paracellular permeability. We identified that the majority of the permeable cyclic penta‐ and hexapeptides possess an N‐methylated cis‐peptide bond, a structural feature that is also present in the orally available peptides cyclosporine A and the tri‐N‐methylated analogue of the Veber–Hirschmann peptide. Based on these observations it appears that the presence of N‐methylated cis‐peptide bonds at certain locations may promote the intestinal permeability of peptides through a suitable conformational preorganization.     

Enantiomeric Cyclic Peptides with Different Caco‐2 Permeability Suggest Carrier‐Mediated Transport

Recently, oral absorption of cyclic hexapeptides was improved by N‐methylation of their backbone amides. However, the number and position of N‐methylations or of solvent exposed NHs did not correlate to intestinal permeability, measured in a Caco‐2 model. In this study, we investigate enantiomeric pairs of three polar and two lipophilic peptides to demonstrate the participation of carrier‐mediated transporters. As expected, all the enantiomeric peptides exhibited identical lipophilicity (logD_7.4) and passive transcellular permeability determined by the parallel artificial membrane permeability assay (PAMPA). However, the enantiomeric polar peptides exhibited different Caco‐2 permeability (P_app) in both directions a–b and b–a. The same trend was observed for one of the lipophilic peptide, whereas the second lipophilic enantiomer pair showed identical Caco‐2 permeability (within the errors). These findings provide the first evidence for the involvement of carrier‐mediated transport for peptides, especially for those of polar nature.     

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.     

Solution Conformations Explain the Chameleonic Behaviour of Macrocyclic Drugs

It has been hypothesised that drugs in the chemical space “beyond the rule of 5” (bRo5) must behave as molecular chameleons to combine otherwise conflicting properties, including aqueous solubility, cell permeability and target binding. Evidence for this has, however, been limited to the cyclic peptide cyclosporine A. Herein, we show that the non-peptidic and macrocyclic drugs roxithromycin, telithromycin and spiramycin behave as molecular chameleons, with rifampicin showing a less pronounced behaviour. In particular roxithromycin, telithromycin and spiramycin display a marked, yet limited flexibility and populate significantly less polar and more compact conformational ensembles in an apolar than in a polar environment. In addition to balancing of membrane permeability and aqueous solubility, this flexibility also allows binding to targets that vary in structure between species. The drugs’ passive cell permeability correlates to their 3D polar surface area and corroborate two theoretical models for permeability, developed for cyclic peptides. We conclude that molecular chameleonicity should be incorporated in the design of orally administered drugs in the bRo5 space.     

Synthesis and conformational studies of gamma-aminoxy peptides

We have synthesized a series of gamma-aminoxy acids, including unsubstituted and gamma4-Ph-, gamma4-alkyl-, and gamma(3,4)-cyclohexyl-substituted systems. Coupling of these monomers to oligomers can be realized using EDCI/HOBt (or HOAt) as the coupling agent. gamma-Aminoxy peptides can form 10-membered-ring intramolecular hydrogen bonds-so-called "gamma N-O turns"-between adjacent residues, the extent of which is controlled by the nature of the side chain of each gamma-aminoxy acid residue, increasing from the unsubstituted gamma-aminoxy peptide to the gamma4-alkyl aminoxy peptides to the gamma4-phenyl- and gamma(3,4)-cyclohexyl-substituted aminoxy peptides. The presence of two consecutive homochiral 10-membered-ring intramolecular hydrogen bonds leads to the formation of a novel helical structure. Theoretical studies on a series of model peptides rationalize very well the experimentally observed conformational features of these gamma-aminoxy peptides.     

Intramolecular Hydrogen Bonding Enables a Zwitterionic Mechanism for Macrocyclic Peptide Formation: Computational Mechanistic Studies of CyClick Chemistry

Macrocyclic peptides have become increasingly important in the pharmaceutical industry. We present a detailed computational investigation of the reaction mechanism of the recently developed “CyClick” chemistry to selectively form imidazolidinone cyclic peptides from linear peptide aldehydes, without using catalysts or directing groups (Angew. Chem. Int. Ed. 2019 , 58, 19073–19080). We conducted computational mechanistic to investigate the effects of intramolecular hydrogen bonds (IMHBs) in promoting a kinetically facile zwitterionic mechanism in “CyClick” of pentapeptide aldehyde AFGPA. Our DFT calculations highlighted the importance of IMHB in pre-organization of the resting state, stabilization of the zwitterion intermediate, and the control of the product stereoselectivity. Furthermore, we have also identified that the low ring strain energy promotes the “CyClick” of hexapeptide aldehyde AAGPFA to form a thermodynamically more stable 15+5 imidazolidinone cyclic peptide product. In contrast, large ring strain energy suppresses “CyClick” reactivity of tetra peptide aldehyde AFPA from forming the 9+5 imidazolidinone cyclic peptide product.     

Influence of hydrogen bonds between sidechains and cooperativity on self‐assembly of the cyclopeptides

The self‐assembly of four cyclic D,L‐octapeptides, [‐(D‐Ala‐Gln)₄‐], [‐(D‐Val‐Gln)₄‐], [‐(D‐Leu‐Gln)₄‐], and [‐(D‐Phe‐Gln)₄‐], was investigated on the theory level in detail. Based on these cyclic peptides, which contain L‐Gln residues and possess C₄ symmetry, a series of oligomers were constructed according to different stacking modes as well as interaction patterns. We employed the semiempirical molecular orbital method AM1 to optimize the structures of all the oligomers, some of which were further studied using density functional method B3PW91/6‐31G to calculate the interaction energies. The studies indicate that when these cyclopeptides aggregate to form oligomers, or even nanotubes, four more hydrogen bonds could form between the sidechains of L‐Gln residues in addition to eight hydrogen bonds formed between the backbones of adjacent two cyclic peptides, a result that would clearly affect the self‐assembling process of cyclic peptides. The main effects can be summarized as follows. First, the dimers of these cyclic peptides with C₄ symmetry are more stable than those with D₄ symmetry due to their additional H‐bonds between Gln sidechains. Second, for the self‐assembly of the cyclopeptides, there is a competition between parallel and antiparallel stacking modes in lower oligomers such as dimers. However, with an increasing degree of oligomerization, energetically there is an increased possibility for the cyclic peptides to take the parallel stacking mode in assembly. Finally, the synergetic effect of weak interactions is the fundamental driving force for cyclic peptides to form stable nanotubes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Unexpectedly fast cis/trans isomerization of Xaa-Pro peptide bonds in disulfide-constrained cyclic peptides

Acyclic dithiol and cyclic disulfide forms of the peptides Ac-Cys-Pro-Xaa-Cys-NH2 (Xaa = Phe, His, Tyr, Gly, and Thr) and Ac-Cys-Gly-Pro-Cys-NH2 and the peptide Ac-Ala-Gly-Pro-Ala-NH2 were synthesized and characterized by mass spectrometry and NMR spectroscopy. Rate constants kct and ktc for cis-to-trans and trans-to-cis isomerization, respectively, across the Cys-Pro or Gly-Pro peptide bonds were determined by magnetization transfer NMR techniques over a range of temperatures, and activation parameters were derived from the temperature dependence of the rate constants. It was found that constraints imposed by the disulfide bond confer an unexpected rate enhancement for cis/trans isomerization, ranging from a factor of 2 to 13. It is proposed that the rate enhancements are a result of an intramolecular catalysis mechanism in which the NH proton of the Pro-Xaa peptide bond hydrogen bonds to the proline nitrogen in the transition state. The peptides Ac-Cys-Pro-Xaa-Cys-NH2 and Ac-Cys-Gly-Pro-Cys-NH2 are model compounds for proline-containing active sites of the thioredoxin superfamily of oxidoreductase enzymes; the results suggest that the backbones of the active sites of the oxidized form of these enzymes may have unusual conformational flexibility.     

Structure prediction of cyclic peptides by molecular dynamics + machine learning

Recent computational methods have made strides in discovering well-structured cyclic peptides that preferentially populate a single conformation. However, many successful cyclic-peptide therapeutics adopt multiple conformations in solution. In fact, the chameleonic properties of some cyclic peptides are likely responsible for their high cell membrane permeability. Thus, we require the ability to predict complete structural ensembles for cyclic peptides, including the majority of cyclic peptides that have broad structural ensembles, to significantly improve our ability to rationally design cyclic-peptide therapeutics. Here, we introduce the idea of using molecular dynamics simulation results to train machine learning models to enable efficient structure prediction for cyclic peptides. Using molecular dynamics simulation results for several hundred cyclic pentapeptides as the training datasets, we developed machine-learning models that can provide molecular dynamics simulation-quality predictions of structural ensembles for all the hundreds of thousands of sequences in the entire sequence space. The prediction for each individual cyclic peptide can be made using less than 1 second of computation time. Even for the most challenging classes of poorly structured cyclic peptides with broad conformational ensembles, our predictions were similar to those one would normally obtain only after running multiple days of explicit-solvent molecular dynamics simulations. The resulting method, termed StrEAMM (Structural Ensembles Achieved by Molecular Dynamics and Machine Learning), is the first technique capable of efficiently predicting complete structural ensembles of cyclic peptides without relying on additional molecular dynamics simulations, constituting a seven-order-of-magnitude improvement in speed while retaining the same accuracy as explicit-solvent simulations.

The StrEAMM method enables predicting the structural ensembles of cyclic peptides that adopt multiple conformations in solution.     

Water clusters on Cu(110): chain versus cyclic structures

Water clusters are assembled and imaged on Cu(110) by using a scanning tunneling microscope. Water molecules are arranged along the Cu row to form "ferroelectric" zigzag chains of trimer to hexamer. The trimer prefers the chain form to a cyclic one in spite of the reduced number of hydrogen bonds, highlighting the crucial role of the water-substrate interaction in the clustering of adsorbed water molecules. On the other hand, the cyclic form with maximal hydrogen bonds becomes more favorable for the tetramer, indicating the crossover from chain to cyclic configurations as the constituent number increases.     

Solubility and Permeability Enhancement of Oleanolic Acid by Solid Dispersion in Poloxamers and γ-CD

Oleanolic acid (OA) is a pentacyclic triterpenoid widely found in the Oleaceae family, and it represents 3.5% of the dry weight of olive leaves. OA has many pharmacological activities, such as hepatoprotection, anti-inflammatory, anti-oxidant, anti-diabetic, anti-tumor, and anti-microbic activities. Its therapeutic application is limited by its poor water solubility, bioavailability, and permeability. In this study, solid dispersions (SDs) were developed to overcome these OA limitations. Solubility studies were conducted to evaluate different hydrophilic polymers, drug-to-polymer ratios, and preparation methods. Poloxamer 188, Poloxamer 407, and γ-CD exhibited the highest increases in terms of OA solubility, regardless of the method of preparation. Binary systems were characterized using differential scanning calorimetry (DSC), X-ray diffraction (XRPD), and Fourier transform infrared spectroscopy (FTIR). In addition, pure compounds and SDs were analyzed using scanning electron microscopy (SEM) in order to observe both the morphology and the particle surface. In vitro dissolution studies were performed for P407, P188, and γ-CD SDs. Preparation using the solvent evaporation method (SEM) produced the highest increase in the dissolution profiles of all three polymers with respect to the OA solution. Finally, the effect of SDs on OA permeability was evaluated with an in vitro parallel artificial membrane permeability assay (PAMPA). The formulation improved passive permeation across the simulated barrier due to OA increased solubility. The dissolution and PAMPA results indicate that the amorphization of OA by SD preparation could be a useful method to enhance its oral absorption, and it is also applicable on an industrial scale.     

Solubility,Permeability, Anti-Inflammatory Action and In Vivo Pharmacokinetic Properties of Several Mechanochemically Obtained Pharmaceutical Solid Dispersions of Nimesulide

Nimesulide (NIM, N-(4-nitro-2-phenoxyphenyl)methanesulfonamide) is a relatively new nonsteroidal anti-inflammatory analgesic drug. It is practically insoluble in water (<0.02 mg/mL). This very poor aqueous solubility of the drug may lead to low bioavailability. The objective of the present study was to investigate the possibility of improving the solubility and the bioavailability of NIM via complexation with polysaccharide arabinogalactan (AG), disodium salt of glycyrrhizic acid (Na₂GA), hydroxypropyl-β-cyclodextrin (HP-β-CD) and MgCO₃. Solid dispersions (SD) have been prepared using a mechanochemical technique. The physical properties of nimesulide SD in solid state were characterized by differential scanning calorimetry and X-ray diffraction studies. The characteristics of the water solutions which form from the obtained solid dispersions were analyzed by reverse phase and gel permeation HPLC. It was shown that solubility increases for all complexes under investigation. These phenomena are obliged by complexation with auxiliary substances, which was shown by ¹H-NMR relaxation methods. The parallel artificial membrane permeability assay (PAMPA) was used for predicting passive intestinal absorption. Results showed that mechanochemically obtained complexes with polysaccharide AG, Na₂GA, and HP-β-CD enhanced permeation of NIM across an artificial membrane compared to that of the pure NIM. The complexes were examined for anti-inflammatory activity on a model of histamine edema. The substances were administered per os to CD-1 mice. As a result, it was found that all investigated complexes dose-dependently reduce the degree of inflammation. The best results were obtained for the complexes of NIM with Na₂GA and HP-β-CD. In noted case the inflammation can be diminished up to 2-fold at equal doses of NIM.     

新型抗菌肽的设计、活性研究及与磷脂相互作用的计算模拟

设计合成了多个具有2个活性序列的线性和环状多肽及具有单个活性序列的短链多肽, 研究了它们的杀菌活性, 发现其杀菌活性顺序为长链肽>环状肽>短链肽, 特别是线性的Linear-KT和Linear-KS对多种革兰氏阴性菌和阳性菌均具有较高的杀菌活性. 采用MTT法考察了Linear-KT和Linear-KS对正常细胞的毒性, 其中Linear-KS表现出较低的细胞毒性, 优于阳性对照多粘菌素B. 利用计算模拟的方法计算了多肽与细菌细胞膜中磷脂酰甘油(DMPG)的相互作用. 结果表明, 多肽和DMPG的结合能也表现出长链肽>环状肽>短链肽的规律, 特别是Linear-KT和Linear-KS具有较高的结合能. 长链肽含有2个活性序列, 可提供多个荷正电的氨基酸与荷负电的磷脂结合, 结合能较大, 杀菌活性较强. 同时, 柔性的结构及Linear-KT和Linear-KS中丝氨酸和苏氨酸的β碳上的羟基可与磷脂上的羰基形成多个氢键, 进一步增大了结合能. 计算模拟的方法为抗菌肽的杀菌活性从理论上提供了一定的依据.     

Hydration-induced changes of structure and vibrational frequencies of methylphosphocholine studied as a model of biomembrane lipids

The chemical characteristics of the polar parts of phospholipids as the main components of biological membranes were investigated by using infrared (IR) spectroscopy and theoretical calculations with water as a probe molecule. The logical key molecule used in this study is methylphosphocholine (MePC) as it is not only a representative model for a polar lipid headgroup but itself has biological significance. Isolated MePC forms a compact (folded) structure which is essentially stabilized by two intramolecular C-H...O type hydrogen bonds. At lower hydration, considerable wavenumber shifts were revealed by IR spectroscopy: the frequencies of the (O-P-O)- stretches were strongly redshifted, whereas methyl and methylene C-H and O-P-O stretches shifted surprisingly to blue. The origin of both red- and blueshifts was rationalized, on the basis of molecular-dynamics and quantum-chemistry calculations. In more detail, the hydration-induced blueshifts of C-H stretches could be shown to arise from several origins: disruption of the intramolecular C-H...O hydrogen bonds, formation of intermolecular C-H...O(water) H-bonds. The stepwise disruption of the intramolecular hydrogen bonds appeared to be the main feature that causes partial unfolding of the compact structure. However, the transition from a folded to extended MePC structure was completed only at high hydration. One might hypothesize that the mechanism of hydration-driven conformational changes as described here for MePC could be transferred to other zwitterions with relevant internal C-H...O hydrogen bonds.     

Folding propensity of cyclohexylether-delta-peptides

[structure: see text] Linear (n = 2-18) and cyclic oligomers (n = 3-8) of a cyclohexylether-delta-amino acid (COA) were prepared in high yield and stereopurity. CD spectra of the linear oligomers were indicative of secondary structure formation. X-ray crystal structures of cyclic COA oligomers revealed hydrophobic packing and internal 5- and 10-membered-ring hydrogen bonds. Ether and amide oxygens reside preferably in an ap orientation. This conformational locking is apparently broken by a C-2 substituent in an asymmetric cyclotrimer, for which a zeolithe-like tubular structure was found.     

Intramolecular co-operative hydrogen bond in calix[<Emphasis Type="Italic">n</Emphasis>]arenes (<Emphasis Type="Italic">n</Emphasis> = 4, 6, 8) bearing bulky substituents

Based on the Fourier transform IR spectroscopy together with the published NMR and X-ray data, it was shown that cyclic co-operative intramolecular hydrogen bond in calix[n]arene (n = 4, 6, 8) molecules is mainly responsible for their conformational state irrespective of the presence or absence of bulky substituents at the upper rim of the molecules. In accordance with the size of a macrocycle (n = 4, 6, 8), the stable conformation, secured by such a hydrogen bond, constitutes a cone, a pinched cone, and a pleated loop, respectively. The new, potentially competing system of hydrogen bonds in calix[6]arenes with 3-carboxymethyl-1-adamantyl substituents does not affect the conformational state of the macrocycle and its H-bonding. Six carboxy groups at the upper rim form in pairs three cyclic dimers, which does not disturb the hydrogen bonds of the hydroxy groups and the conformation of the macrocycle. In addition, the cavity of the molecule is considerably enlarged. The removal or rearrangement of the guest molecules in the solid calixarene by heating up to 180 °C only slightly affects the conformational state of macrocycles bearing bulky substituents, whereas in calixarenes devoid of such substituents, the similar procedure leads to somewhat of a distortion of the macrocycles (judging from the IR spectral indications of hydrogen bonding). Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1062–1068, June, 2007.     

侧链间氢键的协同效应对环状多肽自组装的影响

运用密度泛函(DFT)B3LYP方法和半经验分子轨道方法AM1对四种环状多肽[-(L-Asn-Ala)4-], [-(L-Asp-Ala)4-], [-(L-Gln-Ala)4-] 和[-(L-Glu-Ala)4-]的单体、平行和反平行二聚体到十聚体进行了理论研究. 结果表明, 四种环状多肽无论以平行还是以反平行的方式聚集, 聚集体中相邻两个环状多肽的侧链之间都能形成氢键. 侧链间氢键的相互作用使得这些环状多肽在组装过程中的结构和能量变化均表现出一定的协同效应, 这种协同效应加强了多肽纳米管的稳定性, 同时对聚集模式的选取起到了决定性作用.