Identification of a novel PHGDH covalent inhibitor by chemical proteomics and phenotypic profiling

The first rate-limiting enzyme of the serine synthesis pathway (SSP), phosphoglycerate dehydrogenase (PHGDH), is hyperactive in multiple tumors, which leads to the activation of SSP and promotes tumorigenesis. However, only a few inhibitors of PHGDH have been discovered to date, especially the covalent inhibitors of PHGDH. Here, we identified withangulatin A (WA), a natural small molecule, as a novel covalent inhibitor of PHGDH. Affinity-based protein profiling identified that WA could directly bind to PHGDH and inactivate the enzyme activity of PHGDH. Biolayer interferometry and LC-MS/MS analysis further demonstrated the selective covalent binding of WA to the cysteine 295 residue (Cys295) of PHGDH. With the covalent modification of Cys295, WA blocked the substrate-binding domain (SBD) of PHGDH and exerted an allosteric effect to induce PHGDH inactivation. Further studies revealed that with the inhibition of PHGDH mediated by WA, the glutathione synthesis was decreased and intracellular levels of reactive oxygen species (ROS) were elevated, leading to the inhibition of tumor proliferation. This study indicates WA as a novel PHGDH covalent inhibitor, which identifies Cys295 as a novel allosteric regulatory site of PHGDH and holds great potential in developing anti-tumor agents for targeting PHGDH.

Discriminating between competing models for the allosteric regulation of oncogenic phosphatase SHP2 by characterizing its active state

The Src-homology 2 domain containing phosphatase 2 (SHP2) plays a critical role in crucial signaling pathways and is involved in oncogenesis and in developmental disorders. Its structure includes two SH2 domains (N-SH2 and C-SH2), and a protein tyrosine phosphatase (PTP) domain. Under basal conditions, SHP2 is auto-inhibited, with the N-SH2 domain blocking the PTP active site. Activation involves a rearrangement of the domains that makes the catalytic site accessible, coupled to the association between the SH2 domains and cognate proteins containing phosphotyrosines. Several aspects of this transition are debated and competing mechanistic models have been proposed. A crystallographic structure of SHP2 in an active state has been reported (PDB code 6crf), but several lines of evidence suggests that it is not fully representative of the conformations populated in solution. To clarify the structural rearrangements involved in SHP2 activation, enhanced sampling simulations of the autoinhibited and active states have been performed, for wild type SHP2 and its pathogenic E76K variant. Our results demonstrate that the crystallographic conformation of the active state is unstable in solution, and multiple interdomain arrangements are populated, thus allowing association to bisphosphorylated sequences. Contrary to a recent proposal, activation is coupled to the conformational changes of the N-SH2 binding site, which is significantly more accessible in the active sate, rather than to the structure of the central β-sheet of the domain. In this coupling, a previously undescribed role for the N-SH2 BG loop emerged.

Identification of SARS-CoV-2 RNA-dependent RNA polymerase inhibitors from the major phytochemicals of Nigella sativa: An in silico approach

The coronavirus disease 2019 (COVID-19), which emerged in December 2019, continues to be a serious health concern worldwide. There is an urgent need to develop effective drugs and vaccines to control the spread of this disease. In the current study, the main phytochemical compounds of Nigella sativa were screened for their binding affinity for the active site of the RNA-dependent RNA polymerase (RdRp) enzyme of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The binding affinity was investigated using molecular docking methods, and the interaction of phytochemicals with the RdRp active site was analyzed and visualized using suitable software. Out of the nine phytochemicals of N. sativa screened in this study, a significant docking score was observed for four compounds, namely α-hederin, dithymoquinone, nigellicine, and nigellidine. Based on the findings of our study, we report that α-hederin, which was found to possess the lowest binding energy (–8.6 kcal/mol) and hence the best binding affinity, is the best inhibitor of RdRp of SARS-CoV-2, among all the compounds screened here. Our results prove that the top four potential phytochemical molecules of N. sativa, especially α-hederin, could be considered for ongoing drug development strategies against SARS-CoV-2. However, further in vitro and in vivo testing are required to confirm the findings of this study.

Neuroprotective derivatives of tacrine that target NMDA receptor and acetyl cholinesterase - Design, synthesis and biological evaluation

The complex and multifactorial nature of neuropsychiatric diseases demands multi-target drugs that can intervene with various sub-pathologies underlying disease progression. Targeting the impairments in cholinergic and glutamatergic neurotransmissions with small molecules has been suggested as one of the potential disease-modifying approaches for Alzheimer’s disease (AD). Tacrine, a potent inhibitor of acetylcholinesterase (AChE) is the first FDA approved drug for the treatment of AD. Tacrine is also a low affinity antagonist of N-methyl-D-aspartate receptor (NMDAR). However, tacrine was withdrawn from its clinical use later due to its hepatotoxicity. With an aim to develop novel high affinity multi-target directed ligands (MTDLs) against AChE and NMDAR, with reduced hepatotoxicity, we performed in silico structure-based modifications on tacrine, chemical synthesis of the derivatives and in vitro validation of their activities. Nineteen such derivatives showed inhibition with IC₅₀ values in the range of 18.53 ± 2.09 – 184.09 ± 19.23 nM against AChE and 0.27 ± 0.05 – 38.84 ± 9.64 μM against NMDAR. Some of the selected compounds also protected rat primary cortical neurons from glutamate induced excitotoxicity. Two of the tacrine derived MTDLs, 201 and 208 exhibited in vivo efficacy in rats by protecting against behavioral impairment induced by administration of the excitotoxic agent, monosodium glutamate. Additionally, several of these synthesized compounds also exhibited promising inhibitory activitiy against butyrylcholinesterase. MTDL-201 was also devoid of hepatotoxicity in vivo. Given the therapeutic potential of MTDLs in disease-modifying therapy, our studies revealed several promising MTDLs among which 201 appears to be a potential candidate for immediate preclinical evaluations.

AncPhore: A versatile tool for anchor pharmacophore steered drug discovery with applications in discovery of new inhibitors targeting metallo-β-lactamases and indoleamine/tryptophan 2,3-dioxygenases

We herein describe AncPhore, a versatile tool for drug discovery, which is characterized by pharmacophore feature analysis and anchor pharmacophore (i.e., most important pharmacophore features) steered molecular fitting and virtual screening. Comparative analyses of numerous protein–ligand complexes using AncPhore revealed that anchor pharmacophore features are biologically important, commonly associated with protein conservative characteristics, and have significant contributions to the binding affinity. Performance evaluation of AncPhore showed that it had substantially improved prediction ability on different types of target proteins including metalloenzymes by considering the specific contributions and diversity of anchor pharmacophore features. To demonstrate the practicability of AncPhore, we screened commercially available chemical compounds and discovered a set of structurally diverse inhibitors for clinically relevant metallo-β-lactamases (MBLs); of them, 4 and 6 manifested potent inhibitory activity to VIM-2, NDM-1 and IMP-1 MBLs. Crystallographic analyses of VIM-2:4 complex revealed the precise inhibition mode of 4 with VIM-2, highly consistent with the defined anchor pharmacophore features. Besides, we also identified new hit compounds by using AncPhore for indoleamine/tryptophan 2,3-dioxygenases (IDO/TDO), another class of clinically relevant metalloenzymes. This work reveals anchor pharmacophore as a valuable concept for target-centered drug discovery and illustrates the potential of AncPhore to efficiently identify new inhibitors for different types of protein targets.

Rewards of divergence in sequences, 3-D structures and dynamics of yeast and human spliceosome SF3b complexes

The evolution of homologous and functionally equivalent multiprotein assemblies is intriguing considering sequence divergence of constituent proteins. Here, we studied the implications of protein sequence divergence on the structure, dynamics and function of homologous yeast and human SF3b spliceosomal subcomplexes. Human and yeast SF3b comprise of 7 and 6 proteins respectively, with all yeast proteins homologous to their human counterparts at moderate sequence identity. SF3b6, an additional component in the human SF3b, interacts with the N-terminal extension of SF3b1 while the yeast homologue Hsh155 lacks the equivalent region. Through detailed homology studies, we show that SF3b6 is absent not only in yeast but in multiple lineages of eukaryotes implying that it is critical in specific organisms. We probed for the potential role of SF3b6 in the spliceosome assembled form through structural and flexibility analyses. By analysing normal modes derived from anisotropic network models of SF3b1, we demonstrate that when SF3b1 is bound to SF3b6, similarities in the magnitude of residue motions (0.86) and inter-residue correlated motions (0.94) with Hsh155 are significantly higher than when SF3b1 is considered in isolation (0.21 and 0.89 respectively). We observed that SF3b6 promotes functionally relevant 'open-to-close' transition in SF3b1 by enhancing concerted residue motions. Such motions are found to occur in the Hsh155 without SF3b6. The presence of SF3b6 influences motions of 16 residues that interact with U2 snRNA/branchpoint duplex and supports the participation of its interface residues in long-range communication in the SF3b1. These results advocate that SF3b6 potentially acts as an allosteric regulator of SF3b1 for BPS selection and might play a role in alternative splicing. Furthermore, we observe variability in the relative orientation of SF3b4 and in the local structure of three β-propeller domains of SF3b3 with reference to their yeast counterparts. Such differences influence the inter-protein interactions of SF3b between these two organisms. Together, our findings highlight features of SF3b evolution and suggests that the human SF3b may have evolved sophisticated mechanisms to fine tune its molecular function.

Identification of alkaloids from Justicia adhatoda as potent SARS CoV-2 main protease inhibitors: An in silico perspective

The COVID-19 pandemic, caused by SARS CoV-2, is responsible for millions of death worldwide. No approved/proper therapeutics is currently available which can effectively combat this outbreak. Several attempts have been undertaken in the search of effective drugs to control the spread of SARS CoV-2 infection. The main protease (Mpro), key component for the cleavage of the viral polyprotein, is considered to be one of the important drug targets for treating COVID-19. Various phytochemicals, including polyphenols and alkaloids, have been proposed as potent inhibitors of Mpro. The alkaloids from leaf extracts of Justicia adhatoda have also been reported to possess anti-viral activity. But whether these alkaloids exhibit any inhibitory effect on SARS CoV-2 Mpro is far from clear. To explore this in detail, we have adopted computational approaches. Justicia adhatoda alkaloids possessing proper drug-likeness properties and two anti-HIV drugs (lopinavir and darunavir; having binding affinity -7.3 to -7.4 kcal/mol) were docked against SARS CoV-2 Mpro to study their binding properties. Only one alkaloid (anisotine) had interaction with both the catalytic residues (His41 and Cys145) of Mpro and exhibited good binding affinity (-7.9 kcal/mol). Molecular dynamic simulations (100 ns) revealed that Mpro-anisotine complex is more stable, conformationally less fluctuated; slightly less compact and marginally expanded than Mpro-darunavir/lopinavir complex. Even the number of intermolecular H-bonds and MM-GBSA analysis suggested that anisotine is a more potent Mpro inhibitor than the two previously recommended antiviral drugs (lopinavir and darunavir) and may evolve as a promising anti-COVID-19 drug if proven in animal experiments and on patients.

Bispecific antibodies with Fab-arms featuring exchanged antigen-binding constant domains

Monoclonal antibodies can acquire the property of engagement of a second antigen via fusion methods or modification of their CDR loops, but also by modification of their constant domains, such as in the mAb² format where a set of mutated amino acid residues in the C_H3 domains enables a high-affinity specific interaction with the second antigen. We tested the possibility of introducing multiple binding sites for the second antigen by replacing the Fab C_H1/C_L domain pair with a pair of antigen-binding C_H3 domains in a model scaffold with trastuzumab variable domains and VEGF-binding C_H3 domains. Such bispecific molecules were produced in a “Fab-like” format and in a full-length antibody format. Novel constructs were of expected molecular composition using mass spectrometry. They were expressed at a high level in standard laboratory conditions, purified as monomers with Protein A and gel filtration and were of high thermostability. Their high-affinity binding to both target antigens was retained. Finally, the Her2/VEGF binding domain-exchanged bispecific antibody was able to mediate a potentiated surface Her2-internalization effect on the Her2-overexpressing cell line SK-BR-3 due to improved level of cross-linking with the endogenously secreted cytokine. To conclude, bispecific antibodies with Fabs featuring exchanged antigen-binding C_H3 domains offer an alternative solution in positioning and valency of antigen binding sites.

•

Fab constant domains can be efficiently exchanged for antigen-binding C_H3 domains.

•

Such mutagenesis results in bispecific antibodies with correct chain pairing.

•

Domain-exchanged bispecific Fab- and IgG-like formats are of favorable biophysical properties.

•

Resulting bispecific antibodies show high-affinity binding to both target antigens.

In silico structure modelling of SARS-CoV-2 Nsp13 helicase and Nsp14 and repurposing of FDA approved antiviral drugs as dual inhibitors

The high mortality rate from the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in humans and the lack of effective therapeutic regime for its treatment necessitates the identification of new antivirals. SARS-CoV-2 relies on non-structural proteins such as Nsp13 helicase and nsp14 which are the key components of the replication-transcription complex (RTC) to complete its infectious life cycle. Therefore, targeting these essential viral proteins with small molecules will most likely to halt the disease pathogenesis. The lack of experimental structures of these proteins deters the process of structure-based identification of their specific inhibitors. In the present study, the in silico models of SARS-CoV-2 nsp13 helicase and nsp14 protein were elucidated using a comparative homology modelling approach. These in silico model structures were validated using various parameters such as Ramachandran plot, Verify 3D score, ERRAT score, knowledge-based energy and Z-score. The in silico models were further used for virtual screening of the Food and Drug Administration (FDA) approved antiviral drugs. Simeprevir (SMV), Paritaprevir (PTV) and Grazoprevir (GZR) were the common leads identified which show higher binding affinity to both nsp13 helicase and nsp14 as compared to the control inhibitors and therefore, they might be potential dual-target inhibitors. The leads also establish a network of hydrogen bonds and hydrophobic interactions with the key residues lining the active site pockets. The present findings suggest that these FDA approved antiviral drugs can be subjected to repurposing against SARS-CoV-2 infection after verifying the in silico results through in vitro and in vivo studies.

Crystallographic and cryogenic electron microscopic structures and enzymatic characterization of sulfur oxygenase reductase from Sulfurisphaera tokodaii

•

Sulfur oxygenase reductase (SOR) was biochemically and structurally characterized.

•

High resolution structures of SOR were determined by crystallography and cryo-EM.

•

Twenty-four identical subunits of SOR form a hollow sphere.

•

Catalytic components exhibited different features in the crystal and cryo-EM structures.

Sulfur oxygenase reductases (SORs) are present in thermophilic and mesophilic archaea and bacteria, and catalyze oxygen-dependent oxygenation and disproportionation of elemental sulfur. SOR has a hollow, spherical homo-24-mer structure and reactions take place at active sites inside the chamber. The crystal structures of SORs from Acidianus species have been reported. However, the states of the active site components (mononuclear iron and cysteines) and the entry and exit paths of the substrate and products are still in dispute. Here, we report the biochemical and structural characterizations of SORs from the thermoacidophilic archaeon Sulfurisphaera tokodaii (StSOR) and present high-resolution structures determined by X-ray crystallography and cryogenic electron microscopy (cryo-EM). The crystal structure of StSOR was determined at 1.73 Å resolution. At the catalytic center, iron is ligated to His86, His90, Glu114, and two water molecules. Three conserved cysteines in the cavity are located 9.5–13 Å from the iron and were observed as free thiol forms. A mutational analysis indicated that the iron and one of the cysteines (Cys31) were essential for both activities. The cryo-EM structure was determined at 2.24 Å resolution using an instrument operating at 200 kV. The two structures determined by different methodologies showed similar main chain traces, but the maps exhibited different features at catalytically important components. A possible role of StSOR in the sulfur metabolism of S. tokodaii (an obligate aerobe) is discussed based on this study. Given the high resolution achieved in this study, StSOR was shown to be a good benchmark sample for cryo-EM.

Structural modeling, mutation analysis, and in vitro expression of usherin, a major protein in inherited retinal degeneration and hearing loss

Usherin is the most common causative protein associated with autosomal recessive retinitis pigmentosa (RP) and Usher syndrome (USH), which are characterized by retinal degeneration alone and in combination with hearing loss, respectively. Usherin is essential for photoreceptor survival and hair cell bundle integrity. However, the molecular mechanism underlying usherin function in normal and disease conditions is unclear. In this study, we investigated structural models of usherin domains and localization of usherin pathogenic small in-frame mutations, mainly homozygous missense mutations. We found that usherin fibronectin III (FN3) domains and most laminin-related domains have a β-sandwich structure. Some FN3 domains are predicted to interact with each other and with laminin-related domains. The usherin protein may bend at some FN3 linker regions. RP- and USH-associated small in-frame mutations are differentially located in usherin domains. Most of them are located at the periphery of β-sandwiches, with some at the interface between interacting domains. The usherin laminin epidermal growth factor repeats adopt a rod-shaped structure, which is maintained by disulfide bonds. Most missense mutations and deletion of exon 13 in this region disrupt the disulfide bonds and may affect local protein folding. Despite low expression of the recombinant entire protein and protein fragments in mammalian cell culture, usherin FN3 fragments are more robustly expressed and secreted than its laminin-related fragments. Our findings provide new insights into the usherin structure and the disease mechanisms caused by pathogenic small in-frame mutations, which will help inform future experimental research on diagnosis, disease mechanisms, and therapeutic approaches.

A hallmark of phospholamban functional divergence is located in the N-terminal phosphorylation domain

Sarcoplasmic reticulum Ca²⁺ pump (SERCA) is a critical component of the Ca²⁺ transport machinery in myocytes. There is clear evidence for regulation of SERCA activity by PLB, whose activity is modulated by phosphorylation of its N-terminal domain (residues 1–25), but there is less clear evidence for the role of this domain in PLB’s functional divergence. It is widely accepted that only sarcolipin (SLN), a protein that shares substantial homology with PLB, uncouples SERCA Ca²⁺ transport from ATP hydrolysis by inducing a structural change of its energy-transduction domain; yet, experimental evidence shows that the transmembrane domain of PLB (residues 26–52, PLB_26–52) partially uncouples SERCA in vitro. These apparently conflicting mechanisms suggest that PLB’s uncoupling activity is encoded in its transmembrane domain, and that it is controlled by the N-terminal phosphorylation domain. To test this hypothesis, we performed molecular dynamics simulations (MDS) of the binary complex between PLB_26–52 and SERCA. Comparison between PLB_26–52 and wild-type PLB (PLB_WT) showed no significant changes in the stability and orientation of the transmembrane helix, indicating that PLB_26–52 forms a native-like complex with SERCA. MDS showed that PLB_26–52 produces key intermolecular contacts and structural changes required for inhibition, in agreement with studies showing that PLB_26–52 inhibits SERCA. However, deletion of the N-terminal phosphorylation domain facilitates an order-to-disorder shift in the energy-transduction domain associated with uncoupling of SERCA, albeit weaker than that induced by SLN. This mechanistic evidence reveals that the N-terminal phosphorylation domain of PLB is a primary contributor to the functional divergence among homologous SERCA regulators.

Role of Structural Features in Oligomerization, Active-Site Integrity and Ligand Binding of Ribose-1,5-Bisphosphate Isomerase

Pentose bisphosphate pathway, exclusively found in archaea, is similar to the pentose phosphate pathway present in bacteria and eukarya. In pentose bisphosphate pathway, the conversion of ribose moieties of nucleosides into 3-phosphoglycerate (3-PGA) involves multiple steps; one of them being the conversion of ribose-1,5-bisphosphate (R15P) to ribulose-1,5-bisphosphate (RuBP) catalyzed by an enzyme ribose-1,5-bisphosphate isomerase (R15Pi). The availability of the three-dimensional structure of R15Pi had facilitated the understanding of various structural and functional aspects of the enzyme. Nevertheless, the structure of R15Pi also left several significant questions unanswered that would aid in understanding the structure-function relationship of the enzyme. Thus, we have taken up a computational approach to further understand the role of various structural features of the enzyme R15Pi. Results obtained from molecular dynamics (MD) simulations aided in understanding the obligation of the enzyme R15Pi to oligomerize and also in deciphering the role of catalytic residue(s) in structural stability. Identification of invariant water molecules of the enzyme R15Pi helped in discerning their significance at the active-site pocket and structurally important regions. Further, molecular docking studies allowed the identification of the amino acid residues essential for holding the substrate (R15P) or product (RuBP) in the vicinity of the active site of the enzyme R15Pi. Interestingly, results of the molecular docking studies assisted in the identification of an “alternate binding site” for the substrate, R15P. Finally, based on these results, we propose a mechanism of “substrate sliding” for the enzyme R15Pi.

•

Hexameric state is obligatory for the functionality of R15Pi.

•

Cysteine is more important catalytic residue than aspartate.

•

Invariant water molecules are conserved at the active site and structurally important areas.

•

The substrate does not directly bind to the active site.

•

The substrate slides from the “initial binding site” to the active site.

Allysine modifications perturb tropoelastin structure and mobility on a local and global scale

Elastin provides elastic tissues with resilience through stretch and recoil cycles, and is primarily made of its extensively cross-linked monomer, tropoelastin. Here, we leverage the recently published full atomistic model of tropoelastin to assess how allysine modifications, which are essential to cross-linking, contribute to the dynamics and structural changes that occur in tropoelastin in the context of elastin assembly. We used replica exchange molecular dynamics to generate structural ensembles of allysine containing tropoelastin. We conducted principal component analysis on these ensembles and found that the molecule departs from the canonical structural ensemble. Furthermore, we showed that, while the canonical scissors-twist movement was retained, new movements emerged that deviated from those of the wild type protein, providing evidence for the involvement of a variety of molecular motions in elastin assembly. Additionally, we highlighted secondary structural changes and linked these perturbations to the longevity of specific salt bridges. We propose a model where allysines in tropoelastin contribute to hierarchical elastin assembly through global and local perturbations to molecular structure and dynamics.

•

converting lysine to allysine by lysyl oxidases is needed to generate crosslinks between tropoelastin molecules in order to make elastin

•

structural changes in the intact tropoelastin molecule ensue where modified tropoelastin molecules structurally depart from the canonical ensemble

•

new molecular motions deviate from those of unmodified tropoelastin

•

persistence times of specific salt bridges contribute to these perturbations

•

allysines in tropoelastin contribute to hierarchical elastin assembly through global and local perturbations to molecular structure and dynamics

Molecular dynamics simulation and binding free energy studies of novel leads belonging to the benzofuran class inhibitors of Mycobacterium tuberculosis Polyketide Synthase 13

In this work, the binding mechanism of new Polyketide Synthase 13 (Pks13) inhibitors has been studied through molecular dynamics simulation and free energy calculations. The drug Tam1 and its analogs, belonging to the benzofuran class, were submitted to 100 ns simulations, and according to the results obtained for root mean square deviation, all the simulations converged from approximately 30 ns. For the analysis of backbone flotation, the root mean square fluctuations were plotted for the Cα atoms; analysis revealed that the greatest fluctuation occurred in the residues that are part of the protein lid domain. The binding free energy value (ΔG_bind) obtained for the Tam16 lead molecule was of -51.43 kcal/mol. When comparing this result with the ΔG_bind values for the remaining analogs, the drug Tam16 was found to be the highest ranked: this result is in agreement with the experimental results obtained by Aggarwal and collaborators, where it was verified that the IC₅₀ for Tam16 is the smallest necessary to inhibit the Pks13 (IC₅₀ = 0.19 μM). The energy decomposition analysis suggested that the residues which most interact with inhibitors are: Ser1636, Tyr1637, Asn1640, Ala1667, Phe1670, and Tyr1674, from which the greatest energy contribution to Phe1670 was particularly notable. For the lead molecule Tam16, a hydrogen bond with the hydroxyl of the phenol not observed in the other analogs induced a more stable molecular structure. Aggarwal and colleagues reported this hydrogen bonding as being responsible for the stability of the molecule, optimizing its physic-chemical, toxicological, and pharmacokinetic properties.

A novel identification approach for discovery of 5-HydroxyTriptamine 2A antagonists: combination of 2D/3D similarity screening, molecular docking and molecular dynamics

5-HydroxyTriptamine 2A antagonists are potential targets for treatment of various cerebrovascular and cardiovascular disorders. In this study, we have developed and performed a unique screening pipeline for filtering ZINC database compounds on the basis of similarities to known antagonists to determine novel small molecule antagonists of 5-HydroxyTriptamine 2A. The screening pipeline is based on 2D similarity, 3D dissimilarity and a combination of 2D/3D similarity. The shortlisted compounds were docked to a 5-HydroxyTriptamine 2A homology-based model, and complexes with low binding energies (287 complexes) were selected for molecular dynamics (MD) simulations in a lipid bilayer. The MD simulations of the shortlisted compounds in complex with 5-HydroxyTriptamine 2A confirmed the stability of the complexes and revealed novel interaction insights. The receptor residues S239, N343, S242, S159, Y370 and D155 predominantly participate in hydrogen bonding. π-π stacking is observed in F339, F340, F234, W151 and W336, whereas hydrophobic interactions are observed amongst V156, F339, F234, V362, V366, F340, V235, I152 and W151. The known and potential antagonists shortlisted by us have similar overlapping molecular interaction patterns. The 287 potential 5-HydroxyTriptamine 2A antagonists may be experimentally verified.

Structural distortions due to missense mutations in human formylglycine-generating enzyme leading to multiple sulfatase deficiency

The major candidate for multiple sulfatase deficiency is a defective formylglycine-generating enzyme (FGE). Though adequately produced, mutations in FGE stall the activation of sulfatases and prevent their activity. Missense mutations, viz. E130D, S155P, A177P, W179S, C218Y, R224W, N259I, P266L, A279V, C336R, R345C, A348P, R349Q and R349W associated with multiple sulfatase deficiency are yet to be computationally studied. Aforementioned mutants were initially screened through ws-SNPs&GO^3D program. Mutant R345C acquired the highest score, and hence was studied in detail. Discrete molecular dynamics explored structural distortions due to amino acid substitution. Therein, comparative analyses of wild type and mutant were carried out. Changes in structural contours were observed between wild type and mutant. Mutant had low conformational fluctuation, high atomic mobility and more compactness than wild type. Moreover, free energy landscape showed mutant to vary in terms of its conformational space as compared to wild type. Subsequently, wild type and mutant were subjected to single-model analyses. Mutant had lesser intra molecular interactions than wild type suggesting variations pertaining to its secondary structure. Furthermore, simulated thermal denaturation showed dissimilar pattern of hydrogen bond dilution. Effects of these variations were observed as changes in elements of secondary structure. Docking studies of mutant revealed less favourable binding energy towards its substrate as compared to wild type. Therefore, theoretical explanations for structural distortions of mutant R345C leading to multiple sulfatase deficiency were revealed. The protocol of the study could be useful to examine the effectiveness of pharmacological chaperones prior to experimental studies.

MoFvAb: Modeling the Fv region of antibodies

Knowledge of the 3-dimensional structure of the antigen-binding region of antibodies enables numerous useful applications regarding the design and development of antibody-based drugs. We present a knowledge-based antibody structure prediction methodology that incorporates concepts that have arisen from an applied antibody engineering environment. The protocol exploits the rich and continuously growing supply of experimentally derived antibody structures available to predict CDR loop conformations and the packing of heavy and light chain quickly and without user intervention. The homology models are refined by a novel antibody-specific approach to adapt and rearrange sidechains based on their chemical environment. The method achieves very competitive all-atom root mean square deviation values in the order of 1.5 Å on different evaluation datasets consisting of both known and previously unpublished antibody crystal structures.

Structural insight into activation of homoserine dehydrogenase from the archaeon Sulfolobus tokodaii via reduction

Homoserine dehydrogenase (HSD; 305 amino acid residues) catalyzes an NAD(P)-dependent reversible reaction between l-homoserine and aspartate 4-semialdehyde and is involved in the aspartate pathway. HSD from the hyperthermophilic archaeon Sulfolobus tokodaii was markedly activated (2.5-fold) by the addition of 0.8 mM dithiothreitol. The crystal structure of the homodimer indicated that the activation was caused by cleavage of the disulfide bond formed between two cysteine residues (C303) in the C-terminal regions of the two subunits.

Structural alterations by five disease-causing mutations in the low-pH conformation of human dihydrolipoamide dehydrogenase (hLADH) analyzed by molecular dynamics - Implications in functional loss and modulation of reactive oxygen species generation by pathogenic hLADH forms

Human dihydrolipoamide dehydrogenase (hLADH) is a flavoenzyme component (E3) of the human alpha-ketoglutarate dehydrogenase complex (α-KGDHc) and few other dehydrogenase complexes. Pathogenic mutations of hLADH cause severe metabolic diseases (atypical forms of E3 deficiency) that often escalate to cardiological or neurological presentations and even premature death; the pathologies are generally accompanied by lactic acidosis. hLADH presents a distinct conformation under acidosis (pH 5.5–6.8) with lower physiological activity and the capacity of generating reactive oxygen species (ROS). It has been shown by our laboratory that selected pathogenic mutations, besides lowering the physiological activity of hLADH, significantly stimulate ROS generation by hLADH, especially at lower pH, which might play a role in the pathogenesis of E3-deficiency in respective cases. Previously, we generated by molecular dynamics (MD) simulation the low-pH hLADH structure and analyzed the structural changes induced in this structure by eight of the pathogenic mutations of hLADH. In the absence of high resolution mutant structures these pieces of information are crucial for the mechanistic investigation of the molecular pathogeneses of the hLADH protein. In the present work we analyzed by molecular dynamics simulation the structural changes induced in the low-pH conformation of hLADH by five pathogenic mutations of hLADH; the structures of these disease-causing mutants of hLADH have never been examined before.

•

5 disease-causing mutants of hLADH were subjected to MD to reveal structural changes.

•

MD simulations were carried out both in vacuum and in water supplemented with ions.

•

Functional regions significantly affected by mutation were identified.

•

Implicated residues are to be targeted in mechanistic studies of hLADH dysfunction.

Analysis of consequences of non-synonymous SNP in feed conversion ratio associated TGF-β receptor type 3 gene in chicken

The recent advances in high throughput sequencing technology accelerate possible ways for the study of genome wide variation in several organisms and associated consequences. In the present study, mutations in TGFBR3 showing significant association with FCR trait in chicken during exome sequencing were further analyzed. Out of four SNPs, one nsSNP p.Val451Leu was found in the coding region of TGFBR3. In silico tools such as SnpSift and PANTHER predicted it as deleterious (0.04) and to be tolerated, respectively, while I-Mutant revealed that protein stability decreased. The TGFBR3 I-TASSER model has a C-score of 0.85, which was validated using PROCHECK. Based on MD simulation, mutant protein structure deviated from native with RMSD 0.08 Å due to change in the H-bonding distances of mutant residue. The docking of TGFBR3 with interacting TGFBR2 inferred that mutant required more global energy. Therefore, the present study will provide useful information about functional SNPs that have an impact on FCR traits.

•

Investigated functional nsSNP p.Val451Leu (rs312979494) in feed conversion ratio (FCR) associated TGFBR3 of chicken

•

Computational tools (SIFT and I-Mutant 2.0) predicted that this nsSNP was deleterious.

•

Mutant structure of TGFBR3 showed high energies and RMS deviations compared to native using MD simulation.

•

Molecular docking of TGFBR3 with interacting protein TGFBR2 showed an increase in global energy of mutant compared to native.

•

We have predicted that functional SNP has an impact on TGFBR3 of chicken and thus can be treated as candidate SNP in screening.

MutT from the fish pathogen Aliivibrio salmonicida is a cold-active nucleotide-pool sanitization enzyme with unexpectedly high thermostability

Upon infection by pathogenic bacteria, production of reactive oxygen species (ROS) is part of the host organism's first line of defence. ROS damage a number of macromolecules, and in order to withstand such a harsh environment, the bacteria need to have well-functioning ROS scavenging and repair systems. Herein, MutT is an important nucleotide-pool sanitization enzyme, which degrades 8-oxo-dGTP and thus prevents it from being incorporated into DNA. In this context, we have performed a comparative biochemical and structural analysis of MutT from the fish pathogen Aliivibrio salmonicida (AsMutT) and the human pathogen Vibrio cholerae (VcMutT), in order to analyse their function as nucleotide sanitization enzymes and also determine possible cold-adapted properties of AsMutT. The biochemical characterisation revealed that both enzymes possess activity towards the 8-oxo-dGTP substrate, and that AsMutT has a higher catalytic efficiency than VcMutT at all temperatures studied. Calculations based on the biochemical data also revealed a lower activation energy (E a) for AsMutT compared to VcMutT, and differential scanning calorimetry experiments showed that AsMutT displayed an unexpected higher melting temperature (T m) value than VcMutT. A comparative analysis of the crystal structure of VcMutT, determined to 2.42 Å resolution, and homology models of AsMutT indicate that three unique Gly residues in loops of VcMutT, and additional long range ion-pairs in AsMutT could explain the difference in temperature stability of the two enzymes. We conclude that AsMutT is a stable, cold-active enzyme with high catalytic efficiency and reduced E a, compared to the mesophilic VcMutT.

Contribution of histone N-terminal tails to the structure and stability of nucleosomes

Histones are the protein components of the nucleosome, which forms the basic architecture of eukaryotic chromatin. Histones H2A, H2B, H3, and H4 are composed of two common regions, the "histone fold" and the "histone tail". Many efforts have been focused on the mechanisms by which the post-translational modifications of histone tails regulate the higher-order chromatin architecture. On the other hand, previous biochemical studies have suggested that histone tails also affect the structure and stability of the nucleosome core particle itself. However, the precise contributions of each histone tail are unclear. In the present study, we determined the crystal structures of four mutant nucleosomes, in which one of the four histones, H2A, H2B, H3, or H4, lacked the N-terminal tail. We found that the deletion of the H2B or H3 N-terminal tail affected histone-DNA interactions and substantially decreased nucleosome stability. These findings provide important information for understanding the complex roles of histone tails in regulating chromatin structure.