The High-Affinity Chymotrypsin Inhibitor Eglin C Poorly Inhibits Human Chymotrypsin-Like Protease: Gln192 and Lys218 Are Key Determinants.

Eglin C, a small protein from the medicinal leech, has been long considered a general high-affinity inhibitor of chymotrypsins and elastases. Here, we demonstrate that eglin C inhibits human chymotrypsin-like protease (CTRL) weaker by several orders of magnitude than other chymotrypsins. In order to identify the underlying structural aspects of this unique deviation, we performed comparative molecular dynamics simulations on experimental and AlphaFold model structures of bovine CTRA and human CTRL. Our results indicate that in CTRL, the primary determinants of the observed weak inhibition are amino-acid positions 192 and 218 (using conventional chymotrypsin numbering), which participate in shaping the S1 substrate-binding pocket and thereby affect the stability of the protease-inhibitor complexes.

Correction: Cserzo et al. The First Quarter Century of the Dense Alignment Surface Transmembrane Prediction Method. Int. J. Mol. Sci. 2023, 24, 14016.

There was an error in the original publication [...].

The First Quarter Century of the Dense Alignment Surface Transmembrane Prediction Method.

The dense alignment surface (DAS) transmembrane (TM) prediction method was first published more than 25 years ago. DAS was the one of the earliest tools to discriminate TM proteins from globular ones and to predict the sequence positions of TM helices in proteins with high accuracy from their amino acid sequence alone. The algorithmic improvements that followed in 2002 (DAS-TMfilter) made it one of the best performing tools among those relying on local sequence information for TM prediction. Since then, many more experimental data about membrane proteins (including thousands of 3D structures of membrane proteins) have accumulated but there has been no significant improvement concerning performance in the area of TM helix prediction tools. Here, we report a new implementation of the DAS-TMfilter prediction web server. We reevaluated the performance of the method using a five-times-larger, updated test dataset. We found that the method performs at essentially the same accuracy as the original even without any change to the parametrization of the program despite the much larger dataset. Thus, the approach captures the physico-chemistry of TM helices well, essentially solving this scientific problem.

Molecular Dynamics Simulation as a Tool to Identify Mutual Synergistic Folding Proteins.

Mutual synergistic folding (MSF) proteins belong to a recently emerged subclass of disordered proteins, which are disordered in their monomeric forms but become ordered in their oligomeric forms. They can be identified by experimental methods following their unfolding, which happens in a single-step cooperative process, without the presence of stable monomeric intermediates. Only a limited number of experimentally validated MSF proteins are accessible. The amino acid composition of MSF proteins shows high similarity to globular ordered proteins, rather than to disordered ones. However, they have some special structural features, which makes it possible to distinguish them from globular proteins. Even in the possession of their oligomeric three-dimensional structure, classification can only be performed based on unfolding experiments, which are frequently absent. In this work, we demonstrate a simple protocol using molecular dynamics simulations, which is able to indicate that a protein structure belongs to the MSF subclass. The presumption of the known atomic resolution quaternary structure is an obvious limitation of the method, and because of its high computational time requirements, it is not suitable for screening large databases; still, it is a valuable in silico tool for identification of MSF proteins.

Assortment of Frontiers in Protein Science.

Recent decades have brought significant changes to the protein structure research field [...].

Origin of Increased Solvent Accessibility of Peptide Bonds in Mutual Synergetic Folding Proteins.

Mutual Synergetic Folding (MSF) proteins belong to a recently discovered class of proteins. These proteins are disordered in their monomeric but ordered in their oligomeric forms. Their amino acid composition is more similar to globular proteins than to disordered ones. Our preceding work shed light on important structural aspects of the structural organization of these proteins, but the background of this behavior is still unknown. We suggest that solvent accessibility is an important factor, especially solvent accessibility of the peptide bonds can be accounted for this phenomenon. The side chains of the amino acids which form a peptide bond have a high local contribution to the shielding of the peptide bond from the solvent. During the oligomerization step, other non-local residues contribute to the shielding. We investigated these local and non-local effects of shielding based on Shannon information entropy calculations. We found that MSF and globular homodimeric proteins have different local contributions resulting from different amino acid pair frequencies. Their non-local distribution is also different because of distinctive inter-subunit contacts.

Analysis of Heterodimeric "Mutual Synergistic Folding"-Complexes.

Several intrinsically disordered proteins (IDPs) are capable to adopt stable structures without interacting with a folded partner. When the folding of all interacting partners happens at the same time, coupled with the interaction in a synergistic manner, the process is called Mutual Synergistic Folding (MSF). These complexes represent a discrete subset of IDPs. Recently, we collected information on their complexes and created the MFIB (Mutual Folding Induced by Binding) database. In a previous study, we compared homodimeric MSF complexes with homodimeric and monomeric globular proteins with similar amino acid sequence lengths. We concluded that MSF homodimers, compared to globular homodimeric proteins, have a greater solvent accessible main-chain surface area on the contact surface of the subunits, which becomes buried during dimerization. The main driving force of the folding is the mutual shielding of the water-accessible backbones, but the formation of further intermolecular interactions can also be relevant. In this paper, we will report analyses of heterodimeric MSF complexes. Our results indicate that the amino acid composition of the heterodimeric MSF monomer subunits slightly diverges from globular monomer proteins, while after dimerization, the amino acid composition of the overall MSF complexes becomes more similar to overall amino acid compositions of globular complexes. We found that inter-subunit interactions are strengthened, and additionally to the shielding of the solvent accessible backbone, other factors might play an important role in the stabilization of the heterodimeric structures, likewise energy gain resulting from the interaction of the two subunits with different amino acid compositions. We suggest that the shielding of the ß-sheet backbones and the formation of a buried structural core along with the general strengthening of inter-subunit interactions together could be the driving forces of MSF protein structural ordering upon dimerization.

Physical Background of the Disordered Nature of "Mutual Synergetic Folding" Proteins.

Intrinsically disordered proteins (IDPs) lack a well-defined 3D structure. Their disordered nature enables them to interact with several other proteins and to fulfil their vital biological roles, in most cases after coupled folding and binding. In this paper, we analyze IDPs involved in a new mechanism, mutual synergistic folding (MSF). These proteins define a new subset of IDPs. Recently we collected information on these complexes and created the Mutual Folding Induced by Binding (MFIB) database. These protein complexes exhibit considerable structural variation, and almost half of them are homodimers, but there is a significant amount of heterodimers and various kinds of oligomers. In order to understand the basic background of the disordered character of the monomers found in MSF complexes, the simplest part of the MFIB database, the homodimers are analyzed here. We conclude that MFIB homodimeric proteins have a larger solvent-accessible main-chain surface area on the contact surface of the subunits, when compared to globular homodimeric proteins. The main driving force of the dimerization is the mutual shielding of the water-accessible backbones and the formation of extra intermolecular interactions.

Identification of potential glutaminyl cyclase inhibitors from lead-like libraries by in silico and in vitro fragment-based screening.

A glutaminyl cyclase (QC) fragment library was in silico selected by disconnection of the structure of known QC inhibitors and by lead-like 2D virtual screening of the same set. The resulting fragment library (204 compounds) was acquired from commercial suppliers and pre-screened by differential scanning fluorimetry followed by functional in vitro assays. In this way, 10 fragment hits were identified ([Formula: see text]5 % hit rate, best inhibitory activity: 16 [Formula: see text]). The in vitro hits were then docked to the active site of QC, and the best scoring compounds were analyzed for binding interactions. Two fragments bound to different regions in a complementary manner, and thus, linking those fragments offered a rational strategy to generate novel QC inhibitors. Based on the structure of the virtual linked fragment, a 77-membered QC target focused library was selected from vendor databases and docked to the active site of QC. A PubChem search confirmed that the best scoring analogues are novel, potential QC inhibitors.

The role of stabilization centers in protein thermal stability.

The definition of stabilization centers was introduced almost two decades ago. They are centers of noncovalent long range interaction clusters, believed to have a role in maintaining the three-dimensional structure of proteins by preventing their decay due to their cooperative long range interactions. Here, this hypothesis is investigated from the viewpoint of thermal stability for the first time, using a large protein thermodynamics database. The positions of amino acids belonging to stabilization centers are correlated with available experimental thermodynamic data on protein thermal stability. Our analysis suggests that stabilization centers, especially solvent exposed ones, do contribute to the thermal stabilization of proteins.

Combination of 2D/3D ligand-based similarity search in rapid virtual screening from multimillion compound repositories. Selection and biological evaluation of potential PDE4 and PDE5 inhibitors.

Rapid in silico selection of target focused libraries from commercial repositories is an attractive and cost effective approach. If structures of active compounds are available rapid 2D similarity search can be performed on multimillion compound databases but the generated library requires further focusing by various 2D/3D chemoinformatics tools. We report here a combination of the 2D approach with a ligand-based 3D method (Screen3D) which applies flexible matching to align reference and target compounds in a dynamic manner and thus to assess their structural and conformational similarity. In the first case study we compared the 2D and 3D similarity scores on an existing dataset derived from the biological evaluation of a PDE5 focused library. Based on the obtained similarity metrices a fusion score was proposed. The fusion score was applied to refine the 2D similarity search in a second case study where we aimed at selecting and evaluating a PDE4B focused library. The application of this fused 2D/3D similarity measure led to an increase of the hit rate from 8.5% (1st round, 47% inhibition at 10 µM) to 28.5% (2nd round at 50% inhibition at 10 µM) and the best two hits had 53 nM inhibitory activities.

Synergistic transport of a fluorescent coumarin probe marks coumarins as pharmacological modulators of Organic anion-transporting polypeptide, OATP3A1.

Organic anion-transporting polypeptide 3A1 (OATP3A1) is a membrane transporter mediating the cellular uptake of various hormones such as estrone-3-sulfate, prostaglandins E1 and E2 and thyroxine. OATP3A1 is widely expressed in the human body and its presence in tissue-blood barriers, neurons and muscle cells marks it as a potential pharmacological target. Herein we demonstrate that an otherwise membrane impermeant, zwitterionic fluorescent coumarin probe, bearing a sulfonate function is a potent substrate of human OATP3A1, thus readily transported into HEK-293-OATP3A1 cells allowing functional investigation and the screen of drug interactions of the OATP3A1 transporter. At the same time, dyes lacking either the sulfonate motif or the coumarin scaffold showed a dramatic decrease in affinity or even a complete loss of transport. Furthermore, we observed a distinct inhibition/activation pattern in the OATP3A1-mediated uptake of closely related fluorescent coumarin derivatives differing only in the presence of the sulfonate moiety. Additionally, we detected a synergistic effect between one of the probes tested and the endogenous OATP substrate estrone-3-sulfate. These data, together with docking results indicate the presence of at least two cooperative substrate binding sites in OATP3A1. Besides providing the first sensitive probe for testing OATP3A1 substrate/inhibitor interactions, our results also help to understand substrate recognition and transport mechanism of the poorly characterized OATP3A1. Moreover, coumarins are good candidates for OATP3A1-targeted drug delivery and as pharmacological modulators of OATP3A1.

Identifying new topoisomerase II poison scaffolds by combining publicly available toxicity data and 2D/3D-based virtual screening.

Molecular descriptor (2D) and three dimensional (3D) shape based similarity methods are widely used in ligand based virtual drug design. In the present study pairwise structure comparisons among a set of 4858 DTP compounds tested in the NCI60 tumor cell line anticancer drug screen were computed using chemical hashed fingerprints and 3D molecule shapes to calculate 2D and 3D similarities, respectively. Additionally, pairwise biological activity similarities were calculated by correlating the 60 element vectors of pGI50 values corresponding to the cytotoxicity of the compounds across the NCI60 panel. Subsequently, we compared the power of 2D and 3D structural similarity metrics to predict the toxicity pattern of compounds. We found that while the positive predictive value and sensitivity of 3D and molecular descriptor based approaches to predict biological activity are similar, a subset of molecule pairs yielded contradictory results. By simultaneously requiring similarity of biological activities and 3D shapes, and dissimilarity of molecular descriptor based comparisons, we identify pairs of scaffold hopping candidates displaying characteristic core structural changes such as heteroatom/heterocycle change and ring closure. Attempts to discover scaffold hopping candidates of mitoxantrone recovered known Topoisomerase II (Top2) inhibitors, and also predicted new, previously unknown chemotypes possessing in vitro Top2 inhibitory activity.

SRide: a server for identifying stabilizing residues in proteins.

Residues expected to play key roles in the stabilization of proteins [stabilizing residues (SRs)] are selected by combining several methods based mainly on the interactions of a given residue with its spatial, rather than its sequential neighborhood and by considering the evolutionary conservation of the residues. A residue is selected as a stabilizing residue if it has high surrounding hydrophobicity, high long-range order, high conservation score and if it belongs to a stabilization center. The definition of all these parameters and the thresholds used to identify the SRs are discussed in detail. The algorithm for identifying SRs was originally developed for TIM-barrel proteins [M. M. Gromiha, G. Pujadas, C. Magyar, S. Selvaraj, and I. Simon (2004), Proteins, 55, 316-329] and is now generalized for all proteins of known 3D structure. SRs could be applied in protein engineering and homology modeling and could also help to explain certain folds with significant stability. The SRide server is located at http://sride.enzim.hu.

Phosphorylation-induced transient intrinsic structure in the kinase-inducible domain of CREB facilitates its recognition by the KIX domain of CBP.

Phosphorylation at Ser-133 of the kinase inducible domain of CREB (KID) triggers its binding to the KIX domain of CBP via a concomitant coil-to-helix transition. The exact role of this key event is still puzzling: it does not switch between disordered and ordered states, nor its direct interactions fully account for selectivity. Hence, we reasoned that phosphorylation may shift the conformational preferences of KID towards a binding-competent state. To this end we investigated the intrinsic structural properties of the unbound KID in phosphorylated and unphosphorylated forms by simulated annealing and molecular dynamics simulations. Although helical populations show subtle differences, phosphorylation reduces the flexibility of the turn segment connecting the two helices in the complexed structure and induces a transient structural element that corresponds to its bound conformation. It is stabilized by the pSer-133-Arg-131 interaction, which is absent from the unphosphorylated KID. Diminishing this coupling decreases the 3.1 kcal/mol contribution of pSer-133 to the binding free energy (DeltaGbind) of the phosphorylated KID to KIX by 1.1 kcal/mol, as computed in reference to Ser-133. In a binding competent form of the S133E KID mutant, the contribution of Glu-133 to DeltaGbind is by 1.5 kcal/mol smaller than that of pSer, suggesting that altered structural properties due to pSer --> Glu replacement impair the binding affinity. Thus, we propose that phoshorylation contributes to selectivity not merely by the direct interactions of the phosphate group with KIX, but also by promoting the formation of a transient structural element in the highly conserved turn segment.

Servers for sequence-structure relationship analysis and prediction.

We describe several algorithms and public servers that were developed to analyze and predict various features of protein structures. These servers provide information about the covalent state of cysteine (CYSREDOX), as well as about residues involved in non-covalent cross links that play an important role in the structural stability of proteins (SCIDE and SCPRED). We also discuss methods and servers developed to identify helical transmembrane proteins from large databases and rough genomic data, including two of the most popular transmembrane prediction methods, DAS and HMMTOP. Several biologically interesting applications of these servers are also presented. The servers are available through http://www.enzim.hu/servers.html.

Access Path to the Ligand Binding Pocket May Play a Role in Xenobiotics Selection by AhR.

Understanding of multidrug binding at the atomic level would facilitate drug design and strategies to modulate drug metabolism, including drug transport, oxidation, and conjugation. Therefore we explored the mechanism of promiscuous binding of small molecules by studying the ligand binding domain, the PAS-B domain of the aryl hydrocarbon receptor (AhR). Because of the low sequence identities of PAS domains to be used for homology modeling, structural features of the widely employed HIF-2α and a more recent suitable template, CLOCK were compared. These structures were used to build AhR PAS-B homology models. We performed molecular dynamics simulations to characterize dynamic properties of the PAS-B domain and the generated conformational ensembles were employed in in silico docking. In order to understand structural and ligand binding features we compared the stability and dynamics of the promiscuous AhR PAS-B to other PAS domains exhibiting specific interactions or no ligand binding function. Our exhaustive in silico binding studies, in which we dock a wide spectrum of ligand molecules to the conformational ensembles, suggest that ligand specificity and selection may be determined not only by the PAS-B domain itself, but also by other parts of AhR and its protein interacting partners. We propose that ligand binding pocket and access channels leading to the pocket play equally important roles in discrimination of endogenous molecules and xenobiotics.

Flexibility of prolyl oligopeptidase: molecular dynamics and molecular framework analysis of the potential substrate pathways.

The flexibility of prolyl oligopeptidase has been investigated using molecular dynamics (MD) and molecular framework approaches to delineate the route of the substrate to the active site. The selectivity of the enzyme is mediated by a seven-bladed beta-propeller that in the crystal structure does not indicate the possible passage for the substrate to the catalytic center. Its open topology however, could allow the blades to move apart and let the substrate into the large central cavity. Flexibility analysis of prolyl oligopeptidase structure using the FIRST (Floppy Inclusion and Rigid Substructure Topology) approach and the atomic fluctuations derived from MD simulations demonstrated the rigidity of the propeller domain, which does not permit the substrate to approach the active site through this domain. Instead, a smaller tunnel at the inter-domain region comprising the highly flexible N-terminal segment of the peptidase domain and a facing hydrophilic loop from the propeller (residues 192-205) was identified by cross-correlation analysis and essential dynamics as the only potential pathway for the substrate. The functional importance of the flexible loop has been also verified by kinetic analysis of the enzyme with a split loop. Catalytic effect of engineered disulfide bridges was rationalized by characterizing the concerted motions of the two domains.

Combination of Pharmacophore Matching, 2D Similarity Search, and In Vitro Biological Assays in the Selection of Potential 5-HT6 Antagonists from Large Commercial Repositories.

Rapid in silico selection of target-focused libraries from commercial repositories is an attractive and cost-effective approach. If structures of active compounds are available, rapid 2D similarity search can be performed on multimillion compound databases, but the generated library requires further focusing. We report here a combination of the 2D approach with pharmacophore matching which was used for selecting 5-HT6 antagonists. In the first screening round, 12 compounds showed >85% antagonist efficacy of the 91 screened. For the second-round (hit validation) screening phase, pharmacophore models were built, applied, and compared with the routine 2D similarity search. Three pharmacophore models were created based on the structure of the reference compounds and the first-round hit compounds. The pharmacophore search resulted in a high hit rate (40%) and led to novel chemotypes, while 2D similarity search had slightly better hit rate (51%), but lacking the novelty. To demonstrate the power of the virtual screening cascade, ligand efficiency indices were also calculated and their steady improvement was confirmed.

Locating the stabilizing residues in (alpha/beta)8 barrel proteins based on hydrophobicity, long-range interactions, and sequence conservation.

In nature, 1 out of every 10 proteins has an (alpha/beta)(8) (TIM)-barrel fold, and in most cases, pairwise comparisons show no sequence similarity between them. Hence, delineating the key residues that induce very different sequences to share a common fold is important for understanding the folding and stability of TIM-barrel domains. In this work, we propose a new consensus approach for locating these stabilizing residues based on long-range interactions, hydrophobicity, and conservation of amino acid residues. We have identified 957 stabilizing residues in 63 proteins from a nonredundant set of 71 TIM-barrel domains. Most of these residues are located in the 8-stranded beta-sheet, with nearly one half of them oriented toward the interior of the barrel and the other half oriented toward the surrounding alpha-helices. Several stabilizing residues are found in the N- and C-terminal loops, whereas very few appear in the alpha-helices that surround the internal beta-sheet. Further, these 957 residues are placed in 434 stabilizing segments of various sizes, and each domain contains 1-10 of these segments. We found that 8 segments per domain is the most abundant one, and two thirds of the proteins have 7-9 stabilizing segments. Finally, we verified the identified residues with experimental temperature factors and found that these residues are among the ones with less mobility in the considered proteins. We suggest that our new protocol serves as a powerful tool to identify the stabilizing residues in TIM-barrel domains, which can be used as potential candidates for studying protein folding and stability by means of protein engineering experiments.