Potential energy function calculation

Figure 1. CASSCF potential energy functions calculated along the PT reaction coordinate for the ground- (circles), and for the lowest tttt (squares) and UTT (triangles) excited singlet states of APN (Qpt=0 -form la, Qpt =1 -form Ib).

Fig. 3. Curves calculated using (8) for a series of increasing a values. The curves were calculated using tr = 0.6 nm and e = 0.4 kj/mol. Note that for a = 0.0 the normal 6-12 Lennard Jones potential energy function is recovered.

Molecular dynamics conceptually involves two phases, namely, the force calculations and the numerical integration of the equations of motion. In the first phase, force interactions among particles based on the negative gradient of the potential energy function U,... 

Spedal emphasis was placed on the calculation of spectroscopic properties and properties of distorted molecules. The potential energy function of CFF is domi-... 

How can we apply molecular dynamics simulations practically. This section gives a brief outline of a typical MD scenario. Imagine that you are interested in the response of a protein to changes in the amino add sequence, i.e., to point mutations. In this case, it is appropriate to divide the analysis into a static and a dynamic part. What we need first is a reference system, because it is advisable to base the interpretation of the calculated data on changes compared with other simulations. By taking this relative point of view, one hopes that possible errors introduced due to the assumptions and simplifications within the potential energy function may cancel out. All kinds of simulations, analyses, etc., should always be carried out for the reference and the model systems, applying the same simulation protocols. 

Molecular Dynamics and Monte Carlo Simulations. At the heart of the method of molecular dynamics is a simulation model consisting of potential energy functions, or force fields. Molecular dynamics calculations represent a deterministic method, ie, one based on the assumption that atoms move according to laws of Newtonian mechanics. Molecular dynamics simulations can be performed for short time-periods, eg, 50—100 picoseconds, to examine localized very high frequency motions, such as bond length distortions, or, over much longer periods of time, eg, 500—2000 ps, in order to derive equiUbrium properties. It is worthwhile to summarize what properties researchers can expect to evaluate by performing molecular simulations ... 

A potential energy function is a mathematical equation that allows for the potential energy, V, of a chemical system to be calculated as a function of its tliree-dimensional (3D) structure, R. The equation includes terms describing the various physical interactions that dictate the structure and properties of a chemical system. The total potential energy of a chemical system with a defined 3D strucmre, V(R)iai, can be separated into terms for the internal, V(/ )i,iBmai, and external, V(/ )extemai, potential energy as described in the following equations. 

Practical calculations always consider differences between two or more similar systems. Suppose we effect a change in the system such that the potential energy function is changed into... 

Eigure 3 represents an illustrative biological application an Asp Asn mutation, carried out either in solution or in complex with a protein [25,26]. The calculation uses a hybrid amino acid with both an Asp and an Asn side chain. Eor convenience, we divide the system into subsystems or blocks [27] Block 1 contains the ligand backbone as well as the solvent and protein (if present) block 2 is the Asp moiety of the hybrid ligand side chain block 3 is the Asn moiety. We effect the mutation by making the Asn side chain gradually appear and the Asp side chain simultaneously disappear. We choose initially the hybrid potential energy function to have the form... 

Fig. 81.—Potential energy associated with bond rotation as a function of angle, (a) Symmetrical potential according to Eq. (23) (b) and (c) potential energy functions with lowest minimum at 0=0 corresponding to the planar zigzag form of a polymethylene chain. These curves were calculated by Taylor. 0...

---