Equilibrated atomic models of outward-facing P-glycoprotein and effect of ATP binding on structural dynamics (Supplementary Information)
|
|
- Byron Owen
- 5 years ago
- Views:
Transcription
1 Equilibrated atomic models of outward-facing P-glycoprotein and effect of ATP binding on structural dynamics (Supplementary Information) Lurong Pan 1 and Stephen G. Aller 2 * 1,2 Department of Pharmacology and Toxicology, University of Alabama at Birmingham, th Street South, Birmingham, AL, USA *Correspondence to sgaller@uab.edu Text S1. Equilibration and sampling size The protein C α root-mean-square deviation (RMSD) revealed that the major conformational change occurred in the first 100 ns of simulation and the protein reached different equilibrium states after 100 ns (Figure S1). The equilibrated RMSDs of minor changes within a 1 Å range indicated that the protein adopted certain predominant meta-stable conformations and oscillated between them. The equilibration findings showed that the time scale of the simulations was sufficiently long to reveal various meta-stable conformations that were independent from the starting structure. Both apo-pgp and MgATP-Pgp were simulated in triplicate to acquire sufficient sampling and to avoid biased conclusions. According to the Ergodic hypothesis, the time average is equivalent to the ensemble average. Therefore, the average of the trajectories from triplicate simulations over time was used in order to calculate the dynamic features of the each model in an ensemble with considerable sampling size. Representative microstates with distinct conformational features were also analyzed. Text S2. Molecular dynamic simulations Atomic coordinates of outward-facing Pgp with/without MgATP 2- (apo-pgp/mgatp-pgp) was first placed in a POPC bilayer. The dimension of the inward-facing crystal structures was (Å) and the dimension of the outward-facing homology structure was (Å). The change of dimension involves the fluctuation of TMDs from inward-facing to outward-facing conformations and potential rotation of the protein around the bilayer center normal (z-axis). Therefore, the lipid bilayer was built as (Å) in order to accommodate all possible dynamic conformations. POPC was used as the sole bilayer-forming lipid in our simulations. The long principal axis of the protein was aligned with the bilayer center normal (z-axis). Previous studies indicated that Trp residues are positioned to form hydrogen bonds between the -NH group on the indole ring and carbonyl group of the glyceride on the lipid while the majority of the indole rings are immersed towards the hydrophobic side of the bilayer 1. In order to achieve a proper lipidprotein interface for both leaflets of the lipid bilayer, residues Trp208, Trp311 and Trp851 were aligned with the glyceride groups of the outer leaflet of membrane whereas residues Trp44, Trp132, Trp228, Trp694 and Trp704 were aligned with the glyceride groups of the inner leaflet of the membrane. After the protein was positioned in the lipid bilayer correctly, the overlapping lipid groups were removed with 0.6 Å cutoff distance from the protein. The membrane-embedded protein was then dissolved in a water box without removing the water molecules inside the bilayer. Eliminating water within the lipid would introduce micro vacuum pockets in the initial structure, and we therefore fully solvated the lipid to allow unbiased adjustment of the system to a proper lipid solvation state. After testing different system sizes, we were able to determine that larger padding
2 distances do not provide any more stability or any apparent advantage over standard padding, allowing us to proceed with a smaller overall simulation box for our planned experiments. In this smaller box, the dimensions adjusted to ~ (Å) after equilibration. The padding distance maintains >12 Å from every direction in equilibrium. Free ions with the concentration of 0.15 M NaCl and M of MgCl 2 were added to mimic physiological ionic condition and to neutralize the total system charge. The total number of atoms in the system is ~209K for each model. Conventional MD simulations were performed for all models to examine the thermodynamic stability and potential structural transformation for each starting conformation. The simulations were performed at isothermal isobaric ensemble (NPT ensemble) at body temperature of 310 and 1 atmosphere pressure. The constant temperature was maintained with Langevin dynamics using a damping coefficient of =1 ps -1. Constant pressure ( bar) was maintained using Nosé- Hoover Langevin piston pressure control 2,3 with the Langevin piston period of 200 fs and langevin Piston decay of 100 fs. The time steps were set to 2 fs for the entire simulation and trajectories at every 5000 steps (10 ps) were saved for analysis. Each system was first subjected to conjugate gradient minimization for steps (20 ps). Then, various groups of harmonic constraints were applied for a multi-step simulation to achieve final equilibrium. The system temperature was increased from 0 to 310K at a rate of 1K/ps with the constraints (force constant k=1kcal/mol/ Å 2 ) on protein, ligands and the P atoms of lipids. We then maintained the same constraints for another 5 ns simulation to achieve the equilibrium of water. The constraints on P atoms of lipids were then removed and simulated for another 5 ns to equilibrate the lipid species. The last step was to keep starting constraint (force constant k=5 kcal/mol/å 2 ) on the protein backbone and ligands only, gradually removing the constraints at a speed of 0.5 kcal/mol/å 2 per 500 ps. This step accomplished the maximum packing between the lipid and the protein. In order to observe potential conformation change under conventional unbiased MD simulations, a 200 ns production phase was performed without any constraint on each model. Due to the unconstrained nature of conventional simulations, nucleotide-bound and nucleotide-free Pgp simulations were each performed in triplicate to avoid conclusions based on spurious events. All atomistic MD simulations were performed using NAMD2.9. All MD simulations utilized CHARMM22 with the CMAP correction of all-atom protein force field 4 and CHARMM36 all-atom lipid force field 5. The combination of CHARMM22 with CMAP correction for protein and CHARMM36 for lipid performed well for a similar recent transmembrane protein simulation analysis 6. The force field of ATP 4- was implemented by CHARMM36 7. For both electrostatic and van der Waals calculations, 8 Å switching distance and a 12 Å cutoff distance were used. The particle mesh Ewald (PME) method 8 was employed to calculate the long-range electrostatic interactions with 1.0 grid spacing. The TIP3P explicit water model 9 was utilized. Secondary structural analyses was based on STRIDE algorithm 10. Electrostatic potential calculation used APBS 11 plugin within PyMOL. PyMOL and VMD were used for modeling and molecular image creation.
3 Figure S1. Root-mean-square deviation (RMSD) of protein C α over time. The protein C α RMSDs reveals that the major conformational change took place in the first 100ns of simulation while the protein reach equilibrium states after 100ns. The equilibrate RMSDs with minor changes within 1 Å range indicated that the protein adopt certain predominant meta-stable conformations and oscillate between multiple meta-stable states.
4 Figure S2. Secondary structure component of the apo-pgp run1 (A), run2 (C), run3(e) and MgATP-Pgp run1 (B), run2 (D), run3 (F) over time. Secondary structures were calculated using the STRIDE algorithm. The secondary structure component (alpha helix in blue, turn in red, coil in green, beta sheet in purple, 3_10 helix in orange and beta bridge in light blue) remains consistent during the entire simulation.
5 Figure S3. Secondary structural components of the homology model, X-ray Crystallography structures and equilibrate trajectories from MD simulations. Secondary structures were calculated using the STRIDE algorithm. Homology in blue is the outward-facing mouse-pgp homology model built in this study. Sav1866 (PDB: 2ONJ) in red is the X-ray crystal structure of the outward-facing conformation of Sav1866, mouse-pgp (PDB: 4M1M-A) in green and mouse- Pgp (PDB:4M1M-B) in purple are inward-facing conformations of mouse-pgp from two slightly different X-ray crystal structures. apo-pgp run1-3 in aqua blue, orange and light blue and MgATP- Pgp run1-3 in pink, light green and light purple are conformations from MD simulations. The conformations from MD simulation use average from the last 50 ns (150ns -200ns). The secondary structures of Pgp are conserved (<5% difference) 1) across species between mouse-pgp and Sav1866 and 2) are conserved between conformations: outward-facing conformations (homology model, Sav1866, PDB:2ONJ), intermediate conformations from MD simulations (apo-pgp run1-3, MgATP-Pgp run1-3) and the inward-facing conformations (mouse-pgp PDB: 4M1M-A,B).
6 Figure S4. Per-residue secondary structure probability over time (dynamic secondary structure). Secondary structures were calculated using the STRIDE algorithm. The measurement of dyanmic secondary structure of each residue are based on the last 50 ns equilibrated trajectories averaged over the triplicate simulations for both apo-pgp (A, B) and MgATP-Pgp models (C, D). Panels A and C: first "half" of mouse Pgp. Panels B and D: second "half" of mouse Pgp. Alpha helix in blue, turn in red, coil in green, beta sheet in purple, 3_10 helix in orange and beta bridge in light blue. The figure shows the percentage of time for each residue adopting a certain secondary structure.
7 Figure S5. ATP-binding pocket specifics. The MgATP2- is drawn in licorice with H in white, N in blue, C in cyan, O in red, P in gold and Mg2+ in pink. The protein is shown in cartoon representation with half1 in pink and half2 in ice blue. (A) shows the overview of an MgATP2- complex binding to the mouse P-pg at site 2. (B) shows and H-bond of ATP(Adenine-NH2)-D160/D801(OE). (C) shows multiple H-bonds between the LSGGQ motif with ATP including ATP(H2 )-Q531/Q1176 (OE), ATP(O3 )-Q531/Q1176 (HE), ATP(Oγ)-S528/S1173(HG), ATP(Oγ)-G529/G1174(HN), ATP(Oγ)-G530/G1175(HN). (D) shows five coordination bonds between Mg2+ and the following atoms: E552/E1197 (OE1, OE2), S430/S1073 (OG), and ATP (Oβ, Oγ). (E) shows the H-bond network between MgATP2- and Walker A motif including ATP(Oγ)-S425/S1068(HG), ATP(Oγ)K429/K1072(HZ), ATP(Oβ)-K429/K1072(HZ), ATP(Oβ)-S430/S1073(HN), ATP(Oα)T431/T1074(HG). (F) shows the aromatic interactions between ATP (adenine ring)-y397/y1040 (ring). The H-bonds and coordination bonds are in the range between Å in the initial structure.
8 Figure S6. Electrostatic potential of two halves of mouse Pgp and Sav1866. The crystal structure of Mouse Pgp (PDB: 4M1M) and crystal structure of Sav1866 (PDB: 2ONJ) are drawn in surface using Van der Waals radius. Vacuum electrostatics is calculated using the APBS plugin in PyMOL. Blue/Red color gradient is used to show the charge distribution from the positive end in blue to the negative end in red. The comparison between core regions of TMDs in the Mouse P-pg (A1, A2) and Sav1866 (A) shows an overall neutral region in Mouse P-pg but highly charged region in Sav1866. On the other hand, the interface between two NBDs in mouse P-pg (B1, B2) are more charged than that of Sav1866 (B).
9 Figure S7. Sequence alignment and homology structure prediction. (A) Sequence alignment of a mouse Pgp (first lines) monomer and a Sav1866 (second lines) dimer containing minimal gaps. The linker region (residues ) of mouse Pgp was removed and two tandem copies of Sav1866 were generated. The two resulting amino acid sequences were aligned using the EMBOSS Needle pairwise alignment tool at the EMBL ( using PAM500 similarity matrix, a gap opening penalty of 50 and a gap extension penalty of 10. The resulting alignment had a good amino acid sequence identity of 30%, a similarity of 63%, and only contained 8 gaps (ovals). Due to poor structural conservation of the first and seventh transmembrane helix (TM1, TM7) upon inspection of the inward-facing conformation (PDB: 4M1M), TM1, TM7 and the first extracellular loop (ECL1) of Pgp (residues 1-106) and ECL4 (residue 692 to 747) were removed (top sequence in rectangles). Likewise, Sav1866 TM1 and ECL1 (residues 1-59) were also removed (bottom sequence in rectangles). Four C-terminal residues of Pgp ( ) were removed. (B) The homology model (in orange cartoon) of mouse-pgp has 0.78 Å α-carbon RMSD from the template structure of Sav1866 (in cyan cartoon).
10 Figure S8. Structural alignment and final homology structure. Missing fragments were truncated in the homology model from previous step (Figure S7)(in silver cartoon) using mouse-pgp crystal structure (PDB: 4M1M) (in magenta cartoon) to correctly insert the gap residues. Structure fragments of residue 30 to 107 (A) and residue 681 to 747 (B) were aligned onto the homology structure without disturbing the integrity of the connecting helices. Fragment 368 to 385 was replaced to model missing residue Phe374 (C), Fragment 1009 to 1032 was replaced to model missing residue 1017 (D). Fragment 390 to 410 was replaced to model missing residue Pro398 (E). Fragment 1034 to 1054 was replaced to model residue Pro1044 (F). Fragment 793 to 808 was replaced to model missing residues Asp801 and Asp802 (G). Fragment 1129 to 1139 was replaced to model missing residue 1137, 1138(H).
11 Figure S9. Philosophy of the homology modeling of mouse-pgp outward-facing conformation. The homology model of the outward-facing conformation (C) of mouse Pgp was generated based on both sequence alignment with Sav1866 outward-facing conformation (A) crystal structure (PDB: 2ONJ) in silver and regional structural alignment with inward-facing conformation (B) of the same mouse-pgp (PDB: 4M1M) in magenta.
12 Figure S10. Problematic homology model generated from MODELLER. (I) Problematic homology model generated from MODELLER using the sequence alignment in Figure S2: fragments with problems are labeled in red. (A): incorrect EHL1 sequence residue , the correct loop sequence is (residue 86-92), the fact mouse Pgp has longer TM1 and TM2 than Sav1866. (B) Incorrect EHL3, sequence residue , the correct loop sequence is residue (C) Incorrect orientation after insertion of gap residue Pro398, the loop should point towards the ATP binding pocket. (D) Incorrect orientation after insertion of gap residue Pro1044, the loop should point towards the ATP binding pocket. (E) Poor helix loop orientation after insertion of gap residue Ser1137 Tyr1138. (F) Original helix broken by the poor modeling of gap at Asp801, Asp802. (G)(H) Poor loop orientation after insertion of Phe 374 and Tyr 1017.
13 Figure S11. Morph calculation between inward-facing conformation (A) (PDB: 4M1M) and outward-facing conformation homology model (E). (B), (C), (D) are morph intermediate structures between the two conformations. Secondary structure remained during the conformational change. Secondary structures are calculated using STRIDE algorithm: alpha helix in purple, 3_10 helix in blue, turn in white, coil in cyan, beta sheet in yellow, beta bridge in tan.
14 References. 1 Schiffer, M., Chang, C. H. & Stevens, F. J. The functions of tryptophan residues in membrane proteins. Protein Eng 5, (1992). 2 Martyna, G. J., Klein, M. L. Constant pressure molecular dynamics algorithms. J Chem Phys 101, 4177 (1994). 3 Feller, S. E., Pastor, R. W., Brooks, B. R. Constant pressure molecular dynamics simulation: The Langevin piston method. J. Chem. Phys. 103, 4613 (1995). 4 MacKerell, A. D. J., Bashord, D., M. All- atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B 102, , doi: /jp003919d (1998). 5 Klauda, J. B. et al. Update of the CHARMM all- atom additive force field for lipids: validation on six lipid types. J Phys Chem B 114, , doi: /jp101759q (2010). 6 Andersson, M. et al. Copper- transporting P- type ATPases use a unique ion- release pathway. Nat Struct Mol Biol doi: /nsmb.2721 (2013). 7 Vanommeslaeghe, K. et al. CHARMM general force field: A force field for drug- like molecules compatible with the CHARMM all- atom additive biological force fields. J Comp Chem 31, , doi: /jcc (2010). 8 Tom Darden, D. Y., Lee Pedersen. Particle mesh Ewald: An N [center- dot] log(n) method for Ewald sums in large systems. J Chem Phys 98, 10089, doi: / (1993). 9 Jorgensen, W. L. C., J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liquid water. J Chem Phys 79, (1983). 10 Frishman, D. & Argos, P. Knowledge- based protein secondary structure assignment. Proteins 23, , doi: /prot (1995). 11 Baker, N. A., Sept, D., Joseph, S., Holst, M. J. & McCammon, J. A. Electrostatics of nanosystems: application to microtubules and the ribosome. PNAS 98, , doi: /pnas (2001).
Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains
Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains Fatemeh Khalili-Araghi, Emad Tajkhorshid, Benoît Roux, and Klaus Schulten Department of Physics, Department of Biochemistry,
More informationSupporting Material for. Microscopic origin of gating current fluctuations in a potassium channel voltage sensor
Supporting Material for Microscopic origin of gating current fluctuations in a potassium channel voltage sensor J. Alfredo Freites, * Eric V. Schow, * Stephen H. White, and Douglas J. Tobias * * Department
More informationImproved Resolution of Tertiary Structure Elasticity in Muscle Protein
Improved Resolution of Tertiary Structure Elasticity in Muscle Protein Jen Hsin and Klaus Schulten* Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois
More informationComputational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data
Computational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data ABSTRACT Keyword Lei Shi 1 Advisor: Gianluigi Veglia 1,2 Department of Chemistry 1, & Biochemistry, Molecular
More informationSupplementary Figure S1. Urea-mediated buffering mechanism of H. pylori. Gastric urea is funneled to a cytoplasmic urease that is presumably attached
Supplementary Figure S1. Urea-mediated buffering mechanism of H. pylori. Gastric urea is funneled to a cytoplasmic urease that is presumably attached to HpUreI. Urea hydrolysis products 2NH 3 and 1CO 2
More informationSUPPLEMENTARY INFORMATION
Figure S1. Secondary structure of CAP (in the camp 2 -bound state) 10. α-helices are shown as cylinders and β- strands as arrows. Labeling of secondary structure is indicated. CDB, DBD and the hinge are
More informationMARTINI simulation details
S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained
More informationSupporting Information
Supporting Information Constant ph molecular dynamics reveals ph-modulated binding of two small-molecule BACE1 inhibitors Christopher R. Ellis 1,, Cheng-Chieh Tsai 1,, Xinjun Hou 2, and Jana Shen 1, 1
More informationEffect of intracellular loop 3 on intrinsic dynamics of human β 2 -adrenergic receptor. Ozcan et al.
Effect of intracellular loop 3 on intrinsic dynamics of human β 2 -adrenergic receptor Ozcan et al. Ozcan et al. BMC Structural Biology 2013, 13:29 Ozcan et al. BMC Structural Biology 2013, 13:29 RESEARCH
More informationSubstrate Binding and Formation of an Occluded State in the Leucine Transporter
1600 Biophysical Journal Volume 94 March 2008 1600 1612 Substrate Binding and Formation of an Occluded State in the Leucine Transporter Leyla Celik,* y Birgit Schiøtt, y and Emad Tajkhorshid* *Department
More informationStructure Investigation of Fam20C, a Golgi Casein Kinase
Structure Investigation of Fam20C, a Golgi Casein Kinase Sharon Grubner National Taiwan University, Dr. Jung-Hsin Lin University of California San Diego, Dr. Rommie Amaro Abstract This research project
More informationT H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y. jgp
S u p p l e m e n ta l m at e r i a l jgp Lee et al., http://www.jgp.org/cgi/content/full/jgp.201411219/dc1 T H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y S u p p l e m e n ta l D I S C U S
More informationMedical Research, Medicinal Chemistry, University of Leuven, Leuven, Belgium.
Supporting Information Towards peptide vaccines against Zika virus: Immunoinformatics combined with molecular dynamics simulations to predict antigenic epitopes of Zika viral proteins Muhammad Usman Mirza
More informationComparing crystal structure of M.HhaI with and without DNA1, 2 (PDBID:1hmy and PDBID:2hmy),
Supporting Information 1. Constructing the starting structure Comparing crystal structure of M.HhaI with and without DNA1, 2 (PDBID:1hmy and PDBID:2hmy), we find that: the RMSD of overall structure and
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationStructural and mechanistic insight into the substrate. binding from the conformational dynamics in apo. and substrate-bound DapE enzyme
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 215 Structural and mechanistic insight into the substrate binding from the conformational
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature17991 Supplementary Discussion Structural comparison with E. coli EmrE The DMT superfamily includes a wide variety of transporters with 4-10 TM segments 1. Since the subfamilies of the
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Structure of human carbamoyl phosphate synthetase: deciphering the on/off switch of human ureagenesis Sergio de Cima, Luis M. Polo, Carmen Díez-Fernández, Ana I. Martínez, Javier
More informationSupplementary Information
Supplementary Information Resveratrol Serves as a Protein-Substrate Interaction Stabilizer in Human SIRT1 Activation Xuben Hou,, David Rooklin, Hao Fang *,,, Yingkai Zhang Department of Medicinal Chemistry
More informationType II Kinase Inhibitors Show an Unexpected Inhibition Mode against Parkinson s Disease-Linked LRRK2 Mutant G2019S.
Type II Kinase Inhibitors Show an Unexpected Inhibition Mode against Parkinson s Disease-Linked LRRK2 Mutant G219S. Min Liu@&*, Samantha A. Bender%*, Gregory D Cuny@, Woody Sherman, Marcie Glicksman@ Soumya
More informationElectro-Mechanical Conductance Modulation of a Nanopore Using a Removable Gate
Electro-Mechanical Conductance Modulation of a Nanopore Using a Removable Gate Shidi Zhao a, Laura Restrepo-Pérez b, Misha Soskine c, Giovanni Maglia c, Chirlmin Joo b, Cees Dekker b and Aleksei Aksimentiev
More informationIntroduction to Comparative Protein Modeling. Chapter 4 Part I
Introduction to Comparative Protein Modeling Chapter 4 Part I 1 Information on Proteins Each modeling study depends on the quality of the known experimental data. Basis of the model Search in the literature
More informationPymol Practial Guide
Pymol Practial Guide Pymol is a powerful visualizor very convenient to work with protein molecules. Its interface may seem complex at first, but you will see that with a little practice is simple and powerful
More informationSupplementary Figure 1 Crystal contacts in COP apo structure (PDB code 3S0R)
Supplementary Figure 1 Crystal contacts in COP apo structure (PDB code 3S0R) Shown in cyan and green are two adjacent tetramers from the crystallographic lattice of COP, forming the only unique inter-tetramer
More informationSUPPLEMENTARY INFORMATION
Supplementary Results DNA binding property of the SRA domain was examined by an electrophoresis mobility shift assay (EMSA) using synthesized 12-bp oligonucleotide duplexes containing unmodified, hemi-methylated,
More informationDocking. GBCB 5874: Problem Solving in GBCB
Docking Benzamidine Docking to Trypsin Relationship to Drug Design Ligand-based design QSAR Pharmacophore modeling Can be done without 3-D structure of protein Receptor/Structure-based design Molecular
More informationSupplemental Information for: Characterizing the Membrane-Bound State of Cytochrome P450 3A4: Structure, Depth of Insertion and Orientation
Supplemental Information for: Characterizing the Membrane-Bound State of Cytochrome P450 3A4: Structure, Depth of Insertion and Orientation Javier L. Baylon, Ivan L. Lenov, Stephen G. Sligar and Emad Tajkhorshid
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature12045 Supplementary Table 1 Data collection and refinement statistics. Native Pt-SAD X-ray source SSRF BL17U SPring-8 BL41XU Wavelength (Å) 0.97947 1.07171 Space group P2 1 2 1 2 1 P2
More informationCAP 5510 Lecture 3 Protein Structures
CAP 5510 Lecture 3 Protein Structures Su-Shing Chen Bioinformatics CISE 8/19/2005 Su-Shing Chen, CISE 1 Protein Conformation 8/19/2005 Su-Shing Chen, CISE 2 Protein Conformational Structures Hydrophobicity
More informationSUPPLEMENTARY FIGURES
SUPPLEMENTARY FIGURES Supplementary Figure 1 Protein sequence alignment of Vibrionaceae with either a 40-residue insertion or a 44-residue insertion. Identical residues are indicated by red background.
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature11524 Supplementary discussion Functional analysis of the sugar porter family (SP) signature motifs. As seen in Fig. 5c, single point mutation of the conserved
More informationWhat makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces: Electronic Supplementary
Electronic Supplementary Material (ESI) for Journal of Materials Chemistry B. This journal is The Royal Society of Chemistry 21 What makes a good graphene-binding peptide? Adsorption of amino acids and
More informationβ1 Structure Prediction and Validation
13 Chapter 2 β1 Structure Prediction and Validation 2.1 Overview Over several years, GPCR prediction methods in the Goddard lab have evolved to keep pace with the changing field of GPCR structure. Despite
More informationNature Structural & Molecular Biology: doi: /nsmb Supplementary Figure 1
Supplementary Figure 1 Chemical structure of LPS and LPS biogenesis in Gram-negative bacteria. a. Chemical structure of LPS. LPS molecule consists of Lipid A, core oligosaccharide and O-antigen. The polar
More informationMolecular dynamics simulations of anti-aggregation effect of ibuprofen. Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov
Biophysical Journal, Volume 98 Supporting Material Molecular dynamics simulations of anti-aggregation effect of ibuprofen Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov Supplemental
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature11085 Supplementary Tables: Supplementary Table 1. Summary of crystallographic and structure refinement data Structure BRIL-NOP receptor Data collection Number of crystals 23 Space group
More informationDon t forget to bring your MD tutorial. Potential Energy (hyper)surface
Don t forget to bring your MD tutorial Lab session starts at 1pm You will have to finish an MD/SMD exercise on α-conotoxin in oxidized and reduced forms Potential Energy (hyper)surface What is Force? Energy
More informationSupporting Information
Supporting Information Ottmann et al. 10.1073/pnas.0907587106 Fig. S1. Primary structure alignment of SBT3 with C5 peptidase from Streptococcus pyogenes. The Matchmaker tool in UCSF Chimera (http:// www.cgl.ucsf.edu/chimera)
More informationComputational Predictions of 1-Octanol/Water Partition Coefficient for Imidazolium based Ionic Liquids.
Computational Predictions of 1-Octanol/Water Partition Coefficient for Imidazolium based Ionic Liquids. Ganesh Kamath,* a Navendu Bhatnagar b, Gary A. Baker a, Sheila N. Baker c and Jeffrey J. Potoff b
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationBuilding a Homology Model of the Transmembrane Domain of the Human Glycine α-1 Receptor
Building a Homology Model of the Transmembrane Domain of the Human Glycine α-1 Receptor Presented by Stephanie Lee Research Mentor: Dr. Rob Coalson Glycine Alpha 1 Receptor (GlyRa1) Member of the superfamily
More informationExamples of Protein Modeling. Protein Modeling. Primary Structure. Protein Structure Description. Protein Sequence Sources. Importing Sequences to MOE
Examples of Protein Modeling Protein Modeling Visualization Examination of an experimental structure to gain insight about a research question Dynamics To examine the dynamics of protein structures To
More informationSupplementary Figures for Tong et al.: Structure and function of the intracellular region of the plexin-b1 transmembrane receptor
Supplementary Figures for Tong et al.: Structure and function of the intracellular region of the plexin-b1 transmembrane receptor Figure S1. Plexin-B1 GAP segments are homologous to RasGAPs. Sequence alignment
More informationChowdhury et al. PoreDesigner for tuning solute selectivity in a robust and highly permeable outer membrane pore.
Chowdhury et al. PoreDesigner for tuning solute selectivity in a robust and highly permeable outer membrane pore. Supplementary Methods Isolation of water wires We isolated water-wire trajectories from
More informationNature Structural & Molecular Biology doi: /nsmb Supplementary Figure 1. CRBN binding assay with thalidomide enantiomers.
Supplementary Figure 1 CRBN binding assay with thalidomide enantiomers. (a) Competitive elution assay using thalidomide-immobilized beads coupled with racemic thalidomide. Beads were washed three times
More informationStructural Characterization of Two Metastable ATP-Bound States of P-Glycoprotein
Structural Characterization of Two Metastable ATP-Bound States of P-Glycoprotein Megan L. O Mara 1,2 *, Alan E. Mark 1,3 1 School of Chemistry and Molecular Biosciences (SCMB), The Institute for Molecular
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NCHEM.1821 Cloud-based simulations on Google Exacycle reveal ligand-modulation of GPCR activation pathways Kai J. Kohlhoff 1,4 *, Diwakar Shukla 1,2, Morgan Lawrenz 2, Gregory R. Bowman 2,
More informationSupporting Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Supporting Information Lipid molecules can induce an opening of membrane-facing
More informationCan a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?
Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,
More informationAmino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information
Amino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information Grzegorz Nawrocki and Marek Cieplak Institute of Physics, Polish Academy of Sciences,
More informationPacking of Secondary Structures
7.88 Lecture Notes - 4 7.24/7.88J/5.48J The Protein Folding and Human Disease Professor Gossard Retrieving, Viewing Protein Structures from the Protein Data Base Helix helix packing Packing of Secondary
More informationNitrogenase MoFe protein from Clostridium pasteurianum at 1.08 Å resolution: comparison with the Azotobacter vinelandii MoFe protein
Acta Cryst. (2015). D71, 274-282, doi:10.1107/s1399004714025243 Supporting information Volume 71 (2015) Supporting information for article: Nitrogenase MoFe protein from Clostridium pasteurianum at 1.08
More informationBiophysics 490M Project
Biophysics 490M Project Dan Han Department of Biochemistry Structure Exploration of aa 3 -type Cytochrome c Oxidase from Rhodobacter sphaeroides I. Introduction: All organisms need energy to live. They
More informationComputational Structural Biology and Molecular Simulation. Introduction to VMD Molecular Visualization and Analysis
Computational Structural Biology and Molecular Simulation Introduction to VMD Molecular Visualization and Analysis Emad Tajkhorshid Department of Biochemistry, Beckman Institute, Center for Computational
More informationNature Structural & Molecular Biology: doi: /nsmb Supplementary Figure 1
Supplementary Figure 1 Crystallization. a, Crystallization constructs of the ET B receptor are shown, with all of the modifications to the human wild-type the ET B receptor indicated. Residues interacting
More informationStructural Insights from Molecular Dynamics. Simulations of Tryptophan 7-Halogenase and
Supporting Information Structural Insights from Molecular Dynamics Simulations of Tryptophan 7-Halogenase and Tryptophan 5-halogenase Jon Ainsley 1, Adrian J. Mulholland 2, Gary W. Black 1, Olivier Sparagano
More informationSupplementary figure 1. Comparison of unbound ogm-csf and ogm-csf as captured in the GIF:GM-CSF complex. Alignment of two copies of unbound ovine
Supplementary figure 1. Comparison of unbound and as captured in the GIF:GM-CSF complex. Alignment of two copies of unbound ovine GM-CSF (slate) with bound GM-CSF in the GIF:GM-CSF complex (GIF: green,
More informationSecondary Structure. Bioch/BIMS 503 Lecture 2. Structure and Function of Proteins. Further Reading. Φ, Ψ angles alone determine protein structure
Bioch/BIMS 503 Lecture 2 Structure and Function of Proteins August 28, 2008 Robert Nakamoto rkn3c@virginia.edu 2-0279 Secondary Structure Φ Ψ angles determine protein structure Φ Ψ angles are restricted
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature10955 Supplementary Figures Supplementary Figure 1. Electron-density maps and crystallographic dimer structures of the motor domain. (a f) Stereo views of the final electron-density maps
More informationProtein structure (and biomolecular structure more generally) CS/CME/BioE/Biophys/BMI 279 Sept. 28 and Oct. 3, 2017 Ron Dror
Protein structure (and biomolecular structure more generally) CS/CME/BioE/Biophys/BMI 279 Sept. 28 and Oct. 3, 2017 Ron Dror Please interrupt if you have questions, and especially if you re confused! Assignment
More informationMolecular Modeling Lecture 7. Homology modeling insertions/deletions manual realignment
Molecular Modeling 2018-- Lecture 7 Homology modeling insertions/deletions manual realignment Homology modeling also called comparative modeling Sequences that have similar sequence have similar structure.
More informationMechanical Proteins. Stretching imunoglobulin and fibronectin. domains of the muscle protein titin. Adhesion Proteins of the Immune System
Mechanical Proteins F C D B A domains of the muscle protein titin E Stretching imunoglobulin and fibronectin G NIH Resource for Macromolecular Modeling and Bioinformatics Theoretical Biophysics Group,
More informationMolecular Dynamics Studies of Human β-glucuronidase
American Journal of Applied Sciences 7 (6): 83-88, 010 ISSN 1546-939 010Science Publications Molecular Dynamics Studies of Human β-lucuronidase Ibrahim Ali Noorbatcha, Ayesha Masrur Khan and Hamzah Mohd
More informationSupplementary information for cloud computing approaches for prediction of ligand binding poses and pathways
Supplementary information for cloud computing approaches for prediction of ligand binding poses and pathways Morgan Lawrenz 1, Diwakar Shukla 1,2 & Vijay S. Pande 1,2 1 Department of Chemistry, Stanford
More informationModeling and Dynamics of the Inward-Facing State of a Na + /Cl 2 Dependent Neurotransmitter Transporter Homologue
Modeling and Dynamics of the Inward-Facing State of a Na + /Cl 2 Dependent Neurotransmitter Transporter Homologue Saher Afshan Shaikh 1, Emad Tajkhorshid 1,2 * 1 Department of Biochemistry and Beckman
More informationStructure and evolution of the spliceosomal peptidyl-prolyl cistrans isomerase Cwc27
Acta Cryst. (2014). D70, doi:10.1107/s1399004714021695 Supporting information Volume 70 (2014) Supporting information for article: Structure and evolution of the spliceosomal peptidyl-prolyl cistrans isomerase
More informationMajor Types of Association of Proteins with Cell Membranes. From Alberts et al
Major Types of Association of Proteins with Cell Membranes From Alberts et al Proteins Are Polymers of Amino Acids Peptide Bond Formation Amino Acid central carbon atom to which are attached amino group
More informationPotential Energy (hyper)surface
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = " d dx U(x) Conformation
More informationThe Molecular Dynamics Method
H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular
More informationSupplementary information
Supplementary information The structural basis of modularity in ECF-type ABC transporters Guus B. Erkens 1,2, Ronnie P-A. Berntsson 1,2, Faizah Fulyani 1,2, Maria Majsnerowska 1,2, Andreja Vujičić-Žagar
More informationSupporting Information
S-1 Supporting Information Flaviviral protease inhibitors identied by fragment-based library docking into a structure generated by molecular dynamics Dariusz Ekonomiuk a, Xun-Cheng Su b, Kiyoshi Ozawa
More informationRoutine access to millisecond timescale events with accelerated molecular dynamics
Routine access to millisecond timescale events with accelerated molecular dynamics Levi C.T. Pierce, Romelia Salomon-Ferrer, Cesar Augusto F. de Oliveira #, J. Andrew McCammon #, Ross C. Walker * SUPPORTING
More informationTransmembrane Domains (TMDs) of ABC transporters
Transmembrane Domains (TMDs) of ABC transporters Most ABC transporters contain heterodimeric TMDs (e.g. HisMQ, MalFG) TMDs show only limited sequence homology (high diversity) High degree of conservation
More informationSUPPLEMENTARY INFORMATION
Fig. 1 Influences of crystal lattice contacts on Pol η structures. a. The dominant lattice contact between two hpol η molecules (silver and gold) in the type 1 crystals. b. A close-up view of the hydrophobic
More informationHands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C. Course Details
Hands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C Emad Tajkhorshid Center for Computational Biology and Biophysics Email: emad@life.uiuc.edu or tajkhors@uiuc.edu
More informationThe Dominant Interaction Between Peptide and Urea is Electrostatic in Nature: A Molecular Dynamics Simulation Study
Dror Tobi 1 Ron Elber 1,2 Devarajan Thirumalai 3 1 Department of Biological Chemistry, The Hebrew University, Jerusalem 91904, Israel 2 Department of Computer Science, Cornell University, Ithaca, NY 14853
More informationModel Mélange. Physical Models of Peptides and Proteins
Model Mélange Physical Models of Peptides and Proteins In the Model Mélange activity, you will visit four different stations each featuring a variety of different physical models of peptides or proteins.
More informationTime-dependence of key H-bond/electrostatic interaction distances in the sirna5-hago2 complexes... Page S14
Supporting Information Probing the Binding Interactions between Chemically Modified sirnas and Human Argonaute 2 Using Microsecond Molecular Dynamics Simulations S. Harikrishna* and P. I. Pradeepkumar*
More informationMolecular modeling. A fragment sequence of 24 residues encompassing the region of interest of WT-
SUPPLEMENTARY DATA Molecular dynamics Molecular modeling. A fragment sequence of 24 residues encompassing the region of interest of WT- KISS1R, i.e. the last intracellular domain (Figure S1a), has been
More informationCharacterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi
Characterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi NCSA Blue Waters Symposium for Petascale Science and Beyond Sunriver, Oregon May 11, 2015 Outline Introduction
More informationMolecular dynamics simulations of a single stranded (ss) DNA
Molecular dynamics simulations of a single stranded (ss) DNA Subhasish Chatterjee 1, Bonnie Gersten 1, Siddarth Thakur 2, Alexander Burin 2 1 Department of Chemistry, Queens College and the Graduate Center
More informationSupporting Information
Supporting Information The Predicted Ensemble of Low-Energy Conformations of Human Somatostatin Receptor Subtype 5 and the Binding of Antagonists Sijia S. Dong, [a] Ravinder Abrol, [a, b] and William A.
More informationOxygen Binding in Hemocyanin
Supporting Information for Quantum Mechanics/Molecular Mechanics Study of Oxygen Binding in Hemocyanin Toru Saito and Walter Thiel* Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470
More informationProtein Dynamics. The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron.
Protein Dynamics The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron. Below is myoglobin hydrated with 350 water molecules. Only a small
More informationSUPPLEMENTARY MATERIAL. Supplementary material and methods:
Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 2015 SUPPLEMENTARY MATERIAL Supplementary material and methods: - Computational
More informationComputational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College
Computational Studies of the Photoreceptor Rhodopsin Scott E. Feller Wabash College Rhodopsin Photocycle Dark-adapted Rhodopsin hn Isomerize retinal Photorhodopsin ~200 fs Bathorhodopsin Meta-II ms timescale
More informationSupplementary Figure 3 a. Structural comparison between the two determined structures for the IL 23:MA12 complex. The overall RMSD between the two
Supplementary Figure 1. Biopanningg and clone enrichment of Alphabody binders against human IL 23. Positive clones in i phage ELISA with optical density (OD) 3 times higher than background are shown for
More informationSupplemental Information for: Mechanical properties of a complete microtubule revealed through molecular dynamics simulation
Supplemental Information for: Mechanical properties of a complete microtubule revealed through molecular dynamics simulation avid Wells and Aleksei Aksimentiev Water placement and histidine protonation
More informationFigure 1. Molecules geometries of 5021 and Each neutral group in CHARMM topology was grouped in dash circle.
Project I Chemistry 8021, Spring 2005/2/23 This document was turned in by a student as a homework paper. 1. Methods First, the cartesian coordinates of 5021 and 8021 molecules (Fig. 1) are generated, in
More informationBasics of protein structure
Today: 1. Projects a. Requirements: i. Critical review of one paper ii. At least one computational result b. Noon, Dec. 3 rd written report and oral presentation are due; submit via email to bphys101@fas.harvard.edu
More informationJournal of Pharmacology and Experimental Therapy-JPET#172536
A NEW NON-PEPTIDIC INHIBITOR OF THE 14-3-3 DOCKING SITE INDUCES APOPTOTIC CELL DEATH IN CHRONIC MYELOID LEUKEMIA SENSITIVE OR RESISTANT TO IMATINIB Manuela Mancini, Valentina Corradi, Sara Petta, Enza
More informationModeling Biological Systems Opportunities for Computer Scientists
Modeling Biological Systems Opportunities for Computer Scientists Filip Jagodzinski RBO Tutorial Series 25 June 2007 Computer Science Robotics & Biology Laboratory Protein: πρώτα, "prota, of Primary Importance
More informationMolecular Dynamics Simulations of the Mammalian Glutamate Transporter EAAT3
Molecular Dynamics Simulations of the Mammalian Glutamate Transporter EAAT3 Germano Heinzelmann, Serdar Kuyucak* School of Physics, University of Sydney, NSW, Australia Abstract Excitatory amino acid transporters
More informationW tot (z) = W (z) k B T ln e i q iφ mp(z i )/k B T, [1]
Supporting Text Simulation Details. Simulations were carried out with the program charmm (1) using the PARAM27 force field with standard lipid (2), protein (3), and TIP3P water model (4) parameters. The
More informationTable S1. Overview of used PDZK1 constructs and their binding affinities to peptides. Related to figure 1.
Table S1. Overview of used PDZK1 constructs and their binding affinities to peptides. Related to figure 1. PDZK1 constru cts Amino acids MW [kda] KD [μm] PEPT2-CT- FITC KD [μm] NHE3-CT- FITC KD [μm] PDZK1-CT-
More informationATP Binding Enables Broad Antibiotic Selectivity of Aminoglycoside Phosphotransferase(3 )-IIIa: An Elastic Network Analysis
doi:10.1016/j.jmb.2011.03.061 J. Mol. Biol. (2011) 409, 450 465 Contents lists available at www.sciencedirect.com Journal of Molecular Biology journal homepage: http://ees.elsevier.com.jmb ATP Binding
More informationApplications of Molecular Dynamics
June 4, 0 Molecular Modeling and Simulation Applications of Molecular Dynamics Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru
More informationSupplementary Figures:
Supplementary Figures: Supplementary Figure 1: The two strings converge to two qualitatively different pathways. A) Models of active (red) and inactive (blue) states used as end points for the string calculations
More informationDynamic hydration valve controlled ion permeability and
1 Dynamic hydration valve controlled ion permeability and stability of NaK channel Rong Shen 1, *, Wanlin Guo 1, * & Wenyu Zhong 1 1 Institute of Nano Science, Nanjing University of Aeronautics and Astronautics,
More informationSupporting Information for. Hydrogen Bonding Structure at Zwitterionic. Lipid/Water Interface
Supporting Information for Hydrogen Bonding Structure at Zwitterionic Lipid/Water Interface Tatsuya Ishiyama,, Daichi Terada, and Akihiro Morita,, Department of Applied Chemistry, Graduate School of Science
More information