Tripartite Split-GFP for High Throughput Screening of Small Molecules: A Powerful Strategy for Targeting Transient/Labile Interactors like E2-E3 Ubiquitination Enzymes.

The search for inhibitors of the Ubiquitin Proteasome System (UPS) is an expanding area, due to the crucial role of UPS enzymes in several diseases. The complexity of the UPS and the multiple protein-protein interactions (PPIs) involved, either between UPS proteins themselves or between UPS components and theirs targets, offer an incredibly wide field for the development of chemical compounds for specifically modulating or inhibiting metabolic pathways. However, numerous UPS PPIs are transient/labile, due the processivity of the system (Ubiquitin [Ub] chain elongation, Ub transfer, etc.). Among the different strategies that can be used either for deciphering UPS PPI or for identifying/characterizing small compounds inhibitors, the split-GFP approach offers several advantages notably for high throughput screening of drugs. Split-GFP is based on the principle of protein-fragment complementation assay (PCA). PCA allows addressing PPIs by coupling each protein of interest (POI) to fragments of a reporter protein whose reconstitution is linked to the interaction of the POI. Here, we review the evolution of the split-GFP approach from bipartite to tripartite Split-GFP and its recent applicability for screening chemical compounds targeting the UPS.

Insights into animal septins using recombinant human septin octamers with distinct SEPT9 isoforms.

Septin GTP-binding proteins contribute essential biological functions that range from the establishment of cell polarity to animal tissue morphogenesis. Human septins in cells form hetero-octameric septin complexes containing the ubiquitously expressed SEPT9 subunit (also known as SEPTIN9). Despite the established role of SEPT9 in mammalian development and human pathophysiology, biochemical and biophysical studies have relied on monomeric SEPT9, thus not recapitulating its native assembly into hetero-octameric complexes. We established a protocol that enabled, for the first time, the isolation of recombinant human septin octamers containing distinct SEPT9 isoforms. A combination of biochemical and biophysical assays confirmed the octameric nature of the isolated complexes in solution. Reconstitution studies showed that octamers with either a long or a short SEPT9 isoform form filament assemblies, and can directly bind and cross-link actin filaments, raising the possibility that septin-decorated actin structures in cells reflect direct actin-septin interactions. Recombinant SEPT9-containing octamers will make it possible to design cell-free assays to dissect the complex interactions of septins with cell membranes and the actin and microtubule cytoskeleton.

Modulation of innate immune signaling by a Coxiella burnetii eukaryotic-like effector protein.

The Q fever agent Coxiella burnetii uses a defect in organelle trafficking/intracellular multiplication (Dot/Icm) type 4b secretion system (T4SS) to silence the host innate immune response during infection. By investigating C. burnetii effector proteins containing eukaryotic-like domains, here we identify NopA (nucleolar protein A), which displays four regulator of chromosome condensation (RCC) repeats, homologous to those found in the eukaryotic Ras-related nuclear protein (Ran) guanine nucleotide exchange factor (GEF) RCC1. Accordingly, NopA is found associated with the chromatin nuclear fraction of cells and uses the RCC-like domain to interact with Ran. Interestingly, NopA triggers an accumulation of Ran-GTP, which accumulates at nucleoli of transfected or infected cells, thus perturbing the nuclear import of transcription factors of the innate immune signaling pathway. Accordingly, qRT-PCR analysis on a panel of cytokines shows that cells exposed to the C. burnetii nopA::Tn or a Dot/Icm-defective dotA::Tn mutant strain present a functional innate immune response, as opposed to cells exposed to wild-type C. burnetii or the corresponding nopA complemented strain. Thus, NopA is an important regulator of the innate immune response allowing Coxiella to behave as a stealth pathogen.

Solution structure of the type I polyketide synthase Pks13 from Mycobacterium tuberculosis.

BACKGROUND: Type I polyketide synthases (PKSs) are multifunctional enzymes responsible for the biosynthesis of a group of diverse natural compounds with biotechnological and pharmaceutical interest called polyketides. The diversity of polyketides is impressive despite the limited set of catalytic domains used by PKSs for biosynthesis, leading to considerable interest in deciphering their structure-function relationships, which is challenging due to high intrinsic flexibility. Among nineteen polyketide synthases encoded by the genome of Mycobacterium tuberculosis, Pks13 is the condensase required for the final condensation step of two long acyl chains in the biosynthetic pathway of mycolic acids, essential components of the cell envelope of Corynebacterineae species. It has been validated as a promising druggable target and knowledge of its structure is essential to speed up drug discovery to fight against tuberculosis. RESULTS: We report here a quasi-atomic model of Pks13 obtained using small-angle X-ray scattering of the entire protein and various molecular subspecies combined with known high-resolution structures of Pks13 domains or structural homologues. As a comparison, the low-resolution structures of two other mycobacterial polyketide synthases, Mas and PpsA from Mycobacterium bovis BCG, are also presented. This study highlights a monomeric and elongated state of the enzyme with the apo- and holo-forms being identical at the resolution probed. Catalytic domains are segregated into two parts, which correspond to the condensation reaction per se and to the release of the product, a pivot for the enzyme flexibility being at the interface. The two acyl carrier protein domains are found at opposite sides of the ketosynthase domain and display distinct characteristics in terms of flexibility. CONCLUSIONS: The Pks13 model reported here provides the first structural information on the molecular mechanism of this complex enzyme and opens up new perspectives to develop inhibitors that target the interactions with its enzymatic partners or between catalytic domains within Pks13 itself.

High-content tripartite split-GFP cell-based assays to screen for modulators of small GTPase activation.

The human Ras superfamily of small GTPases controls essential cellular processes such as gene expression and cell proliferation. As their deregulation is widely associated with human cancer, small GTPases and their regulatory proteins have become increasingly attractive for the development of novel therapeutics. Classical methods to monitor GTPase activation include pulldown assays that limit the analysis of GTP-bound form of proteins from cell lysates. Alternatively, live-cell FRET biosensors may be used to study GTPase activation dynamics in response to stimuli, but these sensors often require further optimization for high-throughput applications. Here, we describe a cell-based approach that is suitable to monitor the modulation of small GTPase activity in a high-content analysis. The assay relies on a genetically encoded tripartite split-GFP (triSFP) system that we integrated in an optimized cellular model to monitor modulation of RhoA and RhoB GTPases. Our results indicate the robust response of the reporter, allowing the interrogation of inhibition and stimulation of Rho activity, and highlight potential applications of this method to discover novel modulators and regulators of small GTPases and related protein-binding domains.

Cytoplasmic p27Kip1 promotes tumorigenesis via suppression of RhoB activity.

The cell cycle inhibitor p27Kip1 is a tumor suppressor via the inhibition of CDK complexes in the nucleus. However, p27 also plays other functions in the cell and may acquire oncogenic roles when located in the cytoplasm. Activation of oncogenic pathways such as Ras or PI3K/AKT causes the relocalization of p27 in the cytoplasm, where it can promote tumorigenesis by unclear mechanisms. Here, we investigated how cytoplasmic p27 participates in the development of non-small cell lung carcinomas. We provide molecular and genetic evidence that the oncogenic role of p27 is mediated, at least in part, by binding to and inhibiting the GTPase RhoB, which normally acts as a tumor suppressor in the lung. Genetically modified mice revealed that RhoB expression is preferentially lost in tumors in which p27 is absent and maintained in tumors expressing wild-type p27 or p27CK- , a mutant that cannot inhibit CDKs. Moreover, although the absence of RhoB promoted tumorigenesis in p27-/- animals, it had no effect in p27CK- knock-in mice, suggesting that cytoplasmic p27 may act as an oncogene, at least in part, by inhibiting the activity of RhoB. Finally, in a cohort of lung cancer patients, we identified a subset of tumors harboring cytoplasmic p27 in which RhoB expression is maintained and these characteristics were strongly associated with decreased patient survival. Thus, monitoring p27 localization and RhoB levels in non-small cell lung carcinoma patients appears to be a powerful prognostic marker for these tumors. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Development and Applications of Superfolder and Split Fluorescent Protein Detection Systems in Biology.

Molecular engineering of the green fluorescent protein (GFP) into a robust and stable variant named Superfolder GFP (sfGFP) has revolutionized the field of biosensor development and the use of fluorescent markers in diverse area of biology. sfGFP-based self-associating bipartite split-FP systems have been widely exploited to monitor soluble expression in vitro, localization, and trafficking of proteins in cellulo. A more recent class of split-FP variants, named « tripartite ¼ split-FP, that rely on the self-assembly of three GFP fragments, is particularly well suited for the detection of protein-protein interactions. In this review, we describe the different steps and evolutions that have led to the diversification of superfolder and split-FP reporter systems, and we report an update of their applications in various areas of biology, from structural biology to cell biology.

Water Distribution within Wild-Type NRas Protein and Q61 Mutants during Unrestrained QM/MM Dynamics.

Point mutations in p21ras are associated with ∼30% of human tumors by disrupting its GTP hydrolysis cycle, which is critical to its molecular switch function in cellular signaling pathways. In this work, we investigate the impact of Gln 61 substitutions in the structure of the p21N-ras active site and particularly focus on water reorganization around GTP, which appears to be crucial to evaluate favorable and unfavorable hydration sites for hydrolysis. The NRas-GTP complex is analyzed using a hybrid quantum mechanics/molecular mechanics approach, treating for the first time to our knowledge transient water molecules at the ab initio level and leading to results that account for the electrostatic coupling between the protein complex and the solvent. We show that for the wild-type protein, water molecules are found around the GTP γ-phosphate group, forming an arch extended from residues 12 to 35. Two density peaks are observed, supporting previous results that suggest the presence of two water molecules in the active site, one in the vicinity of residue 35 and a second one stabilized by hydrogen bonds formed with nitrogen backbone atoms of residues 12 and 60. The structural changes observed in NRas Gln 61 mutants result in the drastic delocalization of water molecules that we discuss. In mutants Q61H and Q61K, for which water distribution is overlocalized next to residue 60, the second density peak supports the hypothesis of a second water molecule. We also conclude that Gly 60 indirectly participates in GTP hydrolysis by correctly positioning transient water molecules in the protein complex and that Gln 61 has an indirect steric effect in stabilizing the preorganized catalytic site.

In vivo interactions of TTDA mutant proteins within TFIIH.

Trichothiodystrophy group A (TTD-A) patients carry a mutation in the transcription factor II H (TFIIH) subunit TTDA. Using a novel in vivo tripartite split-GFP system, we show that TTDA interacts with the TFIIH subunit p52 and the p52-TTDA-GFP product is incorporated into TFIIH. p52-TTDA-GFP is able to bind DNA and is recruited to UV-damaged DNA. Furthermore, we show that two patient-mutated TTDA proteins can interact with p52, are able to bind to the DNA and can localize to damaged DNA. Our findings give new insights into the behavior of TTDA within the context of a living cell and thereby shed light on the complex phenotype of TTD-A patients.

4'-Phosphopantetheinyl transferase PptT, a new drug target required for Mycobacterium tuberculosis growth and persistence in vivo.

The cell envelope of Mycobacterium tuberculosis, the causative agent of tuberculosis in humans, contains lipids with unusual structures. These lipids play a key role in both virulence and resistance to the various hostile environments encountered by the bacteria during infection. They are synthesized by complex enzymatic systems, including type-I polyketide synthases and type-I and -II fatty acid synthases, which require a post-translational modification to become active. This modification consists of the covalent attachment of the 4'-phosphopantetheine moiety of Coenzyme A catalyzed by phosphopantetheinyl transferases (PPTases). PptT, one of the two PPTases produced by mycobacteria, is involved in post-translational modification of various type-I polyketide synthases required for the formation of both mycolic acids and lipid virulence factors in mycobacteria. Here we identify PptT as a new target for anti-tuberculosis drugs; we address all the critical issues of target validation to demonstrate that PptT can be used to search for new drugs. We confirm that PptT is essential for the growth of M. bovis BCG in vitro and show that it is required for persistence of M. bovis BCG in both infected macrophages and immunodeficient mice. We generated a conditional expression mutant of M. tuberculosis, in which the expression of the pptT gene is tightly regulated by tetracycline derivatives. We used this construct to demonstrate that PptT is required for the replication and survival of the tubercle bacillus during the acute and chronic phases of infection in mice. Finally, we developed a robust and miniaturized assay based on scintillation proximity assay technology to search for inhibitors of PPTases, and especially of PptT, by high-throughput screening. Our various findings indicate that PptT meets the key criteria for being a therapeutic target for the treatment of mycobacterial infections.

Glucocorticoid receptor-dependent therapeutic efficacy of tauroursodeoxycholic acid in preclinical models of spinocerebellar ataxia type 3.

Spinocerebellar ataxia type 3 (SCA3) is an adult-onset neurodegenerative disease caused by a polyglutamine expansion in the ataxin-3 (ATXN3) gene. No effective treatment is available for this disorder, other than symptom-directed approaches. Bile acids have shown therapeutic efficacy in neurodegenerative disease models. Here, we pinpointed tauroursodeoxycholic acid (TUDCA) as an efficient therapeutic, improving the motor and neuropathological phenotype of SCA3 nematode and mouse models. Surprisingly, transcriptomic and functional in vivo data showed that TUDCA acts in neuronal tissue through the glucocorticoid receptor (GR), but independently of its canonical receptor, the farnesoid X receptor (FXR). TUDCA was predicted to bind to the GR, in a similar fashion to corticosteroid molecules. GR levels were decreased in disease-affected brain regions, likely due to increased protein degradation as a consequence of ATXN3 dysfunction being restored by TUDCA treatment. Analysis of a SCA3 clinical cohort showed intriguing correlations between the peripheral expression of GR and the predicted age at disease onset in presymptomatic subjects and FKBP5 expression with disease progression, suggesting this pathway as a potential source of biomarkers for future study. We have established a novel in vivo mechanism for the neuroprotective effects of TUDCA in SCA3 and propose this readily available drug for clinical trials in SCA3 patients.

Detection of soluble co-factor dependent protein expression in vivo: application to the 4'-phosphopantetheinyl transferase PptT from Mycobacterium tuberculosis.

The need for early-on diagnostic tools to assess the folding and solubility of expressed protein constructs in vivo is of great interest when dealing with recalcitrant proteins. In this paper, we took advantage of the picomolar sensitivity of the bipartite GFP1-10/GFP11 system to investigate the solubility of the Mycobacterium tuberculosis 4'-phosphopantetheinyl transferase PptT, an enzyme essential for the viability of the tubercle bacillus. In vivo and in vitro complementation assays clearly showed the improved solubility of the full-length PptT compared to its N- and C-terminally truncated counterparts. However, initial attempts to purify the full-length enzyme overexpressed in Escherichia coli cells were hampered by aggregation issues overtime that caused the protein to precipitate within hours. The fact that the naturally occurring Coenzyme A and Mg(2+), essentials for PptT to carry out its function, could play a role in stabilizing the enzyme was confirmed using DSF experiments. In vitro activity assays were performed using the ACP substrate from the type I polyketide synthase PpsC from M. tuberculosis, a 2188 amino-acid enzyme that plays a major role in the virulence and pathogenicity of this microbial pathogen. We selected the most soluble and compact ACP fragment (2042-2188), identified by genetic selection of in-frame fragments from random library experiments, to monitor the transfer of the P-pant moiety from Coenzyme A onto a conserved serine residue of this ACP domain.

Experimental mapping of soluble protein domains using a hierarchical approach.

Exploring the function and 3D space of large multidomain protein targets often requires sophisticated experimentation to obtain the targets in a form suitable for structure determination. Screening methods capable of selecting well-expressed, soluble fragments from DNA libraries exist, but require the use of automation to maximize chances of picking a few good candidates. Here, we describe the use of an insertion dihydrofolate reductase (DHFR) vector to select in-frame fragments and a split-GFP assay technology to filter-out constructs that express insoluble protein fragments. With the incorporation of an IPCR step to create high density, focused sublibraries of fragments, this cost-effective method can be performed manually with no a priori knowledge of domain boundaries while permitting single amino acid resolution boundary mapping. We used it on the well-characterized p85α subunit of the phosphoinositide-3-kinase to demonstrate the robustness and efficiency of our methodology. We then successfully tested it onto the polyketide synthase PpsC from Mycobacterium tuberculosis, a potential drug target involved in the biosynthesis of complex lipids in the cell envelope. X-ray quality crystals from the acyl-transferase (AT), dehydratase (DH) and enoyl-reductase (ER) domains have been obtained.

Visualizing the subcellular localization of RHOB-GTP and GTPase-Effector complexes using a split-GFP/nanobody labelling assay.

Small GTPases are highly regulated proteins that control essential signaling pathways through the activity of their effector proteins. Among the RHOA subfamily, RHOB regulates peculiar functions that could be associated with the control of the endocytic trafficking of signaling proteins. Here, we used an optimized assay based on tripartite split-GFP complementation to localize GTPase-effector complexes with high-resolution. The detection of RHOB interaction with the Rhotekin Rho binding domain (RBD) that specifically recognizes the active GTP-bound GTPase, is performed in vitro by the concomitant addition of recombinant GFP1-9 and a GFP nanobody. Analysis of RHOB-RBD complexes localization profiles combined with immunostaining and live cell imaging indicated a serum-dependent reorganization of the endosomal and membrane pool of active RHOB. We further applied this technology to the detection of RHO-effector complexes that highlighted their subcellular localization with high resolution among the different cellular compartments.

Human septins organize as octamer-based filaments and mediate actin-membrane anchoring in cells.

Septins are cytoskeletal proteins conserved from algae and protists to mammals. A unique feature of septins is their presence as heteromeric complexes that polymerize into filaments in solution and on lipid membranes. Although animal septins associate extensively with actin-based structures in cells, whether septins organize as filaments in cells and if septin organization impacts septin function is not known. Customizing a tripartite split-GFP complementation assay, we show that all septins decorating actin stress fibers are octamer-containing filaments. Depleting octamers or preventing septins from polymerizing leads to a loss of stress fibers and reduced cell stiffness. Super-resolution microscopy revealed septin fibers with widths compatible with their organization as paired septin filaments. Nanometer-resolved distance measurements and single-protein tracking further showed that septin filaments are membrane bound and largely immobilized. Finally, reconstitution assays showed that septin filaments mediate actin-membrane anchoring. We propose that septin organization as octamer-based filaments is essential for septin function in anchoring and stabilizing actin filaments at the plasma membrane.

The Brucella TIR-like protein TcpB interacts with the death domain of MyD88.

The pathogen Brucella melitensis secretes a Toll/interleukin-1 receptor (TIR) domain containing protein that abrogates host innate immune responses. In this study, we have characterized the biochemical interactions of Brucella TIR-like protein TcpB with host innate immune adaptor proteins. Using protein-fragment complementation assays based on Gaussia luciferase and green fluorescent protein, we find that TcpB interacts directly with MyD88 and that this interaction is significantly stronger than the interaction of TcpB with TIRAP, the only other adaptor protein that detectably interacts with TcpB. Surprisingly, the TcpB-MyD88 interaction depends on the death domain (DD) of MyD88, and TcpB does not interact with the isolated TIR domain of MyD88. TcpB disrupts MyD88(DD)-MyD88(DD), MyD88(DD)-MyD88(TIR) and MyD88(DD)-MyD88 interactions but not MyD88-MyD88 or MyD88(TIR)-MyD88(TIR) interactions. Structural models consistent with these results suggest how TcpB might inhibit TLR signaling by targeting MyD88 via a DD-TIR domain interface.

Tripartite split-GFP assay to identify selective intracellular nanobody that suppresses GTPase RHOA subfamily downstream signaling.

Strategies based on intracellular expression of artificial binding domains present several advantages over manipulating nucleic acid expression or the use of small molecule inhibitors. Intracellularly-functional nanobodies can be considered as promising macrodrugs to study key signaling pathways by interfering with protein-protein interactions. With the aim of studying the RAS-related small GTPase RHOA family, we previously isolated, from a synthetic phage display library, nanobodies selective towards the GTP-bound conformation of RHOA subfamily proteins that lack selectivity between the highly conserved RHOA-like and RAC subfamilies of GTPases. To identify RHOA/ROCK pathway inhibitory intracellular nanobodies, we implemented a stringent, subtractive phage display selection towards RHOA-GTP followed by a phenotypic screen based on F-actin fiber loss. Intracellular interaction and intracellular selectivity between RHOA and RAC1 proteins was demonstrated by adapting the sensitive intracellular protein-protein interaction reporter based on the tripartite split-GFP method. This strategy led us to identify a functional intracellular nanobody, hereafter named RH28, that does not cross-react with the close RAC subfamily and blocks/disrupts the RHOA/ROCK signaling pathway in several cell lines without further engineering or functionalization. We confirmed these results by showing, using SPR assays, the high specificity of the RH28 nanobody towards the GTP-bound conformation of RHOA subfamily GTPases. In the metastatic melanoma cell line WM266-4, RH28 expression triggered an elongated cellular phenotype associated with a loss of cellular contraction properties, demonstrating the efficient intracellular blocking of RHOA/B/C proteins downstream interactions without the need of manipulating endogenous gene expression. This work paves the way for future therapeutic strategies based on protein-protein interaction disruption with intracellular antibodies.

Vγ9Vδ2 T Cells Activation Through Phosphoantigens Can Be Impaired by a RHOB Rerouting in Lung Cancer.

Vγ9Vδ2 T cells are known to be efficient anti-tumor effectors activated through phosphoantigens (PAg) that are naturally expressed by tumor cells or induced by amino bisphosphonates treatment. This PAg-activation which is TCR and butyrophilin BTN3A dependent can be modulated by NKG2D ligands, immune checkpoint ligands, adhesion molecules, and costimulatory molecules. This could explain the immune-resistance observed in certain clinical trials based on Vγ9Vδ2 T cells therapies. In NSCLC, encouraging responses were obtained with zoledronate administrations for 50% of patients. According to the in vivo results, we showed that the in vitro Vγ9Vδ2 T cell reactivity depends on the NSCLC cell line considered. If the PAg-pretreated KRAS mutated A549 is highly recognized and killed by Vγ9Vδ2 T cells, the EGFR mutated PC9 remains resistant to these killers despite a pre-treatment either with zoledronate or with exogenous BrHPP. The immune resistance of PC9 was shown not to be due to immune checkpoint ligands able to counterbalance NKG2D ligands or adhesion molecules such as ICAM-1 highly expressed by PC9. RHOB has been shown to be involved in the Vγ9Vδ2 TCR signaling against these NSCLC cell lines, in this study we therefore focused on its intracellular behavior. In comparison to a uniform distribution of RHOB in endosomes and at the plasma membrane in A549, the presence of large endosomal clusters of RHOB was visualized by a split-GFP system, suggesting that RHOB rerouting in the PC9 tumor cell could impair the reactivity of the immune response.

Engineering and characterization of a superfolder green fluorescent protein.

Existing variants of green fluorescent protein (GFP) often misfold when expressed as fusions with other proteins. We have generated a robustly folded version of GFP, called 'superfolder' GFP, that folds well even when fused to poorly folded polypeptides. Compared to 'folding reporter' GFP, a folding-enhanced GFP containing the 'cycle-3' mutations and the 'enhanced GFP' mutations F64L and S65T, superfolder GFP shows improved tolerance of circular permutation, greater resistance to chemical denaturants and improved folding kinetics. The fluorescence of Escherichia coli cells expressing each of eighteen proteins from Pyrobaculum aerophilum as fusions with superfolder GFP was proportional to total protein expression. In contrast, fluorescence of folding reporter GFP fusion proteins was strongly correlated with the productive folding yield of the passenger protein. X-ray crystallographic structural analyses helped explain the enhanced folding of superfolder GFP relative to folding reporter GFP.

Hybrid QM/MM vs Pure MM Molecular Dynamics for Evaluating Water Distribution within p21N-ras and the Resulting GTP Electronic Density.

p21ras protein activity, regulated by GTP hydrolysis, constitutes an active field of research for the development of cancer targeted therapies that would concern ∼30% of human tumors to which specific mutations have been associated. Indeed, the catalyzing mechanisms provided by the protein environment during GTP hydrolysis and how they are impaired by specific mutations remain to be fully elucidated. In this article, we present results from molecular mechanics (MM) molecular dynamics (MD) simulations and density functional theory (DFT) calculations carried out for wild-type p21 N-ras and six Gln 61 mutants. In the first part, we present the water distribution within the active site of the wild-type protein according to MM MD. Significant differences are observed when comparing the results to the previous distribution assessed through quantum mechanics/molecular mechanics (QM/MM) MD. Such method-dependent results highlight the importance of accounting for the electrostatic coupling between the protein complex and the solvent molecules in identifying hydration sites. In the second part, we present the results from DFT calculations performed to determine the electronic distribution of the GTP ligand, considering the wild-type active site arrangement according to both classical and hybrid approaches. Only in the QM/MM-based configuration is the ligand electronic density similar to that of a GDP-like state observed experimentally. For this reason, in the last set of calculations carried out for p21 N-ras Gln 61 mutants, only the active site structural conformations obtained through hybrid MD are considered. Through the analysis of the GTP electronic density, we conclude that the wild-type active site arrangement according to QM/MM MD is closer to a catalytically efficient conformation of the protein than the arrangement according to MM MD. Hence, water distribution according to the hybrid approach must correspond to the optimal placement of solvent in the active site. Within all of the studied Gln 61 substituted proteins, p21ras major catalyzing effect, which consists of stabilizing a more GDP-like state, is lost.