3dSS: 3D structural superposition

Semaphorin heterodimerization in cis regulates membrane targeting and neocortical wiring

Disruption of neocortical circuitry and architecture in humans causes numerous neurodevelopmental disorders. Neocortical cytoarchitecture is orchestrated by various transcription factors such as Satb2 that control target genes during strict time windows. In humans, mutations of SATB2 cause SATB2 Associated Syndrome (SAS), a multisymptomatic syndrome involving epilepsy, intellectual disability, speech delay, and craniofacial defects. Here we show that Satb2 controls neuronal migration and callosal axonal outgrowth during murine neocortical development by inducing the expression of the GPI-anchored protein, Semaphorin 7A (Sema7A). We find that Sema7A exerts this biological activity by heterodimerizing in cis with the transmembrane semaphorin, Sema4D. We could also observe that heterodimerization with Sema7A promotes targeting of Sema4D to the plasma membrane in vitro. Finally, we report an epilepsy-associated de novo mutation in Sema4D (Q497P) that inhibits normal glycosylation and plasma membrane localization of Sema4D-associated complexes. These results suggest that neuronal use of semaphorins during neocortical development is heteromeric, and a greater signaling complexity exists than was previously thought.

Structural and thermodynamic properties of conserved water molecules in Mpro native: A combined approach by MD simulation and Grid Inhomogeneous Solvation Theory

The new viral strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are continuously rising, becoming more virulent, and transmissible. Therefore, the development of new antiviral drugs is essential. Due to its significant role in the viral life cycle of SARS-CoV-2, the main protease (Mpro) enzyme is a leading target for antiviral drug design. The Mpro monomer consists of domain DI, DII, and DI-DII interface. Twenty-one conserved water molecules (W4-W24) are occupied at these domains according to multiple crystal structure analyses. The crystal and MD structures reveal the presence of eight conserved water sites in domain DI, DII and remaining in the DI-DII interface. Grid-based inhomogeneous fluid solvation theory (GIST) was employed on MD structures of Mpro native to predict structural and thermodynamic properties of each conserved water site for focusing to identify the specific conserved water molecules that can easily be displaced by proposed ligands. Finally, MD water W13 is emerged as a promising candidate for water mimic drug design due to its low mean interaction energy, loose binding character with the protein, and its involvement in a water-mediated H-bond with catalytic His41 via the interaction Thr25(OG)---W13---W---His41(NE2). In this context, water occupancy, relative interaction energy, entropy, and topologies of W13 are thermodynamically acceptable for the water displacement method. Therefore, the strategic use of W13's geometrical position in the DI domain may be implemented for drug discovery against COVID disease by designing new ligands with appropriately oriented chemical groups to mimic its structural, electronic, and thermodynamic properties.

Distinct characteristics of putative archaeal 5-methylcytosine RNA methyltransferases unveil their substrate specificities and evolutionary ancestries

5-Methylcytosine methyltransferases (m5C MTases) are known to be involved in the modification of RNA. Although these enzymes have been relatively well characterized in bacteria and eukarya, a complete understanding of the archaeal counterparts is lacking. In this study, the identification and characterization of archaeal RNA m5C MTases were performed. As a case study, a hyperthermophilic archaeon, Pyrococcus horikoshii OT3, which possesses five putative RNA m5C MTases, was chosen. Among the five putative RNA m5C MTases, two proteins (PH0851 and PH1991) have been characterized as homologs of a bacterial rRNA MTase (RsmB) and eukaryal tRNA MTase (NSUN6), respectively. The in-depth characterization of the remaining three putative RNA m5C MTases (PH1078, PH1374, and PH1537) in this study suggests the presence of the signature architecture and catalytic residues plausibly involved in the binding of their cognate RNA substrates. Additionally, the results also suggest the existence of two RsmB-like proteins (PH0851 and PH1078) belonging to the same subfamily IV of m5C RNA MTase. However, the proteins PH1374 and PH1537 belong to the same subfamily V but bind to different substrates, rRNA and tRNA, respectively. The findings further indicate that archaeal RNA m5C MTases link those from bacteria and eukarya.Communicated by Ramaswamy H. Sarma.

Proteomic analysis of Cryptostegia grandiflora latex, purification, characterization, and biological activity of two osmotin isoforms

Although latex fluids are found in >20,000 plant species, the biochemical composition and biological function of their proteins are still poorly explored. Thus, this work aimed to conduct a proteomic analysis of Cryptostegia grandiflora latex (CgLP) for subsequent purification and characterization of an antifungal protein. After 2D-SDS-PAGE and mass spectrometry, 27 proteins were identified in CgLP, including a polygalacturonase inhibitor, cysteine peptidases, pathogenesis-related proteins (PR-4), and osmotins. Then, two osmotin isoforms (CgOsm) were purified, and a unique N-terminal sequence was determined (1ATFDIRSNCPYTVWAAAVPGGGRRLDRGQTWTINVAPGTA40). The PCR products revealed a cDNA sequence of 609 nucleotides for CgOsm, which encoded a polypeptide with 203 amino acid residues. The structure of CgOsm has features of typical osmotin or thaumatin-like proteins (TLPs), such as 16 conserved Cys residues, REDDD and FF motifs, an acidic cleft, and three main domains. Atomic force microscopy (AFM) and bioinformatics suggested that CgOsm is associated with three chain units. This result was interesting since the literature describes osmotins and TLPs as monomers. AFM also showed that Fusarium falciforme spores treated with CgOsm were drastically damaged. Therefore, it is speculated that CgOsm forms pores in the membrane of these cells, causing the leakage of cytoplasmic content.

Functional characterization of a hypothetical protein (TTHA1873) from Thermus thermophilus

Thermus thermophilus is an extremely thermophilic organism that thrives at a temperature of 65°C. T. thermophilus genome has ~2218 genes, out of which 66% (1482 genes) have been annotated, and the remaining 34% (736 genes) are assigned as hypothetical proteins. In this work, biochemical and biophysical experiments were performed to characterize the hypothetical protein TTHA1873 from T. thermophilus. The hypothetical protein TTHA1873 acts as a nuclease, which indiscreetly cuts methylated and non-methylated DNA in divalent metal ions and relaxes the plasmid DNA in the presence of ATP. The chelation of metal ions with EDTA inhibits its activity. These results suggest that the hypothetical protein TTHA1873 would be a CRISPR-associated protein with non-specific DNase activity and ATP-dependent DNA-relaxing activity.

Bioinformatics Approaches Applied to the Discovery of Antifungal Peptides

Antifungal peptides (AFPs) comprise a group of substances with a broad spectrum of activities and complex action mechanisms. They develop in nature via an evolutionary process resulting from the interactions between hosts and pathogens. The AFP database is experimentally verified and curated from research articles, patents, and public databases. In this review, we compile information about the primary databases and bioinformatics tools that have been used in the discovery of AFPs during the last 15 years. We focus on the classification and prediction of AFPs using different physicochemical properties, such as polarity, hydrophobicity, hydrophilicity, mass, acidic, basic, and isoelectric indices, and other structural properties. Another method for discovering AFPs is the implementation of a peptidomic approach and bioinformatics filtering, which gave rise to a new family of peptides that exhibit a broad spectrum of antimicrobial activity against Candida albicans with low hemolytic effects. The application of machine intelligence in the sphere of biological sciences has led to the development of automated tools. The progress made in this area has also paved the way for producing new drugs more quickly and effectively. However, we also identified that further advancements are still needed to complete the AFP libraries.

Structural analysis and ensemble docking revealed the binding modes of selected progesterone receptor modulators

Uterine fibroids (UF) are benign smooth muscle neoplasm of uterus that have a significant impact on a woman's quality of life as they perturb hormonal homeostasis resulting in heavy menstrual bleeding, impaired fertility, pregnancy complications and loss. UF can be surgically removed through invasive procedures, but their recurrence rate is often high. Progesterone receptor (PR) has an imperative role in UF management. Mifepristone, ulipristal acetate (UPA) and asoprisnil (ASO) are some selective progesterone receptor modulators (SPRMs), acts on PR, but due to their side effects in long term use, they were withdrawn from the market. Hence, there is a dire need for novel, highly efficient with least side effects, therapeutics for the treatment of UF. To contribute toward the drug discovery for UF, we made an extensive structural comparison of reported PR crystal structures, also elucidated the binding modes of four existing SPRMs against human PR through ensemble docking approach. Our studies revealed the presence of 5 highly repeating water molecules that has an important role in ligand binding and structural stability. Our ensemble docking and MD simulation revealed that studied ligands have preferential selectivity toward the specific conformation of PR. It is anticipated that our study will be a useful resource to all the drug discovery scientists who are engaged in the identification of lead molecules against UF.

Expression and functional analysis of various structural domains of tobacco topoisomerase II: To understand the mechanistic insights of plant type II topoisomerases

In contrast to bacterial, yeast and animal systems, topoisomerases (topo) from plants have not been well studied. In this report, we generated four truncated topoisomerase II (Topo II) cDNA fragments encoding different functional domains of Nicotiana tabacum topo II (NtTopoII). Each of these recombinant polypeptides was expressed alone or in combination in temperature-sensitive topoisomerase II yeast mutants. Recombinant NtTopoII with truncated polypeptides fails to target the yeast nuclei and does not rescue the temperature-sensitive phenotype. In contrast complementation was achieved with the full-length NtTopoII, which localized to the yeast nucleus. These observations suggested the presence of a potent nuclear localization signal (NLS) in the extreme C-terminal 314 amino acid residues of NtTopoII that functioned effectively in the heterologous yeast system. Biochemical characterization of purified recombinant full-length and the partial NtTopoII polypeptides revealed that the ATP-binding and hydrolysis region of NtTopoIIwas located at 413 amino acid N-terminal region and this ATPase domain is functional both when it is expressed as a separate polypeptide or as part of the holoenzyme. The present findings also revealed that all NtTopoII truncated polypeptides were detrimental for in vitro supercoiled DNA relaxation and/or DNA nicking and ligation activity. Further, we discuss the possible disruption of coordinated macromolecular interface movements and the dimer interactions in truncated NtTopoII that are required for functional topoisomerase activity.

In silico testing of flavonoids as potential inhibitors of protease and helicase domains of dengue and Zika viruses

Background

Dengue and Zika are two major vector-borne diseases. Dengue causes up to 25,000 deaths and nearly a 100 million cases worldwide per year, while the incidence of Zika has increased in recent years. Although Zika has been associated to fetal microcephaly and Guillain-Barré syndrome both it and dengue have common clinical symptoms such as severe headache, retroocular pain, muscle and join pain, nausea, vomiting, and rash. Currently, vaccines have been designed and antivirals have been identified for these diseases but there still need for more options for treatment. Our group previously obtained some fractions from medicinal plants that blocked dengue virus (DENV) infection in vitro. In the present work, we explored the possible targets by molecular docking a group of molecules contained in the plant fractions against DENV and Zika virus (ZIKV) NS3-helicase (NS3-hel) and NS3-protease (NS3-pro) structures. Finally, the best ligands were evaluated by molecular dynamic simulations.

Methods

To establish if these molecules could act as wide spectrum inhibitors, we used structures from four DENV serotypes and from ZIKV. ADFR 1.2 rc1 software was used for docking analysis; subsequently molecular dynamics analysis was carried out using AMBER20.

Results

Docking suggested that 3,5-dicaffeoylquinic acid (DCA01), quercetin 3-rutinoside (QNR05) and quercetin 3,7-diglucoside (QND10) can tightly bind to both NS3-hel and NS3-pro. However, after a molecular dynamics analysis, tight binding was not maintained for NS3-hel. In contrast, NS3-pro from two dengue serotypes, DENV3 and DENV4, retained both QNR05 and QND10 which converged near the catalytic site. After the molecular dynamics analysis, both ligands presented a stable trajectory over time, in contrast to DCA01. These findings allowed us to work on the design of a molecule called MOD10, using the QND10 skeleton to improve the interaction in the active site of the NS3-pro domain, which was verified through molecular dynamics simulation, turning out to be better than QNR05 and QND10, both in interaction and in the trajectory.

Discussion

Our results suggests that NS3-hel RNA empty binding site is not a good target for drug design as the binding site located through docking is too big. However, our results indicate that QNR05 and QND10 could block NS3-pro activity in DENV and ZIKV. In the interaction with these molecules, the sub-pocket-2 remained unoccupied in NS3-pro, leaving opportunity for improvement and drug design using the quercetin scaffold. The analysis of the NS3-pro in complex with MOD10 show a molecule that exerts contact with sub-pockets S1, S1', S2 and S3, increasing its affinity and apparent stability on NS3-pro.

Thermally stable and acidic pH tolerant mutant phytases with high catalytic efficiency from Yersinia intermedia for potential application in feed industries

Heat- and pH-stable phytase efficiently hydrolyzes phytic acid. In this research, heat- and pH-stable mutant phytases, T⁸³R, L²⁸⁷R, and T⁸³R/L²⁸⁷R were generated by site-directed mutagenesis from Yersinia intermedia. After the induction and expression of recombinant wild-type and mutant phytases in E. coli BL21, the enzymes were purified using nickel sepharose affinity chromatography, and characterized kinetically and thermodynamically using spectroscopy methods. The mutants showed optimum activity at pH 5.15 and 55-61 °C. The catalytic efficiencies of T⁸³R, L²⁸⁷R, T⁸³R/L²⁸⁷R, and wild-type phytases were calculated to be 2941, 29346, 4906, and 6917 mmol/L^-1s^-1, respectively. Moreover, after the incubation of T⁸³R, L²⁸⁷R, wild-type, and T⁸³R/ L²⁸⁷R phytases at 100 °C for 1 h, the enzymes retained 22, 5, 4, and 2% of their initial activities, respectively. In addition, T⁸³R, T⁸³R/L²⁸⁷R, L²⁸⁷R, and wild-type phytases retained 82, 44, 16 as well as 11% of their initial activities after 1 h at pH 5.15, respectively. Among these mutants, T⁸³R mutant showed 18% increase in thermal stability, 71% increase in pH stability, and +0.103 KJ/mole increase in ΔΔG, while the catalytic efficiency and ΔΔG value of L²⁸⁷R mutant increased by 4 times and +0.0903 KJ/mole, respectively. Thus, the mutants have the potential to be used in feed industries to increase the bioavailability of minerals while decreasing soil and water pollution.

Investigation of structural analogs of hydroxychloroquine for SARS-CoV-2 main protease (Mpro): A computational drug discovery study

The main protease (Mpro) is the key enzyme of nCOVID-19 and plays a decisive role that makes it an attractive drug target. Multiple analysis of crystal structures reveals the presence of W1, W2, and W3 water locations in the active site pocket of Mpro; W1 and W2 are unstable and are weakly bonded with protein in comparison to W3 of Mpro-native. So, we adopt the water displacement method to occupy W1 or W2 sites by triggering HCQ or its analogs to inactivate the enzyme. Virtual screening is employed to find out best analogs of HCQ, molecular docking is used for water displacement from catalytic region of Mpro, and finally, MD simulations are conducted for validation of these findings. The docking study reveals that W1 and W2 are occupied by respective atoms of ZINC28706440 whereas W2 by HCQ and indacaterol. Finally, MD results demonstrate (i) HCQ occupies W1 and W2 positions, but its analogs (indacaterol and ZINC28706440) are inadequate to retain either W1 or W2 (ii) His41 and Asp187 are stabilized by W3 in Mpro-native and His41, Cys145 and HCQ by W7 in ZINC28706440, and W4, W5, and W6 make water mediated bridge between indacaterol with His41. The structural, dynamical, and thermodynamic (WFP and J value) profiling parameters suggest that W3, W4, and W7 are prominent in their corresponding positions in comparison with W5 and W6. The final results conclude that ZINC28706440 may act as a best analog of HCQ with acceptable physico-chemical and toxicological scores and may further be synthesized for experimental validation.

Conserved features of the MlaD domain aid the trafficking of hydrophobic molecules

In Gram-negative bacteria, the maintenance of lipid asymmetry (Mla) system is involved in the transport of phospholipids between the inner (IM) and outer membrane. The Mla system utilizes a unique IM-associated periplasmic solute-binding protein, MlaD, which possesses a conserved domain, MlaD domain. While proteins carrying the MlaD domain are known to be primarily involved in the trafficking of hydrophobic molecules, not much is known about this domain itself. Thus, in this study, the characterization of the MlaD domain employing bioinformatics analysis is reported. The profiling of the MlaD domain of different architectures reveals the abundance of glycine and hydrophobic residues and the lack of cysteine residues. The domain possesses a conserved N-terminal region and a well-preserved glycine residue that constitutes a consensus motif across different architectures. Phylogenetic analysis shows that the MlaD domain archetypes are evolutionarily closer and marked by the conservation of a functionally crucial pore loop located at the C-terminal region. The study also establishes the critical role of the domain-associated permeases and the driving forces governing the transport of hydrophobic molecules. This sheds sufficient light on the structure-function-evolutionary relationship of MlaD domain. The hexameric interface analysis reveals that the MlaD domain itself is not a sole player in the oligomerization of the proteins. Further, an operonic and interactome map analysis reveals that the Mla and the Mce systems are dependent on the structural homologs of the nuclear transport factor 2 superfamily.

Identification and characterization of metal uptake ABC transporters in Mycobacterium tuberculosis unveil their ligand specificity

Mycobacterium tuberculosis, one of the major threats to mankind, requires micronutrients like metal ions for their survival and pathogenicity inside the host system. Intracellular pathogens such as M. tuberculosis have co-evolved to combat the nutritional immunity developed by the host. It has developed eminent mechanisms to sequester essential metal ions from the host system. One such prominent mechanism to scavenge metal ions to thrive in the host cell involves ATP-binding cassette (ABC) transporters, which transport metal ions (in free and/or complex forms) across the cell membrane. This study employs a high-throughput data mining analysis to identify open reading frames (ORFs) encoding metal uptake ABC transporters in M. tuberculosis H37Rv. In total, 19 ORFs resulting in seven ABC transport systems and two P-type ATPases were identified, which are potentially involved in the uptake of different metal ions. The results also suggest the existence of a subunit sharing mechanism in M. tuberculosis where the transmembrane and nucleotide binding domains are shared among different ABC transport systems indicating the import of multiple substrates via a single ABC transporter. Thus, this study reflects an overview of the repertoire of metal-specific ABC transport systems in M. tuberculosis H37Rv, providing potential therapeutic targets for the future.

PilB from Streptococcus sanguinis is a bimodular type IV pilin with a direct role in adhesion

Type IV pili (T4P) are functionally versatile filaments widespread in prokaryotes, composed of type IV pilins and assembled by conserved multiprotein machineries. It remains unclear how such rather simple filaments can be so versatile. Our structure/function analysis of PilB, a minor pilin of Streptococcus sanguinis T4P, offers an elegant explanation for this paradox. We show that PilB is a modular pilin with a bulky module “grafted” onto a small pilin module, which directly mediates adhesion of S. sanguinis to host cells/proteins. This evolutionary tinkering strategy appears to be prevalent in bacteria since a global analysis reveals that modular pilins are widespread and exhibit an astonishing variety of architectures.

Type IV pili (T4P) are functionally versatile filamentous nanomachines, nearly ubiquitous in prokaryotes. They are predominantly polymers of one major pilin but also contain minor pilins whose functions are often poorly defined and likely to be diverse. Here, we show that the minor pilin PilB from the T4P of Streptococcus sanguinis displays an unusual bimodular three-dimensional structure with a bulky von Willebrand factor A–like (vWA) module “grafted” onto a small pilin module via a short loop. Structural modeling suggests that PilB is only compatible with a localization at the tip of T4P. By performing a detailed functional analysis, we found that 1) the vWA module contains a canonical metal ion–dependent adhesion site, preferentially binding Mg²⁺ and Mn²⁺, 2) abolishing metal binding has no impact on the structure of PilB or piliation, 3) metal binding is important for S. sanguinis T4P–mediated twitching motility and adhesion to eukaryotic cells, and 4) the vWA module shows an intrinsic binding ability to several host proteins. These findings reveal an elegant yet simple evolutionary tinkering strategy to increase T4P functional versatility by grafting a functional module onto a pilin for presentation by the filaments. This strategy appears to have been extensively used by bacteria, in which modular pilins are widespread and exhibit an astonishing variety of architectures.

Identification of the Active Principle Conferring Anti-Inflammatory and Antinociceptive Properties in Bamboo Plant

Early plants began colonizing earth about 450 million years ago. During the process of coevolution, their metabolic cellular pathways produced a myriad of natural chemicals, many of which remain uncharacterized biologically. Popular preparations containing some of these molecules have been used medicinally for thousands of years. In Brazilian folk medicine, plant extracts from the bamboo plant Guadua paniculata Munro have been used for the treatment of infections and pain. However, the chemical basis of these therapeutic effects has not yet been identified. Here, we performed protein biochemistry and downstream pharmacological assays to determine the mechanisms underlying the anti-inflammatory and antinociceptive effects of an aqueous extract of the G. paniculata rhizome, which we termed AqGP. The anti-inflammatory and antinociceptive effects of AqGP were assessed in mice. We identified and purified a protein (AgGP), with an amino acid sequence similar to that of thaumatins (~20 kDa), capable of repressing inflammation through downregulation of neutrophil recruitment and of decreasing hyperalgesia in mice. In conclusion, we have identified the molecule and the molecular mechanism responsible for the anti-inflammatory and antinociceptive properties of a plant commonly used in Brazilian folk medicine.

Structural insights into histone chaperone Asf1 and its characterization from Plasmodium falciparum

Asf1 is a highly conserved histone chaperone that regulates tightly coupled nucleosome assembly/disassembly process. We observed that Plasmodium falciparum Asf1 (PfAsf1) is ubiquitously expressed in different stages of the life cycle of the parasite. To gain further insight into its biological activity, we solved the structure of N-terminal histone chaperone domain of PfAsf1 (1-159 amino acids) by X-ray crystallography to a resolution of 2.4 Å. The structure is composed of two beta-sheet to form a beta-sandwich, which resembles an immunoglobulin-like fold. The surface-charge distribution of PfAsf1 is distinct from yAsf1 and hAsf1 although the core-structure shows significant similarity. The crystal-structure indicated that PfAsf1 may exist in a dimeric-state which was further confirmed by solution cross-linking experiment. PfAsf1 was found to specifically interact with Plasmodium histone H3 and H4 and was able to deposit H3/H4 dimer onto DNA-template to form disomes, showing its characteristic histone chaperone activity. We mapped the critical residues of PfAsf1 involved in histone H3/H4 interaction and confirmed by site-directed mutagenesis. Further analysis indicates that histone interacting surface of Asf1 is highly conserved while the dimerization interface is variable. Our results identify the role of PfAsf1 as a mediator of chromatin assembly in Plasmodium falciparum, which is the causative agent of malignant malaria in humans.

Characterization of Three Osmotin-Like Proteins from Plumeria rubra and Prospection for Adiponectin Peptidomimetics

BACKGROUND

Osmotin-Like Proteins (OLPs) have been purified and characterized from different plant tissues, including latex fluids. Besides its defensive role, tobacco osmotin seems to induce adiponectin-like physiological effects, acting as an agonist. However, molecular information about this agonistic effect on adiponectin receptors has been poorly exploited and other osmotins have not been investigated yet.

OBJECTIVE AND METHODS

The present study involved the characterization of three OLPs from Plumeria rubra latex and molecular docking studies to evaluate the interaction between them and adiponectin receptors (AdipoR1 and AdipoR2).

RESULTS

P. rubra Osmotin-Like Proteins (PrOLPs) exhibited molecular masses from 21 to 25 kDa and isoelectric points ranging from 4.4 to 7.7. The proteins have 16 cysteine residues, which are involved in eight disulfide bonds, conserved in the same positions as other plant OLPs. The threedimensional (3D) models exhibited the three typical domains of OLPs, and molecular docking analysis showed that two PrOLP peptides interacted with two adiponectin receptors similarly to tobacco osmotin peptide.

CONCLUSION

As observed for tobacco osmotin, the latex osmotins of P. rubra exhibited compatible interactions with adiponectin receptors. Therefore, these plant defense proteins (without known counterparts in humans) are potential tools to study modulation of glucose metabolism in type II diabetes, where adiponectin plays a pivotal role in homeostasis.

Deciphering Molecular Virulence Mechanism of Mycobacterium tuberculosis Dop isopeptidase Based on Its Sequence-Structure-Function Linkage

The pupylation pathway marks proteins for prokaryotic ubiquitin-like protein (Pup)-proteasomal degradation and survival strategy of mycobacteria inside of the host macrophages. Deamidase of Pup (Dop) plays a central role in the pupylation pathway. It is still a matter of investigation to know the function of Dop in virulence of mycobacterial lineage. Hence, the present study was intended to describe the sequence-structure-function-virulence link of Dop for understanding the molecular virulence mechanism of Mycobacterium tuberculosis H37Rv (Mtb). Phylogenetic analysis of this study indicated that Dop has extensively diverged across the proteasome-harboring bacteria. The functional part of Dop was converged across the pathogenic mycobacterial lineage. The genome-wide analysis pointed out that the pupylation gene locus was identical to each other, but its genome neighborhood differed from species to species. Molecular modeling and dynamic studies proved that the predicted structure of Mtb Dop was energetically stable and low conformational freedom. Moreover, evolutionary constraints in Mtb Dop were intensively analyzed for inferring its sequence-structure-function relationships for the full virulence of Mtb. It indicated that evolutionary optimization was extensively required to stabilize its local structural environment at the side chains of mutable residues. The sequence-structure-function-virulence link of Dop might have retained in Mtb by reordering hydrophobic and hydrogen bonding patterns in the local structural environment. Thus, the results of our study provide a quest to understand the molecular virulence and pathogenesis mechanisms of Mtb during the infection process.

Molecular dynamics analysis of conserved water mediated inter-domain recognition of His667-Trp669 in human ceruloplasmin

The human ceruloplasmin (hCP) is the copper containing ferroxidase enzyme with multifunctional activities (NO-oxidase, NO2-synthase,oxidation of neurotransmitters including antioxidants). Therefore, it is of interest to probe the multi-domain hCP using moleculardynamics simulation. Results explain the role played by several conserved water centers in the intra and inter-domain recognition throughH-bond interaction with the interacting residues. We observed seventeen conserved water centers in the inter-domain recognition. Weshow that five invariant water centers W13, W14, W18, W23 and W26 connect the Domain 5 to Domain 4 (D5…W…W4). We also show thatfive other water centers W19, W20, W27, W30 and W31 connects the Domain 5 to Domain 6 (D5…W…W6) that is unique in the simulatedform. The W7 and W32 water centers are involved in the D1…W…W6 recognition. This is important for the water-mediated interaction ofGlu1032 to the trinuclear copper cluster present at the interface between these domains. The involvement of W10 water center in theD3…W10…D4 recognition through Gln552…W10…His667 H-bond interaction is critical in the complexation of CP with myeloperoxidase(Mpo). These observations provide insights to the molecular recognition of hCP with other biomolecules in the system.

MRPC (Missing Regions in Polypeptide Chains): a knowledgebase

Missing regions in protein crystal structures are those regions that cannot be resolved, mainly owing to poor electron density (if the three-dimensional structure was solved using X-ray crystallography). These missing regions are known to have high B factors and could represent loops with a possibility of being part of an active site of the protein molecule. Thus, they are likely to provide valuable information and play a crucial role in the design of inhibitors and drugs and in protein structure analysis. In view of this, an online database, Missing Regions in Polypeptide Chains (MRPC), has been developed which provides information about the missing regions in protein structures available in the Protein Data Bank. In addition, the new database has an option for users to obtain the above data for non-homologous protein structures (25 and 90%). A user-friendly graphical interface with various options has been incorporated, with a provision to view the three-dimensional structure of the protein along with the missing regions using JSmol. The MRPC database is updated regularly (currently once every three months) and can be accessed freely at the URL http://cluster.physics.iisc.ac.in/mrpc.

Putative role of conserved water molecules in the hydration and inter-domain recognition of mono nuclear copper centers in O2-bound human ceruloplasmin: A comparative study between X-ray and MD simulated structures

Human Ceruloplasmin (hCP) is an unique multicopper oxidase which involves in different biological functions e.g., iron metabolism, copper transportation, biogenic amine oxidation ,and its malfunction causes Wilson's and Menkes diseases. MD- simulation studies of O2- bound solvated structure have revealed the role of several conserved/ semi-conserved water molecules in the hydration of type-I copper centers and their involvement to recognition dynamics of these metal centers. In O2- bound structure, hydration potentiality of CuRS (Cu1106) type-I copper center is observed to be unique, where two water molecules (W1-W2) are interacting with the metal sites, which was not found in X-ray structures of hCP. Generally, in the interdomain recognition of CuCys-His to CuRS, CuRS to CuPR and CuPR to CuCys-His centers, the copper bound His-residue of one domain interacts with the Glu-residue of other complementary domain through conserved/ semi-conserved (W3 to W5) water- mediated hydrogen bonds (Cu-His...W...Glu), however direct salt-bridge (Cu-His...Glu) interaction were observed in the X- ray structures. The MD- simulated and X- ray structures have indicated some possibilities on the Cu-His...W...Glu ↔ Cu-His...Glu transition during the interdomain recognition of type-I copper centers, which may have some importance in biology and chemistry of ceruloplasmin.

Designating ligand specificities to metal uptake ABC transporters in Thermus thermophilus HB8

Acquisition of different metal ions by metal uptake ABC transporters of Thermus thermophilus HB8 for accomplishing its various cellular functions.

RepEx: A web server to extract sequence repeats from protein and DNA sequences

Evolution builds up new genetic material from existing ones, not in random, but in highly ordered and eloquent patterns. Most of these sequence repeats are revelatory of valuable information contributing to areas of disease research and function of macromolecules, to name a few. In the age of next generation genome sequencing, rapid and efficient extraction of all unbiased sequence repeats from macromolecules is the need of the hour. In view of this reckoning, an online web-based computing server, RepEx, has been developed to extract and display all possible repeats for DNA and protein sequences. Apart from exact or identical repeats, the server has been designed adeptly to identify and extract degenerate, inverted, everted and mirror repeats from both DNA and protein sequences. The server has striking output displays, featuring interactive graphs and comprehensive output files. In addition, RepEx has been accoutered with an easy-to-use interface and search filters to facilitate a user-defined query or search and is freely available and accessible via the World Wide Web at http://bioserver2.physics.iisc.ac.in/RepEx/.

Elucidation of conformational diversity of druggable enzymes and classification of chemical modulators based on inhibitor-bound structures

Communicated by Ramaswamy H. Sarma.

Insights into product release dynamics through structural analyses of thymidylate kinase

Several studies on enzyme catalysis have pointed out that the product release event could be a rate limiting step. In this study, we have compared the release event of two products, Adenosine di-phosphate (ADP) and Thymidine di-phosphate (TDP) from the active-site of human and Thermus thermophilus thymidine mono-phosphate kinase (TMPK), referred to as hTMPK and ttTMPK, respectively. TMPK catalyses the conversion of Thymidine mono-phosphate (TMP) to TDP using ATP as phosphoryl donor in the presence of Mg²⁺ ion. Most of the earlier studies on this enzyme have focused on understanding substrate binding and catalysis, but the critical product release event remains elusive. Competitive binding experiments of the substrates and the products using ttTMPK apo crystals have indicated that the substrate (TMP) can replace the bound product (TDP), even in the presence of an ADP molecule. Further, the existing random accelerated molecular dynamics (RAMD) simulation program was modified to study the release of both the products simultaneously from the active site. The RAMD simulations on product-bound structures of both ttTMPK and hTMPK, revealed that while several exit patterns of the products are permissible, the sequential exit mode is the most preferred pattern for both ttTMPK and hTMPK enzymes. Additionally, the product release from the hTMPK was found to be faster and more directional as compared to ttTMPK. Structural investigation revealed that the critical changes in the residue composition in the LID-region of ttTMPK and hTMPK have an effect on the product release and can be attributed to the observed differences during product release event. Understanding of these dissimilarities is of considerable utility in designing potent inhibitors or prodrugs that can distinguish between eukaryotic and prokaryotic homologues of thymidylate kinase.

Putative role of invariant water molecules in the X-ray structures of family G fungal endoxylanases

Fungal endo-1,4-beta-xylanases (EC3.2.1.8), because of their widespread industrial applications have become one of the most researched industrial enzymes in recent times. Despite its significance, the role of conserved water molecules in the catalytic activities and structural stability of these enzymes from the fungi have not been studied to a great extent. Our computational structural bioinformatics and MD simulation studies have identified the existence of seven invariant water molecules (IW1- IW7) and reveals the stereo-chemical and electronic consequences of those conserved water molecules in G-xylanase enzyme from eight different fungi. The buried water molecules IW1 and IW2 may have decisive roles in catalysis and may also be associated with ligand binding process of the enzyme, whereas IW3, IW4, IW5, IW6 and IW7 may be involved in stabilizing the important (H144/R145) residues through H-bonds. Possibly they are also involved in the stabilization of secondary structures and anchor to maintain its stereo-chemical architecture. Moreover, a distorted 'W' shaped signature geometry that is observed at the surface of the enzyme can be used to identify the hydrophilic centers in the electron density map of other unknown members of the family G xylanases. The results from this computational investigation could be of interest to a large number of researchers working with the xylanases.

Recognition dynamics of trinuclear copper cluster and associated histidine residues through conserved or semi-conserved water molecules in human Ceruloplasmin: The involvement of aspartic and glutamic acid gates

Human Ceruloplasmin belongs to the family of multi-copper oxidases and it is involved in different physiological processes, copper ion transport, iron metabolism, iron homeostasis, and biogenic amine metabolism. MD-simulation studies have indicated the higher hydrophilic susceptibility of the trinuclear copper cluster in native CP compared to its oxygen bound form. The copper (T2/T3) atom Cu3047 of the cluster, which is close to T1 copper center Cu3052 (~13 Å) has a higher affinity for water molecules compared to other copper centers. The water molecules of W3, W4, W5, W9, and W12 conserved water sites are coordinated to Cu3047, where W3, W9, and W12 centers are found to play some crucial role in the stabilization of native trinuclear copper cluster. The hydrogen bonding interaction of Asp169, Glu112, Asp995, and Glu1032 residues with the copper-bound conserved water molecules (W3, W4, W5, W10, and W11) in native CP is observed to be unique. The conformational flexibility of Asp169 and Glu112 and their association with the copper-bound water molecules, but the absence of such interaction in O₂-bound simulated structure of the enzyme is indicating some plausible rational on the role of these acidic residues in the gating of O₂ molecule in the native trinuclear Cu cluster of CP. The simulation results may shade some new light on the biochemistry/chemistry of CP, specially on the hydration dynamics of the trinuclear copper cluster.

Identification and characterization of two germin-like proteins with oxalate oxidase activity from Calotropis procera latex

Germin-like proteins (GLPs) have been identified in several plant tissues. However, only one work describes GLP in latex fluids. Therefore, the goal of this study was to investigate GLPs in latex and get new insights concerning the structural and functional aspects of these proteins. Two complete sequences with high identity (>50%) with other GLPs, termed CpGLP1 and CpGLP2, were obtained and consecutively presented 216 and 206 amino acid residues, corresponding to molecular masses of 22.7 and 21.7kDa, pI 6.8 and 6.5. The three-dimensional models revealed overall folding similar to those reported for other plant GLPs. Both deduced sequences were grouped into the GER 2 subfamily. Molecular docking studies indicated a putative binding site consisting of three highly conserved histidines and a glutamate residue, which interacted with oxalate. This interaction was later supported by enzymatic assays. Superoxide dismutase (common activity in GLPs) was not detected for CpGLP1 and CpGLP2 by zymogram. The two proteins were detected in the latex, but not in non-germinated or germinated seeds and calli. These results give additional support that germin-like proteins are broadly distributed in plants and they are tissue-specific. This particularity deserves further studies to better understand their functions in latex.

New biochemical insight of conserved water molecules at catalytic and structural Zn2+ ions in human matrix metalloproteinase-I: a study by MD-simulation

Human matrix metalloproteinase (MMP)-1 or collagenase-1 plays a significant role in embryonic development, tissue remodeling, and is also involved in several diseases like arthritis, metastasis, etc. Molecular dynamics simulation studies on hMMP-1 X-ray structures (PDB Id. 1CGE, 1CGF, 1CGL, 1HFC, and 2TCL) suggest that the three conserved water molecules (W_H/1, W_I, W_S) are coordinated with catalytic zinc (Zn_C), and one water molecule (W) is associated at structural zinc ion (Zn_S). Transition of the coordination geometry around Zn_C from tetrahedral to octahedral and tetrahedral to trigonal bipyramidal at Zn_S are also observed during the dynamics. Recognition of two zinc ions through water mediated bridges (Zn_C - W_H (W₁)…W₂….H₁₈₃ - Zn_S) and stabilization of secondary coordination zone around the metal ions indicates the possibility of Zn_C…Zn_S coupled catalytic mechanism in hMMP-I. This study not only reveals a functionally important role of conserved water molecules in hMMP-I but also highlights the involvement of other non catalytic residues, such as S172 and D170 in the catalytic mechanism. The results obtained in this study could be relevant for importance of conserved water mediated recognition site of the sequence residue id. 202(RWTNNFREY)210, interaction of W(tryptophan)203 to zinc bound histidine, their influence on the water molecules that are involved in bridging between Zn_C and Zn_S, and structure-based design of specific hMMP inhibitors. Graphical abstract Water mediated recognition of structural and catalytic zinc ions of hMMP-1 structure (MD simulatated conformation).

Intergenomic evolution and metabolic cross-talk between rumen and thermophilic autotrophic methanogenic archaea

Methanobrevibacter ruminantium M1 (MRU) is a rumen methanogenic archaean that can be able to utilize formate and CO₂/H₂ as growth substrates. Extensive analysis on the evolutionary genomic contexts considered herein to unravel its intergenomic relationship and metabolic adjustment acquired from the genomic content of Methanothermobacter thermautotrophicus ΔH. We demonstrated its intergenomic distance, genome function, synteny homologs and gene families, origin of replication, and methanogenesis to reveal the evolutionary relationships between Methanobrevibacter and Methanothermobacter. Comparison of the phylogenetic and metabolic markers was suggested for its archaeal metabolic core lineage that might have evolved from Methanothermobacter. Orthologous genes involved in its hydrogenotrophic methanogenesis might be acquired from intergenomic ancestry of Methanothermobacter via Methanobacterium formicicum. Formate dehydrogenase (fdhAB) coding gene cluster and carbon monoxide dehydrogenase (cooF) coding gene might have evolved from duplication events within Methanobrevibacter-Methanothermobacter lineage, and fdhCD gene cluster acquired from bacterial origins. Genome-wide metabolic survey found the existence of four novel pathways viz. l-tyrosine catabolism, mevalonate pathway II, acyl-carrier protein metabolism II and glutathione redox reactions II in MRU. Finding of these pathways suggested that MRU has shown a metabolic potential to tolerate molecular oxygen, antimicrobial metabolite biosynthesis and atypical lipid composition in cell wall, which was acquainted by metabolic cross-talk with mammalian bacterial origins. We conclude that coevolution of genomic contents between Methanobrevibacter and Methanothermobacter provides a clue to understand the metabolic adaptation of MRU in the rumen at different environmental niches.

Computational exploration of the binding mode of heme-dependent stimulators into the active catalytic domain of soluble guanylate cyclase

Soluble guanylate cyclase (sGC), the main target of nitric oxide (NO), has been proven to have a significant role in coronary artery disease, pulmonary hypertension, erectile dysfunction, and myocardial infarction. One of its agonists, BAY 41-2272 (Riociguat), has been recently approved for treatment of pulmonary arterial hypertension (PHA), while some others are in clinical phases of development. However, the location of the binding sites for the two known types of agonists, heme-dependent stimulators and heme-independent activators, is a matter of debate, particularly for the first group where both a location on the regulatory (H-NOX) and on the catalytic domain have been suggested by different authors. Here, we address its potential location on the catalytic domain, the unique well characterized at the structural level, by an "in silico" approach. Homology models of the catalytic domain of sGC in "inactive" or "active" conformations were constructed using the structure of previously described crystals of the catalytic domains of "inactive" sGCs (2WZ1, 3ET6) and of "active" adenylate cyclase (1CJU). Each model was submitted to six independent molecular dynamics simulations of about 1 μs. Docking of YC-1, a classic heme-dependent stimulator, to all frames of representative trajectories of "inactive" and "active" conformations, followed by calculation of absolute binding free energies with the linear interaction energy (LIE) method, revealed a potential high-affinity binding site on the "active" structure. The site, located between the pseudo-symmetric and the catalytic site just over the loop β2 -β3 , does not overlap with the forskolin binding site on adenylate cyclases. Proteins 2016; 84:1534-1548. © 2016 Wiley Periodicals, Inc.

A polyalanine peptide derived from polar fish with anti-infectious activities

Due to the growing concern about antibiotic-resistant microbial infections, increasing support has been given to new drug discovery programs. A promising alternative to counter bacterial infections includes the antimicrobial peptides (AMPs), which have emerged as model molecules for rational design strategies. Here we focused on the study of Pa-MAP 1.9, a rationally designed AMP derived from the polar fish Pleuronectes americanus. Pa-MAP 1.9 was active against Gram-negative planktonic bacteria and biofilms, without being cytotoxic to mammalian cells. By using AFM, leakage assays, CD spectroscopy and in silico tools, we found that Pa-MAP 1.9 may be acting both on intracellular targets and on the bacterial surface, also being more efficient at interacting with anionic LUVs mimicking Gram-negative bacterial surface, where this peptide adopts α-helical conformations, than cholesterol-enriched LUVs mimicking mammalian cells. Thus, as bacteria present varied physiological features that favor antibiotic-resistance, Pa-MAP 1.9 could be a promising candidate in the development of tools against infections caused by pathogenic bacteria.

New insight into the architecture of oxy-anion pocket in unliganded conformation of GAT domains: A MD-simulation study

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N-terminal "glutaminase" (GAT) and C-terminal "synthetase" domain. The enzyme is identified as a potential target for anti-cancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H-bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)-substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non-catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site.

The putative role of some conserved water molecules in the structure and function of human transthyretin

Human transthyretin (hTTR) is a multifunctional protein that is involved in several neurodegenerative diseases. Besides the transportation of thyroxin and vitamin A, it is also involved in the proteolysis of apolipoprotein A1 and Aβ peptide. Extensive analyses of 32 high-resolution X-ray and neutron diffraction structures of hTTR followed by molecular-dynamics simulation studies using a set of 15 selected structures affirmed the presence of 44 conserved water molecules in its dimeric structure. They are found to play several important roles in the structure and function of the protein. Eight water molecules stabilize the dimeric structure through an extensive hydrogen-bonding network. The absence of some of these water molecules in highly acidic conditions (pH ≤ 4.0) severely affects the interfacial hydrogen-bond network, which may destabilize the native tetrameric structure, leading to its dissociation. Three pairs of conserved water molecules contribute to maintaining the geometry of the ligand-binding cavities. Some other water molecules control the orientation and dynamics of different structural elements of hTTR. This systematic study of the location, absence, networking and interactions of the conserved water molecules may shed some light on various structural and functional aspects of the protein. The present study may also provide some rational clues about the conserved water-mediated architecture and stability of hTTR.