Reaction Mechanism of Guanosine Triphosphate Hydrolysis by the Vision-Related Protein Complex Arl3-RP2

Molecular dynamics simulations for the structure-based drug design: targeting small-GTPases proteins

INTRODUCTION

Molecular Dynamics (MD) simulations can support mechanism-based drug design. Indeed, MD simulations by capturing biomolecule motions at finite temperatures can reveal hidden binding sites, accurately predict drug-binding poses, and estimate the thermodynamics and kinetics, crucial information for drug discovery campaigns. Small-Guanosine Triphosphate Phosphohydrolases (GTPases) regulate a cascade of signaling events, that affect most cellular processes. Their deregulation is linked to several diseases, making them appealing drug targets. The broad roles of small-GTPases in cellular processes and the recent approval of a covalent KRas inhibitor as an anticancer agent renewed the interest in targeting small-GTPase with small molecules.

AREA COVERED

This review emphasizes the role of MD simulations in elucidating small-GTPase mechanisms, assessing the impact of cancer-related variants, and discovering novel inhibitors.

EXPERT OPINION

The application of MD simulations to small-GTPases exemplifies the role of MD simulations in the structure-based drug design process for challenging biomolecular targets. Furthermore, AI and machine learning-enhanced MD simulations, coupled with the upcoming power of quantum computing, are promising instruments to target elusive small-GTPases mutations and splice variants. This powerful synergy will aid in developing innovative therapeutic strategies associated to small-GTPases deregulation, which could potentially be used for personalized therapies and in a tissue-agnostic manner to treat tumors with mutations in small-GTPases.

Mechanism of Guanosine Triphosphate Hydrolysis by the Visual Proteins Arl3-RP2: Free Energy Reaction Profiles Computed with Ab Initio Type QM/MM Potentials

We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme-substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.

Amide-imide tautomerization in the glutamine side chain in enzymatic and photochemical reactions in proteins

Amide-imide tautomerization presents a pervasive class of chemical transformations in organic chemistry of natural compounds. In this Perspective, we describe two distinctively different protein systems, in which the amide-imide tautomerization in the glutamine side chain takes place in enzymatic or photochemical reactions. First, hydrolysis of guanosine triphosphate (GTP) catalyzed by the Ras-GAP protein complex suggests the occurrence of the imide tautomer of glutamine in reaction intermediates. Second, photoexcitation of flavin-binding protein domains (BLUFs) initiates a chain of reactions in the chromophore-binding pocket, including amide-imide tautomerization of glutamine. Mechanisms of these reactions at the atomic level have been revealed in quantum mechanics/molecular mechanics (QM/MM) simulations. To reinforce conclusions on the critical role of amide-imide tautomerization of glutamine in these reactions we describe results of new quantum chemistry and QM/MM calculations for relevant molecular model systems. We reexamine results of the recent IR spectroscopy studies of BLUF domains, which provide experimental evidences of Gln tautomerization in proteins. We also propose to validate the glutamine-assisted mechanism of enzymatic GTP hydrolysis by using IR spectroscopy in a proper range of wavenumbers.

Elucidation of Single Hydrogen Bonds in GTPases via Experimental and Theoretical Infrared Spectroscopy

Time-resolved Fourier transform infrared (FTIR) spectroscopy is a powerful tool to elucidate label-free the reaction mechanisms of proteins. After assignment of the absorption bands to individual groups of the protein, the order of events during the reaction mechanism can be monitored and rate constants can be obtained. Additionally, structural information is encoded into infrared spectra and can be decoded by combining the experimental data with biomolecular simulations. We have determined recently the infrared vibrations of GTP and guanosine diphosphate (GDP) bound to Gα_i1, a ubiquitous GTPase. These vibrations are highly sensitive for the environment of the phosphate groups and thereby for the binding mode the GTPase adopts to enable fast hydrolysis of GTP. In this study we calculated these infrared vibrations from biomolecular simulations to transfer the spectral information into a computational model that provides structural information far beyond crystal structure resolution. Conformational ensembles were generated using 15 snapshots of several 100 ns molecular-mechanics/molecular-dynamics (MM-MD) simulations, followed by quantum-mechanics/molecular-mechanics (QM/MM) minimization and normal mode analysis. In comparison with other approaches, no time-consuming QM/MM-MD simulation was necessary. We carefully benchmarked the simulation systems by deletion of single hydrogen bonds between the GTPase and GTP through several Gα_i1 point mutants. The missing hydrogen bonds lead to blue-shifts of the corresponding absorption bands. These band shifts for α-GTP (Gα_i1-T48A), γ-GTP (Gα_i1-R178S), and for both β-GTP/γ-GTP (Gα_i1-K46A, Gα_i1-D200E) were found in agreement in the experimental and the theoretical spectra. We applied our approach to open questions regarding Gα_i1: we show that the GDP state of Gα_i1 carries a Mg²⁺, which is not found in x-ray structures. Further, the catalytic role of K46, a central residue of the P-loop, and the protonation state of the GTP are elucidated.

How to gather useful and valuable information from protein binding measurements using Langmuir lipid monolayers

This review presents data on the influence of various experimental parameters on the binding of proteins onto Langmuir lipid monolayers. The users of the Langmuir methodology are often unaware of the importance of choosing appropriate experimental conditions to validate the data acquired with this method. The protein Retinitis pigmentosa 2 (RP2) has been used throughout this review to illustrate the influence of these experimental parameters on the data gathered with Langmuir monolayers. The methods detailed in this review include the determination of protein binding parameters from the measurement of adsorption isotherms, infrared spectra of the protein in solution and in monolayers, ellipsometric isotherms and fluorescence micrographs.