De novo prediction of the structures of M. tuberculosis membrane proteins

HigA2 (Rv2021c) Is a Transcriptional Regulator with Multiple Regulatory Targets in Mycobacterium tuberculosis

Toxin-antitoxin (TA) systems are the major mechanism for persister formation in Mycobacterium tuberculosis (Mtb). Previous studies found that HigBA2 (Rv2022c-Rv2021c), a predicted type II TA system of Mtb, could be activated for transcription in response to multiple stresses such as anti-tuberculosis drugs, nutrient starvation, endure hypoxia, acidic pH, etc. In this study, we determined the binding site of HigA2 (Rv2021c), which is located in the coding region of the upstream gene higB2 (Rv2022c), and the conserved recognition motif of HigA2 was characterized via oligonucleotide mutation. Eight binding sites of HigA2 were further found in the Mtb genome according to the conserved motif. RT-PCR showed that HigA2 can regulate the transcription level of all eight of these genes and three adjacent downstream genes. DNA pull-down experiments showed that twelve functional regulators sense external regulatory signals and may regulate the transcription of the HigBA2 system. Of these, Rv0903c, Rv0744c, Rv0474, Rv3124, Rv2603c, and Rv3583c may be involved in the regulation of external stress signals. In general, we identified the downstream target genes and possible upstream regulatory genes of HigA2, which paved the way for the illustration of the persistence establishment mechanism in Mtb.

Structural analysis of a peptide fragment of transmembrane transporter protein bilitranslocase

Using a combination of genomic and post-genomic approaches is rapidly altering the number of identified human influx carriers. A transmembrane protein bilitranslocase (TCDB 2.A.65) has long attracted attention because of its function as an organic anion carrier. It has also been identified as a potential membrane transporter for cellular uptake of several drugs and due to its implication in drug uptake, it is extremely important to advance the knowledge about its structure. However, at present, only the primary structure of bilitranslocase is known. In our work, transmembrane subunits of bilitranslocase were predicted by a previously developed chemometrics model and the stability of these polypeptide chains were studied by molecular dynamics (MD) simulation. Furthermore, sodium dodecyl sulfate (SDS) micelles were used as a model of cell membrane and herein we present a high-resolution 3D structure of an 18 amino acid residues long peptide corresponding to the third transmembrane part of bilitranslocase obtained by use of multidimensional NMR spectroscopy. It has been experimentally confirmed that one of the transmembrane segments of bilitranslocase has alpha helical structure with hydrophilic amino acid residues oriented towards one side, thus capable of forming a channel in the membrane.

Structural modelling and dynamics of proteins for insights into drug interactions

Proteins are the workhorses of biomolecules and their function is affected by their structure and their structural rearrangements during ligand entry, ligand binding and protein-protein interactions. Hence, the knowledge of protein structure and, importantly, the dynamic behaviour of the structure are critical for understanding how the protein performs its function. The predictions of the structure and the dynamic behaviour can be performed by combinations of structure modelling and molecular dynamics simulations. The simulations also need to be sensitive to the constraints of the environment in which the protein resides. Standard computational methods now exist in this field to support the experimental effort of solving protein structures. This review presents a comprehensive overview of the basis of the calculations and the well-established computational methods used to generate and understand protein structure and function and the study of their dynamic behaviour with the reference to lung-related targets.

FoldGPCR: structure prediction protocol for the transmembrane domain of G protein-coupled receptors from class A

Building reliable structural models of G protein-coupled receptors (GPCRs) is a difficult task because of the paucity of suitable templates, low sequence identity, and the wide variety of ligand specificities within the superfamily. Template-based modeling is known to be the most successful method for protein structure prediction. However, refinement of homology models within 1-3 A C alpha RMSD of the native structure remains a major challenge. Here, we address this problem by developing a novel protocol (foldGPCR) for modeling the transmembrane (TM) region of GPCRs in complex with a ligand, aimed to accurately model the structural divergence between the template and target in the TM helices. The protocol is based on predicted conserved inter-residue contacts between the template and target, and exploits an all-atom implicit membrane force field. The placement of the ligand in the binding pocket is guided by biochemical data. The foldGPCR protocol is implemented by a stepwise hierarchical approach, in which the TM helical bundle and the ligand are assembled by simulated annealing trials in the first step, and the receptor-ligand complex is refined with replica exchange sampling in the second step. The protocol is applied to model the human beta(2)-adrenergic receptor (beta(2)AR) bound to carazolol, using contacts derived from the template structure of bovine rhodopsin. Comparison with the X-ray crystal structure of the beta(2)AR shows that our protocol is particularly successful in accurately capturing helix backbone irregularities and helix-helix packing interactions that distinguish rhodopsin from beta(2)AR.

Identification of murine T-cell epitopes on low-molecular-mass secretory proteins (CFP11, CFP17, and TB18.5) of Mycobacterium tuberculosis

The low-molecular-mass secretory proteins of Mycobacterium tuberculosis have been shown to be major T-cell antigens during infection with the pathogenic bacterium. In this study, we determined murine T-cell epitopes on three low-molecular-mass proteins, CFP11 (Rv2433c), CFP17 (Rv1827), and TB18.5 (Rv0164) using DNA immunization of inbred mice. We analyzed interferon-gamma production from immune splenocytes in response to overlapping peptides covering these proteins. We identified two CD8+ T-cell epitopes on CFP11 and CFP17, one in BALB/c mice and the other in C57BL/6 mice, respectively. On TB18.5, we identified a CD8+ T-cell epitope in BALB/c mice and a CD4+ T-cell epitope in C57BL/6 mice. With the aid of computer algorithms, we could identify the minimal CD8+ T-cell epitopes. These T-cell epitopes are feasible for analysis of the role of antigen-specific T cells during M. tuberculosis infection.

Contribution of charged and polar residues for the formation of the E1-E2 heterodimer from Hepatitis C Virus

The transmembrane domains of the envelope glycoprotein E1 and E2 have crucial multifunctional roles in the biogenesis of hepatitis C virus. We have performed molecular dynamics simulations to investigate a structural model of the transmembrane segments of the E1-E2 heterodimer. The simulations support the key role of the Lys370-Asp728 ion pair for mediating the E1-E2 heterodimerization. In comparison to these two residues, the simulation results also reveal the differential effect of the conserved Arg730 residue that has been observed in experimental studies. Furthermore, we discovered the formation of inter-helical hydrogen bonds via Asn367 that stabilize dimer formation. Simulations of single and double mutants further demonstrate the importance of the ion-pair and polar interactions between the interacting helix monomers. The conformation of the E1 fragment in the simulation of the E1-E2 heterodimer is in close agreement with an NMR structure of the E1 transmembrane segment. The proposed model of the E1-E2 heterodimer supports the postulated cooperative insertion of both helices by the translocon complex into the bilayer.

The energy landscape for the interaction of the family 1 carbohydrate-binding module and the cellulose surface is altered by hydrolyzed glycosidic bonds

A multiscale simulation model is used to construct potential and free energy surfaces for the carbohydrate-binding module [CBM] from an industrially important cellulase, Trichoderma reesei cellobiohydrolase I, on the hydrophobic face of a coarse-grained cellulose Ibeta polymorph. We predict from computation that the CBM alone exhibits regions of stability on the hydrophobic face of cellulose every 5 and 10 A, corresponding to a glucose unit and a cellobiose unit, respectively. In addition, we predict a new role for the CBM: specifically, that in the presence of hydrolyzed cellulose chain ends, the CBM exerts a thermodynamic driving force to translate away from the free cellulose chain ends. This suggests that the CBM is not only required for binding to cellulose, as has been known for two decades, but also that it has evolved to both assist the enzyme in recognizing a cellulose chain end and exert a driving force on the enzyme during processive hydrolysis of cellulose.