Elmar Krieger

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 Proteins 2002;47:393–402.

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. Proteins 2009.

Making Optimal Use of Empirical Energy Functions: Force-Field Parameterization in Crystal Space

Todays energy functions are not able yet to distinguish reliably between correct and almost correct protein models. Improving these near‐native models is currently a major bottle‐neck in homology modeling or experimental structure determination at low resolution. Increasingly accurate energy functions are required to complete the “last mile of the protein folding problem,” for example during a molecular dynamics simulation. We present a new approach to reach this goal. For 50 high resolution X‐ray structures, the complete unit cell was reconstructed, including disordered water molecules, counter ions, and hydrogen atoms. Simulations were then run at the pH at which the crystal was solved, while force‐field parameters were iteratively adjusted so that the damage done to the structures was minimal. Starting with initial parameters from the AMBER force field, the optimization procedure converged at a new force field called YAMBER (Yet Another Model Building and Energy Refinement force field), which is shown to do significantly less damage to X‐ray structures, often move homology models in the right direction, and occasionally make them look like experimental structures. Application of YAMBER during the CASP5 structure prediction experiment yielded a model for target 176 that was ranked first among 150 submissions. Due to its compatibility with the well‐established AMBER format, YAMBER can be used by almost any molecular dynamics program. The parameters are freely available from www.yasara.org/yamber. Proteins 2004.

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.

YASARA View—molecular graphics for all devices—from smartphones to workstations

Summary: Todays 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.

Identification of different binding sites in the dendritic cell-specific receptor DC-SIGN for intercellular adhesion molecule 3 and HIV-1

The novel dendritic cell (DC)-specific human immunodeficiency virus type 1 (HIV-1) receptor DC-SIGN plays a key role in the dissemination of HIV-1 by DC. DC-SIGN is thought to capture HIV-1 at mucosal sites of entry, facilitating transport to lymphoid tissues, where DC-SIGN efficiently transmits HIV-1 to T cells. DC-SIGN is also important in the initiation of immune responses by regulating DC-T cell interactions through intercellular adhesion molecule 3 (ICAM-3). We have characterized the mechanism of ligand binding by DC-SIGN and identified the crucial amino acids involved in this process. Strikingly, the HIV-1 gp120 binding site in DC-SIGN is different from that of ICAM-3, consistent with the observation that glycosylation of gp120, in contrast to ICAM-3, is not crucial to the interaction with DC-SIGN. A specific mutation in DC-SIGN abrogated ICAM-3 binding, whereas the HIV-1 gp120 interaction was unaffected. This DC-SIGN mutant captured HIV-1 and infected T cells intrans as efficiently as wild-type DC-SIGN, demonstrating that ICAM-3 binding is not necessary for HIV-1 transmission. This study provides a basis for the design of drugs that inhibit or alter interactions of DC-SIGN with gp120 but not with ICAM-3 or vice versa and that have a therapeutic value in immunological diseases and/or HIV-1 infections.

Models@Home: distributed computing in bioinformatics using a screensaver based approach

MOTIVATION Due to the steadily growing computational demands in bioinformatics and related scientific disciplines, one is forced to make optimal use of the available resources. A straightforward solution is to build a network of idle computers and let each of them work on a small piece of a scientific challenge, as done by Seti@Home (http://setiathome.berkeley.edu), the worlds largest distributed computing project. RESULTS We developed a generally applicable distributed computing solution that uses a screensaver system similar to Seti@Home. The software exploits the coarse-grained nature of typical bioinformatics projects. Three major considerations for the design were: (1) often, many different programs are needed, while the time is lacking to parallelize them. Models@Home can run any program in parallel without modifications to the source code; (2) in contrast to the Seti project, bioinformatics applications are normally more sensitive to lost jobs. Models@Home therefore includes stringent control over job scheduling; (3) to allow use in heterogeneous environments, Linux and Windows based workstations can be combined with dedicated PCs to build a homogeneous cluster. We present three practical applications of Models@Home, running the modeling programs WHAT IF and YASARA on 30 PCs: force field parameterization, molecular dynamics docking, and database maintenance.

A mutation in the gamma actin 1 (ACTG1) gene causes autosomal dominant hearing loss (DFNA20/26).

Linkage analysis in a multigenerational family with autosomal dominant hearing loss yielded a chromosomal localisation of the underlying genetic defect in the DFNA20/26 locus at 17q25-qter. The 6-cM critical region harboured the γ-1-actin (ACTG1) gene, which was considered an attractive candidate gene because actins are important structural elements of the inner ear hair cells. In this study, a Thr278Ile mutation was identified in helix 9 of the modelled protein structure. The alteration of residue Thr278 is predicted to have a small but significant effect on the γ 1 actin structure owing to its close proximity to a methionine residue at position 313 in helix 11. Met313 has no space in the structure to move away. Moreover, the Thr278 residue is highly conserved throughout eukaryotic evolution. Using a known actin structure the mutation could be predicted to impair actin polymerisation. These findings strongly suggest that the Thr278Ile mutation in ACTG1 represents the first disease causing germline mutation in a cytoplasmic actin isoform.

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

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

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°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. (HEPATOLOGY 2009;50:1783–1795.)