BMC Structural Biology最新文献

Antibody, the primary effector molecule of the immune system, evolves after initial encounter with the antigen from a precursor form to a mature one to effectively deal with the antigen. Antibodies of a lineage diverge through antigen-directed isolated pathways of maturation to exhibit distinct recognition potential. In the context of evolution in immune recognition, diversity of antigen cannot be ignored. While there are reports on antibody lineage, structural perspective with respect to diverse recognition potential in a lineage has never been studied. Hence, it is crucial to evaluate how maturation leads to topological tailoring within a lineage enabling them to interact with significantly distinct antigens.

A data-driven approach was undertaken for the study. Global experimental mouse and human antibody-antigen complex structures from PDB were compiled into a coherent database of germline-linked antibodies bound with distinct antigens. Structural analysis of all lineages showed variations in CDRs of both H and L chains. Observations of conformational adaptation made from analysis of static structures were further evaluated by characterizing dynamics of interaction in two lineages, mouse V_H1–84 and human V_H5–51. Sequence and structure analysis of the lineages explained that somatic mutations altered the geometries of individual antibodies with common structural constraints in some CDRs. Additionally, conformational landscape obtained from molecular dynamics simulations revealed that incoming pathogen led to further conformational divergence in the paratope (as observed across datasets) even while maintaining similar overall backbone topology. MM-GB/SA analysis showed binding energies to be in physiological range. Results of the study are coherent with experimental observations.

The findings of this study highlight basic structural principles shaping the molecular evolution of a lineage for significantly diverse antigens. Antibodies of a lineage follow different developmental pathways while preserving the imprint of the germline. From the study, it can be generalized that structural diversification of the paratope is an outcome of natural selection of a conformation from an available ensemble, which is further optimized for antigen interaction. The study establishes that starting from a common lineage, antibodies can mature to recognize a wide range of antigens. This hypothesis can be further tested and validated experimentally.

In the backdrop of challenge to obtain a protein structure under the known limitations of both experimental and theoretical techniques, the need of a fast as well as accurate protein structure evaluation method still exists to substantially reduce a huge gap between number of known sequences and structures. Among currently practiced theoretical techniques, homology modelling backed by molecular dynamics based optimization appears to be the most popular one. However it suffers from contradictory indications of different validation parameters generated from a set of protein models which are predicted against a particular target protein. For example, in one model Ramachandran Score may be quite high making it acceptable, whereas, its potential energy may not be very low making it unacceptable and vice versa. Towards resolving this problem, the main objective of this study was fixed as to utilize a simple experimentally derived output, Surface Roughness Index of concerned protein of unknown structure as an intervening agent that could be obtained using ordinary microscopic images of heat denatured aggregates of the same protein.

It was intriguing to observe that direct experimental knowledge of the concerned protein, however simple it may be, might give insight on acceptability of its particular structural model out of a confusion set of models generated from database driven comparative technique for structure prediction. The result obtained from a widely varying structural class of proteins indicated that speed of protein structure evaluation can be further enhanced without compromising with accuracy by recruiting simple experimental output.

In this work, a semi-empirical methodological approach was provided for improving protein structure evaluation. It showed that, once structure models of a protein were obtained through homology technique, the problem of selection of a best model out of a confusion set of Pareto-optimal structures could be resolved by employing a structure agent directly obtainable through experiment with the same protein as experimental ingredient. Overall, in the backdrop of getting a reasonably accurate protein structure of pathogens causing epidemics or biological warfare, such approach could be of use as a plausible solution for fast drug design.

Aryl hydrocarbon receptor (AhR) ligands may act as potential carcinogens or anti-tumor agents. Understanding how some of the residues in AhR ligand binding domain (AhRLBD) modulate their interactions with ligands would be useful in assessing their divergent roles including toxic and beneficial effects. To this end, we have analysed the nature of AhRLBD interactions with 2,3,7,8-tetrachlorodibenzo-ρ-dioxin (TCDD), 6-formylindolo[3,2-b]carbazole (FICZ), indole-3-carbinol (I3C) and its degradation product, 3,3′-diindolylmethane (DIM), Resveratrol (RES) and its analogue, Piceatannol (PTL) using molecular modeling approach followed by molecular dynamic simulations.

Results showed that each of the AhR ligands, TCDD, FICZ, I3C, DIM, RES and PTL affect the local and global conformations of AhRLBD.

The data presented in this study provide a structural understanding of AhR with its ligands and set the basis for its functions in several pathways and their related diseases.

Earlier based on bioinformatics analyses, we had predicted the Mycobacterium tuberculosis (M.tb) proteins; Rv1555 and Rv1554, among the potential new tuberculosis drug targets. According to the ‘TB-drugome’ the Rv1555 protein is ‘druggable’ with sildenafil (Viagra), tadalafil (Cialis) and vardenafil (Levitra) drugs. In the present work, we intended to understand via computer modeling studies, how the above drugs are likely to inhibit the M.tb protein’s function.

The three-dimensional computer models for M.tb proteins; Rv1555 and Rv1554 constructed on the template of equivalent membrane anchor subunits of the homologous E.coli quinol fumarate reductase respiratory protein complex, followed by drug docking analyses, suggested that the binding of above drugs interferes with quinol binding sites. Also, we experimentally observed the in-vitro growth inhibition of E.coli bacteria containing the homologous M.tb protein sequences with sildenafil and tadalafil drugs.

The predicted binding sites of the drugs is likely to affect the above M.tb proteins function as quinol binding is known to be essential for electron transfer function during anaerobic respiration in the homologous E.coli protein complex. Therefore, sildenafil and related drugs currently used in the treatment of male erectile dysfunction targeting the human phosphodiesterase 5 enzyme may be evaluated for their plausible role as repurposed drugs to treat human tuberculosis.

Aquaporins (AQPs) facilitate the passage of small neutral polar molecules across membranes of the cell. In animals there are four distinct AQP subfamilies, whereof AQP8 homologues constitute one of the smallest subfamilies with just one member in man. AQP8 conducts water, ammonia, urea, glycerol and H₂O₂ through various membranes of animal cells. This passive channel has been connected to a number of phenomena, such as volume change of mitochondria, ammonia neurotoxicity, and mitochondrial dysfunction related to oxidative stress. Currently, there is no experimentally determined structure of an AQP8, hence the structural understanding of this subfamily is limited. The recently solved structure of the plant AQP, AtTIP2;1, which has structural and functional features in common with AQP8s, has opened up for construction of homology models that are likely to be more accurate than previous models.

Here we present homology models of seven vertebrate AQP8s. Modeling based on the AtTIP2;1 structure alone resulted in reasonable models except for the pore being blocked by a phenylalanine that is not present in AtTIP2;1. To achieve an open pore, these models were supplemented with models based on the bacterial water specific AQP, EcAqpZ, creating a chimeric monomeric model for each AQP8 isoform. The selectivity filter (also named the aromatic/arginine region), which defines the permeant substrate profile, comprises five amino acid residues in AtTIP2;1, including a histidine coming from loop C. Compared to AtTIP2;1, the selectivity filters of modelled AQP8s only deviates in that they are slightly more narrow and more hydrophobic due to a phenylalanine replacing the histidine from loop C. Interestingly, the models do not exclude the existence of a side pore beneath loop C similar to that described in the structure of AtTIP2;1.

Our models concur that AQP8s are likely to have an AtTIP2;1-like selectivity filter. The detailed description of the expected configuration of residues in the selectivity filters of AQP8s provides an excellent starting point for planning of as well as rationalizing the outcome of mutational studies. Our strategy to compile hybrid models based on several templates may prove useful also for other AQPs for which structural information is limited.

The 5′-methylthioadenosine phosphorylase (MTAP), an enzyme involved in purine and polyamine metabolism and in the methionine salvage pathway, is considered as a potential drug target against cancer and trypanosomiasis. In fact, Trypanosoma and Leishmania parasites lack de novo purine pathways and rely on purine salvage pathways to meet their requirements. Herein, we propose the first comprehensive bioinformatic and structural characterization of the putative Leishmania infantum MTAP (LiMTAP), using a comparative computational approach.

Sequence analysis showed that LiMTAP shared higher identity rates with the Trypanosoma brucei (TbMTAP) and the human (huMTAP) homologs as compared to the human purine nucleoside phosphorylase (huPNP). Motifs search using MEME identified more common patterns and higher relatedness of the parasite proteins to the huMTAP than to the huPNP. The 3D structures of LiMTAP and TbMTAP were predicted by homology modeling and compared to the crystal structure of the huMTAP. These models presented conserved secondary structures compared to the huMTAP, with a similar topology corresponding to the Rossmann fold. This confirmed that both LiMTAP and TbMTAP are members of the NP-I family. In comparison to the huMTAP, the 3D model of LiMTAP showed an additional α-helix, at the C terminal extremity. One peptide located in this specific region was used to generate a specific antibody to LiMTAP. In comparison with the active site (AS) of huMTAP, the parasite ASs presented significant differences in the shape and the electrostatic potentials (EPs). Molecular docking of 5′-methylthioadenosine (MTA) and 5′-hydroxyethylthio-adenosine (HETA) on the ASs on the three proteins predicted differential binding modes and interactions when comparing the parasite proteins to the human orthologue.

This study highlighted significant structural peculiarities, corresponding to functionally relevant sequence divergence in LiMTAP, making of it a potential drug target against Leishmania.

Membrane proteins are difficult targets for structure prediction due to the limited structural data deposited in Protein Data Bank. Most computational methods for membrane protein structure prediction are based on the comparative modeling. There are only few de novo methods targeting that distinct protein family. In this work an example of such de novo method was used to structurally and functionally characterize two representatives of distinct membrane proteins families of solute carrier transporters and G protein-coupled receptors. The well-known Rosetta program and one of its protocols named Broker was used in two test cases. The first case was de novo structure prediction of three N-terminal transmembrane helices of the human concentrative nucleoside transporter 3 (hCNT3) homotrimer belonging to the solute carrier 28 family of transporters (SLC28). The second case concerned the large scale refinement of transmembrane helices of a homology model of the corticotropin-releasing factor receptor 1 (CRFR1) belonging to the G protein-coupled receptors family.

The inward-facing model of the hCNT3 homotrimer was used to propose the functional impact of its single nucleotide polymorphisms. Additionally, the 100?ns molecular dynamics simulation of the unliganded hCNT3 model confirmed its validity and revealed mobility of the selected binding site and homotrimer interface residues. The large scale refinement of transmembrane helices of the CRFR1 homology model resulted in the significant improvement of its accuracy with respect to the crystal structure of CRFR1, especially in the binding site area. Consequently, the antagonist CP-376395 could be docked with Autodock VINA to the CRFR1 model without any steric clashes.

The presented work demonstrated that Rosetta Broker can be a versatile tool for solving various issues referring to protein biology. Two distinct examples of de novo membrane protein structure prediction presented here provided important insights into three major areas of protein biology. Namely, the dynamics of the inward-facing hCNT3 homotrimer system, the structural changes of the CRFR1 receptor upon the antagonist binding and finally, the role of single nucleotide polymorphisms in both, hCNT3 and CRFR1 proteins, were investigated.

T7 lysozyme (T7L), also known as N-acetylmuramoyl-L-alanine amidase, is a T7 bacteriophage gene product. It involves two functions: It can cut amide bonds in the bacterial cell wall and interacts with T7 RNA polymerase (T7RNAP) as a part of transcription inhibition. In this study, with the help of molecular dynamics (MD) calculations and computational interaction studies, we investigated the effect of varying pH conditions on conformational flexibilities of T7L and their influence on T7RNAP -T7L interactions.

From the MD studies of the T7L at three different pH strengths viz. 5, neutral and 7.9 it was observed that T7L structure at pH?5 exhibited less stable nature with more residue level fluctuations, decrease of secondary structural elements and less compactness as compared to its counterparts: neutral pH and pH?7.9. The T-pad analysis of the MD trajectories identified local fluctuations in few residues that influenced the conformational differences in three pH strengths. From the docking of the minimum energy representative structures of T7L at different pH strengths (obtained from the free energy landscape analysis) with T7RNAP structures at same pH strengths, we saw strong interaction patterns at pH?7.9 and pH?5. The MD analysis of these complexes also confirmed the observations of docking study. From the combined in silico studies, it was observed that there are conformational changes in N-terminal and near helix 1 of T7L at different pH strengths, which are involved in the T7RNAP interaction, thereby varying the interaction pattern.

Since T7L has been used for developing novel therapeutics and T7RNAP one of the most biologically useful protein in both in-vitro and in vivo experiments, this in silico study of pH dependent conformational differences in T7L and the differential interaction with T7RNAP at different pH can provide a significant insight into the structural investigations on T7L and T7RNAP in varying pH environments.

The bHLH transcription factor TWIST1 plays a key role in the embryonic development and in tumorigenesis. Some loss-of-function mutations of the TWIST1 gene have been shown to cause an autosomal dominant craniosynostosis, known as the Saethre-Chotzen syndrome (SCS). Although the functional impacts of many TWIST1 mutations have been experimentally reported, little is known on the molecular mechanisms underlying their loss-of-function. In a previous study, we highlighted the predictive value of in silico molecular dynamics (MD) simulations in deciphering the molecular function of TWIST1 residues.

Here, since the substitution of the arginine 154 amino acid by a glycine residue (R154G) is responsible for the SCS phenotype and the substitution of arginine 154 by a proline experimentally decreases the dimerizing ability of TWIST1, we investigated the molecular impact of this point mutation using MD approaches. Consistently, MD simulations highlighted a clear decrease in the stability of the α-helix during the dimerization of the mutated R154P TWIST1/E12 dimer compared to the wild-type TE complex, which was further confirmed in vitro using immunoassays.

Our study demonstrates that MD simulations provide a structural explanation for the loss-of-function associated with the SCS TWIST1 mutation and provides a proof of concept of the predictive value of these MD simulations. This in silico methodology could be used to determine reliable pharmacophore sites, leading to the application of docking approaches in order to identify specific inhibitors of TWIST1 complexes.

Formamidopyrimidine-DNA glycosylase (Fpg) removes abundant pre-mutagenic 8-oxoguanine (oxoG) bases from DNA through nucleophilic attack of its N-terminal proline at C1′ of the damaged nucleotide. Since oxoG efficiently pairs with both C and A, Fpg must excise oxoG from pairs with C but not with A, otherwise a mutation occurs. The crystal structures of several Fpg–DNA complexes have been solved, yet no structure with A opposite the lesion is available.

Here we use molecular dynamic simulation to model interactions in the pre-catalytic complex of Lactococcus lactis Fpg with DNA containing oxoG opposite C or A, the latter in either syn or anti conformation. The catalytic dyad, Pro1–Glu2, was modeled in all four possible protonation states. Only one transition was observed in the experimental reaction rate pH dependence plots, and Glu2 kept the same set of interactions regardless of its protonation state, suggesting that it does not limit the reaction rate. The adenine base opposite oxoG was highly distorting for the adjacent nucleotides: in the more stable syn models it formed non-canonical bonds with out-of-register nucleotides in both the damaged and the complementary strand, whereas in the anti models the adenine either formed non-canonical bonds or was expelled into the major groove. The side chains of Arg109 and Phe111 that Fpg inserts into DNA to maintain its kinked conformation tended to withdraw from their positions if A was opposite to the lesion. The region showing the largest differences in the dynamics between oxoG:C and oxoG:A substrates was unexpectedly remote from the active site, located near the linker joining the two domains of Fpg. This region was also highly conserved among 124 analyzed Fpg sequences. Three sites trapping water molecules through multiple bonds were identified on the protein–DNA interface, apparently helping to maintain enzyme-induced DNA distortion and participating in oxoG recognition.

Overall, the discrimination against A opposite to the lesion seems to be due to incorrect DNA distortion around the lesion-containing base pair and, possibly, to gross movement of protein domains connected by the linker.