- •Foreword
- •Preface
- •Contents
- •Introduction
- •Oren M. Becker
- •Alexander D. MacKerell, Jr.
- •Masakatsu Watanabe*
- •III. SCOPE OF THE BOOK
- •IV. TOWARD A NEW ERA
- •REFERENCES
- •Atomistic Models and Force Fields
- •Alexander D. MacKerell, Jr.
- •II. POTENTIAL ENERGY FUNCTIONS
- •D. Alternatives to the Potential Energy Function
- •III. EMPIRICAL FORCE FIELDS
- •A. From Potential Energy Functions to Force Fields
- •B. Overview of Available Force Fields
- •C. Free Energy Force Fields
- •D. Applicability of Force Fields
- •IV. DEVELOPMENT OF EMPIRICAL FORCE FIELDS
- •B. Optimization Procedures Used in Empirical Force Fields
- •D. Use of Quantum Mechanical Results as Target Data
- •VI. CONCLUSION
- •REFERENCES
- •Dynamics Methods
- •Oren M. Becker
- •Masakatsu Watanabe*
- •II. TYPES OF MOTIONS
- •IV. NEWTONIAN MOLECULAR DYNAMICS
- •A. Newton’s Equation of Motion
- •C. Molecular Dynamics: Computational Algorithms
- •A. Assigning Initial Values
- •B. Selecting the Integration Time Step
- •C. Stability of Integration
- •VI. ANALYSIS OF DYNAMIC TRAJECTORIES
- •B. Averages and Fluctuations
- •C. Correlation Functions
- •D. Potential of Mean Force
- •VII. OTHER MD SIMULATION APPROACHES
- •A. Stochastic Dynamics
- •B. Brownian Dynamics
- •VIII. ADVANCED SIMULATION TECHNIQUES
- •A. Constrained Dynamics
- •C. Other Approaches and Future Direction
- •REFERENCES
- •Conformational Analysis
- •Oren M. Becker
- •II. CONFORMATION SAMPLING
- •A. High Temperature Molecular Dynamics
- •B. Monte Carlo Simulations
- •C. Genetic Algorithms
- •D. Other Search Methods
- •III. CONFORMATION OPTIMIZATION
- •A. Minimization
- •B. Simulated Annealing
- •IV. CONFORMATIONAL ANALYSIS
- •A. Similarity Measures
- •B. Cluster Analysis
- •C. Principal Component Analysis
- •REFERENCES
- •Thomas A. Darden
- •II. CONTINUUM BOUNDARY CONDITIONS
- •III. FINITE BOUNDARY CONDITIONS
- •IV. PERIODIC BOUNDARY CONDITIONS
- •REFERENCES
- •Internal Coordinate Simulation Method
- •Alexey K. Mazur
- •II. INTERNAL AND CARTESIAN COORDINATES
- •III. PRINCIPLES OF MODELING WITH INTERNAL COORDINATES
- •B. Energy Gradients
- •IV. INTERNAL COORDINATE MOLECULAR DYNAMICS
- •A. Main Problems and Historical Perspective
- •B. Dynamics of Molecular Trees
- •C. Simulation of Flexible Rings
- •A. Time Step Limitations
- •B. Standard Geometry Versus Unconstrained Simulations
- •VI. CONCLUDING REMARKS
- •REFERENCES
- •Implicit Solvent Models
- •II. BASIC FORMULATION OF IMPLICIT SOLVENT
- •A. The Potential of Mean Force
- •III. DECOMPOSITION OF THE FREE ENERGY
- •A. Nonpolar Free Energy Contribution
- •B. Electrostatic Free Energy Contribution
- •IV. CLASSICAL CONTINUUM ELECTROSTATICS
- •A. The Poisson Equation for Macroscopic Media
- •B. Electrostatic Forces and Analytic Gradients
- •C. Treatment of Ionic Strength
- •A. Statistical Mechanical Integral Equations
- •VI. SUMMARY
- •REFERENCES
- •Steven Hayward
- •II. NORMAL MODE ANALYSIS IN CARTESIAN COORDINATE SPACE
- •B. Normal Mode Analysis in Dihedral Angle Space
- •C. Approximate Methods
- •IV. NORMAL MODE REFINEMENT
- •C. Validity of the Concept of a Normal Mode Important Subspace
- •A. The Solvent Effect
- •B. Anharmonicity and Normal Mode Analysis
- •VI. CONCLUSIONS
- •ACKNOWLEDGMENT
- •REFERENCES
- •Free Energy Calculations
- •Thomas Simonson
- •II. GENERAL BACKGROUND
- •A. Thermodynamic Cycles for Solvation and Binding
- •B. Thermodynamic Perturbation Theory
- •D. Other Thermodynamic Functions
- •E. Free Energy Component Analysis
- •III. STANDARD BINDING FREE ENERGIES
- •IV. CONFORMATIONAL FREE ENERGIES
- •A. Conformational Restraints or Umbrella Sampling
- •B. Weighted Histogram Analysis Method
- •C. Conformational Constraints
- •A. Dielectric Reaction Field Approaches
- •B. Lattice Summation Methods
- •VI. IMPROVING SAMPLING
- •A. Multisubstate Approaches
- •B. Umbrella Sampling
- •C. Moving Along
- •VII. PERSPECTIVES
- •REFERENCES
- •John E. Straub
- •B. Phenomenological Rate Equations
- •II. TRANSITION STATE THEORY
- •A. Building the TST Rate Constant
- •B. Some Details
- •C. Computing the TST Rate Constant
- •III. CORRECTIONS TO TRANSITION STATE THEORY
- •A. Computing Using the Reactive Flux Method
- •B. How Dynamic Recrossings Lower the Rate Constant
- •IV. FINDING GOOD REACTION COORDINATES
- •A. Variational Methods for Computing Reaction Paths
- •B. Choice of a Differential Cost Function
- •C. Diffusional Paths
- •VI. HOW TO CONSTRUCT A REACTION PATH
- •A. The Use of Constraints and Restraints
- •B. Variationally Optimizing the Cost Function
- •VII. FOCAL METHODS FOR REFINING TRANSITION STATES
- •VIII. HEURISTIC METHODS
- •IX. SUMMARY
- •ACKNOWLEDGMENT
- •REFERENCES
- •Paul D. Lyne
- •Owen A. Walsh
- •II. BACKGROUND
- •III. APPLICATIONS
- •A. Triosephosphate Isomerase
- •B. Bovine Protein Tyrosine Phosphate
- •C. Citrate Synthase
- •IV. CONCLUSIONS
- •ACKNOWLEDGMENT
- •REFERENCES
- •Jeremy C. Smith
- •III. SCATTERING BY CRYSTALS
- •IV. NEUTRON SCATTERING
- •A. Coherent Inelastic Neutron Scattering
- •B. Incoherent Neutron Scattering
- •REFERENCES
- •Michael Nilges
- •II. EXPERIMENTAL DATA
- •A. Deriving Conformational Restraints from NMR Data
- •B. Distance Restraints
- •C. The Hybrid Energy Approach
- •III. MINIMIZATION PROCEDURES
- •A. Metric Matrix Distance Geometry
- •B. Molecular Dynamics Simulated Annealing
- •C. Folding Random Structures by Simulated Annealing
- •IV. AUTOMATED INTERPRETATION OF NOE SPECTRA
- •B. Automated Assignment of Ambiguities in the NOE Data
- •C. Iterative Explicit NOE Assignment
- •D. Symmetrical Oligomers
- •VI. INFLUENCE OF INTERNAL DYNAMICS ON THE
- •EXPERIMENTAL DATA
- •VII. STRUCTURE QUALITY AND ENERGY PARAMETERS
- •VIII. RECENT APPLICATIONS
- •REFERENCES
- •II. STEPS IN COMPARATIVE MODELING
- •C. Model Building
- •D. Loop Modeling
- •E. Side Chain Modeling
- •III. AB INITIO PROTEIN STRUCTURE MODELING METHODS
- •IV. ERRORS IN COMPARATIVE MODELS
- •VI. APPLICATIONS OF COMPARATIVE MODELING
- •VII. COMPARATIVE MODELING IN STRUCTURAL GENOMICS
- •VIII. CONCLUSION
- •ACKNOWLEDGMENTS
- •REFERENCES
- •Roland L. Dunbrack, Jr.
- •II. BAYESIAN STATISTICS
- •A. Bayesian Probability Theory
- •B. Bayesian Parameter Estimation
- •C. Frequentist Probability Theory
- •D. Bayesian Methods Are Superior to Frequentist Methods
- •F. Simulation via Markov Chain Monte Carlo Methods
- •III. APPLICATIONS IN MOLECULAR BIOLOGY
- •B. Bayesian Sequence Alignment
- •IV. APPLICATIONS IN STRUCTURAL BIOLOGY
- •A. Secondary Structure and Surface Accessibility
- •ACKNOWLEDGMENTS
- •REFERENCES
- •Computer Aided Drug Design
- •Alexander Tropsha and Weifan Zheng
- •IV. SUMMARY AND CONCLUSIONS
- •REFERENCES
- •Oren M. Becker
- •II. SIMPLE MODELS
- •III. LATTICE MODELS
- •B. Mapping Atomistic Energy Landscapes
- •C. Mapping Atomistic Free Energy Landscapes
- •VI. SUMMARY
- •REFERENCES
- •Toshiko Ichiye
- •II. ELECTRON TRANSFER PROPERTIES
- •B. Potential Energy Parameters
- •IV. REDOX POTENTIALS
- •A. Calculation of the Energy Change of the Redox Site
- •B. Calculation of the Energy Changes of the Protein
- •B. Calculation of Differences in the Energy Change of the Protein
- •VI. ELECTRON TRANSFER RATES
- •A. Theory
- •B. Application
- •REFERENCES
- •Fumio Hirata and Hirofumi Sato
- •Shigeki Kato
- •A. Continuum Model
- •B. Simulations
- •C. Reference Interaction Site Model
- •A. Molecular Polarization in Neat Water*
- •B. Autoionization of Water*
- •C. Solvatochromism*
- •F. Tautomerization in Formamide*
- •IV. SUMMARY AND PROSPECTS
- •ACKNOWLEDGMENTS
- •REFERENCES
- •Nucleic Acid Simulations
- •Alexander D. MacKerell, Jr.
- •Lennart Nilsson
- •D. DNA Phase Transitions
- •III. METHODOLOGICAL CONSIDERATIONS
- •A. Atomistic Models
- •B. Alternative Models
- •IV. PRACTICAL CONSIDERATIONS
- •A. Starting Structures
- •C. Production MD Simulation
- •D. Convergence of MD Simulations
- •WEB SITES OF INTEREST
- •REFERENCES
- •Membrane Simulations
- •Douglas J. Tobias
- •II. MOLECULAR DYNAMICS SIMULATIONS OF MEMBRANES
- •B. Force Fields
- •C. Ensembles
- •D. Time Scales
- •III. LIPID BILAYER STRUCTURE
- •A. Overall Bilayer Structure
- •C. Solvation of the Lipid Polar Groups
- •IV. MOLECULAR DYNAMICS IN MEMBRANES
- •A. Overview of Dynamic Processes in Membranes
- •B. Qualitative Picture on the 100 ps Time Scale
- •C. Incoherent Neutron Scattering Measurements of Lipid Dynamics
- •F. Hydrocarbon Chain Dynamics
- •ACKNOWLEDGMENTS
- •REFERENCES
- •Appendix: Useful Internet Resources
- •B. Molecular Modeling and Simulation Packages
- •Index
21
Membrane Simulations
Douglas J. Tobias
University of California at Irvine, Irvine, California
I.INTRODUCTION
Biological membranes hold cells together and divide them into compartments. They are also home to the lipid-soluble proteins that comprise roughly one-third of the genome and perform a variety of functions, including energy production and storage, signal transduction, and the formation of channels for the transport of substances into and out of cells and between their compartments. In 1972, Singer and Nicolson [1] assimilated available data into a model for the molecular organization of biological membranes that is now taken for granted. According to this ‘‘fluid mosaic’’ model, a lipid bilayer forms the membrane matrix in which proteins are embedded. This matrix is fluid at physiological temperatures, and both lipids and proteins are free to diffuse in the plane of the membrane.
Biological membranes contain a complex mixture of several different types of lipids. Membrane lipids are amphiphilic molecules, with both polar and nonpolar substituents. A common membrane lipid that has been the subject of numerous experimental and theoretical studies, dipalmitoylphosphatidylcholine (DPPC), is diagrammed in Figure 1. Like other membrane lipids containing two chains, DPPC is classified as an insoluble, swelling amphiphile [2]. These lipids exhibit a variety of liquid crystalline phases that may be lamellar, cubic, or hexagonal, depending on the temperature and water content. Lamellar phases consist of stacks of lipid bilayers separated by layers of water (Fig. 2). The structures and mechanical properties of lipid bilayers depend on the amount of water present. As water is added to a lamellar phase, the spacing between the bilayers increases until the ‘‘fully hydrated’’ state is reached. At this point, which in DPPC is about 28 water molecules per lipid, additional water molecules do not go between the bilayers but rather go into a bulk water phase that is distinct from the lamellar bilayer water phase.
Hydrated bilayers containing one or more lipid components are commonly employed as models for biological membranes. These model systems exhibit a multiplicity of structural ‘‘phases’’ that are not observed in biological membranes. In the state that is analogous to fluid biological membranes, the ‘‘liquid crystal’’ or Lα bilayer phase present above the main bilayer phase transition temperature, Tm, the lipid hydrocarbon chains are conformationally disordered and fluid (‘‘melted’’), and the lipids diffuse in the plane of the bilayer. At temperatures well below Tm, hydrated bilayers exist in the ‘‘gel,’’ or Lβ, state in which the mostly all-trans chains are collectively tilted and pack in a regular two-dimensional
465
466 |
Tobias |
Figure 1 Chemical structure and space-filling representation of a phosphatidylcholine, DPPC. Different parts of the molecule are referred to by the labels at the left; together the choline and phosphate are referred to as the headgroup, which is zwitterionic. In the space-filling model, H atoms are white, O and P gray, and C black. (From Ref. 55.)
array. Because the biologically relevant liquid crystalline phase is highly disordered, it is not amenable to atomic resolution structure determination. However, electron density and scattering length density profiles as a function of depth in the membrane have been mapped out by X-ray and neutron diffraction [3–6], and deuterium NMR has been used to characterize the average order in the hydrocarbon chains [7]. The dynamics of lipids and water molecules in membranes on time scales presently accessible to simulations have been characterized by various experimental techniques, including neutron scattering [8,9] and NMR [8,10].
In recent years molecular dynamics (MD) simulations have emerged as a useful tool for filling in some of the details on the molecular structure and dynamics of membranes that are not available from experiments and for providing feedback on models used to interpret experimental data. Computer hardware and simulation methodology have matured to the point where stable simulations of membranes can now be routinely performed for at least a few nanoseconds, and many of the results of membrane simulations agree remarkably well with available experimental data [11]. With the help of simulations, Singer and Nicolson’s cartoon of a membrane drawn as a flat bilayer of lollypop-like lipids is evolving into images such as the one shown in Figure 2.
This chapter describes some of the technical aspects of simulating membranes, presents results that illustrate the novel insight into membrane structure and dynamics that can be provided by simulation, and discusses the correspondence of the emerging atomic scale picture with the results of NMR and X-ray and neutron scattering experiments. We restrict our attention to pure lipid bilayers, which have largely been the focus of the field of membrane simulations to date (recent reviews can be found in Refs. 11 and 12), with an emphasis on dynamics, which has been somewhat neglected compared to structure. However, we should mention that although the machinery for simulating membranes has been developed and tested on pure lipid bilayers, the number of applications to more
Membrane Simulations |
467 |
Figure 2 Snapshot from an MD simulation of a multilamellar liquid crystalline phase DPPC bilayer. Water molecules are colored white, lipid polar groups gray, and lipid hydrocarbon chains black. The central simulation cell containing 64 DPPC and 1792 water molecules, outlined in the upper left portion of the figure, is shown along with seven replicas generated by the periodic boundary conditions. (From Ref. 55.)
complicated systems has been growing rapidly. A by no means exhaustive list includes the transport of small molecules (water [13], benzene [14], anesthetics [15,16]) and ions [17] within and across bilayers, cholesterol effects on bilayers [18], and membrane-bound peptides and proteins (surface-bound tripeptides [19,20], transmembrane polyalanine [21], an amphipathic helical peptide [22], an amphipathic segment from the human corticotrophin releasing factor [23], the lytic peptide mellitin [24], the antimicrobial peptide dermaseptin B [25], bacteriophage Pf1 coat protein [26], the transmembrane domain of ErbB-2 tyrosine kinase receptor [27], and bacteriorhodopsin [28–30]), including pores and ion channels (gramicidin [31–34], E. coli OmpF porin [35], alamethicin [36], and the influenza A M2 channel [37]).
II. MOLECULAR DYNAMICS SIMULATIONS OF MEMBRANES
A.System Size and Construction
Algorithms and computer hardware have limited the size of most of the all-atom membrane simulations performed to date to 50–100 lipids plus water. However, with the increasing availability of parallel computers, larger systems containing up to 1000 lipids are starting to be simulated [75], and systems containing more than 100 lipids will soon be routine.
468 |
Tobias |
In the typical setup, the lipids are arranged in a bilayer, with water molecules on both sides, in a central simulation cell, or ‘‘box,’’ which is then replicated by using threedimensional periodic boundary conditions to produce an infinite multilamellar system (Fig. 2). It is important to note that the size of the central cell places an upper bound on the wavelength of fluctuations that can be supported by the system.
Because there are no crystal structures of fully hydrated lipid bilayers to start from, the initial configurations for the first membrane simulations had to be constructed ‘‘by hand.’’ Many of the early simulations were based on the crystal structure of dimyristoylphosphatidylcholine dihydrate [38], in which the DMPC molecules form a bilayer with slightly tilted all-trans hydrocarbon chains. The unit cell was replicated to give the desired number of molecules, and the lipid center-of-mass positions were scaled to give the desired area per lipid. Then the system was hydrated and equilibrated to allow hydrocarbon chains to melt and the water molecules to penetrate through the headgroup regions. A less used approach that reduced the equilibration time was to pack together lipids randomly chosen from a library of configurations generated by simulating a single lipid in a mean field [39]. Woolf and Roux extended this approach by using prehydrated lipids [31]. Presently, most simulations are initiated from the end of a previous simulation, and lipid bilayer coordinates are generally provided upon request by authors of simulation papers. Incorporating solutes into bilayers can be tedious. The simplest approach, which consists of simply deleting enough lipids to create ample space, works well in some cases (e.g., for inserting cholesterol molecules or transmembrane helices), but more sophisticated approaches are more generally useful. The latter include using radial forces to create spherical or cylindrical holes in the bilayer [36], and growing groups of soft spheres at locations that create a cavity in the membrane with the desired size and shape [26].
B. Force Fields
The force fields used in atomistic membrane simulations have the same form as the analytical, empirical molecular mechanics potentials used in classical simulations of proteins and nucleic acids (see Chapter 2). These potential functions contain harmonic terms for deforming bonds and bond angles, periodic and harmonic potentials for torsions, and van der Waals and Coulomb nonbonded interactions. The atomic charges are often obtained from quantum chemical calculations on lipid fragments. Most of the other potential parameters are taken from existing force fields for proteins and nucleic acids, which are generally built up from parametrizations of model compounds. Details on biomolecular force field parametrization are given in Refs. 40 and 41, and in Chapter 2. The quality of lipid force fields may be evaluated with simulations of crystals of lipid fragments [40,42,43] and by checking their ability to reproduce well-established experimental results in simulations of hydrated lipid bilayers (see Section III.A). The level of accuracy that we have been able to achieve is typified by the densities of phospholipid and cholesterol crystals from constantpressure MD simulations with fully flexible unit cells, plotted in Figure 3 versus the experimental values.
Calculation of the energies and forces due to the long-range Coulomb interactions between charged atoms is a major problem in simulations of biological molecules (see Chapter 5). In an isolated system the number of these interactions is proportional to N 2, where N is the number of charged atoms, and the evaluation of the electrostatic interactions quickly becomes intractable as the system size is increased. Moreover, when periodic
Membrane Simulations |
469 |
Figure 3 Comparison of the densities (in g/cm3) of model compounds for membrane lipids computed from constant-pressure MD simulations with the corresponding experimental values. The model compounds include solid octane and tricosane, liquid butane, octane, tetradecane, and eicosane, and the glycerylphosphorylcholine, cyclopentylphosphorylcholine monohydrate, dilaurolyglycerol, anhydrous cholesterol, cholesterol monohydrate, and cholesterol acetate crystals. (Models from Refs. 18, 42, and 43).
boundary conditions are used, the direct sum over all the Coulomb interactions in the periodic system is conditionally convergent. There are two common solutions to these problems. The first, which was used in most biomolecular simulations until recently, is to employ a spherical truncation (‘‘cutoff’’) scheme where the electrostatic interactions
˚
are smoothly switched off at a cutoff distance (typically around 10 A). When used in conjunction with a neighbor list, the computational work is proportional to N. An alternative approach, which is essentially exact for crystals, is to include all of the electrostatic interactions by using the Ewald method (or a related technique) for summing long-range interactions in a rapidly convergent fashion in an infinite periodic system. The disadvantages of the Ewald method are that it might enhance the artificial periodicity of a finite system replicated by periodic boundary conditions and that more computational effort is necessary compared to cutoffs (the straightforward implementation scales as N 3/2). The latter deterred most biomolecular simulators from using the Ewald method until recently, when the particle mesh Ewald method (PME) [44] was introduced. In PME, the reciprocal space part of the Ewald sum is evaluated using fast Fourier transforms, resulting in an overall N log N scaling. The potential deleterious effects of the strict periodicity that results when lattice sums are used have not been assessed in membrane simulations. However, it has been explicitly demonstrated that the use of spherical truncation introduces serious artifacts into a variety of structural, thermodynamic, and transport properties of interfacial systems and lipid bilayers [45], even when large cutoff radii are used.