Predicting Conformational Transitions: Conformational Flooding
!["Conformational Booding" lowers free energy barriers of CT's and thus should accelerate the transitions. The figure shows a cut through the free energy landscape F(ci) (bold line) along a particular conformational coordinate ci in the vicinity of a CS (well). To the right, a free energy barrier separates the substate from another one (not shown). Inclusion of the artificial Booding potential Vg into the Hamiltonian of the system reduces the barrier height by an amount delta ΔF (thin line).](/627560/original-1640260542.jpg?t=eyJ3aWR0aCI6MjQ2LCJvYmpfaWQiOjYyNzU2MH0%3D--f880d5cd6612c074accf8193c3a99fc012d4853a)
We present a method to predict complex structural (conformational) transitions in irregular or disordered macromolecular systems, such as proteins or glasses, at the atomic level. Our method aims at rare events, which currently cannot be predicted with traditional molecular dynamics (MD) simulations, since these currently are limited to time scales shorter than a few nanoseconds.
Given an initial conformation of the system, our method identifies one or more product states, which may be separated from the initial state by free energy barriers that are large on the scale of thermal energy. It also provides an approximate reaction path, which can be used to determine barrier heights or reaction rates with the usual techniques. The method employs an artificial potential that destabilizes the initial conformation and, thereby, lowers free energy barriers of structural transitions. As a result, transitions are accelerated and may be observed in MD simulations. An analytical estimate for the acceleration factor is given. The method is applied to two test systems, an argon microcluster and a simplified protein model. By these studies we demonstrated that our method is capable of shortening mean transition times from 0.5 μs (argon cluster) and 1.4 ns (protein model) to a few picoseconds. These results suggest that our method is particularly well suited to study biochemically relevant conformational motions in proteins at a microsecond time scale.