Structures and mechanisms of actin ATP hydrolysis.

The major cytoskeleton protein actin undergoes cyclic transitions between the monomeric G-form and the filamentous F-form, which drive organelle transport and cell motility. This mechanical work is driven by the ATPase activity at the catalytic site in the F-form. For deeper understanding of the actin cellular functions, the reaction mechanism must be elucidated. Here, we show that a single actin molecule is trapped in the F-form by fragmin domain-1 binding and present their crystal structures in the ATP analog-, ADP-Pi-, and ADP-bound forms, at 1.15-Å resolutions. The G-to-F conformational transition shifts the side chains of Gln137 and His161, which relocate four water molecules including W1 (attacking water) and W2 (helping water) to facilitate the hydrolysis. By applying quantum mechanics/molecular mechanics calculations to the structures, we have revealed a consistent and comprehensive reaction path of ATP hydrolysis by the F-form actin. The reaction path consists of four steps: 1) W1 and W2 rotations; 2) PG-O3B bond cleavage; 3) four concomitant events: W1-PO3- formation, OH- and proton cleavage, nucleophilic attack by the OH- against PG, and the abstracted proton transfer; and 4) proton relocation that stabilizes the ADP-Pi-bound F-form actin. The mechanism explains the slow rate of ATP hydrolysis by actin and the irreversibility of the hydrolysis reaction. While the catalytic strategy of actin ATP hydrolysis is essentially the same as those of motor proteins like myosin, the process after the hydrolysis is distinct and discussed in terms of Pi release, F-form destabilization, and global conformational changes.

Elastic network model reveals distinct flexibilities of capping proteins bound to CARMIL and twinfilin-tail.

Capping protein (CP) binds to the barbed end of an actin-filament and inhibits its elongation. CARMIL binds CP and dissociates it from the barbed end of the actin-filament. The binding of CARMIL peptide alters the flexibility of CP, which is considered to facilitate the dissociation. Twinfilin also binds to CP through its C-terminal tail. The complex structures of the CP/twinfilin-tail (TW-tail) peptide indicate that the binding sites of CARMIL and TW-tail overlap. However, TW-tail binding does not facilitate the dissociation of CP from the barbed end. We extensively investigated the flexibilities of CP in the CP/TW-tail or CP/CARMIL complexes using an elastic network model and concluded that TW-tail binding does not alter the flexibility of CP. Our extensive analysis also highlighted that the strong contacts of peptides with the two domains of CP, that is, the CP-L and CP-S domains, are key to changing the flexibilities of CP. CARMIL peptides can interact strongly with both of the domains, while TW-tail peptides exclusively interact with the CP-S domain because the binding site of TW-tail on CP relatively shifts to the CP-S domain compared with that of CP/CARMIL. This result supports our hypothesis that the dissociation of CP from the barbed end is regulated by the flexibility of CP.

MICAN-SQ: a sequential protein structure alignment program that is applicable to monomers and all types of oligomers.

Motivation: Protein structure alignment is a significant tool to understand evolutionary processes and physicochemical properties of proteins. Important targets of structure alignment are not only monomeric but also oligomeric proteins that sometimes include domain swapping or fusions. Although various protein structural alignment programs have been developed, no method is applicable to any protein pair regardless of the number of chain components and oligomeric states with retaining sequential restrictions: structurally equivalent regions must be aligned in the same order along protein sequences. Results: In this paper, we introduced a new sequential protein structural alignment algorithm MICAN-SQ, which is applicable to protein structures in all oligomeric states. In particular, MICAN-SQ allows the complicated structural alignments of proteins with domain swapping or fusion regions. To validate MICAN-SQ, alignment accuracies were evaluated using curated alignments of monomers and examples of domain swapping, and compared with those of pre-existing protein structural alignment programs. The results of this study show that MICAN-SQ has superior accuracy and robustness in comparison with previous programs and offers limited computational times. We also demonstrate that MICAN-SQ correctly aligns very large complexes and fused proteins. The present computations warrant the consideration of MICAN-SQ for studies of evolutionary and physicochemical properties of monomeric structures and all oligomer types. Availability and implementation: The MICAN program was implemented in C. The source code and executable file can be freely downloaded from http://www.tbp.cse.nagoya-u.ac.jp/MICAN/. Supplementary information: Supplementary data are available at Bioinformatics online.

Structural changes of homodimers in the PDB.

Protein complexes are involved in various biological phenomena. These complexes are intrinsically flexible, and structural changes are essential to their functions. To perform a large-scale automated analysis of the structural changes of complexes, we combined two original methods. An application, SCPC, compares two structures of protein complexes and decides the match of binding mode. Another application, Motion Tree, identifies rigid-body motions in various sizes and magnitude from the two structural complexes with the same binding mode. This approach was applied to all available homodimers in the Protein Data Bank (PDB). We defined two complex-specific motions: interface motion and subunit-spanning motion. In the former, each subunit of a complex constitutes a rigid body, and the relative movement between subunits occurs at the interface. In the latter, structural parts from distinct subunits constitute a rigid body, providing the relative movement spanning subunits. All structural changes were classified and examined. It was revealed that the complex-specific motions were common in the homodimers, detected in around 40% of families. The dimeric interfaces were likely to be small and flat for interface motion, while large and rugged for subunit-spanning motion. Interface motion was accompanied by a drastic change in contacts at the interface, while the change in the subunit-spanning motion was moderate. These results indicate that the interface properties of homodimers correlated with the type of complex-specific motion. The study demonstrates that the pipeline of SCPC and Motion Tree is useful for the massive analysis of structural change of protein complexes.

Discovery of Cryoprotective Activity in Human Genome-Derived Intrinsically Disordered Proteins.

Intrinsically disordered proteins (IDPs) are an emerging phenomenon. They may have a high degree of flexibility in their polypeptide chains, which lack a stable 3D structure. Although several biological functions of IDPs have been proposed, their general function is not known. The only finding related to their function is the genetically conserved YSK2 motif present in plant dehydrins. These proteins were shown to be IDPs with the YSK2 motif serving as a core region for the dehydrins' cryoprotective activity. Here we examined the cryoprotective activity of randomly selected IDPs toward the model enzyme lactate dehydrogenase (LDH). All five IDPs that were examined were in the range of 35-45 amino acid residues in length and were equally potent at a concentration of 50 µg/mL, whereas folded proteins, the PSD-95/Dlg/ZO-1 (PDZ) domain, and lysozymes had no potency. We further examined their cryoprotective activity toward glutathione S-transferase as an example of the other enzyme, and toward enhanced green fluorescent protein as a non-enzyme protein example. We further examined the lyophilization protective activity of the peptides toward LDH, which revealed that some IDPs showed a higher activity than that of bovine serum albumin (BSA). Based on these observations, we propose that cryoprotection is a general feature of IDPs. Our findings may become a clue to various industrial applications of IDPs in the future.

Comprehensive analysis of motions in molecular dynamics trajectories of the actin capping protein and its inhibitor complexes.

The actin capping protein (CP) binds to actin filaments to block further elongation. The capping activity is inhibited by proteins V-1 and CARMIL interacting with CP via steric and allosteric mechanisms, respectively. The crystal structures of free CP, CP/V-1, and CP/CARMIL complexes suggest that the binding of CARMIL alters the flexibility of CP rather than the overall structure of CP, and this is an allosteric inhibition mechanism. Here, we performed molecular dynamics (MD) simulations of CP in the free form, and in complex with CARMIL or V-1. The resulting trajectories were analyzed exhaustively using Motion Tree, which identifies various rigid-body motions ranging from small local motions to large domain motions. After enumerating all the motions, CP flexibilities with different ligands were characterized by a list of frequencies for 20 dominant rigid-body motions, some of which were not identified in previous studies. The comparative analysis highlights the influence of the binding of the CARMIL peptide to CP flexibility. In free CP and the CP/V-1 complex, domain motions around a large crevice between the N-stalk and the CP-S domain occur frequently. The CARMIL peptide binds the crevice and suppresses the motions effectively. In addition, the binding of the CARMIL peptide enhances and alters local motions around the pocket that participates in V-1 binding. These newly identified motions are likely to suppress the binding of V-1 to CP. The observed changes in CP motion provide insights that describe the mechanism of allosteric regulation by CARMIL through modulating CP flexibility. Proteins 2016; 84:948-956. © 2016 Wiley Periodicals, Inc.

Interface property responsible for effective interactions of protean segments: Intrinsically disordered regions that undergo disorder-to-order transitions upon binding.

Proteins that lack a well-defined conformation under native conditions are referred to as intrinsically disordered proteins. When interacting with partner proteins, short regions in disordered proteins can undergo disorder-to-order transitions upon binding; these regions are called protean segments (ProSs). It has been indicated that interactions of ProSs are effective: the number of contacts per residue of ProS interface is large. To reveal the properties of ProS interface that are responsible for the interaction efficiency, we classified the interface into core, rim and support, and analyzed them based on the relative accessible surface area (rASA). Despite the effective interactions, the ProS interface is mainly composed of rim residues, rather than core. The ProS rim is more effective than the rim of heterodimers, because the average rASAs of ProS rim, which is significantly large in the monomeric state, provides a large area to be used for the interactions. The amino acid composition of ProSs correlated well with those of heterodimers in both the core and rim. Therefore, the composition cannot explain why the rASAs of the ProS rim are large in the monomeric state. The balance between a small core and a large rim, and the large solvent exposure of the rim in the monomeric state, are the key to the disorder-to-order transition and the effective interactions of ProSs.

IDEAL in 2014 illustrates interaction networks composed of intrinsically disordered proteins and their binding partners.

IDEAL (Intrinsically Disordered proteins with Extensive Annotations and Literature, http://www.ideal.force.cs.is.nagoya-u.ac.jp/IDEAL/) is a collection of intrinsically disordered proteins (IDPs) that cannot adopt stable globular structures under physiological conditions. Since its previous publication in 2012, the number of entries in IDEAL has almost tripled (120 to 340). In addition to the increase in quantity, the quality of IDEAL has been significantly improved. The new IDEAL incorporates the interactions of IDPs and their binding partners more explicitly, and illustrates the protein-protein interaction (PPI) networks and the structures of protein complexes. Redundant experimental data are arranged based on the clustering of Protein Data Bank entries, and similar sequences with the same binding mode are grouped. As a result, the new IDEAL presents more concise and informative experimental data. Nuclear magnetic resonance (NMR) disorder is annotated in a systematic manner, by identifying the regions with large deviations among the NMR models. The ordered/disordered and new domain predictions by DICHOT are available, as well as the domain assignments by HMMER. Some examples of the PPI networks and the highly deviated regions derived from NMR models will be described, together with other advances. These enhancements will facilitate deeper understanding of IDPs, in terms of their flexibility, plasticity and promiscuity.

Hierarchical domain-motion analysis of conformational changes in sarcoplasmic reticulum Ca²âº-ATPase.

Sarco(endo)plasmic reticulum Ca(2+)-ATPase transports two Ca(2+) per ATP-hydrolyzed across biological membranes against a large concentration gradient by undergoing large conformational changes. Structural studies with X-ray crystallography revealed functional roles of coupled motions between the cytoplasmic domains and the transmembrane helices in individual reaction steps. Here, we employed "Motion Tree (MT)," a tree diagram that describes a conformational change between two structures, and applied it to representative Ca(2+) -ATPase structures. MT provides information of coupled rigid-body motions of the ATPase in individual reaction steps. Fourteen rigid structural units, "common rigid domains (CRDs)" are identified from seven MTs throughout the whole enzymatic reaction cycle. CRDs likely act as not only the structural units, but also the functional units. Some of the functional importance has been newly revealed by the analysis. Stability of each CRD is examined on the morphing trajectories that cover seven conformational transitions. We confirmed that the large conformational changes are realized by the motions only in the flexible regions that connect CRDs. The Ca(2+) -ATPase efficiently utilizes its intrinsic flexibility and rigidity to response different switches like ligand binding/dissociation or ATP hydrolysis. The analysis detects functional motions without extensive biological knowledge of experts, suggesting its general applicability to domain movements in other membrane proteins to deepen the understanding of protein structure and function.

A Method for Systematic Assessment of Intrinsically Disordered Protein Regions by NMR.

Intrinsically disordered proteins (IDPs) that lack stable conformations and are highly flexible have attracted the attention of biologists. Therefore, the development of a systematic method to identify polypeptide regions that are unstructured in solution is important. We have designed an "indirect/reflected" detection system for evaluating the physicochemical properties of IDPs using nuclear magnetic resonance (NMR). This approach employs a "chimeric membrane protein"-based method using the thermostable membrane protein PH0471. This protein contains two domains, a transmembrane helical region and a C-terminal OB (oligonucleotide/oligosaccharide binding)-fold domain (named NfeDC domain), connected by a flexible linker. NMR signals of the OB-fold domain of detergent-solubilized PH0471 are observed because of the flexibility of the linker region. In this study, the linker region was substituted with target IDPs. Fifty-three candidates were selected using the prediction tool POODLE and 35 expression vectors were constructed. Subsequently, we obtained 15N-labeled chimeric PH0471 proteins with 25 IDPs as linkers. The NMR spectra allowed us to classify IDPs into three categories: flexible, moderately flexible, and inflexible. The inflexible IDPs contain membrane-associating or aggregation-prone sequences. This is the first attempt to use an indirect/reflected NMR method to evaluate IDPs and can verify the predictions derived from our computational tools.

PSCDB: a database for protein structural change upon ligand binding.

Proteins are flexible molecules that undergo structural changes to function. The Protein Data Bank contains multiple entries for identical proteins determined under different conditions, e.g. with and without a ligand molecule, which provides important information for understanding the structural changes related to protein functions. We gathered 839 protein structural pairs of ligand-free and ligand-bound states from monomeric or homo-dimeric proteins, and constructed the Protein Structural Change DataBase (PSCDB). In the database, we focused on whether the motions were coupled with ligand binding. As a result, the protein structural changes were classified into seven classes, i.e. coupled domain motion (59 structural changes), independent domain motion (70), coupled local motion (125), independent local motion (135), burying ligand motion (104), no significant motion (311) and other type motion (35). PSCDB provides lists of each class. On each entry page, users can view detailed information about the motion, accompanied by a morphing animation of the structural changes. PSCDB is available at http://idp1.force.cs.is.nagoya-u.ac.jp/pscdb/.

IDEAL: Intrinsically Disordered proteins with Extensive Annotations and Literature.

IDEAL, Intrinsically Disordered proteins with Extensive Annotations and Literature (http://www.ideal.force.cs.is.nagoya-u.ac.jp/IDEAL/), is a collection of knowledge on experimentally verified intrinsically disordered proteins. IDEAL contains manual annotations by curators on intrinsically disordered regions, interaction regions to other molecules, post-translational modification sites, references and structural domain assignments. In particular, IDEAL explicitly describes protean segments that can be transformed from a disordered state to an ordered state. Since in most cases they can act as molecular recognition elements upon binding of partner proteins, IDEAL provides a data resource for functional regions of intrinsically disordered proteins. The information in IDEAL is provided on a user-friendly graphical view and in a computer-friendly XML format.

Dual-wield NTPases: A novel protein family mined from AlphaFold DB.

AlphaFold protein structure database (AlphaFold DB) archives a vast number of predicted models. We conducted systematic data mining against AlphaFold DB and discovered an uncharacterized P-loop NTPase family. The structure of the protein family was surprisingly novel, showing an atypical topology for P-loop NTPases, noticeable twofold symmetry, and two pairs of independent putative active sites. Our findings show that structural data mining is a powerful approach to identifying undiscovered protein families.

An assignment of intrinsically disordered regions of proteins based on NMR structures.

Intrinsically disordered proteins (IDPs) do not adopt stable three-dimensional structures in physiological conditions, yet these proteins play crucial roles in biological phenomena. In most cases, intrinsic disorder manifests itself in segments or domains of an IDP, called intrinsically disordered regions (IDRs), but fully disordered IDPs also exist. Although IDRs can be detected as missing residues in protein structures determined by X-ray crystallography, no protocol has been developed to identify IDRs from structures obtained by Nuclear Magnetic Resonance (NMR). Here, we propose a computational method to assign IDRs based on NMR structures. We compared missing residues of X-ray structures with residue-wise deviations of NMR structures for identical proteins, and derived a threshold deviation that gives the best correlation of ordered and disordered regions of both structures. The obtained threshold of 3.2Å was applied to proteins whose structures were only determined by NMR, and the resulting IDRs were analyzed and compared to those of X-ray structures with no NMR counterpart in terms of sequence length, IDR fraction, protein function, cellular location, and amino acid composition, all of which suggest distinct characteristics. The structural knowledge of IDPs is still inadequate compared with that of structured proteins. Our method can collect and utilize IDRs from structures determined by NMR, potentially enhancing the understanding of IDPs.

Substrate-shielding and hydrolytic reaction in hydrolases.

In general, transferases undergo large structural changes and sequester substrate molecules, to shield them from water. By contrast, hydrolases exhibit only small structural changes, and expose substrate molecules to water. However, some hydrolases deeply bury their substrates within the proteins. To clarify the relationship between substrate-shielding and enzymatic functions, we investigated 70 representative hydrolase structures, and examined the relative accessible surface areas of their substrates. As compared to the hydrolases employing the single displacement reaction, the hydrolases employing the double displacement reaction bury the substrate within the proteins. The exo hydrolases display significantly more substrate-shielding from water than the endo hydrolases. It suggests that the substrate-shielding is related to the chemical reaction mechanism of the hydrolases and the substrate specificity.

SCPC: a method to structurally compare protein complexes.

MOTIVATION: Protein-protein interactions play vital functional roles in various biological phenomena. Physical contacts between proteins have been revealed using experimental approaches that have solved the structures of protein complexes at atomic resolution. To examine the huge number of protein complexes available in the Protein Data Bank, an efficient automated method that compares protein complexes is required. RESULTS: We have developed Structural Comparison of Protein Complexes (SCPC), a novel method to structurally compare protein complexes. SCPC compares the spatial arrangements of subunits in a complex with those in another complex using secondary structure elements. Similar substructures are detected in two protein complexes and the similarity is scored. SCPC was applied to dimers, homo-oligomers and haemoglobins. SCPC properly estimated structural similarities between the dimers examined as well as an existing method, MM-align. Conserved substructures were detected in a homo-tetramer and a homo-hexamer composed of homologous proteins. Classification of quaternary structures of haemoglobins using SCPC was consistent with the conventional classification. The results demonstrate that SCPC is a valuable tool to investigate the structures of protein complexes. AVAILABILITY: SCPC is available at http://idp1.force.cs.is.nagoya-u.ac.jp/scpc/. CONTACT: rkoike@is.nagoya-u.ac.jp SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Two distinct mechanisms for actin capping protein regulation--steric and allosteric inhibition.

The actin capping protein (CP) tightly binds to the barbed end of actin filaments, thus playing a key role in actin-based lamellipodial dynamics. V-1 and CARMIL proteins directly bind to CP and inhibit the filament capping activity of CP. V-1 completely inhibits CP from interacting with the barbed end, whereas CARMIL proteins act on the barbed end-bound CP and facilitate its dissociation from the filament (called uncapping activity). Previous studies have revealed the striking functional differences between the two regulators. However, the molecular mechanisms describing how these proteins inhibit CP remains poorly understood. Here we present the crystal structures of CP complexed with V-1 and with peptides derived from the CP-binding motif of CARMIL proteins (CARMIL, CD2AP, and CKIP-1). V-1 directly interacts with the primary actin binding surface of CP, the C-terminal region of the alpha-subunit. Unexpectedly, the structures clearly revealed the conformational flexibility of CP, which can be attributed to a twisting movement between the two domains. CARMIL peptides in an extended conformation interact simultaneously with the two CP domains. In contrast to V-1, the peptides do not directly compete with the barbed end for the binding surface on CP. Biochemical assays revealed that the peptides suppress the interaction between CP and V-1, despite the two inhibitors not competing for the same binding site on CP. Furthermore, a computational analysis using the elastic network model indicates that the interaction of the peptides alters the intrinsic fluctuations of CP. Our results demonstrate that V-1 completely sequesters CP from the barbed end by simple steric hindrance. By contrast, CARMIL proteins allosterically inhibit CP, which appears to be a prerequisite for the uncapping activity. Our data suggest that CARMIL proteins down-regulate CP by affecting its conformational dynamics. This conceptually new mechanism of CP inhibition provides a structural basis for the regulation of the barbed end elongation in cells.

SAHG, a comprehensive database of predicted structures of all human proteins.

Most proteins from higher organisms are known to be multi-domain proteins and contain substantial numbers of intrinsically disordered (ID) regions. To analyse such protein sequences, those from human for instance, we developed a special protein-structure-prediction pipeline and accumulated the products in the Structure Atlas of Human Genome (SAHG) database at http://bird.cbrc.jp/sahg. With the pipeline, human proteins were examined by local alignment methods (BLAST, PSI-BLAST and Smith-Waterman profile-profile alignment), global-local alignment methods (FORTE) and prediction tools for ID regions (POODLE-S) and homology modeling (MODELLER). Conformational changes of protein models upon ligand-binding were predicted by simultaneous modeling using templates of apo and holo forms. When there were no suitable templates for holo forms and the apo models were accurate, we prepared holo models using prediction methods for ligand-binding (eF-seek) and conformational change (the elastic network model and the linear response theory). Models are displayed as animated images. As of July 2010, SAHG contains 42,581 protein-domain models in approximately 24,900 unique human protein sequences from the RefSeq database. Annotation of models with functional information and links to other databases such as EzCatDB, InterPro or HPRD are also provided to facilitate understanding the protein structure-function relationships.

Accidental interaction between PDZ domains and diclofenac revealed by NMR-assisted virtual screening.

In silico approaches have become indispensable for drug discovery as well as drug repositioning and adverse effect prediction. We have developed the eF-seek program to predict protein-ligand interactions based on the surface structure of proteins using a clique search algorithm. We have also developed a special protein structure prediction pipeline and accumulated predicted 3D models in the Structural Atlas of the Human Genome (SAHG) database. Using this database, genome-wide prediction of non-peptide ligands for proteins in the human genome was performed, and a subset of predicted interactions including 14 PDZ domains was then confirmed by NMR titration. Surprisingly, diclofenac, a non-steroidal anti-inflammatory drug, was found to be a non-peptide PDZ domain ligand, which bound to 5 of 15 tested PDZ domains. The critical residues for the PDZ-diclofenac interaction were also determined. Pharmacological implications of the accidental PDZ-diclofenac interaction are further discussed.

Computational Methods to Predict Intrinsically Disordered Regions and Functional Regions in Them.

Intrinsically disordered regions (IDRs) are protein regions that do not adopt fixed tertiary structures. Since these regions lack ordered three-dimensional structures, they should be excluded from the target portions of homology modeling. IDRs can be predicted from the amino acid sequences, because their amino acid compositions are different from that of the structured domains. This chapter provides a review of the prediction methods of IDRs and a case study of IDR prediction.