Structural consequences of site-directed mutagenesis in flexible protein domains: NMR characterization of the L(55,56)S mutant of RhoGDI.

The guanine dissociation inhibitor RhoGDI consists of a folded C-terminal domain and a highly flexible N-terminal region, both of which are essential for biological activity, that is, inhibition of GDP dissociation from Rho GTPases, and regulation of their partitioning between membrane and cytosol. It was shown previously that the double mutation L55S/L56S in the flexible region of RhoGDI drastically decreases its affinity for Rac1. In the present work we study the effect of this double mutation on the conformational and dynamic properties of RhoGDI, and describe the weak interaction of the mutant with Rac1 using chemical shift mapping. We show that the helical content of the region 45-56 of RhoGDI is greatly reduced upon mutation, thus increasing the entropic penalty for the immobilization of the helix, and contributing to the loss of binding. In contrast to wild-type RhoGDI, no interaction with Rac1 could be identified for amino-acid residues of the flexible domain of the mutant RhoGDI and only very weak binding was observed for the folded domain of the mutant. The origins of the effect of the L55S/L56S mutation on the binding constant (decreased by at least three orders of magnitude relative to wild-type) are discussed with particular reference to the flexibility of this part of the protein.

Structure-activity relationships in flexible protein domains: regulation of rho GTPases by RhoGDI and D4 GDI.

The guanine dissociation inhibitors RhoGDI and D4GDI inhibit guanosine 5'-diphosphate dissociation from Rho GTPases, keeping these small GTPases in an inactive state. The GDIs are made up of two domains: a flexible N-terminal domain of about 70 amino acid residues and a folded 134-residue C-terminal domain. Here, we characterize the conformation of the N-terminal regions of both RhoGDI and D4GDI using a series of NMR experiments which include (15)N relaxation and amide solvent accessibility measurements. In each protein, two regions with tendencies to form helices are identified: residues 36 to 58 and 9 to 20 in RhoGDI, and residues 36 to 57 and 20 to 25 in D4GDI. To examine the functional roles of the N-terminal domain of RhoGDI, in vitro and in vivo functional assays have been carried out with N-terminally truncated proteins. These studies show that the first 30 amino acid residues are not required for inhibition of GDP dissociation but appear to be important for GTP hydrolysis, whilst removal of the first 41 residues completely abolish the ability of RhoGDI to inhibit GDP dissociation. The combination of structural and functional studies allows us to explain why RhoGDI and D4GDI are able to interact in similar ways with the guanosine 5'-diphosphate-bound GTPase, but differ in their ability to regulate GTP-bound forms; these functional differences are attributed to the conformational differences of the N-terminal domains of the guanosine 5'-diphosphate dissociation inhibitors. Therefore, the two transient helices, appear to be associated with different biological effects of RhoGDI, providing a clear example of structure-activity relationships in a flexible protein domain.

Mapping the binding site for the GTP-binding protein Rac-1 on its inhibitor RhoGDI-1.

BACKGROUND: Members of the Rho family of small GTP-binding proteins, such as Rho, Rac and Cdc42, have a role in a wide range of cell responses. These proteins function as molecular switches by virtue of a conformational change between the GTP-bound (active) and GDP-bound (inactive) forms. In addition, most members of the Rho and Rac subfamilies cycle between the cytosol and membrane. The cytosolic guanine nucleotide dissociation inhibitors, RhoGDIs, regulate both the GDP/GTP exchange cycle and the membrane association/dissociation cycle. RESULTS: We have used NMR spectroscopy and site-directed mutagenesis to identify the regions of human RhoGDI-1 that are involved in binding Rac-1. The results emphasise the importance of the flexible regions of both proteins in the interaction. At least one specific region (residues 46-57) of the flexible N-terminal domain of RhoGDI, which has a tendency to form an amphipathic helix in the free protein, makes a major contribution to the binding energy of the complex. In addition, the primary site of Rac-1 binding on the folded domain of RhoGDI involves the beta4-beta5 and beta6-beta7 loops, with a slight movement of the 3(10) helix accompanying the interaction. This binding site is on the same face of the protein as the binding site for the isoprenyl group of post-translationally modified Rac-1, but is distinct from this site. CONCLUSIONS: Isoprenylated Rac-1 appears to interact with three distinct sites on RhoGDI. The isoprenyl group attached to the C terminus of Rac-1 binds in a pocket in the folded domain of RhoGDI. This is distinct from the major site on this domain occupied by Rac-1 itself, which involves two loops at the opposite end to the isoprenyl-binding site. It is probable that the flexible C-terminal region of Rac-1 extends from the site at which Rac-1 contacts the folded domain of RhoGDI to allow the isoprenyl group to bind in the pocket at the other end of the RhoGDI molecule. Finally, the flexible N terminus of RhoGDI-1, and particularly residues 48-58, makes a specific interaction with Rac-1 which contributes substantially to the binding affinity.

Stabilization of proteins by enhancement of inter-residue hydrophobic contacts: lessons of T4 lysozyme and barnase.

Although the hydrophobic interactions are considered as the main contributors to the protein stability, not much examples of protein stabilization by rational increasing of this type of interactions still can be found in literature. This is partly due to the lack of proper theoretical "measure" of hydrophobic interactions and their changes upon mutations. In the present paper the molecular hydrophobicity potential approach is used to assess how the changes in type and the strength of inter-residue contacts upon single amino acid mutations are correlated with the changes in thermodynamic stability of T4 lysozyme and barnase mutants, and which factors affect these correlations. Mutations changing unfavorable hydrophilic-to-hydrophobic contacts into favorable hydrophobic were found to enhance the thermodynamic stability in more than 81 % of cases, if these mutations do not create steric bumps and do not involve proline residues and hydrogen-bonded side-chains. Mutations increasing hydrophobic contributions (according to molecular hydrophobicity potential formalism) lead to increase of thermodynamic stability in more than 94% of cases for certain type of mutations (i.e., mutations not involving charged residues, Pro and residues with side-chain hydrogen bonds, when these mutations do not introduce steric bumps and do not involve strongly exposed residues and residues situated at helix N- and C-cap positions). For this type of mutations the correlation was found between the change in hydrophobic contributions of mutated residues deltaCphob and thermodynamic parameters deltaTm (change in melting temperature) and deltadeltaG (change in free energy of unfolding). Although the correlation coefficients were larger if the experimental structures of mutants were used for the calculations (correlation coefficients r(exp) deltaC,deltaT = .85 and r(exp) deltaC,deltadeltaG = 0.87) than if the modeled structures were used instead (r(mod) deltaC,deltaT = 0.74 and r(mod)deltaC,deltadeltaG = 0.76), the modelled structures of mutants in the vast majority of cases can be used for qualitative predictition of the protein stabilization. Basing on the analysis of mutations increasing hydrophobic contributions in T4 lysozyme the substitution matrix was derived, which can be used to decide which new residue should be put instead the old one to increase the stability of protein. The estimation shows that the number of potential mutation sites for enhancement of hydrophobic interactions in T4 lysozyme is quite large, and only approximately 10 per cent of them were studied thus far. Basing on the current analysis of T4 lysozyme and barnase mutations the algorithm for increasing of protein stability via increasing of hydrophobic interactions for the proteins with known spatial structure is proposed.

Recognizing misfolded and distorted protein structures by the assumption-based similarity score.

A new similarity score (sigma-score) is proposed which is able to find the correct protein structure among the very close alternatives and to distinguish between correct and deliberately misfolded structures. This score is based on the general principle 'similar likes similar', and it favors hydrophobic and hydrophilic contacts, and disfavors hydrophobic-to-hydrophilic contacts in proteins. The values of sigma-scores calculated for the high-resolution protein structures from the representative set are compared with those of alternatives: (i) very close alternatives which are only slightly distorted by conformational energy minimization in vacuo; (ii) alternatives with subsequently growing distortions, generated by molecular dynamics simulations in vacuo; (iii) structures derived by molecular dynamics simulation in solvent at 300 K; (iv) deliberately misfolded protein models. In nearly all tested cases the similarity score can successfully distinguish between experimental structure and its alternatives, even if the root mean square displacement of all heavy atoms is less than 1 A. The confidence interval of the similarity score was estimated using the high-resolution X-ray structures of domain pairs related by non-crystallographic symmetry. The similarity score can be used for the evaluation of the general quality of the protein models, choosing the correct structures among the very close alternatives, characterization of models simulating folding/unfolding, etc.

A new method to characterize hydrophobic organization of proteins: application to rational protein engineering of barnase.

We present a new algorithm for characterization of protein spatial structure basing on the molecular hydrophobicity potential approach. The method is illustrated by the analysis of three-dimensional structure of barnase and barnase-barstar complex. Current approach enables identification of amino acid residues situated in unfavorable environment (these residues may be "active" for binding), and to map quantitatively hydrophobic, hydrophilic and unfavorable hydrophobic-hydrophilic intra- and inter-molecular contacts involving backbone and side-chain segments of amino acid residues. Calculation of individual contributions of amino acid residues to such contacts permits identification of structurally-important residues. The contact plots obtained with molecular hydrophobicity potential calculations, provide easy rules to choose sites for mutations, which can increase a strength of intra- or inter-molecular hydrophobic interactions. The unfavorable hydrophobic-hydrophilic contact can be mutated to favorable hydrophobic, and already existing weak hydrophobic contact can be strengthen by increasing hydrophobicity of residues in contact. Basing on the analysis of the contact plots, we suggest several mutations of barnase which are supposed to increase intramolecular hydrophobic interactions, and thus might lead to increased stability of the protein. Part of these mutations was studied previously experimentally, and indeed stabilized barnase. The other of predicted mutations were not studied experimentally yet. Several new mutations of barnase and barstar are also proposed to enhance the hydrophobic interactions on their binding interface.

[The rational evolution of scorpion toxins]. / Ratsional'naya 'evolyutsiya toksinov iz yada skorpionov

A theoretical method for the rational design of a "universal" scorpion toxin with a wider spectrum of specificity for K+ channels and a more stable alpha/beta-folding than in its natural homologues is described. On the basis of the analysis of molecular hydrophobic potentials (MHP) of the protein spatial structures, structural features for a family of five short scorpion toxins were revealed. The analysis of the maps of two-dimensional intramolecular MHP contacts allowed the identification of amino acid residues responsible for the folding of the protein and/or for the manifestation of its specific function. The theoretically predicted structure-function roles of the residues were compared with experimental data on the mutagenesis of charybdotoxin. Based on the results of MHP calculations and with the theory of protein molecular evolution used as an additional criterion for the selection of mutations, the amino acid sequence and the spatial structure of a "universal" scorpion toxin were determined.

Amino acid residue: is it structural or functional?

A new approach is suggested for delineating the structural and functional amino acid residues in proteins with known three-dimensional structure, basing on the involvement of residues in intramolecular hydrophobic and hydrophilic interactions and additional information about the conservativity of the residues. The approach is applied to the families of homologous neurotoxins and cardiotoxins. The results obtained concerning the role of amino acid residues in both families of toxins accord well with the similarity of their fold, but different mechanisms of action. Current approach can be used for detailed characterization of protein spatial structures, as well as for rational protein engineering.

[Refinement of the spatial structure of neurotoxin II from Naja naja oxiana venom]. / Utochnenie prostranstvennoi struktury neirotoksina II iz iada Naja naja oxiana.

A set of 19 conformations of the neurotoxin II from Naja naja oxiana was determined by conformational energy minimization using constraints derived from experimental 1H NMR data. The pairwise average root-mean-square deviations were 0.86 A for the backbone heavy atoms and 1.48 A for all heavy atoms of these conformations. A model of the neurotoxin II dimer is proposed to account for the relatively slow deuterium exchange rates of the Val45 and Leu51 amide protons, which are exposed to the solvent in the calculated conformations of monomeric neurotoxin II. Both the monomeric and dimeric models of neurotoxin II may be useful for detailed studies of the functional, hydrophobic, and electrostatic properties of this molecule.

Detailed assessment of spatial hydrophobic and electrostatic properties of 2D NMR-derived models of neurotoxin II.

2D NMR-derived spatial structures of neurotoxin II (NtII) and several homologous toxins in solution were assessed by comparison with their own amino acid sequences using a three-dimensional (3D) profile method. 3D profiles of all the toxin models match the sequences well and, therefore, the method of 3D profile was demonstrated to work correctly for these well-resolved NMR structures in aqueous solution. At the same time, the profile window plots reveal low scores in the bottom tip of loop II (residues 22-34 in NtII) and in beta-strand of loop III (residues 49-52). Some residues in the first poor-scoring region are of functional importance being involved in binding with nicotinic acetylcholine receptor (AChR). Furthermore, the second segment participates in intermolecular hydrogen bonding upon dimerization of postsynaptic neurotoxins in solution resulting in increasing of the 3D-1D score for residues at the interface between monomers. Therefore, the 3D profile method can be useful for detection functionally-important regions in well-resolved protein structures.

Two-dimensional 1H-NMR study of the spatial structure of neurotoxin II from Naja naja oxiana.

The spatial structure of neurotoxin II from the venom of the central Asian cobra Naja naja oxiana was determined by two-dimensional 1H-NMR techniques and computational analysis. Nearly complete proton resonance assignments for 61 amino acid residues have been made using two-dimensional (2D) homonuclear total correlated spectroscopy, 2D homonuclear double-quantum-filtered correlated spectroscopy and 2D homonuclear NOE spectroscopy (NOESY) experiments. The cross-peak volumes in NOESY spectra spin-spin coupling constants of vicinal protons NH-C alpha H and C alpha H-C beta H and the observation of slow deuterium exchange of amide protons were used to define local structure and a set of constraints for distance geometry program DIANA. The average root-mean-square deviations are 53 pm for backbone heavy atoms and 118 pm for all heavy atoms of 19 final neurotoxin II conformations. The spatial structure is characterized by a short double-stranded (residues 1-5 and 13-17) and a triple-stranded (residues 22-30, 33-41 and 50-54) antiparallel beta-sheets.

The divalent cation-binding sites of gramicidin A transmembrane ion-channel.

The conductance of the gramicidin A single channels in glycerolmonooleate membranes is strongly reduced in the presence of Mn2+ cations. The nmr experiments were performed for N-terminal to N-terminal gramicidin A dimer formed by two right-handed single-stranded helixes incorporated into the sodium dodecyl sulfate micelles in the presence of Mn2+ ions. Dependence of the nonselective spin-lattice relaxation rates of the gramicidin A protons on Mn2+ concentration was analyzed to determine coordinates of the divalent cation binding sites. It is inferred that Mn2+ ions are bound at the channel mouths at distances of 6.4, 8.6, and 8.8 A (+/- 2 A) from the oxygen atoms of exposed carbonyl groups of D-Leu 12, 14, and 10, respectively. The bounded Mn2+ retains its hydrate shell, the size of which (approximately 6 A) exceeds the inner pore diameter (approximately 4 A). That makes the gramicidin A channel impermeable for divalent cations.