The Bio3D packages for structural bioinformatics.

Bio3D is a family of R packages for the analysis of biomolecular sequence, structure, and dynamics. Major functionality includes biomolecular database searching and retrieval, sequence and structure conservation analysis, ensemble normal mode analysis, protein structure and correlation network analysis, principal component, and related multivariate analysis methods. Here, we review recent package developments, including a new underlying segregation into separate packages for distinct analysis, and introduce a new method for structure analysis named ensemble difference distance matrix analysis (eDDM). The eDDM approach calculates and compares atomic distance matrices across large sets of homologous atomic structures to help identify the residue wise determinants underlying specific functional processes. An eDDM workflow is detailed along with an example application to a large protein family. As a new member of the Bio3D family, the Bio3D-eddm package supports both experimental and theoretical simulation-generated structures, is integrated with other methods for dissecting sequence-structure-function relationships, and can be used in a highly automated and reproducible manner. Bio3D is distributed as an integrated set of platform independent open source R packages available from: http://thegrantlab.org/bio3d/.

Comparative Protein Structure Analysis with Bio3D-Web.

Bio3D-web is an online application for the interactive analysis of sequence-structure-dynamics relationships in user-defined protein structure sets. Major functionality includes structure database searching, sequence and structure conservation assessment, inter-conformer relationship mapping and clustering with principal component analysis (PCA), and flexibility prediction and comparison with ensemble normal mode analysis (eNMA). Collectively these methods allow users to start with a single sequence or structure and characterize the structural, conformational, and internal dynamic properties of homologous proteins for which there are high-resolution structures available. Functionality is also provided for the generation of custom PDF, Word, and HTML analysis reports detailing all user-specified analysis settings and corresponding results. Bio3D-web is available at http://thegrantlab.org/bio3d/webapps , as a Docker image https://hub.docker.com/r/bio3d/bio3d-web/ , or downloadable source code https://bitbucket.org/Grantlab/bio3d-web .

A Pharmacological Chaperone Therapy for Acute Intermittent Porphyria.

Mutations in hydroxymethylbilane synthase (HMBS) cause acute intermittent porphyria (AIP), an autosomal dominant disease where typically only one HMBS allele is mutated. In AIP, the accumulation of porphyrin precursors triggers life-threatening neurovisceral attacks and at long-term, entails an increased risk of hepatocellular carcinoma, kidney failure, and hypertension. Today, the only cure is liver transplantation, and a need for effective mechanism-based therapies, such as pharmacological chaperones, is prevailing. These are small molecules that specifically stabilize a target protein. They may be developed into an oral treatment, which could work curatively during acute attacks, but also prophylactically in asymptomatic HMBS mutant carriers. With the use of a 10,000 compound library, we identified four binders that further increased the initially very high thermal stability of wild-type HMBS and protected the enzyme from trypsin digestion. The best hit and a selected analog increased steady-state levels and total HMBS activity in human hepatoma cells overexpressing HMBS, and in an Hmbs-deficient mouse model with a low-expressed wild-type-like allele, compared to untreated controls. Moreover, the concentration of porphyrin precursors decreased in liver of mice treated with the best hit. Our findings demonstrate the great potential of these hits for the development of a pharmacological chaperone-based corrective treatment of AIP by enhancing wild-type HMBS function independently of the patients' specific mutation.

Structure of full-length human phenylalanine hydroxylase in complex with tetrahydrobiopterin.

Phenylalanine hydroxylase (PAH) is a key enzyme in the catabolism of phenylalanine, and mutations in this enzyme cause phenylketonuria (PKU), a genetic disorder that leads to brain damage and mental retardation if untreated. Some patients benefit from supplementation with a synthetic formulation of the cofactor tetrahydrobiopterin (BH4) that partly acts as a pharmacological chaperone. Here we present structures of full-length human PAH (hPAH) both unbound and complexed with BH4 in the precatalytic state. Crystal structures, solved at 3.18-Å resolution, show the interactions between the cofactor and PAH, explaining the negative regulation exerted by BH4 BH4 forms several H-bonds with the N-terminal autoregulatory tail but is far from the catalytic FeII Upon BH4 binding a polar and salt-bridge interaction network links the three PAH domains, explaining the stability conferred by BH4 Importantly, BH4 binding modulates the interaction between subunits, providing information about PAH allostery. Moreover, we also show that the cryo-EM structure of hPAH in absence of BH4 reveals a highly dynamic conformation for the tetramers. Structural analyses of the hPAH:BH4 subunits revealed that the substrate-induced movement of Tyr138 into the active site could be coupled to the displacement of BH4 from the precatalytic toward the active conformation, a molecular mechanism that was supported by site-directed mutagenesis and targeted molecular dynamics simulations. Finally, comparison of the rat and human PAH structures show that hPAH is more dynamic, which is related to amino acid substitutions that enhance the flexibility of hPAH and may increase the susceptibility to PKU-associated mutations.

Validation of electron microscopy initial models via small angle X-ray scattering curves.

MOTIVATION: Cryo electron microscopy (EM) is currently one of the main tools to reveal the structural information of biological macromolecules. The re-construction of three-dimensional (3D) maps is typically carried out following an iterative process that requires an initial estimation of the 3D map to be refined in subsequent steps. Therefore, its determination is key in the quality of the final results, and there are cases in which it is still an open issue in single particle analysis (SPA). Small angle X-ray scattering (SAXS) is a well-known technique applied to structural biology. It is useful from small nanostructures up to macromolecular ensembles for its ability to obtain low resolution information of the biological sample measuring its X-ray scattering curve. These curves, together with further analysis, are able to yield information on the sizes, shapes and structures of the analyzed particles. RESULTS: In this paper, we show how the low resolution structural information revealed by SAXS is very useful for the validation of EM initial 3D models in SPA, helping the following refinement process to obtain more accurate 3D structures. For this purpose, we approximate the initial map by pseudo-atoms and predict the SAXS curve expected for this pseudo-atomic structure. The match between the predicted and experimental SAXS curves is considered as a good sign of the correctness of the EM initial map. AVAILABILITY AND IMPLEMENTATION: The algorithm is freely available as part of the Scipion 1.2 software at http://scipion.i2pc.es/.

Conformational stabilization as a strategy to prevent nucleophosmin mislocalization in leukemia.

Nucleophosmin (NPM) is a nucleolar protein involved in ribosome assembly and cell homeostasis. Mutations in the C-terminal domain of NPM that impair native folding and localization are associated with acute myeloid leukemia (AML). We have performed a high-throughput screening searching for compounds that stabilize the C-terminal domain. We identified three hit compounds which show the ability to increase the thermal stability of both the C-terminal domain as well as full-length NPM. The best hit also seemed to favor folding of an AML-like mutant. Computational pocket identification and molecular docking support a stabilization mechanism based on binding of the phenyl/benzene group of the compounds to a particular hydrophobic pocket and additional polar interactions with solvent-accessible residues. Since these results indicate a chaperoning potential of our candidate hits, we tested their effect on the subcellular localization of AML-like mutants. Two compounds partially alleviated the aggregation and restored nucleolar localization of misfolded mutants. The identified hits appear promising as pharmacological chaperones aimed at therapies for AML based on conformational stabilization of NPM.

Investigating Protein Sequence-structure-dynamics Relationships with Bio3D-web.

We demonstrate the usage of Bio3D-web for the interactive analysis of biomolecular structure data. The Bio3D-web application provides online functionality for: (1) The identification of related protein structure sets to user specified thresholds of similarity; (2) Their multiple alignment and structure superposition; (3) Sequence and structure conservation analysis; (4) Inter-conformer relationship mapping with principal component analysis, and (5) comparison of predicted internal dynamics via ensemble normal mode analysis. This integrated functionality provides a complete online workflow for investigating sequence-structure-dynamic relationships within protein families and superfamilies.

Stable preparations of tyrosine hydroxylase provide the solution structure of the full-length enzyme.

Tyrosine hydroxylase (TH) catalyzes the rate-limiting step in the biosynthesis of catecholamine neurotransmitters. TH is a highly complex enzyme at mechanistic, structural, and regulatory levels, and the preparation of kinetically and conformationally stable enzyme for structural characterization has been challenging. Here, we report on improved protocols for purification of recombinant human TH isoform 1 (TH1), which provide large amounts of pure, stable, active TH1 with an intact N-terminus. TH1 purified through fusion with a His-tagged maltose-binding protein on amylose resin was representative of the iron-bound functional enzyme, showing high activity and stabilization by the natural feedback inhibitor dopamine. TH1 purified through fusion with a His-tagged ZZ domain on TALON is remarkably stable, as it was partially inhibited by resin-derived cobalt. This more stable enzyme preparation provided high-quality small-angle X-ray scattering (SAXS) data and reliable structural models of full-length tetrameric TH1. The SAXS-derived model reveals an elongated conformation (Dmax = 20 nm) for TH1, different arrangement of the catalytic domains compared with the crystal structure of truncated forms, and an N-terminal region with an unstructured tail that hosts the phosphorylation sites and a separated Ala-rich helical motif that may have a role in regulation of TH by interacting with binding partners.

Online interactive analysis of protein structure ensembles with Bio3D-web.

Bio3D-web is an online application for analyzing the sequence, structure and conformational heterogeneity of protein families. Major functionality is provided for identifying protein structure sets for analysis, their alignment and refined structure superposition, sequence and structure conservation analysis, mapping and clustering of conformations and the quantitative comparison of their predicted structural dynamics. AVAILABILITY: Bio3D-web is based on the Bio3D and Shiny R packages. All major browsers are supported and full source code is available under a GPL2 license from http://thegrantlab.org/bio3d-web CONTACT: bjgrant@umich.edu or lars.skjarven@uib.no.

Rapid Characterization of Allosteric Networks with Ensemble Normal Mode Analysis.

Allosteric regulation is a primary means of controlling protein function. By definition, allostery involves the propagation of structural dynamic changes between distal protein sites that yields a functional change. Gaining improved knowledge of these fundamental mechanisms is important for understanding many biomolecular processes and for guiding protein engineering and drug design efforts. In this work we compare and contrast a range of normal mode analysis (NMA) approaches together with network analysis for the prediction of structural dynamics and allosteric sites. Application to heterotrimeric G proteins, hemoglobin, and caspase 7 indicates that atomistic elastic network models provide improved predictions of experimental allosteric mutation sites. Results for G proteins also display an improved consistency with those derived from more computationally demanding MD simulations. Application of this approach across available experimental structures for a given protein family in a unified manner, that we refer to as ensemble NMA, yields the best overall predictive performance. We propose that this atomistic ensemble NMA approach represents an efficient and powerful tool for guiding the exploration of coupled motions and allosteric mechanisms in cases where multiple structures are available and where MD may prove prohibitively expensive.

Dynamic Coupling and Allosteric Networks in the α Subunit of Heterotrimeric G Proteins.

G protein α subunits cycle between active and inactive conformations to regulate a multitude of intracellular signaling cascades. Important structural transitions occurring during this cycle have been characterized from extensive crystallographic studies. However, the link between observed conformations and the allosteric regulation of binding events at distal sites critical for signaling through G proteins remain unclear. Here we describe molecular dynamics simulations, bioinformatics analysis, and experimental mutagenesis that identifies residues involved in mediating the allosteric coupling of receptor, nucleotide, and helical domain interfaces of Gαi. Most notably, we predict and characterize novel allosteric decoupling mutants, which display enhanced helical domain opening, increased rates of nucleotide exchange, and constitutive activity in the absence of receptor activation. Collectively, our results provide a framework for explaining how binding events and mutations can alter internal dynamic couplings critical for G protein function.

Discovery of a Specific Inhibitor of Pyomelanin Synthesis in Legionella pneumophila.

Phenylalanine hydroxylase catalyzes the first step in the synthesis of pyomelanin, a pigment that aids in the acquisition of essential iron in certain bacteria. In this work, we present the development and application of a drug discovery protocol by targeting this enzyme in Legionella pneumophila, the major causative agent of Legionnaires' disease. We employ a combination of high-throughput screening to identify small-molecule binders, enzymatic activity measurements to identify inhibitors in vitro, and the verification of the inhibitory effect in vivo. The most potent inhibitor shows an IC50 value in the low micromolar range and successfully abolishes the synthesis of pyomelanin in L. pneumophila cultures at 10 µM. Thus, this compound represents a novel and effective tool for investigating the role of pyomelanin in the biology and pathogenicity of this organism. Altogether, the results demonstrate a successful pathway for drug development focusing on binding specificity in the initial high-throughput screening steps.

Dynamics, flexibility, and allostery in molecular chaperonins.

The chaperonins are a family of molecular chaperones present in all three kingdoms of life. They are classified into Group I and Group II. Group I consists of the bacterial variants (GroEL) and the eukaryotic ones from mitochondria and chloroplasts (Hsp60), while Group II consists of the archaeal (thermosomes) and eukaryotic cytosolic variants (CCT or TRiC). Both groups assemble into a dual ring structure, with each ring providing a protective folding chamber for nascent and denatured proteins. Their functional cycle is powered by ATP binding and hydrolysis, which drives a series of structural rearrangements that enable encapsulation and subsequent release of the substrate protein. Chaperonins have elaborate allosteric mechanisms to regulate their functional cycle. Long-range negative cooperativity between the two rings ensures alternation of the folding chambers. Positive intra-ring cooperativity, which facilitates concerted conformational transitions within the protein subunits of one ring, has only been demonstrated for Group I chaperonins. In this review, we describe our present understanding of the underlying mechanisms and the structure-function relationships in these complex protein systems with a particular focus on the structural dynamics, allostery, and associated conformational rearrangements.

WEBnm@ v2.0: Web server and services for comparing protein flexibility.

BACKGROUND: Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics-function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins. RESULTS: We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma . CONCLUSION: WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.

Integrating protein structural dynamics and evolutionary analysis with Bio3D.

BACKGROUND: Popular bioinformatics approaches for studying protein functional dynamics include comparisons of crystallographic structures, molecular dynamics simulations and normal mode analysis. However, determining how observed displacements and predicted motions from these traditionally separate analyses relate to each other, as well as to the evolution of sequence, structure and function within large protein families, remains a considerable challenge. This is in part due to the general lack of tools that integrate information of molecular structure, dynamics and evolution. RESULTS: Here, we describe the integration of new methodologies for evolutionary sequence, structure and simulation analysis into the Bio3D package. This major update includes unique high-throughput normal mode analysis for examining and contrasting the dynamics of related proteins with non-identical sequences and structures, as well as new methods for quantifying dynamical couplings and their residue-wise dissection from correlation network analysis. These new methodologies are integrated with major biomolecular databases as well as established methods for evolutionary sequence and comparative structural analysis. New functionality for directly comparing results derived from normal modes, molecular dynamics and principal component analysis of heterogeneous experimental structure distributions is also included. We demonstrate these integrated capabilities with example applications to dihydrofolate reductase and heterotrimeric G-protein families along with a discussion of the mechanistic insight provided in each case. CONCLUSIONS: The integration of structural dynamics and evolutionary analysis in Bio3D enables researchers to go beyond a prediction of single protein dynamics to investigate dynamical features across large protein families. The Bio3D package is distributed with full source code and extensive documentation as a platform independent R package under a GPL2 license from http://thegrantlab.org/bio3d/ .

The structure of the box C/D enzyme reveals regulation of RNA methylation.

Post-transcriptional modifications are essential to the cell life cycle, as they affect both pre-ribosomal RNA processing and ribosome assembly. The box C/D ribonucleoprotein enzyme that methylates ribosomal RNA at the 2'-O-ribose uses a multitude of guide RNAs as templates for the recognition of rRNA target sites. Two methylation guide sequences are combined on each guide RNA, the significance of which has remained unclear. Here we use a powerful combination of NMR spectroscopy and small-angle neutron scattering to solve the structure of the 390 kDa archaeal RNP enzyme bound to substrate RNA. We show that the two methylation guide sequences are located in different environments in the complex and that the methylation of physiological substrates targeted by the same guide RNA occurs sequentially. This structure provides a means for differential control of methylation levels at the two sites and at the same time offers an unexpected regulatory mechanism for rRNA folding.

Accounting for conformational variability in protein-ligand docking with NMR-guided rescoring.

A key component to success in structure-based drug design is reliable information on protein-ligand interactions. Recent development in NMR techniques has accelerated this process by overcoming some of the limitations of X-ray crystallography and computational protein-ligand docking. In this work we present a new scoring protocol based on NMR-derived interligand INPHARMA NOEs to guide the selection of computationally generated docking modes. We demonstrate the performance in a range of scenarios, encompassing traditionally difficult cases such as docking to homology models and ligand dependent domain rearrangements. Ambiguities associated with sparse experimental information are lifted by searching a consensus solution based on simultaneously fitting multiple ligand pairs. This study provides a previously unexplored integration between molecular modeling and experimental data, in which interligand NOEs represent the key element in the rescoring algorithm. The presented protocol should be widely applicable for protein-ligand docking also in a different context from drug design and highlights the important role of NMR-based approaches to describe intermolecular ligand-receptor interactions.

The peripheral binding of 14-3-3γ to membranes involves isoform-specific histidine residues.

Mammalian 14-3-3 protein scaffolds include seven conserved isoforms that bind numerous phosphorylated protein partners and regulate many cellular processes. Some 14-3-3-isoforms, notably γ, have elevated affinity for membranes, which might contribute to modulate the subcellular localization of the partners and substantiate the importance of investigating molecular mechanisms of membrane interaction. By applying surface plasmon resonance we here show that the binding to phospholipid bilayers is stimulated when 14-3-3γ is complexed with its partner, a peptide corresponding to the Ser19-phosphorylated N-terminal region of tyrosine hydroxylase. Moreover, membrane interaction is dependent on salts of kosmotropic ions, which also stabilize 14-3-3γ. Electrostatic analysis of available crystal structures of γ and of the non-membrane-binding ζ-isoform, complemented with molecular dynamics simulations, indicate that the electrostatic potential distribution of phosphopeptide-bound 14-3-3γ is optimal for interaction with the membrane through amphipathic helices at the N-terminal dimerization region. In addition, His158, and especially His195, both specific to 14-3-3γ and located at the convex lateral side, appeared to be pivotal for the ligand induced membrane interaction, as corroborated by site-directed mutagenesis. The participation of these histidine residues might be associated to their increased protonation upon membrane binding. Overall, these results reveal membrane-targeting motifs and give insights on mechanisms that furnish the 14-3-3γ scaffold with the capacity for tuned shuffling from soluble to membrane-bound states.

A dynamic model of long-range conformational adaptations triggered by nucleotide binding in GroEL-GroES.

The molecular chaperone, GroEL, essential for correct protein folding in E. coli, is composed of 14 identical subunits organized in two interacting rings, each providing a folding chamber for non-native substrate proteins. The oligomeric assembly shows positive cooperativity within each ring and negative cooperativity between the rings. Although it is well known that ATP and long-range allosteric interactions drive the functional cycle of GroEL, an atomic resolution view of how ligand binding modulates conformational adaptations over long distances remains a major challenge. Moreover, little is known on the relation between equilibrium dynamics at physiological temperatures and the allosteric transitions in GroEL. Here we present multiple all-atom molecular dynamics simulations of the GroEL-GroES assemblies at different stages of the functional cycle. Combined with an extensive analysis of the complete set of experimentally available structures, principal component analysis and conformer plots, we provide an explicit evaluation of the accessible conformational space of unliganded GroEL. Our results suggest the presence of pre-existing conformers at the equatorial domain level, and a shift of the conformational ensemble upon ATP-binding. At the inter-ring interface the simulations capture a remarkable offset motion of helix D triggered by ATP-binding to the folding active ring. The reorientation of helix D, previously only observed upon GroES association, correlates with a change of the internal dynamics in the opposite ring. This work contributes to the understanding of the molecular mechanisms in GroEL and highlights the ability of all-atom MD simulations to model long-range structural changes and allosteric events in large systems.

Dynamics, flexibility and ligand-induced conformational changes in biological macromolecules: a computational approach.

Biomolecules possess important dynamical properties that enable them to adapt and alternate their conformation as a response to environmental stimuli. Recent advancements in computational resources and methodology allow a higher capability to mimic in vitro conditions and open up the possibility of studying large systems over longer timescales. Here, we describe commonly used computational approaches for studying the dynamic properties of proteins. We review a selected set of simulation studies on ligand-induced changes in the chaperonin GroEL-GroES, a molecular folding machine, maltose-binding protein, a prototypical member of the periplasmic binding proteins, and the bacterial ribosomal A-site, focusing on aminoglycoside antibiotic recognition. We also discuss a recent quantitative reconstruction of the binding process of benzamidine and trypsin. These studies contribute to the understanding and further development of the medicinal regulation of large biomolecular systems.