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). PyCoM facilitates the development of statistical analyses of residue coevolution patterns using filters on biological and structural annotations from UniProtKB/Swiss-Prot, with simple access to PyCoMdb for both novice and advanced users, supporting Jupyter Notebooks, Python scripts, and a web API access. The resource is open source and will help in generating data-driven computational models and methods to study and understand protein structures, stability, function, and design. AVAILABILITY AND IMPLEMENTATION: PyCoM code is freely available from https://github.com/scdantu/pycom and PyCoMdb and the Jupyter Notebook tutorials are freely available from https://pycom.brunel.ac.uk.

Cleaving DNA with DNA: Cooperative Tuning of Structure and Reactivity Driven by Copper Ions.

A copper-dependent self-cleaving DNA (DNAzyme or deoyxyribozyme) previously isolated by in vitro selection has been analyzed by a combination of Molecular Dynamics (MD) simulations and advanced Electron Paramagnetic Resonance (Electron Spin Resonance) EPR/ESR spectroscopy, providing insights on the structural and mechanistic features of the cleavage reaction. The modeled 46-nucleotide deoxyribozyme in MD simulations forms duplex and triplex sub-structures that flank a highly conserved catalytic core. The DNA self-cleaving construct can also form a bimolecular complex that has a distinct substrate and enzyme domains. The highly dynamic structure combined with an oxidative site-specific cleavage of the substrate are two key-aspects to elucidate. By combining EPR/ESR spectroscopy with selectively isotopically labeled nucleotides it has been possible to overcome the major drawback related to the "metal-soup" scenario, also known as "super-stoichiometric" ratios of cofactors versus substrate, conventionally required for the DNA cleavage reaction within those nucleic acids-based enzymes. The focus on the endogenous paramagnetic center (Cu2+) here described paves the way for analysis on mixtures where several different cofactors are involved. Furthermore, the insertion of cleavage reaction within more complex architectures is now a realistic perspective towards the applicability of EPR/ESR spectroscopic studies.

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. This toolkit provides access to uniform, random, stratified, weighted sampling, and bootstrapping sampling methods. Sampling can be performed under the constraint of preserving the original distribution of relevant geometrical properties. Possible applications include simulations post-processing, noise reduction, and structures selection for ensemble docking. AVAILABILITY AND IMPLEMENTATION: MDSubSampler is freely available at https://github.com/alepandini/MDSubSampler, along with guidance on installation and tutorials on how it can be used.

Myricetin protects pancreatic ß-cells from human islet amyloid polypeptide (hIAPP) induced cytotoxicity and restores islet function.

The aberrant misfolding and self-assembly of human islet amyloid polypeptide (hIAPP)-a hormone that is co-secreted with insulin from pancreatic ß-cells-into toxic oligomers, protofibrils and fibrils has been observed in type 2 diabetes mellitus (T2DM). The formation of these insoluble aggregates has been linked with the death and dysfunction of ß-cells. Therefore, hIAPP aggregation has been identified as a therapeutic target for T2DM management. Several natural products are now being investigated for their potential to inhibit hIAPP aggregation and/or disaggregate preformed aggregates. In this study, we attempt to identify the anti-amyloidogenic potential of Myricetin (MYR)- a polyphenolic flavanoid, commonly found in fruits (like Syzygium cumini). Our results from biophysical studies indicated that MYR supplementation inhibits hIAPP aggregation and disaggregates preformed fibrils into non-toxic species. This protection was accompanied by inhibition of oxidative stress, reduction in lipid peroxidation and the associated membrane damage and restoration of mitochondrial membrane potential in INS-1E cells. MYR supplementation also reversed the loss of functionality in hIAPP exposed pancreatic islets via restoration of glucose-stimulated insulin secretion. Molecular dynamics simulation studies suggested that MYR molecules interact with the hIAPP pentameric fibril model at the amyloidogenic core region and thus prevents aggregation and distort the fibrils.

The 'hidden side' of spin labelled oligonucleotides: Molecular dynamics study focusing on the EPR-silent components of base pairing.

Nitroxide labels are combined with nucleic acid structures and are studied using electron paramagnetic resonance experiments (EPR). As X-ray/NMR structures are unavailable with the nitroxide labels, detailed residue level information, down to atomic resolution, about the effect of these nitroxide labels on local RNA structures is currently lacking. This information is critical to evaluate the choice of spin label. In this study, we compare and contrast the effect of TEMPO-based (NT) and rigid spin (Ç) labels (in both 2'-O methylated and not-methylated forms) on RNA duplexes. We also investigate sequence- dependent effects of NT label on RNA duplex along with the more complex G-quadruplex RNA. Distances measured from molecular dynamics simulations between the two spin labels are in agreement with the EPR experimental data. To understand the effect of labelled oligonucleotides on the structure, we studied the local base pair geometries and global structure in comparison with the unlabelled structures. Based on the structural analysis, we can conclude that TEMPO-based and Ç labels do not significantly perturb the base pair arrangements of the native oligonucleotide. When experimental structures for the spin labelled DNA/RNA molecules are not available, general framework offered by the current study can be used to provide information critical to the choice of spin labels to facilitate future EPR studies.

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. In the absence of experimental crystal structures, we employed advanced molecular modeling and simulations in conjunction with biochemical characterization to design a workflow to study five Group H Arabidopsis thaliana (76E1, 76E2, 76E4, 76E5, 76D1) UGTs. Based on our rational structural manipulation and analysis, we identified key amino acids (P129 in 76D1; D374 in 76E2; K275 in 76E4), which when mutated improved donor substrate recognition than wildtype UGTs. Molecular dynamics simulations and deep learning analysis identified structural differences, which drive substrate preferences. The design of these UGTs with broader substrate specificity may play important role in biotechnological and industrial applications. These findings can also serve as basis to study other plant UGTs and thereby advancing UGT enzyme engineering.

A trade-off for covalent and intercalation binding modes: a case study for Copper (II) ions and singly modified DNA nucleoside.

Selective binding to nucleic acids and, more generally, to biopolymers, very often requires at a minimum the presence of specific functionalities and precise spatial arrangement. DNA can fold into defined 3D structures upon binding to metal centers and/or lanthanides. Binding efficiency can be boosted by modified nucleosides incorporated into DNA sequences. In this work the high selectivity of modified nucleosides towards copper (II) ions, when used in the monomeric form, is unexpectedly and drastically reduced upon being covalently attached to the DNA sequence in single-site scenario. Surprisingly, such selectivity is partially retained upon non-covalent (i.e. intercalation) mixture formed by native DNA duplex and a nucleoside in the monomeric form. Exploiting the electron spin properties of such different and rich binding mode scenarios, 1D/2D pulsed EPR experiments have been used and tailored to differentiate among the different modes. An unusual correlation of dispersion of hyperfine couplings and strength of the binding mode(s) is described.

Conformational flexibility of histone variant CENP-ACse4 is regulated by histone H4: A mechanism to stabilize soluble Cse4.

The histone variant CENP-ACse4 is a core component of the specialized nucleosome at the centromere in budding yeast and is required for genomic integrity. Accordingly, the levels of Cse4 in cells are tightly regulated, primarily by ubiquitin-mediated proteolysis. However, structural transitions in Cse4 that regulate its centromeric localization and interaction with regulatory components are poorly understood. Using time-resolved fluorescence, NMR, and molecular dynamics simulations, we show here that soluble Cse4 can exist in a "closed" conformation, inaccessible to various regulatory components. We further determined that binding of its obligate partner, histone H4, alters the interdomain interaction within Cse4, enabling an "open" state that is susceptible to proteolysis. This dynamic model allows kinetochore formation only in the presence of H4, as the Cse4 N terminus, which is required for interaction with other centromeric components, is unavailable in the absence of H4. The specific requirement of H4 binding for the conformational regulation of Cse4 suggests a structure-based regulatory mechanism for Cse4 localization. Our data suggested a novel structural transition-based mechanism where conformational flexibility of the Cse4 N terminus can control Cse4 levels in the yeast cell and prevent Cse4 from interacting with kinetochore components at ectopic locations for formation of premature kinetochore assembly.

Conformational dynamics of Peb4 exhibit "mother's arms" chain model: a molecular dynamics study.

Peb4 from Campylobacter jejuni is an intertwined dimeric, periplasmic holdase, which also exhibits peptidyl prolyl cis/trans isomerase (PPIase) activity. Peb4 gene deletion alters the outer membrane protein profile and impairs cellular adhesion and biofilm formation for C. jejuni. Earlier crystallographic study has proposed that the PPIase domains are flexible and might form a cradle for holding the substrate and these aspects of Peb4 were explored using sub-microsecond molecular dynamics simulations in solution environment. Our simulations have revealed that PPIase domains are highly flexible and undergo a large structural change where they move apart from each other by 8 nm starting at .5 nm. Further, this large conformational change renders Peb4 as a compact protein with crossed-over conformation, forms a central cavity, which can "cradle" the target substrate. As reported for other chaperone proteins, flexibility of linker region connecting the chaperone and PPIase domains is key to forming the "crossed-over" conformation. The conformational transition of the Peb4 protein from the X-ray structure to the crossed-over conformation follows the "mother's arms" chain model proposed for the FkpA chaperone protein. Our results offer insights into how Peb4 and similar chaperones can use the conformational heterogeneity at their disposal to perform its much-revered biological function.

Structural insights into dynamics of RecU-HJ complex formation elucidates key role of NTR and stalk region toward formation of reactive state.

Holliday junction (HJ) resolving enzyme RecU is involved in DNA repair and recombination. We have determined the crystal structure of inactive mutant (D88N) of RecU from Bacillus subtilis in complex with a 12 base palindromic DNA fragment at a resolution of 3.2 Å. This structure shows the stalk region and the essential N-terminal region (NTR) previously unseen in our DNA unbound structure. The flexible nature of the NTR in solution was confirmed using SAXS. Thermofluor studies performed to assess the stability of RecU in complex with the arms of an HJ indicate that it confers stability. Further, we performed molecular dynamics (MD) simulations of wild type and an NTR deletion variant of RecU, with and without HJ. The NTR is observed to be highly flexible in simulations of the unbound RecU, in agreement with SAXS observations. These simulations revealed domain dynamics of RecU and their role in the formation of complex with HJ. The MD simulations also elucidate key roles of the NTR, stalk region, and breathing motion of RecU in the formation of the reactive state.

Mutation of Arg191 in FtsZ Impairs Cytokinetic Abscission of Bacillus subtilis Cells.

FtsZ monomers assemble to form a dynamic Z-ring at the midcell position in bacteria that coordinates bacterial cell division. Antibacterial agents plumbagin and SB-RA-2001 were found to bind to FtsZ and to inhibit Z-ring formation in bacteria. Docking analysis indicated similar binding regions for these two inhibitors on FtsZ, and residue R191 was involved in the binding interaction with both compounds. In this work, the importance of R191 in FtsZ assembly and in bacterial cell division was analyzed. R191A-FtsZ exhibited significantly poorer polymerization ability. Further, the mutant FtsZ could poison the assembly of wild-type FtsZ (WT-FtsZ). The expression of R191A-FtsZ in Bacillus subtilis strain PL2084 perturbed Z-ring formation and produced filamentous cells, indicating that the mutation hindered the division of these cells. The results suggested that the R191A mutation is a dominant negative mutation of FtsZ. Molecular dynamics simulations of R191A-FtsZ and WT-FtsZ revealed a kink in helices H5 and H7 in the active site of R191A-FtsZ compared to that of WT-FtsZ, which is required for FtsZ assembly. The findings suggested that R191 is an important residue for FtsZ assembly, which can be targeted for the design of FtsZ inhibitors.

Backbone and side chain assignments of human cell cycle regulatory protein S-phase kinase-associated protein 1.

Ubiquitination of proteins is required to regulate several cellular mechanisms in cells. Skp1-Cullin-1-F-box (SCF), the largest family of the RING E3 ligases, recognizes and carries out the poly-ubiquitination of many substrate proteins. SCF E3 ligase is a multi-component protein complex, and the human S-phase kinase-associated protein 1 (Skp1) is the adapter protein, which binds and presents the substrate binding protein F-box (FBP) to the rest of the E3 ligase. Several crystallographic studies have solved the partial structure of Skp1 in complex with various FBPs, but there is no structure of standalone Skp1. Understanding the conformational and structural properties of Skp1 with and without FBPs is required to understand the complete mechanism of poly-ubiquitination. Here, we report ~90 % backbone and 64 % side chain (1)H, (13)C, (15)N assignments of Skp1 protein using various double and triple resonance NMR experiments.