Homology modeling with low sequence identity

Three-dimensional models of antigens with serodiagnostic potential for leprosy: An in silico study

BACKGROUND

Leprosy is a disease caused by Mycobacterium leprae (M. leprae), an intracellular pathogen that has tropism and affects skin and nervous system cells. The disease has two forms of presentation: Paucibacillary and multibacillary, with different clinical and immunological manifestations. Unlike what occurs in the multibacillary form , the diagnostic tests for the paucibacillary form are nonspecific and not very sensitive, allowing the existence of infected individuals without treatment, which contributes to the spread of the pathogen in the population. To mitigate this contamination, more sensitive diagnostic tests capable of detecting paucibacillary patients are needed.

AIM

To predict the three-dimensional structure models of M. leprae antigens with serodiagnostic potential for leprosy.

METHODS

In this in silico study, satisfactory templates were selected in the Protein Data Bank (PDB) using Basic Local Alignment Search Tool to predict the structural templates of ML2038, ML0286, ML0050, and 85B antigens by comparative modeling. The templates were selected according to general criteria such as sequence identity, coverage, X-ray resolution, Global Model Quality Estimate value and phylogenetic relationship; Clustal X 2.1 software was used in this analysis. Molecular modeling was completed using the software Modeller 9v13. Visualization of the models was made using ViewerLite 4.2 and PyMol software, and analysis of the quality of the predicted models was performed using the QMEAN score and Z-score. Finally, the three-dimensional moels were validated using the MolProbity and Verify 3D platforms.

RESULTS

The three-dimensional structure models of ML2038, ML0286, ML0050, and 85B antigens of M. leprae were predicted using the templates PDB: 3UOI (90.51% identity), PDB: 3EKL (87.46% identity), PDB: 3FAV (40.00% identity), and PDB: 1F0N (85.21% identity), respectively. The QMEAN and Z-score values indicated the good quality of the structure models. These data refer to the monomeric units of antigens, since some of these antigens have quaternary structure. The validation of the models was performed with the final three-dimensional structure - monomer (ML0050 and 85B antigens) and quaternary structures (ML2038 and ML0286). The majority of amino acid residues were observed in favorable and allowed regions in the Ramachandran plot, indicating correct positioning of the side chain and absence of steric impediment. The MolProbity score value and Verify 3D results of all models indicated a satisfactory prediction.

CONCLUSION

The polarized immune response against M. leprae creates a problem in leprosy detection. The selection of immunodominant epitopes is essential for the development of more sensitive serodiagnostic tests, for this it is important to know the three-dimensional structure of the antigens, which can be predicted with bioinformatics tools.

Solution structure and dynamics of the mitochondrial-targeted GTPase-activating protein (GAP) VopE by an integrated NMR/SAXS approach

The bacterial pathogen Vibrio cholerae use a type III secretion system to inject effector proteins into a host cell. Recently, a putative Toxic GTPase Activating Protein (ToxGAP) called Vibrio outer protein E (VopE) was identified as a T3SS substrate and virulence factor that affected host mitochondrial dynamics and immune response. However, biophysical and structural characterization has been absent. Here, we describe solution NMR structure of the putative GTPase-activating protein (GAP) domain (73-204) of VopE. Using size exclusion chromatography coupled with small-angle x-ray scattering and residual dipolar coupling data, we restrained the MD process to efficiently determine the overall fold and improve the quality of the output calculated structures. Comparing the structure of VopE with other ToxGAP's revealed a similar overall fold with several features unique to VopE. Specifically, the "Bulge 1," α1 helix, and noteworthy "backside linker" elements on the N-terminus are dissimilar to the other ToxGAP's. By using NMR relaxation dispersion experiments, we demonstrate that these regions undergo motions on a > 6 s-1 timescale. Based on the disposition of these mobile regions relative to the putative catalytic arginine residue, we hypothesize that the protein may undergo structural changes to bind cognate GTPases.

Homology Modeling and Molecular Dynamics Simulations of Trypanosoma cruzi Phosphodiesterase b1

Cyclic nucleotide phosphodiesterases have been implicated in the proliferation, differentiation and osmotic regulation of trypanosomatids; in some trypanosomatid species, they have been validated as molecular targets for the development of new therapeutic agents. Because the experimental structure of Trypanosoma cruzi PDEb1 (TcrPDEb1) has not been solved so far, an homology model of the target was created using the structure of Trypanosoma brucei PDEb1 (TbrPDEb1) as a template. The model was refined by extensive enhanced sampling molecular dynamics simulations, and representative snapshots were extracted from the trajectory by combined clustering analysis. This structural ensemble was used to develop a structure-based docking model of the target. The docking accuracy of the model was validated by redocking and cross-docking experiments using all available crystal structures of TbrPDEb1, whereas the scoring accuracy was validated through a retrospective screen, using a carefully curated dataset of compounds assayed against TbrPDEb1 and/or TcrPDEb1. Considering the results from in silico validations, the model may be applied in prospective virtual screening campaigns to identify novel hits, as well as to guide the rational design of potent and selective inhibitors targeting this enzyme.

Trop-2 cleavage by ADAM10 is an activator switch for cancer growth and metastasis

Trop-2 is a transmembrane signal transducer that can induce cancer growth. Using antibody targeting and N-terminal Edman degradation, we show here that Trop-2 undergoes cleavage in the first thyroglobulin domain loop of its extracellular region, between residues R87 and T88. Molecular modeling indicated that this cleavage induces a profound rearrangement of the Trop-2 structure, which suggested a deep impact on its biological function. No Trop-2 cleavage was detected in normal human tissues, whereas most tumors showed Trop-2 cleavage, including skin, ovary, colon, and breast cancers. Coimmunoprecipitation and mass spectrometry analysis revealed that ADAM10 physically interacts with Trop-2. Immunofluorescence/confocal time-lapse microscopy revealed that the two molecules broadly colocalize at the cell membrane. We show that ADAM10 inhibitors, siRNAs and shRNAs abolish the processing of Trop-2, which indicates that ADAM10 is an effector protease. Proteolysis of Trop-2 at R87-T88 triggered cancer cell growth both in vitro and in vivo. A corresponding role was shown for metastatic spreading of colon cancer, as the R87A-T88A Trop-2 mutant abolished xenotransplant metastatic dissemination. Activatory proteolysis of Trop-2 was recapitulated in primary human breast cancers. Together with the prognostic impact of Trop-2 and ADAM10 on cancers of the skin, ovary, colon, lung, and pancreas, these data indicate a driving role of this activatory cleavage of Trop-2 on malignant progression of tumors.

G protein-coupled estrogen receptor-1: homology modeling approaches and application in screening new GPER-1 modulators

G protein-coupled receptors (GPCRs) belong to the largest family of protein targets comprising over 800 members in which at least 500 members are the therapeutic targets. Among the GPCRs, G protein-coupled estrogen receptor-1 (GPER-1) has shown to have the ability in estrogen signaling. As GPER-1 plays a critical role in several physiological responses, GPER-1 has been considered as a potential therapeutic target to treat estrogen-based cancers and other non-communicable diseases. However, the progress in the understanding of GPER-1 structure and function is relatively slow due to the availability of a only a few selective GPER-1 modulators. As with many GPCRs, the X-ray crystal structure of GPER-1 is yet to be resolved and thus has led the researchers to search for new GPER-1 modulators using homology models of GPER-1. In this review, we aim to summarize various approaches used in the generation of GPER-1 homology model and their applications that have resulted in new GPER-1 ligands.

Structural and energetic basis for the inhibitory selectivity of both catalytic domains of dimeric HDAC6

HDAC6 is a protein involved in cancer, neurodegenerative disease and inflammatory disorders. To date, the full three-dimensional (3D) structure of human HDAC6 has not been elucidated; however, there are some experimental 3D structural homologs to HDAC6 that can be used as templates. In this work, we utilized molecular modeling procedures to model both of the catalytic domains of HDAC6 connected by the linker region where DMB region is placed. Once the 3D structure of human HDAC6 was obtained, it was structurally evaluated and submitted to docking and molecular dynamic (MD) simulations along with Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) method to explore the stability and the binding free energy properties of the HDAC6-ligand complexes. In addition, its structural and energetic behavior was explored with each one of the catalytic domains in the molecular recognition of six selective HDAC6 inhibitors, HPOB, CAY10603, Nexturastat, Rocilinostat, Tubacin and Tubastatin A for DD2, and with the so-called 9-peptide which is DD1-HDAC6 selective substrate. The use of the whole system (DD1-DMB-DD2) showed a tendency toward the ligand affinity of DD2, CAY10603> Tubacin > Rocilinostat > Nexturastat > HPOB > Tubastatin > 9-peptide, which is in line with experimental reports. However, 9-peptide showed a higher affinity for DD1, which agrees with experimental reports elsewhere. Principal component analysis provided important information about the structural changes linked to the molecular recognition process, whereas per-residue decomposition analysis revealed the energetic contribution of the key residues in the molecular binding and structural characteristics that could assist in drug design.

The impact of Thr91 mutation on c-Src resistance to UM-164: molecular dynamics study revealed a new opportunity for drug design

The emergence of a drug resistant non-receptor tyrosine kinase (c-Src) in triple-negative breast cancer (TNBC) remains a prime concern in relation to the burden of TNBC among people living with breast cancer and drug development. Thr91 mutation was found to induce a complete loss of protein conformation required for drug fitness. Herein, we provide the first account of the molecular impact of the Thr91 mutation on c-Src resistance to experimental drug UM-164 using various computational approaches, namely molecular dynamics simulation, principal component analysis (PCA), dynamic cross-correlation matrices (DCCM) analysis, hydrogen bond occupancy, thermodynamics calculation, ligand-residue interaction and residue interaction networks (RINs). Findings from this study revealed that Thr91 mutation leads to a steric conflict between UM-164 and the side chain of methionine (Met91); this mutation distorts the UM-164 optimum orientation on the conformational space of mutant c-Src compared to the wild-type; decreases hydrogen bond formation between the residues in the mutant protein structure; decreases the UM-164 binding energy in the mutant by -13.416 kcal mol^-1; reduces the residue correlation in the mutant protein structure; induces a change in the overall protein structure conformation from an inactive to active conformation; and distorts the ligand atomic interaction network and the residue interaction network. This report provides important insights that will assist in the further design of novel dual kinase inhibitors to minimise the chances of drug resistance in triple negative breast cancer.

Characterization and molecular modeling of polygalacturonase isoforms from Saccharomyces cerevisiae

Earlier, low-temperature-active polygalacturonase isoforms from Saccharomyces cerevisiae PVK4 were isolated and purified. Substrate specificity of polygalacturonase isoforms indicated high affinity for pectins and very low enzyme activity towards non-pectic polysaccharides. To characterize the polygalacturonase isoforms, biochemical, spectral, and in silico approaches were used. The apparent K_m and V_max values for hydrolysis of pectin and galacturonic acid were 0.31 mg/ml and 3.15 mmol min/mg, respectively. Interestingly, the polygalacturonase isoforms were found to be more stable at optimal pH and temperature of 4.5 and 40 °C, respectively. These isoforms were reacted with different metal ions; Cd²⁺ and Ni²⁺ severely inhibited the enzyme activity, while Mg²⁺, Zn²⁺, Cd²⁺, Fe²⁺ Cu²⁺, and Ni²⁺ inhibited to a lesser extent, which clearly demonstrated that variations in enzyme activity were due to their differential binding affinity of metal ions. Furthermore, decrease in the viscosity of polygalacturonic acid and citrus pectin by these isoforms was approximately four and six times higher than the rate of release of reducing sugars. This indicates that polygalacturonase isoforms have an endo-mode of action. In addition to the above, thermostability of purified polygalacturonase isoforms was studied by circular dichroism and fluorescence spectroscopy. Circular dichroism showed 18% alpha helix and 57% beta sheets at pH 5, while at pH 7, 8, and 9 there was an increase of random coil. Fluorescence studies revealed small conformational changes, which were observed at 30-50 °C, while unfolding transition region was noticed between 60 and 70 °C. The purified enzyme fractions were analyzed by MALDI-TOF MS. Finally, 3D model structures for isoenzymes of polygalacturonase of S. cerevisiae were generated and validated as good quality models, which are also suitable for molecular interaction studies.

A Numbering System for MFS Transporter Proteins

The Major Facilitator Superfamily (MFS) is one of the largest classes of secondary active transporters and is widely expressed in many domains of life. It is characterized by a common 12-transmembrane helix motif that allows the selective transport of a vast range of diverse substrates across the membrane. MFS transporters play a central role in many physiological processes and are increasingly recognized as potential drug targets. Despite intensive efforts, there are still only a handful of crystal structures and therefore homology modeling is likely to be a necessary process for providing models to interpret experiments for many years to come. However, the diversity of sequences and the multiple conformational states these proteins can exist in makes the process significantly more complicated, especially for sequences for which there is very little sequence identity to known templates. Inspired by the approach adopted many years ago for GPCRs, we have analyzed the large number of MFS sequences now available alongside the current structural information to propose a series of conserved contact points that can provide additional guidance for the homology modeling process. To enable cross-comparison across MFS models we also present a numbering scheme that can be used to provide a point of reference within each of the 12 transmembrane regions.

Macromolecular Modelling and Docking Simulations for the Discovery of Selective GPER Ligands

Estrogens influence multiple physiological processes and are implicated in many diseases as well. Cellular responses to estrogens are mainly mediated by the estrogen receptors (ER)α and ERβ, which act as ligand-activated transcription factors. Recently, a member of the G protein-coupled receptor (GPCR) superfamily, namely GPER/GPR30, has been identified as a further mediator of estrogen signalling in different pathophysiological conditions, including cancer. Today, computational methods are commonly used in all areas of health science research. Among these methods, virtual ligand screening has become an established technique for hit discovery and optimization. The absence of an established three-dimensional structure of GPER promoted studies of structure-based drug design in order to build reliable molecular models of this receptor. Here, we discuss the results obtained through the structure-based virtual ligand screening for GPER, which allowed the identification and synthesis of different selective agonist and antagonist moieties. These compounds led significant advances in our understanding of the GPER function at the cellular, tissue, and organismal levels. In particular, selective GPER ligands were critical toward the evaluation of the role elicited by this receptor in several pathophysiological conditions, including cancer. Considering that structure-based approaches are fundamental in drug discovery, future research breakthroughs with the aid of computer-aided molecular design and chemo-bioinformatics could generate a new class of drugs that, acting through GPER, would be useful in a variety of diseases as well as in innovative anticancer strategies.

Progress and Current Challenges in Modeling Large RNAs

Recent breakthroughs in next-generation sequencing technologies have led to the discovery of several classes of non-coding RNAs (ncRNAs). It is now apparent that RNA molecules are not only just carriers of genetic information but also key players in many cellular processes. While there has been a rapid increase in the number of ncRNA sequences deposited in various databases over the past decade, the biological functions of these ncRNAs are largely not well understood. Similar to proteins, RNA molecules carry out a function by forming specific three-dimensional structures. Understanding the function of a particular RNA therefore requires a detailed knowledge of its structure. However, determining experimental structures of RNA is extremely challenging. In fact, RNA-only structures represent just 1% of the total structures deposited in the PDB. Thus, computational methods that predict three-dimensional RNA structures are in high demand. Computational models can provide valuable insights into structure-function relationships in ncRNAs and can aid in the development of functional hypotheses and experimental designs. In recent years, a set of diverse RNA structure prediction tools have become available, which differ in computational time, input data and accuracy. This review discusses the recent progress and challenges in RNA structure prediction methods.

In silico analysis of missense mutations in LPAR6 reveals abnormal phospholipid signaling pathway leading to hypotrichosis

Autosomal recessive hypotrichosis is a rare genetic irreversible hair loss disorder characterized by sparse scalp hair, sparse to absent eyebrows and eyelashes, and sparse axillary and body hair. The study, presented here, established genetic linkage in four families showing similar phenotypes to lysophosphatidic acid receptor 6 (LPAR6) gene on chromosome 13q14.11-q21.32. Subsequently, sequence analysis of the gene revealed two previously reported missense mutations including p.D63V in affected members of one and p.I188F in three other families. Molecular modeling and docking analysis was performed to investigate binding of a ligand oleoyl-L-alpha-lysophosphatidic acid (LPA) to modeled protein structures of normal and mutated (D63V, G146R, I188F, N248Y, S3T, L277P) LPAR6 receptors. The mutant receptors showed a complete shift in orientation of LPA at the binding site. In addition, hydropathy analysis revealed a significant change in the membrane spanning topology of LPAR6 helical segments. The present study further substantiated involvement of LPAR6-LPA signaling in the pathogenesis of hypotrichosis/woolly hair and provided additional insight into the molecular mechanism of hair development.

Mutually exclusive binding of APPL(PH) to BAR domain and Reptin regulates β-catenin dependent transcriptional events

Reptin functions in a wide range of biological processes including chromatin remodelling, nucleolar organization and transcriptional regulation of WNT signalling. As β-catenin dependent transcriptional repression and activation events involve binding of Reptin and histone deacetylase 1 to APPL endocytic proteins, this complex has become an important target to identify molecules governing endocytic processes and WNT signalling. Here, we describe the structural basis of APPL binding to Reptin to explore their mode of binding in context with APPL1/APPL2 dimerization. There is an evidence that both PH and BAR domains of APPL proteins exhibit alternately conserved regions involved in hetero-dimerization process and our in-silico data also corroborate this fact. Moreover, APPL2(PH) domain binds to the BAR domain region encompassing a nuclear localization signal. We conclude that APPL(PH) binding to BAR domain and Reptin is mutually exclusive which regulates the nucleocytoplasmic shuttling of Reptin. Furthermore, Reptin is unable to bind with membrane-associated APPL proteins. These observations were further expanded by experimental approaches where we identified a novel point mutation D316N lying in the APPL1(PH) domain which resulted in a significantly reduced binding with Reptin. By luciferase assays, we observed that overexpression of APPL1(D316N) and APPL1(WT) stimulated β-catenin/TCF dependent transcriptional activity in a similar manner which suggested that binding of Reptin to APPL1 is not necessary for β-catenin dependent target gene expression. Overall, our data attempt to highlight a comparative role of APPL proteins in controlling β-catenin dependent transcription mechanism which may improve our understanding of gene regulation.

Analysis of the REJ Module of Polycystin-1 Using Molecular Modeling and Force-Spectroscopy Techniques

Polycystin-1 is a large transmembrane protein, which, when mutated, causes autosomal dominant polycystic kidney disease, one of the most common life-threatening genetic diseases that is a leading cause of kidney failure. The REJ (receptor for egg lelly) module is a major component of PC1 ectodomain that extends to about 1000 amino acids. Many missense disease-causing mutations map to this module; however, very little is known about the structure or function of this region. We used a combination of homology molecular modeling, protein engineering, steered molecular dynamics (SMD) simulations, and single-molecule force spectroscopy (SMFS) to analyze the conformation and mechanical stability of the first ~420 amino acids of REJ. Homology molecular modeling analysis revealed that this region may contain structural elements that have an FNIII-like structure, which we named REJd1, REJd2, REJd3, and REJd4. We found that REJd1 has a higher mechanical stability than REJd2 (~190 pN and 60 pN, resp.). Our data suggest that the putative domains REJd3 and REJd4 likely do not form mechanically stable folds. Our experimental approach opens a new way to systematically study the effects of disease-causing mutations on the structure and mechanical properties of the REJ module of PC1.

In silico analysis of the amido phosphoribosyltransferase inhibition by PY873, PY899 and a derivative of isophthalic acid

Selectively decreasing the availability of precursors for the de novo biosynthesis of purine nucleotides is a valid approach towards seeking a cure for leukaemia. Nucleotides and deoxynucleotides are required by living cells for syntheses of RNA, DNA, and cofactors such as NADP(+), FAD(+), coenzyme A and ATP. Nucleotides contain purine and pyrimidine bases, which can be synthesized through salvage pathway as well. Amido phosphoribosyltransferase (APRT), also known as glutamine phosphoribosylpyrophosphate amidotransferase (GPAT), is an enzyme that in humans is encoded by the PPAT (phosphoribosyl pyrophosphate amidotransferase) gene. APRT catalyzes the first committed step of the de novo pathway using its substrate, phosphoribosyl pyrophosphate (PRPP). As APRT is inhibited by many folate analogues, therefore, in this study we focused on the inhibitory effects of three folate analogues on APRT activity. This is extension of our previous wet lab work to analyze and dissect molecular interaction and inhibition mechanism using molecular modeling and docking tools in the current study. Comparative molecular docking studies were carried out for three diamino folate derivatives employing a model of the human enzyme that was built using the 3D structure of Bacillus subtilis APRT (PDB ID; 1GPH) as the template. Binding orientation of interactome indicates that all compounds having nominal cluster RMSD in same active site's deep narrow polar fissure. On the basis of comparative conformational analysis, electrostatic interaction, binding free energy and binding orientation of interactome, we support the possibility that these molecules could behave as APRT inhibitors and therefore may block purine de novo biosynthesis. Consequently, we suggest that PY899 is the most active biological compound that would be a more potent inhibitor for APRT inhibition than PY873 and DIA, which also confirms previous wet lab report.

Bioinformatic analysis of pathogenic missense mutations of activin receptor like kinase 1 ectodomain

Activin A receptor, type II-like kinase 1 (also called ALK1), is a serine-threonine kinase predominantly expressed on endothelial cells surface. Mutations in its ACVRL1 encoding gene (12q11-14) cause type 2 Hereditary Haemorrhagic Telangiectasia (HHT2), an autosomal dominant multisystem vascular dysplasia. The study of the structural effects of mutations is crucial to understand their pathogenic mechanism. However, while an X-ray structure of ALK1 intracellular domain has recently become available (PDB ID: 3MY0), structure determination of ALK1 ectodomain (ALK1_EC) has been elusive so far. We here describe the building of a homology model for ALK1_EC, followed by an extensive bioinformatic analysis, based on a set of 38 methods, of the effect of missense mutations at the sequence and structural level. ALK1_EC potential interaction mode with its ligand BMP9 was then predicted combining modelling and docking data. The calculated model of the ALK1_EC allowed mapping and a preliminary characterization of HHT2 associated mutations. Major structural changes and loss of stability of the protein were predicted for several mutations, while others were found to interfere mainly with binding to BMP9 or other interactors, like Endoglin (CD105), whose encoding ENG gene (9q34) mutations are known to cause type 1 HHT. This study gives a preliminary insight into the potential structure of ALK1_EC and into the structural effects of HHT2 associated mutations, which can be useful to predict the potential effect of each single mutation, to devise new biological experiments and to interpret the biological significance of new mutations, private mutations, or non-synonymous polymorphisms.

Prediction of the Human EP1 Receptor Binding Site by Homology Modeling and Molecular Dynamics Simulation

The prostanoid receptor EP1 is a G-protein-coupled receptor (GPCR) known to be involved in a variety of pathological disorders such as pain, fever and inflammation. These receptors are important drug targets, but design of subtype specific agonists and antagonists has been partially hampered by the absence of three-dimensional structures for these receptors. To understand the molecular interactions of the PGE2, an endogen ligand, with the EP1 receptor, a homology model of the human EP1 receptor (hEP1R) with all connecting loops was constructed from the 2.6 Å resolution crystal structure (PDB code: 1L9H) of bovine rhodopsin. The initial model generated by MODELLER was subjected to molecular dynamics simulation to assess quality of the model. Also, a step by step ligand-supported model refinement was performed, including initial docking of PGE2 and iloprost in the putative binding site, followed by several rounds of energy minimizations and molecular dynamics simulations. Docking studies were performed for PGE2 and some other related compounds in the active site of the final hEP1 receptor model. The docking enabled us to identify key molecular interactions supported by the mutagenesis data. Also, the correlation of r²=0.81 was observed between the Ki values and the docking scores of 15 prostanoid compounds. The results obtained in this study may provide new insights toward understanding the active site conformation of the hEP1 receptor and can be used for the structure-based design of novel specific ligands.

Homology modeling of human CCR5 and analysis of its binding properties through molecular docking and molecular dynamics simulation

In this study, homology modeling, molecular docking and molecular dynamics simulation were performed to explore structural features and binding mechanism of some inhibitors of chemokine receptor type 5 (CCR5), and to construct a model for designing new CCR5 inhibitors for preventing HIV attachment to the host cell. A homology modeling procedure was employed to construct a 3D model of CCR5. For this procedure, the X-ray crystal structure of bovine rhodopsin (1F88A) at 2.80Å resolution was used as template. After inserting the constructed model into a hydrated lipid bilayer, a 20ns molecular dynamics (MD) simulation was performed on the whole system. After reaching the equilibrium, twenty-four CCR5 inhibitors were docked in the active site of the obtained model. The binding models of the investigated antagonists indicate the mechanism of binding of the studied compounds to the CCR5 obviously. Moreover, 3D pictures of inhibitor-protein complex provided precious data regarding the binding orientation of each antagonist into the active site of this protein. One additional 20 ns MD simulation was performed on the initial structure of the CCR5-ligand 21 complex, resulted from the previous docking calculations, embedded in a hydrated POPE bilayer to explore the effects of the presence of lipid bilayer in the vicinity of CCR5-ligand complex. This article is part of a Special Issue entitled Protein translocation across or insertion into membranes.

Structural and functional characterization of HMG-COA reductase from Artemisia annua

Plants synthesize a great variety of isoprenoid products that are required not only for normal growth and development but also for their adaptive responses to environmental challenges. However, despite the remarkable diversity in the structure and function of plant isoprenoids, they all originate from a single metabolic precursor, mevalonic acid. The synthesis of mevalonic acid is catalysed by the enzyme, 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMG- CoA reductase). The analysis of the amino acid sequence of HMG-CoA reductase from Artemisia annua L. plant showed that it belongs to class I HMG-CoA reductase family. The three dimensional structure of HMG-CoA reductase of Artemisia annua has been generated from amino acid sequence using homology modelling with backbone structure of human HMG-CoA reductase as template. The model was generated using the SWISS MODEL SERVER. The generated 3-D structure of HMG-CoA reductase was evaluated at various web interfaced servers to checks the stereo interfaced quality of the structure in terms of bonds, bond angles, dihedral angles and non-bonded atom-atom distances, structural as well as functional domains etc. The generated model was visualized using the RASMOL. Structural analysis of HMG-CoA reductase from Artemisia annua L. plant hypothesize that the N and C-terminals are positioned in cytosol by the two membrane spanning helices and the C-terminals domain shows similarity to the human HMG-CoA reductase enzyme indicating that they both had potential catalytic similarities.

Analysis of the oxidase activity induced by CCl(4) and H(2)O(2) in different recombinant myoglobins

Hemoproteins may present several functions due to their prosthetic groups. After a long time, well-studied proteins such as myoglobin have surprised us with new functions. Myoglobin is a hemoprotein which has some well described and unexpected functions within the organism. Oxidase activity in standard myoglobins has been described and this activity was attributed to a covalent linkage between heme and some amino acid residues such as histidine, when myoglobins are treated with alkyl halides, and tyrosine, and when myoglobins are treated with H(2)O(2). We have found that the oxidase activity, due to H(2)O(2) treatment, can appear in different myoglobins, which presents no key residue, such as Tyr 103, for the oxidase activity previously described in the literature.

Molecular modeling and characterization of the B. thuringiensis and B. thuringiensis LDC-9 cytolytic proteins

The Cyt toxins are able to lyse a wide range of cell types in vitro, unlike the Cry delta-endotoxins. It exerts its activity by the formation of pores within target cell membranes. The structural information available for Cyt2Aa (PDB id: 1CBY) consists of a single domain in which two outer layers of alpha-helix wrap around a mixed beta-sheet. Beta-barrel was suggested as a possible structure of the pores. Hence, this study seeks to investigate the structural properties of other Cytolytic proteins by predicting the three-dimensional (3D) model using Cyt2Aa as template. The predicted models are expected to be significantly more accurate as all the Cyt proteins showed significant similarity with the template (PDB id: 1CBY). The refined homology models revealed similar secondary structures (alpha-helices and beta-sheets) and tertiary features as Cyt2Aa. The variation in the loop regions of the tertiary structure accounts for the differential toxicity.

Computational studies of the binding site of alpha1A-adrenoceptor antagonists

Aimed at achieving a good understanding of the 3-dimensional structures of human alpha1A-adrenoceptor (alpha1A-AR), we have successfully developed its homology model based on the crystal structure of beta2-AR. Subsequent structural refinements were performed to mimic the receptor's natural membrane environment by using molecular mechanics (MM) and molecular dynamics (MD) simulations in the GBSW implicit membrane model. Through molecular docking and further simulations, possible binding modes of subtype-selective alpha1A-AR antagonists, Silodosin, RWJ-69736 and (+)SNAP-7915, were examined. Results of the modeling and docking studies are qualitatively consistent with available experimental data from mutagenesis studies. The homology model built should be very useful for designing more potent subtype-selective alpha1A-AR antagonists and for guiding further mutagenesis studies.

A CD4 homologue in sea bass (Dicentrarchus labrax): molecular characterization and structural analysis

CD4 is a transmembrane glycoprotein fundamental for cell-mediated immunity. Its action as a T cell co-receptor increases the avidity of association between a T cell and an antigen-presenting cell by interacting with portions of the complex between MHC class II and TR molecules. In this paper we report the cDNA cloning, expression and structural analysis of a CD4 homologue from sea bass (Dicentrarchus labrax). The sea bass CD4 cDNA consists of 2071 bp that translates in one reading frame to give the entire molecule containing 480 amino acids. The analysis of the sequence shows the presence of four putative Ig-like domains and that some fundamental structural features, like a disulphide bond in domain D2 and the CXC signalling motif in the cytoplasmic tail, are conserved from sea bass to mammals. Real-time PCR analysis showed that very high levels of CD4 mRNA transcripts are present in thymus, followed by gut and gills. In vitro stimulation of head kidney leukocytes with LPS and PHA-L gave an increase of CD4 mRNA levels after 4h and a decrease after 24h. Homology modelling has been applied to create a 3D model of sea bass CD4 and to investigate its interaction with sea bass MHC-II. The analysis of the 3D complex between sea bass CD4 and sea bass MHC-II suggests that the absence of a disulfide bond in the CD4 D1 domain could make this molecule more flexible, inducing a different conformation and affecting the binding and the way of interaction between CD4 and MHC-II. Our results will add new insights into the sea bass T cell immune responses and will help in the identification of T cell subsets in teleost fishes to better understand the evolution of cell-mediated immunity from fish to mammals.

Dissecting the structural determinants of the interaction between the human cytomegalovirus UL18 protein and the CD85j immune receptor

UL18 is a glycoprotein encoded by the human cytomegalovirus genome and is thought to play a pivotal role during human cytomegalovirus infection, although its exact function is still a matter of debate. UL18 shares structural similarity with MHC class I and binds the receptor CD85j on immune cells. Besides UL18, CD85j binds MHC class I molecules. The binding properties of CD85j to MHC class I molecules have been thoroughly studied. Conversely, very little information is available on the CD85j/UL18 complex, namely that UL18 binds CD85j through its alpha3 domain with an affinity that is approximately 1000-fold higher than the MHC class I affinity for CD85j. Deeper knowledge of features of the UL18/CD85j complex would help to disclose the function of UL18 when it binds to CD85j. In this study we first demonstrated that the UL18alpha3 domain is not sufficient per se for binding and that beta2-microglobulin is necessary for UL18-CD85j interaction. We then dissected structural determinants of binding UL18 to CD85j. To this end, we constructed a three-dimensional model of the complex. The model was used to design mutants in selected regions of the putative interaction interface, the effects of which were measured on binding. Six regions in both the alpha2 and alpha3 domains and specific amino acids within them were identified that are potentially involved in the UL18-CD85j interaction. The higher affinity of UL18 to CD85j, compared with MHC class I, seems to be due not to additional interaction regions but to an overall better fit of the two molecules.

Molecular modelling and comparative structural account of aspartyl beta-semialdehyde dehydrogenase of Mycobacterium tuberculosis (H37Rv)

Aspartyl beta-semialdehyde dehydrogenase (ASADH) is an important enzyme, occupying the first branch position of the biosynthetic pathway of the aspartate family of amino acids in bacteria, fungi and higher plants. It catalyses reversible dephosphorylation of L: -beta-aspartyl phosphate (betaAP) to L: -aspartate-beta-semialdehyde (ASA), a key intermediate in the biosynthesis of diaminopimelic acid (DAP)-an essential component of cross linkages in bacterial cell walls. Since the aspartate pathway is unique to plants and bacteria, and ASADH is the key enzyme in this pathway, it becomes an attractive target for antimicrobial agent development. Therefore, with the objective of deducing comparative structural models, we have described a molecular model emphasizing the uniqueness of ASADH from Mycobacterium tuberculosis (H37Rv) that should generate insights into the structural distinctiveness of this protein as compared to structurally resolved ASADH from other bacterial species. We find that mtASADH exhibits structural features common to bacterial ASADH, while other structural motifs are not present. Structural analysis of various domains in mtASADH reveals structural conservation among all bacterial ASADH proteins. The results suggest that the probable mechanism of action of the mtASADH enzyme might be same as that of other bacterial ASADH. Analysis of the structure of mtASADH will shed light on its mechanism of action and may help in designing suitable antagonists against this enzyme that could control the growth of Mycobacterium tuberculosis.

An evaluation of automated homology modelling methods at low target template sequence similarity

MOTIVATION

There are two main areas of difficulty in homology modelling that are particularly important when sequence identity between target and template falls below 50%: sequence alignment and loop building. These problems become magnified with automatic modelling processes, as there is no human input to correct mistakes. As such we have benchmarked several stand-alone strategies that could be implemented in a workflow for automated high-throughput homology modelling. These include three new sequence-structure alignment programs: 3D-Coffee, Staccato and SAlign, plus five homology modelling programs and their respective loop building methods: Builder, Nest, Modeller, SegMod/ENCAD and Swiss-Model. The SABmark database provided 123 targets with at least five templates from the same SCOP family and sequence identities </=50%.

RESULTS

When using Modeller as the common modelling program, 3D-Coffee outperforms Staccato and SAlign using both multiple templates and the best single template, and across the sequence identity range 20-50%. The mean model RMSD generated from 3D-Coffee using multiple templates is 15 and 28% (or using single templates, 3 and 13%) better than those generated by Staccato and Salign, respectively. 3D-Coffee gives equivalent modelling accuracy from multiple and single templates, but Staccato and SAlign are more successful with single templates, their quality deteriorating as additional lower sequence identity templates are added. Evaluating the different homology modelling programs, on average Modeller performs marginally better in overall modelling than the others tested. However, on average Nest produces the best loops with an 8% improvement by mean RMSD compared to the loops generated by Builder.

A novel structure-based virtual screening model for the hERG channel blockers

The hERG potassium channel is a key effector of cardiac repolarization and the blockade of this channel could cause arrhythmia. Thus, hERG channel blockade plays an important role for the potential pro-arrhythmic liability. In this report, binding of blockers to the hERG potassium channel is investigated using a combination of homology modeling, molecular docking, and molecular simulations, where blockade activities are evaluated using the linear regression model of GoldScore fitness. This structure-based virtual screening model is able to estimate the pIC(50) value of a wide range of ligands for the hERG potassium channel. The docked poses for ligands are also consistent with published mutation. Therefore, this model for the prediction of hERG channel blockade has the potential to provide cost-effective virtual screening tools for the evaluation of the cardiac liability of new chemical entities.

Theoretical model of the three-dimensional structure of a disease resistance gene homolog encoding resistance protein in Vigna mungo

Plant disease resistance (R) genes, the key players of innate immunity system in plants encode 'R' proteins. 'R' protein recognizes product of avirulance gene from the pathogen and activate downstream signaling responses leading to disease resistance. No three dimensional (3D) structural information of any 'R' proteins is available as yet. We have reported a 'R' gene homolog, the 'VMYR1', encoding 'R' protein in Vigna mungo. Here, we describe the homology modeling of the 'VMYR1' protein. The model was created by using the 3D structure of an ATP-binding cassette transporter protein from Vibrio cholerae as a template. The strategy for homology modeling was based on the high structural conservation in the superfamily of P-loop containing nucleoside triphosphate hydrolase in which target and template proteins belong. This is the first report of theoretical model structure of any 'R' proteins.

Determinants of substrate specificity in KdcA, a thiamin diphosphate-dependent decarboxylase

Thiamin diphosphate-dependent decarboxylases catalyze the non-oxidative decarboxylation of 2-keto carboxylic acids. Although they display relatively low sequence similarity, and broadly different range of substrates, these enzymes show a common homotetrameric structure. Here we describe a kinetic characterization of the substrate spectrum of a recently identified member of this class, the branched chain 2-keto acid decarboxylase (KdcA) from Lactococcus lactis. In order to understand the structural basis for KdcA substrate recognition we developed a homology model of its structure. Ser286, Phe381, Val461 and Met358 were identified as residues that appeared to shape the substrate binding pocket. Subsequently, site-directed mutagenesis was carried out on these residues with a view to converting KdcA into a pyruvate decarboxylase. The results show that the mutations all lowered the Km value for pyruvate and both the S286Y and F381W variants also had greatly increased values of k(cat) with pyruvate as a substrate.

Molecular modeling and characterization of Vibrio cholerae transcription regulator HlyU

Background

The SmtB/ArsR family of prokaryotic metal-regulatory transcriptional repressors represses the expression of operons linked to stress-inducing concentrations of heavy metal ions, while derepression results from direct binding of metal ions by these 'metal-sensor' proteins. The HlyU protein from Vibrio cholerae is the positive regulator of haemolysin gene, it also plays important role in the regulation of expression of the virulence genes. Despite the understanding of biochemical properties, its structure and relationship to other protein families remain unknown.

Results

We find that HlyU exhibits structural features common to the SmtB/ArsR family of transcriptional repressors. Analysis of the modeled structure of HlyU reveals that it does not have the key metal-sensing residues which are unique to the SmtB/ArsR family of repressors, yet the tertiary structure is very similar to the family members. HlyU is the only member that has a positive control on transcription, while all the other members in the family are repressors. An evolutionary analysis with other SmtB/ArsR family members suggests that during evolution HlyU probably occurred by gene duplication and mutational events that led to the emergence of this protein from ancestral transcriptional repressor by the loss of the metal-binding sites.

Conclusion

The study indicates that the same protein family can contain both the positive regulator of transcription and repressors – the exact function being controlled by the absence or the presence of metal-binding sites.

Molecular architecture of the glucose 1-phosphate site in ADP-glucose pyrophosphorylases

ADP-Glc pyrophosphorylase (PPase), a key regulatory enzyme in the biosynthetic pathway of starch and bacterial glycogen, catalyzes the synthesis of ADP-Glc from Glc-1-P and ATP. A homology model of the three-dimensional structure of the Escherichia coli enzyme complexed with ADP-Glc has been generated to study the substrate-binding site in detail. A set of amino acids in the model has been identified to be in close proximity to the glucose moiety of the ADP-Glc ligand. The role of these amino acids (Glu(194), Ser(212), Tyr(216), Asp(239), Phe(240), Trp(274), and Asp(276)) was studied by site-directed mutagenesis through the characterization of the kinetic properties and thermal stability of the designed mutants. All purified alanine mutants had 1 or 2 orders of magnitude lower apparent affinity for Glc-1-P compared with the wild type, indicating that the selected set of amino acids plays an important role in their interaction with the substrate. These amino acids, which are conserved within the ADP-Glc PPase family, were replaced with other residues to investigate the effect of size, hydrophobicity, polarity, aromaticity, or charge on the affinity for Glc-1-P. In this study, the architecture of the Glc-1-P-binding site is characterized. The model overlaps with the Glc-1-P site of other PPases such as Pseudomonas aeruginosa dTDP-Glc PPase and Salmonella typhi CDP-Glc PPase. Therefore, the data reported here may have implications for other members of the nucleotide-diphosphoglucose PPase family.

Characterization and kinetics studies of water buffalo (Bubalus bubalis) myoglobin

The colour of buffalo (Bubalus bubalis L.) meat is darker than bovine meat. Since meat colour depends on the concentration of myoglobin (Mb) and its oxidation state, we have determined the main structural and functional properties of buffalo Mb. Buffalo Mb was purified from longissimus dorsi muscles and its molecular mass determined by ESI Q-TOF mass spectrometry. The molecular mass 17,034.50 was 86.20 Da higher than the bovine Mb. This was confirmed by analysing its primary structure, using a combined approach based on Edman degradation and MALDI-TOF mass spectrometry. Comparing the amino acid sequences of both Mbs, we found three amino acid differences out of 153 amino acid residues. One is a conservative substitution (D(bov)141E(buf)), and the other two (A(bov)19T(buf) and A(bov)117D(buf)) are nonconservative. These amino acid substitutions are unlikely to cause structural changes because they are located far from the heme binding pocket, as revealed by the 3D structure of buffalo Mb elaborated by homology modelling. Stability analyses show no difference with the bovine Mb for helix E and only minor differences in the stability values for helices A and G. Moreover, autoxidation rates of purified buffalo and bovine myoglobins at 37 degrees C, pH 7.2, were almost identical, 0.052+/-0.001 h(-1) and 0.054+/-0.002 h(-1), respectively, as were their oxygen-binding Kd values, 3.7+/-0.1 microM and 3.5+/-0.1 microM, respectively. The percent of MetMb values were almost identical. The results presented here suggest that the darker buffalo meat depends on factors other than the oxidation rate of its Mb, as, for example, the Mb content (0.393+/-0.005 g/100 g of tissue) and consequently MetMb, which are almost twice as high as bovine meat (Mb: 0.209+/-0.003 g/100 g of tissue).

FlgM anti-sigma factors: identification of novel members of the family, evolutionary analysis, homology modeling, and analysis of sequence-structure-function relationships

FlgM proteins, also known as Anti-sigma-28 factor (sigma28), are negative regulators of flagellin synthesis. Recently, a three-dimensional structure of the Aquifex aeolicus sigma28/FlgM complex (PDB code: 1rp3) was determined by X-ray crystallography at 2.3 A resolution. Furthermore, experimental data on bacterial FlgM, including site-directed mutagenesis and structural characterization by NMR are also available. However, an interpretation of the sequence-structure-function relationships combining X-ray and NMR data with the evolutionary information extracted from the increasing number of FlgM-related sequences annotated in databases is not available. In the present study, we combined database sequence searches and sequence-analysis tools to update the multiple sequence alignment of a previously characterized cluster of orthologs (COG2747) and the PFAM classification of protein domains (PF04316) for the FlgM family. A phylogenetic analysis of 77 protein sequences revealed the presence of at least three major sequence clades within the FlgM family. Besides, we predicted functional residues using a SequenceSpace method. We also generated homology models for Bacillus subtilis and Salmonella typhimurium FlgM proteins, for which sequence-structure-function relationship data are available, and used the docking program ClusPro to hypothesize about the dimer association between FlgM proteins. In conclusion, the analysis presented in this work will be useful in designing new experiments to understand better protein-protein interactions between FglM, sigma factors, and putative molecules from the flagellar export apparatus. Electronic Supplementary Material is available in the online version of this article at http://link.springer.de/

The role of molecular modelling in biomedical research

The synergy between experimental and computational biology has greatly benefited both fields, providing invaluable information in many different areas of the life sciences. This minireview will focus on one specific aspect of computational biology, molecular modelling, and describe a few examples highlighting the effectiveness of protein structural analysis and modelling in providing relevant information about systems of biomedical interest.

The CD8alpha from sea bass (Dicentrarchus labrax L.): Cloning, expression and 3D modelling

In this paper we describe the cloning, expression and structural study by modelling techniques of the CD8alpha from sea bass (Dicentrarchus labrax L.). The sea bass CD8alpha cDNA is comprised of 1490 bp and is translated in one reading frame to give a protein of 217 amino acids, with a predicted 26 amino acids signal peptide, a 88 bp 5'-UTR and a 748 bp 3'-UTR. A multiple alignment of CD8alpha from sea bass with other known CD8alpha sequences shows the conservation of most amino acid residues involved in the peculiar structural domains found within CD8alpha's. Cysteine residues that are involved in disulfide bonding to form the V domain are conserved. In contrast, an extra cysteine residue found in most mammals in this region is not present in sea bass. The transmembrane and cytoplasmic regions are the most conserved regions within the molecule in the alignment analysis. However, the motif (CXCP) that is thought to be responsible for binding p56lck is missing in the sea bass sequence. Phylogenetic analysis conducted using amino acid sequences showed that sea bass CD8alpha grouped with other known teleost sequences and that three different clusters were formed by the mammalian, avian and fish CD8alpha sequences. The thymus was the tissue with the highest CD8alpha expression, followed by gut, gills, peripheral blood leukocytes and spleen. Lower CD8alpha mRNA levels were found in head kidney, liver and brain. It was possible to create a partial 3D model using the human and mouse structures as template. The CD8alpha 11-120 amino acid region was taken into consideration and the best obtained 3D model shows the presence of ten beta-strands, involving about 50% of the sequence. The global structure was defined as an immunoglobulin-like beta-sandwich made of two anti-parallel sheets. Two cysteines were present in this region and they were at a suitable distance to form an S-S bond as seen in the template human and mouse structures.

Relationship between multiple sequence alignments and quality of protein comparative models

Comparative modeling is the method of choice, whenever applicable, for protein structure prediction, not only because of its higher accuracy compared to alternative methods, but also because it is possible to estimate a priori the quality of the models that it can produce, thereby allowing the usefulness of a model for a given application to be assessed beforehand. By and large, the quality of a comparative model depends on two factors: the extent of structural divergence between the target and the template and the quality of the sequence alignment between the two protein sequences. The latter is usually derived from a multiple sequence alignment (MSA) of as many proteins of the family as possible, and its accuracy depends on the number and similarity distribution of the sequences of the protein family. Here we describe a method to evaluate the expected difficulty, and by extension accuracy, of a comparative model on the basis of the MSA used to build it. The parameter that we derive is used to compare the results obtained in the last two editions of the Critical Assessment of Methods for Structure Prediction (CASP) experiment as a function of the difficulty of the modeling exercise. Our analysis demonstrates that the improvement in the scope and quality of comparative models between the two experiments is largely due to the increased number of available protein sequences and to the consequent increased chance that a large and appropriately spaced set of protein sequences homologous to the proteins of interest is available.

Unravelling the structure of the pneumococcal autolytic lysozyme

The LytC lysozyme of Streptococcus pneumoniae forms part of the autolytic system of this important pathogen. This enzyme is composed of a C-terminal CM (catalytic module), belonging to the GH25 family of glycosyl hydrolases, and an N-terminal CBM (choline-binding module), made of eleven homologous repeats, that specifically recognizes the choline residues that are present in pneumococcal teichoic and lipoteichoic acids. This arrangement inverts the general assembly pattern of the major pneumococcal autolysin, LytA, and the lytic enzymes encoded by pneumococcal bacteriophages that place the CBM (made of six repeats) at the C-terminus. In the present paper, a three-dimensional model of LytC built by homology modelling of each module and consistent with spectroscopic and hydrodynamic studies is shown. In addition, the putative catalytic-pair residues are identified. Despite the inversion in the modular arrangement, LytC and the bacteriophage-encoded Cpl-1 lysozyme most probably adopt a similar global fold. However, the distinct choline-binding ability and their substrate-binding surfaces may reflect a divergent evolution directed by the different roles played by them in the host (LytC) or in the bacteriophage (Cpl-1). The tight binding of LytC to the pneumococcal envelope, mediated by the acquisition of additional choline-binding repeats, could facilitate the regulation of the potentially suicidal activity of this autolysin. In contrast, a looser attachment of Cpl-1 to the cell wall and the establishment of more favourable interactions between its highly negatively charged catalytic surface and the positively charged chains of pneumococcal murein could enhance the lytic activity of the parasite-encoded enzyme and therefore liberation of the phage progeny.

A geometric invariant-based framework for the analysis of protein conformational space

MOTIVATION

Characterization of the restricted nature of the protein local conformational space has remained a challenge, thereby necessitating a computationally expensive conformational search in protein modeling. Moreover, owing to the lack of unilateral structural descriptors, conventional data mining techniques, such as clustering and classification, have not been applied in protein structure analysis.

RESULTS

We first map the local conformations in a fixed dimensional space by using a carefully selected suite of geometric invariants (GIs) and then reduce the number of dimensions via principal component analysis (PCA). Distribution of the conformations in the space spanned by the first four PCs is visualized as a set of conditional bivariate probability distribution plots, where the peaks correspond to the preferred conformations. The locations of the different canonical structures in the PC-space have been interpreted in the context of the weights of the GIs to the first four PCs. Clustering of the available conformations reveals that the number of preferred local conformations is several orders of magnitude smaller than that suggested previously.

SUPPLEMENTARY INFORMATION

www.it.iitb.ac.in/~ashish/bioinfo2005/.

Amino acid 36 in the human immunodeficiency virus type 1 gp41 ectodomain controls fusogenic activity: implications for the molecular mechanism of viral escape from a fusion inhibitor

We have previously described a human immunodeficiency virus type 1 (HIV-1) proviral clone, pL2, derived from defective viral particles with higher fusogenicity than the prototypic NL4-3 virus. In this study, we attempted to determine the region that confers the enhanced fusion activity by creating envelope recombinants between pL2 and pNL4-3, as well as point mutants based on pNL4-3. The results indicate that amino acid 36 of gp41 is key for the fusogenic activity and infectivity enhancement and that glycine 36 (36G) of gp41 in pL2 is conserved in nearly all HIV-1 isolates except for pNL4-3. The mutation 36G-->D in a primary-isolate-derived Env decreased syncytium-forming activity and infectivity. The assays for cell-cell fusion and viral binding suggested that the enhanced fusion mediated by the 36D-->G mutation is not due to increased binding efficiency but is directly due to actual enhancement of viral fusion activity. Interestingly, this amino acid position is exactly equivalent to that at which the mutation of HIV-1 isolates that have escaped from a fusion inhibitor, enfuvirtide (T-20), has been frequently observed. The correlation between these previous findings and our findings was suggested by structural analysis. Our finding, therefore, has implications for a molecular basis of the viral escape from this drug.

Modelling of HLA-DQ2 and its interaction with gluten peptides to explain molecular recognition in celiac disease

Celiac disease (CD) is sustained by abnormal intestinal mucosal T-cell response to gluten and it is strongly associated with HLA class II molecules encoded by DQA1*0501/DQB1*02 (DQ2) or DQA1*03/DQB1*0302 (DQ8). The in vitro stimulatory activity of gliadin increases after treatment with tissue transglutaminase (tTG) which catalyses the deamidation of specific residues of glutamine to glutamate that can serve as anchors for binding to DQ2 as well as to DQ8 molecules. We modelled the three-dimensional structure of the DQ2 dimer protein, the most frequent in celiac patients, by using a homology modelling strategy, and deposited the model in the Protein Data Bank (PDB). Then, we simulated the interactions of DQ2 with different gluten peptides and the deamidation of specific peptide glutamines in the known p4, p6, p7 and p9 anchor positions, as well as in p1 and p5 positions, and other substitutions for which experimental effects on binding are available by previous experimental studies. By evaluating the energy of interaction and the H-bond interactions, we were able to distinguish what substitutions improve the interaction peptide-DQ2, in agreement with previously published experimental data. By analysing the peptide-DQ2 complex at the atom level, we observed that these glutamate side chains can interact with specific positively charged amino acids of DQ2, absent in other HLA alleles not related to celiac disease. The simulation was also extended to other peptides, related to celiac disease but for which no experimental data exists about the effects of glutamine deamidation. Our results give an interpretation at the molecular level of previously reported binding experimental data and open a new window to gain further insights about peptide recognition in celiac disease.