Derivation of rules for comparative protein modeling from a database of protein structure alignments

Nitro-benzylideneoxymorphone, a bifunctional mu and delta opioid receptor ligand with high mu opioid receptor efficacy

Introduction: There is a major societal need for analgesics with less tolerance, dependence, and abuse liability. Preclinical rodent studies suggest that bifunctional ligands with both mu (MOPr) and delta (DOPr) opioid peptide receptor activity may produce analgesia with reduced tolerance and other side effects. This study explores the structure-activity relationships (SAR) of our previously reported MOPr/DOPr lead, benzylideneoxymorphone (BOM) with C7-methylene-substituted analogs. Methods: Analogs were synthesized and tested in vitro for opioid receptor binding and efficacy. One compound, nitro-BOM (NBOM, 12) was evaluated for antinociceptive effects in the warm water tail withdrawal assay in C57BL/6 mice. Acute and chronic antinociception was determined, as was toxicologic effects on chronic administration. Molecular modeling experiments were performed using the Site Identification by Ligand Competitive Saturation (SILCS) method. Results: NBOM was found to be a potent MOPr agonist/DOPr partial agonist that produces high-efficacy antinociception. Antinociceptive tolerance was observed, as was weight loss; this toxicity was only observed with NBOM and not with BOM. Modeling supports the hypothesis that the increased MOPr efficacy of NBOM is due to the substituted benzylidene ring occupying a nonpolar region within the MOPr agonist state. Discussion: Though antinociceptive tolerance and non-specific toxicity was observed on repeated administration, NBOM provides an important new tool for understanding MOPr/DOPr pharmacology.

In silico analysis of single nucleotide polymorphism (rs34377097) of TBXA2R gene and pollen induced bronchial asthma susceptibility in West Bengal population, India

Introduction

Prevalence of asthma is increasing steadily among general population in developing countries over past two decades. One of the causative agents of broncho-constriction in asthma is thromboxane A2 receptor (TBXA2R). However few studies of TBXA2R polymorphism were performed so far. The present study aimed to assess potential association of TBXA2R rs34377097 polymorphism causing missense substitution of Arginine to Leucine (R60L) among 482 patients diagnosed with pollen-induced asthma and 122 control participants from West Bengal, India. Also we performed in-silico analysis of mutated TBXA2R protein (R60L) using homology modeling.

Methods

Clinical parameters like Forced expiratory volume in 1 second (FEV₁), FEV₁/Forced vital capacity (FVC) and Peak expiratory flow rate (PEFR) were assessed using spirometry. Patients' sensitivity was measured by skin prick test (SPT) against 16 pollen allergens. Polymerase chain reaction-based Restriction fragment length polymorphism was done for genotyping. Structural model of wild type and homology model of polymorphic TBXA2R was generated using AlphaFold2 and MODELLER respectively. Electrostatic surface potential was calculated using APBS plugin in PyMol.

Results

Genotype frequencies differed significantly between the study groups (P=0.03). There was no significant deviation from Hardy-Weinberg equilibrium in control population (χ2=1.56). Asthmatic patients have significantly higher frequency of rs34377097TT genotype than control subjects (P=0.03). SPT of patients showed maximum sensitivity in A. indica (87.68%) followed by C. nusifera (83.29%) and C. pulcherima (74.94%). Significant difference existed for pollen sensitivity in adolescent and young adult (P=0.01) and between young and old adult (P=0.0003). Significant negative correlation was found between FEV1/FVC ratio and intensity of SPT reactions (P<0.0001). Significant association of FEV₁, FEV₁/FVC and PEFR was observed with pollen-induced asthma. Furthermore, risk allele T was found to be clinically correlated with lower FEV₁/FVC ratio (P=0.015) in patients. Our data showed R60L polymorphism, which was conserved across mammals, significantly reduced positive electrostatic charge of polymorphic protein in cytoplasmic domain thus altered downstream pathway and induced asthma response.

Discussion

The present in-silico study is the first one to report association of TBXA2R rs34377097 polymorphism in an Indian population. It may be used as prognostic marker of clinical response to asthma in West Bengal and possible target of therapeutics in future.

Mechanistic insights to lactic and formic acid toxicity on benthic oligochaete worm Tubifex tubifex

Lactic and formic acid are two commonly found monocarboxylic organic acids. Lactic acid is discharged into the water bodies as acidic industrial effluent from the food, cosmetic, chemical, and pharmaceutical industries, whereas formic acid is discharged from various paper, leather tanning, and textile processing industries. The present study investigated the toxicity of both organic acids upon the benthic oligochaete worm Tubifex tubifex. The 96-h median lethal concentration (LC₅₀) values for lactic and formic acid are determined as 143.81 mg/l and 57.99 mg/l respectively. The effects of two sublethal concentrations (10% and 30% of 96 h LC₅₀) of these acids on differential expression of oxidative stress enzymes are investigated. The comparative analysis of acute toxicity demonstrates that formic acid exposure is more detrimental to T. tubifex than lactic acid. The in silico structural analysis predicts that formic acid can interact with cytochrome c oxidase of the electron transport system and thereby inhibits its functionality and induces reactive oxygen species production. Integrated biomarker response (IBR) analysis illustrates that overall oxidative stress of formic acid to T. tubifex is significantly higher than that of lactic acid, which supports the structural analysis. It is concluded from this study that toxicokinetic-toxicodynamic and species sensitivity distributions studies are helpful for ecological risk management of environmental toxicants.

Expression of Mutant Glycine Receptors in Xenopus Oocytes Using Canonical and Non-Canonical Amino Acids Reveals Distinct Roles of Conserved Proline Residues

Pentameric ligand-gated ion channels (pLGIC) play important roles in fast neuronal signal transmission. Functional receptors are pentamers, with each subunit having an extracellular domain (ECD), a transmembrane domain (TMD) and an intracellular domain. The binding of the agonist to the ECD induces a structural change that is transduced to the TMD to open the channel. Molecular details of this process are emerging, but a comprehensive understanding is still lacking. Proline (Pro) is one amino acid that has attracted much interest; its unusual features generate bends in loops and kinks and bulges in helices, which can be essential for function in some pLGICs. Here, we explore the roles of four conserved Pros in the glycine receptor (GlyR), creating substitutions with canonical and noncanonical amino acids, characterizing them using two electrode voltage clamp electrophysiology in Xenopus oocytes, and interpreting changes in receptor parameters using structural data from the open and closed states of the receptor. The data reveal that for efficient function, the Pro in the α1β1 loop is needed to create a turn and to be the correct size and shape to interact with nearby residues; the peptide bond of the Pro in the Cys-loop requires the cis conformation; and the Pros in loop A and M1 allow efficient function because of their reduced hydrogen bonding capacity. These data are broadly consistent with data from other pLGICs, and therefore likely represent the important features of these Pros in all members of the family.

Force Field Parameters for Fe2+4S2-4 Clusters of Dihydropyrimidine Dehydrogenase, the 5-Fluorouracil Cancer Drug Deactivation Protein: A Step towards In Silico Pharmacogenomics Studies

The dimeric dihydropyrimidine dehydrogenase (DPD), metalloenzyme, an adjunct anti-cancer drug target, contains highly specialized 4 × Fe²⁺₄S²⁻₄ clusters per chain. These clusters facilitate the catalysis of the rate-limiting step in the pyrimidine degradation pathway through a harmonized electron transfer cascade that triggers a redox catabolic reaction. In the process, the bulk of the administered 5-fluorouracil (5-FU) cancer drug is inactivated, while a small proportion is activated to nucleic acid antimetabolites. The occurrence of missense mutations in DPD protein within the general population, including those of African descent, has adverse toxicity effects due to altered 5-FU metabolism. Thus, deciphering mutation effects on protein structure and function is vital, especially for precision medicine purposes. We previously proposed combining molecular dynamics (MD) and dynamic residue network (DRN) analysis to decipher the molecular mechanisms of missense mutations in other proteins. However, the presence of Fe²⁺₄S²⁻₄ clusters in DPD poses a challenge for such in silico studies. The existing AMBER force field parameters cannot accurately describe the Fe²⁺ center coordination exhibited by this enzyme. Therefore, this study aimed to derive AMBER force field parameters for DPD enzyme Fe²⁺ centers, using the original Seminario method and the collation features Visual Force Field Derivation Toolkit as a supportive approach. All-atom MD simulations were performed to validate the results. Both approaches generated similar force field parameters, which accurately described the human DPD protein Fe²⁺₄S²⁻₄ cluster architecture. This information is crucial and opens new avenues for in silico cancer pharmacogenomics and drug discovery related research on 5-FU drug efficacy and toxicity issues.

An integrated clinical and molecular study of a cohort of Turkish patients with Marfan syndrome harboring known and novel FBN1 variants

Marfan syndrome (MFS) is an autosomal dominant genetic condition that mainly affects connective tissue in many parts of the body. Cardinal manifestations involve the ocular, skeletal, and cardiovascular systems. The diagnosis of MFS relies on the revised Ghent criteria, outlined by international expert opinion to facilitate accurate recognition of this syndrome as well as to improve patient management and counseling. However, it may not always be possible to make a definitive diagnosis according to these criteria in each patient and thus molecular confirmation is necessary in subjects with suspected MFS. This debilitating, if not fatal, disorder is caused by mutations in FBN1, which encodes a major constitutive element of extracellular microfibrils. Here, we present a detailed clinical and molecular analysis of 76 Turkish patients with definitive or suspected MFS diagnosed at our center between 2014 and 2019. We were able to identify a total of 51 different FBN1 variants in our cohort, 31 of which have previously been reported in the relevant scientific literature. The remaining 20 variants have not been documented to date. In one patient, we detected a large deletion including the entire FBN1 gene using the array CGH approach. Currently, there are very few studies on the genotype-phenotype correlation of patients with MFS, and no clear genotype-phenotype maps for MFS have been constructed so far, except for some cases. We believe that our findings will make a rich and peculiar contribution to the elusive genotype-phenotype relationship in MFS, especially in this large and populous ethnic group.

CALR-ETdb, the database of calreticulin variants diversity in essential thrombocythemia

Essential thrombocythemia (ET) is a blood cancer defined by a strong increase of platelet numbers. A quarter of patients suffering from ET show mutations in the last exon of calreticulin (CALR) gene. Two variants named type 1 and type 2 represent 85% of these patients. However, a large number of other variants have been determined. In this study, we have compiled variants taken from COSMIC database and literature leading to 155 different variants. This large number of variants allowed redefining 5 new classes extending the classification of type 1-like and type 2-like to a finer description. These analyses showed that last class, named E, corresponding to more than 10% of CALR variants seemed not attached to ET. Structural properties analyzed showed that CALR variants associated to ET have common features. All the compiled and refined information had been included into a freely dedicated database CALR-ETdb (https://www.dsimb.inserm.fr/CALR-ET).

Protein Structure Modeling with MODELLER

Genome sequencing projects have resulted in a rapid increase in the number of known protein sequences. In contrast, only about one-hundredth of these sequences have been characterized at atomic resolution using experimental structure determination methods. Computational protein structure modeling techniques have the potential to bridge this sequence-structure gap. In the following chapter, we present an example that illustrates the use of MODELLER to construct a comparative model for a protein with unknown structure. Automation of a similar protocol has resulted in models of useful accuracy for domains in more than half of all known protein sequences.

Proline Residues Contribute to Efficient GABAp Receptor Function

GABAp receptors are homomeric pentameric ligand-gated ion channels (pLGICs) and are useful for probing the molecular details of the mechanism of action in this important protein family. Here, we explore the role of proline (Pro) residues by creating mutant receptors, expressing them in HEK293 cells, and using fluorescent membrane potential sensitive dye to monitor receptor activity. The data revealed that 3 of the Pro-to-alanine substitutions resulted in nonfunctional receptors (one in the Cys-loop, one in loop A and one in the β2-β3 loop), 7 resulted in increased EC₅₀ values, and the remaining 13 resulted in receptors with properties similar to WT receptors. Further exploration of the nonfunctional mutants using expression in Xenopus laevis oocytes and whole-cell voltage-clamp electrophysiology, incorporating both canonical and noncanonical amino acids, revealed that the Pro in the Cys-loop had a preference for analogues with a high intrinsic cis bias, the Pro in loop A required a ring, and the β2-β3 loop Pro contributes to expression. The data from the Cys-loop Pro are consistent with those from other pLGICs, while those of the loop A Pro and some of the other Pros surprisingly differ. Thus, overall, the data clarify the roles of many of the Pros in the GABAp receptor and also suggest that caution must be applied in using data from one receptor to understand molecular details of all pLGICs.

Many Proline Residues in the Extracellular Domain Contribute to Glycine Receptor Function

Prolines in signaling proteins are of particular interest because they have a range of unique properties that may be critical for function. Here we show that many proline residues in the extracellular domain (ECD) of the glycine receptor are involved in the correct functioning of this ligand-gated ion channel. We explore their role by creating mutant receptors, expressing them in cells, and using fluorescent membrane potential sensitive dye to monitor receptor activity. We then interpret the changes in receptor parameters using structural information from the open and closed states of the receptor. The data reveal that substitution with alanine of ten of the 13 Pro residues in the ECD alters the function of the receptor: one substitution ablates function, six cause a decrease in the EC₅₀, and three cause an increase. Only three of these mutants result in EC₅₀ values similar to WT. The nonfunctional mutant, Pro30Ala, was further probed in oocytes, and the data suggest a role in both expression and function. Examination of the locations of sensitive Pro residues in the receptor and identification of potential interactions with nearby residues reveal how these residues could contribute to the correct functioning of this typical pentameric ligand-gated ion channel.

Characterizing the binding and function of TARP γ8-selective AMPA receptor modulators

α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid(AMPA)-type glutamate receptors (AMPARs) are the predominant excitatory neurotransmitter receptors in the brain, where they mediate synaptic transmission and plasticity. Excessive AMPAR activation leads to diseases such as epilepsy. AMPAR properties are modulated by auxiliary proteins and foremost by the transmembrane AMPAR regulatory proteins (TARPs). These distribute in unique expression patterns across the brain, rendering AMPAR/TARP complexes promising targets for region-specific therapeutic intervention. TARP γ8 is predominantly expressed in the forebrain and is enriched in the hippocampus, a region associated with temporal lobe epilepsy. Recent high-throughput medicinal chemistry screens have identified multiple promising compounds that selectively target AMPARs associated with γ8 and hold promise for epilepsy treatment. However, how these modulators target the receptor complex is currently unknown. Here, we use a combination of ligand docking, molecular dynamics simulations, and electrophysiology to address this question. We identify a conserved oxindole isostere, shared between three structurally diverse modulators (LY-3130481, JNJ-55511118, and JNJ-61432059) as the major module engaging γ8 by an H-bond to Asn-172 (γ8). The remaining variable region of each molecule likely targets the receptor complex in ligand-selective modes. Functional data reveal parallels in the underlying modulatory action of two prominent compounds. This work will aid development of refined AMPAR epilepsy therapeutics and facilitate to uncover the mechanisms by which TARPs modulate the receptor.

Glutantβase: a database for improving the rational design of glucose-tolerant β-glucosidases

Β-glucosidases are key enzymes used in second-generation biofuel production. They act in the last step of the lignocellulose saccharification, converting cellobiose in glucose. However, most of the β-glucosidases are inhibited by high glucose concentrations, which turns it a limiting step for industrial production. Thus, β-glucosidases have been targeted by several studies aiming to understand the mechanism of glucose tolerance, pH and thermal resistance for constructing more efficient enzymes. In this paper, we present a database of β-glucosidase structures, called Glutantβase. Our database includes 3842 GH1 β-glucosidase sequences collected from UniProt. We modeled the sequences by comparison and predicted important features in the 3D-structure of each enzyme. Glutantβase provides information about catalytic and conserved amino acids, residues of the coevolution network, protein secondary structure, and residues located in the channel that guides to the active site. We also analyzed the impact of beneficial mutations reported in the literature, predicted in analogous positions, for similar enzymes. We suggested these mutations based on six previously described mutants that showed high catalytic activity, glucose tolerance, or thermostability (A404V, E96K, H184F, H228T, L441F, and V174C). Then, we used molecular docking to verify the impact of the suggested mutations in the affinity of protein and ligands (substrate and product). Our results suggest that only mutations based on the H228T mutant can reduce the affinity for glucose (product) and increase affinity for cellobiose (substrate), which indicates an increment in the resistance to product inhibition and agrees with computational and experimental results previously reported in the literature. More resistant β-glucosidases are essential to saccharification in industrial applications. However, thermostable and glucose-tolerant β-glucosidases are rare, and their glucose tolerance mechanisms appear to be related to multiple and complex factors. We gather here, a set of information, and made predictions aiming to provide a tool for supporting the rational design of more efficient β-glucosidases. We hope that Glutantβase can help improve second-generation biofuel production. Glutantβase is available at http://bioinfo.dcc.ufmg.br/glutantbase .

Interaction of LysM BON family protein domain with carbapenems: A putative mechanism of carbapenem resistance

Carbapenem resistance in Gram-negative pathogens has become a global concern for health workers worldwide. In one of our earlier studies, a Klebsiella pneumoniae-carbapenemase-2 producing strain was induced with meropenem to explore differentially expressed proteins under induced and uninduced conditions. There is, LysM domain BON family protein, was found over 12-fold expressed under the induced state. A hypothesis was proposed that LysM domain protein might have an affinity towards carbapenem antibiotics making them unavailable to bind with their target. Hence, we initiated a study to understand the binding mode of carbapenem with LysM domain protein. MICs of imipenem and meropenem against LysM clone were increased by several folds as compared to NP-6 clinical strain as well as DH5 α (PET-28a KPC-2) clone. This study further revealed a strong binding of both antibiotics to LysM domain protein. Molecular simulation studies of LysM domain protein with meropenem and imipenem for 80 ns has also showed stable structure. We concluded that overexpressed LysM domain under induced condition interacted with carbapenems, leading to enhanced resistance as proved by high MIC values. Hence, the study proved the proposed hypothesis that the LysM domain plays a significant role in the putative mechanism of antibiotics resistance.

Negative allosteric modulation of GluN1/GluN3 NMDA receptors

NMDA receptors are ligand-gated ion channels that mediate excitatory neurotransmission. Most native NMDA receptors are tetrameric assemblies of two glycine-binding GluN1 and two glutamate-binding GluN2 subunits. Co-assembly of the glycine-binding GluN1 with glycine-binding GluN3 subunits (GluN3A-B) creates glycine activated receptors that possess strikingly different functional and pharmacological properties compared to GluN1/GluN2 NMDA receptors. The role of GluN1/GluN3 receptors in neuronal function remains unknown, in part due to lack of pharmacological tools with which to explore their physiological roles. We have identified the negative allosteric modulator EU1180-438, which is selective for GluN1/GluN3 receptors over GluN1/GluN2 NMDA receptors, AMPA, and kainate receptors. EU1180-438 is also inactive at GABA, glycine, and P2X receptors, but displays inhibition of some nicotinic acetylcholine receptors. Furthermore, we demonstrate that EU1180-438 produces robust inhibition of glycine-activated current responses mediated by native GluN1/GluN3A receptors in hippocampal CA1 pyramidal neurons. EU1180-438 is a non-competitive antagonist with activity that is independent of membrane potential (i.e. voltage-independent), glycine concentration, and extracellular pH. Non-stationary fluctuation analysis of neuronal current responses provided an estimated weighted mean unitary conductance of 6.1 pS for GluN1/GluN3A channels, and showed that EU1180-438 has no effect on conductance. Site-directed mutagenesis suggests that structural determinants of EU1180-438 activity reside near a short pre-M1 helix that lies parallel to the plane of the membrane below the agonist binding domain. These findings demonstrate that structural differences between GluN3 and other glutamate receptor subunits can be exploited to generate subunit-selective ligands with utility in exploring the roles GluN3 in neuronal function.

Loxoscelism: Advances and Challenges in the Design of Antibody Fragments with Therapeutic Potential

Envenoming due to Loxosceles spider bites still remains a neglected disease of particular medical concern in the Americas. To date, there is no consensus for the treatment of envenomed patients, yet horse polyclonal antivenoms are usually infused to patients with identified severe medical conditions. It is widely known that venom proteins in the 30–35 kDa range with sphingomyelinase D (SMasesD) activity, reproduce most of the toxic effects observed in loxoscelism. Hence, we believe that monoclonal antibody fragments targeting such toxins might pose an alternative safe and effective treatment. In the present study, starting from the monoclonal antibody LimAb7, previously shown to target SMasesD from the venom of L. intermedia and neutralize its dermonecrotic activity, we designed humanized antibody V-domains, then produced and purified as recombinant single-chain antibody fragments (scFvs). These molecules were characterized in terms of humanness, structural stability, antigen-binding activity, and venom-neutralizing potential. Throughout this process, we identified some blocking points that can impact the Abs antigen-binding activity and neutralizing capacity. In silico analysis of the antigen/antibody amino acid interactions also contributed to a better understanding of the antibody’s neutralization mechanism and led to reformatting the humanized antibody fragment which, ultimately, recovered the functional characteristics for efficient in vitro venom neutralization.

Influence of the Structural Accuracy of Homology Models on Their Applicability to Docking-Based Virtual Screening: The β2 Adrenergic Receptor as a Case Study

How accurate do structures of the β₂ adrenergic receptor (β₂AR) need to be to effectively serve as platforms for docking-based virtual screening campaigns? To answer this research question, here, we targeted through controlled virtual screening experiments 23 homology models of the β₂AR endowed with different levels of structural accuracy. Subsequently, we studied the correlation between virtual screening performance and structural accuracy of the targeted models. Moreover, we studied the correlation between virtual screening performance and template/target receptor sequence identity. Our study demonstrates that docking-based virtual screening campaigns targeting homology models of the β₂AR, in the majority of the cases, yielded results that exceeded random expectations in terms of area under the receiver operating characteristic curve (ROC AUC). Moreover, with the most effective scoring method, over one-third and one-quarter of the models yielded results that exceeded random expectation also in terms of enrichment factors (EF1, EF5, and EF10) and BEDROC (α = 160.9), respectively. Not surprisingly, we found a detectable linear correlation between virtual screening performance and structural accuracy of the ligand-binding cavity. We also found a detectable linear correlation between virtual screening performance and structural accuracy of the second extracellular loop (EL2). Finally, our data indicate that, although there is no detectable linear correlation between virtual screening performance and template/β₂AR sequence identity, models built on the basis of templates that show high sequence identity with the β₂AR, especially within the ligand-biding cavity, performed consistently well. Conversely, models with lower sequence identity displayed performance levels that ranged from very good to random, with no apparent correlation with the sequence identity itself.

Investigating targets for neuropharmacological intervention by molecular dynamics simulations

Medical research has identified over 500 brain disorders. Among these, there are still only very few neuropathologies whose causes are fully understood and, consequently, very few drugs whose mechanism of action is known. No FDA drug has been identified for major neurodegenerative diseases, such as Alzheimer's and Parkinson's. We still lack effective treatments and strategies for modulating progression or even early neurodegenerative disease onset diagnostic tools. A great support toward the highly needed identification of neuroactive drugs comes from computer simulation methods and, in particular, from molecular dynamics (MD). This provides insight into structure-function relationship of a target and predicts structure, dynamics and energetics of ligand/target complexes under biologically relevant conditions like temperature and physiological saline concentration. Here, we present examples of the predictive power of MD for neuroactive ligands/target complexes. This brief survey from our own research shows the usefulness of partnerships between academia and industry, and from joint efforts between experimental and theoretical groups.

A Modifying Autoantigen in Graves' Disease

The TSH receptor (TSHR) is the major autoantigen in Graves' disease (GD). Bioinformatic analyses predict the existence of several human TSHR isoforms from alternative splicing, which can lead to the coexpression of multiple receptor forms. The most abundant of these is TSHRv1.3. In silico modeling of TSHRv1.3 demonstrated the structural integrity of this truncated receptor isoform and its potential binding of TSH. Tissue profiling revealed wide expression of TSHRv1.3, with a predominant presence in thyroid, bone marrow, thymus, and adipose tissue. To gain insight into the role of this v1.3 receptor isoform in thyroid pathophysiology, we cloned the entire open reading frame into a mammalian expression vector. Immunoprecipitation studies demonstrated that both TSHR-stimulating antibody and human TSH could bind v1.3. Furthermore, TSHRv1.3 inhibited the stimulatory effect of TSH and TSHR-Ab MS-1 antibody on TSHR-induced cAMP generation in a dose-dependent manner. To confirm the antigenicity of v1.3, we used a peptide ELISA against two different epitopes. Of 13 GD samples, 11 (84.6%) were positive for a carboxy terminal peptide and 10 (76.9%) were positive with a junction region peptide. To demonstrate that intracellular v1.3 could serve as an autoantigen and modulate disease, we used double-transfected Chinese hamster ovary cells that expressed both green fluorescent protein (GFP)-tagged TSHRv1.3 and full-length TSHR. We then induced cell stress and apoptosis using a TSHR monoclonal antibody and observed the culture supernatant contained v1.3-GFP protein, demonstrating the release of the intracellular receptor variant by this mechanism.

Reduced Activity of Nitrate Reductase Under Heavy Metal Cadmium Stress in Rice: An in silico Answer

Cadmium is a well known toxic heavy metal, which has various detrimental effects on plant system. In plants an important enzyme involved in the production of nitric oxide, nitrate reductase, is also affected by cadmium toxicity. According to many studies cadmium has an inhibitory effect on nitrate reductase activity. Similar effect of cadmium was found in our study where an inhibitory effect of cadmium on nitrate reductase activity was noted. However, the mechanism behind this inhibition has not been explored. With the help of homology, 3-D structure of rice-nitrate reductase is modeled in this study. Its binding with nitrate, nitrite and cadmium metal in silico has been explored. The bonds formed between the enzyme-substrate complex, enzyme-cadmium and differences in interactions in presence of cadmium has been studied in detail. The present study should help in understanding the modeled structure of rice-nitrate reductase in 3-D which may in turn guide enzyme related studies in silico. The present study also provides an insight as to how cadmium interacts with nitrate reductase to alter the enzyme activity.

The S6 gate in regulatory Kv6 subunits restricts heteromeric K+ channel stoichiometry

Atypical substitutions in the S6 activation gate sequence distinguish “regulatory” Kv subunits, which cannot homotetramerize due to T1 self-incompatibility. Pisupati et al. show that such substitutions in Kv6 work together with self-incompatibility to restrict Kv2:Kv6 heteromeric stoichiometry to 3:1.

The Shaker-like family of voltage-gated K⁺ channels comprises four functionally independent gene subfamilies, Shaker (Kv1), Shab (Kv2), Shaw (Kv3), and Shal (Kv4), each of which regulates distinct aspects of neuronal excitability. Subfamily-specific assembly of tetrameric channels is mediated by the N-terminal T1 domain and segregates Kv1–4, allowing multiple channel types to function independently in the same cell. Typical Shaker-like Kv subunits can form functional channels as homotetramers, but a group of mammalian Kv2-related genes (Kv5.1, Kv6s, Kv8s, and Kv9s) encodes subunits that have a “silent” or “regulatory” phenotype characterized by T1 self-incompatibility. These channels are unable to form homotetramers, but instead heteromerize with Kv2.1 or Kv2.2 to diversify the functional properties of these delayed rectifiers. While T1 self-incompatibility predicts that these heterotetramers could contain up to two regulatory (R) subunits, experiments show a predominance of 3:1R stoichiometry in which heteromeric channels contain a single regulatory subunit. Substitution of the self-compatible Kv2.1 T1 domain into the regulatory subunit Kv6.4 does not alter the stoichiometry of Kv2.1:Kv6.4 heteromers. Here, to identify other channel structures that might be responsible for favoring the 3:1R stoichiometry, we compare the sequences of mammalian regulatory subunits to independently evolved regulatory subunits from cnidarians. The most widespread feature of regulatory subunits is the presence of atypical substitutions in the highly conserved consensus sequence of the intracellular S6 activation gate of the pore. We show that two amino acid substitutions in the S6 gate of the regulatory subunit Kv6.4 restrict the functional stoichiometry of Kv2.1:Kv6.4 to 3:1R by limiting the formation and function of 2:2R heteromers. We propose a two-step model for the evolution of the asymmetric 3:1R stoichiometry, which begins with evolution of self-incompatibility to establish the regulatory phenotype, followed by drift of the activation gate consensus sequence under relaxed selection to limit stoichiometry to 3:1R.

Chiral analogues of (+)-cyclazosin as potent α1B-adrenoceptor selective antagonist

(+)-Cyclazosin [(+)-1] is one of most selective antagonists of the α_1B-adrenoceptor subtype (selectivity ratios, α_1B/α_1A = 13, α_1B/α_1D = 38-39). To improve the selectivity, we synthesized and pharmacologically studied the blocking activity against α₁-adrenoceptors of several homochiral analogues of (+)-cyclazosin featuring different substituents on the carbonyl or amine groups, namely (-)-2, (+)-3, (-)-4-(-)-8, (+)-9. Moreover, we studied the activity of some their opposite enantiomers, namely (-)-1, (-)-3, (+)-6, and (-)-9, to evaluate the influence of stereochemistry on selectivity. The benzyloxycarbonyl and methyl (4aS,8aR) analogues (+)-3 and (-)-6 improved in a significant way the α_1B selectivity of the progenitor compound: 4 and 14 time vs. the α_1D subtype and 35 and 77 times vs. the α_1A subtype, respectively. The study confirmed the importance of the hydrophobic cis-octahydroquinoxaline moiety of these molecules for the establishment of interactions with the α₁-adrenoceptors as well that of their (4aS,8aR) stereochemistry to grant selectivity for the α_1B subtype. Hypotheses on the mode of interaction of these compounds were advanced on the basis of molecular modeling studies performed on compound (+)-3.

Broadly conserved Na+-binding site in the N-lobe of prokaryotic multidrug MATE transporters

Multidrug and toxic-compound extrusion (MATE) proteins comprise an important but largely uncharacterized family of secondary-active transporters. In both eukaryotes and prokaryotes, these transporters protect the cell by catalyzing the efflux of a broad range of cytotoxic compounds, including human-made antibiotics and anticancer drugs. MATEs are thus potential pharmacological targets against drug-resistant pathogenic bacteria and tumor cells. The activity of MATEs is powered by transmembrane electrochemical ion gradients, but their molecular mechanism and ion specificity are not understood, in part because high-quality structural information is limited. Here, we use computational methods to study PfMATE, from Pyrococcus furiosus, whose structure is the best resolved to date. Analysis of available crystallographic data and additional molecular dynamics simulations unequivocally reveal an occupied Na⁺-binding site in the N-lobe of this transporter, which had not been previously recognized. We find this site to be selective against K⁺ and broadly conserved among prokaryotic MATEs, including homologs known to be Na⁺-dependent such as NorM-VC, VmrA, and ClbM, for which the location of the Na⁺ site had been debated. We note, however, that the chemical makeup of the proposed Na⁺ site indicates it is weakly specific against H⁺, explaining why MATEs featuring this Na⁺-binding motif may be solely driven by H⁺ in laboratory conditions. We further posit that the concurrent coupling to H⁺ and Na⁺ gradients observed for some Na⁺-driven MATEs owes to a second H⁺-binding site, within the C-lobe. In summary, our study provides insights into the structural basis for the complex ion dependency of MATE transporters.

Molecular Modeling of Multidrug Properties of Resistance Nodulation Division (RND) Transporters

Efflux pumps of the resistance nodulation division (RND) superfamily are among the major contributors to intrinsic and acquired multidrug resistance in Gram-negative bacteria. Structural information on AcrAB-TolC and MexAB-OprM, major efflux pumps of Escherichia coli and Pseudomonas aeruginosa respectively, boosted intensive research aimed at understanding the molecular mechanisms ruling the active extrusion processes. In particular, several studies were devoted to the understanding of the determinants behind the extraordinary broad specificity of the RND transporters AcrB and MexB. In this chapter, we discuss the ever-growing role computational methods have been playing in deciphering key structural and dynamical features of these transporters and of their interaction with substrates and inhibitors. We further discuss and illustrate examples from our lab of how molecular docking, homology modeling, all-atom molecular dynamics simulations and in silico free energy estimations can all together give precious insights into the processes of recognition and extrusion of substrates, as well as on the possible inhibition strategies.

The Bioactive Protein-Ligand Conformation of GluN2C-Selective Positive Allosteric Modulators Bound to the NMDA Receptor

N-methyl-d-aspartate (NMDA) receptors are ligand-gated, cation-selective channels that mediate a slow component of excitatory synaptic transmission. Subunit-selective positive allosteric modulators of NMDA receptor function have therapeutically relevant effects on multiple processes in the brain. A series of pyrrolidinones, such as PYD-106, that selectively potentiate NMDA receptors that contain the GluN2C subunit have structural determinants of activity that reside between the GluN2C amino terminal domain and the GluN2C agonist binding domain, suggesting a unique site of action. Here we use molecular biology and homology modeling to identify residues that line a candidate binding pocket for GluN2C-selective pyrrolidinones. We also show that occupancy of only one site in diheteromeric receptors is required for potentiation. Both GluN2A and GluN2B can dominate the sensitivity of triheteromeric receptors to eliminate the actions of pyrrolidinones, thus rendering this series uniquely sensitive to subunit stoichiometry. We experimentally identified NMR-derived conformers in solution, which combined with molecular modeling allows the prediction of the bioactive binding pose for this series of GluN2C-selective positive allosteric modulators of NMDA receptors. These data advance our understanding of the site and nature of the ligand-protein interaction for GluN2C-selective positive allosteric modulators for NMDA receptors.

Asymmetry of inverted-topology repeats in the AE1 anion exchanger suggests an elevator-like mechanism

Anion exchanger 1 catalyzes the transmembrane antiport of chloride and bicarbonate ions through a mechanism that has remained unclear. By modeling its inward-facing state and comparing it with the known outward-facing form, Ficici et al. hypothesize that this transporter features an elevator-like mechanism.

The membrane transporter anion exchanger 1 (AE1), or band 3, is a key component in the processes of carbon-dioxide transport in the blood and urinary acidification in the renal collecting duct. In both erythrocytes and the basolateral membrane of the collecting-duct α-intercalated cells, the role of AE1 is to catalyze a one-for-one exchange of chloride for bicarbonate. After decades of biochemical and functional studies, the structure of the transmembrane region of AE1, which catalyzes the anion-exchange reaction, has finally been determined. Each protomer of the AE1 dimer comprises two repeats with inverted transmembrane topologies, but the structures of these repeats differ. This asymmetry causes the putative substrate-binding site to be exposed only to the extracellular space, consistent with the expectation that anion exchange occurs via an alternating-access mechanism. Here, we hypothesize that the unknown, inward-facing conformation results from inversion of this asymmetry, and we propose a model of this state constructed using repeat-swap homology modeling. By comparing this inward-facing model with the outward-facing experimental structure, we predict that the mechanism of AE1 involves an elevator-like motion of the substrate-binding domain relative to the nearly stationary dimerization domain and to the membrane plane. This hypothesis is in qualitative agreement with a wide range of biochemical and functional data, which we review in detail, and suggests new avenues of experimentation.

KLK14 interactions with HAI-1 and HAI-2 serine protease inhibitors: A molecular dynamics and relative free-energy calculations study

Kallikrein 14 (KLK14) is a serine protease linked to several pathologies including prostate cancer and positively correlates with Gleason score. Though KLK14 functioning in cancer is poorly understood, it has been implicated in HGF/Met signaling, given that KLK14 proteolytically inhibits HGF activator-inhibitor 1 (HAI-1), which strongly inhibits pro-HGF activators, thereby contributing to tumor progression. In this work, KLK14 binding to either hepatocyte growth factor activator inhibitor type-1 (HAI-1) or type-2 (HAI-2) was essayed using homology modeling, molecular dynamic simulations and free-energy calculations through MM/PBSA and MM/GBSA. KLK14 was successfully modeled. Calculated free energies suggested higher binding affinity for the KLK14/HAI-1 interaction than for KLK14/HAI-2. This difference in binding affinity is largely explained by the higher stability of the hydrogen-bond networks in KLK14/HAI-1 along the simulation trajectory. A key arginine residue in both HAI-1 and HAI-2 is responsible for their interaction with the S1 pocket in KLK14. Additionally, MM/GBSA free-energy decomposition postulates that KLK14 Asp174 and Trp196 are hotspots for binding HAI-1 and HAI-2.

Constraining the timing of the Great Oxidation Event within the Rubisco phylogenetic tree

Ribulose 1,5-bisphosphate (RuBP) carboxylase/oxygenase (RuBisCO, or Rubisco) catalyzes a key reaction by which inorganic carbon is converted into organic carbon in the metabolism of many aerobic and anaerobic organisms. Across the broader Rubisco protein family, homologs exhibit diverse biochemical characteristics and metabolic functions, but the evolutionary origins of this diversity are unclear. Evidence of the timing of Rubisco family emergence and diversification of its different forms has been obscured by a meager paleontological record of early Earth biota, their subcellular physiology and metabolic components. Here, we use computational models to reconstruct a Rubisco family phylogenetic tree, ancestral amino acid sequences at branching points on the tree, and protein structures for several key ancestors. Analysis of historic substitutions with respect to their structural locations shows that there were distinct periods of amino acid substitution enrichment above background levels near and within its oxygen-sensitive active site and subunit interfaces over the divergence between Form III (associated with anoxia) and Form I (associated with oxia) groups in its evolutionary history. One possible interpretation is that these periods of substitutional enrichment are coincident with oxidative stress exerted by the rise of oxygenic photosynthesis in the Precambrian era. Our interpretation implies that the periods of Rubisco substitutional enrichment inferred near the transition from anaerobic Form III to aerobic Form I ancestral sequences predate the acquisition of Rubisco by fully derived cyanobacterial (i.e., dual photosystem-bearing, oxygen-evolving) clades. The partitioning of extant lineages at high clade levels within our Rubisco phylogeny indicates that horizontal transfer of Rubisco is a relatively infrequent event. Therefore, it is possible that the mutational enrichment periods between the Form III and Form I common ancestral sequences correspond to the adaptation of key oxygen-sensitive components of Rubisco prior to, or coincident with, the Great Oxidation Event.

Identification of novel inhibitors of the translationally controlled tumor protein (TCTP): insights from molecular dynamics

The translationally controlled tumor protein (TCTP) is a highly conserved multifunctional protein, preferentially expressed in mitotically active tissues and is a potential biomarker and a therapeutic target for lung cancers. An understanding of the biology of this molecule and model systems for the screening of drugs is still awaited. In the absence of complete crystal structure, NMR structures as templates were used for homology modeling and MD optimization of both Dictyostelium discoideum and human TCTPs, which was followed by pocket-site prediction, ligand screening and docking. Rescoring of TCTP-ligand complexes was done using MD and MM-PBSA approaches. D. discoideum TCTP was expressed under a constitutive promoter and the endogenous RNA in multicellular structures formed was localized by in situ hybridization. Based on the interactions and binding energy scores, two novel compounds were identified as the best potential inhibitors that could be further used for the development of drug candidates. Inhibition of cell proliferation was observed in the strain overexpressing Dictyostelium TCTP and in situ hybridization results show them to be localized in the prestalk (dying cell population) cells. D. discoideum and human TCTPs share similar dynamic behaviors; overexpression of Dictyostelium TCTP inhibits cell proliferation. D. discoideum could be used as a model system for understanding the biology of this molecule and also for drug screening.

Parallel homodimer structures of the extracellular domains of the voltage-gated sodium channel β4 subunit explain its role in cell-cell adhesion

Voltage-gated sodium channels (VGSCs) are transmembrane proteins required for the generation of action potentials in excitable cells and essential for propagating electrical impulses along nerve cells. VGSCs are complexes of a pore-forming α subunit and auxiliary β subunits, designated as β1/β1B-β4 (encoded by SCN1B-4B, respectively), which also function in cell-cell adhesion. We previously reported the structural basis for the trans homophilic interaction of the β4 subunit, which contributes to its adhesive function. Here, using crystallographic and biochemical analyses, we show that the β4 extracellular domains directly interact with each other in a parallel manner that involves an intermolecular disulfide bond between the unpaired Cys residues (Cys⁵⁸) in the loop connecting strands B and C and intermolecular hydrophobic and hydrogen-bonding interactions of the N-terminal segments (Ser³⁰-Val³⁵). Under reducing conditions, an N-terminally deleted β4 mutant exhibited decreased cell adhesion compared with the wild type, indicating that the β4 cis dimer contributes to the trans homophilic interaction of β4 in cell-cell adhesion. Furthermore, this mutant exhibited increased association with the α subunit, indicating that the cis dimerization of β4 affects α-β4 complex formation. These observations provide the structural basis for the parallel dimer formation of β4 in VGSCs and reveal its mechanism in cell-cell adhesion.

Functional Evaluation of a De Novo GRIN2A Mutation Identified in a Patient with Profound Global Developmental Delay and Refractory Epilepsy

The N-methyl-d-aspartate receptor (NMDAR), a ligand-gated ionotropic glutamate receptor, plays important roles in normal brain development and a wide range of neurologic disorders, including epilepsy. Here, we evaluate for the first time the functional properties of a de novo GRIN2A missense mutation (p.M817V) in the pre-M4 linker in a child with profound global developmental delay and refractory epilepsy. Electrophysiologic recordings revealed that the mutant GluN2A(M817V)-containing receptors showed enhanced agonist potency, reduced sensitivity to endogenous negative inhibitors (Mg²⁺, proton, and zinc), prolonged synaptic-like response time course, increased single-channel mean open time, and increased channel open probability. These results suggest that the gain-of-function M817V mutation causes overactivation of NMDAR and drives neuronal hyperexcitability, which may contribute to the patient's observed epileptic phenotype. Molecular modeling of the closed channel conformation reveals that this mutation weakens the interaction between GluN2 transmembrane helix M4 and two GluN1 transmembrane helices, and increases atomic fluctuation or movement of the pre-M1 region of GluN1 subunit, suggesting a mechanism by which channel function is enhanced. The functional changes of this mutation on agonist potency occur when the mutation is introduced into all other GluN2 subunits, suggesting a conserved role of this residue in control of NMDAR function through interactions of membrane spanning GluN2 and GluN1 helices. A number of NMDAR-targeted drugs including U.S. Food and Drug Association-approved NMDAR channel blockers were evaluated for their ability to inhibit receptors containing GluN2A(M817V) as a first step to exploring the potential for rescue pharmacology and personalized medicine.

Water transport through the intestinal epithelial barrier under different osmotic conditions is dependent on LI-cadherin trans-interaction

In the intestine water has to be reabsorbed from the chymus across the intestinal epithelium. The osmolarity within the lumen is subjected to high variations meaning that water transport often has to take place against osmotic gradients. It has been hypothesized that LI-cadherin is important in this process by keeping the intercellular cleft narrow facilitating the buildup of an osmotic gradient allowing water reabsorption. LI-cadherin is exceptional among the cadherin superfamily with respect to its localization along the lateral plasma membrane of epithelial cells being excluded from adherens junction. Furthermore it has 7 but not 5 extracellular cadherin repeats (EC1-EC7) and a small cytosolic domain. In this study we identified the peptide VAALD as an inhibitor of LI-cadherin trans-interaction by modeling the structure of LI-cadherin and comparison with the known adhesive interfaces of E-cadherin. This inhibitory peptide was used to measure LI-cadherin dependency of water transport through a monolayer of epithelial CACO2 cells under various osmotic conditions. If LI-cadherin trans-interaction was inhibited by use of the peptide, water transport from the luminal to the basolateral side was impaired and even reversed in the case of hypertonic conditions whereas no effect could be observed at isotonic conditions. These data are in line with a recently published model predicting LI-cadherin to keep the width of the lateral intercellular cleft small. In this narrow cleft a high osmolarity can be achieved due to ion pumps yielding a standing osmotic gradient allowing water absorption from the gut even if the faeces is highly hypertonic.

Protein Structure Modeling with MODELLER

Genome sequencing projects have resulted in a rapid increase in the number of known protein sequences. In contrast, only about one-hundredth of these sequences have been characterized at atomic resolution using experimental structure determination methods. Computational protein structure modeling techniques have the potential to bridge this sequence-structure gap. In the following chapter, we present an example that illustrates the use of MODELLER to construct a comparative model for a protein with unknown structure. Automation of a similar protocol has resulted in models of useful accuracy for domains in more than half of all known protein sequences.

The Diversity of Yellow-Related Proteins in Sand Flies (Diptera: Psychodidae)

Yellow-related proteins (YRPs) present in sand fly saliva act as affinity binders of bioamines, and help the fly to complete a bloodmeal by scavenging the physiological signals of damaged cells. They are also the main antigens in sand fly saliva and their recombinant form is used as a marker of host exposure to sand flies. Moreover, several salivary proteins and plasmids coding these proteins induce strong immune response in hosts bitten by sand flies and are being used to design protecting vaccines against Leishmania parasites. In this study, thirty two 3D models of different yellow-related proteins from thirteen sand fly species of two genera were constructed based on the known protein structure from Lutzomyia longipalpis. We also studied evolutionary relationships among species based on protein sequences as well as sequence and structural variability of their ligand-binding site. All of these 33 sand fly YRPs shared a similar structure, including a unique tunnel that connects the ligand-binding site with the solvent by two independent paths. However, intraspecific modifications found among these proteins affects the charges of the entrances to the tunnel, the length of the tunnel and its hydrophobicity. We suggest that these structural and sequential differences influence the ligand-binding abilities of these proteins and provide sand flies with a greater number of YRP paralogs with more nuanced answers to bioamines. All these characteristics allow us to better evaluate these proteins with respect to their potential use as part of anti-Leishmania vaccines or as an antigen to measure host exposure to sand flies.

Comparative Protein Structure Modeling Using MODELLER

Comparative protein structure modeling predicts the three-dimensional structure of a given protein sequence (target) based primarily on its alignment to one or more proteins of known structure (templates). The prediction process consists of fold assignment, target-template alignment, model building, and model evaluation. This unit describes how to calculate comparative models using the program MODELLER and how to use the ModBase database of such models, and discusses all four steps of comparative modeling, frequently observed errors, and some applications. Modeling lactate dehydrogenase from Trichomonas vaginalis (TvLDH) is described as an example. The download and installation of the MODELLER software is also described. © 2016 by John Wiley & Sons, Inc.

Homology modeling of a Class A GPCR in the inactive conformation: A quantitative analysis of the correlation between model/template sequence identity and model accuracy

With the present work we quantitatively studied the modellability of the inactive state of Class A G protein-coupled receptors (GPCRs). Specifically, we constructed models of one of the Class A GPCRs for which structures solved in the inactive state are available, namely the β₂ AR, using as templates each of the other class members for which structures solved in the inactive state are also available. Our results showed a detectable linear correlation between model accuracy and model/template sequence identity. This suggests that the likely accuracy of the homology models that can be built for a given receptor can be generally forecasted on the basis of the available templates. We also probed whether sequence alignments that allow for the presence of gaps within the transmembrane domains to account for structural irregularities afford better models than the classical alignment procedures that do not allow for the presence of gaps within such domains. As our results indicated, although the overall differences are very subtle, the inclusion of internal gaps within the transmembrane domains has a noticeable a beneficial effect on the local structural accuracy of the domain in question.

Comparative Protein Structure Modeling Using MODELLER

Comparative protein structure modeling predicts the three-dimensional structure of a given protein sequence (target) based primarily on its alignment to one or more proteins of known structure (templates). The prediction process consists of fold assignment, target-template alignment, model building, and model evaluation. This unit describes how to calculate comparative models using the program MODELLER and how to use the ModBase database of such models, and discusses all four steps of comparative modeling, frequently observed errors, and some applications. Modeling lactate dehydrogenase from Trichomonas vaginalis (TvLDH) is described as an example. The download and installation of the MODELLER software is also described. © 2016 by John Wiley & Sons, Inc.

Structure-based site-directed photo-crosslinking analyses of multimeric cell-adhesive interactions of voltage-gated sodium channel β subunits

The β1, β2, and β4 subunits of voltage-gated sodium channels reportedly function as cell adhesion molecules. The present crystallographic analysis of the β4 extracellular domain revealed an antiparallel arrangement of the β4 molecules in the crystal lattice. The interface between the two antiparallel β4 molecules is asymmetric, and results in a multimeric assembly. Structure-based mutagenesis and site-directed photo-crosslinking analyses of the β4-mediated cell-cell adhesion revealed that the interface between the antiparallel β4 molecules corresponds to that in the trans homophilic interaction for the multimeric assembly of β4 in cell-cell adhesion. This trans interaction mode is also employed in the β1-mediated cell-cell adhesion. Moreover, the β1 gene mutations associated with generalized epilepsy with febrile seizures plus (GEFS+) impaired the β1-mediated cell-cell adhesion, which should underlie the GEFS+ pathogenesis. Thus, the structural basis for the β-subunit-mediated cell-cell adhesion has been established.

Structural studies of the protein endostatin in fusion with BAX BH3 death domain, a hybrid that presents enhanced antitumoral activity

Endostatin (ES) is an antiangiogenic protein that exhibits antitumor activity in animal models. However, the activity observed in animals was not observed in human clinical trials. ES-BAX is a fusion protein composed of two functional domains: ES, which presents specificity and is internalized by activated endothelial cells and the proapoptotic BH3 domain of the protein BAX, a peptide inductor of cellular death when internalized. We have previously shown (Chura-Chambi et al., Cell Death Dis, 5, e1371, 2014) that ES-BAX presents improved antitumor activity in relation to wild-type ES. Secondary and tertiary structures of ES-BAX are similar to ES, as indicated by homology-modeling studies and molecular dynamics simulations. Tryptophan intrinsic fluorescence and circular dichroism spectroscopy corroborate these data. ¹⁵ N HSQC NMR indicates that ES-BAX is structured, but some ES residues have suffered chemical shift perturbations, suggesting that the BH3 peptide interacts with some parts of the ES protein. ES and ES-BAX present similar stability to thermal denaturation. The production of stable hybrid proteins can be a new approach to the development of therapeutic agents presenting specificity for tumoral endothelium and improved antitumor effect.

Three conserved MyD88-recruiting TLR residues exert different effects on the human TLR4 signaling pathway

Stimulation of Toll-like receptor (TLR) 4 leads to the activation of both MyD88-dependent and MyD88-independent pathways through the recruitment of adaptors TIRAP/MyD88 and TRIF/TRAM, respectively. However, the molecular basis of the TLR4 Toll/interleukin-1 receptor (TIR) domain in recruiting these downstream adaptors is still not entirely clear. Here, we identify three amino acid residues (714P in the BB loop, 696L in the αA helix and 721N in the αB sheet) conserved in all MyD88-recruited TLRs, but not the TLR3 TIR domain, as being critical for TLR4 responsiveness to LPS. These results were based on the substitution of each residue with a residue of the opposite type (hydrophilic/hydrophobic). However, the responsiveness of the TLR4 mutants to LPS was only partially decreased when each residue was replaced with a residue having the same hydrophilicity/hydrophobicity. This result is likely associated with an alteration in the BB-loop conformation of each TLR4 mutant and its ability to recruit the downstream adaptor TRAM. Thus, we identified three amino acids essential for TLR4 signaling, and their replacement with a residue of the same or opposite hydrophilicity/hydrophobicity greatly affected TLR4 signaling. This study furthers our understanding of the molecular mechanism by which the TLR4 TIR domain modulates TLR4 signaling and also provides new insight for the design of antisepsis therapy.

Comparative Modeling, Molecular Docking, and Revealing of Potential Binding Pockets of RASSF2; a Candidate Cancer Gene

RASSF2, potential tumor suppressor gene, acts as a KRAS-specific effectors protein and may promote apoptosis and cell cycle arrest. It stabilizes STK3/MST2 by protecting it from proteasomal degradation. RASSF2 plays a significant role against the inhibition of cancer. MODELLER (9v15) and online servers (I-Tasser, SwissModel, 3D-JigSaw, ModWeb) were utilized to generate 3D structures of the RASSF2 based on homology modeling. A comparison between models predicted by MODELLER (9v15) and Web servers had been checked through utilized evaluation tools. The most potent model for RASSF2 was analyzed and selected for molecular docking studies. The binding pockets were revealed for binding studies through Site Hound. AutoDock Vina and AutoDock4 were utilized for molecular docking, and the attempt of this experiment was to identify the ligands for RASSF2. The selected compounds may act as regulators and regulate the normal activity of RASSF2. It was also analyzed and observed that the selected compounds showed least binding energy and high-affinity binding in predicted top binding domain. The determination of protein function is based on accurate identification of binding sites in protein structures. The binding site is known, and it may allow the ligand type and protein function to be determined by performing in silico and experimental procedures. The detection, comparison, and analysis of binding pockets are pivotal to drug discovery. It proposed that predicted structure is reliable for the structural insights and functional studies. The predicted binding pockets may lead to further analysis (drug discovery), used against cancer study.

Development and Validation of a Small Single-domain Antibody That Effectively Inhibits Matrix Metalloproteinase 8

A detrimental role for matrix metalloproteinase 8 (MMP8) has been identified in several pathological conditions, e.g., lethal hepatitis and the systemic inflammatory response syndrome. Since matrix MMP8-deficient mice are protected in the above-mentioned diseases, specific MMP8 inhibitors could be of clinical value. However, targeting a specific matrix metalloproteinase remains challenging due to the strong structural homology of matrix metalloproteinases, which form a family of 25 members in mammals. Single-domain antibodies, called nanobodies, offer a range of possibilities toward therapy since they are easy to generate, express, produce, and modify, e.g., by linkage to nanobodies directed against other target molecules. Hence, we generated small MMP8-binding nanobodies, and established a proof-of-principle for developing nanobodies that inhibit matrix metalloproteinase activity. Also, we demonstrated for the first time the possibility of expressing nanobodies systemically by in vivo electroporation of the muscle and its relevance as a potential therapy in inflammatory diseases.

A frequent hypofunctional IRAK2 variant is associated with reduced spontaneous hepatitis C virus clearance

Patients carrying very rare loss-of-function mutations in interleukin-1 receptor-associated kinase 4 (IRAK4), a critical signaling mediator in Toll-like receptor signaling, are severely immunodeficient, highlighting the paramount role of IRAK kinases in innate immunity. We discovered a comparatively frequent coding variant of the enigmatic human IRAK2, L392V (rs3844283), which is found homozygously in ∼15% of Caucasians, to be associated with a reduced ability to induce interferon-alpha in primary human plasmacytoid dendritic cells in response to hepatitis C virus (HCV). Cytokine production in response to purified Toll-like receptor agonists was also impaired. Additionally, rs3844283 was epidemiologically associated with a chronic course of HCV infection in two independent HCV cohorts and emerged as an independent predictor of chronic HCV disease. Mechanistically, IRAK2 L392V showed intact binding to, but impaired ubiquitination of, tumor necrosis factor receptor-associated factor 6, a vital step in signal transduction.

CONCLUSION

Our study highlights IRAK2 and its genetic variants as critical factors and potentially novel biomarkers for human antiviral innate immunity.

Investigating the Conformational Structure and Potential Site Interactions of SOD Inhibitors on Ec-SOD in Marine Mud Crab Scylla serrata: A Molecular Modeling Approach

Superoxide dismutases (SODs) act as a first line of the enzymatic antioxidant defense system to control cellular superoxide anion toxicity. Previously, several inhibitors have been widely identified and catalogued for inhibition of SOD activity; however, still the information about the mechanism of interaction and points toward the inhibitor interactions in structures of SODs in general and in extracellular (Ec)-SOD in particular is still in naive. In the present research, we present an insight to elucidate the molecular basis of interactions of SOD inhibitors with Ec-SOD in mud crab Scylla serrata using molecular modeling and docking approaches. Different inhibitors of SOD such as hydrogen peroxide [Formula: see text], potassium cyanide, sodium dodecyl sulfate (SDS), [Formula: see text]-mercaptoethanol and dithiocarbamate were screened to understand the potential sites that may act as sites for cleavage or blocking in the protein. SOD-SDS and [Formula: see text] complex interactions indicate residues Pro72 and Asp102 of the predicted crab Ec-SOD as common targets. The GOLD result indicates that Pro72, Asp102 and Thr103 are commonly acting as the site of interaction in Ec-SOD of S. serrata with SOD inhibitors. For the first time, the results of this study provide an insight into the structural properties of Ec-SOD of S. serrata and define the possible involvements between the amino acids present in its active sites, i.e., in the regions from 70 to 84 and from 101 to 103 and different inhibitors.

In Silico Approach to Support that p-Nitrophenol Monooxygenase from Arthrobacter sp. Strain JS443 Catalyzes the Initial Two Sequential Monooxygenations

p-Nitrophenol (PNP), used primarily for manufacturing pesticides and dyes, has been recognized as a priority environmental pollutant. It is therefore important to reduce the input of this toxicant into the environment and to establish approaches for its removal from the contaminated sites. PNP monooxygenase, a novel enzyme from Gram-positive bacteria like Arthrobacter sp. and Bacillus sp., that comprises two components, a flavoprotein reductase and an oxygenase, catalyzes the initial two sequential monooxygenations to convert PNP to trihydroxybenzene. Accurate and reliable prediction of this enzyme-substrate interactions and binding affinity are of vital importance in understanding these catalytic mechanisms of the two sequential reactions. As crystal structure of the enzyme has not yet been published, we built a homology model for PNP monooxygenase using crystallized chlorophenol 4-monooxygenase from Burkholderia cepacia AC1100 (3HWC) as the template. The model was assessed for its reliability using PROCHECK, ERRAT and ProSA. Molecular docking of the physiological substrates, PNP and 4-nitrocatechol (4-NC), was carried out using Glide v5.7 implemented in Maestro v9.2, and the binding energies were calculated to substantiate the prediction. Docking complexes formed by molecular level interactions of PNP monooxygenase-PNP/4-NC without or with the cofactors, FAD and NADH, showed good correlation with the established experimental evidence that the two-component PNP monooxygenase catalyzes both the hydroxylation of PNP and the oxidative release of nitrite from 4-NC in B. sphaericus JS905. Furthermore, molecular dynamics simulations performed for docking complexes using Desmond v3.0 showed stable nature of the interactions as well.

Functional conservation of an ancestral Pellino protein in helminth species

The immune system of H. sapiens has innate signaling pathways that arose in ancestral species. This is exemplified by the discovery of the Toll-like receptor (TLR) pathway using free-living model organisms such as Drosophila melanogaster. The TLR pathway is ubiquitous and controls sensitivity to pathogen-associated molecular patterns (PAMPs) in eukaryotes. There is, however, a marked absence of this pathway from the plathyhelminthes, with the exception of the Pellino protein family, which is present in a number of species from this phylum. Helminth Pellino proteins are conserved having high similarity, both at the sequence and predicted structural protein level, with that of human Pellino proteins. Pellino from a model helminth, Schistosoma mansoni Pellino (SmPellino), was shown to bind and poly-ubiquitinate human IRAK-1, displaying E3 ligase activity consistent with its human counterparts. When transfected into human cells SmPellino is functional, interacting with signaling proteins and modulating mammalian signaling pathways. Strict conservation of a protein family in species lacking its niche signalling pathway is rare and provides a platform to examine the ancestral functions of Pellino proteins that may translate into novel mechanisms of immune regulation in humans.