Simplified quality assessment for small-molecule ligands in the Protein Data Bank

Halogen Bonding on Water─A Drop in the Ocean?

Halogen bonding is a valuable interaction in drug design, offering an unconventional way to influence affinity and selectivity by leveraging the halogen atoms' ability to form directional bonds. The present study evaluates halogen-water interactions within protein binding sites, demonstrating that targeting a water molecule via halogen bonding can in specific cases contribute beneficially to ligand binding. In solving and examining the crystal structure of 2-cyclopentyl-7-iodo-1H-indole-3-carbonitrile bound to DYRK1a kinase, we identified a notable iodine-water interaction, where water accepts a halogen bond with good geometric and energetic features. This starting point triggered further investigations into the prevalence of such interactions across various halogen-bearing ligands (chlorine, bromine, iodine) in the PDB. Using QM calculations (MP2/TZVPP), we highlight the versatility and potential benefits of such halogen-water interactions, particularly when the water molecule is a stable part of the binding site's structured environment. While the interaction energies with water are lower compared to other typical halogen bond acceptors, we deem this different binding strength essential for reducing desolvation costs. We suggest that "interstitial" water molecules, as stable parts of the binding site engaging in multiple strong interactions, could be prime targets for halogen bonding. Further systematic studies, combining high-resolution crystal structures and quantum chemistry, are required to scrutinize whether halogen bonding on water is more than a "drop in the ocean".

Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Computational exploration of SLC14A1 genetic variants through structure modeling, protein-ligand docking, and molecular dynamics simulation

The urea transporter UT-B1, encoded by the SLC14A1 gene, has been hypothesized to be a significant protein whose deficiency and dysfunction contribute to the pathogenesis of bladder cancer and many other diseases. Several studies reported the association of genetic alterations in the SLC14A1 (UT-B1) gene with bladder carcinogenesis, suggesting a need for thorough characterization of the UT-B1 protein's coding and non-coding variants. This study used various computational techniques to investigate the commonly occurring germ-line missense and non-coding SNPs (ncSNPs) of the SLC14A1 gene (UT-B1) for their structural, functional, and molecular implications for disease susceptibility and dysfunctionality. SLC14A1 missense variants, primarily identified from the ENSEMBL genome browser, were screened through twelve functionality prediction tools leading to two variants D280Y (predicted detrimental by maximum tools) and D280N (high global MAF) for rs1058396. Subsequently, the ConSurf and NetSurf tools revealed the D280 residue to be in a variable site and exposed on the protein surface. According to I-Mutant2.0 and MUpro, both variants are predicted to cause a significant effect on protein stability. Analysis of molecular docking anticipated these two variants to decrease the binding affinity of UT-B1 protein for the examined ligands to a significant extent. Molecular dynamics also disclosed the possible destabilization of the UT-B1 protein due to single nucleotide polymorphism compared to wild-type protein which may result in impaired protein function. Furthermore, several non-coding SNPs were estimated to affect transcription factor binding and regulation of SLC14A1 gene expression. Additionally, two ncSNPs were found to affect miRNA-based post-transcriptional regulation by creating new seed regions for miRNA binding. This comprehensive in-silico study of SLC14A1 gene variants may serve as a springboard for future large-scale investigations examining SLC14A1 polymorphisms.

PoseBusters: AI-based docking methods fail to generate physically valid poses or generalise to novel sequences

The last few years have seen the development of numerous deep learning-based protein-ligand docking methods. They offer huge promise in terms of speed and accuracy. However, despite claims of state-of-the-art performance in terms of crystallographic root-mean-square deviation (RMSD), upon closer inspection, it has become apparent that they often produce physically implausible molecular structures. It is therefore not sufficient to evaluate these methods solely by RMSD to a native binding mode. It is vital, particularly for deep learning-based methods, that they are also evaluated on steric and energetic criteria. We present PoseBusters, a Python package that performs a series of standard quality checks using the well-established cheminformatics toolkit RDKit. The PoseBusters test suite validates chemical and geometric consistency of a ligand including its stereochemistry, and the physical plausibility of intra- and intermolecular measurements such as the planarity of aromatic rings, standard bond lengths, and protein-ligand clashes. Only methods that both pass these checks and predict native-like binding modes should be classed as having "state-of-the-art" performance. We use PoseBusters to compare five deep learning-based docking methods (DeepDock, DiffDock, EquiBind, TankBind, and Uni-Mol) and two well-established standard docking methods (AutoDock Vina and CCDC Gold) with and without an additional post-prediction energy minimisation step using a molecular mechanics force field. We show that both in terms of physical plausibility and the ability to generalise to examples that are distinct from the training data, no deep learning-based method yet outperforms classical docking tools. In addition, we find that molecular mechanics force fields contain docking-relevant physics missing from deep-learning methods. PoseBusters allows practitioners to assess docking and molecular generation methods and may inspire new inductive biases still required to improve deep learning-based methods, which will help drive the development of more accurate and more realistic predictions.

Integrated virtual screening and molecular dynamics simulation approaches revealed potential natural inhibitors for DNMT1 as therapeutic solution for triple negative breast cancer

Triple negative breast cancers (TNBC) are clinically heterogeneous but mostly aggressive malignancies devoid of expression of the estrogen, progesterone, and HER2 (ERBB2 or NEU) receptors. It accounts for 15-20% of all cases. Altered epigenetic regulation including DNA hypermethylation by DNA methyltransferase 1 (DNMT1) has been implicated as one of the causes of TNBC tumorigenesis. The antitumor effect of DNMT1 has also been explored in TNBC that currently lacks targeted therapies. However, the actual treatment for TNBC is yet to be discovered. This study is attributed to the identification of novel drug targets against TNBC. A comprehensive docking and simulation analysis was performed to optimize promising new compounds by estimating their binding affinity to the target protein. Molecular dynamics simulation of 500 ns well complemented the binding affinity of the compound and revealed strong stability of predicted compounds at the docked site. Calculation of binding free energies using MMPBSA and MMGBSA validated the strong binding affinity between compound and binding pockets of DNMT1. In a nutshell, our study uncovered that Beta-Mangostin, Gancaonin Z, 5-hydroxysophoranone, Sophoraflavanone L, and Dorsmanin H showed maximum binding affinity with the active sites of DNMT1. Furthermore, all of these compounds depict maximum drug-like properties. Therefore, the proposed compounds can be a potential candidate for patients with TNBC, but, experimental validation is needed to ensure their safety.Communicated by Ramaswamy H. Sarma.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Tautomer-Specific Deacylation and Ω-Loop Flexibility Explain the Carbapenem-Hydrolyzing Broad-Spectrum Activity of the KPC-2 β-Lactamase

KPC-2 (Klebsiella pneumoniae carbapenemase-2) is a globally disseminated serine-β-lactamase (SBL) responsible for extensive β-lactam antibiotic resistance in Gram-negative pathogens. SBLs inactivate β-lactams via a mechanism involving a hydrolytically labile covalent acyl-enzyme intermediate. Carbapenems, the most potent β-lactams, evade the activity of many SBLs by forming long-lived inhibitory acyl-enzymes; however, carbapenemases such as KPC-2 efficiently deacylate carbapenem acyl-enzymes. We present high-resolution (1.25-1.4 Å) crystal structures of KPC-2 acyl-enzymes with representative penicillins (ampicillin), cephalosporins (cefalothin), and carbapenems (imipenem, meropenem, and ertapenem) obtained utilizing an isosteric deacylation-deficient mutant (E166Q). The mobility of the Ω-loop (residues 165-170) negatively correlates with antibiotic turnover rates (kcat), highlighting the role of this region in positioning catalytic residues for efficient hydrolysis of different β-lactams. Carbapenem-derived acyl-enzyme structures reveal the predominance of the Δ1-(2R) imine rather than the Δ2 enamine tautomer. Quantum mechanics/molecular mechanics molecular dynamics simulations of KPC-2:meropenem acyl-enzyme deacylation used an adaptive string method to differentiate the reactivity of the two isomers. These identify the Δ1-(2R) isomer as having a significantly (7 kcal/mol) higher barrier than the Δ2 tautomer for the (rate-determining) formation of the tetrahedral deacylation intermediate. Deacylation is therefore likely to proceed predominantly from the Δ2, rather than the Δ1-(2R) acyl-enzyme, facilitated by tautomer-specific differences in hydrogen-bonding networks involving the carbapenem C-3 carboxylate and the deacylating water and stabilization by protonated N-4, accumulating a negative charge on the Δ2 enamine-derived oxyanion. Taken together, our data show how the flexible Ω-loop helps confer broad-spectrum activity upon KPC-2, while carbapenemase activity stems from efficient deacylation of the Δ2-enamine acyl-enzyme tautomer.

Apo2ph4: A Versatile Workflow for the Generation of Receptor-based Pharmacophore Models for Virtual Screening

Pharmacophore models are widely used as efficient virtual screening (VS) filters for the target-directed enrichment of large compound libraries. However, the generation of pharmacophore models that have the power to discriminate between active and inactive molecules traditionally requires structural information about ligand-target complexes or at the very least knowledge of one active ligand. The fact that the discovery of the first known active ligand of a newly investigated target represents a major hurdle at the beginning of every drug discovery project underscores the need for methods that are able to derive high-quality pharmacophore models even without the prior knowledge of any active ligand structures. In this work, we introduce a novel workflow, called apo2ph4, that enables the rapid derivation of pharmacophore models solely from the three-dimensional structure of the target receptor. The utility of this workflow is demonstrated retrospectively for the generation of a pharmacophore model for the M2 muscarinic acetylcholine receptor. Furthermore, in order to show the general applicability of apo2ph4, the workflow was employed for all 15 targets of the recently published LIT-PCBA dataset. Pharmacophore-based VS runs using the apo2ph4-derived models achieved a significant enrichment of actives for 13 targets. In the last presented example, a pharmacophore model derived from the etomidate site of the α1β2γ2 GABAA receptor was used in VS campaigns. Subsequent in vitro testing of selected hits revealed that 19 out of 20 (95%) tested compounds were able to significantly enhance GABA currents, which impressively demonstrates the applicability of apo2ph4 for real-world drug design projects.

RCSB Protein Data bank: Tools for visualizing and understanding biological macromolecules in 3D

Now in its 52nd year of continuous operations, the Protein Data Bank (PDB) is the premiere open-access global archive housing three-dimensional (3D) biomolecular structure data. It is jointly managed by the Worldwide Protein Data Bank (wwPDB) partnership. The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB) is funded by the National Science Foundation, National Institutes of Health, and US Department of Energy and serves as the US data center for the wwPDB. RCSB PDB is also responsible for the security of PDB data in its role as wwPDB-designated Archive Keeper. Every year, RCSB PDB serves tens of thousands of depositors of 3D macromolecular structure data (coming from macromolecular crystallography, nuclear magnetic resonance spectroscopy, electron microscopy, and micro-electron diffraction). The RCSB PDB research-focused web portal (RCSB.org) makes PDB data available at no charge and without usage restrictions to many millions of PDB data consumers around the world. The RCSB PDB training, outreach, and education web portal (PDB101.RCSB.org) serves nearly 700 K educators, students, and members of the public worldwide. This invited Tools Issue contribution describes how RCSB PDB (i) is organized; (ii) works with wwPDB partners to process new depositions; (iii) serves as the wwPDB-designated Archive Keeper; (iv) enables exploration and 3D visualization of PDB data via RCSB.org; and (v) supports training, outreach, and education via PDB101.RCSB.org. New tools and features at RCSB.org are presented using examples drawn from high-resolution structural studies of proteins relevant to treatment of human cancers by targeting immune checkpoints.

Protein Data Bank: A Comprehensive Review of 3D Structure Holdings and Worldwide Utilization by Researchers, Educators, and Students

The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB), funded by the United States National Science Foundation, National Institutes of Health, and Department of Energy, supports structural biologists and Protein Data Bank (PDB) data users around the world. The RCSB PDB, a founding member of the Worldwide Protein Data Bank (wwPDB) partnership, serves as the US data center for the global PDB archive housing experimentally-determined three-dimensional (3D) structure data for biological macromolecules. As the wwPDB-designated Archive Keeper, RCSB PDB is also responsible for the security of PDB data and weekly update of the archive. RCSB PDB serves tens of thousands of data depositors (using macromolecular crystallography, nuclear magnetic resonance spectroscopy, electron microscopy, and micro-electron diffraction) annually working on all permanently inhabited continents. RCSB PDB makes PDB data available from its research-focused web portal at no charge and without usage restrictions to many millions of PDB data consumers around the globe. It also provides educators, students, and the general public with an introduction to the PDB and related training materials through its outreach and education-focused web portal. This review article describes growth of the PDB, examines evolution of experimental methods for structure determination viewed through the lens of the PDB archive, and provides a detailed accounting of PDB archival holdings and their utilization by researchers, educators, and students worldwide.

Analysis of protein structures containing HEPES and MES molecules

X‐ray crystallography is the main experimental method behind ligand–macromolecule complexes found in the Protein Data Bank (PDB). Applying bioinformatics methods to such structural data can fuel drug discovery, albeit under the condition that the information is correct. Regrettably, a small number of structures in the PDB are of suboptimal quality due to incorrectly identified and modeled ligands in protein–ligand complexes. In this paper, we combine a theoretical‐graph approach, nuclear density estimates, bioinformatics methods, and prior chemical knowledge to analyze two non‐physiological ligands, HEPES and MES, that are frequent components of crystallization and purifications buffers. Our analysis includes quantum mechanics calculations and Cambridge Structure Database (CSD) queries to define the ideal conformation of these ligands, geometry analysis of PDB deposits regarding several quality factors, and a search for homologous structures to identify other small molecules that could bind in place of the parasitic ligand. Our results highlight the need for careful refinement of macromolecule–ligand complexes and better validation tools that integrate results from all relevant resources.

PDB Code(s): 3K4L, 3PYI, 5T6L, 6BB0, 1PJX, 3O4P, 6WCF, 3DKE, 3E10, 6G38, 4E8R, 4Z91, 3E9F, 1MOS, 1MOQ, 2ESB, 1VHR, 4P66 and 6NNI;