Four small puzzles that Rosetta doesn't solve

A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on deceptively small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents extensive Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular "puzzles" should serve as useful model systems for developers wishing to make foundational improvements to this powerful modeling suite.Comment: Published in PLoS One as a manuscript for the RosettaCon 2010 Special Collectio

A New Heterobinuclear FeIIICuII Complex with a Single Terminal FeIII–O(phenolate) Bond. Relevance to Purple Acid Phosphatases and Nucleases

A novel heterobinuclear mixed valence complex [Fe^IIICu^II(BPBPMP)(OAc)_2]ClO_4, 1, with the unsymmetrical N_5O_2 donor ligand 2-bis[{(2-pyridylmethyl)aminomethyl}-6-{(2-hydroxybenzyl)(2-pyridylmethyl)} aminomethyl]-4-methylphenol (H_2BPBPMP) has been synthesized and characterized. A combination of data from mass spectrometry, potentiometric titrations, X-ray absorption and electron paramagnetic resonance spectroscopy, as well as kinetics measurements indicates that in ethanol/water solutions an [Fe^III-(nu)OH-Cu^IIOH_2]+ species is generated which is the likely catalyst for 2,4-bis(dinitrophenyl)phosphate and DNA hydrolysis. Insofar as the data are consistent with the presence of an Fe_III-bound hydroxide acting as a nucleophile during catalysis, 1 presents a suitable mimic for the hydrolytic enzyme purple acid phosphatase. Notably, 1 is significantly more reactive than its isostructural homologues with different metal composition (Fe^IIIM^II, where M^II is Zn^II, Mn^II, Ni^II,or Fe^II). Of particular interest is the observation that cleavage of double-stranded plasmid DNA occurs even at very low concentrations of 1 (2.5 nuM), under physiological conditions (optimum pH of 7.0), with a rate enhancement of 2.7 x 10^7 over the uncatalyzed reaction. Thus, 1 is one of the most effective model complexes to date, mimicking the function of nucleases

Application of the PM6 semi-empirical method to modeling proteins enhances docking accuracy of AutoDock

<p>Abstract</p> <p>Background</p> <p>Molecular docking methods are commonly used for predicting binding modes and energies of ligands to proteins. For accurate complex geometry and binding energy estimation, an appropriate method for calculating partial charges is essential. AutoDockTools software, the interface for preparing input files for one of the most widely used docking programs AutoDock 4, utilizes the Gasteiger partial charge calculation method for both protein and ligand charge calculation. However, it has already been shown that more accurate partial charge calculation - and as a consequence, more accurate docking- can be achieved by using quantum chemical methods. For docking calculations quantum chemical partial charge calculation as a routine was only used for ligands so far. The newly developed Mozyme function of MOPAC2009 allows fast partial charge calculation of proteins by quantum mechanical semi-empirical methods. Thus, in the current study, the effect of semi-empirical quantum-mechanical partial charge calculation on docking accuracy could be investigated.</p> <p>Results</p> <p>The docking accuracy of AutoDock 4 using the original AutoDock scoring function was investigated on a set of 53 protein ligand complexes using Gasteiger and PM6 partial charge calculation methods. This has enabled us to compare the effect of the partial charge calculation method on docking accuracy utilizing AutoDock 4 software. Our results showed that the docking accuracy in regard to complex geometry (docking result defined as accurate when the RMSD of the first rank docking result complex is within 2 Å of the experimentally determined X-ray structure) significantly increased when partial charges of the ligands and proteins were calculated with the semi-empirical PM6 method.</p> <p>Out of the 53 complexes analyzed in the course of our study, the geometry of 42 complexes were accurately calculated using PM6 partial charges, while the use of Gasteiger charges resulted in only 28 accurate geometries. The binding affinity estimation was not influenced by the partial charge calculation method - for more accurate binding affinity prediction development of a new scoring function for AutoDock is needed.</p> <p>Conclusion</p> <p>Our results demonstrate that the accuracy of determination of complex geometry using AutoDock 4 for docking calculation greatly increases with the use of quantum chemical partial charge calculation on both the ligands and proteins.</p

Identification and molecular modeling of a novel, plant-like, human purple acid phosphatase

Purple acid phosphatases are a family of binuclear metallohydrolases that have been identified in plants, animals and fungi. Only one isoform of ~35 kDa has been isolated from animals, where it is associated with bone resorption and microbial killing through its phosphatase activity, and hydroxyl radical production, respectively. Using the sensitive PSI-BLAST search method, sequences representing new purple acid phosphatase-like proteins have been identified in mammals, insects and nematodes. These new putative isoforms are closely related to the ~55 kDa purple acid phosphatase characterized from plants. Secondary structure prediction of the new human isoform further confirms its similarity to a purple acid phosphatase from the red kidney bean. A structural model for the human enzyme was constructed based on the red kidney bean purple acid phosphatase structure. This model shows that the catalytic centre observed in other purple acid phosphatases is also present in this new isoform. These observations suggest that the sequences identified in this study represent a novel subfamily of plant-like purple acid phosphatases in animals and humans

Three-dimensional structure of an Fv from a human IgM immunoglobulin

An IgM(κ) immunoglobulin from a patient (Pot) with Waldenstrom's macroglobulinemia was hydrolyzed with pepsin to release a fragment consisting of the 'variable' (V) domains of the light and heavy chains plus eight residue 'tails' from the 'constant' (C) domains. The crystal structure of this fragment was determined at 2.3 Å resolution by molecular replacement and crystallographic refinement methods. When examined separately, the light chain component closely resembles another human κ chain (Rei) in both the β-pleated sheet regions and the 'hypervariable' loops. The conserved pleated sheets in the heavy chain are similar to those in the human Kol IgG1 protein, but the third hypervariable loop in particular is different from that in any immunoglobulin structure described to date. As in the Kol protein, this loop blocks the access to any internal active site along the light-heavy chain interface. Unlike the Kol protein, however, the loop does not protrude beyond the boundaries of a conventional antigen combining site. Instead, it forms a very compact structure, which fills almost all residual space between the domains. This is an example of one dominant complementarity-determining region (CDR) essentially negating the diversity possible with five other CDRs in the two chains. Ordered water molecules are associated with light chain constituents along the interface, but not with CDR3 of the heavy chain. In screening exercises the Pot IgM failed to bind a wide variety of peptides. Together, the results suggest that ligand binding can only occur on external surfaces of the protein. These surfaces carry a limited number of side chains usually assigned to CDRs in more typical antibodies. © 1992

Identification of a non-purple tartrate-resistant acid phosphatase: an evolutionary link to Ser/Thr protein phosphatases?

Background: Tartrate-resistant acid phosphatases (TRAcPs), also known as purple acid phosphatases (PAPs), are a family of binuclear metallohydrolases that have been identified in plants, animals and fungi. The human enzyme is a major histochemical marker for the diagnosis of bone-related diseases. TRAcPs can occur as a small form possessing only the ~35 kDa catalytic domain, or a larger ~55 kDa form possessing both a catalytic domain and an additional N-terminal domain of unknown function. Due to its role in bone resorption the 35 kDa TRAcP has become a promising target for the development of anti-osteoporotic chemotherapeutics. Findings: A new human gene product encoding a metallohydrolase distantly related to the ~55 kDa plant TRAcP was identified and characterised. The gene product is found in a number of animal species, and is present in all tissues sampled by the RIKEN mouse transcriptome project. Construction of a homology model illustrated that six of the seven metal-coordinating ligands in the active site are identical to that observed in the TRAcP family. However, the tyrosine ligand associated with the charge transfer transition and purple color of TRAcPs is replaced by a histidine. Conclusion: The gene product identified here may represent an evolutionary link between TRAcPs and Ser/Thr protein phosphatases. Its biological function is currently unknown but is unlikely to be associated with bone metabolism

An electron transfer path connects subunits of a mycobacterial respiratory supercomplex

We report a 3.5-angstrom-resolution cryo-electron microscopy structure of a respiratory supercomplex isolated from Mycobacterium smegmatis. It comprises a complex III dimer flanked on either side by individual complex IV subunits. Complex III and IV associate so that electrons can be transferred from quinol in complex III to the oxygen reduction center in complex IV by way of a bridging cytochrome subunit. We observed a superoxide dismutase-like subunit at the periplasmic face, which may be responsible for detoxification of superoxide formed by complex III. The structure reveals features of an established drug target and provides a foundation for the development of treatments for human tuberculosis