Dissecting protein architecture with communication blocks and communicating segment pairs

Cysteine hyperoxidation rewires communication pathways in the nucleosome and destabilizes the dyad

Gene activity is tightly controlled by reversible chemical modifications called epigenetic marks, which are of various types and modulate gene accessibility without affecting the DNA sequence. Despite an increasing body of evidence demonstrating the role of oxidative-type modifications of histones in gene expression regulation, there remains a complete absence of structural data at the atomistic level to understand the molecular mechanisms behind their regulatory action. Owing to μs time-scale MD simulations and protein communication networks analysis, we describe the impact of histone H3 hyperoxidation (i.e., S-sulfonylation) on the nucleosome core particle dynamics. Our results reveal the atomic-scale details of the intrinsic structural networks within the canonical histone core and their perturbation by hyperoxidation of the histone H3 C110. We show that this modification involves local rearrangements of the communication networks and destabilizes the dyad, and that one modification is enough to induce a maximal structural signature. Our results suggest that cysteine hyperoxidation in the nucleosome core particle might favor its disassembly.

Computational and biochemical analysis of type IV pilus dynamics and stability

Type IV pili (T4P) are distinctive dynamic filaments at the surface of many bacteria that can rapidly extend and retract and withstand strong forces. T4P are important virulence factors in many human pathogens, including Enterohemorrhagic Escherichia coli (EHEC). The structure of the EHEC T4P has been determined by integrating nuclear magnetic resonance (NMR) and cryo-electron microscopy data. To better understand pilus assembly, stability, and function, we performed a total of 108 ms all-atom molecular dynamics simulations of wild-type and mutant T4P. Extensive characterization of the conformational landscape of T4P in different conditions of temperature, pH, and ionic strength is complemented with targeted mutagenesis and biochemical analyses. Our simulations and NMR experiments reveal a conserved set of residues defining a calcium-binding site at the interface between three pilin subunits. Calcium binding enhances T4P stability ex vivo and in vitro, supporting the role of this binding site as a potential pocket for drug design.

Structural determination of Streptococcus pyogenes M1 protein interactions with human immunoglobulin G using integrative structural biology

Streptococcus pyogenes (Group A streptococcus; GAS) is an important human pathogen responsible for mild to severe, life-threatening infections. GAS expresses a wide range of virulence factors, including the M family proteins. The M proteins allow the bacteria to evade parts of the human immune defenses by triggering the formation of a dense coat of plasma proteins surrounding the bacteria, including IgGs. However, the molecular level details of the M1-IgG interaction have remained unclear. Here, we characterized the structure and dynamics of this interaction interface in human plasma on the surface of live bacteria using integrative structural biology, combining cross-linking mass spectrometry and molecular dynamics (MD) simulations. We show that the primary interaction is formed between the S-domain of M1 and the conserved IgG Fc-domain. In addition, we show evidence for a so far uncharacterized interaction between the A-domain and the IgG Fc-domain. Both these interactions mimic the protein G-IgG interface of group C and G streptococcus. These findings underline a conserved scavenging mechanism used by GAS surface proteins that block the IgG-receptor (FcγR) to inhibit phagocytic killing. We additionally show that we can capture Fab-bound IgGs in a complex background and identify XLs between the constant region of the Fab-domain and certain regions of the M1 protein engaged in the Fab-mediated binding. Our results elucidate the M1-IgG interaction network involved in inhibition of phagocytosis and reveal important M1 peptides that can be further investigated as future vaccine targets.

Streptococcus pyogenes is a human specific pathogen causing both mild and invasive infections. It employs sophisticated mechanisms to evade and circumvent parts of the host’s immune defenses, in part via its major surface associated virulence factor, the family of M proteins. Of these, the M1 protein is the most prevalent serotype. The M1 protein creates a dense coat-like structure with multiple host proteins on the bacterial surface to disguise itself from opsonizing antibodies. It specifically interacts in a non-immune way with human immunoglobulin G (IgG) Fc-domains to disarm their receptor binding site. The molecular level details of this interaction have not been characterized. Here, we describe these interactions from minimally perturbed samples of human plasma adsorbed onto living bacteria using an integrative structural biology approach including cross-linking mass spectrometry, molecular modeling, and molecular dynamics simulations. We identify two distinct M1-peptides that bind IgGs and reveal the stability of these interactions. We show that both peptides block the Fc-receptor binding sites through capturing IgGs via their Fc-domains. These results highlight the importance of describing novel pathogen-derived peptides mediating host immune evasion as potential vaccine targets in future studies.

Predicting substitutions to modulate disorder and stability in coiled-coils

Background

Coiled-coils are described as stable structural motifs, where two or more helices wind around each other. However, coiled-coils are associated with local mobility and intrinsic disorder. Intrinsically disordered regions in proteins are characterized by lack of stable secondary and tertiary structure under physiological conditions in vitro. They are increasingly recognized as important for protein function. However, characterizing their behaviour in solution and determining precisely the extent of disorder of a protein region remains challenging, both experimentally and computationally.

Results

In this work, we propose a computational framework to quantify the extent of disorder within a coiled-coil in solution and to help design substitutions modulating such disorder. Our method relies on the analysis of conformational ensembles generated by relatively short all-atom Molecular Dynamics (MD) simulations. We apply it to the phosphoprotein multimerisation domains (PMD) of Measles virus (MeV) and Nipah virus (NiV), both forming tetrameric left-handed coiled-coils. We show that our method can help quantify the extent of disorder of the C-terminus region of MeV and NiV PMDs from MD simulations of a few tens of nanoseconds, and without requiring an extensive exploration of the conformational space. Moreover, this study provided a conceptual framework for the rational design of substitutions aimed at modulating the stability of the coiled-coils. By assessing the impact of four substitutions known to destabilize coiled-coils, we derive a set of rules to control MeV PMD structural stability and cohesiveness. We therefore design two contrasting substitutions, one increasing the stability of the tetramer and the other increasing its flexibility.

Conclusions

Our method can be considered as a platform to reason about how to design substitutions aimed at regulating flexibility and stability.

Folding with a protein's native shortcut network

A complex network approach to protein folding is proposed, wherein a protein's contact map is reconceptualized as a network of shortcut edges, and folding is steered by a structural characteristic of this network. Shortcut networks are generated by a known message passing algorithm operating on protein residue networks. It is found that the shortcut networks of native structures (SCN0s) are relevant graph objects with which to study protein folding at a formal level. The logarithm form of their contact order (SCN0_lnCO) correlates significantly with folding rate of two-state and nontwo-state proteins. The clustering coefficient of SCN0s (C_SCN0 ) correlates significantly with folding rate, transition-state placement and stability of two-state folders. Reasonable folding pathways for several model proteins are produced when C_SCN0 is used to combine protein segments incrementally to form the native structure. The folding bias captured by C_SCN0 is detectable in non-native structures, as evidenced by Molecular Dynamics simulation generated configurations for the fast folding Villin-headpiece peptide. These results support the use of shortcut networks to investigate the role protein geometry plays in the folding of both small and large globular proteins, and have implications for the design of multibody interaction schemes in folding models. One facet of this geometry is the set of native shortcut triangles, whose attributes are found to be well-suited to identify dehydrated intraprotein areas in tight turns, or at the interface of different secondary structure elements.

Use of Computational Biochemistry for Elucidating Molecular Mechanisms of Nitric Oxide Synthase

Nitric oxide (NO) is an essential signaling molecule in the regulation of multiple cellular processes. It is endogenously synthesized by NO synthase (NOS) as the product of L-arginine oxidation to L-citrulline, requiring NADPH, molecular oxygen, and a pterin cofactor. Two NOS isoforms are constitutively present in cells, nNOS and eNOS, and a third is inducible (iNOS). Despite their biological relevance, the details of their complex structural features and reactivity mechanisms are still unclear. In this review, we summarized the contribution of computational biochemistry to research on NOS molecular mechanisms. We described in detail its use in studying aspects of structure, dynamics and reactivity. We also focus on the numerous outstanding questions in the field that could benefit from more extensive computational investigations.

Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications

Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.

"Infostery" analysis of short molecular dynamics simulations identifies highly sensitive residues and predicts deleterious mutations

Characterizing a protein mutational landscape is a very challenging problem in Biology. Many disease-associated mutations do not seem to produce any effect on the global shape nor motions of the protein. Here, we use relatively short all-atom biomolecular simulations to predict mutational outcomes and we quantitatively assess the predictions on several hundreds of mutants. We perform simulations of the wild type and 175 mutants of PSD95’s third PDZ domain in complex with its cognate ligand. By recording residue displacements correlations and interactions, we identify “communication pathways” and quantify them to predict the severity of the mutations. Moreover, we show that by exploiting simulations of the wild type, one can detect 80% of the positions highly sensitive to mutations with a precision of 89%. Importantly, our analysis describes the role of these positions in the inter-residue communication and dynamical architecture of the complex. We assess our approach on three different systems using data from deep mutational scanning experiments and high-throughput exome sequencing. We refer to our analysis as “infostery”, from “info” - information - and “steric” - arrangement of residues in space. We provide a fully automated tool, COMMA2 (www.lcqb.upmc.fr/COMMA2), that can be used to guide medicinal research by selecting important positions/mutations.

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

Proteins are highly dynamic entities attaining a myriad of different conformations. Protein side chains change their states during dynamics, causing clashes that are propagated at distal sites. A convenient formalism to analyze protein dynamics is based on network theory using Protein Structure Networks (PSNs). Despite their broad applicability, few efforts have been devoted to benchmarking PSN methods and to provide the community with best practices. In many applications, it is convenient to use the centers of mass of the side chains as nodes. It becomes thus critical to evaluate the minimal distance cutoff between the centers of mass which will provide stable network properties. Moreover, when the PSN is derived from a structural ensemble collected by molecular dynamics (MD), the impact of the MD force field has to be evaluated. We selected a dataset of proteins with different fold and size and assessed the two fundamental properties of the PSN, i.e. hubs and connected components. We identified an optimal cutoff of 5 Å that is robust to changes in the force field and the proteins. Our study builds solid foundations for the harmonization and standardization of the PSN approach.

Performance of machine-learning scoring functions in structure-based virtual screening

Classical scoring functions have reached a plateau in their performance in virtual screening and binding affinity prediction. Recently, machine-learning scoring functions trained on protein-ligand complexes have shown great promise in small tailored studies. They have also raised controversy, specifically concerning model overfitting and applicability to novel targets. Here we provide a new ready-to-use scoring function (RF-Score-VS) trained on 15 426 active and 893 897 inactive molecules docked to a set of 102 targets. We use the full DUD-E data sets along with three docking tools, five classical and three machine-learning scoring functions for model building and performance assessment. Our results show RF-Score-VS can substantially improve virtual screening performance: RF-Score-VS top 1% provides 55.6% hit rate, whereas that of Vina only 16.2% (for smaller percent the difference is even more encouraging: RF-Score-VS top 0.1% achieves 88.6% hit rate for 27.5% using Vina). In addition, RF-Score-VS provides much better prediction of measured binding affinity than Vina (Pearson correlation of 0.56 and −0.18, respectively). Lastly, we test RF-Score-VS on an independent test set from the DEKOIS benchmark and observed comparable results. We provide full data sets to facilitate further research in this area (http://github.com/oddt/rfscorevs) as well as ready-to-use RF-Score-VS (http://github.com/oddt/rfscorevs_binary).

The molecular mechanism of hPPARα activation

MD simulations were performed to explore the molecular mechanism of hPPARα activation. 11 key residues favouring binding ligands and the movements of helices and loops playing important roles in inducing the active conformation change of hPPARα were discovered.

The Molecular Mechanism of P2Y1 Receptor Activation

Human purinergic G protein-coupled receptor P2Y1 (P2Y1 R) is activated by adenosine 5'-diphosphate (ADP) to induce platelet activation and thereby serves as an important antithrombotic drug target. Crystal structures of P2Y1 R revealed that one ligand (MRS2500) binds to the extracellular vestibule of this GPCR, whereas another (BPTU) occupies the surface between transmembrane (TM) helices TM2 and TM3. We introduced a total of 20 μs all-atom long-timescale molecular dynamic (MD) simulations to inquire why two molecules in completely different locations both serve as antagonists while ADP activates the receptor. Our results indicate that BPTU acts as an antagonist by stabilizing extracellular helix bundles leading to an increase of the lipid order, whereas MRS2500 blocks signaling by occupying the ligand binding site. Both antagonists stabilize an ionic lock within the receptor. However, binding of ADP breaks this ionic lock, forming a continuous water channel that leads to P2Y1 R activation.