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

NMR tools to detect protein allostery

Allostery is a fundamental mechanism of cellular homeostasis by intra-protein communication between distinct functional sites. It is an internal process of proteins to steer interactions not only with each other but also with other biomolecules such as ligands, lipids, and nucleic acids. In addition, allosteric regulation is particularly important in enzymatic activities. A major challenge in structural and molecular biology today is unraveling allosteric sites in proteins, to elucidate the detailed mechanism of allostery and the development of allosteric drugs. Here we summarize the recently developed tools and approaches which enable the elucidation of regulatory hotspots and correlated motion in biomolecules, focusing primarily on solution-state nuclear magnetic resonance spectroscopy (NMR). These tools open an avenue towards a rational understanding of the mechanism of allostery and provide essential information for the design of allosteric drugs.

Atomistic mechanism of coupling between cytosolic sensor domain and selectivity filter in TREK K2P channels

The two-pore domain potassium (K2P) channels TREK-1 and TREK-2 link neuronal excitability to a variety of stimuli including mechanical force, lipids, temperature and phosphorylation. This regulation involves the C-terminus as a polymodal stimulus sensor and the selectivity filter (SF) as channel gate. Using crystallographic up- and down-state structures of TREK-2 as a template for full atomistic molecular dynamics (MD) simulations, we reveal that the SF in down-state undergoes inactivation via conformational changes, while the up-state structure maintains a stable and conductive SF. This suggests an atomistic mechanism for the low channel activity previously assigned to the down state, but not evident from the crystal structure. Furthermore, experimentally by using (de-)phosphorylation mimics and chemically attaching lipid tethers to the proximal C-terminus (pCt), we confirm the hypothesis that moving the pCt towards the membrane induces the up-state. Based on MD simulations, we propose two gating pathways by which movement of the pCt controls the stability (i.e., conductivity) of the filter gate. Together, these findings provide atomistic insights into the SF gating mechanism and the physiological regulation of TREK channels by phosphorylation.

Key interaction networks: Identifying evolutionarily conserved non-covalent interaction networks across protein families

Protein structure (and thus function) is dictated by non-covalent interaction networks. These can be highly evolutionarily conserved across protein families, the members of which can diverge in sequence and evolutionary history. Here we present KIN, a tool to identify and analyze conserved non-covalent interaction networks across evolutionarily related groups of proteins. KIN is available for download under a GNU General Public License, version 2, from https://www.github.com/kamerlinlab/KIN. KIN can operate on experimentally determined structures, predicted structures, or molecular dynamics trajectories, providing insight into both conserved and missing interactions across evolutionarily related proteins. This provides useful insight both into protein evolution, as well as a tool that can be exploited for protein engineering efforts. As a showcase system, we demonstrate applications of this tool to understanding the evolutionary-relevant conserved interaction networks across the class A β-lactamases.

Mechanism of Mutation-Induced Effects on the Catalytic Function of TEV Protease: A Molecular Dynamics Study

Tobacco etch virus protease (TEVp) is wildly exploited for various biotechnological applications. These applications take advantage of TEVp's ability to cleave specific substrate sequences to study protein function and interactions. A major limitation of this enzyme is its relatively slow catalytic rate. In this study, MD simulations were conducted on TEV enzymes and known highly active mutants (eTEV and uTEV3) to explore the relationship between mutation, conformation, and catalytic function. The results suggest that mutations distant from the active site can influence the substrate-binding pocket through interaction networks. MD analysis of eTEV demonstrates that, by stabilizing the orientation of the substrate at the catalytic site, mutations that appropriately enlarge the substrate-binding pocket will be beneficial for Kcat, enhancing the catalytic efficiency of the enzyme. On the contrary, mutations in uTEV3 reduced the flexibility of the active pocket and increased the hydrogen bonding between the substrate and enzyme, resulting in higher affinity. At the same time, the MD simulation demonstrates that mutations outside of the active site residues could affect the dynamic movement of the binding pocket by altering residue networks and communication pathways, thereby having a profound impact on reactivity. These findings not only provide a molecular mechanistic explanation for the excellent mutants, but also serve as a guiding framework for rational computational design.

Allosteric Activation of RhoA Complexed with p115-RhoGEF Deciphered by Conformational Dynamics

The Ras homologue family member A (RhoA) is a member of the Rho family, a subgroup of the Ras superfamily. RhoA interacts with the 115 kDa guanine nucleotide exchange factor (p115-RhoGEF), which assists in activation and binding with downstream effectors. Here, we use molecular dynamics (MD) simulations and essential dynamics analysis of the inactive RhoA-GDP and active RhoA-GTP, when bound to p115-RhoGEF to decipher the mechanism of RhoA activation at the structural level. We observe that inactive RhoA-GDP maintains its position near the catalytic site on the Dbl homology (DH) domain of p115-RhoGEF through the interaction of its Switch I region with the DH domain. We further show that the active RhoA-GTP is engaged in more interactions with the p115-RhoGEF membrane-bound Pleckstrin homology (PH) domain as compared to RhoA-GDP. We hypothesize that the role of the interactions between the active RhoA-GTP and the PH domain is to help release it from the DH domain upon activation. Our results support this premise, and our simulations uncover the beginning of this process and provide structural details. They also point to allosteric communication pathways that take part in RhoA activation to promote and strengthen the interaction between the active RhoA-GTP and the PH domain. Allosteric regulation also occurs among other members of the Rho superfamily. Collectively, we suggest that in the activation process, the role of the RhoA-GTP interaction with the PH domain is to release RhoA-GTP from the DH domain after activation, making it available to downstream effectors.

Allosteric communication and signal transduction in proteins

Allostery is one of the cornerstones of biological function, as it plays a fundamental role in regulating protein activity. The modelling of allostery has gradually moved from a conformation-based framework, linked to structural changes, to dynamics-based allostery, whereby the effects of ligand binding propagate via signal transduction from the allosteric site to other regions of the protein via inter-residue interactions. Characterising such allosteric signalling pathways, which do not necessarily lead to conformational changes, has been pursued experimentally and complemented by computational analysis of protein networks to detect subtle dynamic propagation paths. Considering allostery from the perspective of signal transduction broadens the understanding of allosteric mechanisms, underscores the importance of protein topology, and can provide insights into allosteric drug design.

Conformational plasticity and allosteric communication networks explain Shelterin protein TPP1 binding to human telomerase

The Shelterin complex protein TPP1 interacts with human telomerase (TERT) by means of the TEL-patch region, controlling telomere homeostasis. Aberrations in the TPP1-TERT heterodimer formation might lead to short telomeres and severe diseases like dyskeratosis congenita and Hoyeraal-Hreidarsson syndrome. In the present study, we provide a thorough characterization of the structural properties of the TPP1's OB-domain by combining data coming from microsecond-long molecular dynamics calculations, time-series analyses, and graph-based networks. Our results show that the TEL-patch conformational freedom is influenced by a network of long-range amino acid communications that together determine the proper TPP1-TERT binding. Furthermore, we reveal that in TPP1 pathological variants Glu169Δ, Lys170Δ and Leu95Gln, the TEL-patch plasticity is reduced, affecting the correct binding to TERT and, in turn, telomere processivity, which eventually leads to accelerated aging of affected cells. Our study provides a structural basis for the design of TPP1-targeting ligands with therapeutic potential against cancer and telomeropathies.

Slc11 Synapomorphy: A Conserved 3D Framework Articulating Carrier Conformation Switch

Transmembrane carriers of the Slc11 family catalyze proton (H+)-dependent uptake of divalent metal ions (Me2+) such as manganese and iron-vital elements coveted during infection. The Slc11 mechanism of high-affinity Me2+ cell import is selective and conserved between prokaryotic (MntH) and eukaryotic (Nramp) homologs, though processes coupling the use of the proton motive force to Me2+ uptake evolved repeatedly. Adding bacterial piracy of Nramp genes spread in distinct environmental niches suggests selective gain of function that may benefit opportunistic pathogens. To better understand Slc11 evolution, Alphafold (AF2)/Colabfold (CF) 3D predictions for bacterial sequences from sister clades of eukaryotic descent (MCb and MCg) were compared using both native and mutant templates. AF2/CF model an array of native MCb intermediates spanning the transition from outwardly open (OO) to inwardly open (IO) carriers. In silico mutagenesis targeting (i) a set of (evolutionarily coupled) sites that may define Slc11 function (putative synapomorphy) and (ii) residues from networked communities evolving during MCb transition indicates that Slc11 synapomorphy primarily instructs a Me2+-selective conformation switch which unlocks carrier inner gate and contributes to Me2+ binding site occlusion and outer gate locking. Inner gate opening apparently proceeds from interaction between transmembrane helix (h) h5, h8 and h1a. MCg1 xenologs revealed marked differences in carrier shape and plasticity, owing partly to an altered intramolecular H+ network. Yet, targeting Slc11 synapomorphy also converted MCg1 IO models to an OO state, apparently mobilizing the same residues to control gates. But MCg1 response to mutagenesis differed, with extensive divergence within this clade correlating with MCb-like modeling properties. Notably, MCg1 divergent epistasis marks the emergence of the genus Bordetella-Achromobacter. Slc11 synapomorphy localizes to the 3D areas that deviate least among MCb and MCg1 models (either IO or OO) implying that it constitutes a 3D network of residues articulating a Me2+-selective carrier conformation switch which is maintained in fast-evolving clades at the cost of divergent epistatic interactions impacting carrier shape and dynamics.

Intrinsic Dynamics of the ClpXP Proteolytic Machine Using Elastic Network Models

ClpXP complex is an ATP-dependent mitochondrial matrix protease that binds, unfolds, translocates, and subsequently degrades specific protein substrates. Its mechanisms of operation are still being debated, and several have been proposed, including the sequential translocation of two residues (SC/2R), six residues (SC/6R), and even long-pass probabilistic models. Therefore, it has been suggested to employ biophysical-computational approaches that can determine the kinetics and thermodynamics of the translocation. In this sense, and based on the apparent inconsistency between structural and functional studies, we propose to apply biophysical approaches based on elastic network models (ENM) to study the intrinsic dynamics of the theoretically most probable hydrolysis mechanism. The proposed models ENM suggest that the ClpP region is decisive for the stabilization of the ClpXP complex, contributing to the flexibility of the residues adjacent to the pore, favoring the increase in pore size and, therefore, with the energy of interaction of its residues with a larger portion of the substrate. It is predicted that the complex may undergo a stable configurational change once assembled and that the deformability of the system once assembled is oriented, to increase the rigidity of the domains of each region (ClpP and ClpX) and to gain flexibility of the pore. Our predictions could suggest under the conditions of this study the mechanism of the interaction of the system, of which the substrate passes through the unfolding of the pore in parallel with a folding of the bottleneck. The variations in the distance calculated by molecular dynamics could allow the passage of a substrate with a size equivalent to ∼3 residues. The theoretical behavior of the pore and the stability and energy of binding to the substrate based on ENM models suggest that in this system, there are thermodynamic, structural, and configurational conditions that allow a possible translocation mechanism that is not strictly sequential.

GSNOR deficiency promotes tumor growth via FAK1 S-nitrosylation

Nitric oxide (NO) production in the tumor microenvironment is a common element in cancer. S-nitrosylation, the post-translational modification of cysteines by NO, is emerging as a key transduction mechanism sustaining tumorigenesis. However, most oncoproteins that are regulated by S-nitrosylation are still unknown. Here we show that S-nitrosoglutathione reductase (GSNOR), the enzyme that deactivates S-nitrosylation, is hypo-expressed in several human malignancies. Using multiple tumor models, we demonstrate that GSNOR deficiency induces S-nitrosylation of focal adhesion kinase 1 (FAK1) at C658. This event enhances FAK1 autophosphorylation and sustains tumorigenicity by providing cancer cells with the ability to survive in suspension (evade anoikis). In line with these results, GSNOR-deficient tumor models are highly susceptible to treatment with FAK1 inhibitors. Altogether, our findings advance our understanding of the oncogenic role of S-nitrosylation, define GSNOR as a tumor suppressor, and point to GSNOR hypo-expression as a therapeutically exploitable vulnerability in cancer.

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

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.

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.

Structural communication between the GTPase Sec4p and its activator Sec2p: Determinants of GEF activity and early deformations to nucleotide release

Ras GTPases are molecular switches that cycle between OFF and ON states depending on the bound nucleotide (i.e. GDP-bound and GTP-bound, respectively).

The Rab GTPase, Sec4p, plays regulatory roles in multiple steps of intracellular vesicle trafficking. Nucleotide release is catalyzed by the Guanine Nucleotide Exchange Factor (GEF) Sec2p.

Here, the integration of structural information with molecular dynamics (MD) simulations addressed a number of questions concerning the intrinsic and stimulated dynamics of Sec2p and Sec4p as well as the chain of structural deformations leading to GEF-assisted activation of the Rab GTPase.

Sec2p holds an intrinsic ability to adopt the conformation found in the crystallographic complexes with Sec4p, thus suggesting that the latter selects and shifts the conformational equilibrium towards a pre-existing bound-like conformation of Sec2p.

The anchoring of Sec4p to a suitable conformation of Sec2p favors the Sec2p-assisted pulling on itself of the α1/switch 1 (SWI) loop and of SWI, which loose any contact with GDP. Those deformations of Sec4p would occur earlier. Formation of the final Sec2p-Sec4p hydrophobic interface, accomplishes later. Disruption of the nucleotide cage would cause firstly loss of interactions with the guanine ring and secondly loss of interactions with the phosphates.

The ease in sampling the energy landscape and adopting a bound-like conformation likely favors the catalyzing ability of GEFs for Ras GTPases.