Selective refinement and selection of near‐native models in protein structure prediction

In recent years in silico protein structure prediction reached a level where fully automated servers can generate large pools of near‐native structures. However, the identification and further refinement of the best structures from the pool of models remain problematic. To address these issues, we have developed (i) a target‐specific selective refinement (SR) protocol; and (ii) molecular dynamics (MD) simulation based ranking (SMDR) method. In SR the all‐atom refinement of structures is accomplished via the Rosetta Relax protocol, subject to specific constraints determined by the size and complexity of the target. The best‐refined models are selected with SMDR by testing their relative stability against gradual heating through all‐atom MD simulations. Through extensive testing we have found that Mufold‐MD, our fully automated protein structure prediction server updated with the SR and SMDR modules consistently outperformed its previous versions. Proteins 2015; 83:1823–1835. © 2015 Wiley Periodicals, Inc.     

Statistical torsion angle potential energy functions for protein structure modeling: A bicubic interpolation approach

A set of grid type knowledge‐based energy functions is introduced for ?–χ₁, ψ–χ₁, ?–ψ, and χ₁–χ₂ torsion angle combinations. Boltzmann distribution is assumed for the torsion angle populations from protein X‐ray structures, and the functions are named as statistical torsion angle potential energy functions. The grid points around periodic boundaries are duplicated to force periodicity, and the remedy relieves the derivative discontinuity problem. The devised functions rapidly improve the quality of model structures. The potential bias in the functions and the usefulness of additional secondary structure information are also investigated. The proposed guiding functions are expected to facilitate protein structure modeling, such as protein structure prediction, protein design, and structure refinement. Proteins 2013. Proteins 2013; 81:1156–1165. © 2013 Wiley Periodicals, Inc.     

On the reliability of peptide nonplanarity seen in ultra‐high resolution crystal structures

Ultra‐high resolution protein crystal structures have been considered as relatively reliable sources for defining details of protein geometry, such as the extent to which the peptide unit deviates from planarity. Chellapa and Rose (Proteins 2015; 83:1687) recently called this into question, reporting that for a dozen representative protein structures determined at ～1 Å resolution, the diffraction data could be equally well fit with models restrained to have highly planar peptides, i.e. having a standard deviation of the ω torsion angles of only ～1° instead of the typically observed value of ～6°. Here, we document both conceptual and practical shortcomings of that study and show that the more tightly restrained models are demonstrably incorrect and do not fit the diffraction data equally well. We emphasize the importance of inspecting electron density maps when investigating the agreement between a model and its experimental data. Overall, this report reinforces that modern standard refinement protocols have been well‐conceived and that ultra‐high resolution protein crystal structures, when evaluated carefully and used with an awareness of their levels of coordinate uncertainty, are powerful sources of information for providing reliable information about the details of protein geometry.     

Accurate geometries for “Mountain pass” regions of the Ramachandran plot using quantum chemical calculations

Unusual local arrangements of protein in Ramachandran space are not well represented by standard geometry tools used in either protein structure refinement using simple harmonic geometry restraints or in protein simulations using molecular mechanics force fields. In contrast, quantum chemical computations using small poly‐peptide molecular models can predict accurate geometries for any well‐defined backbone Ramachandran orientation. For conformations along transition regions—ϕ from −60 to 60°—a very good agreement with representative high‐resolution experimental X‐ray (≤1.5 Å) protein structures is obtained for both backbone C⁻¹‐N‐C_α angle and the nonbonded O⁻¹…C distance, while “standard geometry” leads to the “clashing” of O…C atoms and Amber FF99SB predicts distances too large by about 0.15 Å. These results confirm that quantum chemistry computations add valuable support for detailed analysis of local structural arrangements in proteins, providing improved or missing data for less understood high‐energy or unusual regions.     

Modeling conformational redox‐switch modulation of human succinic semialdehyde dehydrogenase

Succinic semialdehyde dehydrogenase (SSADH) converts succinic semialdehyde (SSA) to succinic acid in the mitochondrial matrix and is involved in the metabolism of the inhibitory neurotransmitter γ‐aminobutyric acid (GABA). The molecular structure of human SSADH revealed the intrinsic regulatory mechanism—redox‐switch modulation—by which large conformational changes are brought about in the catalytic loop through disulfide bonding. The crystal structures revealed two SSADH conformations, and computational modeling of transformation between them can provide substantial insights into detailed dynamic redox modulation. On the basis of these two clear crystal structures, we modeled the conformational motion between these structures in silico. For that purpose, we proposed and used a geometry‐based coarse‐grained mathematical model of long‐range protein motion and the related modeling algorithm. The algorithm is based on solving the special optimization problem, which is similar to the classical Monge–Kantorovich mass transportation problem. The modeled transformation was supported by another morphing method based on a completely different framework. The result of the modeling facilitates better interpretation and understanding of the SSADH biological role. Proteins 2015; 83:2217–2229. © 2015 Wiley Periodicals, Inc.     

ω‐Turn: A novel β‐turn mimic in globular proteins stabilized by main‐chain to side‐chain CH···O interaction

Mimicry of structural motifs is a common feature in proteins. The 10‐membered hydrogen‐bonded ring involving the main‐chain C?O in a β‐turn can be formed using a side‐chain carbonyl group leading to Asx‐turn. We show that the N? H component of hydrogen bond can be replaced by a C^γ‐H group in the side chain, culminating in a nonconventional C? H···O interaction. Because of its shape this β‐turn mimic is designated as ω‐turn, which is found to occur ～three times per 100 residues. Three residues (i to i + 2) constitute the turn with the C? H···O interaction occurring between the terminal residues, constraining the torsion angles ?_{i + 1}, ψ_{i + 1}, ?_{i + 2} and χ′_{1(i + 2)} (using the interacting C^γ atom). Based on these angles there are two types of ω‐turns, each of which can be further divided into two groups. C^β‐branched side‐chains, and Met and Gln have high propensities to occur at i + 2; for the last two residues the carbonyl oxygen may participate in an additional interaction involving the S and amino group, respectively. With Cys occupying the i + 1 position, such turns are found in the metal‐binding sites. N‐linked glycosylation occurs at the consensus pattern Asn‐Xaa‐Ser/Thr; with Thr at i + 2, the sequence can adopt the secondary structure of a ω‐turn, which may be the recognition site for protein modification. Location between two β‐strands is the most common occurrence in protein tertiary structure, and being generally exposed ω‐turn may constitute the antigenic determinant site. It is a stable scaffold and may be used in protein engineering and peptide design. Proteins 2015; 83:203–214. © 2014 Wiley Periodicals, Inc.     

Crystal structure of the PepSY‐containing domain of the YpeB protein involved in germination of bacillus spores

The crystal structure of the C‐terminal domain of the Bacillus megaterium YpeB protein has been solved by X‐ray crystallography to 1.80‐Å resolution. The full‐length protein is essential in stabilising the SleB cortex lytic enzyme in Bacillus spores, and may have a role in regulating SleB activity during spore germination. The YpeB‐C crystal structure comprises three tandemly repeated PepSY domains, which are aligned to form an extended laterally compressed molecule. A predominantly positively charged region located in the second PepSY domain may provide a site for protein interactions that are important in stabilising SleB and YpeB within the spore. Proteins 2015; 83:1914–1921. © 2015 Wiley Periodicals, Inc.     

The effect of calcium on the conformation of cobalamin transporter BtuB

BtuB is a β‐barrel membrane protein that facilitates transport of cobalamin (vitamin B12) from the extracellular medium across the outer membrane of Escherichia coli. It is thought that binding of B12 to BtuB alters the conformation of its periplasm‐exposed N‐terminal residues (the TonB box), which enables subsequent binding of a TonB protein and leads to eventual uptake of B12 into the cytoplasm. Structural studies determined the location of the B12 binding site at the top of the BtuB's β‐barrel, surrounded by extracellular loops. However, the structure of the loops was found to depend on the method used to obtain the protein crystals, which—among other factors—differed in calcium concentration. Experimentally, calcium concentration was found to modulate the binding of the B12 substrate to BtuB. In this study, we investigate the effect of calcium ions on the conformation of the extracellular loops of BtuB and their possible role in B12 binding. Using all‐atom molecular dynamics, we simulate conformational fluctuations of several X‐ray structures of BtuB in the presence and absence of calcium ions. These simulations demonstrate that calcium ions can stabilize the conformation of loops 3–4, 5–6, and 15–16, and thereby prevent occlusion of the binding site. Furthermore, binding of calcium ions to extracellular loops of BtuB was found to enhance correlated motions in the BtuB structure, which is expected to promote signal transduction. Finally, we characterize conformation dynamics of the TonB box in different X‐ray structures and find an interesting correlation between the stability of the TonB box structure and calcium binding. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Heteromerous Bistricyclic Aromatic Enes: Dynamic Stereochemistry of Xanthylidene‐Anthrones

The heteromerous bistricyclic aromatic ene (BAE) 2,2'‐dimethyl‐10‐(9H‐xanthylidene)‐9(10H)‐anthrone (DMXA) was synthesized by a condensation of 10,10‐dichloro‐2‐methylxanthene with 2‐methylanthrone. X‐ray crystallography of (E)‐DMXA and xanthylidene‐anthrone (XA) indicated that the molecules adopt anti‐folded conformations with folding dihedral angles of 44°/44° and 39°/41°, respectively . The crystal structure of anti‐folded (E)‐DMXA does not indicate any xanthenylium–anthracenolate push–pull effect. E,Z‐diastereomerization of DMXA was studied by ¹H‐NMR coalescence‐temperature measurements at different magnetic field strengths and by kinetic equilibration experiments . Free energy of activation for this process was 81.5 (±1.3) kJ/mol. B3LYP/6‐311+G(d,p) calculations showed that anti‐folded conformers of XA, (E)‐DMXA, bianthrone (AA), and dixanthylene (XX) were global minima. The twisted conformers of XA, AA, and XX were local minima (ΔG₂₉₈ = 16, 18, and 24 kJ/mol) with a substantial dipolar xanthenylium–anthracenolate dipolar contribution for XA. Chirality 27:919–928, 2015. © 2015 Wiley Periodicals, Inc.     

Synthetic peptides mimicking the binding site of human acetylcholinesterase for its inhibitor fasciculin 2

Molecules capable of mimicking protein binding and/or functional sites present useful tools for a range of biomedical applications, including the inhibition of protein–ligand interactions. Such mimics of protein binding sites can currently be generated through structure‐based design and chemical synthesis. Computational protein design could be further used to optimize protein binding site mimetics through rationally designed mutations that improve intermolecular interactions or peptide stability. Here, as a model for the study, we chose an interaction between human acetylcholinesterase (hAChE) and its inhibitor fasciculin‐2 (Fas) because the structure and function of this complex is well understood. Structure‐based design of mimics of the hAChE binding site for Fas yielded a peptide that binds to Fas at micromolar concentrations. Replacement of hAChE residues known to be essential for its interaction with Fas with alanine, in this peptide, resulted in almost complete loss of binding to Fas. Computational optimization of the hAChE mimetic peptide yielded a variant with slightly improved affinity to Fas, indicating that more rounds of computational optimization will be required to obtain peptide variants with greatly improved affinity for Fas. CD spectra in the absence and presence of Fas point to conformational changes in the peptide upon binding to Fas. Furthermore, binding of the optimized hAChE mimetic peptide to Fas could be inhibited by hAChE, providing evidence for a hAChE‐specific peptide–Fas interaction. Copyright © 2015 European Peptide Society and John Wiley & Sons, Ltd.     

Structure of the terminal PCP domain of the non‐ribosomal peptide synthetase in teicoplanin biosynthesis

The biosynthesis of the glycopeptide antibiotics, of which teicoplanin and vancomycin are representative members, relies on the combination of non‐ribosomal peptide synthesis and modification of the peptide by cytochrome P450 (Oxy) enzymes while the peptide remains bound to the peptide synthesis machinery. We have structurally characterized the final peptidyl carrier protein domain of the teicoplanin non‐ribosomal peptide synthetase machinery: this domain is believed to mediate the interactions with tailoring Oxy enzymes in addition to its function as a shuttle for intermediates between multiple non‐ribosomal peptide synthetase domains. Using solution state NMR, we have determined structures of this PCP domain in two states, the apo and the post‐translationally modified holo state, both of which conform to a four‐helix bundle assembly. The structures exhibit the same general fold as the majority of known carrier protein structures, in spite of the complex biosynthetic role that PCP domains from the final non‐ribosomal peptide synthetase module must play in glycopeptide antibiotic biosynthesis. These structures thus support the hypothesis that it is subtle rearrangements, rather than dramatic conformational changes, which govern carrier protein interactions and selectivity during non‐ribosomal peptide synthesis. Proteins 2015; 83:711–721. © 2015 Wiley Periodicals, Inc.     

Solution NMR and molecular dynamics reveal a persistent alpha helix within the dynamic region of PsbQ from photosystem II of higher plants

The extrinsic proteins of photosystem II of higher plants and green algae PsbO, PsbP, PsbQ, and PsbR are essential for stable oxygen production in the oxygen evolving center. In the available X‐ray crystallographic structure of higher plant PsbQ residues S14‐Y33 are missing. Building on the backbone NMR assignment of PsbQ, which includes this “missing link”, we report the extended resonance assignment including side chain atoms. Based on nuclear Overhauser effect spectra a high resolution solution structure of PsbQ with a backbone RMSD of 0.81 Å was obtained from torsion angle dynamics. Within the N‐terminal residues 1–45 the solution structure deviates significantly from the X‐ray crystallographic one, while the four‐helix bundle core found previously is confirmed. A short α‐helix is observed in the solution structure at the location where a β‐strand had been proposed in the earlier crystallographic study. NMR relaxation data and unrestrained molecular dynamics simulations corroborate that the N‐terminal region behaves as a flexible tail with a persistent short local helical secondary structure, while no indications of forming a β‐strand are found. Proteins 2015; 83:1677–1686. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Structural characterization of AtmS13, a putative sugar aminotransferase involved in indolocarbazole AT2433 aminopentose biosynthesis

AT2433 from Actinomadura melliaura is an indolocarbazole antitumor antibiotic structurally distinguished by its unique aminodideoxypentose‐containing disaccharide moiety. The corresponding sugar nucleotide‐based biosynthetic pathway for this unusual sugar derives from comparative genomics where AtmS13 has been suggested as the contributing sugar aminotransferase (SAT). Determination of the AtmS13 X‐ray structure at 1.50‐Å resolution reveals it as a member of the aspartate aminotransferase fold type I (AAT‐I). Structural comparisons of AtmS13 with homologous SATs that act upon similar substrates implicate potential active site residues that contribute to distinctions in sugar C5 (hexose vs. pentose) and/or sugar C2 (deoxy vs. hydroxyl) substrate specificity. Proteins 2015; 83:1547–1554. © 2015 Wiley Periodicals, Inc.     

Chiral Resolution and Absolute Configuration of the Enantiomers of the Psychoactive “Designer Drug” 3,4‐Methylenedioxypyrovalerone

Illicit rac‐MDPV (3,4‐methylenedioxypyrovalerone), manufactured in clandestine labs, has become widely abused for its cocaine‐like stimulant properties. It has recently been found as one of the toxic materials in the so‐called “bath salts,” producing, among other effects, psychosis and tachycardia in humans when introduced by any of the several routes of administration (e.g., intravenous, oral, etc.). The considerable toxicity of this “designer drug” probably resides in one of the enantiomers of the racemate. In order to obtain a sufficient amount of the enantiomers of rac‐MDPV to determine their activity, we improved the known synthesis of rac‐MDPV and found chemical resolving agents, (+)‐ and (–)‐2’‐bromotetranilic acid, that gave the MDPV enantiomers in >96% enantiomeric excess as determined by ¹H nuclear magnetic resonance and chiral high‐performance liquid chromatography. The absolute stereochemistry of these enantiomers was determined by single‐crystal X‐ray diffraction studies. Chirality 27:287‐293, 2015. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.     

Accuracy issues involved in modeling in vivo protein structures using PM7

Using the semiempirical method PM7, an attempt has been made to quantify the error in prediction of the in vivo structure of proteins relative to X‐ray structures. Three important contributory factors are the experimental limitations of X‐ray structures, the difference between the crystal and solution environments, and the errors due to PM7. The geometries of 19 proteins from the Protein Data Bank that had small R values, that is, high accuracy structures, were optimized and the resulting drop in heat of formation was calculated. Analysis of the changes showed that about 10% of this decrease in heat of formation was caused by faults in PM7, the balance being attributable to the X‐ray structure and the difference between the crystal and solution environments. A previously unknown fault in PM7 was revealed during tests to validate the geometries generated using PM7. Clashscores generated by the Molprobity molecular mechanics structure validation program showed that PM7 was predicting unrealistically close contacts between nonbonding atoms in regions where the local geometry is dominated by very weak noncovalent interactions. The origin of this fault was traced to an underestimation of the core‐core repulsion between atoms at distances smaller than the equilibrium distance. Proteins 2015; 83:1427–1435. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Solid Tryptophan as a Pseudoracemate: Physicochemical and Crystallographic Characterization

The crystalline nature of solid tryptophan has been characterized by X‐ray single crystal and powder diffraction analyses, differential scanning calorimetry, as well as measurement of solid–liquid equilibrium in water/isopropanol solution. Both the thermodynamic and crystallographic investigations have demonstrated unambiguously that solid tryptophan crystallizes in the form of a pseudoracemate (i.e., solid solution) with maximum melting over the entire enantiomeric composition range. Comparative single‐crystal X‐ray studies show that the crystal structures of racemic and enantiomeric tryptophan give very similar solid‐state packing geometries dictated by hydrogen bonding interactions. Our results indicate that the insignificant difference between homochiral and heterochiral interactions accounts for the formation of a pseudoracemate for this system. Chirality 27:88–94, 2015. © 2014 Wiley Periodicals, Inc.     

Correlation between omega and psi dihedral angles in protein structures

The planarity of the peptide group is one of the fundamental features of protein structure that is described in every chemistry and biochemistry textbook. By surveying a dataset of 163 atomic resolution protein structures we here identify the stereochemical conditions that favor significant deformations of peptide bond planarity. In particular, we demonstrate that the values of the omega dihedral angle are strictly correlated to the values of the adjacent psi angle. This trend is also observed in highly strained states such as those occurring in vicinal disulfide bridges. These findings provide direct evidence for the mutual influence of the geometrical parameters that describe the protein structure.     

Structure–function studies of Escherichia coli RnlA reveal a novel toxin structure involved in bacteriophage resistance

Escherichia coli RnlA–RnlB is a newly identified toxin–antitoxin (TA) system that plays a role in bacteriophage resistance. RnlA functions as a toxin with mRNA endoribonuclease activity and the cognate antitoxin RnlB inhibits RnlA toxicity in E. coli cells. Interestingly, T4 phage encodes the antitoxin Dmd, which acts against RnlA to promote its own propagation, suggesting that RnlA‐Dmd represents a novel TA system. Here, we have determined the crystal structure of RnlA refined to 2.10 Å. RnlA is composed of three independent domains: NTD (N ‐t erminal d omain), NRD (N r epeated d omain) and DBD (D md‐b inding d omain), which is an organization not previously observed among known toxin structures. Small‐angle X‐ray scattering (SAXS) analysis revealed that RnlA forms a dimer in solution via interactions between the DBDs from both monomers. The in vitro and in vivo functional studies showed that among the three domains, only the DBD is responsible for recognition and inhibition by Dmd and subcellular location of RnlA. In particular, the helix located at the C‐terminus of DBD plays a vital role in binding Dmd. Our comprehensive studies reveal the key region responsible for RnlA toxicity and provide novel insights into its structure–function relationship.     

Green‐lighting green fluorescent protein: Faster and more efficient folding by eliminating a cis–trans peptide isomerization event

Wild‐type green fluorescent protein (GFP) folds on a time scale of minutes. The slow step in folding is a cis–trans peptide bond isomerization. The only conserved cis‐peptide bond in the native GFP structure, at P89, was remodeled by the insertion of two residues, followed by iterative energy minimization and side chain design. The engineered GFP was synthesized and found to fold faster and more efficiently than its template protein, recovering 50% more of its fluorescence upon refolding. The slow phase of folding is faster and smaller in amplitude, and hysteresis in refolding has been eliminated. The elimination of a previously reported kinetically trapped state in refolding suggests that X‐P89 is trans in the trapped state. A 2.55 Å resolution crystal structure revealed that the new variant contains only trans‐peptide bonds, as designed. This is the first instance of a computationally remodeled fluorescent protein that folds faster and more efficiently than wild type.