MDN: A Web Portal for Network Analysis of Molecular Dynamics Simulations

SenseNet, a tool for analysis of protein structure networks obtained from molecular dynamics simulations

Computational methods play a key role for investigating allosteric mechanisms in proteins, with the potential of generating valuable insights for innovative drug design. Here we present the SenseNet (“Structure ENSEmble NETworks”) framework for analysis of protein structure networks, which differs from established network models by focusing on interaction timelines obtained by molecular dynamics simulations. This approach is evaluated by predicting allosteric residues reported by NMR experiments in the PDZ2 domain of hPTP1e, a reference system for which previous computational predictions have shown considerable variance. We applied two models based on the mutual information between interaction timelines to estimate the conformational influence of each residue on its local environment. In terms of accuracy our prediction model is comparable to the top performing model published for this system, but by contrast benefits from its independence from NMR structures. Our results are complementary to experimental data and the consensus of previous predictions, demonstrating the potential of our new analysis tool SenseNet. Biochemical interpretation of our model suggests that allosteric residues in the PDZ2 domain form two distinct clusters of contiguous sidechain surfaces. SenseNet is provided as a plugin for the network analysis software Cytoscape, allowing for ease of future application and contributing to a system of compatible tools bridging the fields of system and structural biology.

PSNtools for standalone and web-based structure network analyses of conformational ensembles

Structure graphs, in which interacting amino acids/nucleotides correspond to linked nodes, represent cutting-edge tools to investigate macromolecular function.

The graph-based approach defined as Protein Structure Network (PSN) was initially implemented in the Wordom software and subsequently in the webPSN server. PSNs are computed either on a molecular dynamics (MD) trajectory (PSN-MD) or on a single structure. In the latter case, information on atomic fluctuations is inferred from the Elastic Network Model-Normal Mode Analysis (ENM-NMA) (PSN-ENM). While Wordom performs both PSN-ENM and PSN-MD analyses but without output post-processing, the webPSN server performs only single-structure PSN-EMN but assisting the user in input setup and output analysis.

Here we release for the first time the standalone software PSNtools, which allows calculation and post-processing of PSN analyses carried out either on single structures or on conformational ensembles. Relevant unique and novel features of PSNtools are either comparisons of two networks or computations of consensus networks on sets of homologous/analogous macromolecular structures or conformational ensembles. Network comparisons and consensus serve to infer differences in functionally different states of the same system or network-based signatures in groups of bio-macromolecules sharing either the same functionality or the same fold.

In addition to the new software, here we release also an updated version of the webPSN server, which allows performing an interactive graphical analysis of PSN-MD, following the upload of the PSNtools output.

PSNtools, the auxiliary binary version of Wordom software, and the WebPSN server are freely available at http://webpsn.hpc.unimo.it/wpsn3.php.

MDM-TASK-web: MD-TASK and MODE-TASK web server for analyzing protein dynamics

The web server, MDM-TASK-web, combines the MD-TASK and MODE-TASK software suites, which are aimed at the coarse-grained analysis of static and all-atom MD-simulated proteins, using a variety of non-conventional approaches, such as dynamic residue network analysis, perturbation-response scanning, dynamic cross-correlation, essential dynamics and normal mode analysis. Altogether, these tools allow for the exploration of protein dynamics at various levels of detail, spanning single residue perturbations and weighted contact network representations, to global residue centrality measurements and the investigation of global protein motion. Typically, following molecular dynamic simulations designed to investigate intrinsic and extrinsic protein perturbations (for instance induced by allosteric and orthosteric ligands, protein binding, temperature, pH and mutations), this selection of tools can be used to further describe protein dynamics. This may lead to the discovery of key residues involved in biological processes, such as drug resistance. The server simplifies the set-up required for running these tools and visualizing their results. Several scripts from the tool suites were updated and new ones were also added and integrated with 2D/3D visualization via the web interface. An embedded work-flow, integrated documentation and visualization tools shorten the number of steps to follow, starting from calculations to result visualization. The Django-powered web server (available at https://mdmtaskweb.rubi.ru.ac.za/) is compatible with all major web browsers. All scripts implemented in the web platform are freely available at https://github.com/RUBi-ZA/MD-TASK/tree/mdm-task-web and https://github.com/RUBi-ZA/MODE-TASK/tree/mdm-task-web.

ANCA: A Web Server for Amino Acid Networks Construction and Analysis

Amino acid network (AAN) models empower us to gain insights into protein structures and functions by describing a protein 3D structure as a graph, where nodes represent residues and edges as amino acid interactions. Here, we present the ANCA, an interactive Web server for Amino Acids Network Construction and Analysis based on a single structure or a set of structures from the Protein Data Bank. The main purpose of ANCA is to provide a portal for three types of an environment-dependent residue contact energy (ERCE)-based network model, including amino acid contact energy network (AACEN), node-weighted amino acid contact energy network (NACEN), and edge-weighted amino acid contact energy network (EACEN). For comparison, the C-alpha distance-based network model is also included, which can be extended to protein–DNA/RNA complexes. Then, the analyses of different types of AANs were performed and compared from node, edge, and network levels. The network and corresponding structure can be visualized directly in the browser. The ANCA enables researchers to investigate diverse concerns in the framework of AAN, such as the interpretation of allosteric regulation and functional residues. The ANCA portal, together with an extensive help, is available at http://sysbio.suda.edu.cn/anca/.

webPSN v2.0: a webserver to infer fingerprints of structural communication in biomacromolecules

A mixed Protein Structure Network (PSN) and Elastic Network Model-Normal Mode Analysis (ENM-NMA)-based strategy (i.e. PSN-ENM) was developed to investigate structural communication in bio-macromolecules. Protein Structure Graphs (PSGs) are computed on a single structure, whereas information on system dynamics is supplied by ENM-NMA. The approach was implemented in a webserver (webPSN), which was significantly updated herein. The webserver now handles both proteins and nucleic acids and relies on an internal upgradable database of network parameters for ions and small molecules in all PDB structures. Apart from the radical restyle of the server and some changes in the calculation setup, other major novelties concern the possibility to: a) compute the differences in nodes, links, and communication pathways between two structures (i.e. network difference) and b) infer links, hubs, communities, and metapaths from consensus networks computed on a number of structures. These new features are useful to identify commonalties and differences between two different functional states of the same system or structural-communication signatures in homologous or analogous systems. The output analysis relies on 3D-representations, interactive tables and graphs, also available for download. Speed and accuracy make this server suitable to comparatively investigate structural communication in large sets of bio-macromolecular systems. URL: http://webpsn.hpc.unimore.it.

NAPS update: network analysis of molecular dynamics data and protein-nucleic acid complexes

Network theory is now a method of choice to gain insights in understanding protein structure, folding and function. In combination with molecular dynamics (MD) simulations, it is an invaluable tool with widespread applications such as analyzing subtle conformational changes and flexibility regions in proteins, dynamic correlation analysis across distant regions for allosteric communications, in drug design to reveal alternative binding pockets for drugs, etc. Updated version of NAPS now facilitates network analysis of the complete repertoire of these biomolecules, i.e., proteins, protein–protein/nucleic acid complexes, MD trajectories, and RNA. Various options provided for analysis of MD trajectories include individual network construction and analysis of intermediate time-steps, comparative analysis of these networks, construction and analysis of average network of the ensemble of trajectories and dynamic cross-correlations. For protein–nucleic acid complexes, networks of the whole complex as well as that of the interface can be constructed and analyzed. For analysis of proteins, protein–protein complexes and MD trajectories, network construction based on inter-residue interaction energies with realistic edge-weights obtained from standard force fields is provided to capture the atomistic details. Updated version of NAPS also provides improved visualization features, interactive plots and bulk execution. URL: http://bioinf.iiit.ac.in/NAPS/

gRINN: a tool for calculation of residue interaction energies and protein energy network analysis of molecular dynamics simulations

Atomistic molecular dynamics (MD) simulations generate a wealth of information related to the dynamics of proteins. If properly analyzed, this information can lead to new insights regarding protein function and assist wet-lab experiments. Aiming to identify interactions between individual amino acid residues and the role played by each in the context of MD simulations, we present a stand-alone software called gRINN (get Residue Interaction eNergies and Networks). gRINN features graphical user interfaces (GUIs) and a command-line interface for generating and analyzing pairwise residue interaction energies and energy correlations from protein MD simulation trajectories. gRINN utilizes the features of NAMD or GROMACS MD simulation packages and automatizes the steps necessary to extract residue-residue interaction energies from user-supplied simulation trajectories, greatly simplifying the analysis for the end-user. A GUI, including an embedded molecular viewer, is provided for visualization of interaction energy time-series, distributions, an interaction energy matrix, interaction energy correlations and a residue correlation matrix. gRINN additionally offers construction and analysis of Protein Energy Networks, providing residue-based metrics such as degrees, betweenness-centralities, closeness centralities as well as shortest path analysis. gRINN is free and open to all users without login requirement at http://grinn.readthedocs.io.

Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications

Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.

An optimal distance cutoff for contact-based Protein Structure Networks using side-chain centers of mass

Proteins are highly dynamic entities attaining a myriad of different conformations. Protein side chains change their states during dynamics, causing clashes that are propagated at distal sites. A convenient formalism to analyze protein dynamics is based on network theory using Protein Structure Networks (PSNs). Despite their broad applicability, few efforts have been devoted to benchmarking PSN methods and to provide the community with best practices. In many applications, it is convenient to use the centers of mass of the side chains as nodes. It becomes thus critical to evaluate the minimal distance cutoff between the centers of mass which will provide stable network properties. Moreover, when the PSN is derived from a structural ensemble collected by molecular dynamics (MD), the impact of the MD force field has to be evaluated. We selected a dataset of proteins with different fold and size and assessed the two fundamental properties of the PSN, i.e. hubs and connected components. We identified an optimal cutoff of 5 Å that is robust to changes in the force field and the proteins. Our study builds solid foundations for the harmonization and standardization of the PSN approach.

The allosteric switching mechanism in bacteriophage MS2

We use all-atom simulations to elucidate the mechanisms underlying conformational switching and allostery within the coat protein of the bacteriophage MS2. Assembly of most icosahedral virus capsids requires that the capsid protein adopts different conformations at precise locations within the capsid. It has been shown that a 19 nucleotide stem loop (TR) from the MS2 genome acts as an allosteric effector, guiding conformational switching of the coat protein during capsid assembly. Since the principal conformational changes occur far from the TR binding site, it is important to understand the molecular mechanism underlying this allosteric communication. To this end, we use all-atom simulations with explicit water combined with a path sampling technique to sample the MS2 coat protein conformational transition, in the presence and absence of TR-binding. The calculations find that TR binding strongly alters the transition free energy profile, leading to a switch in the favored conformation. We discuss changes in molecular interactions responsible for this shift. We then identify networks of amino acids with correlated motions to reveal the mechanism by which effects of TR binding span the protein. We find that TR binding strongly affects residues located at the 5-fold and quasi-sixfold interfaces in the assembled capsid, suggesting a mechanism by which the TR binding could direct formation of the native capsid geometry. The analysis predicts amino acids whose substitution by mutagenesis could alter populations of the conformational substates or their transition rates.

Efficient preparation and analysis of membrane and membrane protein systems

Molecular dynamics (MD) simulations have become a highly important technique to consider lipid membrane systems, and quite often they provide considerable added value to laboratory experiments. Rapid development of both software and hardware has enabled the increase of time and size scales reachable by MD simulations to match those attainable by several accurate experimental techniques. However, until recently, the quality and maturity of software tools available for building membrane models for simulations as well as analyzing the results of these simulations have seriously lagged behind. Here, we discuss the recent developments of such tools from the end-users' point of view. In particular, we review the software that can be employed to build lipid bilayers and other related structures with or without embedded membrane proteins to be employed in MD simulations. Additionally, we provide a brief critical insight into force fields and MD packages commonly used for membrane and membrane protein simulations. Finally, we list analysis tools that can be used to study the properties of membrane and membrane protein systems. In all these points we comment on the respective compatibility of the covered tools. We also share our opinion on the current state of the available software. We briefly discuss the most commonly employed tools and platforms on which new software can be built. We conclude the review by providing a few ideas and guidelines on how the development of tools can be further boosted to catch up with the rapid pace at which the field of membrane simulation progresses. This includes improving the compatibility between software tools and promoting the openness of the codes on which these applications rely. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.