Atomic-level protein structure refinement using fragment-guided molecular dynamics conformation sampling

Heat Shock Protein 90 Associates with the Per-Arnt-Sim Domain of Heme-free Soluble Guanylate Cyclase: IMplications for Enzyme Maturation

Heat shock protein 90 (hsp90) drives heme insertion into the β1 subunit of soluble guanylate cyclase (sGC) β1, which enables it to associate with a partner sGCα1 subunit and mature into a nitric oxide (NO)-responsive active form. We utilized fluorescence polarization measurements and hydrogen-deuterium exchange mass spectrometry to define molecular interactions between the specific human isoforms hsp90β and apo-sGCβ1. hsp90β and its isolated M domain, but not its isolated N and C domains, bind with low micromolar affinity to a heme-free, truncated version of sGCβ1 (sGCβ1(1-359)-H105F). Surprisingly, hsp90β and its M domain bound to the Per-Arnt-Sim (PAS) domain of apo-sGC-β1(1-359), which lies adjacent to its heme-binding (H-NOX) domain. The interaction specifically involved solvent-exposed regions in the hsp90β M domain that are largely distinct from sites utilized by other hsp90 clients. The interaction strongly protected two regions of the sGCβ1 PAS domain and caused local structural relaxation in other regions, including a PAS dimerization interface and a segment in the H-NOX domain. Our results suggest a means by which the hsp90β interaction could prevent apo-sGCβ1 from associating with its partner sGCα1 subunit while enabling structural changes to assist heme insertion into the H-NOX domain. This mechanism would parallel that in other clients like the aryl hydrocarbon receptor and HIF1α, which also interact with hsp90 through their PAS domains to control protein partner and small ligand binding interactions.

The Analysis of Pendolino (peo) Mutants Reveals Differences in the Fusigenic Potential among Drosophila Telomeres

Drosophila telomeres are sequence-independent structures that are maintained by transposition to chromosome ends of three specialized retroelements (HeT-A, TART and TAHRE; collectively designated as HTT) rather than telomerase activity. Fly telomeres are protected by the terminin complex (HOAP-HipHop-Moi-Ver) that localizes and functions exclusively at telomeres and by non-terminin proteins that do not serve telomere-specific functions. Although all Drosophila telomeres terminate with HTT arrays and are capped by terminin, they differ in the type of subtelomeric chromatin; the Y, XR, and 4L HTT are juxtaposed to constitutive heterochromatin, while the XL, 2L, 2R, 3L and 3R HTT are linked to the TAS repetitive sequences; the 4R HTT is associated with a chromatin that has features common to both euchromatin and heterochromatin. Here we show that mutations in pendolino (peo) cause telomeric fusions (TFs). The analysis of several peo mutant combinations showed that these TFs preferentially involve the Y, XR and 4th chromosome telomeres, a TF pattern never observed in the other 10 telomere-capping mutants so far characterized. peo encodes a non-terminin protein homologous to the E2 variant ubiquitin-conjugating enzymes. The Peo protein directly interacts with the terminin components, but peo mutations do not affect telomeric localization of HOAP, Moi, Ver and HP1a, suggesting that the peo-dependent telomere fusion phenotype is not due to loss of terminin from chromosome ends. peo mutants are also defective in DNA replication and PCNA recruitment. However, our results suggest that general defects in DNA replication are unable to induce TFs in Drosophila cells. We thus hypothesize that DNA replication in Peo-depleted cells results in specific fusigenic lesions concentrated in heterochromatin-associated telomeres. Alternatively, it is possible that Peo plays a dual function being independently required for DNA replication and telomere capping.

Telomeres are specialized structures that protect chromosome ends from incomplete replication, degradation and end-to-end fusion. Abnormalities in telomere structure or maintenance can promote a variety of human diseases including premature aging and cancer. Although all human telomeres contain the same DNA sequences, they differ from each other in the subtelomeric regions or subtelomeres. Recent work has shown that human subtelomeres control telomere replication and that abnormalities in these structures can lead to localized chromosome instability and disease. However, the relationships between subtelomeres and telomeres are currently poorly understood. Here, we have addressed this problem using the fruit fly Drosophila melanogaster as model system. Drosophila subtelomers are very different from each other as they contain different types of chromatin. We have found that mutations in a gene we called pendolino (peo) cause telomeric fusions (TFs) and that these fusions preferentially involve the telomeres associated with a tightly packed form of chromatin called heterochromatin. Interestingly, none of the 10 mutants with TFs so far described in Drosophila shows the pattern of TFs observed in peo mutants. Thus, our data provide the first demonstration that subtelomeres can affect telomere fusion. We believe that these results will stimulate further studies on the role of subtelomeres in the maintenance of genome stability.

AKTIP/Ft1, a New Shelterin-Interacting Factor Required for Telomere Maintenance

Telomeres are nucleoprotein complexes that protect the ends of linear chromosomes from incomplete replication, degradation and detection as DNA breaks. Mammalian telomeres are protected by shelterin, a multiprotein complex that binds the TTAGGG telomeric repeats and recruits a series of additional factors that are essential for telomere function. Although many shelterin-associated proteins have been so far identified, the inventory of shelterin-interacting factors required for telomere maintenance is still largely incomplete. Here, we characterize AKTIP/Ft1 (human AKTIP and mouse Ft1 are orthologous), a novel mammalian shelterin-bound factor identified on the basis of its homology with the Drosophila telomere protein Pendolino. AKTIP/Ft1 shares homology with the E2 variant ubiquitin-conjugating (UEV) enzymes and has been previously implicated in the control of apoptosis and in vesicle trafficking. RNAi-mediated depletion of AKTIP results in formation of telomere dysfunction foci (TIFs). Consistent with these results, AKTIP interacts with telomeric DNA and binds the shelterin components TRF1 and TRF2 both in vivo and in vitro. Analysis of AKTIP- depleted human primary fibroblasts showed that they are defective in PCNA recruiting and arrest in the S phase due to the activation of the intra S checkpoint. Accordingly, AKTIP physically interacts with PCNA and the RPA70 DNA replication factor. Ft1-depleted p53-/- MEFs did not arrest in the S phase but displayed significant increases in multiple telomeric signals (MTS) and sister telomere associations (STAs), two hallmarks of defective telomere replication. In addition, we found an epistatic relation for MST formation between Ft1 and TRF1, which has been previously shown to be required for replication fork progression through telomeric DNA. Ch-IP experiments further suggested that in AKTIP-depleted cells undergoing the S phase, TRF1 is less tightly bound to telomeric DNA than in controls. Thus, our results collectively suggest that AKTIP/Ft1 works in concert with TRF1 to facilitate telomeric DNA replication.

Bioinformatical analysis of the sequences, structures and functions of fungal polyketide synthase product template domains

The product template (PT) domains, specifically in fungal non-reducing polyketide synthases (NR-PKSs), mediate the regioselective cyclization of polyketides dominating the final structures. However, up to date, the systematic knowledge about PT domains has been insufficient. In present study, the relationships between sequences, structures and functions of the PT domains were analyzed with 661 NR-PKS sequences. Based on the phylogenetic analysis, the PT domains were classified into prominent eight groups (I–VIII) corresponding with the representative compounds and cyclization regioselectivity (C2-C7, C4-C9, and C6-C11). Most of the cavity lining residue (CLR) sites in all groups were common, while the regional CLR site mutations resulted in the appearance of finger-like regions with different orientation. The cavity volumes and shapes, even the catalytic dyad positions of PT domains in different groups were corresponding with characteristic cyclization regioselectivity and compound sizes. The conservative residues in PT sequences were responsible for the cyclization functions and the evolution of the key residues resulted in the differentiations of cyclization functions. The above findings may help to better understand the cyclization mechanisms of PT domains and even predict the structural types of the aromatic polyketide products.

I-TASSER server: new development for protein structure and function predictions

The I-TASSER server (http://zhanglab.ccmb.med.umich.edu/I-TASSER) is an online resource for automated protein structure prediction and structure-based function annotation. In I-TASSER, structural templates are first recognized from the PDB using multiple threading alignment approaches. Full-length structure models are then constructed by iterative fragment assembly simulations. The functional insights are finally derived by matching the predicted structure models with known proteins in the function databases. Although the server has been widely used for various biological and biomedical investigations, numerous comments and suggestions have been reported from the user community. In this article, we summarize recent developments on the I-TASSER server, which were designed to address the requirements from the user community and to increase the accuracy of modeling predictions. Focuses have been made on the introduction of new methods for atomic-level structure refinement, local structure quality estimation and biological function annotations. We expect that these new developments will improve the quality of the I-TASSER server and further facilitate its use by the community for high-resolution structure and function prediction.

Significant Refinement of Protein Structure Models Using a Residue-Specific Force Field

An important application of all-atom explicit-solvent molecular dynamics (MD) simulations is the refinement of protein structures from low-resolution experiments or template-based modeling. A critical requirement is that the native structure is stable with the force field. We have applied a recently developed residue-specific force field, RSFF1, to a set of 30 refinement targets from recent CASP experiments. Starting from their experimental structures, 1.0 μs unrestrained simulations at 298 K retain most of the native structures quite well except for a few flexible terminals and long internal loops. Starting from each homology model, a 150 ns MD simulation at 380 K generates the best RMSD improvement of 0.85 Å on average. The structural improvements roughly correlate with the RMSD of the initial homology models, indicating possible consistent structure refinement. Finally, targets TR614 and TR624 have been subjected to long-time replica-exchange MD simulations. Significant structural improvements are generated, with RMSD of 1.91 and 1.36 Å with respect to their crystal structures. Thus, it is possible to achieve realistic refinement of protein structure models to near-experimental accuracy, using accurate force field with sufficient conformational sampling.

Looking into laminin receptor: critical discussion regarding the non-integrin 37/67-kDa laminin receptor/RPSA protein

The 37/67-kDa laminin receptor (LAMR/RPSA) was originally identified as a 67-kDa binding protein for laminin, an extracellular matrix glycoprotein that provides cellular adhesion to the basement membrane. LAMR has evolutionary origins, however, as a 37-kDa RPS2 family ribosomal component. Expressed in all domains of life, RPS2 proteins have been shown to have remarkably diverse physiological roles that vary across species. Contributing to laminin binding, ribosome biogenesis, cytoskeletal organization, and nuclear functions, this protein governs critical cellular processes including growth, survival, migration, protein synthesis, development, and differentiation. Unsurprisingly given its purview, LAMR has been associated with metastatic cancer, neurodegenerative disease and developmental abnormalities. Functioning in a receptor capacity, this protein also confers susceptibility to bacterial and viral infection. LAMR is clearly a molecule of consequence in human disease, directly mediating pathological events that make it a prime target for therapeutic interventions. Despite decades of research, there are still a large number of open questions regarding the cellular biology of LAMR, the nature of its ability to bind laminin, the function of its intrinsically disordered C-terminal region and its conversion from 37 to 67 kDa. This review attempts to convey an in-depth description of the complexity surrounding this multifaceted protein across functional, structural and pathological aspects.

Cell penetrable human scFv specific to middle domain of matrix protein-1 protects mice from lethal influenza

A new anti-influenza remedy that can tolerate the virus antigenic variation is needed. Influenza virus matrix protein-1 (M1) is highly conserved and pivotal for the virus replication cycle: virus uncoating, assembly and budding. An agent that blocks the M1 functions should be an effective anti-influenza agent. In this study, human scFv that bound to recombinant M1 middle domain (MD) and native M1 of A/H5N1 was produced. Phage mimotope search and computerized molecular docking revealed that the scFv bound to the MD conformational epitope formed by juxtaposed helices 7 and 9 of the M1. The scFv was linked molecularly to a cell penetrable peptide, penetratin (PEN). The PEN-scFv (transbody), when used to treat the cells pre-infected with the heterologous clade/subclade A/H5N1 reduced the viral mRNA intracellularly and in the cell culture fluids. The transbody mitigated symptom severity and lung histopathology of the H5N1 infected mice and caused reduction of virus antigen in the tissues as well as extricated the animals from the lethal challenge in a dose dependent manner. The transbody specific to the M1 MD, either alone or in combination with the cognate human scFvs specific to other influenza virus proteins, should be an effective, safe and mutation tolerable anti-influenza agent.

Protein structure prediction provides comparable performance to crystallographic structures in docking-based virtual screening

Structure based virtual screening has largely been limited to protein targets for which either an experimental structure is available or a strongly homologous template exists so that a high-resolution model can be constructed. The performance of state of the art protein structure predictions in virtual screening in systems where only weakly homologous templates are available is largely untested. Using the challenging DUD database of structural decoys, we show here that even using templates with only weak sequence homology (<30% sequence identity) structural models can be constructed by I-TASSER which achieve comparable enrichment rates to using the experimental bound crystal structure in the majority of the cases studied. For 65% of the targets, the I-TASSER models, which are constructed essentially in the apo conformations, reached 70% of the virtual screening performance of using the holo-crystal structures. A correlation was observed between the success of I-TASSER in modeling the global fold and local structures in the binding pockets of the proteins versus the relative success in virtual screening. The virtual screening performance can be further improved by the recognition of chemical features of the ligand compounds. These results suggest that the combination of structure-based docking and advanced protein structure modeling methods should be a valuable approach to the large-scale drug screening and discovery studies, especially for the proteins lacking crystallographic structures.

Protein surface charge of trypsinogen changes its activation pattern

Background

Trypsinogen is the inactive precursor of trypsin, a serine protease that cleaves proteins and peptides after arginine and lysine residues. In this study, human trypsinogen was used as a model protein to study the influence of electrostatic forces on protein–protein interactions. Trypsinogen is active only after its eight-amino-acid-long activation peptide has been cleaved off by another protease, enteropeptidase. Trypsinogen can also be autoactivated without the involvement of enteropeptidase. This autoactivation process can occur if a trypsinogen molecule is activated by another trypsin molecule and therefore is based on a protein–protein interaction.

Results

Based on a rational protein design based on autoactivation-defective guinea pig trypsinogen, several amino acid residues, all located far away from the active site, were changed to modify the surface charge of human trypsinogen. The influence of the surface charge on the activation pattern of trypsinogen was investigated. The autoactivation properties of mutant trypsinogen were characterized in comparison to the recombinant wild-type enzyme. Surface-charged trypsinogen showed practically no autoactivation compared to the wild-type but could still be activated by enteropeptidase to the fully active trypsin. The kinetic parameters of surface-charged trypsinogen were comparable to the recombinant wild-type enzyme.

Conclusion

The variant with a modified surface charge compared to the wild-type enzyme showed a complete different activation pattern. Our study provides an example how directed modification of the protein surface charge can be utilized for the regulation of functional protein–protein interactions, as shown here for human trypsinogen.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-014-0109-5) contains supplementary material, which is available to authorized users.

Comparative genomics of closely related Salmonella enterica serovar Typhi strains reveals genome dynamics and the acquisition of novel pathogenic elements

Background

Typhoid fever is an infectious disease of global importance that is caused by Salmonella enterica subsp. enterica serovar Typhi (S. Typhi). This disease causes an estimated 200,000 deaths per year and remains a serious global health threat. S. Typhi is strictly a human pathogen, and some recovered individuals become long-term carriers who continue to shed the bacteria in their faeces, thus becoming main reservoirs of infection.

Results

A comparative genomics analysis combined with a phylogenomic analysis revealed that the strains from the outbreak and carrier were closely related with microvariations and possibly derived from a common ancestor. Additionally, the comparative genomics analysis with all of the other completely sequenced S. Typhi genomes revealed that strains BL196 and CR0044 exhibit unusual genomic variations despite S. Typhi being generally regarded as highly clonal. The two genomes shared distinct chromosomal architectures and uncommon genome features; notably, the presence of a ~10 kb novel genomic island containing uncharacterised virulence-related genes, and zot in particular. Variations were also detected in the T6SS system and genes that were related to SPI-10, insertion sequences, CRISPRs and nsSNPs among the studied genomes. Interestingly, the carrier strain CR0044 harboured far more genetic polymorphisms (83% mutant nsSNPs) compared with the closely related BL196 outbreak strain. Notably, the two highly related virulence-determinant genes, rpoS and tviE, were mutated in strains BL196 and CR0044, respectively, which revealed that the mutation in rpoS is stabilising, while that in tviE is destabilising. These microvariations provide novel insight into the optimisation of genes by the pathogens. However, the sporadic strain was found to be far more conserved compared with the others.

Conclusions

The uncommon genomic variations in the two closely related BL196 and CR0044 strains suggests that S. Typhi is more diverse than previously thought. Our study has demonstrated that the pathogen is continually acquiring new genes through horizontal gene transfer in the process of host adaptation, providing novel insight into its unusual genomic dynamics. The understanding of these strains and virulence factors, and particularly the strain that is associated with the large outbreak and the less studied asymptomatic Typhi carrier in the population, will have important impact on disease control.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1007) contains supplementary material, which is available to authorized users.

Advances in GPCR modeling evaluated by the GPCR Dock 2013 assessment: meeting new challenges

Despite tremendous successes of GPCR crystallography, the receptors with available structures represent only a small fraction of human GPCRs. An important role of the modeling community is to maximize structural insights for the remaining receptors and complexes. The community-wide GPCR Dock assessment was established to stimulate and monitor the progress in molecular modeling and ligand docking for GPCRs. The four targets in the present third assessment round presented new and diverse challenges for modelers, including prediction of allosteric ligand interaction and activation states in 5-hydroxytryptamine receptors 1B and 2B, and modeling by extremely distant homology for smoothened receptor. Forty-four modeling groups participated in the assessment. State-of-the-art modeling approaches achieved close-to-experimental accuracy for small rigid orthosteric ligands and models built by close homology, and they correctly predicted protein fold for distant homology targets. Predictions of long loops and GPCR activation states remain unsolved problems.

Evaluation of unrestrained replica-exchange simulations using dynamic walkers in temperature space for protein structure refinement

A central problem of computational structural biology is the refinement of modeled protein structures taken from either comparative modeling or knowledge-based methods. Simulations are commonly used to achieve higher resolution of the structures at the all-atom level, yet methodologies that consistently yield accurate results remain elusive. In this work, we provide an assessment of an adaptive temperature-based replica exchange simulation method where the temperature clients dynamically walk in temperature space to enrich their population and exchanges near steep energetic barriers. This approach is compared to earlier work of applying the conventional method of static temperature clients to refine a dataset of conformational decoys. Our results show that, while an adaptive method has many theoretical advantages over a static distribution of client temperatures, only limited improvement was gained from this strategy in excursions of the downhill refinement regime leading to an increase in the fraction of native contacts. To illustrate the sampling differences between the two simulation methods, energy landscapes are presented along with their temperature client profiles.

Ophiophagus hannah venom: proteome, components bound by Naja kaouthia antivenin and neutralization by N. kaouthia neurotoxin-specific human ScFv

Venomous snakebites are an important health problem in tropical and subtropical countries. King cobra (Ophiophagushannah) is the largest venomous snake found in South and Southeast Asia. In this study, the O. hannah venom proteome and the venom components cross-reactive to N. kaouthia monospecific antivenin were studied. O. hannah venom consisted of 14 different protein families, including three finger toxins, phospholipases, cysteine-rich secretory proteins, cobra venom factor, muscarinic toxin, L-amino acid oxidase, hypothetical proteins, low cysteine protein, phosphodiesterase, proteases, vespryn toxin, Kunitz, growth factor activators and others (coagulation factor, endonuclease, 5’-nucleotidase). N. kaouthia antivenin recognized several functionally different O. hannah venom proteins and mediated paratherapeutic efficacy by rescuing the O. hannah envenomed mice from lethality. An engineered human ScFv specific to N. kaouthia long neurotoxin (NkLN-HuScFv) cross-neutralized the O. hannah venom and extricated the O. hannah envenomed mice from death in a dose escalation manner. Homology modeling and molecular docking revealed that NkLN-HuScFv interacted with residues in loops 2 and 3 of the neurotoxins of both snake species, which are important for neuronal acetylcholine receptor binding. The data of this study are useful for snakebite treatment when and where the polyspecific antivenin is not available. Because the supply of horse-derived antivenin is limited and the preparation may cause some adverse effects in recipients, a cocktail of recombinant human ScFvs for various toxic venom components shared by different venomous snakes, exemplified by the in vitro produced NkLN-HuScFv in this study, should contribute to a possible future route for an improved alternative to the antivenins.

Protein structure refinement of CASP target proteins using GNEIMO torsional dynamics method

A longstanding challenge in using computational methods for protein structure prediction is the refinement of low-resolution structural models derived from comparative modeling methods into highly accurate atomistic models useful for detailed structural studies. Previously, we have developed and demonstrated the utility of the internal coordinate molecular dynamics (MD) technique, generalized Newton-Euler inverse mass operator (GNEIMO), for refinement of small proteins. Using GNEIMO, the high-frequency degrees of freedom are frozen and the protein is modeled as a collection of rigid clusters connected by torsional hinges. This physical model allows larger integration time steps and focuses the conformational search in the low frequency torsional degrees of freedom. Here, we have applied GNEIMO with temperature replica exchange to refine low-resolution protein models of 30 proteins taken from the continuous assessment of structure prediction (CASP) competition. We have shown that GNEIMO torsional MD method leads to refinement of up to 1.3 Å in the root-mean-square deviation in coordinates for 30 CASP target proteins without using any experimental data as restraints in performing the GNEIMO simulations. This is in contrast with the unconstrained all-atom Cartesian MD method performed under the same conditions, where refinement requires the use of restraints during the simulations.

Physics-based protein structure refinement through multiple molecular dynamics trajectories and structure averaging

We used molecular dynamics (MD) simulations for structure refinement of Critical Assessment of Techniques for Protein Structure Prediction 10 (CASP10) targets. Refinement was achieved by selecting structures from the MD-based ensembles followed by structural averaging. The overall performance of this method in CASP10 is described, and specific aspects are analyzed in detail to provide insight into key components. In particular, the use of different restraint types, sampling from multiple short simulations versus a single long simulation, the success of a quality assessment criterion, the application of scoring versus averaging, and the impact of a final refinement step are discussed in detail.

Histatin 5 binds to Porphyromonas gingivalis hemagglutinin B (HagB) and alters HagB-induced chemokine responses

Histatins are human salivary gland peptides with anti-microbial and anti-inflammatory activities. In this study, we hypothesized that histatin 5 binds to Porphyromonas gingivalis hemagglutinin B (HagB) and attenuates HagB-induced chemokine responses in human myeloid dendritic cells. Histatin 5 bound to immobilized HagB in a surface plasmon resonance (SPR) spectroscopy-based biosensor system. SPR spectroscopy kinetic and equilibrium analyses, protein microarray studies, and I-TASSER structural modeling studies all demonstrated two histatin 5 binding sites on HagB. One site had a stronger affinity with a K_D1 of 1.9 μM and one site had a weaker affinity with a K_D2 of 60.0 μM. Binding has biological implications and predictive modeling studies and exposure of dendritic cells both demonstrated that 20.0 μM histatin 5 attenuated (p < 0.05) 0.02 μM HagB-induced CCL3/MIP-1α, CCL4/MIP-1β, and TNFα responses. Thus histatin 5 is capable of attenuating chemokine responses, which may help control oral inflammation.

Evaluation of predictions in the CASP10 model refinement category

Here we report on the assessment results of the third experiment to evaluate the state of the art in protein model refinement, where participants were invited to improve the accuracy of initial protein models for 27 targets. Using an array of complementary evaluation measures, we find that five groups performed better than the naïve (null) method—a marked improvement over CASP9, although only three were significantly better. The leading groups also demonstrated the ability to consistently improve both backbone and side chain positioning, while other groups reliably enhanced other aspects of protein physicality. The top-ranked group succeeded in improving the backbone conformation in almost 90% of targets, suggesting a strategy that for the first time in CASP refinement is successful in a clear majority of cases. A number of issues remain unsolved: the majority of groups still fail to improve the quality of the starting models; even successful groups are only able to make modest improvements; and no prediction is more similar to the native structure than to the starting model. Successful refinement attempts also often go unrecognized, as suggested by the relatively larger improvements when predictions not submitted as model 1 are also considered. Proteins 2014; 82(Suppl 2):98–111.

Targeting lipid esterases in mycobacteria grown under different physiological conditions using activity-based profiling with tetrahydrolipstatin (THL)

Tetrahydrolipstatin (THL) is bactericidal but its precise target spectrum is poorly characterized. Here, we used a THL analog and activity-based protein profiling to identify target proteins after enrichment from whole cell lysates of Mycobacterium bovis Bacillus Calmette-Guérin cultured under replicating and non-replicating conditions. THL targets α/β-hydrolases, including many lipid esterases (LipD, G, H, I, M, N, O, V, W, and TesA). Target protein concentrations and total esterase activity correlated inversely with cellular triacylglycerol upon entry into and exit from non-replicating conditions. Cellular overexpression of lipH and tesA led to decreased THL susceptibility thus providing functional validation. Our results define the target spectrum of THL in a biological species with particularly diverse lipid metabolic pathways. We furthermore derive a conceptual approach that demonstrates the use of such THL probes for the characterization of substrate recognition by lipases and related enzymes.

Princeton_TIGRESS: protein geometry refinement using simulations and support vector machines

Protein structure refinement aims to perform a set of operations given a predicted structure to improve model quality and accuracy with respect to the native in a blind fashion. Despite the numerous computational approaches to the protein refinement problem reported in the previous three CASPs, an overwhelming majority of methods degrade models rather than improve them. We initially developed a method tested using blind predictions during CASP10 which was officially ranked in 5th place among all methods in the refinement category. Here, we present Princeton_TIGRESS, which when benchmarked on all CASP 7,8,9, and 10 refinement targets, simultaneously increased GDT_TS 76% of the time with an average improvement of 0.83 GDT_TS points per structure. The method was additionally benchmarked on models produced by top performing three-dimensional structure prediction servers during CASP10. The robustness of the Princeton_TIGRESS protocol was also tested for different random seeds. We make the Princeton_TIGRESS refinement protocol freely available as a web server at http://atlas.princeton.edu/refinement. Using this protocol, one can consistently refine a prediction to help bridge the gap between a predicted structure and the actual native structure.

Protein folding and de novo protein design for biotechnological applications

In the postgenomic era, the medical/biological fields are advancing faster than ever. However, before the power of full-genome sequencing can be fully realized, the connection between amino acid sequence and protein structure, known as the protein folding problem, needs to be elucidated. The protein folding problem remains elusive, with significant difficulties still arising when modeling amino acid sequences lacking an identifiable template. Understanding protein folding will allow for unforeseen advances in protein design; often referred to as the inverse protein folding problem. Despite challenges in protein folding, de novo protein design has recently demonstrated significant success via computational techniques. We review advances and challenges in protein structure prediction and de novo protein design, and highlight their interplay in successful biotechnological applications.

An evolution-based approach to De Novo protein design and case study on Mycobacterium tuberculosis

Computational protein design is a reverse procedure of protein folding and structure prediction, where constructing structures from evolutionarily related proteins has been demonstrated to be the most reliable method for protein 3-dimensional structure prediction. Following this spirit, we developed a novel method to design new protein sequences based on evolutionarily related protein families. For a given target structure, a set of proteins having similar fold are identified from the PDB library by structural alignments. A structural profile is then constructed from the protein templates and used to guide the conformational search of amino acid sequence space, where physicochemical packing is accommodated by single-sequence based solvation, torsion angle, and secondary structure predictions. The method was tested on a computational folding experiment based on a large set of 87 protein structures covering different fold classes, which showed that the evolution-based design significantly enhances the foldability and biological functionality of the designed sequences compared to the traditional physics-based force field methods. Without using homologous proteins, the designed sequences can be folded with an average root-mean-square-deviation of 2.1 Å to the target. As a case study, the method is extended to redesign all 243 structurally resolved proteins in the pathogenic bacteria Mycobacterium tuberculosis, which is the second leading cause of death from infectious disease. On a smaller scale, five sequences were randomly selected from the design pool and subjected to experimental validation. The results showed that all the designed proteins are soluble with distinct secondary structure and three have well ordered tertiary structure, as demonstrated by circular dichroism and NMR spectroscopy. Together, these results demonstrate a new avenue in computational protein design that uses knowledge of evolutionary conservation from protein structural families to engineer new protein molecules of improved fold stability and biological functionality.

The goal of computational protein design is to create new protein sequences of desirable structure and biological function. Most protein design methods are developed to search for sequences with the lowest free-energy based on physics-based force fields following Anfinsen's thermodynamic hypothesis. A major obstacle of such approaches is the inaccuracy of the force-field design, which cannot accurately describe atomic interactions or correctly recognize protein folds. We propose a novel method which uses evolutionary information, in the form of sequence profiles from structure families, to guide the sequence design. Since sequence profiles are generally more accurate than physics-based potentials in protein fold recognition, a unique advantage lies on that it targets the design procedure to a family of protein sequence profiles to enhance the robustness of designed sequences. The method was tested on 87 proteins and the designed sequences can be folded by I-TASSER to models with an average RMSD 2.1 Å. As a case study of large-scale application, the method is extended to redesign all structurally resolved proteins in the human pathogenic bacteria, Mycobacterium tuberculosis. Five sequences varying in fold and sizes were characterized by circular dichroism and NMR spectroscopy experiments and three were shown to have ordered tertiary structure.

Molecular dynamics simulations on the Tre1 G protein-coupled receptor: exploring the role of the arginine of the NRY motif in Tre1 structure

BACKGROUND

The arginine of the D/E/NRY motif in Rhodopsin family G protein-coupled receptors (GPCRs) is conserved in 96% of these proteins. In some GPCRs, this arginine in transmembrane 3 can form a salt bridge with an aspartic acid or glutamic acid in transmembrane 6. The Drosophila melanogaster GPCR Trapped in endoderm-1 (Tre1) is required for normal primordial germ cell migration. In a mutant form of the protein, Tre1sctt, eight amino acids RYILIACH are missing, resulting in a severe disruption of primordial germ cell development. The impact of the loss of these amino acids on Tre1 structure is unknown. Since the missing amino acids in Tre1sctt include the arginine that is part of the D/E/NRY motif in Tre1, molecular dynamics simulations were performed to explore the hypothesis that these amino acids are involved in salt bridge formation and help maintain Tre1 structure.

RESULTS

Structural predictions of wild type Tre1 (Tre1+) and Tre1sctt were subjected to over 250 ns of molecular dynamics simulations. The ability of the model systems to form a salt bridge between the arginine of the D/E/NRY motif and an aspartic acid residue in transmembrane 6 was analyzed. The results indicate that a stable salt bridge can form in the Tre1+ systems and a weak salt bridge or no salt bridge, using an alternative arginine, is likely in the Tre1sctt systems.

CONCLUSIONS

The weak salt bridge or lack of a salt bridge in the Tre1sctt systems could be one possible explanation for the disrupted function of Tre1sctt in primordial germ cell migration. These results provide a framework for studying the importance of the arginine of the D/E/NRY motif in the structure and function of other GPCRs that are involved in cell migration, such as CXCR4 in the mouse, zebrafish, and chicken.

Interplay of I-TASSER and QUARK for template-based and ab initio protein structure prediction in CASP10

We develop and test a new pipeline in CASP10 to predict protein structures based on an interplay of I-TASSER and QUARK for both free-modeling (FM) and template-based modeling (TBM) targets. The most noteworthy observation is that sorting through the threading template pool using the QUARK-based ab initio models as probes allows the detection of distant-homology templates which might be ignored by the traditional sequence profile-based threading alignment algorithms. Further template assembly refinement by I-TASSER resulted in successful folding of two medium-sized FM targets with >150 residues. For TBM, the multiple threading alignments from LOMETS are, for the first time, incorporated into the ab initio QUARK simulations, which were further refined by I-TASSER assembly refinement. Compared with the traditional threading assembly refinement procedures, the inclusion of the threading-constrained ab initio folding models can consistently improve the quality of the full-length models as assessed by the GDT-HA and hydrogen-bonding scores. Despite the success, significant challenges still exist in domain boundary prediction and consistent folding of medium-size proteins (especially beta-proteins) for nonhomologous targets. Further developments of sensitive fold-recognition and ab initio folding methods are critical for solving these problems.

i3Drefine software for protein 3D structure refinement and its assessment in CASP10

Protein structure refinement refers to the process of improving the qualities of protein structures during structure modeling processes to bring them closer to their native states. Structure refinement has been drawing increasing attention in the community-wide Critical Assessment of techniques for Protein Structure prediction (CASP) experiments since its addition in 8^th CASP experiment. During the 9^th and recently concluded 10^th CASP experiments, a consistent growth in number of refinement targets and participating groups has been witnessed. Yet, protein structure refinement still remains a largely unsolved problem with majority of participating groups in CASP refinement category failed to consistently improve the quality of structures issued for refinement. In order to alleviate this need, we developed a completely automated and computationally efficient protein 3D structure refinement method, i3Drefine, based on an iterative and highly convergent energy minimization algorithm with a powerful all-atom composite physics and knowledge-based force fields and hydrogen bonding (HB) network optimization technique. In the recent community-wide blind experiment, CASP10, i3Drefine (as ‘MULTICOM-CONSTRUCT’) was ranked as the best method in the server section as per the official assessment of CASP10 experiment. Here we provide the community with free access to i3Drefine software and systematically analyse the performance of i3Drefine in strict blind mode on the refinement targets issued in CASP10 refinement category and compare with other state-of-the-art refinement methods participating in CASP10. Our analysis demonstrates that i3Drefine is only fully-automated server participating in CASP10 exhibiting consistent improvement over the initial structures in both global and local structural quality metrics. Executable version of i3Drefine is freely available at http://protein.rnet.missouri.edu/i3drefine/.

Mutation of a Short Variable Region in HCpro Protein of Potato virus A Affects Interactions with a Microtubule-Associated Protein and Induces Necrotic Responses in Tobacco

Helper component proteinase (HCpro) is a multifunctional protein of potyviruses (genus Potyvirus). HCpro of Potato virus A (PVA) interacts with the microtubule-associated protein HIP2 in host cells, and depletion of HIP2 reduces virus accumulation. This study shows that HCpro of Potato virus Y and Tobacco etch virus also interact with HIP2. The C-proximal portion of PVA HCpro determines the interaction with HIP2 and was found to contain a stretch of six residues comprising a highly variable region (HVR) in potyviruses. Mutations in HVR reduced PVA accumulation in tobacco plants and induced necrotic symptoms novel to PVA. Microarray and quantitative reverse transcription polymerase chain reaction analyses revealed induction of many defense-related genes including ethylene- and jasmonic acid-inducible pathways in systemically infected leaves at necrosis onset. Salicylic acid-mediated signaling was dispensable for the response. Genes related to microtubule functions were down-regulated. Structural modeling of HCpro suggested that all mutations in HVR caused conformational changes in adjacent regions containing functionally important motifs conserved in potyviruses. Those mutations, which also caused conformational changes in HVR, led to the greatest reduction of fitness. Our results implicate HVR in the regulation of HCpro conformation and virus-host interactions and suggest that mutation of HVR induces host defense.

GalaxyRefine: Protein structure refinement driven by side-chain repacking

The quality of model structures generated by contemporary protein structure prediction methods strongly depends on the degree of similarity between the target and available template structures. Therefore, the importance of improving template-based model structures beyond the accuracy available from template information has been emphasized in the structure prediction community. The GalaxyRefine web server, freely available at http://galaxy.seoklab.org/refine, is based on a refinement method that has been successfully tested in CASP10. The method first rebuilds side chains and performs side-chain repacking and subsequent overall structure relaxation by molecular dynamics simulation. According to the CASP10 assessment, this method showed the best performance in improving the local structure quality. The method can improve both global and local structure quality on average, when used for refining the models generated by state-of-the-art protein structure prediction servers.

Innovations in proteomic profiling of cancers: alternative splice variants as a new class of cancer biomarker candidates and bridging of proteomics with structural biology

Alternative splicing allows a single gene to generate multiple RNA transcripts which can be translated into functionally diverse protein isoforms. Current knowledge of splicing is derived mainly from RNA transcripts, with very little known about the expression level, 3D structures, and functional differences of the proteins. Splicing is a remarkable phenomenon of molecular and biological evolution. Studies which simply report up-regulation or down-regulation of protein or mRNA expression are confounded by the effects of mixtures of these isoforms. Besides understanding the net biological effects of the mixtures, we may be able to develop biomarker tests based on the observable differential expression of particular splice variants or combinations of splice variants in specific disease states. Here we review our work on differential expression of splice variant proteins in cancers and the feasibility of integrating proteomic analysis with structure-based conformational predictions of the differences between such isoforms.

GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

MOTIVATION

Molecular simulation has historically been a low-throughput technique, but faster computers and increasing amounts of genomic and structural data are changing this by enabling large-scale automated simulation of, for instance, many conformers or mutants of biomolecules with or without a range of ligands. At the same time, advances in performance and scaling now make it possible to model complex biomolecular interaction and function in a manner directly testable by experiment. These applications share a need for fast and efficient software that can be deployed on massive scale in clusters, web servers, distributed computing or cloud resources.

RESULTS

Here, we present a range of new simulation algorithms and features developed during the past 4 years, leading up to the GROMACS 4.5 software package. The software now automatically handles wide classes of biomolecules, such as proteins, nucleic acids and lipids, and comes with all commonly used force fields for these molecules built-in. GROMACS supports several implicit solvent models, as well as new free-energy algorithms, and the software now uses multithreading for efficient parallelization even on low-end systems, including windows-based workstations. Together with hand-tuned assembly kernels and state-of-the-art parallelization, this provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations.

AVAILABILITY

GROMACS is an open source and free software available from http://www.gromacs.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential recognition of Old World and New World arenavirus envelope glycoproteins by subtilisin kexin isozyme 1 (SKI-1)/site 1 protease (S1P)

The arenaviruses are an important family of emerging viruses that includes several causative agents of severe hemorrhagic fevers in humans that represent serious public health problems. A crucial step of the arenavirus life cycle is maturation of the envelope glycoprotein precursor (GPC) by the cellular subtilisin kexin isozyme 1 (SKI-1)/site 1 protease (S1P). Comparison of the currently known sequences of arenavirus GPCs revealed the presence of a highly conserved aromatic residue at position P7 relative to the SKI-1/S1P cleavage side in Old World and clade C New World arenaviruses but not in New World viruses of clades A and B or cellular substrates of SKI-1/S1P. Using a combination of molecular modeling and structure-function analysis, we found that residue Y285 of SKI-1/S1P, distal from the catalytic triad, is implicated in the molecular recognition of the aromatic "signature residue" at P7 in the GPC of Old World Lassa virus. Using a quantitative biochemical approach, we show that Y285 of SKI-1/S1P is crucial for the efficient processing of peptides derived from Old World and clade C New World arenavirus GPCs but not of those from clade A and B New World arenavirus GPCs. The data suggest that during coevolution with their mammalian hosts, GPCs of Old World and clade C New World viruses expanded the molecular contacts with SKI-1/S1P beyond the classical four-amino-acid recognition sequences and currently occupy an extended binding pocket.

Molecular architecture of a sodium channel S6 helix: radial tuning of the voltage-gated sodium channel 1.7 activation gate

BACKGROUND

In-frame deletion mutation (Del-L955) in NaV1.7 sodium channel from a kindred with erythromelalgia hyperpolarizes activation.

RESULTS

Del-L955 twists the S6 helix, displacing the Phe960 activation gate. Replacement of Phe960 at the correct helical position depolarizes activation.

CONCLUSION

Radial tuning of the activation gate is critical to the activation of NaV1.7 channel.

SIGNIFICANCE

Structural modeling guided electrophysiology reveals the functional importance of radial tuning of the S6 segment. Voltage-gated sodium (NaV) channels are membrane proteins that consist of 24 transmembrane segments organized into four homologous domains and are essential for action potential generation and propagation. Although the S6 helices of NaV channels line the ion-conducting pore and participate in channel activation, their functional architecture is incompletely understood. Our recent studies show that a naturally occurring in-frame deletion mutation (Del-L955) of NaV1.7 channel, identified in individuals with a severe inherited pain syndrome (inherited erythromelalgia) causes a substantial hyperpolarizing shift of channel activation. Here we took advantage of this deletion mutation to understand the role of the S6 helix in the channel activation. Based on the recently published structure of a bacterial NaV channel (NaVAb), we modeled the WT and Del-L955 channel. Our structural model showed that Del-L955 twists the DII/S6 helix, shifting location and radial orientation of the activation gate residue (Phe(960)). Hypothesizing that these structural changes produce the shift of channel activation of Del-L955 channels, we restored a phenylalanine in wild-type orientation by mutating Ser(961) (Del-L955/S961F), correcting activation by ∼10 mV. Correction of the displaced Phe(960) (F960S) together with introduction of the rescuing activation gate residue (S961F) produced an additional ∼6-mV restoration of activation of the mutant channel. A simple point mutation in the absence of a twist (L955A) did not produce a radial shift and did not hyperpolarize activation. Our results demonstrate the functional importance of radial tuning of the sodium channel S6 helix for the channel activation.

Structural modelling and mutant cycle analysis predict pharmacoresponsiveness of a Na(V)1.7 mutant channel

Sodium channel Na(V)1.7 is critical for human pain signalling. Gain-of-function mutations produce pain syndromes including inherited erythromelalgia, which is usually resistant to pharmacotherapy, but carbamazepine normalizes activation of Na(V)1.7-V400M mutant channels from a family with carbamazepine-responsive inherited erythromelalgia. Here we show that structural modelling and thermodynamic analysis predict pharmacoresponsiveness of another mutant channel (S241T) that is located 159 amino acids distant from V400M. Structural modelling reveals that Na(v)1.7-S241T is ~2.4 Å apart from V400M in the folded channel, and thermodynamic analysis demonstrates energetic coupling of V400M and S241T during activation. Atomic proximity and energetic coupling are paralleled by pharmacological coupling, as carbamazepine (30 μM) depolarizes S214T activation, as previously reported for V400M. Pharmacoresponsiveness of S241T to carbamazepine was further evident at a cellular level, where carbamazepine normalized the hyperexcitability of dorsal root ganglion neurons expressing S241T. We suggest that this approach might identify variants that confer enhanced pharmacoresponsiveness on a variety of channels.

Three dimensional structure of the anthrax toxin translocon-lethal factor complex by cryo-electron microscopy

We have visualized by cryo-electron microscopy (cryo-EM) the complex of the anthrax protective antigen (PA) translocon and the N-terminal domain of anthrax lethal factor (LF(N) inserted into a nanodisc model lipid bilayer. We have determined the structure of this complex at a nominal resolution of 16 Å by single-particle analysis and three-dimensional reconstruction. Consistent with our previous analysis of negatively stained unliganded PA, the translocon comprises a globular structure (cap) separated from the nanodisc bilayer by a narrow stalk that terminates in a transmembrane channel (incompletely distinguished in this reconstruction). The globular cap is larger than the unliganded PA pore, probably due to distortions introduced in the previous negatively stained structures. The cap exhibits larger, more distinct radial protrusions, previously identified with PA domain three, fitted by elements of the NMFF PA prepore crystal structure. The presence of LF(N), though not distinguished due to the seven-fold averaging used in the reconstruction, contributes to the distinct protrusions on the cap rim volume distal to the membrane. Furthermore, the lumen of the cap region is less resolved than the unliganded negatively stained PA, due to the low contrast obtained in our images of this specimen. Presence of the LF(N) extended helix and N terminal unstructured regions may also contribute to this additional internal density within the interior of the cap. Initial NMFF fitting of the cryoEM-defined PA pore cap region positions the Phe clamp region of the PA pore translocon directly above an internal vestibule, consistent with its role in toxin translocation.

Integrated analysis of cytochrome P450 gene superfamily in the red flour beetle, Tribolium castaneum

Background

The functional and evolutionary diversification of insect cytochrome P450s (CYPs) shaped the success of insects. CYPs constitute one of the largest and oldest gene superfamilies that are found in virtually all aerobic organisms. Because of the availability of whole genome sequence and well functioning RNA interference (RNAi), the red flour beetle, Tribolium castaneum serves as an ideal insect model for conducting functional genomics studies. Although several T. castaneum CYPs had been functionally investigated in our previous studies, the roles of the majority of CYPs remain largely unknown. Here, we comprehensively analyzed the phylogenetic relationship of all T. castaneum CYPs with genes in other insect species, investigated the CYP6BQ gene cluster organization, function and evolution, as well as examined the mitochondrial CYPs gene expression patterns and intron-exon organization.

Results

A total 143 CYPs were identified and classified into 26 families and 59 subfamilies. The phylogenetic trees of CYPs among insects across taxa provided evolutionary insight for the genetic distance and function. The percentage of singleton (33.3%) in T. castaneum CYPs is much less than those in Drosophila melanogaster (52.5%) and Bombyx mori (51.2%). Most members in the largest CYP6BQ gene cluster may make contribution to deltamethrin resistance in QTC279 strain. T. castaneum genome encodes nine mitochondrial CYPs, among them CYP12H1 is only expressed in the final instar larval stage. The intron-exon organizations of these mitochondrial CYPs are highly diverse.

Conclusion

Our studies provide a platform to understand the evolution and functions of T. castaneum CYP gene superfamily which will help reveal the strategies employed by insects to cope with their environment.

Genetic determinants of Potato virus Y required to overcome or trigger hypersensitive resistance to PVY strain group O controlled by the gene Ny in potato

Potato virus Y (PVY) (genus Potyvirus) is the most economically damaging and widely distributed virus in potato. Spread of PVY in the field is controlled by growing resistant cultivars. The dominant potato gene Ny(tbr) for hypersensitive resistance (HR) controls ordinary PVY strains (PVY(O)) but is overcome by PVY(N) strains. Studies with infectious PVY chimeras and mutants indicated that the viral determinants necessary and sufficient to overcome Ny(tbr) reside within the helper component proteinase (HC-Pro) (residues 227 to 327). Specifically, eight residues and the modeled three-dimensional conformation of this HC-Pro region distinguish PVY(N) from PVY(O) strains. According to the model, the conserved IGN and CCCT motifs implicated in potyvirus replication and movement, respectively, are situated in a coiled structure and an α-helix, respectively, within this region in PVY(O); however, their locations are reversed in PVY(N). Two residues (R269 and K270) are crucial for the predicted PVY(O)-specific HC-Pro conformation. Two viral chimeras triggered Ny(tbr) and induced veinal necrosis in tobacco, which is novel for PVY. One chimera belonged to strain group PVY(E). Our results suggest a structure-function relationship in recognition of PVY(O) HC-Pro by Ny(tbr), reveal HC-Pro amino acid signatures specific to PVY(O) and PVY(N), and facilitate identification of PVY strains overcoming Ny(tbr).

Protein Structure Refinement through Structure Selection and Averaging from Molecular Dynamics Ensembles

A molecular dynamics (MD) simulation based protocol for structure refinement of template-based model predictions is described. The protocol involves the application of restraints, ensemble averaging of selected subsets, interpolation between initial and refined structures, and assessment of refinement success. It is found that sub-microsecond MD-based sampling when combined with ensemble averaging can produce moderate but consistent refinement for most systems in the CASP targets considered here.

Mining tertiary structural motifs for assessment of designability

The observation of a limited secondary-structural alphabet in native proteins, with significant sequence preferences, has profoundly influenced the fields of protein design and structure prediction (Simons, Kooperberg, Huang, & Baker, 1997; Verschueren et al., 2011). In the era of structural genomics, as the size of the structural dataset continues to grow rapidly, it is becoming possible to extend this analysis to tertiary structural motifs and their sequences. For a hypothetical tertiary motif, the rate of its utilization in natural proteins may be used to assess its designability-the ease with which the motif can be realized with natural amino acids. This requires a structural similarity search methodology, which rather than looking for global topological agreement (more appropriate for categorization of full proteins or domains), identifies detailed geometric matches. In this chapter, we introduce such a method, called MaDCaT, and demonstrate its use by assessing the designability landscapes of two tertiary structural motifs. We also show that such analysis can establish structure/sequence links by providing the sequence constraints necessary to encode designable motifs. As logical extension of their secondary-structure counterparts, tertiary structural preferences will likely prove extremely useful in de novo protein design and structure prediction.

3Drefine: consistent protein structure refinement by optimizing hydrogen bonding network and atomic-level energy minimization

One of the major limitations of computational protein structure prediction is the deviation of predicted models from their experimentally derived true, native structures. The limitations often hinder the possibility of applying computational protein structure prediction methods in biochemical assignment and drug design that are very sensitive to structural details. Refinement of these low-resolution predicted models to high-resolution structures close to the native state, however, has proven to be extremely challenging. Thus, protein structure refinement remains a largely unsolved problem. Critical assessment of techniques for protein structure prediction (CASP) specifically indicated that most predictors participating in the refinement category still did not consistently improve model quality. Here, we propose a two-step refinement protocol, called 3Drefine, to consistently bring the initial model closer to the native structure. The first step is based on optimization of hydrogen bonding (HB) network and the second step applies atomic-level energy minimization on the optimized model using a composite physics and knowledge-based force fields. The approach has been evaluated on the CASP benchmark data and it exhibits consistent improvement over the initial structure in both global and local structural quality measures. 3Drefine method is also computationally inexpensive, consuming only few minutes of CPU time to refine a protein of typical length (300 residues). 3Drefine web server is freely available at http://sysbio.rnet.missouri.edu/3Drefine/.

Homology modeling of dopamine D2 and D3 receptors: molecular dynamics refinement and docking evaluation

Dopamine (DA) receptors, a class of G-protein coupled receptors (GPCRs), have been targeted for drug development for the treatment of neurological, psychiatric and ocular disorders. The lack of structural information about GPCRs and their ligand complexes has prompted the development of homology models of these proteins aimed at structure-based drug design. Crystal structure of human dopamine D₃ (hD₃) receptor has been recently solved. Based on the hD₃ receptor crystal structure we generated dopamine D₂ and D₃ receptor models and refined them with molecular dynamics (MD) protocol. Refined structures, obtained from the MD simulations in membrane environment, were subsequently used in molecular docking studies in order to investigate potential sites of interaction. The structure of hD₃ and hD_2L receptors was differentiated by means of MD simulations and D₃ selective ligands were discriminated, in terms of binding energy, by docking calculation. Robust correlation of computed and experimental K_i was obtained for hD₃ and hD_2L receptor ligands. In conclusion, the present computational approach seems suitable to build and refine structure models of homologous dopamine receptors that may be of value for structure-based drug discovery of selective dopaminergic ligands.

A novel Drosophila SOD2 mutant demonstrates a role for mitochondrial ROS in neurodevelopment and disease

Reactive oxygen species (ROS) play essential roles in cell signaling, survival, and homeostasis. Aberrant ROS lead to disease and contribute to the aging process. Numerous enzymes and vigilant antioxidant pathways are required to regulate ROS for normal cellular health. Mitochondria are a major source of ROS, and mechanisms to prevent elevated ROS during oxidative phosphorylation require super oxide dismutase (SOD) activity. SOD2, also known as MnSOD, is targeted to mitochondria and is instrumental in regulating ROS by conversion of superoxides to hydrogen peroxide, which is further broken down into H(2)O and oxygen. Here, we describe the identification of a novel mutation within the mitochondrial SOD2 enzyme in Drosophila that results in adults with an extremely shortened life span, sensitivity to hyperoxia, and neuropathology. Additional studies demonstrate that this novel mutant, SOD2(bewildered), exhibits abnormal brain morphology, suggesting a critical role for this protein in neurodevelopment. We investigated the basis of this neurodevelopmental defect and discovered an increase in aberrant axonal that could underlie the aberrant neurodevelopment and brain morphology defects. This novel allele, SOD2(bewildered), provides a unique opportunity to study the effects of increased mitochondrial ROS on neural development, axonal targeting, and neural cell degeneration in vivo.

What is the best reference state for designing statistical atomic potentials in protein structure prediction?

Many statistical potentials were developed in last two decades for protein folding and protein structure recognition. The major difference of these potentials is on the selection of reference states to offset sampling bias. However, since these potentials used different databases and parameter cutoffs, it is difficult to judge what the best reference states are by examining the original programs. In this study, we aim to address this issue and evaluate the reference states by a unified database and programming environment. We constructed distance-specific atomic potentials using six widely-used reference states based on 1022 high-resolution protein structures, which are applied to rank modeling in six sets of structure decoys. The reference state on random-walk chain outperforms others in three decoy sets while those using ideal-gas, quasi-chemical approximation and averaging sample stand out in one set separately. Nevertheless, the performance of the potentials relies on the origin of decoy generations and no reference state can clearly outperform others in all decoy sets. Further analysis reveals that the statistical potentials have a contradiction between the universality and pertinence, and optimal reference states should be extracted based on specific application environments and decoy spaces.

Progress at last

Refinement of protein structures from a correct topology to atomically detailed resolution has proven remarkably difficult. Jian et al. (in this issue of Structure) illustrate a significant advance in this task by carefully incorporating into the refinement process many body interactions extracted from fragment statistics.

A protocol for computer-based protein structure and function prediction

Genome sequencing projects have ciphered millions of protein sequence, which require knowledge of their structure and function to improve the understanding of their biological role. Although experimental methods can provide detailed information for a small fraction of these proteins, computational modeling is needed for the majority of protein molecules which are experimentally uncharacterized. The I-TASSER server is an on-line workbench for high-resolution modeling of protein structure and function. Given a protein sequence, a typical output from the I-TASSER server includes secondary structure prediction, predicted solvent accessibility of each residue, homologous template proteins detected by threading and structure alignments, up to five full-length tertiary structural models, and structure-based functional annotations for enzyme classification, Gene Ontology terms and protein-ligand binding sites. All the predictions are tagged with a confidence score which tells how accurate the predictions are without knowing the experimental data. To facilitate the special requests of end users, the server provides channels to accept user-specified inter-residue distance and contact maps to interactively change the I-TASSER modeling; it also allows users to specify any proteins as template, or to exclude any template proteins during the structure assembly simulations. The structural information could be collected by the users based on experimental evidences or biological insights with the purpose of improving the quality of I-TASSER predictions. The server was evaluated as the best programs for protein structure and function predictions in the recent community-wide CASP experiments. There are currently >20,000 registered scientists from over 100 countries who are using the on-line I-TASSER server.