Binding pocket dynamics along the recovery stroke of human ß-cardiac myosin.

The druggability of small-molecule binding sites can be significantly affected by protein motions and conformational changes. Ligand binding, protein dynamics and protein function have been shown to be closely interconnected in myosins. The breakthrough discovery of omecamtiv mecarbil (OM) has led to an increased interest in small molecules that can target myosin and modulate its function for therapeutic purposes (myosin modulators). In this work, we use a combination of computational methods, including steered molecular dynamics, umbrella sampling and binding pocket tracking tools, to follow the evolution of the OM binding site during the recovery stroke transition of human ß-cardiac myosin. We found that steering two internal coordinates of the motor domain can recapture the main features of the transition and in particular the rearrangements of the binding site, which shows significant changes in size, shape and composition. Possible intermediate conformations were also identified, in remarkable agreement with experimental findings. The differences in the binding site properties observed along the transition can be exploited for the future development of conformation-selective myosin modulators.

Comparative study of binding pocket structure and dynamics in cardiac and skeletal myosin.

The development of small molecule myosin modulators has seen an increased effort in recent years due to their possible use in the treatment of cardiac and skeletal myopathies. Omecamtiv mecarbil (OM) is the first-in-class cardiac myotrope and the first to enter clinical trials. Its selectivity toward slow/beta-cardiac myosin lies at the heart of its function; however, little is known about the underlying reasons for selectivity to this isoform as opposed to other closely related ones such as fast-type skeletal myosins. In this work, we compared the structure and dynamics of the OM binding site in cardiac and in fasttype IIa skeletal myosin to identify possible reasons for OM selectivity. We found that the different shape, size, and composition of the binding pocket in skeletal myosin directly affects the binding mode and related affinity of OM, which is potentially a result of weaker interactions and less optimal molecular recognition. Moreover, we identified a side pocket adjacent to the OM binding site that shows increased accessibility in skeletal myosin compared with the cardiac isoform. These findings could pave the way to the development of skeletal-selective compounds that can target this region of the protein and potentially be used to treat congenital myopathies where muscle weakness is related to myosin loss of function.

Structure and functional analysis of the Legionella pneumophila chitinase ChiA reveals a novel mechanism of metal-dependent mucin degradation.

Chitinases are important enzymes that contribute to the generation of carbon and nitrogen from chitin, a long chain polymer of N-acetylglucosamine that is abundant in insects, fungi, invertebrates and fish. Although mammals do not produce chitin, chitinases have been identified in bacteria that are key virulence factors in severe respiratory, gastrointestinal and urinary diseases. However, it is unclear how these enzymes are able to carry out this dual function. Legionella pneumophila is the causative agent of Legionnaires' disease, an often-fatal pneumonia and its chitinase ChiA is essential for the survival of L. pneumophila in the lung. Here we report the first atomic resolution insight into the pathogenic mechanism of a bacterial chitinase. We derive an experimental model of intact ChiA and show how its N-terminal region targets ChiA to the bacterial surface after its secretion. We provide the first evidence that L. pneumophila can bind mucins on its surface, but this is not dependent on ChiA. This demonstrates that additional peripheral mucin binding proteins are also expressed in L. pneumophila. We also show that the ChiA C-terminal chitinase domain has novel Zn2+-dependent peptidase activity against mammalian mucin-like proteins, namely MUC5AC and the C1-esterase inhibitor, and that ChiA promotes bacterial penetration of mucin gels. Our findings suggest that ChiA can facilitate passage of L. pneumophila through the alveolar mucosa, can modulate the host complement system and that ChiA may be a promising target for vaccine development.

A pH-regulated quality control cycle for surveillance of secretory protein assembly.

To warrant the quality of the secretory proteome, stringent control systems operate at the endoplasmic reticulum (ER)-Golgi interface, preventing the release of nonnative products. Incompletely assembled oligomeric proteins that are deemed correctly folded must rely on additional quality control mechanisms dedicated to proper assembly. Here we unveil how ERp44 cycles between cisGolgi and ER in a pH-regulated manner, patrolling assembly of disulfide-linked oligomers such as IgM and adiponectin. At neutral, ER-equivalent pH, the ERp44 carboxy-terminal tail occludes the substrate-binding site. At the lower pH of the cisGolgi, conformational rearrangements of this peptide, likely involving protonation of ERp44's active cysteine, simultaneously unmask the substrate binding site and -RDEL motif, allowing capture of orphan secretory protein subunits and ER retrieval via KDEL receptors. The ERp44 assembly control cycle couples secretion fidelity and efficiency downstream of the calnexin/calreticulin and BiP-dependent quality control cycles.

Heart Failure Drug Modifies the Intrinsic Dynamics of the Pre-Power Stroke State of Cardiac Myosin.

Omecamtiv mecarbil (OM), currently investigated for the treatment of heart failure, is the first example of a new class of drugs (cardiac myotropes) that can modify muscle contractility by directly targeting sarcomeric proteins. Using atomistic molecular dynamics simulations, we show that the binding of OM to the pre-power stroke state of cardiac myosin inhibits the functional motions of the protein and potentially affects Pi release from the nucleotide binding site. We also show that the changes in myosin ATPase activity induced by a set of OM analogues can be predicted from their relative affinity to the pre-power stroke state compared to the near rigor one, indicating that conformational selectivity plays an important role in determining the activity of these compounds.

In silico identification of rescue sites by double force scanning.

Motivation: A deleterious amino acid change in a protein can be compensated by a second-site rescue mutation. These compensatory mechanisms can be mimicked by drugs. In particular, the location of rescue mutations can be used to identify protein regions that can be targeted by small molecules to reactivate a damaged mutant. Results: We present the first general computational method to detect rescue sites. By mimicking the effect of mutations through the application of forces, the double force scanning (DFS) method identifies the second-site residues that make the protein structure most resilient to the effect of pathogenic mutations. We tested DFS predictions against two datasets containing experimentally validated and putative evolutionary-related rescue sites. A remarkably good agreement was found between predictions and experimental data. Indeed, almost half of the rescue sites in p53 was correctly predicted by DFS, with 65% of remaining sites in contact with DFS predictions. Similar results were found for other proteins in the evolutionary dataset. Availability and implementation: The DFS code is available under GPL at https://fornililab.github.io/dfs/. Supplementary information: Supplementary data are available at Bioinformatics online.

Allosteric modulation of cardiac myosin dynamics by omecamtiv mecarbil.

New promising avenues for the pharmacological treatment of skeletal and heart muscle diseases rely on direct sarcomeric modulators, which are molecules that can directly bind to sarcomeric proteins and either inhibit or enhance their activity. A recent breakthrough has been the discovery of the myosin activator omecamtiv mecarbil (OM), which has been shown to increase the power output of the cardiac muscle and is currently in clinical trials for the treatment of heart failure. While the overall effect of OM on the mechano-chemical cycle of myosin is to increase the fraction of myosin molecules in the sarcomere that are strongly bound to actin, the molecular basis of its action is still not completely clear. We present here a Molecular Dynamics study of the motor domain of human cardiac myosin bound to OM, where the effects of the drug on the dynamical properties of the protein are investigated for the first time with atomistic resolution. We found that OM has a double effect on myosin dynamics, inducing a) an increased coupling of the motions of the converter and lever arm subdomains to the rest of the protein and b) a rewiring of the network of dynamic correlations, which produces preferential communication pathways between the OM binding site and distant functional regions. The location of the residues responsible for these effects suggests possible strategies for the future development of improved drugs and the targeting of specific cardiomyopathy-related mutations.

GSATools: analysis of allosteric communication and functional local motions using a structural alphabet.

MOTIVATION: GSATools is a free software package to analyze conformational ensembles and to detect functional motions in proteins by means of a structural alphabet. The software integrates with the widely used GROMACS simulation package and can generate a range of graphical outputs. Three applications can be supported: (i) investigation of the conformational variability of local structures; (ii) detection of allosteric communication; and (iii) identification of local regions that are critical for global functional motions. These analyses provide insights into the dynamics of proteins and allow for targeted design of functional mutants in theoretical and experimental studies. AVAILABILITY: The C source code of the GSATools, along with a set of pre-compiled binaries, is freely available under GNU General Public License from http://mathbio.nimr.mrc.ac.uk/wiki/GSATools.

Elucidating the Activation Mechanism of the Proton-sensing GPR68 Receptor.

GPR68 is a proton-sensing G-protein Coupled Receptor (GPCR) involved in a variety of physiological processes and disorders including neoplastic pathologies. While GPR68 and few other GPCRs have been shown to be activated by a decrease in the extracellular pH, the molecular mechanism of their activation remains largely unknown. In this work, we used a combined computational and in vitro approach to provide new insight into the activation mechanism of the receptor. Molecular Dynamics simulations of GPR68 were used to model the changes in residue interactions and motions triggered by pH. Global and local rearrangements consistent with partial activation were observed upon protonation of the inactive state. Selected extracellular histidine and transmembrane acidic residues were found to have significantly upshifted pKa values during the simulations, consistently with their previously hypothesised role in activation through changes in protonation state. Moreover, a novel pairing between histidine and acidic residues in the extracellular region was highlighted by both sequence analyses and simulation data and tested through site-directed mutagenesis. At last, we identified a previously unknown hydrophobic lock in the extracellular region that might stabilise the inactive conformation and regulate the transition to the active state.

N4BP1 functions as a dimerization-dependent linear ubiquitin reader which regulates TNF signalling.

Signalling through TNFR1 modulates proinflammatory gene transcription and programmed cell death, and its impairment causes autoimmune diseases and cancer. NEDD4-binding protein 1 (N4BP1) is a critical suppressor of proinflammatory cytokine production that acts as a regulator of innate immune signalling and inflammation. However, our current understanding about the molecular properties that enable N4BP1 to exert its suppressive potential remain limited. Here, we show that N4BP1 is a novel linear ubiquitin reader that negatively regulates NFκB signalling by its unique dimerization-dependent ubiquitin-binding module that we named LUBIN. Dimeric N4BP1 strategically positions two non-selective ubiquitin-binding domains to ensure preferential recognition of linear ubiquitin. Under proinflammatory conditions, N4BP1 is recruited to the nascent TNFR1 signalling complex, where it regulates duration of proinflammatory signalling in LUBIN-dependent manner. N4BP1 deficiency accelerates TNFα-induced cell death by increasing complex II assembly. Under proapoptotic conditions, caspase-8 mediates proteolytic processing of N4BP1, resulting in rapid degradation of N4BP1 by the 26 S proteasome, and acceleration of apoptosis. In summary, our findings demonstrate that N4BP1 dimerization creates a novel type of ubiquitin reader that selectively recognises linear ubiquitin which enables the timely and coordinated regulation of TNFR1-mediated inflammation and cell death.

The Legionella collagen-like protein employs a distinct binding mechanism for the recognition of host glycosaminoglycans.

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual trimer arrangement with a positively charged external surface and negatively charged solvent exposed internal cavity. Through molecular dynamics simulations, we show how the glycosaminoglycan chondroitin-4-sulphate associates with the Lcl-CTD surface via distinct binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate-binding mechanism.

Protein-protein interaction networks studies and importance of 3D structure knowledge.

Protein-protein interaction networks (PPINs) are a powerful tool to study biological processes in living cells. In this review, we present the progress of PPIN studies from abstract to more detailed representations. We will focus on 3D interactome networks, which offer detailed information at the atomic level. This information can be exploited in understanding not only the underlying cellular mechanisms, but also how human variants and disease-causing mutations affect protein functions and complexes' stability. Recent studies have used structural information on PPINs to also understand the molecular mechanisms of binding partner selection. We will address the challenges in generating 3D PPINs due to the restricted number of solved protein structures. Finally, some of the current use of 3D PPINs will be discussed, highlighting their contribution to the studies in genotype-phenotype relationships and in the optimization of targeted studies to design novel chemical compounds for medical treatments.

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics.

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.

Abnormal myosin post-translational modifications and ATP turnover time associated with human congenital myopathy-related RYR1 mutations.

AIM: Conditions related to mutations in the gene encoding the skeletal muscle ryanodine receptor 1 (RYR1) are genetic muscle disorders and include congenital myopathies with permanent weakness, as well as episodic phenotypes such as rhabdomyolysis/myalgia. Although RYR1 dysfunction is the primary mechanism in RYR1-related disorders, other downstream pathogenic events are less well understood and may include a secondary remodeling of major contractile proteins. Hence, in the present study, we aimed to investigate whether congenital myopathy-related RYR1 mutations alter the regulation of the most abundant contractile protein, myosin. METHODS: We used skeletal muscle tissues from five patients with RYR1-related congenital myopathy and compared those with five controls and five patients with RYR1-related rhabdomyolysis/myalgia. We then defined post-translational modifications on myosin heavy chains (MyHCs) using LC/MS. In parallel, we determined myosin relaxed states using Mant-ATP chase experiments and performed molecular dynamics (MD) simulations. RESULTS: LC/MS revealed two additional phosphorylations (Thr1309-P and Ser1362-P) and one acetylation (Lys1410-Ac) on the ß/slow MyHC of patients with congenital myopathy. This method also identified six acetylations that were lacking on MyHC type IIa of these patients (Lys35-Ac, Lys663-Ac, Lys763-Ac, Lys1171-Ac, Lys1360-Ac, and Lys1733-Ac). MD simulations suggest that modifying myosin Ser1362 impacts the protein structure and dynamics. Finally, Mant-ATP chase experiments showed a faster ATP turnover time of myosin heads in the disordered-relaxed conformation. CONCLUSIONS: Altogether, our results suggest that RYR1 mutations have secondary negative consequences on myosin structure and function, likely contributing to the congenital myopathic phenotype.

The Legionella collagen-like protein employs a unique binding mechanism for the recognition of host glycosaminoglycans.

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans (GAGs) on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual dynamic trimer arrangement with a positively charged external surface and a negatively charged solvent exposed internal cavity. Through Molecular Dynamics (MD) simulations, we show how the GAG chondroitin-4-sulphate associates with the Lcl-CTD surface via unique binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate binding mechanism.

Targeting the HER3 pseudokinase domain with small molecule inhibitors.

HER3 is a potent oncogenic growth factor receptor belonging to the human epidermal growth factor (HER/EGFR) family of receptor tyrosine kinases. In contrast to other EGFR family members, HER3 is a pseudokinase, lacking functional kinase activity. As such, efforts to develop small molecule tyrosine kinase inhibitors against this family member have been limited. In response to HER3-specific growth factors such as neuregulin (NRG, also known as heregulin or HRG), HER3 must couple with catalytically active family members, including its preferred partner HER2. Dimerization of the intracellular HER2:HER3 kinase domains is a critical part of the activation mechanism and HER3 plays a specialized role as an allosteric activator of the active HER2 kinase partner. Intriguingly, many pseudokinases retain functionally important nucleotide binding capacity, despite loss of kinase activity. We demonstrated that occupation of the nucleotide pocket of the pseudokinase HER3 retains functional importance for growth factor signaling through oncogenic HER2:HER3 heterodimers. Mutation of the HER3 nucleotide pocket both disrupts signaling and disrupts HER2:HER3 dimerization. Conversely, ATP competitive drugs which bind to HER3, but not HER2, can stabilize HER2:HER3 dimers, induce signaling and promote cell growth in breast cancer models. This indicates a nucleotide-dependent conformational role for the HER3 kinase domain. Critically, our recent proof-of-concept work demonstrated that HER3-directed small molecule inhibitors can also disrupt HER2:HER3 dimerization and signaling, supporting the prospect that HER3 can be a direct drug target despite its lack of intrinsic activity. In this chapter we will describe methods for identifying and validating small molecule inhibitors against the HER3 pseudokinase.

Reconstruction of ARNT PAS-B Unfolding Pathways by Steered Molecular Dynamics and Artificial Neural Networks.

Several experimental studies indicated that large conformational changes, including partial domain unfolding, have a role in the functional mechanisms of the basic helix loop helix Per/ARNT/SIM (bHLH-PAS) transcription factors. Recently, single-molecule atomic force microscopy (AFM) revealed two distinct pathways for the mechanical unfolding of the ARNT PAS-B. In this work we used steered molecular dynamics simulations to gain new insights into this process at an atomistic level. To reconstruct and classify pathways sampled in multiple simulations, we designed an original approach based on the use of self-organizing maps (SOMs). This led us to identify two types of unfolding pathways for the ARNT PAS-B, which are in good agreement with the AFM findings. Analysis of average forces mapped on the SOM revealed a stable conformation of the PAS-B along one pathway, which represents a possible structural model for the intermediate state detected by AFM. The approach here proposed will facilitate the study of other signal transmission mechanisms involving the folding/unfolding of PAS domains.

In silico phosphorylation of the autoinhibited form of p47(phox): insights into the mechanism of activation.

Activation of the multicomponent enzyme NADPH oxidase requires the interaction between the tandem SH3 domain of the cytosolic subunit p47(phox) and the cytoplasmic tail of membrane-bound p22(phox). In the resting state, p47(phox) exists in an autoinhibited conformation stabilized by intramolecular contacts between the SH3 domains and an adjacent polybasic region. Phosphorylation of three serine residues, Ser(303), Ser(304), and Ser(328) within this polybasic region has been shown to be sufficient for the disruption of the intramolecular interactions thereby inducing an active state of p47(phox). This active conformation is accessible to the cytoplasmic tail of p22(phox) and initiates the formation of the membrane-bound functional enzyme complex. Molecular dynamics simulations reveal insights in the mechanism of activation of the autoinhibited form of p47(phox) by in silico phosphorylation, of the three serine residues, Ser(303), Ser(304), and Ser(328). The simulations highlight the major collective coordinates generating the opening and the closing of the two SH3 domains and the residues that cause the unmasking of the p22(phox) binding site.

Structural alphabets derived from attractors in conformational space.

BACKGROUND: The hierarchical and partially redundant nature of protein structures justifies the definition of frequently occurring conformations of short fragments as 'states'. Collections of selected representatives for these states define Structural Alphabets, describing the most typical local conformations within protein structures. These alphabets form a bridge between the string-oriented methods of sequence analysis and the coordinate-oriented methods of protein structure analysis. RESULTS: A Structural Alphabet has been derived by clustering all four-residue fragments of a high-resolution subset of the protein data bank and extracting the high-density states as representative conformational states. Each fragment is uniquely defined by a set of three independent angles corresponding to its degrees of freedom, capturing in simple and intuitive terms the properties of the conformational space. The fragments of the Structural Alphabet are equivalent to the conformational attractors and therefore yield a most informative encoding of proteins. Proteins can be reconstructed within the experimental uncertainty in structure determination and ensembles of structures can be encoded with accuracy and robustness. CONCLUSIONS: The density-based Structural Alphabet provides a novel tool to describe local conformations and it is specifically suitable for application in studies of protein dynamics.

Energy landscapes associated with macromolecular conformational changes from endpoint structures.

Conformational changes modulate macromolecular function by promoting the specific binding of ligands (such as in antigen recognition) or the stabilization of transition states in enzymatic reactions. However, quantitative characterization of the energetics underlying dynamic structural interconversions is still challenging and lacks a unified method. Here, we introduce a novel in silico approach based on the combined use of essential dynamics sampling and nonequilibrium free-energy calculations to obtain quantitative data on conformational energy landscapes. This technique allows the unbiased investigation of highly complex rearrangements, and does not require the crucial definition of user-defined collective variables. We show that free-energy values derived from profiles connecting the unliganded and ligand-bound X-ray structures of a bacterial nucleoside hydrolase match the experimental binding constant. This approach also provides first evidence for a rate-limiting character of the conformational transition in this enzyme, and an unexpected role of the protonation state of a single residue in regulating substrate binding and product release.