Evolutionary sequence and structural basis for the distinct conformational landscapes of Tyr and Ser/Thr kinases.

Protein kinases are molecular machines with rich sequence variation that distinguishes the two main evolutionary branches - tyrosine kinases (TKs) from serine/threonine kinases (STKs). Using a sequence co-variation Potts statistical energy model we previously concluded that TK catalytic domains are more likely than STKs to adopt an inactive conformation with the activation loop in an autoinhibitory "folded" conformation, due to intrinsic sequence effects. Here we investigated the structural basis for this phenomenon by integrating the sequence-based model with structure-based molecular dynamics (MD) to determine the effects of mutations on the free energy difference between active and inactive conformations, using a novel thermodynamic cycle involving many (n=108) protein-mutation free energy perturbation (FEP) simulations in the active and inactive conformations. The sequence and structure-based results are consistent and support the hypothesis that the inactive conformation "DFG-out Activation Loop Folded", is a functional regulatory state that has been stabilized in TKs relative to STKs over the course of their evolution via the accumulation of residue substitutions in the activation loop and catalytic loop that facilitate distinct substrate binding modes in trans and additional modes of regulation in cis for TKs.

All-atom adaptively biased path optimization of Src kinase conformational inactivation: Switched electrostatic network in the concerted motion of αC helix and the activation loop.

A method to optimize a conformational pathway through a space of well-chosen reduced variables is employed to advance our understanding of protein conformational equilibrium. The adaptively biased path optimization strategy utilizes unrestricted, enhanced sampling in the region of a path in the reduced-variable space to identify a broad path between two stable end-states. Application to the inactivation transition of the Src tyrosine kinase catalytic domain reveals new insight into this well studied conformational equilibrium. The mechanistic description gained from identifying the motions and structural features along the path includes details of the switched electrostatic network found to underpin the transition. The free energy barrier along the path results from rotation of a helix, αC, that is tightly correlated with motions in the activation loop (A-loop) as well as distal regions in the C-lobe. Path profiles of the reduced variables clearly demonstrate the strongly correlated motions. The exchange of electrostatic interactions among residues in the network is key to these interdependent motions. In addition, the increased resolution from an all-atom model in defining the path shows multiple components for the A-loop motion and that different parts of the A-loop contribute throughout the length of the path.

Ubiquitin C-Terminal Hydrolase L1: Biochemical and Cellular Characterization of a Covalent Cyanopyrrolidine-Based Inhibitor.

The deubiquitinase (DUB) ubiquitin C-terminal hydrolase L1 (UCHL1) is expressed primarily in the central nervous system under normal physiological conditions. However, UCHL1 is overexpressed in various aggressive forms of cancer with strong evidence supporting UCHL1 as an oncogene in lung, glioma, and blood cancers. In particular, the level of UCHL1 expression in these cancers correlates with increased invasiveness and metastatic behavior, as well as poor patient prognosis. Although UCHL1 is considered an oncogene with potential as a therapeutic target, there remains a significant lack of useful small-molecule probes to pharmacologically validate in vivo targeting of the enzyme. Herein, we describe the characterization of a new covalent cyanopyrrolidine-based UCHL1 inhibitory scaffold in biochemical and cellular studies to better understand the utility of this inhibitor in elucidating the role of UCHL1 in cancer biology.

Substrate binding to Src: A new perspective on tyrosine kinase substrate recognition from NMR and molecular dynamics.

Most signal transduction pathways in humans are regulated by protein kinases through phosphorylation of their protein substrates. Typical eukaryotic protein kinases are of two major types: those that phosphorylate-specific sequences containing tyrosine (~90 kinases) and those that phosphorylate either serine or threonine (~395 kinases). The highly conserved catalytic domain of protein kinases comprises a smaller N lobe and a larger C lobe separated by a cleft region lined by the activation loop. Prior studies find that protein tyrosine kinases recognize peptide substrates by binding the polypeptide chain along the C-lobe on one side of the activation loop, while serine/threonine kinases bind their substrates in the cleft and on the side of the activation loop opposite to that of the tyrosine kinases. Substrate binding structural studies have been limited to four families of the tyrosine kinase group, and did not include Src tyrosine kinases. We examined peptide-substrate binding to Src using paramagnetic-relaxation-enhancement NMR combined with molecular dynamics simulations. The results suggest Src tyrosine kinase can bind substrate positioning residues C-terminal to the phosphoacceptor residue in an orientation similar to serine/threonine kinases, and unlike other tyrosine kinases. Mutagenesis corroborates this new perspective on tyrosine kinase substrate recognition. Rather than an evolutionary split between tyrosine and serine/threonine kinases, a change in substrate recognition may have occurred within the TK group of the human kinome. Protein tyrosine kinases have long been therapeutic targets, but many marketed drugs have deleterious off-target effects. More accurate knowledge of substrate interactions of tyrosine kinases has the potential for improving drug selectivity.

NmrLineGuru: Standalone and User-Friendly GUIs for Fast 1D NMR Lineshape Simulation and Analysis of Multi-State Equilibrium Binding Models.

The ability of high-resolution NMR spectroscopy to readout the response of molecular interactions at multiple atomic sites presents a unique capability to define thermodynamic equilibrium constants and kinetic rate constants for complex, multiple-step biological interactions. Nonetheless, the extraction of the relevant equilibrium binding and rate constants requires the appropriate analysis of not only a readout that follows the equilibrium concentrations of typical binding titration curves, but also the lineshapes of NMR spectra. To best take advantage of NMR data for characterizing molecular interactions, we developed NmrLineGuru, a software tool with a user-friendly graphical user interface (GUI) to model two-state, three-state, and four-state binding processes. Application of NmrLineGuru is through stand-alone GUIs, with no dependency on other software and no scripted input. NMR spectra can be fitted or simulated starting with user-specified input parameters and a chosen kinetic model. The ability to both simulate and fit NMR spectra provides the user the opportunity to not only determine the binding parameters that best reproduce the measured NMR spectra for the selected kinetic model, but to also query the possibility that alternative models agree with the data. NmrLineGuru is shown to provide an accurate, quantitative analysis of complex molecular interactions.

Application of a Substrate-Mediated Selection with c-Src Tyrosine Kinase to a DNA-Encoded Chemical Library.

As aberrant activity of protein kinases is observed in many disease states, these enzymes are common targets for therapeutics and detection of activity levels. The development of non-natural protein kinase substrates offers an approach to protein substrate competitive inhibitors, a class of kinase inhibitors with promise for improved specificity. Also, kinase activity detection approaches would benefit from substrates with improved activity and specificity. Here, we apply a substrate-mediated selection to a peptidomimetic DNA-encoded chemical library for enrichment of molecules that can be phosphorylated by the protein tyrosine kinase, c-Src. Several substrates were identified and characterized for activity. A lead compound (SrcDEL10) showed both the ability to serve as a substrate and to promote ATP hydrolysis by the kinase. In inhibition assays, compounds displayed IC50's ranging from of 8-100 µM. NMR analysis of SrcDEL10 bound to the c-Src:ATP complex was conducted to characterize the binding mode. An ester derivative of the lead compound demonstrated cellular activity with inhibition of Src-dependent signaling in cell culture. Together, the results show the potential for substrate-mediated selections of DNA-encoded libraries to discover molecules with functions other than simple protein binding and offer a new discovery method for development of synthetic tyrosine kinase substrates.

Protein Conformational Transitions from All-Atom Adaptively Biased Path Optimization.

Simulation methods are valuable for elucidating protein conformational transitions between functionally diverse states given that transition pathways are difficult to capture experimentally. Nonetheless, specific computational algorithms are required because of the high free energy barriers between these different protein conformational states. Adaptively biased path optimization (ABPO) is an unrestrained, transition-path optimization method that works in a reduced-variable space to construct an adaptive biasing potential to aid convergence. ABPO was previously applied using a coarse-grained Go̅-model to study conformational activation of Lyn, a Src family tyrosine kinase. How effectively ABPO can be applied at the higher resolution of an all-atom model to explore protein conformational transitions is not yet known. Here, we report the all-atom conformational transition paths of three protein systems constructed using the ABPO methodology. Two systems, triose phosphate isomerase and dihydrofolate reductase, undergo local flipping of a short loop that promotes ligand binding. The third system, estrogen receptor α ligand binding domain, has a helix that adopts different conformations when the protein is bound to an agonist or an antagonist. For each protein, distance-based or torsion-angle reduced variables were identified from unbiased trajectories. ABPO was computed in this reduced variable space to obtain the transition path between the two states. The all-atom ABPO is shown to successfully converge an optimal transition path for each of the three systems.

Entropic allostery dominates the phosphorylation-dependent regulation of Syk tyrosine kinase release from immunoreceptor tyrosine-based activation motifs.

Spleen tyrosine kinase (Syk) is an essential player in immune signaling through its ability to couple multiple classes of membrane immunoreceptors to intracellular signaling pathways. Ligand binding leads to the recruitment of Syk to a phosphorylated cytoplasmic region of the receptors called ITAM. Syk binds to ITAM with high-affinity (nanomolar Kd ) via its tandem pair of SH2 domains. The affinity between Syk and ITAM is allosterically regulated by phosphorylation at Y130 in a linker connecting the tandem SH2 domains; when Y130 is phosphorylated, the binding affinity decreases (micromolar Kd ). Previous equilibrium binding studies attribute the increase in the binding free energy to an intra-molecular binding (isomerization) step of the tandem SH2 and ITAM, but a physical basis for the increased free energy is unknown. Here, we provide evidence that Y130 phosphorylation imposes an entropy penalty to isomerization, but surprisingly, has negligible effect on the SH2 binding interactions with ITAM and thus on the binding enthalpy. An analysis of NMR chemical shift differences characterized conformational effects of ITAM binding, and binding thermodynamics were measured from isothermal titration calorimetry. Together the data support a previously unknown mechanism for the basis of regulating protein-protein interactions through protein phosphorylation. The decreased affinity for Syk association with immune receptor ITAMs by Y130 phosphorylation is an allosteric mechanism driven by an increased entropy penalty, likely contributed by conformational disorder in the SH2-SH2 inter-domain structure, while SH2-ITAM binding contacts are not affected, and binding enthalpy is unchanged.

A minimization principle for transition paths of maximum flux for collective variables.

Considered is the construction of transition paths of conformational changes for proteins and other macromolecules, using methods that do not require the generation of dynamics trajectories. Special attention is given to the use of a reduced set of collective variables for describing such paths. A favored way to define transition paths is to seek channels through the transition state having cross sections with a high reactive flux (density of last hitting points of reactive trajectories). Given here is a formula for reactive flux that is independent of the parameterization of "collective variable space." This formula is needed for the principal curve of the reactive flux (as in the revised finite temperature string method) and for the maximum flux transition (MaxFlux) path. Additionally, a resistance functional is derived for narrow tubes, which when minimized yields a MaxFlux path. A strategy for minimization is outlined in the spirit of the string method. Finally, alternative approaches based on determining trajectories of high probability are considered, and it is observed that they yield paths that depend on the parameterization of collective variable space, except in the case of zero temperature, where such a path coincides with a MaxFlux path.

A simple method for NMR t1 noise suppression.

t1 noise appears as random or semi-random spurious streaks along the indirect t1 (F1) dimension of a 2D or nD NMR spectrum. It can significantly downgrade spectral quality, especially for spectra with strong diagonal signals such as NOESY, because useful and weak cross-peaks can be easily buried under t1 noise. One of the significant contributing factors to t1 noise is unwanted and semi-random F2 signal modulation during t1 acquisition. As such, t1 noise from different acquisitions is unlikely to correlate with each other strongly. In the case of NOESY, co-addition of multiple spectra significantly reduces t1 noise compared with conventional acquisition with the same amount of total acquisition time and resolution.

Analysis of Multidomain Protein Dynamics.

Proteins with a modular architecture of multiple domains connected by linkers often exhibit diversity in the relative positions of domains, while the domain tertiary structure remains unchanged. The biological function of these modular proteins, or the regulation of their activity, depends on the variation in domain orientation and separation. Accordingly, careful characterization of interdomain motion and correlated fluctuations of multidomain systems is relevant for understanding the functional behavior of modular proteins. Molecular dynamics (MD) simulations provides a powerful approach to study these motions in atomic detail. Nevertheless, the common procedure for analyzing fluctuations from MD simulations after rigid-body alignment fails for multidomain proteins; it greatly overestimates correlated positional fluctuations in the presence of relative domain motion. We show here that expressing the atomic motions of a multidomain protein as a combination of displacement within the domain reference frame and motion of the relative domains correctly separates the internal motions to allow a useful description of correlated fluctuations. We illustrate the methodology of separating the domain fluctuations and local fluctuations by application to the tandem SH2 domains of human Syk protein kinase and by characterizing an effect of phosphorylation on the dynamics. Correlated motions are assessed from a distance covariance rather than the more common vector-coordinate covariance. The approach makes it possible to calculate the proper correlations in fluctuations internal to a domain as well as between domains.

Insights into the allosteric regulation of Syk association with receptor ITAM, a multi-state equilibrium.

The phosphorylation of interdomain A (IA), a linker region between tandem SH2 domains of Syk tyrosine kinase, regulates the binding affinity for association of Syk with doubly-phosphorylated ITAM regions of the B cell receptor. The mechanism of this allosteric regulation has been suggested to be a switch from the high-affinity bifunctional binding, mediated through both SH2 domains binding two phosphotyrosine residues of ITAM, to a substantially lower-affinity binding of only one SH2 domain. IA phosphorylation triggers the switch by inducing disorder in IA and weakening the SH2-SH2 interaction. The postulated switch to a single-SH2-domain binding mode is examined using NMR to monitor site-specific binding to each SH2 domain of Syk variants engineered to have IA regions that differ in conformational flexibility. The combined analysis of titration curves and NMR line-shapes provides sufficient information to determine the energetics of inter-molecular binding at each SH2 site along with an intra-molecular binding or isomerization step. A less favorable isomerization equilibrium associated with the changes in the SH2-SH2 conformational ensemble and IA flexibility accounts for the inhibition of Syk association with membrane ITAM regions when IA is phosphorylated, and refutes the proposed switch to single-SH2-domain binding. Syk localizes in the cell through its SH2 interactions, and this basis for allosteric regulation of ITAM association proposes for the first time a phosphorylation-dependent model to regulate Syk binding to alternate receptors and other signaling proteins that differ either in the number of residues separating ITAM phosphotyrosines or by having only one phosphotyrosine, a half ITAM.

Evaluating the dynamics and electrostatic interactions of folded proteins in implicit solvents.

Three implicit solvent models, namely GBMVII, FACTS, and SCPISM, were evaluated for their abilities to emulate an explicit solvent environment by comparing the simulated conformational ensembles, dynamics, and electrostatic interactions of the Src SH2 domain and the Lyn kinase domain. This assessment in terms of structural features in folded proteins expands upon the use of hydration energy as a metric for comparison. All-against-all rms coordinate deviation, average positional fluctuations, and ion-pair distance distribution were used to compare the implicit solvent models with the TIP3P explicit solvent model. Our study shows that the Src SH2 domains solvated with TIP3P, GBMVII, and FACTS sample similar global conformations. Additionally, the Src SH2 ion-pair distance distributions of solvent-exposed side chains corresponding to TIP3P, GBMVII, and FACTS do not differ substantially, indicating that GBMVII and FACTS are capable of modeling these electrostatic interactions. The ion-pair distance distributions of SCPISM are distinct from others, demonstrating that these electrostatic interactions are not adequately reproduced with the SCPISM model. On the other hand, for the Lyn kinase domain, a non-globular protein with bilobal structure and a large concavity on the surface, implicit solvent does not accurately model solvation to faithfully reproduce partially buried electrostatic interactions and lobe-lobe conformations. Our work reveals that local structure and dynamics of small, globular proteins are modeled well using FACTS and GBMVII. Nonetheless, global conformations and electrostatic interactions in concavities of multi-lobal proteins resulting from simulations with implicit solvent models do not match those obtained from explicit water simulations.

Effects of impaired membrane interactions on α-synuclein aggregation and neurotoxicity.

The post-mortem brains of individuals with Parkinson's disease (PD) and other synucleinopathy disorders are characterized by the presence of aggregated forms of the presynaptic protein α-synuclein (aSyn). Understanding the molecular mechanism of aSyn aggregation is essential for the development of neuroprotective strategies to treat these diseases. In this study, we examined how interactions between aSyn and phospholipid vesicles influence the protein's aggregation and toxicity to dopaminergic neurons. Two-dimensional NMR data revealed that two familial aSyn mutants, A30P and G51D, populated an exposed, membrane-bound conformer in which the central hydrophobic region was dissociated from the bilayer to a greater extent than in the case of wild-type aSyn. A30P and G51D had a greater propensity to undergo membrane-induced aggregation and elicited greater toxicity to primary dopaminergic neurons compared to the wild-type protein. In contrast, the non-familial aSyn mutant A29E exhibited a weak propensity to aggregate in the presence of phospholipid vesicles or to elicit neurotoxicity, despite adopting a relatively exposed membrane-bound conformation. Our findings suggest that the aggregation of exposed, membrane-bound aSyn conformers plays a key role in the protein's neurotoxicity in PD and other synucleinopathy disorders.

Relative Binding Enthalpies from Molecular Dynamics Simulations Using a Direct Method.

The potential for reliably predicting relative binding enthalpies, ΔΔE, from a direct method utilizing molecular dynamics is examined for a system of three phosphotyrosyl peptides binding to a protein receptor, the Src SH2 domain. The binding enthalpies were calculated from the potential energy differences between the bound and the unbound end-states of each peptide from equilibrium simulations in explicit water. The statistical uncertainties in the ensemble-mean energy values from multiple, independent simulations were obtained using a bootstrap method. Simulations were initiated with different starting coordinates as well as different velocities. Statistical uncertainties in ΔΔE are 2 to 3 kcal/mol based on calculations from 40, 10 ns trajectories for each system (three SH2-peptide complexes or unbound peptides). Uncertainties in relative component energies, comprising solute-solute, solute-solvent and solvent-solvent interactions, are considerably larger. Energy values were estimated from an unweighted ensemble averaging of multiple trajectories with the a priori assumption that all trajectories are equally likely. Distributions in energy-rmsd space indicate that the trajectories sample the same basin and the difference in mean energy values between trajectories is due to sampling of alternative local regions of this superbasin. The direct estimate of relative binding enthalpies is concluded to be a reasonable approach for well-ordered systems with ΔΔE values greater than ∼3 kcal/mol, although the approach would benefit from future work to determine properly distributed starting points that would enable efficient sampling of conformational space using multiple trajectories.

Prediction of protein relative enthalpic stability from molecular dynamics simulations of the folded and unfolded states.

For proteins of known structure, the relative enthalpic stability with respect to wild-type, ΔΔH(U), can be estimated by direct computation of the folded and unfolded state energies. We propose a model by which the change in stability upon mutation can be predicted from all-atom molecular dynamics simulations for the folded state and a peptide-based model for the unfolded state. The unfolding enthalpies are expressed in terms of environmental and hydration-solvent reorganization contributions that readily allow a residue-specific analysis of ΔΔH(U). The method is applied to estimate the relative enthalpic stability of variants with buried charged groups in T4 lysozyme. The predicted relative stabilities are in good agreement with experimental data. Environmental factors are observed to contribute more than hydration to the overall ΔΔH(U). The residue-specific analysis finds that the effects of burying charge are both localized and long-range. The enthalpy for hydration-solvent reorganization varies considerably among different amino-acid types, but because the variant folded state structures are similar to those of the wild-type, the hydration-solvent reorganization contribution to ΔΔH(U) is localized at the mutation site, in contrast to environmental contributions. Overall, mutation of apolar and polar amino acids to charged amino acids are destabilizing, but the reasons are complex and differ from site to site.

Unrestrained computation of free energy along a path.

We apply the adaptive biasing potential (ABP) method to optimize the principal curve defining a conformational transition between two known end states and to subsequently compute the one-dimensional potential of mean force as a function of arc length along the principal curve. This approach allows the use of the ABP method in a collective variable space of arbitrary dimension and offers several advantages over line-search methods. First, configurations are neither generated along an initial path for the transition nor equilibrated during evolution of the path. Second, and most importantly, the powerful sampling provided by the ABP serves to accelerate the dynamics during the optimization and computation of the free energy. Finally, the free energy is formulated as a potential of mean force that captures changes in the reaction channel along the principal curve, in contrast to the free energy profile evaluated from the local free-energy gradient in restrained path optimization methods. We first demonstrate the ABP formulation of path optimization using a two-dimensional potential surface and then with a more complex system of Src protein tyrosine kinase. The method is shown to be efficient and robust in the case of rugged, free-energy landscapes.

Long-distance correlations of rhinovirus capsid dynamics contribute to uncoating and antiviral activity.

Human rhinovirus (HRV) and other members of the enterovirus genus bind small-molecule antiviral compounds in a cavity buried within the viral capsid protein VP1. These compounds block the release of the viral protein VP4 and RNA from inside the capsid during the uncoating process. In addition, the antiviral compounds prevent "breathing" motions, the transient externalization of the N-terminal regions of VP1 and VP4 from the inside of intact viral capsid. The site for externalization of VP1/VP4 or release of RNA is likely between protomers, distant to the binding cavity for antiviral compounds. Molecular dynamics simulations were conducted to explore how the antiviral compound, WIN 52084, alters properties of the HRV 14 capsid through long-distance effect. We developed an approach to analyze capsid dynamics in terms of correlated radial motion and the shortest paths of correlated motions. In the absence of WIN, correlated radial motion is observed between residues separated by as much as 85 Å, a remarkably long distance. The most frequently populated path segments of the network were localized near the fivefold symmetry axis and included those connecting the N termini of VP1 and VP4 with other regions, in particular near twofold symmetry axes, of the capsid. The results provide evidence that the virus capsid exhibits concerted long-range dynamics, which have not been previously recognized. Moreover, the presence of WIN destroys this radial correlation network, suggesting that the underlying motions contribute to a mechanistic basis for the initial steps of VP1 and VP4 externalization and uncoating.