Alexander S. Rose

The RCSB protein data bank: integrative view of protein, gene and 3D structural information

The Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB, http://rcsb.org), the US data center for the global PDB archive, makes PDB data freely available to all users, from structural biologists to computational biologists and beyond. New tools and resources have been added to the RCSB PDB web portal in support of a ‘Structural View of Biology.’ Recent developments have improved the User experience, including the high-speed NGL Viewer that provides 3D molecular visualization in any web browser, improved support for data file download and enhanced organization of website pages for query, reporting and individual structure exploration. Structure validation information is now visible for all archival entries. PDB data have been integrated with external biological resources, including chromosomal position within the human genome; protein modifications; and metabolic pathways. PDB-101 educational materials have been reorganized into a searchable website and expanded to include new features such as the Geis Digital Archive.

NGL Viewer: a web application for molecular visualization

The NGL Viewer (http://proteinformatics.charite.de/ngl) is a web application for the visualization of macromolecular structures. By fully adopting capabilities of modern web browsers, such as WebGL, for molecular graphics, the viewer can interactively display large molecular complexes and is also unaffected by the retirement of third-party plug-ins like Flash and Java Applets. Generally, the web application offers comprehensive molecular visualization through a graphical user interface so that life scientists can easily access and profit from available structural data. It supports common structural file-formats (e.g. PDB, mmCIF) and a variety of molecular representations (e.g. ‘cartoon, spacefill, licorice’). Moreover, the viewer can be embedded in other web sites to provide specialized visualizations of entries in structural databases or results of structure-related calculations.

Crystal structure of a common GPCR-binding interface for G protein and arrestin.

G-protein-coupled receptors (GPCRs) transmit extracellular signals to activate intracellular heterotrimeric G proteins (Gαβγ) and arrestins. For G protein signalling, the Gα C-terminus (GαCT) binds to a cytoplasmic crevice of the receptor that opens upon activation. A consensus motif is shared among GαCT from the Gi/Gt family and the ‘finger loop’ region (ArrFL1–4) of all four arrestins. Here we present a 2.75 Å crystal structure of ArrFL-1, a peptide analogue of the finger loop of rod photoreceptor arrestin, in complex with the prototypical GPCR rhodopsin. Functional binding of ArrFL to the receptor was confirmed by ultraviolet-visible absorption spectroscopy, competitive binding assays and Fourier transform infrared spectroscopy. For both GαCT and ArrFL, binding to the receptor crevice induces a similar reverse turn structure, although significant structural differences are seen at the rim of the binding crevice. Our results reflect both the common receptor-binding interface and the divergent biological functions of G proteins and arrestins.

Position of Transmembrane Helix 6 Determines Receptor G Protein Coupling Specificity

G protein coupled receptors (GPCRs) transmit extracellular signals into the cell by binding and activating different intracellular signaling proteins, such as G proteins (Gαβγ, families Gi, Gs, Gq, G12/13) or arrestins. To address the issue of Gs vs Gi coupling specificity, we carried out molecular dynamics simulations of lipid-embedded active β2-adrenoceptor (β2AR*) in complex with C-terminal peptides derived from the key interaction site of Gα (GαCT) as surrogate of Gαβγ. We find that GiαCT and GsαCT exploit distinct cytoplasmic receptor conformations that coexist in the uncomplexed β2AR*. The slim GiαCT stabilizes a β2AR* conformation, not accessible to the bulkier GsαCT, which requires a larger TM6 outward tilt for binding. Our results suggest that the TM6 conformational heterogeneity regulates the catalytic activity of β2AR* toward Gi or Gs.

Precision vs Flexibility in GPCR signaling

The G protein coupled receptor (GPCR) rhodopsin activates the heterotrimeric G protein transducin (Gt) to transmit the light signal into retinal rod cells. The rhodopsin activity is virtually zero in the dark and jumps by more than one billion fold after photon capture. Such perfect switching implies both high fidelity and speed of rhodopsin/Gt coupling. We employed Fourier transform infrared (FTIR) spectroscopy and supporting all-atom molecular dynamics (MD) simulations to study the conformational diversity of rhodopsin in membrane environment and extend the static picture provided by the available crystal structures. The FTIR results show how the equilibria of inactive and active protein states of the receptor (so-called metarhodopsin states) are regulated by the highly conserved E(D)RY and Yx7K(R) motives. The MD data identify an intrinsically unstructured cytoplasmic loop region connecting transmembrane helices 5 and 6 (CL3) and show how each protein state is split into conformational substates. The C-termini of the Gtγ- and Gtα-subunits (GαCT and GγCT), prepared as synthetic peptides, are likely to bind sequentially and at different sites of the active receptor. The peptides have different effects on the receptor conformation. While GγCT stabilizes the active states but preserves CL3 flexibility, GαCT selectively stabilizes a single conformational substate with largely helical CL3, as it is found in crystal structures. Based on these results we propose a mechanism for the fast and precise signal transfer from rhodopsin to Gt, which assumes a stepwise and mutual reduction of their conformational space. The mechanism relies on conserved amino acids and may therefore underlie GPCR/G protein coupling in general.

RHYTHM—a server to predict the orientation of transmembrane helices in channels and membrane-coils

RHYTHM is a web server that predicts buried versus exposed residues of helical membrane proteins. Starting from a given protein sequence, secondary and tertiary structure information is calculated by RHYTHM within only a few seconds. The prediction applies structural information from a growing data base of precalculated packing files and evolutionary information from sequence patterns conserved in a representative dataset of membrane proteins (‘Pfam-domains’). The program uses two types of position specific matrices to account for the different geometries of packing in channels and transporters (‘channels’) or other membrane proteins (‘membrane-coils’). The output provides information on the secondary structure and topology of the protein and specifically on the contact type of each residue and its conservation. This information can be downloaded as a graphical file for illustration, a text file for analysis and statistics and a PyMOL file for modeling purposes. The server can be freely accessed at: URL: http://proteinformatics.de/rhythm

Web-based molecular graphics for large complexes

The interactive visualization of very large macromolecular complexes on the web is becoming a challenging problem as experimental techniques advance at an unprecedented rate and deliver structures of increasing size. We have tackled this problem by introducing the binary and compressed Macromolecular Transmission Format (MMTF) to reduce network transfer and parsing time, and by developing NGL, a highly memory-efficient and scalable WebGL-based viewer. MMTF offers over 75% compression over the standard mmCIF format, is over an order of magnitude faster to parse, and contains additional information (e.g., bond information). NGL renders molecular complexes with millions of atoms interactively on desktop computers and smartphones alike, making it a tool of choice for web-based molecular visualization in research and education.

DeepSite: protein-binding site predictor using 3D-convolutional neural networks.

Motivation: An important step in structure‐based drug design consists in the prediction of druggable binding sites. Several algorithms for detecting binding cavities, those likely to bind to a small drug compound, have been developed over the years by clever exploitation of geometric, chemical and evolutionary features of the protein. Results: Here we present a novel knowledge‐based approach that uses state‐of‐the‐art convolutional neural networks, where the algorithm is learned by examples. In total, 7622 proteins from the scPDB database of binding sites have been evaluated using both a distance and a volumetric overlap approach. Our machine‐learning based method demonstrates superior performance to two other competitive algorithmic strategies. Availability and implementation: DeepSite is freely available at www.playmolecule.org. Users can submit either a PDB ID or PDB file for pocket detection to our NVIDIA GPU‐equipped servers through a WebGL graphical interface. Contact: gianni.defabritiis@upf.edu Supplementary information: Supplementary data are available at Bioinformatics online.

MP:PD—a data base of internal packing densities, internal packing defects and internal waters of helical membrane proteins

The membrane protein packing database (MP:PD) (http://proteinformatics.charite.de/mppd) is a database of helical membrane proteins featuring internal atomic packing densities, cavities and waters. Membrane proteins are not tightly packed but contain a considerable number of internal cavities that differ in volume, polarity and solvent accessibility as well as in their filling with internal water. Internal cavities are supposed to be regions of high physical compressibility. By serving as mobile hydrogen bonding donors or acceptors, internal waters likely facilitate transition between different functional states. Despite these distinct functional roles, internal cavities of helical membrane proteins are not well characterized, mainly because most internal waters are not resolved by crystal structure analysis. Here we combined various computational biophysical techniques to characterize internal cavities, reassign positions of internal waters and calculate internal packing densities of all available helical membrane protein structures and stored them in MP:PD. The database can be searched using keywords and entries can be downloaded. Each entry can be visualized in Provi, a Jmol-based protein viewer that provides an integrated display of low energy waters alongside membrane planes, internal packing density, hydrophobic cavities and hydrogen bonds.

SL2: an interactive webtool for modeling of missing segments in proteins

SuperLooper2 (SL2) (http://proteinformatics.charite.de/sl2) is the updated version of our previous web-server SuperLooper, a fragment based tool for the prediction and interactive placement of loop structures into globular and helical membrane proteins. In comparison to our previous version, SL2 benefits from both a considerably enlarged database of fragments derived from high-resolution 3D protein structures of globular and helical membrane proteins, and the integration of a new protein viewer. The database, now with double the content, significantly improved the coverage of fragment conformations and prediction quality. The employment of the NGL viewer for visualization of the protein under investigation and interactive selection of appropriate loops makes SL2 independent of third-party plug-ins and additional installations.