Amide-to-ester substitution as a stable alternative to N-methylation for increasing membrane permeability in cyclic peptides.

Naturally occurring peptides with high membrane permeability often have ester bonds on their backbones. However, the impact of amide-to-ester substitutions on the membrane permeability of peptides has not been directly evaluated. Here we report the effect of amide-to-ester substitutions on the membrane permeability and conformational ensemble of cyclic peptides related to membrane permeation. Amide-to-ester substitutions are shown to improve the membrane permeability of dipeptides and a model cyclic hexapeptide. NMR-based conformational analysis and enhanced sampling molecular dynamics simulations suggest that the conformational transition of the cyclic hexapeptide upon membrane permeation is differently influenced by an amide-to-ester substitution and an amide N-methylation. The effect of amide-to-ester substitution on membrane permeability of other cyclic hexapeptides, cyclic octapeptides, and a cyclic nonapeptide is also investigated to examine the scope of the substitution. Appropriate utilization of amide-to-ester substitution based on our results will facilitate the development of membrane-permeable peptides.

CycPeptMPDB: A Comprehensive Database of Membrane Permeability of Cyclic Peptides.

Recently, cyclic peptides have been considered breakthrough drugs because they can interact with "undruggable" targets such as intracellular protein-protein interactions. Membrane permeability is an essential indicator of oral bioavailability and intracellular targeting, and the development of membrane-permeable peptides is a bottleneck in cyclic peptide drug discovery. Although many experimental data on membrane permeability of cyclic peptides have been reported, a comprehensive database is not yet available. A comprehensive membrane permeability database is essential for developing computational methods for cyclic peptide drug design. In this study, we constructed CycPeptMPDB, the first web-accessible database of cyclic peptide membrane permeability. We collected information on a total of 7334 cyclic peptides, including the structure and experimentally measured membrane permeability, from 45 published papers and 2 patents from pharmaceutical companies. To unambiguously represent cyclic peptides larger than small molecules, we used the hierarchical editing language for macromolecules notation to generate a uniform sequence representation of peptides. In addition to data storage, CycPeptMPDB provides several supporting functions such as online data visualization, data analysis, and downloading. CycPeptMPDB is expected to be a valuable platform to support membrane permeability research on cyclic peptides. CycPeptMPDB can be freely accessed at http://cycpeptmpdb.com.

Realization of the structural fluctuation of biomolecules in solution: Generalized Langevin mode analysis.

A new theoretical method, referred to as Generalized Langevin Mode Analysis (GLMA), is proposed to analyze the mode of structural fluctuations of a biomolecule in solution. The method combines the two theories in the statistical mechanics, or the Generalized Langevin theory and the RISM/3D-RISM theory, to calculate the second derivative, or the Hessian matrix, of the free energy surface of a biomolecule in aqueous solution, which consists of the intramolecular interaction among atoms in the biomolecule and the solvation free energy. The method is applied to calculate the wave-number spectrum of an alanine dipeptide in water for which the optical heterodyne-detected Raman-induced spectroscopy (RIKES) spectrum is available to compare with. The theoretical analysis reproduced the main features of the experimental spectrum with respect to the peak positions of the four bands around ~90 cm-1 , ~240 cm-1 , ~370 cm-1 , and 400 cm-1 , observed in the experimental spectrum, in spite that the physics involved in the two spectrum was not exactly the same: the experimental spectrum includes the contributions from the dipeptide and the water molecules interacting with the solute, while the theoretical one is just concerned with the solute molecule, influenced by solvation. Two major discrepancies between the theoretical and experimental spectra, one in the band intensity around ~100 cm-1 , and the other in the peak positions around ~370 cm-1 , are discussed in terms of the fluctuation mode of water molecules interacting with the dipeptide, which is not taken explicitly into account in the theoretical analysis.

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.

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.

Differential mode of cholesterol inclusion with 2-hydroxypropyl-cyclodextrins increases safety margin in treatment of Niemann-Pick disease type C.

BACKGROUND AND PURPOSE: Niemann-Pick disease type C (NPC) is a lysosomal storage disorder with disrupted intracellular cholesterol trafficking. A cyclic heptasaccharide, 2-hydroxypropyl-ß-cyclodextrin (HP-ß-CD), is a cholesterol solubilizer that is being developed to treat NPC, but its ototoxicity and pulmonary toxicity remain important issues. We have characterized 2-hydroxypropyl-γ-cyclodextrin (HP-γ-CD), a cyclic octasaccharide with a larger cavity than HP-ß-CD, as a candidate drug to treat NPC. However, the molecular target of HP-γ-CD with respect to NPC and its potential for clinical application are still unclear. EXPERIMENTAL APPROACH: We investigated the mode of interaction between HP-γ-CD and cholesterol by phase-solubility analysis, proton NMR spectroscopy and molecular dynamics simulations. We then evaluated the therapeutic effects of HP-γ-CD compared with HP-ß-CD using cellular and murine NPC models. Mouse auditory and pulmonary function tests were also conducted. KEY RESULTS: HP-γ-CD solely formed a 1:1 inclusion complex with cholesterol with an affinity similar to that of HP-ß-CD. In vitro, HP-γ-CD and HP-ß-CD amelioration of NPC-related manifestations was almost equivalent at lower concentrations. However, at higher concentrations, the cholesterol inclusion mode of HP-ß-CD shifted to the highly soluble 2:1 complex whereas that of HP-γ-CD maintained solely the 1:1 complex. The constant lower cholesterol solubilizing ability of HP-γ-CD conferred it with significantly reduced toxicity compared with HP-ß-CD, but equal efficacy in treating a mouse model of NPC. CONCLUSIONS AND IMPLICATIONS: HP-γ-CD can serve as a fine-tuned cholesterol solubilizer for the treatment of NPC with a wider safety margin than HP-ß-CD in terms of ototoxicity and pulmonary toxicity.

Computational Screening of a Functional Cyclodextrin Derivative for Suppressing a Side Effect of Doxorubicin.

The binding affinity of the beta-cyclodextrin (ß-CyD) derivatives with Doxorubicin (Dox) is evaluated by means of the 3D-RISM/KH theory combined with the molecular dynamics simulation in order to screen the compounds for suppressing a side-effect of the cancer drug. A protocol revised for the external and conformational entropies of the host-guest system is employed to calculate the binding free energy. It is found that the direct interactions of CyD with Dox and the desolvation free-energies of the both compounds largely cancel out to leave moderate contributions to the affinity, which are comparable to those from the entropies. The results shed light on the entropy terms for determining the binding affinity, although the external-entropy terms are essentially constant over all the compounds examined and do not affect the screening. The theoretical result is compared with the experimental data of the association constant for a CyD derivative which was predicted to be the best compound by the preliminary calculation without the entropy terms.

Molecular Recognition and Self-Organization in Life Phenomena Studied by a Statistical Mechanics of Molecular Liquids, the RISM/3D-RISM Theory.

There are two molecular processes that are essential for living bodies to maintain their life: the molecular recognition, and the self-organization or self-assembly. Binding of a substrate by an enzyme is an example of the molecular recognition, while the protein folding is a good example of the self-organization process. The two processes are further governed by the other two physicochemical processes: solvation and the structural fluctuation. In the present article, the studies concerning the two molecular processes carried out by Hirata and his coworkers, based on the statistical mechanics of molecular liquids or the RISM/3D-RISM theory, are reviewed.

New Protocol for Predicting the Ligand-Binding Site and Mode Based on the 3D-RISM/KH Theory.

An efficient algorithm to find the binding position and mode of small ligands bound at an active site of protein is proposed based on the spatial distribution function (SDF) obtained from the three-dimensional reference interaction site model (3D-RISM) theory with the Kovalenko-Hirata (KH) closure relation. The ligand examined includes hydrophobic, acidic, and basic molecules and zwitterions. Eighteen different types of proteins, which serve as targets for those ligands, are selected to examine the robustness of the algorithm. An imaginary atom, referred to as an "anchor site", is defined at the center of geometry of a ligand molecule that serves as a center for searching the binding position and mode of the ligand molecule in the translational and rotational spaces. The probable binding sites (PBSs) are identified based on the SDFs of the ligand molecules around the protein, and the PBS is ranked according to the peak height of SDF. The deviations from the mean height of the peak values of SDFs for 50 PBSs are analyzed based on the z-score, which is a measure of prominence of the site. The PBS found at the closest distance from the anchor site of the crystal structure is referred to as the "nearest site". The orientation of the ligand molecule at each PBS is explored by changing the Euler angles, and the most probable binding mode is determined based on the superposition approximation. The binding position of ligand molecules is successfully predicted as one of the distinct peaks in SDF of the anchor site, with a few exceptions. The binding mode of the ligand molecule predicted based on the superposition approximation is consistent with the X-ray crystal structure in nine systems, a half of the systems investigated. The significance of the results is discussed in detail. An application of the new protocol to fragment-based drug discovery is suggested.

Predicting the Binding Mode of 2-Hydroxypropyl-ß-cyclodextrin to Cholesterol by Means of the MD Simulation and the 3D-RISM-KH Theory.

It has been found that a cyclodextrin derivative, 2-hydroxypropyl-ß-cyclodextrin (HPßCD), has reasonable therapeutic effect on Niemann-Pick disease type C, which is caused by abnormal accumulation of unesterified cholesterol and glycolipids in the lysosomes and shortage of esterified cholesterol in other cellular compartments. We study the binding affinity and mode of HPßCD with cholesterol to elucidate the possible mechanism of HPßCD for removing cholesterol from the lysosomes. The dominant binding mode of HPßCD with cholesterol is found based on the molecular dynamics simulation and a statistical mechanics theory of liquids, or the three-dimensional reference interaction site model theory with Kovalenko-Hirata closure relation. We examine the two types of complexes between HPßCD and cholesterol, namely, one-to-one (1:1) and two-to-one (2:1). It is predicted that the 1:1 complex makes two or three types of stable binding mode in solution, in which the ßCD ring tends to be located at the edge of the steroid skeleton. For the 2:1 complex, there are four different types of the complex conceivable, depending on the orientation between the two HPßCDs: head-to-head (HH), head-to-tail (HT), tail-to-head (TH), and tail-to-tail (TT). The HT and HH cyclodextrin dimers show higher affinity to cholesterol compared to the other dimers and to all the binding modes of 1:1 complexes. The physical reason why the HT and HH dimers have higher affinity compared to the other complexes is discussed based on the consistency with the 1:1 complex. On the one hand, in case of the HT and HH dimers, the position of each CD in the dimer along the cholesterol chain comes right on or close to one of the positions where a single CD makes a stable complex. On the other hand, one of the CD molecules is located on unstable region along the cholesterol chain, for the case of TH and TT dimers.

Perspective: Structural fluctuation of protein and Anfinsen's thermodynamic hypothesis.

The thermodynamics hypothesis, casually referred to as "Anfinsen's dogma," is described theoretically in terms of a concept of the structural fluctuation of protein or the first moment (average structure) and the second moment (variance and covariance) of the structural distribution. The new theoretical concept views the unfolding and refolding processes of protein as a shift of the structural distribution induced by a thermodynamic perturbation, with the variance-covariance matrix varying. Based on the theoretical concept, a method to characterize the mechanism of folding (or unfolding) is proposed. The transition state, if any, between two stable states is interpreted as a gap in the distribution, which is created due to an extensive reorganization of hydrogen bonds among back-bone atoms of protein and with water molecules in the course of conformational change. Further perspective to applying the theory to the computer-aided drug design, and to the material science, is briefly discussed.

A Systematic Analysis of the Binding Affinity between the Pim-1 Kinase and Its Inhibitors Based on the MM/3D-RISM/KH Method.

A systematic study of the binding affinities of 16 lead compounds targeting the Pim-1 kinase based on the 3D-RISM/KH theory and MD simulations is reported. The results show a correlation coefficient R = 0.69 between the theoretical and experimental values of the binding free energy. This demonstrates that the method is applicable to the problem of compound screening and lead optimization, for which relative values of the free energy among the compounds have significance. We elucidate the contribution of the ligand fragments to the binding free energy. Our results indicate that the interactions between the residues and the triazolo[4,3-b]pyridazine scaffold as well as the phenyl ring of the ligand molecule make significant contributions to stabilization of the complex. Using the 3D-RISM/KH theory, we further analyze the probability distribution of a ligand fragment around the protein-ligand complex in which the substituent around the phenyl ring is removed from the ligand. The results demonstrate that the 3D-RISM/KH theory is capable of predicting the position of substitution on a ligand that has a higher affinity to a target protein.

Predicting the binding free energy of the inclusion process of 2-hydroxypropyl-ß-cyclodextrin and small molecules by means of the MM/3D-RISM method.

A protocol to calculate the binding free energy of a host-guest system is proposed based on the MM/3D-RISM method, taking cyclodextrin derivatives and their ligands as model systems. The protocol involves the procedure to identify the most probable binding mode (MPBM) of receptors and ligands by means of the umbrella sampling method. The binding free energies calculated by the MM/3D-RISM method for the complexes of the seven ligands with the MPBM of the cyclodextrin, and with the fluctuated structures around it, are in agreement with the corresponding experimental data in a semi-quantitative manner. It suggests that the protocol proposed here is promising for predicting the binding affinity of a small ligand to a relatively rigid receptor such as cyclodextrin.

Topological and sequence information predict that foldons organize a partially overlapped and hierarchical structure.

It has been suggested that proteins have substructures, called foldons, which can cooperatively fold into the native structure. However, several prior investigations define foldons in various ways, citing different foldon characteristics, thereby making the concept of a foldon ambiguous. In this study, we perform a Go model simulation and analyze the characteristics of substructures that cooperatively fold into the native-like structure. Although some results do not agree well with the experimental evidence due to the simplicity of our coarse-grained model, our results strongly suggest that cooperatively folding units sometimes organize a partially overlapped and hierarchical structure. This view makes us easy to interpret some different proposal about the foldon as a difference of the hierarchical structure. On the basis of this finding, we present a new method to assign foldons and their hierarchy, using structural and sequence information. The results show that the foldons assigned by our method correspond to the intermediate structures identified by some experimental techniques. The new method makes it easy to predict whether a protein folds sequentially into the native structure or whether some foldons fold into the native structure in parallel.

Implication of the cause of differences in 3D structures of proteins with high sequence identity based on analyses of amino acid sequences and 3D structures.

BACKGROUND: Proteins that share a high sequence homology while exhibiting drastically different 3D structures are investigated in this study. Recently, artificial proteins related to the sequences of the GA and IgG binding GB domains of human serum albumin have been designed. These artificial proteins, referred to as GA and GB, share 98% amino acid sequence identity but exhibit different 3D structures, namely, a 3α bundle versus a 4ß + α structure. Discriminating between their 3D structures based on their amino acid sequences is a very difficult problem. In the present work, in addition to using bioinformatics techniques, an analysis based on inter-residue average distance statistics is used to address this problem. RESULTS: It was hard to distinguish which structure a given sequence would take only with the results of ordinary analyses like BLAST and conservation analyses. However, in addition to these analyses, with the analysis based on the inter-residue average distance statistics and our sequence tendency analysis, we could infer which part would play an important role in its structural formation. CONCLUSIONS: The results suggest possible determinants of the different 3D structures for sequences with high sequence identity. The possibility of discriminating between the 3D structures based on the given sequences is also discussed.

Analyses of the folding properties of ferredoxin-like fold proteins by means of a coarse-grained Go model: relationship between the free energy profiles and folding cores.

The folding mechanisms of proteins with multi-state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on-going research issues. In this study, we investigate ferredoxin-like fold proteins which have a simple topology and multi-state transitions. We analyze the folding processes by means of a coarse-grained Go model. We are able to reproduce the differences in the folding mechanisms between U1A, which has a high-free-energy intermediate state, and ADA2h and S6, which fold into the native structure through two-state transitions. The folding pathways of U1A, ADA2h, S6, and the S6 circular permutant, S6_p54-55, are reproduced and compared with experimental observations. We show that the ferredoxin-like fold contains two common regions consisting folding cores as predicted in other studies and that U1A produces an intermediate state due to the distinct cooperative folding of each core. However, because one of the cores of S6 loses its cooperativity and the two cores of ADA2h are tightly coupled, these proteins fold into the native structure through a two-state mechanism.

Incorporating into a Cα Go model the effects of geometrical restriction on Cα atoms caused by side chain orientations.

Coarse-grained Go models have been widely used for studying protein-folding mechanisms. Despite the simplicity of the model, these can reproduce the essential features of the folding process of a protein. However, it is also known that side chains significantly contribute to the folding mechanism. Hence, it is desirable to incorporate the side chain effects into a coarse-grained Go model. In this study, to distinguish the effects of side chain orientation and to understand how these effects contribute to folding mechanisms, we incorporate into a Cα Go model not only heterogeneous contact energies but also geometrical restraints around two Cα atoms in contact with each other. We confirm that the heterogeneity of contact energies governs the folding pathway of a protein and that the geometric constraints attributed to side chains reproduce cooperative transitions in folding.