Protein flexibility, not disorder, is intrinsic to molecular recognition

AlphaFold2 modeling and molecular dynamics simulations of an intrinsically disordered protein

We use AlphaFold2 (AF2) to model the monomer and dimer structures of an intrinsically disordered protein (IDP), Nvjp-1, assisted by molecular dynamics (MD) simulations. We observe relatively rigid dimeric structures of Nvjp-1 when compared with the monomer structures. We suggest that protein conformations from multiple AF2 models and those from MD trajectories exhibit a coherent trend: the conformations of an IDP are deviated from each other and the conformations of a well-folded protein are consistent with each other. We use a residue-residue interaction network (RIN) derived from the contact map which show that the residue-residue interactions in Nvjp-1 are mainly transient; however, those in a well-folded protein are mainly persistent. Despite the variation in 3D shapes, we show that the AF2 models of both disordered and ordered proteins exhibit highly consistent profiles of the pLDDT (predicted local distance difference test) scores. These results indicate a potential protocol to justify the IDPs based on multiple AF2 models and MD simulations.

Exploring the Role of Globular Domain Locations on an Intrinsically Disordered Region of p53: A Molecular Dynamics Investigation

The pre-tetramerization loop (PTL) of the human tumor suppressor protein p53 is an intrinsically disordered region (IDR) necessary for the tetramerization process, and its flexibility contributes to the essential conformational changes needed. Although the IDR can be accurately simulated in the traditional manner of molecular dynamics (MD) with the end-to-end distance (EEdist) unhindered, we sought to explore the effects of restraining the EEdist to the values predicted by electron microscopy (EM) and other distances. Simulating the PTL trajectory with a restrained EEdist , we found an increased agreement of nuclear magnetic resonance (NMR) chemical shifts with experiments. Additionally, we observed a plethora of secondary structures and contacts that only appear when the trajectory is restrained. Our findings expand the understanding of the tetramerization of p53 and provide insight into how mutations could make the protein impotent. In particular, our findings demonstrate the importance of restraining the EEdist in studying IDRs and how their conformations change under different conditions. Our results provide a better understanding of the PTL and the conformational dynamics of IDRs in general, which are useful for further studies regarding mutations and their effects on the activity of p53.

Non-canonical helical transitions and conformational switching are associated with characteristic flexibility and disorder indices in TRP and Kv channels

Structural evidence and much experimental data have demonstrated the presence of non-canonical helical substructures (π and 3₁₀) in regions of great functional relevance both in TRP as in Kv channels. Through an exhaustive compositional analysis of the sequences underlying these substructures, we find that each of them is associated with characteristic local flexibility profiles, which in turn are implicated in significant conformational rearrangements and interactions with specific ligands. We found that α-to-π helical transitions are associated with patterns of local rigidity whereas α-to-3₁₀ transitions are mainly leagued with high local flexibility profiles. We also study the relationship between flexibility and protein disorder in the transmembrane domain of these proteins. By contrasting these two parameters, we located regions showing a sort of structural discrepancy between these similar but not identical protein attributes. Notably, these regions are presumably implicated in important conformational rearrangements during the gating in those channels. In that sense, finding these regions where flexibility and disorder are not proportional allows us to detect regions with potential functional dynamism. From this point of view, we highlighted some conformational rearrangements that occur during ligand binding events, the compaction, and refolding of the outer pore loops in several TRP channels, as well as the well-known S4 motion in Kv channels.

In-cell NMR: Why and how?

NMR spectroscopy has been applied to cells and tissues analysis since its beginnings, as early as 1950. We have attempted to gather here in a didactic fashion the broad diversity of data and ideas that emerged from NMR investigations on living cells. Covering a large proportion of the periodic table, NMR spectroscopy permits scrutiny of a great variety of atomic nuclei in all living organisms non-invasively. It has thus provided quantitative information on cellular atoms and their chemical environment, dynamics, or interactions. We will show that NMR studies have generated valuable knowledge on a vast array of cellular molecules and events, from water, salts, metabolites, cell walls, proteins, nucleic acids, drugs and drug targets, to pH, redox equilibria and chemical reactions. The characterization of such a multitude of objects at the atomic scale has thus shaped our mental representation of cellular life at multiple levels, together with major techniques like mass-spectrometry or microscopies. NMR studies on cells has accompanied the developments of MRI and metabolomics, and various subfields have flourished, coined with appealing names: fluxomics, foodomics, MRI and MRS (i.e. imaging and localized spectroscopy of living tissues, respectively), whole-cell NMR, on-cell ligand-based NMR, systems NMR, cellular structural biology, in-cell NMR… All these have not grown separately, but rather by reinforcing each other like a braided trunk. Hence, we try here to provide an analytical account of a large ensemble of intricately linked approaches, whose integration has been and will be key to their success. We present extensive overviews, firstly on the various types of information provided by NMR in a cellular environment (the "why", oriented towards a broad readership), and secondly on the employed NMR techniques and setups (the "how", where we discuss the past, current and future methods). Each subsection is constructed as a historical anthology, showing how the intrinsic properties of NMR spectroscopy and its developments structured the accessible knowledge on cellular phenomena. Using this systematic approach, we sought i) to make this review accessible to the broadest audience and ii) to highlight some early techniques that may find renewed interest. Finally, we present a brief discussion on what may be potential and desirable developments in the context of integrative studies in biology.

AlphaFold2 models indicate that protein sequence determines both structure and dynamics

AlphaFold 2 (AF2) has placed Molecular Biology in a new era where we can visualize, analyze and interpret the structures and functions of all proteins solely from their primary sequences. We performed AF2 structure predictions for various protein systems, including globular proteins, a multi-domain protein, an intrinsically disordered protein (IDP), a randomized protein, two larger proteins (> 1000 AA), a heterodimer and a homodimer protein complex. Our results show that along with the three dimensional (3D) structures, AF2 also decodes protein sequences into residue flexibilities via both the predicted local distance difference test (pLDDT) scores of the models, and the predicted aligned error (PAE) maps. We show that PAE maps from AF2 are correlated with the distance variation (DV) matrices from molecular dynamics (MD) simulations, which reveals that the PAE maps can predict the dynamical nature of protein residues. Here, we introduce the AF2-scores, which are simply derived from pLDDT scores and are in the range of [0, 1]. We found that for most protein models, including large proteins and protein complexes, the AF2-scores are highly correlated with the root mean square fluctuations (RMSF) calculated from MD simulations. However, for an IDP and a randomized protein, the AF2-scores do not correlate with the RMSF from MD, especially for the IDP. Our results indicate that the protein structures predicted by AF2 also convey information of the residue flexibility, i.e., protein dynamics.

Mechanism of GAPDH Redox Signaling by H2O2 Activation of a Two-Cysteine Switch

Oxidation of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) by reactive oxygen species such as H2O2 activate pleiotropic signaling pathways is associated with pathophysiological cell fate decisions. Oxidized GAPDH binds chaperone proteins with translocation of the complex to the nucleus and mitochondria initiating autophagy and cellular apoptosis. In this study, we establish the mechanism by which H2O2-oxidized GAPDH subunits undergo a subunit conformational rearrangement. H2O2 oxidizes both the catalytic cysteine and a vicinal cysteine (four residues downstream) to their respective sulfenic acids. A 'two-cysteine switch' is activated, whereby the sulfenic acids irreversibly condense to an intrachain thiosulfinic ester resulting in a major metastable subunit conformational rearrangement. All four subunits of the homotetramer are uniformly and independently oxidized by H2O2, and the oxidized homotetramer is stabilized at low temperatures. Over time, subunits unfold forming disulfide-linked aggregates with the catalytic cysteine oxidized to a sulfinic acid, resulting from thiosulfinic ester hydrolysis via the highly reactive thiosulfonic ester intermediate. Molecular Dynamic Simulations provide additional mechanistic insights linking GAPDH subunit oxidation with generating a putative signaling conformer. The low-temperature stability of the H2O2-oxidized subunit conformer provides an operable framework to study mechanisms associated with gain-of-function activities of oxidized GAPDH to identify novel targets for the treatment of neurodegenerative diseases.

Nonspecific Binding-Fundamental Concepts and Consequences for Biosensing Applications

Nature achieves differentiation of specific and nonspecific binding in molecular interactions through precise control of biomolecules in space and time. Artificial systems such as biosensors that rely on distinguishing specific molecular binding events in a sea of nonspecific interactions have struggled to overcome this issue. Despite the numerous technological advancements in biosensor technologies, nonspecific binding has remained a critical bottleneck due to the lack of a fundamental understanding of the phenomenon. To date, the identity, cause, and influence of nonspecific binding remain topics of debate within the scientific community. In this review, we discuss the evolution of the concept of nonspecific binding over the past five decades based upon the thermodynamic, intermolecular, and structural perspectives to provide classification frameworks for biomolecular interactions. Further, we introduce various theoretical models that predict the expected behavior of biosensors in physiologically relevant environments to calculate the theoretical detection limit and to optimize sensor performance. We conclude by discussing existing practical approaches to tackle the nonspecific binding challenge in vitro for biosensing platforms and how we can both address and harness nonspecific interactions for in vivo systems.

MEDUSA: Prediction of Protein Flexibility from Sequence

Information on the protein flexibility is essential to understand crucial molecular mechanisms such as protein stability, interactions with other molecules and protein functions in general. B-factor obtained in the X-ray crystallography experiments is the most common flexibility descriptor available for the majority of the resolved protein structures. Since the gap between the number of the resolved protein structures and available protein sequences is continuously growing, it is important to provide computational tools for protein flexibility prediction from amino acid sequence. In the current study, we report a Deep Learning based protein flexibility prediction tool MEDUSA (https://www.dsimb.inserm.fr/MEDUSA). MEDUSA uses evolutionary information extracted from protein homologous sequences and amino acid physico-chemical properties as input for a convolutional neural network to assign a flexibility class to each protein sequence position. Trained on a non-redundant dataset of X-ray structures, MEDUSA provides flexibility prediction in two, three and five classes. MEDUSA is freely available as a web-server providing a clear visualization of the prediction results as well as a standalone utility (https://github.com/DSIMB/medusa). Analysis of the MEDUSA output allows a user to identify the potentially highly deformable protein regions and general dynamic properties of the protein.

Liquid-Liquid Phase Separation by Intrinsically Disordered Protein Regions of Viruses: Roles in Viral Life Cycle and Control of Virus-Host Interactions

Intrinsically disordered proteins (IDPs) are unable to adopt a unique 3D structure under physiological conditions and thus exist as highly dynamic conformational ensembles. IDPs are ubiquitous and widely spread in the protein realm. In the last decade, compelling experimental evidence has been gathered, pointing to the ability of IDPs and intrinsically disordered regions (IDRs) to undergo liquid-liquid phase separation (LLPS), a phenomenon driving the formation of membrane-less organelles (MLOs). These biological condensates play a critical role in the spatio-temporal organization of the cell, where they exert a multitude of key biological functions, ranging from transcriptional regulation and silencing to control of signal transduction networks. After introducing IDPs and LLPS, we herein survey available data on LLPS by IDPs/IDRs of viral origin and discuss their functional implications. We distinguish LLPS associated with viral replication and trafficking of viral components, from the LLPS-mediated interference of viruses with host cell functions. We discuss emerging evidence on the ability of plant virus proteins to interfere with the regulation of MLOs of the host and propose that bacteriophages can interfere with bacterial LLPS, as well. We conclude by discussing how LLPS could be targeted to treat phase separation-associated diseases, including viral infections.

MODELLING OF 3D-STRUCTURES OF THE RARE MELANOCORTIN-1-RECEPTOR MUTATIONS ASSOCIATED TO MELANISM IN THE BANANAQUIT

Melanism in plumage color is often associated to the single nucleotide polymorphism of the melanocortin-1-receptor (MC1R). Despite the striking association between the substitution of a Glutamic-acid by for a Lysine at position 92 on the MC1R protein and a completely black plumage, an in-depth understanding of the effect of missense mutations on the conformational change and behavior of the MC1R in the lipid bilayer caused by the absence of a crystal structure is lacking. We examine the structural basis for receptor activation using DNA sequences from the GenBank to perform in silicoprotein homology-based modeling. Our tridimensional model shows that the Alanine for a 179-Threoninesubstitution is a structural complement of the charge-reversing effect associated to the substitution of a Glutamic-acid by for a Lysine at position 92 on the MC1R. We proposed the possibility of gradual evolution in stability and electrostatic properties of the MC1R by the sequential accumulation of these two rare substitutions. These two rare substitutions further perturb physical-chemical properties that may be necessary folding requirements of the constitutively active MC1R forms without altering of ligand binding affinity. The computational coarse-grained molecular dynamics of the MC1R binding affinities to the melanocyte-stimulating hormone predicted the disparity in ligand binding amongalleles. We speculate that the disparity in structural constraints and ligand binding among the alleles within heterozygous individuals may contribute as a mechanism to the plumage color variation in the Coereba flaveola.

Understanding the Impacts of Conformational Dynamics on the Regulation of Protein-Protein Association by a Multiscale Simulation Method

Complexes formed among diverse proteins carry out versatile functions in nearly all physiological processes. Association rates which measure how fast proteins form various complexes are of fundamental importance to characterize their functions. The association rates are not only determined by the energetic features at binding interfaces of a protein complex but also influenced by the intrinsic conformational dynamics of each protein in the complex. Unfortunately, how this conformational effect regulates protein association has never been calibrated on a systematic level. To tackle this problem, we developed a multiscale strategy to incorporate the information on protein conformational variations from Langevin dynamic simulations into a kinetic Monte Carlo algorithm of protein-protein association. By systematically testing this approach against a large-scale benchmark set, we found the association of a protein complex with a relatively rigid structure tends to be reduced by its conformational fluctuations. With specific examples, we further show that higher degrees of structural flexibility in various protein complexes can facilitate the searching and formation of intermolecular interactions and thereby accelerate their associations. In general, the integration of conformational dynamics can improve the correlation between experimentally measured association rates and computationally derived association probabilities. Finally, we analyzed the statistical distributions of different secondary structural types on protein-protein binding interfaces and their preference to the change of association rates. Our study, to the best of our knowledge, is the first computational method that systematically estimates the impacts of protein conformational dynamics on protein-protein association. It throws lights on the molecular mechanisms of how protein-protein recognition is kinetically modulated.

Disorder and cysteines in proteins: A design for orchestration of conformational see-saw and modulatory functions

Being responsible for more than 90% of cellular functions, protein molecules are workhorses in all the life forms. In order to cater for such a high demand, proteins have evolved to adopt diverse structures that allow them to perform myriad of functions. Beginning with the genetically directed amino acid sequence, the classical understanding of protein function involves adoption of hierarchically complex yet ordered structures. However, advances made over the last two decades have revealed that inasmuch as 50% of eukaryotic proteome exists as partially or fully disordered structures. Significance of such intrinsically disordered proteins (IDPs) is further realized from their ability to exhibit multifunctionality, a feature attributable to their conformational plasticity. Among the coded amino acids, cysteines are considered to be "order-promoting" due to their ability to form inter- or intramolecular disulfide bonds, which confer robust thermal stability to the protein structure in oxidizing conditions. The co-existence of order-promoting cysteines with disorder-promoting sequences seems counter-intuitive yet many proteins have evolved to contain such sequences. In this chapter, we review some of the known cysteine-containing protein domains categorized based on the number of cysteines they possess. We show that many protein domains contain disordered sequences interspersed with cysteines. We show that a positive correlation exists between the degree of cysteines and disorder within the sequences that flank them. Furthermore, based on the computational platform, IUPred2A, we show that cysteine-rich sequences display significant disorder in the reduced but not the oxidized form, increasing the potential for such sequences to function in a redox-sensitive manner. Overall, this chapter provides insights into an exquisite evolutionary design wherein disordered sequences with interspersed cysteines enable potential modulatory protein functions under stress and environmental conditions, which thus far remained largely inconspicuous.

Life in Phases: Intra- and Inter- Molecular Phase Transitions in Protein Solutions

Proteins, these evolutionarily-edited biological polymers, are able to undergo intramolecular and intermolecular phase transitions. Spontaneous intramolecular phase transitions define the folding of globular proteins, whereas binding-induced, intra- and inter- molecular phase transitions play a crucial role in the functionality of many intrinsically-disordered proteins. On the other hand, intermolecular phase transitions are the behind-the-scenes players in a diverse set of macrosystemic phenomena taking place in protein solutions, such as new phase nucleation in bulk, on the interface, and on the impurities, protein crystallization, protein aggregation, the formation of amyloid fibrils, and intermolecular liquid-liquid or liquid-gel phase transitions associated with the biogenesis of membraneless organelles in the cells. This review is dedicated to the systematic analysis of the phase behavior of protein molecules and their ensembles, and provides a description of the major physical principles governing intramolecular and intermolecular phase transitions in protein solutions.

Perspective: the essential role of NMR in the discovery and characterization of intrinsically disordered proteins

The 2019 ISMAR Prize recognized NMR studies of disordered proteins. Here we provide a highly personal perspective on the discovery of intrinsically disordered proteins and the development and application of NMR methods to characterize their conformational ensembles, dynamics, and interactions.

Oncoviruses Can Drive Cancer by Rewiring Signaling Pathways Through Interface Mimicry

Oncoviruses rewire host pathways to subvert host immunity and promote their survival and proliferation. However, exactly how is challenging to understand. Here, by employing the first and to date only interface-based host-microbe interaction (HMI) prediction method, we explore a pivotal strategy oncoviruses use to drive cancer: mimicking binding surfaces-interfaces-of human proteins. We show that oncoviruses can target key human network proteins and transform cells by acquisition of cancer hallmarks. Experimental large-scale mapping of HMIs is difficult and individual HMIs do not permit in-depth grasp of tumorigenic virulence mechanisms. Our computational approach is tractable and 3D structural HMI models can help elucidate pathogenesis mechanisms and facilitate drug design. We observe that many host proteins are unique targets for certain oncoviruses, whereas others are common to several, suggesting similar infectious strategies. A rough estimation of our false discovery rate based on the tissue expression of oncovirus-targeted human proteins is 25%.

Prediction of Host-Pathogen Interactions for Helicobacter pylori by Interface Mimicry and Implications to Gastric Cancer

There is a strong correlation between some pathogens and certain cancer types. One example is Helicobacter pylori and gastric cancer. Exactly how they contribute to host tumorigenesis is, however, a mystery. Pathogens often interact with the host through proteins. To subvert defense, they may mimic host proteins at the sequence, structure, motif, or interface levels. Interface similarity permits pathogen proteins to compete with those of the host for a target protein and thereby alter the host signaling. Detection of host-pathogen interactions (HPIs) and mapping the re-wired superorganism HPI network-with structural details-can provide unprecedented clues to the underlying mechanisms and help therapeutics. Here, we describe the first computational approach exploiting solely interface mimicry to model potential HPIs. Interface mimicry can identify more HPIs than sequence or complete structural similarity since it appears more common than the other mimicry types. We illustrate the usefulness of this concept by modeling HPIs of H. pylori to understand how they modulate host immunity, persist lifelong, and contribute to tumorigenesis. H. pylori proteins interfere with multiple host pathways as they target several host hub proteins. Our results help illuminate the structural basis of resistance to apoptosis, immune evasion, and loss of cell junctions seen in H. pylori-infected host cells.

Biophysical exploration of protein-flavonol recognition: effects of molecular properties and conformational flexibility

The current work explores the biomolecular recognition of a series of flavonols by a protein and then uncovers the influences of the structural features of flavonols and the protein's own characteristics, e.g. the dynamics and flexibility, on the bioavailability of flavonols by using the pivotal biomacromolecule hemoglobin as a model. The experimental results revealed that flavonol may lead to a notable decrease in the steady-state fluorescence intensity of the β-37 Trp residue, and in the meantime the R-T transition of the protein transpired. Such noncovalent recognition forms the ground-state adduct, with an association intensity of 3.991 × 10(4) M(-1) in the reaction process, which has already been authenticated by the detailed analysis of time-resolved fluorescence and UV/vis absorption spectra. Furthermore, flavonol can form hydrogen bonds and π-conjugation effects with several amino acid residues on the polypeptide chain, for example, Trp-37, Arg-40, Asp-99 and Asn-102, and this event would induce self-regulation of the compact, regular conformation of the protein to a certain extent, which explicitly corroborates the results of circular dichroism. According to the study of molecular docking and structure-activity relationships, we could see that the recognition capacities of the protein-flavonols are inversely interrelated with the C log P values of the flavonol molecules. Moreover, the properties of the substituents in the structural B-ring unit of flavonols, i.e. polarity, position and number, will also prominently affect the degree of affinity and bioavailability of the protein-flavonol complexes. The analytical results of molecular dynamics (MD) simulation testified that the discussions of the structure-activity relationships are entirely logical, and the conformations of the amino acid residues forming noncovalent interactions tend to be stable in the MD simulation, as further elucidated from the dynamics data. Plainly, molecular recognition of the protein-flavonols might noticeably cause relatively large changes in protein flexibility, and then manifest different recognition strengths and corresponding biological activities. This issue will be carefully validated by the interpretation of root-mean-square fluctuation.

Common structural features of cholesterol binding sites in crystallized soluble proteins

Cholesterol-protein interactions are essential for the architectural organization of cell membranes and for lipid metabolism. While cholesterol-sensing motifs in transmembrane proteins have been identified, little is known about cholesterol recognition by soluble proteins. We reviewed the structural characteristics of binding sites for cholesterol and cholesterol sulfate from crystallographic structures available in the Protein Data Bank. This analysis unveiled key features of cholesterol-binding sites that are present in either all or the majority of sites: i) the cholesterol molecule is generally positioned between protein domains that have an organized secondary structure; ii) the cholesterol hydroxyl/sulfo group is often partnered by Asn, Gln, and/or Tyr, while the hydrophobic part of cholesterol interacts with Leu, Ile, Val, and/or Phe; iii) cholesterol hydrogen-bonding partners are often found on α-helices, while amino acids that interact with cholesterol's hydrophobic core have a slight preference for β-strands and secondary structure-lacking protein areas; iv) the steroid's C21 and C26 constitute the "hot spots" most often seen for steroid-protein hydrophobic interactions; v) common "cold spots" are C8-C10, C13, and C17, at which contacts with the proteins were not detected. Several common features we identified for soluble protein-steroid interaction appear evolutionarily conserved.

Electron Microscopy Structural Insights into CPAP Oligomeric Behavior: A Plausible Assembly Process of a Supramolecular Scaffold of the Centrosome

Centrosomal P4.1-associated protein (CPAP) is a cell cycle regulated protein fundamental for centrosome assembly and centriole elongation. In humans, the region between residues 897–1338 of CPAP mediates interactions with other proteins and includes a homodimerization domain. CPAP mutations cause primary autosomal recessive microcephaly and Seckel syndrome. Despite of the biological/clinical relevance of CPAP, its mechanistic behavior remains unclear and its C-terminus (the G-box/TCP domain) is the only part whose structure has been solved. This situation is perhaps due in part to the challenges that represent obtaining the protein in a soluble, homogeneous state for structural studies. Our work constitutes a systematic structural analysis on multiple oligomers of HsCPAP^897−1338, using single-particle electron microscopy (EM) of negatively stained (NS) samples. Based on image classification into clearly different regular 3D maps (putatively corresponding to dimers and tetramers) and direct observation of individual images representing other complexes of HsCPAP^897−1338 (i.e., putative flexible monomers and higher-order multimers), we report a dynamic oligomeric behavior of this protein, where different homo-oligomers coexist in variable proportions. We propose that dimerization of the putative homodimer forms a putative tetramer which could be the structural unit for the scaffold that either tethers the pericentriolar material to centrioles or promotes procentriole elongation. A coarse fitting of atomic models into the NS 3D maps at resolutions around 20 Å is performed only to complement our experimental data, allowing us to hypothesize on the oligomeric composition of the different complexes. In this way, the current EM work represents an initial step toward the structural characterization of different oligomers of CPAP, suggesting further insights to understand how this protein works, contributing to the elucidation of control mechanisms for centriole biogenesis.

Conformational diversity analysis reveals three functional mechanisms in proteins

Protein motions are a key feature to understand biological function. Recently, a large-scale analysis of protein conformational diversity showed a positively skewed distribution with a peak at 0.5 Å C-alpha root-mean-square-deviation (RMSD). To understand this distribution in terms of structure-function relationships, we studied a well curated and large dataset of ~5,000 proteins with experimentally determined conformational diversity. We searched for global behaviour patterns studying how structure-based features change among the available conformer population for each protein. This procedure allowed us to describe the RMSD distribution in terms of three main protein classes sharing given properties. The largest of these protein subsets (~60%), which we call "rigid" (average RMSD = 0.83 Å), has no disordered regions, shows low conformational diversity, the largest tunnels and smaller and buried cavities. The two additional subsets contain disordered regions, but with differential sequence composition and behaviour. Partially disordered proteins have on average 67% of their conformers with disordered regions, average RMSD = 1.1 Å, the highest number of hinges and the longest disordered regions. In contrast, malleable proteins have on average only 25% of disordered conformers and average RMSD = 1.3 Å, flexible cavities affected in size by the presence of disordered regions and show the highest diversity of cognate ligands. Proteins in each set are mostly non-homologous to each other, share no given fold class, nor functional similarity but do share features derived from their conformer population. These shared features could represent conformational mechanisms related with biological functions.

Protein function machinery: from basic structural units to modulation of activity

Contemporary protein structure is a result of the trade off between the laws of physics and the evolutionary selection. The polymer nature of proteins played a decisive role in establishing the basic structural and functional units of soluble proteins. We discuss how these elementary building blocks work in the hierarchy of protein domain structure, co-translational folding, as well as in enzymatic activity and molecular interactions. Next, we consider modulators of the protein function, such as intermolecular interactions, disorder-to-order transitions, and allosteric signaling, acting via interference with the protein's structural dynamics. We also discuss the post-translational modifications, which is a complementary intricate mechanism evolved for regulation of protein functions and interactions. In conclusion, we assess an anticipated contribution of discussed topics to the future advancements in the field.

A structural perspective of RNA recognition by intrinsically disordered proteins

Protein-RNA recognition is essential for gene expression and its regulation, which is indispensable for the survival of the living organism at one hand, on the other hand, misregulation of this recognition may lead to their extinction. Polymorphic conformation of both the interacting partners is a characteristic feature of such molecular recognition that promotes the assembly. Many RNA binding proteins (RBP) or regions in them are found to be intrinsically disordered, and this property helps them to play a central role in the regulatory processes. Sequence composition and the length of the flexible linkers between RNA binding domains in RBPs are crucial in making significant contacts with its partner RNA. Polymorphic conformations of RBPs can provide thermodynamic advantage to its binding partner while acting as a chaperone. Prolonged extensions of the disordered regions in RBPs also contribute to the stability of the large cellular machines including ribosome and viral assemblies. The involvement of these disordered regions in most of the significant cellular processes makes RBPs highly associated with various human diseases that arise due to their misregulation.

In-silico structural analysis of E509K mutation in LARGE and T192M mutation in Alpha Dystroglycan in the inhibition of glycosylation of Alpha Dystroglycan by LARGE

Impaired glycosylation of cellular receptor Alpha Dystroglycan (α-DG) leads to dystroglycanopathy. Glycoprotein α-DG is the receptor protein in the Dystrophin Associated Protein Complex (DAPC), a macromolecular gathering on muscle cell membrane to form a bridge between extracellular matrix (ECM) and cellular actin cytoskeleton. Proper glycosylation of α-DG is mediated by the glycosylating enzyme LARGE. Mutations either in α-DG or in LARGE lead to improper glycosylations of α-DG thereby hampering the formation of final Laminin binding form α-DG resulting in dystroglycanopathy. In our current work, we explored the structural changes associated with the presence of mutations in α-DG as well as in the enzyme LARGE. We further extended our research to understand the effect of the mutations onto protein-enzyme interactions. Moreover, since LARGE transfers the sugar moiety (glucuronic acid; GlcA) onto α-DG, we tried to analyze what effect the mutation in LARGE confers on this enzyme ligand interaction. This work for the first time addressed the molecular changes occurring in the structures α-DG, LARGE and their interactions and shed lights on the as yet poorly understood mechanism behind the dystroglycanopathy onset.

Destabilizing an interacting motif strengthens the association of a designed ankyrin repeat protein with tubulin

Affinity maturation by random mutagenesis and selection is an established technique to make binding molecules more suitable for applications in biomedical research, diagnostics and therapy. Here we identified an unexpected novel mechanism of affinity increase upon in vitro evolution of a tubulin-specific designed ankyrin repeat protein (DARPin). Structural analysis indicated that in the progenitor DARPin the C-terminal capping repeat (C-cap) undergoes a 25° rotation to avoid a clash with tubulin upon binding. Additionally, the C-cap appears to be involved in electrostatic repulsion with tubulin. Biochemical and structural characterizations demonstrated that the evolved mutants achieved a gain in affinity through destabilization of the C-cap, which relieves the need of a DARPin conformational change upon tubulin binding and removes unfavorable interactions in the complex. Therefore, this specific case of an order-to-disorder transition led to a 100-fold tighter complex with a subnanomolar equilibrium dissociation constant, remarkably associated with a 30% decrease of the binding surface.

Statistical analysis of structural determinants for protein-DNA-binding specificity

DNA-binding proteins play critical roles in biological processes including gene expression, DNA packaging and DNA repair. They bind to DNA target sequences with different degrees of binding specificity, ranging from highly specific (HS) to nonspecific (NS). Alterations of DNA-binding specificity, due to either genetic variation or somatic mutations, can lead to various diseases. In this study, a comparative analysis of protein-DNA complex structures was carried out to investigate the structural features that contribute to binding specificity. Protein-DNA complexes were grouped into three general classes based on degrees of binding specificity: HS, multispecific (MS), and NS. Our results show a clear trend of structural features among the three classes, including amino acid binding propensities, simple and complex hydrogen bonds, major/minor groove and base contacts, and DNA shape. We found that aspartate is enriched in HS DNA binding proteins and predominately binds to a cytosine through a single hydrogen bond or two consecutive cytosines through bidentate hydrogen bonds. Aromatic residues, histidine and tyrosine, are highly enriched in the HS and MS groups and may contribute to specific binding through different mechanisms. To further investigate the role of protein flexibility in specific protein-DNA recognition, we analyzed the conformational changes between the bound and unbound states of DNA-binding proteins and structural variations. The results indicate that HS and MS DNA-binding domains have larger conformational changes upon DNA-binding and larger degree of flexibility in both bound and unbound states. Proteins 2016; 84:1147-1161. © 2016 Wiley Periodicals, Inc.

A structural perspective of RNA recognition by intrinsically disordered proteins

Protein-RNA recognition is essential for gene expression and its regulation, which is indispensable for the survival of the living organism at one hand, on the other hand, misregulation of this recognition may lead to their extinction. Polymorphic conformation of both the interacting partners is a characteristic feature of such molecular recognition that promotes the assembly. Many RNA binding proteins (RBP) or regions in them are found to be intrinsically disordered, and this property helps them to play a central role in the regulatory processes. Sequence composition and the length of the flexible linkers between RNA binding domains in RBPs are crucial in making significant contacts with its partner RNA. Polymorphic conformations of RBPs can provide thermodynamic advantage to its binding partner while acting as a chaperone. Prolonged extensions of the disordered regions in RBPs also contribute to the stability of the large cellular machines including ribosome and viral assemblies. The involvement of these disordered regions in most of the significant cellular processes makes RBPs highly associated with various human diseases that arise due to their misregulation.

An approach to improve the resolution of helical filaments with a large axial rise and flexible subunits

Single particle analysis is widely used for three-dimensional reconstruction of helical filaments. Near-atomic resolution has been obtained for several well-ordered filaments. However, it is still a challenge to achieve high resolution for filaments with flexible subunits and a large axial rise per subunit relative to pixel size. Here, we describe an approach that improves the resolution in such cases. In filaments with a large axial rise, many segments must be shifted a long distance along the filament axis to match with a reference projection, potentially causing loss of alignment accuracy and hence resolution. In our study of myosin filaments, we overcame this problem by pre-determining the axial positions of myosin head crowns within segments to decrease the alignment error. In addition, homogeneous, well-ordered segments were selected from the raw data set by checking the assigned azimuthal rotation angle of segments in each filament against those expected for perfect helical symmetry. These procedures improved the resolution of the filament reconstruction from 30 Å to 13 Å. This approach could be useful in other helical filaments with a large axial rise and/or flexible subunits.

Multidimensional persistence in biomolecular data

Persistent homology has emerged as a popular technique for the topological simplification of big data, including biomolecular data. Multidimensional persistence bears considerable promise to bridge the gap between geometry and topology. However, its practical and robust construction has been a challenge. We introduce two families of multidimensional persistence, namely pseudomultidimensional persistence and multiscale multidimensional persistence. The former is generated via the repeated applications of persistent homology filtration to high-dimensional data, such as results from molecular dynamics or partial differential equations. The latter is constructed via isotropic and anisotropic scales that create new simiplicial complexes and associated topological spaces. The utility, robustness, and efficiency of the proposed topological methods are demonstrated via protein folding, protein flexibility analysis, the topological denoising of cryoelectron microscopy data, and the scale dependence of nanoparticles. Topological transition between partial folded and unfolded proteins has been observed in multidimensional persistence. The separation between noise topological signatures and molecular topological fingerprints is achieved by the Laplace-Beltrami flow. The multiscale multidimensional persistent homology reveals relative local features in Betti-0 invariants and the relatively global characteristics of Betti-1 and Betti-2 invariants.

Looking below the surface of nicotinic acetylcholine receptors

The amino acid sequences of nicotinic acetylcholine receptors (nAChRs) from diverse species can be compared across extracellular, transmembrane, and intracellular domains. The intracellular domains are most divergent among subtypes, yet relatively consistent among species. The diversity indicates that each nAChR subtype has a unique language for communication with its host cell. The conservation across species also suggests that the intracellular domains have defining functional roles for each subtype. Secondary structure prediction indicates two relatively conserved alpha helices within the intracellular domains of all nAChRs. Among all subtypes, the intracellular domain of α7 nAChR is one of the most well conserved, and α7 nAChRs have effects in non-neuronal cells independent of generating ion currents, making it likely that the α7 intracellular domain directly mediates signal transduction. There are potential phosphorylation and protein-binding sites in the α7 intracellular domain, which are conserved and may be the basis for α7-mediated signal transduction.

Identification of Atg3 as an intrinsically disordered polypeptide yields insights into the molecular dynamics of autophagy-related proteins in yeast

The mechanism of autophagy relies on complex cell signaling and regulatory processes. Each cell contains many proteins that lack a rigid 3-dimensional structure under physiological conditions. These dynamic proteins, called intrinsically disordered proteins (IDPs) and protein regions (IDPRs), are predominantly involved in cell signaling and regulation. Yet, very little is known about their presence among proteins of the core autophagy machinery. In this work, we characterized the autophagy protein Atg3 from yeast and human along with 2 variants to show that Atg3 is an IDPRs-containing protein and that disorder/order predicted for these proteins from their amino acid sequence corresponds to their experimental characteristics. Based on this consensus, we applied the same prediction methods to all known Atg proteins from Saccharomyces cerevisiae. The data presented here provide an insight into the structural dynamics of each Atg protein. They also show that intrinsic disorder at various levels has to be taken into consideration for about half of the Atg proteins. This work should become a useful tool that will facilitate and encourage exploration of protein intrinsic disorder in autophagy.

The statistical conformation of a highly flexible protein: small-angle X-ray scattering of S. aureus protein A

Staphylococcal protein A (SpA) is a multidomain protein consisting of five globular IgG binding domains separated by a conserved six- to nine-residue flexible linker. We collected SAXS data on the N-terminal protein-binding half of SpA (SpA-N) and constructs consisting of one to five domain modules in order to determine statistical conformation of this important S. aureus virulence factor. We fit the SAXS data to a scattering function based on a new polymer physics model, which provides an analytical description of the SpA-N statistical conformation. We describe a protocol for systematically determining the appropriate level of modeling to fit a SAXS data set based on goodness of fit and whether the addition of parameters improves it. In the case of SpA-N, the analytical polymer physics description provides a depiction of the statistical conformation of a flexible protein that, while lacking atomistic detail, properly reflects the information content of the data.

Investigation of arc repressor DNA-binding specificity by comparative molecular dynamics simulations

Transcription factors regulate gene expression through binding to specific DNA sequences. How transcription factors achieve high binding specificity is still not well understood. In this paper, we investigated the role of protein flexibility in protein-DNA-binding specificity by comparative molecular dynamics (MD) simulations. Protein flexibility has been considered as a key factor in molecular recognition, which is intrinsically a dynamic process involving fine structural fitting between binding components. In this study, we performed comparative MD simulations on wild-type and F10V mutant P22 Arc repressor in both free and complex conformations. The F10V mutant has lower DNA-binding specificity though both the bound and unbound main-chain structures between the wild-type and F10V mutant Arc are highly similar. We found that the DNA-binding motif of wild-type Arc is structurally more flexible than the F10V mutant in the unbound state, especially for the six DNA base-contacting residues in each dimer. We demonstrated that the flexible side chains of wild-type Arc lead to a higher DNA-binding specificity through forming more hydrogen bonds with DNA bases upon binding. Our simulations also showed a possible conformational selection mechanism for Arc-DNA binding. These results indicate the important roles of protein flexibility and dynamic properties in protein-DNA-binding specificity.

Structural basis of protein-protein interactions

Regulated interactions between proteins govern signaling pathways within and between cells. Structural studies on protein complexes formed reversibly and/or transiently illustrate the remarkable diversity of interactions, both in terms of interfacial size and nature. In recent years, "domain-peptide" interactions have gained much greater recognition and may be viewed as both pre-translational and posttranslational-dependent functional switches. Our understanding of the multistep regulation of auto-inhibited multidomain proteins has also grown. Their activity may be understood as the "combinatorial" output of multiple input signals, including phosphorylation, location, and mechanical force. The prospects for bridging the gap between the new "systems biology" data and the traditional "reductionist" data are also discussed.

Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species

Background

The identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures.

Results

The current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified.

Conclusions

The results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.

Bovine serum albumin with glycated carboxyl groups shows membrane-perturbing activities

The aim of the present study aimed to investigate whether glycated bovine serum albumin (BSA) showed novel activities on the lipid-water interface. Mannosylated BSA (Man-BSA) was prepared by modification of the carboxyl groups with p-aminophenyl α-d-mannopyranoside. In contrast to BSA, Man-BSA notably induced membrane permeability of egg yolk phosphatidylcholine (EYPC)/egg yolk sphingomyelin (EYSM)/cholesterol (Chol) and EYPC/EYSM vesicles. Noticeably, Man-BSA induced the fusion of EYPC/EYSM/Chol vesicles, but not of EYPC/EYSM vesicles. Although BSA and Man-BSA showed similar binding affinity for lipid vesicles, the lipid-bound conformation of Man-BSA was distinct from that of BSA. Moreover, Man-BSA adopted distinct structure upon binding with the EYPC/EYSM/Chol and EYPC/EYSM vesicles. Man-BSA could induce the fusion of EYPC/EYSM/Chol vesicles with K562 and MCF-7 cells, while Man-BSA greatly induced the leakage of Chol-depleted K562 and MCF-7 cells. The modified BSA prepared by conjugating carboxyl groups with p-aminophenyl α-d-glucopyranoside also showed membrane-perturbing activities. Collectively, our data indicate that conjugation of carboxyl groups with monosaccharide generates functional BSA with membrane-perturbing activities on the lipid-water interface.

Mechanism of Tau-promoted microtubule assembly as probed by NMR spectroscopy

Determining the molecular mechanism of the neuronal Tau protein in the tubulin heterodimer assembly has been a challenge owing to the dynamic character of the complex and the large size of microtubules. We use here defined constructs comprising one or two tubulin heterodimers to characterize their association with a functional fragment of Tau, named TauF4. TauF4 binds with high affinities to the tubulin heterodimer complexes, but NMR spectroscopy shows that it remains highly dynamic, partly because of the interaction with the acidic C-terminal tails of the tubulin monomers. When bound to a single tubulin heterodimer, TauF4 is characterized by an overhanging peptide corresponding to the first of the four microtubule binding repeats of Tau. This peptide becomes immobilized in the complex with two longitudinally associated tubulin heterodimers. The longitudinal associations are favored by the fragment and contribute to Tau's functional role in microtubule assembly.

Protein flexibility facilitates quaternary structure assembly and evolution

The flexibility of individual proteins aids their evolutionary recruitment into complexes with increasing numbers of distinct subunits.

The intrinsic flexibility of proteins allows them to undergo large conformational fluctuations in solution or upon interaction with other molecules. Proteins also commonly assemble into complexes with diverse quaternary structure arrangements. Here we investigate how the flexibility of individual protein chains influences the assembly and evolution of protein complexes. We find that flexibility appears to be particularly conducive to the formation of heterologous (i.e., asymmetric) intersubunit interfaces. This leads to a strong association between subunit flexibility and homomeric complexes with cyclic and asymmetric quaternary structure topologies. Similarly, we also observe that the more nonhomologous subunits that assemble together within a complex, the more flexible those subunits tend to be. Importantly, these findings suggest that subunit flexibility should be closely related to the evolutionary history of a complex. We confirm this by showing that evolutionarily more recent subunits are generally more flexible than evolutionarily older subunits. Finally, we investigate the very different explorations of quaternary structure space that have occurred in different evolutionary lineages. In particular, the increased flexibility of eukaryotic proteins appears to enable the assembly of heteromeric complexes with more unique components.

Proteins often interact with other proteins and assemble into complexes. Here we show that the flexibility of individual proteins is important for their recruitment to complexes, as it facilitates the formation of asymmetric interfaces between different subunits. The role of flexibility becomes increasingly important as a greater number of distinct proteins are packed together within a single complex: the more distinct subunits, the more flexible those subunits need to be. A consequence of this is that, when a protein complex gains a new subunit during evolution, the newer subunit will tend to be more flexible than the older subunits. This suggests that we may be able to partially reconstruct the evolutionary history of a protein complex by considering the flexibility of its subunits. We also find that the types of protein complexes an organism forms are closely related to the flexibility of its proteins, with eukaryotic species, and particularly animals, using their increased flexibility to assemble complexes involving more distinct components.

Protein modeling: what happened to the "protein structure gap"?

Computational modeling of three-dimensional macromolecular structures and complexes from their sequence has been a long-standing vision in structural biology. Over the last 2 decades, a paradigm shift has occurred: starting from a large "structure knowledge gap" between the huge number of protein sequences and small number of known structures, today, some form of structural information, either experimental or template-based models, is available for the majority of amino acids encoded by common model organism genomes. With the scientific focus of interest moving toward larger macromolecular complexes and dynamic networks of interactions, the integration of computational modeling methods with low-resolution experimental techniques allows the study of large and complex molecular machines. One of the open challenges for computational modeling and prediction techniques is to convey the underlying assumptions, as well as the expected accuracy and structural variability of a specific model, which is crucial to understanding its limitations.

Toward in silico biomolecular manipulation through static modes: atomic scale characterization of HIV-1 protease flexibility

Probing biomolecular flexibility with atomic-scale resolution is a challenging task in current computational biology for fundamental understanding and prediction of biomolecular interactions and associated functions. This paper makes use of the static mode method to study HIV-1 protease considered as a model system to investigate the full biomolecular flexibility at the atomic scale, the screening of active site biomechanical properties, the blind prediction of allosteric sites, and the design of multisite strategies to target deformations of interest. Relying on this single calculation run of static modes, we demonstrate that in silico predictive design of an infinite set of complex excitation fields is reachable, thanks to the storage of the static modes in a data bank that can be used to mimic single or multiatom contact and efficiently anticipate conformational changes arising from external stimuli. All along this article, we compare our results to data previously published and propose a guideline for efficient, predictive, and custom in silico experiments.

Intrinsically disordered proteins and intrinsically disordered protein regions

Intrinsically disordered proteins (IDPs) and IDP regions fail to form a stable structure, yet they exhibit biological activities. Their mobile flexibility and structural instability are encoded by their amino acid sequences. They recognize proteins, nucleic acids, and other types of partners; they accelerate interactions and chemical reactions between bound partners; and they help accommodate posttranslational modifications, alternative splicing, protein fusions, and insertions or deletions. Overall, IDP-associated biological activities complement those of structured proteins. Recently, there has been an explosion of studies on IDP regions and their functions, yet the discovery and investigation of these proteins have a long, mostly ignored history. Along with recent discoveries, we present several early examples and the mechanisms by which IDPs contribute to function, which we hope will encourage comprehensive discussion of IDPs and IDP regions in biochemistry textbooks. Finally, we propose future directions for IDP research.

Conserved RNA helicase FRH acts nonenzymatically to support the intrinsically disordered neurospora clock protein FRQ

Protein conformation dictates a great deal of protein function. A class of naturally unstructured proteins, termed intrinsically disordered proteins (IDPs), demonstrates that flexibility in structure can be as important mechanistically as rigid structure. At the core of the circadian transcription/translation feedback loop in Neurospora crassa is the protein FREQUENCY (FRQ), shown here shown to share many characteristics of IDPs. FRQ in turn binds to FREQUENCY-Interacting RNA Helicase (FRH), whose clock function has been assumed to relate to its predicted helicase function. However, mutational analyses reveal that the helicase function of FRH is not essential for the clock, and a region of FRH distinct from the helicase region is essential for stabilizing FRQ against rapid degradation via a pathway distinct from its typical ubiquitin-mediated turnover. These data lead to the hypothesis that FRQ is an IDP and that FRH acts nonenzymatically, stabilizing FRQ to enable proper clock circuitry/function.

Modeling protein association mechanisms and kinetics

Substantial advances have been made in modeling protein association mechanisms and in calculating association rate constants (ka). We now have a clear understanding of the physical factors underlying the wide range of experimental ka values. Half of the association problem, where ka is limited by diffusion, is perhaps solved, and for the other half, where conformational changes become rate-limiting, a number of promising methods are being developed for ka calculations. Notably, the binding kinetics of disordered proteins are receiving growing attention, with 'dock-and-coalesce' emerging as a general mechanism. Progress too has been made in the modeling of protein association kinetics under conditions mimicking the heterogeneous, crowded environments of cells, an endeavor that should ultimately lead to a better understanding of cellular functions.