In Silico Modeling and Conformational Mobility of DNA Duplexes

The review is focused on issues of transferability of the context-sensitive conformational characteristics of DNA estimated from crystallographic structural data on the DNA in aqueous solution. The state of the art in molecular dynamics of charged biopolymers in aqueous solution is covered. Elaboration of expedient force fields and algorithms of calculating long-range electrostatic interactions and solving combined equations of atomic motion have made it possible to generate stable nanosecond trajectories of thermal atomic motion of the biopolymer in aqueous solution in the presence of counterions and salt ions over reasonable time. Tools for analyzing the atomic statistical trajectories of DNA duplexes in aqueous solution to infer context-dependent conformational dynamic characteristics are discussed together with advances in simulating the mechanisms of global axial bend in DNA duplexes. These techniques allow one to consecutively analyze relationships between the contextual composition of the duplex and the basic modes of essential motions, their amplitude and extent of fluctuation. Development of satisfactory methods for estimating the free energy of biopolymer conformations in solution permits qualitative assessment of the conformational thermodynamic stability of biopolymers and their complexes.     

Competition by Estrogens for Catecholamine Receptor Binding In Vitro

Abstract: We have examined the ability of various steroids to compete for high-affinity binding of ³H-labeled ligands to catecholamine receptors in membranes prepared from rat cerebral cortex, striatum, and anterior pituitary. Ligands employed were: [³H]WB4101, [³H]prazosin, [³H]yohimbine, and [³H]clonidine (alpha-noradrenergic); [³H]dihydroalprenolol (beta-noradrenergic); [³H]spiperone and [³H]ADTN (dopaminergic). Only the 17β estrogens were effective and only binding of [³H]spiperone and [³H]ADTN in striatum and [³H]WB4101 and [³H]prazosin in cerebral cortex was reduced. Thus putative dopaminergic and alpha₁-noradrenergic sites alone appear to recognize estrogens. A slight competitive effect on [³H]spiperone binding to anterior pituitary membranes was also observed. Among the 17β estrogens tested, the most effective in all cases was the catechol estrogen 2-hydroxyestradiol (2-OHE₂). The ability of 2-OHE₂ (IC₅₀= 20–30 μM) to inhibit ligand binding to alpha₁ receptors was comparable to that of norepinephrine (IC₅₀= 10–20 μM), whereas for dopamine receptors in striatum and pituitary 2-OHE₂ was an order of magnitude less effective than dopamine (IC₃₀= 12 μM) in reducing binding of ³H ligands. Estradiol-17β and 2-hydroxyestrone were also able to inhibit binding, but the order of steroid potency was different for alpha₁ and dopaminergic receptors. Progesterone, testosterone, and corticosterone were without effect in all cases. These results show that there is specificity of steroid interactions with catecholamine receptors in the brain, both in terms of steroid structure and receptor type. The possible relevance of these interactions to neuroendocrine function is discussed.     

Crystal Structure of the Mineralocorticoid Receptor DNA Binding Domain in Complex with DNA

The steroid hormone receptors regulate important physiological functions such as reproduction, metabolism, immunity, and electrolyte balance. Mutations within steroid receptors result in endocrine disorders and can often drive cancer formation and progression. Despite the conserved three-dimensional structure shared among members of the steroid receptor family and their overlapping DNA binding preference, activation of individual steroid receptors drive unique effects on gene expression. Here, we present the first structure of the human mineralocorticoid receptor DNA binding domain, in complex with a canonical DNA response element. The overall structure is similar to the glucocorticoid receptor DNA binding domain, but small changes in the mode of DNA binding and lever arm conformation may begin to explain the differential effects on gene regulation by the mineralocorticoid and glucocorticoid receptors. In addition, we explore the structural effects of mineralocorticoid receptor DNA binding domain mutations found in type I pseudohypoaldosteronism and multiple types of cancer.     

In Silico Analysis Reveals Sequential Interactions and Protein Conformational Changes during the Binding of Chemokine CXCL-8 to Its Receptor CXCR1

Chemokine CXCL-8 plays a central role in human immune response by binding to and activate its cognate receptor CXCR1, a member of the G-protein coupled receptor (GPCR) family. The full-length structure of CXCR1 is modeled by combining the structures of previous NMR experiments with those from homology modeling. Molecular docking is performed to search favorable binding sites of monomeric and dimeric CXCL-8 with CXCR1 and a mutated form of it. The receptor-ligand complex is embedded into a lipid bilayer and used in multi ns molecular dynamics (MD) simulations. A multi-steps binding mode is proposed: (i) the N-loop of CXCL-8 initially binds to the N-terminal domain of receptor CXCR1 driven predominantly by electrostatic interactions; (ii) hydrophobic interactions allow the N-terminal Glu-Leu-Arg (ELR) motif of CXCL-8 to move closer to the extracellular loops of CXCR1; (iii) electrostatic interactions finally dominate the interaction between the N-terminal ELR motif of CXCL-8 and the EC-loops of CXCR1. Mutation of CXCR1 abrogates this mode of binding. The detailed binding process may help to facilitate the discovery of agonists and antagonists for rational drug design.     

Antagonistic TNF Receptor One-Specific Antibody (ATROSAB): Receptor Binding and In Vitro Bioactivity

Background

Selective inhibition of TNFR1 signaling holds the potential to greatly reduce the pro-inflammatory activity of TNF, while leaving TNFR2 untouched, thus allowing for cell survival and tissue homeostasis. ATROSAB is a humanized antagonistic anti-TNFR1 antibody developed for the treatment of inflammatory diseases.

Methodology/Principal Findings

The epitope of ATROSAB resides in the N-terminal region of TNFR1 covering parts of CRD1 and CRD2. By site-directed mutagenesis, we identified Arg68 and His69 of TNFR1 as important residues for ATROSAB binding. ATROSAB inhibited binding of ¹²⁵I-labeled TNF to HT1080 in the subnanomolar range. Furthermore, ATROSAB inhibited release of IL-6 and IL-8 from HeLa and HT1080 cells, respectively, induced by TNF or lymphotoxin alpha (LTα). Different from an agonistic antibody (Htr-9), which binds to a region close to the ATROSAB epitope but elicits strong TNFR1 activation, ATROSAB showed a negligible induction of IL-6 and IL-8 production over a broad concentration range. We further verified that ATROSAB, comprising mutations within the Fc region known to abrogate complement fixation and antibody-mediated cellular effector functions, indeed lacks binding activity for C1q, FcγRI (CD64), FcγRIIB (CD32b), and FcγRIII (CD16) disabling ADCC and CDC.

Conlusions/Significance

The data corroborate ATROSAB’s unique function as a TNFR1-selective antagonist efficiently blocking both TNF and LTα action. In agreement with recent studies of TNFR1 complex formation and activation, we suggest a model of the underlying mechanism of TNFR1 inhibition by ATROSAB.     

Structure of PP4397 Reveals the Molecular Basis for Different c-di-GMP Binding Modes by Pilz Domain Proteins

Cyclic diguanylate (c-di-GMP) is a global regulator that modulates pathogen virulence and biofilm formation in bacteria. Although a bioinformatic study revealed that PilZ domain proteins are the long-sought c-di-GMP binding proteins, the mechanism by which c-di-GMP regulates them is uncertain. Pseudomonas putida PP4397 is one such protein that contains YcgR-N and PilZ domains and the apo-PP4397 structure was solved earlier by the Joint Center for Structural Genomics. We determined the crystal structure of holo-PP4397 and found that two intercalated c-di-GMPs fit into the junction of its YcgR-N and PilZ domains. Moreover, c-di-GMP binding induces PP4397 to undergo a dimer-to-monomer transition. Interestingly, another PilZ domain protein, VCA0042, binds to a single molecule of c-di-GMP, and both its apo and holo forms are dimeric. Mutational studies and the additional crystal structure of holo-VCA0042 (L135R) showed that the Arg122 residue of PP4397 is crucial for the recognition of two molecules of c-di-GMP. Thus, PilZ domain proteins exhibit different c-di-GMP binding stoichiometry and quaternary structure, and these differences are expected to play a role in generating diverse forms of c-di-GMP-mediated regulation.     

Molecular Dynamics Simulation in Solvent of the Estrogen Receptor Protein DNA Binding Domain in Complex with a Non-Consensus Estrogen Response Element DNA Sequence

Abstract

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the estrogen receptor (ER) protein DNA binding domain (DBD) amino acids and DNA of a non-consensus estrogen response element (ERE) consisting of 29 nucleotide base pairs. This ERE nucleotide sequence occurs naturally upstream of the Xenopus laevis Vitelligenin AI gene. The ER DBD is encoded by three exons. Namely, exons 2 and 3 which encode the two zinc binding motifs and a sequence of exon 4 which encodes a predicted alpha helix. We generated a computer model of the ER DBD using atomic coordinates derived from the average of 30 nuclear magnetic resonance (NMR) spectroscopy coordinate sets. Amino acids on the carboxyl end of the ER DBD were disordered in both X-ray crystallography and NMR determinations and no coordinates were reported. This disordered region includes 10 amino acids of a predicted alpha helix encoded in exon 4 at the exon 3/4 splice junction. These amino acids are known to be important in DNA binding and are also believed to function as a nuclear translocation signal sequence for the ER protein. We generated a computer model of the predicted alpha helix consisting of the 10 amino acids encoded in exon 4 and attached this helix to the carboxyl end of the ER DBD at the exon 3/4 splice junction site. We docked the ER DBD model within the DNA major groove halfsites of the 29 base pair non-consensus ERE and flanking nucleotides. We constructed a solvated model with the ER DBD/ERE complex surrounded by a ten Angstrom water layer and conducted molecular dynamics simulations. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the ER DBD DNA recognition helix formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the ERE DNA right major groove halfsite. Amino acids of the ER DBD exon 4 encoded predicted alpha helix formed direct and water mediated H-bonds with base and backbone sites of their cognate codon-anticodon nucleotides within the minor grooves flanking the ERE DNA major groove halfsites. These interactions together induced bending of the DNA into the protein.     

Molecular Basis of Glucagon-like Peptide 1 Docking to Its Intact Receptor Studied with Carboxyl-terminal Photolabile Probes

The glucagon-like peptide 1 (GLP1) receptor is a member of Family B G protein-coupled receptors and represents an important drug target for type 2 diabetes. Despite recent solution of the structure of the amino-terminal domain of this receptor and that of several close family members, understanding of the molecular basis of natural ligand GLP1 binding to its intact receptor remains limited. The goal of this study was to explore spatial approximations between specific receptor residues within the carboxyl terminus of GLP1 and its receptor as normally docked. Therefore, we developed and characterized two high affinity, full-agonist photolabile GLP1 probes having sites for covalent attachment in positions 24 and 35. Both probes labeled the receptor specifically and saturably. Subsequent peptide mapping using chemical and proteinase cleavages of purified wild-type and mutant GLP1 receptor identified that the Arg¹³¹–Lys¹³⁶ segment at the juxtamembrane region of the receptor amino terminus contained the site of labeling for the position 24 probe, and the specific receptor residue labeled by this probe was identified as Glu¹³³ by radiochemical sequencing. Similarly, nearby residue Glu¹²⁵ within the same region of the receptor amino-terminal domain was identified as the site of labeling by the position 35 probe. These data represent the first direct demonstration of spatial approximation between GLP1 and its intact receptor as docked, providing two important constraints for the modeling of this interaction. This should expand our understanding of the molecular basis of natural agonist ligand binding to the GLP1 receptor and may be relevant to other family members.     

In Silico,In Vitro and In Vivo Analysis of Binding Affinity between N and C-Domains of Clostridium perfringens Alpha Toxin

Clostridium perfringens alpha toxin/phospholipase C (CP-PLC) is one of the most potent bacterial toxins known to cause soft tissue infections like gas gangrene in humans and animals. It is the first bacterial toxin demonstrated to be an enzyme with phospholipase, sphingomyelinase and lecithinase activities. The toxin is comprised of an enzymatic N-domain and a binding C-domain interconnected by a flexible linker. The N-domain alone is non-toxic to mammalian cells, but incubation with C-domain restores the toxicity, the mechanism of which is still not elucidated. The objectives of the current study were to investigate the formation of a stable N and C-domain complex, to determine possible interactions between the two domains in silico and to characterize the in vitro and in vivo correlates of the interaction. To establish the existence of a stable N and C-domain hybrid, in vitro pull down assay and dot-Far Western blotting assays were employed, where it was clearly revealed that the two domains bound to each other to form an intermediate. Using bioinformatics tools like MetaPPISP, PatchDock and FireDock, we predicted that the two domains may interact with each other through electrostatic interactions between at least six pairs of amino acids. This N and C-domains interacted with each other in 1:1 ratio and the hybrid lysed mouse erythrocytes in a slower kinetics when compared with wild type native Cp-PLC. BALB/c mice when challenged with N and C-domain hybrid demonstrated severe myonecrosis at the site of injection while no death was observed. Our results provide further insight into better understanding the mechanism for the toxicity of Cp-PLC N and C-domain mixture.     

Substitution of Glutamate Residue by Lysine in the Dimerization Domain Affects DNA Binding Ability of HapR by Inducing Structural Deformity in the DNA Binding Domain

HapR has been given the status of a high cell density master regulatory protein in Vibrio cholerae. Though many facts are known regarding its structural and functional aspects, much still can be learnt from natural variants of the wild type protein. This work aims at investigating the nature of functional inertness of a HapR natural variant harboring a substitution of a conserved glutamate residue at position 117 which participates in forming a salt bridge by lysine (HapR_V2G-E¹¹⁷K). Experimental evidence presented here reveals the inability of this variant to interact with various cognate promoters by in vitro gel shift assay. Furthermore, the elution profiles of HapR_V2G-E¹¹⁷K protein along with the wild type functional HapR_V2G in size-exclusion chromatography as well as circular dichroism spectra did not reflect any significant differences in its structure, thereby indicating the intactness of dimer in the variant protein. To gain further insight into the global shape of the proteins, small angle X-ray scattering analysis (SAXS) was performed. Intriguingly, increased radius of gyration of HapR_V2G-E¹¹⁷K of 27.5 Å in comparison to the wild type protein from SAXS data analyses implied a significant alteration in the global shape of the dimeric HapR_V2G-E¹¹⁷K protein. Structure reconstruction brought forth that the DNA binding domains were substantially “parted away” in this variant. Taken together, our data illustrates that substitution of the conserved glutamate residue by lysine in the dimerization domain induces separation of the two DNA binding domains from their native-like positioning without altering the dimeric status of HapR variant.     

Toward the Discovery of Vaccine Adjuvants: Coupling In Silico Screening and In Vitro Analysis of Antagonist Binding to Human and Mouse CCR4 Receptors

Background

Adjuvants enhance or modify an immune response that is made to an antigen. An antagonist of the chemokine CCR4 receptor can display adjuvant-like properties by diminishing the ability of CD4+CD25+ regulatory T cells (Tregs) to down-regulate immune responses.

Methodology

Here, we have used protein modelling to create a plausible chemokine receptor model with the aim of using virtual screening to identify potential small molecule chemokine antagonists. A combination of homology modelling and molecular docking was used to create a model of the CCR4 receptor in order to investigate potential lead compounds that display antagonistic properties. Three-dimensional structure-based virtual screening of the CCR4 receptor identified 116 small molecules that were calculated to have a high affinity for the receptor; these were tested experimentally for CCR4 antagonism. Fifteen of these small molecules were shown to inhibit specifically CCR4-mediated cell migration, including that of CCR4⁺ Tregs.

Significance

Our CCR4 antagonists act as adjuvants augmenting human T cell proliferation in an in vitro immune response model and compound SP50 increases T cell and antibody responses in vivo when combined with vaccine antigens of Mycobacterium tuberculosis and Plasmodium yoelii in mice.     

Conformational Properties of the TATA-Box Binding Sequence of DNA

Abstract

Nanosecond scale molecular dynamics simulations in water demonstrate that the DNA oligomer, GCGTATATAAAACGC, which contains a target site for the TATA-box binding protein (TBP), has an intrinsic preference to adopt an A-like conformation in the region of the TATA-box and undergoes bending related to that seen within in the TBP complex. This result is obtained from two independent simulations using different starting structures. In line with earlier suggestions of Guzikevich-Guerstein and Shakked, these simulations imply that an A-DNA conformation may be an important intermediate step in forming the strongly distorted DNA observed within the crystallographically determined complex with TBP. These results also support modeling studies by Lebrun et al. which suggest that the TBP binding mechanism can be broken down into a backbone transition to an A-like form coupled with a mechanical distortion which locally stretches and unwinds the DNA.