Christina Kiel

Insights into Protein-Protein Binding by Binding Free Energy Calculation and Free Energy Decomposition for the Ras-Raf and Ras-RalGDS Complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These hotspots are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.

Engineering Signal Transduction Pathways

Cells respond to their environment by sensing signals and translating them into changes in gene expression. In recent years, synthetic networks have been designed in both prokaryotic and eukaryotic systems to create new functionalities and for specific applications. In this review, we discuss the challenges associated with engineering signal transduction pathways. Furthermore, we address advantages and disadvantages of engineering signaling pathways in prokaryotic and eukaryotic cells, highlighting recent examples, and discuss how progress in synthetic biology might impact biotechnology and biomedicine.

Analysis of disease-linked Rhodopsin mutations based on structure, function, and protein stability calculations

Retinitis pigmentosa (RP) refers to a heterogeneous group of inherited diseases that result in progressive retinal degeneration, characterized by visual field constriction and night blindness. A total of 103 mutations in rhodopsin are linked to RP to date, and the phenotypes range from severe to asymptomatic. To study the relation between phenotype and rhodopsin stability in disease mutants, we used a structure-based approach. For 12 of the mutants located at the protein-lipid interphase, we used the von Heijne water-membrane transfer scale, and we find that 9 of the mutations could affect membrane insertion. For 91 mutants, we used the protein design algorithm FoldX. The 3 asymptomatic mutations had no significant reduced stability, 2 were unsuitable for FoldX analysis since the structure was incorrect in this region, 63 mutations had a significant change in protein stability (>1.6xa0kcal/mol), and 23 mutations had energy change values under the prediction error threshold (<1.6xa0kcal/mol). Out of these 23, the disease-causing effect could be explained by the involvement in other functions (e.g., glycosylation motifs, the interface with arrestin and transducin, and the cilia-binding motif) for 19 mutants. The remaining 4 mutants were probably incorrectly associated with RP or have functionalities not discovered yet. For destabilizing mutations where clinical data were available, we found a highly significant correlation between FoldX energy changes and the average age of night blindness and between FoldX energy changes and daytime vision loss onset. Our detailed structural, functional, and energetic analysis provides a complete picture of the rhodopsin mutations and can guide mutation-specific therapies.

Structural fingerprints of the Ras-GTPase activating proteins neurofibromin and p120GAP

Ras specific GTPase activating proteins (GAPs), neurofibromin and p120GAP, bind GTP bound Ras and efficiently complement its active site. Here we present comparative data from mutations and fluorescence-based assays of the catalytic domains of both RasGAPs and interpret them using the crystal structures. Three prominent regions in RasGAPs, the arginine-finger loop, the phenylalanine-leucine-arginine (FLR) region and alpha7/variable loop contain structural fingerprints governing the GAP function. The finger loop is crucial for the stabilization of the transition state of the GTPase reaction. This function is controlled by residues proximal to the catalytic arginine, which are strikingly different between the two RasGAPs. These residues specifically determine the orientation and therefore the positioning of the arginine finger in the Ras active site. The invariant FLR region, a hallmark for RasGAPs, indirectly contributes to GTPase stimulation by forming a scaffold, which stabilizes Ras switch regions. We show that a long hydrophobic side-chain in the FLR region is crucial for this function. The alpha7/variable loop uses several conserved residues including two lysine residues, which are involved in numerous interactions with the switch I region of Ras. This region determines the specificity of the Ras-RasGAP interaction.

Leucine-Rich Repeat Kinase 2 Binds to Neuronal Vesicles through Protein Interactions Mediated by Its C-Terminal WD40 Domain

ABSTRACT Mutations in the leucine-rich repeat kinase 2 gene (LRRK2) are associated with familial and sporadic Parkinsons disease (PD). LRRK2 is a complex protein that consists of multiple domains, including predicted C-terminal WD40 repeats. In this study, we analyzed functional and molecular features conferred by the WD40 domain. Electron microscopic analysis of the purified LRRK2 C-terminal domain revealed doughnut-shaped particles, providing experimental evidence for its WD40 fold. We demonstrate that LRRK2 WD40 binds and sequesters synaptic vesicles via interaction with vesicle-associated proteins. In fact, a domain-based pulldown approach combined with mass spectrometric analysis identified LRRK2 as being part of a highly specific protein network involved in synaptic vesicle trafficking. In addition, we found that a C-terminal sequence variant associated with an increased risk of developing PD, G2385R, correlates with a reduced binding affinity of LRRK2 WD40 to synaptic vesicles. Our data demonstrate a critical role of the WD40 domain within LRRK2 function.

Structural and functional protein network analyses predict novel signaling functions for rhodopsin

Orchestration of signaling, photoreceptor structural integrity, and maintenance needed for mammalian vision remain enigmatic. By integrating three proteomic data sets, literature mining, computational analyses, and structural information, we have generated a multiscale signal transduction network linked to the visual G protein‐coupled receptor (GPCR) rhodopsin, the major protein component of rod outer segments. This network was complemented by domain decomposition of protein–protein interactions and then qualified for mutually exclusive or mutually compatible interactions and ternary complex formation using structural data. The resulting information not only offers a comprehensive view of signal transduction induced by this GPCR but also suggests novel signaling routes to cytoskeleton dynamics and vesicular trafficking, predicting an important level of regulation through small GTPases. Further, it demonstrates a specific disease susceptibility of the core visual pathway due to the uniqueness of its components present mainly in the eye. As a comprehensive multiscale network, it can serve as a basis to elucidate the physiological principles of photoreceptor function, identify potential disease‐associated genes and proteins, and guide the development of therapies that target specific branches of the signaling pathway.

Challenges ahead in signal transduction: MAPK as an example

Signal transduction is the process of converting one kind of signal or stimulus into another, with the goal of processing external or internal signals into diverse functions. In recent years, progress has been made in identifying signaling pathways, investigating cross-talk and feedbacks in signaling modules, analyzing cell-type specific signaling, and unraveling spatial-temporal aspects, such as receptor clustering into micro-domains, dynamic localization, and mathematical modeling. Here, we used the receptor/MAPK signaling system as an example, and we discuss the current and remaining challenges: the role of scaffolds and signaling machines, the importance of concentration and competition, the use of structural information, the integration of large-scale proteomic datasets, and the future challenge of integrating alternative splicing into signaling pathways.

Protein design in biological networks: from manipulating the input to modifying the output

Protein engineering has been an invaluable tool for the deciphering of protein folding and function and in the understanding of biological signaling networks. From an applied point of view it has been of paramount importance in biotechnological and biopharmaceutical products and applications. Traditionally, the protein engineering tools of choice were classical rational design, or directed evolution-based methods. In recent years, a third tool has matured: computational protein design (CPD). In this review, we summarize the underlying principles of CPD and discuss its application for understanding and modifying biological systems. Three main applications of the use of protein design will be highlighted and reviewed: artificially rewiring of signal transduction networks, prediction and generation of large-scale in silico interaction networks and using protein design to manipulate gene expression.

Improved Binding of Raf to Ras·GDP Is Correlated with Biological Activity

The GTP-binding protein Ras plays a central role in the regulation of various cellular processes, acting as a molecular switch that triggers signaling cascades. Only Ras bound to GTP is able to interact strongly with effector proteins like Raf kinase, phosphatidylinositol 3-kinase, and RalGDS, whereas in the GDP-bound state, the stability of the complex is strongly decreased, and signaling is interrupted. To determine whether this process is only controlled by the stability of the complex, we used computer-aided protein design to improve the interaction between Ras and effector. We challenged the Ras·Raf complex in this study because Raf among all effectors shows the highest Ras affinity and the fastest association kinetics. The proposed mutations were characterized as to their changes in dynamics and binding strength. We demonstrate that Ras-Raf interaction can only be improved at the cost of a loss in specificity of Ras·GTP versus Ras·GDP. As shown by NMR spectroscopy, the Raf mutation A85K leads to a shift of Ras switch I in the GTP-bound as well as in the GDP-bound state, thereby increasing the complex stability. In a luciferase-based reporter gene assay, Raf A85K is associated with higher signaling activity, which appears to be a mere matter of Ras-Raf affinity.

Association Rate Constants of Ras-Effector Interactions Are Evolutionarily Conserved

Evolutionary conservation of protein interaction properties has been shown to be a valuable indication for functional importance. Here we use homology interface modeling of 10 Ras-effector complexes by selecting ortholog proteins from 12 organisms representing the major eukaryotic branches, except plants. We find that with increasing divergence time the sequence similarity decreases with respect to the human protein, but the affinities and association rate constants are conserved as predicted by the protein design algorithm, FoldX. In parallel we have done computer simulations on a minimal network based on Ras-effector interactions, and our results indicate that in the absence of negative feedback, changes in kinetics that result in similar binding constants have strong consequences on network behavior. This, together with the previous results, suggests an important biological role, not only for equilibrium binding constants but also for kinetics in signaling processes involving Ras-effector interactions. Our findings are important to take into consideration in system biology approaches and simulations of biological networks.