Coupling of retinal, protein, and water dynamics in squid rhodopsin

The light-induced isomerization of the retinal from 11-cis to all-trans triggers changes in the conformation of visual rhodopsins that lead to the formation of the activated state, which is ready to interact with the G protein. To begin to understand how changes in the structure and dynamics of the retinal are transmitted to the protein, we performed molecular dynamics simulations of squid rhodopsin with 11-cis and all-trans retinal, and with two different force fields for describing the retinal molecule. The results indicate that structural rearrangements in the binding pocket, albeit small, propagate toward the cytoplasmic side of the protein, and affect the dynamics of internal water molecules. The sensitivity of the active-site interactions on the retinal force-field parameters highlights the coupling between the retinal molecule and its immediate protein environment.     

Early steps of the intramolecular signal transduction in rhodopsin explored by molecular dynamics simulations

We present molecular dynamics simulations of bovine rhodopsin in a membrane mimetic environment based on the recently refined X-ray structure of the pigment. The interactions between the protonated Schiff base and the protein moiety are explored both with the chromophore in the dark-adapted 11-cis and in the photoisomerized all-trans form. Comparison of simulations with Glu181 in different protonation states strongly suggests that this loop residue located close to the 11-cis bond bears a negative charge. Restrained molecular dynamics simulations also provide evidence that the protein tightly confines the absolute conformation of the retinal around the C12-C13 bond to a positive helicity. 11-cis to all-trans isomerization leads to an internally strained chromophore, which relaxes after a few nanoseconds by a switching of the ionone ring to an essentially planar all-trans conformation. This structural transition of the retinal induces in turn significant conformational changes of the protein backbone, especially in helix VI. Our results suggest a possible molecular mechanism for the early steps of intramolecular signal transduction in a prototypical G-protein-coupled receptor.     

Molecular dynamics simulations of the Apo-, Holo-, and acyl-forms of Escherichia coli acyl carrier protein

Acyl carrier protein (ACP) is an essential co-factor protein in fatty acid biosynthesis that shuttles covalently bound fatty acyl intermediates in its hydrophobic pocket to various enzyme partners. To characterize acyl chain-ACP interactions and their influence on enzyme interactions, we performed 19 molecular dynamics (MD) simulations of Escherichia coli apo-, holo-, and acyl-ACPs. The simulations were started with the acyl chain in either a solvent-exposed or a buried conformation. All four short-chain (< or = C10) and one long-chain (C16) unbiased acyl-ACP MD simulation show the transition of the solvent-exposed acyl chain into the hydrophobic pocket of ACP, revealing its pathway of acyl chain binding. Although the acyl chain resides inside the pocket, Thr-39 and Glu-60 at the entrance stabilize the phosphopantetheine linker through hydrogen bonding. Comparisons of the different ACP forms indicate that the loop region between helices II and III and the prosthetic linker may aid in substrate recognition by enzymes of fatty acid synthase systems. The MD simulations consistently show that the hydrophobic binding pocket of ACP is best suited to accommodate an octanoyl group and is capable of adjusting in size to accommodate chain lengths as long as decanoic acid. The simulations also reveal a second, novel binding mode of the acyl chains inside the hydrophobic binding pocket directed toward helix I. This study provides a detailed dynamic picture of acyl-ACPs that is in excellent agreement with available experimental data and, thereby, provides a new understanding of enzyme-ACP interactions.     

Retinal dynamics during light activation of rhodopsin revealed by solid-state NMR spectroscopy

Rhodopsin is a canonical member of class A of the G protein-coupled receptors (GPCRs) that are implicated in many of the drug interventions in humans and are of great pharmaceutical interest. The molecular mechanism of rhodopsin activation remains unknown as atomistic structural information for the active metarhodopsin II state is currently lacking. Solid-state ²H NMR constitutes a powerful approach to study atomic-level dynamics of membrane proteins. In the present application, we describe how information is obtained about interactions of the retinal cofactor with rhodopsin that change with light activation of the photoreceptor. The retinal methyl groups play an important role in rhodopsin function by directing conformational changes upon transition into the active state. Site-specific ²H labels have been introduced into the methyl groups of retinal and solid-state ²H NMR methods applied to obtain order parameters and correlation times that quantify the mobility of the cofactor in the inactive dark state, as well as the cryotrapped metarhodopsin I and metarhodopsin II states. Analysis of the angular-dependent ²H NMR line shapes for selectively deuterated methyl groups of rhodopsin in aligned membranes enables determination of the average ligand conformation within the binding pocket. The relaxation data suggest that the β-ionone ring is not expelled from its hydrophobic pocket in the transition from the pre-activated metarhodopsin I to the active metarhodopsin II state. Rather, the major structural changes of the retinal cofactor occur already at the metarhodopsin I state in the activation process. The metarhodopsin I to metarhodopsin II transition involves mainly conformational changes of the protein within the membrane lipid bilayer rather than the ligand. The dynamics of the retinylidene methyl groups upon isomerization are explained by an activation mechanism involving cooperative rearrangements of extracellular loop E2 together with transmembrane helices H5 and H6. These activating movements are triggered by steric clashes of the isomerized all-trans retinal with the β4 strand of the E2 loop and the side chains of Glu¹²² and Trp²⁶⁵ within the binding pocket. The solid-state ²H NMR data are discussed with regard to the pathway of the energy flow in the receptor activation mechanism.     

Retinal Flip in Rhodopsin Activation?

Rhodopsin is a well-characterized structural model of a G protein-coupled receptor. Photoisomerization of the covalently bound retinal triggers activation. Surprisingly, the x-ray crystal structure of the active Meta-II state has a 180° rotation about the long-axis of the retinal polyene chain. Unbiased microsecond-timescale all-atom molecular dynamics simulations show that the retinal cofactor can flip back to the orientation observed in the inactive state of rhodopsin under conditions favoring the Meta-I state. Our results provide, to our knowledge, the first evidence from molecular dynamics simulations showing how rotation of the retinal ligand within its binding pocket can occur in the activation mechanism of rhodopsin.     

Nanoseconds molecular dynamics simulation of primary mechanical energy transfer steps in F1-ATP synthase

The mitochondrial membrane protein FoF1-ATP synthase synthesizes adenosine triphosphate (ATP), the universal currency of energy in the cell. This process involves mechanochemical energy transfer from a rotating asymmetric gamma-'stalk' to the three active sites of the F1 unit, which drives the bound ATP out of the binding pocket. Here, the primary structural changes associated with this energy transfer in F1-ATP synthase were studied with multi-nanosecond molecular dynamics simulations. By forced rotation of the gamma-stalk that mimics the effect of proton motive Fo-rotation during ATP synthesis, a time-resolved atomic model for the structural changes in the F1 part in terms of propagating conformational motions is obtained. For these, different time scales are found, which allows the separation of nanosecond from microsecond conformational motions. In the simulations, rotation of the gamma-stalk lowers the ATP affinity of the betaTP binding pocket and triggers fast, spontaneous closure of the empty betaE subunit. The simulations explain several mutation studies and the reduced hydrolysis rate of gamma-depleted F1-ATPase.     

Modeling of Arylamide Helix Mimetics in the p53 Peptide Binding Site of hDM2 Suggests Parallel and Anti-Parallel Conformations Are Both Stable

The design of novel α-helix mimetic inhibitors of protein-protein interactions is of interest to pharmaceuticals and chemical genetics researchers as these inhibitors provide a chemical scaffold presenting side chains in the same geometry as an α-helix. This conformational arrangement allows the design of high affinity inhibitors mimicking known peptide sequences binding specific protein substrates. We show that GAFF and AutoDock potentials do not properly capture the conformational preferences of α-helix mimetics based on arylamide oligomers and identify alternate parameters matching solution NMR data and suitable for molecular dynamics simulation of arylamide compounds. Results from both docking and molecular dynamics simulations are consistent with the arylamides binding in the p53 peptide binding pocket. Simulations of arylamides in the p53 binding pocket of hDM2 are consistent with binding, exhibiting similar structural dynamics in the pocket as simulations of known hDM2 binders Nutlin-2 and a benzodiazepinedione compound. Arylamide conformations converge towards the same region of the binding pocket on the 20 ns time scale, and most, though not all dihedrals in the binding pocket are well sampled on this timescale. We show that there are two putative classes of binding modes for arylamide compounds supported equally by the modeling evidence. In the first, the arylamide compound lies parallel to the observed p53 helix. In the second class, not previously identified or proposed, the arylamide compound lies anti-parallel to the p53 helix.     

Molecular dynamics simulations of retinal in rhodopsin: from the dark-adapted state towards lumirhodopsin

The formation of photointermediates and conformational changes observed in the retinal chromophore of bilayer-embedded rhodopsin during the early steps of the protein activation have been studied by molecular dynamics (MD) simulation. In particular, the lysine-bound retinal has been examined, focusing on its conformation in the dark-adapted state (10 ns) and on the early steps after the isomerization of the 11-cis bond to trans (up to 10 ns). The parametrization for the chromophore is based on a recent quantum study [Sugihara, M., Buss, V., Entel, P., Elstner, M., and Frauenheim, T. (2002) Biochemistry 41, 15259-15266] and shows good conformational agreement with recent experimental results. The isomerization, induced by switching the function governing the dihedral angle for the C11=C12 bond, was repeated with several different starting conformations. From the repeated simulations, it is shown that the retinal model exhibits a conserved activation pattern. The conformational changes are sequential and propagate outward from the C11=C12 bond, starting with isomerization of the C11=C12 bond, then a rotation of methyl group C20, and followed by increased fluctuations at the beta-ionone ring. The dynamics of these changes suggest that they are linked with photointermediates observed by spectroscopy. The exact moment when these events occur after the isomerization is modulated by the starting conformation, suggesting that retinal isomerizes through multiple pathways that are slightly different. The amplitudes of the structural fluctuations observed for the protein in the dark-adapted state and after isomerization of the retinal are similar, suggesting a subtle mechanism for the transmission of information from the chromophore to the protein.     

Exploration of the conformational space of myosin recovery stroke via molecular dynamics

Muscle contractions are driven by cyclic conformational changes of myosin, whose molecular mechanisms of operation are being elucidated by recent advances in crystallographic studies and single molecule experiments. To complement such structural studies and consider the energetics of the conformational changes of myosin head, umbrella sampling molecular dynamics (MD) simulations were performed with the all-atom model of the scallop myosin sub-fragment 1 (S1) with a bound ATP in solution in explicit water using the crystallographic near-rigor and transition state conformations as two references. The constraints on RMSD reaction coordinates used for the umbrella sampling were found to steer the conformational changes efficiently, and relatively close correlations have been observed between the set of characteristic structural changes including the lever arm rotation and the closing of the nucleotide binding pocket. The lever arm angle and key residue interaction distances in the nucleotide binding pocket and the relay helix show gradual changes along the recovery stroke reaction coordinate, consistent with previous crystallographic and computational minimum energy studies. Thermal fluctuations, however, appear to make the switch-2 coordination of ATP more flexible than suggested by crystal structures. The local solvation environment of the fluorescence probe, Trp 507 (scallop numbering), also appears highly mobile in the presence of thermal fluctuations.     

Molecular dynamics simulation of dark-adapted rhodopsin and free opsin

Molecular dynamics simulation was carried out for the rhodopsin protein to investigate its conformational changes in respect to inclusion of 11-cis retinal chromophore. Molecular dynamics calculations were performed within the time frame 3000 ps. Totally, 3 X 10(6) configurations ofrhodopsin and free opsin were analyzed and compared. It has been shown that the 11-cis retinal rearrangement (adaptation) in opsin strongly affects the surrounding amino acid residues of protein binding pocket and the protein cytoplasmic region. The extracellular part, however, shows comparatively little changes. On basis of the simulation results obtained we propose a molecular mechanism for the rhodopsin protein function as a G-protein-coupled receptor in the state of darkness. We discuss the role of the retinal chromophore as a ligand-antagonist stabilizing the inactive conformation of rhodopsin.     

Spontaneous conformational changes in the E. coli GroEL subunit from all-atom molecular dynamics simulations

The Escherichia coli chaperonin GroEL is a complex of identical subunit proteins (57 kDa each) arranged in a back-to-back stacking of two heptameric rings. Its hallmarks include nested positive intra-ring and negative inter-ring cooperativity in adenosine trisphosphate (ATP) binding and the ability to mediate the folding of newly transcribed and/or denatured substrate proteins. We performed unbiased molecular dynamics simulations of the GroEL subunit protein in explicit water both with and without the nucleotide KMgATP to understand better the details of the structural transitions that enable these behaviors. Placing KMgATP in the equatorial domain binding pocket of a t state subunit, which corresponds to a low ATP-affinity state, produced a short-lived (6 ns) state that spontaneously transitioned to the high ATP-affinity r state. The important feature of this transition is a large-scale rotation of the intermediate domain's helix M to close the ATP binding pocket. Pivoting of helix M is accompanied by counterclockwise rotation and slight deformation of the apical domain, important for lowering the affinity for substrate protein. Aligning simulation conformations into model heptamer rings demonstrates that the t-->r transition in one subunit is not sterically hindered by t state neighbors, but requires breakage of Arg(197)-Glu(386) intersubunit salt bridges, which are important for inter-ring positive cooperativity. Lowest-frequency quasi-harmonic modes of vibration computed pre- and post-transition clearly show that natural vibrations facilitate the transition. Finally, we propose a novel mechanism for inter-ring cooperativity in ATP binding inspired by the observation of spontaneous insertion of the side chain of Ala(480) into the empty nucleotide pocket.     

All-atom structural investigation of kinesin-microtubule complex constrained by high-quality cryo-electron-microscopy maps

In this study, we have performed a comprehensive structural investigation of three major biochemical states of a kinesin complexed with microtubule under the constraint of high-quality cryo-electron-microscopy (EM) maps. In addition to the ADP and ATP state which were captured by X-ray crystallography, we have also modeled the nucleotide-free or APO state for which no crystal structure is available. We have combined flexible fitting of EM maps with regular molecular dynamics simulations, hydrogen-bond analysis, and free energy calculation. Our APO-state models feature a subdomain rotation involving loop L2 and α6 helix of kinesin, and local structural changes in active site similar to a related motor protein, myosin. We have identified a list of hydrogen bonds involving key residues in the active site and the binding interface between kinesin and microtubule. Some of these hydrogen bonds may play an important role in coupling microtubule binding to ATPase activities in kinesin. We have validated our models by calculating the binding free energy between kinesin and microtubule, which quantitatively accounts for the observation of strong binding in the APO and ATP state and weak binding in the ADP state. This study will offer promising targets for future mutational and functional studies to investigate the mechanism of kinesin motors.     

Extended molecular dynamics simulation of the carbon monoxide migration in sperm whale myoglobin

We report the results of an extended molecular dynamics simulation on the migration of photodissociated carbon monoxide in wild-type sperm whale myoglobin. Our results allow following one possible ligand migration dynamics from the distal pocket to the Xe1 cavity via a path involving the other xenon binding cavities and momentarily two additional packing defects along the pathway. Comparison with recent time resolved structural data obtained by Laue crystallography with subnanosecond to millisecond resolution shows a more than satisfactory agreement. In fact, according to time resolved crystallography, CO, after photolysis, can occupy the Xe1 and Xe4 cavities. However, no information on the trajectory of the ligand from the distal pocket to the Xe1 is available. Our results clearly show one possible path within the protein. In addition, although our data refer to a single trajectory, the local dynamics of the ligand in each cavity is sufficiently equilibrated to obtain local structural and thermodynamic information not accessible to crystallography. In particular, we show that the CO motion and the protein fluctuations are strictly correlated: free energy calculations of the migration between adjacent cavities show that the migration is not a simple diffusion but is kinetically or thermodynamically driven by the collective motions of the protein; conversely, the protein fluctuations are influenced by the ligand in such a way that the opening/closure of the passage between adjacent cavities is strictly correlated to the presence of CO in its proximity. The compatibility between time resolved crystallographic experiments and molecular dynamics simulations paves the way to a deeper understanding of the role of internal dynamics and packing defects in the control of ligand binding in heme proteins.     

Modulation of TRPV2 by endogenous and exogenous ligands: A computational study

Transient receptor potential vanilloid (TRPV) channels play various important roles in human physiology. As membrane proteins, these channels are modulated by their endogenous lipid environment as the recent wealth of structural studies has revealed functional and structural lipid binding sites. Additionally, it has been shown that exogenous ligands can exchange with some of these lipids to alter channel gating. Here, we used molecular dynamics simulations to examine how one member of the TRPV family, TRPV2, interacts with endogenous lipids and the pharmacological modulator cannabidiol (CBD). By computationally reconstituting TRPV2 into a typical plasma membrane environment, which includes phospholipids, cholesterol, and phosphatidylinositol (PIP) in the inner leaflet, we showed that most of the interacting surface lipids are phospholipids without strong specificity for headgroup types. Intriguingly, we observed that the C-terminal membrane proximal region of the channel binds preferentially to PIP lipids. We also modelled two structural lipids in the simulation: one in the vanilloid pocket and the other in the voltage sensor-like domain (VSLD) pocket. The simulation shows that the VSLD lipid dampens the fluctuation of the VSLD residues, while the vanilloid lipid exhibits heterogeneity both in its binding pose and in its influence on protein dynamics. Addition of CBD to our simulation system led to an open selectivity filter and a structural rearrangement that includes a clockwise rotation of the ankyrin repeat domains, TRP helix, and VSLD. Together, these results reveal the interplay between endogenous lipids and an exogenous ligand and their effect on TRPV2 stability and channel gating.     

Thermodynamics of nucleotide binding to actomyosin V and VI: a positive heat capacity change accompanies strong ADP binding

We have measured the energetics of ATP and ADP binding to single-headed actomyosin V and VI from the temperature dependence of the rate and equilibrium binding constants. Nucleotide binding to actomyosin V and VI can be modeled as two-step binding mechanisms involving the formation of collision complexes followed by isomerization to states with high nucleotide affinity. Formation of the actomyosin VI-ATP collision complex is much weaker and slower than for actomyosin V. A three-step binding mechanism where actomyosin VI isomerizes between two conformations, one competent to bind ATP and one not, followed by rapid ATP binding best accounts for the data. ADP binds to actomyosin V more tightly than actomyosin VI. At 25 degrees C, the strong ADP-binding equilibria are comparable for actomyosin V and VI, and the different overall ADP affinities arise from differences in the ADP collision complex affinity. The actomyosin-ADP isomerization leading to strong ADP binding is entropy driven at >15 degrees C and occurs with a large, positive change in heat capacity (DeltaC(P) degrees ) for both actomyosin V and VI. Sucrose slows ADP binding and dissociation from actomyosin V and VI but not the overall equilibrium constants for strong ADP binding, indicating that solvent viscosity dampens ADP-dependent kinetic transitions, presumably a tail swing that occurs with ADP binding and release. We favor a mechanism where strong ADP binding increases the dynamics and flexibility of the actomyosin complex. The heat capacity (DeltaC(P) degrees ) and entropy (DeltaS degrees ) changes are greater for actomyosin VI than actomyosin V, suggesting different extents of ADP-induced structural rearrangement.     

A role for transmembrane domains V and VI in ligand binding and maturation of the angiotensin II AT1 receptor

Several studies have proposed that angiotensin II (Ang II) binds to its receptor AT1 through interactions with residues in helices V and VI, suggesting that the distance between these helices is crucial for ligand binding. Based on a 3D model of AT1 in which the C-terminus of Ang II is docked, we identified the hydrophobic residues of TM V and VI pointing towards the external face of the helices, which may play a role in the structure of the binding pocket and in the structural integrity of the receptor. We performed a systematic mutagenesis study of these residues and examined the binding, localization, maturation, and dimerization of the mutated receptors. We found that mutations of hydrophobic residues to alanine in helix V do not alter binding, whereas mutations to glutamate lead to loss of binding without a loss in cell surface expression, suggesting that the external face of helix V may not directly participate in binding, but may rather contribute to the structure of the binding pocket. In contrast, mutations of hydrophobic residues to glutamate in helix VI lead to a loss in cell surface expression, suggesting that the external surface of helix VI plays a structural role and ensures correct folding of the receptor.     

The signaling pathway of rhodopsin

The signal-transduction mechanism of rhodopsin was studied by molecular dynamics (MD) simulations of the high-resolution, inactive structure in an explicit membrane environment. The simulations were employed to calculate equal-time correlations of the fluctuating interaction energy of residue pairs. The resulting interaction-correlation matrix was used to determine a network that couples retinal to the cytoplasmic interface, where transducin binds. Two highly conserved motifs, D(E)RY and NPxxY, were found to have strong interaction correlation with retinal. MD simulations with restraints on each transmembrane helix indicated that the major signal-transduction pathway involves the interdigitating side chains of helices VI and VII. The functional roles of specific residues were elucidated by the calculated effect of retinal isomerization from 11-cis to all-trans on the residue-residue interaction pattern. It is suggested that Glu134 may act as a "signal amplifier" and that Asp83 may introduce a threshold to prevent background noise from activating rhodopsin.     

Kinetic folding mechanism of an integral membrane protein examined by pulsed oxidative labeling and mass spectrometry

We report the application of pulsed oxidative labeling for deciphering the folding mechanism of a membrane protein. SDS-denatured bacteriorhodopsin (BR) was refolded by mixing with bicelles in the presence of free retinal. At various time points (20 ms to 1 day), the protein was exposed to a microsecond ·OH pulse that induces oxidative modifications at solvent-accessible methionine side chains. The extent of labeling was determined by mass spectrometry. These measurements were complemented by stopped-flow spectroscopy. Major time-dependent changes in solvent accessibility were detected for M20 (helix A) and M118 (helix D). Our kinetic data indicate a sequential folding mechanism, consistent with models previously suggested by others on the basis of optical data. Yet, ·OH labeling provides additional structural insights. An initial folding intermediate I(1) gets populated within 20 ms, concomitantly with formation of helix A. Subsequent structural consolidation leads to a transient species I(2). Noncovalent retinal binding to I(2) induces folding of helix D, thereby generating an intermediate I(R). In the absence of retinal, the latter transition does not take place. Hence, formation of helix D depends on retinal binding, whereas this is not the case for helix A. As the cofactor settles deeper into its binding pocket, a final transient species I(R) is generated. This intermediate converts into native BR within minutes by formation of the retinal-K216 Schiff base linkage. The combination of pulsed covalent labeling and optical spectroscopy employed here should also be suitable for exploring the folding mechanisms of other membrane proteins.     

Molecular dynamics study of the nature and origin of retinal's twisted structure in bacteriorhodopsin

The planarity of the polyene chain of the retinal chromophore in bacteriorhodopsin is studied using molecular dynamics simulation techniques and applying different force-field parameters and starting crystal structures. The largest deviations from a planar structure are observed for the C(13)==C(14) and C(15)==N(16) double bonds in the retinal Schiff base structure. The other dihedral angles along the polyene chain of the chromophore, although having lower torsional barriers in some cases, do not significantly deviate from the planar structure. The results of the simulations of different mutants of the pigment show that, among the studied amino acids of the binding pocket, the side chain of Trp-86 has the largest impact on the planarity of retinal, and the mutation of this amino acid to alanine leads to chromophore planarity. Deletion of the methyl C(20), removal of a water molecule hydrogen-bonded to H(15), or mutation of other amino acids to alanine did not show any significant influence on the distortion of the chromophore. The results from the present study suggest the importance of the bulky residue of Trp-86 in the isomerization process, in both ground and excited states of the chromophore, and in fine-tuning of the pK(a) of the retinal protonated Schiff base in bacteriorhodopsin. The dark adaptation of the pigment and the last step of the bacteriorhodopsin photocycle imply low barriers against the rotation of the double bonds in the Schiff base region. The twisted double bonds found in the present study are consistent with the proposed mechanism of these ground state isomerization events.     

Influence of GTP/GDP and magnesium ion on the solvated structure of the protein FtsZ: a molecular dynamics study

We present here a structural analysis of ten extensive all-atom molecular dynamics simulations of the monomeric protein FtsZ in various binding states. Since the polymerization and GTPase activities of FtsZ depend on the nature of a bound nucleotide as well as on the presence of a magnesium ion, we studied the structural differences between the average conformations of the following five systems: FtsZ-Apo, FtsZ-GTP, FtsZ-GDP, FtsZ-GTP-Mg, and FtsZ-GDP-Mg. The in silico solvated average structure of FtsZ-Apo significantly differs from the crystallographic structure 1W59 of FtsZ which was crystallized in a dimeric form without nucleotide and magnesium. The simulated Apo form of the protein also clearly differs from the FtsZ structures when it is bound to its ligand, the most important discrepancies being located in the loops surrounding the nucleotide binding pocket. The three average structures of FtsZ-GTP, FtsZ-GDP, and FtsZ-GTP-Mg are overall similar, except for the loop T7 located at the opposite side of the binding pocket and whose conformation in FtsZ-GDP notably differs from the one in FtsZ-GTP and FtsZ-GTP-Mg. The presence of a magnesium ion in the binding pocket has no impact on the FtsZ conformation when it is bound to GTP. In contrast, when the protein is bound to GDP, the divalent cation causes a translation of the nucleotide outwards the pocket, inducing a significant conformational change of the loop H6-H7 and the top of helix H7.