ClpS2 is a small protein under development as a probe for selectively recognizing N‐terminal amino acids of N‐degron peptide fragments. To understand the structural basis of ClpS2 specificity for an N‐terminal amino acid, all atom molecular dynamics (MD) simulations were conducted using the sequence of a bench‐stable mutant of ClpS2, called PROSS. We predicted that a single amino acid leucine to asparagine

Leucine rich repeats (LRRs) are present in over 430 000 proteins from viruses to eukaryotes. The LRRs are 20-30 residues long and occur in tandem. Individual LRRs are separated into a highly conserved segment with the consensus of LxxLxLxxNxL or LxxLxLxxNxxL (HCS) and a variable segment (VS). In LRRs parallel stacking of short β-strands (at positions 3-5 in HCS) form a super helix arrangement called

Conversion of 10‐hydroxygeraniol to 10‐oxogeranial is a crucial step in iridoid biosynthesis. This reaction is catalyzed by a zinc‐dependent alcohol dehydrogenase, 10‐hydroxygeraniol dehydrogenase, belonging to the family of medium‐chain dehydrogenase/reductase (MDR). Here, we report the crystal structures of a novel 10‐hydroxygeraniol dehydrogenase from Catharanthus roseus in its apo and nicotinamide

Integration of template‐based modeling, global sampling and precise scoring is crucial for the development of molecular docking programs with improved accuracy. We combined template‐based modeling and ab‐initio docking protocol as hybrid docking strategy called CoDock for the docking and scoring experiments of the seventh CAPRI edition. For CAPRI rounds 38‐45, we obtained acceptable or better models

Targets in the protein docking experiment CAPRI (Critical Assessment of Predicted Interactions) generally present new challenges and contribute to new developments in methodology. In rounds 38 to 45 of CAPRI, most targets could be effectively predicted using template‐based methods. However, the server ClusPro required structures rather than sequences as input, and hence we had to generate and dock

Protein docking is essential for structural characterization of protein interactions. Besides providing the structure of protein complexes, modeling of proteins and their complexes is important for understanding the fundamental principles and specific aspects of protein interactions. The accuracy of protein modeling, in general, is still less than that of the experimental approaches. Thus, it is important

Structural information about protein‐protein interactions, often missing at the interactome scale, is important for mechanistic understanding of cells and rational discovery of therapeutics. Protein docking provides a computational alternative for such information. However, ranking near‐native docked models high among a large number of candidates, often known as the scoring problem, remains a critical

There have been several studies suggesting that protein structures solved by NMR spectroscopy and x-ray crystallography show significant differences. To understand the origin of these differences, we assembled a database of high-quality protein structures solved by both methods. We also find significant differences between NMR and crystal structures-in the root-mean-square deviations of the Cα atomic

Internal structure similarity in proteins can be observed at the domain and at subdomain levels. From an evolutionary perspective, structurally similar elements may arise divergently by gene duplication and fusion events but may also be the product of convergent evolution under physicochemical constraints. The characterization of proteins that contain repeated structural elements has implications for

Kinesin dimer walks processively along a microtubule (MT) protofilament in a hand-over-hand manner, transiting alternately between one-head-bound (1HB) and two-heads-bound (2HB) states. In 1HB state, one head bound by adenosine diphosphate (ADP) is detached from MT and the other head is bound to MT. Here, using all-atom molecular dynamics simulations we determined the position and orientation of the

Cytochromes P450 are versatile heme-based enzymes responsible for vital life processes. Of these, P450cam (substrate camphor) has been most studied. Despite this, precise mechanisms of the key O─O cleavage step remain partly elusive to date; effects observed in various enzyme mutants remain partly unexplained. We have carried out extended (to 1000 ns) MM-MD and follow-on quantum mechanics/molecular

Musashi-2 (MSI2) belongs to Musashi family of RNA binding proteins (RBP). Like Musashi-1 (MSI1), it is overexpressed in a variety of cancers and is a promising therapeutic target. Both MSI proteins contain two N-terminal RNA recognition motifs and play roles in posttranscriptional regulation of target mRNAs. Previously, we have identified several inhibitors of MSI1, all of which bind to MSI2 as well

Hydrostatic pressure has a vital role in the biological adaptation of the piezophiles, organisms that live under high hydrostatic pressure. However, the mechanisms by which piezophiles are able to adapt their proteins to high hydrostatic pressure is not well understood. One proposed hypothesis is that the volume changes of unfolding (ΔVTot ) for proteins from piezophiles is distinct from those of nonpiezophilic

FimH is a bacterial adhesin protein located at the tip of Escherichia coli fimbria that functions to adhere bacteria to host cells. Thus, FimH is a critical factor in bacterial infections such as urinary tract infections and is of interest in drug development. It is also involved in vaccine development and as a model for understanding shear-enhanced catch bond cell adhesion. To date, over 60 structures

Ste24 enzymes, a family of eukaryotic integral membrane proteins, are zinc metalloproteases (ZMPs) originally characterized as "CAAX proteases" targeting prenylated substrates, including a-factor mating pheromone in yeast and prelamin A in humans. Recently, Ste24 was shown to also cleave nonprenylated substrates. Reduced activity of the human ortholog, HsSte24, is linked to multiple disease states

In this paper, using Word2vec, a widely-used natural language processing method, we demonstrate that protein domains may have a learnable implicit semantic "meaning" in the context of their functional contributions to the multi-domain proteins in which they are found. Word2vec is a group of models which can be used to produce semantically meaningful embeddings of words or tokens in a fixed-dimension

The analysis of amino acid coevolution has emerged as a practical method for protein structural modeling by providing structural contact information from alignments of amino acid sequences. In parallel, chemical cross-linking/mass spectrometry (XLMS) has gained attention as a universally applicable method for obtaining low-resolution distance constraints to model the quaternary arrangements of proteins

Immune checkpoint blockade of signaling pathways such as PD-1/PD-L1 has recently opened up a new avenue for highly efficient immunotherapeutic strategies to treat cancer. Since tumor microenvironments are characterized by lower pH (5.5-7.0), pH-dependent protein-ligand interactions can be exploited as efficient means to regulate drug affinity and specificity for a variety of malignancies. In this article

Nostoc sp. PCC 7120 are filamentous cyanobacteria capable of both oxygenic photosynthesis and nitrogen fixation, with the latter taking place in specialized cells known as heterocysts that terminally differentiate from vegetative cells under conditions of nitrogen starvation. Cyanobacteria have existed on earth for more than 2 billion years and are thought to be responsible for oxygenation of the earth's

The nondiscriminating aspartyl‐tRNA synthetase (ND‐AspRS), found in many archaea and bacteria, covalently attaches aspartic acid to tRNAAsp and tRNAAsn generating a correctly charged Asp‐tRNAAsp and an erroneous Asp‐tRNAAsn. This relaxed tRNA specificity is governed by interactions between the tRNA and the enzyme. In an effort to assess the contributions of the anticodon‐binding domain to tRNA specificity

We investigated the structure and Brownian rotational motion of the PEST region (201-268) from human c-Myc oncoprotein, whose overexpression/dysregulation is associated with various types of cancer. The 77-residue PEST fragment revealed a large Stokes radius (~3.1 nm) and CD spectrum highlighting abundance of disordered structure. Changes in structure/dynamics at two specific sites in PEST degron were

The 3D structure of a protein is essential to understand protein dynamics. If experimentally determined structure is unavailable, comparative models could be used to infer dynamics. However, the effectiveness of comparative models, compared to experimental structures, in inferring dynamics is not clear. To address this, we compared dynamics features of ~800 comparative models with their crystal structures

The multidrug resistance (MDR) system effectively expels antibiotics out of bacteria causing serious issues during bacterial infection. In addition to drug, indole, a common metabolic waste of bacteria, is expelled by MDR system of gram-negative bacteria for their survival. Experimental results suggest that AcrB, one of the key components of MDR system, undergoes large scale conformation changes during

We report docking performance on the six targets of Critical Assessment of PRedicted Interactions (CAPRI) rounds 39-45 that involved heteromeric protein-protein interactions and had the solved structures released since the rounds were held. Our general strategy involved protein-protein docking using ZDOCK, reranking using IRAD, and structural refinement using Rosetta. In addition, we made extensive

Comparative docking is based on experimentally determined structures of protein-protein complexes (templates), following the paradigm that proteins with similar sequences and/or structures form similar complexes. Modeling utilizing structure similarity of target monomers to template complexes significantly expands structural coverage of the interactome. Template-based docking by structure alignment

Salt-bridges play a unique role in the structural and functional stability of proteins, especially under harsh environments. How these salt-bridges contribute to the overall thermodynamic stability of protein structure and function across different domains of life is elusive still date. To address the issue, statistical analyses on the energies of salt-bridges, involved in proteins' structure and function

Protein-protein docking plays an important role in the computational prediction of the complex structure between two proteins. For years, a variety of docking algorithms have been developed, as witnessed by the critical assessment of prediction interactions (CAPRI) experiments. However, despite their successes, many docking algorithms often require a series of manual operations like modeling structures

Proteins play important roles in living organisms, and their function is directly linked with their structure. Due to the growing gap between the number of proteins being discovered and their functional characterization (in particular as a result of experimental limitations), reliable prediction of protein function through computational means has become crucial. This paper reviews the machine learning

Bromodomains (BrDs), a conserved structural module in chromatin-associated proteins, are well known for recognizing ε-N-acetyl lysine residues on histones. One of the most relevant BrDs is BRD4, a tandem BrD containing protein (BrD1 and BrD2) that plays a critical role in numerous diseases including cancer. Growing evidence shows that the two BrDs of BRD4 have different biological functions; hence

Macromolecules are characterized by distinctive arrangement of hydrogen bonds. Different patterns of hydrogen bonds give rise to distinct and stable structural motifs. An analysis of 4114 non-redundant protein chains reveals the existence of a three-residue, (i - 1) to (i + 1), structural motif, having two hydrogen-bonded five-membered pseudo rings (the first, an NH···OC involving the first residue

Lysosomal acid lipase (LAL) plays an important role in lipid metabolism by performing hydrolysis of triglycerides and cholesteryl esters in the lysosome. Based upon characteristics of LAL purified from human liver, it has been proposed that LAL is a proprotein with a 55 residue propeptide that may be essential for proper folding, intracellular transport, or enzymatic function. However, the biological

Small ubiquitin-related modifiers (SUMO1 and SUMO2) are ubiquitin family proteins, structurally similar to ubiquitin, differing in terms of their amino acid sequence and functions. Therefore, they provide a great platform for investigating sequence-structure-stability-function relationship. Here, we used chemical denaturation in comparing the folding-unfolding pathways of the SUMO proteins with their

Protein engineering and synthetic biology stand to benefit immensely from recent advances in silico tools for structural and functional analyses of proteins. In the context of designing novel proteins, current in silico tools inform the user on individual parameters of a query protein, with output scores/metrics unique to each parameter. In reality, proteins feature multiple "parts"/functions and modification

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not

The cross-strand disulfides (CSDs) found in β-hairpin antimicrobial peptides (β-AMPs) show a unique disulfide geometry that is characterized by unusual torsion angles and a short Cα-Cα distance. While the sequence and disulfide bond connectivity of disulfide-rich peptides is well studied, much less is known about the disulfide geometry found in CSDs and their role in the stability of β-AMPs. To address

Recognition of antigenic peptides bound to major histocompatibility complex (MHC) proteins by αβ T cell receptors (TCRs) is a hallmark of T cell mediated immunity. Recent data suggest that variations in TCR binding geometry may influence T cell signaling, which could help explain outliers in relationships between physical parameters such as TCR-pMHC binding affinity and T cell function. Traditionally

Smoothened (SMO) antagonist Vismodegib effectively inhibits the Hedgehog pathway in proliferating cancer cells. In early stage of treatment, Vismodegib exhibited promising outcomes to regress the tumors cells, but ultimately relapsed due to the drug resistive mutations in SMO mostly occurring before (primary mutations G497W) or after (acquired mutations D473H/Y) anti-SMO therapy. This study investigates

The selectivity filter (SF) of bacterial voltage-gated sodium channels consists of four glutamate residues arranged in a C4 symmetry. The protonation state population of this tetrad is unclear. To address this question, we simulate the pore domain of bacterial voltage-gated sodium channel of Magnetococcus sp. (Nav Ms) through constant pH methodology in explicit solvent and free energy perturbation

The relative orientation of the two variable domains, VH and VL , influences the shape of the antigen binding site, that is, the paratope, and is essential to understand antigen specificity. ABangle characterizes the VH -VL orientation by using five angles and a distance and compares it to other known structures. Molecular dynamics simulations of antibody variable domains (Fvs) reveal fluctuations

We present the seventh report on the performance of methods for predicting the atomic resolution structures of protein complexes offered as targets to the community-wide initiative on the Critical Assessment of Predicted Interactions. Performance was evaluated on the basis of 36 114 models of protein complexes submitted by 57 groups-including 13 automatic servers-in prediction rounds held during the

Our information-driven docking approach HADDOCK has demonstrated a sustained performance since the start of its participation to CAPRI. This is due, in part, to its ability to integrate data into the modeling process, and to the robustness of its scoring function. We participated in CAPRI both as server and manual predictors. In CAPRI rounds 38-45, we have used various strategies depending on the available

Peptide-protein docking is challenging due to the considerable conformational freedom of the peptide. CAPRI rounds 38-45 included two peptide-protein interactions, both characterized by a peptide forming an additional beta strand of a beta sheet in the receptor. Using the Rosetta FlexPepDock peptide docking protocol we generated top-performing, high-accuracy models for targets 134 and 135, involving

Designing protein sequences that fold to a given three-dimensional (3D) structure has long been a challenging problem in computational structural biology with significant theoretical and practical implications. In this study, we first formulated this problem as predicting the residue type given the 3D structural environment around the C α atom of a residue, which is repeated for each residue of a protein

Hsp16.3, a molecular chaperone, plays a vital role in the growth and survival of Mycobacterium tuberculosis inside the host. We previously reported that deletion of three amino acid residues (142 STN144 ) from C-terminal extension (CTE) of Hsp16.3 triggers its structural perturbation and increases its chaperone activity, which reaches its apex upon the deletion of its entire CTE (141 RSTN144 ). Thus

Perdeuteration in neutron crystallography is an effective method for determining the positions of hydrogen atoms in proteins. However, there is shortage of evidence that the high-resolution details of perdeuterated proteins are consistent with those of the nondeuterated proteins. In this study, we determined the X-ray structure of perdeuterated high-potential iron-sulfur protein (HiPIP) at a high resolution

Understanding the reaction mechanism of CRISPR-associated protein 9 (Cas9) is crucial for the application of programmable gene editing. Despite the availability of the structures of Cas9 in apo- and substrate-bound forms, the catalytically active structure is still unclear. Our first attempt to explore the catalytic mechanism of Cas9 HNH domain has been based on the reasonable assumption that we are

Tropomyosin (Tpm) is an extended α-helical coiled-coil homodimer that regulates actinomyosin interactions in muscle. Molecular simulations of four Tpms, two from the vertebrate class Mammalia (rat and pig), and two from the invertebrate class Malacostraca (shrimp and lobster), showed that despite extensive sequence and structural homology across metazoans, dynamic behavior-particularly long-range structural

The concept of consensus in multiple sequence alignments (MSAs) has been used to design and engineer proteins previously with some success. However, consensus design implicitly assumes that all amino acid positions function independently, whereas in reality, the amino acids in a protein interact with each other and work cooperatively to produce the optimum structure required for its function. Correlation

Protein phosphorylation is one of the essential posttranslation modifications playing a vital role in the regulation of many fundamental cellular processes. We propose a LightGBM-based computational approach that uses evolutionary, geometric, sequence environment, and amino acid-specific features to decipher phosphate binding sites from a protein sequence. Our method, while compared with other existing

Incorporating the dynamic nature of biomolecules in the modeling of their complexes is a challenge, especially when the extent and direction of the conformational changes taking place upon binding is unknown. Estimating whether the binding of a biomolecule to its partner(s) occurs in a conformational state accessible to its unbound form ("conformational selection") and/or the binding process induces

The backbone of a protein is typically represented as either a C α -polyline, a three-dimensional (3D) polyline that passes through the C α atoms, or a tuple of ϕ,ψ pairs while its fold is usually assigned using the 3D topological arrangement of the secondary structure elements (SSEs). It is tricky to obtain the SSE composition for a protein from the C α -polyline representation while its 3D SSE arrangement

Single-stranded DNA-binding proteins (SSBs) are required for all known DNA metabolic events such as DNA replication, recombination and repair. While a wealth of structural and functional data is available on the essential human SSB, hSSB1 (NABP2/OBFC2B), the close homolog hSSB2 (NABP1/OBFC2A) remains relatively uncharacterized. Both SSBs possess a well-structured OB (oligonucleotide/oligosaccharide-binding)

G-protein coupled glucagon receptors (GCGRs) play an important role in glucose homeostasis and pathophysiology of Type-II Diabetes Mellitus (T2DM). The allosteric pocket located at the trans-membrane domain of GCGR consists of hydrophobic (TM5) and hydrophilic (TM7) units. Hydrophobic interactions with the amino acid residues present at TM5, found to facilitate the favorable orientation of antagonist

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple

Several properties of amino acid sequences corresponding to proteins that are known to fold are compared to those of randomly generated sequences and to sequences of intrinsically disordered proteins in order to find properties that distinguish folding sequences from the rest. The properties studied included helix and sheet propensities from secondary structure prediction, adjacency correlations, directionality

This work explores how phosphorylation of an unstructured protein region in inhibitor-2 (I2) regulates protein phosphatase-1 (PP1) enzyme activity using molecular dynamics (MD). Free I2 is largely unstructured; however, when bound to PP1, three segments adopt a stable structure. In particular, an I2 helix (i-helix) blocks the PP1 active site and inhibits phosphatase activity. I2 phosphorylation in

The assembly of protein actin into double-helical filaments promotes many eukaryotic cellular processes that are regulated by actin-binding proteins (ABPs). Actin filaments can adopt multiple conformations, known as structural polymorphism, which possibly influences the interaction between filaments and ABPs. Gelsolin is a Ca2+ -regulated ABP that severs and caps actin filaments. Gelsolin binding modulates

Describing the whole story of protein folding is currently the main enigmatic problem in molecular bioinformatics study. Protein folding mechanisms have been intensively investigated with experimental as well as simulation techniques. Since a protein folds into its specific 3D structure from a unique amino acid sequence, it is interesting to extract as much information as possible from the amino acid