Wavelet Analysis of Protein Motion

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

Structural view on the role of the TRD loop in regulating DNMT3A activity: a molecular dynamics study

DNA methyltransferase 3A (DNMT3A) has been regarded as a potential epigenetic target for the development of cancer therapeutics. A number of DNMT3A inhibitors have been reported, but most of them do not have good potency, high selectivity and/or low cytotoxicity. It has been suggested that a non-conserved region around the target recognition domain (TRD) loop is implicated in the DNMT3A activity under the allosteric regulation of the ATRX-DNMT3-DNMT3L (ADD) domain, but the molecular mechanism of the regulation of the TRD loop on the DNMT3A activity needs to be elucidated. In this study, based on the reported crystal structures, the dynamics of the TRD loop in different multimerization with/without the bound guest molecule, namely the ADD domain or the DNA molecule, was investigated using conventional molecular dynamics (MD) and umbrella sampling simulations. The simulation results illustrate that the TRD loop exhibits relatively higher flexibility than the other components in the whole catalytic domain (CD), which could be well stabilized into different local minima through the binding with either the ADD domain or the DNA molecule by forming tight hydrogen-bond and salt-bridge networks involving distinct residues. Moreover, the movement of the TRD loop away from the catalytic loop upon activation could be triggered simply by the detachment of the ADD domain, but not necessarily induced by the ADD domain relocation on the CD. All these dynamic structural details could be a supplement to the previously reported crystal structure, which underlines the importance of the structural flexibility for the critical residues in the TRD loop, arousing more interest in the rational design of novel DNMT3A inhibitors targeting this region.

Self-assembly Processes in Hydrated Montmorillonite by FTIR Investigations

Experimental findings obtained by FTIR and Raman spectroscopies on montmorillonite-water mixtures at three concentration values are presented. To get some insight into the hydrogen bond network of water within the montmorillonite network, FTIR and Raman spectra have been collected as a function of time and then analyzed following two complementary approaches: An analysis of the intramolecular OH stretching mode in the spectral range of 2700–3900 cm⁻¹ in terms of two Gaussian components, and an analysis of the same OH stretching mode by wavelet cross-correlation. The FTIR and Raman investigations have been carried as a function of time for a montmorillonite-water weight composition (wt%) of 20–80%, 25–75%, and 35–65%, until the dehydrated state where the samples appear as a homogeneous rigid layer of clay. In particular, for both the FTIR and Raman spectra, the decomposition of the OH stretching band into a “closed” and an “open” contribution and the spectral wavelet analysis allow us to extract quantitative information on the time behavior of the system water content. It emerges that, the total water contribution inside the montmorillonite structure decreases as a function of time. However, the relative weight of the ordered water contribution diminishes more rapidly while the relative weight of the disordered water contribution increases, indicating that a residual water content, characterized by a highly structural disorder, rests entrapped in the montmorillonite layer structure for a longer time. From the present study, it can be inferred that the montmorillonite dehydration process promotes the layer self-assembly.

Visualizing Protein Folding and Unfolding

Protein folding/unfolding is a complicated process that defies high-resolution characterization by experimental methods. As an alternative, atomistic molecular dynamics simulations are now routinely employed to elucidate and magnify the accompanying conformational changes and the role of solvent in the folding process. However, the level of detail necessary to map the process at high spatial-temporal resolution provides an overwhelming amount of data. As more and better tools are developed for analysis of these large data sets and validation of the simulations, one is still left with the problem of visualizing the results in ways that provide insight into the folding/unfolding process. While viewing and interrogating static crystal structures has become commonplace, more and different approaches are required for dynamic, interconverting, unfolding, and refolding proteins. Here we review a variety of approaches, ranging from straightforward to complex and unintuitive for multiscale analysis and visualization of protein folding and unfolding.

Using Wavelet Analysis To Assist in Identification of Significant Events in Molecular Dynamics Simulations

Long time scale molecular dynamics (MD) simulations of biological systems are becoming increasingly commonplace due to the availability of both large-scale computational resources and significant advances in the underlying simulation methodologies. Therefore, it is useful to investigate and develop data mining and analysis techniques to quickly and efficiently extract the biologically relevant information from the incredible amount of generated data. Wavelet analysis (WA) is a technique that can quickly reveal significant motions during an MD simulation. Here, the application of WA on well-converged long time scale (tens of μs) simulations of a DNA helix is described. We show how WA combined with a simple clustering method can be used to identify both the physical and temporal locations of events with significant motion in MD trajectories. We also show that WA can not only distinguish and quantify the locations and time scales of significant motions, but by changing the maximum time scale of WA a more complete characterization of these motions can be obtained. This allows motions of different time scales to be identified or ignored as desired.

A comparison of multiscale methods for the analysis of molecular dynamics simulations

Molecular dynamics (MD) is the only technique available for obtaining dynamic protein data at atomic spatial resolution and picosecond or finer temporal resolution. In recent years, the cost of computational resources has decreased exponentially while the number of known protein structures, many of which are not characterized biochemically, has increased rapidly. These events have led to an increase in the use of MD in biological research, both to examine phenomena that cannot be resolved experimentally and to generate hypotheses that direct further experimental research. In fact, several databases of MD simulations have arisen in recent years. MD simulations, and especially MD simulation databases, contain massive amounts of data, yet interesting phenomena often occur over very short time periods and on the scale of only a few atoms. Analysis of such data must balance these fine-detail events with the global picture they create. Here, we address the multiscale nature of the problem by comparing several MD analysis methods to show their strengths and weaknesses at various scales using the wild-type and R282W mutant forms of the DNA-binding domain of protein p53. By leveraging these techniques together, we are able to pinpoint fine-detail and big picture differences between the protein's variants. Our analyses indicate that the R282W mutation of p53 destabilizes the L1 loop and loosens the H2 helix conformation, but the loosened L1 loop can be rescued by residue H115, preventing the R282W mutation from completely destabilizing the protein or abolishing activity.