- •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
18
Simulations of Electron
Transfer Proteins
Toshiko Ichiye
Washington State University, Pullman, Washington
I.INTRODUCTION
Electron transfer proteins serve key roles as electron carriers in a wide variety of processes in all organisms, including the major energy-transducing functions of photosynthesis and respiration, other metabolic functions such as nitrogen fixation, and biosynthesis [1]. One of the most intriguing questions about these proteins is that of how the protein itself influences the electron transfer properties of its redox site. In this chapter, computational methods used to understand the donor–acceptor energetic interactions of electron transfer proteins at a molecular level are described. The focus is on the electron transfer metalloproteins, which are the blue copper proteins, the iron-sulfur proteins, and the cytochromes. However, many of the issues discussed are important for redox-active proteins and/or metalloproteins in general.
II. ELECTRON TRANSFER PROPERTIES
The properties of electron transfer proteins that are discussed here specifically affect the electron transfer reaction and not the association or binding of the reactants. A brief overview of these properties is given here; more detailed discussions may be found elsewhere (e.g., Ref. 1). The process of electron transfer is a very simple chemical reaction, i.e., the transfer of an electron from the donor redox site to the acceptor redox site.
k |
(R1) |
D A → D A |
Thus, it is characterized by ∆G°, the overall free energy of the reaction, and k, the reaction rate.
An electron transfer reaction may be separated into two half-reactions or redox cou-
ples so that the free energy, |
∆ |
G |
° |
, can be separated into |
∆ |
° |
and |
∆ |
° |
the free energies |
|||||||
|
|
|
GA |
|
GD |
||||||||||||
of reduction of the donor (D) and the acceptor (A), respectively, by |
|
||||||||||||||||
∆ |
G |
° |
|
∆ |
° |
∆ |
° |
|
|
|
|
|
|
|
|
|
(1) |
|
|
|
GA |
|
GD |
|
|
|
|
|
|
|
|
|
393
394 Ichiye
Therefore, ∆G° is related to the redox potentials E °i of the donor and acceptor via the Nernst equation,
∆Gi nFEi |
(2) |
where n is the number of electrons transferred and F 23.06 kcal/(mol V) is Faraday’s constant. Thus, an understanding of the redox potentials of electron transfer proteins in an electron transport chain is essential, because the redox potentials determine the direction of favorable flow of the electrons. One caution: There may be differences between the redox potentials of the isolated reactants and those of the reactants in an electron transfer complex [2].
The electron transfer rates in biological systems differ from those between small transition metal complexes in solution because the electron transfer is generally long-
˚
range, often greater than 10 A [1]. For long-range transfer (the nonadiabatic limit), the rate constant is
k |
H AB2 |
|
π |
exp RT∆G‡ |
(3) |
|
λRT |
where HAB is the electronic (or tunneling) matrix element, λ is the environmental reorganization energy, and ∆G‡ is the activation energy for the reaction. HAB is a measure of the electronic coupling between the reactants and the products at the transition state, and λ is a measure of the energy required to change the environment from the equilibrium state around the reactants to what would be the equilibrium state around the products while the reactants remain in their initial state. According to Marcus theory [3], ∆G‡ is given by
∆G‡ (λ ∆G°)2/4λ |
(4) |
where ∆G° is the free energy described above. Thus, understanding λ is important, because it entails understanding how the protein keeps the activation energy ∆G‡ low enough to promote the long-range transfer.
The calculation of E°i and λ from computational methods is the focus here. Generally, the energetics of these quantities are separated into contributions from the inner and outer shells. For transfer between small molecules, the inner shell generally is defined as the entire solutes A and D, and the outer shell is generally defined as only the solvent. However, in a more practical approach for proteins, the inner shell is defined as only the redox site, which consists of the metal plus its ligands no further than atoms of the side chains that are directly coordinated to the metal, and the outer shell is defined as the rest of the protein plus the surrounding solvent. Thus
∆G °i nF(IP SHE) ∆G iin ∆G iout ∆G isite ∆G iout |
(5) |
where IP, the ionization potential, is the negative of the Franck–Condon energy required to add an electron to the species i [4]; ∆SHE, the standard hydrogen electrode correction, is 4.43 V [5]; ∆Giin, the inner shell relaxation energy, is the change in energy for species i when the metal ligand geometry (i.e., the internal coordinates) has the oxidized versus reduced values; and ∆G iout, the solvent or outer shell relaxation energy, is the change upon reduction of the solvation polarization energy [3] of species i. In the second equality, ∆G sitei is the change upon reduction in the intrinsic energy of the redox site i. Similarly,
λ λin λout |
(6) |
Simulations of Electron Transfer Proteins |
395 |
where λin, the inner shell reorganization energy, is the energy required to change the metal ligand geometry of the reactants from the initial to the final values while the charges of the reactants remain in the initial state and λout, the outer shell reorganization energy, is the energy required to change the solvent polarization from an equilibrium configuration around the reactants to one around the products, again while the charges of the reactants remain in the initial state.
The simplest approach to calculating both ∆Gi and λ is to assume that the energetics of the redox site and the outer shell can be determined independently (i.e., the energies are uncoupled). Thus, as a first approximation, the change in the energy of the redox site can be calculated from quantum mechanical calculations of analogs, and the change in the outer shell energy can be calculated from classical calculations. The coupling between the redox site and outer shell energies is through the potential energy parameters, especially the partial charges, of the redox site used in the classical calculations. Because environmental effects due to the protein and/or solvent may influence the electronic structure, a higher order approximation is to use mixed quantum–classical methods (see Chapter 11) such as are used in a study of electronic tunneling pathways in ruthenium-modified myoglobin [6]; however, such calculations are not yet routine.
III.CALCULATION TECHNIQUES FOR ELECTRON TRANSFER PROTEINS
Computer simulations of electron transfer proteins often entail a variety of calculation techniques: electronic structure calculations, molecular mechanics, and electrostatic calculations. In this section, general considerations for calculations of metalloproteins are outlined; in subsequent sections, details for studying specific redox properties are given. Quantum chemistry electronic structure calculations of the redox site are important in the calculation of the energetics of the redox site and in obtaining parameters and are discussed in Sections III.A and III.B. Both molecular mechanics and electrostatic calculations of the protein are important in understanding the outer shell energetics and are discussed in Section III.C, with a focus on molecular mechanics.
A.Quantum Chemistry of the Redox Site
Quantum mechanical electronic structure calculations of the redox site are often necessary for even classical calculations of electron transfer proteins, such as when parameters are lacking for a metal redox site. Although semiempirical methods have been used for these types of studies, here the focus is on ab initio methods. Currently, density functional theory (DFT) methods are the ab initio quality electronic structure calculations of choice for metal sites because Hartree–Fock (HF) methods are insufficient and the inclusion of configuration interactions (CI), which are important for these systems, with high level basis sets leads to intractable calculations (i.e., the computational dependence on the number of electrons n is n 3 for DFT, n 4 for MP2, and n 5 for CI techniques). In addition, DFT gives excellent results for transition metals, including bond dissociation energies, bond lengths, conformational analysis, ionization potentials, and electron affinities [7,8]. An excellent review covers the application of these methods to transition metal complexes, including redox sites within proteins and analogs of protein redox sites [9].
Whether DFT or HF methods are used, there are several issues regarding the details
396 |
Ichiye |
of the calculation that will affect various aspects of the quality of the calculation; these issues are mentioned here and are discussed more fully with respect to the specific calculation later. Comparisons throughout this chapter refer to unrestricted HF calculations of the [1Fe] or rubredoxin-type site by our group, which use a full core on the iron and geometry optimization unless otherwise specified, and to spin-unrestricted DFT calculations by Noodleman and coworkers [9], which use the local density approximation for exchange and correlation with nonlocal corrections as a perturbation or as part of the SCF procedure, frozen core orbitals, and experimental geometries unless otherwise specified.
One issue concerns whether or not the geometry of the redox site is optimized. For the iron-sulfur sites, geometry optimization has given longer FeES bond lengths than are
˚
found in experiment—both an HF calculation of a [1Fe] site (0.1–0.2 A) [10] and a DFT
˚
calculation of a [2Fe-2S] site (0.5–0.2 A) [11]. For the blue copper sites, geometry optimi-
˚
zation has given metal–ligand distances that differ by up to 0.2 A from experimental averages. However, the latter come from crystal structures of the blue copper proteins, which are at much lower resolution than analog structures and have a large range of metal–
˚
ligand bonds (0.2–0.5 A). Another approach is to use experimental geometries from high resolution crystal structures of analogs of the redox site in the proper oxidation state. Less preferably crystal structures of the protein may be used, although they are lower resolution and may exist in only a single oxidation state [12]. The use of analogs is a good approximation for iron-sulfur proteins, because the redox site structures are relatively invariant even for proteins with very different folds and/or redox potentials [13]. However, care must be taken for blue copper proteins, because the redox site structure can vary significantly in different proteins [4].
Another issue is whether the calculation is of the redox site in vacuum or with environmental effects. If the energetics of the metal site are independent of its environment, the calculation of an analog in vacuum is a good approximation. The influence of the protein on the electronic structure of the redox site apparently is rather small for Fe-S proteins, from both experimental [13] and theoretical [10] studies. However, the protein plays a much larger role in the blue copper proteins, as is seen in both experiment [4] and theory [14]. A presumably more accurate but more computationally intensive process entails adding a reaction field that represents the environment. In one approach, a dielectric continuum reaction field (possibly including fixed charges and a low dielectric region representing a protein environment) is incorporated into the self-consistent field (SCF) of the entire system to give a self-consistent reaction field. Examples include DFT calculations of the manganese superoxide dismutase (Mn-SOD) active site using oxidized geometries only in aqueous solution [12] and the [2Fe-2S] site in two different proteins [11]. Other approaches include using QM/MM methods (see Chapter 11).
B. Potential Energy Parameters
The first major obstacle in studying electron transfer and/or metalloproteins is often the lack of potential energy parameters for metal sites in proteins. Although parameters for hemes existed in some of the earliest parameter sets because of the numerous studies of myoglobin [15], hemoglobin [16], and cytochrome [17], there is a dearth of parameters for other metal sites. Parameters for iron-sulfur sites have been recently developed [18–21] based on spectroscopic data for the force constants, crystallographic data for the equilibrium values, and quantum mechanical calculations for the partial charges and for the van der Waals parameters (see Chapter 2). Parameters for other sites have also been developed [22–25].
The force constants and equilibrium values for the internal coordinates should be
Simulations of Electron Transfer Proteins |
397 |
obtained from experimental data given the computational intensity of obtaining electronic structure calculations of sufficient quality for these quantities. If possible, the equilibrium values should be obtained from analogs of the redox site, because generally the structures are at much higher resolution than those of proteins and it is possible that the protein strains the analog away from its equilibrium internal coordinates.
The partial charges for the metal sites are difficult to obtain via experiment, so quantum mechanical calculations provide the best source. Because the electrostatics of the redox site are potentially a crucial factor in calculating electron transfer properties, it is important to understand the sensitivity of the results to the values of the partial charges. The energetics of fixed structures are relatively insensitive to the values of the partial charges as long as the total charge is correct [26,27]. However, the partial charges become important for energy minimization and molecular dynamics simulations because perturbations of the protein may occur [28]. Either DFT or more traditional ab initio calculations can be used to provide wave functions, which are then used to generate single-point ‘‘partial charges.’’ Because Mulliken population analyses tend to be dependent on the basis set, the preferred method of generating the partial charges is to fit the molecular electrostatic potential determined by the nuclear interactions and the electron charge density while constraining the total charge and desired higher moments to match those determined from the electronic wave function using programs such as CHELPG [29]. However, difficulties may be encountered in obtaining good values for a coordinated metal site, because such fitting methods are relatively insensitive to the values of the partial charges of buried atoms such as the metal and can thus be more conformation-dependent than the Mulliken population analysis indicates [10]. Thus, Mulliken charges may be a better indicator of the conformational dependence of the partial charges, whereas fitted charges may be a better indicator of the magnitude of the partial charges. To obtain the best charges, fitted charges should be determined for conformations with a high degree of symmetry, which are likely to be the most accurate. Mulliken population analyses should then be compared for the high symmetry and any other desired conformation to determine if the desired conformation is likely to have the same charge distribution. The details of the electronic structure calculation may also have a significant effect on the values of the partial charges. For instance, the use of an effective core versus a full core on the metal in unrestricted HF calculations of the [1Fe] site lead to 10–15% changes in the CHELPG charge of the Fe [10]. Moreover, the use of geometry for the wrong oxidation state in the same study can lead to 5–10% changes in the CHELPG charge of the Fe. This study also indicated that overestimation of the bond lengths tends to reduce the CHELPG charge on the Fe but increase the Mulliken charge. For partial charges, the treatment of the environment becomes an issue, because the condensed phase environment tends to enhance polarization so that approaches such as using a dielectric continuum to represent either a solvent or protein environment should be considered. The effect of adding environmental effects via a dielectric continuum plus protein charges in a DFT calculation indicates less than a 1% change in charge on the Fe in the [2Fe-2S] site [11], but up to an almost 20% change in CHELPG charge on the Mn of the Mn-SOD site [12]. In determining partial charges of other organic molecules, the tendency for HF to overestimate polarization for the molecule in vacuum had been viewed as an approximate means of determining the increase in polarization as the molecule is transferred to the condensed phase [30]. Experimental verification is essential, because in studies of the [1Fe] site of rubredoxin, DFT calculations [27] and unrestricted HF calculations [10] differ by almost 50% on the Fe even though the identical fitting method was used to obtain the charges.
The van der Waals parameters of the metal in a liganded site may be significantly
398 |
Ichiye |
different from that of the free ion in solution. Experimental determination of these parameters has not been forthcoming, so quantum chemical calculations again provide the best means. In cases in which there are several ligands such as when there is octahedral coordination, the metal may be inaccessible due to the van der Waals spheres of the ligands, so it is not necessary to determine separate van der Waals parameters for the metal as has been assumed for the Fe in the heme group in the CHARMM potential [31,32]. However, for the tetrahedral coordination found in the Fe-S sites, the blue copper sites, and the zinc sites, the metal is very accessible, and thus good metal van der Waals parameters are essential. Although DFT methods are generally preferable for metal sites, they are less well developed for treating the interactions between molecules, so HF methods have been used [21,33].
C.Molecular Mechanics and Electrostatics Calculations of the Protein
Molecular mechanics and electrostatics calculations have both played an important role in studying electron transfer proteins. Molecular mechanics calculations of these proteins use the same techniques (molecular dynamics, energy minimization) as for other proteins, although special consideration must be made in simulation conditions.
The treatment of long-range electrostatics and dielectric effects is an important issue in calculations of the outer shell energetics because of the importance of the electrostatic contribution to the redox properties. Two basic approaches have been used in atomistic studies of biological molecules: (1) treating the protein as a collection of charges embedded in a low dielectric continuum and the solvent as a high dielectric continuum, possibly with an ionic strength associated with it, and solving the Poisson or Poisson–Boltzmann equation as appropriate and (2) treating both the protein and solvent atoms explicitly as a collection of partial charges, with both allowed to reorient to mimic the orientational dielectric response, in a standard molecular mechanics calculation. The molecular mechanics techniques discussed here are energy minimization and molecular dynamics in standard implementations such as in CHARMM and DISCOVER; a review of the Protein Dipoles Langevin Dipoles (PDLD) method in the program POLARIS as applied to iron-sulfur proteins was published in 1996 [34].
Poisson and Poisson–Boltzmann calculations involve calculating the electrostatics due to partial charges of atoms in the protein, generally assuming that the protein itself has a low dielectric constant εp and the surrounding solvent has a high dielectric constant εw ( 80 for water) such as in the programs DelPhi [35,36], UHBD [37,38], and MEAD. Poisson–Boltzmann calculations include nonzero salt concentration, whereas Poisson calculations are for the special case of zero salt concentration. Among the strengths of this approach are that the long-range contributions (i.e., the Born solvation energy) and salt effects are calculated accurately. One of the weaknesses is that a single εp is assumed for the entire protein, whereas dielectric relaxation is a molecular phenomenon so the dielectric response varies within a protein. It is not even clear if there is a good average value for proteins, and values ranging from εp 2 to εp 10 are used. Another potential problem arises if specific interactions of the solvent with the redox site are important, as appears to be the case for rubredoxin from calculations [19] and resonance Raman experiments (J. Sanders-Loehr, personal communication) and also for ferredoxin from resonance Raman experiments [13]. This is an important consideration in general, because the redox sites are close to the surface in many electron transfer proteins. In addition, most Poisson-type calculations do not allow relaxation of the protein or dynamic effects.