mauede at alice.it
2010-Mar-01 05:16 UTC
[R] Is package "dr" appropriate for reducing the dimensionality of molecules conformational space ?
I anticipate lacking of prior experience with dimensionality reduction problems. Some scientists concerned with drug discovery performed several steered Molecular Dynamics simulations of the alanine-dipeptide molecule dragged by a radial force from an equilibrium conformation to another different equilibrium conformation. They sampled at regular intervals 7 dihedral angles, 5 bending angles, and 4 atom pair distances all along the trajectory from the initial to the final conformation. Likewise they also calculated the work done on the molecule by the applied force. When a real system is simulated (this is a toy model) macro-molecules conformational space is huge. Not all the sampled variables change coherently during the chemical-physical reaction (conformation change) taking place. A dimensionality reduction is necessary. We need to identify the subset of variables that are necessary and sufficient to describe such a reaction. Such variables must exhibit the highest correlation with the work curves. The objective is to correlate the work curves with specific structural descriptors (variables like angles, distances, and so on) and automatically discriminate between descriptors that correlate with work curves and descriptors which are not involved in the reactive pathway. I thought I could perform a dimensionality reduction first (experimenting with some common non-linear methods) and then perform regressionto estimate the correlation (dependence) of the survived variables with the work curves. Then i came across R package "dr" which apparently does it all in one step. Unluckily the examples that come with dr documentation are far away from the research field I am involved in. I know dimensionality reduction methods apply to many different fields. I wonder whether dr can help me with finding a reduced representation of molecule conformational changes. Data samples are available upon request. Thank you in advance. Maura tutti i telefonini TIM! [[alternative HTML version deleted]]
Apparently Analagous Threads
- Generative Topographic Map
- WMTSA "wavCWTPeaks" : Error in if (time.start < times.range[1]) time.start <- times.range[1]
- load.molecules function showing error messages
- Molecule-like Notation for Arrays -- anyone interested?
- looking for Variable selections models, techniques, methods