Molecular modeling energy minimization

Model optimization is a further refinement of the secondary and tertiary structure. At a minimum, a molecular mechanics energy minimization is done. Often, molecular dynamics or simulated annealing are used. These are frequently chosen to search the region of conformational space relatively close to the starting structure. For marginal cases, this step is very important and larger simulations should be run. 

The following definition of computational chemistry was published in 1985 (6) quantitative modeling of chemical behavior on a computer by the formalisms of theoretical chemistry. Some quantum theoreticians naturally would like to see computational chemistry as a subset of their field (7). However, today the number of scientists employed as computational chemists well exceeds the number employed as theoreticians (8). A recent textbook author (9) views computational chemistry as encompassing not only quantum mechanics, but also molecular mechanics, [energy] minimization, simulations, conformational analysis, and other computer-based methods for understanding and predicting the behavior of molecular systems. ... 

Molecular modeling using either Monte-Carlo simulations or molecular dynamics is used to apply molecular mechanics energy minimizations to very complex systems [348]. In complex flexible molecules such as proteins or nucleic adds, the number of variable parameters, i.e., bond torsion angles, is such that the global search for energy minima becomes impossible The same problem occurs with theoretical calculations of water structure in aqueous solutions or in heavily hydrated crystals. 

The size, location, and structure of platinum clusters in H-mordenite were modeled by molecular mechanics energy minimization and molecular dynamics simulation techniques [96G1]. It was suggested that the relative stability of monoatomic platinum sites in aluminosilicate mordenites is related to the specific aluminum insertion in T sites of the framework structure. The structural features of the platinum cluster confined to the 12-ring main channel are almost independent of the total Pt content and strongly dependent upon the surrounding zeolite structural field. 

Despite the demands presented by such a calculation, a number of researchers have used ab initio models to treat the electronic and nuclear degrees of freedom for the quantum motif in molecular mechanics, energy minimization studies. Examples of this include the self-consistant reaction field methods developed by Tapia and coworkers [42-44], which represent only the quantum motif explicitly and use continuum models for the environmental effects (classical and boundary regions), and the methods implemented by Kollman and coworkers [45] in their studies of condensed phase (chemical and biochemical) reaction mechanisms. In both of these implementations the expectation value of the quantum motif Hamiltonian, defined in Eqs. (11) and (14) above, is treated at the Hartree Fock level with relatively small basis sets. 

The interpretation and ideas that the chemist gets from viewing a structure will depend on the quality of the computations that led to that representation. This point underscores the importance of understanding the concepts and limitations of the calculations that produced a particular model. For molecular mechanics energy minimizations or for molecular simulations, the quality of the force field and other parameters controlling the calculations underlies the quality of the computed results. In any approach to obtaining three-dimensional molecular structures, the quality of the model should be paramount in the mind of the user. 

Molecular modeling by minimization of the steric energy starts from an assumed structure in which the rigid bonds are replaced by flexible springs of characteristic... 

Features. GEMM is written in a host-independent manner and it has been run with an Apollo, a VAX, and a MicroVAX II as a host. GEMM can currently perform the following operations perform molecular dynamics, perform energy minimizations, compute the energy and forces for a structure, and update the nonbond list (nonbond lists are usually automatic for the other operations). In addition, a wide variety of I/O sequences are possible, such as what is needed for interactive modelling work. 

Fig. 2 Molecular structures of compounds 1 and 2. A 2D structures of 1 and 2, and B their energy-minimized 3D structures, which used the molecular modeling by Guyen et al. (2004) [11]...

On this basis the porosity and surface composition of a number of silicas and zeolites were varied systematically to maximize retention of the isothizolinone structures. For the sake of clarity, data is represented here for only four silicas (Table 1) and three zeolites (Table 2). Silicas 1 and 3 differ in their pore dimensions, these being ca. 20 A and 180 A respectively. Silicas 2 and 4, their counterparts, have been calcined to optimise the number and distribution of isolated silanol sites. Zeolites 1 and 2 are the Na- and H- forms of zeolite-Y respectively. Zeolite 3 is the H-Y zeolite after subjecting to steam calcination, thereby substantially increasing the proportion of Si Al in the structure. The minimum pore dimensions of these materials were around 15 A, selected on the basis that energy-minimized structures obtained by molecular modelling predict the widest dimension of the bulkiest biocide (OIT) to be ca. 13 A, thereby allowing entry to the pore network. 

---