Motivation: Biologically important proteins are often large, multidomain proteins, which are difficult to characterize by high-throughput experimental methods. Efficient domain/boundary predictions are thus increasingly required in diverse area of proteomics research for computationally dissecting proteins into readily analyzable domains. Results: We constructed a support vector machine (SVM)-based domain linker predictor, DROP (Domain linker pRediction using OPtimal features), which was trained with 25 optimal features. The optimal combination of features was identified from a set of 3000 features using a random forest algorithm complemented with a stepwise feature selection. DROP demonstrated a prediction sensitivity and precision of 41.3 and 49.4%, respectively. These values were over 19.9% higher than those of control SVM predictors trained with non-optimized features, strongly suggesting the efficiency of our feature selection method. In addition, the mean NDO-Score of DROP for predicting novel domains in seven CASP8 FM multidomain proteins was 0.760, which was higher than any of the 12 published CASP8 DP servers. Overall, these results indicate that the SVM prediction of domain linkers can be improved by identifying optimal features that best distinguish linker from non-linker regions. Availability: DROP is available at http://tuat.ac.jp/~domserv/DROP.html Contacts:toh-hiroyuki@aist.go.jp; ykuroda@cc.tuat.ac.jp Supplementary information: Supplementary data are available at Bioinformatics online.

NMR studies can provide unique information about protein conformations in solution. In CASP14, three reference structures provided by solution NMR methods were available (T1027, T1029, and T1055), as well as a fourth data set of NMR‐derived contacts for an integral membrane protein (T1088). For the three targets with NMR‐based structures, the best prediction results ranged from very good (GDT_TS = 0.90, for T1055) to poor (GDT_TS = 0.47, for T1029). We explored the basis of these results by comparing all CASP14 prediction models against experimental NMR data. For T1027, NMR data reveal extensive internal dynamics, presenting a unique challenge for protein structure prediction methods. The analysis of T1029 motivated exploration of a novel method of “inverse structure determination,” in which an AlphaFold2 model was used to guide NMR data analysis. NMR data provided to CASP predictor groups for target T1088, a 238‐residue integral membrane porin, was also used to assess several NMR‐assisted prediction methods. Most groups involved in this exercise generated similar beta‐barrel models, with good agreement with the experimental data. However, as was also observed in CASP13, some pure prediction groups that did not use any NMR data generated models for T1088 that better fit the NMR data than the models generated using these experimental data. These results demonstrate the remarkable power of modern methods to predict structures of proteins with accuracies rivaling solution NMR structures, and that it is now possible to reliably use prediction models to guide and complement experimental NMR data analysis.

Gaussia luciferase (GLuc) is a small luciferase (18.2 kDa; 168 residues) and is thus attracting much attention as a reporter protein, but the lack of structural information is hampering further application. Here, we report the first solution structure of a fully active, recombinant GLuc determined by heteronuclear multidimensional NMR. We obtained a natively folded GLuc by bacterial expression and efficient refolding using a Solubility Enhancement Petide (SEP) tag. Almost perfect assignments of GLuc’s ¹H, ¹³C and ¹⁵N backbone signals were obtained. GLuc structure was determined using CYANA, which automatically identified over 2500 NOEs of which > 570 were long-range. GLuc is an all-alpha-helix protein made of nine helices. The region spanning residues 10–18, 36–81, 96–145 and containing eight out of the nine helices was determined with a C_α-atom RMSD of 1.39 Å ± 0.39 Å. The structure of GLuc is novel and unique. Two homologous sequential repeats form two anti-parallel bundles made by 4 helices and tied together by three disulfide bonds. The N-terminal helix 1 is grabbed by these 4 helices. Further, we found a hydrophobic cavity where several residues responsible for bioluminescence were identified in previous mutational studies, and we thus hypothesize that this is a catalytic cavity, where the hydrophobic coelenterazine binds and the bioluminescence reaction takes place.

Protein stabilization is difficult to rationalize, but the detailed thermodynamic and structural analysis of a series of carefully designed mutants may provide experimental insights into the mechanisms underlying stabilization. Here, we report a systematic structural and thermodynamic analysis of bovine pancreatic trypsin inhibitor (BPTI) variants that are significantly stabilized through a single amino acid substitution at residue 38, which is located in a loop mostly exposed on the protein surface. Differential scanning calorimetry indicated that the BPTI‐[5,55]Gly¹⁴ variants with a single mutation at position 38 were stabilized in an enthalpy‐driven manner and that the magnitude of the stabilization increased as the hydrophobicity of residue 38 increased. This increase in the thermal stability of BPTI was unexpected because a hydrophobic residue on a protein surface is usually destabilizing. To identify the structural determinants of this stabilization, we determined the crystal structures of six BPTI‐[5,55]Gly¹⁴ variants (Gly¹⁴Gly³⁸, Gly¹⁴Ala³⁸, Gly¹⁴Val³⁸, Gly¹⁴Leu³⁸, Gly¹⁴Ile³⁸, and Gly¹⁴Lys³⁸) at high resolutions and showed that they retain essentially the same structure as the wild‐type BPTI. A more detailed examination of their structures indicated that the extent of thermal stabilization correlated with both improved local packing and increased hydration around the substitution sites. In particular, the number of water molecules near residue 38 increased upon mutation to a hydrophobic residue suggesting that improved hydration contributed to the enthalpy‐driven stabilization. Increasing a protein's thermal stability by the placement of a hydrophobic amino acid on the protein surface is a novel and unexpected phenomenon, and its exact nature is worth further examination, as it may provide a generic method for stabilizing proteins in an enthalpy‐driven manner.Database: The coordinates and structure factors of Gly¹⁴Gly³⁸, Gly¹⁴Ile³⁸, Gly¹⁴Leu³⁸, and Gly¹⁴Lys³⁸ variants of BPTI‐[5,55] are deposited in the Protein Data Bank under the PDB entry codes 5XX3, 5XX5, 5XX2, and 5XX4, respectively. We previously reported the structures of Gly¹⁴Ala³⁸ (2ZJX) and Gly¹⁴Val³⁸ (2ZVX).