Development of a TSR-Based Method for Protein 3-D Structural Comparison With Its Applications to Protein Classification and Motif Discovery

Development of the TSR-based computational method to investigate spike and monoclonal antibody interactions

Introduction

Monoclonal antibody (mAb) drug treatments have proven effective in reducing COVID-19-related hospitalizations or fatalities, particularly among high-risk patients. Numerous experimental studies have explored the structures of spike proteins and their complexes with ACE2 or mAbs. These 3D structures provide crucial insights into the interactions between spike proteins and ACE2 or mAb, forming a basis for the development of diagnostic tools and therapeutics. However, the field of computational biology has faced substantial challenges due to the lack of methods for precise protein structural comparisons and accurate prediction of molecular interactions. In our previous studies, we introduced the Triangular Spatial Relationship (TSR)-based algorithm, which represents a protein's 3D structure using a vector of integers (keys). These earlier studies, however, were limited to individual proteins.

Purpose

This study introduces new extensions of the TSR-based algorithm, enhancing its ability to study interactions between two molecules. We apply these extensions to gain a mechanistic understanding of spike - mAb interactions.

Method

We expanded the basic TSR method in three novel ways: (1) TSR keys encompassing all atoms, (2) cross keys for interactions between two molecules, and (3) intra-residual keys for amino acids. This TSR-based representation of 3D structures offers a unique advantage by simplifying the search for similar substructures within structural datasets.

Results

The study's key findings include: (i) The method effectively quantified and interpreted conformational changes and steric effects using the newly introduced TSR keys. (ii) Six clusters for CDRH3 and three clusters for CDRL3 were identified using all-atom keys. (iii) We constructed the TSR-STRSUM (TSR-STRucture SUbstitution Matrix), a matrix that represents pairwise similarities between amino acid structures, providing valuable applications in protein sequence and structure comparison. (iv) Intra-residual keys revealed two distinct Tyr clusters characterized by specific triangle geometries.

Conclusion

This study presents an advanced computational approach that not only quantifies and interprets conformational changes in protein backbones, entire structures, or individual amino acids, but also facilitates the search for substructures induced by molecular binding across protein datasets. In some instances, a direct correlation between structures and functions was successfully established.

Development of a TSR-based method for understanding structural relationships of cofactors and local environments in photosystem I

BACKGROUND

All chemical forms of energy and oxygen on Earth are generated via photosynthesis where light energy is converted into redox energy by two photosystems (PS I and PS II). There is an increasing number of PS I 3D structures deposited in the Protein Data Bank (PDB). The Triangular Spatial Relationship (TSR)-based algorithm converts 3D structures into integers (TSR keys). A comprehensive study was conducted, by taking advantage of the PS I 3D structures and the TSR-based algorithm, to answer three questions: (i) Are electron cofactors including P700, A-1 and A0, which are chemically identical chlorophylls, structurally different? (ii) There are two electron transfer chains (A and B branches) in PS I. Are the cofactors on both branches structurally different? (iii) Are the amino acids in cofactor binding sites structurally different from those not in cofactor binding sites?

RESULTS

The key contributions and important findings include: (i) a novel TSR-based method for representing 3D structures of pigments as well as for quantifying pigment structures was developed; (ii) the results revealed that the redox cofactor, P700, are structurally conserved and different from other redox factors. Similar situations were also observed for both A-1 and A0; (iii) the results demonstrated structural differences between A and B branches for the redox cofactors P700, A-1, A0 and A1 as well as their cofactor binding sites; (iv) the tryptophan residues close to A0 and A1 are structurally conserved; (v) The TSR-based method outperforms the Root Mean Square Deviation (RMSD) and the Ultrafast Shape Recognition (USR) methods.

CONCLUSIONS

The structural analyses of redox cofactors and their binding sites provide a foundation for understanding the unique chemical and physical properties of each redox cofactor in PS I, which are essential for modulating the rate and direction of energy and electron transfers.

Development of a novel representation of drug 3D structures and enhancement of the TSR-based method for probing drug and target interactions

Understanding the mechanisms underlying interactions between drugs and target proteins is critical for drug discovery. In our earlier studies, we introduced the Triangular Spatial Relationship (TSR)-based algorithm, which enables the representation of a protein's 3D structure as a vector of integers (TSR keys). These TSR keys correspond to substructures of the 3D structure of a protein and are computed based on the triangles constructed by all possible triples of Cα atoms within the protein. In this study, we report on a new TSR-based algorithm for probing drug and target interactions. Specifically, we have extended the previous algorithm in three novel directions: TSR keys for representing the 3D structure of a drug or a ligand, cross TSR keys between drugs and their targets and intra-residual TSR keys for phosphorylated amino acids. The outcomes illustrate the key contributions as follows: (i) The TSR-based method, which uses the TSR keys as features, is unique in its capability to interpret hierarchical relationships of drugs as well as drug - target complexes using common and specific TSR keys. (ii) The method can distinguish not only the binding sites from the rest of the protein structures, but also the binding sites of primary targets from those of off-targets. (iii) The method has the potential to correlate the 3D structures of drugs with their functions. (iv) Representation of 3D structures by TSR keys has its unique advantage in terms of ease of making searching for similar substructures across structure datasets easier. In summary, this study presents a novel computational methodology, with significant advantages, for providing insights into the mechanism underlying drug and target interactions.

Introducing mirror-image discrimination capability to the TSR-based method for capturing stereo geometry and understanding hierarchical structure relationships of protein receptor family

We have developed a Triangular Spatial Relationship (TSR)-based computational method for protein structure comparison and motif discovery that is both sequence and structure alignment-free. A protein 3D structure is modeled by all possible triangles that are constructed with every three Cα atoms of amino acids as vertices. Every triangle is represented using an integer (a key). The keys are calculated by a rule-based formula which is a function of a representative length, a representative angle, and the vertex labels associated with amino acids. A 3D structure is thereby represented by a vector of integers (TSR keys). Global or local structure comparisons are achieved by computing all keys or a set of keys, respectively. Many enzymatic reactions and notable marketed drugs are highly stereospecific. Thus, in this paper, we propose a modified key calculation formula by including a mechanism for discriminating mirror-image keys to capture stereo geometry. We assign a positive or a negative sign to the integers representing mirror-image keys. Applying the new key calculation function provides the ability to further discriminate mirror-image keys that were previously considered identical. As the result, applying the mirror-image discrimination capability (i) significantly increases the number of distinct keys; (ii) decreases the number of common keys; (iii) decreases structural similarity; (iv) increases the opportunity to identify specific keys for each type of the receptors. The specific keys identified in this study for the cases of without (not applying) and with (applying) mirror-image discrimination can be considered as the structure signatures that exclusively belong to a certain type of receptors. Applying mirror-image discrimination introduces stereospecificity to keys for allowing more precise modeling of ligand - target interactions. The development of mirror-image TSR keys of Cα atom, in conjunction with the integration of Cα TSR keys with all-atom TSR keys for amino acids and drugs, will lead to a new and promising computational method for aiding drug design and discovery.

A study of a hierarchical structure of proteins and ligand binding sites of receptors using the triangular spatial relationship-based structure comparison method and development of a size-filtering feature designed for comparing different sizes of protein structures

The presence of receptors and the specific binding of the ligands determine nearly all cellular responses. Binding of a ligand to its receptor causes conformational changes of the receptor that triggers the subsequent signaling cascade. Therefore, systematically studying structures of receptors will provide insight into their functions. We have developed the triangular spatial relationship (TSR)-based method where all possible triangles are constructed with C_α atoms of a protein as vertices. Every triangle is represented by an integer denoted as a "key" computed through the TSR algorithm. A structure is thereby represented by a vector of integers. In this study, we have first defined substructures using different types of keys. Second, using different types of keys represents a new way to interpret structure hierarchical relations and differences between structures and sequences. Third, we demonstrate the effects of sequence similarity as well as sample size on the structure-based classifications. Fourth, we show identification of structure motifs, and the motifs containing multiple triangles connected by either an edge or a vertex are mapped to the ligand binding sites of the receptors. The structure motifs are valuable resources for the researchers in the field of signal transduction. Next, we propose amino-acid scoring matrices that capture "evolutionary closeness" information based on BLOSUM62 matrix, and present the development of a new visualization method where keys are organized according to evolutionary closeness and shown in a 2D image. This new visualization opens a window for developing tools with the aim of identification of specific and common substructures by scanning pixels and neighboring pixels. Finally, we report a new algorithm called as size filtering that is designed to improve structure comparison of large proteins with small proteins. Collectively, we provide an in-depth interpretation of structure relations through the detailed analyses of different types of keys and their associated key occurrence frequencies, geometries, and labels. In summary, we consider this study as a new computational platform where keys are served as a bridge to connect sequence and structure as well as structure and function for a deep understanding of sequence, structure, and function relationships of the protein family.