Protein dynamics with off-lattice Monte Carlo moves

A Monte Carlo method for dynamics simulation of all-atom protein models is introduced, to reach long times not accessible to conventional molecular dynamics. The considered degrees of freedom are the dihedrals at C$_\alpha$-atoms. Two Monte Carlo moves are used: single rotations about torsion axes, and cooperative rotations in windows of amide planes, changing the conformation globally and locally, respectively. For local moves Jacobians are used to obtain an unbiased distribution of dihedrals. A molecular dynamics energy function adapted to the protein model is employed. A polypeptide is folded into native-like structures by local but not by global moves.Comment: 10 pages, 4 Postscript figures, uses epsf.sty and a4.sty; scheduled tentatively for Phys.Rev.E issue of 1 March 199

Prediction of conformational B-cell epitopes from 3D structures by random forests with a distance-based feature

<p>Abstract</p> <p>Background</p> <p>Antigen-antibody interactions are key events in immune system, which provide important clues to the immune processes and responses. In Antigen-antibody interactions, the specific sites on the antigens that are directly bound by the B-cell produced antibodies are well known as B-cell epitopes. The identification of epitopes is a hot topic in bioinformatics because of their potential use in the epitope-based drug design. Although most B-cell epitopes are discontinuous (or conformational), insufficient effort has been put into the conformational epitope prediction, and the performance of existing methods is far from satisfaction.</p> <p>Results</p> <p>In order to develop the high-accuracy model, we focus on some possible aspects concerning the prediction performance, including the impact of interior residues, different contributions of adjacent residues, and the imbalanced data which contain much more non-epitope residues than epitope residues. In order to address above issues, we take following strategies. Firstly, a concept of 'thick surface patch' instead of 'surface patch' is introduced to describe the local spatial context of each surface residue, which considers the impact of interior residue. The comparison between the thick surface patch and the surface patch shows that interior residues contribute to the recognition of epitopes. Secondly, statistical significance of the distance distribution difference between non-epitope patches and epitope patches is observed, thus an adjacent residue distance feature is presented, which reflects the unequal contributions of adjacent residues to the location of binding sites. Thirdly, a bootstrapping and voting procedure is adopted to deal with the imbalanced dataset. Based on the above ideas, we propose a new method to identify the B-cell conformational epitopes from 3D structures by combining conventional features and the proposed feature, and the random forest (RF) algorithm is used as the classification engine. The experiments show that our method can predict conformational B-cell epitopes with high accuracy. Evaluated by leave-one-out cross validation (LOOCV), our method achieves the mean AUC value of 0.633 for the benchmark bound dataset, and the mean AUC value of 0.654 for the benchmark unbound dataset. When compared with the state-of-the-art prediction models in the independent test, our method demonstrates comparable or better performance.</p> <p>Conclusions</p> <p>Our method is demonstrated to be effective for the prediction of conformational epitopes. Based on the study, we develop a tool to predict the conformational epitopes from 3D structures, available at <url>http://code.google.com/p/my-project-bpredictor/downloads/list</url>.</p

Comparison of the structure and activity of glycosylated and asglycosylated human carboxylesterase 1

Human Carboxylesterase 1 (hCES1) is the key liver microsomal enzyme responsible for detoxification and metabolism of a variety of clinical drugs. To analyse the role of the single N-linked glycan on the structure and activity of the enzyme, authentically glycosylated and aglycosylated hCES1, generated by mutating asparagine 79 to glutamine, were produced in human embryonic kidney cells. Purified enzymes were shown to be predominantly trimeric in solution by analytical ultracentrifugation. The purified aglycosylated enzyme was found to be more active than glycosylated hCES1 and analysis of enzyme kinetics revealed that both enzymes exhibit positive cooperativity. Crystal structures of hCES1 a catalytically inactive mutant (S221A) and the aglycosylated enzyme were determined in the absence of any ligand or substrate to high resolutions (1.86 Å, 1.48 Å and 2.01 Å, respectively). Superposition of all three structures showed only minor conformational differences with a root mean square deviations of around 0.5 Å over all Cα positions. Comparison of the active sites of these un-liganded enzymes with the structures of hCES1-ligand complexes showed that side-chains of the catalytic triad were pre-disposed for substrate binding. Overall the results indicate that preventing N-glycosylation of hCES1 does not significantly affect the structure or activity of the enzyme

Disordered Structural Ensembles of Vasopressin and Oxytocin and Their Mutants

Vasopressin and oxytocin are intrinsically disordered cyclic nonapeptides belonging to a family of neurohypophysial hormones. Although unique in their functions, these peptides differ only by two residues and both feature a tocin ring formed by the disulfide bridge between first and sixth cysteine residues. This sequence and structural similarity are experimentally linked to oxytocin agonism at vasopressin receptors and vasopressin antagonism at oxytocin receptors. Yet single- or double-residue mutations in both peptides have been shown to have drastic impacts on their activities at either receptor, and possibly the ability to bind to their neurophysin carrier protein. In this study we perform molecular dynamics simulations of the unbound native and mutant sequences of the oxytocin and vasopressin hormones to characterize their structural ensembles. We classify the subpopulations of these structural ensembles on the basis of the distributions of radius of gyration and secondary structure and hydrogen-bonding features of the canonical tocin ring and disordered tail region. We then relate the structural changes observed in the unbound form of the different hormone sequences to experimental information about peptide receptor binding, and more indirectly, carrier protein binding affinity, receptor activity, and protease degradation. This study supports the hypothesis that the structural characteristics of the unbound form of an IDP can be used to predict structural or functional preferences of its functional bound form

Structurally encoded intraclass differences in EphA clusters drive distinct cell responses

Functional outcomes of ephrin binding to Eph receptors (Ephs) range from cell repulsion to adhesion. Here we used cell collapse and stripe assays, showing contrasting effects of human ephrinA5 binding to EphA2 and EphA4. Despite equivalent ligand binding affinities, EphA4 triggered greater cell collapse, whereas EphA2-expressing cells adhered better to ephrinA5-coated surfaces. Chimeric receptors showed that the ectodomain is a major determinant of cell response. We report crystal structures of EphA4 ectodomain alone and in complexes with ephrinB3 and ephrinA5. These revealed closed clusters with a dimeric or circular arrangement in the crystal lattice, contrasting with extended arrays previously observed for EphA2 ectodomain. Localization microscopy showed that ligand-stimulated EphA4 induces smaller clusters than does EphA2. Mutant Ephs link these characteristics to interactions observed in the crystal lattices, suggesting a mechanism by which distinctive ectodomain surfaces determine clustering, and thereby signaling, properties. © 2013 Nature America, Inc. All rights reserved

Analog methods: computation and simulation

First ed., 1954, by Walter W. Soroka, published under title: Analog methods in computation and simulation.Mode of access: Internet