Protein dynamics at Eph receptor-ligand interfaces as revealed by crystallography, NMR and MD simulations

An investigation of binding interactions of tumor-targeted peptide conjugated polyphenols with the kinase domain of ephrin B4 and B2 receptors

Recent studies have shown that Ephrin receptors may be upregulated in several types of cancers including breast, ovarian and endometrial cancers, making them a target for drug design. In this work, we have utilized a target-hopping approach to design new natural product-peptide conjugates and examined their interactions with the kinase-binding domain of EphB4 and EphB2 receptors. The peptide sequences were generated through point mutations of the known EphB4 antagonist peptide TNYLFSPNGPIA. Their anticancer properties and secondary structures were analyzed computationally. Conjugates of most optimum of peptides were then designed by binding the N-terminal of the peptides with the free carboxyl group of the polyphenols sinapate, gallate and coumarate, which are known for their inherent anticancer properties. To investigate if these conjugates have a potential to bind to the kinase domain, we carried out docking studies and MMGBSA free energy calculations of the trajectories based on the molecular dynamics simulations, with both the apo and the ATP bound kinase domains of both receptors. In most cases binding interactions occurred within the catalytic loop region, while in some cases the conjugates were found to spread out across the N-lobe and the DFG motif region. The conjugates were further tested for prediction of pharmacokinetic properties using ADME studies. Our results indicated that the conjugates were lipophilic and MDCK permeable with no CYP interactions. These findings provide an insight into the molecular interactions of these peptides and conjugates with the kinase domain of the EphB4 and EphB2 receptor. As a proof of concept, we synthesized and carried out SPR analysis with two of the conjugates (gallate-TNYLFSPNGPIA and sinapate-TNYLFSPNGPIA). Results indicated that the conjugates showed higher binding with the EphB4 receptor and minimal binding to EphB2 receptor. Sinapate-TNYLFSPNGPIA showed inhibitory activity against EphB4. These studies reveal that some of the conjugates may be developed for further investigation into in vitro and in vivo studies and potential development as therapeutics.

NMR-Based Methods for Protein Analysis

Nuclear magnetic resonance (NMR) spectroscopy is a well-established method for analyzing protein structure, interaction, and dynamics at atomic resolution and in various sample states including solution state, solid state, and membranous environment. Thanks to rapid NMR methodology development, the past decade has witnessed a growing number of protein NMR studies in complex systems ranging from membrane mimetics to living cells, which pushes the research frontier further toward physiological environments and offers unique insights in elucidating protein functional mechanisms. In particular, in-cell NMR has become a method of choice for bridging the huge gap between structural biology and cell biology. Herein, we review the recent developments and applications of NMR methods for protein analysis in close-to-physiological environments, with special emphasis on in-cell protein structural determination and the analysis of protein dynamics, both difficult to be accessed by traditional methods.

Exploring Protein Conformational Diversity

The native state of proteins is composed of conformers in dynamical equilibrium. In this chapter, different issues related to conformational diversity are explored using a curated and experimentally based database called CoDNaS (Conformational Diversity in the Native State). This database is a collection of redundant structures for the same sequence. CoDNaS estimates the degree of conformational diversity using different global and local structural similarity measures. It allows the user to explore how structural differences among conformers change as a function of several structural features providing further biological information. This chapter explores the measurement of conformational diversity and its relationship with sequence divergence. Also, it discusses how proteins with high conformational diversity could affect homology modeling techniques.

Structurally- and dynamically-driven allostery of the chymotrypsin-like proteases of SARS, Dengue and Zika viruses

Coronavirus 3C-like and Flavivirus NS2B-NS3 proteases utilize the chymotrypsin fold to harbor their catalytic machineries but also contain additional domains/co-factors. Over the past decade, we aimed to decipher how the extra domains/co-factors mediate the catalytic machineries of SARS 3C-like, Dengue and Zika NS2B-NS3 proteases by characterizing their folding, structures, dynamics and inhibition with NMR, X-ray crystallography and MD simulations, and the results revealed: 1) the chymotrypsin fold of the SARS 3C-like protease can independently fold, while, by contrast, those of Dengue and Zika proteases lack the intrinsic capacity to fold without co-factors. 2) Mutations on the extra domain of SARS 3C-like protease can transform the active catalytic machinery into the inactive collapsed state by structurally-driven allostery. 3) Amazingly, even without detectable structural changes, mutations on the extra domain are sufficient to either inactivate or enhance the catalytic machinery of SARS 3C-like protease by dynamically-driven allostery. 4) Global networks of correlated motions have been identified: for SARS 3C-like protease, N214A inactivates the catalytic machinery by decoupling the network, while STI/A and STIF/A enhance by altering the patterns of the network. The global networks of Dengue and Zika proteases are coordinated by their NS2B-cofactors. 5) Natural products were identified to allosterically inhibit Zika and Dengue proteases through binding a pocket on the back of the active site. Therefore, by introducing extra domains/cofactors, nature develops diverse strategies to regulate the catalytic machinery embedded on the chymotrypsin fold through folding, structurally- and dynamically-driven allostery, all of which might be exploited to develop antiviral drugs.

Non-equilibrium hydrogen exchange for determination of H-bond strength and water accessibility in solid proteins

We demonstrate measurement of non-equilibrium backbone amide hydrogen-deuterium exchange rates (HDX) for solid proteins. The target of this study are the slowly exchanging residues in solid samples, which are associated with stable secondary-structural elements of proteins. These hydrogen exchange processes escape methods measuring equilibrium exchange rates of faster processes. The method was applied to a micro-crystalline preparation of the SH3 domain of chicken α-spectrin. Therefore, from a 100% back-exchanged micro-crystalline protein preparation, the supernatant buffer was exchanged by a partially deuterated buffer to reach a final protonation level of approximately 20% before packing the sample in a 1.3 mm rotor. Tracking of the HN peak intensities for 2 weeks reports on site-specific hydrogen bond strength and also likely reflects water accessibility in a qualitative manner. H/D exchange can be directly determined for hydrogen-bonded amides using ¹H detection under fast magic angle spinning. This approach complements existing methods and provides the means to elucidate interesting site-specific characteristics for protein functionality in the solid state.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Dynamic principle for designing antagonistic/agonistic molecules for EphA4 receptor, the only known ALS modifier

Additional to involvement in diverse physiological and pathological processes such as axon regeneration, synaptic plasticity, and cancers, EphA4 receptor has been recently identified as the only amyotrophic lateral sclerosis (ALS) modifier. Previously, we found that two small molecules bind the same EphA4 channel at almost equivalent affinities but mysteriously trigger opposite signaling outputs: one activated but another inhibited. Here, we determined the solution structure of the 181-residue EphA4 LBD, which represents the first for 16 Eph receptors. Further NMR dynamic studies deciphered that the agonistic and antagonistic effects of two small molecules are dynamically driven, which are achieved by oppositely modulating EphA4 dynamics. Consequently, in design of drugs to target EphA4, the dynamic requirement also needs to be satisfied in addition to the classic criteria. For example, to increase the survival of ALS patients by inhibiting EphA4, the drugs must enhance, or at least not suppress, the EphA4 dynamics.

Development and structural analysis of a nanomolar cyclic peptide antagonist for the EphA4 receptor

The EphA4 receptor is highly expressed in the nervous system, and recent findings suggest that its signaling activity hinders neural repair and exacerbates certain neurodegenerative processes. EphA4 has also been implicated in cancer progression. Thus, EphA4 inhibitors represent potential therapeutic leads and useful research tools to elucidate the role of EphA4 in physiology and disease. Here, we report the structure of a cyclic peptide antagonist, APY, in complex with the EphA4 ligand-binding domain (LBD), which represents the first structure of a cyclic peptide bound to a receptor tyrosine kinase. The structure shows that the dodecameric APY efficiently occupies the ephrin ligand-binding pocket of EphA4 and promotes a "closed" conformation of the surrounding loops. Structure-guided relaxation of the strained APY β-turn and amidation of the C terminus to allow an additional intrapeptide hydrogen bond yielded APY-βAla8.am, an improved APY derivative that binds to EphA4 with nanomolar affinity. APY-βAla8.am potently inhibits ephrin-induced EphA4 activation in cells and EphA4-dependent neuronal growth cone collapse, while retaining high selectivity for EphA4. The two crystal structures of APY and APY-βAla8.am bound to EphA4, in conjunction with secondary phage display screens, highlighted peptide residues that are essential for EphA4 binding as well as residues that can be modified. Thus, the APY scaffold represents an exciting prototype, particularly since cyclic peptides have potentially favorable metabolic stability and are emerging as an important class of molecules for disruption of protein-protein interactions.

Investigation of the interactions between the EphB2 receptor and SNEW peptide variants

EphB2 interacts with cell surface-bound ephrin ligands to relay bidirectional signals. Overexpression of the EphB2 receptor protein has been linked to different types of cancer. The SNEW (SNEWIQPRLPQH) peptide binds with high selectivity and moderate affinity to EphB2, inhibiting Eph-ephrin interactions by competing with ephrin ligands for the EphB2 high-affinity pocket. We used rigorous free energy perturbation (FEP) calculations to re-evaluate the binding interactions of SNEW peptide with the EphB2 receptor, followed by experimental testing of the computational results. Our results provide insight into dynamic interactions of EphB2 with SNEW peptide. While the first four residues of the SNEW peptide are already highly optimized, change of the C-terminal end of the peptide has the potential to improve SNEW-binding affinity. We identified a PXSPY motif that can be similarly aligned with several other EphB2-binding peptides.

Promiscuous and specific recognition among ephrins and Eph receptors

Eph-ephrin interactions control the signal transduction between cells and play an important role in carcinogenesis and other diseases. The interactions between Eph receptors and ephrins of the same subclass are promiscuous; there are cross-interactions between some subclasses, but not all. To understand how Eph-ephrin interactions can be both promiscuous and specific, we investigated sixteen energy landscapes of four Eph receptors (A2, A4, B2, and B4) interacting with four ephrin ligands (A1, A2, A5, and B2). We generated conformational ensembles and recognition energy landscapes starting from separated Eph and ephrin molecules and proceeding up to the formation of Eph-ephrin complexes. Analysis of the Eph-ephrin recognition trajectories and the co-evolution entropy of 400 ligand binding domains of Eph receptor and 241 ephrin ligands identified conserved residues during the recognition process. Our study correctly predicted the promiscuity and specificity of the interactions and provided insights into their recognition. The dynamic conformational changes during Eph-ephrin recognition can be described by progressive conformational selection and population shift events, with two dynamic salt bridges between EphB4 and ephrin-B2 contributing to the specific recognition. EphA3 cancer-related mutations lowered the binding energies. The specificity is not only controlled by the final stage of the interaction across the protein-protein interface, but also has large contributions from binding kinetics with the help of dynamic intermediates along the pathway from the separated Eph and ephrin to the Eph-ephrin complex.

Tie2 and Eph receptor tyrosine kinase activation and signaling

The Eph and Tie cell surface receptors mediate a variety of signaling events during development and in the adult organism. As other receptor tyrosine kinases, they are activated on binding of extracellular ligands and their catalytic activity is tightly regulated on multiple levels. The Eph and Tie receptors display some unique characteristics, including the requirement of ligand-induced receptor clustering for efficient signaling. Interestingly, both Ephs and Ties can mediate different, even opposite, biological effects depending on the specific ligand eliciting the response and on the cellular context. Here we discuss the structural features of these receptors, their interactions with various ligands, as well as functional implications for downstream signaling initiation. The Eph/ephrin structures are already well reviewed and we only provide a brief overview on the initial binding events. We go into more detail discussing the Tie-angiopoietin structures and recognition.

Dynamic conformational switching in the chemokine ligand is essential for G-protein-coupled receptor activation

Chemokines mediate diverse functions from organogenesis to mobilizing leucocytes, and are unusual agonists for class-A GPCRs (G-protein-coupled receptors) because of their large size and multi-domain structure. The current model for receptor activation, which involves interactions between chemokine N-loop and receptor N-terminal residues (Site-I) and between chemokine N-terminal and receptor extracellular loop/transmembrane residues (Site-II), fails to describe differences in ligand/receptor selectivity and the activation of multiple signalling pathways. In the present study, we show in neutrophil-activating chemokine CXCL8 that the highly conserved GP (glycine-proline) motif located distal to both N-terminal and N-loop residues couples Site-I and Site-II interactions. GP mutants showed large differences from native-like to complete loss of function that could not be correlated with the specific mutation, receptor affinity or subtype, or a specific signalling pathway. NMR studies indicated that the GP motif does not influence Site-I interactions, but molecular dynamics simulations suggested that this motif dictates substates of the CXCL8 conformational ensemble. We conclude that the GP motif enables diverse receptor functions by controlling cross-talk between Site-I and Site-II, and further propose that the repertoire of chemokine functions is best described by a conformational ensemble model in which a network of long-range coupled indirect interactions mediate receptor activity.

Specialized Dynamical Properties of Promiscuous Residues Revealed by Simulated Conformational Ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein–protein interaction prediction and design methods.

Unique structure and dynamics of the EphA5 ligand binding domain mediate its binding specificity as revealed by X-ray crystallography, NMR and MD simulations

The 16 EphA and EphB receptors represent the largest family of receptor tyrosine kinases, and their interactions with 9 ephrin-A and ephrin-B ligands initiate bidirectional signals controlling many physiological and pathological processes. Most interactions occur between receptor and ephrins of the same class, and only EphA4 can bind all A and B ephrins. To understand the structural and dynamic principles that enable Eph receptors to utilize the same jellyroll β-sandwich fold to bind ephrins, the VAPB-MSP domain, peptides and small molecules, we have used crystallography, NMR and molecular dynamics (MD) simulations to determine the first structure and dynamics of the EphA5 ligand-binding domain (LBD), which only binds ephrin-A ligands. Unexpectedly, despite being unbound, the high affinity ephrin-binding pocket of EphA5 resembles that of other Eph receptors bound to ephrins, with a helical conformation over the J–K loop and an open pocket. The openness of the pocket is further supported by NMR hydrogen/deuterium exchange data and MD simulations. Additionally, the EphA5 LBD undergoes significant picosecond-nanosecond conformational exchanges over the loops, as revealed by NMR and MD simulations, but lacks global conformational exchanges on the microsecond-millisecond time scale. This is markedly different from the EphA4 LBD, which shares 74% sequence identity and 87% homology. Consequently, the unbound EphA5 LBD appears to comprise an ensemble of open conformations that have only small variations over the loops and appear ready to bind ephrin-A ligands. These findings show how two proteins with high sequence homology and structural similarity are still able to achieve distinctive binding specificities through different dynamics, which may represent a general mechanism whereby the same protein fold can serve for different functions. Our findings also suggest that a promising strategy to design agonists/antagonists with high affinity and selectivity might be to target specific dynamic states of the Eph receptor LBDs.

The underappreciated role of allostery in the cellular network

Allosteric propagation results in communication between distinct sites in the protein structure; it also encodes specific effects on cellular pathways, and in this way it shapes cellular response. One example of long-range effects is binding of morphogens to cell surface receptors, which initiates a cascade of protein interactions that leads to genome activation and specific cellular action. Allosteric propagation results from combinations of multiple factors, takes place through dynamic shifts of conformational ensembles, and affects the equilibria of macromolecular interactions. Here, we (a) emphasize the well-known yet still underappreciated role of allostery in conveying explicit signals across large multimolecular assemblies and distances to specify cellular action; (b) stress the need for quantitation of the allosteric effects; and finally, (c) propose that each specific combination of allosteric effectors along the pathway spells a distinct function. The challenges are colossal; the inspiring reward will be predicting function, misfunction, and outcomes of drug regimes.

Distinctive binding of three antagonistic peptides to the ephrin-binding pocket of the EphA4 receptor

The EphA4 receptor tyrosine kinase interacts with ephrin ligands to regulate many processes, ranging from axon guidance and nerve regeneration to cancer malignancy. Thus antagonists that inhibit ephrin binding to EphA4 could be useful for a variety of research and therapeutic applications. In the present study we characterize the binding features of three antagonistic peptides (KYL, APY and VTM) that selectively target EphA4 among the Eph receptors. Isothermal titration calorimetry analysis demonstrated that all three peptides bind to the ephrin-binding domain of EphA4 with low micromolar affinity. Furthermore, the effects of a series of EphA4 mutations suggest that the peptides interact in different ways with the ephrin-binding pocket of EphA4. Chemical-shift changes observed by NMR spectroscopy upon binding of the KYL peptide involve many EphA4 residues, consistent with extensive interactions and possibly receptor conformational changes. Additionally, systematic replacement of each of the 12 amino acids of KYL and VTM identify the residues critical for EphA4, binding. The peptides exhibit a long half-life in cell culture medium which, with their substantial binding affinity and selectivity for EphA4, makes them excellent research tools to modulate EphA4 function.

Structural and functional characterization of monomeric EphrinA1 binding site to EphA2 receptor

The EphA2 receptor is overexpressed in glioblastoma multiforme and has been to shown to contribute to cell transformation, tumor initiation, progression, and maintenance. EphrinA1 (eA1) is a preferred ligand for the receptor. Treatment with monomeric eA1, the form of eA1 found in the extracellular environment, causes receptor phosphorylation, internalization, and down-regulation with subsequent anti-tumor effects. Here, we investigated the structure-function relationship of a monomeric eA1 focusing on its G-H loop ((108)FQRFTPFTLGKEFKE(123)G), a highly conserved region among eAs that mediates binding to their receptors. Alanine substitution mutants of the G-H loop amino acids were transfected into U-251 MG glioblastoma multiforme cells, and functional activity of each mutant in conditioned media was assessed by EphA2 down-regulation, ERK and AKT activation and cellular response assays. Alanine substitutions at positions Pro-113 Thr-115, Gly-117, Glu-122, and also Gln-109 enhanced the EphA2 receptor down-regulation and decreased p-ERK and p-AKT. Substitution mutants of eA1 at positions Phe-108, Arg-110, Phe-111, Thr-112, Phe-114, Leu-116, Lys-118, Glu-119, and Phe-120 had a deleterious effect on EphA2 down-regulation when compared with eA1-WT. Mutants at positions Phe-108, Lys-18, Lys-121, Gly-123 retained similar properties to eA1-WT. Recombinant eA1-R110A, -T115A, -G117A, and -F120A have been found to exhibit the same characteristics as the ligands contained in the conditioned media mainly due to the differences in their binding to the receptor. Thus, we have identified variants of eA1 that possess either superagonistic or antagonistic properties. These new findings will be important in the understanding of the receptor/ligand interactions and in further design of anti-cancer therapies targeting the eA/EphA system.

Protein dynamics and conformational selection in bidirectional signal transduction

Protein conformational dynamics simultaneously allow promiscuity and specificity in binding. The multiple conformations of the free EphA4 ligand-binding domain observed in two new EphA4 crystal structures provide a unique insight into the conformational dynamics of EphA4 and its signaling pathways. The heterogeneous ensemble and loop dynamics explain how the EphA4 receptor is able to bind multiple A- and B-ephrin ligands and small molecules via conformational selection, which helps to fine-tune cellular signal response in both receptor and ligand cells.

See research article http://www.biomedcentral.com/2046-1682/5/2