Molecular dynamics simulations reveal differences in the conformational stability of FtsZs derived from Staphylococcus aureus and Bacillus subtilis.

FtsZ is highly conserved among bacteria and plays an essential role in bacterial cell division. The tense conformation of FtsZ bound to GTP assembles into a straight filament via head-to-tail associations, and then the upper subunit of FtsZ hydrolyzes GTP bound to the lower FtsZ subunit. The subunit with GDP bound disassembles accompanied by a conformational change in the subunit from the tense to relaxed conformation. Although crystal structures of FtsZ derived from several bacterial species have been determined, the conformational change from the relaxed to tense conformation has only been observed in Staphylococcus aureus FtsZ (SaFtsZ). Recent cryo-electron microscopy analyses revealed the three-dimensional reconstruction of the protofilament, in which tense molecules assemble via head-to-tail associations. However, the lower resolution of the protofilament suggested that the flexibility of the FtsZ protomers between the relaxed and tense conformations caused them to form in less-strict alignments. Furthermore, this flexibility may also prevent FtsZs other than SaFtsZ from crystalizing in the tense conformation, suggesting that the flexibility of bacterial FtsZs differs. In this study, molecular dynamics simulations were performed using SaFtsZ and Bacillus subtilis FtsZ in several situations, which suggested that different features of the FtsZs affect their conformational stability.

Crystal-packing analysis of translation initiation factor 2 reveals new details of its function.

Eukaryotic and archaeal translation initiation factor 2 in complex with GTP delivers the initiator methionyl-tRNA to the small ribosomal subunit. Over the past 20 years, thanks to the efforts of various research groups, including ours, this factor from the archaeon Sulfolobus solfataricus and its individual subunits have been crystallized in ten different space groups. Analysis of the molecular packing in these crystals makes it possible to better understand the roles of functionally significant switches and other elements of the nucleotide-binding pocket during the function of the factor as well as the influence of external effects on its transition between active and inactive states.

Effects of bound nucleotides on the secondary structure, thermal stability, and phosphorylation of Rab3A.

Rab3A is a member of the Rab GTPase family involved in synaptic vesicle trafficking. Recent evidence has demonstrated that Rab3A is phosphorylated by leucine-rich repeat kinase 2 (LRRK2) that is implicated in both familial and sporadic forms of Parkinson's disease (PD), and an abnormal increase in Rab3A phosphorylation has been proposed as a cause of PD. Despite the potential importance of Rab3A in PD pathogenesis, its structural information is limited and the effects of bound nucleotides on its biophysical and biochemical properties remain unclear. Here, we show that GDP-bound Rab3A is preferentially phosphorylated by LRRK2 compared with GTP-bound Rab3A. The secondary structure of Rab3A, measured by circular dichroism (CD) spectroscopy, revealed that Rab3A is resistant to heat-induced denaturation at pH 7.4 or 9.0 regardless of the nucleotides bound. In contrast, Rab3A underwent heat-induced denaturation at pH 5.0 at a lower temperature in its GDP-bound form than in its GTP-bound form. The unfolding temperature of Rab3A was studied by differential scanning fluorimetry, which showed a significantly higher unfolding temperature in GTP-bound Rab3A than in GDP-bound Rab3A, with the highest at pH 7.4. These results suggest that Rab3A has unusual thermal stability under physiologically relevant conditions and that bound nucleotides influence both thermal stability and phosphorylation by LRRK2.

Molecular Mechanism of Phosphorylation-Mediated Impacts on the Conformation Dynamics of GTP-Bound KRAS Probed by GaMD Trajectory-Based Deep Learning.

The phosphorylation of different sites produces a significant effect on the conformational dynamics of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations were combined with deep learning (DL) to explore the molecular mechanism of the phosphorylation-mediated effect on conformational dynamics of the GTP-bound KRAS. The DL finds that the switch domains are involved in obvious differences in conformation contacts and suggests that the switch domains play a key role in the function of KRAS. The analyses of free energy landscapes (FELs) reveal that the phosphorylation of pY32, pY64, and pY137 leads to more disordered states of the switch domains than the wild-type (WT) KRAS and induces conformational transformations between the closed and open states. The results from principal component analysis (PCA) indicate that principal motions PC1 and PC2 are responsible for the closed and open states of the phosphorylated KRAS. Interaction networks were analyzed and the results verify that the phosphorylation alters interactions of GTP and magnesium ion Mg2+ with the switch domains. It is concluded that the phosphorylation pY32, pY64, and pY137 tune the activity of KRAS through changing conformational dynamics and interactions of the switch domains. We anticipated that this work could provide theoretical aids for deeply understanding the function of KRAS.

Structural insights into the GTP-driven monomerization and activation of a bacterial LRRK2 homolog using allosteric nanobodies.

Roco proteins entered the limelight after mutations in human LRRK2 were identified as a major cause of familial Parkinson's disease. LRRK2 is a large and complex protein combining a GTPase and protein kinase activity, and disease mutations increase the kinase activity, while presumably decreasing the GTPase activity. Although a cross-communication between both catalytic activities has been suggested, the underlying mechanisms and the regulatory role of the GTPase domain remain unknown. Several structures of LRRK2 have been reported, but structures of Roco proteins in their activated GTP-bound state are lacking. Here, we use single-particle cryo-electron microscopy to solve the structure of a bacterial Roco protein (CtRoco) in its GTP-bound state, aided by two conformation-specific nanobodies: NbRoco1 and NbRoco2. This structure presents CtRoco in an active monomeric state, featuring a very large GTP-induced conformational change using the LRR-Roc linker as a hinge. Furthermore, this structure shows how NbRoco1 and NbRoco2 collaborate to activate CtRoco in an allosteric way. Altogether, our data provide important new insights into the activation mechanism of Roco proteins, with relevance to LRRK2 regulation, and suggest new routes for the allosteric modulation of their GTPase activity.

Conformational Modulation of Tissue Transglutaminase via Active Site Thiol Alkylating Agents: Size Does Not Matter.

TG2 is a unique member of the transglutaminase family as it undergoes a dramatic conformational change, allowing its mutually exclusive function as either a cross-linking enzyme or a G-protein. The enzyme's dysregulated activity has been implicated in a variety of pathologies (e.g., celiac disease, fibrosis, cancer), leading to the development of a wide range of inhibitors. Our group has primarily focused on the development of peptidomimetic targeted covalent inhibitors, the nature and size of which were thought to be important features to abolish TG2's conformational dynamism and ultimately inhibit both its activities. However, we recently demonstrated that the enzyme was unable to bind guanosine triphosphate (GTP) when catalytically inactivated by small molecule inhibitors. In this study, we designed a library of models targeting covalent inhibitors of progressively smaller sizes (15 to 4 atoms in length). We evaluated their ability to inactivate TG2 by measuring their respective kinetic parameters kinact and KI. Their impact on the enzyme's ability to bind GTP was then evaluated and subsequently correlated to the conformational state of the enzyme, as determined via native PAGE and capillary electrophoresis. All irreversible inhibitors evaluated herein locked TG2 in its open conformation and precluded GTP binding. Therefore, we conclude that steric bulk and structural complexity are not necessary factors to consider when designing TG2 inhibitors to abolish G-protein activity.

Assessing the mechanism of fast-cycling cancer-associated mutations of Rac1 small Rho GTPase.

Rho-GTPases proteins function as molecular switches alternating from an active to an inactive state upon Guanosine triphosphate (GTP) binding and hydrolysis to Guanosine diphosphate (GDP). Among them, Rac subfamily regulates cell dynamics, being overexpressed in distinct cancer types. Notably, these proteins are object of frequent cancer-associated mutations at Pro29 (P29S, P29L, and P29Q). To assess the impact of these mutations on Rac1 structure and function, we performed extensive all-atom molecular dynamics simulations on wild-type (wt) and oncogenic isoforms of this protein in GDP- and GTP-bound states. Our results unprecedentedly elucidate that P29Q/S-induced structural and dynamical perturbations of Rac1 core domain weaken the binding of the catalytic site Mg2+ ion, and reduce the GDP residence time within protein, enhancing the GDP/GTP exchange rate and Rac1 activity. This broadens our knowledge of the role of cancer-associated mutations on small GTPases mechanism supplying valuable information for future drug discovery efforts targeting specific Rac1 isoforms.

Enhancement of intrinsic guanine fluorescence by protonation in DNA of various structures.

Understanding the diversity of DNA structure and functions in biology requires tools to study this biomolecule selectively and thoroughly. Fluorescence methods are powerful technique for non-invasive research. Due to the low quantum yield, the intrinsic fluorescence of nucleotides has not been considered for use in the detection and differentiation of nucleic acid bases. Here, we have studied the influence of protonation of nucleotides on their fluorescence properties. We show that protonation of ATP and GTP leads to enhanced intrinsic fluorescence. Fluorescence enhancement at acidic pH has been observed for double-stranded DNA and single-stranded oligonucleotides. The formation of G4 secondary structures apparently protected certain nucleotides from protonation, resulting in less pronounced fluorescence enhancement. Furthermore, acid-induced depurination under protonation was less noticeable in G4 structures than in double-stranded and single-stranded DNA. We show that changes in the intrinsic fluorescence of guanine can be used as a sensitive sensor for changes in the structure of the DNA and for the protonation of specific nucleotides.

Molecular dynamics investigation of DNA fragments bound to the anti-HIV protein SAMHD1 reveals alterations in allosteric communications.

The sterile alpha motif and histidine-aspartate domain-containing protein 1 (or SAMHD1), a human dNTP-triphosphohydrolase, contributes to HIV-1 restriction in select terminally differentiated cells of the immune system. While the prevailing hypothesis is that the catalytically active form of the protein is an allosterically triggered tetramer, whose HIV-1 restriction properties are attributed to its dNTP - triphosphohydrolase activity, it is also known to bind to ssRNA and ssDNA oligomers. A complete picture of the structure-function relationship of the enzyme is still elusive and the function corresponding to its nucleic acid binding ability is debated. In this in silico study, we investigate the stability, preference and allosteric effects of DNA oligomers bound to SAMHD1. In particular, we compare the binding of DNA and RNA oligomers of the same sequence and also consider the binding of DNA fragments with phosphorothioate bonds in the backbone. The results are compared with the canonical form with the monomers connected by GTP/dATP crossbridges. The simulations indicate that SAMHD1 dimers preferably bind to DNA and RNA oligomers compared to GTP/dATP. However, allosteric communication channels are altered in the nucleic acid acid bound complexes compared to the canonical form. All results are consistent with the hypothesis that the DNA bound form of the protein correspond to an unproductive off-pathway state where the protein is sequestered and not available for dNTP hydrolysis.

Structure of the GDP-bound state of the SRP GTPase FlhF.

The GTPase FlhF, a signal recognition particle (SRP)-type enzyme, is pivotal for spatial-numerical control and bacterial flagella assembly across diverse species, including pathogens. This study presents the X-ray structure of FlhF in its GDP-bound state at a resolution of 2.28 Å. The structure exhibits the classical N- and G-domain fold, consistent with related SRP GTPases such as Ffh and FtsY. Comparative analysis with GTP-loaded FlhF elucidates the conformational changes associated with GTP hydrolysis. These topological reconfigurations are similarly evident in Ffh and FtsY, and play a pivotal role in regulating the functions of these hydrolases.

Native architecture of a human GBP1 defense complex for cell-autonomous immunity to infection.

All living organisms deploy cell-autonomous defenses to combat infection. In plants and animals, large supramolecular complexes often activate immune proteins for protection. In this work, we resolved the native structure of a massive host-defense complex that polymerizes 30,000 guanylate-binding proteins (GBPs) over the surface of gram-negative bacteria inside human cells. Construction of this giant nanomachine took several minutes and remained stable for hours, required guanosine triphosphate hydrolysis, and recruited four GBPs plus caspase-4 and Gasdermin D as a cytokine and cell death immune signaling platform. Cryo-electron tomography suggests that GBP1 can adopt an extended conformation for bacterial membrane insertion to establish this platform, triggering lipopolysaccharide release that activated coassembled caspase-4. Our "open conformer" model provides a dynamic view into how the human GBP1 defense complex mobilizes innate immunity to infection.

Unraveling the Significance of Mg2+ Dependency and Nucleotide Binding Specificity of H-RAS.

RAS is a small GTPase and acts as a binary molecular switch; the transition from its active to inactive state plays a crucial role in various cell signaling processes. Molecular dynamics simulations at the atomistic level suggest that the absence of cofactor Mg2+ ion generally leads to pronounced structural changes in the Switch-I than Switch-II regions and assists GTP binding. The presence of the Mg2+ ion also restricts the rotation of ϒ phosphate and enhances the hydrolysis rate of GTP. Further, the simulations reveal that the stability of the protein is almost uncompromised when Mg2+ is replaced with Zn2+ and not the Ca2+ ion. The specificity of H-RAS to GTP was evaluated by substituting with ATP and CTP, which indicates that the binding pocket tolerates purine bases over pyrimidine bases. However, the D119 residue specifically interacts with the guanine base and serves as one of the primary interactions that leads to the selectivity of GTP over ATP. The ring displacement of 32Y serves as gate dynamics in H-RAS which are important for its interaction with GAP for the nucleotide exchange and is restricted in the presence of ATP. Finally, the point mutations 61, 16, and 32 influence the structural changes, specifically in the Switch-II region, which are expected to impact the GTP hydrolysis and thus are termed oncogenic mutations.

Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics.

The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.

Crystal structure of GTP-dependent dephospho-coenzyme A kinase from the hyperthermophilic archaeon, Thermococcus kodakarensis.

The biosynthesis pathways of coenzyme A (CoA) in most archaea involve several unique enzymes including dephospho-CoA kinase (DPCK) that converts dephospho-CoA to CoA in the final step of CoA biosynthesis in all domains of life. The archaeal DPCK is unrelated to the analogous bacterial and eukaryotic enzymes and shows no significant sequence similarity to any proteins with known structures. Unusually, the archaeal DPCK utilizes GTP as the phosphate donor although the analogous bacterial and eukaryotic enzymes are ATP-dependent kinases. Here, we report the crystal structure of DPCK and its complex with GTP and a magnesium ion from the archaeal hyperthermophile Thermococcus kodakarensis. The crystal structure demonstrates why GTP is the preferred substrate of this kinase. We also report the activity analyses of site-directed mutants of crucial residues determined based on sequence conservation and the crystal structure. From these results, the key residues involved in the reaction of phosphoryl transfer and the possible dephospho-CoA binding site are inferred.

Escherichia coli FtsZ molecular dynamics simulations.

Earlier molecular dynamics studies of the FtsZ protein revealed that the protein has high intrinsic flexibility which the crystal structures cannot reveal. However, the input structure in these simulation studies was based on the available crystal structure data and therefore, the effect of the C-terminal Intrinsically Disordered Region (IDR) of FtsZ could not be observed in any of these studies. Recent investigations have revealed that the C-terminal IDR is crucial for FtsZ assembly in vitro and Z ring formation in vivo. Therefore, in this study, we simulated FtsZ with the IDR. Simulations of the FtsZ monomer in different nucleotide bound forms (without nucleotide, GTP, GDP) were performed. In the conformations of FtsZ monomer with GTP, GTP binds variably with the protein. Such a variable interaction with the monomer has not been observed in any previous simulation studies of FtsZ and not observed in crystal structures. We found that central helix bends towards the C-terminal domain in the GTP bound form, hence, making way for polymerization. A nucleotide dependent shift/rotation of the C-terminal domain was observed in simulation time averaged structures.Communicated by Ramaswamy H. Sarma.

Real-time monitoring of the reaction of KRAS G12C mutant specific covalent inhibitor by in vitro and in-cell NMR spectroscopy.

KRAS mutations are major drivers of various cancers. Recently, allele-specific inhibitors of the KRAS G12C mutant were developed that covalently modify the thiol of Cys12, thereby trapping KRAS in an inactive GDP-bound state. To study the mechanism of action of the covalent inhibitors in both in vitro and intracellular environments, we used real-time NMR to simultaneously observe GTP hydrolysis and inhibitor binding. In vitro NMR experiments showed that the rate constant of ARS-853 modification is identical to that of GTP hydrolysis, indicating that GTP hydrolysis is the rate-limiting step for ARS-853 modification. In-cell NMR analysis revealed that the ARS-853 reaction proceeds significantly faster than that in vitro, reflecting acceleration of GTP hydrolysis by endogenous GTPase proteins. This study demonstrated that the KRAS covalent inhibitor is as effective in the cell as in vitro and that in-cell NMR is a valuable validation tool for assessing the pharmacological properties of the drug in the intracellular context.

Homogeneous luminescent quantitation of cellular guanosine and adenosine triphosphates (GTP and ATP) using QT-LucGTP&ATP assay.

Guanosine triphosphate (GTP) and adenosine triphosphate (ATP) are essential nucleic acid building blocks and serve as energy molecules for a wide range of cellular reactions. Cellular GTP concentration fluctuates independently of ATP and is significantly elevated in numerous cancers, contributing to malignancy. Quantitative measurement of ATP and GTP has become increasingly important to elucidate how concentration changes regulate cell function. Liquid chromatography-coupled mass spectrometry (LC-MS) and capillary electrophoresis-coupled MS (CE-MS) are powerful methods widely used for the identification and quantification of biological metabolites. However, these methods have limitations related to specialized instrumentation and expertise, low throughput, and high costs. Here, we introduce a novel quantitative method for GTP concentration monitoring (GTP-quenching resonance energy transfer (QRET)) in homogenous cellular extracts. CE-MS analysis along with pharmacological control of cellular GTP levels shows that GTP-QRET possesses high dynamic range and accuracy. Furthermore, we combined GTP-QRET with luciferase-based ATP detection, leading to a new technology, termed QT-LucGTP&ATP, enabling high-throughput compatible dual monitoring of cellular GTP and ATP in a homogenous fashion. Collectively, GTP-QRET and QT-LucGTP&ATP offer a unique, high-throughput opportunity to explore cellular energy metabolism, serving as a powerful platform for the development of novel therapeutics and extending its usability across a range of disciplines.

Deoxyguanosine-Linked Bifunctional Inhibitor of SAMHD1 dNTPase Activity and Nucleic Acid Binding.

Sterile alpha motif histidine-aspartate domain protein 1 (SAMHD1) is a deoxynucleotide triphosphohydrolase that exists in monomeric, dimeric, and tetrameric forms. It is activated by GTP binding to an A1 allosteric site on each monomer subunit, which induces dimerization, a prerequisite for dNTP-induced tetramerization. SAMHD1 is a validated drug target stemming from its inactivation of many anticancer nucleoside drugs leading to drug resistance. The enzyme also possesses a single-strand nucleic acid binding function that promotes RNA and DNA homeostasis by several mechanisms. To discover small molecule inhibitors of SAMHD1, we screened a custom ∼69 000-compound library for dNTPase inhibitors. Surprisingly, this effort yielded no viable hits and indicated that exceptional barriers for discovery of small molecule inhibitors existed. We then took a rational fragment-based inhibitor design approach using a deoxyguanosine (dG) A1 site targeting fragment. A targeted chemical library was synthesized by coupling a 5'-phosphoryl propylamine dG fragment (dGpC3NH2) to 376 carboxylic acids (RCOOH). Direct screening of the products (dGpC3NHCO-R) yielded nine initial hits, one of which (R = 3-(3'-bromo-[1,1'-biphenyl]), 5a) was investigated extensively. Amide 5a is a competitive inhibitor against GTP binding to the A1 site and induces inactive dimers that are deficient in tetramerization. Surprisingly, 5a also prevented ssDNA and ssRNA binding, demonstrating that the dNTPase and nucleic acid binding functions of SAMHD1 can be disrupted by a single small molecule. A structure of the SAMHD1-5a complex indicates that the biphenyl fragment impedes a conformational change in the C-terminal lobe that is required for tetramerization.

Exploration of natural compounds against the human mpox virus DNA-dependent RNA polymerase in silico.

BACKGROUND: Last year, the human monkeypox virus (hMPXV) emerged as an alarming threat to the community, with a detectable outbreak outside the African continent for the first time. According to The American Centers for Disease Control and Prevention (CDC), the virus is reported globally, with 86,746 confirmed cases (until April 08, 2023). DNA-dependent RNA polymerase (DdRp) is an essential protein for viral replication; hence it is a promising drug target for developing antiviral drugs against DNA viruses. Therefore, this study was conducted to search for natural compounds that could provide scaffolds for RNA polymerase inhibitors. METHODS: In this study, the DdRp structure of hMPXV was modeled and used to screen the natural compounds database (COCONUT). The virtual screening revealed 15 compounds able to tightly bind to the active site of the DdRp (binding energies less than -7.0 kcal/mol) compared to the physiological nucleotide, guanosine triphosphate (GTP). Molecular dynamics simulation was then performed on the top four hits and compared to GTP RESULTS: The results revealed the potential of four compounds (comp289, comp295, comp441, and comp449) in binding the hMPXV DdRp active site with a comparable binding affinity (-17.06 ± 2.96, -11.6 ± 5.34, -14.85 ± 2.66, and -10.79 ± 4.49 kcal/mol) with GTP (-21.03 ± 7.55 kcal/mol) CONCLUSION: These findings may also pave the way for developing new hMPXV inhibitors based on natural product scaffolds. These results need further experimental validation but promising as it was validated by unbiased all-atom MD simulations and binding free energy calculations.

Difference in Catalytic Loop Repositioning Leads to GMP Variation between Two Human GBP Homologues.

Interferon-gamma-inducible human large GTPases, hGBP1 and hGBP2, have a distinctive feature of hydrolyzing GTP to GDP and GMP through successive phosphate cleavages. In hGBP1, GMP is the major product, which is essential for its anti-pathogenic activities. However, its close homologue hGBP2 produces significantly less GMP, despite having a similar active site architecture. The molecular basis for less GMP formation and catalytic residue(s) in hGBP2 are not fully explored. To address these issues, we performed systematic biochemical, biophysical, and microsecond simulation studies. Our data suggest that the less GMP formation in hGBP2 is due to the lack of H-bond formation between the W79 side-chain (located near the active site) and main-chain carbonyl of K76 (present in the catalytic loop) in the substrate-bound hGBP2. The absence of this H-bond could not redirect the catalytic loop toward the beta phosphate after the cleavage of gamma-phosphate, a step essential for enhanced GMP formation. Furthermore, based on the mutational and structural analyses, this study for the first time indicates that the same residue, T75, mediates both phosphate cleavages in hGBP2 and hGBP1. This suggests the conservation of the catalytic residue in hGBP homologues. These findings emphasize the indispensable role of correct catalytic loop repositioning for efficient beta phosphate cleavage. This led us to propose a new substrate hydrolysis mechanism by hGBP1 and hGBP2, which may also be helpful to understand the GTP hydrolysis in other hGBP homologues. Overall, the study could provide insight into how these two close homologues play crucial roles in host-mediated immunity through different mechanisms.