Viii Heuristic Methods
An alternative approach to the construction of a smooth reaction pathway with a well-refined transition state is the ''conjugate peak refinement'' CPR method of Fischer and Karplus 55 . As in the global methods, the path is optimized as a whole and self-consis-tently. However, all points along the path are not treated equally. The computational effort is always directed at bringing the highest energy segment of the path closer to the valley of the energy surface. Starting from some initial...
References Bxk
1. HA Scheraga. Chem Rev 71 195-217, 1971. 2. KD Gibson, HA Scheraga. In RH Sarma, MH Sarma, eds. Structure and Expression, Vol 1, From Proteins to Ribosomes. New York Adenine Press, 1988, pp 67-94. 3. N Go, T Noguti, T Nishikawa. Proc Natl Acad Sci USA 80 3696-3700, 1983. 4. RA Abagyan. In WF van Gunsteren, PK Weiner, AJ Wilkinson, eds. Computer Simulation of Biomolecular Systems. Dordrecht Kluwer, 1998, pp 363-394. 5. R Lavery. Adv Comput Biol 1 69-145, 1994. 6. N Go, HA Scheraga. J Chem Phys...
Iv Automated Interpretation Of Noe Spectra
The methods discussed in this section extend the original concept of deriving structures from experimental NMR data in two ways. First, during the structure calculation, part of the assignment problem is solved automatically. This allows specification of the NOE data in a form closer to the raw data, which makes the refinement similar to X-ray refinement. Second, the quality of the data is assessed. The methods have been recently reviewed in more detail 64,67 . A. Recognition of Incorrect...
Ii Internal And Cartesian Coordinates
The term internal coordinates'' usually refers to bond lengths, valence angles and dihedrals. They completely define relative atomic positions thus giving an alternative to the Cartesian coordinate description of molecular structures. Dihedrals corresponding to rotations around single bonds are most important because all other internal coordinates are usually considered fixed at their standard values, and the representation thus obtained is referred to as the standard geometry approximation 6 ....
Iii Normal Mode Analysis Of Large Biological Molecules
On the face of it, a normal mode analysis of many biological molecules is a daunting task. A normal mode analysis of the protein citrate synthase, a homodimer of 874 residues, would involve the diagonalization of a matrix whose order is in the tens of thousands. Out of the proteins whose structures have been solved, this is by no means outstandingly large, but a full-scale diagonalization of the Hessian of citrate synthase is still not feasible. However, methods exist that can reduce the size...
H 1
temperature factors 9 . The experimental and theoretical temperature factors are presented in Table 1. The values of the mean-square displacements are in excellent agreement. Any perturbation from ideal space-group symmetry in a crystal will give rise to diffuse scattering. The X-ray diffuse scattering intensity IDU at some point hkl in reciprocal space can be written as where Fn is the structure factor of the nth unit cell and the sum Xm runs over the relative position vectors Rm between the...
J
Note that Eq. 40d , in fact, is redundant, because the other three equations form a closed set. Nonetheless, if we solve the equations of motion for s as well, we can use the following, HNos6, as a diagnostic tool, because this quantity has to be conserved during the simulation even though HNose is no longer a Hamiltonian where g 3N in this real-variable formulation. Theoretically, static quantities are independent of the value chosen for the parameter Q. In practice, however, we observe that...
Vi Summary
A statistical mechanical formulation of implicit solvent representations provides a robust theoretical framework for understanding the influence of solvation biomolecular systems. A decomposition of the free energy in terms of nonpolar and electrostatic contributions, AW AW np AW elec , is central to many approximate treatments. An attractive and widely used treatment consists in representing the nonpolar contribution A W np by a SASA surface tension term with Eq. 15 and the electrostatic...
E Estimating Standard Errors
Computer simulation is an experimental science to the extent that calculated dynamic properties are subject to systematic and statistical errors. Sources of systematic error consist of size dependence, poor equilibration, non-bond interaction cutoff, etc. These should, of course, be estimated and eliminated where possible. It is also essential to obtain an estimate of the statistical significance of the results. Simulation averages are taken over runs of finite length, and this is the main...
References Drx
1. TL Blundell, BL Sibanda, MJE Sternberg, JM Thornton. Knowledge-based prediction of protein structures and the design of novel molecules. Nature 326 347-352, 1987. 2. J Greer. Comparative modelling methods Application to the family of the mammalian serine proteases. Proteins 7 317-334, 1990. 3. MS Johnson, N Srinivasan, R Sowdhamini, TL Blundell. Knowledge-based protein modelling. CRC Crit Rev Biochem Mol Biol 29 1-68, 1994. 4. J Bajorath, R Stenkamp, A Aruffo. Knowledge-based model building...
References Jgb
1. T Havel, K W thrich. Bull Math Biol 46 673-698, 1984. 2. R Kaptein, ERP Zuiderweg, RM Scheek, WA van Gunsteren. J Mol Biol 182 179-182, 1985. 3. W Braun, N Go. J Mol Biol 186 611-626, 1985. 4. AT Brunger, GM Clore, AM Gronenborn, M Karplus. Proc Natl Acad Sci USA 83 38013805, 1986. 5. G Lipari, A Szabo, R Levy. Nature 300 197-198, 1982. 6. ET Olejniczak, CM Dobson, M Karplus, RM Levy. J Am Chem Soc 106 1923-1930, 1984. 7. R Br schweiler, B Roux, M Blackledge, C Griesinger, M Karplus, R...
References Jbg
1. A Radzicka, R Wolfenden. Science 267 90-93, 1995. 2. AR Fersht. Structure and Mechanism in Protein Science A Guide to Enzyme Catalysis and Protein Folding. New York WH Freeman, 1999. 3. K Moffat, R Henderson. Curr Opin Struct Biol 5 656, 1995. 4. GC Schatz. J Phys Chem 100 12839, 1996. 5. M Head-Gordon. J Phys Chem 100 13213-13226, 1996. 6. A Warshel, M Levitt. J Mol Biol 103 227-249, 1976. 7. CL Brooks III, M Karplus, BM Pettitt. Proteins A Theoretical Perspective of Dynamics, Structure and...
Ii Transition State Theory
The original microscopic rate theory is the transition state theory TST 10 12 . This theory is based on two fundamental assumptions about the system dynamics. 1 There is a transition state dividing surface that separates the short-time intrastate dynamics from the long-time interstate dynamics. 2 Once the reactant gains sufficient energy in its reaction coordinate and crosses the transition state the system will lose energy and become deactivated product. That is, the reaction dynamics is...
B Optimization Procedures Used in Empirical Force Fields
Knowledge of the approaches and target data used in the optimization of an empirical force field aids in the selection of the appropriate force field for a given study and acts as the basis for extending a force field to allow for its use with new compounds see below . In this section some of the general considerations that are involved during the development of a force field are presented, followed by a more detailed description of the parameter optimization procedure. Presented in Table 1 is...
A High Temperature Molecular Dynamics
Molecular dynamics simulations, which were discussed in Chapter 3, are among the most useful methods for sampling molecular conformational space. As the simulation proceeds, the classical trajectory that is traced is in fact a subset of the molecular conformations available to the molecule at that energy for microcanonical simulations or temperature for canonical simulations . Assuming that the ergodic hypothesis holds see Chapter 3 , an infinitely long MD trajectory will cover all of...
Y
d , . d ri Fi r or, dp, d . , , d r F r dr dp where , is the Poisson bracket, H is the Hamiltonian of the system, r, and p, are the position and conjugate momentum at coordinate , rt is the time derivative of r, and Fi is the force acting on the ith degree of freedom. The state of the system, r, at time At is then given by r r At , p At U At r r 0 , p 0 49 where U t is composed of a classical time evolution operator, eiLAt, and r could be any arbitrary function that depends on all coordinates...
References 1
1. JD Watson, FHC Crick. Nature 171 737, 1953. 2. F Sanger. Annu Rev Biochem 57 1, 1988. 3. JC Kendrew, G Bodo, MH Dintzis, RG Parrish, H Wyckoff, DC Phillips. Nature 181 622, 1958. 4. MF Perutz, MG Rossmann, AF Cullis, H Muirhead, G Will, ACT North. Nature 185 416, 1960. 5. CCF Blake, DF Koenig, GA Mair, ACT North, DC Phillips, VR Sarma. Nature 206 757, 1965. 6. FC Bernstein, TF Koetzle, GJB Williams, DF Meyer Jr, MD Brice, JR Rodgers, O Kennard, T Shimanouchi, M Tasumi. J Mol Biol 112 535,...
I Introduction
The first hints of the chemical basis of life were noted approximately 150 years ago. Leading up to this initial awareness were a series of insights that living organisms comprise a hierarchy of structures organs, which are composed of individual cells, which are themselves formed of organelles of different chemical compositions, and so on. From this realization and the observation that nonviable extracts from organisms such as yeast could by themselves catalyze chemical reactions, it became...
I Introduction 1
Central to the success of any computational approach to the study of chemical systems is the quality of the mathematical model used to calculate the energy of the system as a function of its structure. For smaller chemical systems studied in the gas phase, quantum mechanical QM approaches are appropriate. The success of these methods was emphasized by the selection of John A. Pople and Walter Kohn as winners of the 1998 Nobel prize in chemistry. These methods, however, are typically limited to...