Comparative Protein Modeling Research Articles

Candidate genes involved in biosynthesis and degradation of the main extracellular matrix polysaccharides of brown algae and their probable evolutionary history

BackgroundBrown algae belong to the Stramenopiles phylum and are phylogenetically distant from plants and other multicellular organisms. This independent evolutionary history has shaped brown algae with numerous metabolic characteristics specific to this group, including the synthesis of peculiar polysaccharides contained in their extracellular matrix (ECM). Alginates and fucose-containing sulphated polysaccharides (FCSPs), the latter including fucans, are the main components of ECMs. However, the metabolic pathways of these polysaccharides remain poorly described due to a lack of genomic data.ResultsAn extensive genomic dataset has been recently released for brown algae and their close sister species, for which we previously performed an expert annotation of key genes involved in ECM-carbohydrate metabolisms. Here we provide a deeper analysis of this set of genes using comparative genomics, phylogenetics analyses, and protein modelling. Two key gene families involved in both the synthesis and degradation of alginate were suggested to have been acquired by the common ancestor of brown algae and their closest sister species Schizocladia ischiensis. Our analysis indicates that this assumption can be extended to additional metabolic steps, and thus to the whole alginate metabolic pathway. The pathway for the biosynthesis of fucans still remains biochemically unresolved and we also investigate putative fucosyltransferase genes that may harbour a fucan synthase activity in brown algae.ConclusionsOur analysis is the first extensive survey of carbohydrate-related enzymes in brown algae, and provides a valuable resource for future research into the glycome and ECM of brown algae. The expansion of specific families related to alginate metabolism may have represented an important prerequisite for the evolution of developmental complexity in brown algae. Our analysis questions the possible occurrence of FCSPs outside brown algae, notably within their closest sister taxon and in other Stramenopiles such as diatoms. Filling this knowledge gap in the future will help determine the origin and evolutionary history of fucan synthesis in eukaryotes.

Plastid-encoded RNA polymerase variation in Pelargonium sect Ciconium

Cyto-Nuclear Incompatibility (CNI), in which there is a mismatch in the interaction between organelles and nucleus, impacts plant species evolution as it has a direct effect on the fitness of plants. It can reduce fertility and/or result in bleached plants devoid of functional chloroplasts. Understanding the processes leading to CNI could help to improve breeding efforts, especially in cases where species with desirable traits need to be crossed into existing cultivars. To better understand the occurrence of CNI and its effects on plant phenotype, we combined near comprehensive crossing series across a clade of species from Pelargonium section Ciconium with comparative genomics and protein modelling for plastid-encoded RNA polymerase (PEP), as the rpo genes encoding PEP subunits were found to be unusually highly divergent, especially in two length-variable regions. Of all plastome-encoded genes, we found these genes to contain more variation than observed across angiosperms and that this underlies structural variation inferred for PEP in P. sect. Ciconium. This variation, resulting in differing physico-chemical properties of the rpo-encoded peptides, provides a possible explanation for the observed CNI, but we cannot directly correlate plastid related CNI phenotypes to rpo genotypes. This suggests that more than one interaction between the nuclear genome and the plastome genes are needed to fully explain the observed patterns.

Comprehensive Genome-Wide Analysis and Expression Pattern Profiling of PLATZ Gene Family Members in Solanum Lycopersicum L. under Multiple Abiotic Stresses.

PLATZ (plant AT-rich sequence and zinc-binding) family proteins with two conserved zinc-dependent DNA-binding motifs are transcription factors specific to the plant kingdom. The functions of PLATZ proteins in growth, development, and adaptation to multiple abiotic stresses have been investigated in various plant species, but their role in tomato has not been explored yet. In the present work, 20 non-redundant Solanum lycopersicum PLATZ (SlPLATZ) genes with three segmentally duplicated gene pairs and four tandemly duplicated gene pairs were identified on eight tomato chromosomes. The comparative modeling and gene ontology (GO) annotations of tomato PLATZ proteins indicated their probable roles in defense response, transcriptional regulation, and protein metabolic processes as well as their binding affinity for various ligands, including nucleic acids, peptides, and zinc. SlPLATZ10 and SlPLATZ17 were only expressed in 1 cm fruits and flowers, respectively, indicating their preferential involvement in the development of these organs. The expression of SlPLATZ1, SlPLATZ12, and SlPLATZ19 was up- or down-regulated following exposure to various abiotic stresses, whereas that of SlPLATZ11 was induced under temperature stresses (i.e., cold and heat stress), revealing their probable function in the abiotic stress tolerance of tomato. Weighted gene co-expression network analysis corroborated the aforementioned findings by spotlighting the co-expression of several stress-associated genes with SlPLATZ genes. Confocal fluorescence microscopy revealed the localization of SlPLATZ−GFP fusion proteins in the nucleus, hinting at their functions as transcription factors. These findings provide a foundation for a better understanding of the structure and function of PLATZ genes and should assist in the selection of potential candidate genes involved in the development and abiotic stress adaptation in tomato.

Bioinformatics analyses of rice calmodulin-lysine N-methyltransferases gene (OsCaM KMT) and characterization of its promoter from Oryza sativa L. sub sp. indica through transgenic approach

• Cloning of OsCaM KMT gene and promoter region from “Swarna”, an indica rice variety. • Presence of OsCAM KMT transcript in meristematic zone of rice plants, with higher expression in root tips suggesting they might be involved in auxin responsiveness. • Promoter analysis of OsCaM KMT::GUS construct depicted strong expression of GUS reporter gene in root tip, lateral roots and young leaf tip region. • Sequencing result indicated single nucleotide difference between the OsCAM KMT CDS from indica and japonica rice variety, resulting in silent mutation. Calmodulin-lysine N-methyltransferases (CaM KMT) belong to a family of protein lysine methyltransferases that catalyzes methylation of non-histone proteins. Lysine present at 116 position in the protein sequence of calmodulin is apparently methylated by CaM KMT in most species, however not characterized in rice. In the present study rice calmodulin-lysine N-methyltransferase ( OsCaM KMT ) gene has been isolated from an indica rice variety ‘Swarna’. Derived amino acid sequences of CaM KMT gene from rice and other plants were used in constructing a phylogenetic tree to show the evolutionary relationship among those species. Motif analysis along with comparative protein modeling was done to find out structural relationships among rice, Arabidopsis, and human CaM KMT. Gene expression analysis revealed differential tissue specific expression of OsCaM KMT . The highest transcript expression was revealed in the root tip followed by flower and new leaf, while, the least expression was found in the old leaf. In silico analysis was performed considering the 2 kilobases (kb) upstream promoter region of OsCaM KMT gene to find out various transcription factor binding sites. Analyses of the putative cis-regulatory elements revealed their involvement in tissue specific and stress dependent responses. For further characterization, the promoter region (2 kb) of OsCaM KMT gene was isolated, cloned, and sequenced. Histochemical GUS detection of transgenic rice depicted strong OsCaM KMT promoter expression in midribs, stomata, and tip of leaf tissues, central cylinder as well as tip parts of the primary root and lateral roots, along with significant expression in pollen grains, ovary, and vascular bundles. OsCaM KMT gene expression was corroborated by the histological GUS detection in OsCaM KMT promoter:: GUS fusion construct. Although several putative microRNAs (miRNAs) has been identified targeting OsCaM KMT coding DNA sequence (CDS) and UTR region, further studies are required to elucidate their tissue-specific regulation.

Identification and prioritization of promising lead molecules from Syzygium aromaticum against Sortase C from Streptococcus pyogenes: an in silico investigation

Sortases are extracellular transpeptidases that play an essential role in the adhesion of secreted proteins to the peptidoglycan layer of the cell wall of Gram-negative bacteria. Sortases are an important drug target protein due to their involvement in synthesizing the peptidoglycan cell wall of Streptococcus pyogenes, and these are not found in Homo sapiens. In this study, initially, we have performed protein sequence analysis to understand the sequential properties of Sortase C. Next, a comparative protein modeling approach was used to predict the three-dimensional model of Sortase C based on the crystal structure of Sortase C from Streptococcus pneumoniae. Virtual screening with an in-house library of phytochemicals from Syzygium aromaticum and molecular docking studies were performed to identify the promising lead molecules. These compounds were also analyzed for their drug-like and pharmacokinetic properties. Subsequently, the protein-ligand complexes of the selected ligands were subjected to molecular dynamics (MD) simulations to investigate their dynamic behavior in physiological conditions. The global and essential dynamics analyses result implied that the Sortase C complexes of the proposed three lead candidates exhibited adequate stability during the MD simulations. Additionally, the three proposed molecules showed favorable MM/PBSA binding free energy values ranging from −13.8 +/- 9.41 to −56.6 +/- 8.82 kcal/mol. After an extensive computational investigation, we have identified three promising lead candidates (CID:13888122, CID:3694932 and CID:102445430) against Sortase C from S. pyogenes. The result obtained from these computational studies can be used to screen and develop the inhibitors against Sortase C from S. pyogenes. Communicated by Ramaswamy H. Sarma

Amalgamation of comparative protein modeling with quantitative structure-retention relationship for prediction of the chromatographic behavior of peptides

Peptide therapeutics plays a prominent role in medical practice. Both peptides and proteins have been used in several disease conditions like diabetes, cancer, bacterial infections etc. The optimization of a peptide library is a time consuming and expensive chore. The tools of computational chemistry offer a way to optimize the properties of peptides. Quantitative Structure Retention (Chromatographic) Relationships (QSRR) is a powerful tool which statistically derives relationships between chromatographic parameters and descriptors that characterize the molecular structure of analytes. In this paper, we show how Comparative Protein ModelingQuantitative Structure Retention Relationship (acronym ComProM-QSRR) can be used to predict the retention time of peptide sequences. This formalism is founded on our earlier published QSAR methodology HomoSAR. ComProM-QSRR can recognize and distinguish the contribution of amino acids at specific positions in the peptide sequences to the retention phenomena through their related physicochemical properties. This study firmly establishes the fact that this approach can be pragmatically used to predict the retention time to all classes of peptides regardless of size or sequence.

Acetylation of conserved lysines fine-tunes mitochondrial malate dehydrogenase activity in land plants.

Plants need to rapidly and flexibly adjust their metabolism to changes of their immediate environment. Since this necessity results from the sessile lifestyle of land plants, key mechanisms for orchestrating central metabolic acclimation are likely to have evolved early. Here, we explore the role of lysine acetylation as a post-translational modification to directly modulate metabolic function. We generated a lysine acetylome of the moss Physcomitrium patens and identified 638 lysine acetylation sites, mostly found in mitochondrial and plastidial proteins. A comparison with available angiosperm data pinpointed lysine acetylation as a conserved regulatory strategy in land plants. Focusing on mitochondrial central metabolism, we functionally analyzed acetylation of mitochondrial malate dehydrogenase (mMDH), which acts as a hub of plant metabolic flexibility. In P. patens mMDH1, we detected a single acetylated lysine located next to one of the four acetylation sites detected in Arabidopsis thaliana mMDH1. We assessed the kinetic behavior of recombinant A. thaliana and P. patens mMDH1 with site-specifically incorporated acetyl-lysines. Acetylation of A. thaliana mMDH1 at K169, K170, and K334 decreases its oxaloacetate reduction activity, while acetylation of P. patens mMDH1 at K172 increases this activity. We found modulation of the malate oxidation activity only in A. thaliana mMDH1, where acetylation of K334 strongly activated it. Comparative homology modeling of MDH proteins revealed that evolutionarily conserved lysines serve as hotspots of acetylation. Our combined analyses indicate lysine acetylation as a common strategy to fine-tune the activity of central metabolic enzymes with likely impact on plant acclimation capacity.

In Silico prediction and Comparative Modeling of proteins in Creeping Fig (Ficus pumila) plant

In the present study, the Ficus pumila have taken to analyze the proteins by their preliminary characters from the database and predicted vital role of different sequences. The F. pumila (Creeping fig) is a prostrate/climbing shrub, experiments proved various active phytochemicals and antioxidant, antimicrobial, antimutagenic, analgesic, anti-inflammatory, antiproliferative, hypoglycemic, hypolipidemic, antihyperprolactinemic, anticholinesterase, nephroprotective properties. In addition to all metabolites, it also constitutes specific proteins that were evaluated through insilico homology modeling. Though it is considered as a poisonous weed, the protein present in this plant is evaluated by physicochemical, phylogeny and amino acid proportions by protparam, Swiss model, SOPMA, Clustal omega tools to describe its structural features and to understand molecular function. The computed theoretical isoelectric point (pI) found to be more than 7 indicates basic nature of proteins. The aliphatic index ranges 67-113 indicates thermal stability of proteins. The predicted Grand average hydropathy(GRAVY) shows possibilities of enhanced interaction of these proteins with water by lowest value. Functional analysis of these proteins was performed by SOSUI server which predicted transmembrane helix and solubility. Secondary structure analysis was carried out by SOPMA revealed that Alpha helix and random coil dominated followed by extended strand, and beta turns among secondary structure elements. The modelling of three-dimensional structure of proteins was performed by Swiss model. The model was validated using protein structure checking tool- VADAR. Particularly, NAD(P)H –quinoneoxidoreductase and Glyceraldehyde-3-phosphate dehydrogenase structures were analysed by phylogenetic analysis to trace relationship and reported. The results suggesting its possible role in cellular and metabolic functions.

Characterization and in planta validation of a CHI4 chitinase from Cajanus platycarpus (Benth.) Maesen for its efficacy against pod borer, Helicoverpa armigera (Hübner).

Pigeonpea, Cajanus cajan is one of the economically important legume food crops and a major source of dietary proteins. Management of pod borer, Helicoverpa armigera has been prominent among crop improvement programs. Lack of resistance sources in the cultivated germplasm and crossing incompatibility with pod borer-resistant wild relatives have prompted biotechnological interventions. Identification and exploitation of genes from pigeonpea wild relatives in host plant resistance towards the pod borer assumes pertinence. Dynamic transcriptome analysis of the wild relative vis a vis cultivated pigeonpea identified a CHI4 chitinase as one of the putative insect resistance genes. The study presents variations in important amino acids in CHI4 chitinases from C. cajan and its wild relative C. platycarpus. Comparative protein modeling and docking analysis of the two proteins demonstrated differences in substrate binding efficacy of the chitinase from C. platycarpus which resulted in a minimum binding energy of -8.7 kcal mol-1 . Furthermore, we successfully evaluated the insecticidal activity of the chitinase from C. platycarpus against H. armigera challenge through heterologous expression in tobacco. Molecular characterization of transgenic plants confirmed that their efficacy against H. armigera was a result of the integration of CHI4 from C. platycarpus. Docking analysis demonstrated effective substrate interaction as a possible reason for efficacy against pod borer in the chitinase from C. platycarpus. This was authenticated by successful overexpression and bioefficacy assessment against H. armigera in tobacco. The CHI4 gene from C. platycarpus can be useful in the mitigation of H. armegira in pigeonpea as well as in other crops. © 2021 Society of Chemical Industry.

Clinicopathological and prognostic significance of GPC3 in human breast cancer and its 3D structure prediction

Glypican-3 (GPC3) is a heparan sulfate proteoglycan that may function as a tumor suppressor in breast cancer (BC). To evaluate the prognostic value of GPC3 in BC, systematic analysis was performed in this study. To evaluate gene alteration during breast carcinogenesis, GPC3 expression was analyzed using the Oncomine, GENT, UALCAN, bcGenExMiner, and UCSC Xena databases. The prognostic role of GPC3 in BC was investigated using KM Plotter and PrognoScan databases. Promoter methylation status and heat map of GPC3 were determined using UALCAN and UCSC Xena. GPC3 expression was significantly downregulated in BC compared to that in normal tissues and was correlated with prognosis. However, estrogen receptor and progesterone receptor status were positively correlated with GPC3 expression, whereas basal-like status, triple-negative breast cancer status, and Scarff, Bloom, and Richardson grade criteria were negatively correlated with GPC3 expression. Further analysis indicated that GPC3 was correlated with Ras-association domain family 6 (RASSF6) expression in BC tissues. GPC3 may thus be considered a significant marker for predicting BC prognosis along with RASSF6. Comparative protein modeling of GPC3 was performed using a FASTA formatted sequence from NCBI, in Swiss model, and GPC3 was visualized using PyMOL tools. The GPC3 active site was identified with CASTp server. The Ramachandran plot in PROCHECK showed significant scores for the protein model in its most favored regions. ProSA server indicated the high accuracy of the protein model. However, large-scale and comprehensive studies are needed to clarify these results.

Screening of promising molecules against MurG as drug target in multi-drug-resistant-Acinetobacter baumannii - insights from comparative protein modeling, molecular docking and molecular dynamics simulation

The UDP-N-acetylglucosamine-N-acetylmuramyl-(pentapeptide) pyrophosphoryl-undecaprenol N-acetylglucosamine transferase (MurG) is located in plasma membrane which plays a crucial role for peptidoglycan biosynthesis in Gram-negative bacteria. Recently, this protein is considered as an important and unique drug target in Acinetobacter baumannii since it plays a key role during the synthesis of peptidoglycan as well as which is not found in Homo sapiens. In this study, initially we performed comparative protein modeling approach to predict the three-dimensional model of MurG based on crystal structure of UDP-N-acetylglucosamine-N-acetylmuramyl-(pentapeptide) pyrophosphoryl-undecaprenol N-acetylglucosamine transferase (PDB ID: 1F0K) from E.coli K12. MurG model has two important functional domains located in N and C- terminus which are separated by a deep cleft. Active site residues are located between two domains and they are Gly20, Arg170, Gly200, Ser201, Gln227, Phe254, Leu275, Thr276, and Glu279 which play essential role for the function of MurG. In order to inhibit the function of MurG, we employed the High Throughput Virtual Screening (HTVS) and docking techniques to identify the promising molecules which will further subjected into screening for computing their drug like and pharmacokinetic properties. From the HTVS, we identified 5279 molecules, among these, 12 were passed the drug-like and pharmacokinetic screening analysis. Based on the interaction analysis in terms of binding affinity, inhibition constant and intermolecular interactions, we selected four molecules for further MD simulation to understand the structural stability of protein-ligand complexes. All the analysis of MD simulation suggested that ZINC09186673 and ZINC09956120 are identified as most promising putative inhibitors for MurG protein in A. baumannii.Communicated by Ramaswamy H. Sarma

InterEvDock2: an expanded server for protein docking using evolutionary and biological information from homology models and multimeric inputs.

Computational protein docking is a powerful strategy to predict structures of protein-protein interactions and provides crucial insights for the functional characterization of macromolecular cross-talks. We previously developed InterEvDock, a server for ab initio protein docking based on rigid-body sampling followed by consensus scoring using physics-based and statistical potentials, including the InterEvScore function specifically developed to incorporate co-evolutionary information in docking. InterEvDock2 is a major evolution of InterEvDock which allows users to submit input sequences – not only structures – and multimeric inputs and to specify constraints for the pairwise docking process based on previous knowledge about the interaction. For this purpose, we added modules in InterEvDock2 for automatic template search and comparative modeling of the input proteins. The InterEvDock2 pipeline was benchmarked on 812 complexes for which unbound homology models of the two partners and co-evolutionary information are available in the PPI4DOCK database. InterEvDock2 identified a correct model among the top 10 consensus in 29% of these cases (compared to 15–24% for individual scoring functions) and at least one correct interface residue among 10 predicted in 91% of these cases. InterEvDock2 is thus a unique protein docking server, designed to be useful for the experimental biology community. The InterEvDock2 web interface is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/InterEvDock2/.

Heterozygous mutations affecting the protein kinase domain of CDK13 cause a syndromic form of developmental delay and intellectual disability

IntroductionRecent evidence has emerged linking mutations in CDK13 to syndromic congenital heart disease. We present here genetic and phenotypic data pertaining to 16 individuals with CDK13 mutations.MethodsPatients were investigated by...

Structural Investigation and In-silico Characterization of Plasmepsins from Plasmodium falciparum

Malaria is the one most important parasitic disease of humans, which affects approximately one hundred countries and threatens half of the world’s population. The Plasmodium aspartic protease called plasmepsins performs a vital role in providing nutrients to the malaria parasite, which make these proteins as an excellent drug target. In this study, we have carried out a comparative protein modeling, active site analysis and structural analysis of all ten plasmepsins from Plasmodium falciparum. In this report we have analyzed in-silico structure modeling and made efforts to characterize plasmepsins structure and further propose its functional information. The phylogenetic analysis and disulfide linkages indicate, plasmepsin I to IV and HAP have similar structure, function property. Whereas, plasmepsin IX to X and plasmepsin VI to VIII belong to a separate cluster. The integral membrane protein plasmepsin V has a functional characterization as compared to the others aspartic proteases from Plasmodium falciparum. The overall study summarizes the need of good model to understand the structure and function activity and to design potent small molecule inhibitors targeting all ten plasmepsins, specifically Plasmepsin V as important target.

Comparative genomic of the BAFF and BAFF-like genes and immune response to bacteria of miiuy croaker (Miichthys miiuy)

As a member of the tumor necrosis factor (TNF) family, B cell activating factor (BAFF), also known as TNF ligand superfamily member 13B (TNF13B), playing a critical role in enhancing immune responses. BAFF is a central cytokine for B-cell survival, proliferation, maturation and immunoglobulin secretion. In the present study, we describe the identification of the miiuy croaker BAFF (designed MmBAFF) and BAFF-like (designed MmBAFF-like) genes. The cDNA of MmBAFF contains an open reading frame (ORF) of 795 nucleotides that are translated into a predicted 264 amino acids. The ORF of MmBAFF-like consists of 705 bases encoding 234 amino acids. Amino acid sequence comparison indicated that MmBAFF and MmBAFF-like possessed the TNF signatures, a predicted transmembrane domain, three conserved cysteine residues and a putative furin protease cleavage site, which were the typical characteristics of TNF gene in mammals and birds. The predicted three-dimensional (3D) structure of the MmBAFF and MmBAFF-like monomer analyzed by comparative protein modeling revealed that they were very similar to human counterpart. Comparative genomic analysis revealed that the locations of MmBAFF and MmBAFF-like genes are conserved among the bony fish. Phylogenetic analysis shows the MmBAFF is most closely related to other teleost BAFFs with the highest similarity to Epinephelus awoara. And BAFF-like cluster get together first to BAFF cluster than three closely related TNF superfamily (TNFSF) members. Real-time quantitative PCR analysis shows the MmBAFF and MmBAFF-like transcripts are expressed in a wide range of tissues with the highest expression in skin and lymphoid tissue spleen. Upon induction by Vibrio anguillarum, their expressions are significantly upregulated in liver, spleen and kidney as compared to phosphate-buffered saline injected control fish. The association of increased BAFF expression after bacterial infection suggests that it plays a potentially important role in immune system of fish.

A Basic Protein Comparative Three-Dimensional Modeling Methodological Workflow Theory and Practice.

When working with proteins and studying its properties, it is crucial to have access to the three-dimensional structure of the molecule. If experimentally solved structures are not available, comparative modeling techniques can be used to generate useful protein models to subsidize structure-based research projects. In recent years, with Bioinformatics becoming the basis for the study of protein structures, there is a crescent need for the exposure of details about the algorithms behind the softwares and servers, as well as a need for protocols to guide in silico predictive experiments. In this article, we explore different steps of the comparative modeling technique, such as template identification, sequence alignment, generation of candidate structures and quality assessment, its peculiarities and theoretical description. We then present a practical step-by-step workflow, to support the Biologist on the in silico generation of protein structures. Finally, we explore further steps on comparative modeling, presenting perspectives to the study of protein structures through Bioinformatics. We trust that this is a thorough guide for beginners that wish to work on the comparative modeling of proteins.

Can natural proteins designed with 'inverted' peptide sequences adopt native-like protein folds?

We have carried out a systematic computational analysis on a representative dataset of proteins of known three-dimensional structure, in order to evaluate whether it would possible to ‘swap’ certain short peptide sequences in naturally occurring proteins with their corresponding ‘inverted’ peptides and generate ‘artificial’ proteins that are predicted to retain native-like protein fold. The analysis of 3,967 representative proteins from the Protein Data Bank revealed 102,677 unique identical inverted peptide sequence pairs that vary in sequence length between 5–12 and 18 amino acid residues. Our analysis illustrates with examples that such ‘artificial’ proteins may be generated by identifying peptides with ‘similar structural environment’ and by using comparative protein modeling and validation studies. Our analysis suggests that natural proteins may be tolerant to accommodating such peptides.

Control of KirBac3.1 Potassium Channel Gating at the Interface between Cytoplasmic Domains

KirBac channels are prokaryotic homologs of mammalian inwardly rectifying potassium (Kir) channels, and recent structures of KirBac3.1 have provided important insights into the structural basis of gating in Kir channels. In this study, we demonstrate that KirBac3.1 channel activity is strongly pH-dependent, and we used x-ray crystallography to determine the structural changes that arise from an activatory mutation (S205L) located in the cytoplasmic domain (CTD). This mutation stabilizes a novel energetically favorable open conformation in which changes at the intersubunit interface in the CTD also alter the electrostatic potential of the inner cytoplasmic cavity. These results provide a structural explanation for the activatory effect of this mutation and provide a greater insight into the role of the CTD in Kir channel gating.

Comparative modeling of proteins: A method for engaging students' interest in bioinformatics tools

The huge increase in data being produced in the genomic era has produced a need to incorporate computers into the research process. Sequence generation, its subsequent storage, interpretation, and analysis are now entirely computer-dependent tasks. Universities from all over the world have been challenged to seek a way of encouraging students to incorporate computational and bioinformatics skills since undergraduation in order to understand biological processes. The aim of this article is to report the experience of awakening students' interest in bioinformatics tools during a course focused on comparative modeling of proteins. The authors start by giving a full description of the course environmental context and students' backgrounds. Then they detail each class and present a general overview of the protein modeling protocol. The positive and negative aspects of the course are also reported, and some of the results generated in class and in projects outside the classroom are discussed. In the last section of the article, general perspectives about the course from students' point of view are given. This work can serve as a guide for professors who teach subjects for which bioinformatics tools are useful and for universities that plan to incorporate bioinformatics into the curriculum.

HomoSAR: Bridging comparative protein modeling with quantitative structural activity relationship to design new peptides

Peptides play significant roles in the biological world. To optimize activity for a specific therapeutic target, peptide library synthesis is inevitable; which is a time consuming and expensive. Computational approaches provide a promising way to simply elucidate the structural basis in the design of new peptides. Earlier, we proposed a novel methodology termed HomoSAR to gain insight into the structure activity relationships underlying peptides. Based on an integrated approach, HomoSAR uses the principles of homology modeling in conjunction with the quantitative structural activity relationship formalism to predict and design new peptide sequences with the optimum activity. In the present study, we establish that the HomoSAR methodology can be universally applied to all classes of peptides irrespective of sequence length by studying HomoSAR on three peptide datasets viz., angiotensin-converting enzyme inhibitory peptides, CAMEL-s antibiotic peptides, and hAmphiphysin-1 SH3 domain binding peptides, using a set of descriptors related to the hydrophobic, steric, and electronic properties of the 20 natural amino acids. Models generated for all three datasets have statistically significant correlation coefficients (r(2)) and predictive r2 (r(pred)2) and cross validated coefficient ( q(LOO)2). The daintiness of this technique lies in its simplicity and ability to extract all the information contained in the peptides to elucidate the underlying structure activity relationships. The difficulties of correlating both sequence diversity and variation in length of the peptides with their biological activity can be addressed. The study has been able to identify the preferred or detrimental nature of amino acids at specific positions in the peptide sequences.