Package: navicat-marc
navicat-marc: modular Analysis of Representative ConformersContentsAboutInstallConceptExamplesCitationAbout↑The code runs on pure python with the following dependencies:numpyscipymatplotlibscikit-learnnetworkxInstall↑You can install marc using pip:pipinstallnavicat_marcAfterwards, you can call marc as:python-mnavicat_marc[-h][-version]-iINPUT[INPUT...][-cC][-mM][-nN][-ewinEWIN][-sfSF][-mine][-yesh][-s][-as][-efileEFILE][-vVERB][-pmPLOTMODE]or simplynavicat_marc[-h][-version]-iINPUT[INPUT...][-cC][-mM][-nN][-ewinEWIN][-sfSF][-mine][-yesh][-s][-as][-efileEFILE][-vVERB][-pmPLOTMODE]Alternatively, you can download the package and execute:pythonsetup.pyinstallAfterwards, you can call marc as:python-mnavicat_marc[-h][-version]-iINPUT[INPUT...][-cC][-mM][-nN][-ewinEWIN][-sfSF][-mine][-yesh][-s][-as][-efileEFILE][-vVERB][-pmPLOTMODE]ornavicat_marc[-h][-version]-iINPUT[INPUT...][-cC][-mM][-nN][-ewinEWIN][-sfSF][-mine][-yesh][-s][-as][-efileEFILE][-vVERB][-pmPLOTMODE]Options can be consulted using the-hflag in either case. The help menu is quite detailed.Note that the main functions are all exposed and called directly in sequential order frommarc.py, in case you want to incorporate them in your own code.Concept↑Several strategies are available for the generation of conformational ensembles. Typically, one then needs to sort the ensemble and proceed with the study of the most energetically favored conformers, which will be the most accesible thermodynamically following a Boltzmann distribution.However, sorting conformers accurately requires high quality energy computations. Accurately determining the energy of every structure may be too computationally demanding. Hence, marc provides a convenient way of accomplishing three goals:Select a handful of conformers that are representative of the diversity of the conformational ensemble using combined metrics.Apply energy cutoffs based on the available energies to remove entire clusters from the space using the-ewinflag and inputting a treshold in kcal/mol.Proceed iteratively, helping the user select non-redundant conformers than can then be refined with a higher level and fed back to marc.The default clustering metric used in marc is the"mix"distance, which measures pairwise similarity based on heavy-atom rmsd times the energy difference times the kernel of the heavy-atom dihedral angles of the system.The logic behind this choice is that rmsd ought to be good except in cases where trivial single bond rotations increase the rmsd without affecting the energy, while the dihedral metric smooths systems that only differ by a few torsions. The possible metrics (to be fed to the-mflag) are"rmsd","erel"(based on the available energies),"da"(based on the most relevant dihedral angle of the molecule),"ewrmsd"(combining geometry and energy) and"mix"(combining geometry, dihedrals and energy).Examples↑The examples subdirectory contains some examples obtained by runningCREST. Any of the xyz files can be run as:navicar_marc-i[FILENAME]Options can be consulted with the-hflag.The input of marc is either a series or xyz files or a single trajectory-like xyz file with many conformers. All structures are expected to be analogous in terms of sorting and molecular topology. Energies per conformer, at any level of theory of your liking but in atomic units, can be provided in atomic units in the title line of each xyz block or file. Alternatively, energies can be provided in a plaintext file whose filename can be passed to theewincommand line argument. Such file must contain the same number of lines as conformers and two numbers per line (separated by blank spaces): an index, and an energy in atomic units. The energy window specified in theewincommand line argument should be in kcal/mol.Note that, by default, marc will select the most representative conformer out of every cluster. If you can provide energy values that you trust strongly, themineflag will ensure that the lowest energy conformer of every cluster is selected.The output of marc arenselected xyz files which will be calledINPUT_selected_n.xyzin the runtime directory. Conformers discarded by theewinthreshold will be printed with therejectedappendix instead. The discarding checks two criteria: if a cluster has an average energy that isminekcal/mol higher than the lowest conformer (plus half a standard deviation), and its lowest energy member is also higher than the threshold, the entire cluster will be discarded.High verbosity levels (-v 1,-v 2, etc.) will print significantly more information while marc runs. To be as automated as possible, reasonable default values are set for most choices, but extreme verbosity can be obtained by raising the value.As a final note, marc does not consider hydrogen atoms for geometry analysis. You can force marc to include them by using theyeshflag.Citation↑Please cite our work with the repository DOI.