- •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
11
Computer Simulation of Biochemical
Reactions with QM–MM Methods
Paul D. Lyne
Biogen, Inc., Cambridge, Massachusetts
Owen A. Walsh
Oxford University, Oxford, England
I.INTRODUCTION
Enzymes are phenomenally proficient molecular machines that are able to catalyze chemical reactions many orders of magnitude faster than the corresponding noncatalyzed reactions in solution [1]. How enzymes achieve these large rate enhancements is a matter of continuing research [2]. In recent years our understanding of the function of enzymes has been greatly increased by advances in structural biology. Current experimental techniques such as X-ray crystallography and nuclear magnetic resonance have produced atomic resolution structures of many enzymes. From the structures of these enzymes it is possible to infer mechanisms from the orientation of amino acids at the active site of the enzyme. In some cases, using the technique of Laue diffraction [3], it is possible to study an enzymatic reaction crystallographically, although this technique is far from being widely applicable. Despite having a detailed three-dimensional view of the active sites of many enzymes, it is often not possible to unambiguously favor one reaction mechanism over another on the basis of the structure. Other experimental methods such as kinetic studies or mutagenesis experiments can augment the information available, but currently it is not possible with experimental techniques to study the entire energy profile of a reaction in the active site of an enzyme. Nor is it possible to identify and quantify on a molecular level the interactions between the substrate and the enzyme as the reaction pathway is traversed.
Computer simulation techniques offer the ability to study the potential energy surfaces of chemical reactions to a high degree of quantitative accuracy [4]. Theoretical studies of chemical reactions in the gas phase are a major field and can provide detailed insights into a variety of processes of fundamental interest in atmospheric and combustion chemistry. In the past decade theoretical methods were extended to the study of reaction processes in mesoscopic systems such as enzymatic reactions in solution, albeit to a more approximate level than the most accurate gas-phase studies.
The best computational approach to the study of chemical reactions uses quantum mechanics; however, in practice the size of the enzyme system precludes the use of tradi-
221
222 |
Lyne and Walsh |
tional quantum mechanical methods [5]. A typical simulation of a solvated enzyme could include approximately 10,000 atoms, which makes the simulation using quantum mechanics intractable even with current computing facilities. A similification of the problem can be achieved by considering only the substrate and selected active site residues while neglecting the remainder of the enzyme and solvent. However, the remainder of the enzyme and solvent, although not having a covalent effect on the reaction, can significantly alter the energetics of the system and thus need to be included in the simulation.
The emergence of hybrid quantum mechanical–molecular mechanical (QM–MM) methods in recent years addresses this problem. Pioneering studies of this type were made by Warshel and Levitt [6]. The method entails the division of the system of interest into a small region that is treated quantum mechanically, with the remainder of the system treated with computationally less expensive classical methods. The quantum region includes all the atoms that are directly involved in the chemical reaction being studied, and the remainder of the system, which is believed to change little during the reaction, is treated with a molecular mechanics force field [7]. The atoms in each system influence the other system through a coupled potential that involves electrostatic and van der Waals interactions [8–18]. Several molecular mechanics programs have been adapted to perform hybrid QM–MM simulations. In the majority of the implementations the quantum region has been treated either by empirical valence bond methods [19] or with a semiempirical method (usually AM1 [20]). These implementations have been applied, for example, to study solvation [12,21], condensed phase spectroscopy [22], conformational flexibility [23], and chemical reactivity in solution [24,25], in enzymes [18,26–36], and in DNA [37].
Although semiempirical methods have the advantage of being computationally inexpensive, they have a number of limitations [38–40]. The major limitations concern their accuracy and reliability. In general, they are less accurate than high level methods, and since they have been parametrized to reproduce the ground-state properties of molecules, they are often not well suited to the study of chemical reactions. A further disadvantage of the semiempirical methods is the limited range of elements for which parameters have been determined.
To overcome these limitations, the hybrid QM–MM potential can employ ad initio or density function methods in the quantum region. Both of these methods can ensure a higher quantitative accuracy, and the density function methods offer a computaitonally less expensive procedure for including electron correlation [5]. Several groups have reported the development of QM–MM programs that employ ab initio [8,10,13,16] or density functional methods [10,41–43].
This chapter presents the implementaiton and applicable of a QM–MM method for studying enzyme-catalyzed reactions. The application of QM–MM methods to study solution-phase reactions has been reviewed elsewhere [44]. Similiarly, empirical valence bond methods, which have been successfully applied to studying enzymatic reactions by Warshel and coworkers [19,45], are not covered in this chapter.
II. BACKGROUND
A. QM–MM Methodology
In the combined QM–MM methodology the system being studies is partitioned into a quantum mechanical region and a molecular mechanical region (Fig. 1). The quantum
Computer Simulation with QM–MM Methods |
223 |
Figure 1 Schematic diagram depicting the partitioning of an enzymatic system into quantum and classical regions. The side chains of a tyrosine and valine are treated quantum mechanically, whereas the remainder of the enzyme and added solvent are treated with a classical force field.
region will normally include the substrate, side chains of residues believed to be involved in the reaction, and any cofactors. The remainder of the protein and solvent is included in the molecular mechanics region. For the QM region, the wave function of the system, Ψ, is a Slater determinant of one-electron molecular orbitals, ψi (or Kohn–Sham orbitals in the case of density functioned theory [46]),
|Ψ(r, Rq, Rc) |ψ1α(1)ψ1β(2) ψNβ(2N) |
(1) |
where α and β refer to spin eigenfunctions, r the coordinates of the electrons, Rq the coordinates of the QM nuclei, and Rc the coordinates of the atoms in the MM region. The total energy of the system is evaluated by solving the Schro¨dinger equation with an effective Hamiltonian for the system:
ˆ |
(2) |
Heff Ψ(r, Rq, Rc) E(Rq, Rc)Ψ(r, Rq, Rc) |
For QM–MM methods it is assumed that the effective Hamiltonian can be partitioned into quantum and classical components by writing [9]
ˆ |
ˆ |
ˆ |
ˆ |
(3) |
Heff HQM HMM HQM–MM |
||||
ˆ |
is the |
pure |
ˆ |
is the classical Hamiltonian, and |
where HQM |
quantum Hamiltonian, HMM |
ˆ |
|
|
|
|
|
HQM–MM is the hybrid QM–MM Hamiltonian. From Eqs. (2) and (3) the energy of the |
|||||
system is given as |
|
|
|
|
|
|
ˆ |
ˆ |
|
|
|
E(Rq, Rc) |
Ψ|HQM |Ψ Ψ|HQM–MM |Ψ |
EMM |
(4) |
||
|
Ψ|Ψ |
|
|||
|
|
|
|
||
E(Rq, Rc) EQM EQM–MM EMM |
|
(5) |
To date the majority of QM–MM applications have employed density functional methods ab initio or semiempirical methods in the quantum region. The energy terms evaluated in these methods are generally similar, but there are specific differences. The relevant equations for the density functional based methods are described first, and this is followed by a description of the specific differences associated with the other methods.
For density functional based QM–MM methods the electronic energy terms depend explicitly on the electron density, ρ(r), of the atoms in the quantum region [46]. The
224 |
Lyne and Walsh |
electron density is determined by solving the one-electron Kohn–Sham equations [46, 47]
ˆ |
(6) |
HDFψi(r) eiψi(r), i 1, . . . , n |
where ei are the eigenfunctions associated with the Kohn–Sham orbitals. The quantum
ˆ
Hamiltonian HDF is given by
ˆ |
¯ 2 |
2 |
|
|
|
Zq |
|
qc |
|
|
|||
h |
|
|
|
|
|
|
|||||||
HDF |
|
|
i |
|
|
|
|
|
|
|
|
||
2mi |
riq |
ric |
|
|
|||||||||
i |
|
|
|
|
i,q |
i,c |
|
|
|||||
|
ZqZq′ |
∫ |
ρ(r′) |
dr′ |
∂EXC |
(7) |
|||||||
Rqq′ |
|
|r r′| |
∂ρ(r) |
||||||||||
q q′ |
|
|
|
|
|
|
|
|
|
|
|
In this equation EXC is the exchange correlation functional [46], qc is the partial charge of an atom in the classical region, Zq is the nuclear charge of an atom in the quantum region, riq is the distance between an electron and quantum atom q, ric is the distance between an electron and a classical atom c; Rqq′ is the distance between two quantum nuclei, and r′ is the coordinate of a second electron. Once the Kohn–Sham equations have been solved, the various energy terms of the DF–MM method are evaluated as
¯ 2 |
|
|
|
|
|
|
|
|
Zq |
|
|
|
|
|
|
EDF |
h |
|
∫ |
|
|
ψi(r) 2ψi (r) dr ∫ |
|
ρ(r) dr |
|
||||||
|
|
|
|rq r| |
|
|||||||||||
|
2m |
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
i |
|
|
q |
|
|
|
|
||
|
1 |
|
∫∫ |
ρ(r)ρ(r′) |
dr dr′ EXC[ρ(r)] |
Zq Zq′ |
(8) |
||||||||
2 |
|
|r r′| |
Rqq′ |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
q q′ |
|
|
|
|
and
EDF–MM ∫ |
qc |
ρ(r)dr |
Zqcq |
Vqc |
(9) |
||
|rc r| |
|
Rqc |
|
||||
c |
q,c |
q,c |
|
where Vqc is the van der Waals interaction energy between the quantum and classical regions. This Lennard-Jones term is given as
|
|
|
|
|
|
|
|
6 |
|
|
|
|
σqc |
12 |
σqc |
|
|
||||
Vqc 4εqc |
|
|
|
|
(10) |
|||||
|
|
Rqc |
|
|
Rqc |
|
|
|
This term is essential to obtain the correct geometry, because there is no Pauli repulsion between quantum and classical atoms. The molecular mechanics energy term, EMM, is calculated with the standard potential energy term from CHARMM [48], AMBER [49], or GROMOS [50], for example.
The implementation of the method using ab initio methods for the quantum region is straightforward. The analogous equations for the electronic Hamiltonian and the corresponding energies in this case are [51]
Computer Simulation with QM–MM Methods |
225 |
|
1 |
|
|
|
|
|
|
|
|
ˆ |
|
qCZα |
|
|
|
|
|
||||
|
|
|
|
|
core |
|
|
Zα Zβ |
|
|
|
|
|
|
|||||||
E |
|
|
PµV (HµV |
FµV) |
|
|
|
|
|
|
|
VαC |
(11) |
||||||||
2 |
Rαβ |
RαC |
|||||||||||||||||||
|
|
|
µ |
V |
|
|
|
α β α |
α,C |
|
α,C |
|
|||||||||
where Hµcorev is defined as |
|
α |
|
|
|
|
|
|
|
|
|
||||||||||
core |
|
|
* |
1 |
2 |
|
* |
C |
|
|
|
|
|||||||||
|
|
|
|
|
Zα |
|
|
qC |
|
|
|
|
|||||||||
Hµv |
∫ φµ |
|
2 |
|
φV dr |
∫ φµ |
|
|Rα r| |
|
|RC r| |
|
φV dr |
(12) |
The indices µ and v refer to the basis set orbitals, φ, and Pµv and Fµv are elements of the density and Fock matrices, respectively.
For QM–MM methods employing semiempirical methods in the quantum region, the implementaiton is not as straightforward. In semiempirical methods such s MNDO, AM1, or PM3 [38], a minimum basis set representation of the valence levels of each atom is employed. Thus, for example, for first row main group elements, the basis set comprises just four atomic orbitals, namely one s and three p orbitals. The remainder of the electrons in the atoms are grouped into core terms. Therefore, the added complication for QM– MM methods involving semiempirical Hamiltonians lies in how to treat the interactions between the core terms on the quantum atoms and the partial charges on atoms in the classical region. The treatment that is currently used follows that suggested by Field et al. [9], where the interaction between the semiempirical cores and the classical atoms are treated in the same manner as pure quantum interactions, with the atoms in the classical regions carrying an s orbital for this purpose. For Am1 or PM3 the interaction between the classical atoms and the cores of the quantum atoms is given as [20]
EcoreQM–MM qcZq(SqSq |ScSc)[1 exp ( αcRqc) exp ( αqRqc)]
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
qcZq |
|
c |
m |
m 2 |
q |
c |
q 2 |
|
||||
|
|
|
|
Ki |
exp [ Li |
(Rqc Mi ) |
] Ki |
exp [ Li |
(Rqc Mi ) ] |
(13) |
|||
R |
qc |
||||||||||||
|
|
i 1 |
|
|
|
|
|
|
|
From this equation it is seen that parameters have been introduced into the QM–MM method, with Kci , Lic, Mic, and αc corresponding to the pseudo s orbital on the classical atom. These parameters can be optimized to reproduce experimental or high level theoretical data. Field et al. [9] performed extensive investigations of the values of these extra parameters and suggested that the parameters Kci , Lci , and Mic (i 1, . . . , 4) can be set to zero and that αc should take a value of 5.0. These are generally the values used in most current QM–MM implementations that employ semiempirical methods in the quantum region.
Finally, the parametrization of the van der Waals part of the QM–MM interaction must be considered. This applies to all QM–MM implementations irrespective of the quantum method being employed. From Eq. (9) it can be seen that each quantum atom needs to have two Lennard-Jones parameters associated with it in order to have a van der Walls interaction with classical atoms. Generally, there are two approaches to this problem. The first is to derive a set of parameters, εq and σq, for each common atom type and then to use this standard set for any study that requires a QM–MM study. This is the most common aproach, and the ‘‘derived’’ Lennard-Jones parameters for the quantum atoms are simply the parameters found in the MM force field for the analogous atom types. For example, a study that employed a QM–MM method implemented in the program CHARMM [48] would use the appropriate Lennard-Jones parameters of the CHARMM force field [52] for the atoms in the quantum region.
226 |
Lyne and Walsh |
The second approach is to derive Lennard-Jones parameters for the quantum atoms that are specific to the problem in hand. This is a less common approach but has been shown to improve the quantitative accuracy of the QM–MM approach in specific cases [53,54]. The disadvantage of this approach, however, is that it is necessary to derive Lennard-Jones parameters for the quantum region for every different study. Since the derivation of Lennard-Jones parameters is not a trivial exercise, this method of finding van der Walls parameters for the QM–MM interaction has not been widely used.
B. The Quantum/Classical Boundary
In the preceding section the essentials of the implementation of QM–MM methods were presented. For studies on systems where there is a discrete natural boundary between the quantum and classical regions, the QM–MM methods can be applied as outlined earlier. An example of such a system might be the study of a chemical reaction in water, such as the SN2 substitution of CH3Cl by Cl . In this case the quantum region could comprise chloroform and the chloride ion, with the solvent being treated by classical methods [24]. For QM–MM studies on the active sites of enzymes there is no natural boundary between the quantum and classical regions. In these situations the quantum region might comprise the enzyme substrate and the side chains of several active-site residues. The active-site residues thus span both the quantum and classical regions, where the side chain is in the quantum region and the main chain atoms are in the classical region. The boundary between the quantum and classical regions falls across a covalent bond between the α- and β-carbons of the residue (see Fig. 2). This is an extra complication for QM–MM implementations, and it is necessary to devise a method that will handle this circumstance.
Because the electrons on atoms in the classical region are not treated explicitly, a method has to be devised that allows the electron density along the QM–MM bonds to be terminated satisfactorily. There are two common approaches that have been proposed to deal with this situation. The first of these is termed the ‘‘link atom approach.’’ In this approach a dummy atom called a link atom is introduced into the quantum system at the location of the boundary between the quantum and classical regions. The link atom serves the purpose of satisfying the truncation of the electron density in the quantum region. The link atom is not seen by the atoms in the classical region. The implementation of the link atom approach depends on the nature of the Hamiltonian used in the quantum region. In its original implementation, Field et al. [55] used the link atom in conjunction with a
Figure 2 A glutamate side chain partitioned into quantum and classical regions. The terminal CH2CO2 group is treated quantum mechanically, and the backbone atoms are treated with the molecular mechanics force field.