Computational studies of the principle of dynamic-change-driven protein interactions

Molecular Dynamics Investigation into the Stability of KRas and CRaf Multimeric Complexes

In the Ras/Raf/MAPK signaling pathway, Ras and Raf proteins interact synergistically to form a tetrameric complex. NMR experiments have demonstrated that Ras dimerizes in solution and binds stably to Raf, forming Ras·Raf complexes. In this study, we constructed the ternary and quaternary complexes of KRas and CRaf based on crystal structures, denoted as (KRas)2·CRaf and (KRas)2·(CRaf)2, respectively. Molecular dynamics (MD) simulations were performed to investigate the stability of these complexes, while hydrogen bonds as well as salt bridges formed at the protein-protein interaction interfaces were analyzed based on simulation trajectories. The results revealed that the KRas·CRaf complex is more stable in explicit solvent compared with the KRas dimer. Formation of the stable quaternary complex (KRas)2·(CRaf)2 might be attributed to the association of two binary KRas·CRaf complexes. Additionally, MD simulations of the KRasG12D·CRaf complex revealed a stable and extended binding site at the KRas-CRaf interaction interface. This binding site was identified as a potential therapeutic target to block abnormal signal transmission in the pathway.

Kinetic and thermodynamic allostery in the Ras protein family

Allostery, the transfer of information between distant parts of a macromolecule, is a fundamental feature of protein function and regulation. However, allosteric mechanisms are usually not explained by protein structure, requiring information on correlated fluctuations uniquely accessible to molecular simulation. Existing work to extract allosteric pathways from molecular dynamics simulations has focused on thermodynamic correlations. Here, we show how kinetic correlations encode complementary information essential to explain observed variations in allosteric regulation. We applied kinetic and thermodynamic correlation analysis on atomistic simulations of H, K, and NRas isoforms in the apo, GTP, and GDP-bound states of Ras protein, with and without complexing to its downstream effector, Raf. We show that switch I and switch II are the primary components of thermodynamic and kinetic allosteric networks, consistent with the key roles of these two motifs. These networks connect the switches to an allosteric loop recently discovered from a crystal structure of HRas. This allosteric loop is inactive in KRas, but is coupled to the hydrolysis arm switch II in NRas and HRas. We find that the mechanism in the latter two isoforms are thermodynamic and kinetic, respectively. Binding of Raf-RBD further activates thermodynamic allostery in HRas and KRas but has limited effect on NRas. These results indicate that kinetic and thermodynamic correlations are both needed to explain protein function and allostery. These two distinct channels of allosteric regulation, and their combinatorial variability, may explain how subtle mutational differences can lead to diverse regulatory profiles among enzymatic proteins.

Allosteric site variants affect GTP hydrolysis on Ras

RAS GTPases are proto-oncoproteins that regulate cell growth, proliferation, and differentiation in response to extracellular signals. The signaling functions of RAS, and other small GTPases, are dependent on their ability to cycle between GDP-bound and GTP-bound states. Structural analyses suggest that GTP hydrolysis catalyzed by HRAS can be regulated by an allosteric site located between helices 3, 4, and loop 7. Here we explore the relationship between intrinsic GTP hydrolysis on HRAS and the position of helix 3 and loop 7 through manipulation of the allosteric site, showing that the two sites are functionally connected. We generated several hydrophobic mutations in the allosteric site of HRAS to promote shifts in helix 3 relative to helix 4. By combining crystallography and enzymology to study these mutants, we show that closure of the allosteric site correlates with increased hydrolysis of GTP on HRAS in solution. Interestingly, binding to the RAS binding domain of RAF kinase (RAF-RBD) inhibits GTP hydrolysis in the mutants. This behavior may be representative of a cluster of mutations found in human tumors, which potentially cooperate with RAF complex formation to stabilize the GTP-bound state of RAS.

A proteome-scale analysis of vertebrate protein amino acid occurrence: Thermoadaptation shows a correlation with protein solvation but less so with dynamics

Despite differences in behaviors and living conditions, vertebrate organisms share the great majority of proteins, often with subtle differences in amino acid sequence. Here, we present a simple way to analyze the difference in amino acid occurrence by comparing highly homologous proteins on a subproteome level between several vertebrate model organisms. Specifically, we use this method to identify a pattern of amino acid conservation as well as a shift in amino acid occurrence between homeotherms (warm-blooded species) and poikilotherms (cold-blooded species). Importantly, this general analysis and a specific example further establish a broad correlation, if not likely connection between the thermal adaptation of protein sequences and two of their physical features: on average a change in their protein dynamics and, even more strongly, in their solvation. For poikilotherms, such as frog and fish, the lower body temperature is expected to increase the protein-protein interaction due to a decrease in protein internal dynamics. In order to counteract the tendency for enhanced binding caused by low temperatures, poikilotherms enhance the solvation of their proteins by favoring polar amino acids. This feature appears to dominate over possible changes in dynamics for some proteins. The results suggest that a general trend for amino acid choice is part of the mechanism for thermoadaptation of vertebrate organisms at the molecular level.

Influence of the Interdomain Interface on Structural and Redox Properties of Multiheme Proteins

Multiheme proteins are important in energy conversion and biogeochemical cycles of nitrogen and sulfur. A diheme cytochrome c4 (c4) was used as a model to elucidate roles of the interdomain interface on properties of iron centers in its hemes A and B. Isolated monoheme domains c4-A and c4-B, together with the full-length diheme c4 and its Met-to-His ligand variants, were characterized by a variety of spectroscopic and stability measurements. In both isolated domains, the heme iron is Met/His-ligated at pH 5.0, as in the full-length c4, but becomes His/His-ligated in c4-B at higher pH. Intradomain contacts in c4-A are minimally affected by the separation of c4-A and c4-B domains, and isolated c4-A is folded. In contrast, the isolated c4-B is partially unfolded, and the interface with c4-A guides folding of this domain. The c4-A and c4-B domains have the propensity to interact even without the polypeptide linker. Thermodynamic cycles have revealed properties of monomeric folded isolated domains, suggesting that ferrous (FeII), but not ferric (FeIII) c4-A and c4-B, is stabilized by the interface. This study illustrates the effects of the interface on tuning structural and redox properties of multiheme proteins and enriches our understanding of redox-dependent complexation.

SAM1 domain of SASH1 harbors distinctive structural heterogeneity

The sterile alpha motif (SAM) domains are among the most versatile protein domains in biology, and the variety of the oligomerization states contribute to their diverse roles in many diseases. A better understanding of the structure and dynamics of various SAM domains will provide a scientific basis for drug development targeting them. Here, we used SEC-MALS, HPLC, NMR, and other biophysical techniques to characterize the structural features and dynamics of the SAM1 domain in SASH1. SASH1 is a scaffold protein belonging to the same family as SASH3. Unlike the dimerization seen in SASH3's SAM domain, our SEC-MALS and SE-HPLC showed that SAM1 exists primarily as a less compact monomer with a minor oligomer. NMR assignment, relaxation, and exchange experiments revealed the presence of both a disordered monomer and a more structured oligomer with multiple timescale exchange regimes in solution. Mutagenesis and SE-HPLC showed that D663A/T664K substitutions in SAM1 increased its oligomerization. In sum, this study is the first to characterize a disordered structure for a SAM domain, provides additional evidence and framework for the diversity of SAM domains, and identifies a region in SAM1 as a potential starting point to further characterize the structural mechanism of oligomerization of the domain.

Hunting for Novel Routes in Anticancer Drug Discovery: Peptides against Sam-Sam Interactions

Among the diverse protein binding modules, Sam (Sterile alpha motif) domains attract attention due to their versatility. They are present in different organisms and play many functions in physiological and pathological processes by binding multiple partners. The EphA2 receptor contains a Sam domain at the C-terminus (EphA2-Sam) that is able to engage protein regulators of receptor stability (including the lipid phosphatase Ship2 and the adaptor Odin). Ship2 and Odin are recruited by EphA2-Sam through heterotypic Sam-Sam interactions. Ship2 decreases EphA2 endocytosis and consequent degradation, producing chiefly pro-oncogenic outcomes in a cellular milieu. Odin, through its Sam domains, contributes to receptor stability by possibly interfering with ubiquitination. As EphA2 is upregulated in many types of tumors, peptide inhibitors of Sam-Sam interactions by hindering receptor stability could function as anticancer therapeutics. This review describes EphA2-Sam and its interactome from a structural and functional perspective. The diverse design strategies that have thus far been employed to obtain peptides targeting EphA2-mediated Sam-Sam interactions are summarized as well. The generated peptides represent good initial lead compounds, but surely many efforts need to be devoted in the close future to improve interaction affinities towards Sam domains and consequently validate their anticancer properties.

Microsecond MD Simulations of the Plexin-B1 RBD: 2. N-H Probability Densities and Conformational Entropy in Ligand-Free, Rac1-Bound, and Dimer RBD

Orientational probability densities, P_eq = exp(-u) (u, local potential), of bond-vectors in proteins provide information on structural flexibility. The related conformational entropy, S_k = -∫P_eq(ln P_eq)dΩ - ln ∫dΩ, provides the entropic contribution to the free energy of the physical/biological process studied. We have developed a new method for deriving P_eq and S_k from MD simulations, using the N-H bond as probe. Recently we used it to study the dimerization of the Rho GTPase binding domain of Plexin-B1 (RBD). Here we use it to study RBD binding to the small GTPase Rac1. In both cases 1 μs MD simulations have been employed. The RBD has the ubiquitin fold with four mostly long loops. L3 is associated with GTPase binding, L4 with RBD dimerization, L2 participates in interdomain interactions, and L1 has not been associated with function. We find that RBD-Rac1 binding renders L1, L3, and L4 more rigid and the turns β₂/α₁ and α₂/β₅ more flexible. By comparison, RBD dimerization renders L4 more rigid, and the α-helices, the β-strands, and L2 more flexible. The rigidity of L1 in RBDRAC is consistent with L1-L3 contacts seen in previous MD simulations. The analysis of the L3-loop reveals two states of distinct flexibility which we associate with involvement in slow conformational exchange processes differing in their rates. Overall, the N-H bonds make an unfavorable entropic contribution of (5.9 ± 0.9) kJ/mol to the free energy of RBD-Rac1 binding; they were found to make a favorably contribution of (-7.0 ± 0.7) kJ/mol to the free energy of RBD dimerization. In summary, the present study provides a new perspective on the impact of Rac1 binding and dimerization on the flexibility characteristics of the RBD. Further studies are stimulated by the results of this work.