The energy landscape of a selective tumor-homing pentapeptide

Recently, a potentially powerful strategy based on phage-display libraries has been presented to target tumors via homing peptides attached to nanoparticles. The Cys-Arg-Glu-Lys-Ala (CREKA) peptide sequence has been identified as a tumor-homing peptide that binds to clotted plasmas proteins present in tumor vessels and interstitium. The aim of this work consists of mapping the conformational profile of CREKA to identify the bioactive conformation. For this purpose, a conformational search procedure based on modified simulated annealing combined with molecular dynamics was applied to three systems that mimic the experimentally used conditions: (i) the free peptide; (ii) the peptide attached to a nanoparticle; and (iii) the peptide inserted in a phage display protein. In addition, the free peptide was simulated in an ionized aqueous solution environment, which mimics the ionic strength of the physiological medium. Accessible minima of all simulated systems reveal a multiple interaction pattern involving the ionized side chains of Arg, Glu, and Lys, which induces a beta-turn motif in the backbone observed in all simulated CREKA systems.     

Engineering strategy to improve peptide analogs: from structure-based computational design to tumor homing

We present a chemical strategy to engineer analogs of the tumor-homing peptide CREKA (Cys-Arg-Glu-Lys-Ala), which binds to fibrin and fibrin-associated clotted plasma proteins in tumor vessels (Simberg et al. in Proc Natl Acad Sci USA 104:932–936, 2007) with improved ability to inhibit tumor growth. Computer modeling using a combination of simulated annealing and molecular dynamics were carried out to design targeted replacements aimed at enhancing the stability of the bioactive conformation of CREKA. Because this conformation presents a pocket-like shape with the charged groups of Arg, Glu and Lys pointing outward, non-proteinogenic amino acids α-methyl and N-methyl derivatives of Arg, Glu and Lys were selected, rationally designed and incorporated into CREKA analogs. The stabilization of the bioactive conformation predicted by the modeling for the different CREKA analogs matched the tumor fluorescence results, with tumor accumulation increasing with stabilization. Here we report the modeling, synthetic procedures, and new biological assays used to test the efficacy and utility of the analogs. Combined, our results show how studies based on multi-disciplinary collaboration can converge and lead to useful biomedical advances.     

A simulation strategy for the atomistic modeling of flexible molecules covalently tethered to rigid surfaces: Application to peptides

A computational strategy to model flexible molecules tethered to a rigid inert surface is presented. The strategy is able to provide uncorrelated relaxed microstructures at the atomistic level. It combines an algorithm to generate molecules tethered to the surface without atomic overlaps, a method to insert solvent molecules and ions in the simulation box, and a powerful relaxation procedure. The reliability of the strategy has been investigated by simulating two different systems: (i) mixed monolayers consisting of binary mixtures of long‐chain alkyl thiols of different lengths adsorbed on a rigid inert surface and (ii) CREKA (Cys‐Arg‐Glu‐Lys‐Ala), a short linear pentapeptide that recognizes clotted plasma proteins and selectively homes to tumors, covalently tethered to a rigid inert surface in aqueous solution. In the first, we examined the segregation of the two species in the monolayers using different long‐chain:short‐chain ratios, whereas in the second, we explored the conformational space of CREKA and ions distribution considering densities of peptides per nm² ranging from 0.03 to 1.67. Results indicate a spontaneous segregation in alkyl thiol monolayers, which enhances when the concentration of longest chains increases. However, the whole conformational profile of CREKA depends on the number of molecules tethered to the surface pointing out the large influence of molecular density on the intermolecular interactions, even though the bioactive conformation was found as the most stable in all cases. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Conformational profile of a proline-arginine hybrid

The intrinsic conformational preferences of a new nonproteinogenic amino acid have been explored by computational methods. This tailored molecule, named ((β)Pro)Arg, is conceived as a replacement for arginine in bioactive peptides when the stabilization of folded turn-like conformations is required. The new residue features a proline skeleton that bears the guanidilated side chain of arginine at the C(β) position of the five-membered pyrrolidine ring, in either a cis or a trans orientation with respect to the carboxylic acid. The conformational profiles of the N-acetyl-N'-methylamide derivatives of the cis and trans isomers of ((β)Pro)Arg have been examined in the gas phase and in solution by B3LYP/6-31+G(d,p) calculations and molecular dynamics simulations. The main conformational features of both isomers represent a balance between geometric restrictions imposed by the five-membered pyrrolidine ring and the ability of the guanidilated side chain to interact with the backbone through hydrogen bonds. Thus, both cis- and trans-((β)Pro)Arg exhibit a preference for the α(L) conformation as a consequence of the interactions established between the guanidinium moiety and the main-chain amide groups.     

A proposed bioactive conformation of Peptide T

The conformational profiles of Peptide T, (5–8)Peptide T, [Abu⁵](4–8)Peptide T and (4–8)Peptide T were computed independently to assess the geometrical characteristics of the bioactive conformation of Peptide T. The conformational profiles of the peptides were computed within the molecular mechanics framework using an effective dielectric constant of 80. The conformational space was thoroughly sampled using an iterative simulated annealing protocol. The bioactive conformation was assessed by pairwise cross comparisons of each of the unique low energy conformations found for each of the different analogs studied. After a putative bioactive conformation was selected, in order to further validate our hypothesis the conformational profile of the potent compound cyclo(Thr-Thr-Asn-Tyr-Thr-Asp) was computed and the putative bioactive conformation was found. The conformation exhibits a pseudo -turn involving the side chain of Thr⁵ and the carbonyl oxygen of Tyr⁷ forming a C12 ring.     

The conformation of cyclo(-D-Pro-Ala4-) as a model for cyclic pentapeptides of the DL4 type

The conformation of the cyclic pentapeptide cyclo(-D-Pro-Ala(4)-) in solution and in the solid state was reinvestigated using modern NMR techniques. To allow unequivocal characterization of hydrogen bonds, relaxation behavior, and intramolecular distances, differently labeled isotopomers were synthesized. The NMR results, supported by extensive MD simulations, demonstrate unambiguously that the preferred conformation previously described by us, but recently questioned, is indeed correct. The validation of the conformational preferences of this cyclic peptide is important given that this system is a template for several bioactive compounds and for controlled "spatial screening" for the search of bioactive conformations.     

Fibrin Association at Hybrid Biointerfaces Made of Clot‐Binding Peptides and Polythiophene

The properties as biointerfaces of electroactive conducting polymer–peptide biocomposites formed by poly(3,4‐ethylenedioxythiophene) (PEDOT) and CREKA or CR(NMe)EKA peptide sequences (where Glu has been replaced by N‐methyl‐Glu in the latter) have been compared. CREKA is a linear pentapeptide that recognizes clotted plasma proteins and selectively homes to tumors, while CR(NMe)EKA is an engineer to improve such properties by altering peptide–fibrin interactions. Differences between PEDOT‐CREKA and PEDOT‐CR(NMe)EKA reflect dissemblance in the organization of the peptides into the polymeric matrix. Both peptides affect fibrinogen thrombin‐catalyzed polymerization causing the immediate formation of fibrin, whereas in the absence of thrombin this phenomenon is only observed for CR(NMe)EKA. Consistently, the fibrin‐adsorption capacity is higher for PEDOT‐CR(NMe)EKA than for PEDOT‐CREKA, even though in both cases adsorbed fibrin exhibits round‐like morphologies rather than the characteristic fibrous structure. PEDOT‐peptide films coated with fibrin are selective in terms of cell adhesion, promoting the attachment of metastatic cells with respect to normal cells.

     

Adsorption of homopolypeptides on gold investigated using atomistic molecular dynamics

We investigate the role of dynamics on adsorption of peptides to gold surfaces using all-atom molecular dynamics simulations in explicit solvent. We choose six homopolypeptides [Ala(10), Ser(10), Thr(10), Arg(10), Lys(10), and Gln(10)], for which experimental surface coverages are not correlated with amino acid level affinities for gold, with the idea that dynamic properties may also play a role. To assess dynamics we determine both conformational movement and flexibility of the peptide within a given conformation. Low conformational movement indicates stability of a given conformation and leads to less adsorption than homopolypeptides with faster conformational movement. Likewise, low flexibility within a given conformation also leads to less adsorption. Neither amino acid affinities nor dynamic considerations alone predict surface coverage; rather both quantities must be considered in peptide adsorption to gold surfaces.     

Conformational preference and potential templates of N-methylated cyclic pentaalanine peptides

Systematic N-methylation of all peptide bonds in the cyclic pentapeptide cyclo(-D-Ala-Ala(4)-) has been performed yielding 30 different N-methylated derivatives, of which only seven displayed a single conformation on the NMR time scale. The conformation of these differentially N-methylated peptides was recently reported by us (J. Am. Chem. Soc. 2006, 128, 15 164-15 172). Here we present the conformational characterization of nine additional N-methylated peptides from the previous library which are not homogeneous but exist as a mixture in which at least one conformation is preferred by over 80 %. The structures of these peptides are investigated employing various 2D-NMR techniques, distance geometry calculations and further refined by molecular dynamics simulations in explicit DMSO. The comparison of the conformation of these nine peptides and the seven conformationally homogeneous peptides allow us to draw conclusions regarding the influence of N-methylation on the peptide backbone of cyclic pentapeptide of the class cyclo(-D-Ala-Ala(4)-). Here we present the different conformational classes of the peptides arising from the definitive pattern of N-methylation which can eventually serve as templates for the design of bioactive peptides.     

Theoretical conformational analysis of opiate peptides Leu-Enkephalin (H-Tyr-Gly-Gly-Phe-Leu-OH) and its two thioamide analogs (H-Tyr-Glyψ[CSNH]Gly-Phe-Leu-OH) and (H-Tyr-Gly-Glyψ[CSNH]Phe-Leu-OH)

In order to find informations on the native structure of the Leu-Enkephalin opiate peptide, the parent peptide and its two thioamide analogs (Thio-Gly2)Leu-Enkephalin and (Thio-Gly3)Leu-Enkephalin were studied by the theoretical method PEPSEA. This comparative conformational analysis showed that the active conformation is a β turn structure centered on Gly3 and Phe4. Moreover, this study showed also that the more active analog (Thio-Gly2)Leu-Enk has a lower tendency to adopt this structure. Consequently, its high activity can only be explained by its long lifetime due to its resistance to enzymatic hydrolysis, following the substitution of the amide linkage by the thioamide one. The weakly active analog (Thio-Gly3)Leu-Enk does not adopt this structure and prefers instead a β turn structure centered on Gly2 and Gly3. This study also confirmed the importance of the distances between the Tyr and Phe residues at positions 1 and 4, and that of the terminal Tyrosine N-H group which must be free of any intramolecular hydrogen bond in order to be available in the molecular recognition process.     

Impact of azaproline on Peptide conformation

The amino acid analog azaproline (azPro) contains a nitrogen atom in place of the C(alpha) of proline. Peptides containing azPro were shown to stabilize the cis-amide conformer for the acyl-azPro bond and prefer type VI beta-turns both in crystals and in organic solvents by NMR. The increased stability for cis-amide conformers was relatively minor with respect to the trans-conformers. Further, their conformational preferences were depended on solvent. To elucidate the impact of azPro substitution on amide cis-trans isomerism and peptide conformation, this paper reports ab initio studies on azPro derivatives and a comparison with their cognate Pro derivatives: 1-acetyl-2-methyl pyrrolidine (1), 1-acetyl-2-methyl pyrazolidine (2), Ac-Pro-NHMe (3), Ac-azPro-NHMe (4), Ac-azPro-NMe(2) (5), Ac-azAzc-NHMe (6), and Ac-azPip-NHMe (7). Conformational preferences were explored at the MP2/6-31+G** level of theory in vacuo. Solvation effects for 1 and 2 were studied implicitly using the polarizable continuum model and explicitly represented by interactions with a single water molecule. An increase in the conformational preference for the cis-amide conformer of azPro was clearly seen. An intramolecular hydrogen bond occurred solely in the trans-amide conformer that reduced the preference for the cis-conformer by 2.2 kcal/mol. The larger ring homolog aza-pipecolic acid (azPip), in which this internal hydrogen bond was diminished, significantly augmented stabilization of the cis-amide conformer. In aqueous solution, the preference for the cis-amide conformers was greatly reduced, mainly as a result of interaction between water and the lone pair of the alpha-nitrogen in the trans-amide conformer that was 3.8 kcal/mol greater than that in the cis-conformer. In the azPro analog, the energy barrier for cis-trans amide isomerization was 6 kcal/mol less than that in the cognate Pro derivative. Because the azPro derivatives can stabilize the cis-amide bond and mimic a type VI beta-turn without incorporation of additional steric bulk, such a simple chemical modification of the peptide backbone provides a useful conformational constraint when incorporated into the structure of selected bioactive peptides. Such modifications can scan receptors for biological recognition of reverse turns containing cis-amide bonds by the incorporation of type VI beta-turn scaffolds with oriented appended side chains.     

Synthesis of a Sialyl Lewis x Mimic with Fixed Carboxylic Acid Group: Chemical Approach toward the Elucidation of the Bioactive Conformation of Sialyl Lewis x

The relation of the solution and bioactive conformation of sialyl Lewis x (sLe(x)) has been addressed by chemical means. To mimic the preferred solution conformation of sLe(x) 1, the more rigid analog 2 has been designed and synthesized. The sialic acid residue of 1 was replaced by a carboxylic acid function which is fixed in the equatorial position of a six membered ring acetal fused to galactose. Due to entropic considerations, an increased biological activity could be expected if the preferred solution conformation and bound form of sLe(x) were similar. Since mimic 2 was found to be inactive in an E-selectin binding assay, the bound form of sLe(x) most probably differs from the prevailing solution conformation.     

Mimicking the Biologically Active Part of the Cyclopeptides Segetalin A and B by “Clipping” of a Linear Tripeptide Derivative by Metal Coordination

A new method for the conformational fixation of bioactive loop-type peptide structures is presented. Hereby, ligand moieties are attached to the termini of a linear peptide sequence. Upon metal complexation, a macrocyclic structure with a loop-type conformation of the peptide is formed. As a representative example, the preparation of a WAG-bridged dicatechol derivative is described which mimics the active part of the natural products Segetalin A and Segetalin B.     

Interaction of a tripeptide with cesium perfluorooctanoate micelles

The interaction of alanyl-phenylalanyl-alanine (Ala-Phe-Ala) with the micelles formed by cesium perfluorooctanoate (CsPFO) in water was studied in the isotropic phase by means of 1H NMR and by molecular dynamics (MD) simulations. Information on the location of the peptide was experimentally obtained from selective variations in Ala-Phe-Ala chemical shifts and from differential line broadening in the presence of the paramagnetic ion Mn2+. The peptide-micelle association constant was estimated analyzing the chemical shift variations of the most sensitive Ala-Phe-Ala resonances with the peptide concentration. MD simulations of Ala-Phe-Ala in the micellar environment confirmed the experimental observations, identifying the hydrogen bonding interactions of the different peptide moieties with the micelle, yielding a binding constant close to the experimental one. NOESY experiments suggest that the peptide in the micellar environment does not adopt a preferred conformation but is mainly unstructured. Details on the conformational behavior of the peptide in the micellar solution observed through MD were consistent with a different conformational equilibrium in the proximity of the micelle. Information on Ala-Phe-Ala dynamics was obtained from 1H T1 data and compared to MD simulation results on the overall tumbling motion.     

Evidence of secondary structure by high-resolution magic angle spinning NMR spectroscopy of a bioactive peptide bound to different solid supports.

The structure of the 19-amino acid peptide epitope, corresponding to the 141-159 sequence of capsid viral protein VP1 of foot-and-mouth disease virus (FMDV), bound to three different resins, namely, polystyrene-MBHA, PEGA, and POEPOP, has been determined by high-resolution magic angle spinning (HRMAS) NMR spectroscopy. A combination of homonuclear and heteronuclear bidimensional experiments was used for the complete peptide resonance assignment and the qualitative characterization of the peptide folding. The influence of the chemicophysical nature of the different polymers on the secondary structure of the covalently attached FMDV peptide was studied in detail. In the case of polystyrene-MBHA and polyacrylamide-PEGA resins, the analysis of the 2D spectra was hampered by missing signals and extensive overlaps, and only a propensity toward a peptide secondary structure could be derived from the assigned NOE correlations. When the FMDV peptide was linked to the polyoxyethylene-based POEPOP resin, it was found to adopt in dimethylformamide a helical conformation encompassing the C-terminal domain from residues 152 to 159. This conformation is very close to that of the free peptide previously analyzed in 2,2,2-trifluoroethanol. Our study clearly demonstrates that a regular helical structure can be adopted by a resin-bound bioactive peptide. Moreover, a change in the folding was observed when the same peptide-POEPOP conjugate was swollen in aqueous solution, displaying the same conformational features as the free peptide in water. The possibility of studying solid-supported ordered secondary structures by the HRMAS NMR technique in a wide range of solvents can be extended either to other biologically relevant peptides and proteins or to new synthetic oligomers.     

Influence of hydrogen bonds on charge distribution and conformation of L-arginine

The molecular recognition of adenosine-5'-triphosphate (ATP) with L-arginine (Arg) through hydrogen bonding interactions has been found using 1H NMR, H-H NOESY, acidity titration and fluorescence spectra techniques. The interactions could influence charge distribution in Arg and induce Arg conformational variation. It is realized that Arg conformation change from a partly folded state to an extended state through the rotation of CC single bonds of Arg side chain during the molecular recognition process.     

De novo Design of Selective Membrane‐Active Peptides by Enzymatic Control of Their Conformational Bias on the Cell Surface

Selectively targeting the membrane‐perturbing potential of peptides towards a distinct cellular phenotype allows one to target distinct populations of cells. We report the de novo design of a new class of peptide whose ability to perturb cellular membranes is coupled to an enzyme‐mediated shift in the folding potential of the peptide into its bioactive conformation. Cells rich in negatively charged surface components that also highly express alkaline phosphatase, for example many cancers, are susceptible to the action of the peptide. The unfolded, inactive peptide is dephosphorylated, shifting its conformational bias towards cell‐surface‐induced folding to form a facially amphiphilic membrane‐active conformer. The fate of the peptide can be further tuned by peptide concentration to affect either lytic or cell‐penetrating properties, which are useful for selective drug delivery. This is a new design strategy to afford peptides that are selective in their membrane‐perturbing activity.     

Importance of ligand bioactive conformation in the discovery of potent indole-diamide inhibitors of the hepatitis C virus NS5B

Significant advances have led to receptor induced-fit and conformational selection models for describing bimolecular recognition, but a more comprehensive view must evolve to also include ligand shape and conformational changes. Here, we describe an example where a ligand's "structural hinge" influences potency by inducing an "L-shape" bioactive conformation, and due to its solvent exposure in the complex, reasonable conformation-activity-relationships can be qualitatively attributed. From a ligand design perspective, this feature was exploited by successful linker hopping to an alternate "structural hinge" that led to a new and promising chemical series which matched the ligand bioactive conformation and the pocket bioactive space. Using a combination of X-ray crystallography, NMR and modeling with support from binding-site resistance mutant studies and photoaffinity labeling experiments, we were able to derive inhibitor-polymerase complexes for various chemical series.     

Reproducing the conformations of protein-bound ligands: a critical evaluation of several popular conformational searching tools

Several programs (Catalyst, Confort, Flo99, MacroModel, and Omega) that are commonly used to generate conformational ensembles have been tested for their ability to reproduce bioactive conformations. The ligands from thirty-two different ligand–protein complexes determined by high-resolution (le2.0 Å) X-ray crystallography have been analyzed. The Low-Mode Conformational Search method (with AMBER^* and the GB/SA hydration model), as implemented in MacroModel, was found to perform better than the other algorithms. The rule-based method Omega, which is orders of magnitude faster than the other methods, also gave reasonable results but were found to be dependent on the input structure. The methods supporting diverse sampling (Catalyst, Confort) performed least well. For the seven ligands in the set having eight or more rotatable bonds, none of the bioactive conformations were ever found, save for one exception (Flo99). These ligands do not bind in a local minimum conformation according to AMBER^*\GB/SA. Taking these last two observations together, it is clear that geometrically similar structures should be collected in order to increase the probability of finding the bioactive conformation among the generated ensembles. Factors influencing bioactive conformational retrieval have been identified and are discussed.