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.

Dual regulatory switch confers tighter control on HtrA2 proteolytic activity.

High-temperature requirement protease A2 (HtrA2), a multitasking serine protease that is involved in critical biological functions and pathogenicity, such as apoptosis and cancer, is a potent therapeutic target. It is established that the C-terminal post-synaptic density protein, Drosophila disc large tumor suppressor, zonula occludens-1 protein (PDZ) domain of HtrA2 plays pivotal role in allosteric modulation, substrate binding and activation, as commonly reported in other members of this family. Interestingly, HtrA2 exhibits an additional level of functional modulation through its unique N-terminus, as is evident from 'inhibitor of apoptosis proteins' binding and cleavage. This phenomenon emphasizes multiple activation mechanisms, which so far remain elusive. Using conformational dynamics, binding kinetics and enzymology studies, we addressed this complex behavior with respect to defining its global mode of regulation and activity. Our findings distinctly demonstrate a novel N-terminal ligand-mediated triggering of an allosteric switch essential for transforming HtrA2 to a proteolytically competent state in a PDZ-independent yet synergistic activation process. Dynamic analyses suggested that it occurs through a series of coordinated structural reorganizations at distal regulatory loops (L3, LD, L1), leading to a population shift towards the relaxed conformer. This precise synergistic coordination among different domains might be physiologically relevant to enable tighter control upon HtrA2 activation for fostering its diverse cellular functions. Understanding this complex rheostatic dual switch mechanism offers an opportunity for targeting various disease conditions with tailored site-specific effector molecules.

Boosting virtual screening enrichments with data fusion: coalescing hits from two-dimensional fingerprints, shape, and docking.

Virtual screening is an effective way to find hits in drug discovery, with approaches ranging from fast information-based similarity methods to more computationally intensive physics-based docking methods. However, the best approach to use for a given project is not clear in advance of the screen. In this work, we show that combining results from multiple methods using a standard score (Z-score) can significantly improve virtual screening enrichments over any of the single screening methods. We show that an augmented Z-score, which considers the best two out of three scores for a given compound, outperforms previously published data fusion algorithms. We use three different virtual screening methods (two-dimensional (2D) fingerprint similarity, shape-based similarity, and docking) and study two different databases (DUD and MDDR). The average enrichment in the top 1% was improved by 9% for DUD and 25% for the MDDR, compared with the top individual method. Improvements of 22% for DUD and 43% for MDDR are seen over the average of the three individual methods. Statistics are presented that show a high significance associated with the findings in this work.

Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments.

Structure-based virtual screening plays an important role in drug discovery and complements other screening approaches. In general, protein crystal structures are prepared prior to docking in order to add hydrogen atoms, optimize hydrogen bonds, remove atomic clashes, and perform other operations that are not part of the x-ray crystal structure refinement process. In addition, ligands must be prepared to create 3-dimensional geometries, assign proper bond orders, and generate accessible tautomer and ionization states prior to virtual screening. While the prerequisite for proper system preparation is generally accepted in the field, an extensive study of the preparation steps and their effect on virtual screening enrichments has not been performed. In this work, we systematically explore each of the steps involved in preparing a system for virtual screening. We first explore a large number of parameters using the Glide validation set of 36 crystal structures and 1,000 decoys. We then apply a subset of protocols to the DUD database. We show that database enrichment is improved with proper preparation and that neglecting certain steps of the preparation process produces a systematic degradation in enrichments, which can be large for some targets. We provide examples illustrating the structural changes introduced by the preparation that impact database enrichment. While the work presented here was performed with the Protein Preparation Wizard and Glide, the insights and guidance are expected to be generalizable to structure-based virtual screening with other docking methods.

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.

Rapid shape-based ligand alignment and virtual screening method based on atom/feature-pair similarities and volume overlap scoring.

Shape-based methods for aligning and scoring ligands have proven to be valuable in the field of computer-aided drug design. Here, we describe a new shape-based flexible ligand superposition and virtual screening method, Phase Shape, which is shown to rapidly produce accurate 3D ligand alignments and efficiently enrich actives in virtual screening. We describe the methodology, which is based on the principle of atom distribution triplets to rapidly define trial alignments, followed by refinement of top alignments to maximize the volume overlap. The method can be run in a shape-only mode or it can include atom types or pharmacophore feature encoding, the latter consistently producing the best results for database screening. We apply Phase Shape to flexibly align molecules that bind to the same target and show that the method consistently produces correct alignments when compared with crystal structures. We then illustrate the effectiveness of the method for identifying active compounds in virtual screening of eleven diverse targets. Multiple parameters are explored, including atom typing, query structure conformation, and the database conformer generation protocol. We show that Phase Shape performs well in database screening calculations when compared with other shape-based methods using a common set of actives and decoys from the literature.

Aromatic-Aromatic Interactions Database, A(2)ID: an analysis of aromatic π-networks in proteins.

The geometrical arrangement of the aromatic rings of phenylalanine, tyrosine, tryptophan and histidine has been analyzed at a database level using the X-ray crystal structure of proteins from PDB in order to find out the aromatic-aromatic (π-π) networks in proteins and to understand how these aromatic rings are connected with each-other in a specific π-π network. A stringent examination of the 7848 proteins indicates that close to 89% of the proteins have occurrence of at least a network of 2π or a higher π-π network. The occurrence of π-π networks in various protein superfamilies based on SCOP, CATH and EC classifiers has also been probed in the present work. In general, we find that multidomain and membrane proteins as well as lyases show a more number of these networks. Analysis of the distribution of angle between planes of two proximal aromatic rings (ϕ) distribution indicates that at a larger cutoff distance (between centroid of two aromatic rings), above 5Å, C-H⋯π interactions (T-shaped orientation) are more prevalent, while π-π interactions (stacked orientation) are more prevalent at a smaller cutoff distance. The connectivity patterns of π-π networks propose strong propensity of finding arrangement of aromatic residues as clusters rather than linear arrangement. We have also made a public domain database "Aromatic-Aromatic Interactions Database" (A(2)ID) comprising of all types of π-π networks and their connectivity pattern present in proteins. It can be accessed by url http://203.199.182.73/gnsmmg/databases/aidb/aidb.html.

Protein ligand interaction database (PLID).

A comprehensive database named, protein ligand interaction database (PLID), is created with 6295 ligands bound to proteins extracted from the protein data bank (PDB). This is by far the most comprehensive database of physico-chemical properties, quantum mechanical descriptors and the residues present in the active site of proteins. It is a publicly available web-based database (via the Internet) at http://203.199.182.73/gnsmmg/databases/plid/.

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/.

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.

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.

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 .

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.

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.

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.

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.

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.