Easy MD Analysis (EMDA) is a Python package based on MDAnalysis created with the aim of providing an easy, yet powerful way to perform analysis of MD simulations.
Several kinds of measures have been implemented in EMDA. Each of them has its own adder function and calculator function. Adder functions are the functions through which the user is able to request a certain kind of measure, while calculators are the engines that perform the measure when the run function is executed.
The available measures are listed below:
- Distance: measures the distance between two sets of atoms
- Angle: measures the angle between three atoms
- Dihedral: measures the dihedral angle between four atoms
- Planar angle: measures the angle between the closest planes to two sets of at least three atoms
- Distance of bridging waters between two sets of atoms: identifies the closest water that is bridging between two sets of atoms and measures the distances to each
- RMSD: measures the RMSD of a set of atoms (or the whole system) in reference of a frame of the structure
- Contacts, both of a group of atoms and of a whole protein: identifies the contacts stablished by a selection in a given radius or the contacts of each residue.
Some analysis can be performed from previous measures and they are stored as Analysis classes. In this case, each analysis has its own analyser. In opposition with measures, analysis are executed when requested. All the analysers functions' names start with the 'analyse_' string.
The available analyis are listed below:
- value: analyses the value of a frame-wise measure (like distance, for instance) and returns frame-wise list containing True if the value is between the given values or False if it is not.
- contacts_frequency: analyses the contacts and returns a dictionary containing the contacts that take place and how many times it takes place (in an absolute or relative number).
- contacts_amounts: analyses the contacts and returns a frame-wise list containing how many contacts a selection (or a residue) stablishes in each frame.
- NACs (near-attack conformations): analyses two or more analysed values (so a frame-wise boolean list) and returns the combination of all the values as a boolean frame-wise list.
Some analysis or measures can be plotted. The plotters functions (named with the plot_ prefix) take the analysis or measures' result and returns a plot depending on the type of data.
The available plotters are listed below:
- values: plots a float-containing frame-wise list. A similar method has been implemented inside the Measure class.
- contacts_frequencies_diff: external plotter (so it is not a method of the EMDA class). It takes two contacts_frequency (analyser) results (or two lists of), compares them so a bar plot is returned containing the contacts that are the most different between the two sets.
pip install EasyMDA- Clone the GitHub repository in your local machine:
git clone https://github.com/MolBioMedUAB/EMDA- Move to the repository direcctory:
cd EMDA- Install using pip:
pip install .EMDA has been design to perform analysis of MD trajectories in three steps: measure, analyse, and plot. Nonetheless, measuring is the most time-consuming task in an analysis. Thus, the code has been structured in such a way that all the requested measures are firstly added using the adders and the ran with the run method. On the other hand, analysis of the measures and plotting tasks are executed at the moment.
An example Jupyter notebook showing how to perform an analysis of a sample trajectory can be found here.
Future features that will be added in the future can be found in the TO-DO file
Additions or fixed bugs in each version can be found in the CHANGELOG file.