BilayerAnalyzer analysis: lipid_collinearity - Lipid-lipid collinearity.

Description

Estimate the lipid-lipid collinearity angles.

This analysis computes the mean lipid-lipid collinearity angle (or order parameter) using the vector represetation of the specified lipids.

This protocol is identified by the analysis key: ‘lipid_collinearity’

Initiated by instance of:

<class 'pybilt.bilayer_analyzer.analysis_protocols.LipidCollinearityProtocol'>

Syntax

lipid_collinearity analysis-ID keyword value
  • lipid_collinearity = analysis-Key - keyword/name for this analysis.

  • analysis-ID = The unique name/ID being assigned to this analysis.

  • keyword value = settings keyword value pairs

    • leaflet (str: ‘upper’, or ‘lower’): Specifies the bilayer leaflet to include in the estimate. Default: ‘upper’

    • style (str: ‘angle’, ‘order’): Specify whether to compute the tilt angle (‘angle’) or the tilt angle order parameter (‘order’). Default: ‘angle’

    • resname_2 (str): Specify the resname of the target lipid type to include in this analysis. Default: ‘first’, the first lipid in the list pulled from the com_frame representation.

    • resname_1 (str): Specify the resname of the reference lipid type to include in this analysis. Default: ‘first’, the first lipid in the list pulled from the com_frame representation.

Examples

Construct analyzer:

analyzer = BilayerAnalyzer(structure='name_of_structure_file',
                           trajectory='name_of_traj_file',
                           selection='resname POPC DOPC')

Add by string - use default settings:

analyzer.add_analysis('lipid_collinearity lipid_collinearity_1') 

Add by string - adjust a setting:

analyzer.add_analysis('lipid_collinearity lipid_collinearity_1 leaflet upper')

Add by list:

analyzer.add_analysis(list(['lipid_collinearity', 'lipid_collinearity_1', dict({'leaflet':'upper'})]))

Add by dict:

analyzer.add_analysis(dict({'analysis_key': 'lipid_collinearity', 'analysis_id': 'lipid_collinearity_1','analysis_settings':dict({'leaflet':'upper'})}))

To remove from analyzer:

analyzer.remove_analysis('lipid_collinearity_1')

Output Info:

Retrieve output after running analyses:

output = analyzer.get_analysis_data('lipid_collinearity_1')

The output is type <type 'numpy.ndarray'>

References

  1. Anton O. Chugunov, Pavel E. Volynsky, Nikolay A. Krylov, Ivan A. Boldyrev, and Roman G. Efremov, Liquid but Durable: Molecular Dynamics Simulations Explain the Unique Properties of Archaeal-Like Membranes, Scientific Reports, 4:7462, 2014, doi:10.1038/srep07462 (https://www.nature.com/articles/srep07462)