Design and synthesis of biologically active peptides: a 'tail' of amino acids can modulate activity of synthetic cyclic peptides

In earlier work, we synthesized a cyclic 9-amino acid peptide (AFPep, cyclo[EKTOVNOGN]) and showed it to be useful for prevention and therapy of breast cancer. In an effort to explore the structure–function relationships of AFPep, we have designed analogs that bear a short ‘tail’ (one or two amino acids) attached to the cyclic peptide distal to its pharmacophore. Analogs that bore a tail of either one or two amino acids, either of which had a hydrophilic moiety in the side chain (e.g., cyclo[EKTOVNOGN]FS) exhibited greatly diminished biological activity (inhibition of estrogen-stimulated uterine growth) relative to AFPep. Analogs that bore a tail of either one or two amino acids which had hydrophobic (aliphatic or aromatic) side chains (e.g., cyclo[EKTOVNOGN]FI) retained (or had enhanced) growth inhibition activity. Combining in the same biological assay a hydrophilic-tailed analog with either AFPep or a hydrophobic-tailed analog resulted in decreased activity relative to that for AFPep or for the hydrophobic-tailed analog alone, suggesting that hydrophilic-tailed analogs are binding to a biologically active receptor. An analog with a disrupted pharmacophore (cyclo[EKTOVGOGN]) exhibited little or no growth inhibition activity. An analog with a hydrophilic tail and a disrupted pharmacophore (cyclo[EKTOVGOGN]FS) exhibited no growth inhibition activity of its own and did not affect the activity of a hydrophobic-tailed analog, but enhanced the growth inhibition activity of AFPep. These results are discussed in the context of a two-receptor model for binding of AFPep and ring-and-tail analogs. We suggest that tails on cyclic peptides may comprise a useful method to enhance diversity of peptide design and specificity of ligand–receptor interactions.     

Limulus polyphemus hemocyanin: 10 A cryo-EM structure, sequence analysis, molecular modelling and rigid-body fitting reveal the interfaces between the eight hexamers

The blue copper protein hemocyanin from the horseshoe crab Limulus polyphemus is among the largest respiratory proteins found in nature (3.5 MDa) and exhibits a highly cooperative oxygen binding. Its 48 subunits are arranged as eight hexamers (1x6mers) that form the native 8x6mer in a nested hierarchy of 2x6mers and 4x6mers. This quaternary structure is established by eight subunit types (termed I, IIA, II, IIIA, IIIB, IV, V, and VI), of which only type II has been sequenced. Crystal structures of the 1x6mer are available, but for the 8x6mer only a 40 A 3D reconstruction exists. Consequently, the structural parameters of the 8x6mer are not firmly established, and the molecular interfaces between the eight hexamers are still to be defined. This, however, is crucial for understanding how allosteric transitions are mediated between the different levels of hierarchy. Here, we show the 10 A structure (FSC(1/2-bit) criterion) of the oxygenated 8x6mer from cryo-electron microscopy (cryo-EM) and single-particle analysis. Moreover, we show its molecular model as obtained by DNA sequencing of subunits II, IIIA, IV and VI, and molecular modelling and rigid-body fitting of all subunit types. Remarkably, the latter enabled us to improve the resolution of the cryo-EM structure from 11 A to the final 10 A. The 10 A structure allows firm assessment of various structural parameters of the 8x6mer, the 4x6mer and the 2x6mer, and reveals a total of 46 inter-hexamer bridges. These group as 11 types of interface: four at the 2x6mer level (II-II, II-IV, V-VI, IV-VI), three form the 4x6mer (V-V, V-VI, VI-IIIB/IV/V), and four are required to assemble the 8x6mer (IIIA-IIIA, IIIA-IIIB, II-IV, IV-IV). The molecular model shows the amino acid residues involved, and reveals that several of the interfaces are intriguingly histidine-rich and likely to transfer allosteric signals between the different levels of the nested hierarchy.     

Helix packing of leucine-rich peptides: a parallel leucine ladder in the structure of Boc-Aib-Leu-Aib-Aib-Leu-Leu-Leu-Aib-Leu-Aib-OMe.

The packing of peptide helices in crystals of the leucine-rich decapeptide Boc-Aib-Leu-Aib-Aib-Leu-Leu-Leu-Aib-Leu-Aib-OMe provides an example of ladder-like leucylleucyl interactions between neighboring molecules. The peptide molecule forms a helix with five 5----1 hydrogen bonds and two 4----1 hydrogen bonds near the C terminus. Three head-to-tail NH ... O = C hydrogen bonds between helices form continuous columns of helices in the crystal. The helicial columns associate in an antiparallel fashion, except for the association of Leu ... Leu side chains, which occurs along the diagonal of the cell where the peptide helices are parallel. The peptide, with formula C56H102N10O13, crystallizes in space group P2(1)2(1)2(1) with Z = 4 and cell parameters a = 16.774(3) A, b = 20.032(3) A and c = 20.117(3) A; overall agreement factor R = 10.7% for 2014 data with magnitude of F(obs) greater than 3 sigma (F); resolution 1.0 A.     

Structure and dynamics of LC8 complexes with KXTQT-motif peptides: swallow and dynein intermediate chain compete for a common site

The dynein light chain LC8 is an integral subunit of the cytoplasmic dynein motor complex that binds directly to and promotes assembly of the dynein intermediate chain (IC). LC8 interacts also with a variety of putative dynein cargo molecules such as Bim, a proapoptotic Bcl2 family protein, which have the KXTQT recognition sequence and neuronal nitric oxide synthase (nNOS), which has the GIQVD fingerprint but shares the same binding grooves at the LC8 dimer interface. The work reported here investigates the interaction of LC8 with IC and a putative cargo, Swallow, which share the KXTQT recognition sequence, and addresses the apparent paradox of how LC8, as part of dynein, mediates binding to cargo. The structures of Drosophila LC8 bound to peptides from IC and Swallow solved by X-ray diffraction show that the IC and Swallow peptides bind in the same grooves at the dimer interface. Differences in flexibility between bound and free LC8 were evaluated from hydrogen isotope exchange experiments using heteronuclear NMR spectroscopy. Peptide binding causes an increase in protection from exchange primarily in residues that interact directly with the peptide, such as the beta-strand intertwined at the interface and the N-terminal end of helix alpha2. There is considerably more protection upon Swallow binding, consistent with tighter binding relative to IC. Comparison with the LC8/nNOS complex shows how both the GIQVD and KXTQT fingerprints are recognized in the same groove. The similar structures of LC8/IC and LC8/Swa and the tighter binding of Swallow call into question the role for LC8 as a cargo adaptor protein, and suggest that binding of LC8 to Swallow serves another function, possibly that of a dimerization engine, which is independent of its role in dynein.     

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.     

Machine learning-enabled discovery and design of membrane-active peptides

Antimicrobial peptides are a class of membrane-active peptides that form a critical component of innate host immunity and possess a diversity of sequence and structure. Machine learning approaches have been profitably employed to efficiently screen sequence space and guide experiment towards promising candidates with high putative activity. In this mini-review, we provide an introduction to antimicrobial peptides and summarize recent advances in machine learning-enabled antimicrobial peptide discovery and design with a focus on a recent work Lee et al. Proc. Natl. Acad. Sci. USA 2016;113(48):13588–13593. This study reports the development of a support vector machine classifier to aid in the design of membrane active peptides. We use this model to discover membrane activity as a multiplexed function in diverse peptide families and provide interpretable understanding of the physicochemical properties and mechanisms governing membrane activity. Experimental validation of the classifier reveals it to have learned membrane activity as a unifying signature of antimicrobial peptides with diverse modes of action. Some of the discriminating rules by which it performs classification are in line with existing “human learned” understanding, but it also unveils new previously unknown determinants and multidimensional couplings governing membrane activity. Integrating machine learning with targeted experimentation can guide both antimicrobial peptide discovery and design and new understanding of the properties and mechanisms underpinning their modes of action.     

Origins of the specificity of inhibitor P218 toward wild-type and mutant PfDHFR: a molecular dynamics analysis

Molecular dynamics simulations were performed to evaluate the origin of the antimalarial effect of the lead compound P218. The simulations of the ligand in the cavities of wild-type, mutant Plasmodium falciparum Dihydrofolate Reductase (PfDHFR) and the human DHFR revealed the differences in the atomic-level interactions and also provided explanation for the specificity of this ligand toward PfDHFR. The binding free energy estimation using Molecular Mechanics Poisson-Boltzmann Surface Area method revealed that P218 has higher binding affinity (~ ?30 to ?35 kcal/mol) toward PfDHFR (both in wild-type and mutant forms) than human DHFR (~ ?22 kcal/mol), corroborating the experimental observations. Intermolecular hydrogen bonding analysis of the trajectories showed that P218 formed two stable hydrogen bonds with human DHFR (Ile7 and Glu30), wild-type and double-mutant PfDHFR’s (Asp54 and Arg122), while it formed three stable hydrogen bonds with quadruple-mutant PfDHFR (Asp54, Arg59, and Arg122). Additionally, P218 binding in PfDHFR is stabilized by hydrogen bonds with residues Ile14 and Ile164. It was found that mutant residues do not reduce the binding affinity of P218 to PfDHFR, in contrast, Cys59Arg mutation strongly favors inhibitor binding to quadruple-mutant PfDHFR. The atomistic-level details explored in this work will be highly useful for the design of non-resistant novel PfDHFR inhibitors as antimalarial agents.     

Structural signatures of enzyme binding pockets from order-independent surface alignment: a study of metalloendopeptidase and NAD binding proteins

Detecting similarities between local binding surfaces can facilitate identification of enzyme binding sites and prediction of enzyme functions, and aid in our understanding of enzyme mechanisms. Constructing a template of local surface characteristics for a specific enzyme function or binding activity is a challenging task, as the size and shape of the binding surfaces of a biochemical function often vary. Here we introduce the concept of signature binding pockets, which captures information on preserved and varied atomic positions at multiresolution levels. For proteins with complex enzyme binding and activity, multiple signatures arise naturally in our model, forming a signature basis set that characterizes this class of proteins. Both signatures and signature basis sets can be automatically constructed by a method called SOLAR (Signature Of Local Active Regions). This method is based on a sequence-order-independent alignment of computed binding surface pockets. SOLAR also provides a structure-based multiple sequence fragment alignment to facilitate the interpretation of computed signatures. By studying a family of evolutionarily related proteins, we show that for metzincin metalloendopeptidase, which has a broad spectrum of substrate binding, signature and basis set pockets can be used to discriminate metzincins from other enzymes, to predict the subclass of metzincins functions, and to identify specific binding surfaces. Studying unrelated proteins that have evolved to bind to the same NAD cofactor, we constructed signatures of NAD binding pockets and used them to predict NAD binding proteins and to locate NAD binding pockets. By measuring preservation ratio and location variation, our method can identify residues and atoms that are important for binding affinity and specificity. In both cases, we show that signatures and signature basis set reveal significant biological insight.     

Allelic variation,fragment length analyses and population genetic models: a case study on Drosophila microsatellites*

The allelic variation of 16 microsatellite loci from selected species of the Drosophila melanogaster and D. obscura group was determined. Intra‐ and interspecific sequence comparisons allowed discrimination of mutations affecting the repetitive microsatellite from those affecting the flanking regions. The hypotheses that slippage needs a minimum number of repeats in order to become efficient with respect to microsatellite variability, and of an increased mutation rate with increased length of the microsatellite are supported by the results of our analyses. There is in particular at the interrupted complex microsatellite locus BICOID in the species of the D. obscura group, extensive variation in the flanking regions in addition to length and sequence variation of the repetitive microsatellite. The allelic variation at this locus can hardly be explained by slippage alone. Estimates of microsatellite variability by fragment length analyses will pick up only a minor fraction of allelic variation at such loci, and conclusions that are based on the stepwise mutation model will not hold.     

Bilayer interaction and localization of cell penetrating peptides with model membranes: A comparative study of a human calcitonin (hCT)-derived peptide with pVEC and pAntp(43-58)

Cell-penetrating peptides (CPPs) are able to translocate problematic therapeutic cargoes across cellular membranes. The exact mechanisms of translocation are still under investigation. However, evidence for endocytic uptake is increasing. We investigated the interactions of CPPs with phospholipid bilayers as first step of translocation. To this purpose, we employed four independent techniques, comprising (i) liposome buffer equilibrium dialysis, (ii) Trp fluorescence quenching, (iii) fluorescence polarization, and (iv) determination of ζ-potentials. Using unilamellar vesicles (LUVs) of different phospholipid composition, we compared weakly cationic human calcitonin (hCT)-derived peptides with the oligocationic CPPs pVEC and penetratin (pAntp). Apparent partition coefficients of hCT-derived peptides in neutral POPC LUVs were dependent on amino acid composition and secondary structure; partitioning in negatively charged POPC/POPG (80:20) LUVs was increased and mainly governed by electrostatic interactions. For hCT(9-32) and its derivatives, D values raised from about 100-200 in POPC to about 1000 to 1500 when negatively charged lipids were present. Localization profiles of CPPs obtained by Trp fluorescence quenching were dependent on the charge density of LUVs. In POPC/POPG, hCT-derived CPPs were located on the bilayer surface, whereas pVEC and pAntp resided deeper in the membrane. In POPG LUVs, an increase of fluorescence polarization was observed for pVEC and pAntp but not for hCT-derived peptides. Generally, we found strong peptide-phospholipid interactions, especially when negatively charged lipids were present.     

Activation of the Ghrelin Receptor is Described by a Privileged Collective Motion: A Model for Constitutive and Agonist-induced Activation of a Sub-class A G-Protein Coupled Receptor (GPCR)

Three homology models of the human ghrelin receptor (GHS-R1a) have been generated from the available X-ray structures of rhodopsin (RHO model), opsin (OPS model) and beta-2 adrenergic receptor (B2 model). The latter was used as a starting point for combined molecular dynamics simulation (MDS) and full atom normal modes analysis (NMA). A low-frequency normal mode (mode 16) perfectly reproduced the intracellular motions observed between B2 and RHO models; in the opposite direction along the same mode, the generated structures are closer to the OPS model, suggesting a direct link with GHS-R1a activation. This was in agreement with motions of the seven transmembranous segments, increase of the solvent accessibility of the 140-ERY-142 sequence, and flip of the Trp276 (C WLP) residue, some features related to GPCRs activation. According to our model, His280 was proposed to stabilize Trp276 in the active state; this was verified by site-directed mutagenesis and biochemical characterization of the resulting H280A and H280S mutants, which were fully functional but sharing an important decrease of their basal activities. Docking performed with short ghrelin derivatives Gly-Ser-Ser _[octa]-Phe-NH ₂ and Gly-Ser-Ser _[octa]-Phe-Leu-NH ₂ allowed the identification of a robust position of these peptides in the active site of the receptor. This model was refined by MDS and validated by docking experiments performed on a set of 55 ghrelin receptor ligands based on the 1,2,4- triazole scaffold. Finally, NMA performed on the obtained peptide-receptor complex suggested stabilization of the Trp276 residue and of the whole receptor in the active state, preventing the motion observed along mode 16 computed for the unbound receptor. Our results show that NMA offers a powerful approach to study the conformational diversity and the activation mechanism of GPCRs.     

Molecular modelling and simulations in cancer research

The complexity of cancer and the vast amount of experimental data available have made computer-aided approaches necessary. Biomolecular modelling techniques are becoming increasingly easier to use, whereas hardware and software are becoming better and cheaper. Cross-talk between theoretical and experimental scientists dealing with cancer-research from a molecular approach, however, is still uncommon. This is in contrast to other fields, such as amyloid-related diseases, where molecular modelling studies are widely acknowledged. The aim of this review paper is therefore to expose some of the more common approaches in molecular modelling to cancer scientists in simple terms, illustrating success stories while also revealing the limitations of computational studies at the molecular level.     

Biphenyl amide p38 kinase inhibitors 3: Improvement of cellular and in vivo activity

The biphenyl amides (BPAs) are a novel series of p38α MAP kinase inhibitor. The optimisation of the series to give compounds that are potent in an in vivo disease model is discussed. SAR is presented and rationalised with reference to the crystallographic binding mode.     

Ascorbyl palmitate interaction with phospholipid monolayers: Electrostatic and rheological preponderancy

Ascorbyl palmitate (ASC₁₆) is an anionic amphiphilic molecule of pharmacological interest due to its antioxidant properties. We found that ASC₁₆ strongly interacted with model membranes. ASC₁₆ penetrated phospholipid monolayers, with a cutoff near the theoretical surface pressure limit. The presence of a lipid film at the interface favored ASC₁₆ insertion compared with a bare air/water surface. The adsorption and penetration time curves showed a biphasic behavior: the first rapid peak evidenced a fast adsorption of charged ASC₁₆ molecules to the interface that promoted a lowering of surface pH, thus partially neutralizing and compacting the film. The second rise represented an approach to the equilibrium between the ASC₁₆ molecules in the subphase and the surface monolayer, whose kinetics depended on the ionization state of the film. Based on the Langmuir dimiristoylphosphatidylcholine + ASC₁₆ monolayer data, we estimated an ASC₁₆ partition coefficient to dimiristoylphosphatidylcholine monolayers of 1.5 × 10⁵ and a ΔG_p = − 6.7 kcal·mol^− 1. The rheological properties of the host membrane were determinant for ASC₁₆ penetration kinetics: a fluid membrane, as provided by cholesterol, disrupted the liquid-condensed ASC₁₆-enriched domains and favored ASC₁₆ penetration. Subphase pH conditions affected ASC₁₆ aggregation in bulk: the smaller structures at acidic pHs showed a faster equilibrium with the surface film than large lamellar ones. Our results revealed that the ASC₁₆ interaction with model membranes has a highly complex regulation. The polymorphism in the ASC₁₆ bulk aggregation added complexity to the equilibrium between the surface and subphase form of ASC₁₆, whose understanding may shed light on the pharmacological function of this drug.     

Conservation of Functionally Important Global Motions in an Enzyme Superfamily across Varying Quaternary Structures

The α‐d‐phosphohexomutase superfamily comprises enzymes involved in carbohydrate metabolism that are found in all kingdoms of life. Recent biophysical studies have shown for the first time that several of these enzymes exist as dimers in solution, prompting an examination of the oligomeric state of all proteins of known structure in the superfamily (11 different proteins; 31 crystal structures) via computational and experimental analyses. We find that these proteins range in quaternary structure from monomers to tetramers, with 6 of the 11 known structures being likely oligomers. The oligomeric state of these proteins not only is associated in some cases with enzyme subgroup (i.e., substrate specificity) but also appears to depend on domain of life, with the two archaeal proteins existing as higher‐order oligomers. Within the oligomers, three distinct interfaces are observed, one of which is found in both archaeal and bacterial proteins. Normal mode analysis shows that the topological arrangement of the oligomers permits domain 4 of each protomer to move independently as required for catalysis. Our analysis suggests that the advantages associated with protein flexibility in this enzyme family are of sufficient importance to be maintained during the evolution of multiple independent oligomers. This study is one of the first showing that global motions may be conserved not only within protein families but also across members of a superfamily with varying oligomeric structures.