ChannelsDB 2.0: a comprehensive database of protein tunnels and pores in AlphaFold era.

ChannelsDB 2.0 is an updated database providing structural information about the position, geometry and physicochemical properties of protein channels-tunnels and pores-within deposited biomacromolecular structures from PDB and AlphaFoldDB databases. The newly deposited information originated from several sources. Firstly, we included data calculated using a popular CAVER tool to complement the data obtained using original MOLE tool for detection and analysis of protein tunnels and pores. Secondly, we added tunnels starting from cofactors within the AlphaFill database to enlarge the scope of the database to protein models based on Uniprot. This has enlarged available channel annotations ∼4.6 times as of 1 September 2023. The database stores information about geometrical features, e.g. length and radius, and physico-chemical properties based on channel-lining amino acids. The stored data are interlinked with the available UniProt mutation annotation data. ChannelsDB 2.0 provides an excellent resource for deep analysis of the role of biomacromolecular tunnels and pores. The database is available free of charge: https://channelsdb2.biodata.ceitec.cz.

αCharges: partial atomic charges for AlphaFold structures in high quality.

The AlphaFold2 prediction algorithm opened up the possibility of exploring proteins' structural space at an unprecedented scale. Currently, >200 million protein structures predicted by this approach are deposited in AlphaFoldDB, covering entire proteomes of multiple organisms, including humans. Predicted structures are, however, stored without detailed functional annotations describing their chemical behaviour. Partial atomic charges, which map electron distribution over a molecule and provide a clue to its chemical reactivity, are an important example of such data. We introduce the web application αCharges: a tool for the quick calculation of partial atomic charges for protein structures from AlphaFoldDB. The charges are calculated by the recent empirical method SQE+qp, parameterised for this class of molecules using robust quantum mechanics charges (B3LYP/6-31G*/NPA) on PROPKA3 protonated structures. The computed partial atomic charges can be downloaded in common data formats or visualised via the powerful Mol* viewer. The αCharges application is freely available at https://alphacharges.ncbr.muni.cz with no login requirement.

Mol* Volumes and Segmentations: visualization and interpretation of cell imaging data alongside macromolecular structure data and biological annotations.

Segmentation helps interpret imaging data in a biological context. With the development of powerful tools for automated segmentation, public repositories for imaging data have added support for sharing and visualizing segmentations, creating the need for interactive web-based visualization of 3D volume segmentations. To address the ongoing challenge of integrating and visualizing multimodal data, we developed Mol* Volumes and Segmentations (Mol*VS), which enables the interactive, web-based visualization of cellular imaging data supported by macromolecular data and biological annotations. Mol*VS is fully integrated into Mol* Viewer, which is already used for visualization by several public repositories. All EMDB and EMPIAR entries with segmentation datasets are accessible via Mol*VS, which supports the visualization of data from a wide range of electron and light microscopy experiments. Additionally, users can run a local instance of Mol*VS to visualize and share custom datasets in generic or application-specific formats including volumes in .ccp4, .mrc, and .map, and segmentations in EMDB-SFF .hff, Amira .am, iMod .mod, and Segger .seg. Mol*VS is open source and freely available at https://molstarvolseg.ncbr.muni.cz/.

Liposomal Copermeation Assay Reveals Unexpected Membrane Interactions of Commonly Prescribed Drugs.

The permeation of small molecules across biological membranes is a crucial process that lies at the heart of life. Permeation is involved not only in the maintenance of homeostasis at the cell level but also in the absorption and biodistribution of pharmacologically active substances throughout the human body. Membranes are formed by phospholipid bilayers that represent an energy barrier for permeating molecules. Crossing this energy barrier is assumed to be a singular event, and permeation has traditionally been described as a first-order kinetic process, proportional only to the concentration gradient of the permeating substance. For a given membrane composition, permeability was believed to be a unique property dependent only on the permeating molecule itself. We provide experimental evidence that this long-held view might not be entirely correct. Liposomes were used in copermeation experiments with a fluorescent probe, where simultaneous permeation of two substances occurred over a single phospholipid bilayer. Using an assay of six commonly prescribed drugs, we have found that the presence of a copermeant can either enhance or suppress the permeation rate of the probe molecule, often more than 2-fold in each direction. This can have significant consequences for the pharmacokinetics and bioavailability of commonly prescribed drugs when used in combination and provide new insight into so-far unexplained drug-drug interactions as well as changing the perspective on how new drug candidates are evaluated and tested.

OverProt: secondary structure consensus for protein families.

SUMMARY: Every protein family has a set of characteristic secondary structures. However, due to individual variations, a single structure is not enough to represent the whole family. OverProt can create a secondary structure consensus, showing the general fold of the family as well as its variation. Our server provides precomputed results for all CATH superfamilies and user-defined computations, visualized by an interactive viewer, which shows the secondary structure element type, length, frequency of occurrence, spatial variability and ß-connectivity. AVAILABILITY AND IMPLEMENTATION: OverProt Server is freely available at https://overprot.ncbr.muni.cz. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Computational Prodrug Design Methodology for Liposome Formulability Enhancement of Small-Molecule APIs.

Encapsulation into liposomes is a formulation strategy that can improve efficacy and reduce side effects of active pharmaceutical ingredients (APIs) that exhibit poor biodistribution or pharmacokinetics when administered alone. However, many APIs are unsuitable for liposomal formulations intended for parenteral administration due to their inherent physicochemical properties─lipid bilayer permeability and water-lipid equilibrium partitioning coefficient. Too high permeability results in premature leakage from liposomes, while too low permeability means the API is not able to pass across biological barriers. There are several options for solving this issue: (i) change of the lipid bilayer composition, (ii) addition of a permeability enhancer, or (iii) modification of the chemical structure of the API to design a prodrug. The latter approach was taken in the present work, and the effect of small changes in the molecular structure of the API on its permeation rate across a lipidic bilayer was systematically explored utilizing computer simulations. An in silico methodology for prodrug design based on the COSMOperm approach has been proposed and applied to four APIs (abiraterone, cytarabine, 5-fluorouracil, and paliperidone). It is shown that the addition of aliphatic hydrocarbon chains via ester or amide bonds can render the molecule more lipophilic and increase its permeability by approximately 1 order of magnitude for each 2 carbon atoms added, while the formation of fructose adducts can provide a more hydrophilic character to the molecule and reduce its lipid partitioning. While partitioning was found to depend only on the size and type of the added group, permeability was found to depend also on the added group location. Overall, it has been shown that both permeability and lipid partitioning coefficient can be systematically shifted into the desired liposome formulability window by appropriate group contributions to the parental drug. This can significantly increase the portfolio of APIs for which liposome or lipid nanoparticle formulations become feasible.

Expanding the squaramide library as mycobacterial ATP synthase inhibitors: Innovative synthetic pathway and biological evaluation.

Mycobacterial ATP synthase is a validated therapeutic target for combating drug-resistant tuberculosis. Inhibition of this enzyme has been featured as an efficient strategy for the development of new antimycobacterial agents against drug-resistant pathogens. In this study, we synthesised and explored two distinct series of squaric acid analogues designed to inhibit mycobacterial ATP synthase. Among the extensive array of compounds investigated, members of the phenyl-substituted sub-library emerged as primary hits. To gain deeper insights into their mechanisms of action, we conducted advanced biological studies, focusing on the compounds displaying a direct binding of a nitrogen heteroatom to the phenyl ring, resulting in the highest potency. Our investigations into spontaneous mutants led to the validation of a single point mutation within the atpB gene (Rv1304), responsible for encoding the ATP synthase subunit a. This genetic alteration sheds light on the molecular basis of resistance to squaramides. Furthermore, we explored the possibility of synergy between squaramides and the reference drug clofazimine using a checkerboard assay, highlighting the promising avenue for enhancing the effectiveness of existing treatments through combined therapeutic approaches. This study contributes to the expansion of investigating squaramides as promising drug candidates in the ongoing battle against drug-resistant tuberculosis.

Mol* Viewer: modern web app for 3D visualization and analysis of large biomolecular structures.

Large biomolecular structures are being determined experimentally on a daily basis using established techniques such as crystallography and electron microscopy. In addition, emerging integrative or hybrid methods (I/HM) are producing structural models of huge macromolecular machines and assemblies, sometimes containing 100s of millions of non-hydrogen atoms. The performance requirements for visualization and analysis tools delivering these data are increasing rapidly. Significant progress in developing online, web-native three-dimensional (3D) visualization tools was previously accomplished with the introduction of the LiteMol suite and NGL Viewers. Thereafter, Mol* development was jointly initiated by PDBe and RCSB PDB to combine and build on the strengths of LiteMol (developed by PDBe) and NGL (developed by RCSB PDB). The web-native Mol* Viewer enables 3D visualization and streaming of macromolecular coordinate and experimental data, together with capabilities for displaying structure quality, functional, or biological context annotations. High-performance graphics and data management allows users to simultaneously visualise up to hundreds of (superimposed) protein structures, stream molecular dynamics simulation trajectories, render cell-level models, or display huge I/HM structures. It is the primary 3D structure viewer used by PDBe and RCSB PDB. It can be easily integrated into third-party services. Mol* Viewer is open source and freely available at https://molstar.org/.

CATH: increased structural coverage of functional space.

CATH (https://www.cathdb.info) identifies domains in protein structures from wwPDB and classifies these into evolutionary superfamilies, thereby providing structural and functional annotations. There are two levels: CATH-B, a daily snapshot of the latest domain structures and superfamily assignments, and CATH+, with additional derived data, such as predicted sequence domains, and functionally coherent sequence subsets (Functional Families or FunFams). The latest CATH+ release, version 4.3, significantly increases coverage of structural and sequence data, with an addition of 65,351 fully-classified domains structures (+15%), providing 500 238 structural domains, and 151 million predicted sequence domains (+59%) assigned to 5481 superfamilies. The FunFam generation pipeline has been re-engineered to cope with the increased influx of data. Three times more sequences are captured in FunFams, with a concomitant increase in functional purity, information content and structural coverage. FunFam expansion increases the structural annotations provided for experimental GO terms (+59%). We also present CATH-FunVar web-pages displaying variations in protein sequences and their proximity to known or predicted functional sites. We present two case studies (1) putative cancer drivers and (2) SARS-CoV-2 proteins. Finally, we have improved links to and from CATH including SCOP, InterPro, Aquaria and 2DProt.

2DProts: database of family-wide protein secondary structure diagrams.

SUMMARY: Secondary structures provide a deep insight into the protein architecture. They can serve for comparison between individual protein family members. The most straightforward way how to deal with protein secondary structure is its visualization using 2D diagrams. Several software tools for the generation of 2D diagrams were developed. Unfortunately, they create 2D diagrams based on only a single protein. Therefore, 2D diagrams of two proteins from one family markedly differ. For this reason, we developed the 2DProts database, which contains secondary structure 2D diagrams for all domains from the CATH and all proteins from PDB databases. These 2D diagrams are generated based on a whole protein family, and they also consider information about the 3D arrangement of secondary structure elements. Moreover, 2DProts database contains multiple 2D diagrams, which provide an overview of a whole protein family's secondary structures. 2DProts is updated weekly and is integrated into CATH. AVAILABILITY AND IMPLEMENTATION: Freely accessible at https://2dprots.ncbr.muni.cz. The web interface was implemented in JavaScript. The database was implemented in Python. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Atomic Charge Calculator II: web-based tool for the calculation of partial atomic charges.

Partial atomic charges serve as a simple model for the electrostatic distribution of a molecule that drives its interactions with its surroundings. Since partial atomic charges are frequently used in computational chemistry, chemoinformatics and bioinformatics, many computational approaches for calculating them have been introduced. The most applicable are fast and reasonably accurate empirical charge calculation approaches. Here, we introduce Atomic Charge Calculator II (ACC II), a web application that enables the calculation of partial atomic charges via all the main empirical approaches and for all types of molecules. ACC II implements 17 empirical charge calculation methods, including the highly cited (QEq, EEM), the recently published (EQeq, EQeq+C), and the old but still often used (PEOE). ACC II enables the fast calculation of charges even for large macromolecular structures. The web server also offers charge visualization, courtesy of the powerful LiteMol viewer. The calculation setup of ACC II is very straightforward and enables the quick calculation of high-quality partial charges. The application is available at https://acc2.ncbr.muni.cz.

ChannelsDB: database of biomacromolecular tunnels and pores.

ChannelsDB (http://ncbr.muni.cz/ChannelsDB) is a database providing information about the positions, geometry and physicochemical properties of channels (pores and tunnels) found within biomacromolecular structures deposited in the Protein Data Bank. Channels were deposited from two sources; from literature using manual deposition and from a software tool automatically detecting tunnels leading to the enzymatic active sites and selected cofactors, and transmembrane pores. The database stores information about geometrical features (e.g. length and radius profile along a channel) and physicochemical properties involving polarity, hydrophobicity, hydropathy, charge and mutability. The stored data are interlinked with available UniProt annotation data mapping known mutation effects to channel-lining residues. All structures with channels are displayed in a clear interactive manner, further facilitating data manipulation and interpretation. As such, ChannelsDB provides an invaluable resource for research related to deciphering the biological function of biomacromolecular channels.

MOLEonline: a web-based tool for analyzing channels, tunnels and pores (2018 update).

MOLEonline is an interactive, web-based application for the detection and characterization of channels (pores and tunnels) within biomacromolecular structures. The updated version of MOLEonline overcomes limitations of the previous version by incorporating the recently developed LiteMol Viewer visualization engine and providing a simple, fully interactive user experience. The application enables two modes of calculation: one is dedicated to the analysis of channels while the other was specifically designed for transmembrane pores. As the application can use both PDB and mmCIF formats, it can be leveraged to analyze a wide spectrum of biomacromolecular structures, e.g. stemming from NMR, X-ray and cryo-EM techniques. The tool is interconnected with other bioinformatics tools (e.g., PDBe, CSA, ChannelsDB, OPM, UniProt) to help both setup and the analysis of acquired results. MOLEonline provides unprecedented analytics for the detection and structural characterization of channels, as well as information about their numerous physicochemical features. Here we present the application of MOLEonline for structural analyses of α-hemolysin and transient receptor potential mucolipin 1 (TRMP1) pores. The MOLEonline application is freely available via the Internet at https://mole.upol.cz.

Tools and data services registry: a community effort to document bioinformatics resources.

Life sciences are yielding huge data sets that underpin scientific discoveries fundamental to improvement in human health, agriculture and the environment. In support of these discoveries, a plethora of databases and tools are deployed, in technically complex and diverse implementations, across a spectrum of scientific disciplines. The corpus of documentation of these resources is fragmented across the Web, with much redundancy, and has lacked a common standard of information. The outcome is that scientists must often struggle to find, understand, compare and use the best resources for the task at hand.Here we present a community-driven curation effort, supported by ELIXIR-the European infrastructure for biological information-that aspires to a comprehensive and consistent registry of information about bioinformatics resources. The sustainable upkeep of this Tools and Data Services Registry is assured by a curation effort driven by and tailored to local needs, and shared amongst a network of engaged partners.As of November 2015, the registry includes 1785 resources, with depositions from 126 individual registrations including 52 institutional providers and 74 individuals. With community support, the registry can become a standard for dissemination of information about bioinformatics resources: we welcome everyone to join us in this common endeavour. The registry is freely available at https://bio.tools.

RH421 binds into the ATP-binding site on the Na+/K+-ATPase.

The Na+/K+-ATPase plays a key role in ion transport across the plasma membrane of all animal cells. The voltage-sensitive styrylpyrimidium dye RH421 has been used in several laboratories for monitoring of Na+/K+-ATPase kinetics. It is known, that RH421 can interact with the enzyme and it can influence its activity at micromolar concentrations, but structural details of this interaction are only poorly understood. Experiments with isolated large cytoplasmic loop (C45) of Na+/K+-ATPase revealed that RH421 can interact with this part of the protein with dissociation constant 1µM. The Trp-to-RH421 FRET performed on six single-tryptophan mutants revealed that RH421 binds directly into the ATP-binding site. This conclusion was further supported by results from molecular docking, site-directed mutagenesis and by competitive experiments using ATP. Experiments with C45/DPPC mixture revealed that RH421 can bind to both C45 and lipids, but only the former interaction was influenced by the presence of ATP.

Molecular insights into the role of a distal F240A mutation that alters CYP1A1 activity towards persistent organic pollutants.

BACKGROUND: Cytochromes P450 are major drug-metabolizing enzymes involved in the biotransformation of diverse xenobiotics and endogenous chemicals. Persistent organic pollutants (POPs) are toxic hydrophobic compounds that cause serious environmental problems because of their poor degradability. This calls for rational design of enzymes capable of catalyzing their biotransformation. Cytochrome P450 1A1 isoforms catalyze the biotransformation of some POPs, and constitute good starting points for the design of biocatalysts with tailored substrate specificity. METHODS: We rationalized the activities of wild type and mutant forms of rat cytochrome P450 1A1 towards 2,3,7,8-tetrachloro-dibenzo-p-dioxin (TCDD) and 3,3',4,4'-tetrachlorobiphenyl (PCB77) using experiments and molecular dynamics simulations. RESULTS: We showed that the enhanced activity of the CYP1A1 mutant towards TCDD was due to more efficient binding of the substrate in the active site even though the mutated site was over 2.5nm away from the catalytic center. Moreover, this mutation reduced activity towards PCB77. GENERAL SIGNIFICANCE: Amino acids that affect substrate access channels can be viable targets for rational enzyme design even if they are located far from the catalytic site.

Novel thidiazuron-derived inhibitors of cytokinin oxidase/dehydrogenase.

KEY MESSAGE: Two new TDZ derivatives (HETDZ and 3FMTDZ) are very potent inhibitors of CKX and are promising candidates for in vivo studies. Cytokinin hormones regulate a wide range of essential processes in plants. Thidiazuron (N-phenyl-N'-1,2,3-thiadiazol-5-yl urea, TDZ), formerly registered as a cotton defoliant, is a well known inhibitor of cytokinin oxidase/dehydrogenase (CKX), an enzyme catalyzing the degradation of cytokinins. TDZ thus increases the lifetime of cytokinins and their effects in plants. We used in silico modeling to design, synthesize and characterize twenty new TDZ derivatives with improved inhibitory properties. Two compounds, namely 1-[1,2,3]thiadiazol-5-yl-3-(3-trifluoromethoxy-phenyl)urea (3FMTDZ) and 1-[2-(2-hydroxyethyl)phenyl]-3-(1,2,3-thiadiazol-5-yl)urea (HETDZ), displayed up to 15-fold lower IC 50 values compared with TDZ for AtCKX2 from Arabidopsis thaliana and ZmCKX1 and ZmCKX4a from Zea mays. Binding modes of 3FMTDZ and HETDZ were analyzed by X-ray crystallography. Crystal structure complexes, solved at 2.0 Å resolution, revealed that HETDZ and 3FMTDZ bound differently in the active site of ZmCKX4a: the thiadiazolyl ring of 3FMTDZ was positioned over the isoalloxazine ring of FAD, whereas that of HETDZ had the opposite orientation, pointing toward the entrance of the active site. The compounds were further tested for cytokinin activity in several cytokinin bioassays. We suggest that the combination of simple synthesis, lowered cytokinin activity, and enhanced inhibitory effects on CKX isoforms, makes 3FMTDZ and HETDZ suitable candidates for in vivo studies.

Exponential repulsion improves structural predictability of molecular docking.

Molecular docking is a powerful tool for theoretical prediction of the preferred conformation and orientation of small molecules within protein active sites. The obtained poses can be used for estimation of binding energies, which indicate the inhibition effect of designed inhibitors, and therefore might be used for in silico drug design. However, the evaluation of ligand binding affinity critically depends on successful prediction of the native binding mode. Contemporary docking methods are often based on scoring functions derived from molecular mechanical potentials. In such potentials, nonbonded interactions are typically represented by electrostatic interactions between atom-centered partial charges and standard 6-12 Lennard-Jones potential. Here, we present implementation and testing of a scoring function based on more physically justified exponential repulsion instead of the standard Lennard-Jones potential. We found that this scoring function significantly improved prediction of the native binding modes in proteins bearing narrow active sites such as serine proteases and kinases. © 2016 Wiley Periodicals, Inc.

The Role of Protein-Protein and Protein-Membrane Interactions on P450 Function.

This symposium summary, sponsored by the ASPET, was held at Experimental Biology 2015 on March 29, 2015, in Boston, Massachusetts. The symposium focused on: 1) the interactions of cytochrome P450s (P450s) with their redox partners; and 2) the role of the lipid membrane in their orientation and stabilization. Two presentations discussed the interactions of P450s with NADPH-P450 reductase (CPR) and cytochrome b5. First, solution nuclear magnetic resonance was used to compare the protein interactions that facilitated either the hydroxylase or lyase activities of CYP17A1. The lyase interaction was stimulated by the presence of b5 and 17α-hydroxypregnenolone, whereas the hydroxylase reaction was predominant in the absence of b5. The role of b5 was also shown in vivo by selective hepatic knockout of b5 from mice expressing CYP3A4 and CYP2D6; the lack of b5 caused a decrease in the clearance of several substrates. The role of the membrane on P450 orientation was examined using computational methods, showing that the proximal region of the P450 molecule faced the aqueous phase. The distal region, containing the substrate-access channel, was associated with the membrane. The interaction of NADPH-P450 reductase (CPR) with the membrane was also described, showing the ability of CPR to "helicopter" above the membrane. Finally, the endoplasmic reticulum (ER) was shown to be heterogeneous, having ordered membrane regions containing cholesterol and more disordered regions. Interestingly, two closely related P450s, CYP1A1 and CYP1A2, resided in different regions of the ER. The structural characteristics of their localization were examined. These studies emphasize the importance of P450 protein organization to their function.

In silico pharmacology: Drug membrane partitioning and crossing.

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; ß-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.