Introduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. Th
|
|
- Thomasina Watson
- 5 years ago
- Views:
Transcription
1 Abstract When conductive, gramicidin monomers are linked by six hydrogen bonds. To understand the details of dissociation and how the channel transits from a state with 6H bonds to ones with 4H bonds or 2H bonds (both believed non-conductive), we generated two-dimensional topographic maps of the total non-bonded inter-monomer interaction energy as functions of the separation distance and rotation angle between monomers. Calculations for four distinct gramicidin structures (1MAG, 1JNO, 1JO3 & 1JO4 pdb files) led to complex, but qualitatively similar, energy maps with three wells (6H, 4H & 2H states), separated by high energy barriers. With rigid body monomer motion of 1MAG fully treated (including tilt and lateral displacement), two distinct dissociative pathways were observed in determinations of both energy and potential of mean force. Along both, initial motion is lateral, leading to a closed pore (nonconducting) structure. At a critical displacement on path 1 the monomers abruptly separate axially with simultaneous lateral relaxation and formation of an open, putatively conducting pore, stabilized by 4H bonds. Then dissociation occurs by further lateral motion, maintaining the 4H bond separation. On path 2 monomer movement is purely lateral, maintaining the 6H bond separation; no open pore forms.
2 Introduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. The assembly is stabilized by six junctional hydrogen bonds; the dimer lifetimes are of the order of seconds [1]. Little is known about the kinetic mechanism for ga dissociation in membranes. Channel gating (closing and opening) is presumed due to the dissociation and association of the monomers with the closing transition triggered by breaking the dimer s stabilizing hydrogen bonds. As the monomers can rotate relative to each other, formation of a dimer with four or two hydrogen bonds is also possible. The channel is believed to be no longer conductive when the monomers separate by ~1.6 Å [2]. To study the detailed kinetics of ga dissociation we calculate the non-bonded interaction energy between the monomers allowing for free spatial and rotational motion. The intermonomer separation distance, the radial displacement and the effect of the tilt angle on the interaction energy are analyzed. The potential of mean force (PMF) along the minimum energy pathway for dimer dissociation is calculated using the free energy perturbation method.
3 Computational Model The non-bonded interaction energy between the ga monomers is calculated using our MCICP (Monte Carlo Ion Channel Proteins) code with the following simplifying assumptions [3]: 1) Bulk water regions are treated as continua with ε = 80. The monomers are immersed in a low dielectric (ε = 1) membrane slab, somewhat overestimating electrostatic effects. The reaction field is treated by the method of images; van der Waals and electrostatic interactions are computed with no cutoff. 2) We use partial charges and van der Waals parameters of the CHARMM22 allhydrogen force field. The crystallographic structure of the ga channel [1] (pdb entry 1MAG resolved using Solid State NMR) was used. The ga structure contains 552 explicit protein atoms including all hydrogens. The two ga monomers are held rigid. They are permitted free 3D relative rotation and translation. The Kinetic Monte Carlo Reaction Path Following (KMCRPF) technique [4] is used to determine the lowest energy pathways and the reaction coordinates for dissociating the dimer into monomers. Monomer B is held fixed and monomer A (Figure 1) allowed translational (three spatial degrees of freedom) and rotational (rotation and tilt angles) motion. To monitor the reaction coordinate only increases in the rotational angle φ are permitted [4].
4 Dissociation Kinetics of the ga Dimer Figure 1 illustrates the molecular model. The lowest energy pathway was determined by rotating monomer A along the φ coordinate from the 6H state to the 4H state, then to the 2H state and finally to the 6H state. The reaction coordinates (separation distance and radial displacement between the monomers, tilt angle of monomer A) corresponding to the lowest energy path were determined simultaneously. Figure 2 presents an energy map as a function of the separation distance and rotation angle between the monomers (the d,φ-map). Radial displacement and tilting of monomer A is forbidden. The map is calculated by the grid method, with separation distance steps of 0.05 Å and angular steps of 0.5º. There are three clear energy wells corresponding to 6H(ydrogen bond), 4H and 2H states. Large energy barriers separate the wells. For the 1MAG structure used here, the energy, separation distance and rotation angle differences between the 6H, 4H and 2H states are tabulated. E [kt] E [kt] d [Å] d [Å] φ [deg] φ [deg] 6H-4H 4H-2H 6H-4H 4H-2H 6H-4H 4H-2H
5 Figure 3 shows the lowest energy contours on the d,φ-map as a function of φ. The results are presented for the four experimental ga structures available, each corresponding to a different conformation and structure. 1MAG and 1JNO are for gramicidin A. 1JO3 is gramicidin B. 1JO4 is gramicidin C. Nonetheless, all have similar lowest energy pathways for dissociation into monomers. The profile shape is conserved. Differences in energy maxima and minima most likely reflect different conformational states of the side chains. These profiles suggest that dissociation kinetics in gramicidin is basically independent of the structural details. Results are unaltered if CO and NH bond bending and stretching is included; the inter-monomer interaction energy does not change (calculations, not shown, were performed for the 1MAG structure). Figure 4 illustrates the effect of the membrane dielectric constant, ε, on the energy profile for axial displacement of monomer A along the Z-axis, starting from the 6H state; rotation is forbidden. Increasing ε greatly alters the 6H energy well depth, with saturation at ε ~10. The 6H energy well (~30 kt) is conserved when all partial charges on the protein atoms are turned off. The inter-monomer van der Waals interaction contributes significantly to the stabilization energy.
6 Figure 5 illustrates the lowest energy profiles corresponding to the reaction coordinate for dimer dissociation. The results of several Monte Carlo (MC) runs are illustrated. Monomer A is permitted radial displacement and tilt. The importance of these additional degrees of freedom is evident from the difference between the MC curves and the d,φ-curve (in black). The energy barriers decreased significantly. The barrier separating 6H and 4H wells is ~40 kt. The well depth difference is ~18 kt. Inter-monomer hydrogen bonds break much more easily if monomer A undergoes screw type motion with simultaneous lateral displacement rather than direct axial separation. At the peak of the barrier between the 6H and 4H states the reaction pathway bifurcates. Path 1 leads to formation of the 4H state. There is switching between hydrogen bonds at the inter-monomer junction (a rotational shift) with an abrupt discontinuity in the separation distance (~1.6 Å, see Figure 2). Two hydrogen bonds break in transiting from the 6H to the 4H state. Along path 2 the monomers dissociate directly from the native 6H state by coupled screw type and lateral motion with neither a separation distance discontinuity nor a rotational shift between hydrogen bonds. The probability for realizing these paths depends on the tilt permitted monomer A, the specific membrane, etc. For φ > 150º (escape from the 4H state) the energy fluctuates between 40 and 60 kt. The barriers between 4H and 2H wells and 2H and 6H wells disappear (compare the MC curves with the d,φ-curve).
7 Figure 6 illustrates the angular dependence of the separation distance between monomers. Zero separation corresponds to the crystal state. For the d,φ-path (no tilt allowed) monomers separate ~4 Å in transition from the 6H state to the 4H state. This arises because of a large electrostatic barrier since for φ ~60º COs and NHs from the individual monomers are in opposition. (HN)VAL 1 (A) abuts (HN)ALA 5 (B) (CO)VAL 1 (A) abuts (CO)ALA 3 (B) (HN)ALA 3 (A) abuts (HN)ALA 3 (B) (CO)ALA 3 (A) abuts (CO)VAL 1 (B) (HN)ALA 5 (A) abuts (HN)VAL 1 (B) The separation distance is very different when radial displacement and tilt are allowed (MC curves). The monomers do not separate for φ 75º and the separation distance fluctuates near zero. This results from strong inter-monomer interaction due to the six native (original) inter-monomer hydrogen bonds, which are not easily broken. Along path 1 the separation then abruptly increases by ~1.6 Å, corresponding to formation of the 4H state. Here there is a rotational shift between the hydrogen bonds at the intermonomer junction. Along path 2 the monomers do not separate; they are displaced radially. Dissociation occurs from the native 6H state with no rotational shift between hydrogen bonds. As the model treated here doesn t account for membrane elasticity, at large radial displacement (> 4 Å) the monomers slip by one another (a non-physical artifact) and the separation distance between them becomes negative.
8 Figure 7 illustrates the radial displacement between the monomers along the dissociative path. The inter-monomer junction is displaced radially by ~3.5 Å at the peak of the barrier separating the 6H and 4H wells; here φ ~75º and the dimer pore is fully occluded. For larger φ along path 1 the monomers radial displacement again approaches zero (forming the 4H state) and the pore reopens fully. Dissociation from the 4H state takes place radially. Beyond φ ~75º along path 2 separation occurs laterally without reforming an open pore. Figure 8 illustrates the tilt angle of monomer A. The tilt angle fluctuates by ~8º on both paths. Depending on the tilt restrictions either path 1 or 2 can be realized. If tilt is forbidden or near zero path 1 occurs frequently (nine out of ten MC runs). When tilt is freely variable direct dissociation from the 6H state (path 2) is most probable. The presence of the membrane greatly modifies monomer behavior [5]. Figure 9 illustrates the potential of mean force (PMF) along the reaction coordinate for dimer dissociation, calculated by free energy perturbation [6]. Unlike computational alchemy, which proceeds via non-physical intermediate states, the PMF is calculated for a physically achievable dissociation process. The reaction coordinate (rotation angle φ) was increased from 0 to 360º in 0.3º steps. In the main the PMF mimics the total energy profiles along both reaction pathways. In the 6H region there is no difference between the PMF and total energy profiles. The PMF barriers separating the 6H and 4H states are lower than the total energy barriers.
9 Conclusions and Observations To overcome a large electrostatic barrier separating the 6H and 4H states the monomers are displaced radially. Axial separation is energetically far more expensive than lateral motion. The two possible reaction pathways are observed for dissociation. On path 1 the 4H state forms (a fully open pore). Here the monomers undergo an abrupt displacement in the axial separation by ~1.6 Å. On path 2 the monomers dissociate laterally and their separation remains near zero. In transiting from the 6H state to the 4H state ionic conduction is interrupted as the monomers undergo a radial shift at the inter-monomer junction by 3 Å with no axial separation. Radial displacement fluctuations in the 6H state may provide a rationale for fast closure events (flickers) at sub-millisecond times [7]. Our results demonstrate that the transition state for channel dissociation (breaking some hydrogen bonds at the dimer junction) may occur by lateral displacement, without axially separating the monomers by 1.6 Å [2]. Introducing membrane elasticity into the model eliminates the non-physical artifact of the monomers slipping by one another at large radial displacements. Membrane elasticity also substantially stabilizes formation of 4H and even 2H states [5]. Along the reaction pathways the PMF mimics the total energy profile. However, the PMF barrier between 6H and 4H states is less than the total energy barrier.
10 Acknowledgements Work supported by a grant from the National Institutes of Health, GM References 1. T.C. Vogt et al., Biochemistry 31: (1992) 2. J.A. Lundbæk & O.S. Andersen, Biophys. J. 76: (1999) 3. P.C. Jordan & G.V. Miloshevsky, Biophys. J. 82:199a (2002) 4. G.V. Miloshevsky & P.C. Jordan, Biophys. J. 84: 412a (2003) 5. M.B. Partenskii et al., Biophys. J. 84: 2520a (2003) 6. A.R. Leach, Molecular modelling: principles and applications, Harlow, England: Longman, K.M. Armstrong & S. Cukierman, Biophys. J. 82: (2002)
11 Y ETA φ d X R θ Z Monomer B FOR Monomer A Figure 1. Two ga monomers (A & B) are illustrated as separate mobile structural elements. The helices are shown in blue and green. Atoms of the FOR and ETA residues are displayed in their conventional colors. Arrows demonstrate the monomers separation distance, radial displacement, tilt angle, and rotation direction.
12 Figure 2. Contour map of the total (electrostatic + vdw) energy as a function of the separation distance d and rotation angle φ between ga monomers. The energy wells corresponding to the states with 6H, 4H and 2H hydrogen bonds are illustrated.
13 Energy difference, (kt) H 4H Angle ϕ, (deg) Figure 3. The energy profiles corresponding to the lowest energy path on the d,φ-map between 6H, 4H and 2H states as a function of the rotation angle φ between monomers. Profiles are shown for four experimental gramicidin structures. 2H 1MAG 1JNO 1JO3 1JO4
14 0 Energy difference, (kt) H ε = 1 ε = 4 no charges Separation Distance d, (Å) Figure 4. Energy profiles as a function of the axial monomer separation distance (movement along the Z-axis directly from the 6H state). Profiles are shown for ε = 1 and 4, and when all partial charges on the protein atoms are excluded (ε ).
15 Energy difference, (kt) 60 path H 4H path 1 (d,ϕ) path MC run 1 MC run 2 MC run 3 MC run 4 MC run 5 MC run Figure 5. Energy profiles corresponding to the lowest energy pathways for dimer dissociation. Two paths are indicated. Path 1 corresponds to formation of the 4H state. Along path 2 monomers dissociate directly from the 6H state by radial displacement. The black curve is the energy profile in d,φ-space (tilt and lateral motion excluded). 2H Angle ϕ, (deg)
16 Separation Distance d, (Å) H (d,ϕ) path MC run 1 MC run 2 MC run 3 path 2 4H path 1 MC run 4 MC run 5 MC run Figure 6. Separation distance between the ga monomers along the lowest energy pathways for paths 1 and 2 and for the d,φ-path. The separation distance discontinuity occurs at 75º on the path to the 4H state. 2H Angle ϕ, (deg)
17 Radial Displacement R, (Å) 10 MC run 4 MC run 5 MC run H path 1 path 2 4H MC run 1 MC run 2 MC run Angle ϕ, (deg) Figure 7. Radial displacement between monomers on the two paths is illustrated. When the 4H state forms (path 1) a conductive ga dimer reforms (R near zero). Along path 2 the monomers dissociate directly from the 6H state.
18 12 Tilt Angle θ, (deg) MC run 1 MC run 2 MC run Angle ϕ, (deg) Figure 8. Tilt angle of monomer A is illustrated. The tilt profiles are similar for paths 1 and 2. The tilt angle fluctuates around 8º.
19 Energy difference, (kt) H 4H Energy, path 1 PMF, path 1 Energy, path 2 PMF, path Angle ϕ, (deg) Figure 9. Total energy and potential of mean force (PMF) profiles along the reaction coordinate (rotation angle φ) for the dimer to monomer reaction. Profiles are illustrated for pathways 1 and 2.
Gating Gramicidin Channels in Lipid Bilayers: Reaction Coordinates and the Mechanism of Dissociation
92 Biophysical Journal Volume 86 January 2004 92 104 Gating Gramicidin Channels in Lipid Bilayers: Reaction Coordinates and the Mechanism of Dissociation Gennady V. Miloshevsky and Peter C. Jordan Department
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationMulti-Scale Hierarchical Structure Prediction of Helical Transmembrane Proteins
Multi-Scale Hierarchical Structure Prediction of Helical Transmembrane Proteins Zhong Chen Dept. of Biochemistry and Molecular Biology University of Georgia, Athens, GA 30602 Email: zc@csbl.bmb.uga.edu
More informationPolypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation
Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation Jakob P. Ulmschneider and William L. Jorgensen J.A.C.S. 2004, 126, 1849-1857 Presented by Laura L. Thomas and
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationDihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769
Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral
More informationAtomic and molecular interaction forces in biology
Atomic and molecular interaction forces in biology 1 Outline Types of interactions relevant to biology Van der Waals interactions H-bond interactions Some properties of water Hydrophobic effect 2 Types
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 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 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 informationProtein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a
Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2004 How do we go from an unfolded polypeptide chain to a compact folded protein? (Folding of thioredoxin, F. Richards) Structure - Function
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationLecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions Van der Waals Interactions
More information3. An Introduction to Molecular Mechanics
3. An Introduction to Molecular Mechanics Introduction When you use Chem3D to draw molecules, the program assigns bond lengths and bond angles based on experimental data. The program does not contain real
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationarxiv: v1 [cond-mat.soft] 22 Oct 2007
Conformational Transitions of Heteropolymers arxiv:0710.4095v1 [cond-mat.soft] 22 Oct 2007 Michael Bachmann and Wolfhard Janke Institut für Theoretische Physik, Universität Leipzig, Augustusplatz 10/11,
More informationLecture 11: Protein Folding & Stability
Structure - Function Protein Folding: What we know Lecture 11: Protein Folding & Stability 1). Amino acid sequence dictates structure. 2). The native structure represents the lowest energy state for a
More informationProtein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall Protein Folding: What we know. Protein Folding
Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2003 Structure - Function Protein Folding: What we know 1). Amino acid sequence dictates structure. 2). The native structure represents
More informationUnfolding CspB by means of biased molecular dynamics
Chapter 4 Unfolding CspB by means of biased molecular dynamics 4.1 Introduction Understanding the mechanism of protein folding has been a major challenge for the last twenty years, as pointed out in the
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 informationChapter 6- An Introduction to Metabolism*
Chapter 6- An Introduction to Metabolism* *Lecture notes are to be used as a study guide only and do not represent the comprehensive information you will need to know for the exams. The Energy of Life
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 informationProteins polymer molecules, folded in complex structures. Konstantin Popov Department of Biochemistry and Biophysics
Proteins polymer molecules, folded in complex structures Konstantin Popov Department of Biochemistry and Biophysics Outline General aspects of polymer theory Size and persistent length of ideal linear
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationThe Potassium Ion Channel: Rahmat Muhammad
The Potassium Ion Channel: 1952-1998 1998 Rahmat Muhammad Ions: Cell volume regulation Electrical impulse formation (e.g. sodium, potassium) Lipid membrane: the dielectric barrier Pro: compartmentalization
More informationPotassium channel gating and structure!
Reading: Potassium channel gating and structure Hille (3rd ed.) chapts 10, 13, 17 Doyle et al. The Structure of the Potassium Channel: Molecular Basis of K1 Conduction and Selectivity. Science 280:70-77
More informationFree Energy Calculations of Gramicidin Dimer Dissociation
pubs.acs.org/jpcb Free Energy Calculations of Gramicidin Dimer Dissociation Surajith N. Wanasundara, Vikram Krishnamurthy,*, and Shin-Ho Chung Department of Electrical and Computer Engineering, University
More information3. An Introduction to Molecular Mechanics
3. An Introduction to Molecular Mechanics Introduction When you use Chem3D to draw molecules, the program assigns bond lengths and bond angles based on experimental data. The program does not contain real
More informationWhat binds to Hb in addition to O 2?
Reading: Ch5; 158-169, 162-166, 169-174 Problems: Ch5 (text); 3,7,8,10 Ch5 (study guide-facts); 1,2,3,4,5,8 Ch5 (study guide-apply); 2,3 Remember Today at 5:30 in CAS-522 is the second chance for the MB
More informationMolecular Mechanics, Dynamics & Docking
Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter
More informationCS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C.
CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring 2006 Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C. Latombe Scribe: Neda Nategh How do you update the energy function during the
More informationLecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions: Van der Waals Interactions
More informationMolecular Basis of K + Conduction and Selectivity
The Structure of the Potassium Channel: Molecular Basis of K + Conduction and Selectivity -Doyle, DA, et al. The structure of the potassium channel: molecular basis of K + conduction and selectivity. Science
More informationThe protein folding problem consists of two parts:
Energetics and kinetics of protein folding The protein folding problem consists of two parts: 1)Creating a stable, well-defined structure that is significantly more stable than all other possible structures.
More information`1AP Biology Study Guide Chapter 2 v Atomic structure is the basis of life s chemistry Ø Living and non- living things are composed of atoms Ø
`1AP Biology Study Guide Chapter 2 v Atomic structure is the basis of life s chemistry Ø Living and non- living things are composed of atoms Ø Element pure substance only one kind of atom Ø Living things
More informationChapter-2 (Page 22-37) Physical and Chemical Properties of Water
Chapter-2 (Page 22-37) Physical and Chemical Properties of Water Introduction About 70% of the mass of the human body is water. Water is central to biochemistry for the following reasons: 1- Biological
More informationLimitations of temperature replica exchange (T-REMD) for protein folding simulations
Limitations of temperature replica exchange (T-REMD) for protein folding simulations Jed W. Pitera, William C. Swope IBM Research pitera@us.ibm.com Anomalies in protein folding kinetic thermodynamic 322K
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationDISCRETE TUTORIAL. Agustí Emperador. Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING:
DISCRETE TUTORIAL Agustí Emperador Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING: STRUCTURAL REFINEMENT OF DOCKING CONFORMATIONS Emperador
More informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
More informationHigh Specificity and Reversibility
Lecture #8 The Cell as a Machine High Specificity and Reversibility In considering the problem of transcription factor binding in the nucleus and the great specificity that is called for to transcribe
More informationBiomolecules are dynamic no single structure is a perfect model
Molecular Dynamics Simulations of Biomolecules References: A. R. Leach Molecular Modeling Principles and Applications Prentice Hall, 2001. M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids",
More informationThere are two types of polysaccharides in cell: glycogen and starch Starch and glycogen are polysaccharides that function to store energy Glycogen Glucose obtained from primary sources either remains soluble
More informationThe potassium channel: Structure, selectivity and diffusion
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 18 8 MAY 2000 The potassium channel: Structure, selectivity and diffusion T. W. Allen a) Protein Dynamics Unit, Department of Chemistry, Australian National
More informationQuiz 2 Morphology of Complex Materials
071003 Quiz 2 Morphology of Complex Materials 1) Explain the following terms: (for states comment on biological activity and relative size of the structure) a) Native State b) Unfolded State c) Denatured
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 informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationAn Introduction to Metabolism
An Introduction to Metabolism I. All of an organism=s chemical reactions taken together is called metabolism. A. Metabolic pathways begin with a specific molecule, which is then altered in a series of
More informationBiophysics II. Hydrophobic Bio-molecules. Key points to be covered. Molecular Interactions in Bio-molecular Structures - van der Waals Interaction
Biophysics II Key points to be covered By A/Prof. Xiang Yang Liu Biophysics & Micro/nanostructures Lab Department of Physics, NUS 1. van der Waals Interaction 2. Hydrogen bond 3. Hydrophilic vs hydrophobic
More informationMolecular Dynamics, Monte Carlo and Docking. Lecture 21. Introduction to Bioinformatics MNW2
Molecular Dynamics, Monte Carlo and Docking Lecture 21 Introduction to Bioinformatics MNW2 Allowed phi-psi angles Red areas are preferred, yellow areas are allowed, and white is avoided 2.3a Hamiltonian
More informationMolecular Mechanics. Yohann Moreau. November 26, 2015
Molecular Mechanics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Mechanics, Label RFCT 2015 November 26, 2015 1 / 29 Introduction A so-called Force-Field
More informationProtein Structure Determination from Pseudocontact Shifts Using ROSETTA
Supporting Information Protein Structure Determination from Pseudocontact Shifts Using ROSETTA Christophe Schmitz, Robert Vernon, Gottfried Otting, David Baker and Thomas Huber Table S0. Biological Magnetic
More informationProtein Structure. W. M. Grogan, Ph.D. OBJECTIVES
Protein Structure W. M. Grogan, Ph.D. OBJECTIVES 1. Describe the structure and characteristic properties of typical proteins. 2. List and describe the four levels of structure found in proteins. 3. Relate
More informationBSc and MSc Degree Examinations
Examination Candidate Number: Desk Number: BSc and MSc Degree Examinations 2018-9 Department : BIOLOGY Title of Exam: Molecular Biology and Biochemistry Part I Time Allowed: 1 hour and 30 minutes Marking
More informationAqueous solutions. Solubility of different compounds in water
Aqueous solutions Solubility of different compounds in water The dissolution of molecules into water (in any solvent actually) causes a volume change of the solution; the size of this volume change is
More informationChemical Kinetics. Goal. Objectives
4 Chemical Kinetics To understand the physical Goal factors that determine reaction rates. bjectives After this chapter, you should be able to: describe the factors that determine reaction rates. identify
More informationSupplementary Information. The Solution Structural Ensembles of RNA Kink-turn Motifs and Their Protein Complexes
Supplementary Information The Solution Structural Ensembles of RNA Kink-turn Motifs and Their Protein Complexes Xuesong Shi, a Lin Huang, b David M. J. Lilley, b Pehr B. Harbury a,c and Daniel Herschlag
More informationThe Molecular Dynamics Method
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 informationProtein Docking by Exploiting Multi-dimensional Energy Funnels
Protein Docking by Exploiting Multi-dimensional Energy Funnels Ioannis Ch. Paschalidis 1 Yang Shen 1 Pirooz Vakili 1 Sandor Vajda 2 1 Center for Information and Systems Engineering and Department of Manufacturing
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationarxiv:cond-mat/ v1 2 Feb 94
cond-mat/9402010 Properties and Origins of Protein Secondary Structure Nicholas D. Socci (1), William S. Bialek (2), and José Nelson Onuchic (1) (1) Department of Physics, University of California at San
More informationChemistry Unit: Chemical Bonding (chapter 7 and 8) Notes
Name: Period: Due Date: 1-18-2019 / 100 Formative pts. Chemistry Unit: Chemical Bonding (chapter 7 and 8) Notes Topic-1: Review: 1. Valence electrons: The electrons in the outermost of an atom Valence
More informationSubject of the Lecture:
Subject of the Lecture: Conceptual basis for the development of force fields. Implementation/validation Water - a worked example Extensions - combining molecular mechanics and quantum mechanics (QM/MM)
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 informationPatrick: An Introduction to Medicinal Chemistry 5e Chapter 04
01) Which of the following statements is not true about receptors? a. Most receptors are proteins situated inside the cell. b. Receptors contain a hollow or cleft on their surface which is known as a binding
More informationMaterial Chemistry KJM 3100/4100. Synthetic Polymers (e.g., Polystyrene, Poly(vinyl chloride), Poly(ethylene oxide))
Material Chemistry KJM 3100/4100 Lecture 1. Soft Materials: Synthetic Polymers (e.g., Polystyrene, Poly(vinyl chloride), Poly(ethylene oxide)) Biopolymers (e.g., Cellulose derivatives, Polysaccharides,
More informationKd = koff/kon = [R][L]/[RL]
Taller de docking y cribado virtual: Uso de herramientas computacionales en el diseño de fármacos Docking program GLIDE El programa de docking GLIDE Sonsoles Martín-Santamaría Shrödinger is a scientific
More informationLec.1 Chemistry Of Water
Lec.1 Chemistry Of Water Biochemistry & Medicine Biochemistry can be defined as the science concerned with the chemical basis of life. Biochemistry can be described as the science concerned with the chemical
More informationDesign of a Novel Globular Protein Fold with Atomic-Level Accuracy
Design of a Novel Globular Protein Fold with Atomic-Level Accuracy Brian Kuhlman, Gautam Dantas, Gregory C. Ireton, Gabriele Varani, Barry L. Stoddard, David Baker Presented by Kate Stafford 4 May 05 Protein
More informationThe Logic of Biological Phenomena
The Logic of Biological Phenomena molecules are lifeless yet, in appropriate complexity and number, molecules compose living things What makes living things distinct? they can grow they can move they can
More informationNanotube AFM Probe Resolution
Influence of Elastic Deformation on Single-Wall Carbon Nanotube AFM Probe Resolution Ian R. Shapiro, Santiago D. Solares, Maria J. Esplandiu, Lawrence A. Wade, William A. Goddard,* and C. Patrick Collier*
More informationDana Alsulaibi. Jaleel G.Sweis. Mamoon Ahram
15 Dana Alsulaibi Jaleel G.Sweis Mamoon Ahram Revision of last lectures: Proteins have four levels of structures. Primary,secondary, tertiary and quaternary. Primary structure is the order of amino acids
More informationRubber elasticity. Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge. February 21, 2009
Rubber elasticity Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge February 21, 2009 A rubber is a material that can undergo large deformations e.g. stretching to five
More informationOrientational degeneracy in the presence of one alignment tensor.
Orientational degeneracy in the presence of one alignment tensor. Rotation about the x, y and z axes can be performed in the aligned mode of the program to examine the four degenerate orientations of two
More informationMolecular Driving Forces
Molecular Driving Forces Statistical Thermodynamics in Chemistry and Biology SUBGfittingen 7 At 216 513 073 / / Ken A. Dill Sarina Bromberg With the assistance of Dirk Stigter on the Electrostatics chapters
More informationParticle-Based Simulation of Bio-Electronic Systems
Particle-Based Simulation of Bio-Electronic Systems Alex Smolyanitsky, and Marco Saraniti Arizona State University Outline Particle-based Brownian dynamics simulations for bioelectronic systems Complex-field
More informationCD Basis Set of Spectra that is used is that derived from comparing the spectra of globular proteins whose secondary structures are known from X-ray
CD Basis Set of Spectra that is used is that derived from comparing the spectra of globular proteins whose secondary structures are known from X-ray crystallography An example of the use of CD Modeling
More informationUNIT 1: BIOCHEMISTRY
UNIT 1: BIOCHEMISTRY UNIT 1: Biochemistry Chapter 6.1: Chemistry of Life I. Atoms, Ions, and Molecules A. Living things consist of atoms of different elements 1. An atom is the smallest basic unit of matter
More information2008 Biowerkzeug Ltd.
2008 Biowerkzeug Ltd. 1 Contents Summary...3 1 Simulation...4 1.1 Setup...4 1.2 Output...4 2 Settings...5 3 Analysis...9 3.1 Setup...9 3.2 Input options...9 3.3 Descriptions...10 Please note that we cannot
More information1) Here we review the various types of interactions that can take place between and among molecules.
Chem 431A-L02-W'05 page 1 of 6 Chem 431A-L02-W'05 Summary of lecture topics discussed in lecture 2-3: 1) Here we review the various types of interactions that can take place between and among molecules.
More information2) Atom manipulation. Xe / Ni(110) Model: Experiment:
2) Atom manipulation D. Eigler & E. Schweizer, Nature 344, 524 (1990) Xe / Ni(110) Model: Experiment: G.Meyer, et al. Applied Physics A 68, 125 (1999) First the tip is approached close to the adsorbate
More informationONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole
ONETEP PB/SA: Application to G-Quadruplex DNA Stability Danny Cole Introduction Historical overview of structure and free energy calculation of complex molecules using molecular mechanics and continuum
More informationM1.(a) (i) giant lattice allow each carbon atom is joined to three others 1
M.(a) (i) giant lattice allow each carbon atom is joined to three others atoms in graphene are covalently bonded max. 2 marks if any reference to wrong type of bonding and covalent bonds are strong or
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 informationChapter 10: Liquids, Solids, and Phase Changes
Chapter 10: Liquids, Solids, and Phase Changes In-chapter exercises: 10.1 10.6, 10.11; End-of-chapter Problems: 10.26, 10.31, 10.32, 10.33, 10.34, 10.35, 10.36, 10.39, 10.40, 10.42, 10.44, 10.45, 10.66,
More informationChapter 4. Glutamic Acid in Solution - Correlations
Chapter 4 Glutamic Acid in Solution - Correlations 4. Introduction Glutamic acid crystallises from aqueous solution, therefore the study of these molecules in an aqueous environment is necessary to understand
More informationPolar Molecules. Textbook pg Molecules in which the charge is not distributed symmetrically among the atoms making up the molecule
Textbook pg. 251-283 Polar Molecules Molecules in which the charge is not distributed symmetrically among the atoms making up the molecule Electronegativity and Polar Molecules Pauling realized that electron
More informationNature Structural and Molecular Biology: doi: /nsmb Supplementary Figure 1. Experimental approach for enhancement of unbiased Fo Fc maps.
Supplementary Figure 1 Experimental approach for enhancement of unbiased Fo Fc maps. a, c, Unbiased Fo-Fc maps of the Tth 70S post-catalysis complex at 2.55 Å resolution with (a) or without (c) bulk solvent
More informationGeneral Biology. The Energy of Life The living cell is a miniature factory where thousands of reactions occur; it converts energy in many ways
Course No: BNG2003 Credits: 3.00 General Biology 5. An Introduction into Cell Metabolism The Energy of Life The living cell is a miniature factory where thousands of reactions occur; it converts energy
More informationI690/B680 Structural Bioinformatics Spring Protein Structure Determination by NMR Spectroscopy
I690/B680 Structural Bioinformatics Spring 2006 Protein Structure Determination by NMR Spectroscopy Suggested Reading (1) Van Holde, Johnson, Ho. Principles of Physical Biochemistry, 2 nd Ed., Prentice
More information6 Shapes of molecules and intermolecular forces Answers to practice questions. OCR Chemistry A. Question Answer Marks Guidance
1 (a) (i) HI, HBr, HCl, HF 1 (a) (ii) CF 4, CH 3 I, CH 2 Br 2, CHCl 2 F 1 (b) (i) CO 2 and HCN: linear H 2 O and SCl 2 : non-linear BF 3 and SO 3 : trigonal planar NH 3 and H 3 O + : pyramidal AlCl 4 and
More informationSolutions and Non-Covalent Binding Forces
Chapter 3 Solutions and Non-Covalent Binding Forces 3.1 Solvent and solution properties Molecules stick together using the following forces: dipole-dipole, dipole-induced dipole, hydrogen bond, van der
More informationMr. Carpenter s Biology Biochemistry. Name Pd
Mr. Carpenter s Biology Biochemistry Name Pd Chapter 2 Vocabulary Atom Element Compound Molecule Ion Cohesion Adhesion Solution Acid Base Carbohydrate Monosaccharide Lipid Protein Amino acid Nucleic acid
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 informationCollision Theory. and I 2
Collision Theory To explain why chemical reactions occur, chemists have proposed a model, known as collision theory, which states that molecules must collide in order to react. These collisions can involve
More informationBiochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,
More informationChiral selection in wrapping, crossover, and braiding of DNA mediated by asymmetric bend-writhe elasticity
http://www.aimspress.com/ AIMS Biophysics, 2(4): 666-694. DOI: 10.3934/biophy.2015.4.666 Received date 28 August 2015, Accepted date 29 October 2015, Published date 06 November 2015 Research article Chiral
More information