Transition conformation

Edman L, Mets Uand Rigler R 1996 Conformational transitions monitored for single molecules in solution Proc. Natl Acad. Sc/. USA 93 6710-15... 

Conformational Transitions of Proteins from Atomistic Simulations... 

The influence of solvent can be incorporated in an implicit fashion to yield so-called langevin modes. Although NMA has been applied to allosteric proteins previously, the predictive power of normal mode analysis is intrinsically limited to the regime of fast structural fluctuations. Slow conformational transitions are dominantly found in the regime of anharmonic protein motion. 

Molecular dynamics simulations ([McCammon and Harvey 1987]) propagate an atomistic system by iteratively solving Newton s equation of motion for each atomic particle. Due to computational constraints, simulations can only be extended to a typical time scale of 1 ns currently, and conformational transitions such as protein domains movements are unlikely to be observed. 

To facilitate conformational transitions in the before-mentioned adenylate kinase, Elamrani and co-workers scaled all atomic masses by a large factor thus allowing the use of a high effective simulation temperature of 2000K ([Elamrani et al. 1996]). To prevent protein unfolding, elements of secondary structure had to be constrained. 

Grubmiiller described a method to induce conformational transitions in proteins and derived rate constants for these ([Grubmiiller 1994]). The method employs subsequent modifications of the original potential function based on a principal component analysis of a short MD simulation. It is discussed in more detail in the chapter of Eichinger et al. in this volume. 

If both starting structure and target structure are known, the method of targeted molecular dynamics simulation can be used to enforce a conformational transition towards the given final structure during a given simulation time ([Schlitter et al. 1994]). 

By applying a pulling force at a portion of the solute molecule in a specific direction (see chapters of Eichinger et al. and Schulten in this volume), conformational transitions can be induced in specific directions. In order to reconstruct information about the underlying potential function governing protein motion, the irreversible work performed on the system by these forces must be discounted ([Balsera et al. 1997]). 

The first term represents the forces due to the electrostatic field, the second describes forces that occur at the boundary between solute and solvent regime due to the change of dielectric constant, and the third term describes ionic forces due to the tendency of the ions in solution to move into regions of lower dielectric. Applications of the so-called PBSD method on small model systems and for the interaction of a stretch of DNA with a protein model have been discussed recently ([Elcock et al. 1997]). This simulation technique guarantees equilibrated solvent at each state of the simulation and may therefore avoid some of the problems mentioned in the previous section. Due to the smaller number of particles, the method may also speed up simulations potentially. Still, to be able to simulate long time scale protein motion, the method might ideally be combined with non-equilibrium techniques to enforce conformational transitions. 

The study of slow protein dynamics is a fascinating field with still many unknowns. We have presented a number of computational techniques that are currently being used to tackle those questions. Most promising for our case seems the development of methods that combine an implicit solvent description with techniques to induce conformational transitions. 

Schlitter et al. 1994] Schlitter, J., Engels, M., Kruger, P. Targeted molecular dynamics A new approach for searching pathways of conformational transitions. J. Mol. Graph. 12 (1994) 84-89... 

As an example for an efficient yet quite accurate approximation, in the first part of our contribution we describe a combination of a structure adapted multipole method with a multiple time step scheme (FAMUSAMM — fast multistep structure adapted multipole method) and evaluate its performance. In the second part we present, as a recent application of this method, an MD study of a ligand-receptor unbinding process enforced by single molecule atomic force microscopy. Through comparison of computed unbinding forces with experimental data we evaluate the quality of the simulations. The third part sketches, as a perspective, one way to drastically extend accessible time scales if one restricts oneself to the study of conformational transitions, which arc ubiquitous in proteins and are the elementary steps of many functional conformational motions. 

Figure 8 shows a one-dimensional sketch of a small fraction of that energy landscape (bold line) including one conformational substate (minimum) as well as, to the right, one out of the typically huge number of barriers separating this local minimum from other ones. Keeping this picture in mind the conformational dynamics of a protein can be characterized as jumps between these local minima. At the MD time scale below nanoseconds only very low barriers can be overcome, so that the studied protein remains in or close to its initial conformational substate and no predictions of slower conformational transitions can be made. 

Fig. 10. Conformational flooding accelerates conformational transitions and makes them accessible for MD simulations. Top left snapshots of the protein backbone of BPTI during a 500 ps-MD simulation. Bottom left a projection of the conformational coordinates contributing most to the atomic motions shows that, on that MD time scale, the system remains in its initial configuration (CS 1). Top right Conformational flooding forces the system into new conformations after crossing high energy barriers (CS 2, CS 3,. . . ). Bottom right The projection visualizes the new conformations they remain stable, even when the applied flooding potentials (dashed contour lines) is switched off.

Mathematical Model of the Nucleic Acids Conformational Transitions with Hysteresis over Hydration-Dehydration Cycle... 

Abstract. A model of the conformational transitions of the nucleic acid molecule during the water adsorption-desorption cycle is proposed. The nucleic acid-water system is considered as an open system. The model describes the transitions between three main conformations of wet nucleic acid samples A-, B- and unordered forms. The analysis of kinetic equations shows the non-trivial bifurcation behaviour of the system which leads to the multistability. This fact allows one to explain the hysteresis phenomena observed experimentally in the nucleic acid-water system. The problem of self-organization in the nucleic acid-water system is of great importance for revealing physical mechanisms of the functioning of nucleic acids and for many specific practical fields. 

For modelling conformational transitions and nonlinear dynamics of NA a phenomenological approach is often used. This allows one not just to describe a phenomenon but also to understand the relationships between the basic physical properties of the system. There is a general algorithm for modelling in the frame of the phenomenological approach determine the dominant motions of the system in the time interval of the process treated and theti write... 

The hydration shell is formed with the increasing of the water content of the sample and the NA transforms from the unordered to A- and then to B form, in the case of DNA and DNA-like polynucleotides and salt concentrations similar to in vivo conditions. The reverse process, dehydration of NA, results in the reverse conformational transitions but they take place at the values of relative humidity (r.h.) less than the forward direction [12]. Thus, there is a conformational hysteresis over the hydration-dehydration loop. The adsorption isotherms of the NAs, i.e. the plots of the number of the adsorbed water molecules versus the r.h. of the sample at constant temperature, also demonstrate the hysteresis phenomena [13]. The hysteresis is i( producible and its value does not decrease for at least a week. 

Taking into account the hydration shell of the NA and the possibility of the water content changing we are forced to consider the water -I- nucleic acid as an open system. In the present study a phenomenological model taking into account the interdependence of hydration and the NA conformation transition processes is offered. In accordance with the algorithm described above we consider two types of the basic processes in the system and thus two time intervals the water adsorption and the conformational transitions of the NA, times of the conformational transitions being much more greater... 

We consider a finite space, which contains the NA sample and is in contact with a bath of water or water vapor. That allows one to maintain the r.h. in the experimental space at a constant level and change it when necessary. Such a scheme corresponds to the real experiments with wet NA samples. A NA molecule is simulated by a sequence of units of the same type. Thus, in the present study, we consider the case of a homogeneous NA or the case where averaging over the unit type is possible. Every unit can be found in the one of three conformational states unordered. A- or B- conformations. The units can reversibly change their conformational state. A unit corresponds to a nucleotide of a real NA. We assume that the NA strands do not diverge during conformational transitions in the wet NA samples [18]. The conformational transitions are considered as cooperative processes that are caused by the unfavorable appearance of an interface between the distinct conformations. 

The processes of the conformational transitions of NA can be illustrated with the scheme ... 

The expressions appearing in the exponents are the free energy change of the NA-water system per unit mole in the U A and A—>B conformational transitions. The terms AF p, introduced to take into account the... 

However, the B.E.T. and modificated B.E.T as well as isotherm of d Arcy and Watt fit the experimental data only in some range of the relative humidities up to about 80-85%. At the same time the adsorption in the interval 90-100% is of great interest for in this interval the A— B conformational transition, which is of biological importance, takes place [17], [18]. This disagreement can be the result of the fact that the adsorbed water molecules can form a regular lattice, structure of which depends on the conformation of the NA. To take into account this fact we assume that the water binding constants depend on the conformational variables of the model, i.e ... 

The expressions for the changes in the free energy of the hydration shell of the NA (per NA units mole) in the U A and A— B conformational transitions should be written in the form ... 

Equations (4) and (9) along with (8) and (7) form the a set of the differential-algebraic equations dependent on X which describes the behaviour of the NA water. system, namely the conformational transitions in... 

The conformational transitions in the presented model take place accord-itig to the all-or-nothing law, i.e. they occur at the certain r.h. value. The same behaviour has been observed, for example, for the helix-coil transition of the model double-stranded structure A(pA)i7-U(pU)i7 [24]. It is worth noting that this structure is homogeneous, the same is supposed in our model. 

---