Re-parameterizations

Unlike the permeability runs, the results showed that, the observed data were not sufficient to distinguish all fifteen values of the porosity. The ill-conditioning of the problem was mainly due to the limited observability and it could be overcome by supplying more information such as additional data or by a re-parameterization of the reservoir model itself (rezoning the reservoir). 

The molecular mechanics method is extremely parameter dependent. A force field equation that has been empirically parameterized for calculating peptides must be used for peptides it cannot be applied to nucleic acids without being re-parameterized for that particular class of molecules. Thankfully, most small organic molecules, with molecular weights less than 800, share similar properties. Therefore, a force field that has been parameterized for one class of drug molecules can usually be transferred to another class of drug molecules. In medicinal chemistry and quantum pharmacology, a number of force fields currently enjoy widespread use. The MM2/MM3/MMX force fields are currently widely used for small molecules, while AMBER and CHARMM are used for macromolecules such as peptides and nucleic acids. 

In fact, the distinction between two-step and direct dynamics is rather fuzzy. The basic issue is what kind and amount of preliminary work is needed before starting a dynamical calculation. Direct ab initio dynamics [90,97-101] requires a minimum of preparation some tests to choose basis sets and other options may suffice. For large systems, however, fully ab initio calculations are impractical, and one has to resort to QM/MM or PCM approaches but then, a host of empirical parameters are introduced, which may need some readjustement to avoid artefacts and to improve the accuracy before starting the dynamical calculations. The same holds for the semiempirical methods in order to represent at best the excited states, one has to re-parameterize the hamiltonian. In particular, our FOMO-SCF-CI method [56-58] differs considerably from the normal SCF or SCF+CIS procedures, so that the standard parameters need to be modified. However, the parameter sets are fairly transferable, and their optimization can be limited to the atoms belonging to the chromophore. In the two-step strategies one fits the ab... 

CM 18.1 to CM 18.3 were assessed in terms of their Cooper statistics, which define an upper limit to predictive performance. In addition, cross-validated Cooper statistics, which provide a more realistic indication of a model s capacity to predict the classifications of independent data, were obtained by applying the threefold cross-validation procedure to the best-sized CTs. In the threefold cross-validation procedure, the data set is randomly divided into three approximately equal parts, the CT is re-parameterized using two thirds of the data, and predicted classifications are made for the remaining third of the data. The cross-validated Cooper statistics are the mean values of the usual Cooper statistics, taken over the three iterations of the cross-validation procedure. The Cooper statistics for CM 18.1 to CM 18.3 are summarized in Table 18.6. 

With the parameters thus obtained, the polarizabilities of six other molecules— not in the training set—containing these atoms were calculated within experimental accuracy. Van Duijnen and Swart [60] re-parameterized the same atomic polarizabilities using restricted Hartree-Fock (RHF)-optimized geometries, extended the training set to 52 molecules, the control set to 18 molecules and the set of atoms with sulfur and the halogens. Also computed molecular polarizabilities were parameterized for enabling comparison with fully quantum-chemical calculations. 

Two-way dynamic parameterisation methods. These involve a dynamic transfer of information between separate classical and quantum calculations, e.g. using successive QM calculations to dynamically re-parameterize the classical potentials the QM atoms of an MM forcefield during a live simulation. 

Next w e show how a simple re-parameterization of the rate constant can reduce the parameter correlation. Consider the mean of the temperatures at which data were collected, and re-parameterize the rate constant as in Chapter 6,. 

Although there are a wide variety of models for simulating protein folding, we have chosen a fairly simple force field and tethered-bead model to demonstrate the optimal histogram methodology. The basic model was first described by Honeycutt and Thirumalai to model the folded states of / -barrel structures [14,15] and has since been re-parameterized by us to model a-helical-type structures as well [16]. The specific sequence studied here... 

All of the aforementioned examples were developed in the context of QM/PCM calculations and would undoubtedly need to be reconsidered, or at least re-parameterized, for classical solutes. 

To improve the accuracy of implicit-solvent potential energy surfaces, non-electrostatic interactions must be included, although such interactions have received only a brief mention here. The smooth, linear-scaling PCM technology that is discussed here is immediately ready for use in MM/PBSA applications [36, 38], as a replacement for finite-difference electrostatics. Other formulas for the non-electrostatic interactions [47] can also be used in PCM calculations, possibly after some re-parameterization. In general these non-electrostatic interaction formulas depend in some way on the cavity surface area, which is smooth and easily calculable by means of the PCM algorithms discussed herein. 

There should be pointed out that the used re-parameterization is not modifying the value of the path integral but is intended to better visualizing of its properties, towards solving it. As such, from expression (3.43) now appears clearer than before that for the systems governed by smooth potentials, the series expansion may now be applied respecting the path fluctuation, here in the second order truncation ... 

Construction curves for boundaries can be extracted from existing surfaces around the boundary surface. In this case, the generation of the surface can maintain continuity with adjacent surfaces automatically. Tangent (Gl) or curvature (G2) continuity at the connection of the resultant surface with the adjacent surfaces may be specified for each boundary curve. The boundary curve can extend beyond its section used as the boundary of a surface. Construction curves are sections of boundary curves (Figure 7-40b). They can be extracted and then re-parameterized, reducing curve data and improving parameterization. 

As such, for total energy the semiempirical AMI (Austin Model 1) and PM3 (re-parameterized AMI with less repulsive nonbonding interactions) were considered among the ab initio //Fand DFTmethods. The results are presented in Tables ill and 3.28. 

Traditionally, it has been generally believed that semiempirical methods are not particularly well suited to hydrogen bonding problems. In order to overcome these limitations, the PM6 method has been re-parameterized and additional empirical terms were added to increase the accuracy of this semiemprical harrultonian for describing weak intermolecular interactions(PM6-DH) (Rezdc et al. 2009) and the second generation corrections scheme PM6-DH2 (Korth et al. 2010). 

The observation that most common united-atom force fields seem to lead to a density that is too high, led Berger et al. to re-parameterize a force field for DPPC starting with the GROMOS bonded parameters, the OPLS nonbonded parameters and three different sets of partial charges. The Lennard-Jones parameters for the hydrocarbon tails were determined from pentadecane simulations. The final choice of these authors results in a density that agrees within % with the experimental density. ... 

The MNDO/d method from Walter Thiel et al. was also introduced in 1992. This method still utilizes a modified point-charge model (multipole expansion as in AMI and PM3) for computing the two-center two-electron integrals (TERIs). MNDO/d is essentially layered onto the older (1977) MNDO model for the elements where Z < 11, Heavier elements were then re-parameterized utilizing an expanded version of... 

---