Dataset Information

MDContactCom: a tool to identify differences of protein molecular dynamics from two MD simulation trajectories in terms of interresidue contacts.

ABSTRACT:

Summary

Comparing results from multiple MD simulations performed under different conditions is essential during the initial stages of analysis. We propose a tool called MD Contact Comparison (MDContactCom) that compares residue-residue contact fluctuations of two MD trajectories, quantifies the differences, identifies sites that exhibit large differences, and visualizes those sites on the protein structure. Using this method, it is possible to identify sites affected by varying simulation conditions and reveal the path of propagation of the effect even when differences between the 3D structure of the molecule and the fluctuation RMSF of each residue is unclear. MDContactCom can monitor differences in complex protein dynamics between two MD trajectories and identify candidate sites to be analyzed in more detail. As such, MDContactCom is a versatile software package for analyzing most MD simulations.

Availability

MDContactCom is freely available for download on GitLab. The software is implemented in Python3. https://gitlab.com/chiemotono/mdcontactcom.

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Motono C

PROVIDER: S-EPMC8696114 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:A lipopolysaccharide (LPS) molecule is a key component of the bacterial outer membrane used to protect the bacterium and to interact with the environment. To gain insight into its function, the study of the LPS conformation and dynamics at the molecular and cellular levels is necessary, but these highly diverse and dynamic membrane-LPS systems are difficult to study. In this work, by using NMR spectroscopy and molecular dynamics (MD) simulations, we determined the conformational preferences of an E. coli O176 O-antigen polysaccharide at the atomic level. Moreover, we analyzed the use of non-uniform sampling (NUS) for the acquisition of high dynamic range spectra, like 1H,1H-NOESY NMR experiments. A comparison of the effective transglycosidic distances derived from conventional uniformly sampled and NUS 1H,1H-NOESY data showed high similarity under equal measuring time conditions. Furthermore, the experimentally derived internuclear distances of the O-antigen polysaccharide with ten repeating units (RUs) showed very good agreement to those calculated from the MD simulations of the same O-antigen polysaccharide in solution. Analysis of the LPS bilayer simulations with five and with ten RUs revealed that, although similar with respect to populated states in solution, the O-antigen in LPS bilayers had more extended chains as a result of spatial limitations due to close packing. Additional MD simulations of O-antigen polysaccharides from E. coli O6 (branched repeating unit) and O91 (negatively charged linear repeating unit) in solution and LPS bilayers were performed and compared to those of O176 (linear polymer). For all three O-antigens, the ensemble of structures present for the polysaccharides in solution were consistent with the results from their 1H,1H-NOESY experiments. In addition, the similarities between the O-antigen on its own and as a constituent of the full LPS in bilayer environment makes it possible to realistically describe the LPS conformation and dynamics from the MD simulations.