Analysis of Molecular Dynamic Simulations Using Wavelet-Based Techniques
A. Abdel-Hafiez
Experimental nuclear physics Department, Nuclear research center, AEA, Cairo, Egypt
Email address:
To cite this article:
A. Abdel-Hafiez. Analysis of Molecular Dynamic Simulations Using Wavelet-Based Techniques. American Journal of Physics and Applications. Vol. 3, No. 4, 2015, pp. 131-137. doi: 10.11648/j.ajpa.20150304.13
Abstract: In this paper, I focus on describe, calculate and analyze of molecular dynamic (MD) simulations using wavelet transform (WT) techniques by analogy with its use in signal and image processing, so that I would like to talk about the theoretical background wavelet transform methods, including what properties they have, their common types, and how to operate them. Secondly, I would introduce the continuous wavelet transform, which is especially well-suited for time course data such as molecular dynamics simulations., the WT permits filtering out the high-frequency noise without completely omitting the high-frequency phenomena whose contribution is crucial in cases where the dynamics is localized in frequency and time. Medical applications could be studied in which biomedical related research requires lots of mathematical and engineering techniques to analyze data. The WT is observed to excel in reconstructing the original signal by a subset of the basis used in the analysis and in identifying the occurrence of rare phenomena by examining the wavelet energies at high-resolution levels.
Keywords: Molecular Dynamics, Wavelet Techniques, Some Applications
1. Introduction
Great progress has been made in applying linear time-invariant techniques in signalprocessing. In such cases the deterministic part of the signal is assumedto be composed of complex exponentials, the solutions to linear timeinvariant differential equations. However, many biomedical signalsdo not meet these assumptions. Thus, the emerging techniques of time-frequency analysis can provide new insights into the natureof biological signals [1-5].
There are several different time frequency analysis methods such as short time Fourier transform, S transform, Wigner transform, Cohen’s distribution and wavelet transforms[6-9]. I would use and focus on wavelet transform since it is used in biomedical signals analysis more commonly than others what it mean it is so sensitive and accurate.
Molecular dynamics (MD) has become a common method for studying the motion of proteins over time, and it is the only available technique for examining continuous fine granularity motion at atomic resolution[10-12]. By numerically integrating Newton’s equations of motion, one can produce a series of snapshots of a protein’s trajectory through time. These snapshots, when saved at sufficiently high resolution, serve as stop-motion photography and provide a great deal of information about how proteins behave.
In this paper, we focus on the analysis of MD simulations using wavelet-based techniques. It is worth noting, however, that any molecular system that evolves over time can be analyzed with these same wavelet techniques. Nrownian dynamics simulations and elastic networks are two examples of systems whose data have a similar structure to MD systems and which could benefit from wavelet analysis as well[13-15].
2. Theory
To extract information from signals and reveal theunderlying dynamics that corresponds to the signals, proper signal processingtechnique is needed. Typically, the process of signal processing transforms atimedomain signal into another domain since the characteristic information embedded within the time domain is not readilyobservable in its original form. Mathematically, this can be achieved by representing the timedomain signal as a series of coefficients, based on a comparisonbetween the signal x(t) and template functions {Y_{n}(t)} [16,17].
The inner productbetween the two functions x(t) and Y_{n}(t) is
The inner product describes an operation of comparing thesimilarity between the signal and the template function, i.e. the degree of closeness between the two functions. This isrealized by observing the similarities between the wavelet transform and othercommonly used techniques, in terms of the choice of thetemplate functions. Anonstationarysignal is shown in Figure1 as an example. The signal consists of four groupsof impulsive signal trains. In these groups, the signals are composed of two major frequencies, 650 and 1500 Hz.
Using the notation of inner product, the Fourier transform of a signal can beexpressed as
Assuming that the signal has finite energy
Short time Fourier transformation ( STFT) can be expressed as
According to the uncertainty principle, the time and frequency resolutions of the STFT techniquecannot be chosen arbitrarily at the same time[18].
Wavelet transform can be expressed as
The definition of continuous wavelet transform
where a shifts time, b modulates the width (not frequency), and is mother wavelt.
It has superposition property. If the continuous wavelet transform of x(t) is X(s, t) and of y(t) is Y(s, t), then the continuous wavelet transform of z(t) = k1x(t) + k2y(t) can be expressed as
Suppose that the continuous wavelet transform of x(t) is , then the transform of x(t-t0) suppose that the continuous wavelet transform of x(t) is , then the continuous wavelet transform ofx(t/a) can be expressed as
This indicates that, when a signal is dilated by a, its correspondingwavelet coefficients are also dilated by a along the scale and time axes.
3. Results and Discussion
The Fourier transform, wavelets give information about the frequency domain of a signal, but, unlike the Fourier transform, which gives only average information about each frequency, wavelets give instantaneous information about how a particular frequency is localized in time. Consequently, one can obtain considerable information about the modes of a particular signal without losing information about when these modes occur or how variable they are (Fig. 2).The wavelet transform shows clearly that the signal is increasing in frequency over time while the Fourier transform shows only that low frequencies are dominant.Figure 3 shows three types of wavelet transforms which are used in analyzing the signals.
3.1. Electrocardiography (ECG)
Electrocardiography (ECG) is a transthoracic interpretation of the electrical activity of the heart over a period of time, as detected by electrodes attached to the outer surface of the skin and recorded by a device external to the body. In short, electrocardiogram is a test that records the electrical activity of the heart. ECG has been used to measure the rate and regularity of heartbeats as well as the size and position of the chambers, the presence of any damage to the heart, and the effects of drugs or devices used to regulate the heart, as shown in Figure 4 and 5.
A typical ECG tracing of the cardiac cycle consists of a P wave, a QRS wave, a T wave, and a U wave. The baseline voltage of the electrocardiogram is known as theisoelectric line. Typically the isoelectric line is measured as the portion of the tracing following the T wave and preceding the next P wave, as shown in Figure 6 and Table 1.
Feature | Description |
RP interval | The interval between an R wave and the next R wave . Normal resting heart rate is between 60 and 100 bpm |
P wave | During normal atrial depolarization, the main electrical vector is directed from the SA node towards the AV node. This turns into the P wave on the ECG. |
PR interval | The PR interval is measured from the beginning of the P wave to the beginning of the QRS wave. The PR interval reflects the time the electrical impulse takes to travel from the sinus node through the AV node and entering the ventricles. The PR interval is therefore a good estimate of AV node function. |
QRS wave | The QRS complex reflects the rapid depolarization of the right and left ventricles. They have a large muscle mass compared to the atria and so the QRS wave usually has a larger amplitude than the Pwave. |
T wave | The T wave represents the repolarization of the ventricles. The interval from the beginning of the QRS wave to the apex of the T wave is referred to as the absolute refractory period. The last half of the T wave is referred to refractory period . |
ST interval | The ST interval is from the J point to the end of the T wave. |
3.2. Edge detection in Medical Images
One of the most important things for the medical images are image compression and image denoising which could be proceed on the following two flowcharts.
This method was developed by StephaneMallat[19,20]. Its principle is based on finding local maxima of wavelet coefficients (Wavelet Transform Modulus Maxima Method), which represent the edges in the image. The method uses only horizontal and vertical coefficients values (LowHigh and High Low coefficients 1) from each level of wavelet decomposition.
Wavelet transform modulus is defined by the equation[21-23]
for each pair of horizontal and vertical coefficients
At the same position in matrices.A matrix of the gradient angles for the same pairs of coefficients is evaluated according to the equation
Points in the original image for which values of M f(u, v, 2^{j} ) are local maxima in one-dimensional neighborhood in the direction of gradient are supposed to be edge pixels. Thesepoints are distributed according the boundary of important structures in the image. To recoveredges, found individual maxima are chained.
3.3. Wavelet Reconstruction
The signal decomposition into detail and approximation coefficients (here we label them as cD and cA) these decomposedsignals contain fewer and fewer coefficients as we move to higher scale. The first scale detail coefficients, for example, are half in number as the original signal. Thisloss in resolution is a problem for localized analysis and that's where the topic ofwavelet reconstruction comes in. The process of reconstruction is used to synthesize signals from detail and approximation coefficients. We could for example retrievethe original signal from a reconstruction of the same coefficients attained by thedecomposition process (see Fig. 11).
Wavelet reconstruction can be used on multi-level transformed signals also toreconstruct the detail each level of the decomposition. Figure 12 is an example thatshows the reconstructed detail of a signal at six different scales. Remember thatthe detail represents the fluctuation in the signal from the average of the signal.With each pass of the wavelet transform (moving downward in Fig. 12) we arelooking at coarser fluctuations corresponding to larger scale regions of the signal. Inthis project we use functions that are part of the Matlab Wavelet Toolbox for bothdecomposition and reconstruction.
4. Conclusion
Wavelet analysis is a powerful tool that can be used to quickly and automatically isolate distinct motions of interest in a protein simulation. Due to their ability to locate subtle changes without being overwhelmed by larger more obvious motions, wavelets represent an ideal method for screening simulations to quickly pinpoint changes or structural rearrangements and for comparing differences in simulations, due to mutation, pH, or temperature changes, for example. Additionally, wavelets can be used to scan large databases of simulations for biochemically relevant events, such as the motion of a catalytic site or of functionally relevant loops.
This work describes how the mathematical formulas come out, what properties they have, and how to select optimal wavelets for biomedical signal processing, we can gain help from the brief categories and descriptions. Electrocardiography, Edge detection in medical images and Wavelet Reconstruction as an application has been studied and explained.
References