Molecular modeling of formate dehydrogenase: the formation of the Michaelis complex

The formation of the reactive enzyme–substrate complex of formate dehydrogenase has been investigated by molecular dynamics techniques accounting for different conformational states of the enzyme. Simulations revealed that the transport of substrate to the active site through the substrate channel proceeds in the open conformation of enzyme due to the crucial role of the Arg284 residue acting as a vehicle. However, formate binding in the active site of the open conformation leads to the formation of a nonproductive enzyme–substrate complex. The productive Michaelis complex is formed only in the closed enzyme conformation after the substrate and coenzyme have bound, when required rigidity of the binding site and reactive formate orientation due to interactions with Arg284, Asn146, Ile122, and His332 residues is attained. Then, the high occupancy (up to 75%) of the reactive substrate–coenzyme conformation is reached, which was demonstrated by hybrid quantum mechanics/molecular mechanics simulations using various semiempirical Hamiltonians.     

The Interaction Between the Novel Intercalator Diethidium Cation and B-Form DNA: a Theoretical Study

Abstract

This research is an effort to further understand the physicochemical interaction between the novel drug molecule diethidium (2,7-diamino 9-[2,7 diamino 10-nN- phenanthridium] 10- nN- phenanthridium) and its biological receptor DNA. The ultimate goal is the elucidation of this novel class of drugs as potential pharmaceutical agents. Understanding the physico- chemical properties of this drug as well as the mechanism by which it interacts with DNA should ultimately allow the rational design of novel anti-cancer or anti-viral drugs.

A novel binding structure for the diethidium cation to B-form DNA is herein described. Molecular modeling on the complex formed between diethidium and a dodecamer of double-stranded B-form DNA, CGCGAATTCGCG, has shown that this complex is indeed fully capable of participating in the formation of a stable intercalation site. It was expected that diethidium would have a mechanism of intercalation significantly different from other classical intercalators because a) Its structure, that of two perpendicular planes, each known to have excellent intercalation properties, is novel b) The linker region length is zero c) The tilt between the two planes of the drug matches the geometry of the space available to this drug in the major groove.

We have studied the complex formed when diethidium enters the central site of the B-DNA dodecamer through the major groove. The complex forms several classes of intercalation structures, which are all stable and vary from “partially intercalate” to “fully intercalated”. Multiple minimizations show the drug to be very mobile within the intercalation site. Further, some structures show organization and concomitant stiffening of the DNA above the intercalation site, with a disorganization and disruption of the regular B-DNA structure immediately below the intercalation site. This particular phenomena may be expected to lead to significantly different physicochemical properties for the diethidium complex with respect to other known intercalators, because this sort of vectorial difference in structure above and below the site of intercalation is unknown in existing intercalators, as far as the authors are aware. In addition, we expect the mechanism of interaction between drug and DNA to be described by “direct ligand transfer”, wherein the drug is transferred from duplex DNA to duplex DNA without re-entering the solvent.¹

This work is the first instance known to the authors of a novel drug entity that was deduced solely by mathematical reasoning ² and described subsequently by computational methods. Evidence that diethidium should interact with its target site DNA differently from other known intercalators is strong.     

Impact of Fc N-glycan sialylation on IgG structure

ABSTRACT

Human IgG antibodies containing terminal alpha 2,6-linked sialic acid on their Fc N-glycans have been shown to reduce antibody-dependent cell-mediated cytotoxicity and possess anti-inflammatory properties. Although terminal sialylation on complex N-glycans can happen via either an alpha 2,3-linkage or an alpha 2,6-linkage, sialic acids on human serum IgG Fc are almost exclusively alpha 2,6-linked. Recombinant IgGs expressed in Chinese hamster ovary (CHO) cells, however, have sialic acids through alpha 2,3-linkages because of the lack of the alpha 2,6-sialyltransferase gene. The impact of different sialylation linkages to the structure of IgG has not been determined. In this work, we investigated the impact of different types of sialylation to the conformational stability of IgG through hydrogen/deuterium exchange (HDX) and limited proteolysis experiments. When human-derived and CHO-expressed IgG1 were analyzed by HDX, sialic acid-containing glycans were found to destabilize the CH2 domain in CHO-expressed IgG, but not human-derived IgG. When structural isomers of sialylated glycans were chromatographically resolved and identified in the limited proteolysis experiment, we found that only alpha 2,3-linked sialic acid on the 6-arm (the major sialylated glycans in CHO-expressed IgG1) destabilizes the CH2 domain, presumably because of the steric effect that decreases the glycan-CH2 domain interaction. The alpha 2,6-linked sialic acid on the 3-arm (the major sialylated glycan in human-derived IgG), and the alpha 2,3-linked sialic acid on the 3-arm, do not have this destabilizing effect.     

New insight into the binding mode of peptides at urotensin‐II receptor by Trp‐constrained analogues of P5U and urantide

Urotensin II (U‐II) is a disulfide bridged peptide hormone identified as the ligand of a G‐protein‐coupled receptor. Human U‐II (H‐Glu‐Thr‐Pro‐Asp‐c[Cys‐Phe‐Trp‐Lys‐Tyr‐Cys]‐Val‐OH) has been described as the most potent vasoconstrictor compound identified to date. We have recently identified both a superagonist of human U‐II termed P5U (H‐Asp‐c[Pen‐Phe‐Trp‐Lys‐Tyr‐Cys]‐Val‐OH) and the compound termed urantide (H‐Asp‐c[Pen‐Phe‐d ‐Trp‐Orn‐Tyr‐Cys]‐Val‐OH), which is the most potent UT receptor peptide antagonist described to date. In the present study, we have synthesized four analogues of P5U and urantide in which the Trp⁷ residue was replaced by the highly constrained l ‐Tpi and d ‐Tpi residues. The replacement of the Trp⁷ by Tpi led to active analogues. Solution NMR analysis allowed improving the knowledge on conformation–activity relationships previously reported on UT receptor ligands. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Mass spectrometry‐based methods to study protein architecture and dynamics

Mass spectrometry is now an indispensable tool in the armamentarium of molecular biophysics, where it is used for tasks ranging from protein sequencing and mapping of post‐translational modifications to studies of higher order structure, conformational dynamics, and interactions of proteins with small molecule ligands and other biopolymers. This mini‐review highlights several popular mass spectrometry‐based tools that are now commonly used for structural studies of proteins beyond their covalent structure with a particular emphasis on hydrogen exchange and direct electrospray ionization mass spectrometry.     

A mass weighted chemical elastic network model elucidates closed form domain motions in proteins

An elastic network model (ENM), usually Cα coarse‐grained one, has been widely used to study protein dynamics as an alternative to classical molecular dynamics simulation. This simple approach dramatically saves the computational cost, but sometimes fails to describe a feasible conformational change due to unrealistically excessive spring connections. To overcome this limitation, we propose a mass‐weighted chemical elastic network model (MWCENM) in which the total mass of each residue is assumed to be concentrated on the representative alpha carbon atom and various stiffness values are precisely assigned according to the types of chemical interactions. We test MWCENM on several well‐known proteins of which both closed and open conformations are available as well as three α‐helix rich proteins. Their normal mode analysis reveals that MWCENM not only generates more plausible conformational changes, especially for closed forms of proteins, but also preserves protein secondary structures thus distinguishing MWCENM from traditional ENMs. In addition, MWCENM also reduces computational burden by using a more sparse stiffness matrix.     

Design,synthesis and evaluation of antimicrobial activity of N-terminal modified Leucocin A analogues

Class IIa bacteriocins are potent antimicrobial peptides produced by lactic acid bacteria to destroy competing microorganisms. The N-terminal domain of these peptides consists of a conserved YGNGV sequence and a disulphide bond. The YGNGV motif is essential for activity, whereas, the two cysteines involved in the disulphide bond can be replaced with hydrophobic residues. The C-terminal region has variable sequences, and folds into a conserved amphipathic α-helical structure. To elucidate the structure–activity relationship in the N-terminal domain of these peptides, three analogues (1–3) of a class IIa bacteriocin, Leucocin A (LeuA), were designed and synthesized by replacing the N-terminal β-sheet residues of the native peptide with shorter β-turn motifs. Such replacement abolished the antibacterial activity in the analogues, however, analogue 1 was able to competitively inhibit the activity of native LeuA. Native LeuA (37-mer) was synthesized using native chemical ligation method in high yield. Solution conformation study using circular dichroism spectroscopy and molecular dynamics simulations suggested that the C-terminal region of analogue 1 adopts helical folding as found in LeuA, while the N-terminal region did not fold into β-sheet conformation. These structure–activity studies highlight the role of proper folding and complete sequence in the activity of class IIa bacteriocins.