Characterization of recombinant chloroperoxidase, and F103A and C29H/C79H/C87H mutants

Mechanistically and structurally chloroperoxidase (CPO) occupies a unique niche among heme containing enzymes. Chloroperoxidase catalyzes a broad range of reactions, such as oxidation of organic substrates, dismutation of hydrogen peroxide, and mono-oxygenation of organic molecules. To expand the synthetic utility of CPO and to appreciate the important interactions that lead to CPO’s exceptional properties, a site-directed mutagenesis study was undertaken. ^ Recombinant CPO and CPO mutants were heterologously expressed in Aspergillus niger. The overall protein structure was almost the same as that of wild type CPO, as determined by UV-vis, NMR and CD spectroscopies. Phenylalanine103, which was proposed to regulate substrate access to the active site by restricting the size of substrates and to control CPO’s enantioselectivity, was mutated to Ala. The ligand binding affinity and most importantly the catalytic activity of F103A was dramatically different from wild type CPO. The mutation essentially eliminated the chlorination and dismutation activities but enhanced, 4-10 fold, the epoxidation, peroxidation, and N-demethylation activities. As expected, the F103A mutant displayed dramatically improved epoxidation activity for larger, more branched styrene derivatives. Furthermore, F103A showed a distinctive enantioselectivity profile: losing enantioselectivity to styrene and cis-β-methylstyrene; having a different configuration preference on α-methylstyrene; showing higher enantioselectivites and conversion rates on larger, more branched substrates. Our results show that F103 acts as a switch box that controls the catalytic activity, substrate specificity, and product enantioselectivity of CPO. Given that no other mutant of CPO has displayed distinct properties, the results with F103A are dramatic. ^ The diverse catalytic activity of CPO has long been attributed to the presence of the proximal thiolate ligand. Surprisingly, a recent report on a C29H mutant suggested otherwise. A new CPO triple mutant C29H/C79H/C87H was prepared, in which all the cysteines were replaced by histidine to eliminate the possibility of cysteine coordinating to the heme. No active form protein was isolated, although, successful transformation and transcription was confirmed. The result suggests that Cys79 and Cys87 are critical to maintaining the structural scaffold of CPO. ^ In vitro biodegradation of nanotubes by CPO were examined by scanning electron microscope method, but little oxidation was observed. ^

Emerging trends at the interface of chemistry and biology: Applications to the design of human therapeutics

This article describes recent developments in the design and implementation of various strategies towards the development of novel therapeutics using first principles from biology and chemistry. Strategies for multi-target therapeutics and network analysis with a focus on cancer and HIV are discussed. Methods for gene and siRNA delivery are presented along with challenges and opportunities for siRNA therapeutics. Advances in protein design methodology and screening are described, with a focus on their application to the design of antibody based therapeutics. Future advances in this area relevant to vaccine design are also mentioned.

Structural analysis of human CCR2b and primate CCR2b by molecular modeling and molecular dynamics simulation

CCR2b, a chemokine receptor for MCP-1, -2, -3, -4, plays an important role in a variety of diseases involving infection, inflammation, and/or injury, as well as being a coreceptor for HIV-1 infection. Two models of human CCR2b (hCCR2b) were generated by h

Structural and functional characterization of the human CCR5 receptor in complex with HIV gp120 envelope glycoprotein and CD4 receptor by molecular modeling studies

The entry of human immunodeficiency virus (HIV) into cells depends on a sequential interaction of the gp120 envelope glycoprotein with the cellular receptors CD4 and members of the chemokine receptor family. The CC chemokine receptor CCR5 is such a receptor for several chemokines and a major coreceptor for the entry of R5 HIV type-1 (HIV-1) into cells. Although many studies focus on the interaction of CCR5 with HIV-1, the corresponding interaction sites in CCR5 and gp120 have not been matched. Here we used an approach combining protein structure modeling, docking and molecular dynamics simulation to build a series of structural models of the CCR5 in complexes with gp120 and CD4. Interactions such as hydrogen bonds, salt bridges and van der Waals contacts between CCR5 and gp120 were investigated. Three snapshots of CCR5-gp120-CD4 models revealed that the initial interactions of CCR5 with gp120 are involved in the negatively charged N-terminus (Nt) region of CCR5 and positively charged bridging sheet region of gp120. Further interactions occurred between extracellular loop2 (ECL2) of CCR5 and the base of V3 loop regions of gp120. These interactions may induce the conformational changes in gp120 and lead to the final entry of HIV into the cell. These results not only strongly support the two-step gp120-CCR5 binding mechanism, but also rationalize extensive biological data about the role of CCR5 in HIV-1 gp120 binding and entry, and may guide efforts to design novel inhibitors.

Searching for novel gene functions in yeast : identification of thousands of novel molecular interactions by protein-fragment complementation assay followed by automated gene function prediction and high-throughput lipidomics

La compréhension de processus biologiques complexes requiert des approches expérimentales et informatiques sophistiquées. Les récents progrès dans le domaine des stratégies génomiques fonctionnelles mettent dorénavant à notre disposition de puissants outils de collecte de données sur l’interconnectivité des gènes, des protéines et des petites molécules, dans le but d’étudier les principes organisationnels de leurs réseaux cellulaires. L’intégration de ces connaissances au sein d’un cadre de référence en biologie systémique permettrait la prédiction de nouvelles fonctions de gènes qui demeurent non caractérisées à ce jour. Afin de réaliser de telles prédictions à l’échelle génomique chez la levure Saccharomyces cerevisiae, nous avons développé une stratégie innovatrice qui combine le criblage interactomique à haut débit des interactions protéines-protéines, la prédiction de la fonction des gènes in silico ainsi que la validation de ces prédictions avec la lipidomique à haut débit. D’abord, nous avons exécuté un dépistage à grande échelle des interactions protéines-protéines à l’aide de la complémentation de fragments protéiques. Cette méthode a permis de déceler des interactions in vivo entre les protéines exprimées par leurs promoteurs naturels. De plus, aucun biais lié aux interactions des membranes n’a pu être mis en évidence avec cette méthode, comparativement aux autres techniques existantes qui décèlent les interactions protéines-protéines. Conséquemment, nous avons découvert plusieurs nouvelles interactions et nous avons augmenté la couverture d’un interactome d’homéostasie lipidique dont la compréhension demeure encore incomplète à ce jour. Par la suite, nous avons appliqué un algorithme d’apprentissage afin d’identifier huit gènes non caractérisés ayant un rôle potentiel dans le métabolisme des lipides. Finalement, nous avons étudié si ces gènes et un groupe de régulateurs transcriptionnels distincts, non préalablement impliqués avec les lipides, avaient un rôle dans l’homéostasie des lipides. Dans ce but, nous avons analysé les lipidomes des délétions mutantes de gènes sélectionnés. Afin d’examiner une grande quantité de souches, nous avons développé une plateforme à haut débit pour le criblage lipidomique à contenu élevé des bibliothèques de levures mutantes. Cette plateforme consiste en la spectrométrie de masse à haute resolution Orbitrap et en un cadre de traitement des données dédié et supportant le phénotypage des lipides de centaines de mutations de Saccharomyces cerevisiae. Les méthodes expérimentales en lipidomiques ont confirmé les prédictions fonctionnelles en démontrant certaines différences au sein des phénotypes métaboliques lipidiques des délétions mutantes ayant une absence des gènes YBR141C et YJR015W, connus pour leur implication dans le métabolisme des lipides. Une altération du phénotype lipidique a également été observé pour une délétion mutante du facteur de transcription KAR4 qui n’avait pas été auparavant lié au métabolisme lipidique. Tous ces résultats démontrent qu’un processus qui intègre l’acquisition de nouvelles interactions moléculaires, la prédiction informatique des fonctions des gènes et une plateforme lipidomique innovatrice à haut débit , constitue un ajout important aux méthodologies existantes en biologie systémique. Les développements en méthodologies génomiques fonctionnelles et en technologies lipidomiques fournissent donc de nouveaux moyens pour étudier les réseaux biologiques des eucaryotes supérieurs, incluant les mammifères. Par conséquent, le stratégie présenté ici détient un potentiel d’application au sein d’organismes plus complexes.

Inhibition of Mycobacterium tuberculosis tyrosine phosphatase PtpA by synthetic chalcones: Kinetics, molecular modeling, toxicity and effect on growth

Tuberculosis (TB) is a major cause of morbidity and mortality throughout the world, and it is estimated that one-third of the world`s population is infected with Mycobacterium tuberculosis. Among a series of tested compounds, we have recently identified five synthetic chalcones which inhibit the activity of M. tuberculosis protein tyrosine phosphatase A (PtpA), an enzyme associated with M. tuberculosis infectivity. Kinetic studies demonstrated that these compounds are reversible competitive inhibitors. In this work we also carried out the analysis of the molecular recognition of these inhibitors on their macromolecular target, PtpA, through molecular modeling. We observed that the predominant determinants responsible for the inhibitory activity of the chalcones are the positions of the two methoxyl groups at the A-ring, that establish hydrogen bonds with the amino acid residues Arg17, His49, and Thr12 in the active site of PtpA, and the substitution of the phenyl ring for a 2-naphthyl group as B-ring, that undergoes p stacking hydrophobic interaction with the Trp48 residue from PtpA. Interestingly, reduction of mycobacterial survival in human macrophages upon inhibitor treatment suggests their potential use as novel therapeutics. The biological activity, synthetic versatility, and low cost are clear advantages of this new class of potential tuberculostatic agents. (C) 2010 Elsevier Ltd. All rights reserved.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

Cytochrome c-promoted cardiolipin oxidation generates singlet molecular oxygen

The interaction of cytochrome c (cyt c) with cardiolipin (CL) induces protein conformational changes that favor peroxidase activity. This process has been correlated with CL oxidation and the induction of cell death. Here we report evidence demonstrating the generation of singlet molecular oxygen [O-2((1)Delta(g))] by a cyt c-CL complex in a model membrane containing CL. The formation of singlet oxygen was directly evidenced by luminescence measurements at 1270 nm and by chemical trapping experiments. Singlet oxygen generation required cyt c-CL binding and occurred at pH values higher than 6, consistent with lipid-protein interactions involving fully deprotonated CL species and positively charged residues in the protein. Moreover, singlet oxygen formation was specifically observed for tetralinoleoyl CL species and was not observed with monounsaturated and saturated CL species. Our results show that there are at least two mechanisms leading to singlet oxygen formation: one with fast kinetics involving the generation of singlet oxygen directly from CL hydroperoxide decomposition and the other involving CL oxidation. The contribution of the first mechanism was clearly evidenced by the detection of labeled singlet oxygen [O-18(2)((1)Delta(g))] from liposomes supplemented with 18-oxygen-labeled CL hydroperoxides. However quantitative analysis showed that singlet oxygen yield from CL hydroperoxides was minor (<5%) and that most of the singlet oxygen is formed from the second mechanism. Based on these data and previous findings we propose a mechanism of singlet oxygen generation through reactions involving peroxyl radicals (Russell mechanism) and excited triplet carbonyl intermediates (energy transfer mechanism).

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Revealing the elusive molecular biology of the vault RNA

Recently several novel and previously reported non-protein-coding RNAs (ncRNAs) have been identified to be upregulated upon Epstein-Barr virus (EBV) infection in human B-lymphocytes. A group of these significantly upregulated ncRNAs are called vault RNAs (vtRNAs). ,b Only about 5% of the total cellular vtRNAs are connected to the vault particle, the largest known ribonucleoprotein particle (RNP) in eukaryotic cells. However the function of this ncRNA family and moreover of the vault particle remains still rather unclear. Our previous findings suggest a link between EBV infection and vtRNA expression. Consequently we are interested which part of the viral genome is responsible for the upregulation and moreover which function the vtRNAs might possess during virus propagation. To address this question we have separately overexpressed specific EBV-encoded, latently expressed proteins in BL2-cells to determine the influence on the vault RNA levels. Thereby we identified one EBV-encoded protein, called Latent Membrane Protein 1 (LMP1), which significantly contributes to the vtRNA upregulation. We used LMP1 mutants to characterize the region of the protein and the responsible pathway for triggering the elevated vtRNA expression. Our results suggest that the NFkB- pathway might be involved in this process. To investigate a possible functional connection between the vtRNA and EBV infection, we have overexpressed vtRNA1-1 in BL41, a cell line usually not expressing this vault RNA. We show that overexpression of vtRNA1-1 leads to a better viral establishment and markedly protects cells from undergoing apoptosis. Knock-down of the major vault protein, the main component of the vault particle, had no effect on EBV infection and apoptosis resistance. Thus these results support the view that the observed phenotype is caused by the vtRNA rather than the vault particle.

Tracking individual membrane proteins and their biochemistry: The power of direct observation

The advent of single molecule fluorescence microscopy has allowed experimental molecular biophysics and biochemistry to transcend traditional ensemble measurements, where the behavior of individual proteins could not be precisely sampled. The recent explosion in popularity of new super-resolution and super-localization techniques coupled with technical advances in optical designs and fast highly sensitive cameras with single photon sensitivity and millisecond time resolution have made it possible to track key motions, reactions, and interactions of individual proteins with high temporal resolution and spatial resolution well beyond the diffraction limit. Within the purview of membrane proteins and ligand gated ion channels (LGICs), these outstanding advances in single molecule microscopy allow for the direct observation of discrete biochemical states and their fluctuation dynamics. Such observations are fundamentally important for understanding molecular-level mechanisms governing these systems. Examples reviewed here include the effects of allostery on the stoichiometry of ligand binding in the presence of fluorescent ligands; the observation of subdomain partitioning of membrane proteins due to microenvironment effects; and the use of single particle tracking experiments to elucidate characteristics of membrane protein diffusion and the direct measurement of thermodynamic properties, which govern the free energy landscape of protein dimerization. The review of such characteristic topics represents a snapshot of efforts to push the boundaries of fluorescence microscopy of membrane proteins to the absolute limit.

Molecular crowding effects of linear polymers in protein solutions

Measurement of protein-polymer second virial coefficients (B-AP) by sedimentation equilibrium studies of carbonic anhydrase and cytochrome c in the presence of dextrans (T10-T80) has revealed an inverse dependence of B-AP upon dextran molecular mass that conforms well with the behaviour predicted for the excluded-volume interaction between a spherical protein solute A and a random-flight representation of the polymeric cosolute P. That model of the protein-polymer interaction is also shown to provide a reasonable description of published gel chromatographic and equilibrium dialysis data on the effect of polymer molecular mass on BAP for human serum albumin in the presence of polyethylene glycols, a contrary finding from analysis of albumin solubility measurements being rejected on theoretical grounds. Inverse dependence upon polymer chainlength is also the predicted excluded-volume effect on the strength of several types of macromolecular equilibria-protein isomerization, protein dimerization, and 1 : 1 complex formation between dissimilar protein reactants. It is therefore concluded that published experimental observations of the reverse dependence, preferential reaction enhancement within DNA replication complexes by larger polyethylene glycols, must reflect the consequences of cosolute chemical interactions that outweigh those of thermodynamic nonideality arising from excluded-volume effects. (c) 2005 Elsevier B.V. All rights reserved.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Addressing accuracy and precision issues in iTRAQ quantitation

iTRAQ (isobaric tags for relative or absolute quantitation) is a mass spectrometry technology that allows quantitative comparison of protein abundance by measuring peak intensities of reporter ions released from iTRAQ-tagged peptides by fragmentation during MS/MS. However, current data analysis techniques for iTRAQ struggle to report reliable relative protein abundance estimates and suffer with problems of precision and accuracy. The precision of the data is affected by variance heterogeneity: low signal data have higher relative variability; however, low abundance peptides dominate data sets. Accuracy is compromised as ratios are compressed toward 1, leading to underestimation of the ratio. This study investigated both issues and proposed a methodology that combines the peptide measurements to give a robust protein estimate even when the data for the protein are sparse or at low intensity. Our data indicated that ratio compression arises from contamination during precursor ion selection, which occurs at a consistent proportion within an experiment and thus results in a linear relationship between expected and observed ratios. We proposed that a correction factor can be calculated from spiked proteins at known ratios. Then we demonstrated that variance heterogeneity is present in iTRAQ data sets irrespective of the analytical packages, LC-MS/MS instrumentation, and iTRAQ labeling kit (4-plex or 8-plex) used. We proposed using an additive-multiplicative error model for peak intensities in MS/MS quantitation and demonstrated that a variance-stabilizing normalization is able to address the error structure and stabilize the variance across the entire intensity range. The resulting uniform variance structure simplifies the downstream analysis. Heterogeneity of variance consistent with an additive-multiplicative model has been reported in other MS-based quantitation including fields outside of proteomics; consequently the variance-stabilizing normalization methodology has the potential to increase the capabilities of MS in quantitation across diverse areas of biology and chemistry.

Molecular phylogeny of the Chinese ranids inferred from nuclear and mitochondrial DNA sequences

The phylogeny of representative species of Chinese ranids was reconstructed using two nuclear (tyrosinase and rhodopsin) and two mitochondrial (12S rRNA, 16S rRNA) DNA fragments. Maximum parsimony, Bayesian, and maximum likelihood analyses were employed.-