Mateusz Banach

Object retrieval systems measure the degree of similarity of the shape of 3D models. They search for the elements of the 3D model databases that resemble the query model. In structural bioinformatics, the query model is a protein tertiary/quaternary structure and the objective is to find similarly shaped molecules in the Protein Data Bank. With the ever-growing size of the PDB, a direct atomic coordinate comparison with all its members is impractical. To overcome this problem, the shape of the molecules can be encoded by fixed-length feature vectors. The distance of a protein to the entire PDB can be measured in this low-dimensional domain in linear time. The state-of-the-art approaches utilize Zernike–Canterakis moments for the shape encoding and supply the retrieval process with geometric data of the input structures. The BioZernike descriptors are a standard utility of the PDB since 2020. However, when trying to calculate the ZC moments locally, the issue of the deficiency of libraries readily available for use in custom programs (i.e., without relying on external binaries) is encountered, in particular programs written in Python. Here, a fast and well-documented Python implementation of the Pozo–Koehl algorithm is presented. In contrast to the more popular algorithm by Novotni and Klein, which is based on the voxelized volume, the PK algorithm produces ZC moments directly from the triangular surface meshes of 3D models. In particular, it can accept the molecular surfaces of proteins as its input. In the presented PK-Zernike library, owing to Numba’s just-in-time compilation, a mesh with 50,000 facets is processed by a single thread in a second at the moment order 20. Since this is the first time the PK algorithm is used in structural bioinformatics, it is employed in a novel, simple, but efficient protein structure retrieval pipeline. The elimination of the outlying chain fragments via a fast PCA-based subroutine improves the discrimination ability, allowing for this pipeline to achieve an 0.961 area under the ROC curve in the BioZernike validation suite (0.997 for the assemblies). The correlation between the results of the proposed approach and of the 3D Surfer program attains values up to 0.99.

In this paper, we present an update to the ellipsoid profile algorithm (EP), a simple technique for the measurement of the globularity of protein structures without the calculation of molecular surfaces. The globularity property is understood in this context as the ability of the molecule to fill a minimum volume enclosing ellipsoid (MVEE) that approximates its assumed globular shape. The more of the interior of this ellipsoid is occupied by the atoms of the protein, the better are its globularity metrics. These metrics are derived from the comparison of the volume of the voxelized representation of the atoms and the volume of all voxels that can fit inside that ellipsoid (a uniform unit Å cube lattice). The so-called ellipsoid profile shows how the globularity changes with the distance from the center. Two of its values, the so-called ellipsoid indexes, are used to classify the structure as globular, semi-globular or non-globular. Here, we enhance the workflow of the EP algorithm via an improved outlier detection subroutine based on principal component analysis. It is capable of robust distinguishing between the dense parts of the molecules and, for example, disordered chain fragments fully exposed to the solvent. The PCA-based method replaces the current approach based on kernel density estimation. The improved EP algorithm was tested on 2124 representatives of domain superfamilies from SCOP 2.08. The second part of this work is dedicated to the survey of globularity of 3594 representatives of biological assemblies from molecules currently deposited in the PDB and analyzed by the 3DComplex database (monomers and complexes up to 60 chains).

In this paper, we show, discuss, and compare the effects of symmetrization in two calculation subroutines of the Fuzzy Oil Drop model, a coarse-grained model of density of hydrophobicity in proteins. In the FOD model, an input structure is enclosed in an axis-aligned ellipsoid called a drop. Two profiles of hydrophobicity are then calculated for its residues: theoretical (based on the 3D Gauss function) and observed (based on pairwise hydrophobic interactions). Condition of the hydrophobic core is revealed by comparing those profiles through relative entropy, while analysis of their local differences allows, in particular, determination of the starting location for the search for protein–protein and protein–ligand interaction areas. Here, we improve the baseline workflow of the FOD model by introducing symmetry to the hydrophobicity profile comparison and ellipsoid bounding procedures. In the first modification (FOD–JS), Kullback–Leibler divergence is enhanced with its Jensen–Shannon variant. In the second modification (FOD-PCA), the molecule is optimally aligned with the axes of the coordinate system via principal component analysis, and the size of its drop is determined by the standard deviation of all its effective atoms, making it less susceptible to structural outliers. Tests on several molecules with various shapes and functions confirm that the proposed modifications improve the accuracy, robustness, speed, and usability of Gauss function-based modeling of the density of hydrophobicity in protein structures.

A computer algorithm for assessment of globularity of protein structures is presented. By enclosing the input protein in a minimum volume ellipsoid (MVEE) and calculating a profile measuring how voxelized space within this shape (cubes on a uniform grid) is occupied by atoms, it is possible to estimate how well the molecule resembles a globule. For any protein to satisfy the proposed globularity criterion, its ellipsoid profile (EP) should first confirm that atoms adequately fill the ellipsoid’s center. This property should then propagate towards the surface of the ellipsoid, although with diminishing importance. It is not required to compute the molecular surface. Globular status (full or partial) is assigned to proteins with values of their ellipsoid profiles, called here the ellipsoid indexes (EI), above certain levels. Due to structural outliers which may considerably distort the measurements, a companion method for their detection and reduction of their influence is also introduced. It is based on kernel density estimation and is shown to work well as an optional input preparation step for MVEE. Finally, the complete workflow is applied to over two thousand representatives of SCOP 2.08 domain superfamilies, surveying the landscape of tertiary structure of proteins from the Protein Data Bank.

Ferredoxin I and II are proteins carrying a specific ligand—an iron-sulfur cluster—which allows transport of electrons. These two classes of ferredoxin in their monomeric and dimeric forms are the object of this work. Characteristic of hydrophobic core in both molecules is analyzed via fuzzy oil drop model (FOD) to show the specificity of their structure enabling the binding of a relatively large ligand and formation of the complex. Structures of FdI and FdII are a promising example for the discussion of influence of hydrophobicity on biological activity but also for an explanation how FOD model can be used as an initial stage adviser (or a scoring function) in the search for locations of ligand binding pockets and protein–protein interaction areas. It is shown that observation of peculiarities in the hydrophobicity distribution present in the molecule (in this case—of a ferredoxin) may provide a promising starting location for computer simulations aimed at the prediction of quaternary structure of proteins.

Protein solubility is based on the compatibility of the specific protein surface with the polar aquatic environment. The exposure of polar residues to the protein surface promotes the protein’s solubility in the polar environment. The aquatic environment also influences the folding process by favoring the centralization of hydrophobic residues with the simultaneous exposure to polar residues. The degree of compatibility of the residue distribution, with the model of the concentration of hydrophobic residues in the center of the molecule, with the simultaneous exposure of polar residues is determined by the sequence of amino acids in the chain. The fuzzy oil drop model enables the quantification of the degree of compatibility of the hydrophobicity distribution observed in the protein to a form fully consistent with the Gaussian 3D function, which expresses an idealized distribution that meets the preferences of the polar water environment. The varied degrees of compatibility of the distribution observed with the idealized one allow the prediction of preferences to interactions with molecules of different polarity, including water molecules in particular. This paper analyzes a set of proteins with different levels of hydrophobicity distribution in the context of the solubility of a given protein and the possibility of complex formation.

The two forms of transthyretin differing slightly in the tertiary structure, despite the presence of five mutations, show radically different properties in terms of susceptibility to the amyloid transformation process. These two forms of transthyretin are the object of analysis. The search for the sources of these differences was carried out by means of a comparative analysis of the structure of these molecules in their native and early intermediate stage forms in the folding process. The criterion for assessing the degree of similarity and differences is the status of the hydrophobic core. The comparison of the level of arrangement of the hydrophobic core and its initial stages is possible thanks to the application of divergence entropy for the early intermediate stage and for the final forms. It was shown that the minimal differences observed in the structure of the hydrophobic core of the forms available in PDB, turned out to be significantly different in the early stage (ES) structure in folding process. The determined values of divergence entropy for both ES forms indicate the presence of the seed of hydrophobic core only in the form resistant to amyloid transformation. In the form of aggressively undergoing amyloid transformation, the structure lacking such a seed is revealed, being a stretched one with a high content of β-type structure. In the discussed case, the active presence of water in the structural transformation of proteins expressed in the fuzzy oil drop model (FOD) is of decisive importance for the generation of the final protein structure. It has been shown that the resistant form tends to generate a centric hydrophobic core with the possibility of creating a globular structure, i.e., a spherical micelle-like form. The aggressively transforming form reveals in the structure of its early intermediate, a tendency to form the ribbon-like micelle as observed in amyloid.

The water environment determines the activity of biological processes. The role of such an environment interpreted in the form of an external field expressed by the 3D Gaussian distribution in the fuzzy oil drop model directs the folding process towards the generation of a centrally located hydrophobic core with the simultaneous exposure of polar residues on the surface. In addition to proteins soluble in the water environment, there is a significant group of membrane proteins that act as receptors or channels, including ion channels in particular. The change of the polar (water) environment into a highly hydrophobic (membrane) environment is quite radical, resulting in a different hydrophobicity distribution within the membrane protein. Modification of the notation of the force field expressing the presence of the hydrophobic environment has been proposed in this work. A modified fuzzy oil drop model with its adaptation to membrane proteins was used to interpret the structure of membrane proteins–mechanosensitive channel. The modified model was also used to describe the so-called negative cases—i.e., for water-soluble proteins with a clear distribution consistent with the fuzzy oil drop model.

In the search for new carriers capable of transporting toxic drugs to a target, particular attention has been devoted to supramolecular systems with a ribbon-like micellar structure of which Congo red is an example. A special promise of the possible use of such systems for directing drugs to a target emerges from their particular affinity to immune complexes and as an independent property, binding many organic compounds including drugs by intercalation. Serum albumin also appeared able to bind micellar particles of such systems. It may protect them against dilution in transport. The mathematical tool, which relies on analysis of the distribution of polarity and hydrophobicity in protein molecules (fuzzy oil drop model), has been used to find the location of binding area in albumin as well as anchorage site for Congo red in heated IgG light chain used as a model presenting immunoglobulin-like structures. Results confirm the suggested formerly binding site of Congo red in V domain of IgG light chain and indicated the cleft between pseudo-symmetric domains of albumin as the area of attachment for the dye.

Research on the protein folding problem differentiates the protein folding process with respect to the duration of this process. The current structure encoded in sequence dogma seems to be clearly justified, especially in the case of proteins referred to as fast-folding, ultra-fast-folding or downhill. In the present work, an attempt to determine the characteristics of this group of proteins using fuzzy oil drop model is undertaken. According to the fuzzy oil drop model, a protein is a specific micelle composed of bi-polar molecules such as amino acids. Protein folding is regarded as a spherical micelle formation process. The presence of covalent peptide bonds between amino acids eliminates the possibility of free mutual arrangement of neighbors. An example would be the construction of co-micelles composed of more than one type of bipolar molecules. In the case of fast folding proteins, the amino acid sequence represents the optimal bipolarity system to generate a spherical micelle. In order to achieve the native form, it is enough to have an external force field provided by the water environment which directs the folding process towards the generation of a centric hydrophobic core. The influence of the external field can be expressed using the 3D Gaussian function which is a mathematical model of the folding process orientation towards the concentration of hydrophobic residues in the center with polar residues exposed on the surface. The set of proteins under study reveals a hydrophobicity distribution compatible with a 3D Gaussian distribution, taken as representing an idealized micelle-like distribution. The structure of the present hydrophobic core is also discussed in relation to the distribution of hydrophobic residues in a partially unfolded form.