Umbrella Sampling and X-ray Crystallographic Analysis Unveil an Arg

Jan 18, 2018 - Phosphopantetheine adenylyltransferase (PPAT) is a rate-limiting enzyme essential for biosynthesis of coenzyme A (CoA), which in turn i...
0 downloads 4 Views 3MB Size
Subscriber access provided by Universitaetsbibliothek | Johann Christian Senckenberg

Article

Umbrella Sampling and X-ray Crystallographic Analysis Unveil an Arg-Asp Gate Facilitating Inhibitor Binding Inside Phosphopantetheine Adenylyltransferase Allosteric Cleft Abhisek Mondal, Rakesh Chatterjee, and Saumen Datta J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.7b09543 • Publication Date (Web): 18 Jan 2018 Downloaded from http://pubs.acs.org on January 18, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Title: Umbrella Sampling and X-ray Crystallographic Analysis Unveil an Arg-Asp Gate Facilitating Inhibitor Binding Inside Phosphopantetheine Adenylyltransferase Allosteric Cleft

A. Mondal, R. Chatterjee and S. Datta*

Structural Biology and Bioinformatics Division, Council of Scientific and Industrial Research-Indian Institute of Chemical Biology, 4 Raja SC Mullick Road, Jadavpur, Kolkata, West Bengal, India. *

Corresponding Author: Email- [email protected] Telephone Number: +913324995896

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Phosphopantetheine adenylyltransferase (PPAT) is a rate limiting enzyme essential for biosynthesis of coenzyme A (CoA) which in turn is responsible to regulate the secretion of exotoxins via type III secretion system in Pseudomonas aeruginosa, causing severe health concerns ranging from nosocomial infections to respiratory failure. Acetyl Coenzyme A (AcCoA) is a newly reported inhibitor of PPAT, believed to regulate the cellular levels of CoA and thereby the pathogenesis. Very little is known so far regarding the mechanistic details of AcCoA binding inside PPAT binding cleft. Herein, we have used extensive umbrella sampling simulations to decipher mechanistic insight regarding the inhibitor accommodation inside the binding cavity. We found out that R90 and D94 residues act like a gate near the binding cavity to accommodate and stabilize the incoming ligand. Mutational models concerning these residues also show considerable difference in AcCoA binding thermodynamics. To substantiate our findings, we have solved the first crystal structure of apo-PPAT from P. aeruginosa which also found to agree with the simulation results. Collectively, these results describe the mechanistic details of accommodation of inhibitor molecule inside PPAT binding cavity and also offer valuable insight for regulating cellular levels of CoA/AcCoA and thus controlling the pathogenicity.

2 ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction: Pseudomonas aeruginosa is a frequent cause of nosocomial pneumonia, having mortality rate of ≤ 60%1. This gram-negative pathogen is a predominant cause of nosocomial bacteremia among immune compromised individuals2. Clinical reports show that chronic colonizing of these bacteria causes eventual respiratory failure and even death in cystic fibrosis patients2-3. World Health Organization published a report on Antimicrobial Resistance in 2014, saying increasing number of multidrug resistant bacteria continues to complicate medical treatment of infected individuals per year. P. aeruginosa infections in immune compromised patients primarily conserved to its own toxin repertoire, in spite of that increasing antibiotic resistance renders the treatment more difficult4-5. As suggested by Osterman and Begley, structural and mechanistic analysis of different conserve pathways, responsible for survival and proliferation of pathogens, would help us to provide new insights in order to develop novel drugs6-7. Pantothenate and Coenzyme A (CoA) biosynthetic pathway which is regulated by pantothenate kinase and phosphopantetheine adenylyltransferase (PPAT) fulfills these criteria. In major metabolic pathways CoA functions as an acyl carrier and also activates carbonyl groups in many reactions. It is also found to be involved in fatty acid biosynthesis, lipopolysaccharide biosynthesis, Krebs cycle6, 8. The biosynthetic precursor of CoA is pantothenate. However, the source of pantothenate varies in different organisms. It is produced de novo in some bacteria and other bacteria, like P. aeruginosa uses extracellular pantetheine for CoA biosynthesis in five different steps. Phosphopantetheine adenylyltransferase (PPAT) catalyzes the rate limiting step to reversibly transfer the adenyl group from adenosine tri-phosphate (ATP) to 4’phosphopanteheine (PhP) to produce 3’-dephosphocoenzyme A (dPCoA) and pyrophosphate (PPi). dPCoA-kinase further converts the dPCoA to produce CoA9. This second rate limiting step in CoA biosynthesis pathway is catalyzed by PPAT and also regulates the cellular CoA content, simultaneously with the catalytic protein pantothenate kinase10. Various virulence pathways are tightly regulated by metabolic pathways which act in a synchronized manner to ensure pathogenesis in correct environmental and cellular condition. Prior studies emphasized that CoA is utilized by Pyruvate Dehydrogenase Complex (PDC) to produce Acetyl-CoA which in turn regulates malate synthase (glcB) and isocitrate lyase (aceA) and thereby modulating the PcrG and PcrV (crucial exotoxin and chaperone for infection) secretion in Type III secretion system1112 . PPAT resembles the characteristic feature of alpha/beta phosphodiesterase superfamily in nucleotide binding mediated by 6 β Rossman fold. This family bears H/TXGH motifs13 which help to stabilize the pentacoordinate transition state during the transition from allosteric to substrate binding state and vice versa. PPAT is found to follow this mechanism and possesses a dinucleotide binding groove consistent with class I aminoacyl-tRNA synthetases which binds with ATP and undergoes 4’-phosphopantehteine mediated nucleophilic attack in an in-line displacement mechanism14-15. Many PPAT co-crystal structures, bound with different ligands, have been solved till date to understand the mechanism of interactions of various substrates or 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ligands with this enzyme14-22. Recent X-ray crystallographic studies of PPAT from P. aeruginosa has reported a new allosteric molecule (AcCoA), which binds in the CoA binding groove and thereby inhibits the enzyme action. Structural characterization revealed that R90 and D94 residues play a crucial role in AcCoA binding in the PPAT allosteric site23. It is proposed earlier that these two residues show a gate or bridge like phenomena using ionic interactions that might allow the AcCoA binding23. Prior structural studies with apo form of M. tuberculosis showed that the ligand accommodation to the binding groove is solely mediated by movement of a loop region. Although in a more recent study, using a series of X-ray crystallographic structures, Chatterjee et al 23 has reported that in P. aeruginosa the mechanism of ligand accommodation happens in a completely different manner. The later study hypothesized that, unlike other structures from different organisms, in P. aeruginosa a ‘gate-like mechanism’ actively controls the entry and exit of ligands inside binding cleft. However, the exact mechanism or existence of such highly choreographed molecular movement still remains poorly understood. In this study, we have attempted to further comprehend the structural details and mechanisms involved in AcCoA binding. COM pulling and Umbrella Sampling simulation are used to decipher how these critical residues allow inhibitor binding from a mechanistic point of view. Additionally, we have also created PPATR90A/D94A double-mutants and performed molecular dynamic simulation to obtain the difference in dissociation free energy landscape. To further corroborate our results, we have solved the apo structure of PPAT for first time from P. aeruginosa using X-ray crystallographic analysis. Our results, clearly shows for the first time, the mechanistic details of the inhibitor binding in the PPAT groove mediated by a gatingmechanism of R90 and D94 residues using enhanced sampling of structures throughout reaction coordinate. The observations described herein are believed to be helpful for targeted drug designing that are able to alter most crucial interactions of the binding site.

Material and method: MD simulations The PPATWT-AcCoA structure used for this study was obtained from RCSB-PDB (PDB ID: 3X1J). This was the only structure reported of PPATWT co-crystallized with its inhibitor AcCoA. The structure had a trimeric asymmetric unit. We used a single chain with lower average Bfactor value for our further analysis. The molecule was parameterized using GROMOS96 43a1 force field and also applied to the system24. Automated Topology Builder web-server was used to build the topology of the acetyl-CoA25. Particle mesh Ewald (PME) algorithm was used for calculation of long-range electrostatics and short-range nonbonded interactions were cut off at 1.4 nm26-27. Periodic boundary condition was applied on the system in all directions. Before starting the pulling simulation, the crystal structure was well equilibrated. The structure was placed in a cubic box (5 nm) with simple point charge (SPC) water28. 100 mM NaCl was added to the system with neutralizing counterions (82 Na and 82 Cl ions). The energy was minimized with steepest descents minimization method and subjected to two phase of equilibration process, 4 ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

with position restraints applied to protein-ligand complex throughout. In the first phase, the structure was simulated for 20 ns with constant volume (NVT) ensemble. Separate temperature cooling bath was used for PPAT-AcCoA and water-ions, maintained at 310 K using Berendsen weak coupling method29. In phase two, another 20 ns of constant pressure (NPT) equilibration was performed. Production MD simulations were conducted for 50 ns after withdrawing all restraints. For production MD procedure, Noose-Hoover thermostat was used for maintaining temperature30-31 and Parrinello-Rahman barostat was used for isotopically regulation of pressure32. To perform all the simulations GROMACS-4.6.2 package was used33. After the equilibration of the structure, the residues interacting with the ligand were identified (depicts the active surface area of the binding cleft) (Figure 1b) and the center of mass (COM) of these residues was used to draw the pulling vector through the COM of ligand molecule towards outside of the cleft. Followed by, AcCoA was pulled along this vector using spring constant of 500 kJ mol-1 nm-2 and pull rate of 0.01 nm ps-1 to achieve complete separation. To check the method of pulling, slower pulling rate was also used which yielded nearly same force-vs-time curves. During the pulling simulation, snapshots were taken to generate staring configurations for umbrella sampling windows34-35. The unbinding vector was defined dynamically for each individual umbrella sampling simulations based on the COM of the ligand and COM of the binding cleft’s opening36-37. We analyzed the pulling trajectories and collected different windows for a robust potential of mean force (PMF) reconstruction with applied harmonic potential for the simulations. We used 35 different trajectories for reconstruction until the PMF appears convergent and total simulation time used was >350 ns. Error estimation was performed using bootstrap method38. Purification, crystallization and X-ray diffraction of apo-PPAT: The phosphopantetheine adenylyltransferase (PPATWT) was cloned and purified according to the previously reported methodologies23. Sparse matrix screening was performed with Crystal screens (Hampton research) using 16.6 mg/ml protein. Oily patches were appeared, after two months, when 2 µl of protein solution was added with 2 µl of Na-formate and incubated against Na-formate. To get diffraction size crystals a gradient of 3.25-3.3 M Na-formate was used and was optimized with or without 10 mM Caacetate. The crystals developed cracks when grown beyond 0.6 mm and had layered morphology. Additive screening of PPAT with 4% Polyvinyl Pyroliddone K-15 showed uneven shaped layered crystals and 1-5% PEG400 yielded showered crystals. Optimal sized crystals obtained in 3.25-3.3 M Na-formate with 10 mM Ca-acetate, were harvested in 9 well trays and were flash frozen in -180º C in cryo stream. All the diffraction data, for optimization purposes, were collected in our home source using Rigaku R-axis IV++ detector with 0.5º oscillation mode. A seed-stock was prepared using 2-3 single crystals and subjected to microseeding matrix screening. Crystals appeared using PEG 4000, 20% isopropanol and HEPES pH 7.0 condition. Gradual evaporation of alcohol made harvesting initial crystals very difficult. Finally, 16-18% 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

PEG 4000, 0.1 M HEPES pH 7.1, 5% isopropanol and 200 mM Na-acetate yielded diffraction quality crystals. 25-30% DMSO was used as cryo-protectant with the reservoir solution. Final data collection was performed in INDUS 2 Beamline PX21 using marCCD X-ray detector with 1º oscillation range at 93 K. Structure Solution: The oscillation files were integrated using HKL-2000 and iMosfilm39, followed by scaling was performed in Scalepack40. Crystals belong to the space group C121. PPAT-AcCoA structure from P. aeruginosa (3X1J) was used as a starting model for structure solution using molecular replacement method in MOLREP41 program. Subsequent model building was performed in Coot0.8.3 and refinement procedure was carried out using PHENIX.refine package42. MolProbity package was used for further stereochemical validation of the structure of PPAT43-44. Schrodinger and Chimera package was used for identifying interaction sites and molecular visualizations45.

RESULTS AND DISCUSSION: Interaction of inhibitor (AcCoA) with binding cavity PPAT is a very crucial enzyme responsible for catalyzing the penultimate rate limiting step in CoA biosynthesis pathway. The relation of this highly conserved metabolic pathway with regulation of bacterial virulence system (e.g. Type III secretion system) makes it an intriguing

Figure 1. Snapshot of Acetyl-CoA(represented in CPK) binding cavity inside PPAT of P. aeruginosa. a. Trimeric PPAT (representing one asymmetric unit) with bound Acetyl-CoA (PDB ID: 3X1J) with each protomer. b. 2D interaction diagram representing the major interactions between AcCoA and PPAT in the binding cleft. These interacting residues were used to calculate the COM of PPAT’s binding site.

subject of study. After all these years of study it came up very recently that Acetyl-CoA also can bind with the PPAT and downregulate the CoA biosynthesis. For our studies, we have obtained 6 ACS Paragon Plus Environment

Page 6 of 28

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the only available crystal structure of PPAT bound with inhibitor Acetyl-CoA (Figure 1a) from RCSB-PDB. Later, we have chosen a single chain with lower average B-factor value for our further mechanistic studies. We performed 50 ns of NVT and NPT simulations on this structure using GROMACS-5.1.4 package and explored the interacting residues around the binding cavity of PPAT. Using Ligplot package46, we have identified all the interacting residues within 3.5 Å cut-off range of the bound inhibitor molecule (Figure 1b). In agreement with the prior report23 the equilibrated structure also shows residues like R90, R87, D94, K41, L73, R87, T129, H17 (neutral state) and T9 involved in electrostatic interactions within that distance cut off. As expected, the non-restrained molecular dynamics results yielded nearly same interaction diagram of allosteric site in accordance with the reported crystal structure (PDB ID: 3X1J). Thermodynamics of Acetyl-CoA binding: The prior study concerned with inhibitor binding also raised an intriguing hypothesis that possibly D94 and R90 residues are acting in a choreographic manner to perform the gating mechanism near binding cavity23. However, further clarification needed the crystal structure of apoPPAT from P. aeruginosa. In this study, we have used umbrella sampling (US) method to understand the energetics by which Acetyl-CoA gets Figure 2. PMF for the dissociation of Acetyl-CoA from the binding site into the binding cavity. of PPAT. The error associated with energy minima is ±0.14 kcal/mol. US gives a reliable way to The reaction coordinate here is represented in ζ (nm). comprehend systems behavior for certain time scales (say, transition phase) which might have otherwise remained inaccessible using both experimental and conventional molecular dynamics (MD) approach. US procedure had been used in earlier studies with reliable results from protein-protein, protein-ligand, lipid bilayerprotein systems46. The pulling force applied on the system perturbs the equilibrium of the system thus making it possible to calculate thermodynamic quantities from US trajectories without large errors. Further, we used weighted histogram analysis method (WHAM) for calculation of free 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

energy47 as it uses simulations performed on number of configurational windows generated from umbrella sampling simulation for calculation of ∆G.

We used COM pulling and US simulation to estimate the binding energetics as well as to understand the mechanism of inhibitor binding in the cavity of PPAT. After equilibrating the PPAT-AcCoA crystal structure (Figure S6), the residues involved in interaction with the ligand (within 3.5 Å cutoff range) were identified (Figure 1b). The unbinding vector was defined dynamically for each individual umbrella sampling simulation based on the COM of these interacting residues (representing the active surface area of the ligand binding cleft) through the COM of the ligand towards the opening of the cleft36-37 (Figure 3a). Starting configurations and velocities for the pulling simulations were taken at 15 ps intervals. The ∆G obtained from this US simulation (Figure 2) is in close agreement with the previously reported experimental observation23, which justifies the authenticity of our sampling procedure. Analysis of the snapshots collected from the trajectory showed that R90 and D94 residues together moved quite a lot to maintain the interaction with the ligand while in binding cavity (Figure 2-inset pictures). After the ligand got completely dissociated the displacement of R90 was found to be 12.5 Å and D94 residue was also found to be 4.8 Å away from the starting position (Video S11). Figure 4 clearly shows the difference in orientation of R90 and D94 as the ligand is moving out of its binding cleft. These observations are in good agreement with the hypothesis derived previously regarding the Figure 3. a. Overview of the vector used for movement of these residues. Therefore, it can now be center-of mass pulling along the egression of the binding cleft. b. View of PPAT binding said that these two residues do play very crucial role cleft after 0.5 nm displacement of AcCoA in insertion and accommodation of ligand inside the along the cleft opening. c. Changes in binding cavity. A large helix movement was also neighboring contacts (represented in black seen after AcCoA pulled out and possibly was to dotted lines) inside binding cleft as the ligand prevent entry of water molecules inside allosteric further egress 1 nm. cleft (discussed later) (Figure 4d).

8 ACS Paragon Plus Environment

Page 8 of 28

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4. Snapshot of Acetyl-CoA at the entry gate of PPAT binding cavity. a. Shows the orientation of critical residues when Acetyl-CoA approaches at the gate. b. 2D interaction diagram of ligand acetyl-CoA with the PPAT while the ligand is residing near the gate of binding site. c. Superposed structures of apoPPAT(off white) with AcCoA-PPAT(sky blue) showing outward movement of R90 and D94 while binding cavity is empty. d. Superimposed structure of simulated-PPAT (after complete dissociation AcCoA, in off-white color) and apoPPAT crystal structure (sky blue) showing nearly same kind of outward movements of R90 and D94.

Thermodynamics of PPATR90A/D94A-AcCoA binding: Using PMF calculation and sampling over a large set of configurational windows, we have found evidences that led our belief towards the understanding of R90-D94 mediated gating mechanism. However, to completely rely on the results produced so far, we needed to check the free-energy profile of the same system, mutating those two critical residues. For this purpose, we took the structure prepared earlier after 100 ns of equilibration. The in-silico mutations intended were performed using Pymol package by replacing R90 and D94 residues to alanines (∆PPAT or PPATR90A/D94A). We used alanine mutations to eliminate any possible formation of ionic interaction of these two positions with ligand molecule. Thus, it would provide us a way to analyze either R90-D94 residues were acting as a gate or we just hit an artefactual mechanism.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The mutant was further equilibrated for 100 ns using both NVT and NPT method until the backbone RMSD reached certain stability. After performing thorough equilibration, the structure was subjected to US procedure. The COM of interacting residues around the binding site of ligand (within 3.5 Å cut-off distance) was used to calculate the force vector through the COM of Figure 5. PMF for AcCoA dissociation from PPATR90A/D94A binding ligand (same as used cleft. R90 and D94 residues are replaced with alanine double mutation. previously). COM Well equilibrated structure is used for PMF reconstruction. The error associated with energy minima is ±0.13 kcal/mol. The reaction pulling was performed with 500 kJ mol-1 nm-2 coordinate here is represented in ζ (nm). and pull rate of 0.01 nm -1 ps to achieve complete separation of ligand from PPATR90A/D94A. During this puling 36 different trajectories were recorded (Figure S9) and further umbrella sampling (US) simulations were performed utilizing over 360 ns of computation time for obtaining a convergent PMF∆PPAT. As evident very clearly from the figure that the PMF∆PPAT shows a considerable difference compared to PMFPPAT. PMF calculation shows that the ∆G of binding happens to be near -4.5 kcal/mol for PPATR90A/D94A-AcCoA. COM pulling started the dissociation around the ξ = -2.6 nm and it continued till the ligand got completely separated around ξ = 0.0 nm. The difference observed in PMF∆PPAT clearly suggests that due to the mutation of the R90 and D94 the ligand failed to maintain its position inside the binding cavity. Based on the results we must admit that the differences in ∆G between PPATWT and mutant are small (-2.5 kcal/mol). However, we would like to consider this significant as the ∆G obtained in case of PPATWT goes in accordance with previously reported experimental evidences23. We gave the same treatment for PPATR90A/D94A structure and thus we have a reason to believe that the ∆∆G, however small, but carries significant relevance towards the focus of this work. The dissociation of ligand from allosteric site in COM pulling simulation was also examined using different pulling rates (0.10 and 0.05 nm ps-1) to ensure the pulling rate (0.01 nm ps-1) used in this article did not induce any artifacts. Figure S7 shows that the ligand dissociation occurred in same way in different pulling rates, yielding same sequence of configurations. The required time to attain the maximum force and the magnitude of this force depends on the rate of pulling as the structural changes are happening at different rates. 10 ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The PPATWT is found in hexameric state in nature. In our previous study23, we have showed that the AcCoA binds to all the protomers and the binding occurs mostly in independent or noncooperative manner. Structural analysis also showed that there is not much difference among individual AcCoA-bound protomers. However, compared to structures reported in other organisms, we found a unique region that possibly assists inhibitor binding/unbinding. We performed the simulation studies using a monomer of P. aeruginosa-PPAT-AcCoA as our primary concern was to decipher the unique ‘gating-mechanism’ of R90-D94 residues upon ligand movement. X-ray crystallographic analysis of apo-PPAT from P. aeruginosa: Based on the extensive simulations performed on both PPAT-AcCoA and PPATR90A/D94A -

Figure 6. Cartoon superposition with Cα traced superposition of the PPAT enzyme from various organisms. a-b. Cartoon representation and its corresponding Cα traced superposition. Showing clearly the loop movement (α3-β2) of M. tuberculosis-PPAT. c-d. 90˚ rotational view of a-b representation. Color code: Black – newly obtained structure of apo-PPAT(P. aeruginosa), PDB ID: 5X6F, red – 5TS2, green – 4RUK, blue – 3X1M, brown – 3X1K, pink – 3X1J, dark blue – 1TFU.

AcCoA structures, quite a vivid picture is formed regarding the mechanistic point of view of ligand accommodation. However, prior studies using X-ray crystallographic structures of PPAT from M. tuberculosis showed that the ligand accommodation is facilitated by movement of a loop region near the binding site. In a recent report Chatterjee et al.23 hypothesized that the mechanism employed by P. aeruginosa is completely different with that of reported earlier in M.

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

tuberculosis. Our results so far corroborate with the later statement clarifying the fact that instead of loop movement P. aeruginosa uses R90 and D94 residue’s gating mechanism (Figure 6). However, to establish such claim we needed to obtain the first crystal structure of full-length poform of PPAT from P. aeruginosa. Obtaining the crystal structure was very difficult as the purified protein always found to be in complex with either AcCoA or CoA. This experimental observation of strong binding also corroborates with the slow dissociation pattern of AcCoA from binding cavity upon application of constant force. The crystals of apo-PPAT always appeared in a layered manner, leading to interfere with data collection strategy (Figure S1, S2, S3). However, eventually we solved the crystal structure in C121 space group (Figure 7, Table S10).

Figure 7. Crystal structure of apo-PPAT from P. aeruginosa. The structure shows a hexameric asymmetric unit. The protomers of the asymmetric unit is shown in cartoon representation.

Analysis of the active site using PPAT-AcCoA and PPAT crystal structure revealed a striking similarity in movement of R90-D94 residues between the obtained results from US simulation and the crystal structure (Figure 4d). Figure 4c shows that in absence of bound AcCoA both R90 and D94 residues are shown to deviate 4.3 Å and 2.8 Å with side chain flip, respectively. E96 residue also found to be displaced by 3.2 Å in the PPAT crystal structure for avoiding any clash due to the movement of D94. However, a little difference in displacements is encountered between COM-pulling results and X-ray crystallographic structure. Solvent content analysis of the crystal structure revealed that in absence of any ligand, a glycerol and an isopropanol molecule are found to sit inside the binding cavity. This observation also indicates that R90 and D94 residue shows a lot of activity near binding cavity to provide ionic interactions and thereby maintaining stability of the binding cavity.

12 ACS Paragon Plus Environment

Page 12 of 28

Page 13 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 8. Crystal structure superposition of apo forms of M. tuberculosis (sky blue) and P. aeruginosa (off white). a. Superposition cartoon diagram of apo structure from both the organism. The loop movement is very clear between α3 and β2. b-d. Cartoon superposition showing the deviations of critical residue near binding cavity of PPAT in P. aeruginosa, b. apoPPAT(off white) – anpnp-PPAT, c. apoPPAT – coA-PPAT(sea green), d. apoPPAT – dPCoA-PPAT(sea green).

To further illustrate the role of R90 and D94 gate, we used all the previously reported structures of PPAT (ligand bound) from P. aeruginosa and aligned the with apo-PPAT crystal structure. Figure 8b shows that D94 and R90 residue moved 6.9 Å and 3.7 Å, respectively to facilitate binding of AMP-PNP inside binding cavity. For providing binding interactions to dPCoA, D94 residue moved 2 Å towards binding cavity and R90 seems to move 2.1 Å closer to the ligand for providing stability (Figure 8d). Nearly similar observations were found in case of CoA binding also. The D94 residue is found to move 4.5 Å towards the binding pocket while R90 residue moved 2.1 Å towards the CoA to stabilize it (Figure 8c). Structural comparison with apo form of PPAT of M. tuberculosis made it evident that such gating mechanism is a characteristic phenomenon of P. aeruginosa for ligand accommodation in binding cavity. Figure 8a shows the loop movement in P. aeruginosa, with D94 and R90 residues, compared to M. tuberculosis in absence of any ligand. Indeed, Chatterjee et al23 argued about the existence of different kind of binding mechanism in P. aeruginosa which found to be regulating cellular level of CoA and AcCoA (Figure 9); affecting various pathogenic pathways like type III secretion system. 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 28

However, our work here have provided the first demonstration regarding how this highly choreographed phenomena facilitate ligand accommodation and provide stability inside the binding cavity.

Figure 9. A schematic representation showing metabolic regulation of T3SS by cellular levels of AcetylCoA and CoA.

Conclusion: We have performed extensive molecular dynamic simulations using various approaches to provide detailed understanding of ligand binding inside PPAT binding cavity. We used explicit solvent MD simulations, COM pulling and umbrella sampling strategy to comprehend the mechanism of binding mediated by the key residues. In agreement with previous studies23, we find that the mode of ligand accommodation inside binding cavity is very different in P. aeruginosa compared to earlier report in M. tuberculosis. We have further deciphered the R90 and D94 residue mediated gating mechanism near the binding cavity that help maintain ligand accommodation and its stability by moving in and out from binding cleft. The PMF obtained from COM pulling MD simulations unraveled the thermodynamic details indicating the rapid movement of R90 and D94 residues to accommodate ligand molecule in binding cavity and providing stability to it. Expanding upon these findings, we obtained the first crystal structure of apo-PPAT from P. aeruginosa to further validate our prior theoretical observations. Corroborating with our extensive molecular dynamic studies, the X-ray crystallographic study also revealed that D94 and R90 residues indeed show a large movement near binding site and hence provide a gating mechanism for ligand accommodation and stability. Further validating these finding, we have shown considerable difference in binding free energy using PMF 14 ACS Paragon Plus Environment

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

calculation of PPATR90A/D94A mutant. Thus, we propose that small molecular compounds, interfering with these particular residues or blocking the gating mechanism may provide a viable way for controlling the cellular level of CoA or Acetyl-CoA. This in turn could be helpful to down regulate the expression level of various pathogenic pathways like type III secretion system’s exotoxins and thereby avoiding the infection.

Supporting Information Supporting information contains figures, table and video files relevant to this manuscript.

Acknowledgements: We thank R.P. Thangavelu (Sr. Principal Scientist, CSIR CMMACS, Bangalore) for providing access to the Supercomputing facility. Financial support was provided by the Department of Science and Technology, Govt. of India (Grant no: SB/SO/BB – 36/2014) and Council of Scientific and Industrial Research, Govt. of India (Grant no: TREAT – BSC0113, UNSEEN – BSC0116). CSIR-SRF fellowship was provided to AM.

References: 1. Roy-Burman, A.; Savel, R. H.; Racine, S.; Swanson, B. L.; Revadigar, N. S.; Fujimoto, J.; Sawa, T.; Frank, D. W.; Wiener-Kronish, J. P. Type III protein secretion is associated with death in lower respiratory and systemic Pseudomonas aeruginosa infections. J. Infect. Dis. 2001, 183, 1767-1774. 2. Kurahashi, K.; Kajikawa, O.; Sawa, T.; Ohara, M.; Gropper, M. A.; Frank, D. W.; Martin, T. R.; Wiener-Kronish, J. P. Pathogenesis of septic shock in Pseudomonas aeruginosa pneumonia. J. Clin. Invest. 1999, 104, 743. 3. Lynch, J. P. Hospital-acquired pneumonia: risk factors, microbiology, and treatment. CHEST 2001, 119, 373S-384S. 4. Otter, J. A.; Yezli, S.; French, G. L. The role played by contaminated surfaces in the transmission of nosocomial pathogens. Infect. Control Hosp. Epidemiol. 2011, 32, 687-699. 5. Sulakvelidze, A.; Alavidze, Z.; Morris, J. G. Bacteriophage therapy. Antimicrob. Agents Chemother. 2001, 45, 649-659. 6. Begley, T. P.; Kinsland, C.; Strauss, E. The biosynthesis of coenzyme A in bacteria. Vitam. Horm. 2001, 61, 157-171. 7. Osterman, A.; Begley, T. A subsystems-based approach to the identification of drug targets in bacterial pathogens. In: Boshoff, H. I.; Barry, C. E. Eds; Systems Biological Approaches in Infectious Diseases. Progress in Drug Research 2007; vol 64. 131-170. 8. Spry, C.; Kirk, K.; Saliba, K. J. Coenzyme A biosynthesis: an antimicrobial drug target. FEMS Microbiol. Rev. 2007, 32, 56-106. 9. Jackowski, S.; Rock, C. Regulation of coenzyme A biosynthesis. J. Bacteriol. 1981, 148, 926-932. 10. Worrall, D.; Tubbs, P. A bifunctional enzyme complex in coenzyme A biosynthesis: purification of pantetheine phosphate adenylyltransferase and dephospho-CoA kinase. Biochem. J. 1983, 215, 153157. 11. Viala, J. P.; Prima, V.; Puppo, R.; Agrebi, R.; Canestrari, M. J.; Lignon, S.; Chauvin, N.; Méresse, S.; Mignot, T.; Lebrun, R. Acylation of the Type 3 Secretion System Translocon Using a Dedicated Acyl Carrier Protein. PLoS Genet. 2017, 13, e1006556. 15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

12. Wilharm, G.; Heider, C. Interrelationship between type three secretion system and metabolism in pathogenic bacteria. Front. Cell. Infect. Microbiol. 2014, 4. 13. Rao, S. T.; Rossmann, M. G. Comparison of super-secondary structures in proteins. J.Mol. Biol. 1973, 76, 241-56. 14. Izard, T. The crystal structures of phosphopantetheine adenylyltransferase with bound substrates reveal the enzyme’s catalytic mechanism. J. Mol. Biol. 2002, 315, 487-495. 15. Izard, T.; Geerlof, A. The crystal structure of a novel bacterial adenylyltransferase reveals half of sites reactivity. EMBO J. 1999, 18, 2021-2030. 16. Edwards, T. E.; Leibly, D. J.; Bhandari, J.; Statnekov, J. B.; Phan, I.; Dieterich, S. H.; Abendroth, J.; Staker, B. L.; Van Voorhis, W. C.; Myler, P. J. Structures of phosphopantetheine adenylyltransferase from Burkholderia pseudomallei. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 2011, 67, 1032-1037. 17. Eom, S. J.; Ahn, H. J.; Kim, H. W.; Baek, S.-H.; Suh, S. W. Crystallization and preliminary X-ray crystallographic studies of phosphopantetheine adenylyltransferase from Helicobacter pylori. Acta Crystallogr. Sect. D Biol. Crystallogr. 2003, 59, 561-562. 18. Izard, T.; Geerlof, A.; Lewendon, A.; Barker, J. J. Cubic crystals of phosphopantetheine adenylyltransferase from Escherichia coli. Acta Crystallogr. Sect. D Biol. Crystallogr. 1999, 55, 1226-1228. 19. Kang, J. Y.; Lee, H. H.; Yoon, H. J.; Kim, H. S.; Suh, S. W. Overexpression, crystallization and preliminary X-ray crystallographic analysis of phosphopantetheine adenylyltransferase from Enterococcus faecalis. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 2006, 62, 1131-1133. 20. Lee, H. H.; Yoon, H.-J.; Kang, J. Y.; Park, J. H.; Kim, D. J.; Choi, K.-H.; Lee, S.-K.; Song, J.; Kim, H.-J.; Suh, S. W. The structure of Staphylococcus aureus phosphopantetheine adenylyltransferase in complex with 3′-phosphoadenosine 5′-phosphosulfate reveals a new ligand-binding mode. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 2009, 65, 987-991. 21. Takahashi, H.; Inagaki, E.; Fujimoto, Y.; Kuroishi, C.; Nodake, Y.; Nakamura, Y.; Arisaka, F.; Yutani, K.; Kuramitsu, S.; Yokoyama, S. Structure and implications for the thermal stability of phosphopantetheine adenylyltransferase from Thermus thermophilus. Acta Crystallogr. Sect. D Biol. Crystallogr. 2004, 60, 97-104. 22. Timofeev, V.; Smirnova, E.; Chupova, L.; Esipov, R.; Kuranova, I. X-ray study of the conformational changes in the molecule of phosphopantetheine adenylyltransferase from Mycobacterium tuberculosis during the catalyzed reaction. Acta Crystallogr. Sect. D Biol. Crystallogr. 2012, 68, 1660-1670. 23. Chatterjee, R.; Mondal, A.; Basu, A.; Datta, S. Transition of phosphopantetheine adenylyltransferase from catalytic to allosteric state is characterized by ternary complex formation in Pseudomonas aeruginosa. Biochim. Biophys. Acta 2016, 1864, 773-786. 24. Pol-Fachin, L.; Fernandes, C. L.; Verli, H. GROMOS96 43a1 performance on the characterization of glycoprotein conformational ensembles through molecular dynamics simulations. Carbohydr. Res. 2009, 344, 491-500. 25. Malde, A. K.; Zuo, L.; Breeze, M.; Stroet, M.; Poger, D.; Nair, P. C.; Oostenbrink, C.; Mark, A. E. An automated force field topology builder (ATB) and repository: version 1.0. J. Chem. Theory Comput. 2011, 7, 4026-4037. 26. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N⋅ log (N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089-10092. 27. Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577-8593. 28. Berendsen, H. J.; Postma, J. P.; van Gunsteren, W. F.; Hermans, J. Interaction models for water in relation to protein hydration. In: Pullman, B. Eds.; Intermolecular forces, The Jerusalem Symposia on Quantum Chemistry and Biochemistry, Springer: 1981; pp 331-342.

16 ACS Paragon Plus Environment

Page 16 of 28

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

29. Berendsen, H. J.; Postma, J. v.; van Gunsteren, W. F.; DiNola, A.; Haak, J. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684-3690. 30. Hoover, W. G. Canonical dynamics: equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695. 31. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511-519. 32. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: a new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182-7190. 33. Hess, B.; Kutzner, C.; Van Der Spoel, D.; Lindahl, E. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 2008, 4, 435-447. 34. Patey, G.; Valleau, J. The free energy of spheres with dipoles: monte carlo with multistage sampling. Chem. Phys. Lett. 1973, 21, 297-300. 35. Torrie, G. M.; Valleau, J. P., Monte Carlo free energy estimates using non-Boltzmann sampling: application to the sub-critical Lennard-Jones fluid. Chem. Phys. Lett. 1974, 28, 578-581. 36. Allen, W. J.; Bevan, D. R. Steered molecular dynamics simulations reveal important mechanisms in reversible monoamine oxidase B inhibition. Biochem. 2011, 50, 6441-6454. 37. Mascarenhas, N. M.; Kästner, J. How maltose influences structural changes to bind to maltosebinding protein: Results from umbrella sampling simulation. Proteins 2013, 81, 185-198. 38. Hub, J. S.; De Groot, B. L.; Van Der Spoel, D. g_wham - A free weighted histogram analysis implementation including robust error and autocorrelation estimates. J. Chem.Theory Comput. 2010, 6, 3713-3720. 39. Battye, T. G.; Kontogiannias, L.; Johnson, O.; Powell, H. R.; Leslie, A. G. iMOSFLM: a new graphical interface for diffraction-image processing with MOSFLM. Acta Crystallogr. D Biol. Crystallogr. 2011, 67, 271-281. 40. Otwinowski, Z.; Minor, W. [20] Processing of X-ray diffraction data collected in oscillation mode. Methods Enzymol. 1997, 276, 307-326. 41. Vagin, A.; Teplyakov, A. MOLREP: an automated program for molecular replacement. J. Appl. Crystallogr.1997, 30, 1022-1025. 42. Afonine, P. V.; Grosse-Kunstleve, R. W.; Echols, N.; Headd, J. J.; Moriarty, N. W.; Mustyakimov, M.; Terwilliger, T. C.; Urzhumtsev, A.; Zwart, P. H.; Adams, P. D. Towards automated crystallographic structure refinement with phenix. refine. Acta Crystallogr. Sect. D Biol. Crystallogr. 2012, 68, 352-367. 43. Chen, V. B.; Arendall, W. B.; Headd, J. J.; Keedy, D. A.; Immormino, R. M.; Kapral, G. J.; Murray, L. W.; Richardson, J. S.; Richardson, D. C. MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr. Sect. D Biol. Crystallogr. 2010, 66, 12-21. 44. Davis, I. W.; Leaver-Fay, A.; Chen, V. B.; Block, J. N.; Kapral, G. J.; Wang, X.; Murray, L. W.; Arendall III, W. B.; Snoeyink, J.; Richardson, J. S. MolProbity: all-atom contacts and structure validation for proteins and nucleic acids. Nucleic Acids Res. 2007, 35, W375-W383. 45. Pettersen, E. F.; Goddard, T. D.; Huang, C. C.; Couch, G. S.; Greenblatt, D. M.; Meng, E. C.; Ferrin, T. E. UCSF Chimera—a visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605-1612. 46. Izrailev, S.; Stepaniants, S.; Isralewitz, B.; Kosztin, D.; Lu, H.; Molnar, F.; Wriggers, W.; Schulten, K. Steered molecular dynamics. In: Deuflhard, P.; Hermans, J.; Leimkuhler, B.; Mark, A.E.; Reich, S.; Skeel, R.D. Eds.; Computational Molecular Dynamics: Challenges, Methods, Ideas. Springer, Berlin, Heidelberg 1999, 4, 39-65. 47. Kumar, S.; Rosenberg, J. M.; Bouzida, D.; Swendsen, R. H.; Kollman, P. A. THE weighted histogram analysis method for free-energy calculations on biomolecules. I. The method. J. Comput. Chem. 1992, 13, 1011-1021.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic Title: Umbrella Sampling and X-ray Crystallographic Analysis Unveil an Arg-Asp Gate Facilitating Inhibitor Binding Inside Phosphopantetheine Adenylyltransferase Allosteric Cleft A. Mondal, R. Chatterjee and S. Datta*

18 ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. Snapshot of Acetyl-CoA(represented in CPK) binding cavity inside PPAT of P. aeruginosa. a. Trimeric PPAT (representing one asymmetric unit) with bound Acetyl-CoA (PDB ID: 3X1J) with each protomer. b. 2D interaction diagram representing the major interactions between AcCoA and PPAT in the binding cleft. These interacting residues were used to calculate the COM of PPAT’s binding site.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. PMF for the dissociation of Acetyl-CoA from the binding site of PPAT. The error associated with energy minima is ±0.14 kcal/mol. The reaction coordinate here is represented in ζ (nm). 38x29mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. a. Overview of the vector used for center-of mass pulling along the egression of the binding cleft. b. View of PPAT binding cleft after 0.5 nm displacement of AcCoA along the cleft opening. c. Changes in neighboring contacts (represented in black dotted lines) inside binding cleft as the ligand further egress 1 nm. 356x905mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Snapshot of Acetyl-CoA at the entry gate of PPAT binding cavity. a. Shows the orientation of critical residues when Acetyl-CoA approaches at the gate. b. 2D interaction diagram of ligand acetyl-CoA with the PPAT while the ligand is residing near the gate of binding site. c. Superposed structures of apoPPAT(off white) with AcCoA-PPAT(sky blue) showing outward movement of R90 and D94 while binding cavity is empty. d. Superimposed structure of simulated-PPAT (after complete dissociation AcCoA, in offwhite color) and apoPPAT crystal structure (sky blue) showing nearly same kind of outward movements of R90 and D94.

ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. PMF for AcCoA dissociation from PPATR90A/D94A binding cleft. R90 and D94 residues are replaced with alanine double mutation. Well equilibrated structure is used for PMF reconstruction. The error associated with energy minima is ±0.13 kcal/mol. The reaction coordinate here is represented in ζ (nm). 42x35mm (600 x 600 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Cartoon superposition with Cα traced superposition of the PPAT enzyme from various organisms. a-b. Cartoon representation and its corresponding Cα traced superposition. Showing clearly the loop movement (α3-β2) of M. tuberculosis-PPAT. c-d. 90˚ rotational view of a-b representation. Color code: Black – newly obtained structure of apo-PPAT(P. aeruginosa), PDB ID: 5X6F, red – 5TS2, green – 4RUK, blue – 3X1M, brown – 3X1K, pink – 3X1J, dark blue – 1TFU. 34x23mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 7. Crystal structure of apo-PPAT from P. aeruginosa. The structure shows a hexameric asymmetric unit. The protomers of the asymmetric unit is shown in cartoon representation.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 8. Crystal structure superposition of apo forms of M. tuberculosis (sky blue) and P. aeruginosa (off white). a. Superposition cartoon diagram of apo structure from both the organism. The loop movement is very clear between α3 and β2. b-d. Cartoon superposition showing the deviations of critical residue near binding cavity of PPAT in P. aeruginosa, b. apoPPAT(off white) – anpnp-PPAT, c. apoPPAT – coA-PPAT(sea green), d. apoPPAT – dPCoA-PPAT(sea green). 37x27mm (600 x 600 DPI)

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 9. A schematic representation showing metabolic regulation of T3SS by cellular levels of Acetyl-CoA and CoA. 118x81mm (600 x 600 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic 42x21mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 28 of 28