CFTI: conformational energy/free energy calculations * Menu: * Constraints:: Note on constrained optimization implementation * CFTINT:: Description and syntax of standard conformational free energy thermodynamic integration * CFTIM:: Description and syntax of multidimensional onformational free energy thermodynamic integration
Constraints: Energy minimization with holonomic constraints has been implemented. There are no special commands for this option. Charlie Brook's TSM module allows for MD simulations with constrained values of selected conformational coordinates - distances, atoms, dihedrals. This has been expanded to also allow energy minimization using several algorithms. The method is an alternative to using harmonic restraints in generating structures of flexible molecules with desired properties, or generating adiabatic profiles. To use this option, simply enter the 'TSM' module and give set of 'FIX' commands to define set of fixed internal coordinates (see perturb.doc for details). Next specify an energy minimization (see minmiz.doc). Algorithms that work: SD, CONJ, POWE (ABNR works also, for reasons unclear to me, KK)
CFTI: standard (one-dimensional) conformational thermodynamic integration Description of method Method expands the capabilities of the TSM module. The TSM module employs the Thermodynamic Perturbation (TP) approach to conformational free energy simulations. The basis of the calculation is a MD simulation with a constrained value of a conformational coordinate. With minimal modifications, the alternative Thermodynamic Integration (TI) method is added on. In the modified code the user has the option of using TP only (as previously) or activating TI, in which case the same simulation and data files are used to give both TP and TI results. [SYNTAX CFTI] All commands are parsed by the TSM command parser, so should be within a 'TSM ... END' block. CFTI command activates the thermodynamic integration calculation the context of use should be the same as for a thermodynamic perturbation run, i.e. some coordinates should be fixed by 'FIX ,,,', saving data to a disk file should be specified by 'SAVI ...', and one perturbation should be defined by 'MOVE ...' The derivative dA/dx is calculated for the coordinate deifned in the 'MOVE ...' statement. This coordinate has to be also fixed with 'FIX ...'; other coordinates may also be fixed if desired. See test case cftigas.inp for details Notes: 1) The formatted data file generated by 'SAVI ...' may be read by both TI postprocessing command (CFTJ) and TP postprocessing (POST). The SAVI 'NWIN' keyword has meaning only for TP, it can be set to an arbitrary value if TI only is to be used. 2) For consistency with TP, the 'BY <real>' part of the 'MOVE' command was retained. The <real> value has meaning for TP only, it can be set to an arbitrary number if TI only is to be used. CFTJ [TEMP <real>] [UICP <int>] [CONT <int>] Command to calculate the conformational free energy derivative dA/dx = <dU/dx> as well as the energy-entropy components: d<U>/dx, -TdS/dx Data is read in from the formatted file generated by the 'SAVI ...' command TEMP - specifies temperature, needed for energy-entropy components UICP - specifies unit with data CONT - defines length of data block for error analysis e.g. if data file has 1000 entries, 'CONT 100' will divide data into 10 blocks and calculate the standard deviation of the mean of the block averages CFTA [FIRSt <int>] [NUNIt <int>] [BEGIn <int>] [STOP <int>] [SKIP <int>] [CONT <int>] [TEMP <real>] Command activates analysis of CFTI-generated trajectory. Trajectory coordinate file(s) should be in consecutive units FIRST, NUNI, BEGIN, STOP, SKIP - define trajectory reading CONT, TEMP - as in CFTJ Examples of usage : see test case cftigas.inp
CFTM: multidimensional conformational thermodynamic integration Description of method This is a new approach. MD simulations are performed with several conformational coordinates simultaneously constrained to fixed values. The partial derivatives of the conformational free energy with respect to all the coordinates in the fixed set are calculated from this one simulation. The free energy gradient may be used in different ways to explore conformational free energy surfaces of flexible molecules. Method expands the capabilities of the TSM module. Only TI calculations possible, no corresponding TP analysis possible. [SYNTAX CFTM] All commands are parsed by the TSM command parser, so should be within a 'TSM ... END' block. CFTM command activates the multidimensional TI method the context of use should be the same as for a thermodynamic perturbation run, i.e. several coordinates should be fixed by 'FIX ,,,', saving data to a disk file should be specified by 'SAVI ...'. and a perturbation should be defined by a 'MOVE ...' statement for each of the fixed coordinates. See test case cftmgas.inp for details Notes: 1) The formatted data file generated by 'SAVI ...' may be read by both TI postprocessing command (CFTJ) and TP postprocessing (POST). The SAVI 'NWIN' keyword has no meaning in CFTM. The CFTM data file format is different than that for CFTI. 2) For consistency with TP, the 'BY <real>' part of the 'MOVE' command was retained. The <real> value has no meaning in CFTM. 'INTE' keyword has to be specified within the 'MOVE' command. CFTC [TEMP <real>] [UICP <int>] [CONT <int>] Command to calculate the conformational free energy derivatives dA/dx_i = <dU/dx_i> as well as the energy-entropy components: d<U>/dx_i, -TdS/dx_i Data is read in from the formatted file generated by the 'SAVI ...' command TEMP - specifies temperature, needed for energy-entropy components UICP - specifies unit with data CONT - defines length of data block for error analysis e.g. if data file has 1000 entries, 'CONT 100' will divide data into 10 blocks and calculate the standard deviation of the mean of the block averages Output includes all individual partial derivatives, and optionally their analysis into groups. The derivative with respect to a path direction is also calculated. CFTB [FIRSt <int>] [NUNIt <int>] [BEGIn <int>] [STOP <int>] [SKIP <int>] [CONT <int>] [TEMP <real>] Command activates analysis of CFTM-generated trajectory. Trajectory coordinate file(s) should be in consecutive units FIRST, NUNI, BEGIN, STOP, SKIP - define trajectory reading CONT, TEMP - as in CFTJ Output is the free energy gradient with respect to the set of fixed coordinates, the derivative along a specified direction (see DIRE) and optionally a group contribution analysis. CFTS [FIRSt <int>] [NUNIt <int>] [BEGIn <int>] [STOP <int>] [SKIP <int>] [CONT <int>] [TEMP <real>] [DUNI <int>] Analogous to CFTB, additionally writes out potential energy and dU/dx_i to a disk file specified by DUNI. NCOR NUMB <int> NUMB specifies the number of internal coordinates involved (=NICP). Used in calculating the path derivative. DIRE LAMB <int> <real, real, ... , real> The LAMB value specifies number of step (progress along reaction path). The following line(s) contain NICP real numbers defining a path vector. The vector will be normalized in side the program. The unit vector will be used to calculate derivatives of dA/dl, d<U>/dl, -TdS/dl along the path from the gradients. Note: the vector components are read in free format CFTG NGRUp <int> <int, int, ..., int> <string,string,...,string> Define groups for group contribution analysis to free energy NGRUP is the number of groups. The following line(s) contain the integer group numbers of the coordinates (LGRUP(J),J=1,NICP) in free format After that follow line(s) with group symbols (i.e. tags that will be used to denote the groups) in (20A4) format (GSYM(J),J=1,NGRUP) Example of usage: The system is a decapeptide, we calculate derivatives with respect to all phi and psi backbone dihedrals (NICP=18). In the 18 'MOVE ...' commands we specify the 9 phi first and the 9 psi at the end. The following will calculate and print out an aggregate of all phi and all psi contributions labelled by the tags 'PHI' and 'PSI': cftg ngrup 2 1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2 PHI PSI cfts Note on sign of derivatives: In both CFTI and CFTM it is possible to obtain a derivative value with incorrect sign by cleverly manipulating the atom selections in the 'MOVE ...' command. A simple way of checking the sign is to run a 1-D test case using both TI and TP postprocessing (see test case cftigas.inp). A general rule is to think about how the coordinate is defined and how motions of fragments influence it. E.g. for a distance between atoms A and B, the coordinate is the length of the vector from A to B. Perturbations (TP) involve actual displacements of A and B along the =vector= from A to B; Derivative calculations (TI) do not involve actual motions of atoms, but rather predictions of how atomic positions will vary with infinitesimal coordinate changes. Moving B along this coordinate by delta > 0 will increase the coordinate, while moving A by delta will decrease the coordinate. I.e. to get correct sign of derivative you have to use the following scheme: either FIX DIST <spec atom A> <spec atom B> MOVE DIST <spec atom A> <spec atom B> BY 1.0 INTE - sele <atom B> end or FIX DIST <spec atom A> <spec atom B> MOVE DIST <spec atom A> <spec atom B> BY 1.0 INTE - sele <atom B> end sele <atom A> end see test cases cftifas.inp and cftmgas.inp for examples.
NIH/DCRT/Laboratory for Structural Biology
FDA/CBER/OVRR Biophysics Laboratory
Modified, updated and generalized by C.L. Brooks, III
The Scripps Research Institute