Insights into stereoselective ring formation in canonical strigolactone: Identification of a dirigent domain-containing enzyme catalyzing orobanchol synthesis.

Strigolactones (SLs) are plant apocarotenoids with diverse roles and structures. Canonical SLs, widespread and characterized by structural variations in their tricyclic lactone (ABC-ring), are classified into two types based on C-ring configurations. The steric C-ring configuration emerges during the BC-ring closure, downstream of the biosynthetic intermediate, carlactonoic acid (CLA). Most plants produce either type of canonical SLs stereoselectively, e.g., tomato (Solanum lycopersicum) yields orobanchol with an α-oriented C-ring. The mechanisms driving SL structural diversification are partially understood, with limited insight into functional implications. Furthermore, the exact molecular mechanism for the stereoselective BC-ring closure reaction is yet to be known. We identified an enzyme, the stereoselective BC-ring-forming factor (SRF), from the dirigent protein (DIR) family, specifically the DIR-f subfamily, whose biochemical function had not been characterized, making it a key enzyme in stereoselective canonical SL biosynthesis with the α-oriented C-ring. We first confirm the precise catalytic function of the tomato cytochrome P450 SlCYP722C, previously shown to be involved in orobanchol biosynthesis [T. Wakabayashi et al., Sci. Adv. 5, eaax9067 (2019)], to convert CLA to 18-oxocarlactonoic acid. We then show that SRF catalyzes the stereoselective BC-ring closure reaction of 18-oxocarlactonoic acid, forming orobanchol. Our methodology combines experimental and computational techniques, including SRF structure prediction and conducting molecular dynamics simulations, suggesting a catalytic mechanism based on the conrotatory 4π-electrocyclic reaction for the stereoselective BC-ring formation in orobanchol. This study sheds light on the molecular basis of how plants produce SLs with specific stereochemistry in a controlled manner.

Chemigenetic indicators based on synthetic chelators and green fluorescent protein.

Molecular fluorescent indicators are versatile tools for dynamic imaging of biological systems. We now report a class of indicators that are based on the chemigenetic combination of a synthetic ion-recognition motif and a protein-based fluorophore. Specifically, we have developed a calcium ion (Ca2+) indicator that is based on genetic insertion of circularly permuted green fluorescent protein into HaloTag protein self-labeled with a ligand containing the Ca2+ chelator 1,2-bis(2-aminophenoxy)ethane-N,N,N',N'-tetraacetic acid. We have demonstrated the versatility of this design by also developing a sodium ion (Na+) indicator using a crown-ether-containing ligand. This approach affords bright and sensitive ion indicators that can be applicable to cell imaging. This design can enable the development of chemigenetic indicators with ion or molecular specificities that have not been realized with fully protein-based indicators.

QM/MM Study of the Catalytic Mechanism and Substrate Specificity of the Aromatic Substrate C-Methyltransferase Fur6.

In the field of medical chemistry and other organic chemistry, introducing a methyl group into a designed position has been difficult to achieve. However, owing to the vigorous developments in the field of enzymology, methyltransferases are considered potential tools for addressing this problem. Within the methyltransferase family, Fur6 catalyzes the methylation of C3 of 1,2,4,5,7-pentahydroxynaphthalene (PHN) using S-adenosyl-l-methionine (SAM) as the methyl donor. Here, we report the catalytic mechanism and substrate specificity of Fur6 based on computational studies. Our molecular dynamics (MD) simulation studies reveal the reactive form of PHN and its interactions with the enzyme. Our hybrid quantum mechanics/molecular mechanics (QM/MM) calculations suggest the reaction pathway of the methyl transfer step in which the energy barrier is 8.6 kcal mol-1. Our free-energy calculations with a polarizable continuum model (PCM) indicate that the final deprotonation step of the methylated intermediate occurs after it is ejected into the water solvent from the active center pocket of Fur6. Additionally, our studies on the protonation states, the highest occupied molecular orbital (HOMOs), and the energy barriers of the methylation reaction for the analogs of PHN demonstrate the mechanism of the specificity to PHN. Our study provides valuable insights into Fur6 chemistry, contributing to a deeper understanding of molecular mechanisms and offering an opportunity to engineer the enzyme to achieve high yields of the desired product(s).

Identification of difructose dianhydride I synthase/hydrolase from an oral bacterium establishes a novel glycoside hydrolase family.

Fructooligosaccharides and their anhydrides are widely used as health-promoting foods and prebiotics. Various enzymes acting on ß-D-fructofuranosyl linkages of natural fructan polymers have been used to produce functional compounds. However, enzymes that hydrolyze and form α-D-fructofuranosyl linkages have been less studied. Here, we identified the BBDE_2040 gene product from Bifidobacterium dentium (α-D-fructofuranosidase and difructose dianhydride I synthase/hydrolase from Bifidobacterium dentium [αFFase1]) as an enzyme with α-D-fructofuranosidase and α-D-arabinofuranosidase activities and an anomer-retaining manner. αFFase1 is not homologous with any known enzymes, suggesting that it is a member of a novel glycoside hydrolase family. When caramelized fructose sugar was incubated with αFFase1, conversions of ß-D-Frup-(2→1)-α-D-Fruf to α-D-Fruf-1,2':2,1'-ß-D-Frup (diheterolevulosan II) and ß-D-Fruf-(2→1)-α-D-Fruf (inulobiose) to α-D-Fruf-1,2':2,1'-ß-D-Fruf (difructose dianhydride I [DFA I]) were observed. The reaction equilibrium between inulobiose and DFA I was biased toward the latter (1:9) to promote the intramolecular dehydrating condensation reaction. Thus, we named this enzyme DFA I synthase/hydrolase. The crystal structures of αFFase1 in complex with ß-D-Fruf and ß-D-Araf were determined at the resolutions of up to 1.76 Å. Modeling of a DFA I molecule in the active site and mutational analysis also identified critical residues for catalysis and substrate binding. The hexameric structure of αFFase1 revealed the connection of the catalytic pocket to a large internal cavity via a channel. Molecular dynamics analysis implied stable binding of DFA I and inulobiose to the active site with surrounding water molecules. Taken together, these results establish DFA I synthase/hydrolase as a member of a new glycoside hydrolase family (GH172).

Cryptic Oxidative Transamination of Hydroxynaphthoquinone in Natural Product Biosynthesis.

Pyridoxal 5'-phosphate (PLP)-dependent enzymes are a group of versatile enzymes that catalyze various reactions, but only a small number of them react with O2. Here, we report an unprecedented PLP-dependent enzyme, NphE, that catalyzes both transamination and two-electron oxidation using O2 as an oxidant. Our intensive analysis reveals that NphE transfers the l-glutamate-derived amine to 1,3,6,8-tetrahydroxynaphthalene-derived mompain to form 8-amino-flaviolin (8-AF) via a highly conjugated quinonoid intermediate that is reactive with O2. During the NphE reaction, O2 is reduced to yield H2O2. An integrated technique involving NphE structure prediction by AlphaFold v2.0 and molecular dynamics simulation suggested the O2-accessible cavity. Our in vivo results demonstrated that 8-AF is a genuine biosynthetic intermediate for the 1,3,6,8-tetrahydroxynaphthalene-derived meroterpenoid naphterpin without an amino group, which was supported by site-directed mutagenesis. This study clearly establishes the NphE reaction product 8-AF as a common intermediate with a cryptic amino group for the biosynthesis of terpenoid-polyketide hybrid natural products.

The α- and ß-Subunit Boundary at the Stem of the Mushroom-Like α3ß3-Type Oxygenase Component of Rieske Non-Heme Iron Oxygenases Is the Rieske-Type Ferredoxin-Binding Site.

Cumene dioxygenase (CumDO) is an initial enzyme in the cumene degradation pathway of Pseudomonas fluorescens IP01 and is a Rieske non-heme iron oxygenase (RO) that comprises two electron transfer components (reductase [CumDO-R] and Rieske-type ferredoxin [CumDO-F]) and one catalytic component (α3ß3-type oxygenase [CumDO-O]). Catalysis is triggered by electrons that are transferred from NAD(P)H to CumDO-O by CumDO-R and CumDO-F. To investigate the binding mode between CumDO-F and CumDO-O and to identify the key CumDO-O amino acid residues for binding, we simulated docking between the CumDO-O crystal structure and predicted model of CumDO-F and identified two potential binding sites: one is at the side-wise site and the other is at the top-wise site in mushroom-like CumDO-O. Then, we performed alanine mutagenesis of 16 surface amino acid residues at two potential binding sites. The results of reduction efficiency analyses using the purified components indicated that CumDO-F bound at the side-wise site of CumDO-O, and K117 of the α-subunit and R65 of the ß-subunit were critical for the interaction. Moreover, these two positively charged residues are well conserved in α3ß3-type oxygenase components of ROs whose electron donors are Rieske-type ferredoxins. Given that these residues were not conserved if the electron donors were different types of ferredoxins or reductases, the side-wise site of the mushroom-like structure is thought to be the common binding site between Rieske-type ferredoxin and α3ß3-type oxygenase components in ROs. IMPORTANCE We clarified the critical amino acid residues of the oxygenase component (Oxy) of Rieske non-heme iron oxygenase (RO) for binding with Rieske-type ferredoxin (Fd). Our results showed that Rieske-type Fd-binding site is commonly located at the stem (side-wise site) of the mushroom-like α3ß3 quaternary structure in many ROs. The resultant binding site was totally different from those reported at the top-wise site of the doughnut-like α3-type Oxy, although α3-type Oxys correspond to the cap (α3 subunit part) of the mushroom-like α3ß3-type Oxys. Critical amino acid residues detected in this study were not conserved if the electron donors of Oxys were different types of Fds or reductases. Altogether, we can suggest that unique binding modes between Oxys and electron donors have evolved, depending on the nature of the electron donors, despite Oxy molecules having shared α3ß3 quaternary structures.

Mutagenesis analysis of GMN motif in Arabidopsis thaliana Mg2+ transporter MRS2-1.

Magnesium is an important nutrient for plants, but much is still unknown about plant Mg2+ transporters. Combining with the structural prediction of AlphaFold2, we used mutagenesis and 28Mg uptake assay to study the highly conserved "GMN" motif of Arabidopsis thaliana MRS2-1 (AtMRS2-1) transporter. We demonstrated that the glycine and methionine in GMN motif are essential for AtMRS2-1 to transport Mg2+.

Structural and Molecular Basis of the Catalytic Mechanism of Geranyl Pyrophosphate C6-Methyltransferase: Creation of an Unprecedented Farnesyl Pyrophosphate C6-Methyltransferase.

Prenyl pyrophosphate methyltransferases enhance the structural diversity of terpenoids. However, the molecular basis of their catalytic mechanisms is poorly understood. In this study, using multiple strategies, we characterized a geranyl pyrophosphate (GPP) C6-methyltransferase, BezA. Biochemical analysis revealed that BezA requires Mg2+ and solely methylates GPP. The crystal structures of BezA and its complex with S-adenosyl homocysteine were solved at 2.10 and 2.56 Å, respectively. Further analyses using site-directed mutagenesis, molecular docking, molecular dynamics simulations, and quantum mechanics/molecular mechanics calculations revealed the molecular basis of the methylation reaction. Importantly, the function of E170 as a catalytic base to complete the methylation reaction was established. We also succeeded in switching the substrate specificity by introducing a W210A substitution, resulting in an unprecedented farnesyl pyrophosphate C6-methyltransferase.

EXPRORER: Rational Cosolvent Set Construction Method for Cosolvent Molecular Dynamics Using Large-Scale Computation.

Cosolvent molecular dynamics (CMD) simulations involve an MD simulation of a protein in the presence of explicit water molecules mixed with cosolvent molecules to perform hotspot detection, binding site identification, and binding energy estimation, while other existing methods (e.g., MixMD, SILCS, and MDmix) utilize small molecules that represent functional groups of compounds. However, the cosolvent selections employed in these methods differ and there are only a few cosolvents that are commonly used in these methods. In this study, we proposed a systematic method for constructing a set of cosolvents for drug discovery, termed the EXtended PRObes set construction by REpresentative Retrieval (EXPRORER). First, we extracted typical substructures from FDA-approved drugs, generated 138 cosolvent structures, and for each cosolvent molecule, we conducted CMD simulations to generate a spatial probability distribution map of cosolvent atoms (PMAP). Analyses of PMAP similarity revealed that a cosolvent pair with a PMAP similarity greater than 0.70-0.75 shared similar structural features. We present a method for the construction of a cosolvent subset that satisfies a similarity threshold for all cosolvents, and we tested the constructed sets for four proteins. To our knowledge, this is the first study to include a systematic proposal for cosolvent set construction, and thus, the EXPRORER cosolvents will provide deeper insights into ligand binding sites of various proteins.

Substrate recognition mechanism of tRNA-targeting ribonuclease, colicin D, and an insight into tRNA cleavage-mediated translation impairment.

Colicin D is a plasmid-encoded bacteriocin that specifically cleaves tRNAArg of sensitive Escherichia coli cells. E. coli has four isoaccepting tRNAArgs; the cleavage occurs at the 3' end of anticodon-loop, leading to translation impairment in the sensitive cells. tRNAs form a common L-shaped structure and have many conserved nucleotides that limit tRNA identity elements. How colicin D selects tRNAArgs from the tRNA pool of sensitive E. coli cells is therefore intriguing. Here, we reveal the recognition mechanism of colicin D via biochemical analyses as well as structural modelling. Colicin D recognizes tRNAArgICG, the most abundant species of E. coli tRNAArgs, at its anticodon-loop and D-arm, and selects it as the most preferred substrate by distinguishing its anticodon-loop sequence from that of others. It has been assumed that translation impairment is caused by a decrease in intact tRNA molecules due to cleavage. However, we found that intracellular levels of intact tRNAArgICG do not determine the viability of sensitive cells after such cleavage; rather, an accumulation of cleaved ones does. Cleaved tRNAArgICG dominant-negatively impairs translation in vitro. Moreover, we revealed that EF-Tu, which is required for the delivery of tRNAs, does not compete with colicin D for binding tRNAArgICG, which is consistent with our structural model. Finally, elevation of cleaved tRNAArgICG level decreases the viability of sensitive cells. These results suggest that cleaved tRNAArgICG transiently occupies ribosomal A-site in an EF-Tu-dependent manner, leading to translation impairment. The strategy should also be applicable to other tRNA-targeting RNases, as they, too, recognize anticodon-loops.Abbreviations: mnm5U: 5-methylaminomethyluridine; mcm5s2U: 5-methoxycarbonylmethyl-2-thiouridine.

Identification, characterization, and structural analyses of a fungal endo-ß-1,2-glucanase reveal a new glycoside hydrolase family.

endo-ß-1,2-Glucanase (SGL) is an enzyme that hydrolyzes ß-1,2-glucans, which play important physiological roles in some bacteria as a cyclic form. To date, no eukaryotic SGL has been identified. We purified an SGL from Talaromyces funiculosus (TfSGL), a soil fungus, to homogeneity and then cloned the complementary DNA encoding the enzyme. TfSGL shows no significant sequence similarity to any known glycoside hydrolase (GH) families, but shows significant similarity to certain eukaryotic proteins with unknown functions. The recombinant TfSGL (TfSGLr) specifically hydrolyzed linear and cyclic ß-1,2-glucans to sophorose (Glc-ß-1,2-Glc) as a main product. TfSGLr hydrolyzed reducing-end-modified ß-1,2-gluco-oligosaccharides to release a sophoroside with the modified moiety. These results indicate that TfSGL is an endo-type enzyme that preferably releases sophorose from the reducing end of substrates. Stereochemical analysis demonstrated that TfSGL is an inverting enzyme. The overall structure of TfSGLr includes an (α/α)6 toroid fold. The substrate-binding mode was revealed by the structure of a Michaelis complex of an inactive TfSGLr mutant with a ß-1,2-glucoheptasaccharide. Mutational analysis and action pattern analysis of ß-1,2-gluco-oligosaccharide derivatives revealed an unprecedented catalytic mechanism for substrate hydrolysis. Glu-262 (general acid) indirectly protonates the anomeric oxygen at subsite -1 via the 3-hydroxy group of the Glc moiety at subsite +2, and Asp-446 (general base) activates the nucleophilic water via another water. TfSGLr is apparently different from a GH144 SGL in the reaction and substrate recognition mechanism based on structural comparison. Overall, we propose that TfSGL and closely-related enzymes can be classified into a new family, GH162.

Understanding the Molecular Mechanism Underlying the High Catalytic Activity of p-Hydroxybenzoate Hydroxylase Mutants for Producing Gallic Acid.

p-Hydroxybenzoate hydroxylase (PHBH) is a flavoprotein monooxygenase that catalyzes the hydroxylation of p-hydroxybenzoate (p-OHB) to 3,4-dihydroxybenzoate (3,4-DOHB). PHBH can bind to other benzoate derivatives in addition to p-OHB; however, hydroxylation does not occur on 3,4-DOHB. Replacement of Tyr385 with Phe forms a mutant, which enables the production of 3,4,5-trihydroxybenzonate (gallic acid) from 3,4-DOHB, although the catalytic activity of the mutant is quite low. In this study, we report how the L199V/Y385F double mutant exhibits activity for producing gallic acid 4.3-fold higher than that of the Y385F single mutant. This improvement in catalytic activity is primarily due to the suppression of a shunt reaction that wastes reduced nicotinamide adenine dinucleotide phosphate by producing H2O2. To further elucidate the molecular mechanism underlying this higher catalytic activity, we performed molecular dynamics simulations and quantum mechanics/molecular mechanics calculations, in addition to determining the crystal structure of the Y385F·3,4-DOHB complex. The simulations showed that the Y385F mutation facilitates the deprotonation of the 4-hydroxy group of 3,4-DOHB, which is necessary for initiating hydroxylation. Moreover, the L199V mutation in addition to the Y385F mutation allows the OH moiety in the peroxide group of C-(4a)-flavin hydroperoxide to come into the proximity of the C5 atom of 3,4-DOHB. Overall, this study provides a consistent explanation for the change in the catalytic activity of PHBH caused by mutations, which will enable us to better design an enzyme with different activities.

Structural and thermodynamic insights into ß-1,2-glucooligosaccharide capture by a solute-binding protein in Listeria innocua.

ß-1,2-Glucans are bacterial carbohydrates that exist in cyclic or linear forms and play an important role in infections and symbioses involving Gram-negative bacteria. Although several ß-1,2-glucan-associated enzymes have been characterized, little is known about how ß-1,2-glucan and its shorter oligosaccharides (Sop n s) are captured and imported into the bacterial cell. Here, we report the biochemical and structural characteristics of the Sop n -binding protein (SO-BP, Lin1841) associated with the ATP-binding cassette (ABC) transporter from the Gram-positive bacterium Listeria innocua Calorimetric analysis revealed that SO-BP specifically binds to Sop n s with a degree of polymerization of 3 or more, with Kd values in the micromolar range. The crystal structures of SO-BP in an unliganded open form and in closed complexes with tri-, tetra-, and pentaoligosaccharides (Sop3-5) were determined to a maximum resolution of 1.6 Å. The binding site displayed shape complementarity to Sop n , which adopted a zigzag conformation. We noted that water-mediated hydrogen bonds and stacking interactions play a pivotal role in the recognition of Sop3-5 by SO-BP, consistent with its binding thermodynamics. Computational free-energy calculations and a mutational analysis confirmed that interactions with the third glucose moiety of Sop n s are significantly responsible for ligand binding. A reduction in unfavorable changes in binding entropy that were in proportion to the lengths of the Sop n s was explained by conformational entropy changes. Phylogenetic and sequence analyses indicated that SO-BP ABC transporter homologs, glycoside hydrolases, and other related proteins are co-localized in the genomes of several bacteria. This study may improve our understanding of bacterial ß-1,2-glucan metabolism and promote the discovery of unidentified ß-1,2-glucan-associated proteins.

Trisaccharide containing α2,3-linked sialic acid is a receptor for mumps virus.

Mumps virus (MuV) remains an important pathogen worldwide, causing epidemic parotitis, orchitis, meningitis, and encephalitis. Here we show that MuV preferentially uses a trisaccharide containing α2,3-linked sialic acid in unbranched sugar chains as a receptor. Crystal structures of the MuV attachment protein hemagglutinin-neuraminidase (MuV-HN) alone and in complex with the α2,3-sialylated trisaccharide revealed that in addition to the interaction between the MuV-HN active site residues and sialic acid, other residues, including an aromatic residue, stabilize the third sugar of the trisaccharide. The importance of the aromatic residue and the third sugar in the MuV-HN-receptor interaction was confirmed by computational energy calculations, isothermal titration calorimetry studies, and glycan-binding assays. Furthermore, MuV-HN was found to bind more efficiently to unbranched α2,3-sialylated sugar chains compared with branched ones. Importantly, the strategically located aromatic residue is conserved among the HN proteins of sialic acid-using paramyxoviruses, and alanine substitution compromised their ability to support cell-cell fusion. These results suggest that not only the terminal sialic acid but also the adjacent sugar moiety contribute to receptor function for mumps and these paramyxoviruses. The distribution of structurally different sialylated glycans in tissues and organs may explain in part MuV's distinct tropism to glandular tissues and the central nervous system. In the crystal structure, the epitopes for neutralizing antibodies are located around the α-helices of MuV-HN that are not well conserved in amino acid sequences among different genotypes of MuV. This may explain the fact that MuV reinfection sometimes occurs.

Comparative analysis of membrane protein structure databases.

BACKGROUND: Membrane proteins play important roles in cell survival and cell communication, as they function as transporters, receptors, anchors and enzymes. They are also potential targets for drugs that block receptors or inhibit enzymes related to diseases. Although the number of known structures of membrane proteins is still small relative to the size of the proteome as a whole, many new membrane protein structures have been determined recently. SCOPE OF THE ARTICLE: We compared and analyzed the widely used membrane protein databases, mpstruc, Orientations of Proteins in Membranes (OPM), and PDBTM, as well as the extended dataset of mpstruc based on sequence similarity, the PDB structures whose classification field indicates that they are "membrane proteins" and the proteins with Structural Classification of Proteins (SCOP) class-f domains. We evaluated the relationships between these databases or datasets based on the overlap in their contents and the degree of consistency in the structural, topological, and functional classifications and in the transmembrane domain assignment. MAJOR CONCLUSIONS: The membrane databases differ from each other in their coverage, and in the criteria that they use for annotation and classification. To ensure the efficient use of these databases, it is important to understand their differences and similarities. The establishment of more detailed and consistent annotations for the sequence, structure, membrane association, and function of membrane proteins is still required. GENERAL SIGNIFICANCE: Considering the recent growth of experimentally determined structures, a broad survey and cumulative analysis of the sum of knowledge as presented in the membrane protein structure databases can be helpful to elucidate structures and functions of membrane proteins. We also aim to provide a framework for future research and classification of membrane proteins.

Multiscale enhanced sampling of glucokinase: Regulation of the enzymatic reaction via a large scale domain motion.

Enhanced sampling yields a comprehensive structural ensemble or a free energy landscape, which is beyond the capability of a conventional molecular dynamics simulation. Our recently developed multiscale enhanced sampling (MSES) method employs a coarse-grained model coupled with the target physical system for the efficient acceleration of the dynamics. MSES has demonstrated applicability to large protein systems in solution, such as intrinsically disordered proteins and protein-protein and protein-ligand interactions. Here, we applied the MSES simulation to an important drug discovery target, glucokinase (GCK), to elucidate the structural basis of the positive cooperativity of the enzymatic reaction at an atomistic resolution. MSES enabled us to compare two sets of the free energy landscapes of GCK, for the glucose-bound and glucose-unbound forms, and thus demonstrated the drastic change of the free energy surface depending on the glucose concentration. In the glucose-bound form, we found two distinct basins separated by a high energy barrier originating from the domain motion and the folding/unfolding of the α13 helix. By contrast, in the glucose-unbound form, a single flat basin extended to the open and super-open states. These features illustrated the two distinct phases achieving the cooperativity, the fast reaction cycle staying in the closed state at a high glucose concentration and the slow cycle primarily in the open/super-open state at a low concentration. The weighted ensemble simulations revealed the kinetics of the structural changes in GCK with the synergetic use of the MSES results; the rate constant of the transition between the closed state and the open/super-open states, kC/O = 1.1 ms-1, is on the same order as the experimental catalytic rate, kcat = 0.22 ms-1. Finally, we discuss the pharmacological activities of GCK activators (small molecular drugs modulating the GCK activity) in terms of the slight changes in the domain motion, depending on their chemical structures as regulators. The present study demonstrated the capability of the enhanced sampling and the associated kinetic calculations for understanding the atomistic structural dynamics of protein systems in physiological environments.

Development of a protein-ligand-binding site prediction method based on interaction energy and sequence conservation.

We present a new method for predicting protein-ligand-binding sites based on protein three-dimensional structure and amino acid conservation. This method involves calculation of the van der Waals interaction energy between a protein and many probes placed on the protein surface and subsequent clustering of the probes with low interaction energies to identify the most energetically favorable locus. In addition, it uses amino acid conservation among homologous proteins. Ligand-binding sites were predicted by combining the interaction energy and the amino acid conservation score. The performance of our prediction method was evaluated using a non-redundant dataset of 348 ligand-bound and ligand-unbound protein structure pairs, constructed by filtering entries in a ligand-binding site structure database, LigASite. Ligand-bound structure prediction (bound prediction) indicated that 74.0 % of predicted ligand-binding sites overlapped with real ligand-binding sites by over 25 % of their volume. Ligand-unbound structure prediction (unbound prediction) indicated that 73.9 % of predicted ligand-binding residues overlapped with real ligand-binding residues. The amino acid conservation score improved the average prediction accuracy by 17.0 and 17.6 points for the bound and unbound predictions, respectively. These results demonstrate the effectiveness of the combined use of the interaction energy and amino acid conservation in the ligand-binding site prediction.

Development of a prediction system for tail-anchored proteins.

BACKGROUND: "Tail-anchored (TA) proteins" is a collective term for transmembrane proteins with a C-terminal transmembrane domain (TMD) and without an N-terminal signal sequence. TA proteins account for approximately 3-5 % of all transmembrane proteins that mediate membrane fusion, regulation of apoptosis, and vesicular transport. The combined use of TMD and signal sequence prediction tools is typically required to predict TA proteins. RESULTS: Here we developed a prediction system named TAPPM that predicted TA proteins solely from target amino acid sequences according to the knowledge of the sequence features of TMDs and the peripheral regions of TA proteins. Manually curated TA proteins were collected from published literature. We constructed hidden markov models of TA proteins as well as three different types of transmembrane proteins with similar structures and compared their likelihoods as TA proteins. CONCLUSIONS: Using the HMM models, we achieved high prediction accuracy; area under the receiver operator curve values reaching 0.963. A command line tool written in Python is available at https://github.com/davecao/tappm_cli .

Energy landscape of all-atom protein-protein interactions revealed by multiscale enhanced sampling.

Protein-protein interactions are regulated by a subtle balance of complicated atomic interactions and solvation at the interface. To understand such an elusive phenomenon, it is necessary to thoroughly survey the large configurational space from the stable complex structure to the dissociated states using the all-atom model in explicit solvent and to delineate the energy landscape of protein-protein interactions. In this study, we carried out a multiscale enhanced sampling (MSES) simulation of the formation of a barnase-barstar complex, which is a protein complex characterized by an extraordinary tight and fast binding, to determine the energy landscape of atomistic protein-protein interactions. The MSES adopts a multicopy and multiscale scheme to enable for the enhanced sampling of the all-atom model of large proteins including explicit solvent. During the 100-ns MSES simulation of the barnase-barstar system, we observed the association-dissociation processes of the atomistic protein complex in solution several times, which contained not only the native complex structure but also fully non-native configurations. The sampled distributions suggest that a large variety of non-native states went downhill to the stable complex structure, like a fast folding on a funnel-like potential. This funnel landscape is attributed to dominant configurations in the early stage of the association process characterized by near-native orientations, which will accelerate the native inter-molecular interactions. These configurations are guided mostly by the shape complementarity between barnase and barstar, and lead to the fast formation of the final complex structure along the downhill energy landscape.