scholarly journals Phosphorylation of interfacial phosphosite leads to increased binding of Rap-Raf complex

2021 ◽  
Author(s):  
Devanand T ◽  
Susmita Ghosh ◽  
Prasanna Venkatraman ◽  
Satyavani Vemparala

The effect of phosphorylation of a serine residue in the Rap protein, residing at the complex interface of Rap-Raf complex is studied using atomistic molecular dynamics simulations. As the phosphosite of interest (SER39) is buried at the interface of the Rap-Raf complex, phosphorylation of only Rap protein was simulated and then complexed with the RBD of Raf for further analysis of complex stability. Our simulations reveal that the phosophorylation increases the binding of complex through strong electrostatic interactions and changes the charge distribution of the interface significantly. This is manifested as an increase in stable salt-bridge interactions between the Rap and Raf of the complex. Network analysis clearly shows that the phosphorylation of SER39 reorganizes the community network to include the entire region of Raf chain, including, Raf L4 loop potentially affecting downstream signalling.

2019 ◽  
Vol 24 (9) ◽  
pp. 928-938 ◽  
Author(s):  
Luca Palazzolo ◽  
Chiara Paravicini ◽  
Tommaso Laurenzi ◽  
Sara Adobati ◽  
Simona Saporiti ◽  
...  

SLC6A14 (ATB0,+) is a sodium- and chloride-dependent neutral and dibasic amino acid transporter that regulates the distribution of amino acids across cell membranes. The transporter is overexpressed in many human cancers characterized by an increased demand for amino acids; as such, it was recently acknowledged as a novel target for cancer therapy. The knowledge on the molecular mechanism of SLC6A14 transport is still limited, but some elegant studies on related transporters report the involvement of the 12 transmembrane α-helices in the transport mechanism, and describe structural rearrangements mediated by electrostatic interactions with some pivotal gating residues. In the present work, we constructed a SLC6A14 model in outward-facing conformation via homology modeling and used molecular dynamics simulations to predict amino acid residues critical for substrate recognition and translocation. We docked the proteinogenic amino acids and other known substrates in the SLC6A14 binding site to study both gating regions and the exposed residues involved in transport. Interestingly, some of these residues correspond to those previously identified in other LeuT-fold transporters; however, we could also identify a novel relevant residue with such function. For the first time, by combined approaches of molecular docking and molecular dynamics simulations, we highlight the potential role of these residues in neutral amino acid transport. This novel information unravels new aspects of the human SLC6A14 structure–function relationship and may have important outcomes for cancer treatment through the design of novel inhibitors of SLC6A14-mediated transport.


2014 ◽  
Vol 169 ◽  
pp. 303-321 ◽  
Author(s):  
Ariane Allain ◽  
Isaure Chauvot de Beauchêne ◽  
Florent Langenfeld ◽  
Yann Guarracino ◽  
Elodie Laine ◽  
...  

Allostery is a universal phenomenon that couples the information induced by a local perturbation (effector) in a protein to spatially distant regulated sites. Such an event can be described in terms of a large scale transmission of information (communication) through a dynamic coupling between structurally rigid (minimally frustrated) and plastic (locally frustrated) clusters of residues. To elaborate a rational description of allosteric coupling, we propose an original approach – MOdular NETwork Analysis (MONETA) – based on the analysis of inter-residue dynamical correlations to localize the propagation of both structural and dynamical effects of a perturbation throughout a protein structure. MONETA uses inter-residue cross-correlations and commute times computed from molecular dynamics simulations and a topological description of a protein to build a modular network representation composed of clusters of residues (dynamic segments) linked together by chains of residues (communication pathways). MONETA provides a brand new direct and simple visualization of protein allosteric communication. A GEPHI module implemented in the MONETA package allows the generation of 2D graphs of the communication network. An interactive PyMOL plugin permits drawing of the communication pathways between chosen protein fragments or residues on a 3D representation. MONETA is a powerful tool for on-the-fly display of communication networks in proteins. We applied MONETA for the analysis of communication pathways (i) between the main regulatory fragments of receptors tyrosine kinases (RTKs), KIT and CSF-1R, in the native and mutated states and (ii) in proteins STAT5 (STAT5a and STAT5b) in the phosphorylated and the unphosphorylated forms. The description of the physical support for allosteric coupling by MONETA allowed a comparison of the mechanisms of (a) constitutive activation induced by equivalent mutations in two RTKs and (b) allosteric regulation in the activated and non-activated STAT5 proteins. Our theoretical prediction based on results obtained with MONETA was validated for KIT by in vitro experiments. MONETA is a versatile analytical and visualization tool entirely devoted to the understanding of the functioning/malfunctioning of allosteric regulation in proteins – a crucial basis to guide the discovery of next-generation allosteric drugs.


2018 ◽  
Author(s):  
Mustapha Carab Ahmed ◽  
Elena Papaleo ◽  
Kresten Lindorff-Larsen

AbstractSalt bridges form between pairs of ionisable residues in close proximity and are important interactions in proteins. While salt bridges are known to be important both for protein stability, recognition and regulation, we still do not have fully accurate predictive models to assess the energetic contributions of salt bridges. Molecular dynamics simulations is one technique that may be used study the complex relationship between structure, solvation and energetics of salt bridges, but the accuracy of such simulations depend on the force field used. We have used NMR data on the B1 domain of protein G (GB1) to benchmark molecular dynamics simulations. Using enhanced sampling simulations, we calculated the free energy of forming a salt bridge for three possible ionic interactions in GB1. The NMR experiments showed that these interactions are either not formed, or only very weakly formed, in solution. In contrast, we show that the stability of the salt bridges is slightly overestimated in simulations of GB1 using six commonly used combinations of force fields and water models. We therefore conclude that further work is needed to refine our ability to model quantitatively the stability of salt bridges through simulations, and that comparisons between experiments and simulations will play a crucial role in furthering our understanding of this important interaction.


2018 ◽  
Author(s):  
Annie M. Westerlund ◽  
Lucie Delemotte

AbstractCalmodulin (CaM) is a calcium sensing protein that regulates the function of a large number of proteins, thus playing a crucial part in many cell signaling path- ways. CaM has the ability to bind more than 300 different target peptides in a Ca2+-dependent manner, mainly through the exposure of hydrophobic residues. How CaM can bind a large number of targets while retaining some selectivity is a fascinating open question.Here, we explore the mechanism of CaM selective promiscuity for selected target proteins. Analyzing enhanced sampling molecular dynamics simulations of Ca2+-bound and Ca2+-free CaM via spectral clustering has allowed us to identify distinct conformational states, characterized by interhelical angles, secondary structure determinants and the solvent exposure of specific residues. We searched for indicators of conformational selection by mapping solvent exposure of residues in these conformational states to contacts in structures of CaM/target peptide complexes. We thereby identified CaM states involved in various binding classes arranged along a depth binding gradient. Binding Ca2+ modifies the accessible hydrophobic surface of the two lobes and allows for deeper binding. Apo CaM indeed shows shallow binding involving predominantly polar and charged residues. Furthermore, binding to the C-terminal lobe of CaM appears selective and involves specific conformational states that can facilitate deep binding to target proteins, while binding to the N-terminal lobe appears to happen through a more flexible mechanism. Thus the long-ranged electrostatic interactions of the charged residues of the N-terminal lobe of CaM may initiate binding, while the short-ranged interactions of hydrophobic residues in the C-terminal lobe of CaM may account for selectivity.This work furthers our understanding of the mechanism of CaM binding and selectivity to different target proteins and paves the way towards a comprehensive model of CaM selectivity.Author summaryCalmodulin is a protein involved in the regulation of a variety of cell signaling pathways. It acts by making usually calcium-insensitive proteins sensitive to changes in the calcium concentration inside the cell. Its two lobes bind calcium and allow the energetically unfavorable exposure of hydrophobic residues to the aqueous environment which can then bind target proteins. The mechanisms behind the simultaneous specificity and variation of target protein binding is yet unknown but will aid understanding of the calcium-signaling and regulation that occur in many of our cellular processes.Here, we used molecular dynamics simulations and data analysis techniques to investigate what effect calcium has on the binding modes of calmodulin. The simulations and analyses allow us to observe and differentiate specific states. One domain of calmodulin is shown to be selective with binding involving short- distance interactions between hydrophobic residues, while the other binds target proteins through a more flexible mechanism involving long-distance electrostatic interactions.


2020 ◽  
Vol 22 (5) ◽  
pp. 2887-2907
Author(s):  
Trevor R. Fisher ◽  
Guobing Zhou ◽  
Yijun Shi ◽  
Liangliang Huang

Molecular dynamics simulations reveal that the maximum contribution of H-bonds between water and glycerol occurs around 30 mol% glycerol. Such a concentration is also where several of the mixture's properties have an observed maxima or minima.


Sign in / Sign up

Export Citation Format

Share Document