Geometry optimization Conformation search, Molecular

Geometry optimization of the proposed mimetic is included as part of the design analysis to ensure the feasibility of the desired molecular conformation. MM and semiempirical quantum mechanical methods have been used most extensively for these purposes. Conformational analysis of the proposed mimetic allows the determination of an energy profile for the molecule under consideration. This has been used by researchers to assess where the desired conformation for the mimetic resides on the molecular potential energy surface. Monte Carlo, MD, and distance geometry-based conformational search techniques have been employed extensively to sample conformational space. Computational methods that attempt to approximate the efifects of aqueous solvation on the conformational profile of the mimetic are being used more frequently as part of these efforts. 

In order for this to work, the force field must be designed to describe inter-molecular forces and vibrations away from equilibrium. If the purpose of the simulation is to search conformation space, a force field designed for geometry optimization is often used. For simulating bulk systems, it is more common to use a force field that has been designed for this purpose, such as the GROMOS or OPLS force fields. 

For a conformation in a relatively deep local minimum, a room temperature molecular dynamics simulation may not overcome the barrier and search other regions of conformational space in reasonable computing time. To overcome barriers, many conformational searches use elevated temperatures (600-1200 K) at constant energy. To search conformational space adequately, run simulations of 0.5-1.0 ps each at high temperature and save the molecular structures after each simulation. Alternatively, take a snapshot of a simulation at about one picosecond intervals to store the structure. Run a geometry optimization on each structure and compare structures to determine unique low-energy conformations. 

You can often use experimental data, such as Nuclear Overhauser Effect (NOE) signals from 2D NMR studies, as restraints. NOE signals give distances between pairs of hydrogens in a molecule. Use these distances to limit distances during a molecular mechanics geometry optimization or molecular dynamics calculation. Information on dihedral angles, deduced from NMR, can also limit a conformational search. 

---