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.

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.

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.

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.

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.

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).

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.

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.

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.

Phenylalanine 368 of multidrug resistance-associated protein 4 (MRP4/ABCC4) plays a crucial role in substrate-specific transport activity.

Multidrug resistance-associated protein 4 (MRP4) is a membrane transporter that mediates the cellular efflux of a wide range of anionic drugs and endogenous molecules. MRP4 transport can influence the pharmacokinetics of drugs and their metabolites, therefore more knowledge about the molecular determinants important for its transport function would be of relevance. Here, we substituted amino acids Phe(368), Trp(995), and Arg(998) with conservative or non-conservative residues, and determined the effect on transport of the model substrates estradiol 17-β-d-glucuronide (E(2)17βG), cyclic guanosine monophosphate (cGMP), methotrexate (MTX), and folic acid into membrane vesicles isolated from baculovirus transduced HEK293 cells overexpressing the mutant MRP4 proteins.

Assignment of protonation states in proteins and ligands: combining pKa prediction with hydrogen bonding network optimization.

Among the many applications of molecular modeling, drug design is probably the one with the highest demands on the accuracy of the underlying structures. During lead optimization, the position of every atom in the binding site should ideally be known with high precision to identify those chemical modifications that are most likely to increase drug affinity. Unfortunately, X-ray crystallography at common resolution yields an electron density map that is too coarse, since the chemical elements and their protonation states cannot be fully resolved.

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.

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.

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.

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.

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.

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.

MYO15A (DFNB3) mutations in Turkish hearing loss families and functional modeling of a novel motor domain mutation.

Myosin XVA is an unconventional myosin which has been implicated in autosomal recessive nonsyndromic hearing impairment (ARNSHI) in humans. In Myo15A mouse models, vestibular dysfunction accompanies the autosomal recessive hearing loss. Genomewide homozygosity mapping and subsequent fine mapping in two Turkish families with ARNSHI revealed significant linkage to a critical interval harboring a known deafness gene MYO15A on chromosome 17p13.

MPP1 links the Usher protein network and the Crumbs protein complex in the retina.

The highly ordered distribution of neurons is an essential feature of a functional mammalian retina. Disruptions in the apico-basal polarity complexes at the outer limiting membrane (OLM) of the retina are associated with retinal patterning defects in vertebrates. We have analyzed the binding repertoire of MPP5/Pals1, a key member of the apico-basal Crumbs polarity complex, that has functionally conserved counterparts in zebrafish (nagie oko) and Drosophila (Stardust).

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.

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.

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.

Interaction of nephrocystin-4 and RPGRIP1 is disrupted by nephronophthisis or Leber congenital amaurosis-associated mutations.

RPGR-interacting protein 1 (RPGRIP1) is a key component of cone and rod photoreceptor cells, where it interacts with RPGR (retinitis pigmentosa GTPase regulator). Mutations in RPGRIP1 lead to autosomal recessive congenital blindness [Leber congenital amaurosis (LCA)]. Most LCA-associated missense mutations in RPGRIP1 are located in a segment that encodes two C2 domains.

Definition of a new information-based per-residue quality parameter.

For biomolecular NMR structures typically only a poor correspondence is observed between statistics derived from the experimental input data and structural quality indicators obtained from the structure ensembles. Here, we investigate the relationship between the amount of available NMR data and structure quality. By generating datasets with a predetermined information content and evaluating the quality of the resulting structure ensembles we show that there is, in contrast to previous findings, a linear relation between the information contained in experimental data and structural quality.

Reconstruction of the complete ouabain-binding pocket of Na,K-ATPase in gastric H,K-ATPase by substitution of only seven amino acids.

Although cardiac glycosides have been used as drugs for more than 2 centuries and their primary target, the sodium pump (Na,K-ATPase), has already been known for 4 decades, their exact binding site is still elusive. In our efforts to define the molecular basis of digitalis glycosides binding we started from the fact that a closely related enzyme, the gastric H,K-ATPase, does not bind glycosides like ouabain. Previously, we showed that a chimera of these two enzymes, in which only the M3-M4 and M5-M6 hairpins were of Na,K-ATPase, bound ouabain with high affinity (Koenderink, J.