YASARA Model-Interactive Molecular Modeling from Two Dimensions to Virtual Realities.

The industry's transition from three-dimensional (3D) glasses to virtual reality (VR) headsets has left modelers stranded without hardware supply, since walking around and waving arms in a virtual world is a great experience but also very tiring when doing time-intensive modeling work. We present a novel software implementation that uses a VR headset while sitting at a desk in front of the normal screen, which is beamed into the virtual reality together with keyboard, mouse, and chair using the headset's cameras and an extra tracker attached to the seat-back. Compared to 3D glasses, this yields a comparably relaxing but much more immersive workplace and provides additional possibilities such as taking molecules into one's hands, standing up, and walking or teleporting through the models. This VR functionality has been combined with a molecular graphics engine based on Vulkan, a next-generation cross-platform application programming interface (API) for GPUs and the successor of the widely used Open Graphics Library (OpenGL). It is built into the YASARA Model program, which includes many features like small and large molecule builders, electron densities, partial surfaces, contact analysis, coordinate manipulation, and animations. Interactive tutorials are provided to guide modelers into VR and familiarize them with the molecular modeling features. YASARA Model is available for Linux, Windows, Android, and MacOS (the latter without VR) with an introductory video at www.YASARA.org/vr.

Assembly of Biomolecular Gigastructures and Visualization with the Vulkan Graphics API.

Building and displaying all-atom models of biomolecular structures with millions or billions of atoms, like virus particles or cells, remain a challenge due to the sheer size of the data, the required levels of automated building, and the visualization limits of today's graphics hardware. Based on concepts introduced with the CellPack program, we report new algorithms to create such large-scale models using an intermediate coarse-grained "pet representation" of biomolecules with 1/10th the normal size. Pet atoms are placed such that they optimally trace the surface of the original molecule with just ∼1/50th the original atom number and are joined with covalent bonds. Molecular dynamics simulations of pet molecules allow for efficient packing optimization, as well as the generation of realistic DNA/RNA conformations. This pet world can be expanded back to the all-atom representation to be explored and visualized with full details. Essential for the efficient interactive visualization of gigastructures is the use of multiple levels of detail (LODs), where distant molecules are drawn with a heavily reduced polygon count. We present a grid-based algorithm to create such LODs for all common molecular graphics styles (including ball-and-sticks, ribbons, and cartoons) that do not require monochrome molecules to hide LOD transitions. As a practical application, we built all-atom models of SARS-CoV-2, HIV, and an entire presynaptic bouton with 1 µm diameter and 3.6 billion atoms, using modular building blocks to significantly reduce GPU memory requirements through instancing. We employ the Vulkan graphics API to maximize performance on consumer grade hardware and describe how to use the mmCIF format to efficiently store such giant models. An implementation is available as part of the YASARA molecular modeling and simulation program from www.YASARA.org. The free YASARA View program can be used to explore the presented models, which can be downloaded from www.YASARA.org/petworld, a Creative Commons platform for sharing giant biomolecular structures.

The drug diazaborine blocks ribosome biogenesis by inhibiting the AAA-ATPase Drg1.

The drug diazaborine is the only known inhibitor of ribosome biogenesis and specifically blocks large subunit formation in eukaryotic cells. However, the target of this drug and the mechanism of inhibition were unknown. Here we identify the AAA-ATPase Drg1 as a target of diazaborine. Inhibitor binding into the second AAA domain of Drg1 requires ATP loading and results in inhibition of ATP hydrolysis in this site. As a consequence the physiological activity of Drg1, i.e. the release of Rlp24 from pre-60S particles, is blocked, and further progression of cytoplasmic preribosome maturation is prevented. Our results identify the first target of an inhibitor of ribosome biogenesis and provide the mechanism of inhibition of a key step in large ribosomal subunit formation.

New ways to boost molecular dynamics simulations.

We describe a set of algorithms that allow to simulate dihydrofolate reductase (DHFR, a common benchmark) with the AMBER all-atom force field at 160 nanoseconds/day on a single Intel Core i7 5960X CPU (no graphics processing unit (GPU), 23,786 atoms, particle mesh Ewald (PME), 8.0 Å cutoff, correct atom masses, reproducible trajectory, CPU with 3.6 GHz, no turbo boost, 8 AVX registers). The new features include a mixed multiple time-step algorithm (reaching 5 fs), a tuned version of LINCS to constrain bond angles, the fusion of pair list creation and force calculation, pressure coupling with a "densostat," and exploitation of new CPU instruction sets like AVX2. The impact of Intel's new transactional memory, atomic instructions, and sloppy pair lists is also analyzed. The algorithms map well to GPUs and can automatically handle most Protein Data Bank (PDB) files including ligands. An implementation is available as part of the YASARA molecular modeling and simulation program from www.YASARA.org.

YASARA View - molecular graphics for all devices - from smartphones to workstations.

SUMMARY: Today's graphics processing units (GPUs) compose the scene from individual triangles. As about 320 triangles are needed to approximate a single sphere-an atom-in a convincing way, visualizing larger proteins with atomic details requires tens of millions of triangles, far too many for smooth interactive frame rates. We describe a new approach to solve this 'molecular graphics problem', which shares the work between GPU and multiple CPU cores, generates high-quality results with perfectly round spheres, shadows and ambient lighting and requires only OpenGL 1.0 functionality, without any pixel shader Z-buffer access (a feature which is missing in most mobile devices). AVAILABILITY AND IMPLEMENTATION: YASARA View, a molecular modeling program built around the visualization algorithm described here, is freely available (including commercial use) for Linux, MacOS, Windows and Android (Intel) from www.YASARA.org. CONTACT: elmar@yasara.org SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A series of PDB related databases for everyday needs.

The Protein Data Bank (PDB) is the world-wide repository of macromolecular structure information. We present a series of databases that run parallel to the PDB. Each database holds one entry, if possible, for each PDB entry. DSSP holds the secondary structure of the proteins. PDBREPORT holds reports on the structure quality and lists errors. HSSP holds a multiple sequence alignment for all proteins. The PDBFINDER holds easy to parse summaries of the PDB file content, augmented with essentials from the other systems. PDB_REDO holds re-refined, and often improved, copies of all structures solved by X-ray. WHY_NOT summarizes why certain files could not be produced. All these systems are updated weekly. The data sets can be used for the analysis of properties of protein structures in areas ranging from structural genomics, to cancer biology and protein design.

CING: an integrated residue-based structure validation program suite.

We present a suite of programs, named CING for Common Interface for NMR Structure Generation that provides for a residue-based, integrated validation of the structural NMR ensemble in conjunction with the experimental restraints and other input data. External validation programs and new internal validation routines compare the NMR-derived models with empirical data, measured chemical shifts, distance- and dihedral restraints and the results are visualized in a dynamic Web 2.0 report. A red-orange-green score is used for residues and restraints to direct the user to those critiques that warrant further investigation. Overall green scores below ~20 % accompanied by red scores over ~50 % are strongly indicative of poorly modelled structures. The publically accessible, secure iCing webserver ( https://nmr.le.ac.uk ) allows individual users to upload the NMR data and run a CING validation analysis.

Folding defects in P-type ATP 8B1 associated with hereditary cholestasis are ameliorated by 4-phenylbutyrate.

UNLABELLED: Deficiency in P-type ATP8B1 is a severe and clinically highly variable hereditary disorder that is primarily characterized by intrahepatic cholestasis. It presents either as a progressive (progressive familial intrahepatic cholestasis type 1 [PFIC1]) or intermittent (benign recurrent intrahepatic cholestasis type 1 [BRIC1]) disease. ATP8B1 deficiency is caused by autosomal recessive mutations in the gene encoding ATP8B1, a putative aminophospholipid-translocating P-type adenosine triphosphatase. The exact pathogenesis of the disease is elusive, and no effective pharmacological therapy is currently available. Here, the molecular consequences of six distinct ATP8B1 missense mutations (p.L127P, p.G308V, p.D454G, p.D554N, p.I661T, and p.G1040R) and one nonsense mutation (p.R1164X) associated with PFIC1 and/or BRIC1 were systematically characterized. Except for the p.L127P mutation, all mutations resulted in markedly reduced ATP8B1 protein expression, whereas messenger RNA expression was unaffected. Five of seven mutations resulted in (partial) retention of ATP8B1 in the endoplasmic reticulum. Reduced protein expression was partially restored by culturing the cells at 30 degrees C and by treatment with proteasomal inhibitors, indicating protein misfolding and subsequent proteosomal degradation. Protein misfolding was corroborated by predicting the consequences of most mutations onto a homology model of ATP8B1. Treatment with 4-phenylbutyrate, a clinically approved pharmacological chaperone, partially restored defects in expression and localization of ATP8B1 substitutions G308V, D454G, D554N, and in particular I661T, which is the most frequently identified mutation in BRIC1. CONCLUSION: A surprisingly large proportion of ATP8B1 mutations resulted in aberrant folding and decreased expression at the plasma membrane. These effects were partially restored by treatment with 4-phenylbutyrate. We propose that treatment with pharmacological chaperones may represent an effective therapeutic strategy to ameliorate the recurrent attacks of cholestasis in patients with intermittent (BRIC1) disease.

Reduced expression of ATP7B affected by Wilson disease-causing mutations is rescued by pharmacological folding chaperones 4-phenylbutyrate and curcumin.

UNLABELLED: Wilson disease (WD) is an autosomal recessive copper overload disorder of the liver and basal ganglia. WD is caused by mutations in the gene encoding ATP7B, a protein localized to the trans-Golgi network that primarily facilitates hepatic copper excretion. Current treatment comprises reduction of circulating copper by zinc supplementation or copper chelation. Despite treatment, a significant number of patients have neurological deterioration. The aim of this study was to investigate the possibility that defects arising from some WD mutations are ameliorated by drug treatment aimed at improvement of protein folding and restoration of protein function. This necessitated systematic characterization of the molecular consequences of distinct ATP7B missense mutations associated with WD. With the exception of p.S1363F, all mutations tested (p.G85V, p.R778L, p.H1069Q, p.C1104F, p.V1262F, p.G1343V, and p.S1363F) resulted in reduced ATP7B protein expression, whereas messenger RNA abundance was unaffected. Retention of mutant ATP7B in the endoplasmic reticulum, increased protein expression, and normalization of localization after culturing cells at 30 degrees C, and homology modeling suggested that these proteins were misfolded. Four distinct mutations exhibited residual copper export capacity, whereas other mutations resulted in complete disruption of copper export by ATP7B. Treatment with pharmacological chaperones 4-phenylbutyrate (4-PBA) and curcumin, a clinically approved compound, partially restored protein expression of most ATP7B mutants. CONCLUSION: These findings might enable novel treatment strategies in WD by directly enhancing the protein expression of mutant ATP7B with residual copper export activity. 1795.).

Homology modelling and spectroscopy, a never-ending love story.

Homology modelling is normally the technique of choice when experimental structure data are not available but three-dimensional coordinates are needed, for example, to aid with detailed interpretation of results of spectroscopic studies. Herein, the state of the art of homology modelling will be described in the light of a series of recent developments, and an overview will be given of the problems and opportunities encountered in this field. The major topic, the accuracy and precision of homology models, will be discussed extensively due to its influence on the reliability of conclusions drawn from the combination of homology models and spectroscopic data. Three real-world examples will illustrate how both homology modelling and spectroscopy can be beneficial for (bio)medical research.

Facilities that make the PDB data collection more powerful.

We describe a series of databases and tools that directly or indirectly support biomedical research on macromolecules, with focus on their applicability in protein structure bioinformatics research. DSSP, that determines secondary structures of proteins, has been updated to work well with extremely large structures in multiple formats. The PDBREPORT database that lists anomalies in protein structures has been remade to remove many small problems. These reports are now available as PDF-formatted files with a computer-readable summary. The VASE software has been added to analyze and visualize HSSP multiple sequence alignments for protein structures. The Lists collection of databases has been extended with a series of databases, most noticeably with a database that gives each protein structure a grade for usefulness in protein structure bioinformatics projects. The PDB-REDO collection of reanalyzed and re-refined protein structures that were solved by X-ray crystallography has been improved by dealing better with sugar residues and with hydrogen bonds, and adding many missing surface loops. All academic software underlying these protein structure bioinformatics applications and databases are now publicly accessible, either directly from the authors or from the GitHub software repository.

Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: Four approaches that performed well in CASP8.

A correct alignment is an essential requirement in homology modeling. Yet in order to bridge the structural gap between template and target, which may not only involve loop rearrangements, but also shifts of secondary structure elements and repacking of core residues, high-resolution refinement methods with full atomic details are needed. Here, we describe four approaches that address this "last mile of the protein folding problem" and have performed well during CASP8, yielding physically realistic models: YASARA, which runs molecular dynamics simulations of models in explicit solvent, using a new partly knowledge-based all atom force field derived from Amber, whose parameters have been optimized to minimize the damage done to protein crystal structures. The LEE-SERVER, which makes extensive use of conformational space annealing to create alignments, to help Modeller build physically realistic models while satisfying input restraints from templates and CHARMM stereochemistry, and to remodel the side-chains. ROSETTA, whose high resolution refinement protocol combines a physically realistic all atom force field with Monte Carlo minimization to allow the large conformational space to be sampled quickly. And finally UNDERTAKER, which creates a pool of candidate models from various templates and then optimizes them with an adaptive genetic algorithm, using a primarily empirical cost function that does not include bond angle, bond length, or other physics-like terms.

Molecular dynamics simulations of the chemokine CCL2 in complex with pull down-derived heparan sulfate hexasaccharides.

BACKGROUND: Binding of chemokines to glycosaminoglycans (GAGs) is a crucial step in leukocyte recruitment to inflamed tissues. METHODS: A disaccharide compositional analysis of the HS dp6 fraction in combination with MS analysis of the CCL2-depleted dp6 fraction was the basis for target GAG ligand structure suggestions. Four experimentally-derived heparan sulfate hexasaccharides, two potentially chemokine-specific and two unspecific, have been docked to CCL2. Subsequent 300 ns molecular dynamics simulations were used to improve the docked complexes. RESULTS: Hexasaccharides with four sulfations and no acetylations are suggested for selective and high affinity chemokine binding. Using the Antithromin-III/heparin complex as positive control for docking, we were able to recover the correct complex structure only if the previously liganded ATIII structure was used as input. Since the liganded structure is not known for a CCL2-GAG complex, we investigated if molecular dynamics simulations could improve initial docking results. We found that all four GAG oligosaccharides ended up in close contact with the known binding residues after about 100 ns simulation time. CONCLUSIONS: A discrimination of specific vs. unspecific CCL2 GAG ligands is not possible by this approach. Long-time molecular dynamics simulations are, however, well suited to capture the delicate enthalpy/entropy balance of GAG binding and improve results obtained from docking. GENERAL SIGNIFICANCE: With the comparison of two methods, MS-based ligand identification and molecular modelling, we have shown the current limitations of our molecular understanding of complex ligand binding which is could be due to the numerical inaccessibility of ligand-induced protein conformational changes.

Functional role of arginine 375 in transmembrane helix 6 of multidrug resistance protein 4 (MRP4/ABCC4).

Multidrug resistance protein (MRP) 4 transports a variety of endogenous and xenobiotic organic anions. MRP4 is widely expressed in the body and specifically localized to the renal apical proximal tubule cell membrane, where it mediates the excretion of these compounds into urine. To characterize the MRP4 substrate-binding site, the amino acids Phe368, Phe369, Glu374, Arg375, and Glu378 of transmembrane helix 6, and Arg998 of helix 12, localized in the intracellular half of the central pore, were mutated into the corresponding amino acids of MRP1 and MRP2. Membrane vesicles isolated from human embryonic kidney 293 cells overexpressing these mutants showed significantly reduced methotrexate (MTX) and cGMP transport activity compared with vesicles that expressed wild-type MRP4. The only exception was substitution of Arg375 with serine, which had no effect on cGMP transport but significantly decreased the affinity of MTX. Substitution of the same amino acid with a positively charged lysine returned the MTX affinity to that of the wild type. Furthermore, MTX inhibition of MRP4-mediated cGMP transport was noncompetitive, and the inhibition constant was increased by introduction of the R375S mutation. A homology model of MRP4 showed that Arg375 and Arg998 face right into the central aqueous pore of MRP4. We conclude that positively charged amino acids in transmembrane helices 6 and 12 contribute to the MRP4 substrate-binding pocket.

GTP-Ras disrupts the intramolecular complex of C1 and RA domains of Nore1.

The novel Ras effector mNore1, capable of inducing apoptosis, is a multidomain protein. It comprises a C1 domain homologous to PKC and an RA domain similar to the Ras effectors AF-6 and RalGDS. Here, we determine the affinity of these two domains to the active forms of Ras and Rap1 using isothermal calorimetric titration. The interaction of Ras/Rap1-GTP with the RA domain of mNore1 is weakened significantly by direct binding of the C1 domain to the RA domain. In order to analyze this observation in atomic detail, we solved the C1 solution structure by NMR. By determining chemical shifts and relaxation rates, we can show an intramolecular complex of C1-RA. GTP-Ras titration and binding to RA disrupts this complex and displaces the C1 domain. Once the C1 domain tumbles freely in solution, a lipid binding interface becomes accessible. Furthermore, we provide evidence of phosphatidylinositol 3-phosphate binding of the free C1 domain.

The human Vps29 retromer component is a metallo-phosphoesterase for a cation-independent mannose 6-phosphate receptor substrate peptide.

The retromer complex is involved in the retrograde transport of the CI-M6PR (cation-independent mannose 6-phosphate receptor) from endosomes to the Golgi. It is a hetero-trimeric complex composed of Vps26 (vacuolar sorting protein 26), Vps29 and Vps35 proteins, which are conserved in eukaryote evolution. Recently, elucidation of the crystal structure of Vps29 revealed that Vps29 contains a metallo-phosphoesterase fold [Wang, Guo, Liang, Fan, Zhu, Zang, Zhu, Li, Teng, Niu et al. (2005) J. Biol. Chem. 280, 22962-22967; Collins, Skinner, Watson, Seaman and Owen (2005) Nat. Struct. Mol. Biol. 12, 594-602]. We demonstrate that recombinant hVps29 (human Vps29) displays in vitro phosphatase activity towards a serine-phosphorylated peptide, containing the acidic-cluster dileucine motif of the cytoplasmatic tail of the CI-M6PR. Efficient dephosphorylation required the additional presence of recombinant hVps26 and hVps35 proteins, which interact with hVps29. Phosphatase activity of hVps29 was greatly decreased by alanine substitutions of active-site residues that are predicted to co-ordinate metal ions. Using inductively coupled plasma MS, we demonstrate that recombinant hVps29 binds zinc. Moreover, hVps29-dependent phosphatase activity is greatly reduced by non-specific and zinc-specific metal ion chelators, which can be completely restored by addition of excess ZnCl2. The binuclear Zn2+ centre and phosphate group were modelled into the hVps29 catalytic site and pKa calculations provided further insight into the molecular mechanisms of Vps29 phosphatase activity. We conclude that the retromer complex displays Vps29-dependent in vitro phosphatase activity towards a serinephosphorylated acidic-cluster dileucine motif that is involved in endosomal trafficking of the CI-M6PR. The potential significance of these findings with respect to regulation of transport of cycling trans-Golgi network proteins is discussed.

A conserved inter-domain communication mechanism regulates the ATPase activity of the AAA-protein Drg1.

AAA-ATPases fulfil essential roles in different cellular pathways and often act in form of hexameric complexes. Interaction with pathway-specific substrate and adaptor proteins recruits them to their targets and modulates their catalytic activity. This substrate dependent regulation of ATP hydrolysis in the AAA-domains is mediated by a non-catalytic N-terminal domain. The exact mechanisms that transmit the signal from the N-domain and coordinate the individual AAA-domains in the hexameric complex are still the topic of intensive research. Here, we present the characterization of a novel mutant variant of the eukaryotic AAA-ATPase Drg1 that shows dysregulation of ATPase activity and altered interaction with Rlp24, its substrate in ribosome biogenesis. This defective regulation is the consequence of amino acid exchanges at the interface between the regulatory N-domain and the adjacent D1 AAA-domain. The effects caused by these mutations strongly resemble those of pathological mutations of the AAA-ATPase p97 which cause the hereditary proteinopathy IBMPFD (inclusion body myopathy associated with Paget's disease of the bone and frontotemporal dementia). Our results therefore suggest well conserved mechanisms of regulation between structurally, but not functionally related members of the AAA-family.

Fast empirical pKa prediction by Ewald summation.

pK(a) calculations for macromolecules are normally performed by solving the Poisson-Boltzmann equation, accounting for the different dielectric constants of solvent and solute, as well as the ionic strength. Despite the large number of successful applications, there are some situations where the current algorithms are not suitable: (1) large scale, high-throughput analysis which requires calculations to be completed within a fraction of a second, e.g. when permanently monitoring pK(a) shifts during a molecular dynamics simulation; (2) prediction of pK(a)s in periodic boundaries, e.g. when reconstructing entire protein crystal unit cells from PDB files, including the correct protonation patterns at experimental pH. Such in silico crystals are needed by 'self-parameterizing' molecular dynamics force fields like YASARA YAMBER, that optimize their parameters while energy-minimizing high-resolution protein crystals. To address both problems, we define an empirical equation that expresses the pK(a) as a function of electrostatic potential, hydrogen bonds and accessible surface area. The electrostatic potential is evaluated by Ewald summation, which captures periodic crystal environments and the uncertainty in atom positions using Gaussian charge densities. The empirical proportionality constants are derived from 217 experimentally determined pK(a)s, and despite its simplicity, this pK(a) calculation method reaches a high overall jack-knifed accuracy, and is fast enough to be used during a molecular dynamics simulation. A reliable null-model to judge pK(a) prediction accuracies is also presented.

MPP5 recruits MPP4 to the CRB1 complex in photoreceptors.

PURPOSE: Mutations in the human Crumbs homologue 1 (CRB1) gene are a frequent cause of Leber congenital amaurosis (LCA) and various forms of retinitis pigmentosa. CRB1 is thought to organize an intracellular protein scaffold in the retina that is involved in photoreceptor polarity. This study was focused on the identification, subcellular localization, and binding characteristics of a novel member of the protein scaffold connected to CRB1. METHODS: To dissect the protein scaffold connected to CRB1, the yeast two-hybrid approach was used to screen for interacting proteins. Glutathione S-transferase (GST) pull-down analysis and immunoprecipitation were used to verify protein-protein interactions. The subcellular localization of the proteins was visualized by immunohistochemistry and confocal microscopy on human retinas and immunoelectron microscopy on mouse retinas. RESULTS: A novel member of the scaffold connected to CRB1, called membrane palmitoylated protein (MPP) subfamily member 4 (MPP4), a membrane-associated guanylate kinase (MAGUK) protein, was identified. MPP4 was found to exist in a complex with CRB1 through direct interaction with the MPP subfamily member MPP5 (PALS1). 3D homology modeling provided evidence for a mechanism that regulates the recruitment of both homo- and heterodimers of MPP4 and -5 proteins to the complex. Localization studies in the retina showed that CRB1, MPP5, and MPP4 colocalize at the outer limiting membrane (OLM). CONCLUSIONS: These data imply that MPP4 and -5 have a role in photoreceptor polarity and, by association with CRB1, pinpoint the cognate genes as functional candidate genes for inherited retinopathies.

Increasing the precision of comparative models with YASARA NOVA--a self-parameterizing force field.

One of the conclusions drawn at the CASP4 meeting in Asilomar was that applying various force fields during refinement of template-based models tends to move predictions in the wrong direction, away from the experimentally determined coordinates. We have derived an all-atom force field aimed at protein and nucleotide optimization in vacuo (NOVA), which has been specifically designed to avoid this problem. NOVA resembles common molecular dynamics force fields but has been automatically parameterized with two major goals: (i) not to make high resolution X-ray structures worse and (ii) to improve homology models built by WHAT IF. Force-field parameters were not required to be physically correct; instead, they were optimized with random Monte Carlo moves in force-field parameter space, each one evaluated by simulated annealing runs of a 50-protein optimization set. Errors inherent to the approximate force-field equation could thus be canceled by errors in force-field parameters. Compared with the optimization set, the force field did equally well on an independent validation set and is shown to move in silico models closer to reality. It can be applied to modeling applications as well as X-ray and NMR structure refinement. A new method to assign force-field parameters based on molecular trees is also presented. A NOVA server is freely accessible at http://www.yasara.com/servers