Free‐energy function based on an all‐atom model for proteins

We have developed a free‐energy function based on an all‐atom model for proteins. It comprises two components, the hydration entropy (HE) and the total dehydration penalty (TDP). Upon a transition to a more compact structure, the number of accessible configurations arising from the translational displacement of water molecules in the system increases, leading to a water‐entropy gain. To fully account for this effect, the HE is calculated using a statistical‐mechanical theory applied to a molecular model for water. The TDP corresponds to the sum of the hydration energy and the protein intramolecular energy when a fully extended structure, which possesses the maximum number of hydrogen bonds with water molecules and no intramolecular hydrogen bonds, is chosen as the standard one. When a donor and an acceptor (e.g., N and O, respectively) are buried in the interior after the break of hydrogen bonds with water molecules, if they form an intramolecular hydrogen bond, no penalty is imposed. When a donor or an acceptor is buried with no intramolecular hydrogen bond formed, an energetic penalty is imposed. We examine all the donors and acceptors for backbone‐backbone, backbone‐side chain, and side chain‐side chain intramolecular hydrogen bonds and calculate the TDP. Our free‐energy function has been tested for three different decoy sets. It is better than any other physics‐based or knowledge‐based potential function in terms of the accuracy in discriminating the native fold from misfolded decoys and the achievement of high Z‐scores. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Synthesis of cyclic herpes simplex virus peptides containing 281–284 epitope of glycoprotein D‐1 in endo‐ or exo‐position

We have prepared two types of cyclopeptides containing the ²⁸¹DPVG²⁸⁴ sequence from the 276–284 region of glycoprotein gD‐1 of the Herpes simplex virus (HSV). The syntheses were performed by solid phase methodology using MBHA or BHA resin and orthogonal protection schemes. Head‐to‐side‐chain cyclization included the N‐terminal part of the epitope, while side‐chain‐to‐side‐chain lactam bridge formation resulted in a peptide containing a C‐terminal cycle. Peptides elongated by Cys at the N‐terminal of the sequence were also prepared. Boc chemistry using Fmoc and OFm orthogonal protection was applied for on‐resin cyclization. Based on the orthogonality of Bzl and cHex esters under a 1 m TMSOTf‐thioanisole/TFA cleavage condition, a new approach for the cyclization on BHA‐resin has also been developed. Preliminary studies on solution conformation of the cyclic peptides by CD spectroscopy indicated the importance of the location and the size of the cycle within the epitope sequence. Copyright © 1999 European Peptide Society and John Wiley & Sons, Ltd.     

‘Click’ chemistry synthesis and capillary electrophoresis study of 1,4‐linked 1,2,3‐triazole AZT‐systemin conjugate

The Cu(I) catalyzed Huisgen 1,3‐dipolar azide‐alkyne cycloaddition (CuAAC) was applied for a nucleoside‐peptide bioconjugation. Systemin (Sys), an 18‐aa plant signaling peptide naturally produced in response to wounding or pathogen attack, was chemically synthesized as its N‐propynoic acid functionalized analog (Prp‐Sys) using the SPPS. Next, CuAAC was applied to conjugate Prp‐Sys with 3′‐azido‐2′,3′‐dideoxythymidine (AZT), a model cargo molecule. 1,4‐Linked 1,2,3‐triazole AZT‐Sys conjugate was designed to characterize the spreading properties and ability to translocate of cargo molecules of systemin. CuAAC allowed the synthesis of the conjugate in a chemoselective and regioselective manner, with high purity and yield. The presence of Cu(I) ions generated in situ drove the CuAAC reaction to completion within a few minutes without any by‐products. Under typical separation conditions of phosphate ‘buffer’ at low pH and uncoated fused bare‐silica capillary, an increasing peak intensity assigned to triazole‐linked AZT‐Sys conjugate was observed using capillary electrophoresis (CE) during CuAAC. CE analysis showed that systemin peptides are stable in tomato leaf extract for up to a few hours. CE‐ESI‐MS revealed that the native Sys and its conjugate with AZT are translocated through the tomato stem and can be directly detected in stem exudates. The results show potential application of systemin as a transporter of low molecular weight cargo molecules in tomato plant and of CE method to characterize a behavior of plant peptides and its analogs. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

Complete structure determination of N‐acetyl‐D‐galactosamine‐binding mistletoe lectin‐3 from Viscum album L. album

The primary structure of the B chain of the N‐acetyl‐D ‐galactosamine‐recognizing mistletoe lectin‐3 (ML‐3B) has been deduced from proteolytic digest peptides of the purified glycoprotein, their HPLC‐separation and Edman degradation and confirmation of the peptide sequences by MALDI‐MS. ML‐3B consists of 262 amino acid residues including 10 cysteine moieties. The structure and linkage of the carbohydrate side chains, connected to two N‐glycosylation sites at positions Asn⁹⁵ and Asn¹³⁵ of the lectin, were determined by a combination of glycosidase treatment and MALDI‐MS of corresponding glycopeptide fragments. The sequence alignment reveals a high homology with other B chains of type‐II RIPs, although there are remarkable differences in the D ‐galactose‐specific mistletoe lectin‐1B chain. The recently published primary structure of the mistletoe lectin‐3A chain 1 and the now available primary sequence of the 3B chain allowed the construction of a preliminary homology model of ML‐3. The model demonstrates, unequivocally, that ML‐3 is a member of the type‐II RIP family with rigid conservation of the enzymatic active site of the A chain and an identical overall protein fold. Specific amino acid residue exchanges and the different glycosylation pattern in comparison with ML‐1 are discussed and related to the properties of the two glycoproteins. The knowledge of the complete primary structure of mistletoe lectin‐3 is a major contribution towards more insight into the mechanism of the biological activity of commercial mistletoe preparations. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Cryo‐EM map interpretation and protein model‐building using iterative map segmentation

A procedure for building protein chains into maps produced by single‐particle electron cryo‐microscopy (cryo‐EM) is described. The procedure is similar to the way an experienced structural biologist might analyze a map, focusing first on secondary structure elements such as helices and sheets, then varying the contour level to identify connections between these elements. Since the high density in a map typically follows the main‐chain of the protein, the main‐chain connection between secondary structure elements can often be identified as the unbranched path between them with the highest minimum value along the path. This chain‐tracing procedure is then combined with finding side‐chain positions based on the presence of density extending away from the main path of the chain, allowing generation of a C_α model. The C_α model is converted to an all‐atom model and is refined against the map. We show that this procedure is as effective as other existing methods for interpretation of cryo‐EM maps and that it is considerably faster and produces models with fewer chain breaks than our previous methods that were based on approaches developed for crystallographic maps.     

Stabilisation of a short α‐helical VIP fragment by side chain to side chain cyclisation: a comparison of common cyclisation motifs by circular dichroism

A model octapeptide segment derived from vasoactive intestinal peptide (VIP) was utilised to investigate the effect of several conventional cyclisation methods on the α‐helical conformation in short peptide fragments. Three of the classical macrocyclisation techniques (i.e. lactamisation, ring‐closing metathesis and Huisgen cycloaddition) were applied, and the conformations of the resulting cyclic peptides, as well as their linear precursors, were compared by CD analysis. The visibly higher folding propensity of the triazole‐tethered peptide after azide‐alkyne CuAAC macrocyclisation illustrates that the secondary structure of a short peptide fragment can differ significantly depending on the chemical strategy used to covalently cross‐link side chain residues in a ‘helical’ fragment. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Preferred side‐chain conformation of arginine residues in a triple‐helical structure

The crystal structure of the triple‐helical peptide (Pro‐Hyp‐Gly)₃‐Pro‐Arg‐Gly‐(Pro‐Hyp‐Gly)₄ (POG3‐PRG‐POG4) was determined at 1.45 Å resolution. POG3‐PRG‐POG4 was designed to permit investigation of the side‐chain conformation of the Arg residues in a triple‐helical structure. Because of the alternative structure of one of three Arg residues, four side‐chain conformations were observed in an asymmetric unit. Among them, three adopt a ttg^?t conformation and the other adopts a tg^?g^?t conformation. A statistical analysis of 80 Arg residues in various triple‐helical peptides showed that, unlike those in globular proteins, they preferentially adopt a tt conformation for χ₁ and χ₂, as observed in POG3‐PRG‐POG4. This conformation permits van der Waals contacts between the side‐chain atoms of Arg and the main‐chain atoms of the adjacent strand in the same molecule. Unlike many other host–guest peptides, in which there is a significant difference between the helical twists in the guest and the host peptides, POG3‐PRG‐POG4 shows a marked difference between the helical twists in the N‐terminal peptide and those in the C‐terminal peptide, separated near the Arg residue. This suggested that the unique side‐chain conformation of the Arg residue affects not only the conformation of the guest peptide, but also the conformation of the peptide away from the Arg residue. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 1000–1009, 2014.     

Insights into the helix‐sense inversion of poly(β‐phenethyl‐L‐aspartate) by two‐dimensional Raman correlation spectroscopy

In this study, the transition process of the helix‐sense inversion of poly(β‐phenethyl‐L‐aspartate) was investigated by Raman scattering and 2‐dimensional correlation spectroscopy. Temperature‐dependent Raman spectra were obtained during the helix‐sense inversion. The results of 2‐dimensional correlation analysis in the spectral regions of 1600‐1800 and 3200‐3400 cm^?1 showed that the intensity changes of the side‐chain ester C═O stretching bands occurred prior to those of amide A and amide I bands in the unwinding process of αR‐helix on heating. The sequential order of the intensity changes for amide A, amide I, and the side‐chain ester C═O stretching bands during the inversion process was determined. It was found that the conformation change of the side chain occurred prior to that of the main chain for the αR‐helix on heating. Thus, we concluded that the transformation of the backbone chain from right‐handed to left‐handed is triggered by the conformational change of the side chains.     

Protein side‐chain modeling with a protein‐dependent optimized rotamer library

Despite years of effort, the problem of predicting the conformations of protein side chains remains a subject of inquiry. This problem has three major issues, namely defining the conformations that a side chain may adopt within a protein, developing a sampling procedure for generating possible side‐chain packings, and defining a scoring function that can rank these possible packings. To solve the former of these issues, most procedures rely on a rotamer library derived from databases of known protein structures. We introduce an alternative method that is free of statistics. We begin with a rotamer library that is based only on stereochemical considerations; this rotamer library is then optimized independently for each protein under study. We show that this optimization step restores the diversity of conformations observed in native proteins. We combine this protein‐dependent rotamer library (PDRL) method with the self‐consistent mean field (SCMF) sampling approach and a physics‐based scoring function into a new side‐chain prediction method, SCMF–PDRL. Using two large test sets of 831 and 378 proteins, respectively, we show that this new method compares favorably with competing methods such as SCAP, OPUS‐Rota, and SCWRL4 for energy‐minimized structures. Proteins 2014; 82:2000–2017. © 2014 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.     

Solvent‐dependent conformation of a regioselective amylose carbamate: Amylose‐2‐acetyl‐3,6‐bis(phenylcarbamate)

Six amylose‐2‐acetyl‐3,6‐bis(phenylcarbamate) (AAPC) samples ranging in weight‐average molar mass M_w from 1.8 × 10⁴ g mol^?1 to 1.1 × 10⁶ g mol^?1 have been prepared from enzymatically synthesized amylose samples. Static light scattering, small‐angle X‐ray scattering, sedimentation equilibrium, and viscosity measurements were made for the samples in 1,4‐dioxane (DIOX), 2‐ethoxyethanol (2EE), and 2‐butanone (MEK) all at 25°C to determine particle scattering functions, z‐average radii of gyration, intrinsic viscosities, as well as M_w. The data were analyzed in terms of the wormlike cylinder model mainly to yield the helix pitch per residue h and the Kuhn segment length λ^?1, which corresponds to twice of the persistence length. The latter parameters (λ^?1) in 2EE (11 nm) and MEK (12 nm) are quite smaller than those for amylose tris(phenylcarbamate) (ATPC) in the same solvent (16 nm in 2EE and 18 nm in MEK) whereas those for AAPC (21 nm) and ATPC (22 nm) in DIOX are essentially the same as each other. This indicates that the chain stiffness of AAPC is more strongly influenced by the solvents since the number of intramolecular H‐bonds of AAPC is more changeable than that for ATPC. © 2012 Wiley Periodicals, Inc. Biopolymers 97:1010–1017, 2012.     

Assessment of protein side‐chain conformation prediction methods in different residue environments

Computational prediction of side‐chain conformation is an important component of protein structure prediction. Accurate side‐chain prediction is crucial for practical applications of protein structure models that need atomic‐detailed resolution such as protein and ligand design. We evaluated the accuracy of eight side‐chain prediction methods in reproducing the side‐chain conformations of experimentally solved structures deposited to the Protein Data Bank. Prediction accuracy was evaluated for a total of four different structural environments (buried, surface, interface, and membrane‐spanning) in three different protein types (monomeric, multimeric, and membrane). Overall, the highest accuracy was observed for buried residues in monomeric and multimeric proteins. Notably, side‐chains at protein interfaces and membrane‐spanning regions were better predicted than surface residues even though the methods did not all use multimeric and membrane proteins for training. Thus, we conclude that the current methods are as practically useful for modeling protein docking interfaces and membrane‐spanning regions as for modeling monomers. Proteins 2014; 82:1971–1984. © 2014 Wiley Periodicals, Inc.     

A three‐dimensional model of mammalian tyrosinase active site accounting for loss of function mutations

Tyrosinases are the first and rate‐limiting enzymes in the synthesis of melanin pigments responsible for colouring hair, skin and eyes. Mutation of tyrosinases often decreases melanin production resulting in albinism, but the effects are not always understood at the molecular level. Homology modelling of mouse tyrosinase based on recently published crystal structures of non‐mammalian tyrosinases provides an active site model accounting for loss‐of‐function mutations. According to the model, the copper‐binding histidines are located in a helix bundle comprising four densely packed helices. A loop containing residues M374, S375 and V377 connects the CuA and CuB centres, with the peptide oxygens of M374 and V377 serving as hydrogen acceptors for the NH‐groups of the imidazole rings of the copper‐binding His367 and His180. Therefore, this loop is essential for the stability of the active site architecture. A double substitution ³⁷⁴MS³⁷⁵→³⁷⁴GG³⁷⁵ or a single M374G mutation lead to a local perturbation of the protein matrix at the active site affecting the orientation of the H367 side chain, that may be unable to bind CuB reliably, resulting in loss of activity. The model also accounts for loss of function in two naturally occurring albino mutations, S380P and V393F. The hydroxyl group in S380 contributes to the correct orientation of M374, and the substitution of V393 for a bulkier phenylalanine sterically impedes correct side chain packing at the active site. Therefore, our model explains the mechanistic necessity for conservation of not only active site histidines but also adjacent amino acids in tyrosinase.     

Syntheses and Antitumor Evaluation of C(6)‐Isobutyl‐ and C(6)‐Isobutenyl‐Substituted Pyrimidines,and Dihydropyrrolo[1,2‐c]pyrimidine‐1,3‐diones

A growing body of evidence supports that pyrimidine derivatives, in which the sugar residues have been replaced by acyclic side chains, might be developed as promising anticancer agents that interfere with tumor cell proliferation, survival, and metastatic formation. In this work, we prepared novel pyrimidines bearing i‐Bu (i.e., 3, 4 , and 7 – 9 ) and isobutenyl (i.e., 5 and 10 ) side chains at C(6) and examined their in vitro effects on tumor cell lines. The dihydropyrrolo[1,2‐c]pyrimidine‐1,3‐diones 6 and 11 were obtained as products of intramolecular cyclization, which occurred during the removal of Bn in 5 or MeO protecting groups in 10 . Fluorination of 3 with diethylaminosulfur trifluoride (DAST) and then dehydrohalogenation of the resulting fluorinated derivative 4 afforded 6‐isobut‐2′‐enyl pyrimidine derivative 5 with a C(2′)C(3′) bond. For the preparation of 6‐isobut‐1′‐en‐1‐yl pyrimidine 10 , a synthetic strategy involving acetylation of the 1,3‐diols was applied. Antitumor evaluation of compounds 3 – 11 showed that 2,4‐dimethoxypyrimidine containing 6‐[(1,3‐dibenzyloxy)‐2‐hydroxy]methyl side chain, 3 , exerted a strong antiproliferative effect on the studied tumor cell lines. Additionally, it was shown that the mechanism of antiproliferative effect of 3 in HeLa cells include early G2/M arrest and apoptosis, as well as a p53‐independent S‐phase arrest upon prolonged treatment.     

OPUS‐SSF: A side‐chain‐inclusive scoring function for ranking protein structural models

We introduce a side‐chain‐inclusive scoring function, named OPUS‐SSF, for ranking protein structural models. The method builds a scoring function based on the native distributions of the coordinate components of certain anchoring points in a local molecular system for peptide segments of 5, 7, 9, and 11 residues in length. Differing from our previous OPUS‐CSF [Xu et al., Protein Sci. 2018; 27: 286–292], which exclusively uses main chain information, OPUS‐SSF employs anchoring points on side chains so that the effect of side chains is taken into account. The performance of OPUS‐SSF was tested on 15 decoy sets containing totally 603 proteins, and 571 of them had their native structures recognized from their decoys. Similar to OPUS‐CSF, OPUS‐SSF does not employ the Boltzmann formula in constructing scoring functions. The results indicate that OPUS‐SSF has achieved a significant improvement on decoy recognition and it should be a very useful tool for protein structural prediction and modeling.     

Syntheses and activities of backbone‐side chain cyclic octapeptide ligands with N‐functionalized phosphotyrosine for the N‐terminal SH2‐domain of the protein tyrosine phosphatase SHP‐1

The protein tyrosine phosphatase SHP‐1 plays an important role in many physiological and pathophysiological processes. This phosphatase is activated through binding of ligands to its SH2‐domains, mainly to the N‐terminal one. Based on a theoretical docking model, backbone‐to‐side chain cyclized octapeptides were designed as ligands. Assembly of such modelled structures required the synthesis of N‐functionalized tyrosine derivatives and their incorporation into the sequence. Because of difficulties encountered in the condensation of N‐protected amino acids to the N‐alkylated tyrosine‐peptide we synthesized and used preformed dipeptide building units. As all attempts to obtain phosphorylated dipeptide units failed, the syntheses had to be performed with a free phenolic function. Use of different N‐alkyl or cycloalkyl residues in the N‐functionalized side chains allowed to investigate the effect of ring size, flexibility and hydrophobicity of formed lactam bridges on stimulatory activity. All tested linear and cyclic octapeptides stimulate the phosphatase activity of SHP‐1. Stimulatory activities of cyclic ligands increase with the chain length of the lactam bridges resulting in increased flexibility and better entropic preformation of the binding conformation. The strong activity of some cyclic octapeptides supports the modelled structure. Copyright © 2010 European Peptide Society and John Wiley & Sons, Ltd.     

Prediction of protein mutation effects based on dehydration and hydrogen bonding – A large‐scale study

Reliable computational prediction of protein side chain conformations and the energetic impact of amino acid mutations are the key aspects for the optimization of biotechnologically relevant enzymatic reactions using structure‐based design. By improving the protein stability, higher yields can be achieved. In addition, tuning the substrate selectivity of an enzymatic reaction by directed mutagenesis can lead to higher turnover rates. This work presents a novel approach to predict the conformation of a side chain mutation along with the energetic effect on the protein structure. The HYDE scoring concept applied here describes the molecular interactions primarily by evaluating the effect of dehydration and hydrogen bonding on molecular structures in aqueous solution. Here, we evaluate its capability of side‐chain conformation prediction in classic remutation experiments. Furthermore, we present a new data set for evaluating “cross‐mutations,” a new experiment that resembles real‐world application scenarios more closely. This data set consists of protein pairs with up to five point mutations. Thus, structural changes are attributed to point mutations only. In the cross‐mutation experiment, the original protein structure is mutated with the aim to predict the structure of the side chain as in the paired mutated structure. The comparison of side chain conformation prediction (“remutation”) showed that the performance of HYDE_protein is qualitatively comparable to state‐of‐the art methods. The ability of HYDE_protein to predict the energetic effect of a mutation is evaluated in the third experiment. Herein, the effect on protein stability is predicted correctly in 70% of the evaluated cases. Proteins 2017; 85:1550–1566. © 2017 Wiley Periodicals, Inc.     

All‐atom chain‐building by optimizing MODELLER energy function using conformational space annealing

We have investigated the effect of rigorous optimization of the MODELLER energy function for possible improvement in protein all‐atom chain‐building. For this we applied the global optimization method called conformational space annealing (CSA) to the standard MODELLER procedure to achieve better energy optimization than what MODELLER provides. The method, which we call MODELLERCSA , is tested on two benchmark sets. The first is the 298 proteins taken from the HOMSTRAD multiple alignment set. By simply optimizing the MODELLER energy function, we observe significant improvement in side‐chain modeling, where MODELLERCSA provides about 10.7% (14.5%) improvement for χ₁ (χ₁ + χ₂) accuracy compared to the standard MODELLER modeling. The improvement of backbone accuracy by MODELLERCSA is shown to be less prominent, and a similar improvement can be achieved by simply generating many standard MODELLER models and selecting lowest energy models. However, the level of side‐chain modeling accuracy by MODELLERCSA could not be matched either by extensive MODELLER strategies, side‐chain remodeling by SCWRL3, or copying unmutated rotamers. The identical procedure was successfully applied to 100 CASP7 template base modeling domains during the prediction season in a blind fashion, and the results are included here for comparison. From this study, we observe a good correlation between the MODELLER energy and the side‐chain accuracy. Our findings indicate that, when a good alignment between a target protein and its templates is provided, thorough optimization of the MODELLER energy function leads to accurate all‐atom models. Proteins 2009. © 2008 Wiley‐Liss, Inc.