Universal peptidomimetics

This paper concerns peptidomimetic scaffolds that can present side chains in conformations resembling those of amino acids in secondary structures without incurring excessive entropic or enthalpic penalties. Compounds of this type are referred to here as minimalist mimics. The core hypothesis of this paper is that small sets of such scaffolds can be designed to analogue local pairs of amino acids (including noncontiguous ones) in any secondary structure; i.e., they are universal peptidomimetics. To illustrate this concept, we designed a set of four peptidomimetic scaffolds. Libraries based on them were made bearing side chains corresponding to many of the protein-derived amino acids. Modeling experiments were performed to give an indication of kinetic and thermodynamic accessibilities of conformations that can mimic secondary structures. Together, peptidomimetics based on these four scaffolds can adopt conformations that resemble almost any combination of local amino acid side chains in any secondary structure. Universal peptidomimetics of this kind are likely to be most useful in the design of libraries for high-throughput screening against diverse targets. Consequently, data arising from submission of these molecules to the NIH Molecular Libraries Small Molecule Repository (MLSMR) are outlined.     

Generate: a program for 3-D structure generation and conformational analysis of peptides and peptidomimetics

The program Generate, aimed at generating 3-D structures for peptides and peptidomimetics, is presented. The algorithm is based on a build-up procedure, using a library of conformations of amino acid residues. This library is built from conformational analysis of amino acids placed in a di- or tripeptide environment to mimic the surroundings of the amino acid in a true peptide, considering different positions of the residue in the peptide chain (peptidyl fragment, NH(+)(3)-terminus or COO(-)-terminus). Cis-trans isomerism in the amide bonds is taken into account by construction of rotamer libraries for different isomers. Water solvation is included through the GB/SA model. New amino acid residues can easily be added to the libraries, making it possible to generate conformations of peptidomimetics.     

A Multifaceted Secondary Structure Mimic Based On Piperidine‐piperidinones

Minimalist secondary structure mimics are typically made to resemble one interface in a protein–protein interaction (PPI), and thus perturb it. We recently proposed suitable chemotypes can be matched with interface regions directly, without regard for secondary structures. Here we describe a modular synthesis of a new chemotype 1 , simulation of its solution‐state conformational ensemble, and correlation of that with ideal secondary structures and real interface regions in PPIs. Scaffold 1 presents amino acid side‐chains that are quite separated from each other, in orientations that closely resemble ideal sheet or helical structures, similar non‐ideal structures at PPI interfaces, and regions of other PPI interfaces where the mimic conformation does not resemble any secondary structure. 68 different PPIs where conformations of 1 matched well were identified. A new method is also presented to determine the relevance of a minimalist mimic crystal structure to its solution conformations. Thus dld ‐ 1 faf crystallized in a conformation that is estimated to be 0.91 kcal mol^?1 above the minimum energy solution state.     

Molecular modelling prediction of ligand binding site flexibility

We have investigated the efficacy of generating multiple sidechain conformations using a rotamer library in order to find the experimentally observed ligand binding site conformation of a protein in the presence of a bound ligand. We made use of a recently published algorithm that performs an exhaustive conformational search using a rotamer library to enumerate all possible sidechain conformations in a binding site. This approach was applied to a dataset of proteins whose structures were determined by X-ray and NMR methods. All chosen proteins had two or more structures, generally involving different bound ligands. By taking one of these structures as a reference, we were able in most cases to successfully reproduce the experimentally determined conformations of the other structures, as well as to suggest alternative low-energy conformations of the binding site. In those few cases where this procedure failed, we observed that the bound ligand had induced a high-energy conformation of the binding site. These results suggest that for most proteins that exhibit limited backbone motion, ligands tend to bind to low energy conformations of their binding sites. Our results also reveal that it is possible in most cases to use a rotamer search-based approach to predict alternative low-energy protein binding site conformations that can be used by different ligands. This opens the possibility of incorporating alternative binding site conformations to improve the efficacy of docking and structure-based drug design algorithms.     

Vibrational spectroscopic studies on fibrinogen adsorption at polystyrene/protein solution interfaces: hydrophobic side chain and secondary structure changes

Structural changes of fibrinogen after adsorption to polystyrene (PS) were examined at the PS/protein solution interface in situ using sum frequency generation (SFG) vibrational spectroscopy and attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR). Different behaviors of hydrophobic side chains and secondary structures of adsorbed fibrinogen molecules have been observed. Our results indicate that upon adsorption, the hydrophobic PS surface induces fast structural changes of fibrinogen molecules by aligning some hydrophobic side chains in fibrinogen so that they face to the surface. Such structural changes of fibrinogen hydrophobic side chains are local changes and do not immediately induce significant changes of the protein secondary structures. Our research also shows that the interactions between adsorbed fibrinogen and the PS surface can induce significant changes of protein secondary structures or global conformations which occur on a much longer time scale.     

NMR-assisted computational studies of peptidomimetic inhibitors bound in the hydrophobic pocket of HIV-1 glycoprotein 41

Due to the inherently flexible nature of a protein–protein interaction surface, it is difficult both to inhibit the association with a small molecule, and to predict how it might bind to the surface. In this study, we have examined small molecules that mediate the interaction between a WWI motif on the C-helix of HIV-1 glycoprotein-41 (gp41) and a deep hydrophobic pocket contained in the interior N-helical trimer. Association between these two components of gp41 leads to virus–cell and cell–cell fusion, which could be abrogated in the presence of an inhibitor that binds tightly in the pocket. We have studied a comprehensive combinatorial library of α-helical peptidomimetics, and found that compounds with strongly hydrophobic side chains had the highest affinity. Computational docking studies produced multiple possible binding modes due to the flexibility of both the binding site and the peptidomimetic compounds. We applied a transferred paramagnetic relaxation enhancement experiment to two selected members of the library, and showed that addition of a few experimental constraints enabled definitive identification of unique binding poses. Computational docking results were extremely sensitive to side chain conformations, and slight variations could preclude observation of the experimentally validated poses. Different receptor structures were required for docking simulations to sample the correct pose for the two compounds. The study demonstrated the sensitivity of predicted poses to receptor structure and indicated the importance of experimental verification when docking to a malleable protein–protein interaction surface.     

Electrophoretic dynamics of large DNA stars in polymer solutions and gels

We have developed a procedure for synthesizing large stable branched DNA structures that enables visualization via fluorescence microscopy. Using this procedure we have synthesized large DNA stars and observed their electrophoretic behavior in polymer solutions and gels. In dilute polyacrylamide solutions, the DNA stars move as random coils and appear to experience only brief collisions with the polymer chains in solution. The effect of polymer solution concentration on the electrophoretic mobility of stars in the dilute regime is found to be in good accord with predictions of the transient entanglement coupling (TEC) model. In semidilute polymer solutions, the star arms extend in the field direction and drag the core through the matrix. The star arms form several U-shaped conformations as they collide and engage with polyacrylamide chains. The U-shaped conformations occasionally evolve into J-shaped conformations as the star arms slide off the matrix chains they engage during electrophoretic migration. In concentrated polymer solutions, the arms of the star extend and form V-shaped structures with the core as the apex. The arms then pull the core through the matrix. These V-shaped conformations are much longer-lived than U-shaped ones and, unlike the latter, do not transform to J-shaped conformations. In polyacrylamide and agarose gels, where matrix entanglements are fixed, DNA stars become trapped when entanglements with matrix molecules prevent the core from being pulled through the matrix by the extended arms. This trapping was observed at all gel concentrations and electric fields studied.     

Importance of chirality and reduced flexibility of protein side chains: a study with square and tetrahedral lattice models

Side chains of amino acid residues are the determining factor that distinguishes proteins from other unstable chain polymers. In simple models they are often represented implicitly (e.g., by spin states) or simplified as one atom. Here we study side chain effects using two-dimensional square lattice and three-dimensional tetrahedral lattice models, with explicitly constructed side chains formed by two atoms of different chirality and flexibility. We distinguish effects due to chirality and effects due to side chain flexibilities, since residues in proteins are L residues, and their side chains adopt different rotameric states. For short chains, we enumerate exhaustively all possible conformations. For long chains, we sample effectively rare events such as compact conformations and obtain complete pictures of ensemble properties of conformations of these models at all compactness region. This is made possible by using sequential Monte Carlo techniques based on chain growth method. Our results show that both chirality and reduced side chain flexibility lower the folding entropy significantly for globally compact conformations, suggesting that they are important properties of residues to ensure fast folding and stable native structure. This corresponds well with our finding that natural amino acid residues have reduced effective flexibility, as evidenced by statistical analysis of rotamer libraries and side chain rotatable bonds. We further develop a method calculating the exact side chain entropy for a given backbone structure. We show that simple rotamer counting underestimates side chain entropy significantly for both extended and near maximally compact conformations. We find that side chain entropy does not always correlate well with main chain packing. With explicit side chains, extended backbones do not have the largest side chain entropy. Among compact backbones with maximum side chain entropy, helical structures emerge as the dominating configurations. Our results suggest that side chain entropy may be an important factor contributing to the formation of alpha helices for compact conformations.     

Metal complexes of chiral pentaazacrowns as conformational templates for beta-turn recognition

Examples of reverse turns as recognition motifs in biological systems can be found in high-resolution crystal structures of antibody-peptide complexes. Development of peptidomimetics is often based on replacing the amide backbone of peptides by sugar rings, steroids, benzodiazepines, or other hetero- and carbocycles. In this approach, the chemical scaffold of the peptide backbone can be replaced while retaining activity as long as the pharmacophoric groups of the peptide side chains stay in relatively the same place; in other words, similar functional groups must overlap in space for interaction with critical receptor sites. This study evaluates the potential of metal complexes of chiral pentaazacrowns (PAC) derived by reduction of cyclic pentapeptides as -turn mimetics. Due to the limited flexibility of the pendant chiral side groups in these metal complexes, one can potentially elicit information about the receptor-bound conformation from their binding affinities. 11 PAC crystal structures with different substitution patterns complexed with 3 different metals (Mn, Fe, Cd) as a prototypical database of potential side-chain orientations. Complexation with different metals induces subtle differences in the conformations of a particular azacrown scaffold. The lack of parameterization of transition metals for force field calculations precludes a thorough theoretical study. Thus, this study utilizes a simple geometrical comparison between the experimental data for crystalline PAC complexes and the side-chain orientations seen in classic -turns. The FOUNDATION program was used to overlap the C-C vectors of the corresponding ideal -turn side-chains to all possible leaving groups of the PAC complexes. When comparing the relative orientations of the chiral side chains, a strong overlap of the bonds (between about 0.1 Å to about 0.5 Å RMS for 3 residues and up to about 1 Å RMS for 4 residues) was observed for many of the molecules. Such metal complexes may lack complete peptidomimetic activity due to the lack of spatial overlap of all four side-chain residues, however, if only three peptide side chains are needed for receptor recognition and/or binding, the metal complexes should show biological activity.     

Thioether‐derived Macrocycle for Peptide Secondary Structure Fixation

Recently, we developed methods to stabilize peptides into various secondary structures, including α‐helix, type III turn and β‐hairpin via proper thioether based macrocyclization. These conformationally constrained peptidomimetics confer enhanced biophysical properties and provide a valuable avenue towards clinically‐relevant therapeutic molecules. In this personal account, thioether‐derived macrocyclization methods developed by our group for stabilization of α‐helix, type‐III β turn and β‐hairpin conformations are discussed.     

Understanding metalloprotein folding using a de novo design strategy

Metal ions play significant roles in most biological systems. Over the past two decades, there has been significant interest in the redesign of existing metal binding sites in proteins/peptides and the introduction of metals into folded proteins/peptides. Recent research has focused on the effects of metal binding on the overall secondary and tertiary conformations of unstructured peptides/proteins. In this context, de novo design of metallopeptides has become a valuable approach for studying the consequence of metal binding. It has been seen that metal ions not only direct folding of partially folded peptides but have at times also been the elixir for properly folding random-coil-like structures in stable secondary conformations. Work in our group has focused on binding of heavy metal ions such as Hg(II) to de novo designed alpha-helical three stranded coiled coil peptides with sequences based on the heptad repeat motif. Removal from or addition of a heptad to the parent 30-residue TRI peptide with the amino acid sequence Ac-G(LKALEEK)(4)G-NH(2) generated peptides whose self-aggregation affinities were seen to be dependent on their lengths. It was noted that adjustment in the position of the thiol from an "a" position in the case of the shorter BabyL9C to a "d" position for BabyL12C resulted in a peptide with low association affinities for itself, weaker binding with Hg(II), and a considerably faster kinetic profile for metal insertion. Similar differences in thermodynamic and kinetic parameters were also noted for the longer TRI peptides. At the same time, metal insertion into the prefolded and longer TRI and Grand peptides has clearly demonstrated that the metal binding is both thermodynamically as well kinetically different from that to unassociated peptides.     

A new method for ligand docking to flexible receptors by dual alanine scanning and refinement (SCARE)

Protein binding sites undergo ligand specific conformational changes upon ligand binding. However, most docking protocols rely on a fixed conformation of the receptor, or on the prior knowledge of multiple conformations representing the variation of the pocket, or on a known bounding box for the ligand. Here we described a general induced fit docking protocol that requires only one initial pocket conformation and identifies most of the correct ligand positions as the lowest score. We expanded a previously used diverse "cross-docking" benchmark to thirty ligand-protein pairs extracted from different crystal structures. The algorithm systematically scans pairs of neighbouring side chains, replaces them by alanines, and docks the ligand to each 'gapped' version of the pocket. All docked positions are scored, refined with original side chains and flexible backbone and re-scored. In the optimal version of the protocol pairs of residues were replaced by alanines and only one best scoring conformation was selected from each 'gapped' pocket for refinement. The optimal SCARE (SCan Alanines and REfine) protocol identifies a near native conformation (under 2 angstroms RMSD) as the lowest rank for 80% of pairs if the docking bounding box is defined by the predicted pocket envelope, and for as many as 90% of the pairs if the bounding box is derived from the known answer with approximately 5 angstroms margin as used in most previous publications. The presented fully automated algorithm takes about 2 h per pose of a single processor time, requires only one pocket structure and no prior knowledge about the binding site location. Furthermore, the results for conformationally conserved pockets do not deteriorate due to substantial increase of the pocket variability.     

Structural mining: self-consistent design on flexible protein-peptide docking and transferable binding affinity potential

A flexible protein-peptide docking method has been designed to consider not only ligand flexibility but also the flexibility of the protein. The method is based on a Monte Carlo annealing process. Simulations with a distance root-mean-square (dRMS) virtual energy function revealed that the flexibility of protein side chains was as important as ligand flexibility for successful protein-peptide docking. On the basis of mean field theory, a transferable potential was designed to evaluate distance-dependent protein-ligand interactions and atomic solvation energies. The potential parameters were developed using a self-consistent process based on only 10 known complex structures. The effectiveness of each intermediate potential was judged on the basis of a Z score, approximating the gap between the energy of the native complex and the average energy of a decoy set. The Z score was determined using experimentally determined native structures and decoys generated by docking with the intermediate potentials. Using 6600 generated decoys and the Z score optimization criterion proposed in this work, the developed potential yielded an acceptable correlation of R(2) = 0.77, with binding free energies determined for known MHC I complexes (Class I Major Histocompatibility protein HLA-A(*)0201) which were not present in the training set. Test docking on 25 complexes further revealed a significant correlation between energy and dRMS, important for identifying native-like conformations. The near-native structures always belonged to one of the conformational classes with lower predicted binding energy. The lowest energy docked conformations are generally associated with near-native conformations, less than 3.0 Angstrom dRMS (and in many cases less than 1.0 Angstrom) from the experimentally determined structures.     

Structure‐Based Design of Inhibitors of Protein–Protein Interactions: Mimicking Peptide Binding Epitopes

Protein–protein interactions (PPIs) are involved at all levels of cellular organization, thus making the development of PPI inhibitors extremely valuable. The identification of selective inhibitors is challenging because of the shallow and extended nature of PPI interfaces. Inhibitors can be obtained by mimicking peptide binding epitopes in their bioactive conformation. For this purpose, several strategies have been evolved to enable a projection of side chain functionalities in analogy to peptide secondary structures, thereby yielding molecules that are generally referred to as peptidomimetics. Herein, we introduce a new classification of peptidomimetics (classes A–D) that enables a clear assignment of available approaches. Based on this classification, the Review summarizes strategies that have been applied for the structure‐based design of PPI inhibitors through stabilizing or mimicking turns, β‐sheets, and helices.     

Effect of molecular conformations on the adsorption behavior of gold-binding peptides

Despite extensive recent reports on combinatorially selected inorganic-binding peptides and their bionanotechnological utility as synthesizers and molecular linkers, there is still only limited knowledge about the molecular mechanisms of peptide binding to solid surfaces. There is, therefore, much work that needs to be carried out in terms of both the fundamentals of solid-binding kinetics of peptides and the effects of peptide primary and secondary structures on their recognition and binding to solid materials. Here we discuss the effects of constraints imposed on FliTrx-selected gold-binding peptide molecular structures upon their quantitative gold-binding affinity. We first selected two novel gold-binding peptide (AuBP) sequences using a FliTrx random peptide display library. These were, then, synthesized in two different forms: cyclic (c), reproducing the original FliTrx gold-binding sequence as displayed on bacterial cells, and linear (l) dodecapeptide gold-binding sequences. All four gold-binding peptides were then analyzed for their adsorption behavior using surface plasmon resonance spectroscopy. The peptides exhibit a range of binding affinities to and adsorption kinetics on gold surfaces, with the equilibrium constant, Keq, varying from 2.5x10(6) to 13.5x10(6) M(-1). Both circular dichroism and molecular mechanics/energy minimization studies reveal that each of the four peptides has various degrees of random coil and polyproline type II molecular conformations in solution. We found that AuBP1 retained its molecular conformation in both the c- and l-forms, and this is reflected in having similar adsorption behavior. On the other hand, the c- and l-forms of AuBP2 have different molecular structures, leading to differences in their gold-binding affinities.     

Secondary and Tertiary Structure of Polysaccharides in Solutions and Gels

Many polysaccharide chains can adopt ordered helical and ribbon-like secondary structures. It seems however that these chains are often so stiff and extended that the cooperative interactions necessary for stability in the solvent environment can only be achieved when inter-chain as well as intra-chain interactions are favorable. Hence we commonly find two-or more-stranded associations of helices, of ribbons, or of helices with ribbons. These can be regarded as tertiary and higher levels of structure. The ordered secondary structure characteristically requires a regular repeating sequence of sugar residues, and the termination of this sequence by insertion of a residue of different type may also terminate the secondary structure and hence the association in which it is involved. This is the mechanism by which native polysaccharides may link up to form three dimensional networks, or gels, in which state they perform their natural roles in maintaining the hydration and integrity of biological tissues. For several polysaccharides there is evidence that the mechanism of biological control over the fine topology and properties of the gel network is mediated by enzymes which modify sugar residues at the polymer level to change the pattern of “interrupting” sugar residues.     

Why not consider a spherical protein? Implications of backbone hydrogen bonding for protein structure and function

The intrinsic ability of protein structures to exhibit the geometric features required for molecular function in the absence of evolution is examined in the context of three systems: the reference set of real, single domain protein structures, a library of computationally generated, compact homopolypeptides, artificial structures with protein-like secondary structural elements, and quasi-spherical random proteins packed at the same density as proteins but lacking backbone secondary structure and hydrogen bonding. Without any evolutionary selection, the library of artificial structures has similar backbone hydrogen bonding, global shape, surface to volume ratio and statistically significant structural matches to real protein global structures. Moreover, these artificial structures have native like ligand binding cavities, and a tiny subset has interfacial geometries consistent with native-like protein-protein interactions and DNA binding. In contrast, the quasi-spherical random proteins, being devoid of secondary structure, have a lower surface to volume ratio and lack ligand binding pockets and intermolecular interaction interfaces. Surprisingly, these quasi-spherical random proteins exhibit protein like distributions of virtual bond angles and almost all have a statistically significant structural match to real protein structures. This implies that it is local chain stiffness, even without backbone hydrogen bonding, and compactness that give rise to the likely completeness of the library solved single domain protein structures. These studies also suggest that the packing of secondary structural elements generates the requisite geometry for intermolecular binding. Thus, backbone hydrogen bonding plays an important role not only in protein structure but also in protein function. Such ability to bind biological molecules is an inherent feature of protein structure; if combined with appropriate protein sequences, it could provide the non-zero background probability for low-level function that evolution requires for selection to occur.     

Constraining Cyclic Peptides To Mimic Protein Structure Motifs

Many proteins exert their biological activities through small exposed surface regions called epitopes that are folded peptides of well‐defined three‐dimensional structures. Short synthetic peptide sequences corresponding to these bioactive protein surfaces do not form thermodynamically stable protein‐like structures in water. However, short peptides can be induced to fold into protein‐like bioactive conformations (strands, helices, turns) by cyclization, in conjunction with the use of other molecular constraints, that helps to fine‐tune three‐dimensional structure. Such constrained cyclic peptides can have protein‐like biological activities and potencies, enabling their uses as biological probes and leads to therapeutics, diagnostics and vaccines. This Review highlights examples of cyclic peptides that mimic three‐dimensional structures of strand, turn or helical segments of peptides and proteins, and identifies some additional restraints incorporated into natural product cyclic peptides and synthetic macrocyclic peptidomimetics that refine peptide structure and confer biological properties.     

Biaryl amino acid templates in place of D-Pro-L-Pro in cyclic beta-hairpin cationic antimicrobial peptidomimetics

The turn-forming D-Pro-L-Pro template has been frequently used to promote regular beta-hairpin conformations in cyclic protein epitope mimetics. Here the use of three isomeric biaryl templates has been studied as alternatives to D-Pro-L-Pro in the preparation of beta-hairpin peptidomimetics. The o,o'- o,m'- and m,m'-isomers of carboxymethyl- and aminomethyl-substituted biaryl templates have been incorporated into novel macrocyclic mimics of the naturally occurring cationic antimicrobial peptide protegrin I. The presence of the o-carboxymethyl-o'-aminomethyl-biaryl template within the macrocyclic peptide resulted in the appearance of slowly interconverting atropisomers. Although none of the resulting mimetics adopted stable beta-hairpin structures in aqueous solution, they all nevertheless retained a significant antimicrobial activity against Gram positive and Gram negative bacteria. These mimetics provide interesting starting points for an optimization program in the search for potent and novel antimicrobial compounds.     

Characterisation by immobilised metal ion affinity chromatographic procedures of the binding behaviour of several synthetic peptides designed to have high affinity for Cu(II) ions.

In this investigation, several peptides containing an increasing number of histidine residues have been designed and synthesised. The peptides involved repeat units of either the pentameric EAEHA or the tetrameric HLLH sequence motifs. Adsorption isotherms for these synthetic peptides and hexahistidine (hexa-His) as a control substance were measured under batch equilibrium binding conditions with an immobilised Cu(II)-iminodiacetic acid (IDA) sorbent. The experimental data were analysed in terms of Langmuirean binding behaviour. In common with previous studies with synthetic peptides, these investigations have demonstrate that the sequential organisation of the histidine side chains in these peptides can affect the selectivity of the coordination interactions with borderline metal ions in immobilised metal ion affinity chromatographic systems. The results also confirm that peptides selected on the basis of their potential to form amphipathic secondary structures with their histidine residues presented on one face of the molecule can exhibit equivalent or higher affinity constants towards copper ions than hexa-His, although they contain fewer histidine residues. These findings are thus relevant to the selection of peptides produced inter alia by combinatorial synthetic procedures to have enhanced binding properties for Cu(II) or Ni(II) ions, or intended for use as peptide tags in the fusion handle approach for the affinity chromatographic purification of recombinant proteins.