Generalized correlation-based dynamical network analysis: a new high-performance approach for identifying allosteric communications in molecular dynamics trajectories

Abstract

Molecular interactions are essential for regulation of cellular processes, from the formation of multiprotein complexes, to the allosteric activation of enzymes. Identifying the essential residues and molecular features that regulate such interactions is paramount for understanding the biochemical process in question, allowing for suppression of a reaction through drug interventions, or optimization of a chemical process using bioengineered molecules. In order to identify important residues and information pathways within molecular complexes, the Dynamical Network Analysis method was developed and has since been broadly applied in the literature. However, in the dawn of exascale computing, this method is generally limited to relatively small biomolecular systems. In this work we provide an evolution of the method, application and interface. All data processing and analysis is conducted through Jupyter notebooks, providing automatic detection of important solvent and ion residues, an optimized and parallel generalized correlation implementation that is linear with respect to the number of nodes in the system, and subsequent community clustering, calculation of betweenness of contacts, and determination optimal paths. Using the popular visualization program VMD, high-quality renderings of the networks over the biomolecular structures can be produced. Our new implementation was employed to investigate three different systems, with up to 2.5 M atoms, namely the OMP-decarboxylase, the Leucyl-tRNA synthetase complexed with its cognate tRNA and adenylate, and the respiratory complex I in a membrane environment. Our enhanced and updated protocol provides the community with an intuitive and interactive interface, which can be easily applied to large macromolecular complexes.1/25

Markov state models and NMR uncover an overlooked allosteric loop in p53

“… 66 The allosteric communication within the DBD can be also further investigated using other recent computational approaches such as the energy decomposition method 67–69 and the dynamical network analysis. 70 …”

The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: the case of the spike A222V mutation

“…The dynamic network analysis was performed using an implementation by Melo et al . (Melo et al 2020), following a standard set of parameters, i.e. a contact cut-off of 4.5 Å and contact persistence of 0.5.…”

How Communication Pathways Bridge Local and Global Conformations in an IgG4 Antibody: a Molecular Dynamics Study

“…51 The number of hydrogen bonds was computed with VMD, 66 setting an angle cutoff of 25°and a maximum donor-acceptor distance of 3Å. Mutual information, generalized correlation coefficient and correlation score were computed with in-house scripts, and network analysis was performed using the scripts produced by Melo et al 67 For a detailed description of the methods used for data analysis, the reader is referred to the Supplementary Information. and RMSD of antigen; MM/PBSA calculations of antigen/pembrolizumab binding; distributions of the generalized correlation coefficient; intra-and inter-domain correlation score; RMSD values of binding site in clusters 0 H , 1 A , and 0 H after antigen removal; antigen-pembrolizumab mutual information.…”