Kinetic Models of Cyclosporin A in Polar and Apolar Environments Reveal Multiple Congruent Conformational States

Design Principles for Balancing Lipophilicity and Permeability in beyond Rule of 5 Space

An ab initio conformational analysis of oral beyond Rule of 5 (bRo5) drugs was complemented with measured permeability and logP(octanol) to derive design principles conferring oral bioavailability. 3D polar surface area (PSA) thresholds for oral bRo5 drugs coincided with those reported for Ro5 space. The majority of oral bRo5 drugs exceeded the Ro5 logP threshold of 5, reflecting a bias for permeability. Above 500 Da molecular weight (MW), oral drugs and highly permeable Novartis compounds occupy a narrow polarity range (topological or TPSA/MW) of 0.1-0.3 Å2 /Da, whose upper half coincides with the lower 90 percentiles of the Novartis logP set. This TPSA/MW range and 3D PSA below 100 Å2 define the "Rule of ~1 /₅" for balancing lipophilicity and permeability. Neutral TPSA, defined as TPSA minus 3D PSA occurs independent of conformation, intramolecular hydrogen bonds (IMHB) and MW, suggesting it is an intrinsic molecular property. Neutral TPSA increased in the lead optimization (LO) campaigns of three first in class de novo designed bRo5 drugs and may be a useful design parameter in bRo5 space.

Lipid bilayer membrane permeability mechanism of the K-Ras(G12D)-inhibitory bicyclic peptide KS-58 elucidated by molecular dynamics simulations

Peptides are mid-size molecules (700-2000 g/mol) and have attracted particular interest as therapeutic modalities as they are superior in controlling protein-protein interactions, a process that is a typical drug target category, compared with small molecules (<500 g/mol). In 2020, we identified KS-58 (1333 g/mol) as a K-Ras(G12D)-inhibitory bicyclic peptide and suggested its cell membrane permeability. However, the membrane permeability mechanism had not been elucidated. In this study, we aim to clarify the mechanism by molecular dynamics (MD) simulations. Initially, we simulated the molecular conformations of KS-58 in water (a polar solvent) and in chloroform (a non-polar solvent). The identified stable conformations were significantly different in each solvent. KS-58 behaves as a chameleon-like molecule as it alters its polar surface area (PSA) depending on the solvent environment. It was also discovered that orientation of Asp's side chain is a critical energy barrier for KS-58 altering its conformation from hydrophilic to lipophilic. Taking these properties into consideration, we simulated its lipid bilayer membrane permeability. KS-58 shifted toward the inside of the lipid bilayer membrane with altering its conformations to lipophilic. When the simulation condition was set in deionized form of that carboxy group of Asp, KS-58 traveled deeper inside the cell membrane. PSA and the depth of the membrane penetration correlated. In vitro data suggested that cell membrane permeability of KS-58 is improved in weakly acidic conditions leading to partial deionization of the carboxy group. Our data provide an example of the molecular properties of mid-size peptides with membrane accessibility and propose an effective metadynamics approach to elucidate such molecular mechanisms by MD simulations.

Oral Absorption of Middle-to-Large Molecules and Its Improvement, with a Focus on New Modality Drugs

To meet unmet medical needs, middle-to-large molecules, including peptides and oligonucleotides, have emerged as new therapeutic modalities. Owing to their middle-to-large molecular sizes, middle-to-large molecules are not suitable for oral absorption, but there are high expectations around orally bioavailable macromolecular drugs, since oral administration is the most convenient dosing route. Therefore, extensive efforts have been made to create bioavailable middle-to-large molecules or develop absorption enhancement technology, from which some successes have recently been reported. For example, Rybelsus® tablets and Mycapssa® capsules, both of which contain absorption enhancers, were approved as oral medications for type 2 diabetes and acromegaly, respectively. The oral administration of Rybelsus and Mycapssa exposes their pharmacologically active peptides with molecular weights greater than 1000, namely, semaglutide and octreotide, respectively, into systemic circulation. Although these two medications represent major achievements in the development of orally absorbable peptide formulations, the oral bioavailability of peptides after taking Rybelsus and Mycapssa is still only around 1%. In this article, we review the approaches and recent advances of orally bioavailable middle-to-large molecules and discuss challenges for improving their oral absorption.

Conformations of Macrocyclic Peptides Sampled by Nuclear Magnetic Resonance: Models for Cell-Permeability

The biological activities and pharmacological properties of peptides and peptide mimetics are determined by their conformational states. Therefore, a detailed understanding of the conformational landscape is crucial for rational drug design. Nuclear magnetic resonance (NMR) is the only method for structure determination in solution. However, it remains challenging to determine the structures of peptides using NMR because of very weak nuclear Overhauser effects (NOEs), the semiquantitative nature of the rotating frame Overhauser effect (ROE), and the low number of NOEs/ROEs in N-methylated peptides. In this study, we introduce a new approach to investigating the structures of modified macrocyclic peptides. We utilize exact NOEs (eNOEs) in viscous solvent mixtures to replicate various cellular environments. eNOEs provide detailed structural information for highly dynamic modified peptides. Structures of high precision were obtained for cyclosporin A, with a backbone atom rmsd of 0.10 Å. Distinct conformational states in different environments were identified for omphalotin A (OmphA), a fungal nematotoxic and multiple backbone N-methylated macrocyclic peptides. A model for cell-permeation is presented for OmphA, based on its structures in polar, apolar, and mixed polarity solvents. During the transition from a polar to an apolar environment, OmphA undergoes a rearrangement of its H-bonding network, accompanied by a cis to trans isomerization of the ω torsion angle within a type VIa β-turn. We hypothesize that the kinetics of these conformational transitions play a crucial role in determining the membrane-permeation capabilities of OmphA.

Controlling the Chameleonic Behavior and Membrane Permeability of Cyclosporine Derivatives via Backbone and Side Chain Modifications

Some macrocycles exhibit enhanced membrane permeability through conformational switching in different environmental polarities, a trait known as chameleonic behavior. In this study, we demonstrate specific backbone and side chain modifications that can control chameleonic behavior and passive membrane permeability using a cyclosporin O (CsO) scaffold. To quantify chameleonic behavior, we used a ratio of the population of the closed conformation obtained in polar solvent and nonpolar solvent for each CsO derivative. We found that β-hydroxylation at position 1 (1 and 3) can encode chameleonicity and improve permeability. However, the conformational stabilization induced by adding an additional transannular H-bond (2 and 5) leads to a much slower rate of membrane permeation. Our CsO scaffold provides a platform for the systematic study of the relationship among conformation, membrane permeability, solubility, and protein binding. This knowledge contributes to the discovery of potent beyond the rule of five (bRo5) macrocycles capable of targeting undruggable targets.

Molecular Design of Cyclic Peptides with Cell Membrane Permeability and Development of MDMX-p53 Inhibitor

Cyclic peptides have been expected to be one of the modalities of intracellular protein-protein interaction (PPI) inhibitors, but they are generally known to have low cell membrane permeability. In this study, we focused on the conformation of cyclic peptides in the cell membrane to determine the requirement for their cell membrane permeability through passive diffusion. Utilizing the requirement, we searched for structures with high affinity for MDMX via computational chemistry and acquired cyclic peptide 19 (Papp = 0.80 × 10-6 cm s-1, IC50 = 0.07 μM).

Chamelogk: A Chromatographic Chameleonicity Quantifier to Design Orally Bioavailable Beyond-Rule-of-5 Drugs

New chemical modalities in drug discovery include molecules belonging to the bRo5 chemical space. Because of their complex and flexible structure, bRo5 compounds often suffer from a poor solubility/permeability profile. Chameleonicity describes the capacity of a molecule to adapt to the environment through conformational changes; the design of molecular chameleons is a medicinal chemistry strategy simultaneously optimizing solubility and permeability. A default method to quantify chameleonicity in early drug discovery is still missing. Here we introduce Chamelogk, an automated, fast, and cheap chromatographic descriptor of chameleonicity. Moreover, we report measurements for 55 Ro5 and bRo5 compounds and validate our method with literature data. Then, selected case studies (macrocycles, nonmacrocyclic compounds, and PROTACs) are used to illustrate the application of Chamelogk in combination with lipophilicity (BRlogD) and polarity (Δ log kwIAM) descriptors. Overall, we show how Chamelogk deserves being included in property-based drug discovery strategies to design oral bioavailable bRo5 compounds.

Cell-permeable chameleonic peptides: Exploiting conformational dynamics in de novo cyclic peptide design

Membrane-traversing peptides offer opportunities for targeting intracellular proteins and oral delivery. Despite progress in understanding the mechanisms underlying membrane traversal in natural cell-permeable peptides, there are still several challenges to designing membrane-traversing peptides with diverse shapes and sizes. Conformational flexibility appears to be a key determinant of membrane permeability of large macrocycles. We review recent developments in the design and validation of chameleonic cyclic peptides, which can switch between alternative conformations to enable improved permeability through cell membranes, while still maintaining reasonable solubility and exposed polar functional groups for target protein binding. Finally, we discuss the principles, strategies, and practical considerations for rational design, discovery, and validation of permeable chameleonic peptides.

Preparation and Surface Characterization of Chitosan-Based Coatings for PET Materials

Poly(ethylene terephthalate)-PET-is one of the most frequently used polymers in biomedical applications. Due to chemical inertness, PET surface modification is necessary to gain specific properties, making the polymer biocompatible. The aim of this paper is to characterize the multi-component films containing chitosan (Ch), phospholipid 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), immunosuppressant cyclosporine A (CsA) and/or antioxidant lauryl gallate (LG) which can be utilized as a very attractive material for developing the PET coatings. Chitosan was employed owing to its antibacterial activity and also its ability to promote cell adhesion and proliferation favorable for tissue engineering and regeneration purposes. Moreover, the Ch film can be additionally modified with other substances of biological importance (DOPC, CsA and LG). The layers of varying compositions were prepared using the Langmuir-Blodgett (LB) technique on the air plasma-activated PET support. Then their nanostructure, molecular distribution, surface chemistry and wettability were determined by atomic force microscopy (AFM), time-of-flight secondary ion mass spectrometry (TOF-SIMS), X-ray photoelectron spectroscopy (XPS), contact angle (CA) measurements and the surface free energy and its components' determination, respectively. The obtained results show clearly the dependence of the surface properties of the films on the molar ratio of components and allow for a better understanding of the coating organization and mechanisms of interactions at the molecular level both inside the films and between the films and the polar/apolar liquids imitating the environment of different properties. The organized layers of this type can be helpful in gaining control over the surface properties of the biomaterial, thus getting rid of the limitations in favor of increased biocompatibility. This is a good basis for further investigations on the correlation of the immune system response to the presence of biomaterial and its physicochemical properties.

Lessons for Oral Bioavailability: How Conformationally Flexible Cyclic Peptides Enter and Cross Lipid Membranes

Cyclic peptides extend the druggable target space due to their size, flexibility, and hydrogen-bonding capacity. However, these properties impact also their passive membrane permeability. As the "journey" through membranes cannot be monitored experimentally, little is known about the underlying process, which hinders rational design. Here, we use molecular simulations to uncover how cyclic peptides permeate a membrane. We show that side chains can act as "molecular anchors", establishing the first contact with the membrane and enabling insertion. Once inside, the peptides are positioned between headgroups and lipid tails─a unique polar/apolar interface. Only one of two distinct orientations at this interface allows for the formation of the permeable "closed" conformation. In the closed conformation, the peptide crosses to the lower leaflet via another "anchoring" and flipping mechanism. Our findings provide atomistic insights into the permeation process of flexible cyclic peptides and reveal design considerations for each step of the process.

Binding of Different Cyclosporin Variants to Micelles Evidenced by NMR and MD Simulations

Peptides play a critical role in the life of organisms, performing completely different functions. The biological activity of some peptides, such as cyclosporins, can be determined by the degree of membrane permeability. Thus, it becomes important to study how the molecule interacts with lipid bilayers. Cyclosporins C, E, H and L were characterised molecular dynamics simulation; NMR spectroscopy studies were also carried out for cyclosporins C and E. The comparison of one- and two-dimensional spectra revealed certain similarities between spatial structures of the studied cyclosporin variants. Upon dissolving in water containing DPC micelles, which serve as model membranes, subtle changes in the NMR spectra appear, but in a different way for different cyclosporins. In order to understand whether observed changes are related to any structural modifications, simulation of the interaction of the peptide with the phospholipid micelle was performed. The onset of the interaction was observed, when the peptide is trapped to the surface of the micelle. Simulations of this kind are also of interest in the light of the well-known membrane permeability of cyclosporin, which is important for its biological action.

Conformational Sampling Deciphers the Chameleonic Properties of a VHL-Based Degrader

Chameleonicity (the capacity of a molecule to adapt its conformations to the environment) may help to identify orally bioavailable drugs in the beyond-Rule-of-5 chemical space. Computational methods to predict the chameleonic behaviour of degraders have not yet been reported and the identification of molecular chameleons still relies on experimental evidence. Therefore, there is a need to tune predictions with experimental data. Here, we employ PROTAC-1 (a passively cell-permeable degrader), for which NMR and physicochemical data prove the chameleonic behaviour, to benchmark the capacity of two conformational sampling algorithms and selection schemes. To characterize the conformational ensembles in both polar and nonpolar environments, we compute three molecular properties proven to be essential for cell permeability: conformer shape (radius of gyration), polarity (3D PSA), and the number of intramolecular hydrogen bonds. Energetic criteria were also considered. Infographics monitored the simultaneous variation of those properties in computed and NMR conformers. Overall, we provide key points for tuning conformational sampling tools to reproduce PROTAC-1 chameleonicity according to NMR evidence. This study is expected to improve the design of PROTAC drugs and the development of computational sustainable strategies to exploit the potential of new modalities in drug discovery.

Simulation Reveals the Chameleonic Behavior of Macrocycles

Conformational analysis is central to the design of bioactive molecules. It is particularly challenging for macrocycles due to noncovalent transannular interactions, steric interactions, and ring strain that are often coupled. Herein, we simulated the conformations of five macrocycles designed to express a progression of increasing complexity in environment-dependent intramolecular interactions and verified the results against NMR measurements in chloroform and dimethyl sulfoxide. Molecular dynamics using an explicit solvent model, but not the Monte Carlo method with implicit solvation, handled both solvents correctly. Refinement of conformations at the ab initio level was fundamental to reproducing the experimental observations─standard state-of-the-art molecular mechanics force fields were insufficient. Our simulations correctly predicted the intramolecular interactions between side chains and the macrocycle and revealed an unprecedented solvent-induced conformational switch of the macrocyclic ring. Our results provide a platform for the rational, prospective design of molecular chameleons that adapt to the properties of the environment.

Ensemble Generation for Linear and Cyclic Peptides Using a Reservoir Replica Exchange Molecular Dynamics Implementation in GROMACS

The profile of shapes presented by a cyclic peptide modulates its therapeutic efficacy and is represented by the ensemble of its sampled conformations. Although some algorithms excel at creating a diverse ensemble of cyclic peptide conformations, they seldom address the entropic contribution of flexible conformations and often have significant practical difficulty producing an ensemble with converged and reliable thermodynamic properties. In this study, an accelerated molecular dynamics (MD) method, namely, reservoir replica exchange MD (R-REMD or Res-REMD), was implemented in GROMACS ver. 4.6.7 and benchmarked on two small cyclic peptide model systems: a cyclized furin cleavage site of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (cyclo-(CGPRRARSG)) and oxytocin (disulfide-bonded CYIQNCPLG). Additionally, we also benchmarked Res-REMD on alanine dipeptide and Trpzip2 to demonstrate its validity and efficiency over REMD. For Trpzip2, Res-REMD coupled with an umbrella-sampling-derived reservoir generated similar folded fractions as regular REMD but on a much faster time scale. For cyclic peptides, Res-REMD appeared to be marginally faster than REMD in ensemble generation. Finally, Res-REMD was more effective in sampling rare events such as trans to cis peptide bond isomerization. We provide a GitHub page with the modified GROMACS source code for running Res-REMD at https://github.com/PlotkinLab/Reservoir-REMD.

Surface Properties of the Polyethylene Terephthalate (PET) Substrate Modified with the Phospholipid-Polypeptide-Antioxidant Films: Design of Functional Biocoatings

Surface properties of polyethylene terephthalate (PET) coated with the ternary monolayers of the phospholipid 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), the immunosuppressant cyclosporine A (CsA), and the antioxidant lauryl gallate (LG) were examined. The films were deposited, by means of the Langmuir-Blodgett (LB) technique, on activated by air low temperature plasma PET plates (PET_air). Their topography and surface chemistry were determined with the help of atomic force microscopy (AFM) and time-of-flight secondary ion mass spectrometry (TOF-SIMS), respectively, while wettability was evaluated by the contact angle measurements. Then, the surface free energy and its components were calculated from the Lifshitz-van der Waals/Acid-Base (LWAB) approach. The AFM imaging showed that the Langmuir monolayers were transferred effectively and yielded smoothing of the PET_air surface. Mass spectrometry confirmed compatibility of the quantitative and qualitative compositions of the monolayers before and after the transfer onto the substrate. Moreover, the molecular arrangement in the LB films and possible mechanisms of DOPC-CsA-LG interactions were determined. The wettability studies provided information on the type and magnitude of the interactions that can occur between the biocoatings and the liquids imitating different environments. It was found that the changes from open to closed conformation of CsA molecules are driven by the hydrophobic environment ensured by the surrounding DOPC and LG molecules. This process is of significance to drug delivery where the CsA molecules can be released directly from the biomaterial surface by passive diffusion. The obtained results showed that the chosen techniques are complementary for the characterization of the molecular organization of multicomponent LB films at the polymer substrate as well as for designing biocompatible coatings with precisely defined wettability.

Development of the force field for cyclosporine A

Membrane permeability of cyclic peptides is an important factor in drug design. To investigate the membrane permeability of cyclic peptides using molecular dynamics (MD) simulations, the accurate force fields for unnatural amino acids present in the cyclic peptides are required. Therefore, we developed the CHARMM force fields of the unnatural amino acids present in cyclosporin A (CsA), a cyclic peptide used as an immune suppressor. Especially for N-methyl amino acids, which contribute to the membrane permeability of cyclic peptides, we developed a grid correction map (CMAP) of the energy surface using the φ and ψ dihedral angles in the main chain of CsA. To validate the developed force field, we performed MD simulations, including the generalized replica exchange with solute tempering method, of CsA in water and chloroform solvents. The conformations of CsA in water and chloroform sampled using the developed force field were consistent with those of the experimental results of the solution nuclear magnetic resonance spectroscopy.

Lipid Composition Is Critical for Accurate Membrane Permeability Prediction of Cyclic Peptides by Molecular Dynamics Simulations

Cyclic peptides have attracted attention as a promising pharmaceutical modality due to their potential to selectively inhibit previously undruggable targets, such as intracellular protein-protein interactions. Poor membrane permeability is the biggest bottleneck hindering successful drug discovery based on cyclic peptides. Therefore, the development of computational methods that can predict membrane permeability and support elucidation of the membrane permeation mechanism of drug candidate peptides is much sought after. In this study, we developed a protocol to simulate the behavior in membrane permeation steps and estimate the membrane permeability of large cyclic peptides with more than or equal to 10 residues. This protocol requires the use of a more realistic membrane model than a single-lipid phospholipid bilayer. To select a membrane model, we first analyzed the effect of cholesterol concentration in the model membrane on the potential of mean force and hydrogen bonding networks along the direction perpendicular to the membrane surface as predicted by molecular dynamics simulations using cyclosporine A. These results suggest that a membrane model with 40 or 50 mol % cholesterol was suitable for predicting the permeation process. Subsequently, two types of membrane models containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 40 and 50 mol % cholesterol were used. To validate the efficiency of our protocol, the membrane permeability of 18 ten-residue peptides was predicted. Correlation coefficients of R > 0.8 between the experimental and calculated permeability values were obtained with both model membranes. The results of this study demonstrate that the lipid membrane is not just a medium but also among the main factors determining the membrane permeability of molecules. The computational protocol proposed in this study and the findings obtained on the effect of membrane model composition will contribute to building a schematic view of the membrane permeation process. Furthermore, the results of this study will eventually aid the elucidation of design rules for peptide drugs with high membrane permeability.

Conformational Effects on the Passive Membrane Permeability of Synthetic Macrocycles

Macrocyclic compounds (MCs) can have complex conformational properties that affect pharmacologically important behaviors such as membrane permeability. We measured the passive permeability of 3600 diverse nonpeptidic MCs and used machine learning to analyze the results. Incorporating selected properties based on the three-dimensional (3D) conformation gave models that predicted permeability with Q² = 0.81. A biased spatial distribution of polar versus nonpolar regions was particularly important for good permeability, consistent with a mechanism in which the initial insertion of nonpolar portions of a MC helps facilitate the subsequent membrane entry of more polar parts. We also examined effects on permeability of 800 substructural elements by comparing matched molecular pairs. Some substitutions were invariably beneficial or invariably deleterious to permeability, while the influence of others was highly contextual. Overall, the work provides insights into how the permeability of MCs is influenced by their 3D conformational properties and suggests design hypotheses for achieving macrocycles with high membrane permeability.

Atomic View of Aqueous Cyclosporine A: Unpacking a Decades-Old Mystery

An atomic view of a main aqueous conformation of cyclosporine A (CycA), an important 11-amino-acid macrocyclic immunosuppressant, is reported. For decades, it has been a grand challenge to determine the conformation of free CycA in an aqueous-like solution given its poor water solubility. Using a combination of X-ray and single-crystal neutron diffraction, we unambiguously resolve a unique conformer (A1) with a novel cis-amide between residues 11 and 1 and two water ligands that stabilize hydrogen bond networks. NMR spectroscopy and titration experiments indicate that the novel conformer is as abundant as the closed conformer in 90/10 (v/v) methanol/water and is the main conformer at 10/90 methanol/water. Five other conformers were also detected in 90/10 methanol/water, one in slow exchange with A1, another one in slow exchange with the closed form and three minor ones, one of which contains two cis amides Abu2-Sar3 and MeBmt1-MeVal11. These conformers help better understand the wide spectrum of membrane permeability observed for CycA analogues and, to some extent, the binding of CycA to protein targets.

A New Amino Acid for Improving Permeability and Solubility in Macrocyclic Peptides through Side Chain-to-Backbone Hydrogen Bonding

Despite the notoriously poor membrane permeability of peptides, many cyclic peptide natural products show high passive membrane permeability and potently inhibit a variety of "undruggable" intracellular targets. A major impediment to the design of cyclic peptides with good permeability is the high desolvation energy associated with the peptide backbone amide NH groups. While several strategies have been proposed to mitigate this deleterious effect, only few studies have used polar side chains to sequester backbone NH groups. We investigated the ability of N,N-pyrrolidinylglutamine (Pye), whose side chain contains a powerful hydrogen-bond-accepting C═O amide group but no hydrogen-bond donors, to sequester exposed backbone NH groups in a series of cyclic hexapeptide diastereomers. Analyses revealed that specific Leu-to-Pye substitutions conferred dramatic improvements in aqueous solubility and permeability in a scaffold- and position-dependent manner. Therefore, this approach offers a complementary tool for improving membrane permeability and solubility in cyclic peptides.

Polar/apolar interfaces modulate the conformational behavior of cyclic peptides with impact on their passive membrane permeability

Cyclic peptides have the potential to vastly extend the scope of druggable proteins and lead to new therapeutics for currently untreatable diseases. However, cyclic peptides often suffer from poor bioavailability. To uncover design principles for permeable cyclic peptides, a promising strategy is to analyze the conformational dynamics of the peptides using molecular dynamics (MD) and Markov state models (MSMs). Previous MD studies have focused on the conformational dynamics in pure aqueous or apolar environments to rationalize membrane permeability. However, during the key steps of the permeation through the membrane, cyclic peptides are exposed to interfaces between polar and apolar regions. Recent studies revealed that these interfaces constitute the free energy minima of the permeation process. Thus, a deeper understanding of the behavior of cyclic peptides at polar/apolar interfaces is desired. Here, we investigate the conformational and kinetic behavior of cyclic decapeptides at a water/chloroform interface using unbiased MD simulations and MSMs. The distinct environments at the interface alter the conformational equilibrium as well as the interconversion kinetics of cyclic peptide conformations. For peptides with low population of the permeable conformation in aqueous solution, the polar/apolar interface facilitates the interconversion to the closed conformation, which is required for membrane permeation. Comparison to unbiased MD simulations with a POPC bilayer reveals that not only the conformations but also the orientations are relevant in a membrane system. These findings allow us to propose a permeability model that includes both ‘prefolding’ and ‘non-prefolding’ cyclic peptides – an extension that can lead to new design considerations for permeable cyclic peptides.

This study uses molecular dynamics and Markov state models to analyse how interfaces interact with cyclic decapeptides and modulate their dynamic and equilibrium properties.

Enhanced sampling without borders: on global biasing functions and how to reweight them

Molecular dynamics (MD) simulations are a powerful tool to follow the time evolution of biomolecular motions in atomistic resolution. However, the high computational demand of these simulations limits the timescales of motions that can be observed. To resolve this issue, so called enhanced sampling techniques are developed, which extend conventional MD algorithms to speed up the simulation process. Here, we focus on techniques that apply global biasing functions. We provide a broad overview of established enhanced sampling methods and promising new advances. As the ultimate goal is to retrieve unbiased information from biased ensembles, we also discuss benefits and limitations of common reweighting schemes. In addition to concisely summarizing critical assumptions and implications, we highlight the general application opportunities as well as uncertainties of global enhanced sampling.

Incorporating NOE-Derived Distances in Conformer Generation of Cyclic Peptides with Distance Geometry

Nuclear magnetic resonance (NMR) data from NOESY (nuclear Overhauser enhancement spectroscopy) and ROESY (rotating frame Overhauser enhancement spectroscopy) experiments can easily be combined with distance geometry (DG) based conformer generators by modifying the molecular distance bounds matrix. In this work, we extend the modern DG based conformer generator ETKDG, which has been shown to reproduce experimental crystal structures from small molecules to large macrocycles well, to include NOE-derived interproton distances. In noeETKDG, the experimentally derived interproton distances are incorporated into the distance bounds matrix as loose upper (or lower) bounds to generate large conformer sets. Various subselection techniques can subsequently be applied to yield a conformer bundle that best reproduces the NOE data. The approach is benchmarked using a set of 24 (mostly) cyclic peptides for which NOE-derived distances as well as reference solution structures obtained by other software are available. With respect to other packages currently available, the advantages of noeETKDG are its speed and that no prior force-field parametrization is required, which is especially useful for peptides with unnatural amino acids. The resulting conformer bundles can be further processed with the use of structural refinement techniques to improve the modeling of the intramolecular nonbonded interactions. The noeETKDG code is released as a fully open-source software package available at www.github.com/rinikerlab/customETKDG.

Cyclosporin A: Conformational Complexity and Chameleonicity

The chameleonic behavior of cyclosporin A (CsA) was investigated through conformational ensembles employing multicanonical molecular dynamics simulations that could sample the cis and trans isomers of N-methylated amino acids; these assessments were conducted in explicit water, dimethyl sulfoxide, acetonitrile, methanol, chloroform, cyclohexane (CHX), and n-hexane (HEX) using AMBER ff03, AMBER10:EHT, AMBER12:EHT, and AMBER14:EHT force fields. The conformational details were discussed employing the free-energy landscapes (FELs) at T = 300 K; it was observed that the experimentally determined structures of CsA were only a part of the conformational space. Comparing the ROESY measurements in CHX-d12 and HEX-d14, the major conformations in those apolar solvents were essentially the same as that in CDCl₃ except for the observation of some sidechain rotamers. The effects of the metal ions on the conformations, including the cis/trans isomerization, were also investigated. Based on the analysis of FELs, it was concluded that the AMBER ff03 force field best described the experimentally derived conformations, indicating that CsA intrinsically formed membrane-permeable conformations and that the metal ions might be the key to the cis/trans isomerization of N-methylated amino acids before binding a partner protein.

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.

Cyclosporin Structure and Permeability: From A to Z and Beyond

Cyclosporins are natural or synthetic undecapeptides with a wide range of actual and potential pharmaceutical applications. Several members of the cyclosporin compound family have remarkably high passive membrane permeabilities that are not well-described by simple structural metrics. Here we review experimental studies of cyclosporin structure and permeability, including cyclosporin-metal complexes. We also discuss models for the conformation-dependent permeability of cyclosporins and similar compounds. Finally, we identify current knowledge gaps in the literature and provide recommendations regarding future avenues of exploration.

Effect of Flexibility, Lipophilicity, and the Location of Polar Residues on the Passive Membrane Permeability of a Series of Cyclic Decapeptides

Cyclic peptides have received increasing attention over the recent years as potential therapeutics for "undruggable" targets. One major obstacle is, however, their often relatively poor bioavailability. Here, we investigate the structure-permeability relationship of 24 cyclic decapeptides that share the same backbone N-methylation pattern but differ in their side chains. The peptides cover a large range of values for passive membrane permeability as well as lipophilicity and solubility. To rationalize the observed differences in permeability, we extracted for each peptide the population of the membrane-permeable conformation in water from extensive explicit-solvent molecular dynamics simulations and used this as a metric for conformational rigidity or "prefolding." The insights from the simulations together with lipophilicity measurements highlight the intricate interplay between polarity/lipophilicity and flexibility/rigidity and the possible compensating effects on permeability. The findings allow us to better understand the structure-permeability relationship of cyclic peptides and extract general guiding principles.

Ensemble completeness in conformer sampling: the case of small macrocycles

In this study we compare the three algorithms for the generation of conformer ensembles Biovia BEST, Schrödinger Prime macrocycle sampling (PMM) and Conformator (CONF) form the University of Hamburg, with ensembles derived for exhaustive molecular dynamics simulations applied to a dataset of 7 small macrocycles in two charge states and three solvents. Ensemble completeness is a prerequisite to allow for the selection of relevant diverse conformers for many applications in computational chemistry. We apply conformation maps using principal component analysis based on ring torsions. Our major finding critical for all applications of conformer ensembles in any computational study is that maps derived from MD with explicit solvent are significantly distinct between macrocycles, charge states and solvents, whereas the maps for post-optimized conformers using implicit solvent models from all generator algorithms are very similar independent of the solvent. We apply three metrics for the quantification of the relative covered ensemble space, namely cluster overlap, variance statistics, and a novel metric, Mahalanobis distance, showing that post-optimized MD ensembles cover a significantly larger conformational space than the generator ensembles, with the ranking PMM > BEST >> CONF. Furthermore, we find that the distributions of 3D polar surface areas are very similar for all macrocycles independent of charge state and solvent, except for the smaller and more strained compound 7, and that there is also no obvious correlation between 3D PSA and intramolecular hydrogen bond count distributions.

Rational Control of Molecular Properties Is Mandatory to Exploit the Potential of PROTACs as Oral Drugs

To obtain new oral drugs in the beyond rule of five space, PROTACs among others, molecular properties should be optimized in early drug discovery. Degraders call for design strategies which focus on intramolecular interaction and chameleonicity. In parallel, tailored revalidation of permeability assessment and prediction methods becomes fundamental in this innovative chemical space.

Large-Scale Membrane Permeability Prediction of Cyclic Peptides Crossing a Lipid Bilayer Based on Enhanced Sampling Molecular Dynamics Simulations

Membrane permeability is a significant obstacle facing the development of cyclic peptide drugs. However, membrane permeation mechanisms are poorly understood. To investigate common features of permeable (and nonpermeable) designs, it is necessary to reproduce the membrane permeation process of cyclic peptides through the lipid bilayer. We simulated the membrane permeation process of 100 six-residue cyclic peptides across the lipid bilayer based on steered molecular dynamics (MD) and replica-exchange umbrella sampling simulations and predicted membrane permeability using the inhomogeneous solubility-diffusion model and a modified version of it. Furthermore, we confirmed the effectiveness of this protocol by predicting the membrane permeability of 56 eight-residue cyclic peptides with diverse chemical structures, including some confidential designs from a pharmaceutical company. As a result, a reasonable correlation between experimentally assessed and calculated membrane permeability of cyclic peptides was observed for the peptide libraries, except for strongly hydrophobic peptides. Our analysis of the MD trajectory demonstrated that most peptides were stabilized in the boundary region between bulk water and membrane and that for most peptides, the process of crossing the center of the membrane is the main obstacle to membrane permeation. The height of this barrier is well correlated with the electrostatic interaction between the peptide and the surrounding media. The structural and energetic features of the representative peptide at each vertical position within the membrane were also analyzed, revealing that peptides permeate the membrane by changing their orientation and conformation according to the surrounding environment.

Interplay among Conformation, Intramolecular Hydrogen Bonds, and Chameleonicity in the Membrane Permeability and Cyclophilin A Binding of Macrocyclic Peptide Cyclosporin O Derivatives

A macrocyclic peptide scaffold with well-established structure-property relationship is desirable for tackling undruggable targets. Here, we adopted a natural macrocycle, cyclosporin O (CsO) and its derivatives (CP1-3), and evaluated the impact of conformation on membrane permeability, cyclophilin A (CypA) binding, and the pharmacokinetic (PK) profile. In nonpolar media, CsO showed a similar conformation to cyclosporin A (CsA), a well-known chameleonic macrocycle, but less chameleonic behavior in a polar environment. The weak chameleonicity of CsO resulted in decreased membrane permeability; however, the more rigid conformation of CsO was not detrimental to its PK profile. CsO exhibited a higher plasma concentration than CsA, which resulted from minimal CypA binding and lower accumulation in red blood cells and moderate oral bioavailability (F = 12%). Our study aids understanding of CsO, a macrocyclic peptide that is less explored than CsA but with greater potential for diversity generation and rational design.

Improvement on Permeability of Cyclic Peptide/Peptidomimetic: Backbone N-Methylation as A Useful Tool

Peptides have a three-dimensional configuration that can adopt particular conformations for binding to proteins, which are well suited to interact with larger contact surface areas on target proteins. However, low cell permeability is a major challenge in the development of peptide-related drugs. In recent years, backbone N-methylation has been a useful tool for manipulating the permeability of cyclic peptides/peptidomimetics. Backbone N-methylation permits the adjustment of molecule’s conformational space. Several pathways are involved in the drug absorption pathway; the relative importance of each N-methylation to total permeation is likely to differ with intrinsic properties of cyclic peptide/peptidomimetic. Recent studies on the permeability of cyclic peptides/peptidomimetics using the backbone N-methylation strategy and synthetic methodologies will be presented in this review.

Permeability prediction in the beyond-Rule-of 5 chemical space: Focus on cyclic hexapeptides

Cyclic peptides (CPs) are gaining more and more relevance in drug discovery. Since one of their main drawbacks is poor permeability, the discovery of new orally available CP drugs requires computational tools that predict CP permeability in very early drug discovery. In this study we used a literature dataset of 62 cyclic hexapeptides to evaluate the performances of a number of in silico tools based on different computational theory to model and rationalize PAMPA and Caco-2 permeability values. In particular, we submitted the dataset to a) online calculators, b) QSPR strategies, c) a physics-based tool, d) a mixed approach and e) a kinetic method. This latter is an emergent strategy in which a few relevant conformations retrieved from a set of molecular dynamics (MD) simulations by the Markov State Model (MSM) are used to establish the compounds permeability. Both free and commercial software were used. Results were compared with a model based on experimental physicochemical descriptors. All the computational approaches but online calculators performed quite well and show that lipophilicity and not polarity is the main determinant of the investigated event. A second major outcome of the study is that the impact of flexibility on the permeability of the considered dataset cannot be unambiguously assessed. Finally, our comparative analysis, which also included not common applied strategies, allowed a sound evaluation of the pros and cons of the applied computational approaches.

Analysis of Molecular Interactions between Components in Phospholipid-Immunosuppressant-Antioxidant Mixed Langmuir Films

The study of Langmuir monolayers incorporating biomimetic and bioactive substances plays an important role today in assessing the properties and quality of the molecular films for potential biomedical applications. Here, miscibility of binary and ternary monolayers of phospholipid (dioleoyl phosphatidylcholine, DOPC), immunosuppressant (cyclosporine A, CsA), and antioxidant (lauryl gallate, LG) of varying molar fractions was analyzed by means of the Langmuir technique coupled with a surface potential (ΔV) module at the air-water interface. The surface pressure-area per molecule (π-A) isotherms provided information on the physical state of the films at a given surface pressure, the monolayer packing and ordering, and the type and strength of intermolecular interactions. Surface potential-area (ΔV-A) isotherms revealed the molecular orientation changes at the interface upon compression. In addition, the apparent dipole moment of the monolayer-forming molecules was determined from the surface potential isotherms. The obtained results indicated that the film compression provoked subsequent changes of CsA conformation and/or orientation, conferring better affinity for the hydrocarbon environment. The mutual interactions between the components were analyzed here in terms of the excess and total Gibbs energy of mixing, whose values depended on the stoichiometry of the mixed films. The strongest attraction, thus the highest thermodynamic stability, was found for a DOPC-CsA-LG mixture with a 1:1:2 molar ratio. Based on these results, a molecular model for the organization of the molecules within the Langmuir film was proposed. Through this model, we elucidated the significant role of LG in improving the miscibility of CsA in the model DOPC membrane and thus in increasing the stability of self-assembled monolayers by noncovalent interactions, such as H-bonds and Lifshitz-van der Waals forces. The above 1:1:2 combination of three components is revealed as the most promising film composition for the modification of implant device surfaces to improve their biocompatibility. Further insight into mechanisms concerning drug-membrane interactions at the molecular level is provided, which results in great importance for biocoating design and development as well as for drug release at target sites.

Preparative separation and purification of Cyclosporin D from fungus Hypoxylon Spp. by improved closed-loop recycling counter-current chromatography

Improved closed-loop recycling counter-current chromatography (CLR CCC) with a two-phase solvent system composed of n-hexane-acetonitrile (1:1, v/v) was developed for separation, purification and preparation of cyclosporin D from the crude extract of fungus Hypoxylon Spp. (sj18). 28 mg cyclosporin D was successfully purified from 300 mg crude extract sample. The purity was 95.2% after five cycles, determined by HPLC. The structure of cyclosporin D was identified and assigned by ¹H NMR, ¹³C NMR and mass spectrometric analyses. In addition, in the study, we show an interesting phenomenon that cyclosporin D can be prepared by the conventional CCC in n-hexane-ethyl acetate-methanol-water solvent system (2.5:1:2.5:1, v/v/v/v), and can also be prepared by the improved closed-loop recycling CCC in n-hexane-acetonitrile solvent system (1:1, v/v), but the efficiency of preparation varies greatly.

Flexibility and Cell Permeability of Cyclic Ras-Inhibitor Peptides Revealed by the Coupled Nosé-Hoover Equation

Quantifying the cell permeability of cyclic peptides is crucial for their rational drug design. However, the reasons remain unclear why a minor chemical modification, such as the difference between Ras inhibitors cyclorasin 9A5 and 9A54, can substantially change a peptide's permeability. To address this question, we performed enhanced sampling simulations of these two 11-mer peptides using the coupled Nosé-Hoover equation (cNH) we recently developed. The present cNH simulations realized temperature fluctuations over a wide range (240-600 K) in a dynamic manner, allowing structural samplings that were well validated by nuclear Overhauser effect measurements. The derived structural ensembles were comprehensively analyzed by all-atom structural clustering, mapping the derived clusters onto principal components (PCs) that characterize the cyclic structure, and calculating cluster-dependent geometric and chemical properties. The planar-open conformation was dominant in aqueous solvent, owing to inclusion of the Trp side chain in the main-chain ring, while the compact-closed conformation, which favors cell permeation due to its compactness and high polarity, was also accessible. Conformation-dependent cell permeability was observed in one of the derived PCs, demonstrating that decreased cell permeability in 9A54 is due to the high free energy barrier separating the two conformations. The origin of the change in free energy surface was determined to be loss of flexibility in the modified residues 2-3, resulting from the increased bulkiness of their side chains. The derived molecular mechanism of cell permeability highlights the significance of complete structural dynamics surveys for accelerating drug development with cyclic peptides.

Markov models from the square root approximation of the Fokker-Planck equation: calculating the grid-dependent flux

Molecular dynamics (MD) are extremely complex, yet understanding the slow components of their dynamics is essential to understanding their macroscopic properties. To achieve this, one models the MD as a stochastic process and analyses the dominant eigenfunctions of the associated Fokker-Planck operator, or of closely related transfer operators. So far, the calculation of the discretized operators requires extensive MD simulations. The square-root approximation of the Fokker-Planck equation is a method to calculate transition rates as a ratio of the Boltzmann densities of neighboring grid cells times a flux, and can in principle be calculated without a simulation. In a previous work we still used MD simulations to determine the flux. Here, we propose several methods to calculate the exact or approximate flux for various grid types, and thus estimate the rate matrix without a simulation. Using model potentials we test computational efficiency of the methods, and the accuracy with which they reproduce the dominant eigenfunctions and eigenvalues. For these model potentials, rate matrices with up to [Formula: see text] states can be obtained within seconds on a single high-performance compute server if regular grids are used.

Volume-scaled common nearest neighbor clustering algorithm with free-energy hierarchy

The combination of Markov state modeling (MSM) and molecular dynamics (MD) simulations has been shown in recent years to be a valuable approach to unravel the slow processes of molecular systems with increasing complexity. While the algorithms for intermediate steps in the MSM workflow such as featurization and dimensionality reduction have been specifically adapted to MD datasets, conventional clustering methods are generally applied to the discretization step. This work adds to recent efforts to develop specialized density-based clustering algorithms for the Boltzmann-weighted data from MD simulations. We introduce the volume-scaled common nearest neighbor (vs-CNN) clustering that is an adapted version of the common nearest neighbor (CNN) algorithm. A major advantage of the proposed algorithm is that the introduced density-based criterion directly links to a free-energy notion via Boltzmann inversion. Such a free-energy perspective allows a straightforward hierarchical scheme to identify conformational clusters at different levels of a generally rugged free-energy landscape of complex molecular systems.

Elucidating Solution Structures of Cyclic Peptides Using Molecular Dynamics Simulations

Protein-protein interactions are vital to biological processes, but the shape and size of their interfaces make them hard to target using small molecules. Cyclic peptides have shown promise as protein-protein interaction modulators, as they can bind protein surfaces with high affinity and specificity. Dozens of cyclic peptides are already FDA approved, and many more are in various stages of development as immunosuppressants, antibiotics, antivirals, or anticancer drugs. However, most cyclic peptide drugs so far have been natural products or derivatives thereof, with de novo design having proven challenging. A key obstacle is structural characterization: cyclic peptides frequently adopt multiple conformations in solution, which are difficult to resolve using techniques like NMR spectroscopy. The lack of solution structural information prevents a thorough understanding of cyclic peptides' sequence-structure-function relationship. Here we review recent development and application of molecular dynamics simulations with enhanced sampling to studying the solution structures of cyclic peptides. We describe novel computational methods capable of sampling cyclic peptides' conformational space and provide examples of computational studies that relate peptides' sequence and structure to biological activity. We demonstrate that molecular dynamics simulations have grown from an explanatory technique to a full-fledged tool for systematic studies at the forefront of cyclic peptide therapeutic design.

Structural Flexibility of Cyclosporine A Is Mediated by Amide Cis-Trans Isomerization and the Chameleonic Roles of Calcium

Falling outside of Lipinski's rule of five, macrocyclic drugs have accessed unique binding sites of their target receptors unreachable by traditional small molecules. Cyclosporin(e) A (CycA), an extensively studied macrocyclic natural product, is an immunosuppressant with undesirable side effects such as electrolytic imbalances. In this work, a comprehensive view on the conformational landscape of CycA, its interactions with Ca²⁺, and host-guest interactions with cyclophilin A (CypA) is reported through exhaustive analyses that combine ion-mobility spectrometry-mass spectrometry (IMS-MS), nuclear magnetic resonance (NMR) spectroscopy, distance-geometry modeling, and NMR-driven molecular dynamics. Our IMS-MS data show that CycA can adopt extremely compact conformations with significantly smaller collisional cross sections than the closed conformation observed in CDCl₃. To adopt these conformations, the macrocyclic ring has to twist and bend via cis-trans isomerization of backbone amides, and thus, we termed this family of structures the "bent" conformation. Furthermore, NMR measurements indicate that the closed conformation exists at 19% in CD₃OD/H₂O and 55% in CD₃CN. However, upon interacting with Ca²⁺, in addition to the bent and previously reported closed conformations of free CycA, the CycA:Ca²⁺ complex is open and has all-trans peptide bonds. Previous NMR studies using calcium perchlorate reported only the closed conformation of CycA (which contains one cis peptide bond). Here, calcium chloride, a more biologically relevant salt, was used, and interestingly, it helps converting the cis-MeLeu9-MeLeu10 peptide bond into a trans bond. Last, we were able to capture the native binding of CycA and CypA to give forth evidence that IMS-MS is able to probe the solution-phase structures of the complexes and that the Ca²⁺:CycA complex may play an essential role in the binding of CycA to CypA.

Computational design of mixed chirality peptide macrocycles with internal symmetry

Cyclic symmetry is frequent in protein and peptide homo‐oligomers, but extremely rare within a single chain, as it is not compatible with free N‐ and C‐termini. Here we describe the computational design of mixed‐chirality peptide macrocycles with rigid structures that feature internal cyclic symmetries or improper rotational symmetries inaccessible to natural proteins. Crystal structures of three C2‐ and C3‐symmetric macrocycles, and of six diverse S2‐symmetric macrocycles, match the computationally‐designed models with backbone heavy‐atom RMSD values of 1 Å or better. Crystal structures of an S4‐symmetric macrocycle (consisting of a sequence and structure segment mirrored at each of three successive repeats) designed to bind zinc reveal a large‐scale zinc‐driven conformational change from an S4‐symmetric apo‐state to a nearly inverted S4‐symmetric holo‐state almost identical to the design model. These symmetric structures provide promising starting points for applications ranging from design of cyclic peptide based metal organic frameworks to creation of high affinity binders of symmetric protein homo‐oligomers. More generally, this work demonstrates the power of computational design for exploring symmetries and structures not found in nature, and for creating synthetic switchable systems.

PDB Code(s): 6UFU, 6UG2, 6UG3, 6UG6, 6UGB, 6UGC, 6UCX, 6UD9, 6UDR, 6UDW, 6UDZ, 6UF4, 6UF7, 6UF8, 6UFA and 6UF9;

Depsipeptide synthesis using a late-stage Ag(i)-promoted macrolactonisation of peptide thioamides

Ag(i)-Triggered activation of the thioamide and attack by the C-terminal carboxylate generates an isoimide intermediate that undergoes an intramolecular acyl transfer to furnish the cyclic depsipeptide.

The importance of intramolecular hydrogen bonds on the translocation of the small drug piracetam through a lipid bilayer

The number of hydrogen bond donors and acceptors is a fundamental molecular descriptor to predict the oral bioavailability of small drug candidates. In fact, the most widely used oral bioavailability rules (such as the Lipinsky's rule-of-five and the Veber rules) make use of this molecular descriptor. It is generally assumed that hydrogen bond donors and acceptors impact on passive diffusion across cell membranes, a fundamental event during drug absorption and distribution. Although the relationship between the number of these motifs and the probability of having good oral bioavailability has been studied and described for more than 20 years, little attention has been given to their spatial distribution in the molecule. In this paper, we used molecular dynamics to describe the effect of intramolecular hydrogen bonding on the passive diffusion of a small drug (piracetam) through a lipid membrane. The results indicated that the formation of an intramolecular hydrogen bond decreases the barrier for translocation by ca. 4 kcal mol-1 and increases the permeability of the tested molecule, partially compensating the desolvation penalty arising from the penetration of the drug into the biological membrane core. This effect was apparent in simulations where the formation of this interaction was prevented with the help of modified potentials, and in simulations with a similar compound to piracetam that was not able to form this intramolecular hydrogen bond due to a larger distance between the hydrogen bond donor and acceptor groups. These results were also supported by coarse-grained methods, which are becoming an important resource for sampling a larger chemical space of molecules, with reduced computational effort. Furthermore, entropy and enthalpy derived profiles were also obtained as the compounds translocated across the membrane, suggesting that, even though the process of formation of internal hydrogen bonds is entropically unfavorable, the enthalpic gain is such that the formation of these interactions is beneficial for the passive diffusion across cell membranes.

Drug-Like Properties in Macrocycles above MW 1000: Backbone Rigidity versus Side-Chain Lipophilicity

Large macrocyclic peptides can achieve surprisingly high membrane permeability, although the properties that govern permeability in this chemical space are only beginning to come into focus. We generated two libraries of cyclic decapeptides with stable cross-β conformations, and found that peptoid substitutions within the β-turns of the macrocycle preserved the rigidity of the parent scaffold, whereas peptoid substitutions in the opposing β-strands led to "chameleonic" species that were rigid in nonpolar media but highly flexible in water. Both rigid and chameleonic compounds showed high permeability over a wide lipophilicity range, with peak permeabilities differing significantly depending on scaffold rigidity. Our findings indicate that modulating lipophilicity can be used to engineer favorable ADME properties into both rigid and flexible macrocyclic peptides, and that scaffold rigidity can be used to tune optimal lipophilicity.

The role of NMR in leveraging dynamics and entropy in drug design

Nuclear magnetic resonance (NMR) spectroscopy has contributed to structure-based drug development (SBDD) in a unique way compared to the other biophysical methods. The potency of a ligand binding to a protein is dictated by the binding free energy, which is an intricate interplay between entropy and enthalpy. In addition to providing the atomic resolution structural information, NMR can help to identify protein-ligand interactions that potentially contribute to the enthalpic component of the free energy. NMR can also illuminate dynamic aspects of the interaction, which correspond to the entropic term of the free energy. The ability of NMR to access both terms in the free energy equation stems from the suite of experiments developed to shed light on various aspects that contribute to both entropy and enthalpy, deepening our understanding of the biological function of macromolecules and assisting to target them in physiological conditions. Here we provide a brief account of the contribution of NMR to SBDD, highlighting hallmark examples and discussing the challenges that demand further method development. In the era of integrated biology, the unique ability of NMR to directly ascertain structural and dynamical aspects of macromolecule and monitor changes in these properties upon engaging a ligand can be combined with computational and other structural and biophysical methods to provide a more complete picture of the energetics of drug engagement with the target. Such efforts can be used to engineer better drugs.

Automated Design of Macrocycles for Therapeutic Applications: From Small Molecules to Peptides and Proteins

Macrocycles and cyclic peptides are increasingly attractive therapeutic modalities as they often have improved affinity, are able to bind to extended protein surfaces, and otherwise have favorable properties. Macrocyclization of a known binder may stabilize its bioactive conformation and improve its metabolic stability, cell permeability, and in certain cases oral bioavailability. Herein, we present implementation and application of an approach that automatically generates, evaluates, and proposes cyclizations utilizing a library of well-established chemical reactions and reagents. Using the three-dimensional (3D) conformation of the linear molecule in complex with a target protein as the starting point, this approach identifies attachment points, generates linkers, evaluates their geometric compatibility, and ranks the resulting molecules with respect to their predicted conformational stability and interactions with the target protein. As we show here with prospective and retrospective case studies, this procedure can be applied for the macrocyclization of small molecules and peptides and even PROteolysis TArgeting Chimeras (PROTACs) and proteins.

Connecting the conformational behavior of cyclic octadepsipeptides with their ionophoric property and membrane permeability

Cyclic octadepsipeptides such as PF1022A and its synthetic derivative emodepside exhibit anthelmintic activity with the latter sold as a commercial drug treatment against gastrointestinal nematodes for animal health use. The structure-permeability relationship of these cyclic depsipeptides that could ultimately provide insights into the compound bioavailability is not yet well understood. The fully N-methylated amide backbone and apolar sidechain residues do not allow for the formation of intramolecular hydrogen bonds, normally observed in the membrane-permeable conformations of cyclic peptides. Hence, any understanding gained on these depsipeptides would serve as a prototype for future design strategies. In previous nuclear magnetic resonance (NMR) studies, two macrocyclic core conformers of emodepside were detected, one with all backbone amides in trans-configuration (hereon referred as the symmetric conformer) and the other with one amide in cis-configuration (hereon referred as the asymmetric conformer). In addition, these depsipeptides were also reported to be ionophores with a preference of potassium over sodium. In this study, we relate the conformational behavior of PF1022A, emodepside, and closely related analogs with their ionophoric characteristic probed using NMR and molecular dynamics (MD) simulations and finally evaluated their passive membrane permeability using PAMPA. We find that the equilibrium between the two core conformers shifts more towards the symmetric conformer upon addition of monovalent cations with selectivity for potassium over sodium. Both the NMR experiments and the theoretical Markov state models based on extensive MD simulations indicate a more rigid backbone for the asymmetric conformation, whereas the symmetric conformation shows greater flexibility. The experimental results further advocate for the symmetric conformation binding the cation. The PAMPA results suggest that the investigated depsipeptides are retained in the membrane, which may be advantageous for the likely target, a membrane-bound potassium channel.