MDavocado: Analysis and Visualization of Protein Motion by Time-Dependent Angular Diagrams.

Extracting meaningful information from atomistic molecular dynamics (MD) simulations of proteins remains a challenging task due to the high-dimensionality and complexity of the data. MD simulations yield trajectories that contain the positions of thousands of atoms in millions of steps. Gaining a comprehensive understanding of local dynamical events across the entire trajectory is often difficult.

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

Signalling through TNFR1 modulates proinflammatory gene transcription and programmed cell death, and its impairment causes autoimmune diseases and cancer. NEDD4-binding protein 1 (N4BP1) is a critical suppressor of proinflammatory cytokine production that acts as a regulator of innate immune signalling and inflammation. However, our current understanding about the molecular properties that enable N4BP1 to exert its suppressive potential remain limited.

PyCoM: a python library for large-scale analysis of residue-residue coevolution data.

Motivation: Computational methods to detect correlated amino acid positions in proteins have become a valuable tool to predict intra- and inter-residue protein contacts, protein structures, and effects of mutation on protein stability and function. While there are many tools and webservers to compute coevolution scoring matrices, there is no central repository of alignments and coevolution matrices for large-scale studies and pattern detection leveraging on biological and structural annotations already available in UniProt.

Results: We present a Python library, PyCoM, which enables users to query and analyze coevolution matrices and sequence alignments of 457 622 proteins, selected from UniProtKB/Swiss-Prot database (length ≤ 500 residues), from a precompiled coevolution matrix database (PyCoMdb).

MDSubSampler: a posteriori sampling of important protein conformations from biomolecular simulations.

Motivation: Molecular dynamics (MD) simulations have become routine tools for the study of protein dynamics and function. Thanks to faster GPU-based algorithms, atomistic and coarse-grained simulations are being used to explore biological functions over the microsecond timescale, yielding terabytes of data spanning multiple trajectories, thereby extracting relevant protein conformations without losing important information is often challenging.

Results: We present MDSubSampler, a Python library and toolkit for a posteriori subsampling of data from multiple trajectories.

Recognizing the Binding Pattern and Dissociation Pathways of the p300 Taz2-p53 TAD2 Complex.

The dynamic association and dissociation between proteins are the basis of cellular signal transduction. This process becomes much more complicated if one or both interaction partners are intrinsically disordered because intrinsically disordered proteins can undergo disorder-to-order transitions upon binding to their partners. p53, a transcription factor with disordered regions, plays significant roles in many cellular signaling pathways.

PathDetect-SOM: A Neural Network Approach for the Identification of Pathways in Ligand Binding Simulations.

Understanding the process of ligand-protein recognition is important to unveil biological mechanisms and to guide drug discovery and design. Enhanced-sampling molecular dynamics is now routinely used to simulate the ligand binding process, resulting in the need for suitable tools for the analysis of large data sets of binding events. Here, we designed, implemented, and tested PathDetect-SOM, a tool based on self-organizing maps to build concise visual models of the ligand binding pathways sampled along single simulations or replicas.

Dynamic Profiling of β-Coronavirus 3CL M Protease Ligand-Binding Sites.

β-coronavirus (CoVs) alone has been responsible for three major global outbreaks in the 21st century. The current crisis has led to an urgent requirement to develop therapeutics. Even though a number of vaccines are available, alternative strategies targeting essential viral components are required as a backup against the emergence of lethal viral variants.

Machine Learning Prediction of Allosteric Drug Activity from Molecular Dynamics.

Allosteric drugs have been attracting increasing interest over the past few years. In this context, it is common practice to use high-throughput screening for the discovery of non-natural allosteric drugs. While the discovery stage is supported by a growing amount of biological information and increasing computing power, major challenges still remain in selecting allosteric ligands and predicting their effect on the target protein's function.

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

Several experimental studies indicated that large conformational changes, including partial domain unfolding, have a role in the functional mechanisms of the basic helix loop helix Per/ARNT/SIM (bHLH-PAS) transcription factors. Recently, single-molecule atomic force microscopy (AFM) revealed two distinct pathways for the mechanical unfolding of the ARNT PAS-B. In this work we used steered molecular dynamics simulations to gain new insights into this process at an atomistic level.

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). Constitutively active D13K-Y106W CheY has been an important tool for motor physiology.

Structure-based enzyme engineering improves donor-substrate recognition of Arabidopsis thaliana glycosyltransferases.

Glycosylation of secondary metabolites involves plant UDP-dependent glycosyltransferases (UGTs). UGTs have shown promise as catalysts in the synthesis of glycosides for medical treatment. However, limited understanding at the molecular level due to insufficient biochemical and structural information has hindered potential applications of most of these UGTs.

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

is a Gram-negative bacterium that is able to replicate within a broad range of aquatic protozoan hosts. 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 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.

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.

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.

In silico identification of rescue sites by double force scanning.

Motivation: A deleterious amino acid change in a protein can be compensated by a second-site rescue mutation. These compensatory mechanisms can be mimicked by drugs. In particular, the location of rescue mutations can be used to identify protein regions that can be targeted by small molecules to reactivate a damaged mutant.

Binding of Myomesin to Obscurin-Like-1 at the Muscle M-Band Provides a Strategy for Isoform-Specific Mechanical Protection.

The sarcomeric cytoskeleton is a network of modular proteins that integrate mechanical and signaling roles. Obscurin, or its homolog obscurin-like-1, bridges the giant ruler titin and the myosin crosslinker myomesin at the M-band. Yet, the molecular mechanisms underlying the physical obscurin(-like-1):myomesin connection, important for mechanical integrity of the M-band, remained elusive.

The Gearbox of the Bacterial Flagellar Motor Switch.

Switching of flagellar motor rotation sense dictates bacterial chemotaxis. Multi-subunit FliM-FliG rotor rings couple signal protein binding in FliM with reversal of a distant FliG C-terminal (FliGC) helix involved in stator contacts. Subunit dynamics were examined in conformer ensembles generated by molecular simulations from the X-ray structures.

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.

Coevolved Mutations Reveal Distinct Architectures for Two Core Proteins in the Bacterial Flagellar Motor.

Switching of bacterial flagellar rotation is caused by large domain movements of the FliG protein triggered by binding of the signal protein CheY to FliM. FliG and FliM form adjacent multi-subunit arrays within the basal body C-ring. The movements alter the interaction of the FliG C-terminal (FliGC) "torque" helix with the stator complexes.

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).

Bridging topological and functional information in protein interaction networks by short loops profiling.

Protein-protein interaction networks (PPINs) have been employed to identify potential novel interconnections between proteins as well as crucial cellular functions. In this study we identify fundamental principles of PPIN topologies by analysing network motifs of short loops, which are small cyclic interactions of between 3 and 6 proteins. We compared 30 PPINs with corresponding randomised null models and examined the occurrence of common biological functions in loops extracted from a cross-validated high-confidence dataset of 622 human protein complexes.

The crystal structure of the human titin:obscurin complex reveals a conserved yet specific muscle M-band zipper module.

M10 is the most C-terminal immunoglobulin (Ig) domain of the giant protein titin and a frequent target of disease-linked mutations. Currently, it is the only known muscle Ig domain able to interact with two alternative ligands-obscurin and obscurin-like-1 (Obsl1)-in different sarcomeric subregions. Obscurin and Obsl1 use their homologous N-terminal Ig domain (O1 in obscurin and OL1 in Obsl1) to bind M10 in a mutually exclusive manner.

BCR-ABL residues interacting with ponatinib are critical to preserve the tumorigenic potential of the oncoprotein.

Patients with chronic myeloid leukemia in whom tyrosine kinase inhibitors (TKIs) fail often present mutations in the BCR-ABL catalytic domain. We noticed a lack of substitutions involving 4 amino acids (E286, M318, I360, and D381) that form hydrogen bonds with ponatinib. We therefore introduced mutations in each of these residues, either preserving or altering their physicochemical properties.