Dataset Information

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations.

ABSTRACT: Molecular dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth has led to the production of a large amount of simulation data that need to be filtered for relevant information to address specific biomedical and biochemical questions. One of the most relevant molecular properties that can be investigated by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as molecular recognition, binding strength, and mechanical and structural stability. Extracting, evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. We have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to investigate biomolecular interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid analysis and data visualization without any additional programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical analysis representation is interactively combined with full mapping of the results on the molecular system through the synergistic connection between PyContact and VMD. We showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X3 receptor. As a second application, we examine the protein-protein interaction network of the mechanically ultrastable cohesin-dockering complex.

SUBMITTER: Scheurer M

PROVIDER: S-EPMC5985026 | biostudies-literature | 2018 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations.

Scheurer Maximilian M Rodenkirch Peter P Siggel Marc M Bernardi Rafael C RC Schulten Klaus K Tajkhorshid Emad E Rudack Till T

Biophysical journal 20180201 3

Molecular dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth has led to the production of a large amount of simulation data that need to be filtered for relevant information to address specific biomedical and biochemical questions. One of the most relevant molecular properties that can be investigated by all-atom MD simulation ...[more]

PMID: 29414703

Similar Datasets

Project description:Anionic lipids and cholesterols (CHOL) are critical to the function of nicotinic acetylcholine receptors (nAChR). We investigated their interactions with an open- and closed-channel alpha4beta2 nAChR by over 10 ns molecular dynamics simulations in a ternary lipid mixture of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl phosphatidic acid (POPA), and CHOL with a ratio of 3:1:1 (Haddadian et al., J. Phys. Chem. B 2008, 112, 13981). On average there were 65 and 74 interfacial lipids around the closed- and open-channel alpha4beta2 nAChR, respectively, in the equilibrated simulation systems. In the open-channel system, 42% of the interfacial POPA had acyl chains partially inserted into intra- or intersubunit cavities, as compared to only 7% in the closed-channel alpha4beta2. No CHOL was found in cavities within single subunits, though some CHOL infiltrated into the gaps between subunits. Because of its smaller headgroup, POPA could access some nonannular sites where POPC could not easily reach due to steric exclusion. Furthermore, POPA acted not only as an acceptor for hydrogen bonding (H bonding) as POPC did, but also as a donor through its hydroxyl group for H bonding with the backbone of the protein. The charged headgroup of POPA allowed the lipid to form stable salt bridges with conserved Arg and Lys residues at the interfaces of the transmembrane (TM) and extracellular (EC) or intracellular (IC) domains of the alpha4beta2. A higher number of salt bridges and hydrogen bonds (H bonds) between POPA and the alpha4beta2 nAChR were found in the open system than in the closed system, suggesting a potential role of POPA in the equilibrium between different channel states. Most interfacial POPA molecules showed lower order parameters than the bulk POPA due to the mixed effect of gauche defects, hydrophobic mismatch, and the lipid orientations near the magic angle. These unique properties enable the interfacial POPA to achieve what POPC cannot with regard to specific interactions with the protein, thereby making POPA essential for the function of nAChR.

Project description:Protein rigidity analysis is an efficient computational method for extracting flexibility information from static, X-ray crystallography protein data. Atoms and bonds are modeled as a mechanical structure and analyzed with a fast graph-based algorithm, producing a decomposition of the flexible molecule into interconnected rigid clusters. The result depends critically on noncovalent atomic interactions, primarily on how hydrogen bonds and hydrophobic interactions are computed and modeled. Ongoing research points to the stringent need for benchmarking rigidity analysis software systems, towards the goal of increasing their accuracy and validating their results, either against each other and against biologically relevant (functional) parameters. We propose two new methods for modeling hydrogen bonds and hydrophobic interactions that more accurately reflect a mechanical model, without being computationally more intensive. We evaluate them using a novel scoring method, based on the B-cubed score from the information retrieval literature, which measures how well two cluster decompositions match.To evaluate the modeling accuracy of KINARI, our pebble-game rigidity analysis system, we use a benchmark data set of 20 proteins, each with multiple distinct conformations deposited in the Protein Data Bank. Cluster decompositions for them were previously determined with the RigidFinder method from Gerstein's lab and validated against experimental data. When KINARI's default tuning parameters are used, an improvement of the B-cubed score over a crude baseline is observed in 30% of this data. With our new modeling options, improvements were observed in over 70% of the proteins in this data set. We investigate the sensitivity of the cluster decomposition score with case studies on pyruvate phosphate dikinase and calmodulin.To substantially improve the accuracy of protein rigidity analysis systems, thorough benchmarking must be performed on all current systems and future extensions. We have measured the gain in performance by comparing different modeling methods for noncovalent interactions. We showed that new criteria for modeling hydrogen bonds and hydrophobic interactions can significantly improve the results. The two new methods proposed here have been implemented and made publicly available in the current version of KINARI (v1.3), together with the benchmarking tools, which can be downloaded from our software's website, http://kinari.cs.umass.edu.

Dataset Information

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations.

Publications

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets