Inhibition of Aurora-A/N-Myc Protein-Protein Interaction Using Peptidomimetics: Understanding the Role of Peptide Cyclization.

Using N-Myc61-89 as a starting template we showcase the systematic use of truncation and maleimide constraining to develop peptidomimetic inhibitors of the N-Myc/Aurora-A protein-protein interaction (PPI); a potential anticancer drug discovery target. The most promising of these - N-Myc73-94-N85C/G89C-mal - is shown to favour a more Aurora-A compliant binding ensemble in comparison to the linear wild-type sequence as observed through fluorescence anisotropy competition assays, circular dichroism (CD) and nuclear magnetic resonance (NMR) experiments. Further in silico investigation of this peptide in its Aurora-A bound state, by molecular dynamics (MD) simulations, imply (i) the bound conformation is more stable as a consequence of the constraint, which likely suppresses dissociation and (ii) the constraint may make further stabilizing interactions with the Aurora-A surface. Taken together this work unveils the first orthosteric N-Myc/Aurora-A inhibitor and provides useful insights on the biophysical properties and thus design of constrained peptides, an attractive therapeutic modality.

α-Helix stabilization by co-operative side chain charge-reinforced interactions to phosphoserine in a basic kinase-substrate motif.

How cellular functions are regulated through protein phosphorylation events that promote or inhibit protein-protein interactions (PPIs) is key to understanding regulatory molecular mechanisms. Whilst phosphorylation can orthosterically or allosterically influence protein recognition, phospho-driven changes in the conformation of recognition motifs are less well explored. We recently discovered that clathrin heavy chain recognizes phosphorylated TACC3 through a helical motif that, in the unphosphorylated protein, is disordered. However, it was unclear whether and how phosphorylation could stabilize a helix in a broader context. In the current manuscript, we address this challenge using poly-Ala-based model peptides and a suite of circular dichroism and nuclear magnetic resonance spectroscopies. We show that phosphorylation of a Ser residue stabilizes the α-helix in the context of an Arg(i-3)pSeri Lys(i+4) triad through charge-reinforced side chain interactions with positive co-operativity, whilst phosphorylation of Thr induces an opposing response. This is significant as it may represent a general method for control of PPIs by phosphorylation; basic kinase-substrate motifs are common with 55 human protein kinases recognizing an Arg at a position -3 from the phosphorylated Ser, whilst the Arg(i-3)Seri Lys(i+4) is a motif found in over 2000 human proteins.

Dynamic ion pair behavior stabilizes single α-helices in proteins.

Ion pairs are key stabilizing interactions between oppositely charged amino acid side chains in proteins. They are often depicted as single conformer salt bridges (hydrogen-bonded ion pairs) in crystal structures, but it is unclear how dynamic they are in solution. Ion pairs are thought to be particularly important in stabilizing single α-helix (SAH) domains in solution. These highly stable domains are rich in charged residues (such as Arg, Lys, and Glu) with potential ion pairs across adjacent turns of the helix. They provide a good model system to investigate how ion pairs can contribute to protein stability. Using NMR spectroscopy, small-angle X-ray light scattering (SAXS), and molecular dynamics simulations, we provide here experimental evidence that ion pairs exist in a SAH in murine myosin 7a (residues 858-935), but that they are not fixed or long lasting. In silico modeling revealed that the ion pairs within this α-helix exhibit dynamic behavior, rapidly forming and breaking and alternating between different partner residues. The low-energy helical state was compatible with a great variety of ion pair combinations. Flexible ion pair formation utilizing a subset of those available at any one time avoided the entropic penalty of fixing side chain conformations, which likely contributed to helix stability overall. These results indicate the dynamic nature of ion pairs in SAHs. More broadly, thermodynamic stability in other proteins is likely to benefit from the dynamic behavior of multi-option solvent-exposed ion pairs.

Effect of external pulling forces on the length distribution of peptides.

BACKGROUND: The distribution of the length of a polypeptide, or that of the distance between any two of its atoms, is an important property as it can be analytically or numerically estimated for a number of polymer models. Importantly, it is directly measurable through a number of different experimental techniques. Length distributions can be straightforwardly assessed from molecular dynamics simulation; however, true convergence through full accurate coverage of the length range is difficult to achieve. METHODS: The application of external constant force combined with the weighted-histogram analysis method (WHAM) is used to enhance sampling of unlikely 'long' or 'short' conformations and obtain the potential of mean force, while also collecting dynamic properties of the chain under variable tension. RESULTS: We demonstrate the utility of constant force to enhance the sampling efficiency and obtain experimentally measurable quantities on a series of short peptides, including charge-rich sequences that are known to be highly helical but whose properties are distinct from those of helical peptides undergoing helix-coil transitions. CONCLUSIONS: Force-enhanced sampling enhances the range and accuracy of the length-based potential of mean force of the peptide, in particular those sequences that contain increased numbers of charged residues. GENERAL SIGNIFICANCE: This approach allows users to simultaneously probe the force-dependent behaviour of peptides directly, enhance the range and accuracy of the length-based PMF of the peptide and also test the convergence of simulations by comparing the overlap of PMF profiles from different constant forces. This article is part of a special issue entitled Recent developments of molecular dynamics.

Differential Effects of Hydrophobic Core Packing Residues for Thermodynamic and Mechanical Stability of a Hyperthermophilic Protein.

Proteins from organisms that have adapted to environmental extremes provide attractive systems to explore and determine the origins of protein stability. Improved hydrophobic core packing and decreased loop-length flexibility can increase the thermodynamic stability of proteins from hyperthermophilic organisms. However, their impact on protein mechanical stability is not known. Here, we use protein engineering, biophysical characterization, single-molecule force spectroscopy (SMFS), and molecular dynamics (MD) simulations to measure the effect of altering hydrophobic core packing on the stability of the cold shock protein TmCSP from the hyperthermophilic bacterium Thermotoga maritima. We make two variants of TmCSP in which a mutation is made to reduce the size of aliphatic groups from buried hydrophobic side chains. In the first, a mutation is introduced in a long loop (TmCSP L40A); in the other, the mutation is introduced on the C-terminal ß-strand (TmCSP V62A). We use MD simulations to confirm that the mutant TmCSP L40A shows the most significant increase in loop flexibility, and mutant TmCSP V62A shows greater disruption to the core packing. We measure the thermodynamic stability (ΔGD-N) of the mutated proteins and show that there is a more significant reduction for TmCSP L40A (ΔΔG = 63%) than TmCSP V62A (ΔΔG = 47%), as might be expected on the basis of the relative reduction in the size of the side chain. By contrast, SMFS measures the mechanical stability (ΔG*) and shows a greater reduction for TmCSP V62A (ΔΔG* = 8.4%) than TmCSP L40A (ΔΔG* = 2.5%). While the impact on the mechanical stability is subtle, the results demonstrate the power of tuning noncovalent interactions to modulate both the thermodynamic and mechanical stability of a protein. Such understanding and control provide the opportunity to design proteins with optimized thermodynamic and mechanical properties.

Tuning protein mechanics through an ionic cluster graft from an extremophilic protein.

Proteins from extremophilic organisms provide excellent model systems to determine the role of non-covalent interactions in defining protein stability and dynamics as well as being attractive targets for the development of robust biomaterials. Hyperthermophilic proteins have a prevalence of salt bridges, relative to their mesophilic homologues, which are thought to be important for enhanced thermal stability. However, the impact of salt bridges on the mechanical properties of proteins is far from understood. Here, a combination of protein engineering, biophysical characterisation, single molecule force spectroscopy (SMFS) and molecular dynamics (MD) simulations directly investigates the role of salt bridges in the mechanical stability of two cold shock proteins; BsCSP from the mesophilic organism Bacillus subtilis and TmCSP from the hyperthermophilic organism Thermotoga maritima. Single molecule force spectroscopy shows that at ambient temperatures TmCSP is mechanically stronger yet, counter-intuitively, its native state can withstand greater deformation before unfolding (i.e. it is mechanically soft) compared with BsCSP. MD simulations were used to identify the location and quantify the population of salt bridges, and reveal that TmCSP contains a larger number of highly occupied salt bridges than BsCSP. To test the hypothesis that salt-bridges endow these mechanical properties on the hyperthermophilic CSP, a charged triple mutant (CTM) variant of BsCSP was generated by grafting an ionic cluster from TmCSP into the BsCSP scaffold. As expected CTM is thermodynamically more stable and mechanically softer than BsCSP. We show that a grafted ionic cluster can increase the mechanical softness of a protein and speculate that it could provide a mechanical recovery mechanism and that it may be a design feature applicable to other proteins.

Stable single α-helices are constant force springs in proteins.

Single α-helix (SAH) domains are rich in charged residues (Arg, Lys, and Glu) and stable in solution over a wide range of pH and salt concentrations. They are found in many different proteins where they bridge two functional domains. To test the idea that their high stability might enable these proteins to resist unfolding along their length, the properties and unfolding behavior of the predicted SAH domain from myosin-10 were characterized. The expressed and purified SAH domain was highly helical, melted non-cooperatively, and was monomeric as shown by circular dichroism and mass spectrometry as expected for a SAH domain. Single molecule force spectroscopy experiments showed that the SAH domain unfolded at very low forces (<30 pN) without a characteristic unfolding peak. Molecular dynamics simulations showed that the SAH domain unfolds progressively as the length is increased and refolds progressively as the length is reduced. This enables the SAH domain to act as a constant force spring in the mechanically dynamic environment of the cell.

Myosin tails and single α-helical domains.

The human genome contains 39 myosin genes, divided up into 12 different classes. The structure, cellular function and biochemical properties of many of these isoforms remain poorly characterized and there is still some controversy as to whether some myosin isoforms are monomers or dimers. Myosin isoforms 6 and 10 contain a stable single α-helical (SAH) domain, situated just after the canonical lever. The SAH domain is stiff enough to be able to lengthen the lever allowing the myosin to take a larger step. In addition, atomic force microscopy and atomistic simulations show that SAH domains unfold at relatively low forces and have a high propensity to refold. These properties are likely to be important for protein function, enabling motors to carry cargo in dense actin networks, and other proteins to remain attached to binding partners in the crowded cell.

Life in extreme environments: single molecule force spectroscopy as a tool to explore proteins from extremophilic organisms.

Extremophiles are organisms which survive and thrive in extreme environments. The proteins from extremophilic single-celled organisms have received considerable attention as they are structurally stable and functionally active under extreme physical and chemical conditions. In this short article, we provide an introduction to extremophiles, the structural adaptations of proteins from extremophilic organisms and the exploitation of these proteins in industrial applications. We provide a review of recent developments which have utilized single molecule force spectroscopy to mechanically manipulate proteins from extremophilic organisms and the information which has been gained about their stability, flexibility and underlying energy landscapes.

Exploring the Conformational Landscape and Stability of Aurora A Using Ion-Mobility Mass Spectrometry and Molecular Modeling.

Protein kinase inhibitors are highly effective in treating diseases driven by aberrant kinase signaling and as chemical tools to help dissect the cellular roles of kinase signaling complexes. Evaluating the effects of binding of small molecule inhibitors on kinase conformational dynamics can assist in understanding both inhibition and resistance mechanisms. Using gas-phase ion-mobility mass spectrometry (IM-MS), we characterize changes in the conformational landscape and stability of the protein kinase Aurora A (Aur A) driven by binding of the physiological activator TPX2 or small molecule inhibition. Aided by molecular modeling, we establish three major conformations, the relative abundances of which were dependent on the Aur A activation status: one highly populated compact conformer similar to that observed in most crystal structures, a second highly populated conformer possessing a more open structure infrequently found in crystal structures, and an additional low-abundance conformer not currently represented in the protein databank. Notably, inhibitor binding induces more compact configurations of Aur A, as adopted by the unbound enzyme, with both IM-MS and modeling revealing inhibitor-mediated stabilization of active Aur A.

Biomembrane force probe investigation of RNA dissociation.

Investigations into the energy pathways of biomolecular interactions by use of dynamic force spectroscopy are limited by the range of loading rates accessible with a single technique. In the work discussed in this paper, this range has been extended for a previously studied system by using the biomembrane force probe (BFP). This work builds on our previous single-molecule atomic force microscopy (AFM) study of the dissociation of a bulge-motif-containing RNA complex. The disparity observed, at high loading rates, between the dissociation of a 12-base pair complex with and without a central three-base pair bulge was not observed at low rates. This suggests that the two species share a similar outer barrier to dissociation and that inclusion of the bulge motif creates an additional barrier at a distance closer to the bound state. Experiments performed in different buffer environments yielded similar results. The results, when combined with those of previous studies, suggest that the shared outer barrier to dissociation is that due to a rearrangement and fraying of the ends of the helix.

Vancomycin dimer formation between analogues of bacterial peptidoglycan surfaces probed by force spectroscopy.

Functionalised thiols presenting peptides found in the peptidoglycan of vancomycin-sensitive and -resistant bacteria were synthesised and used to form self-assembled monolayers (SAMs) on gold surfaces. This model bacterial cell-wall surface mimic was used to study binding interactions with vancomycin. Force spectroscopy, using the atomic force microscope (AFM), was used to investigate the specific rupture of interfacial vancomycin dimer complexes formed between pairs of vancomycin molecules bound to peptide-coated AFM probe and substrate surfaces. Clear adhesive contacts were observed between the vancomycin-sensitive peptide surfaces when vancomycin was present in solution, and the adhesion force demonstrated a clear dependence on antibiotic concentration.

Protein mechanics probed using simple molecular models.

BACKGROUND: Single-molecule experimental techniques such as optical tweezers or atomic force microscopy are a direct probe of the mechanical unfolding/folding of individual proteins. They are also a means to investigate free energy landscapes. Protein force spectroscopy alone provides limited information; theoretical models relate measurements to thermodynamic and kinetic properties of the protein, but do not reveal atomic level information. By building a molecular model of the protein and probing its properties through numerical simulation, one can gauge the response to an external force for individual interatomic interactions and determine structures along the unfolding pathway. In combination, single-molecule force probes and molecular simulations contribute to uncover the rich behavior of proteins when subjected to mechanical force. SCOPE OF REVIEW: We focus on how simplified protein models have been instrumental in showing how general properties of the free energy landscape of a protein relate to its response to mechanical perturbations. We discuss the role of simple protein models to explore the complexity of free energy landscapes and highlight important conceptual issues that more chemically accurate models with all-atom representations of proteins and solvent cannot easily address. MAJOR CONCLUSIONS: Native-centric, coarse-grained models, despite simplifications in chemical detail compared to all-atom models, can reproduce and interpret experimental results. They also highlight instances where the theoretical framework used to interpret single-molecule data is too simple. However, these simple models are not able to reproduce experimental findings where non-native contacts are involved. GENERAL SIGNIFICANCE: Mechanical forces are ubiquitous in the cell and it is increasingly clear that the way a protein responds to mechanical perturbation is important.

Helical Polyampholyte Sequences Have Unique Thermodynamic Properties.

Helices are the most common structural pattern observed in structured proteins. Polypeptide sequences that form helices in isolation have been identified and extensively studied. These are generally rich in alanine, the amino acid with strongest helical propensity. Insertion of charged or polar amino acids has been shown to be necessary to make alanine-rich peptides soluble and sometimes even increase the helicity of the peptides. More recently sequences that contain mostly charged residues (E-R/K rich) have been found in naturally occurring proteins that are highly helical, soluble, and extended regardless their length. Artificial sequences composed mostly or exclusively of charged amino acids have been designed that are also highly helical, depending on the specific pattern of oppositely charged residues. Here we explore the thermodynamic properties of a number of 16-residue long peptides with varying helical propensity by performing equilibrium simulations over a broad range of temperatures. We observe quantitative differences in the peptides' helical propensities that can be related to qualitative differences in the free energy landscape, depending on the ampholytic patterns in the sequence. The results provide hints on how the specific physical properties of naturally occurring long sequences with similar patterns of charged residues may relate to their biological function.

Determining Stable Single Alpha Helical (SAH) Domain Properties by Circular Dichroism and Atomic Force Microscopy.

Stable, single α-helical (SAH) domains exist in a number of unconventional myosin isoforms, as well as other proteins. These domains are formed from sequences rich in charged residues (Arg, Lys, and Glu), they can be hundreds of residues long, and in isolation they can tolerate significant changes in pH and salt concentration without loss in helicity. Here we describe methods for the preparation and purification of SAH domains and SAH domain-containing constructs, using the myosin 10 SAH domain as an example. We go on to describe the use of circular dichroism spectroscopy and force spectroscopy with the atomic force microscope for the elucidation of structural and mechanical properties of these unusual helical species.

A1603P and K1617del, Mutations in ß-Cardiac Myosin Heavy Chain that Cause Laing Early-Onset Distal Myopathy, Affect Secondary Structure and Filament Formation In Vitro and In Vivo.

Over 20 mutations in ß-cardiac myosin heavy chain (ß-MHC), expressed in cardiac and slow muscle fibers, cause Laing early-onset distal myopathy (MPD-1), a skeletal muscle myopathy. Most of these mutations are in the coiled-coil tail and commonly involve a mutation to a proline or a single-residue deletion, both of which are predicted to strongly affect the secondary structure of the coiled coil. To test this, we characterized the effects of two MPD-1 causing mutations: A1603P and K1617del in vitro and in cells. Both mutations affected secondary structure, decreasing the helical content of 15 heptad and light meromyosin constructs. Both mutations also severely disrupted the ability of glutathione S-transferase-light meromyosin fusion proteins to form minifilaments in vitro, as demonstrated by negative stain electron microscopy. Mutant eGFP-tagged ß-MHC accumulated abnormally into the M-line of sarcomeres in cultured skeletal muscle myotubes. Incorporation of eGFP-tagged ß-MHC into sarcomeres in adult rat cardiomyocytes was reduced. Molecular dynamics simulations using a composite structure of part of the coiled coil demonstrated that both mutations affected the structure, with the mutation to proline (A1603P) having a smaller effect compared to K1617del. Taken together, it seems likely that the MPD-1 mutations destabilize the coiled coil, resulting in aberrant myosin packing in thick filaments in muscle sarcomeres, providing a potential mechanism for the disease.

Nonexponential Kinetics of Loop Formation in Proteins and Peptides: A Signature of Rugged Free Energy Landscapes?

The kinetics of loop formation, i.e., the occurrence of contact between two atoms of a polypeptide, remains the focus of continuing interest. One of the reasons is that contact formation is the elementary event underlying processes such as folding and binding. More importantly, it is experimentally measurable and can be predicted theoretically for ideal polymers. Deviations from single exponential kinetics have sometimes been interpreted as a signature of rugged, protein-like, free energy landscapes. Here we present simulations, with different atomistic models, of short peptides with varied structural propensity, and of a structured protein. Results show exponential contact formation kinetics (or relaxation) at long times, and a power law relaxation at very short times. At intermediate times, a deviation from either power law or simple exponential kinetics is observed that appears to be characteristic of polypeptides with either specific or nonspecific attractive interactions but disappears if attractive interactions are absent. Our results agree with recent experimental measurements on peptides and proteins and offer a comprehensive interpretation for them.

Characterization of long and stable de novo single alpha-helix domains provides novel insight into their stability.

Naturally-occurring single α-helices (SAHs), are rich in Arg (R), Glu (E) and Lys (K) residues, and stabilized by multiple salt bridges. Understanding how salt bridges promote their stability is challenging as SAHs are long and their sequences highly variable. Thus, we designed and tested simple de novo 98-residue polypeptides containing 7-residue repeats (AEEEXXX, where X is K or R) expected to promote salt-bridge formation between Glu and Lys/Arg. Lys-rich sequences (EK3 (AEEEKKK) and EK2R1 (AEEEKRK)) both form SAHs, of which EK2R1 is more helical and thermo-stable suggesting Arg increases stability. Substituting Lys with Arg (or vice versa) in the naturally-occurring myosin-6 SAH similarly increased (or decreased) its stability. However, Arg-rich de novo sequences (ER3 (AEEERRR) and EK1R2 (AEEEKRR)) aggregated. Combining a PDB analysis with molecular modelling provides a rational explanation, demonstrating that Glu and Arg form salt bridges more commonly, utilize a wider range of rotamer conformations, and are more dynamic than Glu-Lys. This promiscuous nature of Arg helps explain the increased propensity of de novo Arg-rich SAHs to aggregate. Importantly, the specific K:R ratio is likely to be important in determining helical stability in de novo and naturally-occurring polypeptides, giving new insight into how single α-helices are stabilized.

Modification of commercially available image analysis software for semi-automated qualitative analysis of axon regeneration and myelination in the rat sciatic nerve.

BACKGROUND: Manual measurement of three routine parameters of axon regeneration and myelination, such as axon counts, axon caliber and G-Ratio, is tedious and time consuming. An automated or semi-automated computer program could improve the efficiency of this process. The concern, however, is that the automated method will lack the accuracy of manual counting and measuring. METHOD: We introduce semi-automated axon analysis software utilizing customized Image Pro Plus software designed to identify, count, and measure axons and their surrounding myelin sheaths. Histologic specimens were subsequently analyzed using either conventional manual techniques (manual counting and computer assisted measurements) or this semi-automated software (software performed counting and measurements with human supervision and adjustment of axon identification). RESULTS: G-Ratios, axon numbers, and axon diameter values did not differ between the two methods though the semi-automated method took approximately 50% less time. COMPARISON WITH EXISTING METHODS: Other described semi-automated nerve analysis tools require complex, original software development. CONCLUSIONS: Commercially available software can be modified to improve nerve histomorphology analysis time while maintaining accuracy.

Nanomechanical detection of antibiotic-mucopeptide binding in a model for superbug drug resistance.

The alarming growth of the antibiotic-resistant superbugs methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant Enterococcus (VRE) is driving the development of new technologies to investigate antibiotics and their modes of action. We report the label-free detection of vancomycin binding to bacterial cell wall precursor analogues (mucopeptides) on cantilever arrays, with 10 nM sensitivity and at clinically relevant concentrations in blood serum. Differential measurements have quantified binding constants for vancomycin-sensitive and vancomycin-resistant mucopeptide analogues. Moreover, by systematically modifying the mucopeptide density we gain new insights into the origin of surface stress. We propose that stress is a product of a local chemical binding factor and a geometrical factor describing the mechanical connectivity of regions activated by local binding in terms of a percolation process. Our findings place BioMEMS devices in a new class of percolative systems. The percolation concept will underpin the design of devices and coatings to significantly lower the drug detection limit and may also have an impact on our understanding of antibiotic drug action in bacteria.