Search

The performance of molecular docking can be improved by comparing the shape similarity of the flexibly sampled poses against the target proteins' inverted binding cavities. The effectiveness of these pseudo-ligands or negative image-based models in docking rescoring is boosted further by performing enrichment-driven optimization. Here, we introduce a novel shape-focused pharmacophore modeling algorithm O-LAP that generates a new class of cavity-filling models by clumping together overlapping atomic content via pairwise distance graph clustering. Top-ranked poses of flexibly docked active ligands were used as the modeling input and multiple alternative clustering settings were benchmark-tested thoroughly with five demanding drug targets using random training/test divisions. In docking rescoring, the O-LAP modeling typically improved massively on the default docking enrichment; furthermore, the results indicate that the clustered models work well in rigid docking. The C+ +/Qt5-based algorithm O-LAP is released under the GNU General Public License v3.0 via GitHub ( https://github.com/jvlehtonen/overlap-toolkit ). SCIENTIFIC CONTRIBUTION: This study introduces O-LAP, a C++/Qt5-based graph clustering software for generating new type of shape-focused pharmacophore models. In the O-LAP modeling, the target protein cavity is filled with flexibly docked active ligands, the overlapping ligand atoms are clustered, and the shape/electrostatic potential of the resulting model is compared against the flexibly sampled molecular docking poses. The O-LAP modeling is shown to ensure high enrichment in both docking rescoring and rigid docking based on comprehensive benchmark-testing., (© 2024. The Author(s).)

Molecular docking is a key method used in virtual screening (VS) campaigns to identify small-molecule ligands for drug discovery targets. While docking provides a tangible way to understand and predict the protein-ligand complex formation, the docking algorithms are often unable to separate active ligands from inactive molecules in practical VS usage. Here, a novel docking and shape-focused pharmacophore VS protocol is demonstrated for facilitating effective hit discovery using retinoic acid receptor-related orphan receptor gamma t (RORγt) as a case study. RORγt is a prospective target for treating inflammatory diseases such as psoriasis and multiple sclerosis. First, a commercial molecular database was flexibly docked. Second, the alternative docking poses were rescored against the shape/electrostatic potential of negative image-based (NIB) models that mirror the target's binding cavity. The compositions of the NIB models were optimized via iterative trimming and benchmarking using a greedy search-driven algorithm or brute force NIB optimization. Third, a pharmacophore point-based filtering was performed to focus the hit identification on the known RORγt activity hotspots. Fourth, free energy binding affinity evaluation was performed on the remaining molecules. Finally, twenty-eight compounds were selected for in vitro testing and eight compounds were determined to be low μM range RORγt inhibitors, thereby showing that the introduced VS protocol generated an effective hit rate of ~29%.

Despite the pivotal role of molecular docking in modern drug discovery, the default docking scoring functions often fail to recognize active ligands in virtual screening campaigns. Negative image-based rescoring improves docking enrichment by comparing the shape/electrostatic potential (ESP) of the flexible docking poses against the target protein's inverted cavity volume. By optimizing these negative image-based (NIB) models using a greedy search, the docking rescoring yield can be improved massively and consistently. Here, a fundamental modification is implemented to this shape-focused pharmacophore modelling approach-actual ligand 3D coordinates are incorporated into the NIB models for the optimization. This hybrid approach, labelled as ligand-enhanced brute-force negative image-based optimization (LBR-NiB), takes the best from both worlds, i.e., the all-roundedness of the NIB models and the difficult to emulate atomic arrangements of actual protein-bound small-molecule ligands. Thorough benchmarking, focused on proinflammatory targets, shows that the LBR-NiB routinely improves the docking enrichment over prior iterations of the R-NiB methodology. This boost can be massive, if the added ligand information provides truly essential binding information that was lacking or completely missing from the cavity-based NIB model. On a practical level, the results indicate that the LBR-NiB typically works well when the added ligand 3D data originates from a high-quality source, such as X-ray crystallography, and, yet, the NIB model compositions can also sometimes be improved by fusing into them, for example, with flexibly docked solvent molecules. In short, the study demonstrates that the protein-bound ligands can be used to improve the shape/ESP features of the negative images for effective docking rescoring use in virtual screening.

Molecular docking is a key in silico method used routinely in modern drug discovery projects. Although docking provides high-quality ligand binding predictions, it regularly fails to separate the active compounds from the inactive ones. In negative image-based rescoring (R-NiB), the shape/electrostatic potential (ESP) of docking poses is compared to the negative image of the protein's ligand binding cavity. While R-NiB often improves the docking yield considerably, the cavity-based models do not reach their full potential without expert editing. Accordingly, a greedy search-driven methodology, brute force negative image-based optimization (BR-NiB), is presented for optimizing the models via iterative editing and benchmarking. Thorough and unbiased training, testing and stringent validation with a multitude of drug targets, and alternative docking software show that BR-NiB ensures excellent docking efficacy. BR-NiB can be considered as a new type of shape-focused pharmacophore modeling, where the optimized models contain only the most vital cavity information needed for effectively filtering docked actives from the inactive or decoy compounds. Finally, the BR-NiB code for performing the automated optimization is provided free-of-charge under MIT license via GitHub (https://github.com/jvlehtonen/brutenib) for boosting the success rates of docking-based virtual screening campaigns.

In this study, we use an optogenetic inhibitor of c-Jun NH 2 -terminal kinase (JNK) in dendritic spine sub-compartments of rat hippocampal neurons. We show that JNK inhibition exerts rapid (within seconds) reorganization of actin in the spine-head. Using real-time Förster resonance energy transfer (FRET) to measure JNK activity, we find that either excitotoxic insult (NMDA) or endocrine stress (corticosterone), activate spine-head JNK causing internalization of AMPARs and spine retraction. Both events are prevented upon optogenetic inhibition of JNK, and rescued by JNK inhibition even 2 h after insult. Moreover, we identify that the fast-acting anti-depressant ketamine reduces JNK activity in hippocampal neurons suggesting that JNK inhibition may be a downstream mediator of its anti-depressant effect. In conclusion, we show that JNK activation plays a role in triggering spine elimination by NMDA or corticosterone stress, whereas inhibition of JNK facilitates regrowth of spines even in the continued presence of glucocorticoid. This identifies that JNK acts locally in the spine-head to promote AMPAR internalization and spine shrinkage following stress, and reveals a protective function for JNK inhibition in preventing spine regression., (Copyright © 2020 Hollos et al.)

We introduce a statistical method for evaluating atomic level 3D interaction patterns of protein-ligand contacts. Such patterns can be used for fast separation of likely ligand and ligand binding site combinations out of all those that are geometrically possible. The practical purpose of this probabilistic method is for molecular docking and scoring, as an essential part of a scoring function. Probabilities of interaction patterns are calculated conditional on structural x-ray data and predefined chemical classification of molecular fragment types. Spatial coordinates of atoms are modeled using a Bayesian statistical framework with parametric 3D probability densities. The parameters are given distributions a priori, which provides the possibility to update the densities of model parameters with new structural data and use the parameter estimates to create a contact hierarchy. The contact preferences can be defined for any spatial area around a specified type of fragment. We compared calculated contact point hierarchies with the number of contact atoms found near the contact point in a reference set of x-ray data, and found that these were in general in a close agreement. Additionally, using substrate binding site in cathechol-O-methyltransferase and 27 small potential binder molecules, it was demonstrated that these probabilities together with auxiliary parameters separate well ligands from decoys (true positive rate 0.75, false positive rate 0). A particularly useful feature of the proposed Bayesian framework is that it also characterizes predictive uncertainty in terms of probabilities, which have an intuitive interpretation from the applied perspective.

Reliable and effective virtual high-throughput screening (vHTS) methods are desperately needed to minimize the expenses involved in drug discovery projects. Here, we present an improvement to the negative image-based (NIB) screening: the shape, the electrostatics, and the solvation state of the target protein's ligand-binding site are included into the vHTS. Additionally, the initial vHTS results are postprocessed with molecular mechanics/generalized Born surface area (MMGBSA) calculations to estimate the favorability of ligand-protein interactions. The results show that docking produces very good early enrichment for phosphodiesterase-5 (PDE-5); however, in general, the NIB and the ligand-based screening performed better with or without the added electrostatics. Furthermore, the postprocessing of the NIB screening results using MMGBSA calculations improved the early enrichment for the PDE-5 considerably, thus, making hit discovery affordable.

Endosialidase (endo-N-acetylneuraminidase) is a tailspike enzyme of bacteriophages specific for human pathogenic Escherichia coli K1, which specifically recognizes and degrades polySia (polysialic acid). polySia is also a polysaccharide of the capsules of other meningitis- and sepsis-causing bacteria, and a post-translational modification of the NCAM (neural cell-adhesion molecule). We have cloned and sequenced three spontaneously mutated endosialidases of the PK1A bacteriophage and one of the PK1E bacteriophage which display lost or residual enzyme activity but retain the binding activity to polySia. Single to triple amino acid substitutions were identified, and back-mutation constructs indicated that single substitutions accounted for only partial reduction of enzymic activity. A homology-based structural model of endosialidase revealed that all substituted amino acid residues localize to the active site of the enzyme. The results reveal the importance of non-catalytic amino acid residues for the enzymatic activity. The results reveal the molecular background for the dissociation of the polySia binding and cleaving activities of endosialidase and for the evolvement of 'host range' mutants of E. coli K1 bacteriophages.

BODIL is a molecular modeling environment geared to help the user to quickly identify key features of proteins critical to molecular recognition, especially (1) in drug discovery applications, and (2) to understand the structural basis for function. The program incorporates state-of-the-art graphics, sequence and structural alignment methods, among other capabilities needed in modern structure-function-drug target research. BODIL has a flexible design that allows on-the-fly incorporation of new modules, has intelligent memory management, and fast multi-view graphics. A beta version of BODIL and an accompanying tutorial are available at http://www.abo.fi/fak/mnf/bkf/research/johnson/bodil.html.

A detailed comparison of the structures of aspartate aminotransferase, alanine race-mase, the beta subunit of tryptophan synthase, D-amino acid aminotransferase and glycogen phosphorylase has revealed more extensive structural similarities among pyridoxal phosphate (PLP)-binding domains in these enzymes than was observed previously. These similarities consist of seven common structural segments of the polypeptide chain, which form an extensive common structural organization of the backbone chain responsible for the appropriate disposition of key residues, some from the aligned fragments and some from variable loops joined to these fragments, interacting with PLPs in these enzymes. This common structural organization contains an analogous hydrophobic minicore formed from four amino acid side chains present in the two most conserved structural elements. In addition, equivalent alpha-beta-alpha-beta supersecondary structures are formed by these seven fragments in three of the five structures: alanine racemase, tryptophan synthase and glycogen phosphorylase. Despite these similarities, it is generally accepted that these proteins do not share a common heritage, but have arisen on five separate occasions. The common and contiguous alpha-beta-alpha-beta structure accounts for only 28 residues and all five enzymes differ greatly in both the orientation of the PLP pyridoxal rings and their contacts with residues close to the common structural elements.

We have developed a generic tool for the automatic identification of regions of local structural similarity in unrelated proteins having different folds, as well as for defining more global similarities that result from homologous protein structures. The computer program GENFIT has evolved from the genetic algorithm-based three-dimensional protein structure comparison program GA_FIT. GENFIT, however, can locate and superimpose regions of local structural homology regardless of their position in a pair of structures, the fold topology, or the chain direction. Furthermore, it is possible to restrict the search to a volume centered about a region of interest (e.g., catalytic site, ligand-binding site) in two protein structures. We present a number of examples to illustrate the function of the program, which is a parallel processing implementation designed for distribution to multiple machines over a local network or to run on a single multiprocessor computer.

Three ATP-dependent enzymes with different folds, cAMP-dependent protein kinase, D-Ala:D-Ala ligase and the alpha-subunit of the alpha2beta2 ribonucleotide reductase, have a similar organization of their ATP-binding sites. The most meaningful similarity was found over 23 structurally equivalent residues in each protein and includes three strands each from their beta-sheets, in addition to a connecting loop. The equivalent secondary structure elements in each of these enzymes donate four amino acids forming key hydrogen bonds responsible for the common orientation of the "AMP" moieties of their ATP-ligands. One lysine residue conserved throughout the three families binds the alpha-phosphate in each protein. The common fragments of structure also position some, but not all, of the equivalent residues involved in hydrophobic contacts with the adenine ring. These examples of convergent evolution reinforce the view that different proteins can fold in different ways to produce similar structures locally, and nature can take advantage of these features when structure and function demand it, as shown here for the common mode of ATP-binding by three unrelated proteins.

Two proteins, D-alanine:D-alanine ligase and cAMP-dependent protein kinase, share a remarkable degree of structural convergence despite having different three-dimensional folds and different enzymatic functions. Here we report that as many as 103 residues from 10 segments form two identical super-secondary structures between which the cofactor ATP is bound. The cofactor, two bound metal cations, and several water molecules form a large network of electrostatic and hydrophobic interactions common to both enzymes, and these are mediated by the similar placement of equivalent amino acids within the common supersecondary structures.