Characterizing protein conformations by correlation analysis of coarse-grained contact matrices

Computational Exploration of Potential CFTR Binding Sites for Type I Corrector Drugs

Cystic fibrosis (CF) is a recessive genetic disease that is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) protein. The recent development of a class of drugs called "correctors", which repair the structure and function of mutant CFTR, has greatly enhanced the life expectancy of CF patients. These correctors target the most common disease causing CFTR mutant F508del and are exemplified by the FDA-approved VX-809. While one binding site of VX-809 to CFTR was recently elucidated by cryo-electron microscopy, four additional binding sites have been proposed in the literature and it has been theorized that VX-809 and structurally similar correctors may engage multiple CFTR binding sites. To explore these five binding sites, ensemble docking was performed on wild-type CFTR and the F508del mutant using a large library of structurally similar corrector drugs, including VX-809 (lumacaftor), VX-661 (tezacaftor), ABBV-2222 (galicaftor), and a host of other structurally related molecules. For wild-type CFTR, we find that only one site, located in membrane spanning domain 1 (MSD1), binds favorably to our ligand library. While this MSD1 site also binds our ligand library for F508del-CFTR, the F508del mutation also opens a binding site in nucleotide binding domain 1 (NBD1), which enables strong binding of our ligand library to this site. This NBD1 site in F508del-CFTR exhibits the strongest overall binding affinity for our library of corrector drugs. This data may serve to better understand the structural changes induced by mutation of CFTR and how correctors bind to the protein. Additionally, it may aid in the design of new, more effective CFTR corrector drugs.

Protein conformational switch discerned via network centrality properties

Network analysis has emerged as a powerful tool for examining structural biology systems. The spatial organization of the components of a biomolecular structure has been rendered as a graph representation and analyses have been performed to deduce the biophysical and mechanistic properties of these components. For proteins, the analysis of protein structure networks (PSNs), especially via network centrality measurements and cluster coefficients, has led to identifying amino acid residues that play key functional roles and classifying amino acid residues in general. Whether these network properties examined in various studies are sensitive to subtle (yet biologically significant) conformational changes remained to be addressed. Here, we focused on four types of network centrality properties (betweenness, closeness, degree, and eigenvector centralities) for conformational changes upon ligand binding of a sensor protein (constitutive androstane receptor) and an allosteric enzyme (ribonucleotide reductase). We found that eigenvector centrality is sensitive and can distinguish salient structural features between protein conformational states while other centrality measures, especially closeness centrality, are less sensitive and rather generic with respect to the structural specificity. We also demonstrated that an ensemble-informed, modified PSN with static edges removed (which we term PSN*) has enhanced sensitivity at discerning structural changes.

Dynamics, a Powerful Component of Current and Future in Silico Approaches for Protein Design and Engineering

Computational prediction has become an indispensable aid in the processes of engineering and designing proteins for various biotechnological applications. With the tremendous progress in more powerful computer hardware and more efficient algorithms, some of in silico tools and methods have started to apply the more realistic description of proteins as their conformational ensembles, making protein dynamics an integral part of their prediction workflows. To help protein engineers to harness benefits of considering dynamics in their designs, we surveyed new tools developed for analyses of conformational ensembles in order to select engineering hotspots and design mutations. Next, we discussed the collective evolution towards more flexible protein design methods, including ensemble-based approaches, knowledge-assisted methods, and provable algorithms. Finally, we highlighted apparent challenges that current approaches are facing and provided our perspectives on their further development.

Using contact statistics to characterize structure transformation of biopolymer ensembles

As a unique subset of functional polymers, many biopolymers have a set of well-defined three-dimensional (3D) structural characteristics that can be described by spatial contacts between monomers. Statistical analysis of the contacts has been extremely productive in characterizing the biopolymer structural ensemble, such as for 3D chromosome structures. Often, native contacts and compartment structures are the focus of the studies, while the generic polymer aspect, such as the overall decaying of contacts with increasing sequence distance, is analyzed separately or preemptively removed. Here, we explore insights that can be gained by performing "compartment analysis" that keeps the distance decay, which we believe is particularly useful for characterizing the structure transformation of biopolymers. We tested contact analysis on several such transformations under physical perturbation or biological processes, including (1) unfolding of proteins induced by thermal denaturation, (2) chromosome conformation transition during the cell cycle, and (3) chromosome unpacking by physicochemical perturbations. Useful score functions were developed to further quantitatively characterize the transformation judging from the contact analysis. We also find that the sinusoidal undertone of eigenvector patterns (the "unwanted," low frequency signal, in contrast to the detailed A/B compartment) that had previously been attributed to biological effects of centromere proximal and distal interactions may in fact reflect a universal feature of polymers that have relatively weaker long-range contacts.

Disorder Mediated Oligomerization of DISC1 Proteins Revealed by Coarse-Grained Molecular Dynamics Simulations

Disrupted-in-schizophrenia-1 (DISC1) is a scaffold protein of significant importance for neuro-development and a prominent candidate protein in the etiology of mental disorders. In this work, we investigate the role of conformational heterogeneity and local structural disorder in the oligomerization pathway of the full-length DISC1 and of two truncation variants. Through extensive coarse-grained molecular dynamics simulations with a predictive energy landscape-based model, we shed light on the interplay of local and global disorder which lead to different oligomerization pathways. We found that both global conformational heterogeneity and local structural disorder play an important role in shaping distinct oligomerization trends of DISC1. This study also sheds light on the differences in oligomerization pathways of the full-length protein compared to the truncated variants produced by a chromosomal translocation associated with schizophrenia. We report that oligomerization of full-length DISC1 sequence works in a nonadditive manner with respect to truncated fragments that do not mirror the conformational landscape or binding affinities of the full-length unit.

Characterizing the 3D structure and dynamics of chromosomes and proteins in a common contact matrix framework

Conformational ensembles of biopolymers, whether proteins or chromosomes, can be described using contact matrices. Principal component analysis (PCA) on the contact data has been used to interrogate both protein and chromosome structures and/or dynamics. However, as these fields have developed separately, variants of PCA have emerged. Previously, a variant we hereby term Implicit-PCA (I-PCA) has been applied to chromosome contact matrices and revealed the spatial segregation of active and inactive chromatin. Separately, Explicit-PCA (E-PCA) has previously been applied to proteins and characterized their correlated structure fluctuations. Here, we swapped analysis methods (I-PCA and E-PCA), applying each to a different biopolymer type (chromosome or protein) than the one for which they were initially developed. We find that applying E-PCA to chromosome distance matrices derived from microscopy data can reveal the dominant motion (concerted fluctuation) of these chromosomes. Further, by applying E-PCA to Hi-C data across the human blood cell lineage, we isolated the aspects of chromosome structure that most strongly differentiate cell types. Conversely, when we applied I-PCA to simulation snapshots of proteins, the major component reported the consensus features of the structure, making this a promising approach for future analysis of semi-structured proteins.

Effector-Binding-Directed Dimerization and Dynamic Communication between Allosteric Sites of Ribonucleotide Reductase

Proteins forming dimers or larger complexes can be strongly influenced by their effector-binding status. We investigated how the effector-binding event is coupled with interface formation via computer simulations, and we quantified the correlation of two types of contact interactions: between the effector and its binding pocket and between protein monomers. This was achieved by connecting the protein dynamics at the monomeric level with the oligomer interface information. We applied this method to ribonucleotide reductase (RNR), an essential enzyme for de novo DNA synthesis. RNR contains two important allosteric sites, the s-site (specificity site) and the a-site (activity site), which bind different effectors. We studied these different binding states with atomistic simulation and used their coarse-grained contact information to analyze the protein dynamics. The results reveal that the effector-protein dynamics at the s-site and dimer interface formation are positively coupled. We further quantify the resonance level between these two events, which can be applied to other similar systems. At the a-site, different effector-binding states (ATP vs dATP) drastically alter the protein dynamics and affect the activity of the enzyme. On the basis of these results, we propose a new mechanism of how the a-site regulates enzyme activation.

CAMERRA: An analysis tool for the computation of conformational dynamics by evaluating residue-residue associations

A computational method which extracts the dominant motions from an ensemble of biomolecular conformations via a correlation analysis of residue-residue contacts is presented. The algorithm first renders the structural information into contact matrices, then constructs the collective modes based on the correlated dynamics of a selected set of dynamic contacts. Associated programs can bridge the results for further visualization using graphics software. The aim of this method is to provide an analysis of conformations of biopolymers from the contact viewpoint. It may assist a systematical uncovering of conformational switching mechanisms existing in proteins and biopolymer systems in general by statistical analysis of simulation snapshots. In contrast to conventional correlation analyses of Cartesian coordinates (such as distance covariance analysis and Cartesian principal component analysis), this program also provides an alternative way to locate essential collective motions in general. Herein, we detail the algorithm in a stepwise manner and comment on the importance of the method as applied to decoding allosteric mechanisms. © 2018 Wiley Periodicals, Inc.