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

Fast, accurate ranking of engineered proteins by target-binding propensity using structure modeling

Deep-learning-based methods for protein structure prediction have achieved unprecedented accuracy, yet their utility in the engineering of protein-based binders remains constrained due to a gap between the ability to predict the structures of candidate proteins and the ability toprioritize proteins by their potential to bind to a target. To bridge this gap, we introduce Automated Pairwise Peptide-Receptor Analysis for Screening Engineered proteins (APPRAISE), a method for predicting the target-binding propensity of engineered proteins. After generating structural models of engineered proteins competing for binding to a target using an established structure prediction tool such as AlphaFold-Multimer or ESMFold, APPRAISE performs a rapid (under 1 CPU second per model) scoring analysis that takes into account biophysical and geometrical constraints. As proof-of-concept cases, we demonstrate that APPRAISE can accurately classify receptor-dependent vs. receptor-independent adeno-associated viral vectors and diverse classes of engineered proteins such as miniproteins targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike, nanobodies targeting a G-protein-coupled receptor, and peptides that specifically bind to transferrin receptor or programmed death-ligand 1 (PD-L1). APPRAISE is accessible through a web-based notebook interface using Google Colaboratory (https://tiny.cc/APPRAISE). With its accuracy, interpretability, and generalizability, APPRAISE promises to expand the utility of protein structure prediction and accelerate protein engineering for biomedical applications.

Importance of Inter-residue Contacts for Understanding Protein Folding and Unfolding Rates, Remote Homology, and Drug Design

Inter-residue interactions in protein structures provide valuable insights into protein folding and stability. Understanding these interactions can be helpful in many crucial applications, including rational design of therapeutic small molecules and biologics, locating functional protein sites, and predicting protein-protein and protein-ligand interactions. The process of developing machine learning models incorporating inter-residue interactions has been improved recently. This review highlights the theoretical models incorporating inter-residue interactions in predicting folding and unfolding rates of proteins. Utilizing contact maps to depict inter-residue interactions aids researchers in developing computer models for detecting remote homologs and interface residues within protein-protein complexes which, in turn, enhances our knowledge of the relationship between sequence and structure of proteins. Further, the application of contact maps derived from inter-residue interactions is highlighted in the field of drug discovery. Overall, this review presents an extensive assessment of the significant models that use inter-residue interactions to investigate folding rates, unfolding rates, remote homology, and drug development, providing potential future advancements in constructing efficient computational models in structural biology.

Strengthening E-cadherin adhesion via antibody-mediated binding stabilization

E-cadherins (Ecads) are a crucial cell-cell adhesion protein with tumor suppression properties. Ecad adhesion can be enhanced by the monoclonal antibody 66E8, which has potential applications in inhibiting cancer metastasis. However, the biophysical mechanisms underlying 66E8-mediated adhesion strengthening are unknown. Here, we use molecular dynamics simulations, site-directed mutagenesis, and single-molecule atomic force microscopy experiments to demonstrate that 66E8 strengthens Ecad binding by stabilizing the primary Ecad adhesive conformation: the strand-swap dimer. By forming electrostatic interactions with Ecad, 66E8 stabilizes the swapped β-strand and its hydrophobic pocket and impedes Ecad conformational changes, which are necessary for rupture of the strand-swap dimer. Our findings identify fundamental mechanistic principles for strengthening of Ecad binding using monoclonal antibodies.

Drug Aggregation of Sparingly-Soluble Ionizable Drugs: Molecular Dynamics Simulations of Papaverine and Prostaglandin F2α

Aggregation in aqueous solution can have important implications on both the in vivo exposure of a drug and its pharmaceutical manufacturability. However, the drug aggregates formed can be very small and, thus, difficult to interrogate experimentally. On the other hand, at higher supersaturations where larger aggregates are supported, the chemical system is inherently metastable and therefore likewise challenging to study from an experimental standpoint. Understanding aggregation behavior is further complicated in the case of ionizable drugs where, unlike neutral compounds, there can be uncertainty in the kinds of drug molecules (i.e., charged, neutral, or both) that become incorporated into various clusters, particularly at pH values near the pKa. In this paper, we apply physics-based all-atom molecular dynamics (MD) simulations to study aggregation in the weakly basic drug papaverine and in the weakly acidic drug prostaglandin F2α. We employ in silico tools to construct simulation workflows and comprehensive cluster analysis protocols to elucidate the size distributions and dynamics of the drug aggregates formed at both an experimentally relevant concentration and at high supersaturation. We build on a previously published treatment [Solubility of sparingly soluble ionizable drugs. Adv. Drug Deliv. Rev. 2007, 59, 568-590, DOI: 10.1016/j.addr.2007.05.008] to translate the predicted aggregate distributions of each ionized drug into corresponding pH-solubility curves that can be compared directly to experiment. Our findings show that the assumption of a single predominant (charged) aggregate can be misleading in interpreting experimental pH-solubility curves, as it does not adequately reflect the rich diversity revealed in our simulations. Beyond not accounting for the distribution of ionized drug-containing clusters actually observed in solution, for both drugs we find evidence that neutral drug molecules can also participate in the aggregation phenomena. Notably, we observe that many drug molecules remain as free monomers in solution even under simulated conditions designed to mimic those where there is significant deviation of the experimental pH-solubility curve from the Henderson-Hasselbalch (H-H) equation, often taken to be a clear signpost of drug aggregation.

SARS-CoV-2 antibodies recognize 23 distinct epitopic sites on the receptor binding domain

The COVID-19 pandemic and SARS-CoV-2 variants have dramatically illustrated the need for a better understanding of antigen (epitope)-antibody (paratope) interactions. To gain insight into the immunogenic characteristics of epitopic sites (ES), we systematically investigated the structures of 340 Abs and 83 nanobodies (Nbs) complexed with the Receptor Binding Domain (RBD) of the SARS-CoV-2 spike protein. We identified 23 distinct ES on the RBD surface and determined the frequencies of amino acid usage in the corresponding CDR paratopes. We describe a clustering method for analysis of ES similarities that reveals binding motifs of the paratopes and that provides insights for vaccine design and therapies for SARS-CoV-2, as well as a broader understanding of the structural basis of Ab-protein antigen (Ag) interactions.

Connected Component Analysis of Dynamical Perturbation Contact Networks

Describing protein dynamical networks through amino acid contacts is a powerful way to analyze complex biomolecular systems. However, due to the size of the systems, identifying the relevant features of protein-weighted graphs can be a difficult task. To address this issue, we present the connected component analysis (CCA) approach that allows for fast, robust, and unbiased analysis of dynamical perturbation contact networks (DPCNs). We first illustrate the CCA method as applied to a prototypical allosteric enzyme, the imidazoleglycerol phosphate synthase (IGPS) enzyme from Thermotoga maritima bacteria. This approach was shown to outperform the clustering methods applied to DPCNs, which could not capture the propagation of the allosteric signal within the protein graph. On the other hand, CCA reduced the DPCN size, providing connected components that nicely describe the allosteric propagation of the signal from the effector to the active sites of the protein. By applying the CCA to the IGPS enzyme in different conditions, i.e., at high temperature and from another organism (yeast IGPS), and to a different enzyme, i.e., a protein kinase, we demonstrated how CCA of DPCNs is an effective and transferable tool that facilitates the analysis of protein-weighted networks.

PyInteraph2 and PyInKnife2 to Analyze Networks in Protein Structural Ensembles

Due to the complex nature of noncovalent interactions and their long-range effects, analyzing protein conformations using network theory can be enlightening. Protein Structure Networks (PSNs) provide a convenient formalism to study protein structures in relation to essential properties such as key residues for structural stability, allosteric communication, and the effects of modifications of the protein. PSNs can be defined according to very different principles, and the available tools have limitations in input formats, supported models, and version control. Other outstanding problems are related to the definition of network cutoffs and the assessment of the stability of the network properties. The protein science community could benefit from a common framework to carry out these analyses and make them easier to reproduce, reuse, and evaluate. We here provide two open-source software packages, PyInteraph2 and PyInKnife2, to implement and analyze PSNs in a reproducible and documented manner. PyInteraph2 interfaces with multiple formats for protein ensembles and incorporates different network models with the possibility of integrating them into a macronetwork and performing various downstream analyses, including hubs, connected components, and several other centrality measures, and visualizes the networks or further analyzes them thanks to compatibility with Cytoscape.PyInKnife2 that supports the network models implemented in PyInteraph2. It employs a jackknife resampling approach to estimate the convergence of network properties and streamline the selection of distance cutoffs. We foresee that the modular structure of the code and the supported version control system will promote the transition to a community-driven effort, boost reproducibility, and establish common protocols in the PSN field. As developers, we will guarantee the introduction of new functionalities and maintenance, assistance, and training of new contributors.

Strengthening E-cadherin adhesion via antibody mediated binding stabilization

E-cadherins (Ecads) are a crucial cell-cell adhesion protein with tumor suppression properties. Ecad adhesion can be enhanced by the monoclonal antibody 66E8, which has potential applications in inhibiting cancer metastasis. However, the biophysical mechanisms underlying 66E8 mediated adhesion strengthening are unknown. Here, we use molecular dynamics simulations, site directed mutagenesis and single molecule atomic force microscopy experiments to demonstrate that 66E8 strengthens Ecad binding by stabilizing the primary Ecad adhesive conformation: the strand-swap dimer. By forming electrostatic interactions with Ecad, 66E8 stabilizes the swapped β-strand and its hydrophobic pocket and impedes Ecad conformational changes, which are necessary for rupture of the strand-swap dimer. Our findings identify fundamental mechanistic principles for strengthening of Ecad binding using monoclonal antibodies.

The ClusPro AbEMap web server for the prediction of antibody epitopes

Antibodies play an important role in the immune system by binding to molecules called antigens at their respective epitopes. These interfaces or epitopes are structural entities determined by the interactions between an antibody and an antigen, making them ideal systems to analyze by using docking programs. Since the advent of high-throughput antibody sequencing, the ability to perform epitope mapping using only the sequence of the antibody has become a high priority. ClusPro, a leading protein-protein docking server, together with its template-based modeling version, ClusPro-TBM, have been re-purposed to map epitopes for specific antibody-antigen interactions by using the Antibody Epitope Mapping server (AbEMap). ClusPro-AbEMap offers three different modes for users depending on the information available on the antibody as follows: (i) X-ray structure, (ii) computational/predicted model of the structure or (iii) only the amino acid sequence. The AbEMap server presents a likelihood score for each antigen residue of being part of the epitope. We provide detailed information on the server's capabilities for the three options and discuss how to obtain the best results. In light of the recent introduction of AlphaFold2 (AF2), we also show how one of the modes allows users to use their AF2-generated antibody models as input. The protocol describes the relative advantages of the server compared to other epitope-mapping tools, its limitations and potential areas of improvement. The server may take 45-90 min depending on the size of the proteins.

Computational Insights into the Dynamic Structural Features and Binding Characteristics of Recombinase UvsX Compared with RecA

RecA family recombinases are the core enzymes in the process of homologous recombination, and their normal operation ensures the stability of the genome and the healthy development of organisms. The UvsX protein from bacteriophage T4 is a member of the RecA family recombinases and plays a central role in T4 phage DNA repair and replication, which provides an important model for the biochemistry and genetics of DNA metabolism. UvsX shares a high degree of structural similarity and function with RecA, which is the most deeply studied member of the RecA family. However, the detailed molecular mechanism of UvsX has not been resolved. In this study, a comprehensive all-atom molecular dynamics simulation of the UvsX protein dimer complex was carried out in order to investigate the conformational and binding properties of UvsX in combination with ATP and DNA, and the simulation of RecA was synchronized with the property comparison learning for UvsX. This study confirmed the highly conserved molecular structure characteristics and catalytic centers of RecA and UvsX, and also discovered differences in regional conformation, volatility and the ability to bind DNA between the two proteins at different temperatures, which would be helpful for the subsequent understanding and application of related recombinases.

Mapping of antibody epitopes based on docking and homology modeling

Antibodies are key proteins produced by the immune system to target pathogen proteins termed antigens via specific binding to surface regions called epitopes. Given an antigen and the sequence of an antibody the knowledge of the epitope is critical for the discovery and development of antibody based therapeutics. In this work, we present a computational protocol that uses template-based modeling and docking to predict epitope residues. This protocol is implemented in three major steps. First, a template-based modeling approach is used to build the antibody structures. We tested several options, including generation of models using AlphaFold2. Second, each antibody model is docked to the antigen using the fast Fourier transform (FFT) based docking program PIPER. Attention is given to optimally selecting the docking energy parameters depending on the input data. In particular, the van der Waals energy terms are reduced for modeled antibodies relative to x-ray structures. Finally, ranking of antigen surface residues is produced. The ranking relies on the docking results, that is, how often the residue appears in the docking poses' interface, and also on the energy favorability of the docking pose in question. The method, called PIPER-Map, has been tested on a widely used antibody-antigen docking benchmark. The results show that PIPER-Map improves upon the existing epitope prediction methods. An interesting observation is that epitope prediction accuracy starting from antibody sequence alone does not significantly differ from that of starting from unbound (i.e., separately crystallized) antibody structure.

psnGPCRdb: The Structure-network Database of G Protein Coupled Receptors

G protein coupled receptors (GPCRs) are critical eukaryotic signal transduction gatekeepers and represent the largest protein superfamily in the human proteome, with more than 800 members. They share seven transmembrane helices organized in an up-down bundle architecture. GPCR-mediated signaling pathways have been linked to numerous human diseases, and GPCRs are the targets of approximately 35% of all drugs currently on the market. Structure network analysis, a graph theory-based approach, represents a cutting-edge tool to deeply understand GPCR function, which strongly relies on communication between the extracellular and intracellular poles of their structure. psnGPCRdb stores the structure networks (i.e., linked nodes, hubs, communities and communication pathways) computed on all updated GPCR structures in the Protein Data Bank, in their isolated states or in complex with extracellular and/or intracellular molecules. The structure communication signatures of a sub-family or family of GPCRs as well as of their small-molecule activators or inhibitors are stored as consensus networks. The database stores also all meaningful structure network-based comparisons (i.e., difference networks) of functionally different states (i.e., inactive or active) of a given receptor sub-type, or of consensus networks representative of a receptor sub-type, type, sub-family or family. Single or consensus GPCR networks hold also information on amino acid conservation. The database allows to graphically analyze 3D structure networks together with interactive data-tables. Ligand-centric networks can be analyzed as well. psnGPCRdb is unique and represents a powerful resource to unravel GPCR function with important implications in cell signaling and drug design. psnGPCRdb is freely available at: http://webpsn.hpc.unimo.it/psngpcr.php.

A correlation coefficient-based feature selection approach for virus-host protein-protein interaction prediction

Prediction of virus-host protein-protein interactions (PPI) is a broad research area where various machine-learning-based classifiers are developed. Transforming biological data into machine-usable features is a preliminary step in constructing these virus-host PPI prediction tools. In this study, we have adopted a virus-host PPI dataset and a reduced amino acids alphabet to create tripeptide features and introduced a correlation coefficient-based feature selection. We applied feature selection across several correlation coefficient metrics and statistically tested their relevance in a structural context. We compared the performance of feature-selection models against that of the baseline virus-host PPI prediction models created using different classification algorithms without the feature selection. We also tested the performance of these baseline models against the previously available tools to ensure their predictive power is acceptable. Here, the Pearson coefficient provides the best performance with respect to the baseline model as measured by AUPR; a drop of 0.003 in AUPR while achieving a 73.3% (from 686 to 183) reduction in the number of tripeptides features for random forest. The results suggest our correlation coefficient-based feature selection approach, while decreasing the computation time and space complexity, has a limited impact on the prediction performance of virus-host PPI prediction tools.

The Importance of Charge Transfer and Solvent Screening in the Interactions of Backbones and Functional Groups in Amino Acid Residues and Nucleotides

Quantum mechanical (QM) calculations at the level of density-functional tight-binding are applied to a protein–DNA complex (PDB: 2o8b) consisting of 3763 atoms, averaging 100 snapshots from molecular dynamics simulations. A detailed comparison of QM and force field (Amber) results is presented. It is shown that, when solvent screening is taken into account, the contributions of the backbones are small, and the binding of nucleotides in the double helix is governed by the base–base interactions. On the other hand, the backbones can make a substantial contribution to the binding of amino acid residues to nucleotides and other residues. The effect of charge transfer on the interactions is also analyzed, revealing that the actual charge of nucleotides and amino acid residues can differ by as much as 6 and 8% from the formal integer charge, respectively. The effect of interactions on topological models (protein -residue networks) is elucidated.

Connecting the dots: A practical evaluation of web-tools for describing protein dynamics as networks

Protein Structure Networks (PSNs) are a well-known mathematical model for estimation and analysis of the three-dimensional protein structure. Investigating the topological architecture of PSNs may help identify the crucial amino acid residues for protein stability and protein-protein interactions, as well as deduce any possible mutational effects. But because proteins go through conformational changes to give rise to essential biological functions, this has to be done dynamically over time. The most effective method to describe protein dynamics is molecular dynamics simulation, with the most popular software programs for manipulating simulations to infer interaction networks being RING, MD-TASK, and NAPS. Here, we compare the computational approaches used by these three tools-all of which are accessible as web servers-to understand the pathogenicity of missense mutations and talk about their potential applications as well as their advantages and disadvantages.

Conformation-Dependent Influences of Hydrophobic Amino Acids in Two In-Register Parallel β-Sheet Amyloids, an α-Synuclein Amyloid and a Local Structural Model of PrPSc

Prions are unconventional pathogens that encode the pathogenic information in conformations of the constituent abnormal isoform of prion protein (PrP^Sc), independently of the nucleotide genome. Therefore, conformational diversity of PrP^Sc underlies the existence of many prion strains and species barriers of prions, although the conformational information is extremely limited. Interestingly, differences between polymorphic or species-specific residues responsible for the species/strain barriers are often caused by conservative replacements between hydrophobic amino acids. This implies that subtle differences among hydrophobic amino acids are significant for PrP^Sc structures. Here we analyzed the influence of different hydrophobic residues on the structures of an in-register parallel β-sheet amyloid of α-synuclein (αSyn) using molecular dynamics (MD) simulation and applied the knowledge from the αSyn amyloid to modeling a local structure of human PrP^Sc encompassing residues 107-143. We found that mutations equivalent to polymorphisms that cause transmission barriers substantially affect the stabilities of the local structures; for example, the G127V mutation, which makes the host resistant to various human prion diseases, greatly destabilized the local structure of the model amyloid. Our study indicates that subtle differences among hydrophobic side chains can considerably affect the interaction network, including hydrogen bonds, and demonstrates specifically how and in what structures hydrophobic residues can exert unique effects on in-register parallel β-sheet amyloids.

Dynamic stability of salt stable cowpea chlorotic mottle virus capsid protein dimers and pentamers of dimers

Intermediates of the self-assembly process of the salt stable cowpea chlorotic mottle virus (ss-CCMV) capsid can be modelled atomistically on realistic computational timescales either by studying oligomers in equilibrium or by focusing on their dissociation instead of their association. Our previous studies showed that among the three possible dimer interfaces in the icosahedral capsid, two are thermodynamically relevant for capsid formation. The aim of the current study is to evaluate the relative structural stabilities of the three different ss-CCMV dimers and to find and understand the conditions that lead to their dissociation. Long timescale molecular dynamics simulations at 300 K of the various dimers and of the pentamer of dimers underscore the importance of large contact surfaces on stabilizing the capsid subunits within an oligomer. Simulations in implicit solvent show that at higher temperature (350 K), the N-terminal tails of the protein units act as tethers, delaying dissociation for all but the most stable interface. The pentamer of dimers is also found to be stable on long timescales at 300 K, with an inherent flexibility of the outer protein chains.

Novel protein contact points among TP53 and minichromosome maintenance complex proteins 2, 3, and 5

OBJECTIVE

Identify protein contact points between TP53 and minichromosome maintenance (MCM) complex proteins 2, 3, and 5 with high resolution allowing for potential novel Cancer drug design.

METHODS

A next-generation sequencing-based protein-protein interaction method developed in our laboratory called AVA-Seq was applied to a gold-standard human protein interaction set. Proteins including TP53, MCM2, MCM3, MCM5, HSP90AA1, PCNA, NOD1, and others were sheared and ligated into the AVA-Seq system. Protein-protein interactions were then identified in both mild and stringent selective conditions.

RESULTS

Known interactions among MCM2, MCM3, and MCM5 were identified with the AVA-Seq system. The interacting regions detected between these three proteins overlap with the structural data of the MCM complex, and novel domains were identified with high resolution determined by multiple overlapping fragments. Fragments of wild type TP53 were shown to interact with MCM2, MCM3, and MCM5, and details on the location of the interactions were provided. Finally, a mini-network of known and novel cancer protein interactions was provided, which could have implications for fundamental changes in multiple cancers.

CONCLUSION

We provide a high-resolution mini-interactome that could direct novel drug targets and implicate possible effects of specific cancer mutations.

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.

Substrate Access Mechanism in a Novel Membrane-Bound Phospholipase A of Pseudomonas aeruginosa Concordant with Specificity and Regioselectivity

PlaF is a cytoplasmic membrane-bound phospholipase A₁ from Pseudomonas aeruginosa that alters the membrane glycerophospholipid (GPL) composition and fosters the virulence of this human pathogen. PlaF activity is regulated by a dimer-to-monomer transition followed by tilting of the monomer in the membrane. However, how substrates reach the active site and how the characteristics of the active site tunnels determine the activity, specificity, and regioselectivity of PlaF for natural GPL substrates have remained elusive. Here, we combined unbiased and biased all-atom molecular dynamics (MD) simulations and configurational free-energy computations to identify access pathways of GPL substrates to the catalytic center of PlaF. Our results map out a distinct tunnel through which substrates access the catalytic center. PlaF variants with bulky tryptophan residues in this tunnel revealed decreased catalysis rates due to tunnel blockage. The MD simulations suggest that GPLs preferably enter the active site with the sn-1 acyl chain first, which agrees with the experimentally demonstrated PLA₁ activity of PlaF. We propose that the acyl chain-length specificity of PlaF is determined by the structural features of the access tunnel, which results in favorable free energy of binding of medium-chain GPLs. The suggested egress route conveys fatty acid (FA) products to the dimerization interface and, thus, contributes to understanding the product feedback regulation of PlaF by FA-triggered dimerization. These findings open up opportunities for developing potential PlaF inhibitors, which may act as antibiotics against P. aeruginosa.

Mapping Function from Dynamics: Future Challenges for Network-Based Models of Protein Structures

Proteins fulfill complex and diverse biological functions through the controlled atomic motions of their structures (functional dynamics). The protein composition is given by its amino-acid sequence, which was assumed to encode the function. However, the discovery of functional sequence variants proved that the functional encoding does not come down to the sequence, otherwise a change in the sequence would mean a change of function. Likewise, the discovery that function is fulfilled by a set of structures and not by a unique structure showed that the functional encoding does not come down to the structure either. That leaves us with the possibility that a set of atomic motions, achievable by different sequences and different structures, encodes a specific function. Thanks to the exponential growth in annual depositions in the Protein Data Bank of protein tridimensional structures at atomic resolutions, network models using the Cartesian coordinates of atoms of a protein structure as input have been used over 20 years to investigate protein features. Combining networks with experimental measures or with Molecular Dynamics (MD) simulations and using typical or ad-hoc network measures is well suited to decipher the link between protein dynamics and function. One perspective is to consider static structures alone as alternatives to address the question and find network measures relevant to dynamics that can be subsequently used for mining and classification of dynamic sequence changes functionally robust, adaptable or faulty. This way the set of dynamics that fulfill a function over a diversity of sequences and structures will be determined.

Mutation Edgotype Drives Fitness Effect in Human

Missense mutations are known to perturb protein-protein interaction networks (known as interactome networks) in different ways. However, it remains unknown how different interactome perturbation patterns (“edgotypes”) impact organismal fitness. Here, we estimate the fitness effect of missense mutations with different interactome perturbation patterns in human, by calculating the fractions of neutral and deleterious mutations that do not disrupt PPIs (“quasi-wild-type”), or disrupt PPIs either by disrupting the binding interface (“edgetic”) or by disrupting overall protein stability (“quasi-null”). We first map pathogenic mutations and common non-pathogenic mutations onto homology-based three-dimensional structural models of proteins and protein-protein interactions in human. Next, we perform structure-based calculations to classify each mutation as either quasi-wild-type, edgetic, or quasi-null. Using our predicted as well as experimentally determined interactome perturbation patterns, we estimate that >∼40% of quasi-wild-type mutations are effectively neutral and the remaining are mostly mildly deleterious, that >∼75% of edgetic mutations are only mildly deleterious, and that up to ∼75% of quasi-null mutations may be strongly detrimental. These estimates are the first such estimates of fitness effect for different network perturbation patterns in any interactome. Our results suggest that while mutations that do not disrupt the interactome tend to be effectively neutral, the majority of human PPIs are under strong purifying selection and the stability of most human proteins is essential to human life.

pyProGA-A PyMOL plugin for protein residue network analysis

The field of protein residue network (PRN) research has brought several useful methods and techniques for structural analysis of proteins and protein complexes. Many of these are ripe and ready to be used by the proteomics community outside of the PRN specialists. In this paper we present software which collects an ensemble of (network) methods tailored towards the analysis of protein-protein interactions (PPI) and/or interactions of proteins with ligands of other type, e.g. nucleic acids, oligosaccharides etc. In parallel, we propose the use of the network differential analysis as a method to identify residues mediating key interactions between proteins. We use a model system, to show that in combination with other, already published methods, also included in pyProGA, it can be used to make such predictions. Such extended repertoire of methods allows to cross-check predictions with other methods as well, as we show here. In addition, the possibility to construct PRN models from various kinds of input is so far a unique asset of our code. One can use structural data as defined in PDB files and/or from data on residue pair interaction energies, either from force-field parameters or fragment molecular orbital (FMO) calculations. pyProGA is a free open-source software available from https://gitlab.com/Vlado_S/pyproga.

The Odd Faces of Oligomers: The Case of TRAF2-C, A Trimeric C-Terminal Domain of TNF Receptor-Associated Factor

TNF Receptor Associated Factor 2 (TRAF2) is a trimeric protein that belongs to the TNF receptor associated factor family (TRAFs). The TRAF2 oligomeric state is crucial for receptor binding and for its interaction with other proteins involved in the TNFR signaling. The monomer-trimer equilibrium of a C- terminal domain truncated form of TRAF2 (TRAF2-C), plays also a relevant role in binding the membrane, causing inward vesiculation. In this study, we have investigated the conformational dynamics of TRAF2-C through circular dichroism, fluorescence, and dynamic light scattering, performing temperature-dependent measurements. The data indicate that the protein retains its oligomeric state and most of its secondary structure, while displaying a significative increase in the heterogeneity of the tyrosines signal, increasing the temperature from ≈15 to ≈35 °C. The peculiar crowding of tyrosine residues (12 out of 18) at the three subunit interfaces and the strong dependence on the trimer concentration indicate that such conformational changes mainly involve the contact areas between each pair of monomers, affecting the oligomeric state. Molecular dynamic simulations in this temperature range suggest that the interfaces heterogeneity is an intrinsic property of the trimer that arises from the continuous, asymmetric approaching and distancing of its subunits. Such dynamics affect the results of molecular docking on the external protein surface using receptor peptides, indicating that the TRAF2-receptor interaction in the solution might not involve three subunits at the same time, as suggested by the static analysis obtainable from the crystal structure. These findings shed new light on the role that the TRAF2 oligomeric state might have in regulating the protein binding activity in vivo.

Binding mode characterization of 13b in the monomeric and dimeric states of SARS-CoV-2 main protease using molecular dynamics simulations

The main protease, M^pro/3CL^pro, plays an essential role in processing polyproteins translated from viral RNA to produce functional viral proteins and therefore serve as an attractive target for discovering COVID-19 therapeutics. The availability of both monomer and dimer crystal bound with a common ligand, '13b' (α-ketoamide inhibitor), opened up opportunities to understand the M^pro mechanism of action. A comparative analysis of both forms of M^pro was carried out to elucidate the binding site architectural differences in the presence and absence of '13b'. Molecular dynamics simulations suggest that the presence of '13b' enhances the stability of M^pro than the unbound APO form. The N- and C- terminals of both the protomers stabilize each other, and making it's interface essential for the active form of M^pro. In comparison to monomer, the relatively high affinity of '13b' is gained in dimer pocket due to the high stability of the pocket by the interaction of S1 residue of chain B with residues F140, E166 and H172 of chain A, which is absent in monomer. The comprehensive essential dynamics, protein structure network analysis and thermodynamic profiling highlight the hot-spots, pivotal in molecular recognition process at protein-ligand and protein-protein interaction levels, cross-validated through computational alanine scanning study. A comparative description of '13b' binding mechanism in both forms illustrates valuable insights into the inhibition mechanism and the selection of critical residues suitable for the structure-based approaches for the identification of more potent M^pro inhibitors.Communicated by Ramaswamy H. Sarma.

Bcl-xL Dynamics under the Lens of Protein Structure Networks

Understanding the finely orchestrated interactions leading to or preventing programmed cell death (apoptosis) is of utmost importance in cancer research because the failure of these systems could eventually lead to the onset of the disease. In this regard, the maintenance of a delicate balance between the promoters and inhibitors of mitochondrial apoptosis is crucial, as demonstrated by the interplay among the Bcl-2 family members. In particular, B-cell lymphoma extra-large (Bcl-x_L) is a target of interest due to the forefront role of its dysfunctions in cancer development. Bcl-x_L prevents apoptosis by binding both the pro-apoptotic BH3-only proteins, like PUMA, and the noncanonical partners, such as p53, at different sites. An allosteric communication between the BH3-only protein binding pocket and the p53 binding site, mediating the release of p53 from Bcl-x_L upon PUMA binding, has been postulated and supported by nuclear magnetic resonance and other biophysical data. The molecular details of this mechanism, especially at the residue level, remain unclear. In this work, we investigated the distal communication between these two sites in Bcl-x_L in its free state and when bound to PUMA. We also evaluated how missense mutations of Bcl-x_L found in cancer samples might impair this communication and therefore the allosteric mechanism. We employed all-atom explicit solvent microsecond molecular dynamics simulations, analyzed through a Protein Structure Network approach and integrated with calculations of changes in free energies upon cancer-related mutations identified by genomics studies. We found a subset of candidate residues responsible for both maintaining protein stability and for conveying structural information between the two binding sites and hypothesized possible communication routes between specific residues at both sites.

A compact vocabulary of paratope-epitope interactions enables predictability of antibody-antigen binding

Antibody-antigen binding relies on the specific interaction of amino acids at the paratope-epitope interface. The predictability of antibody-antigen binding is a prerequisite for de novo antibody and (neo-)epitope design. A fundamental premise for the predictability of antibody-antigen binding is the existence of paratope-epitope interaction motifs that are universally shared among antibody-antigen structures. In a dataset of non-redundant antibody-antigen structures, we identify structural interaction motifs, which together compose a commonly shared structure-based vocabulary of paratope-epitope interactions. We show that this vocabulary enables the machine learnability of antibody-antigen binding on the paratope-epitope level using generative machine learning. The vocabulary (1) is compact, less than 10⁴ motifs; (2) distinct from non-immune protein-protein interactions; and (3) mediates specific oligo- and polyreactive interactions between paratope-epitope pairs. Our work leverages combined structure- and sequence-based learning to demonstrate that machine-learning-driven predictive paratope and epitope engineering is feasible.

A graph-based algorithm for detecting rigid domains in protein structures

Background

Conformational transitions are implicated in the biological function of many proteins. Structural changes in proteins can be described approximately as the relative movement of rigid domains against each other. Despite previous efforts, there is a need to develop new domain segmentation algorithms that are capable of analysing the entire structure database efficiently and do not require the choice of protein-dependent tuning parameters such as the number of rigid domains.

Results

We develop a graph-based method for detecting rigid domains in proteins. Structural information from multiple conformational states is represented by a graph whose nodes correspond to amino acids. Graph clustering algorithms allow us to reduce the graph and run the Viterbi algorithm on the associated line graph to obtain a segmentation of the input structures into rigid domains. In contrast to many alternative methods, our approach does not require knowledge about the number of rigid domains. Moreover, we identified default values for the algorithmic parameters that are suitable for a large number of conformational ensembles. We test our algorithm on examples from the DynDom database and illustrate our method on various challenging systems whose structural transitions have been studied extensively.

Conclusions

The results strongly suggest that our graph-based algorithm forms a novel framework to characterize structural transitions in proteins via detecting their rigid domains. The web server is available at http://azifi.tz.agrar.uni-goettingen.de/webservice/.

Disclosing Allostery Through Protein Contact Networks

Proteins are located in the twilight zone between chemistry and biology, where a peculiar kind of complexity starts. Proteins are the smallest 'devices' showing a sensible adaptation to their environment by the production of appropriate behavior when facing a specific stimulus. This fact qualifies (from the 'effector' side) proteins as nanomachines working as catalysts, motors, or switches. However (from the sensor side), the need to single out the 'specific stimulus' out of thermal noise qualifies proteins as information processing devices. Allostery corresponds to the modification of the configuration (in a broad sense) of the protein molecule in response to a specific stimulus in a non-strictly local way, thereby connecting the sensor and effector sides of the nanomachine. This is why the 'disclosing' of allostery phenomenon is at the very heart of protein function; in this chapter, we will demonstrate how a network-based representation of protein structure in terms of nodes (aminoacid residues) and edges (effective contacts between residues) is the natural language for getting rid of allosteric phenomena and, more in general, of protein structure/function relationships.

FOXP1 Interacts with MyoD to Repress its Transcription and Myoblast Conversion

Forkhead transcription factors (TFs) often dimerize outside their extensive family, whereas bHLH transcription factors typically dimerize with E12/E47. Based on structural similarities, we predicted that a member of the former, Forkhead Box P1 (FOXP1), might heterodimerize with a member of the latter, MYOD1 (MyoD). Data shown here support this hypothesis and further demonstrate the specificity of this forkhead/myogenic interaction among other myogenic regulatory factors. We found that FOXP1-MyoD heterodimerization compromises the ability of MyoD to bind to E-boxes and to transactivate E box- containing promoters. We observed that FOXP1 is required for the full ability of MyoD to convert fibroblasts into myotubules. We provide a model in which FOXP1 displaces ID and E12/E47 to repress MyoD during the proliferative phase of myoblast differentiation. These data identify FOXP1 as a hitherto unsuspected transcriptional repressor of MyoD. We suggest that isolation of paired E-box and forkhead sites within 1 turn helical spacings provides potential for cooperative interactions among heretofore distinct classes of transcription factors.

Investigating Conformational Dynamics and Allostery in the p53 DNA-Binding Domain Using Molecular Simulations

The p53 tumor suppressor is a multifaceted context-dependent protein, which is involved in multiple cellular pathways, with the ability to either keep the cells alive or to kill them through mechanisms such as apoptosis. To complicate this picture, cancer cells that express mutant p53 becomes addicted to the mutant activity, so that the mutant variant features a myriad of gain-of-function activities, opening different venues for therapy. This makes essential to think outside the box and apply new approaches to the study of p53 structure-(mis)function relationship to find new critical components of its pathway or to understand how known parts are interconnected, compete, or cooperate. In this context, I will here illustrate how to integrate different computational methods to the identification of possible allosteric effects transmitted from the DNA binding interface of p53 to regions for cofactor recruitment. The protocol can be extended to any other cases of study. Indeed, it does not necessarily apply only to the study of DNA-induced effects, but more broadly to the investigation of long-range effects induced by a biological partner that binds to a biomolecule of interest.

Surveying the Side-Chain Network Approach to Protein Structure and Dynamics: The SARS-CoV-2 Spike Protein as an Illustrative Case

Network theory-based approaches provide valuable insights into the variations in global structural connectivity between different dynamical states of proteins. Our objective is to review network-based analyses to elucidate such variations, especially in the context of subtle conformational changes. We present technical details of the construction and analyses of protein structure networks, encompassing both the non-covalent connectivity and dynamics. We examine the selection of optimal criteria for connectivity based on the physical concept of percolation. We highlight the advantages of using side-chain-based network metrics in contrast to backbone measurements. As an illustrative example, we apply the described network approach to investigate the global conformational changes between the closed and partially open states of the SARS-CoV-2 spike protein. These conformational changes in the spike protein is crucial for coronavirus entry and fusion into human cells. Our analysis reveals global structural reorientations between the two states of the spike protein despite small changes between the two states at the backbone level. We also observe some differences at strategic locations in the structures, correlating with their functions, asserting the advantages of the side-chain network analysis. Finally, we present a view of allostery as a subtle synergistic-global change between the ligand and the receptor, the incorporation of which would enhance drug design strategies.

Topology Results on Adjacent Amino Acid Networks of Oligomeric Proteins

In this chapter, we focus on topology measurements of the adjacent amino acid networks for a data set of oligomeric proteins and some of its subnetworks. The aim is to present many mathematical tools in order to understand the structures of proteins implicitly coded in such networks and subnetworks. We mainly investigate four important networks by computing the number of connected components, the degree distribution, and assortativity measures. We compare each result in order to prove that the four networks have quite independent topologies.

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.

Identification of potential allosteric binding sites in cathepsin K based on intramolecular communication

Network theory methods and molecular dynamics (MD) simulations are accepted tools to study allosteric regulation. Indeed, dynamic networks built upon correlation analysis of MD trajectories provide detailed information about communication paths between distant sites. In this context, we aimed to understand whether the efficiency of intramolecular communication could be used to predict the allosteric potential of a given site. To this end, we performed MD simulations and network theory analyses in cathepsin K (catK), whose allosteric sites are well defined. To obtain a quantitative measure of the efficiency of communication, we designed a new protocol that enables the comparison between properties related to ensembles of communication paths obtained from different sites. Further, we applied our strategy to evaluate the allosteric potential of different catK cavities not yet considered for drug design. Our predictions of the allosteric potential based on intramolecular communication correlate well with previous catK experimental and theoretical data. We also discuss the possibility of applying our approach to other proteins from the same family.

Understanding the structural details of APOBEC3-DNA interactions using graph-based representations

Human APOBEC3 (A3; apolipoprotein B mRNA editing catalytic polypeptide-like 3) is a family of seven enzymes involved in generating mutations in nascent reverse transcripts of many retroviruses, as well as the human genome in a range of cancer types. The structural details of the interaction between A3 proteins and DNA molecules are only available for a few family members. Here we use homology modelling techniques to address the difference in structural coverage of human A3 enzymes interacting with different DNA substrates. A3-DNA interfaces are represented as residue networks ("graphs"), based on which features at these interfaces are compared and quantified. We demonstrate that graph-based representations are effective in highlighting structural features of A3-DNA interfaces. By large-scale in silico mutagenesis of the bound DNA chain, we predicted the preference of substrate DNA sequence for multiple A3 domains. These data suggested that computational modelling approaches could contribute in the exploration of the structural basis for sequence specificity in A3 substrate selection, and demonstrated the utility of graph-based approaches in evaluating a large number of structural models generated in silico.

•

APOBEC3(A3)-DNA structures have been resolved with modified deaminase domains.

•

Structural modelling of interaction between wild-type A3 domains and DNA substrates.

•

Graph-based representations reveal structural differences across A3-DNA interfaces.

•

Using in silico mutagenesis we compared substrate preference of multiple A3 domains.

•

Graph-based approaches can efficiently compare a large number of structural models.

ProSNEx: a web-based application for exploration and analysis of protein structures using network formalism

ProSNEx (Protein Structure Network Explorer) is a web service for construction and analysis of Protein Structure Networks (PSNs) alongside amino acid flexibility, sequence conservation and annotation features. ProSNEx constructs a PSN by adding nodes to represent residues and edges between these nodes using user-specified interaction distance cutoffs for either carbon-alpha, carbon-beta or atom-pair contact networks. Different types of weighted networks can also be constructed by using either (i) the residue-residue interaction energies in the format returned by gRINN, resulting in a Protein Energy Network (PEN); (ii) the dynamical cross correlations from a coarse-grained Normal Mode Analysis (NMA) of the protein structure; (iii) interaction strength. Upon construction of the network, common network metrics (such as node centralities) as well as shortest paths between nodes and k-cliques are calculated. Moreover, additional features of each residue in the form of conservation scores and mutation/natural variant information are included in the analysis. By this way, tool offers an enhanced and direct comparison of network-based residue metrics with other types of biological information. ProSNEx is free and open to all users without login requirement at http://prosnex-tool.com.

Long-range Regulation of Partially Folded Amyloidogenic Peptides

Neurodegeneration involves abnormal aggregation of intrinsically disordered amyloidogenic peptides (IDPs), usually mediated by hydrophobic protein-protein interactions. There is mounting evidence that formation of α-helical intermediates is an early event during self-assembly of amyloid-β42 (Aβ42) and α-synuclein (αS) IDPs in Alzheimer’s and Parkinson’s disease pathogenesis, respectively. However, the driving force behind on-pathway molecular assembly of partially folded helical monomers into helical oligomers assembly remains unknown. Here, we employ extensive molecular dynamics simulations to sample the helical conformational sub-spaces of monomeric peptides of both Aβ42 and αS. Our computed free energies, population shifts, and dynamic cross-correlation network analyses reveal a common feature of long-range intra-peptide modulation of partial helical folds of the amyloidogenic central hydrophobic domains via concerted coupling with their charged terminal tails (N-terminus of Aβ42 and C-terminus of αS). The absence of such inter-domain fluctuations in both fully helical and completely unfolded (disordered) states suggests that long-range coupling regulates the dynamicity of partially folded helices, in both Aβ42 and αS peptides. The inter-domain coupling suggests a form of intra-molecular allosteric regulation of the aggregation trigger in partially folded helical monomers. This approach could be applied to study the broad range of amyloidogenic peptides, which could provide a new path to curbing pathogenic aggregation of partially folded conformers into oligomers, by inhibition of sites far from the hydrophobic core.

Non-symmetrical structural behavior of a symmetric protein: the case of homo-trimeric TRAF2 (tumor necrosis factor-receptor associated factor 2)

The oligomeric state of TRAF2 (tumor necrosis factor-receptor associated factor 2), a TNF (tumor necrosis factor) receptor-associated factor, is crucial for membrane binding and probably plays a fundamental role in regulating the protein function in vivo. In this study we have combined molecular dynamics with the protein contact network approach to characterize the interaction of the three identical subunits of TRAF2. The average structure obtained after a 225 ns simulation reveals that two clusters of different size are formed, one of which includes almost completely two subunits, while the third monomer appears to be more independent. This picture is also confirmed by the estimated average number of inter-subunit contacts and by the comparison of side chains mobility in each monomer. The analysis of equilibrium pressure-induced dissociation measurements supports such findings, indicating that the dimeric-monomeric (2 + 1) might be prevalent with respect to the trimeric configuration, especially in the case of more diluted samples. These findings suggest that the formation of monomeric species, which is crucial for the formation of intra-luminal vesicles, might depend on preferential asymmetric interactions among the three subunits.Communicated by Ramaswamy H. Sarma.

Characterizing Protein-Protein Interactions with the Fragment Molecular Orbital Method

Proteins are vital components of living systems, serving as building blocks, molecular machines, enzymes, receptors, ion channels, sensors, and transporters. Protein-protein interactions (PPIs) are a key part of their function. There are more than 645,000 reported disease-relevant PPIs in the human interactome, but drugs have been developed for only 2% of these targets. The advances in PPI-focused drug discovery are highly dependent on the availability of structural data and accurate computational tools for analysis of this data. Quantum mechanical approaches are often too expensive computationally, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed and the ability to reveal key interactions that would otherwise be hard to detect. FMO provides essential information for PPI drug discovery, namely, identification of key interactions formed between residues of two proteins, including their strength (in kcal/mol) and their chemical nature (electrostatic or hydrophobic). In this chapter, we have demonstrated how three different FMO-based approaches (pair interaction energy analysis (PIE analysis), subsystem analysis (SA) and analysis of protein residue networks (PRNs)) have been applied to study PPI in three protein-protein complexes.

Alterations of the interactome of Bcl-2 proteins in breast cancer at the transcriptional, mutational and structural level

Apoptosis is an essential defensive mechanism against tumorigenesis. Proteins of the B-cell lymphoma-2 (Bcl-2) family regulate programmed cell death by the mitochondrial apoptosis pathway. In response to intracellular stress, the apoptotic balance is governed by interactions of three distinct subgroups of proteins; the activator/sensitizer BH3 (Bcl-2 homology 3)-only proteins, the pro-survival, and the pro-apoptotic executioner proteins. Changes in expression levels, stability, and functional impairment of pro-survival proteins can lead to an imbalance in tissue homeostasis. Their overexpression or hyperactivation can result in oncogenic effects. Pro-survival Bcl-2 family members carry out their function by binding the BH3 short linear motif of pro-apoptotic proteins in a modular way, creating a complex network of protein-protein interactions. Their dysfunction enables cancer cells to evade cell death. The critical role of Bcl-2 proteins in homeostasis and tumorigenesis, coupled with mounting insight in their structural properties, make them therapeutic targets of interest. A better understanding of gene expression, mutational profile, and molecular mechanisms of pro-survival Bcl-2 proteins in different cancer types, could help to clarify their role in cancer development and may guide advancement in drug discovery. Here, we shed light on the pro-survival Bcl-2 proteins in breast cancer using different bioinformatic approaches, linking -omics with structural data. We analyzed the changes in the expression of the Bcl-2 proteins and their BH3-containing interactors in breast cancer samples. We then studied, at the structural level, a selection of interactions, accounting for effects induced by mutations found in the breast cancer samples. We find two complexes between the up-regulated Bcl2A1 and two down-regulated BH3-only candidates (i.e., Hrk and Nr4a1) as targets associated with reduced apoptosis in breast cancer samples for future experimental validation. Furthermore, we predict L99R, M75R as damaging mutations altering protein stability, and Y120C as a possible allosteric mutation from an exposed surface to the BH3-binding site.

Structural topology defines protective CD8+ T cell epitopes in the HIV proteome

Mutationally constrained epitopes of variable pathogens represent promising targets for vaccine design but are not reliably identified by sequence conservation. In this study, we employed structure-based network analysis, which applies network theory to HIV protein structure data to quantitate the topological importance of individual amino acid residues. Mutation of residues at important network positions disproportionately impaired viral replication and occurred with high frequency in epitopes presented by protective human leukocyte antigen (HLA) class I alleles. Moreover, CD8⁺ T cell targeting of highly networked epitopes distinguished individuals who naturally control HIV, even in the absence of protective HLA alleles. This approach thereby provides a mechanistic basis for immune control and a means to identify CD8⁺ T cell epitopes of topological importance for rational immunogen design, including a T cell-based HIV vaccine.

Go-Kit: A Tool To Enable Energy Landscape Exploration of Proteins

Coarse-grained Go̅-like models, based on the principle of minimal frustration, provide valuable insight into fundamental questions in the field of protein folding and dynamics. In conjunction with commonly used molecular dynamics (MD) simulations, energy landscape exploration methods like discrete path sampling (DPS) with Go̅-like models can provide quantitative details of the thermodynamics and kinetics of proteins. Here we present Go-kit, a software that facilitates the setup of MD and DPS simulations of several flavors of Go̅-like models. Go-kit is designed for use with MD (GROMACS) and DPS (PATHSAMPLE) simulation engines that are open source. The Go-kit code is written in python2.7 and is also open source. A case study for the ribosomal protein S6 is discussed to illustrate the utility of the software, which is available at https://github.com/gokit1/gokit .

The protonation state of an evolutionarily conserved histidine modulates domainswapping stability of FoxP1

Forkhead box P (FoxP) proteins are members of the versatile Fox transcription factors, which control the timing and expression of multiple genes for eukaryotic cell homeostasis. Compared to other Fox proteins, they can form domain-swapped dimers through their DNA-binding –forkhead– domains, enabling spatial reorganization of distant chromosome elements by tethering two DNA molecules together. Yet, domain swapping stability and DNA binding affinity varies between different FoxP proteins. Experimental evidence suggests that the protonation state of a histidine residue conserved in all Fox proteins is responsible for pH-dependent modulation of these interactions. Here, we explore the consequences of the protonation state of another histidine (H59), only conserved within FoxM/O/P subfamilies, on folding and dimerization of the forkhead domain of human FoxP1. Dimer dissociation kinetics and equilibrium unfolding experiments demonstrate that protonation of H59 leads to destabilization of the domain-swapped dimer due to an increase in free energy difference between the monomeric and transition states. This pH–dependence is abolished when H59 is mutated to alanine. Furthermore, anisotropy measurements and molecular dynamics evidence that H59 has a direct impact in the local stability of helix H3. Altogether, our results highlight the relevance of H59 in domain swapping and folding stability of FoxP1.

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.

Analyzing Biomolecular Ensembles

Several techniques are available to generate conformational ensembles of proteins and other biomolecules either experimentally or computationally. These methods produce a large amount of data that need to be analyzed to identify structure-dynamics-function relationship. In this chapter, we will cover different tools to unveil the information hidden in conformational ensemble data and to guide toward the rationalization of the data. We included routinely used approaches such as dimensionality reduction, as well as new methods inspired by high-order statistics and graph theory.

Protein Residue Networks from Energetic and Geometric Data: Are They Identical?

Protein residue networks (PRN) from energetic and geometric data are probably not identical. PRNs constructed from ab initio pair interaction energies are analyzed for the first time and compared to PRN based on center of mass separation. We use modern, previously unused algorithms such as global and local efficiencies to quantitatively confirm that both types of PRNs do exhibit small-world character. The main novelty finding is that interaction energy-based PRNs preserve small-world character even when clustered. A node hierarchy independent of the cutoff energy used for the edge creation is characteristic for them. Efficiency centrality identifies hubs responsible for such behavior. The interaction energy-based PRNs seem to comply with the scale-free network model with respect to efficiency centrality distribution as opposed to distance based PRNs. Community detection is introduced into protein network research as an extension beyond cluster analysis to study tertiary and quaternary structures.