Specialized Dynamical Properties of Promiscuous Residues Revealed by Simulated Conformational Ensembles

ThermoPCD: a database of molecular dynamics trajectories of antibody-antigen complexes at physiologic and fever-range temperatures

Progression of various cancers and autoimmune diseases is associated with changes in systemic or local tissue temperatures, which may impact current therapies. The role of fever and acute inflammation-range temperatures on the stability and activity of antibodies relevant for cancers and autoimmunity is unknown. To produce molecular dynamics (MD) trajectories of immune complexes at relevant temperatures, we used the Research Collaboratory for Structural Bioinformatics (RCSB) database to identify 50 antibody:antigen complexes of interest, in addition to single antibodies and antigens, and deployed Groningen Machine for Chemical Simulations (GROMACS) to prepare and run the structures at different temperatures for 100-500 ns, in single or multiple random seeds. MD trajectories are freely available. Processed data include Protein Data Bank outputs for all files obtained every 50 ns, and free binding energy calculations for some of the immune complexes. Protocols for using the data are also available. Individual datasets contain unique DOIs. We created a web interface, ThermoPCD, as a platform to explore the data. The outputs of ThermoPCD allow the users to relate thermally-dependent changes in epitopes:paratopes interfaces to their free binding energies, or against own experimentally derived binding affinities. ThermoPCD is a free to use database of immune complexes' trajectories at different temperatures that does not require registration and allows for all the data to be available for download. Database URL: https://sites.google.com/view/thermopcd/home.

ATLAS: protein flexibility description from atomistic molecular dynamics simulations

Dynamical behaviour is one of the most crucial protein characteristics. Despite the advances in the field of protein structure resolution and prediction, analysis and prediction of protein dynamic properties remains a major challenge, mostly due to the low accessibility of data and its diversity and heterogeneity. To address this issue, we present ATLAS, a database of standardised all-atom molecular dynamics simulations, accompanied by their analysis in the form of interactive diagrams and trajectory visualisation. ATLAS offers a large-scale view and valuable insights on protein dynamics for a large and representative set of proteins, by combining data obtained through molecular dynamics simulations with information extracted from experimental structures. Users can easily analyse dynamic properties of functional protein regions, such as domain limits (hinge positions) and residues involved in interaction with other biological molecules. Additionally, the database enables exploration of proteins with uncommon dynamic properties conditioned by their environment such as chameleon subsequences and Dual Personality Fragments. The ATLAS database is freely available at https://www.dsimb.inserm.fr/ATLAS.

Dynamozones are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease

Proteins are not static but are flexible molecules that can adopt many different conformations. The HIV-1 protease is an important target for the development of therapies to treat AIDS, due to its critical role in the viral life cycle. We investigated several dynamics studies on the HIV-1 protease families to illustrate the significance of examining the dynamic behaviors and molecular motions for an entire understanding of their dynamics-structure-function relationships. Using computer simulations and principal component analysis approaches, the dynamics data obtained revealed that: (i) The flap regions are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease; (ii) There are dynamic structural regions in some proteins that contribute to the biological function and allostery of proteins via appropriate flexibility. These regions are a clear sign of the evolution of conformational dynamics of proteins, which we call dynamozones. The flap regions are one of the most important dynamozones members that are critical for HIV-1 protease function. Due to the existence of other members of dynamozones in different proteins, we propose to consider dynamozones as a footprint of the evolution of the conformational dynamics of proteins.

Modeling the Dynamics of Protein–Protein Interfaces, How and Why?

Protein–protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein–protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein–protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein–protein interface.

Pathogenic missense protein variants affect different functional pathways and proteomic features than healthy population variants

Missense variants are present amongst the healthy population, but some of them are causative of human diseases. A classification of variants associated with “healthy” or “diseased” states is therefore not always straightforward. A deeper understanding of the nature of missense variants in health and disease, the cellular processes they may affect, and the general molecular principles which underlie these differences is essential to offer mechanistic explanations of the true impact of pathogenic variants. Here, we have formalised a statistical framework which enables robust probabilistic quantification of variant enrichment across full-length proteins, their domains, and 3D structure-defined regions. Using this framework, we validate and extend previously reported trends of variant enrichment in different protein structural regions (surface/core/interface). By examining the association of variant enrichment with available functional pathways and transcriptomic and proteomic (protein half-life, thermal stability, abundance) data, we have mined a rich set of molecular features which distinguish between pathogenic and population variants: Pathogenic variants mainly affect proteins involved in cell proliferation and nucleotide processing and are enriched in more abundant proteins. Additionally, rare population variants display features closer to common than pathogenic variants. We validate the association between these molecular features and variant pathogenicity by comparing against existing in silico variant impact annotations. This study provides molecular details into how different proteins exhibit resilience and/or sensitivity towards missense variants and provides the rationale to prioritise variant-enriched proteins and protein domains for therapeutic targeting and development. The ZoomVar database, which we created for this study, is available at fraternalilab.kcl.ac.uk/ZoomVar. It allows users to programmatically annotate missense variants with protein structural information and to calculate variant enrichment in different protein structural regions.

How do can one improve the classification of genetic variants as harmful or harmless? This study uses a robust statistical analysis to exploit the interplay between protein structure, proteomic measurements and functional pathways to enable better discrimination between missense variants in health and disease.

Allosteric Priming of E. coli CheY by the Flagellar Motor Protein FliM

Phosphorylation of Escherichia coli CheY protein transduces chemoreceptor stimulation to a highly cooperative flagellar motor response. CheY binds to the N-terminal peptide of the FliM motor protein (FliM_N). Constitutively active D13K-Y106W CheY has been an important tool for motor physiology. The crystal structures of CheY and CheY ⋅ FliM_N with and without D13K-Y106W have shown FliM_N-bound CheY contains features of both active and inactive states. We used molecular dynamics (MD) simulations to characterize the CheY conformational landscape accessed by FliM_N and D13K-Y106W. Mutual information measures identified the central features of the long-range CheY allosteric network between D57 phosphorylation site and the FliM_N interface, namely the closure of the α4-β4 hinge and inward rotation of Y- or W106 with W58. We used hydroxy-radical foot printing with mass spectroscopy (XFMS) to track the solvent accessibility of these and other side chains. The solution XFMS oxidation rate correlated with the solvent-accessible area of the crystal structures. The protection of allosteric relay side chains reported by XFMS confirmed the intermediate conformation of the native CheY ⋅ FliM_N complex, the inactive state of free D13K-Y106W CheY, and the MD-based network architecture. We extended the MD analysis to determine temporal coupling and energetics during activation. Coupled aromatic residue rotation was a graded rather than a binary switch, with Y- or W106 side-chain burial correlated with increased FliM_N affinity. Activation entrained CheY fold stabilization to FliM_N affinity. The CheY network could be partitioned into four dynamically coordinated sectors. Residue substitutions mapped to sectors around D57 or the FliM_N interface according to phenotype. FliM_N increased sector size and interactions. These sectors fused between the substituted K13-W106 residues to organize a tightly packed core and novel surfaces that may bind additional sites to explain the cooperative motor response. The community maps provide a more complete description of CheY priming than proposed thus far.

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.

Protein-Protein Interactions Mediated by Intrinsically Disordered Protein Regions Are Enriched in Missense Mutations

Because proteins are fundamental to most biological processes, many genetic diseases can be traced back to single nucleotide variants (SNVs) that cause changes in protein sequences. However, not all SNVs that result in amino acid substitutions cause disease as each residue is under different structural and functional constraints. Influential studies have shown that protein–protein interaction interfaces are enriched in disease-associated SNVs and depleted in SNVs that are common in the general population. These studies focus primarily on folded (globular) protein domains and overlook the prevalent class of protein interactions mediated by intrinsically disordered regions (IDRs). Therefore, we investigated the enrichment patterns of missense mutation-causing SNVs that are associated with disease and cancer, as well as those present in the healthy population, in structures of IDR-mediated interactions with comparisons to classical globular interactions. When comparing the different categories of interaction interfaces, division of the interface regions into solvent-exposed rim residues and buried core residues reveal distinctive enrichment patterns for the various types of missense mutations. Most notably, we demonstrate a strong enrichment at the interface core of interacting IDRs in disease mutations and its depletion in neutral ones, which supports the view that the disruption of IDR interactions is a mechanism underlying many diseases. Intriguingly, we also found an asymmetry across the IDR interaction interface in the enrichment of certain missense mutation types, which may hint at an increased variant tolerance and urges further investigations of IDR interactions.

Mining the PDB for Tractable Cases Where X-ray Crystallography Combined with Fragment Screens Can Be Used to Systematically Design Protein-Protein Inhibitors: Two Test Cases Illustrated by IL1β-IL1R and p38α-TAB1 Complexes

Nowadays, it is possible to combine X-ray crystallography and fragment screening in a medium throughput fashion to chemically probe the surfaces used by proteins to interact and use the outcome of the screens to systematically design protein-protein inhibitors. To prove it, we first performed a bioinformatics analysis of the Protein Data Bank protein complexes, which revealed over 400 cases where the crystal lattice of the target in the free form is such that large portions of the interacting surfaces are free from lattice contacts and therefore accessible to fragments during soaks. Among the tractable complexes identified, we then performed single fragment crystal screens on two particular interesting cases: the Il1β-ILR and p38α-TAB1 complexes. The result of the screens showed that fragments tend to bind in clusters, highlighting the small-molecule hotspots on the surface of the target protein. In most of the cases, the hotspots overlapped with the binding sites of the interacting proteins.

Structure, Dynamics and Cellular Insight Into Novel Substrates of the Legionella pneumophila Type II Secretion System

Legionella pneumophila is a Gram-negative bacterium that is able to replicate within a broad range of aquatic protozoan hosts. L. pneumophila is also an opportunistic human pathogen that can infect macrophages and epithelia in the lung and lead to Legionnaires’ disease. The type II secretion system is a key virulence factor of L. pneumophila and is used to promote bacterial growth at low temperatures, regulate biofilm formation, modulate host responses to infection, facilitate bacterial penetration of mucin gels and is necessary for intracellular growth during the initial stages of infection. The L. pneumophila type II secretion system exports at least 25 substrates out of the bacterium and several of these, including NttA to NttG, contain unique amino acid sequences that are generally not observed outside of the Legionella genus. NttA, NttC, and NttD are required for infection of several amoebal species but it is unclear what influence other novel substrates have within their host. In this study, we show that NttE is required for optimal infection of Acanthamoeba castellanii and Vermamoeba vermiformis amoeba and is essential for the typical colony morphology of L. pneumophila. In addition, we report the atomic structures of NttA, NttC, and NttE and through a combined biophysical and biochemical hypothesis driven approach we propose novel functions for these substrates during infection. This work lays the foundation for future studies into the mechanistic understanding of novel type II substrate functions and how these relate to L. pneumophila ecology and disease.

Functional cross-talk between allosteric effects of activating and inhibiting ligands underlies PKM2 regulation

Several enzymes can simultaneously interact with multiple intracellular metabolites, however, how the allosteric effects of distinct ligands are integrated to coordinately control enzymatic activity remains poorly understood. We addressed this question using, as a model system, the glycolytic enzyme pyruvate kinase M2 (PKM2). We show that the PKM2 activator fructose 1,6-bisphosphate (FBP) alone promotes tetramerisation and increases PKM2 activity, but addition of the inhibitor L-phenylalanine (Phe) prevents maximal activation of FBP-bound PKM2 tetramers. We developed a method, AlloHubMat, that uses eigenvalue decomposition of mutual information derived from molecular dynamics trajectories to identify residues that mediate FBP-induced allostery. Experimental mutagenesis of these residues identified PKM2 variants in which activation by FBP remains intact but cannot be attenuated by Phe. Our findings reveal residues involved in FBP-induced allostery that enable the integration of allosteric input from Phe and provide a paradigm for the coordinate regulation of enzymatic activity by simultaneous allosteric inputs.

Conformational coupling by trans-phosphorylation in calcium calmodulin dependent kinase II

The calcium calmodulin-dependent protein kinase II (CaMKII) is a dodecameric holoenzyme important for encoding memory. Its activation, triggered by binding of calcium-calmodulin, persists autonomously after calmodulin dissociation. One (receiver) kinase captures and subsequently phosphorylates the regulatory domain peptide of a donor kinase forming a chained dimer as the first stage of autonomous activation. Protein dynamics simulations examined the conformational changes triggered by dimer formation and phosphorylation, aimed to provide a molecular rationale for human mutations that result in learning disabilities. Ensembles generated from X-ray crystal structures were characterized by network centrality and community analysis. Mutual information related collective motions to local fragment dynamics encoded with a structural alphabet. Implicit solvent tCONCOORD conformational ensembles revealed the dynamic architecture of inactive kinase domains was co-opted in the activated dimer but the network hub shifted from the nucleotide binding cleft to the captured peptide. Explicit solvent molecular dynamics (MD) showed nucleotide and substrate binding determinants formed coupled nodes in long-range signal relays between regulatory peptides in the dimer. Strain in the extended captured peptide was balanced by reduced flexibility of the receiver kinase C-lobe core. The relays were organized around a hydrophobic patch between the captured peptide and a key binding helix. The human mutations aligned along the relays. Thus, these mutations could disrupt the allosteric network alternatively, or in addition, to altered binding affinities. Non-binding protein sectors distant from the binding sites mediated the allosteric signalling; providing possible targets for inhibitor design. Phosphorylation of the peptide modulated the dielectric of its binding pocket to strengthen the patch, non-binding sectors, domain interface and temporal correlations between parallel relays. These results provide the molecular details underlying the reported positive kinase cooperativity to enrich the discussion on how autonomous activation by phosphorylation leads to long-term behavioural effects.

Protein kinases play central roles in intracellular signalling. Auto-phosphorylation by bound nucleotide typically precedes phosphate transfer to multiple substrates. Protein conformational changes are central to kinase function, altering binding affinities to change cellular location and shunt from one signal pathway to another. In the brain, the multi-subunit kinase, CaMKII is activated by calcium-calmodulin upon calcium jumps produced by synaptic stimulation. Auto-transphosphorylation of a regulatory peptide enables the kinase to remain activated and mediate long-term behavioural effects after return to basal calcium levels. A database of mutated residues responsible for these effects is difficult to reconcile solely with impaired nucleotide or substrate binding. Therefore, we have computationally generated interaction networks to map the conformational plasticity of the kinase domains where most mutations localize. The network generated from the atomic structure of a phosphorylated dimer resolves protein sectors based on their collective motions. The sectors link nucleotide and substrate binding sites in self-reinforcing relays between regulatory peptides. The self-reinforcement is strengthened by phosphorylation consistent with the reported positive cooperativity of kinase activity with calcium-calmodulin concentration. The network gives a better match with the mutations and, in addition, reveals target sites for drug development.

Ligand-induced perturbation of the HIF-2α:ARNT dimer dynamics

Hypoxia inducible factors (HIFs) are transcription factors belonging to the basic helix−loop−helix PER-ARNT-SIM (bHLH-PAS) protein family with a role in sensing oxygen levels in the cell. Under hypoxia, the HIF-α degradation pathway is blocked and dimerization with the aryl hydrocarbon receptor nuclear translocator (ARNT) makes HIF-α transcriptionally active. Due to the common hypoxic environment of tumors, inhibition of this mechanism by destabilization of HIF-α:ARNT dimerization has been proposed as a promising therapeutic strategy. Following the discovery of a druggable cavity within the PAS-B domain of HIF-2α, research efforts have been directed to identify artificial ligands that can impair heterodimerization. Although the crystallographic structures of the HIF-2α:ARNT complex have elucidated the dimer architecture and the 0X3-inhibitor placement within the HIF-2α PAS-B, unveiling the inhibition mechanism requires investigation of how ligand-induced perturbations could dynamically propagate through the structure and affect dimerization. To this end, we compared evolutionary features, intrinsic dynamics and energetic properties of the dimerization interfaces of HIF-2α:ARNT in both the apo and holo forms. Residue conservation analysis highlighted inter-domain connecting elements that have a role in dimerization. Analysis of domain contributions to the dimerization energy demonstrated the importance of bHLH and PAS-A of both partners and of HIF-2α PAS-B domain in dimer stabilization. Among quaternary structure oscillations revealed by Molecular Dynamics simulations, the hinge-bending motion of the ARNT PAS-B domain around the flexible PAS-A/PAS-B linker supports a general model for ARNT dimerization in different heterodimers. Comparison of the HIF-2α:ARNT dynamics in the apo and 0X3-bound forms indicated a model of inhibition where the HIF-2α-PAS-B interfaces are destabilised as a result of water-bridged ligand-protein interactions and these local effects allosterically propagate to perturb the correlated motions of the domains and inter-domain communication. These findings will guide the design of improved inhibitors to contrast cell survival in tumor masses.

A low oxygen condition, called hypoxia, often occurs in tumor masses and generally correlates with worse prognosis. Cells in a tumor react to low oxygen levels with a metabolism modification induced by the activation of hypoxia inducible factors (HIFs) through dimerization with a partner protein and binding to a DNA target. Disrupting this protein-protein interaction could be a potential therapeutic strategy, but directly interfering with dimer formation can be troublesome because of the difficulty to design drugs that bind to protein interfaces. However, ligands that bind internal protein cavities can indirectly perturb the interfaces reducing dimers stability. Albeit protein crystallography had offered a detailed static picture of a HIF dimer bound to candidate inhibitors, it is not able to describe either the perturbation caused by binding or the molecular mechanism of dimer destabilization. Here we exploit molecular dynamics to identify the crucial interfaces in the HIF dimer stabilization and, by comparing the results obtained in the bound and unbound forms, we reveal the mechanism of ligand inhibition at atomic detail. All these findings will guide toward the design of improved dimerization inhibitors, to contrast cell survival in tumor masses.

Genetic variants and protein-protein interactions: a multidimensional network-centric view

We review recent progress in the mapping of genetic variants to proteins, in the context of their interactions, as measured from experiments and/or computational predictions. Such variants can impact on the molecular mechanisms underlying an interaction and its stability. We highlight recent work which relies on the effective use of protein-protein interaction networks (PPINs), integrated with 3D structural information, for evaluating disease-associated variants. Furthermore, we discuss how the integration of multiple layers of biological information, in the context of PPINs, can improve the interpretation of genetic variants and inspire new therapeutic strategies.

Investigating Ebola virus pathogenicity using molecular dynamics

BACKGROUND

Ebolaviruses have been known to cause deadly disease in humans for 40 years and have recently been demonstrated in West Africa to be able to cause large outbreaks. Four Ebolavirus species cause severe disease associated with high mortality in humans. Reston viruses are the only Ebolaviruses that do not cause disease in humans. Conserved amino acid changes in the Reston virus protein VP24 compared to VP24 of other Ebolaviruses have been suggested to alter VP24 binding to host cell karyopherins resulting in impaired inhibition of interferon signalling, which may explain the difference in human pathogenicity. Here we used protein structural analysis and molecular dynamics to further elucidate the interaction between VP24 and KPNA5.

RESULTS

As a control experiment, we compared the interaction of wild-type and R137A-mutant (known to affect KPNA5 binding) Ebola virus VP24 with KPNA5. Results confirmed that the R137A mutation weakens direct VP24-KPNA5 binding and enables water molecules to penetrate at the interface. Similarly, Reston virus VP24 displayed a weaker interaction with KPNA5 than Ebola virus VP24, which is likely to reduce the ability of Reston virus VP24 to prevent host cell interferon signalling.

CONCLUSION

Our results provide novel molecular detail on the interaction of Reston virus VP24 and Ebola virus VP24 with human KPNA5. The results indicate a weaker interaction of Reston virus VP24 with KPNA5 than Ebola virus VP24, which is probably associated with a decreased ability to interfere with the host cell interferon response. Hence, our study provides further evidence that VP24 is a key player in determining Ebolavirus pathogenicity.

Structural and evolutionary analyses reveal determinants of DNA binding specificities of nucleoid-associated proteins HU and IHF

Nucleoid-associated proteins (NAPs) are chromosome-organizing factors, which affect the transcriptional landscape of a bacterial cell. HU is an NAP, which binds to DNA with a broad specificity while homologous IHF (Integration Host Factor), binds DNA with moderately higher specificity. Specificity and differential binding affinity of HU/IHF proteins towards their target binding sites play a crucial role in their regulatory dynamics. Decades of biochemical and genomic studies have been carried out for HU and IHF like proteins. Yet, questions related to their DNA binding specificity, and differential ability to bend DNA thus affecting the binding site length remained unanswered. In addition, the problem has not been investigated from an evolutionary perspective. Our phylogenetic analysis revealed three major clades belonging to HU, IHFα and IHFβ like proteins with reference to E. coli. We carried out a comparative analysis of three-dimensional structures of HU/IHF proteins to gain insight into the structural basis of clade division. The present study revealed three major features which contribute to differential DNA binding specificity of HU/IHF proteins, (I) conformational restriction of DNA binding residues due to salt-bridge formation, (II) the enrichment of alanine in the DNA binding site increasing conformational space of flexible side chains in its vicinity and (III) nature of DNA binding residue (Arg to Lys bias in different clades) which interacts differentially to DNA bases. We observed an extended electropositive surface at the DNA draping site for IHF clade proteins compared to HU, which stabilizes the DNA bend. Differences in the dimer stabilization strategies between HU and IHF were also observed. Our analysis reveals a comprehensive evolutionary picture, which rationalizes the origin of multi-specificity of HU/IHF proteins using sequence and structure-based determinants, which could also be applied to understand differences in binding specificities of other nucleic acid binding proteins.

Promiscuous antibodies characterised by their physico-chemical properties: From sequence to structure and back

Human B cells produce antibodies, which bind to their cognate antigen based on distinct molecular properties of the antibody CDR loop. We have analysed a set of 10 antibodies showing a clear difference in their binding properties to a panel of antigens, resulting in two subsets of antibodies with a distinct binding phenotype. We call the observed binding multiplicity ‘promiscuous’ and selected physico-chemical CDRH3 characteristics and conformational preferences may characterise these promiscuous antibodies. To classify CDRH3 physico-chemical properties playing a role in their binding properties, we used statistical analyses of the sequences annotated by Kidera factors. To characterise structure-function requirements for antigen binding multiplicity we employed Molecular Modelling and Monte Carlo based coarse-grained simulations. The ability to predict the molecular causes of promiscuous, multi-binding behaviour would greatly improve the efficiency of the therapeutic antibody discovery process.

Deciphering Dimerization Modes of PAS Domains: Computational and Experimental Analyses of the AhR:ARNT Complex Reveal New Insights Into the Mechanisms of AhR Transformation

The Aryl hydrocarbon Receptor (AhR) is a transcription factor that mediates the biochemical response to xenobiotics and the toxic effects of a number of environmental contaminants, including dioxins. Recently, endogenous regulatory roles for the AhR in normal physiology and development have also been reported, thus extending the interest in understanding its molecular mechanisms of activation. Since dimerization with the AhR Nuclear Translocator (ARNT) protein, occurring through the Helix-Loop-Helix (HLH) and PER-ARNT-SIM (PAS) domains, is needed to convert the AhR into its transcriptionally active form, deciphering the AhR:ARNT dimerization mode would provide insights into the mechanisms of AhR transformation. Here we present homology models of the murine AhR:ARNT PAS domain dimer developed using recently available X-ray structures of other bHLH-PAS protein dimers. Due to the different reciprocal orientation and interaction surfaces in the different template dimers, two alternative models were developed for both the PAS-A and PAS-B dimers and they were characterized by combining a number of computational evaluations. Both well-established hot spot prediction methods and new approaches to analyze individual residue and residue-pairwise contributions to the MM-GBSA binding free energies were adopted to predict residues critical for dimer stabilization. On this basis, a mutagenesis strategy for both the murine AhR and ARNT proteins was designed and ligand-dependent DNA binding ability of the AhR:ARNT heterodimer mutants was evaluated. While functional analysis disfavored the HIF2α:ARNT heterodimer-based PAS-B model, most mutants derived from the CLOCK:BMAL1-based AhR:ARNT dimer models of both the PAS-A and the PAS-B dramatically decreased the levels of DNA binding, suggesting this latter model as the most suitable for describing AhR:ARNT dimerization. These novel results open new research directions focused at elucidating basic molecular mechanisms underlying the functional activity of the AhR.

Computational modeling combined with experimental validation may give insight into structural and functional properties of protein systems. The basic Helix-Loop-Helix PER-ARNT-SIM (bHLH-PAS) proteins show conserved functional domains despite the broad range of functions exerted by the different systems. Within this protein family, the Aryl hydrocarbon Receptor (AhR) is known to mediate the toxic effects of a number of environmental contaminants, including dioxins and dioxin-like chemicals, and it also exerts other biochemical and physiological effects. Despite the absence of experimentally determined structures, theoretical models of the AhR PAS domains developed on the basis of homologous systems have allowed understanding of some aspects of the molecular mechanisms underlying its function. In this work we present alternative structural models of the transcriptionally active complex of AhR with the AhR Nuclear Translocator (ARNT) protein. Computational characterization of the modeled protein-protein interaction interfaces guided the design of mutagenesis experiments, and evaluation of the DNA binding ability of the resulting AhR:ARNT dimer mutants allowed validation of the models and selection of the most reliable one. These findings open new research directions for understanding the molecular mechanisms underlying the functional activity of the AhR.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

The Phylogenetic Signature Underlying ATP Synthase c-Ring Compliance

The proton-driven ATP synthase (FOF1) is comprised of two rotary, stepping motors (FO and F1) coupled by an elastic power transmission. The elastic compliance resides in the rotor module that includes the membrane-embedded FO c-ring. Proton transport by FO is firmly coupled to the rotation of the c-ring relative to other FO subunits (ab2). It drives ATP synthesis. We used a computational method to investigate the contribution of the c-ring to the total elastic compliance. We performed principal component analysis of conformational ensembles built using distance constraints from the bovine mitochondrial c-ring x-ray structure. Angular rotary twist, the dominant ring motion, was estimated to show that the c-ring accounted in part for the measured compliance. Ring rotation was entrained to rotation of the external helix within each hairpin-shaped c-subunit in the ring. Ensembles of monomer and dimers extracted from complete c-rings showed that the coupling between collective ring and the individual subunit motions was independent of the size of the c-ring, which varies between organisms. Molecular determinants were identified by covariance analysis of residue coevolution and structural-alphabet-based local dynamics correlations. The residue coevolution gave a readout of subunit architecture. The dynamic couplings revealed that the hinge for both ring and subunit helix rotations was constructed from the proton-binding site and the adjacent glycine motif (IB-GGGG) in the midmembrane plane. IB-GGGG motifs were linked by long-range couplings across the ring, while intrasubunit couplings connected the motif to the conserved cytoplasmic loop and adjacent segments. The correlation with principal collective motions shows that the couplings underlie both ring rotary and bending motions. Noncontact couplings between IB-GGGG motifs matched the coevolution signal as well as contact couplings. The residue coevolution reflects the physiological importance of the dynamics that may link proton transfer to ring compliance.

PinSnps: structural and functional analysis of SNPs in the context of protein interaction networks

Summary: We present a practical computational pipeline to readily perform data analyses of protein–protein interaction networks by using genetic and functional information mapped onto protein structures. We provide a 3D representation of the available protein structure and its regions (surface, interface, core and disordered) for the selected genetic variants and/or SNPs, and a prediction of the mutants’ impact on the protein as measured by a range of methods. We have mapped in total 2587 genetic disorder-related SNPs from OMIM, 587 873 cancer-related variants from COSMIC, and 1 484 045 SNPs from dbSNP. All result data can be downloaded by the user together with an R-script to compute the enrichment of SNPs/variants in selected structural regions.

Availability and Implementation: PinSnps is available as open-access service at http://fraternalilab.kcl.ac.uk/PinSnps/

Contact:franca.fraternali@kcl.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Using Local States To Drive the Sampling of Global Conformations in Proteins

Conformational changes associated with protein function often occur beyond the time scale currently accessible to unbiased molecular dynamics (MD) simulations, so that different approaches have been developed to accelerate their sampling. Here we investigate how the knowledge of backbone conformations preferentially adopted by protein fragments, as contained in precalculated libraries known as structural alphabets (SA), can be used to explore the landscape of protein conformations in MD simulations. We find that (a) enhancing the sampling of native local states in both metadynamics and steered MD simulations allows the recovery of global folded states in small proteins; (b) folded states can still be recovered when the amount of information on the native local states is reduced by using a low-resolution version of the SA, where states are clustered into macrostates; and (c) sequences of SA states derived from collections of structural motifs can be used to sample alternative conformations of preselected protein regions. The present findings have potential impact on several applications, ranging from protein model refinement to protein folding and design.

Anatomy of protein disorder, flexibility and disease-related mutations

Integration of protein structural information with human genetic variation and pathogenic mutations is essential to understand molecular mechanisms associated with the effects of polymorphisms on protein interactions and cellular processes. We investigate occurrences of non-synonymous SNPs in ordered and disordered protein regions by systematic mapping of common variants and disease-related SNPs onto these regions. We show that common variants accumulate in disordered regions; conversely pathogenic variants are significantly depleted in disordered regions. These different occurrences of pathogenic and common SNPs can be attributed to a negative selection on random mutations in structurally highly constrained regions. New approaches in the study of quantitative effects of pathogenic-related mutations should effectively account for all the possible contexts and relative functional constraints in which the sequence variation occurs.

Protein flexibility in the light of structural alphabets

Protein structures are valuable tools to understand protein function. Nonetheless, proteins are often considered as rigid macromolecules while their structures exhibit specific flexibility, which is essential to complete their functions. Analyses of protein structures and dynamics are often performed with a simplified three-state description, i.e., the classical secondary structures. More precise and complete description of protein backbone conformation can be obtained using libraries of small protein fragments that are able to approximate every part of protein structures. These libraries, called structural alphabets (SAs), have been widely used in structure analysis field, from definition of ligand binding sites to superimposition of protein structures. SAs are also well suited to analyze the dynamics of protein structures. Here, we review innovative approaches that investigate protein flexibility based on SAs description. Coupled to various sources of experimental data (e.g., B-factor) and computational methodology (e.g., Molecular Dynamic simulation), SAs turn out to be powerful tools to analyze protein dynamics, e.g., to examine allosteric mechanisms in large set of structures in complexes, to identify order/disorder transition. SAs were also shown to be quite efficient to predict protein flexibility from amino-acid sequence. Finally, in this review, we exemplify the interest of SAs for studying flexibility with different cases of proteins implicated in pathologies and diseases.

Evolution of conformational dynamics determines the conversion of a promiscuous generalist into a specialist enzyme

β-Lactamases are produced by many modern bacteria as a mechanism of resistance toward β-lactam antibiotics, the most common antibiotics in use. β-Lactamases, however, are ancient enzymes that originated billions of years ago. Recently, proteins corresponding to 2- to 3-Gy-old Precambrian nodes in the evolution of Class A β-lactamases have been prepared and shown to be moderately efficient promiscuous catalysts, able to degrade a variety of antibiotics with catalytic efficiency levels similar to those of an average modern enzyme. Remarkably, there are few structural differences (in particular at the active-site regions) between the resurrected enzymes and a penicillin-specialist modern β-lactamase. Here, we propose that the ancestral promiscuity originates from conformational dynamics. We investigate the differences in conformational dynamics of the ancient and extant β-lactamases through MD simulations and quantify the contribution of each position to functionally related dynamics through Dynamic Flexibility Index. The modern TEM-1 lactamase shows a comparatively rigid active-site region, likely reflecting adaptation for efficient degradation of a specific substrate (penicillin), whereas enhanced deformability at the active-site neighborhood in the ancestral resurrected proteins likely accounts for the binding and subsequent degradation of antibiotic molecules of different size and shape. Clustering of the conformational dynamics on the basis of Principal Component Analysis is in agreement with the functional divergence, as the ancient β-lactamases cluster together, separated from their modern descendant. Finally, our analysis leads to testable predictions, as sites of potential relevance for the evolution of dynamics are identified and mutations at those sites are expected to alter substrate-specificity.