Ing minimized output conformations, a series of MD simulations had been performed for 20 picoseconds (ps) at ten K, with an integration step of 0.five femtoseconds (fs), to further reinforce the search for minimum-energy states. The relative stabilities of each conformation, obtained from the ten K MD final frame, had been used to construct relaxed power contour plots (Fig. S7) describing the conformation of every single glycosidic linkage.Components and Methods Application and NomenclatureAll saccharide and PAPS topologies have been generated using the PRODRG server [29]. Posterior manipulation in the structures was performed applying MOLDEN [30], to draw ligand structures and manage ab initio-derived files, VMD [31], for visualization of the trajectories, and PyMOL, to create the NST mutants [32]. Molecular dynamics (MD), including each calculations and analyses, was performed making use of GROMACS simulation suite, version four.5.1 and GROMOS96 43a1 force field [33]. The relative orientation of a pair of contiguous carbohydrate residues is described by two or three torsional angles at the glycosidic linkage. For an (1R4) linkage, the W and Y have been defined as shown in the equations 1 and 2: w O5{C1{OX {CX y C1{OX {CX {C(X {1)Atomic D1 Receptor drug Charge CalculationThe calculation of atomic charges suitable for the PAPS MD simulations was performed as previously described for other sulfated compounds [34,35]. Briefly, the PAPS was submitted to full geometry optimization at the HF/3-21G level using GAMESS [36]. Subsequently, the minimal energy conformations were submitted to single-point calculations at the HF/6-31G levelPLOS ONE | plosone.orgDocking ProceduresAutoDock4.2 was used as grid-based docking procedure [41]. Although heparan sulfate docking validation against crystal structures such as IL9 and have been performed elsewhere [42],Molecular Dynamics of N-Sulfotransferase Activitywe performed a docking using the 1.8 A resolution structure of 3O-sulfotransferase bound to a heptasaccharide substrate using Autdock (PDBiD 1T8T). The obtained reference RMSD was 0.49 for the lowest scoring energy population (212 Kcal/mol able S1). The crystal structure of the sulfotransferase domain of human heparan sulfate N-deacetylase/N-sulfotransferase 1 bound to PAP obtained from Brookhaven Protein Data Bank (PDB ID code: 1NST) [25] was used in the docking experiments. Missing side chain atoms and the Kollman united atom partial charges to the PAPS molecule were included [43,44]. Concerning carbohydrate structures, the Lowdin atomic charges, as previously calculated for sulfated saccharides [39,40], were employed and all torsion angles were considered flexible. The grid maps, calculated using AutoGrid, were chosen to be large enough to include the active site, as well as a significant portion of the surrounding surface. The dimensions of the grids were thus 50 A 650 A 640 A, with 0.3 A spacing between the grid points. Docking of the disaccharide to 1NST was Kinesin review carried out using the empirical free energy function and the Lamarckian genetic algorithm, applying a standard protocol with an initial population of 500 randomly placed individuals, a maximum number of 2.56108 energy evaluations, a 0.02 mutation rate, a 0.80 crossover rate, and an elitism value of 1, where the average of the worst energy was calculated over a window of the previous 10 generations. One hundred independent docking runs were carried out for the disaccharide. Results were clustered according to the 0.5 A root-mean-square deviation (RMSD).
