Modeling the mutational effects on calmodulin structure: prediction of alteration in the amino acid interactions

Characterization and functional analysis of the calmodulin-binding domain of Rac1 GTPase

Rac1, a member of the Rho family of small GTPases, has been shown to promote formation of lamellipodia at the leading edge of motile cells and affect cell migration. We previously demonstrated that calmodulin can bind to a region in the C-terminal of Rac1 and that this interaction is important in the activation of platelet Rac1. Now, we have analyzed amino acid residue(s) in the Rac1-calmodulin binding domain that are essential for the interaction and assessed their functional contribution in Rac1 activation. The results demonstrated that region 151–164 in Rac1 is essential for calmodulin binding. Within the 151–164 region, positively-charged amino acids K153 and R163 were mutated to alanine to study impact on calmodulin binding. Mutant form of Rac1 (K153A) demonstrated significantly reduced binding to calmodulin while the double mutant K153A/R163A demonstrated complete lack of binding to calmodulin. Thrombin or EGF resulted in activation of Rac1 in CHRF-288-11 or HeLa cells respectively and W7 inhibited this activation. Immunoprecipitation studies demonstrated that higher amount of CaM was associated with Rac1 during EGF dependent activation. In cells expressing mutant forms of Rac1 (K153A or K153A/R163A), activation induced by EGF was significantly decreased in comparison to wild type or the R163A forms of Rac1. The lack of Rac1 activation in mutant forms was not due to an inability of GDP-GTP exchange or a change in subcelllular distribution. Moreover, Rac1 activation was decreased in cells where endogenous level of calmodulin was reduced using shRNA knockdown and increased in cells where calmodulin was overexpressed. Docking analysis and modeling demonstrated that K153 in Rac1 interacts with Q41 in calmodulin. These results suggest an important role for calmodulin in the activation of Rac1 and thus, in cytoskeleton reorganization and cell migration.

Modeling the 3D structure of wheat subtilisin/chymotrypsin inhibitor (WSCI). Probing the reactive site with two susceptible proteinases by time-course analysis and molecular dynamics simulations

Comparative modeling and time-course hydrolysis experiments have been applied to investigate two enzyme-inhibitor complexes formed between the wheat subtilisin-chymotrypsin inhibitor (WSCI) and two susceptible proteinases. WSCI represents the first case of a wheat protein inhibitor active against animal chymotrypsins and bacterial subtilisins. The model was created using as template structure that of the CI-2A inhibitor from barley (PDB code: 2CI2), which shares 87% sequence identity with WSCI. Under these conditions of high similarity, the comparative modeling approach can be successfully applied. We predicted the WSCI 3D model and used it to investigate enzyme-inhibitor complex systems. Experimental observations indicated that chymotrypsin, but not subtilisin, in addition to cleavage at the primary reactive site Met48-Glu49, is able to hydrolyze a second peptide bond between Phe58 and Val59. Here, we report on cleavage of the peptide bond at the inhibitor's reactive site (Met48-Glu49) determined using time-course hydrolysis experiments; the same event was investigated for both subtilisin/WSCI and chymotrypsin/WSCI complexes using molecular dynamics simulations. The molecular details of the initial inhibitor-enzyme interactions, as well as of the changes observed during the simulations, allow us to speculate on the different fates of the two WSCI-proteinase complexes.