Interplay between structural rigidity and electrostatic interactions in the ligand binding domain of GluR2

Using molecular dynamics (MD) simulations, computational protein modifications, and a novel theoretical methodology that determines structural rigidity/flexibility (the FIRST algorithm), we investigate how molecular structure and dynamics of the glutamate receptor ligand binding domain (GluR2 S1S2) facilitate its conformational transition. S1S2 is a two-lobe protein, which undergoes a cleft closure conformational transition upon binding an agonist in the cleft between the two lobes; hence it is expected that the mechanism of this conformational transition can be characterized as a hinge-type. However, in the rigidity analysis one lobe of the protein is identified as a single rigid cluster while the other one is structurally flexible, inconsistent with a presumed mechanical hinge mechanism. Instead, we characterize the cleft-closing transition as a load and lock mechanism. We find that when two cross-cleft hydrogen bonds are disrupted the protein undergoes a rapid cleft opening transition. At the same time, the dynamical behavior of the cleft in the presence of the glutamate ligand is only weakly affected by the S652 peptide bond in its flipped conformation observed in the crystal structure. The residue E705 plays significant role in stabilization of the closed conformation via electrostatic interactions. The presence of the E705-K730 salt bridge seems to correlate strongly withthe cleft opening transition in the MD simulations.     

Flexibility of prolyl oligopeptidase: molecular dynamics and molecular framework analysis of the potential substrate pathways

The flexibility of prolyl oligopeptidase has been investigated using molecular dynamics (MD) and molecular framework approaches to delineate the route of the substrate to the active site. The selectivity of the enzyme is mediated by a seven-bladed beta-propeller that in the crystal structure does not indicate the possible passage for the substrate to the catalytic center. Its open topology however, could allow the blades to move apart and let the substrate into the large central cavity. Flexibility analysis of prolyl oligopeptidase structure using the FIRST (Floppy Inclusion and Rigid Substructure Topology) approach and the atomic fluctuations derived from MD simulations demonstrated the rigidity of the propeller domain, which does not permit the substrate to approach the active site through this domain. Instead, a smaller tunnel at the inter-domain region comprising the highly flexible N-terminal segment of the peptidase domain and a facing hydrophilic loop from the propeller (residues 192-205) was identified by cross-correlation analysis and essential dynamics as the only potential pathway for the substrate. The functional importance of the flexible loop has been also verified by kinetic analysis of the enzyme with a split loop. Catalytic effect of engineered disulfide bridges was rationalized by characterizing the concerted motions of the two domains.     

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.     

Molecular dynamics simulations of the Ras:Raf and Rap:Raf complexes

The protein Raf is an immediate downstream target of Ras in the MAP kinase signalling pathway. The complex of Ras with the Ras-binding domain (RBD) of Raf has been modelled by homology to the (E30D,K31E)-Rap1A:RBD complex, and both have been subjected to multiple molecular dynamics simulations in solution. While both complexes are stable, several rearrangements occur in the Ras:RBD simulations: the RBD loop 100-109 moves closer to Ras, Arg73 in the RBD moves towards Ras to form a salt bridge with Ras-Asp33, and Loop 4 of the Ras switch II region shifts upwards toward the RBD. The Ras:RBD interactions (including the RBD-Arg73 interaction) are consistent with available NMR and mutagenesis data on the Ras: RBD complex in solution. The Ras switch II region does not interact directly with the RBD, although indirect interactions exist through the effector domain and bridging water molecules. No large-scale RBD motion is seen in the Ras:RBD complex, compared to the Rap:RBD complex, to suggest an allosteric activation of Raf by Ras. This may be because the Raf kinase domain (whose structure is unknown) is not included in the model.     

What Makes Ras an Efficient Molecular Switch: A Computational, Biophysical, and Structural Study of Ras-GDP Interactions with Mutants of Raf

Ras is a small GTP-binding protein that is an essential molecular switch for a wide variety of signaling pathways including the control of cell proliferation, cell cycle progression and apoptosis. In the GTP-bound state, Ras can interact with its effectors, triggering various signaling cascades in the cell. In the GDP-bound state, Ras looses its ability to bind to known effectors. The interaction of the GTP-bound Ras (Ras^GTP) with its effectors has been studied intensively. However, very little is known about the much weaker interaction between the GDP-bound Ras (Ras^GDP) and Ras effectors. We investigated the factors underlying the nucleotide-dependent differences in Ras interactions with one of its effectors, Raf kinase. Using computational protein design, we generated mutants of the Ras-binding domain of Raf kinase (Raf) that stabilize the complex with Ras^GDP. Most of our designed mutations narrow the gap between the affinity of Raf for Ras^GTP and Ras^GDP, producing the desired shift in binding specificity towards Ras^GDP. A combination of our best designed mutation, N71R, with another mutation, A85K, yielded a Raf mutant with a 100-fold improvement in affinity towards Ras^GDP. The Raf A85K and Raf N71R/A85K mutants were used to obtain the first high-resolution structures of Ras^GDP bound to its effector. Surprisingly, these structures reveal that the loop on Ras previously termed the switch I region in the Ras^GDP·Raf mutant complex is found in a conformation similar to that of Ras^GTP and not Ras^GDP. Moreover, the structures indicate an increased mobility of the switch I region. This greater flexibility compared to the same loop in Ras^GTP is likely to explain the natural low affinity of Raf and other Ras effectors to Ras^GDP. Our findings demonstrate that an accurate balance between a rigid, high-affinity conformation and conformational flexibility is required to create an efficient and stringent molecular switch.     

The Scaffold Protein Shoc2/SUR-8 Accelerates the Interaction of Ras and Raf

Shoc2/SUR-8 positively regulates Ras/ERK MAP kinase signaling by serving as a scaffold for Ras and Raf. Here, we examined the role of Shoc2 in the spatio-temporal regulation of Ras by using a fluorescence resonance energy transfer (FRET)-based biosensor, together with computational modeling. In epidermal growth factor-stimulated HeLa cells, RNA-mediated Shoc2 knockdown reduced the phosphorylation of MEK and ERK with half-maximal inhibition, but not the activation of Ras. For the live monitoring of Ras binding to Raf, we utilized a FRET biosensor wherein Ras and the Ras-binding domain of Raf were connected tandemly and sandwiched with acceptor and donor fluorescent proteins for the FRET measurement. With this biosensor, we found that Shoc2 was required for the rapid interaction of Ras with Raf upon epidermal growth factor stimulation. To decipher the molecular mechanisms underlying the kinetics, we developed two computational models that might account for the action of Shoc2 in the Ras-ERK signaling. One of these models, the Shoc2 accelerator model, provided a reasonable explanation of the experimental observations. In this Shoc2 accelerator model, Shoc2 accelerated both the association and dissociation of Ras-Raf interaction. We propose that Shoc2 regulates the spatio-temporal patterns of the Ras-ERK signaling pathway primarily by accelerating the Ras-Raf interaction.     

Comparative molecular dynamics of mesophilic and psychrophilic protein homologues studied by 1.2 ns simulations

It is well established that the dynamic motion of proteins plays an important functional role, and that the adaptation of a protein molecule to its environment requires optimization of internal non-covalent interactions and protein-solvent interactions. Serine proteinases in general, and trypsin in particular has been used as a model system in exploring possible structural features for cold adaptation. In this study, a 500 p.s. and a 1200 p.s. molecular dynamics (MD) simulation at 300 K of both anionic salmon trypsin and cationic bovine trypsin are analyzed in terms of molecular flexibility, internal non-covalent interactions and protein-solvent interactions. The present MD simulations do not indicate any increased flexibility of the cold adapted enzyme on an overall basis. However, the apparent higher flexibility and deformability of the active site of anionic salmon trypsin may lower the activation energy for ligand binding and for catalysis, and might be a reason for the increased binding affinity and catalytic efficiency compared to cationic bovine trypsin.     

Identification of putative, stable binding regions through flexibility analysis of HIV-1 gp120

The acquired-immunodeficiency syndrome has evolved into a major worldwide epidemic. Significant effort has been made in the development of antiviral therapies. A new strategy for vaccine and drug design that complements the existing cocktail therapy is to target entry of the human immunodeficiency virus (HIV). Such an approach provides the advantage of interfering with multiple intermediates in this multi-step process. The extraordinary conformational flexibility, glycosylation, and strain variations of viral glycoprotein gp120 cause general viral evasion of humoral immune response and thus complicate the development of an effective vaccine. Especially difficult to define are the conformation of gp120 before CD4 engagement as well as the relative orientations of the V1/V2 and V3 loops with respect to the inner and outer domains. In this study, we used Floppy Inclusion and Rigid Substructure Topography (FIRST), a program based on graph theory, to analyze the flexibility and rigidity of all known HIV-1 gp120 structures. A flexibility index is used to describe and compare the spatial distribution of protein flexibility and rigidity of these structures in isolation and in complex with CD4, CD4-mimics, and neutralizing antibodies. Using this flexibility analysis, we identified a universal rigid region (the alpha2 helix) as well as the consensus largest rigid cluster involving a beta-sheet located on the coreceptor binding face. Both of these regions may serve as stable targets for vaccine design and drug discovery. Detailed comparisons of the changes in flexibility based on strain variations, stabilizing mutations, binding features of CD4 mimics, and impact of b12 binding are reported.     

Analyzing the flexibility of RNA structures by constraint counting

RNA requires conformational dynamics to undergo its diverse functional roles. Here, a new topological network representation of RNA structures is presented that allows analyzing RNA flexibility/rigidity based on constraint counting. The method extends the FIRST approach, which identifies flexible and rigid regions in atomic detail in a single, static, three-dimensional molecular framework. Initially, the network rigidity of a canonical A-form RNA is analyzed by counting on constraints of network elements of increasing size. These considerations demonstrate that it is the inclusion of hydrophobic contacts into the RNA topological network that is crucial for an accurate flexibility prediction. The counting also explains why a protein-based parameterization results in overly rigid RNA structures. The new network representation is then validated on a tRNA^ASP structure and all NMR-derived ensembles of RNA structures currently available in the Protein Data Bank (with chain length ≥40). The flexibility predictions demonstrate good agreement with experimental mobility data, and the results are superior compared to predictions based on two previously used network representations. Encouragingly, this holds for flexibility predictions as well as mobility predictions obtained by constrained geometric simulations on these networks. Potential applications of the approach to analyzing the flexibility of DNA and RNA/protein complexes are discussed.     

Modeling correlated main-chain motions in proteins for flexible molecular recognition

We describe a new method for modeling protein and ligand main-chain flexibility, and show its ability to model flexible molecular recognition. The goal is to sample the full conformational space, including large-scale motions that typically cannot be reached in molecular dynamics simulations due to the computational intensity, as well as conformations that have not been observed yet by crystallography or NMR. A secondary goal is to assess the degree of flexibility consistent with protein-ligand recognition. Flexibility analysis of the target protein is performed using the graph-theoretic algorithm FIRST, which also identifies coupled networks of covalent and noncovalent bonds within the protein. The available conformations of the flexible regions are then explored with ROCK by random-walk sampling of the rotatable bonds. ROCK explores correlated motions by only sampling dihedral angles that preserve the coupled bond networks in the protein and generates conformers with good stereochemistry, without using a computationally expensive potential function. A representative set of the conformational ensemble generated this way can be used as targets for docking with SLIDE, which handles the flexibility of protein and ligand side-chains. The realism of this protein main-chain conformational sampling is assessed by comparison with time-resolved NMR studies of cyclophilin A motions. ROCK is also effective for modeling the flexibility of large cyclic and polycyclic ligands, as demonstrated for cyclosporin and zearalenol. The use of this combined approach to perform docking with main-chain flexibility is illustrated for the cyclophilin A-cyclosporin complex and the estrogen receptor in complex with zearalenol, while addressing the question of how much flexibility is allowed without hindering molecular recognition.     

Critical binding and regulatory interactions between Ras and Raf occur through a small, stable N-terminal domain of Raf and specific Ras effector residues.

Genetic and biochemical evidence suggests that the Ras protooncogene product regulates the activation of the Raf kinase pathway, leading to the proposal that Raf is a direct mitogenic effector of activated Ras. Here we report the use of a novel competition assay to measure in vitro the relative affinity of the c-Raf-1 regulatory region for Ras-GTP, Ras-GDP, and 10 oncogenic and effector mutant Ras proteins. c-Raf-1 associates with normal Ras and the oncogenic V12 and L61 forms of Ras with equal affinity. The moderately transforming mutant Ras[E30K31] also bound to the c-Raf-1 regulatory region with normal affinity. Transformation-defective Ras effector mutants Ras[N33], Ras[S35], and Ras[N38] bound poorly. In contrast, the transformation defective Ras[G26I27] and Ras[E45] mutants bound to the c-Raf-1 regulatory region with nearly wild-type affinity. A stable, high-affinity Ras-binding region of c-Raf-1 was mapped to a 99-amino-acid subfragment of the first 257 residues. The smallest Ras-binding region identified consisted of N-terminal residues 51 to 131, although stable expression of the domain and high-affinity binding were improved by the presence of residues 132 to 149. Deletion of the Raf zinc finger region did not reduce Ras-binding affinity, while removal of the first 50 amino acids greatly increased affinity. Phosphorylation of Raf[1-149] by protein kinase A on serine 43 resulted in significant inhibiton of Ras binding. demonstrating that the mechanism of cyclic AMP downregulation results through structural changes occurring exclusively in this small Ras-binding domain.     

Protein flexibility predictions using graph theory

Techniques from graph theory are applied to analyze the bond networks in proteins and identify the flexible and rigid regions. The bond network consists of distance constraints defined by the covalent and hydrogen bonds and salt bridges in the protein, identified by geometric and energetic criteria. We use an algorithm that counts the degrees of freedom within this constraint network and that identifies all the rigid and flexible substructures in the protein, including overconstrained regions (with more crosslinking bonds than are needed to rigidify the region) and underconstrained or flexible regions, in which dihedral bond rotations can occur. The number of extra constraints or remaining degrees of bond-rotational freedom within a substructure quantifies its relative rigidity/flexibility and provides a flexibility index for each bond in the structure. This novel computational procedure, first used in the analysis of glassy materials, is approximately a million times faster than molecular dynamics simulations and captures the essential conformational flexibility of the protein main and side-chains from analysis of a single, static three-dimensional structure. This approach is demonstrated by comparison with experimental measures of flexibility for three proteins in which hinge and loop motion are essential for biological function: HIV protease, adenylate kinase, and dihydrofolate reductase.     

Association of yeast adenylyl cyclase with cyclase-associated protein CAP forms a second Ras-binding site which mediates its Ras-dependent activation

Posttranslational modification, in particular farnesylation, of Ras is crucial for activation of Saccharomyces cerevisiae adenylyl cyclase (CYR1). Based on the previous observation that association of CYR1 with cyclase-associated protein (CAP) is essential for its activation by posttranslationally modified Ras, we postulated that the associated CAP might contribute to the formation of a Ras-binding site of CYR1, which mediates CYR1 activation, other than the primary Ras-binding site, the leucine-rich repeat domain. Here, we observed a posttranslational modification-dependent association of Ras with a complex between CAP and CYR1 C-terminal region. When CAP mutants defective in Ras signaling but retaining the CYR1-binding activity were isolated by screening of a pool of randomly mutagenized CAP, CYR1 complexed with two of the obtained three mutants failed to be activated efficiently by modified Ras and exhibited a severely impaired ability to bind Ras, providing a genetic evidence for the importance of the physical association with Ras at the second Ras-binding site. On the other hand, CYR1, complexed with the other CAP mutant, failed to be activated by Ras but exhibited a greatly enhanced binding to Ras. Conversely, a Ras mutant E31K, which exhibits a greatly enhanced binding to the CYR1-CAP complex, failed to activate CYR1 efficiently. Thus, the strength of interaction at the second Ras-binding site appears to be a critical determinant of CYR1 regulation by Ras: too-weak and too-strong interactions are both detrimental to CYR1 activation. These results, taken together with those obtained with mammalian Raf, suggest the importance of the second Ras-binding site in effector regulation.     

Long time molecular dynamic simulation on the agonist binding and activation of the β2-adrenergic receptor

β₂-Adrenergic receptor (β₂AR) plays a critical role in mediating the effects of catecholamine hormones. Due to the flexibility of the structure of its active state, study of agonist–β2AR is usually performed by molecular dynamic (MD) simulation. In this study we show the representative characteristics of agonist binding and activation on β₂AR by MD simulation. The binding site and the conformational changes in the specific regions of β2AR are reasonable which confirmed the conclusion that agonist–β2AR reached its active-like state. We have studied the disruption of non-covalent intramolecular interactions, including the conserved DRY motif, the rotamer toggle switch and the ionic lock, the cytoplasmic ends of transmembranes 5 and 6, and some water-mediated hydrogen bond network regions. We conclude that agonist induced β₂AR to its active conformation, or at least the active-like conformation.     

Comprehensive analysis of motions in molecular dynamics trajectories of the actin capping protein and its inhibitor complexes

The actin capping protein (CP) binds to actin filaments to block further elongation. The capping activity is inhibited by proteins V‐1 and CARMIL interacting with CP via steric and allosteric mechanisms, respectively. The crystal structures of free CP, CP/V‐1, and CP/CARMIL complexes suggest that the binding of CARMIL alters the flexibility of CP rather than the overall structure of CP, and this is an allosteric inhibition mechanism. Here, we performed molecular dynamics (MD) simulations of CP in the free form, and in complex with CARMIL or V‐1. The resulting trajectories were analyzed exhaustively using Motion Tree, which identifies various rigid‐body motions ranging from small local motions to large domain motions. After enumerating all the motions, CP flexibilities with different ligands were characterized by a list of frequencies for 20 dominant rigid‐body motions, some of which were not identified in previous studies. The comparative analysis highlights the influence of the binding of the CARMIL peptide to CP flexibility. In free CP and the CP/V‐1 complex, domain motions around a large crevice between the N‐stalk and the CP‐S domain occur frequently. The CARMIL peptide binds the crevice and suppresses the motions effectively. In addition, the binding of the CARMIL peptide enhances and alters local motions around the pocket that participates in V‐1 binding. These newly identified motions are likely to suppress the binding of V‐1 to CP. The observed changes in CP motion provide insights that describe the mechanism of allosteric regulation by CARMIL through modulating CP flexibility. Proteins 2016; 84:948–956. © 2016 Wiley Periodicals, Inc.     

Molecular trigger for pre-transfer editing pathway in Valyl-tRNA synthetase: A molecular dynamics simulation study

Pre-transfer editing pathway in Valyl-tRNA synthetase (ValRS) is a very important process to maintain the high fidelity of protein synthesis. However, molecular basis for this pathway remains unclear. Here we employed molecular dynamics (MD) simulation to study two complexes, ValRS·tRNA^val·Val-AMP (complex V) and ValRS·tRNA^val·Thr-AMP (complex T), and compared their simulation trajectories, in order to understand how the pre-transfer editing pathway is triggered by the noncognate substrate Thr-AMP. The MD simulations showed that the binding of Thr-AMP to ValRS led to different motions from those in complex V: clockwise rotation of the editing domain along the hinge region, and strong motions in the catalytic domain, especially in KMSKS loop. We found that the changed motion of Trp495 induced by Thr-AMP serves as a signal to discriminate Thr-AMP from Val-AMP, and the rigid ⁴⁹¹ILFL⁴⁹⁴ segment then propagates this signal from Trp495 to Asp490 and induces dissociation of the salt-bridge Asp490-Arg346 and formation of the salt-bridge Glu189-Lys533. The change in salt-bridges alters the motion of KMSKS loop and the editing domain, and eventually triggers the pre-transfer editing pathway. This study provides a model for the molecular trigger of the pre-transfer editing pathway in ValRS, and is useful for further exploring this process.     

Modeling the intact form of the alpha 1-proteinase inhibitor

The structure of the intact form of the serpin alpha 1-proteinase inhibitor has been modeled based on the assumption that the central strand s4A of the six-stranded beta-sheet A of the cleaved inhibitor is not incorporated into the sheet of intact alpha 1-proteinase inhibitor. This strand was removed from its position in the center of the sheet by suitable rotations about the backbone dihedrals of Lys343 using molecular graphics. The resulting structure was then annealed using molecular dynamics (MD) while applying progressive distance restraints to the reactive peptide bond (Met358-Ser359) for 50 ps. During this time, the disrupted beta-sheet reformed to create a five-stranded beta-sheet with strands 3 and 5 in a parallel arrangement. This change and accompanying structural rearrangements are largely confirmed by the X-ray structure of plakalbumin, whose structure reflects the overall structure of intact serpins. The successful modeling experiment demonstrates the utility of MD for making gross structural predictions based on related structures. The binding loop of the intact form is modeled to allow docking with serine proteinases, in particular thrombin, which most highly constrains the possible conformations of the binding loop.     

Stability of the beta-sheet of the WW domain: A molecular dynamics simulation study.

The WW domain consists of approximately 40 residues, has no disulfide bridges, and forms a three-stranded antiparallel beta-sheet that is monomeric in solution. It thus provides a model system for studying beta-sheet stability in native proteins. We performed molecular dynamics simulations of two WW domains, YAP65 and FBP28, with very different stability characteristics, in order to explore the initial unfolding of the beta-sheet. The less stable YAP domain is much more sensitive to simulation conditions than the FBP domain. Under standard simulation conditions in water (with or without charge-balancing counterions) at 300 K, the beta-sheet of the YAP WW domain disintegrated at early stages of the simulations. Disintegration commenced with the breakage of a hydrogen bond between the second and third strands of the beta-sheet due to an anticorrelated transition of the Tyr-28 psi and Phe-29 phi angles. Electrostatic interactions play a role in this event, and the YAP WW domain structure is more stable when simulated with a complete explicit model of the surrounding ionic strength. Other factors affecting stability of the beta-sheet are side-chain packing, the conformational entropy of the flexible chain termini, and the binding of cognate peptide.     

Homology modeling and molecular dynamics simulations of lymphotactin

We have modeled the structure of human lymphotactin (hLpnt), by homology modeling and molecular dynamics simulations. This chemokine is unique in having a single disulfide bond and a long C-terminal tail. Because other structural classes of chemokines have two pairs of Cys residues, compared to one in Lpnt, and because it has been shown that both disulfide bonds are required for stability and function, the question arises how the Lpnt maintains its structural integrity. The initial structure of hLpnt was constructed by homology modeling. The first 63 residues in the monomer of hLpnt were modeled using the structure of the human CC chemokine, RANTES, whose sequence appeared most similar. The structure of the long C-terminal tail, missing in RANTES, was taken from the human muscle fatty-acid binding protein. In a Protein Data Bank search, this protein was found to contain a sequence that was most homologous to the long tail. Consequently, the modeled hLpnt C-terminal tail consisted of both alpha-helical and beta-motifs. The complete model of the hLpnt monomer consisted of two alpha-helices located above the five-stranded beta-sheet. Molecular dynamics simulations of the solvated initial model have indicated that the stability of the predicted fold is related to the geometry of Pro78. The five-stranded beta-sheet appeared to be preserved only when Pro78 was modeled in the cis conformation. Simulations were also performed both for the C-terminal truncated forms of the hLpnt that contained one or two (CC chemokine-like) disulfide bonds, and for the chicken Lpnt (cLpnt). Our MD simulations indicated that the turn region (T30-G34) in hLpnt is important for the interactions with the receptor, and that the long C-terminal region stabilizes both the turn (T30-G34) and the five-stranded beta-sheet. The major conclusion from our theoretical studies is that the lack of one disulfide bond and the extension of the C-terminus in hLptn are mutually complementary. It is very likely that removal of two Cys residues sufficiently destabilizes the structure of a chemokine molecule, particularly the core beta-sheet, to abolish its biological function. However, this situation is rectified by the long C-terminal segment. The role of this long region is most likely to stabilize the first beta-turn region and alpha-helix H1, explaining how this chemokine can function with a single disulfide bond.