Molecular Property Diagnostic Suite for COVID-19 (MPDSCOVID-19): an open-source disease-specific drug discovery portal

Molecular Property Diagnostic Suite (MPDS) was conceived and developed as an open-source disease-specific web portal based on Galaxy. MPDSCOVID-19 was developed for COVID-19 as a one-stop solution for drug discovery research. Galaxy platforms enable the creation of customized workflows connecting various modules in the web server. The architecture of MPDSCOVID-19 effectively employs Galaxy v22.04 features, which are ported on CentOS 7.8 and Python 3.7. MPDSCOVID-19 provides significant updates and the addition of several new tools updated after six years. Tools developed by our group in Perl/Python and open-source tools are collated and integrated into MPDSCOVID-19 using XML scripts. Our MPDS suite aims to facilitate transparent and open innovation. This approach significantly helps bring inclusiveness in the community while promoting free access and participation in software development.

Availability & Implementation

The MPDSCOVID-19 portal can be accessed at https://mpds.neist.res.in:8085/.

Analyzing the cation-aromatic interactions in proteins: Cation-aromatic database V2.0

The cation-aromatic database (CAD) is a comprehensive repository of cation-aromatic motifs found in experimentally determined protein structures, first reported in 2007 [Proteins, 2007, 67, 1179]. The present article is an update of CAD that contains information of approximately 27.26 million cation-aromatic motifs. CAD uses three distance parameters (r, d1, and d2) to determine the position of the cation relative to the centroid of the aromatic residue and classifies the motifs as cation-π or cation-σ interactions. As of June 2023, about 193 936 protein structures were retrieved from Protein Data Bank, and this resulted in the identification of an impressive number of 27 255 817 cation-aromatic motifs. Among these motifs, spherical motifs constituted 94.09%, while cylindrical motifs made up the remaining 5.91%. When considering the interaction of metal ions with aromatic residues, 965 564 motifs are identified. Remarkably, 82.08% of these motifs involved the binding of metal ions to the amino acid HIS. Moreover, the analysis of binding preferences between cations and aromatic residues revealed that the HIS-HIS, PHE-ARG, and TRP-ARG pairs exhibited a preferential geometry. The motif pair HIS-HIS was the most prevalent, accounting for 19.87% of the total, closely followed by TYR-LYS at 10.17%. Conversely, the motif pair TRP-HIS had the lowest occurrence, representing only 4.20% of the total. The data generated help in revealing the characteristics and biological functions of cation-aromatic interactions in biological molecules. The updated version of CAD (Cation-Aromatic Database V2.0) can be accessed at https://acds.neist.res.in/cadv2.

Analyzing the aromatic-aromatic interactions in proteins: A2ID 2.0

The Aromatic-Aromatic Interactions Database (A2ID) is a comprehensive repository dedicated to documenting aromatic-aromatic (π-π) networks observed in experimentally determined protein structures. The first version of A2ID was reported in 2011 [Int J Biol Macromol, 2011, 48, 540]. It has undergone a series of significant updates, leading to its current version, which focuses on the identification and analysis of 3,444,619 π-π networks from proteins. The geometrical parameters such as centroid-centroid distances (r) and interplanar angles (ϕ) were used to identify and characterize π-π networks. It was observed that among the 84,500 proteins with at least one aromatic π-π network, about 92.50 % of the instances are found to be either 2π (77.34 %) or 3π (15.23 %) networks. The analysis of interacting amino acid pairs in 2π networks indicated a dominance of PHE residues followed by TYR. The updated version of A2ID incorporates analysis of π-π networks based on SCOP2 and ECOD classifiers, in addition to the existing SCOP, CATH, and EC classifications. This expanded scope allows researchers to explore the characteristics and functional implications of π-π networks in protein structures from multiple perspectives. The current version of A2ID along with its extensive dataset and detailed geometric information is publicly accessible using https://acds.neist.res.in/a2idv2.

Interactions between Hofmeister anions and the binding pocket of a protein

This paper uses the binding pocket of human carbonic anhydrase II (HCAII, EC 4.2.1.1) as a tool to examine the properties of Hofmeister anions that determine (i) where, and how strongly, they associate with concavities on the surfaces of proteins and (ii) how, upon binding, they alter the structure of water within those concavities. Results from X-ray crystallography and isothermal titration calorimetry show that most anions associate with the binding pocket of HCAII by forming inner-sphere ion pairs with the Zn(2+) cofactor. In these ion pairs, the free energy of anion-Zn(2+) association is inversely proportional to the free energetic cost of anion dehydration; this relationship is consistent with the mechanism of ion pair formation suggested by the "law of matching water affinities". Iodide and bromide anions also associate with a hydrophobic declivity in the wall of the binding pocket. Molecular dynamics simulations suggest that anions, upon associating with Zn(2+), trigger rearrangements of water that extend up to 8 Å away from their surfaces. These findings expand the range of interactions previously thought to occur between ions and proteins by suggesting that (i) weakly hydrated anions can bind complementarily shaped hydrophobic declivities, and that (ii) ion-induced rearrangements of water within protein concavities can (in contrast with similar rearrangements in bulk water) extend well beyond the first hydration shells of the ions that trigger them. This study paints a picture of Hofmeister anions as a set of structurally varied ligands that differ in size, shape, and affinity for water and, thus, in their ability to bind to—and to alter the charge and hydration structure of—polar, nonpolar, and topographically complex concavities on the surfaces of proteins.

Allosteric regulation of serine protease HtrA2 through novel non-canonical substrate binding pocket

HtrA2, a trimeric proapoptotic serine protease is involved in several diseases including cancer and neurodegenerative disorders. Its unique ability to mediate apoptosis via multiple pathways makes it an important therapeutic target. In HtrA2, C-terminal PDZ domain upon substrate binding regulates its functions through coordinated conformational changes the mechanism of which is yet to be elucidated. Although allostery has been found in some of its homologs, it has not been characterized in HtrA2 so far. Here, with an in silico and biochemical approach we have shown that allostery does regulate HtrA2 activity. Our studies identified a novel non-canonical selective binding pocket in HtrA2 which initiates signal propagation to the distal active site through a complex allosteric mechanism. This non-classical binding pocket is unique among HtrA family proteins and thus unfolds a novel mechanism of regulation of HtrA2 activity and hence apoptosis.

From subtle to substantial: role of metal ions on pi-pi interactions

Quantum chemistry calculations reveal that the subtle pi-pi interactions, usually in the range 2-4 kcal/mol, will become substantially significant, from 6 to 17 kcal/mol, in the presence of metal ion. The metal ions have higher affinity toward a pi-pi dimer compared to a single pi-moiety. Considering the widespread occurrence of cation-pi-pi motifs in chemistry and biology, as evident from the database analysis, we propose that the two key noncovalent forces, which govern the macromolecular structure, cation-pi and pi-pi, work in concert.

Cation-aromatic database

Cation-aromatic database (CAD) is a publicly available web-based database that aims to provide further understanding of interaction between a cation and the pi interactions. A tool to identify the interactions in a user-given protein is also added to the database. CAD is freely accessible via the Internet at http://203.199.182.73/gnsmmg/databases/cad/.

Active site acidic residues and structural analysis of modelled human aromatase: a potential drug target for breast cancer

This study sheds new light on the role of acidic residues present in the active site cavity of human aromatase. Eight acidic residues (E129, D222, E245, E302, D309, E379, D380 and D476) lining the cavity are identified and studied using comparative modeling, docking, molecular dynamics as well as statistical techniques. The structural environment of these acidic residues is studied to assess the stability of the corresponding carboxylate anions. Results indicate that the environment of the residues E245, E302 and D222 is most suitable for carboxylate ion formation in the uncomplexed form. However, the stability of D309, D222 and D476 anions is seen to increase on complexation to steroidal substrates. In particular, the interaction between D309 and T310, which assists proton transfer, is found to be formed following androgen/nor-androgen complexation. The residue D309 is found to be clamped in the presence of substrate which is not observed in the case of the other residues although they exhibit changes in properties following substrate binding. Information entropic analysis indicates that the residues D309, D222 and D476 have more conformational flexibility compared to E302 and E245 prior to substrate binding. Interaction similar to that between D476 and D309, which is expected to assist androgen aromatization, is proposed between E302 and E245. The inhibition of aromatase activity by 4-hydroxy androstenedione (formestane) is attributed to a critical hydrogen bond formation between the hydroxy moiety and T310/D309 as well as the large distance from D476. The results corroborate well with earlier site directed mutagenesis studies.

Homology modeling of membrane proteins: A critical assessment

Evaluation and validation of homology modeling protocols are indispensable for membrane proteins as experimental determination of their three-dimensional structure is an arduous task. The prediction ability of Modeller, MOE, InsightII-Homology and Swiss-PdbViewer (SPV) with different sequence alignments CLUSTALW, BLAST and 3D-JIGSAW have been assessed. The sequence identity of the target and template was chosen to be in the range of 25-35%. Validation protocols to assess the structure, fold and stereochemical quality, are employed by comparing with experimental structures. Two different ranking schemes are suggested to evaluate the performance of each methodology based on the validation scores. While unambiguous preference for any given procedure did not surface, statistically Modeller and the sequence alignment technique, 3D-JIGSAW, gave best results amongst the chosen protocols. The present study helps in selecting the right protocols when modeling membrane proteins, which form a major class of drug targets.

DSTHO: database of siRNAs targeted at human oncogenes: a statistical analysis

Existing treatments of human cancer, which is characterized by abnormal proliferation of cells often lead to fatal outcomes. Sequence selective silencing of oncogene expression using siRNA technology is emerging as a potential solution for cancer treatment. The exclusive selectivity and easy application to virtually any therapeutic target including intracellular factors and transcription factors renders siRNA oligonucleotide applications very promising. However, synthesis of siRNA having sufficient knockdown efficiency is laborious and cost intensive. The database is designed in order to aid the synthesis of siRNAs, which target human oncogenes (OsiRNAs). It provides OsiRNAs of known efficacy from previous experiments with links to published literature and theoretically pre-generated putative target sequences. In addition, links to available theoretical tools, databases and literature corresponding to siRNAs in general are also provided. The links to literature provide information about role of siRNA in therapeutics, chemical properties and transfection methods. Statistical analysis of mono-, di- and tri- mers located in OsiRNAs of known efficacies is performed to identify positional preferences and screen specific motifs. This analysis aids the design and synthesis of effective siRNAs, which particularly target human oncogenes. The database can be accessed at .

Proton binding sites and conformational analysis of H+K(+)-ATPase

It is proposed that the hydronium ion, H3O+, binds to the E1 conformation of the alpha-subunit of gastric proton pump. The H3O+ binding cavities are characterized parametrically based on valence, sequence, geometry, and size considerations from comparative modeling. The cavities have scope for accommodating monovalent cations of different ionic radii. The H3O+ transport is proposed to be aided by arenes which are arranged regularly along the pump starting from N-domain through the transmembrane region. Step-by-step structural changes accompanying H3O+ occlusion are studied in detail. The observations corroborate well with earlier experimental studies.

Structural and active site analysis of plasmepsins of Plasmodium falciparum: Potential anti-malarial targets

Comparative protein modeling, active site analysis and binding site specificity for the homologous series of plasmepsins (PM's), present in food vacuole of Plasmodium falciparum, are carried out. Four loops (L1, L2, L3 and L4), which show maximum structural deviations irrespective of type of inhibitor, have been identified. Comparison of the crystal structures of ligand complexes reveal that residues belonging to these loops have negligible coulomb and VDW interactions with the inhibitor but play major role in determining the openness of the binding cavity. The coulomb and VDW interactions between the PMII subsite pockets and inhibitors, which play a major role in determining the inhibition constants, are delineated. Besides small displacements, the catalytic residues D32 of PMII undergoes rotation around the Cgamma-Cbeta single bond to assist catalysis whereas side chain conformational deviations are not observed in D214 on plasmepsin activation. The mutant S79D of PMII (and the corresponding residues of PMI and PMIV) which helps in recognizing and cleaving substrates containing lysine at P1 position is surrounded by highly polar atmosphere stabilized by lysine. However, in PMIII significantly lower polar atmosphere around the mutant A78S/A78D is observed. Large buried side chain area of residues located at M15 and I289 of PMII (and corresponding residues of PMI and PMIV) corroborates well with increase in specificity constant for hydrophobic substrates.

Characterization of calcium and magnesium binding domains of human 5-lipoxygenase

Two calcium binding sites, separated by about 9.3A, present in the loops that connect the beta-sheets of N-terminal domain contain the ligating residues F14, A15, G16, D79, and D18, D19, L76, respectively. Magnesium is found to bind in regions, which are marginally different owing to the disparity in the ionic radii of Ca2+ and Mg2+. The entropy analysis on the loops of 5-lipoxygenase, implementing the wormlike chain model, explains that the N-terminal beta-barrel is well suited to accommodate calcium binding sites. The large buried side chain area of W102 (compared to W13 and W75) and comparatively smaller fraction of side chain exposed to polar atoms corroborate the calcium induced higher affinity to phosphatidylcholine (PC). However, W80 lying in close proximity of the calcium binding sites is expected to have considerable PC affinity but negligible calcium induced effect on PC binding.

Structural and conformational changes concomitant with the E1–E2 transition in H+K+-ATPase: a comparative protein modeling study

Comparative modeling studies on conserved regions of the gastric H(+)K(+)-ATPase reveal that the E1-E2 conformational transition induces significant tertiary structural changes while conserving the secondary structure. The residues 516-530 of the cytoplasmic domain and TM10 within the transmembrane (TM) regions undergo maximum tertiary structural changes. The luminal regions exhibit comparatively lesser tertiary structural deviations. Residues 249-304 show maximum secondary structural deviation in the conformational transition. The Cys-815 and Cys-323 residues involved in inhibitor binding are found to have smaller buried side chain areas in the E1 conformation compared to E2. Retention of activity correlates well with the buried side chain area when selected amino acid residues in TM6 are mutated using modeling techniques with bulkier amino acid residues. Conformational specificity for ion binding is corroborated with the fraction of side chains exposed to polar atoms of the residues E345, D826, V340, A341, V343, and E822.

Trehalose prevents myoglobin collapse and preserves its internal mobility

A quantitative model, which involves diffusion on a temperature-dependent potential, is utilized to analyze the time-dependence of geminate CO recombination to sperm whale myoglobin in a trehalose glass and the accompanying spectral shifts. Most of the recombination is inhomogeneous. This is due to higher geminate reactivity rather than slower protein relaxation. A fraction of the hemes undergoes relaxation with a concomitant increase in the barrier height for recombination. The activation energy for conformational diffusion (relaxation) is considerably lower than in glycerol/water. "Protein collapse", manifested in glycerol/water by a decrease in the equilibrium conformational separation between the bound and deoxy states, is completely prevented in trehalose. We postulate that the high internal viscosity in glycerol/water is due to dehydration of the heme pocket. Trehalose prevents the escape of the few vital internal water molecules and thus preserves the internal lability of the protein. This might be important in understanding the ability of trehalose to protect against the adverse effects of dehydration.